title: Finite-size scaling of typicality-based estimates

abstract: According to the concept of typicality, an ensemble average can be accurately approximated by an expectation value with respect to a single pure state drawn at random from a high-dimensional Hilbert space. This random-vector approximation, or trace estimator, provides a powerful approach to, e.g., thermodynamic quantities for systems with large Hilbert-space sizes, which usually cannot be treated exactly, analytically or numerically. Here, we discuss the finite-size scaling of the accuracy of such trace estimators from two perspectives. First, we study the full probability distribution of random-vector expectation values and, second, the full temperature dependence of the standard deviation. With the help of numerical examples, we find pronounced Gaussian probability distributions and the expected decrease of the standard deviation with system size, at least above certain system-specific temperatures. Below and in particular for temperatures smaller than the excitation gap, simple rules are not available.


title: Combating False Negatives in Adversarial Imitation Learning

abstract: In adversarial imitation learning, a discriminator is trained to differentiate agent episodes from expert demonstrations representing the desired behavior. However, as the trained policy learns to be more successful, the negative examples (the ones produced by the agent) become increasingly similar to expert ones. Despite the fact that the task is successfully accomplished in some of the agent's trajectories, the discriminator is trained to output low values for them. We hypothesize that this inconsistent training signal for the discriminator can impede its learning, and consequently leads to worse overall performance of the agent. We show experimental evidence for this hypothesis and that the 'False Negatives' (i.e. successful agent episodes) significantly hinder adversarial imitation learning, which is the first contribution of this paper. Then, we propose a method to alleviate the impact of false negatives and test it on the BabyAI environment. This method consistently improves sample efficiency over the baselines by at least an order of magnitude.


title: Fast Generating A Large Number of Gumbel-Max Variables

abstract: The well-known Gumbel-Max Trick for sampling elements from a categorical distribution (or more generally a nonnegative vector) and its variants have been widely used in areas such as machine learning and information retrieval. To sample a random element $i$ (or a Gumbel-Max variable $i$) in proportion to its positive weight $v_i$, the Gumbel-Max Trick first computes a Gumbel random variable $g_i$ for each positive weight element $i$, and then samples the element $i$ with the largest value of $g_i+\ln v_i$. Recently, applications including similarity estimation and graph embedding require to generate $k$ independent Gumbel-Max variables from high dimensional vectors. However, it is computationally expensive for a large $k$ (e.g., hundreds or even thousands) when using the traditional Gumbel-Max Trick. To solve this problem, we propose a novel algorithm, \emph{FastGM}, that reduces the time complexity from $O(kn^+)$ to $O(k \ln k + n^+)$, where $n^+$ is the number of positive elements in the vector of interest. Instead of computing $k$ independent Gumbel random variables directly, we find that there exists a technique to generate these variables in descending order. Using this technique, our method FastGM computes variables $g_i+\ln v_i$ for all positive elements $i$ in descending order. As a result, FastGM significantly reduces the computation time because we can stop the procedure of Gumbel random variables computing for many elements especially for those with small weights. Experiments on a variety of real-world datasets show that FastGM is orders of magnitude faster than state-of-the-art methods without sacrificing accuracy and incurring additional expenses.


title: Approximation by periodic multivariate quasi-projection operators

abstract: Approximation properties of periodic quasi-projection operators with matrix dilations are studied. Such operators are generated by a sequence of functions $\varphi_j$ and a sequence of distributions/functions $\widetilde{\varphi}_j$. Error estimates for sampling-type quasi-projection operators are obtained under the periodic Strang-Fix conditions for $\varphi_j$ and the compatibility conditions for $\varphi_j$ and $\widetilde{\varphi}_j$. These estimates are given in terms of the Fourier coefficients of approximated functions and provide analogs of some known non-periodic results. Under some additional assumptions error estimates are given in other terms in particular using the best approximation. A number of examples are provided.


title: Closed geodesics on reversible Finsler 2-spheres

abstract: We extend two celebrated theorems on closed geodesics of Riemannian 2-spheres to the larger class of reversible Finsler 2-spheres: Lusternik-Schnirelmann's theorem asserting the existence of three simple closed geodesics, and Bangert-Franks-Hingston's theorem asserting the existence of infinitely many closed geodesics. In order to prove the first theorem, we employ the generalization of Grayson's curve shortening flow developed by Angenent-Oaks.


title: Using Machine Learning for Model Physics: an Overview

abstract: In the overview, a generic mathematical object (mapping) is introduced, and its relation to model physics parameterization is explained. Machine learning (ML) tools that can be used to emulate and/or approximate mappings are introduced. Applications of ML to emulate existing parameterizations, to develop new parameterizations, to ensure physical constraints, and control the accuracy of developed applications are described. Some ML approaches that allow developers to go beyond the standard parameterization paradigm are discussed.


title: WaveTTS: Tacotron-based TTS with Joint Time-Frequency Domain Loss

abstract: Tacotron-based text-to-speech (TTS) systems directly synthesize speech from text input. Such frameworks typically consist of a feature prediction network that maps character sequences to frequency-domain acoustic features, followed by a waveform reconstruction algorithm or a neural vocoder that generates the time-domain waveform from acoustic features. As the loss function is usually calculated only for frequency-domain acoustic features, that doesn't directly control the quality of the generated time-domain waveform. To address this problem, we propose a new training scheme for Tacotron-based TTS, referred to as WaveTTS, that has 2 loss functions: 1) time-domain loss, denoted as the waveform loss, that measures the distortion between the natural and generated waveform; and 2) frequency-domain loss, that measures the Mel-scale acoustic feature loss between the natural and generated acoustic features. WaveTTS ensures both the quality of the acoustic features and the resulting speech waveform. To our best knowledge, this is the first implementation of Tacotron with joint time-frequency domain loss. Experimental results show that the proposed framework outperforms the baselines and achieves high-quality synthesized speech.


title: A Time Delay Dynamical Model for Outbreak of 2019-nCoV and the Parameter Identification

abstract: In this paper, we propose a novel dynamical system with time delay to describe the outbreak of 2019-nCoV in China. One typical feature of this epidemic is that it can spread in latent period, which is therefore described by the time delay process in the differential equations. The accumulated numbers of classified populations are employed as variables, which is consistent with the official data and facilitates the parameter identification. The numerical methods for the prediction of outbreak of 2019-nCoV and parameter identification are provided, and the numerical results show that the novel dynamic system can well predict the outbreak trend so far. Based on the numerical simulations, we suggest that the transmission of individuals should be greatly controlled with high isolation rate by the government.


title: Dynamic Frame Structure for Next Generation Wireless Networks

abstract: In this paper, we devise a novel radio resource block (RB) structure named dynamic resource block structure (D-RBS) which can handle low latency traffics and large fluctuations in data rates by exploiting smart time and frequency duplexing. In our framework, the main resource block with a predefined bandwidth and time duration is divided into several small blocks with the same bandwidth and time duration. Depending on the service requirements, e.g., data rate and latency, the users are assigned to some these small blocks which could be noncontiguous both in frequency and time. This is in contrast to the previously introduced static resource block structure (S-RBS) where the size of each RB is predetermined and fixed. We provide resource allocation frameworks for this RB structure and formulate the optimization problems whose solutions are obtained by alternate search method (ASM) based on successive convex approximation approach (SCA). We provide a global optimal solution by exploiting the monotonic optimization method. By simulation we study the performance of our proposed scheme with S-RBS scheme and show it has 26% gain compared to the S-RBS scheme.


title: Long-distance transmission of quantum key distribution coexisting with classical optical communication over weakly-coupled few-mode fiber

abstract: Quantum key distribution (QKD) is one of the most practical applications in quantum information processing, which can generate information-theoretical secure keys between remote parties. With the help of the wavelength-division multiplexing technique, QKD has been integrated with the classical optical communication networks. The wavelength-division multiplexing can be further improved by the mode-wavelength dual multiplexing technique with few-mode fiber (FMF), which has additional modal isolation and large effective core area of mode, and particularly is practical in fabrication and splicing technology compared with the multi-core fiber. Here, we present for the first time a QKD implementation coexisting with classical optical communication over weakly-coupled FMF using all-fiber mode-selective couplers. The co-propagation of QKD with one 100 Gbps classical data channel at -2.60 dBm launched power is achieved over 86 km FMF with 1.3 kbps real-time secure key generation. Compared with single-mode fiber, the average Raman noise in FMF is reduced by 86% at the same fiber-input power. Our work implements an important approach to the integration between QKD and classical optical communication and previews the compatibility of quantum communications with the next-generation mode division multiplexing networks


title: Choice Set Optimization Under Discrete Choice Models of Group Decisions

abstract: The way that people make choices or exhibit preferences can be strongly affected by the set of available alternatives, often called the choice set. Furthermore, there are usually heterogeneous preferences, either at an individual level within small groups or within sub-populations of large groups. Given the availability of choice data, there are now many models that capture this behavior in order to make effective predictions--however, there is little work in understanding how directly changing the choice set can be used to influence the preferences of a collection of decision-makers. Here, we use discrete choice modeling to develop an optimization framework of such interventions for several problems of group influence, namely maximizing agreement or disagreement and promoting a particular choice. We show that these problems are NP-hard in general, but imposing restrictions reveals a fundamental boundary: promoting a choice can be easier than encouraging consensus or sowing discord. We design approximation algorithms for the hard problems and show that they work well on real-world choice data.


title: Gap opening in the spectrum of some Dirac-like pseudo-differential operators

abstract: In this paper, we study the opening of a spectral gap for a class of 2-dimensional periodic Hamiltonians which include those modelling multilayer graphene. The kinetic part of the Hamiltonian is given by $\sigma \cdot F(-i\nabla)$, where $\sigma$ denotes the Pauli matrices and $F$ is a sufficiently regular vector-valued function which equals 0 at the origin and grows at infinity. Its spectrum is the whole real line. We prove that a gap appears for perturbations in a certain class of periodic matrix-valued potentials depending on $F$, and we study how this gap depends on different parameters.


