title: Primordial power spectrum and cosmology from black-box galaxy surveys

abstract: We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computation using a novel effective likelihood, and upon the linearisation of black-box models around an expansion point. Consequently, we obtain simple "filter equations" for an effective posterior of the primordial power spectrum, and a straightforward scheme for cosmological parameter inference. We demonstrate that the workload is computationally tractable, fixed a priori, and perfectly parallel. As a proof of concept, we apply our framework to a realistic synthetic galaxy survey, with a data model accounting for physical structure formation and incomplete and noisy galaxy observations. In doing so, we show that the use of non-linear numerical models allows the galaxy power spectrum to be safely fitted up to at least $k_\mathrm{max} = 0.5$ $h$/Mpc, outperforming state-of-the-art backward-modelling techniques by a factor of $\sim 5$ in the number of modes used. The result is an unbiased inference of the primordial matter power spectrum across the entire range of scales considered, including a high-fidelity reconstruction of baryon acoustic oscillations. It translates into an unbiased and robust inference of cosmological parameters. Our results pave the path towards easy applications of likelihood-free simulation-based inference in cosmology. We have made our code pySELFI and our data products publicly available at http://pyselfi.florent-leclercq.eu.


title: Quark mass dependence of $\gamma^{*}\pi\rightarrow\pi\pi$

abstract: Usually the simulation of scattering processes in lattice QCD is carried out at unphysically high values of the quark masses. Hence, a method to extrapolate data obtained in lattice calculations to physical masses is needed to allow for comparison between theory and experiment. To obtain a sound extrapolation, dispersion relations and chiral perturbation theory can be invoked. While a simple combined approach known as the inverse amplitude method allows for a successful extrapolation of $\pi\pi\to\pi\pi$ data, a more complicated framework is needed for inelastic processes such as $\gamma^{*}\pi\to\pi\pi$. By employing a well-established dispersive description, the extrapolation can be performed for $\gamma^{*}\pi\to\pi\pi$ both for on-shell as well as virtual photons, the decay $\gamma^*\rightarrow\pi\pi\pi$ is also within the range of applicability. This particular process is interesting due to both its contribution to the anomalous magnetic moment of the muon and its connection to the axial anomaly.


title: Detectable molecular features above hydrocarbon haze via transmission spectroscopy with JWST: Case studies of GJ 1214b-, GJ 436b-, HD 97658b-, and Kepler-51b-like planets

abstract: Some of the exoplanets so far observed show featureless or flat transmission spectra, possibly indicating the existence of clouds and/or haze in their atmospheres. Thanks to its large aperture size and broad wavelength coverage, JWST is expected to enable detailed investigation of exoplanet atmospheres, which could provide important constraints on the atmospheric composition obscured by clouds/haze. Here, we use four warm ($\lesssim 1000$ K) planets suitable for atmospheric characterization via transmission spectroscopy, GJ 1214b, GJ 436b, HD 97658b, and Kepler-51b, as examples to explore molecular absorption features detectable by JWST even in the existence of hydrocarbon haze in the atmospheres. We simulate photochemistry, the growth of hydrocarbon haze particles, and transmission spectra for the atmospheres of these four planets. Among the planetary parameters considered, super-Earths with hazy, relatively hydrogen-rich atmospheres are mostly expected to produce detectable molecular absorption features such as a quite prominent $\mathrm{CH_4}$ feature at 3.3 ${\rm \mu}$m even for the extreme case of the most efficient production of photochemical haze. For a planet that has extremely low gravity, such as Kepler-51b, haze particles grow significantly large in the upper atmosphere due to the small sedimentation velocity, resulting in the featureless or flat transmission spectrum in a wide wavelength range. This investigation shows that the transmission spectra with muted features measured by HST in most cases do not preclude strong features at the longer wavelengths accessible by JWST.


title: Dust attenuation, dust emission, and dust temperature in galaxies at z>=5: a view from the FIRE-2 simulations

