title: Universality and individuality in neural dynamics across large populations of recurrent networks

abstract: Task-based modeling with recurrent neural networks (RNNs) has emerged as a popular way to infer the computational function of different brain regions. These models are quantitatively assessed by comparing the low-dimensional neural representations of the model with the brain, for example using canonical correlation analysis (CCA). However, the nature of the detailed neurobiological inferences one can draw from such efforts remains elusive. For example, to what extent does training neural networks to solve common tasks uniquely determine the network dynamics, independent of modeling architectural choices? Or alternatively, are the learned dynamics highly sensitive to different model choices? Knowing the answer to these questions has strong implications for whether and how we should use task-based RNN modeling to understand brain dynamics. To address these foundational questions, we study populations of thousands of networks, with commonly used RNN architectures, trained to solve neuroscientifically motivated tasks and characterize their nonlinear dynamics. We find the geometry of the RNN representations can be highly sensitive to different network architectures, yielding a cautionary tale for measures of similarity that rely representational geometry, such as CCA. Moreover, we find that while the geometry of neural dynamics can vary greatly across architectures, the underlying computational scaffold---the topological structure of fixed points, transitions between them, limit cycles, and linearized dynamics---often appears universal across all architectures.


title: Multiple ferromagnetic transitions and structural distortion in the van-der-Waals ferromagnet VI$_3$ at ambient and finite pressures

abstract: We present a combined study of zero-field $^{51}$V and $^{127}$I NMR at ambient pressure and specific heat and magnetization measurements under pressure up to 2.08 GPa on bulk single crystals of the van-der-Waals ferromagnet VI$_3$. At ambient pressure, our results consistently demonstrate that VI$_3$ undergoes a structural transition at $T_s \approx $78 K, followed by two subsequent ferromagnetic transitions at $T_{FM1} \approx $50 K and $T_{FM2} \approx $36 K upon cooling. At lowest temperature ($T < T_{FM2}$), two magnetically-ordered V sites exist, whereas only one magnetically-ordered V site is observed for $T_{FM1} < T\,< T_{FM2}$. Whereas $T_{FM1}$ is almost unaffected by external pressure, $T_{FM2}$ is highly responsive to pressure and merges with the $T_{FM1}$ line at $p \approx 0.6 $GPa. At even higher pressures ($p \approx $1.25\,GPa), the $T_{FM2}$ line merges with the structural transition at $T_s$ which becomes moderately suppressed with $p$ for $p < 1.25$ GPa. Taken together, our data point towards a complex magnetic structure and an interesting interplay of magnetic and structural degrees of freedom in VI$_3$.


title: Parton distributions from nonlocal chiral SU(3) effective theory: Flavor asymmetries

abstract: Using recently derived results for one-loop hadronic splitting functions from a nonlocal implementation of chiral effective theory, we study the contributions from pseudoscalar meson loops to flavor asymmetries in the proton. Constraining the parameters of the regulating functions by inclusive production of $n$, $\Delta^{++}$, $\Lambda$ and $\Sigma^{*+}$ baryons in $pp$ collisions, we compute the shape of the light antiquark asymmetry $\bar{d}-\bar{u}$ in the proton and the strange asymmetry $s-\bar{s}$ in the nucleon sea. With these constraints, the magnitude of the $\bar{d}-\bar{u}$ asymmetry is found to be compatible with that extracted from the Fermilab E866 Drell-Yan measurement, with no indication of a sign change at large values of $x$, and an integrated value in the range $\langle \bar d-\bar u \rangle \approx 0.09-0.17$. The $s-\bar s$ asymmetry is predicted to be positive at $x > 0$, with compensating negative contributions at $x=0$, and an integrated $x$-weighted moment in the range $\langle x (s-\bar s) \rangle \approx (0.9-2.5) \times 10^{-3}$.


title: Roots of generalised Hermite polynomials when both parameters are large

abstract: We study the roots of the generalised Hermite polynomials $H_{m,n}$ when both $m$ and $n$ are large. We prove that the roots, when appropriately rescaled, densely fill a bounded quadrilateral region, called the elliptic region, and organise themselves on a deformed rectangular lattice, as was numerically observed by Clarkson. We describe the elliptic region and the deformed lattice in terms of elliptic integrals and their degenerations. Keywords: Generalised Hermite polynomials; roots asymptotics; Painleve IV; Boutroux Curves; Tritronquee solution.


title: LightGuider: Guiding Interactive Lighting Design using Suggestions, Provenance, and Quality Visualization

abstract: LightGuider is a novel guidance-based approach to interactive lighting design, which typically consists of interleaved 3D modeling operations and light transport simulations. Rather than having designers use a trial-and-error approach to match their illumination constraints and aesthetic goals, LightGuider supports the process by simulating potential next modeling steps that can deliver the most significant improvements. LightGuider takes predefined quality criteria and the current focus of the designer into account to visualize suggestions for lighting-design improvements via a specialized provenance tree. This provenance tree integrates snapshot visualizations of how well a design meets the given quality criteria weighted by the designer's preferences. This integration facilitates the analysis of quality improvements over the course of a modeling workflow as well as the comparison of alternative design solutions. We evaluate our approach with three lighting designers to illustrate its usefulness.


title: Underlying topological Dirac nodal line mechanism of anomalously large electron-phonon coupling strength on Be (0001) surface

abstract: Beryllium was recently discovered to harbor a Dirac nodal line (DNL) in its bulk phase and the DNL-induced non-trivial drumhead-like surface states (DNSSs) on its (0001) surface, rationalizing several already-existing historic puzzles [Phys. Rev. Lett., \textbf{117}, 096401 (2016)]. However, to date the underlying mechanism, as to why its (0001) surface exhibits an anomalously large electron-phonon coupling effect ($\lambda_{e-ph}^s$ $\approx$ 1.0), remains unresolved. Here, by means of first-principles calculations we have evidenced that the coupling of the DNSSs with the phononic states mainly contributes to its novel surface \emph{e-ph} enhancement. Besides that the experimentally observed $\lambda_{e-ph}^s$ and the main Eliashberg coupling function (ECF) peaks have been reproduced well, we have decomposed the ECF, $\alpha^{2}$$F$(\emph{k},\textbf{\emph{q}};\emph{v}), and the \emph{e-ph} coupling strength $\lambda(\emph{k},\textbf{\emph{q}};\emph{v})$ as a function of each electron momentum (\emph{k}), each phonon momentum (\textbf{\emph{q}}) and each phonon mode ($v$), evidencing the robust connection between the DNSSs and both $\alpha^{2}$$F$(\emph{k},\textbf{\emph{q}};\emph{v}) and $\lambda(\emph{k},\textbf{\emph{q}};\emph{v})$. The results reveal the strong \emph{e-ph} coupling between the DNSSs and the phonon modes, which contributes over 80$\%$ of the $\lambda_{e-ph}^s$ coefficient on the Be (0001) surface. It highlights that the anomalously large \emph{e-ph} coefficient on the Be (0001) surface can be attributed to the presence of its DNL-induced DNSSs, clarifying the long-term debated mechanism.


title: Parametric Oscillation and amplification with gate controlled capacitor-within-capacitor

abstract: A Capacitor-within-Capacitor (CWC) is a nested electronic element that has two components: the cell (e.g., the outer capacitor) and the gate (e.g., the inner capacitor). The designations for the gate and cell may be interchanged. Here, an analysis and experimentation on a diode-interfaced parametric oscillations and amplification are described. By replacing the diode with a junction, made of doped nano-graphene films, we demonstrated a new structure whose doping may be electronically and chemically controlled.


title: D-GAN: Deep Generative Adversarial Nets for Spatio-Temporal Prediction

abstract: Spatio-temporal (ST) data for urban applications, such as taxi demand, traffic flow, regional rainfall is inherently stochastic and unpredictable. Recently, deep learning based ST prediction models are proposed to learn the ST characteristics of data. However, it is still very challenging (1) to adequately learn the complex and non-linear ST relationships; (2) to model the high variations in the ST data volumes as it is inherently dynamic, changing over time (i.e., irregular) and highly influenced by many external factors, such as adverse weather, accidents, traffic control, PoI, etc.; and (3) as there can be many complicated external factors that can affect the accuracy and it is impossible to list them explicitly. To handle the aforementioned issues, in this paper, we propose a novel deep generative adversarial network based model (named, D-GAN) for more accurate ST prediction by implicitly learning ST feature representations in an unsupervised manner. D-GAN adopts a GAN-based structure and jointly learns generation and variational inference of data. More specifically, D-GAN consists of two major parts: (1) a deep ST feature learning network to model the ST correlations and semantic variations, and underlying factors of variations and irregularity in the data through the implicit distribution modelling; (2) a fusion module to incorporate external factors for reaching a better inference. To the best our knowledge, no prior work studies ST prediction problem via deep implicit generative model and in an unsupervised manner. Extensive experiments performed on two real-world datasets show that D-GAN achieves more accurate results than traditional as well as deep learning based ST prediction methods.


title: Rearrangement of 2D aggregates of droplets under compression: signatures of the energy landscape from crystal to glass

abstract: We study signatures of the energy landscape's evolution through the crystal-to-glass transition by compressing 2D finite aggregates of oil droplets. Droplets of two distinct sizes are used to compose small aggregates in an aqueous environment. Aggregates range from perfectly ordered monodisperse single crystals to disordered bidisperse glasses. The aggregates are compressed between two parallel boundaries, with one acting as a force sensor. The compression force provides a signature of the aggregate composition and gives insight into the energy landscape. In particular, crystals dissipate all the stored energy through single catastrophic fracture events whereas the glassy aggregates break step-by-step. Remarkably, the yielding properties of the 2D aggregates are strongly impacted by even a small amount of disorder.


title: Eigenfunctions of the Fourier Transform with specified zeros

abstract: We give a unified description of the modular and quasi-modular functions used in Viazovska's proof of the best packing bounds in dimension 8 and the proof by Cohn, Kumar, Miller, Radchenko, and Viazovska of the best packing bound in dimension 24. We show that necessarily modular forms have to be used to obtain these results. We extend these constructions to arbitrary dimensions divisible by 4.


title: An Algorithm and Estimates for the Erd\H{o}s-Selfridge Function (work in progress)

