title: Constraints on Spin-Independent Nucleus Scattering with sub-GeV Weakly Interacting Massive Particle Dark Matter from the CDEX-1B Experiment at the China Jin-Ping Laboratory

abstract: We report results on the searches of weakly interacting massive particles (WIMPs) with sub-GeV masses ($m_{\chi}$) via WIMP-nucleus spin-independent scattering with Migdal effect incorporated. Analysis on time-integrated (TI) and annual modulation (AM) effects on CDEX-1B data are performed, with 737.1 kg$\cdot$day exposure and 160 eVee threshold for TI analysis, and 1107.5 kg$\cdot$day exposure and 250 eVee threshold for AM analysis. The sensitive windows in $m_{\chi}$ are expanded by an order of magnitude to lower DM masses with Migdal effect incorporated. New limits on $\sigma_{\chi N}^{\rm SI}$ at 90\% confidence level are derived as $2\times$10$^{-32}\sim7\times$10$^{-35}$ $\rm cm^2$ for TI analysis at $m_{\chi}\sim$ 50$-$180 MeV/$c^2$, and $3\times$10$^{-32}\sim9\times$10$^{-38}$ $\rm cm^2$ for AM analysis at $m_{\chi}\sim$75 MeV/$c^2-$3.0 GeV/$c^2$.


title: Compactification of Extensive Game Structures and Backward Dominance Procedure

abstract: We study the relationship between invariant transformations on extensive game structures and backward dominance procedure (BD), a generalization of the classical backward induction introduced in Perea (2014). We show that behavioral equivalence with unambiguous orderings of information sets, a critical property that guarantees BD's applicability, can be characterized by the classical Coalescing and a modified Interchange/Simultanizing in Battigalli et al. (2020). We also give conditions on transformations that improve BD's efficiency. In addition, we discuss the relationship between transformations and Bonanno (2014)'s generalized backward induction.


title: Analogue Black Hole Spectroscopy; or, how to listen to dumb holes

abstract: Spectroscopy is a fundamental tool in science which consists in studying the response of a system as a function of frequency. Among its many applications in Physics, Biology, Chemistry and other fields, the possibility of identifying objects and structures through their emission spectra is remarkable and incredibly useful. In this paper we apply the spectroscopy idea to a numerically simulated hydrodynamical flow, with the goal of developing a new, non-invasive flow measurement technique. Our focus lies on an irrotational draining vortex, which can be seen, under specific conditions, as the analogue of a rotating black hole (historically named a dumb hole). This paper is a development of a recent experiment that suggests that irrotational vortices and rotating black holes share a common relaxation process, known as the ringdown phase. We apply techniques borrowed from black hole physics to identify vortex flows from their characteristic spectrum emitted during this ringdown phase. We believe that this technique is a new facet of the fluid-gravity analogy and constitutes a promising way to investigate experimentally vortex flows in fluids and superfluids alike.


title: Web Test Dependency Detection

abstract: E2E web test suites are prone to test dependencies due to the heterogeneous multi-tiered nature of modern web apps, which makes it difficult for developers to create isolated program states for each test case. In this paper, we present the first approach for detecting and validating test dependencies present in E2E web test suites. Our approach employs string analysis to extract an approximated set of dependencies from the test code. It then filters potential false dependencies through natural language processing of test names. Finally, it validates all dependencies, and uses a novel recovery algorithm to ensure no true dependencies are missed in the final test dependency graph. Our approach is implemented in a tool called TEDD and evaluated on the test suites of six open-source web apps. Our results show that TEDD can correctly detect and validate test dependencies up to 72% faster than the baseline with the original test ordering in which the graph contains all possible dependencies. The test dependency graphs produced by TEDD enable test execution parallelization, with a speed-up factor of up to 7x.


title: Deep Learning Methods for Mean Field Control Problems with Delay

abstract: We consider a general class of mean field control problems described by stochastic delayed differential equations of McKean-Vlasov type. Two numerical algorithms are provided based on deep learning techniques, one is to directly parameterize the optimal control using neural networks, the other is based on numerically solving the McKean-Vlasov forward anticipated backward stochastic differential equation (MV-FABSDE) system. In addition, we establish a necessary and sufficient stochastic maximum principle for this class of mean field control problems with delay based on the differential calculus on function of measures, as well as existence and uniqueness results for the associated MV-FABSDE system.


title: Ground state of the mass-critical inhomogeneous nonlinear Schrodinger functional

abstract: We study the ground state problem of the nonlinear Schrodinger functional with a mass-critical inhomogeneous nonlinear term. We provide the optimal condition for the existence of ground states and show that in the critical focusing regime there is a universal blow-up profile given by the unique optimizer of a Gagliardo-Nirenberg interpolation inequality.


title: Information-Theoretic Considerations in Batch Reinforcement Learning

abstract: Value-function approximation methods that operate in batch mode have foundational importance to reinforcement learning (RL). Finite sample guarantees for these methods often crucially rely on two types of assumptions: (1) mild distribution shift, and (2) representation conditions that are stronger than realizability. However, the necessity ("why do we need them?") and the naturalness ("when do they hold?") of such assumptions have largely eluded the literature. In this paper, we revisit these assumptions and provide theoretical results towards answering the above questions, and make steps towards a deeper understanding of value-function approximation.


title: Measuring dark matter-neutrino relative velocity on cosmological scales

abstract: We present a new method to measure neutrino masses using the dark matter-neutrino relative velocity. The relative motion between dark matter and neutrinos results in a parity-odd bispectrum which can be measured from cross-correlation of different cosmic fields. This new method is not affected by most systematics which are parity even and not limited by the knowledge of optical depth to the cosmic microwave background. We estimate the detectability of the relative velocity effect and find that the minimal sum of neutrino masses could be detected at high significance with upcoming surveys.


title: Inverse source problems for degenerate time-fractional PDE

abstract: In this paper, we investigate two inverse source problems for degenerate time-fractional partial differential equation in rectangular domains. The first problem involves a space-degenerate partial differential equation and the second one involves a time-degenerate partial differential equation. Solutions to both problem are expressed in series expansions. For the first problem, we obtained solutions in the form of Fourier-Legendre series. Convergence and uniqueness of solutions have been discussed. Solutions to the second problem are expressed in the form of Fourier-Sine series and they involve a generalized Mittag- Leffler type function. Moreover, we have established a new estimate for this generalized Mittag-Leffler type function. The obtained results are illustrated by providing example solutions using certain given data at the initial and final time.


title: Finding beam loss locations in a linac with oscillating dipole correctors

abstract: The paper proposes a method of finding the beam loss locations in a linac. If the beam is scraped at an aperture limitation, moving its centroid with two dipole correctors located upstream and oscillating in sync produces a line at the corresponding frequency in spectra of current-sensitive devices downstream of the loss point. The phase of this signal contains information about the location of the beam loss. Similar lines appear also in the position signals of Beam Position Monitors (BPMs). The phases of the BPM position lines change monotonically (within each 2{\pi}) along the linac and can be used a reference system. The phase of the loss signal compared with this reference system pinpoints the beam loss location, assuming that longitudinal coordinates of the BPMs are known. If the correctors deflection amplitudes and the phase offset between their waveforms are chosen optimally and well calibrated, the same measurement provides values of the \b{eta}-function and the betatron phase advance at the BPM locations. Optics measurements of this type can be made parasitically, with negligible effect on the emittance, if a long measurement time is acceptable.


title: Matching for the Israeli "Mechinot" Gap-Year Programs: Handling Rich Diversity Requirements

abstract: We describe our experience with designing and running a matching market for the Israeli "Mechinot" gap-year programs. The main conceptual challenge in the design of this market was the rich set of diversity considerations, which necessitated the development of an appropriate preference-specification language along with corresponding choice-function semantics, which we also theoretically analyze. Our contribution extends the existing toolbox for two-sided matching with soft constraints. This market was run for the first time in January 2018 and matched 1,607 candidates (out of a total of 3,120 candidates) to 35 different programs, has been run twice more since, and has been adopted by the Joint Council of the "Mechinot" gap-year programs for the foreseeable future.


title: Quantum Generalized Linear Models

abstract: Generalized linear models (GLM) are link function based statistical models. Many supervised learning algorithms are extensions of GLMs and have link functions built into the algorithm to model different outcome distributions. There are two major drawbacks when using this approach in applications using real world datasets. One is that none of the link functions available in the popular packages is a good fit for the data. Second, it is computationally inefficient and impractical to test all the possible distributions to find the optimum one. In addition, many GLMs and their machine learning extensions struggle on problems of overdispersion in Tweedie distributions. In this paper we propose a quantum extension to GLM that overcomes these drawbacks. A quantum gate with non-Gaussian transformation can be used to continuously deform the outcome distribution from known results. In doing so, we eliminate the need for a link function. Further, by using an algorithm that superposes all possible distributions to collapse to fit a dataset, we optimize the model in a computationally efficient way. We provide an initial proof-of-concept by testing this approach on both a simulation of overdispersed data and then on a benchmark dataset, which is quite overdispersed, and achieved state of the art results. This is a game changer in several applied fields, such as part failure modeling, medical research, actuarial science, finance and many other fields where Tweedie regression and overdispersion are ubiquitous.


title: Mass and spin measurements for the neutron star 4U1608-52 through the relativistic precession model

abstract: We present a systematic analysis of the complete set of observations of the neutron star low-mass X-ray binary 4U1608-52 obtained by the Rossi X-ray Timing Explorer's Proportional Counter Array. We study the spectral and fast-time variability properties of the source in order to determine the mass and spin of the neutron star via the relativistic precession model, and find 24 observations containing usable sets of the necessary three quasi-periodic oscillations (triplets) with which to accomplish this task, along with a further 7 observations containing two of the three quasi-periodic oscillations each. We calculate the spin and mass of the source for each of the triplets, and find that they give physically realistic estimates clustering in the spin range $0.19 < a < 0.35$ and mass range $2.15 < M/\textrm{M}_{\odot} < 2.6$. Neutron stars present environments for studying matter under the most extreme conditions of pressure and density; as their equation of state is not yet known, accurate measurements of their mass and spin will eventually allow for the discrimination between various models. We discuss the implications of our findings in the context of equation of state predictions, physically allowed spin ranges, emission proximity to the innermost stable circular orbit and possible model inaccuracies.