title: An Experimental Study of Formula Embeddings for Automated Theorem Proving in First-Order Logic

abstract: Automated theorem proving in first-order logic is an active research area which is successfully supported by machine learning. While there have been various proposals for encoding logical formulas into numerical vectors -- from simple strings to more involved graph-based embeddings -- little is known about how these different encodings compare. In this paper, we study and experimentally compare pattern-based embeddings that are applied in current systems with popular graph-based encodings, most of which have not been considered in the theorem proving context before. Our experiments show that the advantages of simpler encoding schemes in terms of runtime are outdone by more complex graph-based embeddings, which yield more efficient search strategies and simpler proofs. To support this, we present a detailed analysis across several dimensions of theorem prover performance beyond just proof completion rate, thus providing empirical evidence to help guide future research on neural-guided theorem proving towards the most promising directions.


title: Hyperon-nuclear interactions from SU(3) chiral effective field theory

abstract: The interaction between hyperons and nucleons has a wide range of applications in strangeness nuclear physics and is a topic of continuing great interest. These interactions are not only important for hyperon-nucleon scattering but also essential as basic input to studies of hyperon-nuclear few- and many-body systems including hypernuclei and neutron star matter. We review the systematic derivation and construction of such baryonic forces from the symmetries of quantum chromodynamics within non-relativistic SU(3) chiral effective field theory. Several applications of the resulting potentials are presented for topics of current interest in strangeness nuclear physics.


title: Condensed Generalized Finite Element Method (CGFEM)

abstract: Generalized or extended finite element methods (GFEM/XFEM) are in general badly conditioned and have numerous additional degrees of freedom (DOF) compared with the FEM because of introduction of enriched functions. In this paper, we develop an approach to establish a subspace of a conventional GFEM/XFEM approximation space using partition of unity (PU) techniques and local least square procedures. The proposed GFEM is referred to as condensed GFEM (CGFEM), which (i) possesses as many DOFs as the preliminary FEM, (ii) enjoys similar approximation properties with the GFEM/XFEM, and (iii) is well-conditioned in a sense that its conditioning is of the same order as that of the FEM. The fundamental approximation properties of CGFEM is proven mathematically. The CGFEM is applied to a problem of high order polynomial approximations and a Poisson crack problem; optimal convergence orders of the former are proven rigorously. The numerical experiments and comparisons with the conventional GFEM/XFEM and FEM are made to verify the theory and effectiveness of CGFEM.


title: A Simple Prediction Model for the Development Trend of 2019-nCov Epidemics Based on Medical Observations

abstract: In order to predict the development trend of the 2019 coronavirus (2019-nCov), we established an prediction model to predict the number of diagnoses case in China except Hubei Province. From January 25 to January 29, 2020, we optimized 6 prediction models, 5 of them based on the number of medical observations to predicts the peak time of confirmed diagnosis will appear on the period of morning of January 29 from 24:00 to February 2 before 5 o'clock 24:00. Then we tracked the data from 24 o'clock on January 29 to 24 o'clock on January 31, and found that the predicted value of the data on the 3rd has a small deviation from the actual value, and the actual value has always remained within the range predicted by the comprehensive prediction model 6. Therefore we discloses this finding and will continue to track whether this pattern can be maintained for longer. We believe that the changes medical observation case number may help to judge the trend of the epidemic situation in advance.


title: Dynamic maintenance policy for systems with repairable components subject to mutually dependent competing failure processes

abstract: In this paper, a repairable multi-component system is studied where all the components can be repaired individually within the system. The whole system is inspected at inspection intervals and the failed components are detected and replaced with a new one, while the other components continue functioning. Replacing components individually within the system makes their initial age to be different at each inspection time. Different initial age of all the components have affect on the system reliability and probability of failure and consequently the optimal inspection time, which is for the whole system not individual components. A dynamic maintenance policy is proposed to find the next inspection time based on the initial age of all the components. Two competing failure processes of degradation and a shock process are considered for each component. In our paper, there is a mutual dependency between the degradation process and shock process. Each incoming shock adds additional abrupt damages on the cumulative degradation path of all the components, and the shock arrival process is affected by the system degradation process. A realistic numerical example is presented to illustrate the proposed reliability and maintenance model.


title: Formulation for renormalon-free perturbative predictions beyond large-$\beta_0$ approximation

abstract: We present a formulation to give renormalon-free predictions consistently with fixed order perturbative results. The formulation has a similarity to Lee's method in that the renormalon-free part consists of two parts: one is given by a series expansion which does not contain renormalons and the other is the real part of the Borel integral of a singular Borel transform. The main novel aspect is to evaluate the latter object using a dispersion relation technique, which was possible only in the large-$\beta_0$ approximation. Here, we introduce an "ambiguity function," which is defined by inverse Mellin transform of the singular Borel transform. With the ambiguity function, we can rewrite the Borel integral in an alternative formula, which allows us to obtain the real part using analytic techniques similarly to the case of the large-$\beta_0$ approximation. We also present detailed studies of renormalization group properties of the formulation. As an example, we apply our formulation to the fixed-order result of the static QCD potential, whose closest renormalon is already visible.


title: Uncertainty Weighted Causal Graphs

abstract: Causality has traditionally been a scientific way to generate knowledge by relating causes to effects. From an imaginery point of view, causal graphs are a helpful tool for representing and infering new causal information. In previous works, we have generated automatically causal graphs associated to a given concept by analyzing sets of documents and extracting and representing the found causal information in that visual way. The retrieved information shows that causality is frequently imperfect rather than exact, feature gathered by the graph. In this work we will attempt to go a step further modelling the uncertainty in the graph through probabilistic improving the management of the imprecision in the quoted graph.


title: Assessing and Minimizing Collisions in Satellite Mega-Constellations

abstract: We aim to provide satellite operators and researchers with an efficient means for evaluating and mitigating collision risk during the design process of mega-constellations. We first establish a baseline for evaluating various techniques for close-encounter prediction and collision-probability calculation (Hoots et al. 1984, Gronchi 2005, JeongAhn and Malhotra 2015) by carrying out brute-force numerical simulations and using a sequence of filters to greatly reduce the computational expense of the algorithm. Next, we estimate conjunction events in the orbital environment following the anticipated deployments of the OneWeb LEO and SpaceX Starlink mega-constellations. As a final step, we investigate Minimum Space Occupancy (MiSO) orbits (Bombardelli et al. 2018), a generalization of the well-known frozen orbits that account for the perturbed-Keplerian dynamics of the Earth-Moon-Sun-satellite system. We evaluate the ability of MiSO configurations of the proposed mega-constellations, as suggested by Bombardelli et al. 2018, to reduce the risk of endogenous (intra-constellation) collisions. The results indicate that the adoption of the MiSO orbital configuration can significantly reduce risk with nearly indistinguishable adjustments to the nominal orbital elements of the constellation satellites.


title: Integral p-adic Hodge filtrations in low dimension and ramification

abstract: Given an integral p-adic variety, we observe that if the integral Hodge--de Rham spectral sequence behaves nicely, then the special fiber knows the Hodge numbers of the generic fiber. Applying recent advancements of integral p-adic Hodge theory, we show that such a nice behavior is guaranteed if the p-adic variety can be lifted to an analogue of second Witt vectors and satisfies some bound on dimension and ramification index. This is a (ramified) mixed characteristic analogue of results due to Deligne--Illusie and Fontaine--Messing. Lastly, we discuss an example illustrating the necessity of the aforementioned lifting condition, which is of independent interest.


title: The KdV soliton crosses a dissipative and dispersive border

abstract: We demonstrate the behavior of the soliton which, while moving in non-dissipative and dispersion-constant medium encounters a finite-width barrier with varying dissipation and/or dispersion; beyond the layer dispersion is constant (but not necessarily of the same value) and dissipation is null. The passed wave either retains the form of a soliton (though of different parameters) or becomes a bi-soliton. And a reflection wave may be negligible or absent. This models a situation similar to a light passing from a humid air to a dry one through the vapour saturation/condensation area. Some rough estimations for a prediction of an output are given using relative decay of the KdV conserved quantities are given. Keywords: KdV- Burgers, non-homogeneous layered media, soliton, bi-soliton, reflection, refraction.


title: On badly approximable numbers

abstract: Motivated by a wonderful paper by Ngoc Ai Van Nguyen, Anthony Po\"els and Damien Roy, where a powerful method was introduced, we prove a criterion for a vector $\pmb{\alpha}\in \mathbb{R}^d$ to be a badly approximable vector. Moreover we construct certain examples which show that a more general version of our criterion is not valid.


title: Integrating Deep Reinforcement Learning with Model-based Path Planners for Automated Driving

abstract: Automated driving in urban settings is challenging. Human participant behavior is difficult to model, and conventional, rule-based Automated Driving Systems (ADSs) tend to fail when they face unmodeled dynamics. On the other hand, the more recent, end-to-end Deep Reinforcement Learning (DRL) based model-free ADSs have shown promising results. However, pure learning-based approaches lack the hard-coded safety measures of model-based controllers. Here we propose a hybrid approach for integrating a path planning pipe into a vision based DRL framework to alleviate the shortcomings of both worlds. In summary, the DRL agent is trained to follow the path planner's waypoints as close as possible. The agent learns this policy by interacting with the environment. The reward function contains two major terms: the penalty of straying away from the path planner and the penalty of having a collision. The latter has precedence in the form of having a significantly greater numerical value. Experimental results show that the proposed method can plan its path and navigate between randomly chosen origin-destination points in CARLA, a dynamic urban simulation environment. Our code is open-source and available online.


title: Parafermionic bases of standard modules for affine Lie algebras

abstract: In this paper we construct combinatorial bases of parafermionic spaces associated with the standard modules of the rectangular highest weights for the untwisted affine Lie algebras. Our construction is a modification of G. Georgiev's construction for the affine Lie algebra $\widehat{\mathfrak sl}(n+1,\mathbb C)$---the constructed parafermionic bases are projections of the quasi-particle bases of the principal subspaces, obtained previously in a series of papers by the first two authors. As a consequence we prove the character formula of A. Kuniba, T. Nakanishi and J. Suzuki for all non-simply-laced untwisted affine Lie algebras.