abstract: Let $p(n)$ denote the smallest prime divisor of the integer $n$. Define the function $g(k)$ to be the smallest integer $>k+1$ such that $p(\binom{g(k)}{k})>k$. So we have $g(2)=6$ and $g(3)=g(4)=7$. In this paper we present the following new results on the Erd\H{o}s-Selfridge function $g(k)$: We present a new algorithm to compute the value of $g(k)$, and use it to both verify previous work and compute new values of $g(k)$, with our current limit being $$ g(323)= 1\ 69829\ 77104\ 46041\ 21145\ 63251\ 22499. $$ We define a new function $\hat{g}(k)$, and under the assumption of our Uniform Distribution Heuristic we show that $$ \log g(k) = \log \hat{g}(k) + O(\log k) $$ with high "probability". We also provide computational evidence to support our claim that $\hat{g}(k)$ estimates $g(k)$ reasonably well in practice. There are several open conjectures on the behavior of $g(k)$ which we are able to prove for $\hat{g}(k)$, namely that $$ 0.525\ldots +o(1) \quad \le \quad \frac{\log \hat{g}(k)}{k/\log k} \quad \le \quad 1+o(1), $$ and that $$ \limsup_{k\rightarrow\infty} \frac{\hat{g}(k+1)}{\hat{g}(k)}=\infty.$$ Let $G(x,k)$ count the number of integers $n\le x$ such that $p(\binom{n}{k})>k$. Unconditionally, we prove that for large $x$, $G(x,k)$ is asymptotic to $x/\hat{g}(k)$. And finally, we show that the running time of our new algorithm is at most $g(k) \exp[ -c (k\log\log k) /(\log k)^2 (1+o(1))]$ for a constant $c>0$.


title: The LAPW method with eigendecomposition based on the Hari--Zimmermann generalized hyperbolic SVD

abstract: In this paper we propose an accurate, highly parallel algorithm for the generalized eigendecomposition of a matrix pair $(H, S)$, given in a factored form $(F^{\ast} J F, G^{\ast} G)$. Matrices $H$ and $S$ are generally complex and Hermitian, and $S$ is positive definite. This type of matrices emerges from the representation of the Hamiltonian of a quantum mechanical system in terms of an overcomplete set of basis functions. This expansion is part of a class of models within the broad field of Density Functional Theory, which is considered the golden standard in condensed matter physics. The overall algorithm consists of four phases, the second and the fourth being optional, where the two last phases are computation of the generalized hyperbolic SVD of a complex matrix pair $(F,G)$, according to a given matrix $J$ defining the hyperbolic scalar product. If $J = I$, then these two phases compute the GSVD in parallel very accurately and efficiently.


title: Dynamics of Intermediate Measurements in Quantum Systems

abstract: We study the dynamics of a quantum system in which an intermediate property $m$ is measured in between initial and final measurements of two different non-commuting properties $a$ and $b$. Since this intermediate measurement must involve an interaction, we use this case to explore in more detail the dynamics of measurement and propose that the physics of this change is described by a unitary transformation parametrised by a random parameter, with the selection of $m$ described by a superposition of these unitaries. We prove a set of conditions which must hold in order for this superposition to remain unitary and then argue that in the case where $U$ is a matrix, it must be a random unitary matrix. We outline a numerical example with matrices and make some predictions based on the conditions we have found which are capable of experimental verification. We finish with a brief discussion of the complex phases which appear in the relations between non-commuting properties and their relevance for the dynamics associated with measurement.


title: First demonstration of ionization cooling by the Muon Ionization Cooling Experiment

abstract: High-brightness muon beams of energy comparable to those produced by state-of-the-art electron, proton and ion accelerators have yet to be realised. Such beams have the potential to carry the search for new phenomena in lepton-antilepton collisions to extremely high energy and also to provide uniquely well-characterised neutrino beams. A muon beam may be created through the decay of pions produced in the interaction of a proton beam with a target. To produce a high-brightness beam from such a source requires that the phase space volume occupied by the muons be reduced (cooled). Ionization cooling is the novel technique by which it is proposed to cool the beam. The Muon Ionization Cooling Experiment collaboration has constructed a section of an ionization cooling cell and used it to provide the first demonstration of ionization cooling. We present these ground-breaking measurements.


title: Synthetic anyons in noninteracting systems

abstract: Synthetic anyons can be implemented in a noninteracting many-body system, by using specially tailored localized (physical) probes, which supply the demanded nontrivial topology in the system. We consider the Hamiltonian for noninteracting electrons in two-dimensions (2D), in a uniform magnetic field, where the probes are external solenoids with a magnetic flux that is a fraction of the flux quantum. The Hamiltonian could also be implemented in an ultracold (fermionic) atomic gas in 2D, in a uniform synthetic magnetic field, where the probes are lasers giving rise to synthetic solenoid gauge potentials. We find analytically and numerically the ground state of this system when only the lowest Landau level states are occupied. It is shown that the ground state is anyonic in the coordinates of the probes. We show that these synthetic anyons cannot be considered as emergent quasiparticles. The fusion rules of synthetic anyons are discussed for different microscopic realizations of the fusion process.


title: The glow of annihilating dark matter in Omega Centauri

abstract: Dark matter (DM) is the most abundant material in the Universe, but has so far been detected only via its gravitational effects. Several theories suggest that pairs of DM particles can annihilate into a flash of light at gamma-ray wavelengths. While gamma-ray emission has been observed from environments where DM is expected to accumulate, such as the centre of our Galaxy, other high energy sources can create a contaminating astrophysical gamma-ray background, thus making DM detection difficult. In principle, dwarf galaxies around the Milky Way are a better place to look -- they contain a greater fraction of DM with no astrophysical gamma-ray background -- but they are too distant for gamma-rays to have been seen. A range of observational evidence suggests that Omega Centauri (omega Cen or NGC 5139), usually classified as the Milky Way's largest globular cluster, is really the core of a captured and stripped dwarf galaxy. Importantly, Omega Cen is ten times closer to us than known dwarfs. Here we show that not only does Omega Cen contain DM with density as high as compact dwarf galaxies, but also that it emits gamma-rays with an energy spectrum matching that expected from the annihilation of DM particles with mass 31$\pm$4 GeV (68\% confidence limit). No astrophysical sources have been found that would otherwise explain Omega Cen's gamma-ray emission, despite deep multi-wavelength searches. We anticipate our results to be the starting point for even deeper radio observations of Omega Cen. If multi-wavelength searches continue to find no astrophysical explanations, this pristine, nearby clump of DM will become the best place to study DM interactions through forces other than gravity.


title: Integrality of matrices, finiteness of matrix semigroups, and dynamics of linear and additive cellular automata

abstract: Let $\mathbb{K}$ be a finite commutative ring, and let $\mathbb{L}$ be a commutative $\mathbb{K}$-algebra. Let $A$ and $B$ be two $n \times n$-matrices over $\mathbb{L}$ that have the same characteristic polynomial. The main result of this paper states that the set $\left\{ A^0,A^1,A^2,\ldots\right\}$ is finite if and only if the set $\left\{ B^0,B^1,B^2,\ldots\right\}$ is finite. We apply this result to Cellular Automata (CA). Indeed, it gives a complete and easy-to-check characterization of sensitivity to initial conditions and equicontinuity for linear CA over the alphabet $\mathbb{K}^n$ for $\mathbb{K} = \mathbb{Z}/m\mathbb{Z}$ i.e., CA in which the local rule is defined by $n\times n$-matrices with elements in $\mathbb{Z}/m\mathbb{Z}$. To prove our main result, we derive an integrality criterion for matrices that is likely of independent interest. Namely, let $\mathbb{K}$ be any commutative ring (not necessarily finite), and let $\mathbb{L}$ be a commutative $\mathbb{K}$-algebra. Consider any $n \times n$-matrix $A$ over $\mathbb{L}$. Then, $A \in \mathbb{L}^{n \times n}$ is integral over $\mathbb{K}$ (that is, there exists a monic polynomial $f \in \mathbb{K}\left[t\right]$ satisfying $f\left(A\right) = 0$) if and only if all coefficients of the characteristic polynomial of $A$ are integral over $\mathbb{K}$. The proof of this fact relies on a strategic use of exterior powers (a trick pioneered by Gert Almkvist). Furthermore, we extend the decidability result concerning sensitivity and equicontinuity to the wider class of additive CA over a finite abelian group. For such CA, we also prove the decidability of injectivity, surjectivity, topological transitivity and all the properties (as, for instance, ergodicity) that are equivalent to the latter.


title: Clustering Higher Order Data: An Application to Pediatric Multi-variable Longitudinal Data

abstract: Physical activity levels are an important predictor of cardiovascular health and increasingly being measured by sensors, like accelerometers. Accelerometers produce rich multivariate data that can inform important clinical decisions related to individual patients and public health. The CHAMPION study, a study of youth with chronic inflammatory conditions, aims to determine the links between heart health, inflammation, physical activity, and fitness. The accelerometer data from CHAMPION is represented as 4-dimensional arrays, and a finite mixture of multidimensional arrays model is developed for clustering. The use of model-based clustering for multidimensional arrays has thus far been limited to two-dimensional arrays, i.e., matrices or order-two tensors, and the work in this paper can also be seen as an approach for clustering D-dimensional arrays for D > 2 or, in other words, for clustering order-D tensors.


title: Orientation control and nonlinear trajectory tracking of colloidal particles using microfluidics

abstract: Suspensions of anisotropic Brownian particles are commonly encountered in a wide array of applications such as drug delivery and manufacturing of fiber-reinforced composites. Technological applications and fundamental studies of small anisotropic particles critically require precise control of particle orientation over defined trajectories and paths. In this work, we demonstrate robust control over the two-dimensional (2D) center-of-mass position and orientation of anisotropic Brownian particles using only fluid flow. We implement a path-following model predictive control scheme to manipulate colloidal particles over defined trajectories in position space, where the speed of movement along the path is a degree of freedom in the controller design. We further explore how the external flow field affects the orientation dynamics of anisotropic particles in steady and transient extensional flow using a combination of experiments and analytical modeling. Overall, this technique offers new avenues for fundamental studies of anisotropic colloidal particles using only fluid flow, without the need for external electric or optical fields.


title: Linear Stability Analysis of a Normal Shock Train in a Constant Area Isolator of a Hypersonic ScramJet