title: Probing Cosmic Axions through Resonant Emission and Absorption in Atomic Systems with Superradiance

abstract: The $\mu$eV-mass axion is one of the most promising candidates for cold dark matter, and remains to be a well-motivated solution to the CP problem of Quantum Chromodynamics (QCD) via the Peccei-Quinn mechanism. In this paper, we propose a novel method to detect the dark-matter axions in our galaxy via the resonant emission $|{\rm e}\rangle \to |{\rm g}\rangle + \gamma + \gamma^{\prime}_{} + a$ (or absorption $a + |{\rm e}\rangle \to |{\rm g}\rangle + \gamma + \gamma^{\prime}_{}$) in an atomic system with superradiance, where $|{\rm e}\rangle$ and $|{\rm g}\rangle$ stand for the excited and ground energy levels of atoms, respectively. A similar process via $|{\rm e}\rangle \to |{\rm g}\rangle + \gamma + a$ (or $a + |{\rm e}\rangle \to |{\rm g}\rangle + \gamma$) is also put forward to probe the axion-electron coupling. For the nominal experimental setup assuming a background-free environment, most of the parameter space for typical QCD axion models can be covered with parahydrogen molecules or ytterbium atoms. However, the background in a realistic experimental setup remains to be a major issue that needs to be solved in future studies. Searching for better atomic or molecular candidates may be required for a bigger signal-to-noise ratio.


title: All Adapted Topologies are Equal

abstract: A number of researchers have introduced topological structures on the set of laws of stochastic processes. A unifying goal of these authors is to strengthen the usual weak topology in order to adequately capture the temporal structure of stochastic processes. Aldous defines an extended weak topology based on the weak convergence of prediction processes. In the economic literature, Hellwig introduced the information topology to study the stability of equilibrium problems. Bion-Nadal and Talay introduce a version of the Wasserstein distance between the laws of diffusion processes. Pflug and Pichler consider the nested distance (and the weak nested topology) to obtain continuity of stochastic multistage programming problems. These distances can be seen as a symmetrization of Lassalle's causal transport problem, but there are also further natural ways to derive a topology from causal transport. Our main result is that all of these seemingly independent approaches define the same topology in finite discrete time. Moreover we show that this 'weak adapted topology' is characterized as the coarsest topology that guarantees continuity of optimal stopping problems for continuous bounded reward functions.


title: Fast hashing with Strong Concentration Bounds

abstract: Previous work on tabulation hashing by Patrascu and Thorup from STOC'11 on simple tabulation and from SODA'13 on twisted tabulation offered Chernoff-style concentration bounds on hash based sums, e.g., the number of balls/keys hashing to a given bin, but under some quite severe restrictions on the expected values of these sums. The basic idea in tabulation hashing is to view a key as consisting of $c=O(1)$ characters, e.g., a 64-bit key as $c=8$ characters of 8-bits. The character domain $\Sigma$ should be small enough that character tables of size $|\Sigma|$ fit in fast cache. The schemes then use $O(1)$ tables of this size, so the space of tabulation hashing is $O(|\Sigma|)$. However, the concentration bounds by Patrascu and Thorup only apply if the expected sums are $\ll |\Sigma|$. To see the problem, consider the very simple case where we use tabulation hashing to throw $n$ balls into $m$ bins and want to analyse the number of balls in a given bin. With their concentration bounds, we are fine if $n=m$, for then the expected value is $1$. However, if $m=2$, as when tossing $n$ unbiased coins, the expected value $n/2$ is $\gg |\Sigma|$ for large data sets, e.g., data sets that do not fit in fast cache. To handle expectations that go beyond the limits of our small space, we need a much more advanced analysis of simple tabulation, plus a new tabulation technique that we call \emph{tabulation-permutation} hashing which is at most twice as slow as simple tabulation. No other hashing scheme of comparable speed offers similar Chernoff-style concentration bounds.


title: GAMA+KiDS: Alignment of galaxies in galaxy groups and its dependence on galaxy scale

abstract: Intrinsic galaxy alignments are a source of bias for weak lensing measurements as well as a tool for understanding galaxy formation and evolution. In this work, we measure the alignment of shapes of satellite galaxies, in galaxy groups, with respect to the brightest group galaxy (BGG), as well as alignments of the BGG shape with the satellite positions, using the highly complete Galaxy And Mass Assembly (GAMA) spectroscopic survey and deep imaging from the Kilo Degree Survey. We control systematic errors with dedicated image simulations and measure accurate shapes using the DEIMOS shape measurement method. We find a significant satellite radial alignment signal, which vanishes at large separations from the BGG. We do not identify any strong trends of the signal with galaxy absolute magnitude or group mass. The alignment signal is dominated by red satellites. We also find that the outer regions of galaxies are aligned more strongly than their inner regions, by varying the radial weight employed during the shape measurement process. This behaviour is evident for both red and blue satellites. BGGs are also found to be aligned with satellite positions, with this alignment being stronger when considering the innermost satellites, using red BGGs and the shape of the outer region of the BGG. Lastly, we measure the global intrinsic alignment signal in the GAMA sample for two different radial weight functions and find no significant difference.


title: Stochastic Gravitational-Wave Background from Axion-Monodromy Oscillons in String Theory During Preheating

abstract: Axion monodromy is a generic phenomenon in string compactifications. We investigate the production of gravitational-wave background from axion monodromy inflation during preheating. By performing lattice simulations using pseudo-spectral methods, we find that significant single-peak stochastic gravitational-wave backgrounds are generated during preheating for asymptotic linear $|\phi|$ and $\phi^{2/3}$ potentials in axion monodromy inflation and cuspy potentials. This may allow us to explore string compactifications through gravitational waves.


title: Gender Classification from Iris Texture Images Using a New Set of Binary Statistical Image Features

abstract: Soft biometric information such as gender can contribute to many applications like as identification and security. This paper explores the use of a Binary Statistical Features (BSIF) algorithm for classifying gender from iris texture images captured with NIR sensors. It uses the same pipeline for iris recognition systems consisting of iris segmentation, normalisation and then classification. Experiments show that applying BSIF is not straightforward since it can create artificial textures causing misclassification. In order to overcome this limitation, a new set of filters was trained from eye images and different sized filters with padding bands were tested on a subject-disjoint database. A Modified-BSIF (MBSIF) method was implemented. The latter achieved better gender classification results (94.6\% and 91.33\% for the left and right eye respectively). These results are competitive with the state of the art in gender classification. In an additional contribution, a novel gender labelled database was created and it will be available upon request.


title: Autonomous engines driven by active matter: Energetics and design principles

abstract: Because of its nonequilibrium character, active matter in a steady state can drive engines that autonomously deliver work against a constant mechanical force or torque. As a generic model for such an engine, we consider systems that contain one or several active components and a single passive one that is asymmetric in its geometrical shape or its interactions. Generally, one expects that such an asymmetry leads to a persistent, directed current in the passive component, which can be used for the extraction of work. We validate this expectation for a minimal model consisting of an active and a passive particle on a one-dimensional lattice. It leads us to identify thermodynamically consistent measures for the efficiency of the conversion of isotropic activity to directed work. For systems with continuous degrees of freedom, work cannot be extracted using a one-dimensional geometry under quite general conditions. In contrast, we put forward two-dimensional shapes of a movable passive obstacle that are best suited for the extraction of work, which we compare with analytical results for an idealised work-extraction mechanism. For a setting with many noninteracting active particles, we use a mean-field approach to calculate the power and the efficiency, which we validate by simulations. Surprisingly, this approach reveals that the interaction with the passive obstacle can mediate cooperativity between otherwise noninteracting active particles, which enhances the extracted power per active particle significantly.


title: The Relationship between Centaurs and Jupiter Family Comets with Implications for K-Pg-type Impacts

abstract: Centaurs--icy bodies orbiting beyond Jupiter and interior to Neptune--are believed to be dynamically related to Jupiter Family Comets (JFCs), which have aphelia near Jupiter's orbit and perihelia in the inner Solar System. Previous dynamical simulations have recreated the Centaur/JFC conversion, but the mechanism behind that process remains poorly described. We have performed a numerical simulation of Centaur analogues that recreates this process, generating a dataset detailing over 2.6 million close planet/planetesimal interactions. We explore scenarios stored within that database and, from those, describe the mechanism by which Centaur objects are converted into JFCs. Because many JFCs have perihelia in the terrestrial planet region, and since Centaurs are constantly resupplied from the Scattered Disk and other reservoirs, the JFCs are an ever-present impact threat.


title: Dynamics of Merging Galaxy Clusters from Simulated Analogs

abstract: Merging galaxy clusters may provide a unique window into the behavior of dark matter and the evolution of member galaxies. To interpret these natural collider experiments we must account for how much time has passed since pericenter passage (TSP), the maximum relative speed of the merging subclusters, merger phase (outbound after first pericenter or returning for second pericenter), and other dynamical parameters that are not directly observable. These quantities are often inferred from staged simulations or analytical timing arguments that include neither substructure, large-scale structure, nor a cosmologically motivated range of impact parameters. We include all these effects by extracting dynamical parameters from analog systems in a cosmological n-body simulation, and we present constraints for 11 observed systems. The TSP and viewing angles we derive are consistent with those of staged hydrodynamical simulations, but we find lower maximum speeds. Compared to the analytical MCMAC method we find lower TSP, and viewing angles that put the separation vector closer to the plane of the sky; we attribute this to the MCMAC assumption of zero pericenter distance. We discuss potential extensions to the basic analog method as well as complementarities between methods.