title: On the structural properties of the bounded control set of a linear control system

abstract: The present paper shows that the bounded control set of a linear system on a connected Lie group $G$ contains all the bounded orbits of the system. As a consequence, its closure is the continuous image of the cartesian product of the set of control functions by the central subgroup associated with the drift of the system.


title: Duel of cosmological screening lengths

abstract: Two distinct perturbative approaches have been recently formulated within General Relativity, arguing for the screening of gravity in the $\Lambda$CDM Universe. In this paper we compare them and show that the offered screening concepts, each characterized by its own interaction range, can peacefully coexist. Accordingly, we advance a united scheme, determining the gravitational potential at all scales, including regions of nonlinear density contrasts, by means of a simple Helmholtz equation with the effective cosmological screening length. In addition, we claim that cosmic structures may not grow at distances above this Yukawa range and confront its current value with dimensions of the largest known objects in the Universe.


title: Intermittent dynamics of antiferromagnetic phase in inhomogeneous iron-based chalcogenide superconductor

abstract: Coexistence of phases, characterized by different electronic degrees of freedom, commonly occurs in layered superconductors. Among them, alkaline intercalated chalcogenides are model systems showing microscale coexistence of paramagnetic (PAR) and antiferromagnetic (AFM) phases, however, temporal behavior of different phases is still unknown. Here, we report the first visualization of the atomic motion in the granular phase of K$_{x}$Fe$_{2-y}$Se$_2$ using X-ray photon correlation spectroscopy. Unlike the PAR phase, the AFM texture reveals an intermittent dynamics with avalanches as in martensites. When cooled down across the superconducting transition temperature T$_c$, the AFM phase goes through an anomalous slowing behavior suggesting a direct relationship between the atomic motions in the AFM phase and the superconductivity. In addition of providing a compelling evidence of avalanche-like dynamics in a layered superconductor, the results provide a basis for new theoretical models to describe quantum states in inhomogeneous solids.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el


title: A millimeter-wave Bell Test using a ferrite parametric amplifier and a homodyne interferometer

abstract: A combined ferrite parametric amplifier and millimeter-wave homodyne interferometer are proposed as an ambient temperature Bell Test. It is shown that the non-linear magnetic susceptibility of the yttrium iron garnet (YIG) ferrite, on account of its narrow line-width Larmor precessional resonance, make it an ideal material for the creation of entangled photons. These can be measured using a homodyne interferometer, as the much larger number of thermally generated photons associated with ambient temperature emission can be screened out. The proposed architecture may enable YIG quantum technology-based sensors to be developed, mimicking in the millimeter-wave band the large number of quantum optical experiments in the near-infrared and visible regions which had been made possible by use of the nonlinear beta barium borate ferroelectric, an analogue of YIG. It is illustrated here how the YIG parametric amplifier can reproduce quantum optical Type I and Type II wave interactions, which can be used to create entangled photons in the millimeter-wave band. It is estimated that when half a cubic centimeter of YIG crystal is placed in a magnetic field of a few Tesla and pumped with 5 Watts of millimeter-wave radiation, approximately 0.5x10^12 entangled millimeter-wave photon pairs per second are generated by the spin-wave interaction. This means an integration time of only a few tens of seconds is needed for a successful Bell Test. A successful demonstration of this will lead to novel architectures of entanglement-based quantum technology room temperature sensors, re-envisioning YIG as a modern quantum material.


title: Simultaneous Left Atrium Anatomy and Scar Segmentations via Deep Learning in Multiview Information with Attention

abstract: Three-dimensional late gadolinium enhanced (LGE) cardiac MR (CMR) of left atrial scar in patients with atrial fibrillation (AF) has recently emerged as a promising technique to stratify patients, to guide ablation therapy and to predict treatment success. This requires a segmentation of the high intensity scar tissue and also a segmentation of the left atrium (LA) anatomy, the latter usually being derived from a separate bright-blood acquisition. Performing both segmentations automatically from a single 3D LGE CMR acquisition would eliminate the need for an additional acquisition and avoid subsequent registration issues. In this paper, we propose a joint segmentation method based on multiview two-task (MVTT) recursive attention model working directly on 3D LGE CMR images to segment the LA (and proximal pulmonary veins) and to delineate the scar on the same dataset. Using our MVTT recursive attention model, both the LA anatomy and scar can be segmented accurately (mean Dice score of 93% for the LA anatomy and 87% for the scar segmentations) and efficiently (~0.27 seconds to simultaneously segment the LA anatomy and scars directly from the 3D LGE CMR dataset with 60-68 2D slices). Compared to conventional unsupervised learning and other state-of-the-art deep learning based methods, the proposed MVTT model achieved excellent results, leading to an automatic generation of a patient-specific anatomical model combined with scar segmentation for patients in AF.


title: Don't Forget to Lock the Front Door! Inferring the Deployment of Source Address Validation of Inbound Traffic

abstract: This paper concerns the problem of the absence of ingress filtering at the network edge, one of the main causes of important network security issues. Numerous network operators do not deploy the best current practice - Source Address Validation (SAV) that aims at mitigating these issues. We perform the first Internet-wide active measurement study to enumerate networks not filtering incoming packets by their source address. The measurement method consists of identifying closed and open DNS resolvers handling requests coming from the outside of the network with the source address from the range assigned inside the network under the test. The proposed method provides the most complete picture of the inbound SAV deployment state at network providers. We reveal that 32 673 Autonomous Systems (ASes) and 197 641 Border Gateway Protocol (BGP) prefixes are vulnerable to spoofing of inbound traffic. Finally, using the data from the Spoofer project and performing an open resolver scan, we compare the filtering policies in both directions.


title: New moduli spaces of one-dimensional sheaves on $\mathbb{P}^3$

abstract: We define a one-dimensional family of "Euler" stability conditions on $\mathbb{P}^n$ which are conjectured to converge to Gieseker stability for coherent sheaves. Here, we focus on ${\mathbb P}^3$, first identifying Euler stability conditions with double-tilt stability conditions, and then we consider moduli of one-dimensional sheaves, proving some asymptotic results, boundedness for walls, and then explicitly computing walls and wall-crossings for sheaves supported on rational curves of degrees $3$ and $4$.


title: Hong-Ou-Mandel interference on a real beam splitter

abstract: Hong-Ou-Mandel (HOM) effect is known to be one of the main phenomena in quantum optics. The effect occurs when two identical single-photon waves enter a 1:1 beam splitter, one in each input port. When the photons are identical, they will extinguish each other. One of the main elements of the HOM interferometer is the beam splitter, which has its own coefficients of reflection $R = 1/2$ and transmission $ T = 1/2 $. In this work, the general mechanism of the interaction of two photons in a beam splitter is considered using an analytical solution, which shows that in the HOM theory of the effect it is necessary to know (including when planning the experiment) not only $ R = 1/2 $ and $ T = 1/2 $, but also their root-mean-square fluctuations $ \Delta R ^ 2, \Delta T ^ 2 $, which arise due to the dependence of $R = R(\omega_1, \omega_2) $ and $ T = T (\omega_1, \omega_2) $ on the frequencies where $\omega_1, \omega_2$ are the frequencies of the first and second photons, respectively. Under certain conditions, specifically when the dependence of the fluctuations $ \Delta R^2 $ and $ \Delta T^2 $ can be neglected and $ R=T=1/2 $ is chosen, the developed theory coincides with previously known results.


title: Deep Reinforcement Learning for Autonomous Driving: A Survey

abstract: With the development of deep representation learning, the domain of reinforcement learning (RL) has become a powerful learning framework now capable of learning complex policies in high dimensional environments. This review summarises deep reinforcement learning (DRL) algorithms and provides a taxonomy of automated driving tasks where (D)RL methods have been employed, while addressing key computational challenges in real world deployment of autonomous driving agents. It also delineates adjacent domains such as behavior cloning, imitation learning, inverse reinforcement learning that are related but are not classical RL algorithms. The role of simulators in training agents, methods to validate, test and robustify existing solutions in RL are discussed.


title: Geometric Methods In Heegaard Theory

abstract: We survey some recent geometric methods for studying Heegaard splittings of 3-manifolds


title: Topological Valley Transport of Gapped Dirac Magnons in Bilayer Ferromagnetic Insulators

abstract: Bilayer Heisenberg ferromagnetic insulators hold degenerate terahertz Dirac magnon modes associated with two opposite valleys of the hexagonal Brillouin zone. We show that this energy degeneracy can be removed by breaking of the inversion symmetry (I), leading to a topological magnon valley current. We show furthermore that this current leads to valley Seebeck effect for magnons and is thereby detectable. We perform calculations in the specific example of bilayer CrBr3, where I can be broken by electrostatic doping.


title: On Sum-Of-Tails Identities

abstract: In this article, a finite analogue of the generalized sum-of-tails identity of Andrews and Freitas is obtained. We derive several interesting results as special cases of this analogue, in particular, a recent identity of Dixit, Eyyyunni, Maji and Sood. We derive a new extension of Abel's lemma with the help of which we obtain a one-parameter generalization of a sum-of-tails identity of Andrews, Garvan and Liang, an identity of Ramanujan as well as two new results - one for Ramanujan's function $\sigma(q)$ and another for the function recently introduced by Andrews and Ballantine. Later we introduce a new generalization $\mathrm{FFW}_{c}(n)$ of a function of Fokkink, Fokkink and Wang and derive an identity for its generating function. This gives, as a special case, a recent representation for the generating function of $\mathrm{spt}(n)$ given by Andrews, Garvan and Liang. We also obtain some weighted partition identities along with new representations for two of Ramanujan's third order mock theta functions through combinatorial techniques.


title: Stability of synchronous states in sparse neuronal networks