abstract: Recent experiments in hypersonic vehicle technology have had some successful outcomes (ex. NASA X-43, Boeing X-51) but also have seen some unexpected failures (DARPA Hypersonic Test Vehicle 2). One important failure mode is the problem of engine unstart. Morgan, et al (2012) has discussed and analyzed a Scramjet isolator using two unit problems, viz. that of an oblique shock impinging on a turbulent boundary layer and that of a normal shock train in a constant-area duct. In order to analyze unstart, we use the numerical data obtained from the LES calculations of the normal shock train in a constant area duct and perform temporal linear stability analysis on it under the assumptions of inviscid semi-parallel flow. This paper describes the methodology and the details of this linear stability analysis. The results obtained from the investigations have been documented and interpreted with their relevance to the physical problem. The report concludes that temporal linear stability analysis is insufficient to analyze the problem of unstart. Finally, the need for a new approach (namely, dynamic mode decomposition (DMD)) has been discussed.


title: Practical trapped-ion protocols for universal qudit-based quantum computing

abstract: The notion of universal quantum computation can be generalized to multi-level qudits, which offer advantages in resource usage and algorithmic efficiencies. Trapped ions, which are pristine and well-controlled quantum systems, offer an ideal platform to develop qudit-based quantum information processing. Previous work has not fully explored the practicality of implementing trapped-ion qudits accounting for known experimental error sources. Here, we describe a universal set of protocols for state preparation, single-qudit gates, a new generalization of the M\o{}lmer-S\o{}rensen gate for two-qudit gates, and a measurement scheme which utilizes shelving to a meta-stable state. We numerically simulate known sources of error from previous trapped ion experiments, and show that there are no fundamental limitations to achieving fidelities above \(99\%\) for three-level qudits encoded in \(^{137}\mathrm{Ba}^+\) ions. Our methods are extensible to higher-dimensional qudits, and our measurement and single-qudit gate protocols can achieve \(99\%\) fidelities for five-level qudits. We identify avenues to further decrease errors in future work. Our results suggest that three-level trapped ion qudits will be a useful technology for quantum information processing.


title: Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

abstract: We present a numerically exact Inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.


title: Multi-component scalar fields and the complete factorization of its equations of motion

abstract: In the paper by Bazeia D. et al., EPL, 119 (2017) 61002, the authors demonstrate the equivalence between the second-order differential equation of motion and a family of first-order differential equations of Bogomolnyi type for the cases of single real and complex scalar field theories with non-canonical dynamics. The goal of this paper is to demonstrate that this equivalence is also valid for a more general classes of real scalar field models. We start the paper by demonstrating the equivalence in a single real scalar model. The first goal is to generalize the equivalence presented in papers by Bazeia et al. to a single real scalar field model without a specific form for its Lagrangian. The second goal is to use the setup presented in the first demonstration to show that this equivalence can be achieved also in a real multi-component scalar field model again without a specific form for its Lagrangian. The main goal of this paper is to show that this equivalence can be achieved in real scalar field scenarios that can be standard, or non-standard, with single, or multi-component, scalar fields.


title: Modeling the non-thermal emission of the gamma Cygni Supernova Remnant up to the highest energies

abstract: The gamma Cygni supernova remnant (SNR) is a middle-aged, Sedov-phase SNR in the Cygnus region. It is a known source of non-thermal emission at radio, X-ray, and gamma-ray energies. Very-high energy (VHE, >100 GeV) gamma-ray emission from gamma Cygni was first detected by the VERITAS observatory and it has since been observed by other experiments. Observations so far indicate that there must be a population of non-thermal particles present in the remnant which produces the observed emission. However, it is not clear what kind of particles (protons/ions or electrons) are accelerated in the remnant, how efficient the acceleration is, and up to which energy particles can be accelerated. Accurate measurements of the VHE gamma-ray spectrum are crucial to investigate particle acceleration above TeV energies. This presentation will focus on multi-wavelength observations of the gamma Cygni SNR and their interpretation. We will present improved measurements of the VHE gamma-ray emission spectrum of gamma Cygni by the High-Altitude Water Cherenkov (HAWC) gamma-ray observatory, and use these results as well as measurements from other instruments to model the underlying particle populations producing this emission. HAWC's excellent sensitivity at TeV energies and above enables us to extend spectral measurements to higher energies and better constrain the maximum acceleration energy.


title: Real-time evolution and quantized charge pumping in magnetic Weyl semimetals

abstract: Real-time evolution and charge pumping in magnetic Weyl semimetals are studied by solving the time-dependent Schr\"{o}dinger equations. In the adiabatic limit of the real-time evolution, we show that the total pumped charge is quantized in the magnetic Weyl semimetals as in the quantum Hall system although the Weyl semimetal has no bulk gap. We examine how the disorder affects the charge pumping. As a result, we show that the quantized pumped charge is robust against the small disorder and find that the pumped charge increases in the intermediate disorder region. We also examine the doping effects on the charge pumping and show that the remnant of the quantized pumped charge at zero doping can be detected. Our results show that the real-time evolution is a useful technique for detecting the topological properties of the systems with no bulk gap and/or disorders.


title: Quantifying coherence with respect to general quantum measurements

abstract: Coherence is a cornerstone of quantum theory and a prerequisite for the advantage of quantum technologies. In recent work, the notion of coherence with respect to a general quantum measurement (POVM) was introduced and embedded into a resource-theoretic framework that generalizes the standard resource theory of coherence. In particular, POVM-incoherent (free) states and operations were established. In this work, we explore features of this framework which arise due to the rich structure of POVMs compared to projective measurements. Moreover, we introduce a rigorous, probabilisitic framework for POVM-based coherence measures and free operations. This leads to the introduction of new, strongly monotonic resource measures that neatly generalize well-known standard coherence measures. Finally, we show that the relative entropy of POVM-coherence is equal to the cryptographic randomness gain, providing an important operational meaning to the concept of coherence with respect to a general measurement.


title: A Survey of TeV emission from Galactic Supernova Remnants with HAWC

abstract: Supernova remnants (SNRs) have long been hypothesized as the main source of Galactic Cosmic Rays up to PeV energies. Some of them have indeed been shown to accelerate protons to TeV energies and above. But which of them are indeed efficient accelerators of protons and nuclei? And up to which energies can they accelerate these particles? Measurements of non-thermal emission, especially in the X-ray and gamma-ray regimes, are essential to answer these questions. The High-Altitude Water Cherenkov (HAWC) observatory, surveying the northern TeV gamma-ray sky, is currently the most sensitive wide field-of-view survey instrument in the VHE (very-high-energy, $>$100 GeV) range. With more than three years of data recorded, HAWC is ideally suited for an unbiased survey of gamma-ray emission from galactic SNRs, particularly at TeV energies and above. In this proceeding, I will give an overview of recent measurements of VHE gamma-ray emission from SNRs with the HAWC Observatory. Combined with data from other wavelengths, these measurements are used to derive information about the underlying particle populations such as the maximum acceleration energy and whether leptonic or hadronic processes are responsible for the emission.


title: Transverse Temperature Interfaces in the Katz-Lebowitz-Spohn Driven Lattice Gas

abstract: We explore the intriguing spatial patterns that emerge in a two-dimensional spatially inhomogeneous Katz-Lebowitz-Spohn (KLS) driven lattice gas with attractive nearest-neighbor interactions. The domain is split into two regions with hopping rates governed by different temperatures T > T_c and T_c, respectively, where T_c indicates the critical temperature for phase ordering, and with the temperature boundaries oriented perpendicular to the drive. In the hotter region, the system behaves like the (totally) asymmetric exclusion processes (TASEP), and experiences particle blockage in front of the interface to the critical region. To explain this particle density accumulation near the interface, we have measured the steady-state current in the KLS model at T > T_c and found it to decay as 1/T. In analogy with TASEP systems containing "slow" bonds, we argue that transport in the high-temperature subsystem is impeded by the lower current in the cooler region, which tends to set the global stationary particle current value. This blockage is induced by the extended particle clusters, growing logarithmically with system size, in the critical region. We observe the density profiles in both high-and low-temperature subsystems to be similar to the well-characterized coexistence and maximal-current phases in (T)ASEP models with open boundary conditions, which are respectively governed by hyperbolic and trigonometric tangent functions. Yet if the lower temperature is set to T_c, we detect marked fluctuation corrections to the mean-field density profiles, e.g., the corresponding critical KLS power law density decay near the interfaces into the cooler region.


title: Learning Multimorbidity Patterns from Electronic Health Records Using Non-negative Matrix Factorisation

abstract: Multimorbidity, or the presence of several medical conditions in the same individual, has been increasing in the population, both in absolute and relative terms. However, multimorbidity remains poorly understood, and the evidence from existing research to describe its burden, determinants and consequences has been limited. Previous studies attempting to understand multimorbidity patterns are often cross-sectional and do not explicitly account for multimorbidity patterns' evolution over time; some of them are based on small datasets and/or use arbitrary and narrow age ranges; and those that employed advanced models, usually lack appropriate benchmarking and validations. In this study, we (1) introduce a novel approach for using Non-negative Matrix Factorisation (NMF) for temporal phenotyping (i.e., simultaneously mining disease clusters and their trajectories); (2) provide quantitative metrics for the evaluation of disease clusters from such studies; and (3) demonstrate how the temporal characteristics of the disease clusters that result from our model can help mine multimorbidity networks and generate new hypotheses for the emergence of various multimorbidity patterns over time. We trained and evaluated our models on one of the world's largest electronic health records (EHR), with 7 million patients, from which over 2 million where relevant to this study.


title: Testing with Fewer Resources: An Adaptive Approach to Performance-Aware Test Case Generation

