title: Combining unsupervised and supervised learning for predicting the final stroke lesion

abstract: Predicting the final ischaemic stroke lesion provides crucial information regarding the volume of salvageable hypoperfused tissue, which helps physicians in the difficult decision-making process of treatment planning and intervention. Treatment selection is influenced by clinical diagnosis, which requires delineating the stroke lesion, as well as characterising cerebral blood flow dynamics using neuroimaging acquisitions. Nonetheless, predicting the final stroke lesion is an intricate task, due to the variability in lesion size, shape, location and the underlying cerebral haemodynamic processes that occur after the ischaemic stroke takes place. Moreover, since elapsed time between stroke and treatment is related to the loss of brain tissue, assessing and predicting the final stroke lesion needs to be performed in a short period of time, which makes the task even more complex. Therefore, there is a need for automatic methods that predict the final stroke lesion and support physicians in the treatment decision process. We propose a fully automatic deep learning method based on unsupervised and supervised learning to predict the final stroke lesion after 90 days. Our aim is to predict the final stroke lesion location and extent, taking into account the underlying cerebral blood flow dynamics that can influence the prediction. To achieve this, we propose a two-branch Restricted Boltzmann Machine, which provides specialized data-driven features from different sets of standard parametric Magnetic Resonance Imaging maps. These data-driven feature maps are then combined with the parametric Magnetic Resonance Imaging maps, and fed to a Convolutional and Recurrent Neural Network architecture. We evaluated our proposal on the publicly available ISLES 2017 testing dataset, reaching a Dice score of 0.38, Hausdorff Distance of 29.21 mm, and Average Symmetric Surface Distance of 5.52 mm.


title: Multi-stage Deep Layer Aggregation for Brain Tumor Segmentation

abstract: Gliomas are among the most aggressive and deadly brain tumors. This paper details the proposed Deep Neural Network architecture for brain tumor segmentation from Magnetic Resonance Images. The architecture consists of a cascade of three Deep Layer Aggregation neural networks, where each stage elaborates the response using the feature maps and the probabilities of the previous stage, and the MRI channels as inputs. The neuroimaging data are part of the publicly available Brain Tumor Segmentation (BraTS) 2020 challenge dataset, where we evaluated our proposal in the BraTS 2020 Validation and Test sets. In the Test set, the experimental results achieved a Dice score of 0.8858, 0.8297 and 0.7900, with an Hausdorff Distance of 5.32 mm, 22.32 mm and 20.44 mm for the whole tumor, core tumor and enhanced tumor, respectively.


title: A New Framework for Inference on Markov Population Models

abstract: In this work we construct a joint Gaussian likelihood for approximate inference on Markov population models. We demonstrate that Markov population models can be approximated by a system of linear stochastic differential equations with time-varying coefficients. We show that the system of stochastic differential equations converges to a set of ordinary differential equations. We derive our proposed joint Gaussian deterministic limiting approximation (JGDLA) model from the limiting system of ordinary differential equations. The results is a method for inference on Markov population models that relies solely on the solution to a system deterministic equations. We show that our method requires no stochastic infill and exhibits improved predictive power in comparison to the Euler-Maruyama scheme on simulated susceptible-infected-recovered data sets. We use the JGDLA to fit a stochastic susceptible-exposed-infected-recovered system to the Princess Diamond COVID-19 cruise ship data set.


title: Low temperature Universalities in amorphous systems: role of microscopic length scales

abstract: We find that a competition between dispersion forces among molecules in solids and their phonon mediated coupling leads to a natural length scale based on molecular parameters and relevant to decipher glass anomalies. For amorphous systems, the length scale is of the medium range topological orders and its ratio with the distance of closest approach between molecules turns out to be a constant. This in turn leads to a material independent, constant value of the ratio ${\gamma_l \over \gamma_t}$, with $\gamma_l$ and $\gamma_t$ as the coupling-strength for two amorphous molecules mediated by longitudinal and transverse phonons (also referred as Meissner-Berret ratio) and thereby provides a theoretical explanation of their experimentally obserevd quantitative universality in \cite{mb}. The above length scales are also related to Ioffe-Regel frequency and boson peak frequency of the vibrational spectrum and indicate that the former is of the same order as the latter for transverse phonon-dynamics.


title: Brown-York Energy in Spacetimes with Horizon

abstract: We obtain a simple relationship between the change in the Brown-York energy inside of a closed two-surface just outside of a horizon of spacetime, and the change in the area of that two-surface.


title: A Provably Efficient Algorithm for Linear Markov Decision Process with Low Switching Cost

abstract: Many real-world applications, such as those in medical domains, recommendation systems, etc, can be formulated as large state space reinforcement learning problems with only a small budget of the number of policy changes, i.e., low switching cost. This paper focuses on the linear Markov Decision Process (MDP) recently studied in [Yang et al 2019, Jin et al 2020] where the linear function approximation is used for generalization on the large state space. We present the first algorithm for linear MDP with a low switching cost. Our algorithm achieves an $\widetilde{O}\left(\sqrt{d^3H^4K}\right)$ regret bound with a near-optimal $O\left(d H\log K\right)$ global switching cost where $d$ is the feature dimension, $H$ is the planning horizon and $K$ is the number of episodes the agent plays. Our regret bound matches the best existing polynomial algorithm by [Jin et al 2020] and our switching cost is exponentially smaller than theirs. When specialized to tabular MDP, our switching cost bound improves those in [Bai et al 2019, Zhang et al 20020]. We complement our positive result with an $\Omega\left(dH/\log d\right)$ global switching cost lower bound for any no-regret algorithm.


title: Volterra model-based control for nonlinear systems via Carleman linearization

abstract: This paper presents detailed insights of embedding Carleman linearization into nonlinear systems for designing Volterra model-based control technique. Volterra series method is a competent mathematical tool, which extends the convolution integral for linear systems to nonlinear systems. First, we utilize the Carleman linearization technique to arrive at the bilinear approximation of the nonlinear system. Secondly, the third-order Volterra model representation is computed from the Carleman bilinearized model. Then, Volterra model-based control strategy is developed. The proposed method is effectuated on the benchmark van de Vusse reactor exhibiting non-minimum phase response. The open and closed-loop simulation results are presented, demonstrating the superiority and practical utility of the proposed method.


title: Smart Car Features using Embedded Systems and IoT

abstract: There has been a tremendous rise in technological advances in the field of automobiles and autonomous vehicles. With the increase in the number of driven vehicles, the safety concerns with the same have also risen. The cases of accidents and life-threatening injuries have skyrocketed. It has become a necessity to provide adequate safety measures in automobiles. This project aims to develop a prototype for a smart vehicle system that provides real-time location of the vehicle on detection of a crash and alert the police station and relatives of the user, it has a panic button feature for a passenger's safety. We also demonstrate a mechanism for cabin monitoring and an interactive interface between a user and a car, where the user can inquire about the temperature, humidity, and other variables inside the car remotely by sending a text message to the GSM module which is present in the car. The GSM module connects to the Arduino, which fetches the readings from sensors attached to it and sends it back to the user through a text message. We show the integration of MQ3 Alcohol sensor with Arduino for drunk driving prevention.


title: Generalizing the Swampland: Embedding $P(X, \varphi)$ Inflationary Theories in a Curved Multi-field Space

abstract: We study the general embedding of a $ P(X, \varphi) $ inflationary theory into a two-field theory with curved field space metric, which was proposed as a possible way to examine the relation between de Sitter Swampland conjecture and \textit{k}-inflation. We show that this embedding method fits into the special type of two-field model in which the heavy field can be integrated out at the full action level. However, this embedding is not exact due to the upper bound of the effective mass of the heavy field. We quantify the deviation between the speed of sound calculated via the $ P(X, \varphi) $ theory and the embedding two-field picture to next leading order terms. We especially focus on the first potential slow roll parameter defined in the two-field picture and obtain an upper bound on it.


title: Nonstandard Expansiveness