abstract: The stability of synchronous states is analysed in the context of two populations of inhibitory and excitatory neurons, characterized by different pulse-widths. The problem is reduced to that of determining the eigenvalues of a suitable class of sparse random matrices, randomness being a consequence of the network structure. A detailed analysis, which includes also the study of finite-amplitude perturbations, is performed in the limit of narrow pulses, finding that the stability depends crucially on the relative pulse-width. This has implications for the overall property of the asynchronous (balanced) regime.


title: Dynamic Set Values for Nonzero Sum Games with Multiple Equilibriums

abstract: Nonzero sum games typically have multiple Nash equilibriums (or no equilibrium), and unlike the zero sum case, they may have different values at different equilibriums. Instead of focusing on the existence of individual equilibriums, we study the set of values over all equilibriums, which we call the set value of the game. The set value is unique by nature and always exists (with possible value $\emptyset$). Similar to the standard value function in control literature, it enjoys many nice properties such as regularity, stability, and more importantly the dynamic programming principle. There are two main features in order to obtain the dynamic programming principle: (i) we must use closed-loop controls (instead of open-loop controls); (ii) we must allow for path dependent controls, even if the problem is in a state dependent (Markovian) setting. We shall consider both discrete and continuous time models with finite time horizon. For the latter we will also provide a duality approach through certain standard PDE (or path dependent PDE), which is quite efficient for numerically computing the set value of the game.


title: The empirical mean position of a branching L\'evy process

abstract: We consider a supercritical branching L\'evy process on the real line. Under mild moment assumptions on the number of offspring and their displacements, we prove a second-order limit theorem on the empirical mean position.


title: Fair Allocation Based Soft Load Shedding

abstract: Renewable sources are taking center stage in electricity generation. Due to the intermittent nature of these renewable resources, the problem of the demand-supply gap arises. To solve this problem, several techniques have been proposed in the literature in terms of cost (adding peaker plants), availability of data (Demand Side Management "DSM"), hardware infrastructure (appliance controlling DSM) and safety (voltage reduction). However, these solutions are not fair in terms of electricity distribution. In many cases, although the available supply may not match the demand in peak hours, however, the total aggregated demand remains less than the total supply for the whole day. Load shedding (complete blackout) is a commonly used solution to deal with the demand-supply gap, which can cause substantial economic losses. To solve the demand-supply gap problem, we propose a solution called Soft Load Shedding (SLS), which assigns electricity quota to each household in a fair way. We measure the fairness of SLS by defining a function for household satisfaction level. We model the household utilities by parametric function and formulate the problem of SLS as a social welfare problem. We also consider revenue generated from the fair allocation as a performance measure. To evaluate our approach, extensive experiments have been performed on both synthetic and real-world datasets, and our model is compared with several baselines to show its effectiveness in terms of fair allocation and revenue generation.


title: Core-shell microgels via precipitation polymerization: computer simulations

abstract: In this work we presented a novel computational model of precipitation polymerization allowing one to obtain core-shell microgels via a realistic cross-linking process based on the experimental procedure. We showed that the cross-linker-monomer reactivity ratios r are responsible for the microgel internal structure. Values of r lower than 1 correspond to the case when alternating sequences occur at the early reaction stages; this leads to the formation of microgels with pronounced core-shell structure. The distribution of dangling ends for small values of r becomes bimodal with two well-distinguished peaks, which correspond to the core (short dangling ends) and corona (long dangling ends) regions. The density profiles confirm the existence of two distinct regions for small r: a densely cross-linked core and a loose corona entirely consisting of dangling ends with no cross-linker. The consumption of the cross-linker in the course in the microgel formation was found to be in a perfect agreement with the predictions of Monte Carlo (MC) model in the sequence space.


title: DropClass and DropAdapt: Dropping classes for deep speaker representation learning

abstract: Many recent works on deep speaker embeddings train their feature extraction networks on large classification tasks, distinguishing between all speakers in a training set. Empirically, this has been shown to produce speaker-discriminative embeddings, even for unseen speakers. However, it is not clear that this is the optimal means of training embeddings that generalize well. This work proposes two approaches to learning embeddings, based on the notion of dropping classes during training. We demonstrate that both approaches can yield performance gains in speaker verification tasks. The first proposed method, DropClass, works via periodically dropping a random subset of classes from the training data and the output layer throughout training, resulting in a feature extractor trained on many different classification tasks. Combined with an additive angular margin loss, this method can yield a 7.9% relative improvement in equal error rate (EER) over a strong baseline on VoxCeleb. The second proposed method, DropAdapt, is a means of adapting a trained model to a set of enrolment speakers in an unsupervised manner. This is performed by fine-tuning a model on only those classes which produce high probability predictions when the enrolment speakers are used as input, again also dropping the relevant rows from the output layer. This method yields a large 13.2% relative improvement in EER on VoxCeleb. The code for this paper has been made publicly available.


title: Accelerated expansion of the universe from the perspective of inhomogeneous cosmology

abstract: We present a special case of the Stephani solution with spherical symmetry while considering different values of spatial curvature. We investigate the dynamics of the universe evolution in our model, build the R--T-regions for the resulting spacetime and analyze the behavior of the deceleration parameter. The singularities of the model are also discussed. The geometry of the spatial part of the obtained solution is explored.


title: Random walks on tori and normal numbers in self similar sets

abstract: We study random walks on a $d$-dimensional torus by affine expanding maps whose linear parts commute. Assuming an irrationality condition on their translation parts, we prove that the Haar measure is the unique stationary measure. We deduce that if $K \subset \mathbb{R}^d$ is an attractor of a finite iterated function system of $n\geq 2$ maps of the form $x \mapsto D^{-r_i} x + t_i \ (i=1, \ldots, n)$, where $D$ is an expanding $d\times d$ integer matrix, and is the same for all the maps, and $r_{i} \in\mathbb{N}$, under an irrationality condition on the translation parts $t_i$, almost every point in $K$ (w.r.t. any Bernoulli measure) has an equidistributed orbit under the map $x\mapsto Dx$ (multiplication mod $\mathbb{Z}^{d}$). In the one-dimensional case, this conclusion amounts to normality to base $D$. Thus for example, almost every point in an irrational dilation of the middle-thirds Cantor set is normal to base 3.


title: Permissioned Blockchain-Based Security for SDN in IoT Cloud Networks

abstract: The advancement in cloud networks has enabled connectivity of both traditional networked elements and new devices from all walks of life, thereby forming the Internet of Things (IoT). In an IoT setting, improving and scaling network components as well as reducing cost is essential to sustain exponential growth. In this domain, software-defined networking (SDN) is revolutionizing the network infrastructure with a new paradigm. SDN splits the control/routing logic from the data transfer/forwarding. This splitting causes many issues in SDN, such as vulnerabilities of DDoS attacks. Many solutions (including blockchain based) have been proposed to overcome these problems. In this work, we offer a blockchain-based solution that is provided in redundant SDN (load-balanced) to service millions of IoT devices. Blockchain is considered as tamper-proof and impossible to corrupt due to the replication of the ledger and consensus for verification and addition to the ledger. Therefore, it is a perfect fit for SDN in IoT Networks. Blockchain technology provides everyone with a working proof of decentralized trust. The experimental results show gain and efficiency with respect to the accuracy, update process, and bandwidth utilization.


title: On the classification of Smale-Barden manifolds with Sasakian structures

abstract: Smale-Barden manifolds $M$ are classified by their second homology $H_2(M,{\mathbb Z})$ and the Barden invariant $i(M)$. It is an important and dificult question to decide when $M$ admits a Sasakian structure in terms of these data. In this work we show methods of doing this. In particular we realize all $M$ with $H_2(M)={\mathbb Z}^k\oplus(\oplus_{i=1}^r{\mathbb Z}_{m_i}^{2g_i})$ and $i=0,\infty$, provided that $k\geq 1$, $m_i\geq 2$, $g_i\geq 1$, $m_i$ are pairwise coprime. Using our methods we also contribute to the problem of the existence of definite Sasakian structures on rational homology spheres. Also, we give a complete solution to the problem of the existence of Sasakian structures on rational homology spheres in the class of semi-regular Sasakian structures.


title: A Classification Framework for TDMA Scheduling techniques in WSNs

abstract: One of the major challenges in wireless sensor networks (WSNs) is the mitigation of collisions due to simultaneous transmissions by multiple nodes over a common channel which are located in a proximity. TDMA-based channel access provides energy-efficient and collision-free transmissions. It is especially suitable for traffic with periodic transmission patterns and guaranteed QoS requirements. For that reason, a large number of TDMA-scheduling algorithms are available in the literature, and consequently, a good number of survey papers on TDMA-scheduling algorithms have been written. In this work, we propose a novel classification framework to categorize the existing TDMA-scheduling algorithms available for WSNs. As against existing survey works, the proposed framework possess certain new dimensions (categories) to classify existing TDMA-scheduling algorithms. Additionally, we introduce a couple of new sub-categories for the existing classes which would help researchers to even differentiate between two TDMA-Scheduling algorithms that are assumed to be similar as per existing classification schemes. Finally, we also discuss few important works in the context of proposed classification scheme.


title: Annual modulations from secular variations: relaxing DAMA?

abstract: The DAMA collaboration reported an annually modulated rate with a phase compatible with a Dark Matter induced signal. We point out that a slowly varying rate can bias or even simulate an annual modulation if data are analyzed in terms of residuals computed by subtracting approximately yearly averages starting from a fixed date, rather than a background continuous in time. In the most extreme case, the amplitude and phase of the annual modulation reported by DAMA could be alternatively interpreted as a decennial growth of the rate. This possibility appears mildly disfavoured by a detailed study of the available data, but cannot be safely excluded. In general, a decreasing or increasing rate could partially reduce or enhance a true annual modulation, respectively. The issue could be clarified by looking at the full time-dependence of the DAMA total rate, not explicitly published so far.


title: Regularizing Reasons for Outfit Evaluation with Gradient Penalty