title: Astro2020 Science White Paper: Probing the Structure of Interstellar Dust from Micron to Kpc Scales with X-ray Imaging

abstract: The X-ray regime is a largely underused resource for constraining interstellar dust grain models and improving our understanding of the physical processes that dictate how grains evolve over their lifetimes. This is mostly due to current detectors' relatively low sensitivity and high background, limiting the targets to the brightest sources. The improved sensitivity of the next generation of X-ray detectors will allow studies of much fainter sources, at much higher angular resolution, expanding our sampled sightlines in both quality and quantity.

categories: astro-ph.HE astro-ph.EP astro-ph.GA astro-ph.SR


title: Developing a large scale population screening tool for the assessment of Parkinson's disease using telephone-quality voice

abstract: Recent studies have demonstrated that analysis of laboratory-quality voice recordings can be used to accurately differentiate people diagnosed with Parkinson's disease (PD) from healthy controls (HC). These findings could help facilitate the development of remote screening and monitoring tools for PD. In this study, we analyzed 2759 telephone-quality voice recordings from 1483 PD and 15321 recordings from 8300 HC participants. To account for variations in phonetic backgrounds, we acquired data from seven countries. We developed a statistical framework for analyzing voice, whereby we computed 307 dysphonia measures that quantify different properties of voice impairment, such as, breathiness, roughness, monopitch, hoarse voice quality, and exaggerated vocal tremor. We used feature selection algorithms to identify robust parsimonious feature subsets, which were used in combination with a Random Forests (RF) classifier to accurately distinguish PD from HC. The best 10-fold cross-validation performance was obtained using Gram-Schmidt Orthogonalization (GSO) and RF, leading to mean sensitivity of 64.90% (standard deviation, SD 2.90%) and mean specificity of 67.96% (SD 2.90%). This large-scale study is a step forward towards assessing the development of a reliable, cost-effective and practical clinical decision support tool for screening the population at large for PD using telephone-quality voice.


title: Gradient-free activation maximization for identifying effective stimuli

abstract: A fundamental question for understanding brain function is what types of stimuli drive neurons to fire. In visual neuroscience, this question has also been posted as characterizing the receptive field of a neuron. The search for effective stimuli has traditionally been based on a combination of insights from previous studies, intuition, and luck. Recently, the same question has emerged in the study of units in convolutional neural networks (ConvNets), and together with this question a family of solutions were developed that are generally referred to as "feature visualization by activation maximization." We sought to bring in tools and techniques developed for studying ConvNets to the study of biological neural networks. However, one key difference that impedes direct translation of tools is that gradients can be obtained from ConvNets using backpropagation, but such gradients are not available from the brain. To circumvent this problem, we developed a method for gradient-free activation maximization by combining a generative neural network with a genetic algorithm. We termed this method XDream (EXtending DeepDream with real-time evolution for activation maximization), and we have shown that this method can reliably create strong stimuli for neurons in the macaque visual cortex (Ponce et al., 2019). In this paper, we describe extensive experiments characterizing the XDream method by using ConvNet units as in silico models of neurons. We show that XDream is applicable across network layers, architectures, and training sets; examine design choices in the algorithm; and provide practical guides for choosing hyperparameters in the algorithm. XDream is an efficient algorithm for uncovering neuronal tuning preferences in black-box networks using a vast and diverse stimulus space.


title: Local metrics of the Gaussian free field

abstract: We introduce the concept of a local metric of the Gaussian free field (GFF) $h$, which is a random metric coupled with $h$ in such a way that it depends locally on $h$ in a certain sense. This definition is a metric analog of the concept of a local set for $h$. We establish general criteria for two local metrics of the same GFF $h$ to be bi-Lipschitz equivalent to each other and for a local metric to be a.s. determined by $h$. Our results are used in subsequent works which prove the existence, uniqueness, and basic properties of the $\gamma$-Liouville quantum gravity (LQG) metric for all $\gamma \in (0,2)$, but no knowledge of LQG is needed to understand this paper.


title: Weak LQG metrics and Liouville first passage percolation

abstract: For $\gamma \in (0,2)$, we define a weak $\gamma$-Liouville quantum gravity (LQG) metric to be a function $h\mapsto D_h$ which takes in an instance of the planar Gaussian free field (GFF) and outputs a metric on the plane satisfying a certain list of natural axioms. We show that these axioms are satisfied for any subsequential limits of Liouville first passage percolation. Such subsequential limits were proven to exist by Ding-Dub\'edat-Dunlap-Falconet (2019). It is also known that these axioms are satisfied for the $\sqrt{8/3}$-LQG metric constructed by Miller and Sheffield (2013-2016). For any weak $\gamma$-LQG metric, we obtain moment bounds for diameters of sets as well as point-to-point, set-to-set, and point-to-set distances. We also show that any such metric is locally bi-H\"older continuous with respect to the Euclidean metric and compute the optimal H\"older exponents in both directions. Finally, we show that LQG geodesics cannot spend a long time near a straight line or the boundary of a metric ball. These results are used in subsequent work by Gwynne and Miller which proves that the weak $\gamma$-LQG metric is unique for each $\gamma \in (0,2)$, which in turn gives the uniqueness of the subsequential limit of Liouville first passage percolation. However, most of our results are new even in the special case when $\gamma=\sqrt{8/3}$.


title: Confluence of geodesics in Liouville quantum gravity for $\gamma \in (0,2)$

abstract: We prove that for any metric which one can associate with a Liouville quantum gravity (LQG) surface for $\gamma \in (0,2)$ satisfying certain natural axioms, its geodesics exhibit the following confluence property. For any fixed point $z$, a.s.\ any two $\gamma$-LQG geodesics started from distinct points other than $z$ must merge into each other and subsequently coincide until they reach $z$. This is analogous to the confluence of geodesics property for the Brownian map proven by Le Gall (2010). Our results apply for the subsequential limits of Liouville first passage percolation and are an important input in the proof of the existence and uniqueness of the LQG metric for all $\gamma\in (0,2)$.


title: The Higgs program and open questions in particle physics and cosmology

abstract: The Higgs program is relevant to many of the open fundamental questions in particle physics and in cosmology. Thus, when discussing future collider experiments, one way of comparing them is by assessing their potential contributions to progress on these questions. We discuss in detail the capabilities of the various proposed experiments in searching for singlet scalars, which are relevant to several of the open questions, and in measuring Higgs decays to fermion pairs, which are relevant to the flavor puzzles. On other interesting questions, we list the most relevant observables within the Higgs program.


title: Existence and uniqueness of the Liouville quantum gravity metric for $\gamma \in (0,2)$

abstract: We show that for each $\gamma \in (0,2)$, there is a unique metric (i.e., distance function) associated with $\gamma$-Liouville quantum gravity (LQG). More precisely, we show that for the whole-plane Gaussian free field (GFF) $h$, there is a unique random metric $D_h$ associated with the Riemannian metric tensor "$e^{\gamma h} (dx^2 + dy^2)$" on $\mathbb C$ which is characterized by a certain list of axioms: it is locally determined by $h$ and it transforms appropriately when either adding a continuous function to $h$ or applying a conformal automorphism of $\mathbb C$ (i.e., a complex affine transformation). Metrics associated with other variants of the GFF can be constructed using local absolute continuity. The $\gamma$-LQG metric can be constructed explicitly as the scaling limit of Liouville first passage percolation (LFPP), the random metric obtained by exponentiating a mollified version of the GFF. Earlier work by Ding, Dub\'edat, Dunlap, and Falconet (2019) showed that LFPP admits non-trivial subsequential limits. This paper shows that the subsequential limit is unique and satisfies our list of axioms. In the case when $\gamma = \sqrt{8/3}$, our metric coincides with the $\sqrt{8/3}$-LQG metric constructed in previous work by Miller and Sheffield, which in turn is equivalent to the Brownian map for a certain variant of the GFF. For general $\gamma \in (0,2)$, we conjecture that our metric is the Gromov-Hausdorff limit of appropriate weighted random planar map models, equipped with their graph distance. We include a substantial list of open problems.


title: Conformal covariance of the Liouville quantum gravity metric for $\gamma \in (0,2)$

abstract: For $\gamma \in (0,2)$, $U\subset \mathbb C$, and an instance $h$ of the Gaussian free field (GFF) on $U$, the $\gamma$-Liouville quantum gravity (LQG) surface associated with $(U,h)$ is formally described by the Riemannian metric tensor $e^{\gamma h} (dx^2 + dy^2)$ on $U$. Previous work by the authors showed that one can define a canonical metric (distance function) $D_h$ on $U$ associated with a $\gamma$-LQG surface. We show that this metric is conformally covariant in the sense that it respects the coordinate change formula for $\gamma$-LQG surfaces. That is, if $U,\widetilde{U}$ are domains, $\phi \colon U \to \widetilde{U}$ is a conformal transformation, $Q=2/\gamma+\gamma/2$, and $\widetilde h = h\circ\phi^{-1} + Q\log|(\phi^{-1})'|$, then $D_h(z,w) = D_{\widetilde{h}}(\phi(z),\phi(w))$ for all $z,w \in U$. This proves that $D_h$ is intrinsic to the quantum surface structure of $(U,h)$, i.e., it does not depend on the particular choice of parameterization.


title: Permanence and universal classification of discrete-time competitive systems via the carrying simplex