abstract: Let $(X,d)$ be a metric space and $f: X \rightarrow X$ be a homeomorphism. We say that a dynamical system $(X,f)$ is \emph{expansive}, with constant of expansivity $c \in \mathbb{R{^+}}$, if for all $x,y \in X$ , $x \neq y$, exists $n \in \mathbb{Z}$, such that $d(f^n(x), f^n(y)) >c$. In this paper we will use the theory of Nonstandard Analysis to study a subfamily of these dynamics, which verify that for all $x,y \in X$, if $x\neq y$ then the set $\lbrace n \in \mathbb{Z} : d(f^n(x), f^n(y) > c \rbrace$ is infinite.


title: A Molecular Dynamics Study of Laser-Excited Gold

abstract: The structural evolution of laser-excited systems of gold has previously been measured through ultrafast MeV electron diffraction. However, there has been a long-standing inability of atomistic simulations to provide a consistent picture of the melt process, concluding in large discrepancies between the predicted threshold energy density for complete melt, as well as the transition between heterogeneous and homogeneous melting. We make use of two-temperature classical molecular dynamics simulations utilizing three highly successful interatomic potentials and reproduce electron diffraction data presented by Mo et al. We recreate the experimental electron diffraction data employing both a constant and temperature-dependent electron-ion equilibration rate. In all cases we are able to match time-resolved electron diffraction data, and find consistency between atomistic simulations and experiments, only by allowing laser energy to be transported away from the interaction region. This additional energy-loss pathway, which scales strongly with laser fluence, we attribute to hot electrons leaving the target on a timescale commensurate with melting.


title: Signal anticipation and delay in excitable media: group delay of the FitzHugh-Nagumo model

abstract: An expression for the group delay of the FitzHugh-Nagumo model in response to low amplitude input is obtained by linearisation of the cubic term of the voltage equation around its stable fixed-point. It is found that a negative group delay exists for low frequencies, indicating that the evolution of slowly fluctuating signals are anticipated by the voltage dynamics. The effects of the group delay for different types of signals are shown numerically for the non-linearised FitzHugh-Nagumo model, and some observations on the signal aspects that are anticipated are stated.


title: Birational geometry of sextic double solids with a compound $A_n$ singularity

abstract: Sextic double solids, double covers of $\mathbb P^3$ branched along a sextic surface, are the lowest degree Gorenstein Fano 3-folds, hence are expected to behave very rigidly in terms of birational geometry. Smooth sextic double solids, and those which are $\mathbb Q$-factorial with ordinary double points, are known to be birationally rigid. In this article, we study sextic double solids with an isolated compound $A_n$ singularity. We prove a sharp bound $n \leq 8$, describe models for each $n$ explicitly and prove that sextic double solids with $n > 3$ are birationally non-rigid.


title: Optical variability power spectrum analysis of blazar sources on intranight timescales

abstract: We report the first results of a systematic investigation to characterize blazar variability power spectral densities (PSDs) at optical frequencies using densely sampled (5--15 minutes integration time), high photometric accuracy ($\lesssim$0.2--0.5\%) R-band intranight light curves, covering timescales ranging from several hours to $\sim$15\,minutes. Our sample consists of 14 optically bright blazars, including nine BL Lacertae objects (BL Lacs) and five flat-spectrum radio quasars (FSRQs) which have shown statistically significant variability during 29 monitoring sessions. We model the intranight PSDs as simple power--laws and derive the best-fit slope along with uncertainty using the `power spectral response' method. Our main results are the following: (1) on 19 out of 29 monitoring sessions, the intranight PSDs show an acceptable fit to simple power-laws at the rejection confidence $\leq$ 90\%; (2) for these 19 instances, the PSD slopes show a large range from 1.4 to 4.0, consistent with statistical characters of red (slope$\sim$2) and black (slope$\geq$3) noise stochastic processes; (3) the average PSD slopes for the BL Lacs and FSRQs are indistinguishable from one another; (4) the normalization of intranight PSDs for individual blazar sources which were monitored on more than one occasion turns out to be consistent with one another with a few exceptions. The average PSD slope, 2.9$\pm$0.3 (1$\sigma$ uncertainty) is steeper than the red-noise type character of variability found on longer timescales (many decades to days), indicative of a cutoff in the variability spectrum on timescales around a few days at the synchrotron frequencies of the emission spectrum.


title: Modularity maximisation for graphons

abstract: Networks are a widely-used tool to investigate the large-scale connectivity structure in complex systems and graphons have been proposed as an infinite size limit of dense networks. The detection of communities or other meso-scale structures is a prominent topic in network science as it allows the identification of functional building blocks in complex systems. When such building blocks may be present in graphons is an open question. In this paper, we define a graphon-modularity and demonstrate that it can be maximised to detect communities in graphons. We then investigate specific synthetic graphons and show that they may show a wide range of different community structures. We also reformulate the graphon-modularity maximisation as a continuous optimisation problem and so prove the optimal community structure or lack thereof for some graphons, something that is usually not possible for networks. Furthermore, we demonstrate that estimating a graphon from network data as an intermediate step can improve the detection of communities, in comparison with exclusively maximising the modularity of the network. While the choice of graphon-estimator may strongly influence the accord between the community structure of a network and its estimated graphon, we find that there is a substantial overlap if an appropriate estimator is used. Our study demonstrates that community detection for graphons is possible and may serve as a privacy-preserving way to cluster network data.

categories: stat.CO cs.LG cs.SI nlin.AO physics.soc-ph


title: Acoustic radiation-free surface phononic crystal resonator for in-liquid low-noise gravimetric detection

abstract: Acoustic wave resonators are promising for gravimetric biosensing. However, they generally suffer from strong acoustic radiation in liquid, which limits their quality factor and increases their frequency noise. This article presents an acoustic-radiation-free gravimetric biosensor based on a locally-resonant surface phononic crystal (SPC) consisting of periodic high aspect ratio electrodes to ad-dress the above issue. The acoustic wave generated in the SPC is slower than the sound wave in water, hence preventing acoustic propagation in the fluid and resulting in energy confinement near the electrode surface. This energy confinement results in a significant quality factor improvement and thus reduces the frequency noise. The proposed SPC resonator is numerically studied by finite element analysis and experimentally implemented by an electroplating based fabrication process. Experimental results show that the SPC resonator exhibits an in-liquid quality factor 15 times higher than a conventional Rayleigh wave resonator with a similar operating frequency. The proposed radiation suppression method using SPC can also be applied in other types of acoustic wave resonators. It can thus serve as a general technique for boosting the in-liquid quality factor and the sensing performance of many acoustic biosensors.


title: Compressible fluids interacting with plates -- regularity and weak-strong uniqueness

abstract: In this paper, we study a nonlinear interaction problem between compressible viscous fluids and plates. For this problem, we introduce relative entropy and relative energy inequality for the finite energy weak solutions (FEWS). First, we prove that for all FEWS, the structure displacement enjoys improved regularity by utilizing the dissipation effects of the fluid onto the structure and that all FEWS satisfy the relative energy inequality. Then, we show that all FEWS enjoy the weak-strong uniqueness property, thus extending the classical result for compressible Navier-Stokes system to this fluid-structure interaction problem.


title: Influence of NLTE effects in Fe I lines on inverted atmosphere II. 6301 A and 6302 A lines formed in 3DNLTE

abstract: This paper forms the second part of our study on how the neglect of NLTE conditions in the formation of Fe I 6301.5 A and the 6302.5 A lines influences the atmosphere obtained by inverting their profiles in LTE. The main cause of NLTE effects is the line opacity deficit due to the excess ionization of the Fe I atoms by the UV photons in the Sun. In the first paper, the above photospheric lines were assumed to have formed in 1DNLTE and the effects of horizontal radiation transfer (RT) were neglected. In the present paper, the iron lines are computed in 3DNLTE. We investigate the influence of horizontal RT on the inverted atmosphere and how it can enhance or reduce the errors due to the neglect of 1DNLTE effects. The iron lines are computed in LTE, 1DNLTE and 3DNLTE. They all are inverted using an LTE inversion code. The atmosphere from the inversion of LTE profiles is taken as the reference model. The test atmospheres from the inversion of 1DNLTE and 3DNLTE profiles are compared with it. The differences between models are analysed and correspondingly attributed to NLTE and 3D effects. The effects of horizontal RT are evident in regions surrounded by strong horizontal gradients in temperature. In some regions, the 3D effects enhance the 1DNLTE effects while in some, they weaken. The errors due to neglecting the 3D effects are less than 5% in temperature while the errors are mostly less than 20% in both velocity and magnetic field strength. These errors are found to survive spatial and spectral degradation. The neglect of horizontal RT is found to introduce errors in the derived atmosphere. How large the errors are depends on how strong the local horizontal gradients are in temperature. Compared to the 1DNLTE effect, the 3D effects are more localised to specific regions in the atmosphere and overall less dominant.


title: Minimizing the number of complete bipartite graphs in a $K_s$-saturated graph

abstract: A graph $G$ is $F$-saturated if it contains no copy of $F$ as a subgraph but the addition of any new edge to $G$ creates a copy of $F$. We prove that for $s \geq 3$ and $t \geq 2$, the minimum number of copies of $K_{1,t}$ in a $K_s$-saturated graph is $\Theta ( n^{t/2})$. More precise results are obtained when $t = 2$ where the problem is related to Moore graphs with diameter 2 and girth 5. We prove that for $s \geq 4$ and $t \geq 3$, the minimum number of copies of $K_{2,t}$ in an $n$-vertex $K_s$-saturated graph is at least $\Omega( n^{t/5 + 8/5})$ and at most $O(n^{t/2 + 3/2})$. These results answer a question of Chakraborti and Loh. General estimates on the number of copies of $K_{a,b}$ in a $K_s$-saturated graph are also obtained, but finding an asymptotic formula remains open.


title: Clark Measures for Rational Inner Functions

abstract: We analyze the fine structure of Clark measures and Clark isometries associated with two-variable rational inner functions on the bidisk. In the degree (n,1) case, we give a complete description of supports and weights for both generic and exceptional Clark measures, characterize when the associated embedding operators are unitary, and give a formula for those embedding operators. We also highlight connections between our results and both the structure of Agler decompositions and study of extreme points for the set of positive pluriharmonic measures on 2-torus.


title: Regularization-based Continual Learning for Anomaly Detection in Discrete Manufacturing

abstract: The early and robust detection of anomalies occurring in discrete manufacturing processes allows operators to prevent harm, e.g. defects in production machinery or products. While current approaches for data-driven anomaly detection provide good results on the exact processes they were trained on, they often lack the ability to flexibly adapt to changes, e.g. in products. Continual learning promises such flexibility, allowing for an automatic adaption of previously learnt knowledge to new tasks. Therefore, this article discusses different continual learning approaches from the group of regularization strategies, which are implemented, evaluated and compared based on a real industrial metal forming dataset.


title: Linear Properties of Electrostatic Wave in Two-Component Fermi Plasma

abstract: We study the two component Fermi plasma. Two components are electrons and ions. Using the Quantum-Hydrodynamic model (QHD), we study the linear properties of electrostatic wave. We derive the linear dispersion relation for the system from dynamical governing equations for the system. We study the dependence of linear-dispersion relation on various parameters of the system.


title: A new approach to achieving high granularity for silicon diode detectors with impact ionization gain

abstract: Low Gain Avalanche Diodes (LGADs) are thin (20-50 $\mu m$)silicon di ode sensors with modest internal gain (typically 5 to 50) and exceptional time resolution (17 $ps$ to 50 $ps$). However, the granularity of such devices is limited to the millimeter scale due to the need to include protection structures at the boundaries of the readout pads to avoid premature breakdown due to large local electric fields. In this paper we present a new approach -- the Deep-Junction LGAD (DJ-LGAD) -- that decouples the high-field gain region from the readout plane. This approach is expected to improve the achievable LGAD granularity to the tens-of-micron scale while maintaining direct charge collection on the segmented electrodes.


title: Fast solution of fully implicit Runge-Kutta and discontinuous Galerkin in time for numerical PDEs, Part I: the linear setting

abstract: Fully implicit Runge-Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). This framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time stepping schemes, and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite-element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally implicit RK methods.


title: On sequences of homomorphisms into measure algebras and the Efimov Problem

abstract: For given Boolean algebras $\mathbb{A}$ and $\mathbb{B}$ we endow the space $\mathcal{H}(\mathbb{A},\mathbb{B})$ of all Boolean homomorphisms from $\mathbb{A}$ to $\mathbb{B}$ with various topologies and study convergence properties of sequences in $\mathcal{H}(\mathbb{A},\mathbb{B})$. We are in particular interested in the situation when $\mathbb{B}$ is a measure algebra as in this case we obtain a natural tool for studying topological convergence properties of sequences of ultrafilters on $\mathbb{A}$ in random extensions of the set-theoretical universe. This appears to have strong connections with Dow and Fremlin's result stating that there are Efimov spaces in the random model. We also investigate relations between topologies on $\mathcal{H}(\mathbb{A},\mathbb{B})$ for a Boolean algebra $\mathbb{B}$ carrying a strictly positive measure and convergence properties of sequences of measures on $\mathbb{A}$.


title: Envelopes for multivariate linear regression with linearly constrained coefficients

abstract: A constrained multivariate linear model is a multivariate linear model with the columns of its coefficient matrix constrained to lie in a known subspace. This class of models includes those typically used to study growth curves and longitudinal data. Envelope methods have been proposed to improve estimation efficiency in the class of unconstrained multivariate linear models, but have not yet been developed for constrained models that we develop in this article. We first compare the standard envelope estimator based on an unconstrained multivariate model with the standard estimator arising from a constrained multivariate model in terms of bias and efficiency. Then, to further improve efficiency, we propose a novel envelope estimator based on a constrained multivariate model. Novel envelope-based testing methods are also proposed. We provide support for our proposals by simulations and by studying the classical dental data and data from the China Health and Nutrition Survey and a study of probiotic capacity to reduced Salmonella infection.


title: A General Deep Reinforcement Learning Framework for Grant-Free NOMA Optimization in mURLLC

abstract: Grant-free non-orthogonal multiple access (GF-NOMA) is a potential technique to support massive Ultra-Reliable and Low-Latency Communication (mURLLC) service. However, the dynamic resource configuration in GF-NOMA systems is challenging due to random traffics and collisions, that are unknown at the base station (BS). Meanwhile, joint consideration of the latency and reliability requirements makes the resource configuration of GF-NOMA for mURLLC more complex. To address this problem, we develop a general learning framework for signature-based GF-NOMA in mURLLC service taking into account the multiple access signature collision, the UE detection, as well as the data decoding procedures for the K-repetition GF and the Proactive GF schemes. The goal of our learning framework is to maximize the long-term average number of successfully served users (UEs) under the latency constraint. We first perform a real-time repetition value configuration based on a double deep Q-Network (DDQN) and then propose a Cooperative Multi-Agent learning technique based on the DQN (CMA-DQN) to optimize the configuration of both the repetition values and the contention-transmission unit (CTU) numbers. Our results show that the number of successfully served UEs under the same latency constraint in our proposed learning framework is up to ten times for the K-repetition scheme, and two times for the Proactive scheme, more than that with fixed repetition values and CTU numbers. In addition, the superior performance of CMA-DQN over the conventional load estimation-based approach (LE-URC) demonstrates its capability in dynamically configuring in long term. Importantly, our general learning framework can be used to optimize the resource configuration problems in all the signature-based GF-NOMA schemes.


title: Some properties of q-Gaussian distributions

abstract: The q-Gaussian is a probability distribution generalizing the Gaussian one. In spite of a q-normal distribution is popular, there is a problem when calculating an expectation value with a corresponding normalized distribution and not a q-normal distribution itself. In this paper, two q-moments types called normalized and unormalized q-moments are introduced in details. Some properties of q-moments are given, and several relationships between them are established, and some results related to q-moments are also obtained. Moreover, we show that these new q-moments may be regarded as a generelazation of the classical case for q = 1. Firstly, we determine the q-moments of q-Gaussian distribution. Especially, we give explicitely the kurtosis parameters. Secondly, we compute the expression of the q-Laplace transform of the q-Gaussian distribution. Finally, we study the distribution of sum of q-independent Gaussian distributions.


title: Advanced and comprehensive research on the dynamics of COVID-19 under mass communication outlets intervention and quarantine strategy: a deterministic and probabilistic approach

abstract: The ongoing Coronavirus disease 2019 (COVID-19) is a major crisis that has significantly affected the healthcare sector and global economies, which made it the main subject of various fields in scientific and technical research. To properly understand and control this new epidemic, mathematical modelling is presented as a very effective tool that can illustrate the mechanisms of its propagation. In this regard, the use of compartmental models is the most prominent approach adopted in the literature to describe the dynamics of COVID-19. Along the same line, we aim during this study to generalize and ameliorate many existing works that consecrated to analyse the behaviour of this epidemic. Precisely, we propose an SQEAIHR epidemic system for Coronavirus. Our constructed model is enriched by taking into account the media intervention and vital dynamics. By the use of the next-generation matrix method, the theoretical basic reproductive number $R_0$ is obtained for COVID-19. Based on some nonstandard and generalized analytical techniques, the local and global stability of the disease-free equilibrium are proven when $R_0 < 1$. Moreover, in the case of $R_0 > 1$, the uniform persistence of COVID-19 model is also shown. In order to better adapt our epidemic model to reality, the randomness factor is taken into account by considering a proportional white noises, which leads to a well-posed stochastic model. Under appropriate conditions, interesting asymptotic properties are proved, namely: extinction and persistence in the mean. The theoretical results show that the dynamics of the perturbed COVID-19 model are determined by parameters that are closely related to the magnitude of the stochastic noise. Finally, we present some numerical illustrations to confirm our theoretical results and to show the impact of media intervention and quarantine strategies.


title: Spin evolution of neutron stars in two modes: implication for millisecond pulsars

abstract: An understanding of spin frequency ($\nu$) evolution of neutron stars in the low-mass X-ray binary (LMXB) phase is essential to explain the observed $\nu$-distribution of millisecond pulsars (MSPs), and to probe the stellar and binary physics, including the possibility of continuous gravitational wave emission. Here, using numerical computations we conclude that $\nu$ can evolve in two distinctly different modes, as $\nu$ may approach a lower spin equilibrium value ($\nu_{\rm eq,per}$) for persistent accretion for a long-term average accretion rate ($\dot{M}_{\rm av}$) greater than a critical limit ($\dot{M}_{\rm av,crit}$), and may approach a higher effective spin equilibrium value ($\nu_{\rm eq,eff}$) for transient accretion for $\dot{M}_{\rm av} < \dot{M}_{\rm av,crit}$. For example, when $\dot{M}_{\rm av}$ falls below $\dot{M}_{\rm av,crit}$ for an initially persistent source, $\nu$ increases considerably due to transient accretion, which is counterintuitive. We also find that, contrary to what was suggested, a fast or sudden decrease of $\dot{M}_{\rm av}$ to zero in the last part of the LMXB phase is not essential for the genesis of spin-powered MSPs, and neutron stars could spin up in this $\dot{M}_{\rm av}$-decreasing phase. Our findings imply that the traditional way of $\nu$-evolution computation is inadequate in most cases, even for initially persistent sources, and may not even correctly estimate whether $\nu$ increases or decreases.


title: Depth resolution optimized sputter depth profiling of a polycrystalline Al layer

abstract: Depth profiling of thin films by inert gas sputtering is one of the most important applications of Auger electron spectroscopy and X-ray photoelectron spectroscopy. Such an analysis monitors the elemental in-depth composition of the thin film system and controls the quality of the manufacturing process. A challenging task of a sputter depth profile measurement is the determination of the interface contaminations. Interface contaminations are only detectable if the depth profile data were measured with a sufficient depth resolution. Interface contaminations are crucial, because they strongly influence the mechanical and electrical properties of the whole thin film system. The depth resolution of unidirectional sputter depth profiling for polycrystalline samples is limited by the sputter yield differences attributed to grains having different crystalline orientations relative to the incoming ion beam. Therefore depth resolution optimized sputter depth profiling of polycrystalline thin films requires dedicated experimental approaches. If the sample is rotated during ion sputtering the depth profile is recorded with a better depth resolution because the ion impact direction relative to the grains orientation varies. If depth profiling can be performed on a single grain only, the depth resolution is improved, too. These approaches are applied to the depth profiling and interface contamination analysis of a 5000 nm thick polycrystalline Al layer. An O contamination at the interface of the Al to the Si substrate was detected utilizing these high depth resolution depth profiling methods.


title: Axion-Like Particles at the ILC Giga-Z

abstract: Axion-Like Particles (ALPs) are a generic, calculable, and well motivated extension of the Standard Model with far reaching phenomenology. ALPs that couple only to hypercharge represent one subset of such models, coupling the ALP to both photons and the $Z$ boson. We examine the current constraints on this class of models with an ALP mass in the 100 MeV to 100 GeV range, paying particular attention to the region between 100 MeV to 10 GeV, a portion of parameter space which is ill constrained by current experiments. We show that the more than $10^{9}$ $Z$ bosons produced in the Giga-Z mode of the future ILC experiment, combined with the highly granular nature of its detectors, will allow for ALPs coupled to hypercharge to be discovered with couplings down to nearly $10^{-5}\,\rm{GeV^{-1}}$ over a range of masses from 0.4 to 50 GeV.


title: Improving DGA-Based Malicious Domain Classifiers for Malware Defense with Adversarial Machine Learning

abstract: Domain Generation Algorithms (DGAs) are used by adversaries to establish Command and Control (C\&C) server communications during cyber attacks. Blacklists of known/identified C\&C domains are often used as one of the defense mechanisms. However, since blacklists are static and generated by signature-based approaches, they can neither keep up nor detect never-seen-before malicious domain names. Due to this shortcoming of blacklist domain checking, machine learning algorithms have been used to address the problem to some extent. However, when training is performed with limited datasets, the algorithms are likely to fail in detecting new DGA variants. To mitigate this weakness, we successfully applied a DGA-based malicious domain classifier using the Long Short-Term Memory (LSTM) method with a novel feature engineering technique. Our model's performance shows a higher level of accuracy compared to a previously reported model from prior research. Additionally, we propose a new method using adversarial machine learning to generate never-before-seen malware-related domain families that can be used to illustrate the shortcomings of machine learning algorithms in this regard. Next, we augment the training dataset with new samples such that it makes training of the machine learning models more effective in detecting never-before-seen malicious domain name variants. Finally, to protect blacklists of malicious domain names from disclosure and tampering, we devise secure data containers that store blacklists and guarantee their protection against adversarial access and modifications.


title: Domain Adaptation for the Segmentation of Confidential Medical Images

abstract: Convolutional neural networks (CNNs) have led to significant improvements in the semantic segmentation of images. When source and target datasets come from different modalities, CNN performance suffers due to domain shift. In such cases data annotation in the target domain becomes necessary to maintain model performance. To circumvent the re-annotation process, unsupervised domain adaptation (UDA) is proposed to adapt a model to new modalities using solely unlabeled target data. Common UDA algorithms require access to source domain data during adaptation, which may not be feasible in medical imaging due to data sharing restrictions. In this work, we develop an algorithm for UDA where the source domain data is inaccessible during target adaptation. Our approach is based on encoding the source domain information into an internal distribution that is used to guide adaptation in the absence of source samples. We demonstrate the effectiveness of our algorithm by comparing it to state-of-the-art medical image semantic segmentation approaches on two medical image semantic segmentation datasets.


title: Topological data analysis distinguishes parameter regimes in the Anderson-Chaplain model of angiogenesis

abstract: Angiogenesis is the process by which blood vessels form from pre-existing vessels. It plays a key role in many biological processes, including embryonic development and wound healing, and contributes to many diseases including cancer and rheumatoid arthritis. The structure of the resulting vessel networks determines their ability to deliver nutrients and remove waste products from biological tissues. Here we simulate the Anderson-Chaplain model of angiogenesis at different parameter values and quantify the vessel architectures of the resulting synthetic data. Specifically, we propose a topological data analysis (TDA) pipeline for systematic analysis of the model. TDA is a vibrant and relatively new field of computational mathematics for studying the shape of data. We compute topological and standard descriptors of model simulations generated by different parameter values. We show that TDA of model simulation data stratifies parameter space into regions with similar vessel morphology. The methodologies proposed here are widely applicable to other synthetic and experimental data including wound healing, development, and plant biology.


title: One-shot Representational Learning for Joint Biometric and Device Authentication

abstract: In this work, we propose a method to simultaneously perform (i) biometric recognition (i.e., identify the individual), and (ii) device recognition, (i.e., identify the device) from a single biometric image, say, a face image, using a one-shot schema. Such a joint recognition scheme can be useful in devices such as smartphones for enhancing security as well as privacy. We propose to automatically learn a joint representation that encapsulates both biometric-specific and sensor-specific features. We evaluate the proposed approach using iris, face and periocular images acquired using near-infrared iris sensors and smartphone cameras. Experiments conducted using 14,451 images from 15 sensors resulted in a rank-1 identification accuracy of upto 99.81% and a verification accuracy of upto 100% at a false match rate of 1%.


title: The autoregressive filter problem for multivariable degree one symmetric polynomials

abstract: The multivariable autoregressive filter problem asks for a polynomial $p(z)=p(z_1, \ldots , z_d)$ without roots in the closed $d$-disk based on prescribed Fourier coefficients of its spectral density function $1/|p(z)|^2$. The conditions derived in this paper for the construction of a degree one symmetric polynomial reveal a major divide between the case of at most two variables vs. the the case of three or more variables. The latter involves multivariable elliptic functions, while the former (due to [J. S. Geronimo and H. J. Woerdeman, Ann. of Math. (2), 160(3):839--906, 2004]) only involve polynomials. The three variable case is treated with more detail, and entails hypergeometric functions. Along the way, we identify a seemingly new relation between $_2F_1(\frac13,\frac23;1;z)$ and $_2F_1(\frac12,\frac12;1;\widetilde{z})$.


title: Analysis of some Epidemic Models in complex networks and some ideas about isolation strategies

abstract: Many models of virus propagation in Computer Networks inspired by {\bf SIS,SIR,}\\ {\bf SEIR}, etc. epidemic disease propagation mathematical models that can be found in the epidemiology field have been proposed in the last two decades. The purpose of these models has been to determine the conditions under which a virus becomes rapidly extinct in a network. The most common models proposed in the field of virus propagation in networks are inspired by SIS-type models or their variants. In such models, the conditions that lead to a rapid extinction of the spread of a computer virus have been calculated and its dependence on some parameters inherent to the mathematical model has been observed. In this article we will try to analyze a particular model proposed in the past and show through simulations the influence that topology has on the dynamics of the spread of a virus in different networks. A consequence of knowing the impact of the topology of a network can serve to propose effective isolation strategies to reduce the spread of a virus through modifications to the original network of contacts. I will talk about this subject at the final section of the present article.


title: Statistical Inference on the Hilbert Sphere with Application to Random Densities

abstract: The infinite-dimensional Hilbert sphere $S^\infty$ has been widely employed to model density functions and shapes, extending the finite-dimensional counterpart. We consider the Fr\'echet mean as an intrinsic summary of the central tendency of data lying on $S^\infty$. To break a path for sound statistical inference, we derive properties of the Fr\'echet mean on $S^\infty$ by establishing its existence and uniqueness as well as a root-$n$ central limit theorem (CLT) for the sample version, overcoming obstructions from infinite-dimensionality and lack of compactness on $S^\infty$. Intrinsic CLTs for the estimated tangent vectors and covariance operator are also obtained. Asymptotic and bootstrap hypothesis tests for the Fr\'echet mean based on projection and norm are then proposed and are shown to be consistent. The proposed two-sample tests are applied to make inference for daily taxi demand patterns over Manhattan modeled as densities, of which the square roots are analyzed on the Hilbert sphere. Numerical properties of the proposed hypothesis tests which utilize the spherical geometry are studied in the real data application and simulations, where we demonstrate that the tests based on the intrinsic geometry compare favorably to those based on an extrinsic or flat geometry.


title: Toward a flow-structure-based wall-modeled large-eddy simulation paradigm

abstract: A promising and cost-effective method for numerical simulation of high Re wall-bounded flows is wall-modeled large-eddy simulation. Most wall models are formulated from the Reynolds-averaged Navier-Stokes equations (RANS). These RANS-based wall models are calibrated using mean turbulence data and make no use of the current vast knowledge on turbulent flow structure. Moreover, RANS-based models are limited to predicting the mean velocity profile, and the mean wall shear stress. Using the knowledge of the wall-normal self-similarity of high-$Re_\tau$ wall-bounded turbulent flows, we present a coupling between a near-wall patch of DNS resolution fixed in inner-units and an outer LES flow field. The near-wall patch captures the near-wall self-sustaining cycle as well as the lower portion of the self-similar hierarchy of eddies. Given that both the near-wall patch and the LES capture separate portions of the self-similar hierarchy of eddies an instantaneous top boundary condition for the patch is formulated. The near-wall model is capable of predicting subgrid-scale quantities such as the wall stress, velocity fluctuations, kinetic energy spectra, and flow structure across the entire near-wall layer.


title: VinVL: Revisiting Visual Representations in Vision-Language Models

abstract: This paper presents a detailed study of improving visual representations for vision language (VL) tasks and develops an improved object detection model to provide object-centric representations of images. Compared to the most widely used \emph{bottom-up and top-down} model \cite{anderson2018bottom}, the new model is bigger, better-designed for VL tasks, and pre-trained on much larger training corpora that combine multiple public annotated object detection datasets. Therefore, it can generate representations of a richer collection of visual objects and concepts. While previous VL research focuses mainly on improving the vision-language fusion model and leaves the object detection model improvement untouched, we show that visual features matter significantly in VL models. In our experiments we feed the visual features generated by the new object detection model into a Transformer-based VL fusion model \oscar \cite{li2020oscar}, and utilize an improved approach \short\ to pre-train the VL model and fine-tune it on a wide range of downstream VL tasks. Our results show that the new visual features significantly improve the performance across all VL tasks, creating new state-of-the-art results on seven public benchmarks. We will release the new object detection model to public.


title: Polarized radiation and the Emergence of Biological Homochirality on Earth and Beyond

abstract: It has been proposed that spin-polarized cosmic radiation can induce asymmetric changes in helical biopolymers that may account for the emergence of biological homochirality. The parity violation in the weak interaction has direct consequences on the transport of polarization in cosmic ray showers. In this paper, we show that muons retain their polarization down to energies at which they can initiate enantioselective mutagenesis. Therefore, muons are most likely to succeed in establishing the connection between broken symmetries in the standard model of particle physics and that found in living organisms. We calculate the radiation doses deposited by primary and secondary cosmic rays at various prime targets for the searches of life in the solar system: Mars, Venus, Titan, icy moons and planetesimals, and discuss the implications for the enantioselective mutagenesis proposed as to be the driver of homochiralization. Earth is unusual in that spin-polarized muons dominate the cosmic radiation at its surface.


title: Context-Aware Safe Reinforcement Learning for Non-Stationary Environments

abstract: Safety is a critical concern when deploying reinforcement learning agents for realistic tasks. Recently, safe reinforcement learning algorithms have been developed to optimize the agent's performance while avoiding violations of safety constraints. However, few studies have addressed the non-stationary disturbances in the environments, which may cause catastrophic outcomes. In this paper, we propose the context-aware safe reinforcement learning (CASRL) method, a meta-learning framework to realize safe adaptation in non-stationary environments. We use a probabilistic latent variable model to achieve fast inference of the posterior environment transition distribution given the context data. Safety constraints are then evaluated with uncertainty-aware trajectory sampling. The high cost of safety violations leads to the rareness of unsafe records in the dataset. We address this issue by enabling prioritized sampling during model training and formulating prior safety constraints with domain knowledge during constrained planning. The algorithm is evaluated in realistic safety-critical environments with non-stationary disturbances. Results show that the proposed algorithm significantly outperforms existing baselines in terms of safety and robustness.


title: A Warped Resolvent Algorithm to Construct Nash Equilibria

abstract: We propose an asynchronous block-iterative decomposition algorithm to solve Nash equilibrium problems involving a mix of nonsmooth and smooth functions acting on linear mixtures of strategies. The methodology relies heavily on monotone operator theory and in particular on warped resolvents.


title: Cutoff phenomenon for the warp-transpose top with random shuffle

abstract: Let $\{G_n\}_1^{\infty}$ be a sequence of non-trivial finite groups. In this paper, we study the properties of a random walk on the complete monomial group $G_n\wr S_n$ generated by the elements of the form $(\text{e},\dots,\text{e},g;\text{id})$ and $(\text{e},\dots,\text{e},g^{-1},\text{e},\dots,\text{e},g;(i,n))$ for $g\in G_n,\;1\leq i< n$. We call this the warp-transpose top with random shuffle on $G_n\wr S_n$. We find the spectrum of the transition probability matrix for this shuffle. We prove that the mixing time for this shuffle is $O\left(n\log n+\frac{1}{2}n\log (|G_n|-1)\right)$. We show that this shuffle exhibits $\ell^2$-cutoff at $n\log n+\frac{1}{2}n\log (|G_n|-1)$ and total variation cutoff at $n\log n$.


title: Multiple ergodic averages for variable polynomials

abstract: In this paper we study multiple ergodic averages for "good" variable polynomials. In particular, under an additional assumption, we show that these averages converge to the expected limit, making progress related to an open problem posted by Frantzikinakis (Problem 10, "Some open problems on multiple ergodic averages. Bulletin of the Hellenic Mathematical Society. 60 (2016), 41-90"). Corresponding averages along prime numbers are studied too. These general convergence results imply various variable extensions of classical recurrence, combinatorial and number theoretical results which are presented as well.


title: RV-GAN: Segmenting Retinal Vascular Structure in Fundus Photographs using a Novel Multi-scale Generative Adversarial Network

abstract: High fidelity segmentation of both macro and microvascular structure of the retina plays a pivotal role in determining degenerative retinal diseases, yet it is a difficult problem. Due to successive resolution loss in the encoding phase combined with the inability to recover this lost information in the decoding phase, autoencoding based segmentation approaches are limited in their ability to extract retinal microvascular structure. We propose RV-GAN, a new multi-scale generative architecture for accurate retinal vessel segmentation to alleviate this. The proposed architecture uses two generators and two multi-scale autoencoding discriminators for better microvessel localization and segmentation. In order to avoid the loss of fidelity suffered by traditional GAN-based segmentation systems, we introduce a novel weighted feature matching loss. This new loss incorporates and prioritizes features from the discriminator's decoder over the encoder. Doing so combined with the fact that the discriminator's decoder attempts to determine real or fake images at the pixel level better preserves macro and microvascular structure. By combining reconstruction and weighted feature matching loss, the proposed architecture achieves an area under the curve (AUC) of 0.9887, 0.9914, and 0.9887 in pixel-wise segmentation of retinal vasculature from three publicly available datasets, namely DRIVE, CHASE-DB1, and STARE, respectively. Additionally, RV-GAN outperforms other architectures in two additional relevant metrics, mean intersection-over-union (Mean-IOU) and structural similarity measure (SSIM).


title: Computing Cliques and Cavities in Networks

abstract: Complex networks contain complete subgraphs such as nodes, edges, triangles, etc., referred to as simplices and cliques of different orders. Notably, cavities consisting of higher-order cliques play an important role in brain functions. Since searching for maximum cliques is an NP-complete problem, we use k-core decomposition to determine the computability of a given network. For a computable network, we design a search method with an implementable algorithm for finding cliques of different orders, obtaining also the Euler characteristic number. Then, we compute the Betti numbers by using the ranks of boundary matrices of adjacent cliques. Furthermore, we design an optimized algorithm for finding cavities of different orders. Finally, we apply the algorithm to the neuronal network of C. elegans with data from one typical dataset, and find all of its cliques and some cavities of different orders, providing a basis for further mathematical analysis and computation of its structure and function.


title: Intersections of Deligne--Lusztig varieties and Springer fibres

abstract: In this paper we prove a direct geometric relation between Deligne--Lusztig varieties and Springer fibres in type $\mathsf{A}$: For any rational unipotent element, the Springer fibre cuts out a unique component of a specific Deligne--Lusztig variety; moreover, this component forms an open dense subset of a component of the Springer fibre. This boils down to a map from the unipotent variety to the Weyl group, and combines several constructions with a combinatorial flavour (like Weyr normal forms, Robinson--Schensted correspondence, and Spaltenstein's and Steinberg's labellings); it also provides a geometric interpretation of a classical dimension formula of unipotent centralisers.


title: On a strengthening of the Blaschke-Leichtweiss theorem

abstract: The Blaschke-Leichtweiss theorem (Abh. Math. Sem. Univ. Hamburg 75: 257-284, 2005) states that the smallest area convex domain of constant width $w$ in the $2$-dimensional spherical space ${\mathbb S}^2$ is the spherical Reuleaux triangle for all $0<w\leq\frac{\pi}{2}$. In this paper we extend this result to the family of wide $r$-disk domains of ${\mathbb S}^2$, where $0<r\leq\frac{\pi}{2}$. Here a wide $r$-disk domain is an intersection of spherical disks of radius $r$ with centers contained in their intersection. This gives a new and elementary proof of the Blaschke-Leichtweiss theorem. Furthermore, we investigate the higher dimensional analogue of wide $r$-disk domains called wide $r$-ball bodies. In particular, we determine their minimum spherical width (resp., inradius) in the spherical $d$-space ${\mathbb S}^d$ for all $d\geq 2$. Also, it is shown that any minimum volume wide $r$-ball body is of constant width $r$ in ${\mathbb S}^d$, $d\geq 2$.


title: Two-dimensional superconducting fluctuations associated with charge density wave stripes in La$_{1.87}$Sr$_{0.13}$Cu$_{0.99}$Fe$_{0.01}$O$_4$

abstract: The presence of a small concentration of in-plane Fe dopants in La$_{1.87}$Sr$_{0.13}$Cu$_{0.99}$Fe$_{0.01}$O$_4$ is known to enhance stripe-like spin and charge density wave (SDW and CDW) order, and suppress the superconducting $T_c$. Here, we show that it also induces highly two-dimensional (2D) superconducting correlations that have been argued to be signatures of a new form of superconducting order, so-called pair-density-wave (PDW) order. In addition, using the resonant soft x-ray scattering, we find that the 2D superconducting fluctuation is strongly associated with the CDW stripe. In particular, the PDW signature first appears when the correlation length of the CDW stripe grows over eight times the lattice unit ($\sim$ 8$a$). These results provide critical conditions for the formation of PDW order.


title: Attentive Tree-structured Network for Monotonicity Reasoning

abstract: Many state-of-art neural models designed for monotonicity reasoning perform poorly on downward inference. To address this shortcoming, we developed an attentive tree-structured neural network. It consists of a tree-based long-short-term-memory network (Tree-LSTM) with soft attention. It is designed to model the syntactic parse tree information from the sentence pair of a reasoning task. A self-attentive aggregator is used for aligning the representations of the premise and the hypothesis. We present our model and evaluate it using the Monotonicity Entailment Dataset (MED). We show and attempt to explain that our model outperforms existing models on MED.


title: Time fractional gradient flows: Theory and numerics

abstract: We develop the theory of fractional gradient flows: an evolution aimed at the minimization of a convex, l.s.c.~energy, with memory effects. This memory is characterized by the fact that the negative of the (sub)gradient of the energy equals the so-called Caputo derivative of the state. We introduce the notion of energy solutions, for which we provide existence, uniqueness and certain regularizing effects. We also consider Lipschitz perturbations of this energy. For these problems we provide an a posteriori error estimate and show its reliability. This estimate depends only on the problem data, and imposes no constraints between consecutive time-steps. On the basis of this estimate we provide an a priori error analysis that makes no assumptions on the smoothness of the solution.


title: An Efficient Transformer Decoder with Compressed Sub-layers

abstract: The large attention-based encoder-decoder network (Transformer) has become prevailing recently due to its effectiveness. But the high computation complexity of its decoder raises the inefficiency issue. By examining the mathematic formulation of the decoder, we show that under some mild conditions, the architecture could be simplified by compressing its sub-layers, the basic building block of Transformer, and achieves a higher parallelism. We thereby propose Compressed Attention Network, whose decoder layer consists of only one sub-layer instead of three. Extensive experiments on 14 WMT machine translation tasks show that our model is 1.42x faster with performance on par with a strong baseline. This strong baseline is already 2x faster than the widely used standard baseline without loss in performance.


title: Centralized and Distributed Age of Information Minimization with non-linear Aging Functions in the Internet of Things

abstract: Resource management in Internet of Things (IoT) systems is a major challenge due to the massive scale and heterogeneity of the IoT system. For instance, most IoT applications require timely delivery of collected information, which is a key challenge for the IoT. In this paper, novel centralized and distributed resource allocation schemes are proposed to enable IoT devices to share limited communication resources and to transmit IoT messages in a timely manner. In the considered system, the timeliness of information is captured using non-linear age of information (AoI) metrics that can naturally quantify the freshness of information. To model the inherent heterogeneity of the IoT system, the non-linear aging functions are defined in terms of IoT device types and message content. To minimize AoI, the proposed resource management schemes allocate the limited communication resources considering AoI. In particular, the proposed centralized scheme enables the base station to learn the device types and to determine aging functions. Moreover, the proposed distributed scheme enables the devices to share the limited communication resources based on available information on other devices and their AoI. The convergence of the proposed distributed scheme is proved, and the effectiveness in reducing the AoI with partial information is analyzed. Simulation results show that the proposed centralized scheme achieves significantly lower average instantaneous AoI when compared to simple centralized allocation without learning, while the proposed distributed scheme achieves significantly lower average instantaneous AoI when compared to random allocation. The results also show that the proposed centralized scheme outperforms the proposed distributed scheme in almost all cases, but the distributed approach is more viable for a massive IoT.


title: On the non-very generic intersections in discriminantal arrangements

abstract: In 1985 Crapo introduced in \cite{Crapo} a new mathematical object that he called $\textit{geometry of circuits}$. Four years later, in 1989, Manin and Schechtman defined in \cite{MS} the same object and called it $\textit{discriminantal arrangement}$, the name by which it is known now a days. Those discriminantal arrangements $\mathcal{B}(n,k,\mathcal{A}^0)$ are builded from an arrangement $\mathcal{A}^0$ of $n$ hyperplanes in general position in a $k$-dimensional space and their combinatorics depends on the arrangement $\mathcal{A}^0$. On this basis, in 1997 Bayer and Brandt (see \cite{BB}) distinguished two different type of arrangements $\mathcal{A}^0$ calling $\textit{very generic}$ the ones for which the intersection lattice of $\mathcal{B}(n,k,\mathcal{A}^0)$ has maximum cardinality and $\textit{non-very generic}$ the others. Results on the combinatorics of $\mathcal{B}(n,k,\mathcal{A}^0)$ in the very generic case already appear in Crapo \cite{Crapo} and in 1997 in Athanasiadis \cite{Atha} while the first known result on non-very generic case is due to Libgober and the first author in 2018. In their paper \cite{LS} they provided a necessary and sufficient condition on $\mathcal{A}^0$ for which the cardinality of rank 2 intersections in $\mathcal{B}(n,k,\mathcal{A}^0)$ is not maximal anymore. In this paper we further develop their result providing a sufficient condition on $\mathcal{A}^0$ for which the cardinality of rank r, $r \geq 2$, intersections in $\mathcal{B}(n,k,\mathcal{A}^0)$ decreases.


title: A Hybrid Attention Mechanism for Weakly-Supervised Temporal Action Localization

abstract: Weakly supervised temporal action localization is a challenging vision task due to the absence of ground-truth temporal locations of actions in the training videos. With only video-level supervision during training, most existing methods rely on a Multiple Instance Learning (MIL) framework to predict the start and end frame of each action category in a video. However, the existing MIL-based approach has a major limitation of only capturing the most discriminative frames of an action, ignoring the full extent of an activity. Moreover, these methods cannot model background activity effectively, which plays an important role in localizing foreground activities. In this paper, we present a novel framework named HAM-Net with a hybrid attention mechanism which includes temporal soft, semi-soft and hard attentions to address these issues. Our temporal soft attention module, guided by an auxiliary background class in the classification module, models the background activity by introducing an "action-ness" score for each video snippet. Moreover, our temporal semi-soft and hard attention modules, calculating two attention scores for each video snippet, help to focus on the less discriminative frames of an action to capture the full action boundary. Our proposed approach outperforms recent state-of-the-art methods by at least 2.2% mAP at IoU threshold 0.5 on the THUMOS14 dataset, and by at least 1.3% mAP at IoU threshold 0.75 on the ActivityNet1.2 dataset. Code can be found at: https://github.com/asrafulashiq/hamnet.


title: Optimal stopping time on discounted semi-Markov processes

abstract: This paper attempts to study the optimal stopping time for semi-Markov processes (SMPs) under the discount optimization criteria with unbounded cost rates. In our work, we introduce an explicit construction of the equivalent semi-Markov decision processes (SMDPs). The equivalence is embodied in the value functions of SMPs and SMDPs, that is, every stopping time of SMPs can induce a policy of SMDPs such that the value functions are equal, and vice versa. The existence of the optimal stopping time of SMPs is proved by this equivalence relation. Next, we give the optimality equation of the value function and develop an effective iterative algorithm for computing it. Moreover, we show that the optimal and {\epsilon}-optimal stopping time can be characterized by the hitting time of the special sets. Finally, to illustrate the validity of our results, an example of a maintenance system is presented in the end.


title: Driven-dissipative control of cold atoms in tilted optical lattices

abstract: We present a sequence of driven-dissipative protocols for controlling cold atoms in tilted optical lattices. These experimentally accessible examples are templates that demonstrate how dissipation can be used to manipulate quantum many-body systems. We consider bosonic atoms trapped in a tilted optical lattice, immersed in a superfluid bath, and excited by coherent Raman lasers. With these ingredients, we are able to controllably transport atoms in the lattice and produce self-healing quantum states: a Mott insulator and the topologically ordered spin-1 AKLT state.


title: Better understanding of the multivariate hypergeometric distribution with implications in design-based survey sampling

abstract: Multivariate hypergeometric distribution arises frequently in elementary statistics and probability courses, for simultaneously studying the occurence law of specified events, when sampling without replacement from a finite population with fixed number of classification. Covariance matrix of this distribution is well known to be identical to its multinomial counterpart multiplied by 1-(n-1)/(N-1), with N and n being population and sample sizes, respectively. It appears to however, have been less discussed in the literature about the meaning of this relationship, especially regarding the specific form of the multiplier. Based on an augmenting argument together with probabilistic symmetry, we present a more transparent understanding for the covariance structure of the multivariate hypergeometric distribution. We discuss implications of these combined techniques and provide a unified description about the relative efficiency for estimating population mean based on simple random sampling, probability proportional-to-size sampling and adaptive cluster sampling, with versus without replacement. We also provide insight into the classic random group method for variance estimation.


title: Study of the ${\Upsilon}(1S)$ ${\to}$ $DP$ decays

abstract: Inspired by the potential prospects of high-luminosity dedicated colliders and the high enthusiasms in searching for new physics in the flavor sector at the intensity frontier, the ${\Upsilon}(1S)$ ${\to}$ $D^{-}{\pi}^{+}$, $\overline{D}^{0}{\pi}^{0}$ and $D_{s}^{-}K^{+}$ weak decays are studied with the perturbative QCD approach. It is found within the standard model that the branching ratios for the concerned processes are tiny, about ${\cal O}(10^{-18})$, and far beyond the detective ability of current experiments unless there exists some significant enhancements from a noval interaction.


title: Percolation and the pandemic

abstract: This paper is dedicated to the memory of Dietrich Stauffer, who was a pioneer in percolation theory and applications of it to problems of society, such as epidemiology. An epidemic is a percolation process gone out of control, that is, going beyond the critical transition threshold $p_c$. Here we discuss how the threshold is related to the basic infectivity of neighbors $R_0$, for trees (Bethe lattice), trees with triangular cliques, and in non-planar lattice percolation with extended-range connectivity. It is shown how having a smaller range of contacts increases the critical value of $R_0$ above the value $R_{0,c}=1$ appropriate for a tree, an infinite-range system or a large completely connected graph.


title: Aspects of quantum states asymmetry for the magnetic dipolar interaction dynamics

abstract: We investigate the asymmetry properties of quantum states in relation to the Hamiltonian responsible for the magnetic dipolar interaction (MDI) dynamics, and we evaluate its relationship to entanglement production. We consider some classes of pure and mixed quantum states of two qubits evolved under MDI and, using the asymmetry measure defined via the Wigner-Yanase skew information, we describe the asymmetry dependence on the Hamiltonian parameters and initial conditions of the system. In addition, we define and calculate the dynamics of the asymmetry of local states, characterizing their temporal and interaction parameters dependence. Finally, because the MDI Hamiltonian has a null eigenvalue, the group generator-based asymmetry measure does not adequately quantify the state susceptibility with respect to the action of the subspace generated by the eigenvectors associated with this eigenvalue. For this reason, we also define and study the group element-based asymmetry measure with relation to the unitary operator associated with the MDI Hamiltonian.


title: Hyponormal dual Toeplitz operators on the orthogonal complement of the Harmonic Bergman space

abstract: In this paper, we mainly study the hyponormality of dual Toeplitz operators on the orthogonal complement of the harmonic Bergman space. First we show that the dual Toeplitz operator with bounded symbol is hyponormal if and only if it is normal. Then we obtain a necessary and sufficient condition for the dual Toeplitz operator $S_\varphi$ with the symbol $\varphi(z) = az^{n_1}\overline{z}^{m_1} + bz^{n_2} \overline{z}^{m_2}$, $(n_1,n_2,m_1,m_2\in \mathbb {N}$ and $a,b \in \mathbb{C})$ to be hyponormal. Finally, we show that the rank of the commutator of two dual Toeplitz operators must be an even number if the commutator has a finite rank.


title: Transcendental pairs of generic extensions

abstract: We isolate a new preservation class of Suslin forcings and prove several associated consistency results in the choiceless theory ZF+DC regarding countable chromatic numbers of various Borel hypergraphs.


title: Global field reconstruction from sparse sensors with Voronoi tessellation-assisted deep learning

abstract: Achieving accurate and robust global situational awareness of a complex time-evolving field from a limited number of sensors has been a longstanding challenge. This reconstruction problem is especially difficult when sensors are sparsely positioned in a seemingly random or unorganized manner, which is often encountered in a range of scientific and engineering problems. Moreover, these sensors can be in motion and can become online or offline over time. The key leverage in addressing this scientific issue is the wealth of data accumulated from the sensors. As a solution to this problem, we propose a data-driven spatial field recovery technique founded on a structured grid-based deep-learning approach for arbitrary positioned sensors of any numbers. It should be noted that the na\"ive use of machine learning becomes prohibitively expensive for global field reconstruction and is furthermore not adaptable to an arbitrary number of sensors. In the present work, we consider the use of Voronoi tessellation to obtain a structured-grid representation from sensor locations enabling the computationally tractable use of convolutional neural networks. One of the central features of the present method is its compatibility with deep-learning based super-resolution reconstruction techniques for structured sensor data that are established for image processing. The proposed reconstruction technique is demonstrated for unsteady wake flow, geophysical data, and three-dimensional turbulence. The current framework is able to handle an arbitrary number of moving sensors, and thereby overcomes a major limitation with existing reconstruction methods. The presented technique opens a new pathway towards the practical use of neural networks for real-time global field estimation.


title: On Koopman Mode Decomposition and Tensor Component Analysis

abstract: Koopman mode decomposition and tensor component analysis (also known as CANDECOMP/PARAFAC or canonical polyadic decomposition) are two popular approaches of decomposing high dimensional data sets into low dimensional modes that capture the most relevant features and/or dynamics. Despite their similar goal, the two methods are largely used by different scientific communities and formulated in distinct mathematical languages. We examine the two together and show that, under a certain (reasonable) condition on the data, the theoretical decomposition given by tensor component analysis is the \textit{same} as that given by Koopman mode decomposition. This provides a "bridge" with which the two communities should be able to more effectively communicate. When this condition is not met, Koopman mode decomposition still provides a tensor decomposition with an \textit{a priori} computable error, providing an alternative to the non-convex optimization that tensor component analysis requires. Our work provides new possibilities for algorithmic approaches to Koopman mode decomposition and tensor component analysis, provides a new perspective on the success of tensor component analysis, and builds upon a growing body of work showing that dynamical systems, and Koopman operator theory in particular, can be useful for problems that have historically made use of optimization theory.


title: Zero-dynamics Attack, Variations, and Countermeasures

abstract: This chapter presents an overview on actuator attacks that exploit zero dynamics, and countermeasures against them. First, zero-dynamics attack is re-introduced based on a canonical representation called normal form. Then it is shown that the target dynamic system is at elevated risk if the associated zero dynamics is unstable. From there on, several questions are raised in series to ensure when the target system is immune to the attack of this kind. The first question is: Is the target system secure from zero-dynamics attack if it does not have any unstable zeros? An answer provided for this question is: No, the target system may still be at risk due to another attack surface emerging in the process of implementation. This is followed by a series of next questions, and in the course of providing answers, variants of the classic zero-dynamics attack are presented, from which the vulnerability of the target system is explored in depth. At the end, countermeasures are proposed to render the attack ineffective. Because it is known that the zero-dynamics in continuous-time systems cannot be modified by feedback, the main idea of the countermeasure is to relocate any unstable zero to a stable region in the stage of digital implementation through modified digital samplers and holders. Adversaries can still attack actuators, but due to the re-located zeros, they are of little use in damaging the target system.


title: Silicon Photonic Microring Based Chip-Scale Accelerator for Delayed Feedback Reservoir Computing

abstract: To perform temporal and sequential machine learning tasks, the use of conventional Recurrent Neural Networks (RNNs) has been dwindling due to the training complexities of RNNs. To this end, accelerators for delayed feedback reservoir computing (DFRC) have attracted attention in lieu of RNNs, due to their simple hardware implementations. A typical implementation of a DFRC accelerator consists of a delay loop and a single nonlinear neuron, together acting as multiple virtual nodes for computing. In prior work, photonic DFRC accelerators have shown an undisputed advantage of fast computation over their electronic counterparts. In this paper, we propose a more energy-efficient chip-scale DFRC accelerator that employs a silicon photonic microring (MR) based nonlinear neuron along with on-chip photonic waveguides-based delayed feedback loop. Our evaluations show that, compared to a well-known photonic DFRC accelerator from prior work, our proposed MR-based DFRC accelerator achieves 35% and 98.7% lower normalized root mean square error (NRMSE), respectively, for the prediction tasks of NARMA10 and Santa Fe time series. In addition, our MR-based DFRC accelerator achieves 58.8% lower symbol error rate (SER) for the Non-Linear Channel Equalization task. Moreover, our MR-based DFRC accelerator has 98% and 93% faster training time, respectively, compared to an electronic and a photonic DFRC accelerators from prior work.


title: Partial regularity for an exponential PDE in crystal surface models

abstract: We study the regularity properties of a weak solution to the boundary value problem for the equation $-\Delta \rho +a u=f$ in a bounded domain $\Omega\subset \mathbb{R}^N$, where $\rho=e^{-\mbox{div}\left(|\nabla u|^{p-2}\nabla u+\beta_0|\nabla u|^{-1}\nabla u\right)}$. This problem is derived from the mathematical modeling of crystal surfaces. It is known that the exponent term can exhibit singularity. In this paper we obtain a partial regularity result for the weak solution. It asserts that there exists an open subset $\Omega_0\subset \Omega$ such that $|\Omega\setminus\Omega_0|=0$ and the exponent term is locally bounded in $\Omega_0$. Furthermore, if $x_0\in \Omega\setminus\Omega_0$, then $\rho$ vanishes of $N+2-\varepsilon$ order at $x_0$ for each $\varepsilon\in(0,2)$. Our results reveal that the exponent term behaves well if it stays away from negative infinity.


title: Existence and the reducibility of the Hilbert scheme of linearly normal curves in $\mathbb{P}^r$ of relatively high degrees

abstract: Let $\mathcal{H}_{d,g,r}$ be the Hilbert scheme parametrizing smooth irreducible and non-degenerate curves of degree $d$ and genus $g$ in $\mathbb{P}^r.$ We denote by $\mathcal{H}^\mathcal{L}_{d,g,r}$ the union of those components of $\mathcal{H}_{d,g,r}$ whose general element is linearly normal. In this article we show that $\mathcal{H}^\mathcal{L}_{d,g,r}$ ($d\ge g+r-3$) is non-empty in a certain optimal range of triples $(d,g,r)$ and is empty outside the range. This settles the existence (or non-emptiness if one prefers) of the Hilbert scheme $\mathcal{H}^\mathcal{L}_{d,g,r}$ of linearly normal curves of degree $d$ and genus $g$ in $\mathbb{P}^r$ for $g+r-3\le d\le g+r$, $r\ge 3$. We also determine all the triples $(d,g,r)$ with $g+r-3\le d\le g+r$ for which $\mathcal{H}^\mathcal{L}_{d,g,r}$ is reducible (or irreducible).


title: Gaussian Happy Numbers

abstract: This paper extends the concept of a $B$-happy number, for $B \geq 2$, from the rational integers, $\mathbb{Z}$, to the Gaussian integers, $\mathbb{Z}[i]$. We investigate the fixed points and cycles of the Gaussian $B$-happy functions, determining them for small values of $B$ and providing a method for computing them for any $B \geq 2$. We discuss heights of Gaussian $B$-happy numbers, proving results concerning the smallest Gaussian $B$-happy numbers of certain heights. Finally, we prove conditions for the existence and non-existence of arbitrarily long arithmetic sequences of Gaussian $B$-happy numbers.


title: Six-channel Image Representation for Cross-domain Object Detection

abstract: Most deep learning models are data-driven and the excellent performance is highly dependent on the abundant and diverse datasets. However, it is very hard to obtain and label the datasets of some specific scenes or applications. If we train the detector using the data from one domain, it cannot perform well on the data from another domain due to domain shift, which is one of the big challenges of most object detection models. To address this issue, some image-to-image translation techniques have been employed to generate some fake data of some specific scenes to train the models. With the advent of Generative Adversarial Networks (GANs), we could realize unsupervised image-to-image translation in both directions from a source to a target domain and from the target to the source domain. In this study, we report a new approach to making use of the generated images. We propose to concatenate the original 3-channel images and their corresponding GAN-generated fake images to form 6-channel representations of the dataset, hoping to address the domain shift problem while exploiting the success of available detection models. The idea of augmented data representation may inspire further study on object detection and other applications.


title: Few-shot Image Classification: Just Use a Library of Pre-trained Feature Extractors and a Simple Classifier

abstract: Recent papers have suggested that transfer learning can outperform sophisticated meta-learning methods for few-shot image classification. We take this hypothesis to its logical conclusion, and suggest the use of an ensemble of high-quality, pre-trained feature extractors for few-shot image classification. We show experimentally that a library of pre-trained feature extractors combined with a simple feed-forward network learned with an L2-regularizer can be an excellent option for solving cross-domain few-shot image classification. Our experimental results suggest that this simpler sample-efficient approach far outperforms several well-established meta-learning algorithms on a variety of few-shot tasks.


title: Learning Neural Networks on SVD Boosted Latent Spaces for Semantic Classification

abstract: The availability of large amounts of data and compelling computation power have made deep learning models much popular for text classification and sentiment analysis. Deep neural networks have achieved competitive performance on the above tasks when trained on naive text representations such as word count, term frequency, and binary matrix embeddings. However, many of the above representations result in the input space having a dimension of the order of the vocabulary size, which is enormous. This leads to a blow-up in the number of parameters to be learned, and the computational cost becomes infeasible when scaling to domains that require retaining a colossal vocabulary. This work proposes using singular value decomposition to transform the high dimensional input space to a lower-dimensional latent space. We show that neural networks trained on this lower-dimensional space are not only able to retain performance while savoring significant reduction in the computational complexity but, in many situations, also outperforms the classical neural networks trained on the native input space.


title: Mapping class groups of surfaces with noncompact boundary components

abstract: We show that the pure mapping class group is uniformly perfect for a certain class of infinite type surfaces with noncompact boundary components. We then combine this result with recent work in the remaining cases to give a complete classification of the perfect and uniformly perfect pure mapping class groups for infinite type surfaces. We also develop a method to cut a general surface into simpler surfaces and extend some mapping class group results to the general case.


title: Estimation of Tempered Stable L\'{e}vy Models of Infinite Variation

abstract: We propose a new method for the estimation of a semiparametric tempered stable L\'{e}vy model. The estimation procedure combines iteratively an approximate semiparametric method of moment estimator, Truncated Realized Quadratic Variations (TRQV), and a newly found small-time high-order approximation for the optimal threshold of the TRQV of tempered stable processes. The method is tested via simulations to estimate the volatility and the Blumenthal-Getoor index of the generalized CGMY model as well as the integrated volatility of a Heston-type model with CGMY jumps. The method outperforms other efficient alternatives proposed in the literature when working with a L\'evy process (i.e., the volatility is constant), or when the index of jump intensity $Y$ is larger than $3/2$ in the presence of stochastic volatility.


title: Airplane-Aided Integrated Next-Generation Networking

abstract: A high-rate yet low-cost air-to-ground (A2G) communication backbone is conceived for integrating the space and terrestrial network by harnessing the opportunistic assistance of the passenger planes or high altitude platforms (HAPs) as mobile base stations (BSs) and millimetre wave communication. The airliners act as the network-provider for the terrestrial users while relying on satellite backhaul. Three different beamforming techniques relying on a large-scale planar array are used for transmission by the airliner/HAP for achieving a high directional gain, hence minimizing the interference among the users. Furthermore, approximate spectral efficiency (SE) and area spectral efficiency (ASE) expressions are derived and quantified for diverse system parameters.


title: ASIST: Annotation-free Synthetic Instance Segmentation and Tracking by Adversarial Simulations

abstract: Background: The quantitative analysis of microscope videos often requires instance segmentation and tracking of cellular and subcellular objects. The traditional method consists of two stages: (1) performing instance object segmentation of each frame, and (2) associating objects frame-by-frame. Recently, pixel-embedding-based deep learning approaches these two steps simultaneously as a single stage holistic solution. In computer vision, annotated training data with consistent segmentation and tracking is resource intensive, the severity of which is multiplied in microscopy imaging due to (1) dense objects (e.g., overlapping or touching), and (2) high dynamics (e.g., irregular motion and mitosis). Adversarial simulations have provided successful solutions to alleviate the lack of such annotations in dynamics scenes in computer vision, such as using simulated environments (e.g., computer games) to train real-world self-driving systems. Methods: In this paper, we propose an annotation-free synthetic instance segmentation and tracking (ASIST) method with adversarial simulation and single-stage pixel-embedding based learning. Contribution: The contribution of this paper is three-fold: (1) the proposed method aggregates adversarial simulations and single-stage pixel-embedding based deep learning; (2) the method is assessed with both the cellular (i.e., HeLa cells) and subcellular (i.e., microvilli) objects; and (3) to the best of our knowledge, this is the first study to explore annotation-free instance segmentation and tracking study for microscope videos. Results: The ASIST method achieved an important step forward, when compared with fully supervised approaches: ASIST shows 7% to 11% higher segmentation, detection and tracking performance on microvilli relative to fully supervised methods, and comparable performance on Hela cell videos.


title: Almost automorphy of minimal sets for $C^1$-smooth strongly monotone skew-product semiflows on Banach spaces

abstract: We focus on the presence of almost automorphy in strongly monotone skew-product semiflows on Banach spaces. Under the $C^1$-smoothness assumption, it is shown that any linearly stable minimal set must be almost automorphic. This extends the celebrated result of Shen and Yi [Mem. Amer. Math. Soc. 136(1998), No. 647] for the classical $C^{1,\alpha}$-smooth systems. Based on this, one can reduce the regularity of the almost periodically forced differential equations and obtain the almost automorphic phenomena in a wider range.


title: Magnetic surface reconstruction in the van-der-Waals antiferromagnet Fe$_{1+x}$Te

abstract: Fe$_{1+x}$Te is a two dimensional van der Waals antiferromagnet that becomes superconducting on anion substitution on the Te site. The parent phase of Fe$_{1+x}$Te is sensitive to the amount of interstitial iron situated between the iron-tellurium layers displaying collinear magnetic order coexisting with low temperature metallic resistivity for small concentrations of interstitial iron $x$ and helical magnetic order for large values of $x$. While this phase diagram has been established through scattering [see for example E. E. Rodriguez $\textit{et al.}$ Phys. Rev. B ${\bf{84}}$, 064403 (2011) and S. R\"ossler $\textit{et al.}$ Phys. Rev. B ${\bf{84}}$, 174506 (2011)], recent scanning tunnelling microscopy measurements [C. Trainer $\textit{et al.}$ Sci. Adv. ${\bf{5}}$, eaav3478 (2019)] have observed a different magnetic structure for small interstitial iron concentrations $x$ with a significant canting of the magnetic moments along the crystallographic $c$ axis of $\theta$=28 $\pm$ 3$^{\circ}$. In this paper, we revisit the magnetic structure of Fe$_{1.09}$Te using spherical neutron polarimetry and scanning tunnelling microscopy to search for this canting in the bulk phase and compare surface and bulk magnetism. The results show that the bulk magnetic structure of Fe$_{1.09}$Te is consistent with collinear in-plane order ($\theta=0$ with an error of $\sim$ 5$^{\circ}$). Comparison with scanning tunnelling microscopy on a series of Fe$_{1+x}$Te samples reveals that the surface exhibits a magnetic surface reconstruction with a canting angle of the spins of $\theta=29.8^{\circ}$. We suggest that this is a consequence of structural relaxation of the surface layer resulting in an out-of-plane magnetocrystalline anisotropy. The magnetism in Fe$_{1+x}$Te displays different properties at the surface when the symmetry constraints of the bulk are removed.


title: A new parsimonious method for classifying Cancer Tissue-of-Origin Based on DNA Methylation 450K data

abstract: DNA methylation is a well-studied genetic modification that regulates gene transcription of Eukaryotes. Its alternations have been recognized as a significant component of cancer development. In this study, we use the DNA methylation 450k data from The Cancer Genome Atlas to evaluate the efficacy of DNA methylation data on cancer classification for 30 cancer types. We propose a new method for gene selection in high dimensional data(over 450 thousand). Variance filtering is first introduced for dimension reduction and Recursive feature elimination (RFE) is then used for feature selection. We address the problem of selecting a small subsets of genes from large number of methylated sites, and our parsimonious model is demonstrated to be efficient, achieving an accuracy over 91%, outperforming other studies which use DNA micro-arrays and RNA-seq Data . The performance of 20 models, which are based on 4 estimators (Random Forest, Decision Tree, Extra Tree and Support Vector Machine) and 5 classifiers (k-Nearest Neighbours, Support Vector Machine, XGboost, Light GBM and Multi-Layer Perceptron), is compared and robustness of the RFE algorithm is examined. Results suggest that the combined model of extra tree plus catboost classifier offers the best performance in cancer identification, with an overall validation accuracy of 91% , 92.3%, 93.3% and 93.5% for 20, 30, 40 and 50 features respectively. The biological functions in cancer development of 50 selected genes is also explored through enrichment analysis and the results show that 12 out of 16 of our top features have already been identified to be specific with cancer and we also propose some more genes to be tested for future studies. Therefore, our method may be utilzed as an auxiliary diagnostic method to determine the actual clinicopathological status of a specific cancer.


title: Harnessing confinement and driving to tune active particle dynamics

abstract: A distinguishing feature of active particles is the nature of the non-equilibrium noise driving their dynamics. Control of these noise properties is, therefore, of both fundamental and applied interest. We demonstrate emergent tuning of the active noise of a granular self-propelled particle by confining it to a quasi one-dimensional channel. We find that this particle, moving like an active Brownian particle (ABP) in two-dimensions, displays run-and-tumble (RTP) characteristics in confinement. We show that the dynamics of the relative orientation co-ordinate of the particle maps to that of a Brownian particle in a periodic potential subject to a constant force, in analogy to the dynamics of a molecular motor. This mapping captures the essential statistical characteristics of the one-dimensional RTP motion. Specifically, our theoretical analysis is in agreement with the empirical distributions of the relative orientation co-ordinate and the run-times (tumble-rates) of the particle. Finally, we explicitly control these emergent run-and-tumble like noise parameters by external driving. Altogether, our work illustrates geometry-induced tuning of the active dynamics of self-propelled units thus suggesting an independent route to harness their internal dynamics.


title: Eigenvalues of stochastic Hamiltonian systems with boundary conditions and its application

abstract: In this paper we solve the eigenvalue problem of stochastic Hamiltonian system with boundary conditions. Firstly, we extend the results in S. Peng \cite{peng} from time-invariant case to time-dependent case, proving the existence of a series of eigenvalues $\{\lambda_m\}$ and construct corresponding eigenfunctions. Moreover, the order of growth for these $\{\lambda_m\}$ are obtained: $\lambda_m\sim m^2$, as $m\rightarrow+\infty$. As applications, we give an explicit estimation formula about the statistic period of solutions of Forward-Backward SDEs. Besides, by a meticulous example we show the subtle situation in time-dependent case that some eigenvalues appear when the solution of the associated Riccati equation does not blow-up, which does not happen in time-invariant case.


title: WiMesh: Leveraging Mesh Networking For Disaster Communication in Poor Regions of the World

abstract: This paper discusses the design, implementation and field trials of WiMesh - a resilient Wireless Mesh Network (WMN) based disaster communication system purpose-built for underdeveloped and rural parts of the world. Mesh networking is a mature area, and the focus of this paper is not on proposing novel models, protocols or other mesh solutions. Instead, the paper focuses on the identification of important design considerations and justifications for several design trade offs in the context of mesh networking for disaster communication in developing countries with very limited resources. These trade-offs are discussed in the context of key desirable traits including security, low cost, low power, size, availability, customization, portability, ease of installation and deployment, and coverage area among others. We discuss at length the design, implementation, and field trial results of the WiMesh system which enables users spread over large geographical regions, to communicate with each other despite the lack of cellular coverage, power, and other communication infrastructure by leveraging multi-hop mesh networking and Wi-Fi equipped handheld devices. Lessons learned along with real-world results are shared for WiMesh deployment in a remote rural mountainous village of Pakistan, and the source code is shared with the research community.


title: StarNet: Gradient-free Training of Deep Generative Models using Determined System of Linear Equations

abstract: In this paper we present an approach for training deep generative models solely based on solving determined systems of linear equations. A network that uses this approach, called a StarNet, has the following desirable properties: 1) training requires no gradient as solution to the system of linear equations is not stochastic, 2) is highly scalable when solving the system of linear equations w.r.t the latent codes, and similarly for the parameters of the model, and 3) it gives desirable least-square bounds for the estimation of latent codes and network parameters within each layer.