abstract: Automated test case generation is an effective technique to yield high-coverage test suites. While the majority of research effort has been devoted to satisfying coverage criteria, a recent trend emerged towards optimizing other non-coverage aspects. In this regard, runtime and memory usage are two essential dimensions: less expensive tests reduce the resource demands for the generation process and later regression testing phases. This study shows that performance-aware test case generation requires solving two main challenges: providing a good approximation of resource usage with minimal overhead and avoiding detrimental effects on both final coverage and fault detection effectiveness. To tackle these challenges, we conceived a set of performance proxies -- inspired by previous work on performance testing -- that provide a reasonable estimation of the test execution costs (i.e., runtime and memory usage). Thus, we propose an adaptive strategy, called aDynaMOSA, which leverages these proxies by extending DynaMOSA, a state-of-the-art evolutionary algorithm in unit testing. Our empirical study -- involving 110 non-trivial Java classes -- reveals that our adaptive approach generates test suite with statistically significant improvements in runtime (-25%) and heap memory consumption (-15%) compared to DynaMOSA. Additionally, aDynaMOSA has comparable results to DynaMOSA over seven different coverage criteria and similar fault detection effectiveness. Our empirical investigation also highlights that the usage of performance proxies (i.e., without the adaptiveness) is not sufficient to generate more performant test cases without compromising the overall coverage.


title: On Approximate Range Mode and Range Selection

