title: Learning Compressed Sentence Representations for On-Device Text Processing

abstract: Vector representations of sentences, trained on massive text corpora, are widely used as generic sentence embeddings across a variety of NLP problems. The learned representations are generally assumed to be continuous and real-valued, giving rise to a large memory footprint and slow retrieval speed, which hinders their applicability to low-resource (memory and computation) platforms, such as mobile devices. In this paper, we propose four different strategies to transform continuous and generic sentence embeddings into a binarized form, while preserving their rich semantic information. The introduced methods are evaluated across a wide range of downstream tasks, where the binarized sentence embeddings are demonstrated to degrade performance by only about 2% relative to their continuous counterparts, while reducing the storage requirement by over 98%. Moreover, with the learned binary representations, the semantic relatedness of two sentences can be evaluated by simply calculating their Hamming distance, which is more computational efficient compared with the inner product operation between continuous embeddings. Detailed analysis and case study further validate the effectiveness of proposed methods.


title: Polarization-dependent conductivity of grain boundaries in BiFeO3 thin films

abstract: Charge transport across the interfaces in complex oxides attracts a lot of attention because it allows creating novel functionalities useful for device applications. In particular, it has been observed that movable domain walls in epitaxial BiFeO3 films possess enhanced conductivity that can be used for read out in ferroelectric-based memories. In this work, the relation between the polarization and conductivity in sol-gel BiFeO3 films with special emphasis on grain boundaries as natural interfaces in polycrystalline ferroelectrics is investigated. The grains exhibit self-organized domain structure in these films, so that the "domain clusters" consisting of several grains with aligned polarization directions are formed. Surprisingly, grain boundaries between these clusters (with antiparallel polarization direction) have significantly higher electrical conductivity in comparison to "inter-cluster" grain boundaries, in which the conductivity was even smaller than in the bulk. As such, polarization-dependent conductivity of the grain boundaries was observed for the first time in ferroelectric thin films. The results are rationalized by thermodynamic modelling combined with finite element simulations of the charge and stress accumulation at the grain boundaries giving major contribution to conductivity. The observed polarization-dependent conductivity of grain boundaries in ferroelectrics opens up a new avenue for exploiting these materials in electronic devices.


title: Gamma-ray emission in radio galaxies under the VLBI scope -- I. Parsec-scale kinematics and high-energy properties of $\gamma$-ray detected TANAMI radio galaxies

abstract: In the framework of the TANAMI multi-wavelength and VLBI monitoring, we study the evolution of the parsec-scale radio emission in radio galaxies in the southern hemisphere and their relationship to the $\gamma$-ray properties. In this first paper, we focus on Fermi-LAT-detected sources. We perform a kinematic analysis for five $\gamma$-ray detected radio galaxies using multi-epoch 8.4 GHz VLBI images, deriving limits on intrinsic jet parameters. We analyzed Fermi-LAT data in order to study possible connections between the $\gamma$-ray properties and the pc-scale jets of Fermi-LAT-detected radio galaxies, both in terms of variability and average properties. We discuss the individual source results and draw preliminary conclusions on sample properties including published VLBI results from the MOJAVE survey, with a total of fifteen sources. We find that the first $\gamma$-ray detection of Pictor A might be associated with the passage of a new VLBI component through the radio core. For the peculiar AGN PKS 0521-36, we detect subluminal parsec-scale jet motions, and we confirm the presence of fast $\gamma$-ray variability in the source down to timescales of 6 hours. We robustly confirm the presence of significant superluminal motion, up to $\beta_{app}\sim$3, in the jet of the TeV radio galaxy PKS 0625-35. Finally, we place a lower limit on the age of the Compact Symmetric Object (CSO) PKS 1718-649. We draw some preliminary conclusions on the relationship between pc-scale jets and $\gamma$-ray emission in radio galaxies. We find that the VLBI core flux density correlates with the $\gamma$-ray flux, as seen in blazars. On the other hand, the $\gamma$-ray luminosity does not show any dependence on the core brightness temperature and core dominance, two indicators of Doppler boosting, suggesting that $\gamma$-ray emission in radio galaxies is not driven by orientation-dependent effects.


title: Optimal designs for estimating individual coefficients in polynomial regression with no intercept

abstract: In a seminal paper \cite{studden1968} characterized $c$-optimal designs in regression models, where the regression functions form a Chebyshev system. He used these results to determine the optimal design for estimating the individual coefficients in a polynomial regression model on the interval $[-1,1]$ explicitly. In this note we identify the optimal design for estimating the individual coefficients in a polynomial regression model with no intercept (here the regression functions do not form a Chebyshev system).


title: Multi-resolution Multi-task Gaussian Processes

abstract: We consider evidence integration from potentially dependent observation processes under varying spatio-temporal sampling resolutions and noise levels. We develop a multi-resolution multi-task (MRGP) framework while allowing for both inter-task and intra-task multi-resolution and multi-fidelity. We develop shallow Gaussian Process (GP) mixtures that approximate the difficult to estimate joint likelihood with a composite one and deep GP constructions that naturally handle biases in the mean. By doing so, we generalize and outperform state of the art GP compositions and offer information-theoretic corrections and efficient variational approximations. We demonstrate the competitiveness of MRGPs on synthetic settings and on the challenging problem of hyper-local estimation of air pollution levels across London from multiple sensing modalities operating at disparate spatio-temporal resolutions.


title: Topological invariants for interface modes

abstract: We consider topologically non-trivial interface Hamiltonians, which find several applications in materials science and geophysical fluid flows. The non-trivial topology manifests itself in the existence of topologically protected, asymmetric edge states at the interface between two two-dimensional half spaces in different topological phases. It is characterized by a quantized interface conductivity. The objective of this paper is to compute such a conductivity and show its stability under perturbations. We present two methods. The first one computes the conductivity using the winding number of branches of absolutely continuous spectrum of the interface Hamiltonian. This calculation is independent of any bulk properties but requires a sufficient understanding of the spectral decomposition of the Hamiltonian. In the fluid flow setting, it also applies in cases where the so-called bulk-interface correspondence fails. The second method establishes a bulk-interface correspondence between the interface conductivity and a so-called bulk-difference invariant. We introduce the bulk-difference invariants characterizing pairs of half spaces. We then relate the interface conductivity to the bulk-difference invariant by means of a Fedosov-H\"ormander formula, which computes the index of a related Fredholm operator. The two methods are used to compute invariants for representative $2\times2$ and $3\times3$ systems of equations that appear in the applications.


title: On the Geramita-Harbourne-Migliore conjecture

abstract: Let $\Sigma$ be a finite collection of linear forms in $\mathbb K[x_0,\ldots,x_n]$, where $\mathbb K$ is a field. Denote ${\rm Supp}(\Sigma)$ to be the set of all nonproportional elements of $\Sigma$, and suppose ${\rm Supp}(\Sigma)$ is generic, meaning that any $n+1$ of its elements are linearly independent. Let $1\leq a\leq |\Sigma|$. In this article we prove the conjecture that the ideal generated by (all) $a$-fold products of linear forms of $\Sigma$ has linear graded free resolution. As a consequence we prove the Geramita-Harbourne-Migliore conjecture concerning the primary decomposition of ordinary powers of defining ideals of star configurations, and we also determine the resurgence of these ideals.


title: Logarithmic integrals, zeta values, and tiered binomial coefficients

abstract: We study logarithmic integrals of the form $\int_0^1 x^i\ln^n(x)\ln^m(1-x)dx$. They are expressed as a rational linear combination of certain rational numbers $(n,m)_i$, which we call tiered binomial coefficients, and products of the zeta values $\zeta(2)$, $\zeta(3)$,\dots. Various properties of the tiered binomial coefficients are established. They involve, amongst others, the binomial transform, truncated multiple zeta and multiple zeta star values, as well as special functions. As an application we discuss the limit law of the number of comparisons of the Quicksort algorithm: we reprove that the moments of the limit law are rational polynomials in the zeta values. A novel expression for the cumulants of the Quicksort limit is also presented.


title: Some algebras that are not silting connected

abstract: We give examples of finite-dimensional algebras $A$ for which the silting objects in $K^b(\mbox{proj-}A)$ are not connected by any sequence of (possibly reducible) silting mutations. The argument is based on the fact that silting mutation preserves invariance under twisting by a fixed algebra automorphism, combined with the existence of spherical modules that are not invariant under such a twist.


title: Foreword to the Focus Issue on Machine Learning in Astronomy and Astrophysics

abstract: Astronomical observations already produce vast amounts of data through a new generation of telescopes that cannot be analyzed manually. Next-generation telescopes such as the Large Synoptic Survey Telescope and the Square Kilometer Array are planned to become operational in this decade and the next, and will increase the data volume by many orders of magnitude. The increased spatial, temporal and spectral resolution afford a powerful magnifying lens on the physical processes that underlie the data but, at the same time, generate unprecedented complexity hard to exploit for knowledge extraction. It is therefore imperative to develop machine intelligence, machine learning (ML) in particular, suitable for processing the amount and variety of astronomical data that will be collected, and capable of answering scientific questions based on the data. Astronomical data exhibit the usual challenges associated with 'big data' such as immense volumes, high dimensionality, missing or highly distorted observations. In addition, astronomical data can exhibit large continuous observational gaps, very low signal-to-noise ratio and the need to distinguish between true missing data and non-detections due to upper limits). There are strict laws of physics behind the data production which can be assimilated into ML mechanisms to improve over general off-the-shelf state-of-the-art methods. Significant progress in the face of these challenges can be achieved only via the new discipline of Astroinformatics: a symbiosis of diverse disciplines, such as ML, probabilistic modeling, astronomy and astrophysics, statistics, distributed computing and natural computation. This editorial summarizes the contents of a soon to appear Focus Issue of the PASP on Machine Learning in Astronomy and Astrophysics (with contributions by 69 authors representing 15 countries, from 6 continents).


title: Precision measurement of the $\Lambda_c^+$, $\Xi_c^+$ and $\Xi_c^0$ baryon lifetimes