abstract: In this paper, we build an outfit evaluation system which provides feedbacks consisting of a judgment with a convincing explanation. The system is trained in a supervised manner which faithfully follows the domain knowledge in fashion. We create the EVALUATION3 dataset which is annotated with judgment, the decisive reason for the judgment, and all corresponding attributes (e.g. print, silhouette, and material \etc.). In the training process, features of all attributes in an outfit are first extracted and then concatenated as the input for the intra-factor compatibility net. Then, the inter-factor compatibility net is used to compute the loss for judgment. We penalize the gradient of judgment loss of so that our Grad-CAM-like reason is regularized to be consistent with the labeled reason. In inference, according to the obtained information of judgment, reason, and attributes, a user-friendly explanation sentence is generated by the pre-defined templates. The experimental results show that the obtained network combines the advantages of high precision and good interpretation.


title: Effect of Analysis Window and Feature Selection on Classification of Hand Movements Using EMG Signal

abstract: Electromyography (EMG) signals have been successfully employed for driving prosthetic limbs of a single or double degree of freedom. This principle works by using the amplitude of the EMG signals to decide between one or two simpler movements. This method underperforms as compare to the contemporary advances done at the mechanical, electronics, and robotics end, and it lacks intuition. Recently, research on myoelectric control based on pattern recognition (PR) shows promising results with the aid of machine learning classifiers. Using the approach termed as, EMG-PR, EMG signals are divided into analysis windows, and features are extracted for each window. These features are then fed to the machine learning classifiers as input. By offering multiple class movements and intuitive control, this method has the potential to power an amputated subject to perform everyday life movements. In this paper, we investigate the effect of the analysis window and feature selection on classification accuracy of different hand and wrist movements using time-domain features. We show that effective data preprocessing and optimum feature selection helps to improve the classification accuracy of hand movements. We use publicly available hand and wrist gesture dataset of $40$ intact subjects for experimentation. Results computed using different classification algorithms show that the proposed preprocessing and features selection outperforms the baseline and achieve up to $98\%$ classification accuracy.


title: Multi-Toeplitz operators associated with regular polydomains

abstract: In this paper we introduce and study the class of weighted multi-Toeplitz operators associated with noncommutative polydomains ${\bf D_f^m}$, ${\bf m}:=(m_1,\ldots, m_k)\in {\bf N}^k$, generated by $k$-tuples ${\bf f}:=(f_1,\ldots, f_k)$ of positive regular free holomorphic functions in a neighborhood of the origin. These operators are acting on the tensor product $F^2(H_{n_1})\otimes \cdots \otimes F^2(H_{n_k})$ of full Fock spaces with $n_i$ generators or, equivalently, they can be viewed as multi-Toeplitz operators acting on tensor products of weighted full Fock spaces. For a large class of polydomains, we show that there are no non-zero compact multi-Toeplitz operators. We characterize the weighted multi-Toeplitz operators in terms of bounded free $k$-pluriharmonic functions on the radial part of ${\bf D_f^m}$ and use the result to obtain an analogue of the Dirichlet extension problem for free $k$-pluriharmonic functions. We show that the weighted multi-Toeplitz operators have noncommutative Fourier representations which can be viewed as noncommutative symbols and can be used to recover the associated operators. We also prove that the weighted multi-Toeplitz operators satisfy a Brown-Halmos type equation associated with the polydomain ${\bf D_f^m}$.


title: Asymptotic spectra of large (grid) graphs with a uniform local structure

abstract: We are concerned with sequences of graphs having a grid geometry, with a uniform local structure in a bounded domain $\Omega\subset {\mathbb R}^d$, $d\ge 1$. We assume $\Omega$ to be Lebesgue measurable with regular boundary and contained, for convenience, in the cube $[0,1]^d$. When $\Omega=[0,1]$, such graphs include the standard Toeplitz graphs and, for $\Omega=[0,1]^d$, the considered class includes $d$-level Toeplitz graphs. In the general case, the underlying sequence of adjacency matrices has a canonical eigenvalue distribution, in the Weyl sense, and we show that we can associate to it a symbol $f$. The knowledge of the symbol and of its basic analytical features provide many informations on the eigenvalue structure, of localization, spectral gap, clustering, and distribution type. Few generalizations are also considered in connection with the notion of generalized locally Toeplitz sequences and applications are discussed, stemming e.g. from the approximation of differential operators via numerical schemes.


title: Full-Blind Delegating Private Quantum Computation

abstract: The delegating private quantum computation (DQC) protocol with the universal quantum gate set $\left\{ {X,Z,H,P,R,CNOT} \right\}$ was firstly proposed by Broadbent \emph{et al.}, and then Tan \emph{et al.} tried to put forward an half-blind DQC protocol (HDQC) with another universal set $\left\{ {H,P,CNOT,T} \right\}$. However, the decryption circuit of \emph{Toffoli} gate (i.e., \emph{T}) is a little redundant, and Tan \emph{et al}.'s protocol exists the information leak. In addition, both of these two protocols just focus on the blindness of data (i.e., the client's input and output), but do not consider the blindness of computation (i.e., the delegated quantum operation). For solving these problems, we propose a full-blind DQC protocol (FDQC) with quantum gate set $\left\{ {H,P,CNOT,T} \right\}$ , where the desirable delegated quantum operation, one of $\left\{ {H,P,CNOT,T} \right\}$ , is replaced by a fixed sequence $\left \{ {H,P,T,CZ,CNOT} \right\}$ to make the computation blind, and the decryption circuit of \emph{Toffoli} gate is also optimized. Analysis shows that our protocol can not only correctly perform any delegated quantum computation, but also holds the characteristics of data blindness and computation blindness.


title: On Smoothness of the Abel Equation Solution in Terms of the Jacoby Series Coefficients

abstract: In this paper we continue the investigation of the Abel equation with the right part belonging to a Lebesgue weighted space. We have improved the previously known result - the uniqueness and existence theorem formulated in terms of the Jacoby series coefficients that gives us an opportunity to find and classify a solution due to an asymptotic of some relation containing the Jacoby coefficients of the right part. The new main results are in the following: The conditions imposed on the parameters, under which the Abel equation has a unique solution represented by the series, are formulated; The relationship between the values of the parameters and the solution smoothness is established. The independence between one of the parameters and the smoothness of the solution is proved.


title: Hurwitz numbers from matrix integrals over Gaussian measure

abstract: We explain how Gaussian integrals over ensemble of complex matrices with source matrices generate Hurwitz numbers of the most general type, namely, Hurwitz numbers with arbitrary orientable or non-orientable base surface and arbitrary profiles at branch points. We use the Feynman diagram approach. The connections with topological theories and also with certain classical and quantum integrable theories in particular with Witten's description of two-dimensional gauge theory are shown.


title: Safe Exploration for Optimizing Contextual Bandits

abstract: Contextual bandit problems are a natural fit for many information retrieval tasks, such as learning to rank, text classification, recommendation, etc. However, existing learning methods for contextual bandit problems have one of two drawbacks: they either do not explore the space of all possible document rankings (i.e., actions) and, thus, may miss the optimal ranking, or they present suboptimal rankings to a user and, thus, may harm the user experience. We introduce a new learning method for contextual bandit problems, Safe Exploration Algorithm (SEA), which overcomes the above drawbacks. SEA starts by using a baseline (or production) ranking system (i.e., policy), which does not harm the user experience and, thus, is safe to execute, but has suboptimal performance and, thus, needs to be improved. Then SEA uses counterfactual learning to learn a new policy based on the behavior of the baseline policy. SEA also uses high-confidence off-policy evaluation to estimate the performance of the newly learned policy. Once the performance of the newly learned policy is at least as good as the performance of the baseline policy, SEA starts using the new policy to execute new actions, allowing it to actively explore favorable regions of the action space. This way, SEA never performs worse than the baseline policy and, thus, does not harm the user experience, while still exploring the action space and, thus, being able to find an optimal policy. Our experiments using text classification and document retrieval confirm the above by comparing SEA (and a boundless variant called BSEA) to online and offline learning methods for contextual bandit problems.


title: Zeolite-inspired 3d printed structures with enhanced mechanical properties

abstract: Specific strength (strength/density) is a crucial factor while designing high load bearing architecture in areas of aerospace and defence. Strength of the material can be enhanced by blending with high strength component or, by compositing with high strength fillers but both the options has limitations such as at certain load, materials fail due to poor filler and matrix interactions. Therefore, researchers are interested in enhancing strength of materials by playing with topology/geometry and therefore nature is best option to mimic for structures whereas, complexity limits nature mimicked structures. In this paper, we have explored Zeolite-inspired structures for load bearing capacity. Zeolite-inspired structure were obtained from molecular dynamics simulation and then fabricated via Fused deposition Modeling. The atomic scale complex topology from simulation is experimentally synthesized using 3D printing. Compressibility of as-fabricated structures was tested in different direction and compared with simulation results. Such complex architecture can be used for ultralight aerospace and automotive parts.


title: WeatherBench: A benchmark dataset for data-driven weather forecasting

abstract: Data-driven approaches, most prominently deep learning, have become powerful tools for prediction in many domains. A natural question to ask is whether data-driven methods could also be used to predict global weather patterns days in advance. First studies show promise but the lack of a common dataset and evaluation metrics make inter-comparison between studies difficult. Here we present a benchmark dataset for data-driven medium-range weather forecasting, a topic of high scientific interest for atmospheric and computer scientists alike. We provide data derived from the ERA5 archive that has been processed to facilitate the use in machine learning models. We propose simple and clear evaluation metrics which will enable a direct comparison between different methods. Further, we provide baseline scores from simple linear regression techniques, deep learning models, as well as purely physical forecasting models. The dataset is publicly available at https://github.com/pangeo-data/WeatherBench and the companion code is reproducible with tutorials for getting started. We hope that this dataset will accelerate research in data-driven weather forecasting.


title: Reduced-order modeling of advection-dominated systems with recurrent neural networks and convolutional autoencoders