abstract: We study the permanence and impermanence for discrete-time Kolmogorov systems admitting a carrying simplex. Sufficient conditions to guarantee permanence and impermanence are provided based on the existence of a carrying simplex. Particularly, for low-dimensional systems, permanence and impermanence can be determined by boundary fixed points. For a class of competitive systems whose fixed points are determined by linear equations, there always exists a carrying simplex. We provide a universal classification via the equivalence relation relative to local dynamics of boundary fixed points for the three-dimensional systems by the index formula on the carrying simplex. There are a total of 33 stable equivalence classes which are described in terms of inequalities on parameters, and we present the phase portraits on their carrying simplices. Moreover, every orbit converges to some fixed point in classes 1-25 and 33; there is always a heteroclinic cycle in class 27; Neimark-Sacker bifurcations may occur in classes 26-31 but cannot occur in class 32. Based on our permanence criteria and the equivalence classification, we obtain the specific conditions on parameters for permanence and impermanence. Only systems in classes 29,31,33 and those in class 27 with a repelling heteroclinic cycle are permanent. Applications to discrete population models including the Leslie-Gower models, Atkinson-Allen models and Ricker models are given.


title: Strong magnetophonon oscillations in extra-large graphene

abstract: Van der Waals materials and their heterostructures offer a versatile platform for studying a variety of quantum transport phenomena due to their unique crystalline properties and the exceptional ability in tuning their electronic spectrum. However, most experiments are limited to devices that have lateral dimensions of only a few micrometres. Here, we perform magnetotransport measurements on graphene/hexagonal boron-nitride Hall bars and show that wider devices reveal additional quantum effects. In devices wider than ten micrometres we observe distinct magnetoresistance oscillations that are caused by resonant scattering of Landau-quantised Dirac electrons by acoustic phonons in graphene. The study allows us to accurately determine graphene's low energy phonon dispersion curves and shows that transverse acoustic modes cause most of phonon scattering. Our work highlights the crucial importance of device width when probing quantum effects and also demonstrates a precise, spectroscopic method for studying electron-phonon interactions in van der Waals heterostructures.


title: Spacecraft Attitude Control with Nonconvex Constraints: An Explicit Reference Governor Approach

abstract: This paper introduces a novel attitude controller for spacecraft subject to actuator saturation and multiple exclusion cone constraints. The proposed solution relies on a two-layer approach where the first layer prestabilizes the system dynamics whereas the second layer enforces constraint satisfaction by suitably manipulating the reference of the prestabilized system. In particular, constraint satisfaction is guaranteed by taking advantage of set invariance properties, whereas asymptotic convergence is achieved by implementing a non-conservative navigation field which is devoid of undesired stagnation points. Multiple numerical examples illustrate the good behavior of the proposed scheme.


title: Weinstein Homotopies

abstract: We give a partial answer to a question asked by Eliashberg in one of his recent papers.


title: The impact of incoming preparation and demographics on performance in Physics I: a multi-institution comparison

abstract: We have studied the impact of incoming preparation and demographic variables on student performance on the final exam in physics 1, the standard introductory, calculus-based mechanics course This was done at three different institutions using multivariable regression analysis to determine the extent to which exam scores can be predicted by a variety of variables that are available to most faculty and departments. We have found that the results are surprisingly consistent across the institutions, with the only two variables that have predictive power being math SAT/ACT scores and concept inventory pre-scores. The importance of both variables is comparable and fairly similar across the institutions. They explain 20 - 30 percent of the variation in students' performance on the final exam. Most notably, the demographic variables (gender, under-represented minority, first generation to attend college) are not significant. In all cases, although there appear to be gaps in exam performance if one considers only the demographic variable, once these two proxies of incoming preparation are included in the model, there is no longer a demographic gap. There is only a preparation gap that applies equally across the entire student population. This work shows that to properly understand differences in student performance across a diverse population, and hence to design more effective instruction, it is important to do statistical analyses that take multiple variables into account. It also illustrates the importance of having measures that are sensitive to both subject specific and more general preparation. The results suggest that better matching of the course design and teaching to the incoming student preparation will likely be the most effective way to eliminate observed performance gaps across demographic groups while also improving the success of all students.


title: Interfacing PDM MEMS microphones with PFM spiking systems: Application for Neuromorphic Auditory Sensors

abstract: In neuromorphic engineering, computation is commonly performed asynchronously, mimicking the way in which nervous systems process information: spike by spike. The Neuromorphic Auditory Sensor (NAS) has been implemented under this principle: applying different spike-based Signal Processing blocks. Computation in the spike domain requires the conversion of signals from analog or digital representation to the spike domain, which could present a speed constraint in many cases. This paper presents a spike-based system to convert audio information from low-power pulse density modulation (PDM) MicroElectroMechanical Systems (MEMS) microphones into rate coded spike frequencies. These spikes represent the input signal of the NAS, avoiding the analog or digital to spike conversion, and therefore improving the time response of the NAS. This conversion has been done in VHDL as an interface for PDM microphones, converting their pulses into temporal distributed spikes following a pulse frequency modulation (PFM) scheme with an accurate Inter-Spike-Interval, known as "PDM to spikes interface" (PSI). This was tested in two scenarios, first as a stand-alone circuit for its characterization, and then integrated with a NAS for verification. The PSI achieves a Total Harmonic Distortion (THD) of -39.51dB and a Signal-to-Noise Ratio (SNR) of 59.12dB, demands less than 1\% of the resources of a Spartan-6 FPGA and has a power consumption below 5mW.


title: Estimation of Tissue Oxygen Saturation from RGB images and Sparse Hyperspectral Signals based on Conditional Generative Adversarial Network

abstract: Purpose: Intra-operative measurement of tissue oxygen saturation (StO2) is important in the detection of ischemia, monitoring perfusion and identifying disease. Hyperspectral imaging (HSI) measures the optical reflectance spectrum of the tissue and uses this information to quantify its composition, including StO2. However, real-time monitoring is difficult due to the capture rate and data processing time. Methods: An endoscopic system based on a multi-fiber probe was previously developed to sparsely capture HSI data (sHSI). These were combined with RGB images, via a deep neural network, to generate high-resolution hypercubes and calculate StO2. To improve accuracy and processing speed, we propose a dual-input conditional generative adversarial network (cGAN), Dual2StO2, to directly estimate StO2 by fusing features from both RGB and sHSI. Results: Validation experiments were carried out on in vivo porcine bowel data, where the ground truth StO2 was generated from the HSI camera. The performance was also compared to our previous super-spectral-resolution network, SSRNet in terms of mean StO2 prediction accuracy and structural similarity metrics. Dual2StO2 was also tested using simulated probe data with varying fiber number. Conclusions: StO2 estimation by Dual2StO2 is visually closer to ground truth in general structure, achieves higher prediction accuracy and faster processing speed than SSRNet. Simulations showed that results improved when a greater number of fibers are used in the probe. Future work will include refinement of the network architecture, hardware optimization based on simulation results, and evaluation of the technique in clinical applications beyond StO2 estimation.


title: Contextual bound states for qudit magic state distillation

abstract: Identifying necessary and sufficient conditions for universal quantum computing is a long-standing open problem for which contextuality is, perhaps, the only promising solution [Howard et al., Nature 510, 351 (2014)]. To justify this conjecture, Howard et al. showed that contextuality is equivalent to Wigner negativity for qudits, and is therefore necessary and possibly sufficient for qudit magic state distillation. Here, we reformulate magic state distillation in the language of discrete phase space, to show that any finite qudit magic state distillation routine contains bound states outside the Wigner polytope -- these states exhibit contextuality but are useless for magic state distillation.


title: Probabilistic Predictive Principal Component Analysis for Spatially-Misaligned and High-Dimensional Air Pollution Data with Missing Observations

abstract: Accurate predictions of pollutant concentrations at new locations are often of interest in air pollution studies on fine particulate matters (PM$_{2.5}$), in which data is usually not measured at all study locations. PM$_{2.5}$ is also a mixture of many different chemical components. Principal component analysis (PCA) can be incorporated to obtain lower-dimensional representative scores of such multi-pollutant data. Spatial prediction can then be used to estimate these scores at new locations. Recently developed predictive PCA modifies the traditional PCA algorithm to obtain scores with spatial structures that can be well predicted at unmeasured locations. However, these approaches require complete data, whereas multi-pollutant data tends to have complex missing patterns in practice. We propose probabilistic versions of predictive PCA which allow for flexible model-based imputation that can account for spatial information and subsequently improve the overall predictive performance.


title: Overcharging a multi black hole system and cosmic censorship

abstract: We study the generalization of the gadenken experiment of overcharging an extremal black hole proposed by Wald in the context of a multi black hole solution. In particular, we attempt to overcharge a system of two extremal black holes via test particle absorption to produce a system involving a black hole and a naked singularity. If such a process is possible, then this would be a potential violation of the cosmic censorship hypothesis. However, we find that, analogous to Wald's result for a single charged black hole, such a test particle which can expose the singularity, would not be able to enter the horizon. This provides an interesting and non-trivial example that supports the validity of the cosmic censorship hypothesis in four-dimensional general relativity.


title: Type Ia supernovae from dark matter core collapse

abstract: Dark matter (DM) which sufficiently heats a local region in a white dwarf will trigger runaway fusion, igniting a type Ia supernova (SN). In a companion paper, this instability was used to constrain DM heavier than $10^{16}$ GeV which ignites SN through the violent interaction of one or two individual DM particles with the stellar medium. Here we study the ignition of supernovae by the formation and self-gravitational collapse of a DM core containing many DM particles. For non-annihilating DM, such a core collapse may lead to a mini black hole that can ignite SN through the emission of Hawking radiation, or possibly as a by-product of accretion. For annihilating DM, core collapse leads to an increasing annihilation rate and can ignite SN through a large number of rapid annihilations. These processes extend the previously derived constraints on DM to masses as low as $10^{5}$ GeV.


title: Sex-Prediction from Periocular Images across Multiple Sensors and Spectra