abstract: We report measurements of the lifetimes of the $\Lambda_c^+$, $\Xi_c^+$ and $\Xi_c^0$ charm baryons using proton-proton collision data at center-of-mass energies of 7 and 8\tev, corresponding to an integrated luminosity of 3.0 fb$^{-1}$, collected by the LHCb experiment. The charm baryons are reconstructed through the decays $\Lambda_c^+\to pK^-\pi^+$, $\Xi_c^+\to pK^-\pi^+$ and $\Xi_c^0\to pK^-K^-\pi^+$, and originate from semimuonic decays of beauty baryons. The lifetimes are measured relative to that of the $D^+$ meson, and are determined to be \begin{align*} \tau_{\Lambda_c^+} &= 203.5\pm1.0\pm1.3\pm1.4~{\rm fs}, \newline \tau_{\Xi_c^+} &= 456.8\pm3.5\pm2.9\pm3.1~{\rm fs}, \newline \tau_{\Xi_c^0} &= 154.5\pm1.7\pm1.6\pm1.0~{\rm fs}, \end{align*} where the uncertainties are statistical, systematic, and due to the uncertainty in the $D^+$ lifetime. The measurements are approximately 3--4 times more precise than the current world average values. The $\Lambda_c^+$ and $\Xi_c^+$ lifetimes are in agreement with previous measurements; however, the $\Xi_c^0$ baryon lifetime is approximately 3.3 standard deviations larger than the world average value.


title: Towards Lakosian Multilingual Software Design Principles

abstract: Large software systems often comprise programs written in different programming languages. In the case when cross-language interoperability is accomplished with a Foreign Function Interface (FFI), for example pybind11, Boost.Python, Emscripten, PyV8, or JNI, among many others, common software engineering tools, such as call-graph analysis, are obstructed by the opacity of the FFI. This complicates debugging and fosters potential inefficiency and security problems. One contributing issue is that there is little rigorous software design advice for multilingual software. In this paper, we present our progress towards a more rigorous design approach to multilingual software. The approach is based on the existing approach to the design of large-scale C++ systems developed by Lakos. The Lakosian approach is one of the few design methodologies to address physical design rather than just logical design. Using the MLSA toolkit developed in prior work for analysis of multilingual software, we focus in on one FFI -- the pybind11 FFI. An extension to the Lakosian C++ design rules is proposed to address multilingual software that uses pybind11. Using a sample of 50 public GitHub repositories that use pybind11, we measure how many repositories would currently satisfy these rules. We conclude with a proposed generalization of the pybind11-based rules for any multilingual software using an FFI interface.


title: Surface passivated and encapsulated ZnO atomic layer by high-$\kappa$ ultrathin MgO layer

abstract: Atomically transparent vertically aligned ZnO-based van der Waals material have been developed by surface passivation and encapsulation with atomic layers of MgO using materials by design; the physical properties investigated. The passivation and encapsulation led to a remarkable improvement in optical and electronic properties. The valence-band offset $\Delta E_v$ between MgO and ZnO, ZnO and MgO/ZnO, and ZnO and MgO/ZnO/MgO heterointerfaces are determined to be 0.37 $\pm$0.02, -0.05$\pm$0.02, and -0.11$\pm$0.02 eV, respectively; the conduction-band offset $\Delta E_c$ is deduced to be 0.97$\pm$0.02, 0.46$\pm$0.02, and 0.59$\pm$0.02 eV indicating straddling type-I in MgO and ZnO, and staggering type-II heterojunction band alignment in ZnO and the various heterostructures. The band-offsets and interfacial charge transfer are used to explain the origin of $n$-type conductivity in the superlattices. Enhanced optical absorption due to carrier confinement in the layers demonstrates that MgO is an excellent high-$\kappa$ dielectric gate oxide for encapsulating ZnO-based optoelectronic devices.


title: Magnetic dipole moments of the spin-$\frac{3}{2}$ doubly heavy baryons

abstract: The magnetic dipole moments of the spin-$\frac{3}{2}$ doubly charmed, bottom and charmed-bottom baryons are obtained by means of the light-cone QCD sum rule. The magnetic dipole moments of these baryons encode essential knowledge of their inner structure and shape deformations. The numerical results are given as, $\mu_{\Xi_{cc}^{*++}} = 2.94 \pm 0.95$, $\mu_{\Xi_{cc}^{*+}} = - 0.67 \pm 0.11$, $\mu_{\Omega_{cc}^{*+}} =- 0.52 \pm 0.07$, $\mu_{\Xi_{bb}^{*0}} = 2.30 \pm 0.55$, $\mu_{\Xi_{bb}^{*-}} = -1.39 \pm 0.32$, $\mu_{\Omega_{bb}^{*-}} = -1.56 \pm 0.33$, $\mu_{\Xi_{bc}^{*+}} = 2.63 \pm 0.82$, $\mu_{\Xi_{bc}^{*0}} = - 0.96 \pm 0.32$ and $\mu_{\Omega_{bc}^{*+}} =- 1.11 \pm 0.33$, respectively.


title: Investigation of Near-Surface Defects of Nanodiamonds by High-Frequency EPR and DFT Calculation

abstract: Nanodiamond (ND) hosting nitrogen-vacancy (NV) centers is a promising platform for quantum sensing applications. Sensitivity of the applications using NV centers in NDs is often limited due to presence of paramagnetic impurity contents near the ND surface. Here, we investigate near-surface paramagnetic impurities in NDs. Using high-frequency (HF) electron paramagnetic resonance spectroscopy, the near-surface paramagnetic impurity within the shell of NDs is probed and its g-value is determined to be 2.0028(3). Furthermore, HF electron-electron double resonance-detected nuclear magnetic resonance spectroscopy and a first principle calculation show that a possible structure of the near-surface impurity is the negatively charged vacancy V-. The identification of the near-surface impurity by the present investigation provides a promising pathway to improve the NV properties in NDs and the NV-based sensing techniques.


title: Lightcone Construction for HACC Cosmological Simulations with LANTERN