abstract: A common strategy for the dimensionality reduction of nonlinear partial differential equations relies on the use of the proper orthogonal decomposition (POD) to identify a reduced subspace and the Galerkin projection for evolving dynamics in this reduced space. However, advection-dominated PDEs are represented poorly by this methodology since the process of truncation discards important interactions between higher-order modes during time evolution. In this study, we demonstrate that an encoding using convolutional autoencoders (CAEs) followed by a reduced-space time evolution by recurrent neural networks overcomes this limitation effectively. We demonstrate that a truncated system of only two latent-space dimensions can reproduce a sharp advecting shock profile for the viscous Burgers equation with very low viscosities, and a six-dimensional latent space can recreate the evolution of the inviscid shallow water equations. Additionally, the proposed framework is extended to a parametric reduced-order model by directly embedding parametric information into the latent space to detect trends in system evolution. Our results show that these advection-dominated systems are more amenable to low-dimensional encoding and time evolution by a CAE and recurrent neural network combination than the POD Galerkin technique.


title: Wavelength transduction from a 3D microwave cavity to telecom using piezoelectric optomechanical crystals

abstract: Microwave to optical transduction has received a great deal of interest from the cavity optomechanics community as a landmark application for electro-optomechanical systems. In this Letter, we demonstrate a novel transducer that combines high-frequency mechanical motion and a microwave cavity for the first time. The system consists of a 3D microwave cavity and a gallium arsenide optomechanical crystal, which has been placed in the microwave electric field maximum. This allows the microwave cavity to actuate the gigahertz-frequency mechanical breathing mode in the optomechanical crystal through the piezoelectric effect, which is then read out using a telecom optical mode. The gallium arsenide optomechanical crystal is a good candidate for low-noise microwave-to-telecom transduction, as it has been previously cooled to the mechanical ground state in a dilution refrigerator. Moreover, the 3D microwave cavity architecture can naturally be extended to couple to superconducting qubits and to create hybrid quantum systems.

categories: cond-mat.mes-hall physics.app-ph physics.optics quant-ph


title: One-thimble regularisation of lattice field theories: is it only a dream?

abstract: Lefschetz thimbles regularisation of (lattice) field theories was put forward as a possible solution to the sign problem. Despite elegant and conceptually simple, it has many subtleties, a major one boiling down to a plain question: how many thimbles should we take into account? In the original formulation, a single thimble dominance hypothesis was put forward: in the thermodynamic limit, universality arguments could support a scenario in which the dominant thimble (associated to the global minimum of the action) captures the physical content of the field theory. We know by now many counterexamples and we have been pursuing multi-thimble simulations ourselves. Still, a single thimble regularisation would be the real breakthrough. We report on ongoing work aiming at a single thimble formulation of lattice field theories, in particular putting forward the proposal of performing Taylor expansions on the dominant thimble.


title: METTEOR: Robust Multi-Traffic Topology Engineering for Commercial Data Center Networks

abstract: Numerous optical circuit switched data center networks have been proposed over the past decade for higher capacity, though commercial adoption of these architectures have been minimal so far. One major challenge commonly facing these architectures is the difficulty of handling bursty traffic with optical circuit switches (OCS) with high switching latency. Prior works generally rely on fast-switching OCS prototypes to better react to traffic changes via frequent reconfigurations. This approach, unfortunately, adds further complexity to the control plane. We propose METTEOR, an easily deployable solution for optical circuit switched data centers, that is designed for the current capabilities of commercial OCSs. Using multiple predicted traffic matrices, METTEOR designs data center topologies that are less sensitive to traffic changes, thus eliminating the need of frequently reconfiguring OCSs upon traffic changes. Results based on extensive evaluations using production traces show that METTEOR increases the percentage of direct-hop traffic by about 80% over a fat tree at comparable cost, and by about 30% over a uniform mesh, at comparable maximum link utilizations. Compared to ideal solutions that reconfigure OCSs on every traffic matrix, METTEOR achieves close-to-optimal bandwidth utilization even with biweekly reconfiguration. This drastically lowers the controller and management complexity needed to perform METTEOR in commercial settings.


title: Deep Learning Based Spatial User Mapping on Extra Large MIMO Arrays

abstract: In an extra-large scale MIMO (XL-MIMO) system, the antenna arrays have a large physical size that goes beyond the dimensions in traditional MIMO systems. Because of this large dimensionality, the optimization of an XL-MIMO system leads to solutions with prohibitive complexity when relying on conventional optimization tools. In this paper, we propose a design based on machine learning for the downlink of a multi-user setting with linear pre-processing, where the goal is to select a limited mapping area per user, i.e. a small portion of the array that contains the beamforming energy to the user. We refer to this selection as spatial user mapping (SUM). Our solution relies on learning using deep convolutional neural networks with a distributed architecture that is built to manage the large system dimension. This architecture contains one network per user where all the networks work in parallel and exploit specific non-stationary properties of the channels along the array. Our results show that, once the parallel networks are trained, they provide the optimal SUM solution in more than $80\%$ of the instances, resulting in a negligible sum-rate loss compared to a system using the optimal SUM solution while providing an insightful approach to rethink these kinds of problems that have no closed-form solution.


title: Phase-space resolved rates in driven multidimensional chemical reactions

abstract: Chemical reactions in multidimensional driven systems are typically described by a time-dependent rank-1 saddle associated with one reaction and several orthogonal coordinates (including the solvent bath). To investigate reactions in such systems, we develop a fast and robust method ---viz., local manifold analysis (LMA)--- for computing the instantaneous decay rate of reactants. Specifically, it computes the instantaneous decay rates along saddle-bound trajectories near the activated complex by exploiting local properties of the stable and unstable manifold associated with the normally hyperbolic invariant manifold (NHIM). The LMA method offers substantial reduction of numerical effort and increased reliability in comparison to direct ensemble integration. It provides an instantaneous flux that can be assigned to every point on the NHIM and which is associated with a trajectory ---regardless of whether it is periodic, quasi-periodic, or chaotic--- that is bound on the NHIM. The time average of these fluxes in the driven system corresponds to the average rate through a given local section containing the corresponding point on the NHIM. We find good agreement between the results of the LMA and direct ensemble integration obtained using numerically constructed, recrossing-free dividing surfaces.


title: Sound Event Detection with Depthwise Separable and Dilated Convolutions

abstract: State-of-the-art sound event detection (SED) methods usually employ a series of convolutional neural networks (CNNs) to extract useful features from the input audio signal, and then recurrent neural networks (RNNs) to model longer temporal context in the extracted features. The number of the channels of the CNNs and size of the weight matrices of the RNNs have a direct effect on the total amount of parameters of the SED method, which is to a couple of millions. Additionally, the usually long sequences that are used as an input to an SED method along with the employment of an RNN, introduce implications like increased training time, difficulty at gradient flow, and impeding the parallelization of the SED method. To tackle all these problems, we propose the replacement of the CNNs with depthwise separable convolutions and the replacement of the RNNs with dilated convolutions. We compare the proposed method to a baseline convolutional neural network on a SED task, and achieve a reduction of the amount of parameters by 85% and average training time per epoch by 78%, and an increase the average frame-wise F1 score and reduction of the average error rate by 4.6% and 3.8%, respectively.


title: Uniform full stability of recovering convolutional perturbation of the Sturm-Liouville operator from the spectrum

abstract: The perturbation of the Sturm-Liouville operator on a finite interval with Dirichlet boundary conditions by a convolution operator is considered. Local stability and global unique solvability of the inverse problem of recovering the convolution kernel from the spectrum, provided that the potential is given a priori, is known. In the present work, we establish uniform full stability of this inverse problem involving a uniform estimate of deviations of the convolution kernel via deviations of the spectrum and the potential within balls of any fixed radii.


title: Neural network approach for the dynamics on the normally hyperbolic invariant manifold of periodically driven systems

abstract: Chemical reactions in multidimensional systems are often described by a rank-1 saddle, whose stable and unstable manifolds intersect in the normally hyperbolic invariant manifold (NHIM). Trajectories started on the NHIM in principle never leave this manifold when propagated forward or backward in time. However, the numerical investigation of the dynamics on the NHIM is difficult because of the instability of the motion. We apply a neural network to describe time-dependent NHIMs and use this network to stabilize the motion on the NHIM for a periodically driven model system with two degrees of freedom. The method allows us to analyze the dynamics on the NHIM via Poincar\'e surfaces of section (PSOS) and to determine the transition state (TS) trajectory as a periodic orbit with the same periodicity as the driving saddle, viz. a fixed point of the PSOS surrounded by near-integrable tori. Based on Transition State Theory and a Floquet analysis of a periodic TS trajectory we compute the rate constant of the reaction with significantly reduced numerical effort compared to the propagation of a large trajectory ensemble.


title: Neural Sign Language Translation by Learning Tokenization

abstract: Sign Language Translation has attained considerable success recently, raising hopes for improved communication with the Deaf. A pre-processing step called tokenization improves the success of translations. Tokens can be learned from sign videos if supervised data is available. However, data annotation at the gloss level is costly, and annotated data is scarce. The paper utilizes Adversarial, Multitask, Transfer Learning to search for semi-supervised tokenization approaches without burden of additional labeling. It provides extensive experiments to compare all the methods in different settings to conduct a deeper analysis. In the case of no additional target annotation besides sentences, the proposed methodology attains 13.25 BLUE-4 and 36.28 ROUGE scores which improves the current state-of-the-art by 4 points in BLUE-4 and 5 points in ROUGE.


title: Multilevel Ensemble Kalman Filtering based on a sample average of independent EnKF estimators

abstract: We introduce a new multilevel ensemble Kalman filter method (MLEnKF) which consists of a hierarchy of independent samples of ensemble Kalman filters (EnKF). This new MLEnKF method is fundamentally different from the preexisting method introduced by Hoel, Law and Tempone in 2016, and it is suitable for extensions towards multi-index Monte Carlo based filtering methods. Robust theoretical analysis and supporting numerical examples show that under appropriate regularity assumptions, the MLEnKF method has better complexity than plain vanilla EnKF in the large-ensemble and fine-resolution limits, for weak approximations of quantities of interest. The method is developed for discrete-time filtering problems with finite-dimensional state space and linear observations polluted by additive Gaussian noise.