abstract: We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of halos up to M_halo~10^12 M_sun (M_star~10^10.5 M_sun) at z>=5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust extinction, dust emission, and dust temperature within these simulated z>=5 galaxies. Our sample forms a tight correlation between infrared excess (IRX=F_IR/F_UV) and ultraviolet (UV)-continuum slope (beta_UV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX-beta_UV relation is mainly determined by the shape of the extinction curve and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (L_IR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given L_IR, the peak wavelength of the dust spectral energy distributions for z>=5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z=0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (LFs) at z=5-10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IRLFs and UV luminosity densities at z=5-10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z>=5 galaxies.


title: Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores

abstract: Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.


title: Lattice models for Non-Fermi Liquids with Tunable Transport Scalings

abstract: A variety of exotic non-fermi liquid (NFL) states have been observed in many condensed matter systems, with different scaling relations between transport coefficients and temperature. The "standard" approach to studying these NFLs is by coupling a Fermi liquid to quantum critical fluctuations, which potentially can drive the system into a NFL. In this work we seek for an alternative understanding of these various NFLs in a unified framework. We first construct two "elementary" randomness-free models with four-fermion interactions only, whose many properties can be analyzed exactly in certain limit just like the Sachdev-Ye-Kitaev (SYK) model. The most important new feature of our models is that, the fermion scaling dimension in the conformal invariant solution in the infrared limit is tunable by charge density. Then based on these elementary models, we propose two versions of lattice models with four fermion interactions which give us non-fermi liquid behaviors with DC resistivity scaling $\varrho \sim T^\alpha$ in a finite temperature window, and $\alpha \in [1, 2)$ depends on the fermion density in the model, which is a rather universal feature observed in many experimental systems.


title: Joint Bayesian analysis of large angular scale CMB temperature anomalies

abstract: Cosmic microwave background measurements show an agreement with the concordance cosmology model except for a few notable anomalies: Power Suppression, the lack of large scale power in the temperature data compared to what is expected in the concordance model, and Cosmic Hemispherical Asymmetry, a dipolar breakdown of statistical isotropy. An expansion of the CMB covariance in Bipolar Spherical Harmonics naturally parametrizes both these large-scale anomalies, allowing us to perform an exhaustive, fully Bayesian joint analysis of the power spectrum and violations of statistical isotropy up to the dipole level. Our analysis sheds light on the scale dependence of the Cosmic Hemispherical Asymmetry. Assuming a scale-dependent dipole modulation model with a two-parameter power law form, we explore the posterior pdf of amplitude $A(l = 16)$ and the power law index $\alpha$ and find the maximum a posteriori values $A_*(l = 16) = 0.064 \pm 0.022$ and $\alpha_* = -0.92 \pm 0.22$. The maximum a posteriori direction associated with the Cosmic Hemispherical Asymmetry is $(l,b) = (247.8^o, -19.6^o)$ in Galactic coordinates, consistent with previous analyses. We evaluate the Bayes factor $B_{SI-DM}$ to compare the Cosmic Hemispherical Asymmetry model with the isotropic model. The data prefer but do not substantially favor the anisotropic model ($B_{SI-DM}=0.4$). We consider several priors and find that this evidence ratio is robust to prior choice. The large-scale power suppression does not soften when jointly inferring both the isotropic power spectrum and the parameters of the asymmetric model, indicating no evidence that these anomalies are coupled.


title: Searching new physics in rare $B$-meson decays into multiple muons

abstract: New heavy vector bosons and light scalars are predicted in a plethora of models of new physics. In particular, in new strongly interacting sectors they play the role of the $\rho$ and $\pi$ mesons in QCD. We show that some of their interactions, for example those required for the explanation of the $B$ anomalies and the $g-2$ of the muon, can be only probed in $B$ meson decays. We highlight new golden channels not yet studied experimentally, including $B^+ \rightarrow K^+(D^+)\mu^+\mu^-\mu^+\mu^-$ and $B^0\rightarrow K^{*0}\mu^+\mu^-\mu^+\mu^-$. Relying on generator level simulations for data taking with the LHCb detector, we determine the reach of this facility to the aforementioned processes. We show that branching ratios as small as $9\times 10^{-12}$ ($3.2\times 10^{-10}$) and $2.7\times 10^{-11}$ can be tested at the $95\, \%$ CL respectively.


title: Holographic Tensor Networks in Full AdS/CFT

abstract: We present a general procedure for constructing tensor networks for geometric states in the Anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence. Given a state in a large-$N$ CFT with a static, semiclassical gravitational dual, our procedure produces a tensor network for the boundary state whose internal geometry matches (a discretization of) the bulk spacetime geometry. By invoking the "holographic entanglement of purification" conjecture, our construction can be made to capture the structure of the bulk spacetime at sub-AdS scales.


title: The dichotomy of dark matter fraction and total mass density slope of galaxies over five dex in mass

abstract: We analyse the mass density distribution in the centres of galaxies across five orders of magnitude in mass range. Using high-quality spiral galaxy rotation curves and infrared photometry from SPARC, we conduct a systematic study of their central dark matter fraction ($f_{\rm DM}$) and their mass density slope ($\alpha$), within their effective radius. We show that lower-mass spiral galaxies are more dark matter dominated and have more shallow mass density slopes when compared with more massive galaxies, which have density profiles closer to isothermal. Low-mass ($M_{\star} \lesssim 10^{10}\, \rm M_{\odot}$) gas-rich spirals span a wide range of \fdm\ values, but systematically lower than in gas-poor systems of similar mass. With increasing galaxy mass, the values of \fdm\ decrease and the density profiles steepen. In the most massive late-type gas-poor galaxies, a possible flattening of these trends is observed. When comparing these results to massive ($M_{\star} \gtrsim 10^{10}M_{\odot}$) elliptical galaxies from SPIDER and to dwarf ellipticals from SMACKED, these trends result to be inverted. Hence, the values of both $f_{\rm DM}$ and $\alpha$, as a function of $M_{\star}$, exhibit a U-shape trend. At a fixed stellar mass, the mass density profiles in dwarf ellipticals are steeper than in spirals. These trends can be understood by stellar feedback from a more prolonged star formation period in spirals, causing a transformation of the initial steep density cusp to a more shallow profile via differential feedback efficiency by supernovae, and by galaxy mergers or AGN feedback in higher-mass galaxies.


title: The Role of Machine Learning in the Next Decade of Cosmology

abstract: In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.


title: Chromatic Adaptation Transform by Spectral Reconstruction (Preprint)

abstract: A color appearance model (CAM) is an advanced colorimetric tool used to predict color appearance under a wide variety of viewing conditions. A chromatic adaptation transform (CAT) is an integral part of a CAM. Its role is to predict "corresponding colors," that is, a pair of colors that have the same color appearance when viewed under different illuminants, after partial or full adaptation to each illuminant. Modern CATs perform well when applied to a limited range of illuminant pairs and a limited range of source (test) colors. However, they can fail if operated outside these ranges. For imaging applications, it is important to have a CAT that can operate on any real color and illuminant pair without failure. This paper proposes a new CAT that does not operate on the standard von Kries model of adaptation. Instead it relies on spectral reconstruction and how these reconstructions behave with respect to different illuminants. It is demonstrated that the proposed CAT is immune to some of the limitations of existing CATs (such as producing colors with negative tristimulus values). The proposed CAT does not use established empirical corresponding-color datasets to optimize performance, as most modern CATs do, yet it performs as well as or better than the most recent CATs when tested against the corresponding-color datasets. This increase in robustness comes at the expense of additional complexity and computational effort. If robustness is of prime importance, then the proposed method may be justifiable.


title: A geodesic Witten diagram description of holographic entanglement entropy and its quantum corrections

abstract: We use the formalism of geodesic Witten diagrams to study the holographic realization of the conformal block expansion for entanglement entropy of two disjoint intervals. The agreement between the Ryu-Takayanagi formula and the identity block contribution has a dual realization as the product of bulk to boundary propagators. Quantum bulk corrections instead arise from stripped higher order diagrams and back-reaction effects; these are also mapped to the structure for $G_N^0$ terms found in \cite{Faulkner:2013ana}, with the former identified as the bulk entanglement entropy across the Ryu-Takayanagi surfaces. An independent derivation of this last statement is provided by implementing a twist-line formalism in the bulk, and additional checks from the computation of mutual information and single interval entanglement entropy. Finally an interesting correspondence is found between the recently proposed holographic entanglement of purification, and an approximated form for certain $1/c$ Renyi entropies corrections.


title: Coloring Big Graphs with AlphaGoZero

abstract: We show that recent innovations in deep reinforcement learning can effectively color very large graphs -- a well-known NP-hard problem with clear commercial applications. Because the Monte Carlo Tree Search with Upper Confidence Bound algorithm used in AlphaGoZero can improve the performance of a given heuristic, our approach allows deep neural networks trained using high performance computing (HPC) technologies to transform computation into improved heuristics with zero prior knowledge. Key to our approach is the introduction of a novel deep neural network architecture (FastColorNet) that has access to the full graph context and requires $O(V)$ time and space to color a graph with $V$ vertices, which enables scaling to very large graphs that arise in real applications like parallel computing, compilers, numerical solvers, and design automation, among others. As a result, we are able to learn new state of the art heuristics for graph coloring.


title: Weighted games of best choice

abstract: The game of best choice (also known as the secretary problem) is a model for sequential decision making with a long history and many variations. The classical setup assumes that the sequence of candidate rankings are uniformly distributed. Given a statistic on the symmetric group, one can instead weight each permutation according to an exponential function in the statistic. We play the game of best choice on the Ewens and Mallows distributions that are obtained in this way from the number of left-to-right maxima and number of inversions in the permutation, respectively. For each of these, we give the optimal strategy and probability of winning. Moreover, we introduce a general class of permutation statistics that always produces games of best choice whose optimal strategies are positional, which simplifies their analysis considerably.


title: Gravitational wave echoes from black hole area quantization

abstract: Gravitational-wave astronomy has the potential to substantially advance our knowledge of the cosmos, from the most powerful astrophysical engines to the initial stages of our universe. Gravitational waves also carry information about the nature of black holes. Here we investigate the potential of gravitational-wave detectors to test a proposal by Bekenstein and Mukhanov that the area of black hole horizons is quantized in units of the Planck area. Our results indicate that this quantization could have a potentially observable effect on the classical gravitational wave signals received by detectors. In particular, we find distorted gravitational-wave "echoes" in the post-merger waveform describing the inspiral and merger of two black holes. These echoes have a specific frequency content that is characteristic of black hole horizon area quantization.


title: Diameter of P.A. random graphs with edge-step functions

abstract: In this work we prove general bounds for the diameter of random graphs generated by a preferential attachment model whose parameter is a function $f:\mathbb{N}\to[0,1]$ that drives the asymptotic proportion between the numbers of vertices and edges. These results are sharp when $f$ is a \textit{regularly varying function at infinity} with strictly negative index of regular variation~$-\gamma$. For this particular class, we prove a characterization for the diameter that depends only on~$-\gamma$. More specifically, we prove that the diameter of such graphs is of order $1/\gamma$ with high probability, although its vertex set order goes to infinity polynomially. Sharp results for the diameter for a wide class of \textit{slowly varying functions} are also obtained.


title: Electromagnetic multipole moments of baryons

abstract: We calculate the charge quadrupole and magnetic octupole moments of baryons using a group theoretical approach based on broken SU(6) spin-flavor symmetry. The latter is an approximate symmetry of the QCD Lagrangian which becomes exact in the large color N_c limit. Spin-flavor symmetry breaking is induced by one-, two-, and three-quark terms in the electromagnetic current operator. Two- and three-quark currents provide the leading contributions for higher multipole moments, despite being of higher order in an 1/N_c expansion. Our formalism leads to relations between N --> N* transition multipole moments and nucleon ground state properties. We compare our results to experimental quadrupole and octupole transition moments extracted from measured helicity amplitudes.


title: One-loop weak corrections to Higgs production

abstract: We compute mixed QCD-weak corrections to inclusive Higgs production at the LHC from the partonic process $g g \rightarrow H q \bar{q}$. We start from the UV- and IR-finite one-loop weak amplitude and consider its interference with the corresponding one-loop QCD amplitude. This contribution is a $\mathcal{O}(\alpha_s\alpha)$ correction to the leading-order gluon-fusion cross section, and was not numerically assessed in previous works. We also compute the cross section from the square of this weak amplitude, suppressed by $\mathcal{O}(\alpha^2)$. Finally, we consider contributions from the partonic process $g q \rightarrow H q$, which are one order lower in $\alpha_s$, as a reference for the size of terms which are not enhanced by the large gluon luminosity. We find that, given the magnitude of the uncertainties on current state-of-the-art predictions for Higgs production, all contributions computed in this work can be safely ignored, both fully inclusively and in the boosted Higgs regime. This result supports the approximate factorisation of QCD and weak corrections to that process.


title: Fast Approximation of Frequent $k$-mers and Applications to Metagenomics

abstract: Estimating the abundances of all $k$-mers in a set of biological sequences is a fundamental and challenging problem with many applications in biological analysis. While several methods have been designed for the exact or approximate solution of this problem, they all require to process the entire dataset, that can be extremely expensive for high-throughput sequencing datasets. While in some applications it is crucial to estimate all $k$-mers and their abundances, in other situations reporting only frequent $k$-mers, that appear with relatively high frequency in a dataset, may suffice. This is the case, for example, in the computation of $k$-mers' abundance-based distances among datasets of reads, commonly used in metagenomic analyses. In this work, we develop, analyze, and test, a sampling-based approach, called SAKEIMA, to approximate the frequent $k$-mers and their frequencies in a high-throughput sequencing dataset while providing rigorous guarantees on the quality of the approximation. SAKEIMA employs an advanced sampling scheme and we show how the characterization of the VC dimension, a core concept from statistical learning theory, of a properly defined set of functions leads to practical bounds on the sample size required for a rigorous approximation. Our experimental evaluation shows that SAKEIMA allows to rigorously approximate frequent $k$-mers by processing only a fraction of a dataset and that the frequencies estimated by SAKEIMA lead to accurate estimates of $k$-mer based distances between high-throughput sequencing datasets. Overall, SAKEIMA is an efficient and rigorous tool to estimate $k$-mers abundances providing significant speed-ups in the analysis of large sequencing datasets.


title: Self-Organization in Spontaneous Movements of Neonates generates Self-specifying Sensory Experiences

abstract: Movement experience and the coordination of perception and action are the basis of developing body awareness, emotion, motivation and cognition and the sense of self. The four limbs play a key role in the developing sense of body ownership, agency and peripersonal space. Neonatal limb movements were investigated respective processes of self-organization and developing body awareness. Application of symbolic dynamics to kinematic data of radial distance trajectories of the hands and feet determined 16 coordination patterns according to the distance of the effectors to the body center. From time series analysis of occurrence and recurrence of these configurations of 25 movement episodes of 14 infants both characteristics of self-organization as well as features assumed to provide neonates with experiences involved in the development of body awareness were found. With increasing age a shift from configurations with proximal to distal positions suggests a role of the proximal-distal dimension in movement development. We conclude that self-organization in spontaneous movements provides neonates with perceptual body- and self-specifying stimuli involved in developing body awareness and postulate the involvement of emotional and cognitive processes.


title: Nonlinear Approximation via Compositions

abstract: Given a function dictionary $\cal D$ and an approximation budget $N\in\mathbb{N}^+$, nonlinear approximation seeks the linear combination of the best $N$ terms $\{T_n\}_{1\le n\le N}\subseteq{\cal D}$ to approximate a given function $f$ with the minimum approximation error\[\varepsilon_{L,f}:=\min_{\{g_n\}\subseteq{\mathbb{R}},\{T_n\}\subseteq{\cal D}}\|f(x)-\sum_{n=1}^N g_n T_n(x)\|.\]Motivated by recent success of deep learning, we propose dictionaries with functions in a form of compositions, i.e.,\[T(x)=T^{(L)}\circ T^{(L-1)}\circ\cdots\circ T^{(1)}(x)\]for all $T\in\cal D$, and implement $T$ using ReLU feed-forward neural networks (FNNs) with $L$ hidden layers. We further quantify the improvement of the best $N$-term approximation rate in terms of $N$ when $L$ is increased from $1$ to $2$ or $3$ to show the power of compositions. In the case when $L>3$, our analysis shows that increasing $L$ cannot improve the approximation rate in terms of $N$. In particular, for any function $f$ on $[0,1]$, regardless of its smoothness and even the continuity, if $f$ can be approximated using a dictionary when $L=1$ with the best $N$-term approximation rate $\varepsilon_{L,f}={\cal O}(N^{-\eta})$, we show that dictionaries with $L=2$ can improve the best $N$-term approximation rate to $\varepsilon_{L,f}={\cal O}(N^{-2\eta})$. We also show that for H\"older continuous functions of order $\alpha$ on $[0,1]^d$, the application of a dictionary with $L=3$ in nonlinear approximation can achieve an essentially tight best $N$-term approximation rate $\varepsilon_{L,f}={\cal O}(N^{-2\alpha/d})$. Finally, we show that dictionaries consisting of wide FNNs with a few hidden layers are more attractive in terms of computational efficiency than dictionaries with narrow and very deep FNNs for approximating H\"older continuous functions if the number of computer cores is larger than $N$ in parallel computing.


title: Ground-state energy estimation of the water molecule on a trapped ion quantum computer

abstract: Quantum computing leverages the quantum resources of superposition and entanglement to efficiently solve computational problems considered intractable for classical computers. Examples include calculating molecular and nuclear structure, simulating strongly-interacting electron systems, and modeling aspects of material function. While substantial theoretical advances have been made in mapping these problems to quantum algorithms, there remains a large gap between the resource requirements for solving such problems and the capabilities of currently available quantum hardware. Bridging this gap will require a co-design approach, where the expression of algorithms is developed in conjunction with the hardware itself to optimize execution. Here, we describe a scalable co-design framework for solving chemistry problems on a trapped ion quantum computer, and apply it to compute the ground-state energy of the water molecule. The robust operation of the trapped ion quantum computer yields energy estimates with errors approaching the chemical accuracy, which is the target threshold necessary for predicting the rates of chemical reaction dynamics.


title: Near Optimal Algorithms for Hard Submodular Programs with Discounted Cooperative Costs

abstract: In this paper, we investigate a class of submodular problems which in general are very hard. These include minimizing a submodular cost function under combinatorial constraints, which include cuts, matchings, paths, etc., optimizing a submodular function under submodular cover and submodular knapsack constraints, and minimizing a ratio of submodular functions. All these problems appear in several real world problems but have hardness factors of $\Omega(\sqrt{n})$ for general submodular cost functions. We show how we can achieve constant approximation factors when we restrict the cost functions to low rank sums of concave over modular functions. A wide variety of machine learning applications are very naturally modeled via this subclass of submodular functions. Our work therefore provides a tighter connection between theory and practice by enabling theoretically satisfying guarantees for a rich class of expressible, natural, and useful submodular cost models. We empirically demonstrate the utility of our models on real world problems of cooperative image matching and sensor placement with cooperative costs.


title: Online Learning with Continuous Ranked Probability Score

abstract: Probabilistic forecasts in the form of probability distributions over future events have become popular in several fields of statistical science. The dissimilarity between a probability forecast and an outcome is measured by a loss function (scoring rule). Popular example of scoring rule for continuous outcomes is the continuous ranked probability score (CRPS). We consider the case where several competing methods produce online predictions in the form of probability distribution functions. In this paper, the problem of combining probabilistic forecasts is considered in the prediction with expert advice framework. We show that CRPS is a mixable loss function and then the time independent upper bound for the regret of the Vovk's aggregating algorithm using CRPS as a loss function can be obtained. We present the results of numerical experiments illustrating the proposed methods.


title: Two dimensional ordering and collective magnetic excitations in the dilute ferromagnetic topological insulator (Bi$_{0.95}$Mn$_{0.05}$)$_{2}$Te$_{3}$

abstract: Employing elastic and inelastic neutron scattering (INS) techniques, we report on detailed microscopic properties of the ferromagnetism in he magnetic topological insulator (Bi$_{0.95}$Mn$_{0.05}$)$_{2}$Te$_{3}$. Neutron diffraction of polycrystalline samples show the ferromagnetic (FM) ordering is long-range within the basal plane, and mainly 2D in character with short-range correlations between layers below $T_{\mathrm{C}} \approx 13$ K. Despite the random distribution of the dliute Mn atoms, we find that the 2D-like magnetic peaks are commensurate with the chemical structure, and the absence of (00L) magnetic peaks denote that the Mn$^{2+}$ magnetic moments are normal to the basal planes. Surprisingly, we observed collective magnetic excitations, in this dilute magnetic system. Despite the dilute nature, the excitations are typical of quasi-2D FM systems, albeit are severely broadened at short wavelengths, likely due to the random spatial distribution of Mn atoms in the Bi planes. Detailed analysis of the INS provide energy scales of the exchange couplings and the single ion anisotropy.


title: Six new supermassive black hole mass determinations from adaptive-optics assisted SINFONI observations

abstract: Different massive black hole mass - host galaxy scaling relations suggest that the growth of massive black holes is entangled with the evolution of their host galaxies. The number of measured black hole masses is still limited, and additional measurements are necessary to understand the underlying physics of this apparent co-evolution. We add six new black hole mass (MBH) measurements of nearby fast rotating early-type galaxies to the known black hole mass sample, namely NGC 584, NGC 2784, NGC 3640, NGC 4570, NGC 4281 and NGC 7049. Our target galaxies have effective velocity dispersions ({\sigma}e) between 170 and 245 km s^(-1), and thus this work provides additional insight into the black hole properties of intermediate-mass early-type galaxies. We combine high-resolution adaptive-optics SINFONI data with large-scale MUSE, VIMOS and SAURON data from ATLAS3D to derive two-dimensional stellar kinematics maps. We then build both Jeans Anisotropic Models and axisymmetric Schwarzschild models to measure the central black hole masses. Our Schwarzschild models provide black hole masses which are consistent with recent MBH-{\sigma}e scaling relations. NGC 3640 has a velocity dispersion dip and NGC 7049 a constant velocity dispersion in the center, but we can clearly constrain their lower black hole mass limit. We conclude our analysis with a test on NGC 4570 taking into account a variable mass-to-light ratio (M/L) when constructing dynamical models. When considering M/L variations linked mostly to radial changes in the stellar metallicity, we find that the dynamically determined black hole mass from NGC 4570 decreases by 30%. Further investigations are needed in the future to account for the impact of radial M/L gradients on dynamical modeling.


title: A Memoization Framework for Scaling Submodular Optimization to Large Scale Problems

abstract: We are motivated by large scale submodular optimization problems, where standard algorithms that treat the submodular functions in the \emph{value oracle model} do not scale. In this paper, we present a model called the \emph{precomputational complexity model}, along with a unifying memoization based framework, which looks at the specific form of the given submodular function. A key ingredient in this framework is the notion of a \emph{precomputed statistic}, which is maintained in the course of the algorithms. We show that we can easily integrate this idea into a large class of submodular optimization problems including constrained and unconstrained submodular maximization, minimization, difference of submodular optimization, optimization with submodular constraints and several other related optimization problems. Moreover, memoization can be integrated in both discrete and continuous relaxation flavors of algorithms for these problems. We demonstrate this idea for several commonly occurring submodular functions, and show how the precomputational model provides significant speedups compared to the value oracle model. Finally, we empirically demonstrate this for large scale machine learning problems of data subset selection and summarization.


title: The optimal packing of eight points in the real projective plane

abstract: How can we arrange $n$ lines through the origin in three-dimensional Euclidean space in a way that maximizes the minimum interior angle between pairs of lines? Conway, Hardin and Sloane (1996) produced line packings for $n \leq 55$ that they conjectured to be within numerical precision of optimal in this sense, but until now only the cases $n \leq 7$ have been solved. In this paper, we resolve the case $n = 8$. Drawing inspiration from recent work on the Tammes problem, we enumerate contact graph candidates for an optimal configuration and eliminate those that violate various combinatorial and geometric necessary conditions. The contact graph of the putatively optimal numerical packing of Conway, Hardin and Sloane is the only graph that survives, and we recover from this graph an exact expression for the minimum distance of eight optimally packed points in the real projective plane.


title: Unmasking Clever Hans Predictors and Assessing What Machines Really Learn

abstract: Current learning machines have successfully solved hard application problems, reaching high accuracy and displaying seemingly "intelligent" behavior. Here we apply recent techniques for explaining decisions of state-of-the-art learning machines and analyze various tasks from computer vision and arcade games. This showcases a spectrum of problem-solving behaviors ranging from naive and short-sighted, to well-informed and strategic. We observe that standard performance evaluation metrics can be oblivious to distinguishing these diverse problem solving behaviors. Furthermore, we propose our semi-automated Spectral Relevance Analysis that provides a practically effective way of characterizing and validating the behavior of nonlinear learning machines. This helps to assess whether a learned model indeed delivers reliably for the problem that it was conceived for. Furthermore, our work intends to add a voice of caution to the ongoing excitement about machine intelligence and pledges to evaluate and judge some of these recent successes in a more nuanced manner.


title: An extension of the Bourgain-Sarnak-Ziegler theorem with modular applications

abstract: We prove an extension of the Bourgain-Sarnak-Ziegler theorem and then apply it to bound certain polynomial exponential sums with modular coefficients.


title: SCOPE: SCUBA-2 Continuum Observations of Pre-protostellar Evolution - Survey Description and Compact Source Catalogue

abstract: We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.


title: Mach-Zehnder interferometer with quantum beamsplitters

abstract: A one-dimensional waveguide enables a single two-level emitter to route the propagation of a single photon, as to provide a quantum mirror or a quantum beamsplitter. Here we present a fully-quantum Mach-Zehnder interferometer (QMZ) for single-photon pulses comprised of two quantum beamsplitters. We theoretically show how nonlinearities of the QMZ due to photon-emitter detunings and to the spectral linewidth of the pulse contribute to the versatility of the device with respect to the classical-beamsplitters scenario. We employ a quantum dynamics framework to obtain analytical expressions for the photodetection probabilities and prove, in the monochromatic regime, the equivalence with a transfer-matrix approach.


title: Obstacle-aware Adaptive Informative Path Planning for UAV-based Target Search

abstract: Target search with unmanned aerial vehicles (UAVs) is relevant problem to many scenarios, e.g., search and rescue (SaR). However, a key challenge is planning paths for maximal search efficiency given flight time constraints. To address this, we propose the Obstacle-aware Adaptive Informative Path Planning (OA-IPP) algorithm for target search in cluttered environments using UAVs. Our approach leverages a layered planning strategy using a Gaussian Process (GP)-based model of target occupancy to generate informative paths in continuous 3D space. Within this framework, we introduce an adaptive replanning scheme which allows us to trade off between information gain, field coverage, sensor performance, and collision avoidance for efficient target detection. Extensive simulations show that our OA-IPP method performs better than state-of-the-art planners, and we demonstrate its application in a realistic urban SaR scenario.


title: The effect of dust composition and shape on radiation-pressure forces and blowout sizes of particles in debris disks

abstract: The light scattered from dust grains in debris disks is typically modeled as compact spheres using Lorenz-Mie theory or as porous spheres by incorporating an effective medium theory. In this work we examine the effect of incorporating a more realistic particle morphology on estimated radiation-pressure blowout sizes. To calculate the scattering and absorption cross sections of irregularly shaped dust grains, we use the discrete dipole approximation. These cross sections are necessary to calculate the $\beta$-ratio, which determines whether dust grains can remain gravitationally bound to their star. We calculate blowout sizes for a range of stellar spectral types corresponding with stars known to host debris disks. As with compact spheres, more luminous stars blow out larger irregularly shaped dust grains. We also find that dust grain composition influences blowout size such that absorptive grains are more readily removed from the disk. Moreover, the difference between blowout sizes calculated assuming spherical particles versus particle morphologies more representative of real dust particles is compositionally dependent as well, with blowout size estimates diverging further for transparent grains. We find that the blowout sizes calculated have a strong dependence on the particle model used, with differences in the blowout size calculated being as large as an order of magnitude for particles of similar porosities.


title: Convergence in uncertain linear systems

abstract: State convergence is essential in several scientific areas, e.g. multi-agent consensus/disagreement, distributed optimization, monotone game theory, multi-agent learning over time-varying networks. This paper is the first on state convergence in both continuous- and discrete-time linear systems affected by polytopic uncertainty. First, we characterize state convergence in linear time invariant systems via equivalent necessary and sufficient conditions. In the presence of uncertainty, we complement the canonical definition of (weak) convergence with a stronger notion of convergence, which requires the existence of a common kernel among the generator matrices of the difference/differential inclusion (strong convergence). We investigate under which conditions the two definitions are equivalent. Then, we characterize weak and strong convergence by means of Lyapunov and LaSalle arguments, (linear) matrix inequalities and separability of the eigenvalues of the generator matrices. We also show that, unlike asymptotic stability, state convergence lacks of duality.


title: On some functional generalizations of the regularity of topological spaces

abstract: We introduce and study some generalizations of regular spaces, which were motivated by studying continuity properties of functions between (regular) topological spaces. In particular, we prove that a first-countable Hausdorff topological space is regular if and only if it does not contain a topological copy of the Gutik hedgehog.


title: Attention is not Explanation

abstract: Attention mechanisms have seen wide adoption in neural NLP models. In addition to improving predictive performance, these are often touted as affording transparency: models equipped with attention provide a distribution over attended-to input units, and this is often presented (at least implicitly) as communicating the relative importance of inputs. However, it is unclear what relationship exists between attention weights and model outputs. In this work, we perform extensive experiments across a variety of NLP tasks that aim to assess the degree to which attention weights provide meaningful `explanations' for predictions. We find that they largely do not. For example, learned attention weights are frequently uncorrelated with gradient-based measures of feature importance, and one can identify very different attention distributions that nonetheless yield equivalent predictions. Our findings show that standard attention modules do not provide meaningful explanations and should not be treated as though they do. Code for all experiments is available at https://github.com/successar/AttentionExplanation.


title: Numerical Approximation of Young Measure Solutions to Parabolic Systems of Forward-Backward Type

abstract: This paper is concerned with the proof of existence and numerical approximation of large-data global-in-time Young measure solutions to initial-boundary-value problems for multidimensional nonlinear parabolic systems of forward-backward type of the form $\partial_t u - \mbox{div}(a(Du)) + Bu = F$, where $B \in \mathbb{R}^{m \times m}$, $Bv \cdot v \geq 0$ for all $v \in \mathbb{R}^m$, $F$ is an $m$-component vector-function defined on a bounded open Lipschitz domain $\Omega \subset \mathbb{R}^n$, and $a$ is a locally Lipschitz mapping of the form $a(A)=K(A)A$, where $K\,:\, \mathbb{R}^{m \times n} \rightarrow \mathbb{R}$. The function $a$ may have a nonstandard growth rate, in the sense that it is permitted to have unequal lower and upper growth rates. Furthermore, $a$ is not assumed to be monotone, nor is it assumed to be the gradient of a potential. Problems of this type arise in mathematical models of the atmospheric boundary layer and fall beyond the scope of monotone operator theory. We develop a numerical algorithm for the approximate solution of problems in this class, and we prove the convergence of the algorithm to a Young measure solution of the system under consideration.


title: On the Mortality Problem: from multiplicative matrix equations to linear recurrence sequences and beyond

abstract: We consider the following variant of the Mortality Problem: given $k\times k$ matrices $A_1, A_2, \dots,A_{t}$, does there exist nonnegative integers $m_1, m_2, \dots,m_t$ such that the product $A_1^{m_1} A_2^{m_2} \cdots A_{t}^{m_{t}}$ is equal to the zero matrix? It is known that this problem is decidable when $t \leq 2$ for matrices over algebraic numbers but becomes undecidable for sufficiently large $t$ and $k$ even for integral matrices. In this paper, we prove the first decidability results for $t>2$. We show as one of our central results that for $t=3$ this problem in any dimension is Turing equivalent to the well-known Skolem problem for linear recurrence sequences. Our proof relies on the Primary Decomposition Theorem for matrices that was not used to show decidability results in matrix semigroups before. As a corollary we obtain that the above problem is decidable for $t=3$ and $k \leq 3$ for matrices over algebraic numbers and for $t=3$ and $k=4$ for matrices over real algebraic numbers. Another consequence is that the set of triples $(m_1,m_2,m_3)$ for which the equation $A_1^{m_1} A_2^{m_2} A_3^{m_3}$ equals the zero matrix is equal to a finite union of direct products of semilinear sets. For $t=4$ we show that the solution set can be non-semilinear, and thus it seems unlikely that there is a direct connection to the Skolem problem. However we prove that the problem is still decidable for upper-triangular $2 \times 2$ rational matrices by employing powerful tools from transcendence theory such as Baker's theorem and S-unit equations.


title: Variational Inference to Measure Model Uncertainty in Deep Neural Networks

abstract: We present a novel approach for training deep neural networks in a Bayesian way. Classical, i.e. non-Bayesian, deep learning has two major drawbacks both originating from the fact that network parameters are considered to be deterministic. First, model uncertainty cannot be measured thus limiting the use of deep learning in many fields of application and second, training of deep neural networks is often hampered by overfitting. The proposed approach uses variational inference to approximate the intractable a posteriori distribution on basis of a normal prior. The variational density is designed in such a way that the a posteriori uncertainty of the network parameters is represented per network layer and depending on the estimated parameter expectation values. This way, only a few additional parameters need to be optimized compared to a non-Bayesian network. We apply this Bayesian approach to train and test the LeNet architecture on the MNIST dataset. Compared to classical deep learning, the test error is reduced by 15%. In addition, the trained model contains information about the parameter uncertainty in each layer. We show that this information can be used to calculate credible intervals for the prediction and to optimize the network architecture for a given training data set.


title: High Flux Passive Imaging with Single-Photon Sensors

abstract: Single-photon avalanche diodes (SPADs) are an emerging technology with a unique capability of capturing individual photons with high timing precision. SPADs are being used in several active imaging systems (e.g., fluorescence lifetime microscopy and LiDAR), albeit mostly limited to low photon flux settings. We propose passive free-running SPAD (PF-SPAD) imaging, an imaging modality that uses SPADs for capturing 2D intensity images with unprecedented dynamic range under ambient lighting, without any active light source. Our key observation is that the precise inter-photon timing measured by a SPAD can be used for estimating scene brightness under ambient lighting conditions, even for very bright scenes. We develop a theoretical model for PF-SPAD imaging, and derive a scene brightness estimator based on the average time of darkness between successive photons detected by a PF-SPAD pixel. Our key insight is that due to the stochastic nature of photon arrivals, this estimator does not suffer from a hard saturation limit. Coupled with high sensitivity at low flux, this enables a PF-SPAD pixel to measure a wide range of scene brightness, from very low to very high, thereby achieving extreme dynamic range. We demonstrate an improvement of over 2 orders of magnitude over conventional sensors by imaging scenes spanning a dynamic range of 1,000,000:1.


title: EvolveGCN: Evolving Graph Convolutional Networks for Dynamic Graphs

abstract: Graph representation learning resurges as a trending research subject owing to the widespread use of deep learning for Euclidean data, which inspire various creative designs of neural networks in the non-Euclidean domain, particularly graphs. With the success of these graph neural networks (GNN) in the static setting, we approach further practical scenarios where the graph dynamically evolves. Existing approaches typically resort to node embeddings and use a recurrent neural network (RNN, broadly speaking) to regulate the embeddings and learn the temporal dynamics. These methods require the knowledge of a node in the full time span (including both training and testing) and are less applicable to the frequent change of the node set. In some extreme scenarios, the node sets at different time steps may completely differ. To resolve this challenge, we propose EvolveGCN, which adapts the graph convolutional network (GCN) model along the temporal dimension without resorting to node embeddings. The proposed approach captures the dynamism of the graph sequence through using an RNN to evolve the GCN parameters. Two architectures are considered for the parameter evolution. We evaluate the proposed approach on tasks including link prediction, edge classification, and node classification. The experimental results indicate a generally higher performance of EvolveGCN compared with related approaches. The code is available at \url{https://github.com/IBM/EvolveGCN}.


title: A Graph Computation based Sequential Power Flow Calculation for Large-Scale ACDC Systems

abstract: This paper proposes a graph computation based sequential power flow calculation method for Line Commutated Converter (LCC) based large-scale AC/DC systems to achieve a high computing performance. Based on the graph theory, the complex AC/DC system is first converted to a graph model and stored in a graph database. Then, the hybrid system is divided into several isolated areas with graph partition algorithm by decoupling AC and DC networks. Thus, the power flow analysis can be executed in parallel for each independent area with the new selected slack buses. Furthermore, for each area, the node-based parallel computing (NPC) and hierarchical parallel computing (HPC) used in graph computation are employed to speed up fast decoupled power flow (FDPF). Comprehensive case studies on the IEEE 300-bus, polished South Carolina 12,000-bus system and a China 11,119-bus system are performed to demonstrate the accuracy and efficiency of the proposed method


title: On the Idiosyncrasies of the Mandarin Chinese Classifier System

abstract: While idiosyncrasies of the Chinese classifier system have been a richly studied topic among linguists (Adams and Conklin, 1973; Erbaugh, 1986; Lakoff, 1986), not much work has been done to quantify them with statistical methods. In this paper, we introduce an information-theoretic approach to measuring idiosyncrasy; we examine how much the uncertainty in Mandarin Chinese classifiers can be reduced by knowing semantic information about the nouns that the classifiers modify. Using the empirical distribution of classifiers from the parsed Chinese Gigaword corpus (Graff et al., 2005), we compute the mutual information (in bits) between the distribution over classifiers and distributions over other linguistic quantities. We investigate whether semantic classes of nouns and adjectives differ in how much they reduce uncertainty in classifier choice, and find that it is not fully idiosyncratic; while there are no obvious trends for the majority of semantic classes, shape nouns reduce uncertainty in classifier choice the most.


title: Learning to See the Wood for the Trees: Deep Laser Localization in Urban and Natural Environments on a CPU

abstract: Localization in challenging, natural environments such as forests or woodlands is an important capability for many applications from guiding a robot navigating along a forest trail to monitoring vegetation growth with handheld sensors. In this work we explore laser-based localization in both urban and natural environments, which is suitable for online applications. We propose a deep learning approach capable of learning meaningful descriptors directly from 3D point clouds by comparing triplets (anchor, positive and negative examples). The approach learns a feature space representation for a set of segmented point clouds that are matched between a current and previous observations. Our learning method is tailored towards loop closure detection resulting in a small model which can be deployed using only a CPU. The proposed learning method would allow the full pipeline to run on robots with limited computational payload such as drones, quadrupeds or UGVs.


title: Multivariate analysis of covariance when standard assumptions are violated

abstract: In applied research, it is often sensible to account for one or several covariates when testing for differences between multivariate means of several groups. However, the "classical" parametric multivariate analysis of covariance (MANCOVA) tests (e.g., Wilks' Lambda) are based on quite restrictive assumptions (homoscedasticity and normality of the errors), which might be difficult to justify in small sample size settings. Furthermore, existing potential remedies (e.g., heteroskedasticity-robust approaches) become inappropriate in cases where the covariance matrices are singular. Nevertheless, such scenarios are frequently encountered in the life sciences and other fields, when for example, in the context of standardized assessments, a summary performance measure as well as its corresponding subscales are analyzed. In the present manuscript, we consider a general MANCOVA model, allowing for potentially heteroskedastic and even singular covariance matrices as well as non-normal errors. We combine heteroskedasticity-consistent covariance matrix estimation methods with our proposed modified MANCOVA ANOVA-type statistic (MANCATS) and apply two different bootstrap approaches. We provide the proofs of the asymptotic validity of the respective testing procedures as well as the results from an extensive simulation study, which indicate that especially the parametric bootstrap version of the MANCATS outperforms its competitors in most scenarios, both in terms of type I error rates and power. These considerations are further illustrated and substantiated by examining real-life data from standardized achievement tests.


title: Transmission of Video Image in Underwater Acoustic Communication

abstract: Acoustic communication is currently considered as the best way to transmit information over long distances under water, since acoustic waves have lower attenuation in water than the other information transmission media. However, current information transmission rates of the underwater acoustic communication is usually limited by low-frequency bandwidth due to the large attenuation which characterized in high-frequency sound waves in water. Spiral acoustic beams with helical wavefront dislocations carrying orbital angular momentum have been applied to acoustic levitation and capture. Recently, harnessing the orthogonality and the infinite dimension of Hilbert space in underwater acoustic communication has been considered promising. Here, we construct a real-time underwater acoustic communication system to transmit video image. The system multiplexes eight orbital angular momentum topology charges ranging from +1 to +8 and encoded with on-off keying modulation format, achieving a spectral efficiency of 8 (bit/s)/Hz. This demonstration suggests that multiplexed information-carrying orbital angular momentum acoustic beams presents tremendous potential for increasing the capacity of underwater communication systems.


title: RotatE: Knowledge Graph Embedding by Relational Rotation in Complex Space

abstract: We study the problem of learning representations of entities and relations in knowledge graphs for predicting missing links. The success of such a task heavily relies on the ability of modeling and inferring the patterns of (or between) the relations. In this paper, we present a new approach for knowledge graph embedding called RotatE, which is able to model and infer various relation patterns including: symmetry/antisymmetry, inversion, and composition. Specifically, the RotatE model defines each relation as a rotation from the source entity to the target entity in the complex vector space. In addition, we propose a novel self-adversarial negative sampling technique for efficiently and effectively training the RotatE model. Experimental results on multiple benchmark knowledge graphs show that the proposed RotatE model is not only scalable, but also able to infer and model various relation patterns and significantly outperform existing state-of-the-art models for link prediction.


title: Competition with Three-Tier Spectrum Access and Spectrum Monitoring

abstract: The Citizens Broadband Radio Service (CBRS) recently adopted in the U.S. enables two tiers of commercial users to share spectrum with a third tier of incumbent users. This sharing can be further assisted by Environmental Sensing Capability operators (ESCs), that monitor the spectrum occupancy to determine when use of the spectrum will not harm incumbents. Two key aspects of this framework that impact how firms may compete are the differences in information provided by different ESCs and the different tiers in which a firm may access the spectrum. We develop a game theoretic model that captures both of these features and analyze it to gain insight into their impact. Specifically, we consider a priority access (PA) tier firm has access to the both licensed band and unlicensed band, and a general authorized access (GAA) tier firm has access only to the unlicensed band. The PA tier and GAA tier firms compete for users. Our analysis reveals that the amount of unlicensed and licensed bandwidth in the CBRS must be chosen judiciously in order to maximize the social welfare. We also show that a limited amount of unlicensed access of PA tier firm is beneficial to the user's surplus as well as to the social welfare.


title: Chemical Abundances of Main-Sequence, Turnoff, Subgiant and red giant Stars from APOGEE spectra II: Atomic Diffusion in M67 Stars

abstract: Chemical abundances for 15 elements (C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are presented for 83 stellar members of the 4 Gyr old solar-metallicity open cluster M67. The sample contains stars spanning a wide range of evolutionary phases, from G dwarfs to red clump stars. The abundances were derived from near-IR ($\lambda$1.5 -- 1.7$\mu$m) high-resolution spectra ($R$ = 22,500) from the SDSS-IV/APOGEE survey. A 1-D LTE abundance analysis was carried out using the APOGEE synthetic spectral libraries, via chi-square minimization of the synthetic and observed spectra with the qASPCAP code. We found significant abundance differences ($\sim$0.05 -- 0.30 dex) between the M67 member stars as a function of the stellar mass (or position on the HR diagram), where the abundance patterns exhibit a general depletion (in [X/H]) in stars at the main-sequence turnoff. The amount of the depletion is different for different elements. We find that atomic diffusion models provide, in general, good agreement with the abundance trends for most chemical species, supporting recent studies indicating that measurable atomic diffusion operates in M67 stars.


title: Differentiable Scene Graphs

abstract: Reasoning about complex visual scenes involves perception of entities and their relations. Scene graphs provide a natural representation for reasoning tasks, by assigning labels to both entities (nodes) and relations (edges). Unfortunately, reasoning systems based on SGs are typically trained in a two-step procedure: First, training a model to predict SGs from images; Then, a separate model is created to reason based on predicted SGs. In many domains, it is preferable to train systems jointly in an end-to-end manner, but SGs are not commonly used as intermediate components in visual reasoning systems because being discrete and sparse, scene-graph representations are non-differentiable and difficult to optimize. Here we propose Differentiable Scene Graphs (DSGs), an image representation that is amenable to differentiable end-to-end optimization, and requires supervision only from the downstream tasks. DSGs provide a dense representation for all regions and pairs of regions, and do not spend modelling capacity on areas of the images that do not contain objects or relations of interest. We evaluate our model on the challenging task of identifying referring relationships (RR) in three benchmark datasets, Visual Genome, VRD and CLEVR. We describe a multi-task objective, and train in an end-to-end manner supervised by the downstream RR task. Using DSGs as an intermediate representation leads to new state-of-the-art performance.


title: Curves with more than one inner Galois point

abstract: Let $\mathcal{C}$ be an irreducible plane curve of $\text{PG}(2,\mathbb{K})$ where $\mathbb{K}$ is an algebraically closed field of characteristic $p\geq 0$. A point $Q\in \mathcal{C}$ is an inner Galois point for $\mathcal{C}$ if the projection $\pi_Q$ from $Q$ is Galois. Assume that $\mathcal{C}$ has two different inner Galois points $Q_1$ and $Q_2$, both simple. Let $G_1$ and $G_2$ be the respective Galois groups. Under the assumption that $G_i$ fixes $Q_i$, for $i=1,2$, we provide a complete classification of $G=\langle G_1,G_2 \rangle$ and we exhibit a curve for each such $G$. Our proof relies on deeper results from group theory.


title: Inflation in the R\'{e}nyi cosmology

abstract: The description of the accelerated phases for the current and primordial cosmos has been obtained in the framework of R\'{e}nyi cosmology. The results show that corrections to the Friedmann equations, originated from differences between R\'{e}nyi and Bekenstein entropies, can describe both the current accelerated and inflationary eras. The effects of these corrections can be appointed to a hypothetical fluid which is classically unstable. Although inflation can be obtained in R\'{e}nyi cosmology as the vacuum solution, the values found for the parameters of the hypothetical fluid which supports the current accelerated and primary inflationary cosmos differ from each other.


title: The local and global dynamics of a general cancer tumor growth model with multiphase structure

abstract: We present a phase-space analysis of a mathematical model of tumor growth with an immune responses. We consider mathematical analysis of the model equations with multipoint initial condition regarding to dissipativity, boundedness of solutions, invariance of non-negativity, local and global stability and the basins of attractions. We derive some features of behavior of one of three-dimensional tumor growth models with dynamics described in terms of densities of three cells populations: tumor cells, healthy host cells and effector immune cells. We found sufficient conditions, under which trajectories from the positive domain of feasible multipoint initial conditions tend to one of equilibrium points. Here, cases of the small tumor mass equilibrium points-the healthy equilibrium point, the "death" equilibrium point have been examined. Biological implications of our results are discussed.


title: A constructive solution to a problem of ranking tournaments

abstract: A tournament is an oriented complete graph. The problem of ranking tournaments was firstly investigated by P. Erd\H{o}s and J. W. Moon. By probabilistic methods, the existence of "unrankable" tournaments was proved. On the other hand, they also mentioned the problem of explicit constructions. However, there seems to be only a few of explicit constructions of such tournaments. In this note, we give a construction of many such tournaments by using skew Hadamard difference sets which have been investigated in combinatorial design theory.


title: Deep MR Fingerprinting with total-variation and low-rank subspace priors

abstract: Deep learning (DL) has recently emerged to address the heavy storage and computation requirements of the baseline dictionary-matching (DM) for Magnetic Resonance Fingerprinting (MRF) reconstruction. Fed with non-iterated back-projected images, the network is unable to fully resolve spatially-correlated corruptions caused from the undersampling artefacts. We propose an accelerated iterative reconstruction to minimize these artefacts before feeding into the network. This is done through a convex regularization that jointly promotes spatio-temporal regularities of the MRF time-series. Except for training, the rest of the parameter estimation pipeline is dictionary-free. We validate the proposed approach on synthetic and in-vivo datasets.


title: Superconductivity at 161 K in Thorium Hydride $ThH_{10}$: Synthesis and Properties

abstract: Here we report targeted high-pressure synthesis of two novel high-$T_C$ hydride superconductors, $P6_3/mmc$-$ThH_9$ and $Fm\bar{3}m$-$ThH_{10}$, with the experimental critical temperatures ($T_C$) of 146 K and 159-161 K and upper critical magnetic fields ($\mu$$H_C$) 38 and 45 Tesla at pressures 170-175 Gigapascals, respectively. Superconductivity was evidenced by the observation of zero resistance and a decrease of $T_C$ under external magnetic field up to 16 Tesla. This is one of the highest critical temperatures that has been achieved experimentally in any compounds, along with such materials as $LaH_{10}$, $H_3S$ and $HgBa_2Ca_xCu_2O_{6+z}$. Our experiments show that $fcc$-$ThH_{10}$ has stabilization pressure of 85 GPa, making this material unique among all known high-$T_C$ metal polyhydrides. Two recently predicted Th-H compounds, $I4/mmm$-$ThH_4$ (> 86 GPa) and $Cmc2_1$-$ThH_6$ (86-104 GPa), were also synthesized. Equations of state of obtained thorium polyhydrides were measured and found to perfectly agree with the theoretical calculations. New phases were examined theoretically and their electronic, phonon, and superconducting properties were calculated.


title: Transverse properties of parabolic subgroups of Garside groups

abstract: Let $G$ be a Garside group endowed with the generating set $\mathcal{S}$ of non-trivial simple elements, and let $H$ be a parabolic subgroup of $G$. We determine a transversal $T$ of $H$ in $G$ such that each $\theta \in T$ is of minimal length in its right-coset, $H \theta$, for the word length with respect to $\mathcal{S}$. We show that there exists a regular language $L$ on $\mathcal{S} \cup \mathcal{S}^{-1}$ and a bijection $\mathrm{ev} : L \to T$ satisfying $\mathrm{lg} (U) = \mathrm{lg}_\mathcal{S}( \mathrm{ev}(U))$ for all $U \in L$. From this we deduce that the coset growth series of $H$ in $G$ is rational. Finally, we show that $G$ has fellow projections on $H$ but does not have bounded projections on $H$.


title: Numerical approximations for the variable coefficient fractional diffusion equations with non-smooth data

abstract: In this article we study the numerical approximation of a variable coefficient fractional diffusion equation. Using a change of variable, the variable coefficient fractional diffusion equation is transformed into a constant coefficient fractional diffusion equation of the same order. The transformed equation retains the desirable stability property of being an elliptic equation. A spectral approximation scheme is proposed and analyzed for the transformed equation, with error estimates for the approximated solution derived. An approximation to the unknown of the variable coefficient fractional diffusion equation is then obtained by post processing the computed approximation to the transformed equation. Error estimates are also presented for the approximation to the unknown of the variable coefficient equation with both smooth and non-smooth diffusivity coefficient and right-hand side. Three numerical experiments are given whose convergence results are in strong agreement with the theoretically derived estimates.


title: Proof-of-principle experimental demonstration of twin-field type quantum key distribution

abstract: The twin-field (TF) quantum key distribution (QKD) protocol and its variants are highly attractive because they can beat the well-known rate-loss limit (i.e., the PLOB bound) for QKD protocols without quantum repeaters. In this paper, we perform a proof-of-principle experimental demonstration of TF-QKD based on the protocol proposed by Curty et al. which removes from the original TF-QKD scheme the need for post-selection on the matching of a global phase, and can deliver nearly an order of magnitude higher secret key rate. Furthermore, we overcome the major difficulty in the practical implementation of TF-QKD, namely, the need to stabilize the phase of the quantum state over kilometers of fiber. A Sagnac loop structure is utilized to ensure excellent phase stability between the different parties. Using decoy states, we demonstrate secret-key generation rates that beat the PLOB bound when the channel loss is above 40 dB.


title: Adaptive Robust Energy Management Strategy for Campus-Based Commercial Buildings Considering Comprehensive Comfort Levels

abstract: Neglecting consumers' comfort always leads to failure or slow-response to demand response request. In this paper, we propose several comprehensive comfort level models for various appliances in campus-based commercial buildings (CBs). The objective of the proposed system is to minimize O\&M costs of campus-based CBs and maximize various comfort levels simultaneously under the worst-case scenarios. Adaptive robust optimization (ARO) is leveraged to handle various uncertainties within the proposed system: (i) demand response signals sending from the distribution system operator (DSO); (ii) arrival state-of-charge (SoC) conditions of plug-in electric vehicles (PEVs); (iii) power outputs of renewable energy sources (RESs); and (iv) load demand of other appliances. Benders decomposition, such as column-and-constraint generation (C\&CG) algorithm, is used to solve the reformulated NP-hard min-max problem. Extensive simulation results demonstrate the effectiveness of the proposed optimal energy management strategy for campus-based CBs in both minimizing O\&M costs and maximizing comprehensive comfort levels.


title: Band-gap engineering in AB(O$_x$S$_{1-x}$)$_{3}$ perovskite oxysulfides: A route to strongly polar materials for photocatalytic water splitting

abstract: Polar heterogeneous photocatalysts were shown to lead to enhanced charge-carrier separation that results in superior activity for example for photocatalytic water splitting. Promising photocatalyst materials such as oxynitrides can be rendered polar by epitaxial strain, which however also increases their band gap, making them unsuitable for visible light absorption. This suggests a trade-off between small band gaps and polar distortions - both being crucial for the catalyst's efficiency. In this paper we investigate, using density functional theory calculations, the suitability of strained AB(O$_x$S$_{1-x}$)$_{3}$ perovskites for photocatalytic water splitting. These materials normally have band gaps too small for water splitting but inducing polar distortions via epitaxial strain can increase the band gap to the suitable range. We find perovskite BaZr$_y$Ti$_{1-y}$O$_2$S compounds to be highly promising for photocatalytic water splitting due to large polar distortions and suitable band gaps.


title: Tameness, powerful images, and large cardinals

abstract: We provide comprehensive, level-by-level characterizations of large cardinals, in the range from weakly compact to strongly compact, by closure properties of powerful images of accessible functors. In the process, we show that these properties are also equivalent to various forms of tameness for abstract elementary classes. This systematizes and extends results of [BU17], [BTR16], [Lie18], and [LR16].


title: Reliable Deep Grade Prediction with Uncertainty Estimation

abstract: Currently, college-going students are taking longer to graduate than their parental generations. Further, in the United States, the six-year graduation rate has been 59% for decades. Improving the educational quality by training better-prepared students who can successfully graduate in a timely manner is critical. Accurately predicting students' grades in future courses has attracted much attention as it can help identify at-risk students early so that personalized feedback can be provided to them on time by advisors. Prior research on students' grade prediction include shallow linear models; however, students' learning is a highly complex process that involves the accumulation of knowledge across a sequence of courses that can not be sufficiently modeled by these linear models. In addition to that, prior approaches focus on prediction accuracy without considering prediction uncertainty, which is essential for advising and decision making. In this work, we present two types of Bayesian deep learning models for grade prediction. The MLP ignores the temporal dynamics of students' knowledge evolution. Hence, we propose RNN for students' performance prediction. To evaluate the performance of the proposed models, we performed extensive experiments on data collected from a large public university. The experimental results show that the proposed models achieve better performance than prior state-of-the-art approaches. Besides more accurate results, Bayesian deep learning models estimate uncertainty associated with the predictions. We explore how uncertainty estimation can be applied towards developing a reliable educational early warning system. In addition to uncertainty, we also develop an approach to explain the prediction results, which is useful for advisors to provide personalized feedback to students.


title: Implicit Kernel Learning

abstract: Kernels are powerful and versatile tools in machine learning and statistics. Although the notion of universal kernels and characteristic kernels has been studied, kernel selection still greatly influences the empirical performance. While learning the kernel in a data driven way has been investigated, in this paper we explore learning the spectral distribution of kernel via implicit generative models parametrized by deep neural networks. We called our method Implicit Kernel Learning (IKL). The proposed framework is simple to train and inference is performed via sampling random Fourier features. We investigate two applications of the proposed IKL as examples, including generative adversarial networks with MMD (MMD GAN) and standard supervised learning. Empirically, MMD GAN with IKL outperforms vanilla predefined kernels on both image and text generation benchmarks; using IKL with Random Kitchen Sinks also leads to substantial improvement over existing state-of-the-art kernel learning algorithms on popular supervised learning benchmarks. Theory and conditions for using IKL in both applications are also studied as well as connections to previous state-of-the-art methods.


title: Sharp Constants of Approximation Theory. III. Certain Polynomial Inequalities of Different Metrics on Convex Sets

abstract: Let $V\subset\R^m$ be a centrally symmetric convex body and let $V^*\subset\R^m$ be its polar. We prove limit relations between the sharp constants in the multivariate Markov-Bernstein-Nikolskii type inequalities for algebraic polynomials on $V^*$ and the corresponding constants for entire functions of exponential type with the spectrum in $V$.


title: Improving students' understanding of rotating frames of reference using videos from different perspectives

abstract: The concepts of the Coriolis and the centrifugal force are essential in various scientific fields and they are standard components of introductory physics lectures. In this paper we explore how students understand and apply concepts of rotating frames of reference in the context of an exemplary lecture demonstration experiment. We found in a $Predict-Observe-Explain$-setting, that after predicting the outcome prior to the demonstration, only one out of five physics students correctly reported the observation of the trajectory of a sphere rolling over a rotating disc. Despite this low score, a detailed analysis of distractors revealed significant conceptual learning during the observation of the experiment. In this context, we identified three main misconceptions and learning difficulties. First, the centrifugal force is only required to describe the trajectory if the object is coupled to the rotating system. Second, inertial forces cause a reaction of an object on which they act. And third, students systematically mix-up the trajectories in the stationary and the rotating frame of reference. Furthermore, we captured students' eye movements during the $Predict$ task and found that physics students with low confidence ratings focused longer on relevant task areas than confident students despite having a comparable score. Consequently, this metric is a helpful tool for the identification of misconceptions using eye tracking. Overall, the results help to understand the complexity of concept learning from demonstration experiments and provide important implications for instructional design of introductions to rotating frames of reference.


title: Compatibility of linear-response theory with the Second Law of Thermodynamics and the emergence of negative entropy production rates

abstract: The reliability of physical theories depends on whether they agree with well established physical laws. In this work, we address the compatibility of the Hamiltonian formulation of linear-response theory with the Second Law of Thermodynamics. In order to do so, we verify three complementary aspects often understood as statements of the Second Law: 1. No dissipation for quasistatic process; 2. Dissipation for finite-time processes; 3. Positive entropy production rate. Our analysis focus on two classes of nonequilibrium isothermal processes: slowly-varying and finite-time but weak ones. For the former, we show that these aspects are easily verified. For the later, we present conditions for the achievement of the first two aspects. We also show that the third one is not always verified, presenting an example based on Brownian motion in which we observe negative values in the entropy production rate. In particular, we compare linear-response and exact results for this example.


title: Room Temperature Quantum Coherence vs. Electron Transfer in a Rhodanine Derivative Chromophore

abstract: Understanding electron transfer in organic molecules is of great interest in quantum materials for light harvesting, energy conversion, and integration of molecules into solar cells. This, however, poses the challenge of designing specific optimal molecular structure for which the processes of ultrafast quantum coherence and electron transport are not so well understood. In this work, we investigate subpicosecond time scale quantum dynamics and electron transfer in an efficient electron acceptor Rhodanine chromophoric complex. We consider an open quantum system approach to model the complex-solvent interaction and compute the crossover from weak to strong dissipation on the reduced system dynamics for both a polar (Methanol) and a non-polar solvent (Toluene). We show that the electron transfer rates are enhanced in the strong chromophore-solvent coupling regime, being the highest transfer rates those found at room temperature. Even though the computed dynamics are highly non-Markovian, and they may exhibit a quantum character up to hundreds of femtoseconds, we show that quantum coherence does not necessarily optimize the electron transfer in the chromophore.


title: Magnetohydrodynamic-Particle-in-Cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution

abstract: The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfv\'en waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($\lambda \ll \rm pc$), low cosmic ray number density ($n_{\rm CR}/n_{\rm ISM} \sim 10^{-9}$), and weak CR anisotropy ($\sim v_A/c$) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic-particle-in-cell (MHD-PIC) method, which has unique advantages to alleviate these issues, we conduct one-dimensional simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the $\delta f$ method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum, equally shared between left and right handed Alfv\'en modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires pitch angles of most CRs to efficiently cross $90^\circ$, and can be captured in simulations with relatively high wave amplitude and/or high spatial resolution. We attribute this crossing to non-linear wave-particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate our methodology will open up opportunities for future investigations that incorporate additional physics.


title: A Near-coplanar Stellar Flyby of the Planet Host Star HD 106906

abstract: We present an investigation into the kinematics of HD 106906 using the newly released Gaia DR2 catalog to search for close encounters with other members of the Scorpius-Centaurus (Sco-Cen) association. HD 106906 is an eccentric spectroscopic binary that hosts both a large asymmetric debris disk extending out to at least 500 au and a directly imaged planetary-mass companion at a projected separation of 738 au. The cause of the asymmetry in the debris disk and the unusually wide separation of the planet is not currently known. Using a combination of Gaia DR2 astrometry and ground-based radial velocities, we explore the hypothesis that a close encounter with another cluster member within the last 15 Myr is responsible for the present configuration of the system. Out of 461 stars analyzed, we identified two candidate perturbers that had a median closest approach (CA) distance within 1 pc of HD 106906: HIP 59716 at $D_{\rm CA}=0.65_{-0.40}^{+0.93}$ pc ($t_{\rm CA}=-3.49_{-1.76}^{+0.90}$ Myr) and HIP 59721 at $D_{\rm CA}=0.71_{-0.11}^{+0.18}$ pc ($t_{\rm CA}=-2.18_{-1.04}^{+0.54}$ Myr), with the two stars likely forming a wide physical binary. The trajectories of both stars relative to HD 106906 are almost coplanar with the inner disk ($\Delta\theta = 5.4\pm1.7$ deg and $4.2_{-1.1}^{+0.9}$ deg). These two stars are the best candidates of the currently known members of Sco-Cen for having a dynamically important close encounter with HD 106906, which may have stabilized the orbit of HD 106906 b in the scenario where the planet formed in the inner system and attained high eccentricity by interaction with the central binary.


title: On Ramsey numbers of hedgehogs

abstract: The hedgehog $H_t$ is a 3-uniform hypergraph on vertices $1,\dots,t+\binom{t}{2}$ such that, for any pair $(i,j)$ with $1\le i<j\le t$, there exists a unique vertex $k>t$ such that $\{i,j,k\}$ is an edge. Conlon, Fox, and R\"odl proved that the two-color Ramsey number of the hedgehog grows polynomially in the number of its vertices, while the four-color Ramsey number grows exponentially in the number of its vertices. They asked whether the two-color Ramsey number of the hedgehog $H_t$ is nearly linear in the number of its vertices. We answer this question affirmatively, proving that $r(H_t) = O(t^2\ln t)$.


title: ROMANet: Fine-Grained Reuse-Driven Off-Chip Memory Access Management and Data Organization for Deep Neural Network Accelerators

abstract: Enabling high energy efficiency is crucial for embedded implementations of deep learning. Several studies have shown that the DRAM-based off-chip memory accesses are one of the most energy-consuming operations in deep neural network (DNN) accelerators, and thereby limit the designs from achieving efficiency gains at the full potential. DRAM access energy varies depending upon the number of accesses required as well as the energy consumed per-access. Therefore, searching for a solution towards the minimum DRAM access energy is an important optimization problem. Towards this, we propose the ROMANet methodology that aims at reducing the number of memory accesses, by searching for the appropriate data partitioning and scheduling for each layer of a network using a design space exploration, based on the knowledge of the available on-chip memory and the data reuse factors. Moreover, ROMANet also targets decreasing the number of DRAM row buffer conflicts and misses, by exploiting the DRAM multi-bank burst feature to improve the energy-per-access. Besides providing the energy benefits, our proposed DRAM data mapping also results in an increased effective DRAM throughput, which is useful for latency-constraint scenarios. Our experimental results show that the ROMANet saves DRAM access energy by 12% for the AlexNet, by 36% for the VGG-16, and by 46% for the MobileNet, while also improving the DRAM throughput by 10%, as compared to the state-of-the-art.


title: Virtual Environments for Rehabilitation of Postural Control Dysfunction

abstract: We developed a novel virtual reality [VR] platform with 3-dimensional sounds to help improve sensory integration and visuomotor processing for postural control and fall prevention in individuals with balance problems related to sensory deficits, such as vestibular dysfunction (disease of the inner ear). The system has scenes that simulate scenario-based environments. We can adjust the intensity of the visual and audio stimuli in the virtual scenes by controlling the user interface (UI) settings. A VR headset (HTC Vive or Oculus Rift) delivers stereo display while providing real-time position and orientation of the participants' head. The 3D game-like scenes make participants feel immersed and gradually exposes them to situations that may induce dizziness, anxiety or imbalance in their daily-living.


title: Prediction of the disease controllability in a complex network using machine learning algorithms

abstract: The application of machine learning (ML) techniques span a vast spectrum ranging from speech, face and character recognition, medical diagnosis, anomaly detection in data to the general classification, prediction, and regression problems. In the present work, we solve the problem of predicting R_0 for disease spreading on complex networks using the regression-based state-of-art ML techniques. R_0 is a metric that determines whether the disease-free epidemic or an endemic state is asymptotically stable and hence indicates the controllability of the disease spread. We predict R_0 , based on training the ML models with structural properties of complex networks, irrespective of the network type. The prediction is possible because: (a) The structure of complex networks plays an essential role in the spreading processes on networks (b) The regression techniques such as Support Vector Regression and Artificial Neural Network Model can be very efficiently used for prediction problems, even for non-linear data. We obtained good accuracy in the prediction of R_0 for the simulated networks as well as real-world networks using these techniques. Moreover, the ML model training is a one-time investment cost in terms of training time and memory, and the trained model can be used for predicting R_0 on unseen/new examples of networks.


title: Population modelling of FRBs from intrinsic properties

abstract: We present a method to estimate the source properties of FRBs from observations by assuming a fixed DM contribution from a MW-like host galaxy, pulse temporal broadening models for turbulent plasma and a flat FRB energy spectrum. We then perform Monte Carlo simulations to constrain the properties of the FRB source, its host galaxy and scattering in the intervening plasma from the observational data of FRBs detected with Parkes. The typical scatter broadening of the intrinsic pulse is found to be considerably small $\lesssim 10^{-2}-1\ {\rm ms}$ from physical models, with the ISM contribution suppressed significantly relative to IGM. The intrinsic width for non-repeating FRBs is broadened by a factor $\sim 2-3$ on average primarily due to dispersive smearing. From the simulations, we find that the host galaxy DM contribution is likely to be comparable to the Galactic contribution and the FRB energy decreases significantly at high frequencies with a negative spectral index. The FRB spatial density is found to increase up to redshift $\sim 2.0$ and then drops significantly at larger distances. We obtain the energy distribution for FRB 121102 with repetition rate $\sim 0.1-0.3\ {\rm hr^{-1}}$ and exponential energy cutoff that is significantly smaller compared to typical FRB energies. We find that the probability of observing none of the other FRBs to be repeating at Parkes is $\sim 0.8-0.9$ with the current follow-up data insufficient to suggest more than one class of FRB progenitors.


title: TensorMap: Lidar-Based Topological Mapping and Localization via Tensor Decompositions

abstract: We propose a technique to develop (and localize in) topological maps from light detection and ranging (Lidar) data. Localizing an autonomous vehicle with respect to a reference map in real-time is crucial for its safe operation. Owing to the rich information provided by Lidar sensors, these are emerging as a promising choice for this task. However, since a Lidar outputs a large amount of data every fraction of a second, it is progressively harder to process the information in real-time. Consequently, current systems have migrated towards faster alternatives at the expense of accuracy. To overcome this inherent trade-off between latency and accuracy, we propose a technique to develop topological maps from Lidar data using the orthogonal Tucker3 tensor decomposition. Our experimental evaluations demonstrate that in addition to achieving a high compression ratio as compared to full data, the proposed technique, $\textit{TensorMap}$, also accurately detects the position of the vehicle in a graph-based representation of a map. We also analyze the robustness of the proposed technique to Gaussian and translational noise, thus initiating explorations into potential applications of tensor decompositions in Lidar data analysis.


title: Chirality from quantum walks without quantum coin

abstract: Quantum walks (QWs) describe the evolution of quantum systems on graphs. An intrinsic degree of freedom---called the coin and represented by a finite-dimensional Hilbert space---is associated to each node. Scalar quantum walks are QWs with a one-dimensional coin. We propose a general strategy allowing one to construct scalar QWs on a broad variety of graphs, which admit embedding in Eulidean spaces, thus having a direct geometric interpretation. After reviewing the technique that allows one to regroup cells of nodes into new nodes, transforming finite spatial blocks into internal degrees of freedom, we prove that no QW with a two-dimensional coin can be derived from an isotropic scalar QW in this way. Finally we show that the Weyl and Dirac QWs can be derived from scalar QWs in spaces of dimension up to three, via our construction.


title: Continual Prediction from EHR Data for Inpatient Acute Kidney Injury

abstract: Acute kidney injury (AKI) commonly occurs in hospitalized patients and can lead to serious medical complications. In order to optimally predict AKI before it develops at any time during a hospital stay, we present a novel framework in which AKI is continually predicted automatically from EHR data over the entire hospital stay instead of at only one particular time. The continual model predicts AKI every time a patients AKI-relevant variable changes in the EHR. Thus the model is not only independent of a particular time for making predictions, but it can also leverage the latest values of all the AKI-relevant patient variables for making predictions. Using data of 44,691 hospital stays of duration longer than 24 hours we evaluated our continual prediction model and compared it with the traditional one-time prediction models. Excluding hospitals stays in which AKI occurred within 24 hours from admission, the one-time prediction model predicting at 24 hours from admission obtained area under ROC curve (AUC) of 0.653 while the continual prediction model obtained AUC of 0.724. The one-time prediction model that predicts at 24 hours obviously cannot predict AKI incidences that occur within 24 hours of admission which when included in the evaluation reduced its AUC to 0.57. In comparison, the continual prediction model had AUC of 0.709. The continual prediction model also did better than all other one-time prediction models predicting at other fixed times. By being able to take into account the latest values of AKI-relevant patient variables and by not being limited to a particular time of prediction, the continual prediction model out-performed one-time prediction models in predicting AKI.


title: Report on the Physics at the HL-LHC and Perspectives for the HE-LHC

abstract: Manuscript of volume 2 of the Yellow Report on the "Workshop on the Physics of HL-LHC and perspectives at HE-LHC" containing the notes produced by the ATLAS and CMS collaborations as input to the Yellow Report volume 1.


title: Ion-liquid based super-capacitors with inner gate diode-like separators

abstract: We demonstrate that the capacitance of ionic-liquid filled supercapacitors is substantially increased by placing a diode-like structure on the separator membrane. We call the structured separator: gate, and demonstrate that the order of a p-n layout with respect to the auxiliary electrode affects the overall cell's capacitance. The smallest ESR and the largest capacitance values are noted when the p-side is facing the auxiliary electrode.