abstract: This document describes the LANTERN (Lightcone generAtion via sNapshoT intERpolatioN) code module, our current approach for generating lightcones from HACC [5] simulation products. In sec. 1, we derive the condition which defines all spacetime events that can be seen by an observer at any given time (we parameterize the surface of the observer's past lightcone). In sec. 2, we describe how we apply this condition to find approximate lightcone-crossing times of objects whose positions and velocities are known only at discrete time snapshots. Higher order effects and corrections are discussed in sec. 3. Sec. 4 describes simulation box replication capabilities that allow generation of arbitrary lightcone sizes.


title: Updated measurement of time-dependent $CP$-violating observables in $B^{0}_{s}\to J/\psi K^+ K^-$ decays

abstract: The decay-time-dependent $CP$ asymmetry in $B^{0}_{s}\to J/\psi K^{+} K^{-}$ decays is measured using proton-proton collision data, corresponding to an integrated luminosity of $1.9\,\mathrm{fb^{-1}}$, collected with the LHCb detector at a centre-of-mass energy of $13\,\mathrm{TeV}$ in 2015 and 2016. Using a sample of approximately 117\,000 signal decays with an invariant $K^{+} K^{-}$ mass in the vicinity of the $\phi(1020)$ resonance, the $CP$-violating phase $\phi_s$ is measured, along with the difference in decay widths of the light and heavy mass eigenstates of the $B^{0}_{s}$-$\bar{B}^{0}_{s}$ system, $\Delta\Gamma_s$. The difference of the average $B^{0}_{s}$ and $B^{0}$ meson decay widths, $\Gamma_s-\Gamma_d$, is determined using in addition a sample of $B^{0} \to J/\psi K^{+} \pi^{-}$ decays. The values obtained are $\phi_s = -0.083\pm0.041\pm0.006\,\mathrm{rad}$, $\Delta\Gamma_s = 0.077 \pm 0.008 \pm 0.003\, \mathrm{ps^{-1}}$ and $\Gamma_s-\Gamma_d = -0.0041 \pm 0.0024 \pm 0.0015\,\mathrm{ps^{-1}}$, where the first uncertainty is statistical and the second systematic. These are the most precise single measurements of these quantities to date and are consistent with expectations based on the Standard Model and with a previous LHCb analysis of this decay using data recorded at centre-of-mass energies 7 and 8 TeV. Finally, the results are combined with recent results from $B^{0}_{s}\to J/\psi \pi^{+} \pi^{-}$ decays obtained using the same dataset as this analysis, and with previous independent LHCb results.


title: The Age-Period-Cohort-Interaction Model for Describing and Investigating Inter-Cohort Deviations and Intra-Cohort Life-Course Dynamics

abstract: Social scientists have frequently sought to understand the distinct effects of age, period, and cohort, but disaggregation of the three dimensions is difficult because cohort = period - age. We argue that this technical difficulty reflects a disconnection between how cohort effect is conceptualized and how it is modeled in the traditional age-period-cohort framework. We propose a new method, called the age-period-cohort-interaction (APC-I) model, that is qualitatively different from previous methods in that it represents Ryder's (1965) theoretical account about the conditions under which cohort differentiation may arise. This APC-I model does not require problematic statistical assumptions and the interpretation is straightforward. It quantifies inter-cohort deviations from the age and period main effects and also permits hypothesis testing about intra-cohort life-course dynamics. We demonstrate how this new model can be used to examine age, period, and cohort patterns in women's labor force participation.


title: Transfinite mean value interpolation over polygons

abstract: Mean value interpolation is a method for fitting a smooth function to piecewise-linear data prescribed on the boundary of a polygon of arbitrary shape, and has applications in computer graphics and curve and surface modelling. The method generalizes to transfinite interpolation, i.e., to any continuous data on the boundary but a mathematical proof that interpolation always holds has so far been missing. The purpose of this note is to complete this gap in the theory.


title: An axiomatic nonparametric production function estimator: Modeling production in Japan's cardboard industry

abstract: We develop a new approach to estimate a production function based on the economic axioms of the Regular Ultra Passum law and convex non-homothetic input isoquants. Central to the development of our estimator is stating the axioms as shape constraints and using shape constrained nonparametric regression methods. We implement this approach using data from the Japanese corrugated cardboard industry from 1997-2007. Using this new approach, we find most productive scale size is a function of the capital-to-labor ratio and the largest firms operate close to the largest most productive scale size associated with a high capital-to-labor ratio. We measure the productivity growth across the panel periods based on the residuals from our axiomatic model. We also decompose productivity into scale, input mix, and unexplained effects to clarify the sources the productivity differences and provide managers guidance to make firms more productive.


title: Frequentist Inference without Repeated Sampling

abstract: Frequentist inference typically is described in terms of hypothetical repeated sampling but there are advantages to an interpretation that uses a single random sample. Contemporary examples are given that indicate probabilities for random phenomena are interpreted as classical probabilities, and this interpretation is applied to statistical inference using urn models. Both classical and limiting relative frequency interpretations can be used to communicate statistical inference, and the effectiveness of each is discussed. Recent descriptions of p-values, confidence intervals, and power are viewed through the lens of classical probability based on a single random sample from the population.


title: Transformation of XML Documents with Prolog

abstract: Transforming XML documents with conventional XML languages, like XSL-T, is disadvantageous because there is too lax abstraction on the target language and it is rather difficult to recognize rule-oriented transformations. Prolog as a programming language of declarative paradigm is especially good for implementation of analysis of formal languages. Prolog seems also to be good for term manipulation, complex schema-transformation and text retrieval. In this report an appropriate model for XML documents is proposed, the basic transformation language for Prolog LTL is defined and the expressiveness power compared with XSL-T is demonstrated, the implementations used throughout are multi paradigmatic.


title: Trepan Reloaded: A Knowledge-driven Approach to Explaining Artificial Neural Networks

abstract: Explainability in Artificial Intelligence has been revived as a topic of active research by the need of conveying safety and trust to users in the `how' and `why' of automated decision-making. Whilst a plethora of approaches have been developed for post-hoc explainability, only a few focus on how to use domain knowledge, and how this influences the understandability of global explanations from the users' perspective. In this paper, we show how ontologies help the understandability of global post-hoc explanations, presented in the form of symbolic models. In particular, we build on Trepan, an algorithm that explains artificial neural networks by means of decision trees, and we extend it to include ontologies modeling domain knowledge in the process of generating explanations. We present the results of a user study that measures the understandability of decision trees using a syntactic complexity measure, and through time and accuracy of responses as well as reported user confidence and understandability. The user study considers domains where explanations are critical, namely, in finance and medicine. The results show that decision trees generated with our algorithm, taking into account domain knowledge, are more understandable than those generated by standard Trepan without the use of ontologies.


title: Topological Formulae for the Zeroth Cohomology of Line Bundles on del Pezzo and Hirzebruch Surfaces

abstract: We show that the zeroth cohomology of effective line bundles on del Pezzo and Hirzebruch surfaces can always be computed in terms of a topological index.


title: Analysis of a Century's Worth of AR Scorpii Photometry from DASCH and ASAS-SN

abstract: AR Scorpii (AR Sco) is a binary star system containing the only known white dwarf (WD) pulsar. Previously reported photometric datasets only provide coverage back to 2005, but we extend the observational baseline for AR Sco back to the beginning of the 20th century by analyzing observations from the Digital Access to a Sky Century at Harvard project (DASCH). We find that the orbital waveform of AR Sco remained constant across that baseline with no significant deviations from its average brightness. This result strongly suggests that the absence of accretion in modern observations is a long-term feature of AR Sco. Additionally, the DASCH light curve provides an opportunity to test an earlier prediction that an obliquity of the WD would result in a precessional period observable in long-term studies of the orbital light curve. The DASCH observations do not indicate the presence of such a period, and we show that previous, inconclusive tests of this hypothesis were insensitive to the existence of a precessional period. Furthermore, the long DASCH baseline enables us to constrain the rate of change of the orbital frequency to $\dot{\nu} \lesssim 3.8\times10^{-20}$ Hz s$^{-1}$, constraining the efficacy of magnetic braking as a mechanism of angular-momentum loss in this system. Finally, we discuss how the combination of the orbital waveform's stability, high amplitude, and short period should make it possible to identify additional WD pulsars in all-sky survey data.


title: Extraction of hierarchical functional connectivity components in human brain using resting-state fMRI

abstract: The study of hierarchy in networks of the human brain has been of significant interest among the researchers as numerous studies have pointed out towards a functional hierarchical organization of the human brain. This paper provides a novel method for the extraction of hierarchical connectivity components in the human brain using resting-state fMRI. The method builds upon prior work of Sparse Connectivity Patterns (SCPs) by introducing a hierarchy of sparse overlapping patterns. The components are estimated by deep factorization of correlation matrices generated from fMRI. The goal of the paper is to extract interpretable hierarchical patterns using correlation matrices where a low rank decomposition is formed by a linear combination of a high rank decomposition. We formulate the decomposition as a non-convex optimization problem and solve it using gradient descent algorithms with adaptive step size. We also provide a method for the warm start of the gradient descent using singular value decomposition. We demonstrate the effectiveness of the developed method on two different real-world datasets by showing that multi-scale hierarchical SCPs are reproducible between sub-samples and are more reproducible as compared to single scale patterns. We also compare our method with existing hierarchical community detection approaches. Our method also provides novel insight into the functional organization of the human brain.


title: Full non-differentiability sets of typical Lipschitz functions

abstract: In this paper we prove that the typical Lipschitz function has no directional derivative at any point of a Borel set $E$ if and only if $E$ is contained in a countable union of closed purely unrectifiable sets.


title: Stability of the Kaczmarz Reconstruction for Stationary Sequences

abstract: The Kaczmarz algorithm is an iterative method to reconstruct an unknown vector $f$ from inner products $\langle f , \varphi_{n} \rangle $. We consider the problem of how additive noise affects the reconstruction under the assumption that $\{ \varphi_{n} \}$ form a stationary sequence. Unlike other reconstruction methods, such as frame reconstructions, the Kaczmarz reconstruction is unstable in the presence of noise. We show, however, that the reconstruction can be stabilized by relaxing the Kaczmarz algorithm; this relaxation corresponds to Abel summation when viewed as a reconstruction on the unit disc. We show, moreover, that for certain noise profiles, such as those that lie in $H^{\infty}(\mathbb{D})$ or certain subspaces of $H^{2}(\mathbb{D})$, the relaxed version of the Kaczmarz algorithm can fully remove the corruption by noise in the inner products. Using the spectral representation of stationary sequences, we show that our relaxed version of the Kaczmarz algorithm also stabilizes the reconstruction of Fourier series expansions in $L^2(\mu)$ when $\mu$ is singular.


title: Anisotropic Tolman VII solution by gravitational decoupling

abstract: Using the gravitational decoupling by the minimal geometric deformation approach, we build an anisotropic version of the well-known Tolman VII solution, determining an exact and physically acceptable interior two-fluid solution that can represent behavior of compact objects. Comparison of the effective density and density of the perfect fluid is demonstrated explicitly. We show that the radial and tangential pressure are different in magnitude giving thus the anisotropy of the modified Tolman VII solution. The dependence of the anisotropy on the coupling constant is also shown.


title: Unification of Template-Expansion and XML-Validation

abstract: The processing of XML documents often includes creation and validation. These two operations are typically performed in two different nodes within a computer network that do not correlate with each other. The process of creation is also called instantiation of a template and can be described by filling a template with data from external repositories. Initial access to arbitrary sources can be formulated as an expression of certain command languages like XPath. Filling means copying invariant element nodes to the target document and unfolding variable parts from a given template. Validation is a descision problem returning true if a given XML document satisfies a schema and false otherwise. The main subject is to find a language that unions the template expansion and the validation. [..].


title: Driving Path Stability in VANETs

abstract: Vehicular Ad Hoc Network has attracted both research and industrial community due to its benefits in facilitating human life and enhancing the security and comfort. However, various issues have been faced in such networks such as information security, routing reliability, dynamic high mobility of vehicles, that influence the stability of communication. To overcome this issue, it is necessary to increase the routing protocols performances, by keeping only the stable path during the communication. The effective solutions that have been investigated in the literature are based on the link prediction to avoid broken links. In this paper, we propose a new solution based on machine learning concept for link prediction, using LR and Support Vector Regression (SVR) which is a variant of the Support Vector Machine (SVM) algorithm. SVR allows predicting the movements of the vehicles in the network which gives us a decision for the link state at a future time. We study the performance of SVR by comparing the generated prediction values against real movement traces of different vehicles in various mobility scenarios, and to show the effectiveness of the proposed method, we calculate the error rate. Finally, we compare this new SVR method with Lagrange interpolation solution.


title: Fine-grained complexity of the graph homomorphism problem for bounded-treewidth graphs

abstract: For graphs $G$ and $H$, a \emph{homomorphism} from $G$ to $H$ is an edge-preserving mapping from the vertex set of $G$ to the vertex set of $H$. For a fixed graph $H$, by \textsc{Hom($H$)} we denote the computational problem which asks whether a given graph $G$ admits a homomorphism to $H$. If $H$ is a complete graph with $k$ vertices, then \textsc{Hom($H$)} is equivalent to the $k$-\textsc{Coloring} problem, so graph homomorphisms can be seen as generalizations of colorings. It is known that \textsc{Hom($H$)} is polynomial-time solvable if $H$ is bipartite or has a vertex with a loop, and NP-complete otherwise [Hell and Ne\v{s}et\v{r}il, JCTB 1990]. In this paper we are interested in the complexity of the problem, parameterized by the treewidth of the input graph $G$. If $G$ has $n$ vertices and is given along with its tree decomposition of width $\mathrm{tw}(G)$, then the problem can be solved in time $|V(H)|^{\mathrm{tw}(G)} \cdot n^{\mathcal{O}(1)}$, using a straightforward dynamic programming. We explore whether this bound can be improved. We show that if $H$ is a \emph{projective core}, then the existence of such a faster algorithm is unlikely: assuming the Strong Exponential Time Hypothesis (SETH), the \textsc{Hom($H$)} problem cannot be solved in time $(|V(H)|-\epsilon)^{\mathrm{tw}(G)} \cdot n^{\mathcal{O}(1)}$, for any $\epsilon > 0$. This result provides a full complexity characterization for a large class of graphs $H$, as almost all graphs are projective cores. We also notice that the naive algorithm can be improved for some graphs $H$, and show a complexity classification for all graphs $H$, assuming two conjectures from algebraic graph theory. In particular, there are no known graphs $H$ which are not covered by our result. In order to prove our results, we bring together some tools and techniques from algebra and from fine-grained complexity.


title: First order covariance inequalities via Stein's method

abstract: We propose probabilistic representations for inverse Stein operators (i.e. solutions to Stein equations) under general conditions; in particular we deduce new simple expressions for the Stein kernel. These representations allow to deduce uniform and non-uniform Stein factors (i.e. bounds on solutions to Stein equations) and lead to new covariance identities expressing the covariance between arbitrary functionals of an arbitrary {univariate} target in terms of a weighted covariance of the derivatives of the functionals. Our weights are explicit, easily computable in most cases, and expressed in terms of objects familiar within the context of Stein's method. Applications of the Cauchy-Schwarz inequality to these weighted covariance identities lead to sharp upper and lower covariance bounds and, in particular, weighted Poincar\'e inequalities. Many examples are given and, in particular, classical variance bounds due to Klaassen, Brascamp and Lieb or Otto and Menz are corollaries. Connections with more recent literature are also detailed.


title: Minimal definable graphs of definable chromatic number at least three

abstract: We show that there is a Borel graph on a standard Borel space of Borel chromatic number three that admits a Borel homomorphism to every analytic graph on a standard Borel space of Borel chromatic number at least three. Moreover, we characterize the Borel graphs on standard Borel spaces of vertex-degree at most two with this property, and show that the analogous result for digraphs fails.


title: Predicting the Voltage Distribution for Low Voltage Networks using Deep Learning

abstract: The energy landscape for the Low-Voltage (LV) networks are beginning to change; changes resulted from the increase penetration of renewables and/or the predicted increase of electric vehicles charging at home. The previously passive `fit-and-forget' approach to LV network management will be inefficient to ensure its effective operations. A more adaptive approach is required that includes the prediction of risk and capacity of the circuits. Many of the proposed methods require full observability of the networks, motivating the installations of smart meters and advance metering infrastructure in many countries. However, the expectation of `perfect data' is unrealistic in operational reality. Smart meter (SM) roll-out can have its issues, which may resulted in low-likelihood of full SM coverage for all LV networks. This, together with privacy requirements that limit the availability of high granularity demand power data have resulted in the low uptake of many of the presented methods. To address this issue, Deep Learning Neural Network is proposed to predict the voltage distribution with partial SM coverage. The results show that SM measurements from key locations are sufficient for effective prediction of voltage distribution.


title: Conformal geodesics on gravitational instantons

abstract: We study the integrability of the conformal geodesic flow (also known as the conformal circle flow) on the $SO(3)$--invariant gravitational instantons. On a hyper--K\"ahler four--manifold the conformal geodesic equations reduce to geodesic equations of a charged particle moving in a constant self--dual magnetic field. In the case of the anti--self--dual Taub NUT instanton we integrate these equations completely by separating the Hamilton--Jacobi equations, and finding a commuting set of first integrals. This gives the first example of an integrable conformal geodesic flow on a four--manifold which is not a symmetric space. In the case of the Eguchi--Hanson we find all conformal geodesics which lie on the three--dimensional orbits of the isometry group. In the non--hyper--K\"ahler case of the Fubini--Study metric on $\CP^2$ we use the first integrals arising from the conformal Killing--Yano tensors to recover the known complete integrability of conformal geodesics.


title: On infinite covariance expansions

abstract: In this paper we provide a probabilistic representation of Lagrange's identity which we use to obtain Papathanasiou-type variance expansions of arbitrary order. Our expansions lead to generalized sequences of weights which depend on an arbitrarily chosen sequence of (non-decreasing) test functions. The expansions hold for arbitrary univariate target distribution under weak assumptions, in particular they hold for continuous and discrete distributions alike. The weights are studied under different sets of assumptions either on the test functions or on the underlying distributions. Many concrete illustrations for standard probability distributions are provided (including Pearson, Ord, Laplace, Rayleigh, Cauchy, and Levy distributions).


title: Consequences of functional equations for pairs of p-adic L-functions

abstract: We prove consequences of functional equations of p-adic L-functions for elliptic curves at supersingular primes p. The results include a relationship between the leading and sub-leading terms (for which we use ideas of Wuthrich and Bianchi), a parity result of orders of vanishing, and invariance of Iwasaswa invariants under conjugate twists of the p-adic L-functions.


title: Fast Computational Convolution Methods For Extended Source Effects In Microlensing Lightcurves

abstract: Extended source effects can be seen in gravitational lensing events when sources cross critical lines. Those events probe the stellar intensity profile and could be used to measure limb darkening coefficients to test stellar model predictions. A data base of accurately measured stellar profiles will be needed to correctly subtract the stellar flux in planetary transient events. The amount of data that is being and will be produced in current and future microlensing surveys, from the space and the ground, requires algorithms that can quickly compute light curves for different source-lens configurations. Based on the convolution method we describe a general formalism to compute those curves for single lenses. We develop approximations in terms of quadratures of elliptic integrals that we integrate by solving the associated first order differential equations. We construct analytic solutions for a limb darkening and, for the first time, for a parabolic profile that are accurate at the $\sim 1-3\%$ and $0.5\%$ level, respectively. These solutions can be computed orders of magnitude faster than other integration routines. They can be implemented in pipelines processing large data sets to extract stellar parameters in real time.


title: Considerations for the Interpretation of Bias Measures of Word Embeddings

abstract: Word embedding spaces are powerful tools for capturing latent semantic relationships between terms in corpora, and have become widely popular for building state-of-the-art natural language processing algorithms. However, studies have shown that societal biases present in text corpora may be incorporated into the word embedding spaces learned from them. Thus, there is an ethical concern that human-like biases contained in the corpora and their derived embedding spaces might be propagated, or even amplified with the usage of the biased embedding spaces in downstream applications. In an attempt to quantify these biases so that they may be better understood and studied, several bias metrics have been proposed. We explore the statistical properties of these proposed measures in the context of their cited applications as well as their supposed utilities. We find that there are caveats to the simple interpretation of these metrics as proposed. We find that the bias metric proposed by Bolukbasi et al. 2016 is highly sensitive to embedding hyper-parameter selection, and that in many cases, the variance due to the selection of some hyper-parameters is greater than the variance in the metric due to corpus selection, while in fewer cases the bias rankings of corpora vary with hyper-parameter selection. In light of these observations, it may be the case that bias estimates should not be thought to directly measure the properties of the underlying corpus, but rather the properties of the specific embedding spaces in question, particularly in the context of hyper-parameter selections used to generate them. Hence, bias metrics of spaces generated with differing hyper-parameters should be compared only with explicit consideration of the embedding-learning algorithms particular configurations.


title: 2D Linear Time-Variant Controller for Human's Intention Detection for Reach-to-Grasp Trajectories in Novel Scenes

abstract: Designing robotic assistance devices for manipulation tasks is challenging. This work is concerned with improving accuracy and usability of semi-autonomous robots, such as human operated manipulators or exoskeletons. The key insight is to develop a system that takes into account context- and user-awareness to take better decisions in how to assist the user. The context-awareness is implemented by enabling the system to automatically generate a set of candidate grasps and reach-to-grasp trajectories in novel, cluttered scenes. The user-awareness is implemented as a linear time-variant feedback controller to facilitate the motion towards the most promising grasp. Our approach is demonstrated in a simple 2D example in which participants are asked to grasp a specific object in a clutter scene. Our approach also reduce the number of controllable dimensions for the user by providing only control on x- and y-axis, while orientation of the end-effector and the pose of its fingers are inferred by the system. The experimental results show the benefits of our approach in terms of accuracy and execution time with respect to a pure manual control.


title: Metrics and Benchmarks for Remote Shared Controllers in Industrial Applications

abstract: Remote manipulation is emerging as one of the key robotics tasks needed in extreme environments. Several researchers have investigated how to add AI components into shared controllers to improve their reliability. Nonetheless, the impact of novel research approaches in real-world applications can have a very slow in-take. We propose a set of benchmarks and metrics to evaluate how the AI components of remote shared control algorithms can improve the effectiveness of such frameworks for real industrial applications. We also present an empirical evaluation of a simple intelligent share controller against a manually operated manipulator in a tele-operated grasping scenario.


title: An Open-World Extension to Knowledge Graph Completion Models

abstract: We present a novel extension to embedding-based knowledge graph completion models which enables them to perform open-world link prediction, i.e. to predict facts for entities unseen in training based on their textual description. Our model combines a regular link prediction model learned from a knowledge graph with word embeddings learned from a textual corpus. After training both independently, we learn a transformation to map the embeddings of an entity's name and description to the graph-based embedding space. In experiments on several datasets including FB20k, DBPedia50k and our new dataset FB15k-237-OWE, we demonstrate competitive results. Particularly, our approach exploits the full knowledge graph structure even when textual descriptions are scarce, does not require a joint training on graph and text, and can be applied to any embedding-based link prediction model, such as TransE, ComplEx and DistMult.


title: Global Convergence of Policy Gradient Methods to (Almost) Locally Optimal Policies

abstract: Policy gradient (PG) methods are a widely used reinforcement learning methodology in many applications such as video games, autonomous driving, and robotics. In spite of its empirical success, a rigorous understanding of the global convergence of PG methods is lacking in the literature. In this work, we close the gap by viewing PG methods from a nonconvex optimization perspective. In particular, we propose a new variant of PG methods for infinite-horizon problems that uses a random rollout horizon for the Monte-Carlo estimation of the policy gradient. This method then yields an unbiased estimate of the policy gradient with bounded variance, which enables the tools from nonconvex optimization to be applied to establish global convergence. Employing this perspective, we first recover the convergence results with rates to the stationary-point policies in the literature. More interestingly, motivated by advances in nonconvex optimization, we modify the proposed PG method by introducing periodically enlarged stepsizes. The modified algorithm is shown to escape saddle points under mild assumptions on the reward and the policy parameterization. Under a further strict saddle points assumption, this result establishes convergence to essentially locally-optimal policies of the underlying problem, and thus bridges the gap in existing literature on the convergence of PG methods. Results from experiments on the inverted pendulum are then provided to corroborate our theory, namely, by slightly reshaping the reward function to satisfy our assumption, unfavorable saddle points can be avoided and better limit points can be attained. Intriguingly, this empirical finding justifies the benefit of reward-reshaping from a nonconvex optimization perspective.

categories: math.OC cs.LG cs.SY eess.SY math.ST stat.TH


title: Endotactic Networks and Toric Differential Inclusions

abstract: An important dynamical property of biological interaction networks is persistence, which intuitively means that "no species goes extinct". It has been conjectured that dynamical system models of weakly reversible networks (i.e., networks for which each reaction is part of a cycle) are persistent. The property of persistence is also related to the well known global attractor conjecture. An approach for the proof of the global attractor conjecture uses an embedding of weakly reversible dynamical systems into toric differential inclusions. We show that the larger class of endotactic dynamical systems can also be embedded into toric differential inclusions. Moreover, we show that, essentially, endotactic networks form the largest class of networks with this property.


title: Complex phase retrieval from subgaussian measurements

abstract: Phase retrieval refers to the problem of reconstructing an unknown vector $x_0 \in \mathbb{C}^n$ or $x_0 \in \mathbb{R}^n $ from $m$ measurements of the form $y_i = \big\vert \langle \xi^{\left(i\right)}, x_0 \rangle \big\vert^2 $, where $ \left\{ \xi^{\left(i\right)} \right\}^m_{i=1} \subset \mathbb{C}^m $ are known measurement vectors. While Gaussian measurements allow for recovery of arbitrary signals provided the number of measurements scales at least linearly in the number of dimensions, it has been shown that ambiguities may arise for certain other classes of measurements $ \left\{ \xi^{\left(i\right)} \right\}^{m}_{i=1}$ such as Bernoulli measurements or Fourier measurements. In this paper, we will prove that even when a subgaussian vector $ \xi^{\left(i\right)} \in \mathbb{C}^m $ does not fulfill a small-ball probability assumption, the PhaseLift method is still able to reconstruct a large class of signals $x_0 \in \mathbb{R}^n$ from the measurements. This extends recent work by Krahmer and Liu from the real-valued to the complex-valued case. However, our proof strategy is quite different and we expect some of the new proof ideas to be useful in several other measurement scenarios as well. We then extend our results $x_0 \in \mathbb{C}^n $ up to an additional assumption which, as we show, is necessary.


title: Inherent Tradeoffs in Learning Fair Representations

abstract: Real-world applications of machine learning tools in high-stakes domains are often regulated to be fair, in the sense that the predicted target should satisfy some quantitative notion of parity with respect to a protected attribute. However, the exact tradeoff between fairness and accuracy is not entirely clear, even for the basic paradigm of classification problems. In this paper, we characterize an inherent tradeoff between statistical parity and accuracy in the classification setting by providing a lower bound on the sum of group-wise errors of any fair classifiers. Our impossibility theorem could be interpreted as a certain uncertainty principle in fairness: if the base rates differ among groups, then any fair classifier satisfying statistical parity has to incur a large error on at least one of the groups. We further extend this result to give a lower bound on the joint error of any (approximately) fair classifiers, from the perspective of learning fair representations. To show that our lower bound is tight, assuming oracle access to Bayes (potentially unfair) classifiers, we also construct an algorithm that returns a randomized classifier that is both optimal (in terms of accuracy) and fair. Interestingly, when the protected attribute can take more than two values, an extension of this lower bound does not admit an analytic solution. Nevertheless, in this case, we show that the lower bound can be efficiently computed by solving a linear program, which we term as the TV-Barycenter problem, a barycenter problem under the TV-distance. On the upside, we prove that if the group-wise Bayes optimal classifiers are close, then learning fair representations leads to an alternative notion of fairness, known as the accuracy parity, which states that the error rates are close between groups. Finally, we also conduct experiments on real-world datasets to confirm our theoretical findings.


title: Experience Replay Optimization

abstract: Experience replay enables reinforcement learning agents to memorize and reuse past experiences, just as humans replay memories for the situation at hand. Contemporary off-policy algorithms either replay past experiences uniformly or utilize a rule-based replay strategy, which may be sub-optimal. In this work, we consider learning a replay policy to optimize the cumulative reward. Replay learning is challenging because the replay memory is noisy and large, and the cumulative reward is unstable. To address these issues, we propose a novel experience replay optimization (ERO) framework which alternately updates two policies: the agent policy, and the replay policy. The agent is updated to maximize the cumulative reward based on the replayed data, while the replay policy is updated to provide the agent with the most useful experiences. The conducted experiments on various continuous control tasks demonstrate the effectiveness of ERO, empirically showing promise in experience replay learning to improve the performance of off-policy reinforcement learning algorithms.


title: Bi-Hamiltonian Systems in (2+1) and Higher Dimensions Defined by Novikov Algebras

abstract: The results from the article [Strachan I.A.B., Szablikowski B.M., Stud. Appl. Math. 133 (2014), 84-117] are extended over consideration of central extensions allowing the introducing of additional independent variables. Algebraic conditions associated to the first-order central extension with respect to additional independent variables are derived. As result $(2+1)$- and, in principle, higher-dimensional multicomponent bi-Hamiltonian systems are constructed. Necessary classification of the central extensions for low-dimensional Novikov algebras is performed and the theory is illustrated by significant $(2+1)$- and $(3+1)$-dimensional examples.


title: Extended range of dipole-dipole interactions in periodically structured photonic media

abstract: The interaction between quantum two-level systems is typically short-range in free space and most photonic environments. Here we show that diminishing momentum isosurfaces with equal frequencies can create a significantly extended range of interaction between distant quantum systems. The extended range is robust and does not rely on a specific location or orientation of the transition dipoles. A general relation between the interaction range and properties of the isosurface is described for structured photonic media. It provides a new way to mediate long-range quantum behavior.


title: Ground state solution for a class of modified nonlinear fourth-order elliptic equation with sign-changing unbounded potential

abstract: We are concerned on the fourth-order elliptic equation \begin{equation}\tag{$P_\lambda$} \left\{ \begin{array}[c]{ll} \Delta^2 u- \Delta u + V(x)u -\lambda \Delta[\rho(u^2)]\rho'(u^2)u= f(u)\, \, \mbox{in} \, \, \mathbb{R}^N, & u\in W^{2,2}(\mathbb{R}^N), \end{array} \right. \end{equation} where $\Delta^2 = \Delta(\Delta)$ is the biharmonic operator, $3\leq N\leq 6$, the radially symmetric potential $V$ may change sign and $\inf_{\mathbb{R}^N}V(x)=-\infty$ is allowed. If $f$ satisfies a type of nonquadracity and monotonicity conditions and $\rho$ is a suitable smooth function, we prove, via variational approach, the existence of a radially symmetric nontrivial ground state solution $u_\lambda$ for problem $(P_\lambda)$ for all $\lambda\geq 0$.


title: Radial Current and Rotation Profile Tailoring in Highly Ionized Linear Plasma Devices

abstract: In a rotating magnetized plasma cylinder with shear, cross-field current can arise from inertial mechanisms and from the cross-field viscosity. Considering these mechanisms, it is possible to calculate the irreducible radial current draw in a cylindrical geometry as a function of the rotation frequency. The resulting expressions raise novel possibilities for tailoring the electric field profile by controlling the density and temperature profiles of a plasma.


title: Key biology you should have learned in physics class: Using ideal-gas mixtures to understand biomolecular machines

abstract: The biological cell exhibits a fantastic range of behaviors, but ultimately these are governed by a handful of physical and chemical principles. Here we explore simple theory, known for decades and based on the simple thermodynamics of mixtures of ideal gases, which illuminates several key functions performed within the cell. Our focus is the free-energy-driven import and export of molecules, such as nutrients and other vital compounds, via transporter proteins. Complementary to a thermodynamic picture is a description of transporters via "mass-action" chemical kinetics, which lends further insights into biological machinery and free energy use. Both thermodynamic and kinetic descriptions can shed light on the fundamental non-equilibrium aspects of transport. On the whole, our biochemical-physics discussion will remain agnostic to chemical details, but we will see how such details ultimately enter a physical description through the example of the cellular fuel ATP.


title: Robust Machine Translation with Domain Sensitive Pseudo-Sources: Baidu-OSU WMT19 MT Robustness Shared Task System Report

abstract: This paper describes the machine translation system developed jointly by Baidu Research and Oregon State University for WMT 2019 Machine Translation Robustness Shared Task. Translation of social media is a very challenging problem, since its style is very different from normal parallel corpora (e.g. News) and also include various types of noises. To make it worse, the amount of social media parallel corpora is extremely limited. In this paper, we use a domain sensitive training method which leverages a large amount of parallel data from popular domains together with a little amount of parallel data from social media. Furthermore, we generate a parallel dataset with pseudo noisy source sentences which are back-translated from monolingual data using a model trained by a similar domain sensitive way. We achieve more than 10 BLEU improvement in both En-Fr and Fr-En translation compared with the baseline methods.


title: Topological electronic structure and intrinsic magnetization in MnBi$_4$Te$_7$: a Bi$_2$Te$_3$-derivative with a periodic Mn sublattice

abstract: Combinations of non-trivial band topology and long-range magnetic order hold promise for realizations of novel spintronic phenomena, such as the quantum anomalous Hall effect and the topological magnetoelectric effect. Following theoretical advances material candidates are emerging. Yet, a compound with a band-inverted electronic structure and an intrinsic net magnetization remains unrealized. MnBi$_2$Te$_4$ is a candidate for the first antiferromagnetic topological insulator and the progenitor of a modular (Bi$_2$Te$_3$)$_n$(MnBi$_2$Te$_4$) series. For $n$ = 1, we confirm a non-stoichiometric composition proximate to MnBi$_4$Te$_7$ and establish an antiferromagnetic state below 13 K followed by a state with net magnetization and ferromagnetic-like hysteresis below 5 K. Angle-resolved photoemission experiments and density-functional calculations reveal a topological surface state on the MnBi$_4$Te$_7$(0001) surface, analogous to the non-magnetic parent compound Bi$_2$Te$_3$. Our results render MnBi$_4$Te$_7$ as a band-inverted material with an intrinsic net magnetization and a complex magnetic phase diagram providing a versatile platform for the realization of different topological phases.


title: Correlation between the Dzyaloshinskii-Moriya interaction and the orbital angular momentum at an oxide / ferromagnet interface

abstract: We report on the Dzyaloshinskii-Moriya (DMI) interaction at the interface between a ferromagnet and an oxide. We demonstrate experimentally that oxides can give rise to DMI. By comparison of systems comprised of Pt/Co90Fe10/Oxide and Cu/Co90Fe10/Oxide, we also show how oxidation of one interface can enhance and add to the total DMI of that generated by the Pt interface. This is due to the fact that the DMI on both interfaces promotes left-handed chirality. Finally, by use of ferromagnetic resonance spectroscopy, we show that the DMI and the spectroscopic splitting factor, which is a measure of the orbital momentum, are correlated. This indicates the importance of hybridization and charge transfer at the oxide interface for the DMI.


title: Universality in Learning from Linear Measurements

abstract: We study the problem of recovering a structured signal from independently and identically drawn linear measurements. A convex penalty function $f(\cdot)$ is considered which penalizes deviations from the desired structure, and signal recovery is performed by minimizing $f(\cdot)$ subject to the linear measurement constraints. The main question of interest is to determine the minimum number of measurements that is necessary and sufficient for the perfect recovery of the unknown signal with high probability. Our main result states that, under some mild conditions on $f(\cdot)$ and on the distribution from which the linear measurements are drawn, the minimum number of measurements required for perfect recovery depends only on the first and second order statistics of the measurement vectors. As a result, the required of number of measurements can be determining by studying measurement vectors that are Gaussian (and have the same mean vector and covariance matrix) for which a rich literature and comprehensive theory exists. As an application, we show that the minimum number of random quadratic measurements (also known as rank-one projections) required to recover a low rank positive semi-definite matrix is $3nr$, where $n$ is the dimension of the matrix and $r$ is its rank. As a consequence, we settle the long standing open question of determining the minimum number of measurements required for perfect signal recovery in phase retrieval using the celebrated PhaseLift algorithm, and show it to be $3n$.


title: Adversarial Self-Paced Learning for Mixture Models of Hawkes Processes

abstract: We propose a novel adversarial learning strategy for mixture models of Hawkes processes, leveraging data augmentation techniques of Hawkes process in the framework of self-paced learning. Instead of learning a mixture model directly from a set of event sequences drawn from different Hawkes processes, the proposed method learns the target model iteratively, which generates "easy" sequences and uses them in an adversarial and self-paced manner. In each iteration, we first generate a set of augmented sequences from original observed sequences. Based on the fact that an easy sample of the target model can be an adversarial sample of a misspecified model, we apply a maximum likelihood estimation with an adversarial self-paced mechanism. In this manner the target model is updated, and the augmented sequences that obey it are employed for the next learning iteration. Experimental results show that the proposed method outperforms traditional methods consistently.


title: A nonlocal game for witnessing quantum networks

abstract: Nonlocal game as a novel witness of the nonlocality of entanglement is of fundamental importance in various fields. The known nonlocal games or equivalent linear Bell inequalities are only useful for Bell networks of single entanglement. Our goal in this paper is to propose a unified method for constructing cooperating games in network scenarios. We first propose an efficient method to construct numerous multipartite games from any graphs. The main idea is the graph representation of entanglement-based quantum networks. We further specify these graphic games with quantum advantages by providing a simple sufficient and necessary condition. The graphic games imply the first linear testing of the nonlocality of general quantum networks consisting of EPR states. It also allows generating new instances going beyond well-known CHSH games. Our result has interesting applications in quantum networks, Bell theory, computational complexity, and theoretical computer science.


title: Sequential Experimental Design for Transductive Linear Bandits

abstract: In this paper we introduce the transductive linear bandit problem: given a set of measurement vectors $\mathcal{X}\subset \mathbb{R}^d$, a set of items $\mathcal{Z}\subset \mathbb{R}^d$, a fixed confidence $\delta$, and an unknown vector $\theta^{\ast}\in \mathbb{R}^d$, the goal is to infer $\text{argmax}_{z\in \mathcal{Z}} z^\top\theta^\ast$ with probability $1-\delta$ by making as few sequentially chosen noisy measurements of the form $x^\top\theta^{\ast}$ as possible. When $\mathcal{X}=\mathcal{Z}$, this setting generalizes linear bandits, and when $\mathcal{X}$ is the standard basis vectors and $\mathcal{Z}\subset \{0,1\}^d$, combinatorial bandits. Such a transductive setting naturally arises when the set of measurement vectors is limited due to factors such as availability or cost. As an example, in drug discovery the compounds and dosages $\mathcal{X}$ a practitioner may be willing to evaluate in the lab in vitro due to cost or safety reasons may differ vastly from those compounds and dosages $\mathcal{Z}$ that can be safely administered to patients in vivo. Alternatively, in recommender systems for books, the set of books $\mathcal{X}$ a user is queried about may be restricted to well known best-sellers even though the goal might be to recommend more esoteric titles $\mathcal{Z}$. In this paper, we provide instance-dependent lower bounds for the transductive setting, an algorithm that matches these up to logarithmic factors, and an evaluation. In particular, we provide the first non-asymptotic algorithm for linear bandits that nearly achieves the information theoretic lower bound.


title: Design of Large Sequential Conformational Change in Mechanical Networks

abstract: From the complex motions of robots to the oxygen binding of hemoglobin, the function of many mechanical systems depends on large, coordinated movements of their components. Such movements arise from a network of physical interactions in the form of links that transmit forces between constituent elements. However, the principled design of specific movements is made difficult by the number and nonlinearity of interactions. Here, we model mechanical systems as linkages of rigid bonds (edges) connected by joints (nodes), and formulate a simple but powerful framework for designing full nonlinear coordinated motions using concepts from dynamical systems theory. We begin with principles for designing finite and infinitesimal motions in small modules, and show that each module is a one-dimensional map between distances across pairs of nodes. Next, we represent the act of combining modules as an iteration of this map, and design networks whose geometries reflect the map's fixed points, limit cycles, and chaos. We use this representation to design different folding sequences from a deployable network and a soliton, to a branched network acting as a mechanical AND gate. Finally, we design large changes in curvature of the entire network, and construct physical networks from laser-cut acrylic, origami, and 3D printed material to demonstrate the framework's potential and versatility for designing the full conformational trajectory of morphing metamaterials and structures.


title: Hierarchical Document Encoder for Parallel Corpus Mining

abstract: We explore using multilingual document embeddings for nearest neighbor mining of parallel data. Three document-level representations are investigated: (i) document embeddings generated by simply averaging multilingual sentence embeddings; (ii) a neural bag-of-words (BoW) document encoding model; (iii) a hierarchical multilingual document encoder (HiDE) that builds on our sentence-level model. The results show document embeddings derived from sentence-level averaging are surprisingly effective for clean datasets, but suggest models trained hierarchically at the document-level are more effective on noisy data. Analysis experiments demonstrate our hierarchical models are very robust to variations in the underlying sentence embedding quality. Using document embeddings trained with HiDE achieves state-of-the-art performance on United Nations (UN) parallel document mining, 94.9% P@1 for en-fr and 97.3% P@1 for en-es.


title: Polarons from first principles, without supercells

abstract: We develop a formalism and a computational method to study polarons in insulators and semi-conductors from first principles. Unlike in standard calculations requiring large supercells, we solve a secular equation involving phonons and electron-phonon matrix elements from density-functional perturbation theory, in a spirit similar to the Bethe-Salpeter equation for excitons. We show that our approach describes seamlessly large and small polarons, and we illustrate its capability by calculating wavefunctions, formation energies, and spectral decomposition of polarons in LiF and Li2O2.


title: How the Avengers assemble: Ecological modelling of effective cast sizes for movies

abstract: The number of characters in a movie is an interesting feature. However, it is non-trivial to measure directly. Naive metrics such as the number of credited characters vary wildly. Here, we show that a metric based on the notion of "ecological diversity" as expressed through a Shannon-entropy based metric can characterise the number of characters in a movie, and is useful in taxonomic classification. We also show how the metric can be generalised using Jensen-Shannon divergence to provide a measure of the similarity of characters appearing in different movies, for instance of use in recommender systems, e.g., Netflix. We apply our measures to the Marvel Cinematic Universe (MCU), and show what they teach us about this highly successful franchise of movies. In particular, these measures provide a useful predictor of "success" for films in the MCU, as well as a natural means to understand the relationships between the stories in the overall film arc.


title: The Structure of the Hadron-Quark Reaction Zone

abstract: Hadron-quark combustion in dense matter is a central topic in the study of phases in compact stars and their high-energy astrophysics. We critically review the literature on hadron-quark combustion, dividing them into a "first wave" that treats the problem as a steady-state burning with or without constraints of mechanical equilibrium, and a "second wave" which uses numerical techniques to resolve the burning front and solves the underlying Partial Differential Equations for the chemistry of the burning front under less restrictive conditions. We detail the inaccuracies that the second wave amends over the first wave, and highlight crucial differences between various approaches in the second wave. We also include results from time-dependent simulations of the reaction zone that include a hadronic EOS, neutrinos, and self-consistent thermodynamics without using parameterized shortcuts.


title: Propagator identities, holographic conformal blocks, and higher-point AdS diagrams

abstract: Conformal blocks are the fundamental, theory-independent building blocks in any CFT, so it is important to understand their holographic representation in the context of AdS/CFT. We describe how to systematically extract the holographic objects which compute higher-point global (scalar) conformal blocks in arbitrary spacetime dimensions, extending the result for the four-point block, known in the literature as a geodesic Witten diagram, to five- and six-point blocks. The main new tools which allow us to obtain such representations are various higher-point propagator identities, which can be interpreted as generalizations of the well-known flat space star-triangle identity, and which compute integrals over products of three bulk-to-bulk and/or bulk-to-boundary propagators in negatively curved spacetime. Using the holographic representation of the higher-point conformal blocks and higher-point propagator identities, we develop geodesic diagram techniques to obtain the explicit direct-channel conformal block decomposition of a broad class of higher-point AdS diagrams in a scalar effective bulk theory, with closed-form expressions for the decomposition coefficients. These methods require only certain elementary manipulations and no bulk integration, and furthermore provide quite trivially a simple algebraic origin of the logarithmic singularities of higher-point tree-level AdS diagrams. We also provide a more compact repackaging in terms of the spectral decomposition of the same diagrams, as well as an independent discussion on the closely related but computationally simpler framework over $p$-adics which admits comparable statements for all previously mentioned results.


title: Tighter monogamy and polygamy relations of multiparty quantum entanglement

abstract: We investigate the tight monogamy and polygamy relations of multiparty entanglement for arbitrary quantum states. By using the power of the bipartite measure of entanglement, we establish a class of tight monogamy relations of multiparty entanglement with larger lower bounds than the existing monogamy relations. We also give a class of tight polygamy relations of multiparty entanglement with smaller upper bounds than the existing polygamy relations, by using the power of the entanglement of assistance. It is shown that these new monogamy and polygamy relations are tighter than the former results.


title: Parameter Enhancement for MELP Speech Codec in Noisy Communication Environment

abstract: In this paper, we propose a deep learning (DL)-based parameter enhancement method for a mixed excitation linear prediction (MELP) speech codec in noisy communication environment. Unlike conventional speech enhancement modules that are designed to obtain clean speech signal by removing noise components before speech codec processing, the proposed method directly enhances codec parameters on either the encoder or decoder side. As the proposed method has been implemented by a small network without any additional processes required in conventional enhancement systems, e.g., time-frequency (T-F) analysis/synthesis modules, its computational complexity is very low. By enhancing the noise-corrupted codec parameters with the proposed DL framework, we achieved an enhancement system that is much simpler and faster than conventional T-F mask-based speech enhancement methods, while the quality of its performance remains similar.


title: Ab initio theory of polarons: formalism and applications

abstract: We develop a theoretical and computational framework to study polarons in semiconductors and insulators from first principles. Our approach provides the formation energy, excitation energy, and wavefunction of both electron and hole polarons, and takes into account the coupling of the electron or hole to all phonons. An important feature of the present method is that it does not require supercell calculations, and relies exclusively on electron band structures, phonon dispersions, and electron-phonon matrix elements obtained from calculations in the crystal unit cell. Starting from the Kohn-Sham (KS) equations of density-functional theory, we formulate the polaron problem as a variational minimization, and we obtain a nonlinear eigenvalue problem in the basis of KS states and phonon eigenmodes. In our formalism the electronic component of the polaron is expressed as a coherent superposition of KS states, in close analogy with the solution of the Bethe-Salpeter equation for the calculation of excitons. We demonstrate the power of the methodology by studying polarons in LiF and Li2O2. We show that our method describes both small and large polarons, and seamlessly captures Frohlich-type polar electron-phonon coupling and non-Frohlich coupling to acoustic and optical phonons. To analyze in quantitative terms the electron-phonon coupling mechanisms leading to the formation of polarons, we introduce spectral decompositions similar to the Eliashberg spectral function. We validate our theory using both analytical results and direct calculations on large supercells. This study constitutes a first step toward complete ab initio many-body calculations of polarons in real materials.


title: Ongoing Vaccine and Monoclonal Antibody HIV Prevention Efficacy Trials and Considerations for Sequel Efficacy Trial Designs

abstract: Four randomized placebo-controlled efficacy trials of a candidate vaccine or passively infused monoclonal antibody for prevention of HIV-1 infection are underway (HVTN 702 in South African men and women; HVTN 705 in sub-Saharan African women; HVTN 703/HPTN 081 in sub-Saharan African women; HVTN 704/HPTN 085 in U.S., Peruvian, Brazilian, and Swiss men or transgender persons who have sex with men). Several challenges are posed to the optimal design of the sequel efficacy trials, including: (1) how to account for the evolving mosaic of effective prevention interventions that may be part of the trial design or standard of prevention; (2) how to define viable and optimal sequel trial designs depending on the primary efficacy results and secondary 'correlates of protection' results of each of the ongoing trials; and (3) how to define the primary objective of sequel efficacy trials if HIV-1 incidence is expected to be very low in all study arms such that a standard trial design has a steep opportunity cost. After summarizing the ongoing trials, I discuss statistical science considerations for sequel efficacy trial designs, both generally and specifically to each trial listed above. One conclusion is that the results of 'correlates of protection' analyses, which ascertain how different host immunological markers and HIV-1 viral features impact HIV-1 risk and prevention efficacy, have an important influence on sequel trial design. This influence is especially relevant for the monoclonal antibody trials because of the focused pre-trial hypothesis that potency and coverage of serum neutralization constitutes a surrogate endpoint for HIV-1 infection... (see manuscript for the full abstract)


title: Optimal Reinsurance and Investment Strategies under Mean-Variance Criteria: Partial and Full Information

abstract: This paper is concerned with an optimal reinsurance and investment problem for an insurance firm under the criterion of mean-variance. The driving Brownian motion and the rate in return of the risky asset price dynamic equation cannot be directly observed. And the short-selling of stocks is prohibited. The problem is formulated as a stochastic linear-quadratic (LQ) optimal control problem where the control variables are constrained. Based on the separation principle and stochastic filtering theory, the partial information problem is solved. Efficient strategies and efficient frontier are presented in closed forms via solutions to two extended stochastic Riccati equations. As a comparison, the efficient strategies and efficient frontier are given by the viscosity solution for the Hamilton-Jacobi-Bellman (HJB) equation in the full information case. Some numerical illustrations are also provided.


title: A novel linear battery energy storage system (BESS) life loss calculation model for BESS-integrated wind farm in scheduled power tracking

abstract: Recently, rapid development of battery technology makes it feasible to integrate renewable generations with battery energy storage system (BESS). The consideration of BESS life loss for different BESS application scenarios is economic imperative. In this paper, a novel linear BESS life loss calculation model for BESS-integrated wind farm in scheduled power tracking is proposed. Firstly, based on the life cycle times-depth of discharge (DOD) relation-curve, the BESS life loss coefficient for unit throughput energy with different state of charge (SOC) can be determined from the life cycle times-DOD relation-curve fitting function directly. Secondly, as unidirectional variation of SOC in a single time step, the BESS life loss can be calculated through integration of the life loss coefficient-SOC relation function. A linear BESS life loss calculation model is established through self-optimal piecewise linearization of the primitive function of the life loss coefficient-SOC relation function. Thirdly, the proposed life loss calculation model is incorporated in the BESS-integrated wind farm scheduled power tracking optimization. Case studies demonstrate that with the proposed method, the BESS life loss item can be incorporated in the optimization model effectively, and the scheduled power tracking cost of the BESS-integrated wind farm can be determined and optimized more comprehensively.


title: Data Interpolating Prediction: Alternative Interpretation of Mixup

abstract: Data augmentation by mixing samples, such as Mixup, has widely been used typically for classification tasks. However, this strategy is not always effective due to the gap between augmented samples for training and original samples for testing. This gap may prevent a classifier from learning the optimal decision boundary and increase the generalization error. To overcome this problem, we propose an alternative framework called Data Interpolating Prediction (DIP). Unlike common data augmentations, we encapsulate the sample-mixing process in the hypothesis class of a classifier so that train and test samples are treated equally. We derive the generalization bound and show that DIP helps to reduce the original Rademacher complexity. Also, we empirically demonstrate that DIP can outperform existing Mixup.


title: Testing SU(3) Flavor Symmetry in Semileptonic and Two-body Nonleptonic Decays of Hyperons

abstract: The semileptonic decays and two-body nonleptonic decays of light baryon octet ($T_8$) and decuplet ($T_{10}$) consisting of light $u,d,s$ quarks are studied with the SU(3) flavor symmetry in this work. We obtain the amplitude relations between different decay modes by the SU(3) irreducible representation approach, and then predict relevant branching ratios by present experimental data within $1 \sigma$ error. We find that the predictions for all branching ratios except $\mathcal{B}(\Xi\rightarrow \Lambda^0\pi)$ and $\mathcal{B}(\Xi^*\rightarrow \Xi\pi)$ are in good agreement with present experimental data, that implies the neglected $C_+$ terms or SU(3) breaking effects might contribute at the order of a few percent in $\Xi\rightarrow \Lambda^0\pi$ and $\Xi^*\rightarrow \Xi\pi$ weak decays. We predict that $\mathcal{B}(\Xi^{-}\rightarrow \Sigma^0\mu^-\bar{\nu}_\mu)=(1.13\pm0.08)\times10^{-6}$, $\mathcal{B}(\Xi^{-}\rightarrow\Lambda^0\mu^-\bar{\nu}_\mu)=(1.58\pm0.04)\times10^{-4}$, $\mathcal{B}(\Omega^-\rightarrow\Xi^0\mu^-\bar{\nu}_\mu)=(3.7\pm1.8)\times10^{-3}$, $\mathcal{B}(\Sigma^-\rightarrow \Sigma^0e^-\bar{\nu}_e)=(1.35\pm0.28)\times10^{-10}$, $\mathcal{B}(\Xi^-\rightarrow \Xi^0e^-\bar{\nu}_e)=(4.2\pm2.4)\times10^{-10}$. We also study $T_{10}\to T_8 P_8$ weak, electromagnetic or strong decays. Some of these decay modes could be observed by the BESIII, LHCb and other experiments in the near future. Due to the very small life times of $\Sigma^0$, $\Xi^{*0,-}$, $\Sigma^{*0,-}$ and $\Delta^{0,-}$, the branching ratios of these baryon weak decays are only at the order of $\mathcal{O}(10^{-20}-10^{-13}$), which are too small to be reached by current experiments. Furthermore, the longitudinal branching ratios of $T_{8A} \to T_{8B} \ell^- \bar{\nu}_\ell~(\ell=\mu,e)$ decays are also given.


title: Stable homotopy, 1-dimensional NCCW complexes, and Property (H)

abstract: In this paper, we show that the homomorphisms between two unital one-dimensional NCCW complexes with the same KK-class are stably homotopic, i.e., with adding on a common homomorphism (with finite dimensional image), they are homotopic. As a consequence, any one-dimensional NCCW complex has the Property (H).


title: A Monaural Speech Enhancement Method for Robust Small-Footprint Keyword Spotting

abstract: Robustness against noise is critical for keyword spotting (KWS) in real-world environments. To improve the robustness, a speech enhancement front-end is involved. Instead of treating the speech enhancement as a separated preprocessing before the KWS system, in this study, a pre-trained speech enhancement front-end and a convolutional neural networks (CNNs) based KWS system are concatenated, where a feature transformation block is used to transform the output from the enhancement front-end into the KWS system's input. The whole model is trained jointly, thus the linguistic and other useful information from the KWS system can be back-propagated to the enhancement front-end to improve its performance. To fit the small-footprint device, a novel convolution recurrent network is proposed, which needs fewer parameters and computation and does not degrade performance. Furthermore, by changing the input features from the power spectrogram to Mel-spectrogram, less computation and better performance are obtained. our experimental results demonstrate that the proposed method significantly improves the KWS system with respect to noise robustness.


title: Improving the robustness of ImageNet classifiers using elements of human visual cognition

abstract: We investigate the robustness properties of image recognition models equipped with two features inspired by human vision, an explicit episodic memory and a shape bias, at the ImageNet scale. As reported in previous work, we show that an explicit episodic memory improves the robustness of image recognition models against small-norm adversarial perturbations under some threat models. It does not, however, improve the robustness against more natural, and typically larger, perturbations. Learning more robust features during training appears to be necessary for robustness in this second sense. We show that features derived from a model that was encouraged to learn global, shape-based representations (Geirhos et al., 2019) do not only improve the robustness against natural perturbations, but when used in conjunction with an episodic memory, they also provide additional robustness against adversarial perturbations. Finally, we address three important design choices for the episodic memory: memory size, dimensionality of the memories and the retrieval method. We show that to make the episodic memory more compact, it is preferable to reduce the number of memories by clustering them, instead of reducing their dimensionality.


title: A parallax distance to 3C 273 through spectroastrometry and reverberation mapping

abstract: Distance measurements for extragalactic objects are a fundamental problem in astronomy and cosmology. In the era of precision cosmology, we urgently need better measurements of cosmological distances to observationally test the increasing $H_{0}$ tension of the Hubble constant measured from different tools. Using spectroastrometry, GRAVITY at The Very Large Telescope Interferometer successfully revealed the structure, kinematics and angular sizes of the broad-line region (BLR) of 3C 273 with an unprecedentedly high spatial resolution. Fortunately, reverberation mapping (RM) of active galactic nuclei (AGNs) reliably provides linear sizes of their BLRs. Here we report a joint analysis of spectroastrometry and RM observations to measure AGN distances. We apply this analysis to 3C 273 observed by both GRAVITY and an RM campaign,and find an angular distance of $551.5_{-78.7}^{+97.3}\, {\rm Mpc}$ and $H_{0}=71.5_{-10.6}^{+11.9}\,{\rm km\,s^{-1}\,Mpc^{-1}}$. Advantages of the analysis are 1) its pure geometrical measurements and 2) it simultaneously yields mass of the central black hole in the BLR. Moreover, we can conveniently repeat measurements of selected AGNs to efficiently reduce the statistical and systematic errors. Future observations of a reasonably sized sample ($\sim 30$ AGNs) will provide distances of the AGNs and hence a new way of measuring $H_{0}$ with a high precision $\left(\lesssim 3\%\right)$ to test the $H_{0}$ tension.