title: Assessment of Amazon Comprehend Medical: Medication Information Extraction

abstract: In November 27, 2018, Amazon Web Services (AWS) released Amazon Comprehend Medical (ACM), a deep learning based system that automatically extracts clinical concepts (which include anatomy, medical conditions, protected health information (PH)I, test names, treatment names, and medical procedures, and medications) from clinical text notes. Uptake and trust in any new data product relies on independent validation across benchmark datasets and tools to establish and confirm expected quality of results. This work focuses on the medication extraction task, and particularly, ACM was evaluated using the official test sets from the 2009 i2b2 Medication Extraction Challenge and 2018 n2c2 Track 2: Adverse Drug Events and Medication Extraction in EHRs. Overall, ACM achieved F-scores of 0.768 and 0.828. These scores ranked the lowest when compared to the three best systems in the respective challenges. To further establish the generalizability of its medication extraction performance, a set of random internal clinical text notes from NYU Langone Medical Center were also included in this work. And in this corpus, ACM garnered an F-score of 0.753.


title: A Relativistic GRW Flash Process With Interaction

abstract: In 2004, I described a relativistic version of the Ghirardi-Rimini-Weber (GRW) model of spontaneous wave function collapse for N non-interacting distinguishable particles. Here I present a generalized version for N interacting distinguishable particles. Presently, I do not know how to set up a similar model for indistinguishable particles or a variable number of particles. The present interacting model is constructed from a given interacting unitary Tomonaga-Schwinger type evolution between spacelike hypersurfaces, into which discrete collapses are inserted. I assume that this unitary evolution is interaction-local (i.e., no interaction at spacelike separation). The model is formulated in terms of Bell's flash ontology but is also compatible with Ghirardi's matter density ontology. It is non-local and satisfies microscopic parameter independence and no-signaling; it also works in curved space-time; in the non-relativistic limit, it reduces to the known non-relativistic GRW model.


title: Caterpillars on three and four leaves are sufficient to reconstruct normal networks

abstract: While every rooted binary phylogenetic tree is determined by its set of displayed rooted triples, such a result does not hold for an arbitrary rooted binary phylogenetic network. In particular, there exist two non-isomorphic rooted binary temporal normal networks that display the same set of rooted triples. Moreover, without any structural constraint on the rooted phylogenetic networks under consideration, similarly negative results have also been established for binets and trinets which are rooted subnetworks on two and three leaves, respectively. Hence, in general, piecing together a rooted phylogenetic network from such a set of small building blocks appears insurmountable. In contrast to these results, in this paper, we show that a rooted binary normal network is determined by its sets of displayed caterpillars (particular type of subtrees) on three and four leaves. The proof is constructive and realises a polynomial-time algorithm that takes the sets of caterpillars on three and four leaves displayed by a rooted binary normal network and, up to isomorphism, reconstructs this network.


title: DeepLocNet: Deep Observation Classification and Ranging Bias Regression for Radio Positioning Systems

abstract: WiFi technology has been used pervasively in fine-grained indoor localization, gesture recognition, and adaptive communication. Achieving better performance in these tasks generally boils down to differentiating Line-Of-Sight (LOS) from Non-Line-Of-Sight (NLOS) signal propagation reliably which generally requires expensive/specialized hardware due to the complex nature of indoor environments. Hence, the development of low-cost accurate positioning systems that exploit available infrastructure is not entirely solved. In this paper, we develop a framework for indoor localization and tracking of ubiquitous mobile devices such as smartphones using on-board sensors. We present a novel deep LOS/NLOS classifier which uses the Received Signal Strength Indicator (RSSI), and can classify the input signal with an accuracy of 85\%. The proposed algorithm can globally localize and track a smartphone (or robot) with a priori unknown location, and with a semi-accurate prior map (error within 0.8 m) of the WiFi Access Points (AP). Through simultaneously solving for the trajectory and the map of access points, we recover a trajectory of the device and corrected locations for the access points. Experimental evaluations of the framework show that localization accuracy is increased by using the trained deep network; furthermore, the system becomes robust to any error in the map of APs.


title: Concentration dependence of diffusion-limited reaction rates and its consequences

abstract: Diffusion-limited association reactions are ubiquitous in nature. They are particularly important for biological reactions, where the reaction rates are often determined by the diffusive transport of the molecules on two-dimensional surfaces, such as the cell membrane. The peculiarities of diffusion on two-dimensional surfaces may lead to nontrivial reaction kinetics, such as concentration dependent rate of association between two molecules. However, traditionally, the kinetics of biomolecular association reactions has been modeled using the law of mass action, which assumes that the rate of reaction is a concentration independent constant. In this paper, using multiscale molecular simulation, we investigate the concentration dependence of diffusion-limited association reactions on 2D surfaces. In particular, we quantify the influence of short-ranged pair interactions on the concentration dependence of the reaction rates and codify it in an empirical law. Using this law in a chemical kinetic model, we find that the the steady state behaviors of simple chemical systems are drastically modified by the presence of concentration dependent rates. In particular, we find that it leads to suppression of intrinsic noise in dimerization reaction and destabilizes robust oscillation in Lotka-Volterra predator-prey systems. In fact, we see a transition from robust to fine-tuned behavior in the latter. In addition, we show that concentration dependent reaction rates arise naturally in stochastic predator-prey systems due to intrinsic noise. We comment on the consequences of these results and discuss their implications in the modeling of complex chemical and biological systems. In particular, we comment on the range of validity of the law of mass action, which is a staple in all theoretical modeling of these systems.

categories: cond-mat.soft cond-mat.stat-mech physics.bio-ph physics.chem-ph q-bio.SC


title: Hyperbolic Secant Varieties of M-Curves

abstract: We relate the geometry of curves to the notion of hyperbolicity in real algebraic geometry. A hyperbolic variety is a real algebraic variety that (in particular) admits a real fibered morphism to a projective space whose dimension is equal to the dimension of the variety. We study hyperbolic varieties with a special interest in the case of hypersurfaces that admit a real algebraic ruling. The central part of the paper is concerned with secant varieties of real algebraic curves where the real locus has the maximal number of connected components, which is determined by the genus of the curve. For elliptic normal curves, we further obtain definite symmetric determinantal representations for the hyperbolic secant hypersurfaces, which implies the existence of symmetric Ulrich sheaves of rank one on these hypersurfaces.


title: Virtual Critical Coupling

abstract: Electromagnetic resonators are a versatile platform to harvest, filter and trap electromagnetic energy, at the basis of many applications from microwaves to optics. Resonators with a large intrinsic quality factor (Q) are highly desirable since they can store a large amount of energy, leading to sharp filtering and low loss. Exciting high-Q cavities with monochromatic signals, however, suffers from poor excitation efficiency, i.e., most of the impinging energy is lost in the form of reflection, since high-Q resonators are weakly coupled to external radiation. Although critical coupling eliminates reflections in steady-state by matching the intrinsic and coupling decay rates, this approach requires the introduction of loss in the resonator, causing dissipation and lowering the overall Q-factor. Here, we extend the notion of critical coupling to high-Q lossless resonators based on tailoring the temporal profile of the excitation wave. Utilizing coupled-mode theory, we demonstrate an effect analogous to critical coupling by mimicking loss with non-monochromatic excitations at complex frequencies. Remarkably, we show that this approach enables unitary excitation efficiency in open systems, even in the limit of extreme quality factors in the regime of quasi-bound states in the continuum.


title: Price of Anarchy in Multiuser Massive MIMO: Coordinated versus Uncoordinated Channel Tracking for High-Rate Internet of Things

abstract: Incorporation of high-rate internet of things (IoT) service into a massive MIMO framework is investigated. It is revealed that massive MIMO possess the inherent potential to offer such service provided it knows the channels for all devices. Our proposed method is to jointly estimate and track the channels of all devices irrespective of their current activity. Using the dynamical model for devices' channels evolution over time, optimal and sub-optimal trackers are developed for coordinated scenario. Furthermore, we introduce a new paradigm where the BS need not know the pilot access patterns of devices in advance which we refer to as uncoordinated setup. After motivating this scenario, we derive the optimal tracker which is intractable. Then, target tracking approaches are applied to address uncertainties in the measurements and derive sub-optimal trackers. Our proposed approaches explicitly address the channel aging problem and will not require downlink paging and uplink access request control channels which can become bottlenecks in crowded scenarios. The fundamental minimum mean square error (MMSE) gap between optimal coordinated and uncoordinated trackers which is defined as price of anarchy is evaluated and upper-bounded. Stability of optimal trackers is also investigated. Finally, performance of various proposed trackers are numerically compared.


title: Motion of finite-size spheres released in a turbulent boundary layer

abstract: Individual magnetic wax spheres with specific gravities of 1.006, 1.054 and 1.152 were released from rest on a smooth wall in water at friction Reynolds numbers, Re_{\tau}=680 and 1320 (d^+ = 58 and 122 viscous units, respectively). Three-dimensional tracking was conducted to understand the effect of turbulence and wall friction on sphere motions. Spheres subjected to sufficient mean shear initially lifted off of the wall before descending back towards it. These lifting spheres translated with the fluid above the wall, undergoing saltation or resuspension, with minimal rotation about any axis. By contrast, spheres that did not lift off upon release mainly slid along the wall. These denser spheres lagged the fluid more significantly due to greater wall friction. As they slid downstream, they began to roll forward after which small repeated lift-off events occurred. These spheres also rotated about both the streamwise and wall-normal axes. In all cases, the sphere trajectories were limited to the buffer and logarithmic regions, and all wall collisions were completely inelastic. Sphere streamwise velocities fluctuated up to 20% from the mean value even after the sphere had attained an approximate terminal velocity. In the plane parallel to the wall, the spheres migrated in the spanwise direction about 12% of the streamwise distance traveled suggesting that spanwise forces are important. The variations in sphere kinematics were likely induced by high and low momentum zones in the boundary layer, vortex shedding in the sphere wakes, and wall friction. The repeated lift-offs of the forward rolling denser sphere were attributed to a Magnus lift.