abstract: In this paper, we provide a comprehensive analysis of periocular-based sex-prediction (commonly referred to as gender classification) using state-of-the-art machine learning techniques. In order to reflect a more challenging scenario where periocular images are likely to be obtained from an unknown source, i.e. sensor, convolutional neural networks are trained on fused sets composed of several near-infrared (NIR) and visible wavelength (VW) image databases. In a cross-sensor scenario within each spectrum an average classification accuracy of approximately 85\% is achieved. When sex-prediction is performed across spectra an average classification accuracy of about 82\% is obtained. Finally, a multi-spectral sex-prediction yields a classification accuracy of 83\% on average. Compared to proposed works, obtained results provide a more realistic estimation of the feasibility to predict a subject's sex from the periocular region.


title: Fast AutoAugment

abstract: Data augmentation is an essential technique for improving generalization ability of deep learning models. Recently, AutoAugment has been proposed as an algorithm to automatically search for augmentation policies from a dataset and has significantly enhanced performances on many image recognition tasks. However, its search method requires thousands of GPU hours even for a relatively small dataset. In this paper, we propose an algorithm called Fast AutoAugment that finds effective augmentation policies via a more efficient search strategy based on density matching. In comparison to AutoAugment, the proposed algorithm speeds up the search time by orders of magnitude while achieves comparable performances on image recognition tasks with various models and datasets including CIFAR-10, CIFAR-100, SVHN, and ImageNet.


title: Hydrodynamics of broken global symmetries in the bulk

abstract: We consider holographic theories at finite temperature in which a continuous global symmetry in the bulk is spontaneously broken. We study the linear response of operators in a regime which is dual to time dependent, long wavelength deformations of solutions generated by the symmetry. By computing the boundary theory retarded Green's function we show the existence of a gapless mode with a diffusive dispersion relation. The diffusive character of the mode is compatible with the absence of a conserved charge from the field theory point of view. We give an analytic expression for the corresponding diffusion constant in terms of thermodynamic data and a new transport coefficient $\sigma_{b}$ which is fixed by the black hole horizon data. After adding a perturbative source on the boundary, we compute the resulting gap $\delta\omega_{g}$ as a simple function of $\sigma_{b}$ and of data of the thermal state.


title: Network Calculus-based Timing Analysis of AFDX networks incorporating multiple TSN/BLS traffic classes

abstract: We propose a formal timing analysis of an extension of the AFDX standard, incorporating the TSN/BLS shaper, to homogenize the avionics communication architecture, and enable the interconnection of different avionics domains with mixed-criticality levels, e.g., current AFDX traffic, Flight Control and In-Flight Entertainment. Existing Network Calculus models are limited to three classes, but applications with heterogeneous traffic require additional classes. Hence, we propose to generalize an existing Network Calculus model to do a worst-case timing analysis of an architecture with multiple BLS on multi-hop networks, to infer real-time bounds. Then, we conduct the performance analysis of such a proposal. First we evaluate the model on a simple 3-classes single-hop network to assess the sensitivity and tightness of the model, and compare it to existing models (CPA and Network Calculus). Secondly, we study a realistic AFDX configuration with six classes and two BLS. Finally, we compute a real use-case to add A350 flight control traffic to the AFDX. Results show the good properties of the generalized Network Calculus model compared to the CPA model and the efficiency of the extended AFDX to noticeably enhance the medium priority level delay bounds, while respecting the higher priority level constraints, in comparison with the current AFDX standard.


title: Evolution of structural, magnetic and transport properties in MnBi2-xSbxTe4

abstract: Here we report the evolution of structural, magnetic and transport properties in MnBi$_{2-x}$Sb$_x$Te$_4$ (0$\leq x \leq$2) single crystals. MnSb$_2$Te$_4$, isostructural to MnBi$_2$Te$_4$, has the lattice parameters of \textit{a}=4.2445(3)$\AA$ and \textit{c}=40.869(5)$\AA$, respectively. With increasing Sb content in MnBi$_{2-x}$Sb$_x$Te$_4$, the \textit{a}-lattice decreases linearly following the Vegards law while the \textit{c}-lattice shows little compositional dependence. The \textit{a}-lattice contraction occurs by reducing Mn-Te-Mn bond angle while Mn-Te bond length remains nearly constant. The anisotropic magnetic properties suggest an antiferromagnetic order below T$_N$=19\,K for MnSb$_2$Te$_4$ with the magnetic moments aligned along the crystallographic \textit{c}-axis. The antiferromagnetic ordering temperature slightly decreases from 24\,K for MnBi$_2$Te$_4$ to 19\,K for MnSb$_2$Te$_4$. More dramatic change was observed for the critical magnetic fields required for the spin-flop transition and moment saturation. With increasing Sb content, both critical fields decrease and in MnSb$_2$Te$_4$ a small field of 3\,kOe is enough to saturate the moment. In high magnetic fields, the saturation moment shows significant suppression from 3.56$\mu_B$/Mn for MnBi$_2$Te$_4$ to 1.57$\mu_B$/Mn for MnSb$_2$Te$_4$. Data analyses suggest that both the interlayer magnetic interaction and single ion anisotropy decrease with increasing Sb content. The partial substitution of Bi by Sb also dramatically affects the transport properties. A crossover from n-type to p-type conducting behavior is observed around x=0.63. Our results show close correlation between structural, magnetic and transport properties in MnBi$_{2-x}$Sb$_x$Te$_4$ and that partial substitution of Bi by Sb is an effective approach to fine tuning both the magnetism and transport properties of MnBi$_{2-x}$Sb$_x$Te$_4$.


title: Learn Stereo, Infer Mono: Siamese Networks for Self-Supervised, Monocular, Depth Estimation

abstract: The field of self-supervised monocular depth estimation has seen huge advancements in recent years. Most methods assume stereo data is available during training but usually under-utilize it and only treat it as a reference signal. We propose a novel self-supervised approach which uses both left and right images equally during training, but can still be used with a single input image at test time, for monocular depth estimation. Our Siamese network architecture consists of two, twin networks, each learns to predict a disparity map from a single image. At test time, however, only one of these networks is used in order to infer depth. We show state-of-the-art results on the standard KITTI Eigen split benchmark as well as being the highest scoring self-supervised method on the new KITTI single view benchmark. To demonstrate the ability of our method to generalize to new data sets, we further provide results on the Make3D benchmark, which was not used during training.


title: Next-Paradigm Programming Languages: What Will They Look Like and What Changes Will They Bring?

abstract: The dream of programming language design is to bring about orders-of-magnitude productivity improvements in software development tasks. Designers can endlessly debate on how this dream can be realized and on how close we are to its realization. Instead, I would like to focus on a question with an answer that can be, surprisingly, clearer: what will be the common principles behind next-paradigm, high-productivity programming languages, and how will they change everyday program development? Based on my decade-plus experience of heavy-duty development in declarative languages, I speculate that certain tenets of high-productivity languages are inevitable. These include, for instance, enormous variations in performance (including automatic transformations that change the asymptotic complexity of algorithms); a radical change in a programmer's workflow, elevating testing from a near-menial task to an act of deep understanding; a change in the need for formal proofs; and more.


title: Green Functions, Sommerfeld Images, and Wormholes

abstract: Electrostatic Green functions for grounded equipotential circular and elliptical rings, and grounded hyperspheres in n-dimension electrostatics, are constructed using Sommerfeld's method. These electrostatic systems are treated geometrically as different radial p-norm wormhole metrics that are deformed to be the Manhattan norm, namely "squashed wormholes". Differential geometry techniques are discussed to show how Riemannian geometry plays a rule in Sommerfeld's method. A comparison is made in terms of strength and position of the image charges for Sommerfeld's method with those for the more conventional Kelvin's method. Both methods are shown to be mathematically equivalent in terms of the corresponding Green functions. However, the two methods provide different physics perspectives, especially when studying different limits of those electrostatic systems. Further studies of ellipsoidal cases are suggested.


title: Observation of Complete Photonic Bandgap in Low Refractive Index Contrast Inverse Rod-Connected Diamond Structured Chalcogenides

abstract: Three-dimensional complete photonic bandgap materials or photonic crystals block light propagation in all directions. The rod-connected diamond structure exhibits the largest photonic bandgap known to date and supports a complete bandgap for the lowest refractive index contrast ratio down to $n_{high}/n_{low} \sim 1.9$. We confirm this threshold by measuring a complete photonic bandgap in the infrared region in Sn--S--O $(n\sim1.9)$ and Ge--Sb--S--O $(n\sim2)$ inverse rod-connected diamond structures. The structures were fabricated using a low-temperature chemical vapor deposition process via a single-inversion technique. This provides a reliable fabrication technique of complete photonic bandgap materials and expands the library of backfilling materials, leading to a wide range of future photonic applications.


title: LDPC Coded Multiuser Shaping for the Gaussian Multiple Access Channel

abstract: The joint design of input constellation and low-density parity-check (LDPC) codes to approach the symmetric capacity of the two-user Gaussian multiple access channel is studied. More specifically, multilevel coding is employed at each user to construct a high-order input constellation and the constellations of the users are jointly designed so as to maximize the multiuser shaping gain. At the receiver, each layer of the multilevel coding is jointly decoded among users, while successive cancellation is employed across layers. The LDPC code employed by each user in each layer is designed using EXIT charts to support joint decoding among users for the prescribed per-layer rate and SNR. Numerical simulations are provided to validate the proposed constellation and LDPC code designs.


title: Dynamic Origin-Destination Matrix Prediction with Line Graph Neural Networks and Kalman Filter

abstract: Modern intelligent transportation systems provide data that allow real-time dynamic demand prediction, which is essential for planning and operations. The main challenge of prediction of dynamic Origin-Destination (O-D) demand matrices is that demands cannot be directly measured by traffic sensors; instead, they have to be inferred from aggregate traffic flow data on traffic links. Specifically, spatial correlation, congestion and time dependent factors need to be considered in general transportation networks. In this paper we propose a novel O-D prediction framework combining heterogeneous prediction in graph neural networks and Kalman filter to recognize spatial and temporal patterns simultaneously. The underlying road network topology is converted into a corresponding line graph in the newly designed Fusion Line Graph Convolutional Networks (FL-GCNs), which provide a general framework of predicting spatial-temporal O-D flows from link information. Data from New Jersey Turnpike network are used to evaluate the proposed model. The results show that our proposed approach yields the best performance under various prediction scenarios. In addition, the advantage of combining deep neural networks and Kalman filter is demonstrated.