title: Gridless Multisnapshot Variational Line Spectral Estimation from Coarsely Quantized Samples

abstract: Due to the increasing demand for low power and higher sampling rates, low resolution quantization for data acquisition has drawn great attention recently. Consequently, line spectral estimation (LSE) with multiple measurement vectors (MMVs) from coarsely quantized samples is of vital importance in cutting edge array signal processing applications such as range estimation and DOA estimation in millimeter wave radar systems. In this paper, we combine the low complexity gridless variational line spectral estimation (VALSE) and expectation propagation (EP) and propose an MVALSE-EP algorithm to estimate the frequencies from coarsely quantized samples. In addition, the Cram\'{e}r Rao bound (CRB) is derived as a benchmark performance of the proposed algorithm, and insights are provided to reveal the effects of system parameters on estimation performance. It is shown that snapshots benefits the frequency estimation, especially in coarsely quantized scenarios. Numerical experiments are conducted to demonstrate the effectiveness of MVALSE-EP, including real data set.


title: The Dark Matter Profiles in the Milky Way

abstract: We investigate the dark matter profile of the Milky Way using the observed rotation curve data out to 100 kpc. The baryonic matter of the Milky Way is divided into bulge, disk and gas components, and each component is modelled using various possible mass profiles available in literature. The arbitrary combination of seven bulge profiles, four disk profiles and two gas profiles results in fifty-six baryon models. These baryon models are combined with one of the four dark matter profiles: Burkert profile, core-modified profile, pseudo-isothermal profile and NFW profile, to fit the observed rotation curve data. Results show that in general the NFW profile fits the data better than the Burkert profile, while the core-modified profile and the pseudo-isothermal profile are essentially ruled out. The best-fitting NFW model has the scale length $r_0=8.1\pm 0.7$ kpc, and the corresponding local density of dark matter is $\rho_{\rm dm}(R=R_\odot)=0.51\pm0.09~{\rm GeV}/{\rm cm}^3$.