title: Network Inference using Sinusoidal Probing

abstract: The aim of this manuscript is to present a non-invasive method to recover the network structure of a dynamical system. We propose to use a controlled probing input and to measure the response of the network, in the spirit of what is done to determine oscillation modes in large electrical networks. For a large class of dynamical systems, we show that this approach is analytically tractable and we confirm our findings by numerical simulations of networks of Kuramoto oscillators. Our approach also allows us to determine the number of agents in the network by probing and measuring a single one of them.


title: Stochastic Navier-Stokes Equations and Related Models

abstract: Regularization by noise for certain classes of fluid dynamic equations, a theme dear to Giuseppe Da Prato (see G. Da Prato and A. Debussche, Ergodicity for the 3D stochastic Navier-Stokes equations, J. Math. Pures Appl., 2003), is reviewed focusing on 3D Navier-Stokes equations and dyadic models of turbulence.


title: Overfitting Can Be Harmless for Basis Pursuit, But Only to a Degree

abstract: Recently, there have been significant interests in studying the so-called "double-descent" of the generalization error of linear regression models under the overparameterized and overfitting regime, with the hope that such analysis may provide the first step towards understanding why overparameterized deep neural networks (DNN) still generalize well. However, to date most of these studies focused on the min $\ell_2$-norm solution that overfits the data. In contrast, in this paper we study the overfitting solution that minimizes the $\ell_1$-norm, which is known as Basis Pursuit (BP) in the compressed sensing literature. Under a sparse true linear regression model with $p$ i.i.d. Gaussian features, we show that for a large range of $p$ up to a limit that grows exponentially with the number of samples $n$, with high probability the model error of BP is upper bounded by a value that decreases with $p$. To the best of our knowledge, this is the first analytical result in the literature establishing the double-descent of overfitting BP for finite $n$ and $p$. Further, our results reveal significant differences between the double-descent of BP and min $\ell_2$-norm solutions. Specifically, the double-descent upper-bound of BP is independent of the signal strength, and for high SNR and sparse models the descent-floor of BP can be much lower and wider than that of min $\ell_2$-norm solutions.


title: Automorphic Schwarzian equations

abstract: This paper concerns the study of the Schwarz differential equation $\{h,\tau \}=s\,E_4(\tau)$ where $E_4$ is the weight 4 Eisenstein series and $s$ is a complex parameter. In particular, we determine all the values of $s$ for which the solutions $h$ are modular functions for a finite index subgroup of $\mbox{SL}_2(\mathbb Z)$. We do so using the theory of equivariant functions on the complex upper-half plane as well as an analysis of the representation theory of $\mbox{SL}_2(\mathbb Z)$. This also leads to the solutions to the Fuchsian differential equation $y''+s\,E_4\,y=0$.


title: Free actions of large groups on complex threefolds

abstract: We classify compact K\"ahler threefolds $X$ with a free group of automorphisms acting freely on $X$.


title: Active Learning for Identification of Linear Dynamical Systems

abstract: We propose an algorithm to actively estimate the parameters of a linear dynamical system. Given complete control over the system's input, our algorithm adaptively chooses the inputs to accelerate estimation. We show a finite time bound quantifying the estimation rate our algorithm attains and prove matching upper and lower bounds which guarantee its asymptotic optimality, up to constants. In addition, we show that this optimal rate is unattainable when using Gaussian noise to excite the system, even with optimally tuned covariance, and analyze several examples where our algorithm provably improves over rates obtained by playing noise. Our analysis critically relies on a novel result quantifying the error in estimating the parameters of a dynamical system when arbitrary periodic inputs are being played. We conclude with numerical examples that illustrate the effectiveness of our algorithm in practice.


title: Excluding a ladder

abstract: A ladder is a $2 \times k$ grid graph. When does a graph class $\mathcal{C}$ exclude some ladder as a minor? We show that this is the case if and only if all graphs $G$ in $\mathcal{C}$ admit a proper vertex coloring with a bounded number of colors such that for every $2$-connected subgraph $H$ of $G$, there is a color that appears exactly once in $H$. This type of vertex coloring is a relaxation of the notion of centered coloring, where for every connected subgraph $H$ of $G$, there must be a color that appears exactly once in $H$. The minimum number of colors in a centered coloring of $G$ is the treedepth of $G$, and it is known that classes of graphs with bounded treedepth are exactly those that exclude a fixed path as a subgraph, or equivalently, as a minor. In this sense, the structure of graphs excluding a fixed ladder as a minor resembles the structure of graphs without long paths. Another similarity is as follows: It is an easy observation that every connected graph with two vertex-disjoint paths of length $k$ has a path of length $k+1$. We show that every $3$-connected graph which contains as a minor a union of sufficiently many vertex-disjoint copies of a $2 \times k$ grid has a $2 \times (k+1)$ grid minor. Our structural results have applications to poset dimension. We show that posets whose cover graphs exclude a fixed ladder as a minor have bounded dimension. This is a new step towards the goal of understanding which graphs are unavoidable as minors in cover graphs of posets with large dimension.


title: Accelerating Cooperative Planning for Automated Vehicles with Learned Heuristics and Monte Carlo Tree Search

abstract: Efficient driving in urban traffic scenarios requires foresight. The observation of other traffic participants and the inference of their possible next actions depending on the own action is considered cooperative prediction and planning. Humans are well equipped with the capability to predict the actions of multiple interacting traffic participants and plan accordingly, without the need to directly communicate with others. Prior work has shown that it is possible to achieve effective cooperative planning without the need for explicit communication. However, the search space for cooperative plans is so large that most of the computational budget is spent on exploring the search space in unpromising regions that are far away from the solution. To accelerate the planning process, we combined learned heuristics with a cooperative planning method to guide the search towards regions with promising actions, yielding better solutions at lower computational costs.


title: DYNOTEARS: Structure Learning from Time-Series Data

abstract: We revisit the structure learning problem for dynamic Bayesian networks and propose a method that simultaneously estimates contemporaneous (intra-slice) and time-lagged (inter-slice) relationships between variables in a time-series. Our approach is score-based, and revolves around minimizing a penalized loss subject to an acyclicity constraint. To solve this problem, we leverage a recent algebraic result characterizing the acyclicity constraint as a smooth equality constraint. The resulting algorithm, which we call DYNOTEARS, outperforms other methods on simulated data, especially in high-dimensions as the number of variables increases. We also apply this algorithm on real datasets from two different domains, finance and molecular biology, and analyze the resulting output. Compared to state-of-the-art methods for learning dynamic Bayesian networks, our method is both scalable and accurate on real data. The simple formulation and competitive performance of our method make it suitable for a variety of problems where one seeks to learn connections between variables across time.


title: Detecting Anomalous Time Series by GAMLSS-Akaike-Weights-Scoring

abstract: An extensible statistical framework for detecting anomalous time series including those with heavy-tailed distributions and non-stationarity in higher-order moments is introduced based on penalized likelihood distributional regression. Specifically, generalized additive models for location, scale, and shape are used to infer sample path representations defined by a parametric distribution with parameters comprised of basis functions. Akaike weights are then applied to each model and time series, yielding a probability measure that can be effectively used to classify and rank anomalous time series. A mathematical exposition is also given to justify the proposed Akaike weight scoring under a suitable model embedding as a way to asymptotically identify anomalous time series. Studies evaluating the methodology on both multiple simulations and real-world datasets also confirm that high accuracy can be obtained detecting many different and complex types of shape anomalies. Both code implementing GAWS for running on a local machine and the datasets referenced in this paper are available online.


title: Exceptional scatteredness in prime degree

abstract: Let $q$ be an odd prime power and $n$ be a positive integer. Let $\ell\in \mathbb F_{q^n}[x]$ be a $q$-linearised $t$-scattered polynomial of linearized degree $r$. Let $d=\max\{t,r\}$ be an odd prime number. In this paper we show that under these assumptions it follows that $\ell=x$. Our technique involves a Galois theoretical characterization of $t$-scattered polynomials combined with the classification of transitive subgroups of the general linear group over the finite field $\mathbb F_q$.


title: How to measure musical preference on Facebook? Evidence from a mixed-method data collection

abstract: More and more digital data is available for social science analysis. This provides new ways of measuring several concepts. But when we start using new data sources, we have to understand how the new data source could be processed and how it could be analysed effectively. It is especially for Facebook data since there is no established gold standard analysis-framework. However, researchers have in-depth knowledge on how to measure different concepts using survey data. Thus, cross-referencing Facebook data with survey data is a reasonable way to support Facebook data analysis at different decision points. In this paper, we present how music preference could be measured by Facebook data and how survey data could support the selection of main indicators. Based on our results, we provide some general suggestions for Facebook data processing and indicator operationalization.


title: On the Erd\H{o}s distance problem

abstract: In this paper, using the method of compression, we recover the lower bound for the Erd\H{o}s unit distance problem and provide an alternative proof to the distinct distance conjecture. In particular, in $\mathbb{R}^k$ for all $k\geq 2$, we have \begin{align} \#\bigg\{||\vec{x_j}-\vec{x_t}||:~||\vec{x_j}-\vec{x_t}||=1,~1\leq t,j \leq n,~\vec{x_j},~\vec{x}_t \in \mathbb{R}^k\bigg\}\gg_k \frac{\sqrt{k}}{2}n^{1+o(1)}.\nonumber \end{align}We also show that \begin{align} \# \bigg\{d_j:d_j=||\vec{x_s}-\vec{y_t}||,~d_j\neq d_i,~1\leq s,t\leq n\bigg\}\gg_k \frac{\sqrt{k}}{2}n^{\frac{2}{k}-o(1)}.\nonumber \end{align}These lower bounds generalizes the lower bounds of the Erd\H{o}s unit distance and the distinct distance problem to higher dimensions.