title: Recurrent strongly continuous operator families on Banach space

abstract: In this paper, we analyze recurrent $C_{0}$-semigroups of bounded operators on Banach spaces. We also introduce the notion of a (uniformly) $C_{0}$-rigid semigroups of bounded operators and give a structural characterization of them. A characterization of a Li-Yorke chaoticity of the translation semigroup $(T(t))_{t\geq 0}$ on weighted spaces of integrable functions and continuous functions in terms of admissible weight function is given. The recurrent $C_0$-semigroups induced by semiflows are characterized on the spaces of integrable functions and of spaces of continuous functions.


title: Tidal deformations of compact objects and gravitational wave emission

abstract: Observations of gravitational wave (GW) signals produced by coalescing binary neutron stars (NS), like the GW event GW170817, can be exploited to constrain the equation of state (EoS) of matter in the stars' inner core. The information on the internal structure and composition of the stars is encoded in their tidal Love numbers, which leave an imprint in the waveform of the GW signal emitted from the binary during the late inspiral phase. We extended the theory of tidal deformations of compact objects by computing the spin-tidal corrections that affect the dynamics and the GW emission of a binary system at the leading post-Newtonian (PN) order and to linear order in the spin. These corrections are divided into two classes: terms due to the coupling between the standard tidal Love numbers and the spins of the objects, and terms depending on the rotational tidal Love numbers. Both enter the GW phase at 6.5PN order. We analysed the impact of the spin-tidal couplings by estimating the parameter bias induced on GW170817-like events, assuming second- and third-generation ground based interferometers. If relatively high spinning ($\chi \gtrsim 0.1$) NS binaries exist in nature, these effects might be observed by the next generation of detectors. Lastly, we proved the feasibility of solving the so-called inverse stellar problem using GW detections, i.e., reconstructing the EoS from the measurement of NS masses and tidal Love numbers. Our results show that few observations of coalescing binary NS by a network of advanced detectors would allow us to put interesting constraints on the phenomenological parameters of a piecewise polytropic representation of the EoS, and to perform a model selection among the realistic EoS proposed in the literature.


title: Signals for Lorentz violation in gravitational waves

abstract: Lorentz violations in gravitational waves are investigated. Plane-wave solutions for arbitrary gauge-invariant violations in linearized gravity are constructed. Signatures of Lorentz violation include dispersion, birefringence, and anisotropies. Modifications to waves from coalescing compact binaries and to strain signals in gravitational-wave detectors are derived.


title: Minkowski Functionals of Convergence Maps and the Lensing Figure of Merit

abstract: Minkowski functionals (MFs) quantify the topological properties of a given field probing its departure from Gaussianity. We investigate their use on lensing convergence maps in order to see whether they can provide further insights on the underlying cosmology with respect to the standard second-order statistics, i.e., cosmic shear tomography. To this end, we first present a method to match theoretical predictions with measured MFs taking care of the shape noise, imperfections in the map reconstruction, and inaccurate description of the nonlinearities in the matter power spectrum and bispectrum. We validate this method against simulated maps reconstructed from shear fields generated by the MICE simulation. We then perform a Fisher matrix analysis to forecast the accuracy on cosmological parameters from a joint MFs and shear tomography analysis. It turns out that MFs are indeed helpful to break the $\Omega_{\rm m}$--$\sigma_8$ degeneracy thus generating a sort of chain reaction leading to an overall increase of the Figure of Merit.


title: An analysis of sparsity preserving pivot strategies for discontinuous Galerkin methods applied to acoustic scattering

abstract: In this paper we discuss and analyze the sparse structure of matrices associated to the interior penalty discontinuous Galerkin (IP-DG) method applied to the Helmholtz equation. It is well-known that LU-factorization causes fill-in and this paper discusses three pivoting strategies: approximate minimal degree (AMD), nested dissection, and reverse Cuthill-McKee, that can reduce fill-in associated to the LU-factorization. Numerical experiments are included that demonstrate the performance of these pivoting strategies. These experiments include both uniform and non-uniform mesh structures, the inclusion of a scattering boundary, and both piecewise linear and quadratic solution spaces.


title: Ricci-flat metrics on vector bundles over flag manifolds

abstract: We construct explicit complete Ricci-flat metrics on the total spaces of certain vector bundles over flag manifolds of the group $SU(n)$, for all K\"ahler classes. These metrics are natural generalizations of the metrics of Candelas-de la Ossa on the conifold, Pando Zayas-Tseytlin on the canonical bundle over $\mathbb{CP}^1\times \mathbb{CP}^1$, as well as the metrics on canonical bundles over flag manifolds, recently constructed by van Coevering.


title: Pushing the Boundaries of View Extrapolation with Multiplane Images

abstract: We explore the problem of view synthesis from a narrow baseline pair of images, and focus on generating high-quality view extrapolations with plausible disocclusions. Our method builds upon prior work in predicting a multiplane image (MPI), which represents scene content as a set of RGB$\alpha$ planes within a reference view frustum and renders novel views by projecting this content into the target viewpoints. We present a theoretical analysis showing how the range of views that can be rendered from an MPI increases linearly with the MPI disparity sampling frequency, as well as a novel MPI prediction procedure that theoretically enables view extrapolations of up to $4\times$ the lateral viewpoint movement allowed by prior work. Our method ameliorates two specific issues that limit the range of views renderable by prior methods: 1) We expand the range of novel views that can be rendered without depth discretization artifacts by using a 3D convolutional network architecture along with a randomized-resolution training procedure to allow our model to predict MPIs with increased disparity sampling frequency. 2) We reduce the repeated texture artifacts seen in disocclusions by enforcing a constraint that the appearance of hidden content at any depth must be drawn from visible content at or behind that depth. Please see our results video at: https://www.youtube.com/watch?v=aJqAaMNL2m4.


title: Similarity of Neural Network Representations Revisited

abstract: Recent work has sought to understand the behavior of neural networks by comparing representations between layers and between different trained models. We examine methods for comparing neural network representations based on canonical correlation analysis (CCA). We show that CCA belongs to a family of statistics for measuring multivariate similarity, but that neither CCA nor any other statistic that is invariant to invertible linear transformation can measure meaningful similarities between representations of higher dimension than the number of data points. We introduce a similarity index that measures the relationship between representational similarity matrices and does not suffer from this limitation. This similarity index is equivalent to centered kernel alignment (CKA) and is also closely connected to CCA. Unlike CCA, CKA can reliably identify correspondences between representations in networks trained from different initializations.


title: ALMA reveals the magnetic field evolution in the high-mass star forming complex G9.62+0.19

abstract: Context. The role of magnetic fields during the formation of high-mass stars is not yet fully understood, and the processes related to the early fragmentation and collapse are largely unexplored today. The high-mass star forming region G9.62+0.19 is a well known source, presenting several cores at different evolutionary stages. Aims. We determine the magnetic field morphology and strength in the high-mass star forming region G9.62+0.19, to investigate its relation to the evolutionary sequence of the cores. Methods. We use Band 7 ALMA observations in full polarisation mode and we analyse the polarised dust emission. We estimate the magnetic field strength via the Davis-Chandrasekhar-Fermi and the Structure Function methods. Results. We resolve several protostellar cores embedded in a bright and dusty filamentary structure. The polarised emission is clearly detected in six regions. Moreover the magnetic field is oriented along the filament and appears perpendicular to the direction of the outflows. We suggest an evolutionary sequence of the magnetic field, and the less evolved hot core exhibits a magnetic field stronger than the more evolved one. We detect linear polarisation from thermal line emission and we tentatively compared linear polarisation vectors from our observations with previous linearly polarised OH masers observations. We also compute the spectral index, the column density and the mass for some of the cores. Conclusions. The high magnetic field strength and the smooth polarised emission indicate that the magnetic field could play an important role for the fragmentation and the collapse process in the star forming region G9.62+019 and that the evolution of the cores can be magnetically regulated. On average, the magnetic field derived by the linear polarised emission from dust, thermal lines and masers is pointing in the same direction and has consistent strength.


title: RadiX-Net: Structured Sparse Matrices for Deep Neural Networks

abstract: The sizes of deep neural networks (DNNs) are rapidly outgrowing the capacity of hardware to store and train them. Research over the past few decades has explored the prospect of sparsifying DNNs before, during, and after training by pruning edges from the underlying topology. The resulting neural network is known as a sparse neural network. More recent work has demonstrated the remarkable result that certain sparse DNNs can train to the same precision as dense DNNs at lower runtime and storage cost. An intriguing class of these sparse DNNs is the X-Nets, which are initialized and trained upon a sparse topology with neither reference to a parent dense DNN nor subsequent pruning. We present an algorithm that deterministically generates RadiX-Nets: sparse DNN topologies that, as a whole, are much more diverse than X-Net topologies, while preserving X-Nets' desired characteristics. We further present a functional-analytic conjecture based on the longstanding observation that sparse neural network topologies can attain the same expressive power as dense counterparts


title: Viability of the $R+e^T$ cosmology

abstract: We propose a new theoretical approach for a cosmological model, which starts from an exponential of the trace of the energy-momentum tensor-dependence on the gravitational action, to be summed to the Ricci scalar. We derive the referred field equations and Friedmann-like equations. We derive the scale factor, Hubble parameter and deceleration parameter, in terms of both time and redshift. In possession of those parameters in terms of the redshift, we confront their predictions with the observational Hubble dataset and the outcomes are pretty satisfactory so that the model can be seen as a new alternative to the cosmological constant problem. We also present the statefinder diagnostic and discuss profoundly the dynamical behaviour of the model.