title: Causal Inference from Possibly Unbalanced Split-Plot Designs: A Randomization-based Perspective

abstract: Split-plot designs find wide applicability in multifactor experiments with randomization restrictions. Practical considerations often warrant the use of unbalanced designs. This paper investigates randomization based causal inference in split-plot designs that are possibly unbalanced. Extension of ideas from the recently studied balanced case yields an expression for the sampling variance of a treatment contrast estimator as well as a conservative estimator of the sampling variance. However, the bias of this variance estimator does not vanish even when the treatment effects are strictly additive. A careful and involved matrix analysis is employed to overcome this difficulty, resulting in a new variance estimator, which becomes unbiased under milder conditions. A construction procedure that generates such an estimator with minimax bias is proposed.


title: Reliable data from low cost ozone sensors in a hierarchical network

abstract: We demonstrate how a hierarchical network comprising a number of compliant reference stations and a much larger number of low-cost sensors can deliver reliable high temporal-resolution ozone data at neighbourhood scales. The framework, demonstrated originally for a smaller scale regional network deployed in the Lower Fraser Valley, BC was tested and refined using two much more extensive networks of gas-sensitive semiconductor-based (GSS) sensors deployed at neighbourhood scales in Los Angeles: one of ~20 and one of ~45 GSS ozone sensors. Of these, ten sensors were co-located with different regulatory measurement stations, allowing a rigorous test of the accuracy of the algorithms used for off-site calibration and adjustment of low cost sensors. The method is based on adjusting the gain and offset of the low-cost sensor to match the first two moments of the probability distribution of the sensor result to that of a proxy: a calibrated independent measurement (usually derived from regulatory monitors) whose probability distribution evaluated over a time that emphasizes diurnal variations is similar to that at the test location. The regulatory measurement station physically closest to the low-cost sensor was a good proxy for most sites. The algorithms developed were successful in detecting and correcting sensor drift, and in identifying locations where geographical features resulted in significantly different patterns of ozone variation due to the relative dominance of different dispersion, emission and chemical processes. The entire network results show very large variations in ozone concentration that take place on short time- and distance scales across the Los-Angeles region. Such patterns were not captured by the more sparsely distributed stations of the existing regulatory network and demonstrate the need for reliable data from dense networks of monitors.