abstract: For any $\epsilon \in (0,1)$, a $(1+\epsilon)$-approximate range mode query asks for the position of an element whose frequency in the query range is at most a factor $(1+\epsilon)$ smaller than the true mode. For this problem, we design an $O(n/\epsilon)$ bit data structure supporting queries in $O(\lg(1/\epsilon))$ time. This is an encoding data structure which does not require access to the input sequence; we prove the space cost is asymptotically optimal for constant $\epsilon$. Our solution improves the previous best result of Greve et al. (Cell Probe Lower Bounds and Approximations for Range Mode, ICALP'10) by reducing the space cost by a factor of $\lg n$ while achieving the same query time. We also design an $O(n)$-word dynamic data structure that answers queries in $O(\lg n /\lg\lg n)$ time and supports insertions and deletions in $O(\lg n)$ time, for any constant $\epsilon \in (0,1)$. This is the first result on dynamic approximate range mode; it can also be used to obtain the first static data structure for approximate 3-sided range mode queries in two dimensions. We also consider approximate range selection. For any $\alpha \in (0,1/2)$, an $\alpha$-approximate range selection query asks for the position of an element whose rank in the query range is in $[k - \alpha s, k + \alpha s]$, where $k$ is a rank given by the query and $s$ is the size of the query range. When $\alpha$ is a constant, we design an $O(n)$-bit encoding data structure that can answer queries in constant time and prove this space cost is asymptotically optimal. The previous best result by Krizanc et al. (Range Mode and Range Median Queries on Lists and Trees, Nordic Journal of Computing, 2005) uses $O(n\lg n)$ bits, or $O(n)$ words, to achieve constant approximation for range median only. Thus we not only improve the space cost, but also provide support for any arbitrary $k$ given at query time.


title: Optical control of synchronous phases in a programmable polariton cell

abstract: We demonstrate deterministic control of the nearest and next-nearest neighbor coupling in the unit cell of a square lattice of microcavity exciton-polariton condensates. We tune the coupling in a continuous and reversible manner by optically imprinting potential barriers of variable height, in the form of spatially localized incoherent exciton reservoirs that modify the particle flow between condensates. By controlling the couplings in a $2\times2$ polariton cluster, we realize ferromagnetic, anti-ferromagnetic and paired ferromagnetic phases. Our approach paves the way towards simulating complex condensed matter phases through precise control of the individual couplings in networks of optical nonlinear oscillators.


title: Partitioning a reflecting stationary set

abstract: We address the question of whether a reflecting stationary set may be partitioned into two or more reflecting stationary subsets, providing various affirmative answers in ZFC. As an application to singular cardinals combinatorics, we infer that it is never the case that there exists a singular cardinal all of whose scales are very good.


title: Systematic benchmarking of planar nitrogen-incorporated ultrananocrystalline diamond field emission electron source: rf conditioning and beam spatio-temporal characteristics

abstract: Planar nitrogen-incorporated ultrananocrystalline diamond, (N)UNCD, has emerged as a unique field emission source attractive for accelerator applications because of its capability to generate high charge beam and handle moderate vacuum conditions. Most importantly, (N)UNCD sources are simple to produce: conventional high aspect ratio isolated emitters are not required to be formed on the surface, and the actual emitter surface roughness is on the order of only 100~nm. Careful reliability assessment of (N)UNCD is required before it may find routine application in accelerator systems. In the present study using an L-band normal conducting single-cell rf gun, a (N)UNCD cathode has been conditioned to $\sim$42~MV/m in a well-controlled manner. It reached a maximum output charge of 15~nC corresponding to an average current of 6~mA during an emission period of 2.5~$\mu$s. Imaging of emission current revealed a large number of isolated emitters (density over 100/cm$^{2}$) distributed on the cathode, which is consistent with previous tests in dc environments. The performance metrics, the emission imaging, and the systematic study of emission properties during rf conditioning in a wide gradient range assert (N)UNCD as an enabling electron source for rf injector designs serving industrial and scientific applications. These studies also improve the fundamental knowledge of the practical conditioning procedure via better understanding of emission mechanisms.


title: Interplay between $\Delta$ Particles and Hyperons in Neutron Stars

abstract: We analyze the effects of including $\Delta(1232)$ isobars in an equation of state (EoS) for cold, $\beta$-stable neutron star matter, employing relativistic nuclear mean field theory. The selected EoS reproduces the properties of nuclear matter and finite nuclei and, in the astrophysical context, allows for the presence of hyperons in neutron stars having masses larger than 2$M_{\odot}$. We find that the composition and structure of neutron stars is critically influenced by the addition of the $\Delta$ isobars, which allows us to constrain their interaction with the meson fields taking into account astrophysical information. Imposing that the EoS is stable and ensures the existence of 2$M_{\odot}$ neutron stars, as well as requiring agreement with data of $\Delta$ excitation in nuclei, we find that, in the absence of other mechanisms stiffening the EoS at high densities, the interaction of the $\Delta$ isobars with the sigma and omega meson fields must be at least 10\% stronger than that of the nucleons. Moreover, the neutron star moment of inertia turns out to be sensitive to the presence of $\Delta$ isobars, whereas the inclusion of $\Delta$ isobars in the EoS allows for smaller stellar radii and for a lower value of the tidal deformability consistent with the analysis of the GW170817 merger event.


title: CraftAssist: A Framework for Dialogue-enabled Interactive Agents

abstract: This paper describes an implementation of a bot assistant in Minecraft, and the tools and platform allowing players to interact with the bot and to record those interactions. The purpose of building such an assistant is to facilitate the study of agents that can complete tasks specified by dialogue, and eventually, to learn from dialogue interactions.


title: Measuring the local non-convexity of real algebraic curves

abstract: The goal of this paper is to measure the non-convexity of compact and smooth connected components of real algebraic plane curves. We study these curves first in a general setting and then in an asymptotic one. In particular, we consider sufficiently small levels of a real bivariate polynomial in a small enough neighbourhood of a strict local minimum at the origin of the real affine plane. We introduce and describe a new combinatorial object, called the Poincare-Reeb graph, whose role is to encode the shape of such curves and to allow us to quantify their non-convexity. Moreover, we prove that in this setting the Poincare-Reeb graph is a plane tree and can be used as a tool to study the asymptotic behaviour of level curves near a strict local minimum. Finally, using the real polar curve, we show that locally the shape of the levels stabilises and that no spiralling phenomena occur near the origin.


title: CityScopeAR: Urban Design and Crowdsourced Engagement Platform

abstract: Processes of urban planning, urban design and architecture are inherently tangible, iterative and collaborative. Nevertheless, the majority of tools in these fields offer virtual environments and single user experience. This paper presents CityScopeAR: a computational-tangible mixed-reality platform designed for collaborative urban design processes. It portrays the evolution of the tool and presents an overview of the history and limitations of notable CAD and TUI platforms. As well, it depicts the development of a distributed networking system between TUIs and CityScopeAR, as a key in design collaboration. It shares the potential advantage of broad and decentralized community-engagement process using such tools. Finally, this paper demonstrates several real-world tests and deployments of CityScopeAR and proposes a path to future integration of AR/MR devices in urban design and public participation.


title: Attitude Control of a Novel Tailsitter: Swiveling Biplane-Quadrotor

abstract: This paper proposes a solution to the attitude tracking problem for a novel quadrotor tailsitter unmanned aerial vehicle called swiveling biplane quadrotor. The proposed vehicle design addresses the lack of yaw control authority in conventional biplane quadrotor tailsitters by proposing a new design wherein two wings with two attached propellers are joined together with a rod through a swivel mechanism. The yaw torque is generated by relative rotation of the thrust vector of each wing. The unique design of this configuration having two rigid bodies interconnected through a rod with zero torsional rigidity makes the vehicle underactuated in the attitude configuration manifold. An output tracking problem is posed which results in a single equivalent rigid body attitude tracking problem with second-order moment dynamics. The proposed controller is uniformly valid for all attitudes and is based on dynamic feedback linearization in a geometric control framework. Almost-global asymptotic stability of the desired equilibrium of the tracking error dynamics is shown. The efficacy of the controller is shown with numerical simulation and flight tests.


title: On Usefulness of the Deep-Learning-Based Bug Localization Models to Practitioners

abstract: Background: Developers spend a significant amount of time and efforts to localize bugs. In the literature, many researchers proposed state-of-the-art bug localization models to help developers localize bugs easily. The practitioners, on the other hand, expect a bug localization tool to meet certain criteria, such as trustworthiness, scalability, and efficiency. The current models are not capable of meeting these criteria, making it harder to adopt these models in practice. Recently, deep-learning-based bug localization models have been proposed in the literature, which show a better performance than the state-of-the-art models. Aim: In this research, we would like to investigate whether deep learning models meet the expectations of practitioners or not. Method: We constructed a Convolution Neural Network and a Simple Logistic model to examine their effectiveness in localizing bugs. We train these models on five open source projects written in Java and compare their performance with the performance of other state-of-the-art models trained on these datasets. Results: Our experiments show that although the deep learning models perform better than classic machine learning models, they meet the adoption criteria set by the practitioners only partially. Conclusions: This work provides evidence that the practitioners should be cautious while using the current state of the art models for production-level use-cases. It also highlights the need for standardization of performance benchmarks to ensure that bug localization models are assessed equitably and realistically.


title: Spectral Analysis of Latent Representations

abstract: We propose a metric, Layer Saturation, defined as the proportion of the number of eigenvalues needed to explain 99% of the variance of the latent representations, for analyzing the learned representations of neural network layers. Saturation is based on spectral analysis and can be computed efficiently, making live analysis of the representations practical during training. We provide an outlook for future applications of this metric by outlining the behaviour of layer saturation in different neural architectures and problems. We further show that saturation is related to the generalization and predictive performance of neural networks.


title: One-loop Holography with Strings in $AdS_4\times\mathbb {CP}^3$

abstract: We compute the one-loop effective action of string configurations embedded in $AdS_4\times\mathbb{CP}^3$ which are dual to $\frac{1}{6}$-BPS latitude Wilson Loops in the ABJM theory. To avoid ambiguities in the string path integral we subtract the $\frac{1}{2}$-BPS case. The one-loop determinants are computed by Fourier-decomposing the two dimensional operators and then using the Gel'fand-Yaglom method. We comment on various aspects related to the regularization procedure, showing the cancellation of a hierarchy of divergences. After taking into account an IR anomaly from a change in topology, we find a precise agreement with the field theory result known from supersymmetric localization.


title: Zermelo's problem: Optimal point-to-point navigation in 2D turbulent flows using Reinforcement Learning

abstract: To find the path that minimizes the time to navigate between two given points in a fluid flow is known as Zermelo's problem. Here, we investigate it by using a Reinforcement Learning (RL) approach for the case of a vessel which has a slip velocity with fixed intensity, Vs , but variable direction and navigating in a 2D turbulent sea. We show that an Actor-Critic RL algorithm is able to find quasi-optimal solutions for both time-independent and chaotically evolving flow configurations. For the frozen case, we also compared the results with strategies obtained analytically from continuous Optimal Navigation (ON) protocols. We show that for our application, ON solutions are unstable for the typical duration of the navigation process, and are therefore not useful in practice. On the other hand, RL solutions are much more robust with respect to small changes in the initial conditions and to external noise, even when V s is much smaller than the maximum flow velocity. Furthermore, we show how the RL approach is able to take advantage of the flow properties in order to reach the target, especially when the steering speed is small.

categories: nlin.CD cs.AI cs.LG cs.SY eess.SY physics.flu-dyn


title: Kernel Mode Decomposition and programmable/interpretable regression networks

abstract: Mode decomposition is a prototypical pattern recognition problem that can be addressed from the (a priori distinct) perspectives of numerical approximation, statistical inference and deep learning. Could its analysis through these combined perspectives be used as a Rosetta stone for deciphering mechanisms at play in deep learning? Motivated by this question we introduce programmable and interpretable regression networks for pattern recognition and address mode decomposition as a prototypical problem. The programming of these networks is achieved by assembling elementary modules decomposing and recomposing kernels and data. These elementary steps are repeated across levels of abstraction and interpreted from the equivalent perspectives of optimal recovery, game theory and Gaussian process regression (GPR). The prototypical mode/kernel decomposition module produces an optimal approximation $(w_1,w_2,\cdots,w_m)$ of an element $(v_1,v_2,\ldots,v_m)$ of a product of Hilbert subspaces of a common Hilbert space from the observation of the sum $v:=v_1+\cdots+v_m$. The prototypical mode/kernel recomposition module performs partial sums of the recovered modes $w_i$ based on the alignment between each recovered mode $w_i$ and the data $v$. We illustrate the proposed framework by programming regression networks approximating the modes $v_i= a_i(t)y_i\big(\theta_i(t)\big)$ of a (possibly noisy) signal $\sum_i v_i$ when the amplitudes $a_i$, instantaneous phases $\theta_i$ and periodic waveforms $y_i$ may all be unknown and show near machine precision recovery under regularity and separation assumptions on the instantaneous amplitudes $a_i$ and frequencies $\dot{\theta}_i$. The structure of some of these networks share intriguing similarities with convolutional neural networks while being interpretable, programmable and amenable to theoretical analysis.


title: Microwave-to-optical conversion using lithium niobate thin-film acoustic resonators

abstract: We demonstrate conversion of up to 4.5 GHz-frequency microwaves to 1500 nm-wavelength light using optomechanical interactions on suspended thin-film lithium niobate. Our method utilizes an interdigital transducer that drives a free-standing 100 $\mu$m-long thin-film acoustic resonator to modulate light travelling in a Mach-Zehnder interferometer or racetrack cavity. Owing to the strong microwave-to-acoustic coupling offered by the transducer in conjunction with the strong photoelastic, piezoelectric, and electro-optic effects of lithium niobate, we achieve a half-wave voltage of $V_\pi$ = 4.6 V and $V_\pi$ = 0.77 V for the Mach-Zehnder interferometer and racetrack resonator, respectively. The acousto-optic racetrack cavity exhibits an optomechancial single-photon coupling strength of 1.1 kHz. Our integrated nanophotonic platform coherently leverages the compelling properties of lithium niobate to achieve microwave-to-optical transduction. To highlight the versatility of our system, we also demonstrate a lossless microwave photonic link, which refers to a 0 dB microwave power transmission over an optical channel.


title: The inner dust shell of Betelgeuse detected by polarimetric aperture-masking interferometry

abstract: Theory surrounding the origin of the dust-laden winds from evolved stars remains mired in controversy. Characterizing the formation loci and the dust distribution within approximately the first stellar radius above the surface is crucial for understanding the physics that underlie the mass-loss phenomenon. By exploiting interferometric polarimetry, we derive the fundamental parameters that govern the dust structure at the wind base of a red supergiant. We present near-infrared aperture-masking observations of Betelgeuse in polarimetric mode obtained with the NACO/SAMPol instrument. We used both parametric models and radiative transfer simulations to predict polarimetric differential visibility data and compared them to SPHERE/ZIMPOL measurements. Using a thin dust shell model, we report the discovery of a dust halo that is located at only 0.5 R$_{\star}$ above the photosphere (i.e. an inner radius of the dust halo of 1.5 R$_{\star}$). By fitting the data under the assumption of Mie scattering, we estimate the grain size and density for various dust species. By extrapolating to the visible wavelengths using radiative transfer simulations, we compare our model with SPHERE/ZIMPOL data and find that models based on dust mixtures that are dominated by forsterite are most favored. Such a close dusty atmosphere has profound implications for the dust formation mechanisms around red supergiants.


title: QPROP with faster calculation of photoelectron spectra

abstract: The calculation of accurate photoelectron spectra (PES) for strong-field laser-atom experiments is a demanding computational task, even in single-active-electron approximation. The QPROP code, published in 2006, has been extended in 2016 in order to provide the possibility to calculate PES using the so-called t-SURFF approach [L. Tao, A. Scrinzi, New J. Phys. 14, 013021 (2012)]. In t-SURFF, the flux through a surface while the laser is on is monitored. Calculating PES from this flux through a surface enclosing a relatively small computational grid is much more efficient than calculating it from the widely spread wavefunction at the end of the laser pulse on a much larger grid. However, the smaller the minimum photoelectron energy of interest is, the more post-propagation after the actual laser pulse is necessary. This drawback of t-SURFF has been overcome by Morales et al. [F. Morales, T. Bredtmann, S. Patchkovskii, J. Phys. B: At. Mol. Opt. Phys. 49, 245001 (2016)] by noticing that the propagation of the wavefunction from the end of the laser pulse to infinity can be performed very efficiently in a single step. In this work, we introduce QPROP 3.0, in which this single-step post-propagation (dubbed i-SURFV) is added. Examples, illustrating the new feature, are discussed. A few other improvements, concerning mainly the parameter files, are also explained.


title: Generalized Lieb-Schultz-Mattis theorem on bosonic symmetry protected topological phases

abstract: We propose and prove a family of generalized Lieb-Schultz-Mattis (LSM) theorems for symmetry protected topological (SPT) phases on boson/spin models in any dimensions. The "conventional" LSM theorem, applicable to e.g. any translation invariant system with an odd number of spin-1/2 particles per unit cell, forbids a symmetric short-range-entangled ground state in such a system. Here we focus on systems with no LSM anomaly, where global/crystalline symmetries and fractional spins within the unit cell ensure that any symmetric SRE ground state must be a nontrivial SPT phase with anomalous boundary excitations. Depending on models, they can be either strong or "higher-order" crystalline SPT phases, characterized by nontrivial surface/hinge/corner states. Furthermore, given the symmetry group and the spatial assignment of fractional spins, we are able to determine all possible SPT phases for a symmetric ground state, using the real space construction for SPT phases based on the spectral sequence of cohomology theory. We provide examples in one, two and three spatial dimensions, and discuss possible physical realization of these SPT phases based on condensation of topological excitations in fractionalized phases.


title: A refined Brill-Noether theory over Hurwitz spaces

abstract: Let $f\colon C \rightarrow \mathbb{P}^1$ be a degree $k$ genus $g$ cover. The stratification of line bundles $L \in \mathrm{Pic}^d(C)$ by the splitting type of $f_*L$ is a refinement of the stratification by Brill-Noether loci $W^r_d(C)$. We prove that for general degree $k$ covers, these strata are smooth of the expected dimension. In particular, this determines the dimensions of all irreducible components of $W^r_d(C)$ for a general $k$-gonal curve (there are often components of different dimensions), extending results of Pflueger and Jensen-Ranganathan. The results here apply over any algebraically closed field.


title: A Wearable Medical Sensor for Provisional Healthcare

abstract: Thispaper presents the design and realization of a context-aware wireless health monitoring system for recording the heartbeat (HR) and respiration (RR) rate based on an indirect measurement approach. The system consists of a contact-less medical sensor as well as a communication infrastructure for handling the transmission and reception of the measured results. The contact-less sensor includes a highly sensitive tri-axial accelerometer, an accurate temperature and air pressure sensor that enable one to inspect patients' health condition by continuously monitoring of two critical signs related to the cardiorespiratory system. The developed system can also be utilized in performing a number of long-term inspection on the heart and lungs while measuring the HR and RR values in addition to calculating the HR and RR ratio, which is denoted by HRR. The obtained results show the potential of the developed system for versatile monitoring applications applied to telemedicine


title: Conformal Three-Point Correlation Functions from the Operator Product Expansion

abstract: We show how to construct embedding space three-point functions for operators in arbitrary Lorentz representations by employing the formalism developed in arXiv:1905.00036 and arXiv:1905.00434. We study tensor structures that intertwine the operators with the derivatives in the OPE and examine properties of OPE coefficients under permutations of operators. Several examples are worked out in detail. We point out that the group theoretic objects used in this work can be applied directly to construct three-point functions without any reference to the OPE.


title: Learning sparsity in reservoir computing through a novel bio-inspired algorithm

abstract: The mushroom body is the key network for the representation of learned olfactory stimuli in Drosophila and insects. The sparse activity of Kenyon cells, the principal neurons in the mushroom body, plays a key role in the learned classification of different odours. In the specific case of the fruit fly, the sparseness of the network is enforced by an inhibitory feedback neuron called APL, and by an intrinsic high firing threshold of the Kenyon cells. In this work we took inspiration from the fruit fly brain to formulate a novel machine learning algorithm that is able to optimize the sparsity level of a reservoir by changing the firing thresholds of the nodes. The sparsity is only applied on the readout layer so as not to change the timescales of the reservoir and to allow the derivation of a one-layer update rule for the firing thresholds. The proposed algorithm is a combination of learning a neuron-specific sparsity threshold via gradient descent and a global sparsity threshold via a Markov chain Monte Carlo method. The proposed model outperforms the standard gradient descent, which is limited to the readout weights of the reservoir, on two example tasks. It demonstrates how the learnt sparse representation can lead to better classification performance, memorization ability and convergence time.


title: Joule meets van der Waals: Mechanical dissipation via image potential states on a topological insulator surface

abstract: Dissipation mechanisms are experimentally studied on topological insulator surfaces of Bi2Te3, where common Joule dissipation was observed to be suppressed due to topologically protected surface states. Thus, a novel type of dissipation mechanism is observed by pendulum AFM, which is related to single electron tunneling resonances into image potential states that are slightly above the Bi2Te3 surface. The application of a magnetic field leads to the break down of the topological protection of the surface states and restores the expected Joule dissipation process. Nanomechanical energy dissipation experienced by the cantilever of pendulum AFM provides a novel source of information on the dissipative nature of the quantum-tunneling phenomena on the topological insulator surface.


title: Probing Image Potential States on Topological Semimetal Antimony Surface

abstract: A point charge near the surface of a topological insulator (TI) with broken time-reversal symmetry is predicted to generate an image magnetic charge in addition to an image electric charge. We use scanning tunneling spectroscopy to study the image potential states (IPS) of the topological semimetal Sb(111) surface. We observe five IPS with discrete energy levels that are well described by a one-dimensional model. The spatial variation of the IPS energies and lifetimes near surface step edges shows the first local signature of resonant interband scattering between IPS, which suggests that image charges too may interact. Our work motivates the exploration of the TI surface geometry necessary to realize and manipulate a magnetic charge.


title: A scaling hypothesis for matrix product states

abstract: We revisit the question of describing critical spin systems and field theories using matrix product states, and formulate a scaling hypothesis in terms of operators, eigenvalues of the transfer matrix, and lattice spacing in the case of field theories. Critical exponents and central charge are determined by optimizing the exponents such as to obtain a data collapse. We benchmark this method by studying critical Ising and Potts models, where we also obtain a scaling ansatz for the correlation length and entanglement entropy. The formulation of those scaling functions turns out to be crucial for studying critical quantum field theories on the lattice. For the case of $\lambda\phi^4$ with mass $\mu^2$ and lattice spacing $a$, we demonstrate a double data collapse for the correlation length $ \delta \xi(\mu,\lambda,D)=\tilde{\xi} \left((\alpha-\alpha_c)(\delta/a)^{-1/\nu}\right)$ with $D$ the bond dimension, $\delta$ the gap between eigenvalues of the transfer matrix, and $\alpha_c=\mu_R^2/\lambda$ the parameter which fixes the critical quantum field theory.

categories: cond-mat.stat-mech hep-lat hep-th quant-ph


title: Stimulated emission depletion microscopy with diamond silicon-vacancy centers

abstract: The spatial resolution and fluorescence signal amplitude in stimulated emission depletion (STED) microscopy is limited by the photostability of available fluorophores. Here, we show that negatively-charged silicon vacancy (SiV) centers in diamond are promising fluorophores for STED microscopy, owing to their photostable, near-infrared emission and favorable photophysical properties. A home-built pulsed STED microscope was used to image shallow implanted SiV centers in bulk diamond at room temperature. The SiV stimulated emission cross section for 765-800 nm light is found to be (4.0 +/- 0.3) x 10^(-17) cm^2, which is approximately 2-4 times larger than that of the negatively-charged diamond nitrogen vacancy center and approaches that of commonly-used organic dye molecules. We performed STED microscopy on isolated SiV centers and observed a lateral full-width-at-half-maximum spot size of 89 +/- 2 nm, limited by the low available STED laser pulse energy (0.4 nJ). For a pulse energy of 5 nJ, the resolution is expected to be ~20 nm. We show that the present microscope can resolve SiV centers separated by <150 nm that cannot be resolved by confocal microscopy.

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


title: A Detection of CMB-Cluster Lensing using Polarization Data from SPTpol

abstract: We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization gradient. Using data from the SPTpol 500 deg$^{2}$ survey at the locations of roughly 18,000 clusters with richness $\lambda \ge 10$ from the Dark Energy Survey (DES) Year-3 full galaxy cluster catalog, we detect lensing at $4.8\sigma$. The mean stacked mass of the selected sample is found to be $(1.43 \pm 0.4)\ \times 10^{14}\ {\rm M_{\odot}}$ which is in good agreement with optical weak lensing based estimates using DES data and CMB-lensing based estimates using SPTpol temperature data. This measurement is a key first step for cluster cosmology with future low-noise CMB surveys, like CMB-S4, for which CMB polarization will be the primary channel for cluster lensing measurements.


title: Coherence manipulation with dephasing-covariant operations

abstract: We characterize the operational capabilities of quantum channels which can neither create nor detect quantum coherence vis-\`a-vis efficiently manipulating coherence as a resource. We study the class of dephasing-covariant operations (DIO), unable to detect the coherence of any input state, as well as introduce an operationally-motivated class of channels $\rho$-DIO which is tailored to a specific input state. We first show that pure-state transformations under DIO are completely governed by majorization, establishing necessary and sufficient conditions for such transformations and adding to the list of operational paradigms where majorization plays a central role. We then show that $\rho$-DIO are strictly more powerful: although they cannot detect the coherence of the input state $\rho$, the operations $\rho$-DIO can distill more coherence than DIO. However, the advantage disappears in the task of coherence dilution as well as generally in the asymptotic limit, where both sets of operations achieve the same rates in all transformations.


title: Parallel Algorithms for Butterfly Computations

abstract: Butterflies are the smallest non-trivial subgraph in bipartite graphs, and therefore having efficient computations for analyzing them is crucial to improving the quality of certain applications on bipartite graphs. In this paper, we design a framework called ParButterfly that contains new parallel algorithms for the following problems on processing butterflies: global counting, per-vertex counting, per-edge counting, tip decomposition (vertex peeling), and wing decomposition (edge peeling). The main component of these algorithms is aggregating wedges incident on subsets of vertices, and our framework supports different methods for wedge aggregation, including sorting, hashing, histogramming, and batching. In addition, ParButterfly supports different ways of ranking the vertices to speed up counting, including side ordering, approximate and exact degree ordering, and approximate and exact complement coreness ordering. For counting, ParButterfly also supports both exact computation as well as approximate computation via graph sparsification. We prove strong theoretical guarantees on the work and span of the algorithms in ParButterfly. We perform a comprehensive evaluation of all of the algorithms in ParButterfly on a collection of real-world bipartite graphs using a 48-core machine. Our counting algorithms obtain significant parallel speedup, outperforming the fastest sequential algorithms by up to 13.6x with a self-relative speedup of up to 38.5x. Compared to general subgraph counting solutions, we are orders of magnitude faster. Our peeling algorithms achieve self-relative speedups of up to 10.7x and outperform the fastest sequential baseline by up to several orders of magnitude.


title: The inverse problem for arboreal Galois representations of index two

abstract: This paper introduces a systematic approach towards the inverse problem for arboreal Galois representations of finite index attached to quadratic polynomials. Let $F$ be a field of characteristic $\neq 2$, $f\in F[x]$ be monic and quadratic and $\rho_f$ be the arboreal Galois representation associated to $f$, taking values in the group $\Omega_{\infty}$ of automorphisms of the infinite binary tree. We give a complete description of the maximal closed subgroups of each closed subgroup of index at most two of $\Omega_{\infty}$ in terms of linear relations modulo squares among certain universal functions evaluated in elements of the critical orbit of $f$. We use such description in order to derive necessary and sufficient criteria for the image of $\rho_f$ to be a given subgroup of index two of $\Omega_\infty$. These depend exclusively on the arithmetic of the critical orbit of $f$. Afterwards, we prove that if $\phi=x^2+t\in\mathbb Q(t)[x]$, then there exist exactly five distinct subgroups of index two of $\Omega_{\infty}$ that can appear as images of $\rho_{\phi_{t_0}}$ for infinitely many $t_0\in\mathbb Q$, where $\phi_{t_0}$ is the specialized polynomial. We show that two of them appear infinitely often, and if Vojta's conjecture over $\mathbb Q$ holds true, then so do the remaining ones. Finally, we give an explicit description of the derived series of each subgroup of index two. Using this, we introduce a sequence of combinatorial invariants for subgroups of index two of $\Omega_\infty$. With a delicate use of these invariants we are able to establish that such subgroups are pairwise non-isomorphic as topological groups, a result of independent interest. This implies, in particular, that the five aforementioned groups are pairwise distinct topological groups, and therefore yield five genuinely different instances of the infinite inverse Galois problem over $\mathbb Q$.


title: Effective Lagrangian for Nambu-Goldstone modes in nonequilibrium open systems

abstract: We develop the effective field theory of diffusive Nambu-Goldstone (NG) modes associated with spontaneous internal symmetry breaking taking place in nonequilibrium open systems. The effective Lagrangian describing semi-classical dynamics of the NG modes is derived and matching conditions for low-energy coefficients are also investigated. Due to new terms peculiar to open systems, the associated NG modes show diffusive gapless behaviors in contrast to the propagating NG mode in closed systems. We demonstrate two typical situations relevant to the condensed matter physics and high-energy physics, where diffusive type-A or type-B NG modes appear.

categories: hep-th cond-mat.quant-gas cond-mat.stat-mech hep-ph nucl-th


title: Lookahead Optimizer: k steps forward, 1 step back

abstract: The vast majority of successful deep neural networks are trained using variants of stochastic gradient descent (SGD) algorithms. Recent attempts to improve SGD can be broadly categorized into two approaches: (1) adaptive learning rate schemes, such as AdaGrad and Adam, and (2) accelerated schemes, such as heavy-ball and Nesterov momentum. In this paper, we propose a new optimization algorithm, Lookahead, that is orthogonal to these previous approaches and iteratively updates two sets of weights. Intuitively, the algorithm chooses a search direction by looking ahead at the sequence of fast weights generated by another optimizer. We show that Lookahead improves the learning stability and lowers the variance of its inner optimizer with negligible computation and memory cost. We empirically demonstrate Lookahead can significantly improve the performance of SGD and Adam, even with their default hyperparameter settings on ImageNet, CIFAR-10/100, neural machine translation, and Penn Treebank.


title: Distributions.jl: Definition and Modeling of Probability Distributions in the JuliaStats Ecosystem

abstract: Random variables and their distributions are a central part in many areas of statistical methods. The Distributions.jl package provides Julia users and developers tools for working with probability distributions, leveraging Julia features for their intuitive and flexible manipulation, while remaining highly efficient through zero-cost abstractions.


title: Medical Imaging with Deep Learning: MIDL 2019 -- Extended Abstract Track

abstract: This compendium gathers all the accepted extended abstracts from the Second International Conference on Medical Imaging with Deep Learning (MIDL 2019), held in London, UK, 8-10 July 2019. Note that only accepted extended abstracts are listed here, the Proceedings of the MIDL 2019 Full Paper Track are published as Volume 102 of the Proceedings of Machine Learning Research (PMLR) http://proceedings.mlr.press/v102/.


title: Generic "Unnecessary" Quantum Critical Points with Minimal Degrees of Freedom

abstract: We explore generic "unnecessary" quantum critical points with minimal degrees of freedom. These quantum critical points can be avoided with strong enough symmetry-allowed deformations of the Hamiltonian, but these deformations are irrelevant perturbations below certain threshold at the quantum critical point. These quantum critical points are hence unnecessary, but also unfine-tuned (generic). The previously known examples of such generic unnecessary quantum critical points involve at least eight Dirac fermions in both two and three spatial dimensions. In this work we seek for examples of generic unnecessary quantum critical points with minimal degrees of freedom. In particular, in three dimensional space, we identify two examples of such generic unnecessary quantum critical points. The first example occurs in a $3d$ interacting topological insulator, and it is described by $two$ $(3+1)d$ massless Dirac fermions in the infrared limit; the second example occurs in a $3d$ topological superconductor, and it is formally described only $one$ $(3+1)d$ massless Dirac fermion.


title: Black hole-neutron star mergers from triples II: the role of metallicity and spin-orbit misalignment

abstract: Observations of black hole-neutron star (BH-NS) mergers via gravitational waves (GWs) are of great interest for their electromagnetic counterparts, such as short gamma-ray bursts, and could provide crucial information on the nature of BHs and the NS crust and magnetosphere. While no event has been confirmed, a recent possible detection of a BH-NS merger event by the LIGO-Virgo collaboration has attracted a lot of attention to these sources. In this second paper of the series, we follow-up our study of the dynamical evolution of triples comprised of an inner BH-NS binary. In particular, we examine how the progenitor metallicity affects the characteristics of the BH-NS mergers in triples. We determine the distributions of masses, orbital parameters and merger times, as a function of the progenitor metallicity and initial triple orbital distributions, and show that the typical eccentricity in the LIGO band is $\sim 10^{-2}-10^{-1}$. We derive a merger rate range of $\Gamma_\mathrm{BH-NS}=1.9\times 10^{-4}-22 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, consistent the LIGO-Virgo upper limit. Finally, we study the expected spin-orbit misalignments of merging BH-NS binaries from this channel, and find that typically the effective spin distribution is peaked at $\chi_{\rm eff}\sim 0$ with significant tails.


title: Logical Segmentation of Source Code

abstract: Many software analysis methods have come to rely on machine learning approaches. Code segmentation - the process of decomposing source code into meaningful blocks - can augment these methods by featurizing code, reducing noise, and limiting the problem space. Traditionally, code segmentation has been done using syntactic cues; current approaches do not intentionally capture logical content. We develop a novel deep learning approach to generate logical code segments regardless of the language or syntactic correctness of the code. Due to the lack of logically segmented source code, we introduce a unique data set construction technique to approximate ground truth for logically segmented code. Logical code segmentation can improve tasks such as automatically commenting code, detecting software vulnerabilities, repairing bugs, labeling code functionality, and synthesizing new code.


title: On special matrices related to Cauchy and Toeplitz matrices

abstract: In this paper, we are going to calculate the determinant of a certain type of square matrices, which are related to the well-known Cauchy and Toeplitz matrices. Then, we will use the results to determine the rank of special non-square matrices.


title: Detection of extensive cosmic ray air showers by measuring radio emission

abstract: So far, cosmic ray air showers have been detected using scintillation counter arrays on the ground widely. And also air Cherenkov detection method, which is limited its observation period in moonless nights, has been adopted. The detection method of radio emission from cosmic ray air showers is not new, but rather old method. Radio emission from cosmic ray air showers has not been detected with the method of self-trigger system. If the detection method of radio emission were available, there is no limit of the observation like Cherenkov counter. The developments of high radio frequency (RF) technology might make the detection of radio emission from extensive air showers possible. Antennas calibrated in laboratory are available. With those antennas, we can directly obtain the absolute intensity of radio frequency from cosmic ray air showers. Main issue to detect radio emission is that signal level is quite low. The thermal noise particularly causes background noise source. Those issues and detection method are discussed. We also describe an antenna detection method to generate self-trigger signal and a method to identify air shower events from natural and artificial noises.


title: Null Decomposition of Unicyclic Graphs

abstract: In this work we obtain basis for the null space of unicyclic graphs. We extend the null decomposition of trees from [11] for unicyclic graphs. As an application, we obtain closed formulas for the independence and matching numbers of unicyclic graphs just using the support of the graph.


title: Real-time UAV Complex Missions Leveraging Self-Adaptive Controller with Elastic Structure

abstract: The expectation of unmanned air vehicles (UAVs) pushes the operation environment to narrow spaces, where the systems may fly very close to an object and perform an interaction. This phase brings the variation in UAV dynamics: thrust and drag coefficient of the propellers might change under different proximity. At the same time, UAVs may need to operate under external disturbances to follow time-based trajectories. Under these challenging conditions, a standard controller approach may not handle all missions with a fixed structure, where there may be a need to adjust its parameters for each different case. With these motivations, practical implementation and evaluation of an autonomous controller applied to a quadrotor UAV are proposed in this work. A self-adaptive controller based on a composite control scheme where a combination of sliding mode control (SMC) and evolving neuro-fuzzy control is used. The parameter vector of the neuro-fuzzy controller is updated adaptively based on the sliding surface of the SMC. The autonomous controller possesses a new elastic structure, where the number of fuzzy rules keeps growing or get pruned based on bias and variance balance. The interaction of the UAV is experimentally evaluated in real time considering the ground effect, ceiling effect and flight through a strong fan-generated wind while following time-based trajectories.


title: Bishop-Phelps-Bollob\'as property for positive operators when the domain is $L_\infty $

abstract: We prove that the class of positive operators from $L_\infty (\mu)$ to $Y$ has the Bishop-Phelps-Bollob\'as property for any positive measure $\mu$, whenever $Y$ is a uniformly monotone Banach lattice with a weak unit. The same result also holds for the pair $(c_0, Y)$ for any uniformly monotone Banach lattice $Y.$ Further we show that these results are optimal in case that $Y$ is strictly monotone.


title: A lattice model for active--passive pedestrian dynamics: a quest for drafting effects

abstract: We study the pedestrian escape from an obscure corridor using a lattice gas model with two species of particles. One species, called passive, performs a symmetric random walk on the lattice, whereas the second species, called active, is subject to a drift guiding the particles towards the exit. The drift mimics the awareness of some pedestrians of the geometry of the corridor and of the location of the exit. We provide numerical evidence that, in spite of the hard core interaction between particles -- namely, there can be at most one particle of any species per site, -- adding a fraction of active particles in the system enhances the evacuation rate of all particles from the corridor. A similar effect is also observed when looking at the outgoing particle flux, when the system is in contact with an external particle reservoir that induces the onset of a steady state. We interpret this phenomenon as a discrete space counterpart of the drafting effect typically observed in a continuum set--up as the aerodynamic drag experienced by pelotons of competing cyclists.


title: Twisted Dirac Operators and the noncommutative residue for manifolds with boundary II

abstract: In this paper, we establish two kinds of Kastler-Kalau-Walze type theorems for Dirac operators and signature operators twisted by a vector bundle with a non-unitary connection on six-dimensional manifolds with boundary.


title: Effects of changing population or density on urban carbon dioxide emissions

abstract: The question of whether urbanization contributes to increasing carbon dioxide emissions has been mainly investigated via scaling relationships with population or population density. However, these approaches overlook the correlations between population and area, and ignore possible interactions between these quantities. Here, we propose a generalized framework that simultaneously considers the effects of population and area along with possible interactions between these urban metrics. Our results significantly improve the description of emissions and reveal the coupled role between population and density on emissions. These models show that variations in emissions associated with proportionate changes in population or density may not only depend on the magnitude of these changes but also on the initial values of these quantities. For US areas, the larger the city, the higher is the impact of changing its population or density on its emissions; but population changes always have a greater effect on emissions than population density.


title: Scale-dependent measure of network centrality from diffusion dynamics

abstract: Classic measures of graph centrality capture distinct aspects of node importance, from the local (e.g., degree) to the global (e.g., closeness). Here we exploit the connection between diffusion and geometry to introduce a multiscale centrality measure. A node is defined to be central if it breaks the metricity of the diffusion as a consequence of the effective boundaries and inhomogeneities in the graph. Our measure is naturally multiscale, as it is computed relative to graph neighbourhoods within the varying time horizon of the diffusion. We find that the centrality of nodes can differ widely at different scales. In particular, our measure correlates with degree (i.e., hubs) at small scales and with closeness (i.e., bridges) at large scales, and also reveals the existence of multi-centric structures in complex networks. By examining centrality across scales, our measure thus provides an evaluation of node importance relative to local and global processes on the network.


title: On the disc reflection spectroscopy of NS LMXB Serpens~X-1: analysis of a recent NuSTAR observation

abstract: We present \nustar{} observation of the atoll type neutron star (NS) low-mass X-ray binary (LMXB) Serpens~X-1 (Ser~X-1) performed on 17 February 2018. We observed Ser~X-1 in a soft X-ray spectral state with 3-79 keV luminosity of $L_\text{X}\sim0.4\times 10^{38}$ erg s$^{-1}$ ($\sim 23\%$ of the Eddington luminosity), assuming a distance of 7.7 kpc. A positive correlation between intensity and hardness ratio suggests that the source was in the banana branch during this observation. The broadband 3-30 keV NuSTAR energy spectrum can be well described either by a three-component continuum model consisting of a disk blackbody, a single temperature blackbody and a power-law or by a two-component continuum model consisting of a disk blackbody and a Comptonization component. A broad iron line $\sim 5-8$ keV and the Compton back-scattering hump peaking at $\sim$10-20 keV band are clearly detected in the X-ray spectrum. These features are best interpreted by a self-consistent relativistic reflection model. Fits with relativistically blurred disc reflection model suggests that the inner disc radius $R_{in}$ is truncated prior to the ISCO at $(1.9-2.5)\;R_{ISCO}$ ($\simeq11.4-15\,R_{g}\: \text{or}\: 26-34$ km) and the accretion disc is viewed at an low inclination of $i\simeq16^\circ-20^\circ$. The disc is likely to be truncated either by a boundary layer or by the magnetosphere. Based on the measured flux and the mass accretion rate, the maximum radial extension for the boundary layer is estimated to be $\sim6.4\:R_{g}$ from the NS surface. The truncated inner disc in association with pressure from a magnetic field sets an upper limit of $B\leq1.9\times10^{9}$ G.


title: A new class of polynomials from to the spectrum of a graph, and its application to bound the $k$-independence number

abstract: The $k$-independence number of a graph is the maximum size of a set of vertices at pairwise distance greater than $k$. A graph is called $k$-partially walk-regular if the number of closed walks of a given length $l\le k$, rooted at a vertex $v$, only depends on $l$. In particular, a distance-regular graph is also $k$-partially walk-regular for any $k$. In this note, we introduce a new family of polynomials obtained from the spectrum of a graph. These polynomials, together with the interlacing technique, allow us to give tight spectral bounds on the $k$-independence number of a $k$-partially walk-regular graph. Together with some examples where the bounds are tight, we also show that the odd graph $O_{\ell}$ with $\ell$ odd has no $1$-perfect code.


title: Extent of occurrence reconstruction using a new data-driven support estimator

abstract: Given a random sample of points from some unknown distribution, we propose a new data-driven method for estimating its probability support S. Under the mild assumption that S is r-convex, the smallest r-convex set which contains the sample points is the natural estimator. The main problem for using this estimator in practice is that r is an unknown geometric characteristic of the set S. A stochastic algorithm is proposed for determining an optimal estimate of r from the data under mild regularity assumptions on the density function. The resulting data-driven reconstruction of S attains the same convergence rates as the convex hull for estimating convex sets, but under a much more flexible smoothness shape condition. The new support estimator will be used for reconstructing the extent of occurrence of an assemblage of invasive plant species in the Azores archipelago.


title: Discovery and spectroscopy of dark matter and dark sectors with microscopic black holes at next generation colliders

abstract: If the length scale of possible extra dimensions is large enough, the effective Planck scale is lowered such that microscopic black holes could be produced in collisions of high-energy particles at colliders. These black holes evaporate through Hawking radiation of a handful of energetic particles drawn from the set of all kinematically and thermally allowed degrees of freedom, including dark matter. Here, we perform the first numerical black hole spectroscopic study of the dark sector. We find that if the next generation of colliders can produce microscopic black holes, then missing momentum signatures can reveal the existence of any new light ($\lesssim 10$ TeV) particle, regardless of the strength of its coupling to the Standard Model, even if there exists no such non-gravitational coupling at all.


title: VOAs and rank-two instanton SCFTs

abstract: We analyze the N=2 superconformal field theories that arise when a pair of D3-branes probe an F-theory singularity from the perspective of the associated vertex operator algebra. We identify these vertex operator algebras for all cases; we find that they have a completely uniform description, parameterized by the dual Coxeter number of the corresponding global symmetry group. We further present free field realizations for these algebras in the style of recent work by three of the authors. These realizations transparently reflect the algebraic structure of the Higgs branches of these theories. We find fourth-order linear modular differential equations for the vacuum characters/Schur indices of these theories, which are again uniform across the full family of theories and parameterized by the dual Coxeter number. We comment briefly on expectations for the still higher-rank cases.


title: Systematic classification of two loop d = 4 Dirac neutrino mass models and the Diracness-dark matter stability connection

abstract: We provide a complete systematic classification of all two-loop realizations of the dimension four operator for Dirac neutrino masses. Our classification is multi-layered, starting first with a classification in terms of all possible distinct two-loop topologies. Then we discuss the possible diagrams for each topology. Model-diagrams originating from each diagram are then considered. The criterion for genuineness is also defined and discussed at length. Finally, as examples, we construct two explicit models which also serve to highlight the intimate connection between the Dirac nature of neutrinos and dark matter stability.


title: Determination of Intrinsic Effective Fields and Microwave Polarizations by High-Resolution Spectroscopy of Single NV Center Spins

abstract: We present high-resolution optically detected magnetic resonance (ODMR) spectroscopy on single nitrogen-vacancy (NV) center spins in diamond at and around zero magnetic field. The experimentally observed transitions depend sensitively on the interplay between the microwave (MW) probing field and the local intrinsic effective field comprising strain and electric fields, which act on the NV spin. Based on a theoretical model of the magnetic dipole transitions and the MW driving field, we extract both the strength and the direction of the transverse component of the effective field. Our results reveal that for the diamond crystal under study, strain is the dominant contribution to the effective field. Our experiments further yield a method for MW polarization analysis in a tunable, linear basis, which we demonstrate on a single NV spin. Our results are of importance to low-field quantum sensing applications using NV spins and form a relevant addition to the ever-growing toolset of spin-based quantum sensing.


title: BAT AGN Spectroscopic Survey -- XVII: The Parsec-scale Jet Properties of the Ultra Hard X-ray Selected Local AGN

abstract: We have performed a very long baseline interferometry (VLBI) survey of local (z < 0.05) ultra hard X-ray (14-195 keV) selected active galactic nuclei (AGN) from the Swift Burst Alert Telescope (BAT) using KVN, KaVA, and VLBA. We first executed fringe surveys of 142 BAT-detected AGN at 15 or 22 GHz. Based on the fringe surveys and archival data, we find 10/279 nearby AGN (~4%) VLBI have 22 GHz flux above 30 mJy. This implies that the X-ray AGN with a bright nuclear jet are not common. Among these 10 radio-bright AGN, we obtained 22 GHz VLBI imaging data of our own for four targets and reprocessed archival data for six targets. We find that, although our 10 AGN observed with VLBI span a wide range of pc-scale morphological types, they lie on a tight linear relation between accretion luminosity and nuclear jet luminosity. Our result suggests that a powerful nuclear radio jet correlates with the accretion disc luminosity. We also probed the fundamental plane of black hole activity at VLBI scales (e.g., few milli-arcsecond). The jet luminosity and size distribution among our sample roughly fit into the proposed AGN evolutionary scenario, finding powerful jets after the blow-out phase based on the Eddington ratio (\lambda_{Edd})-hydrogen column density (N_{H}) relation. In addition, we find some hints of gas inflow or galaxy-galaxy merger in the majority of our sample. This implies that gas supply via tidal interactions in galactic scale may help the central AGN to launch a powerful parsec-scale jet.


title: CO and Fine-Structure Lines Reveal Low Metallicity in a Stellar-Mass-Rich Galaxy at z ~ 1?

abstract: We present detections of the CO(4-3) and [C I] 609 $\mu$m spectral lines, as well as the dust continuum at 480.5 GHz (rest-frame), in 3C 368, a Fanaroff-Riley class II (FR-II) galaxy at redshift (z) 1.131. 3C 368 has a large stellar mass, ~ 3.6 x 10$^{11}$ M$_\odot$, and is undergoing an episode of vigorous star formation, at a rate of ~ 350 M$_\odot$/yr, and active galactic nucleus (AGN) activity, with radio-emitting lobes extended over ~ 73 kpc. Our observations allow us to inventory the molecular-gas reservoirs in 3C 368 by applying three independent methods: (1) using the CO(4-3)-line luminosity, excitation state of the gas, and an $\alpha_{CO}$ conversion factor, (2) scaling from the [C I]-line luminosity, and (3) adopting a gas-to-dust conversion factor. We also present gas-phase metallicity estimates in this source, both using far-infrared (FIR) fine-structure lines together with radio free-free continuum emission and independently employing the optical [O III] 5007 A and [O II] 3727 A lines (R$_{23}$ method). Both methods agree on a sub-solar gas-phase metallicity of ~ 0.3 Z$_\odot$. Intriguingly, comparing the molecular-gas mass estimated using this sub-solar metallicity, M$_{gas}$ ~ 6.4 x 10$^{10}$ M$_\odot$, to dust-mass estimates from multi-component spectral energy distribution (SED) modeling, M$_{dust}$ ~ 1.4 x 10$^8$ M$_\odot$, yields a gas-to-dust ratio within ~ 15% of the accepted value for a metallicity of 0.3 Z$_\odot$. The derived gas-mass puts 3C 368 on par with other galaxies at z ~ 1 in terms of specific star-formation rate and gas fraction. However, it does not explain how a galaxy can amass such a large stellar population while maintaining such a low gas-phase metallicity. Perhaps 3C 368 has recently undergone a merger, accreting pristine molecular gas from an external source.