title: Optimal Storage Arbitrage under Net Metering using Linear Programming

abstract: We formulate the optimal energy arbitrage problem for a piecewise linear cost function for energy storage devices using linear programming (LP). The LP formulation is based on the equivalent minimization of the epigraph. This formulation considers ramping and capacity constraints, charging and discharging efficiency losses of the storage, inelastic consumer load and local renewable generation in presence of net-metering which facilitates selling of energy to the grid and incentivizes consumers to install renewable generation and energy storage. We consider the case where the consumer loads, electricity prices, and renewable generations at different instances are uncertain. These uncertain quantities are predicted using an Auto-Regressive Moving Average (ARMA) model and used in a model predictive control (MPC) framework to obtain the arbitrage decision at each instance. In numerical results we present the sensitivity analysis of storage performing arbitrage with varying ramping batteries and different ratio of selling and buying price of electricity.


title: Variational Bayesian Inference for Mixed Logit Models with Unobserved Inter- and Intra-Individual Heterogeneity

abstract: Variational Bayes (VB), a method originating from machine learning, enables fast and scalable estimation of complex probabilistic models. Thus far, applications of VB in discrete choice analysis have been limited to mixed logit models with unobserved inter-individual taste heterogeneity. However, such a model formulation may be too restrictive in panel data settings, since tastes may vary both between individuals as well as across choice tasks encountered by the same individual. In this paper, we derive a VB method for posterior inference in mixed logit models with unobserved inter- and intra-individual heterogeneity. In a simulation study, we benchmark the performance of the proposed VB method against maximum simulated likelihood (MSL) and Markov chain Monte Carlo (MCMC) methods in terms of parameter recovery, predictive accuracy and computational efficiency. The simulation study shows that VB can be a fast, scalable and accurate alternative to MSL and MCMC estimation, especially in applications in which fast predictions are paramount. VB is observed to be between 2.8 and 17.7 times faster than the two competing methods, while affording comparable or superior accuracy. Besides, the simulation study demonstrates that a parallelised implementation of the MSL estimator with analytical gradients is a viable alternative to MCMC in terms of both estimation accuracy and computational efficiency, as the MSL estimator is observed to be between 0.9 and 2.1 times faster than MCMC.


title: Restricted Connection Orthogonal Matching Pursuit For Sparse Subspace Clustering

abstract: Sparse Subspace Clustering (SSC) is one of the most popular methods for clustering data points into their underlying subspaces. However, SSC may suffer from heavy computational burden. Orthogonal Matching Pursuit applied on SSC accelerates the computation but the trade-off is the loss of clustering accuracy. In this paper, we propose a noise-robust algorithm, Restricted Connection Orthogonal Matching Pursuit for Sparse Subspace Clustering (RCOMP-SSC), to improve the clustering accuracy and maintain the low computational time by restricting the number of connections of each data point during the iteration of OMP. Also, we develop a framework of control matrix to realize RCOMP-SCC. And the framework is scalable for other data point selection strategies. Our analysis and experiments on synthetic data and two real-world databases (EYaleB & Usps) demonstrate the superiority of our algorithm compared with other clustering methods in terms of accuracy and computational time.


title: A Novel Trend Symbolic Aggregate Approximation for Time Series

abstract: Symbolic Aggregate approximation (SAX) is a classical symbolic approach in many time series data mining applications. However, SAX only reflects the segment mean value feature and misses important information in a segment, namely the trend of the value change in the segment. Such a miss may cause a wrong classification in some cases, since the SAX representation cannot distinguish different time series with similar average values but different trends. In this paper, we present Trend Feature Symbolic Aggregate approximation (TFSAX) to solve this problem. First, we utilize Piecewise Aggregate Approximation (PAA) approach to reduce dimensionality and discretize the mean value of each segment by SAX. Second, extract trend feature in each segment by using trend distance factor and trend shape factor. Then, design multi-resolution symbolic mapping rules to discretize trend information into symbols. We also propose a modified distance measure by integrating the SAX distance with a weighted trend distance. We show that our distance measure has a tighter lower bound to the Euclidean distance than that of the original SAX. The experimental results on diverse time series data sets demonstrate that our proposed representation significantly outperforms the original SAX representation and an improved SAX representation for classification.


title: Time-series Insights into the Process of Passing or Failing Online University Courses using Neural-Induced Interpretable Student States

abstract: This paper addresses a key challenge in Educational Data Mining, namely to model student behavioral trajectories in order to provide a means for identifying students most at-risk, with the goal of providing supportive interventions. While many forms of data including clickstream data or data from sensors have been used extensively in time series models for such purposes, in this paper we explore the use of textual data, which is sometimes available in the records of students at large, online universities. We propose a time series model that constructs an evolving student state representation using both clickstream data and a signal extracted from the textual notes recorded by human mentors assigned to each student. We explore how the addition of this textual data improves both the predictive power of student states for the purpose of identifying students at risk for course failure as well as for providing interpretable insights about student course engagement processes.


title: On the characterization of Danielewski surfaces by their automorphism group

abstract: In this note we show that if the automorphism group of a normal affine surface $S$ is isomorphic to the automorphism group of a Danielewski surface, then $S$ is isomorphic to a Danielewski surface.


title: An ADMM Based Framework for AutoML Pipeline Configuration

abstract: We study the AutoML problem of automatically configuring machine learning pipelines by jointly selecting algorithms and their appropriate hyper-parameters for all steps in supervised learning pipelines. This black-box (gradient-free) optimization with mixed integer & continuous variables is a challenging problem. We propose a novel AutoML scheme by leveraging the alternating direction method of multipliers (ADMM). The proposed framework is able to (i) decompose the optimization problem into easier sub-problems that have a reduced number of variables and circumvent the challenge of mixed variable categories, and (ii) incorporate black-box constraints along-side the black-box optimization objective. We empirically evaluate the flexibility (in utilizing existing AutoML techniques), effectiveness (against open source AutoML toolkits),and unique capability (of executing AutoML with practically motivated black-box constraints) of our proposed scheme on a collection of binary classification data sets from UCI ML& OpenML repositories. We observe that on an average our framework provides significant gains in comparison to other AutoML frameworks (Auto-sklearn & TPOT), highlighting the practical advantages of this framework.


title: Stochastic ordering results in parallel and series systems with Gumble distributed random variables

abstract: The stochastic comparisons of parallel and series system are worthy of study. In this paper, we present some stochastic comparisons of parallel and series systems having independent components from Gumble distribution with two parameters (one location and one shape). Here, we first put a condition for the likelihood ratio ordering of the parallel systems and second we use the concept of vector majorization technique to compare the systems by the reversed hazard rate ordering, the hazard rate ordering, the dispersive ordering, and the less uncertainty ordering with respect to the location parameter.


title: Large-scale Maps of the Cosmic Infrared Background from Planck

abstract: The cosmic infrared background (CIB) is a powerful probe of large-scale structure across a very large redshift range, and consists of unresolved redshifted infrared emission from dusty galaxies. It can be used to study the astrophysics of galaxies, the star formation history of the universe, and the connection between dark and luminous matter. It can furthermore be used as a tracer of the large-scale structure and thus assist in de-lensing of the cosmic microwave background. The major difficulty in its use lies in obtaining accurate and unbiased large-scale CIB images that are cleaned of the contamination by Galactic dust. We used data on neutral atomic hydrogen from the recently-released HI4PI Survey to create template maps of Galactic dust, allowing us to remove this component from the Planck intensity maps from 353 to 857 GHz for approximately $25\%$ of the sky. This allows us to constrain the CIB power spectrum down to $\ell\gtrsim 70$. We present these CIB maps and the various processing and validation steps that we have performed to ensure their quality, as well as a comparison with previous studies. All our data products are made publicly available at https://doi.org/10.7910/DVN/8A1SR3, thereby enabling the community to investigate a wide range of questions related to the universe's large-scale structure.


title: Black Hole and Neutron Star Binary Mergers in Triple Systems: II. Merger Eccentricity and Spin-Orbit Misalignment

abstract: We study the dynamical signatures of black hole (BH) and neutron star (NS) binary mergers via Lidov-Kozai oscillations induced by tertiary companions in hierarchical triple systems. For each type of binaries (BH-BH and BH-NS), we explore a wide range of binary/triple parameters that lead to binary mergers, and determine the distributions of eccentricity ($e_\mathrm{m}$) and spin-orbit misalignment angle ($\theta_\rm {sl}^\rm f$) when the binary enters the LIGO/VIRGO band. We use the double-averaged (over both orbits) and single-averaged (over the inner orbit) secular equations, as well as N-body integration, to evolve systems with different hierarchy levels, including the leading-order post-Newtonian effect, spin-orbit coupling and gravitational radiation. We find that for merging BH-BH binaries with comparable masses, about $7\%$ have $e_\mathrm{m}>0.1$ and $0.7\%$ have $e_{\rm m}>0.9$. The majority of the mergers have significant eccentricities in the LISA band. The BH spin evolution and $\theta_\rm {sl}^\rm f$ are correlated with the orbital evolution and $e_{\rm m}$. Mergers with $e_{\rm m} \lesssim 10^{-3}$ have a distribution of $\theta_\rm {sl}^\rm f$ that peaks around $90^\circ$ (and thus favoring a projected binary spin parameter $\chi_{\rm eff}\sim 0$), while mergers with larger $e_{\rm m}$ have a more isotropic spin-orbit misalignments. For typical BH-NS binaries, strong octuple effects lead to more mergers with non-negligible $e_\mathrm{m}$ (with $\sim 18\%$ having $e_\mathrm{m}>0.1$ and $2.5\%$ having $e_\mathrm{m}>0.9$), and the final BH spin axis tends to be randomly orientated. Measurements or constraints on eccentric mergers and $\theta_\rm {sl}^\rm f$ from LIGO/VIRGO and LISA would provide useful diagnostics on the dynamical formation of merging BH or NS binaries in triples.