title: Van der Waals stacked multilayer in-plane graphene/hexagonal boron nitride heterostructure: its interfacial thermal transport properties

abstract: Combining both vertical and in-plane two-dimensional (2D) heterostructures opens up the possibility to create an unprecedented architecture using 2D atomic layer building blocks. The thermal transport properties of such mixed heterostructures, critical to various applications in nanoelectronics, however, have not been thoroughly explored. Herein, we construct two configurations of multilayer in-plane graphene/hexagonal boron nitride (Gr/h-BN) heterostructures (i.e. mixed heterostructures) via weak van der Waals (vdW) interactions and systematically investigate the dependence of their interfacial thermal conductance (ITC) on the number of layers using non-equilibrium molecular dynamics (NEMD) simulations. The computational results show that the ITC of two configurations of multilayer in-plane Gr/h-BN heterostructures (MIGHHs) decrease with increasing layer number n and both saturate at n = 3. And surprisingly, we find that the MIGHH is more advantageous to interfacial thermal transport than the monolayer in-plane Gr/h-BN heterostructure, which is in strong contrast to the commonly held notion that the multilayer structures of Gr and h-BN suppress the phonon transmission. The underlying physical mechanisms for these puzzling phenomena are probed through the analyses of heat flux, temperature jump, stress concentration factor, overlap of phonon vibrational spectra and phonon participation ratio. In particular, by changing the stacking angle of MIGHH, a higher ITC can be obtained due to the thermal rectification behavior. Furthermore, we find that the ITC in MIGHH can be well-regulated by controlling the coupling strength between layers. Our findings here are of significance for understanding the interfacial thermal transport behaviors of multilayer in-plane Gr/h-BN heterostructure, and are expected to attract extensive interest in exploring its new physics and applications.