title: Identification of RR Lyrae stars in multiband, sparsely-sampled data from the Dark Energy Survey using template fitting and Random Forest classification

abstract: Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (median 4 observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~31% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets.


title: "GammaMaP" - A Mathematica Package for Clifford Algebras, Gamma Matrices and Spinors

abstract: We present a Mathematica package for doing computations with gamma matrices, spinors, tensors and other objects, in any dimension and signature. The approach we use is based on defining the commutation relations of the relevant matrices, and is thus general and flexible. As examples, we reproduce torsion conditions for $AdS_3$ compactification of type IIB supergravity from literature, and check the vanishing of supersymmetry variations of 10-dimensional supersymmetric Yang-Mills theory.


title: Fluorescent rings in star-free dark matter haloes

abstract: Photoheating of the gas in low-mass dark matter (DM) haloes prevents baryons from cooling, leaving the haloes free of stars. Gas in these 'dark' haloes remains exposed to the ultraviolet background (UVB), and so is expected to emit via fluorescent recombination lines. We present a set of radiative transfer simulations, which model dark haloes as spherical gas clouds in hydrostatic equilibrium with a DM halo potential, and in thermal equilibrium with the UVB at redshift z = 0. We use these simulations to predict surface brightnesses in H-alpha, which we show to have a characteristic ring-shaped morphology for haloes in a narrow mass range between 10^9.5 and 10^9.6 M_sun. We explore how this emission depends on physical parameters such as the DM density profile and the UVB spectrum. We predict the abundance of fluorescent haloes on the sky, and discuss possible strategies for their detection. We demonstrate how detailed observations of fluorescent rings can be used to infer the properties of the haloes which host them, such as their density profiles and the mass-concentration relation, as well as to directly measure the UVB amplitude.


title: The [CII]--SFR correlation in dwarf galaxies across cosmic time

abstract: Current galaxy observations suggest that a roughly linear correlation exists between the [CII] emission and the star formation rate, either as spatially-resolved or integrated quantities. Observationally, this correlation seems to be independent of metallicity, but the very large scatter does not allow to properly assess whether this is true. On the other hand, theoretical models tend to suggest a metallicity dependence of the correlation. In this study, we investigate the metallicity evolution of the correlation via a high-resolution zoom-in cosmological simulation of a dwarf galaxy employing state-of-the-art sub-grid modelling for gas cooling, star formation, and stellar feedback, and that self-consistently evolves the abundances of metal elements out of equilibrium. Our results suggest that the correlation should evolve with metallicity, in agreement with theoretical predictions, but also that this evolution can be hardly detected in observations, because of the large scatter. We also find that most of the [CII] emission is associated with neutral gas at low-intermediate densities, whereas the highest emissivity is produced by the densest regions around star-forming regions.


title: Clues to the nature of dark matter from first galaxies

abstract: We use thirty-eight high-resolution simulations of galaxy formation between redshift 10 and 5 to study the impact of a 3 keV warm dark matter (WDM) candidate on the high-redshift Universe. We focus our attention on the stellar mass function and the global star formation rate and consider the consequences for reionization, namely the neutral hydrogen fraction evolution and the electron scattering optical depth. We find that three different effects contribute to differentiate warm and cold dark matter (CDM) predictions: WDM suppresses the number of haloes with mass less than few $10^9$ M$_{\odot}$; at a fixed halo mass, WDM produces fewer stars than CDM; and finally at halo masses below $10^9$ M$_{\odot}$, WDM has a larger fraction of dark haloes than CDM post-reionization. These three effects combine to produce a lower stellar mass function in WDM for galaxies with stellar masses at and below $\sim 10^7$ M$_{\odot}$. For $z > 7$, the global star formation density is lower by a factor of two in the WDM scenario, and for a fixed escape fraction, the fraction of neutral hydrogen is higher by 0.3 at $z \sim 6$. This latter quantity can be partially reconciled with CDM and observations only by increasing the escape fraction from 23 per cent to 34 per cent. Overall, our study shows that galaxy formation simulations at high redshift are a key tool to differentiate between dark matter candidates given a model for baryonic physics.


title: FIREBall-2: advancing TRL while doing proof-of-concept astrophysics on a suborbital platform

abstract: Here we discuss advances in UV technology over the last decade, with an emphasis on photon counting, low noise, high efficiency detectors in sub-orbital programs. We focus on the use of innovative UV detectors in a NASA astrophysics balloon telescope, FIREBall-2, which successfully flew in the Fall of 2018. The FIREBall-2 telescope is designed to make observations of distant galaxies to understand more about how they evolve by looking for diffuse hydrogen in the galactic halo. The payload utilizes a 1.0-meter class telescope with an ultraviolet multi-object spectrograph and is a joint collaboration between Caltech, JPL, LAM, CNES, Columbia, the University of Arizona, and NASA. The improved detector technology that was tested on FIREBall-2 can be applied to any UV mission. We discuss the results of the flight and detector performance. We will also discuss the utility of sub-orbital platforms (both balloon payloads and rockets) for testing new technologies and proof-of-concept scientific ideas


title: New Methods for Conformal Correlation Functions

abstract: The most general operator product expansion in conformal field theory is obtained using the embedding space formalism and a new uplift for general quasi-primary operators. The uplift introduced here, based on quasi-primary operators with spinor indices only and standard projection operators, allows a unified treatment of all quasi-primary operators irrespective of their Lorentz group irreducible representations. This unified treatment works at the level of the operator product expansion and hence applies to all correlation functions. A very useful differential operator appearing in the operator product expansion is established and its action on appropriate products of embedding space coordinates is explicitly computed. This computation leads to tensorial generalizations of the usual Exton function for all correlation functions. Several important identities and contiguous relations are also demonstrated for these new tensorial functions. From the operator product expansion all correlation functions for all quasi-primary operators, irrespective of their Lorentz group irreducible representations, can be computed recursively in a systematic way. The resulting answer can be expressed in terms of tensor structures that carry all the Lorentz group information and linear combinations of the new tensorial functions. Finally, a summary of the well-defined rules allowing the computation of all correlation functions constructively is presented.


title: Anomaly indicators for topological orders with $U(1)$ and time-reversal symmetry

abstract: We study anomalies in time-reversal ($\mathbb{Z}_2^T$) and $U(1)$ symmetric topological orders. In this context, an anomalous topological order is one that cannot be realized in a strictly $(2+1)$-D system but can be realized on the surface of a $(3+1)$-D symmetry-protected topological (SPT) phase. To detect these anomalies we propose several anomaly indicators -- functions that take as input the algebraic data of a symmetric topological order and that output a number indicating the presence or absence of an anomaly. We construct such indicators for both structures of the full symmetry group, i.e. $U(1)\rtimes\mathbb{Z}_2^T$ and $U(1)\times\mathbb{Z}_2^T$, and for both bosonic and fermionic topological orders. In all cases we conjecture that our indicators are complete in the sense that the anomalies they detect are in one-to-one correspondence with the known classification of $(3+1)$-D SPT phases with the same symmetry. We also show that one of our indicators for bosonic topological orders has a mathematical interpretation as a partition function for the bulk $(3+1)$-D SPT phase on a particular manifold and in the presence of a particular background gauge field for the $U(1)$ symmetry.


title: The bias of cosmic voids in the presence of massive neutrinos

abstract: Cosmic voids offer an extraordinary opportunity to study the effects of massive neutrinos on cosmological scales. Because they are freely streaming, neutrinos can penetrate the interior of voids more easily than cold dark matter or baryons, which makes their relative contribution to the mass budget in voids much higher than elsewhere in the Universe. In simulations it has recently been shown how various characteristics of voids in the matter distribution are affected by neutrinos, such as their abundance, density profiles, dynamics, and clustering properties. However, the tracers used to identify voids in observations (e.g., galaxies or halos) are affected by neutrinos as well, and isolating the unique neutrino signatures inherent to voids becomes more difficult. In this paper we make use of the DEMNUni suite of simulations to investigate the clustering bias of voids in Fourier space as a function of their core density and compensation. We find a clear dependence on the sum of neutrino masses that remains significant even for void statistics extracted from halos. In particular, we observe that the amplitude of the linear void bias increases with neutrino mass for voids defined in dark matter, whereas this trend gets reversed and slightly attenuated when measuring the relative void-halo bias using voids identified in the halo distribution. Finally, we argue how the original behaviour can be restored when considering observations of the total matter distribution (e.g. via weak lensing), and comment on scale-dependent effects in the void bias that may provide additional information on neutrinos in the future.


title: Bubble mapping with the Square Kilometer Array -- I. Detecting galaxies with Euclid, JWST, WFIRST and ELT within ionized bubbles in the intergalactic medium at z>6

abstract: The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.


title: Unitary Extension of Exotic Massive 3D Gravity from Bi-gravity

abstract: We obtain a new 3D gravity model from two copies of parity-odd Einstein-Cartan theories. Using Hamiltonian analysis, we demonstrate that the only local degrees of freedom are two massive spin-2 modes. Unitarity of the model in anti-de Sitter and Minkowski backgrounds can be satisfied for vast choices of the parameters without fine-tuning. The recent "exotic massive 3D gravity" model arises as a limiting case of the new model. We also show that there exist trajectories on the parameter space of the new model which cross the boundary between unitary and non-unitary regions. At the crossing point, one massive graviton decouples resulting in a unitary model with just one bulk degree of freedom but two positive central charges at odds with the usual expectation that the critical model has at least one vanishing central charge. Given the fact that a suitable non-relativistic version of bi-gravity has been used as an effective theory for gapped spin-2 fractional quantum Hall states, our model may have interesting applications in condensed matter physics.