title: The existence of optimal control for continuous-time Markov decision processes in random environments

abstract: In this work, we investigate the optimal control problem for continuous-time Markov decision processes with the random impact of the environment. We provide conditions to show the existence of optimal controls under finite-horizon criteria. Under appropriate conditions, the value function is continuous and satisfies the dynamic programming principle. These results are established by introducing some restriction on the regularity of the optimal controls and by developing a new compactification method for continuous-time Markov decision processes, which is originally used to solve the optimal control problem for jump-diffusion processes.


title: Cryptanalysis of Khatoon et al.'s ECC-based Authentication Protocol for Healthcare Systems

abstract: Telecare medical information systems are gaining rapid popularity in terms of providing the delivery of online health-related services such as online remote health profile access for patients and doctors. Due to being installed entirely on Internet, these systems are exposed to various security and privacy threats. Hence, establishing a secure key agreement and authentication process between the patients and the medical servers is an important challenge. Recently, Khatoon et.al proposed an ECC-based unlink-able authentication and key agreement method for healthcare related application in smart city. In this article, we provide a descriptive analysis on their proposed scheme and prove that Khatoon et al.'s scheme is vulnerable to known-session-specific temporary information attack and is not able to provide perfect forward secrecy.


title: The existence of optimal feedback controls for stochastic dynamical systems with regime-switching

abstract: In this work we provide explicit conditions on the existence of optimal feedback controls for stochastic processes with regime-switching. We use the compactification method which needs less regularity conditions on the coefficients of the studied stochastic systems. Two different kinds of controls have been considered: one is the control on the coefficients of the diffusion processes, another is the control on the transition rate matrices of the continuous-time Markov chains. Moreover, the dynamic programming principle is established after showing the continuity of the value function.