title: On Semi-parametric Bernstein-von Mises Theorems for BART

abstract: Few methods in Bayesian non-parametric statistics/ machine learning have received as much attention as Bayesian Additive Regression Trees (BART). While BART is now routinely performed for prediction tasks, its theoretical properties began to be understood only very recently. In this work, we continue the theoretical investigation of BART initiated by Rockova and van der Pas (2017). In particular, we study the Bernstein-von Mises (BvM) phenomenon (i.e. asymptotic normality) for smooth linear functionals of the regression surface within the framework of non-parametric regression with fixed covariates. As with other adaptive priors, the BvM phenomenon may fail when the regularities of the functional and the truth are not compatible. To overcome the curse of adaptivity under hierarchical priors, we induce a self-similarity assumption to ensure convergence towards a single Gaussian distribution as opposed to a Gaussian mixture. Similar qualitative restrictions on the functional parameter are known to be necessary for adaptive inference. Many machine learning methods lack coherent probabilistic mechanisms for gauging uncertainty. BART readily provides such quantification via posterior credible sets. The BvM theorem implies that the credible sets are also confidence regions with the same asymptotic coverage. This paper presents the first asymptotic normality result for BART priors, providing another piece of evidence that BART is a valid tool from a frequentist point of view.


title: On the canonical ideal of the Ehrhart ring of the chain polytope of a poset

abstract: Let P be a poset, O(P) the order polytope of P and C(P) the chain polytope of P. In this paper, we study the canonical ideal of the Ehrhart ring K[C(P)] of C(P) over a field K and characterize the level (resp. anticanonical level) property of K[C(P)] by a combinatorial structure of P. In particular, we show that if K[C(P)] is level (resp. anticanonical level), then so is K[O(P)]. We exhibit examples which show the converse does not hold. Moreover, we show that the symbolic powers of the canonical ideal of K[C(P)] are identical with ordinary ones and degrees of the generators of the canonical and anticanonical ideals are consecutive integers.


title: Notes on the $D=11$ pure spinor superparticle

abstract: The $D=11$ pure spinor superparticle has been shown to describe linearized $D=11$ supergravity in a manifestly covariant way. A number of authors have proposed that its correlation functions be used to compute amplitudes. The use of the scalar structure of the eleven-dimensional pure spinor top cohomology introduces a natural measure for computing such correlation functions. This prescription requires the construction of ghost number one and zero vertex operators. In these notes, we construct explicitly a ghost number one vertex operator but show the incompatibiliy of a ghost number zero vertex operator satisfying a standard descent equation for $D=11$ supergravity.


title: Search for Ultra-High-Energy Neutrinos with the Telescope Array Surface Detector

abstract: We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function.


title: Modelling radiation damage to pixel sensors in the ATLAS detector

abstract: Silicon pixel detectors are at the core of the current and planned upgrade of the ATLAS experiment at the LHC. Given their close proximity to the interaction point, these detectors will be exposed to an unprecedented amount of radiation over their lifetime. The current pixel detector will receive damage from non-ionizing radiation in excess of $10^{15}$ 1 MeV ${n}_{eq}/{cm}^2$, while the pixel detector designed for the high-luminosity LHC must cope with an order of magnitude larger fluence. This paper presents a digitization model incorporating effects of radiation damage to the pixel sensors. The model is described in detail and predictions for the charge collection efficiency and Lorentz angle are compared with collision data collected between 2015 and 2017 ($\leq 10^{15}$ 1 MeV ${n}_{eq}/{cm}^2$).


title: NGC 300 ULX1: spin evolution, super-Eddington accretion and outflows

abstract: NGC300 ULX1 is an ultra-luminous X-ray pulsar, showing an unprecedented spin evolution, from about 126 s to less than 20 s in only 4 years, consistent with steady mass accretion rate. Following its discovery we have been monitoring the system with Swift/XRT and NICER to further study its properties. We found that even though the observed flux of the system dropped by a factor of $\gtrsim$20, the spin-up rate remained almost constant. A possible explanation is that the decrease in the observed flux is a result of increased absorption of obscuring material due to outflows or a precessing accretion disk.


title: New composite supernova remnant toward HESS J1844-030?

abstract: AIMS: HESS J1844-030 is a newly confirmed TeV source in the direction of the X-ray pulsar wind nebula (PWN) candidate G29.4+0.1 and the complex radio source G29.37+0.1, which is likely formed by the superposition of a background radio galaxy and a Galactic supernova remnant (SNR). We investigate the possible connection between the SNR, the PWN G29.4+0.1, and HESS J1844-030 to shed light on the astrophysical origin of the TeV emission. METHODS: We performed an imaging and spectral study of the X-ray emission from the PWN G29.4+0.1 using archival observations obtained with the Chandra and XMM-Newton telescopes. Public radio continuum and HI data were used to derive distance constraints for the SNR that is linked to G29.37+0.1 and to investigate the interstellar medium where it is expanding. We applied a simple model of the evolution of a PWN inside an SNR to analyze the association between G29.4+0.1 and the radio emission from G29.37+0.1. We compared the spectral properties of the system with the population of TeV PWNe to investigate if HESS J1844-030 is the very high energy counterpart of the X-ray PWN G29.4+0.1. RESULTS: We conclude that G29.4+0.1 is a PWN and that a point source embedded on it is the powering pulsar. The HI data revealed that the SNR linked to G29.37+0.1 is a Galactic source at 6.5 kpc and expanding in a nonuniform medium. From the analysis of the pulsar motion and the pressure balance at the boundary of X-ray emission, we conclude that G29.4+0.1 could be a PWN that is located inside its host remnant, forming a new composite SNR. Based on the magnetic field of the PWN obtained from the X-ray luminosity, we found that the population of electrons producing synchrotron radiation in the keV band can also produce IC photons in the TeV band. This suggests that HESS J1844-030 could be the very high energy counterpart of G29.4+0.1.


title: Calculating energy derivatives for quantum chemistry on a quantum computer

abstract: Modeling chemical reactions and complicated molecular systems has been proposed as the `killer application' of a future quantum computer. Accurate calculations of derivatives of molecular eigenenergies are essential towards this end, allowing for geometry optimization, transition state searches, predictions of the response to an applied electric or magnetic field, and molecular dynamics simulations. In this work, we survey methods to calculate energy derivatives, and present two new methods: one based on quantum phase estimation, the other on a low-order response approximation. We calculate asymptotic error bounds and approximate computational scalings for the methods presented. Implementing these methods, we perform the world's first geometry optimization on an experimental quantum processor, estimating the equilibrium bond length of the dihydrogen molecule to within 0.014 Angstrom of the full configuration interaction value. Within the same experiment, we estimate the polarizability of the H2 molecule, finding agreement at the equilibrium bond length to within 0.06 a.u. (2% relative error).


title: Interactive Image Generation Using Scene Graphs

abstract: Recent years have witnessed some exciting developments in the domain of generating images from scene-based text descriptions. These approaches have primarily focused on generating images from a static text description and are limited to generating images in a single pass. They are unable to generate an image interactively based on an incrementally additive text description (something that is more intuitive and similar to the way we describe an image). We propose a method to generate an image incrementally based on a sequence of graphs of scene descriptions (scene-graphs). We propose a recurrent network architecture that preserves the image content generated in previous steps and modifies the cumulative image as per the newly provided scene information. Our model utilizes Graph Convolutional Networks (GCN) to cater to variable-sized scene graphs along with Generative Adversarial image translation networks to generate realistic multi-object images without needing any intermediate supervision during training. We experiment with Coco-Stuff dataset which has multi-object images along with annotations describing the visual scene and show that our model significantly outperforms other approaches on the same dataset in generating visually consistent images for incrementally growing scene graphs.


title: A catalog of known Galactic K-M stars of class I, candidate RSGs, in Gaia DR2

abstract: We investigate individual distances and luminosities of a sample of 889 nearby candidate red supergiants with reliable parallaxes (plx/plxerr > 4 and RUWE < 2.7) from Gaia DR2. The sample was extracted from the historical compilation of spectroscopically derived spectral types by Skiff (2014), and consists of K-M stars that are listed with class I at least once. The sample includes well-known red supergiants from Humphreys (1978), Elias et al. (1985), Jura and Kleinmann (1990), and Levesque et al. (2005). Infrared and optical measurements from the 2MASS, CIO, MSX, WISE, MIPSGAL, GLIMPSE, and NOMAD catalogs allow us to estimate the stellar bolometric magnitudes. We analyze the stars in the luminosity versus effective temperature plane and confirm that 43 sources are highly-probably red supergiants with Mbol<-7.1 mag. 43% of the sample is made of stars with masses > 7 Msun. Another 30% of the sample consists of giant stars.


title: PhysPort use and growth: Supporting physics teaching with research-based resources since 2011

abstract: PhysPort (www.physport.org) has become the go-to place for physics faculty to learn to apply research-based teaching and assessment in their classrooms. Usage has doubled every two years since the site was released in 2011, and 20% of all U.S. physics faculty are now verified educators on PhysPort. The lead author conceived of PhysPort in 2007 after meeting many physics faculty interested in incorporating results of physics education research (PER) in their classrooms but with no idea where to start. At the time, most PER results were buried in research journals not accessible to ordinary physics educators, and there was no central place to learn about results and implications for classroom practice. Throughout its development, PhysPort has been based on user research: we interview physics faculty about their needs, design the site based on those needs, and conduct usability testing to see how we meet those needs. PhysPort is a joint product of the American Association of Physics Teachers (AAPT) and Kansas State University (K-State), with contributions from many universities, funded by the National Science Foundation (NSF), and a library within the ComPADRE.org digital library. The site (originally called the PER User's Guide) was first funded by an NSF grant in 2009 and released in November 2011. This article presents an overview of resources on PhysPort, discussion of research and development of the site, and data on the continuing growth of site usage.


title: Research Note: An Open Source Bluespec Compiler

abstract: In this Research Note we report on an open-source compiler for the Bluespec hardware description language.


title: Approximate Bayesian computation with the Wasserstein distance

abstract: A growing number of generative statistical models do not permit the numerical evaluation of their likelihood functions. Approximate Bayesian computation (ABC) has become a popular approach to overcome this issue, in which one simulates synthetic data sets given parameters and compares summaries of these data sets with the corresponding observed values. We propose to avoid the use of summaries and the ensuing loss of information by instead using the Wasserstein distance between the empirical distributions of the observed and synthetic data. This generalizes the well-known approach of using order statistics within ABC to arbitrary dimensions. We describe how recently developed approximations of the Wasserstein distance allow the method to scale to realistic data sizes, and propose a new distance based on the Hilbert space-filling curve. We provide a theoretical study of the proposed method, describing consistency as the threshold goes to zero while the observations are kept fixed, and concentration properties as the number of observations grows. Various extensions to time series data are discussed. The approach is illustrated on various examples, including univariate and multivariate g-and-k distributions, a toggle switch model from systems biology, a queueing model, and a L\'evy-driven stochastic volatility model.


title: Arbitrarily large iterative tomographic reconstruction on multiple GPUs using the TIGRE toolbox

abstract: Tomographic image sizes keep increasing over time and while the GPUs that compute the tomographic reconstruction are also increasing in memory size, they are not doing so fast enough to reconstruct the largest datasets. This problem is often solved by reconstructing data in large clusters of GPUs with enough devices to fit the measured X-ray projections and reconstructed volume. Often this requires tens of GPUs, which is a very economically expensive solution. Access to single-node machines designed to reconstruct using just one or a few GPUs is more common in the field, but current software does not allow iterative reconstruction of volumes that do not fit in those GPUs. In this work, we propose a strategy to execute efficiently the required operations for iterative reconstruction for arbitrarily large images with any number of GPUs with arbitrarily small memories in a single node. Strategies for both the forward and backprojection operators are presented, along with two regularization approaches that are easily generalized to other projection types or regularizers. The proposed improvement also accelerates reconstruction of smaller images on single or multiple GPUs, providing faster code for time-critical medical applications. The resulting algorithm has been added to the TIGRE toolbox, a repository for iterative reconstruction algorithms for general CT, but this memory-saving and problem-splitting strategy can be easily adapted for any other GPU-based CT code.


title: MHD simulations of the in-situ generation of kink and sausage waves in the solar corona by collision of dense plasma clumps

abstract: MHD waves are ubiquitous in the solar corona where the highly structured magnetic fields provide efficient wave guides for their propagation. While MHD waves have been observed originating from lower layers of the solar atmosphere, recent studies have shown that some can be generated in-situ by the collision of dense counter propagating flows. In this theoretical study, we analyse the mechanism that triggers the propagation of kink and sausage modes in the solar corona, following the collision of counter-propagating flows and how the properties of the flows affect the properties of the generated waves. In order to study in detail such a mechanism we run a series of ideal 2D and 3D MHD simulations where we vary the properties of the counter-propagating flows and by means of a simple technique to estimate the amplitudes of the kink and sausage modes, we investigate their role on the generation and propagation of the MHD waves. We find that the amplitude of the waves is largely dependent on the kinetic energy of the flows and that the onset of kink or sausage modes depends on the asymmetries between the colliding blobs. Moreover, the initial wavelength of the MHD waves is associated with the magnetic configuration resulting from the collision of the flows. We also find that genuine 3D systems respond with smaller wave amplitudes. In this study, we present a parameter space description of the mechanism that leads to the generation of MHD waves from the collision of flows in the corona. Future observations of these waves can be used to understand the properties of the plasma and magnetic field of the solar corona.


title: Probing neutrino transition magnetic moments with coherent elastic neutrino-nucleus scattering

abstract: We explore the potential of current and next generation of coherent elastic neutrino-nucleus scattering (CE$\nu$NS) experiments in probing neutrino electromagnetic interactions. On the basis of a thorough statistical analysis, we determine the sensitivities on each component of the Majorana neutrino transition magnetic moment (TMM), $\left \vert \Lambda_i \right \vert$, that follow from low-energy neutrino-nucleus experiments. We derive the sensitivity to neutrino TMM from the first CE$\nu$NS measurement by the COHERENT experiment, at the Spallation Neutron Source. We also present results for the next phases of COHERENT using HPGe, LAr and NaI[Tl] detectors and for reactor neutrino experiments such as CONUS, CONNIE, MINER, TEXONO and RED100. The role of the CP violating phases in each case is also briefly discussed. We conclude that future CE$\nu$NS experiments with low-threshold capabilities can improve current TMM limits obtained from Borexino data.


title: The Hierarchy of Excitation Lifetimes in Two-Dimensional Fermi Gases

abstract: Momentum-conserving quasiparticle collisions in two-dimensional Fermi gases give rise to a large family of exceptionally long-lived excitation modes. The lifetimes of these modes exceed by a factor $(T_F/T)^2\gg 1$ the conventional Landau Fermi-liquid lifetimes $\tau\sim T_F/T^2$. The long-lived modes have a distinct angular structure, taking the form of $\cos m\theta$ and $\sin m\theta$ with odd $m$ values for a circular Fermi surface, with relaxation rate dependence on $m$ of the form $m^4\log m$, valid at not-too-large $m$. In contrast, the even-$m$ harmonics feature conventional lifetimes with a weak $m$ dependence. The long-time dynamics, governed by the long-lived modes, takes the form of angular (super)diffusion over the Fermi surface. Altogether, this leads to unusual long-time memory effects, defining an intriguing transport regime that lies between the conventional ballistic and hydrodynamic regimes.


abstract: Efficiency is crucial to the online recommender systems. Representing users and items as binary vectors for Collaborative Filtering (CF) can achieve fast user-item affinity computation in the Hamming space, in recent years, we have witnessed an emerging research effort in exploiting binary hashing techniques for CF methods. However, CF with binary codes naturally suffers from low accuracy due to limited representation capability in each bit, which impedes it from modeling complex structure of the data. In this work, we attempt to improve the efficiency without hurting the model performance by utilizing both the accuracy of real-valued vectors and the efficiency of binary codes to represent users/items. In particular, we propose the Compositional Coding for Collaborative Filtering (CCCF) framework, which not only gains better recommendation efficiency than the state-of-the-art binarized CF approaches but also achieves even higher accuracy than the real-valued CF method. Specifically, CCCF innovatively represents each user/item with a set of binary vectors, which are associated with a sparse real-value weight vector. Each value of the weight vector encodes the importance of the corresponding binary vector to the user/item. The continuous weight vectors greatly enhances the representation capability of binary codes, and its sparsity guarantees the processing speed. Furthermore, an integer weight approximation scheme is proposed to further accelerate the speed. Based on the CCCF framework, we design an efficient discrete optimization algorithm to learn its parameters. Extensive experiments on three real-world datasets show that our method outperforms the state-of-the-art binarized CF methods (even achieves better performance than the real-valued CF method) by a large margin in terms of both recommendation accuracy and efficiency.


title: $\alpha$-attractor dark energy in view of next-generation cosmological surveys

abstract: The $\alpha$-attractor inflationary models are nowadays favored by CMB Planck observations. Their similarity with canonical quintessence models motivates the exploration of a common framework that explains both inflation and dark energy. We study the expected constraints that next-generation cosmological experiments will be able to impose for the dark energy $\alpha$-attractor model. We systematically account for the constraining power of SNIa from WFIRST, BAO from DESI and WFIRST, galaxy clustering and shear from LSST and Stage-4 CMB experiments. We assume a tensor-to-scalar ratio, $10^{-3} < r < 10^{-2}$, which permits to explore the wide regime sufficiently close, but distinct, to a cosmological constant, without need of fine tunning the initial value of the field. We find that the combination S4CMB + LSST + SNIa will achieve the best results, improving the FoM by almost an order of magnitude; respect to the S4CMB + BAO + SNIa case. We find this is also true for the FoM of the $w_0 - w_a$ parameters. Therefore, future surveys will be uniquely able to probe models connecting early and late cosmic acceleration.


title: Sharp asymptotics for Fredholm Pfaffians related to interacting particle systems and random matrices

abstract: It has been known since the pioneering paper of Mark Kac, that the asymptotics of Fredholm determinants can be studied using probabilistic methods. We demonstrate the efficacy of Kac' approach by studying the Fredholm Pfaffian describing the statistics of both non-Hermitian random matrices and annihilating Brownian motions. Namely, we establish the following two results. Firstly, let $\sqrt{N}+\lambda_{max}$ be the largest real eigenvalue of a random $N\times N$ matrix with independent $N(0,1)$ entries (the `real Ginibre matrix'). Consider the limiting $N\rightarrow \infty$ distribution $\mathbb{P}[\lambda_{max}<-L]$ of the shifted maximal real eigenvalue $\lambda_{max}$. Then \[ \lim_{L\rightarrow \infty} e^{\frac{1}{2\sqrt{2\pi}}\zeta\left(\frac{3}{2}\right)L} \mathbb{P}\left(\lambda_{max}<-L\right) =e^{C_e}, \] where $\zeta$ is the Riemann zeta-function and \[ C_e=\frac{1}{2}\log 2+\frac{1}{4\pi}\sum_{n=1}^{\infty}\frac{1}{n} \left(-\pi+\sum_{m=1}^{n-1}\frac{1}{\sqrt{m(n-m)}}\right). \] Secondly, let $X_t^{(max)}$ be the position of the rightmost particle at time $t$ for a system of annihilating Brownian motions (ABM's) started from every point of $\mathbb{R}_{-}$. Then \[ \lim_{L\rightarrow \infty} e^{\frac{1}{2\sqrt{2\pi}}\zeta\left(\frac{3}{2}\right)L} \mathbb{P}\left(\frac{X_{t}^{(max)}}{\sqrt{4t}}<-L\right) =e^{C_e}. \] These statements are a sharp counterpart of our previous results improved by computing the terms of order $L^{0}$ in the asymptotic expansion of the corresponding Fredholm Pfaffian.


title: Statistical analysis of interplanetary magnetic field path lengths from solar energetic electron events observed by WIND

abstract: We calculate the interplanetary magnetic field path lengths traveled by electrons in solar electron events detected by the WIND 3DP instrument from $1994$ to $2016$. The velocity dispersion analysis method is applied for electrons at energies of $\sim$ $27$ keV to $310$ keV. Previous velocity dispersion analyses employ the onset times, which are often affected by instrumental effects and the pre-existing background flux, leading to large uncertainties. We propose a new method here. Instead of using the peak or onset time, we apply the velocity dispersion analysis to the times that correspond to the rising phase of the fluxes that are a fraction, $\eta$, of the peak flux. We perform statistical analysis on selected events whose calculated path lengths have uncertainties smaller than $0.1$ AU. The mean and standard deviation, ($\mu$, $\sigma$), of the calculated path lengths corresponding to $\eta=$ $3/4$, $1/2$, and $1/3$ of the peak flux is ($1.17$ AU, $0.17$ AU), ($1.11$ AU, $0.14$ AU), and ($1.06$ AU, $0.15$ AU). The distribution of the calculated path lengths is also well fitted by a Gaussian distribution for the $\eta=3/4$ and $1/3$ cases. These results suggest that in these electron events the interplanetary magnetic field topology is close to the nominal Parker spiral with little field line meandering. Our results have important implications for particles' perpendicular diffusion.


title: Measuring Spin of the Remnant Black Hole from Maximum Amplitude

abstract: Gravitational waves emitted during the merger of two black holes carry information about the remnant black hole, namely its mass and spin. This information is typically found from the ringdown radiation as the black hole settles to a final state. We find that the remnant black hole spin is already known at the peak amplitude of the gravitational wave strain. Using this knowledge, we present a new method for measuring the final spin that is template independent, using only the chirp mass, the instantaneous frequency of the strain and its derivative at maximum amplitude, all template independent.


title: Enhanced stability of the focus obtained by wavefront optimization in dynamical scattering media

abstract: Focusing scattered light using wavefront shaping provides interesting perspectives to image deep in opaque samples, as e.g. in nonlinear fluorescence microscopy. Applying these technics to in vivo imaging remains challenging due to the short decorrelation time of the speckle in depth, as focusing and imaging has to be achieved within the order of the decorrelation time. In this paper, we experimentally study the focus lifetime after focusing through dynamical scattering media, when iterative wavefront optimization and speckle decorrelation occur over the same timescale. We show experimental situations with heterogeneous stability of the scattering sequences, where the focus presents significantly higher stability than the surrounding speckle.


title: Super-pancyclic hypergraphs and bipartite graphs

abstract: We find Dirac-type sufficient conditions for a hypergraph $\mathcal H$ with few edges to be hamiltonian. We also show that these conditions provide that $\mathcal H$ is {\em super-pancyclic}, i.e., for each $A \subseteq V(\mathcal H)$ with $|A| \geq 3$, $\mathcal H$ contains a Berge cycle with vertex set $A$. We mostly use the language of bipartite graphs, because every bipartite graph is the incidence graph of a multihypergraph. In particular, we extend some results of Jackson on the existence of long cycles in bipartite graphs where the vertices in one part have high minimum degree. Furthermore, we prove a conjecture of Jackson from 1981 on long cycles in 2-connected bipartite graphs.


title: Nonlinear Polariton Fluids in a Flatband Reveal Discrete Gap Solitons

abstract: Phase frustration in periodic lattices is responsible for the formation of dispersionless flat bands. The absence of any kinetic energy scale makes flat band physics critically sensitive to perturbations and interactions. We report here on the experimental investigation of the nonlinear dynamics of cavity polaritons in the gapped flat band of a one-dimensional Lieb lattice. We observe the formation of gap solitons with quantized size and very abrupt edges, signature of the frozen propagation of switching fronts. This type of gap solitons belongs to the class of truncated Bloch waves, and had only been observed in closed systems up to now. Here the driven-dissipative character of the system gives rise to a complex multistability of the nonlinear domains generated in the flat band. These results open up interesting perspective regarding more complex 2D lattices and the generation of correlated photon phases.


title: Monte Carlo Co-Ordinate Ascent Variational Inference

abstract: In Variational Inference (VI), coordinate-ascent and gradient-based approaches are two major types of algorithms for approximating difficult-to-compute probability densities. In real-world implementations of complex models, Monte Carlo methods are widely used to estimate expectations in coordinate-ascent approaches and gradients in derivative-driven ones. We discuss a Monte Carlo Co-ordinate Ascent VI (MC-CAVI) algorithm that makes use of Markov chain Monte Carlo (MCMC) methods in the calculation of expectations required within Co-ordinate Ascent VI (CAVI). We show that, under regularity conditions, an MC-CAVI recursion will get arbitrarily close to a maximiser of the evidence lower bound (ELBO) with any given high probability. In numerical examples, the performance of MC-CAVI algorithm is compared with that of MCMC and -- as a representative of derivative-based VI methods -- of Black Box VI (BBVI). We discuss and demonstrate MC-CAVI's suitability for models with hard constraints in simulated and real examples. We compare MC-CAVI's performance with that of MCMC in an important complex model used in Nuclear Magnetic Resonance (NMR) spectroscopy data analysis -- BBVI is nearly impossible to be employed in this setting due to the hard constraints involved in the model.


title: Deep Learning for TDD and FDD Massive MIMO: Mapping Channels in Space and Frequency

abstract: Can we map the channels at one set of antennas and one frequency band to the channels at another set of antennas---possibly at a different location and a different frequency band? If this channel-to-channel mapping is possible, we can expect dramatic gains for massive MIMO systems. For example, in FDD massive MIMO, the uplink channels can be mapped to the downlink channels or the downlink channels at one subset of antennas can be mapped to the downlink channels at all the other antennas. This can significantly reduce (or even eliminate) the downlink training/feedback overhead. In the context of cell-free/distributed massive MIMO systems, this channel mapping can be leveraged to reduce the fronthaul signaling overhead as only the channels at a subset of the distributed terminals need to be fed to the central unit which can map them to the channels at all the other terminals. This mapping can also find interesting applications in mmWave beam prediction, MIMO radar, and massive MIMO based positioning. In this paper, we introduce the new concept of channel mapping in space and frequency, where the channels at one set of antennas and one frequency band are mapped to the channels at another set of antennas and frequency band. First, we prove that this channel-to-channel mapping function exists under the condition that the mapping from the candidate user positions to the channels at the first set of antennas is bijective; a condition that can be achieved with high probability in several practical MIMO communication scenarios. Then, we note that the channel-to-channel mapping function, even if it exists, is typically unknown and very hard to characterize analytically as it heavily depends on the various elements of the surrounding environment. With this motivation, we propose to leverage the powerful learning capabilities of deep neural networks ....


title: Semidirect Products and Wreath Products of Localities

abstract: We develop a theory of semidirect products of partial groups and localities. Our concepts generalize the notions of direct products of partial groups and localities, and of semidirect products of groups.


title: Diffractive dijet production in impact parameter dependent saturation models

abstract: We study coherent diffractive dijet production in electron-hadron and electron-nucleus collisions within the dipole picture. We provide semi-analytic results for the differential cross section and elliptic anisotropy in the angle between dijet transverse momentum and hadron recoil momentum. We demonstrate the direct relation between angular moments of the dipole amplitude in coordinate space and angular moments of the diffractive dijet cross sections. To perform explicit calculations we employ two different saturation models, extended to include the target geometry. In the limit of large photon virtuality or quark masses, we find fully analytic results that allow direct insight into how the differential cross section and elliptic anisotropy depend on the saturation scale, target geometry, and kinematic variables. We further provide numerical results for more general kinematics in collisions at a future electron-ion collider, and study the effects of approaching the saturated regime on diffractive dijet observables.


title: Higgs Boson Studies at Future Particle Colliders

abstract: This document aims to provide an assessment of the potential of future colliding beam facilities to perform Higgs boson studies. The analysis builds on the submissions made by the proponents of future colliders to the European Strategy Update process, and takes as its point of departure the results expected at the completion of the HL-LHC program. This report presents quantitative results on many aspects of Higgs physics for future collider projects of sufficient maturity using uniform methodologies. A first version of this report was prepared for the purposes of discussion at the Open Symposium in Granada (13-16/05/2019). Comments and feedback received led to the consideration of additional run scenarios as well as a refined analysis of the impact of electroweak measurements on the Higgs coupling extraction.


title: Non-Relativistic and Relativistic Equations for the Kratzer Potential plus a Dipole in 2D Systems

abstract: In this work, we study the wave equations in 2D Euclidian space for a new non-central potential consisting of a Kratzer term and a dipole term. For Schrodinger equation, we obtain the analytical expressions of the energies and the wave functions of the system. For Klein-Gordon and Dirac equations, we do the study in both spin and pseudo-spin symmetries to get the eigenfunctions and the eingenvalues. Then we study the dependence of energies on the radial moment D_{r} and the angular moment D_{theta}. We find that the D_{theta} term tends to dissociate the system, and thus counteracts the Coulomb binding effect, and that the D_{r} term can either amplify or decrease this effect according to its sign.


title: AutoRegressive Planet Search: Feasibility Study for Irregular Time Series

abstract: Sensitive signal processing methods are needed to detect transiting planets from ground-based photometric surveys. Caceres et al. (2019) show that the AutoRegressive Planet Search (ARPS) method --- a combination of autoregressive integrated moving average (ARIMA) parametric modeling, a new Transit Comb Filter (TCF) periodogram, and machine learning classification --- is effective when applied to evenly spaced light curves from space-based missions. We investigate here whether ARIMA and TCF will be effective for ground-based survey light curves that are often sparsely sampled with high noise levels from atmospheric and instrumental conditions. The ARPS procedure is applied to selected light curves with strong planetary signals from the Kepler mission that have been altered to simulate the conditions of ground-based exoplanet surveys. Typical irregular cadence patterns are used from the HATSouth survey. We also evaluate recovery of known planets from HATSouth. Simulations test transit signal recovery as a function of cadence pattern and duration, stellar magnitude, planet orbital period and transit depth. Detection rates improve for shorter periods and deeper transits. The study predicts that the ARPS methodology will detect planets with $\gtrsim 0.1$\% transit depth and periods $\lesssim 40$ days in HATSouth stars brighter than $\sim$15 mag. ARPS methodology is therefore promising for planet discovery from ground-based exoplanet surveys with sufficiently dense cadence patterns.


title: Learning Interpretable Features via Adversarially Robust Optimization

abstract: Neural networks are proven to be remarkably successful for classification and diagnosis in medical applications. However, the ambiguity in the decision-making process and the interpretability of the learned features is a matter of concern. In this work, we propose a method for improving the feature interpretability of neural network classifiers. Initially, we propose a baseline convolutional neural network with state of the art performance in terms of accuracy and weakly supervised localization. Subsequently, the loss is modified to integrate robustness to adversarial examples into the training process. In this work, feature interpretability is quantified via evaluating the weakly supervised localization using the ground truth bounding boxes. Interpretability is also visually assessed using class activation maps and saliency maps. The method is applied to NIH ChestX-ray14, the largest publicly available chest x-rays dataset. We demonstrate that the adversarially robust optimization paradigm improves feature interpretability both quantitatively and visually.


title: The Return of the WIMP: Missing Energy Signals and the Galactic Center Excess

abstract: In a recent work, we emphasized that an excess in tri-lepton events plus missing energy observed by the ATLAS experiment at the LHC could be interpreted as a signal of low energy supersymmetry. In such a scenario the lightest neutralino mass is approximately $m_\chi \simeq 60$ GeV and the direct Dark Matter detection cross section is naturally below the current bound. In this work we present simple extensions of this scenario that lead to an explanation of the gamma ray excess at the center of the galaxy observed by Fermi-LAT, as well as the anti-proton excess observed by AMS-02. These extensions include the addition of a small CP violating phase in the neutralino sector or the addition of a light CP-odd Higgs scalar. Our study is of special relevance in view of a recent analysis that casts doubt on the previously accepted preference for mili-second pulsars as the origin of the galactic center excess.


title: Homogenous $In_{x}Ga_{1-x}N$ alloys on ZnO substrates: A new approach for high performance thermoelectric materials

abstract: High performance thermoelectric materials for wide-range temperature applications still remains a challenge. In this study, we have produced high-quality homogeneous $In_{0.32}Ga_{0.68}N$ on ZnO substrates, with no phase separation at high Indium content, using metal organic chemical vapor deposition for thermoelectric applications. A record high room temperature figure of merit zT is obtained of 0.86, which is five times larger than that of SiGe, the current state of the art high temperature thermoelectric material. These materials are shown to have a nearly perfect doping concentration to maximize zT regardless of the scattering mechanism. This almost one order of magnitude increase in zT is due to large electrical conductivities from oxygen co-doping as well as low thermal conductivities from alloy scattering. The maximum power factor reached was $77.98x10^{-4} W/mK^{2}$ at 300K for $In_{0.32}Ga_{0.68}N$ alloys at a carrier concentration $~6.25x10^{20} cm^{-3}$. This work indicates that $In_{x}Ga_{1-x}N$ alloys have great potential for thermoelectric applications especially at a high temperature range.


title: QSO2 outflow characterization using data obtained with OSIRIS at the Gran Telescopio Canarias

abstract: Ionized outflows are ubiquitous in non radio-loud obscured quasars (QSO2s) at different redshifts. We search for large-scale ionized outflows associated with six optically selected QSO2 (five non-radio-loud and one radio-loud) at $z\sim$ 0.2-0.5, targeting objects with extended radio structures. We have obtained OSIRIS/GTC optical long slit spectroscopy data for these six QSO2 with the slit located along the radio axis. We traced the gas kinematics with the [OIII]$\lambda$4959,5007 lines to investigate ionized outflows and characterize the dynamical state of the host galaxies. This second study has been complemented with previously published FORS2/VLT spectroscopic data of 13 more QSO2 at similar z. We identify ionized outflows in four out of the six QSO2 observed with the GTC. The outflows are spatially unresolved in two QSO2 and compact in a third (radial size of R=0.8$\pm$0.3 kpc). Of particular interest is the radio-quiet QSO2 SDSS 0741+3020 at z=0.47, associated with a giant $\sim$ 112 kpc nebula. An ionized outflow probably induced by the radio structures has been detected along the axis defined by the central $\sim$1\arcsec\ radio structure, extending up to at least $\sim$4 kpc from the active galactic nucleus (AGN). Turbulent gas ($\sigma\sim$130 km s$^{-1}$) has also been detected across the giant gas nebula up to $\sim$40 kpc from the AGN. This turbulence may have been induced by outflows triggered by the interaction between a so-far undetected large-scale radio source and the nebula. Regarding the dynamical state of the host galaxies, we find that the majority of the QSO2 show v/$\sigma<$ 1, implying that they are dominated by random motions (so-called dispersion-dominated systems). Most (17 of 19) fall in the area of the E/S0 galaxies in the dynamical diagram v/$\sigma$ versus $\sigma$. None are consistent with spiral or disk galaxies.


title: Resummed prediction for Higgs boson production through $b\bar{b}$ annihilation at N$^3$LL

abstract: We present an accurate theoretical prediction for the production of Higgs boson through bottom quark annihilation at the LHC up to next-to-next-to-next-to leading order (N$^3$LO) plus next-to-next-to-next-to-leading logarithmic (N$^3$LL) accuracy. We determine the third order perturbative Quantum Chromodynamics (QCD) correction to the process dependent constant in the resummed expression using the three loop bottom quark form factor and third order quark soft distribution function. Thanks to the recent computation of N$^3$LO corrections to this production cross-section from all the partonic channels, an accurate matching can be obtained for a consistent predictions at N$^3$LO+N$^3$LL accuracy in QCD. We have studied in detail the impact of resummed threshold contributions to inclusive cross-sections at various centre-of-mass energies and also discussed their sensitivity to renormalization and factorization scales at next-to-next-to leading order (NNLO) matched with next-to-next-to leading logarithm (NNLL). At N$^3$LO+N$^3$LL, we predict the cross-section for different centre-of-mass energies using the recently available results in \cite{Duhr:2019kwi} as well as study the renormalization scale dependence at the same order.


title: Anatomy of the $tthh$ Physics at HL-LHC

abstract: The $tthh$ production at colliders contain rich information on the nature of Higgs boson. In this article, we systematically studied its physics at High-Luminosity Large Hadron Collider (HL-LHC), using exclusive channels with multiple ($\geq 5$) $b$-jets and one lepton ($5b1\ell$), multiple ($\geq 5$) $b$-jets and opposite-sign di-lepton ($5b2\ell$), same-sign di-lepton (SS2$\ell$), multiple leptons (multi-$\ell$), and di-tau resonance ($\tau\tau$). The scenarios analyzed include: (1) the $tthh$ production in Standard Model; (2) the $tthh$ production mediated by anomalous cubic Higgs self-coupling and $tthh$ contact interaction; (3) heavy Higgs ($H$) production with $tt H \to tthh$; and (4) pair production of fermionic top partners ($T$) with $T T \to tthh$. To address the complication of event topologies and the mess of combinatorial backgrounds, a tool of Boosted-Decision-Tree was applied in the analyses. The $5b1\ell$ and SS2$\ell$ analyses define the two most promising channels, resulting in slightly different sensitivities. For non-resonant $tthh$ production, a combination of these exclusive analyses allows for its measurment in the SM with a statistical significance $\sim 0.9\sigma$ (with $S/B > 1 \%$), and may assist partially breaking the sensitivity degeneracy w.r.t. the cubic Higgs self-coupling, a difficulty usually thought to exist in gluon fusion di-Higgs analysis at HL-LHC. These sensitivities were also projected to future hadron colliders at 27 TeV and 100 TeV. For resonant $tthh$ productions, the heavy Higgs boson in type II Two-Higgs-Doublet-Model could be efficiently searched for between the mass thresholds $2 m_h < m_H < 2 m_t$ and even beyond that, for relatively small $\tan\beta$, while the fermionic top partners in composite Higgs models could be probed for up to $\sim 1.5$ TeV and $\sim 1.7$ TeV, for Br$(T\to th)=25\%$ and $50\%$, respectively.


title: The Vickrey Auction with a Single Duplicate Bidder Approximates the Optimal Revenue

abstract: Bulow and Klemperer's well-known result states that, in a single-item auction where the $n$ bidders' values are independently and identically drawn from a regular distribution, the Vickrey auction with one additional bidder (a duplicate) extracts at least as much revenue as the optimal auction without the duplicate. Hartline and Roughgarden, in their influential 2009 paper, removed the requirement that the distributions be identical, at the cost of allowing the Vickrey auction to recruit $n$ duplicates, one from each distribution, and relaxing its revenue advantage to a $2$-approximation. In this work we restore Bulow and Klemperer's number of duplicates in Hartline and Roughgarden's more general setting with a worse approximation ratio. We show that recruiting a duplicate from one of the distributions suffices for the Vickrey auction to $10$-approximate the optimal revenue. We also show that in a $k$-items unit demand auction, recruiting $k$ duplicates suffices for the VCG auction to $O(1)$-approximate the optimal revenue. As another result, we tighten the analysis for Hartline and Roughgarden's Vickrey auction with $n$ duplicates for the case with two bidders in the auction. We show that in this case the Vickrey auction with two duplicates obtains at least $3/4$ of the optimal revenue. This is tight by meeting a lower bound by Hartline and Roughgarden. En route, we obtain a transparent analysis of their $2$-approximation for $n$~bidders, via a natural connection to Ronen's lookahead auction.


title: Pair-density-wave superconductivity of faces, edges and vertices in systems with imbalanced fermions

abstract: We describe boundary effects in superconducting systems with Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconducting instability, using Bogoliubov-de-Gennes and Ginzburg-Landau (GL) formalisms. First, we show that in dimensions larger than one the standard GL functional formalism for FFLO superconductors is unbounded from below. This is demonstrated by finding solutions with zero Laplacian terms near boundaries. We generalize the GL formalism for these systems by retaining higher order terms. Next, we demonstrate that a cuboid sample of a superconductor with imbalanced fermions at a mean-field level has a sequence of the phase transitions. At low temperatures it forms Larkin-Ovchinnikov state in the bulk but has a different modulation pattern close to the boundaries. When temperature is increased the first phase transition occurs when the bulk of the material becomes normal while the faces remain superconducting. The second transition occurs at higher temperature where the system retains superconductivity on the edges. The third transition is associated with the loss of edge superconductivity while retaining superconducting gap in the vertices. We obtain the same sequence of phase transition by numerically solving the Bogoliubov-de Gennes model.


title: Discrete gauge groups in certain F-theory models in six dimensions

abstract: We construct six-dimensional (6D) F-theory models in which discrete $\mathbb{Z}_5, \mathbb{Z}_4, \mathbb{Z}_3,$ and $\mathbb{Z}_2$ gauge symmetries arise. We demonstrate that a special family of "Fano 3-folds" is a useful tool for constructing the aforementioned models. The geometry of Fano 3-folds in the constructions of models can be useful for understanding discrete gauge symmetries in 6D F-theory compactifications. We argue that the constructions of the aforementioned models are applicable to Calabi-Yau genus-one fibrations over any base space, except models with a discrete $\mathbb{Z}_5$ gauge group. We construct 6D F-theory models with a discrete $\mathbb{Z}_5$ gauge group over the del Pezzo surfaces, as well as over $\mathbb{P}^1\times\mathbb{P}^1$ and $\mathbb{P}^2$. We also discuss some applications to four-dimensional F-theory models with discrete gauge symmetries.


title: The Effect of Network Width on Stochastic Gradient Descent and Generalization: an Empirical Study

abstract: We investigate how the final parameters found by stochastic gradient descent are influenced by over-parameterization. We generate families of models by increasing the number of channels in a base network, and then perform a large hyper-parameter search to study how the test error depends on learning rate, batch size, and network width. We find that the optimal SGD hyper-parameters are determined by a "normalized noise scale," which is a function of the batch size, learning rate, and initialization conditions. In the absence of batch normalization, the optimal normalized noise scale is directly proportional to width. Wider networks, with their higher optimal noise scale, also achieve higher test accuracy. These observations hold for MLPs, ConvNets, and ResNets, and for two different parameterization schemes ("Standard" and "NTK"). We observe a similar trend with batch normalization for ResNets. Surprisingly, since the largest stable learning rate is bounded, the largest batch size consistent with the optimal normalized noise scale decreases as the width increases.


title: Holographic Shannon Entropy: The Outer Entropy of Entanglement Wedges

abstract: We introduce a simple geometrical construction similar to covariant holographic entanglement entropy but with the addition of a new term proportional to boundary region volume. This new procedure has properties strongly resembling classical Shannon entropy of probability distributions rather than von Neumann entropy, so we call the quantity holographic Shannon entropy. The holographic Shannon entropy of a region $A$ is divergent in AdS/CFT, but upon regulation, it appears to be equal to the outer entropy of the entanglement wedge of $A$ with the entanglement wedge of the complement of $A$ held fixed. The construction is unambiguous when applied to compact surfaces with convex shape in general spacetimes obeying the null curvature condition (of which AdS with a regulator is a special case). In this context we prove that holographic Shannon entropy obeys monotonicity, a key property of Shannon entropy, as well as all known balanced inequalities of dynamical holographic entanglement entropy. In the static case, we explain why there must exist some classical probability distribution on random variables locally distributed on the boundary with the property that the Shannon entropies of all marginals are exactly reproduced by the holographic Shannon entropy formula.


title: Splitting hairs with transcendental entire functions

abstract: In recent years, there has been significant progress in the understanding of the dynamics of transcendental entire functions with bounded postsingular set. In particular, for certain classes of such functions, a complete description of their topological dynamics in terms of a simpler model has been given inspired by methods from polynomial dynamics. In this paper, and for the first time, we give analogous results in cases when the postsingular set is unbounded. More specifically, we show that if $f$ is of finite order, has bounded criticality on its Julia set $J(f)$, and its singular set consists of finitely many critical values that escape to infinity and satisfy a certain separation condition, then $J(f)$ is a collection of dynamic rays or hairs, that split at critical points, together with their corresponding landing points. In fact, our result holds for a much larger class of functions with bounded singular set. Moreover, this result is a consequence of a significantly more general one: we provide a topological model for the action of $f$ on its Julia set.


title: Non-Gaussianity in the Weak Lensing Correlation Function Likelihood -- Implications for Cosmological Parameter Biases

abstract: We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in one-dimensional marginal distributions of shear two-point correlation functions in simulated weak lensing data. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in one-dimensional posteriors of $\Omega_{\rm m}$ and $\sigma_{\rm 8}$ due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data ($100$ deg$^{2}$). We also present a systematic approach to constructing approximate multivariate likelihoods with one-dimensional parametric functions by assuming independence or more flexible non-parametric multivariate methods after decorrelating the data points using principal component analysis (PCA). While the use of PCA does not modify the non-Gaussianity of the multivariate likelihood, we find empirically that the one-dimensional marginal sampling distributions of the PCA components exhibit less skewness and kurtosis than the original shear correlation functions.Modeling the likelihood with marginal parametric functions based on the assumption of independence between PCA components thus gives a lower limit for the biases. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of $\sim$5.


title: Low-dimensional de Sitter quantum gravity

abstract: We study aspects of Jackiw-Teitelboim (JT) quantum gravity in two-dimensional nearly de Sitter (dS) spacetime, as well as pure de Sitter quantum gravity in three dimensions. These are each theories of boundary modes, which include a reparameterization field on each connected component of the boundary as well as topological degrees of freedom. In two dimensions, the boundary theory is closely related to the Schwarzian path integral, and in three dimensions to the quantization of coadjoint orbits of the Virasoro group. Using these boundary theories we compute loop corrections to the wavefunction of the universe, and investigate gravitational contributions to scattering. Along the way, we show that JT gravity in dS$_2$ is an analytic continuation of JT gravity in Euclidean AdS$_2$, and that pure gravity in dS$_3$ is a continuation of pure gravity in Euclidean AdS$_3$. We define a genus expansion for de Sitter JT gravity by summing over higher genus generalizations of surfaces used in the Hartle-Hawking construction. Assuming a conjecture regarding the volumes of moduli spaces of such surfaces, we find that the de Sitter genus expansion is the continuation of the recently discovered AdS genus expansion. Then both may be understood as coming from the genus expansion of the same double-scaled matrix model, which would provide a non-perturbative completion of de Sitter JT gravity.


title: Revisiting the Temperature of the Diffuse ISM with CHESS Sounding Rocket Observations

abstract: Measuring the temperature and abundance patterns of clouds in the interstellar medium (ISM) provides an observational basis for models of the physical conditions within the clouds, which play an important role in studies of star and planet formation. The Colorado High-resolution Echelle Stellar Spectrograph (CHESS) is a far ultraviolet rocket-borne instrument designed to study the atomic-to-molecular transitions within diffuse molecular and translucent cloud regions. The final two flights of the instrument observed $\beta^{1}$ Scorpii ($\beta$ Sco) and $\gamma$ Arae. We present flight results of interstellar molecular hydrogen (H$_{\rm 2}$) excitation on the sightlines, including measurements of the column densities and temperatures. These results are compared to previous values that were measured using the damping wings of low J$^{\prime \prime}$ H$_{\rm 2}$ absorption features (Savage et al. 1977). For $\beta$ Sco, we find that the derived column density of the J$^{\prime \prime}$ = 1 rotational level differs by a factor of 2-3 when compared to the previous observations. We discuss the discrepancies between the two measurements and show that the source of the difference is due to the opacity of higher rotational levels contributing to the J$^{\prime \prime}$ = 1 absorption wing, increasing the inferred column density in the previous work. We extend this analysis to 9 $Copernicus$ and 13 $FUSE$ spectra to explore the interdependence of the column densities of different rotational levels and how the H$_{\rm 2}$ kinetic temperature is influenced by these relationships. We find a revised average gas kinetic temperature of the diffuse molecular ISM of T$_{01}$ = 68 $\pm$ 13 K, 12% lower than the value found previously.


title: Super-resolution limit of the ESPRIT algorithm

abstract: The problem of imaging point objects can be formulated as estimation of an unknown atomic measure from its $M+1$ consecutive noisy Fourier coefficients. The standard resolution of this inverse problem is $1/M$ and super-resolution refers to the capability of resolving atoms at a higher resolution. When any two atoms are less than $1/M$ apart, this recovery problem is highly challenging and many existing algorithms either cannot deal with this situation or require restrictive assumptions on the sign of the measure. ESPRIT is an efficient method that does not depend on the sign of the measure. This paper provides an explicit error bound on the support matching distance of ESPRIT in terms of the minimum singular value of Vandermonde matrices. When the support consists of multiple well-separated clumps and noise is sufficiently small, the support error by ESPRIT scales like ${\rm SRF}^{2\lambda+2} \times {\rm Noise}$, where the Super-Resolution Factor (${\rm SRF}$) governs the difficulty of the problem and $\lambda$ is the cardinality of the largest clump. {If the support contains one clump of closely spaced atoms, the min-max error is ${\rm SRF}^{2\lambda +2} \times {\rm Noise}/M$. Our error bound matches the min-max rate up to a factor of $M$ in the small noise regime. Our results therefore establishes the near-optimality of ESPRIT,} and our theory is validated by numerical experiments.


title: Thermo-magnetic corrections to $\pi$-$\pi$ Scattering Lengths in the Linear Sigma Model

abstract: In this article, we extend our previous study of the $\pi$-$\pi$ scattering lengths under the presence of an external magnetic field, including finite temperature effects. The novelty of this work is precisely the introduction of temperature into the discussion, and its interplay with the magnetic field. As in the previous article, we base our analysis in the linear sigma-model, and our calculations are exact within this context. Although the effects are comparatively small, it is interesting to remark that magnetic field and temperature display opposite effects over the scattering lengths.


title: Relation between gravitational mass and baryonic mass for non-rotating and rapidly rotating neutron stars

abstract: With a selected sample of neutron star (NS) equation-of-states (EOSs) that are consistent with the current observations and have a range of maximum masses, we investigate the relations between NS gravitational mass $M_g$ and baryonic mass $M_b$, and the relations between the maximum NS mass supported through uniform rotation ($M_{\rm max}$) and that of nonrotating NSs ($M_{\rm TOV}$). We find that if one intends to apply an EOS-independent quadratic, universal transformation formula ($M_b=M_g+A\times M_{g}^2$) to all EOSs, the best fit $A$ value is 0.080 for non-rotating NSs only and 0.073 when different spin periods are considered. The residual error of the transformation is as large as $\sim0.1M_{\odot}$. For different EOSs, we find that the parameter $A$ for non-rotating NSs is proportional to $R_{1.4}^{-1}$ (where $R_{1.4}$ is NS radius for 1.4$M_\odot$ in unit of km). For a particular EOS, if one adopts the best-fit parameters for different spin periods, the residual error of the transformation is smaller, which is of the order of 0.01$M_\odot$ for the quadratic form and less than 0.01$M_\odot$ for the cubic form ($M_b=M_g+A_1\times M_{g}^2+A_2\times M_{g}^3$). We also find a very tight and general correlation between the normalized mass gain due to spin $\Delta m\equiv(M_{\rm max}-M_{\rm TOV})/M_{\rm TOV}$ and the spin period normalized to the Keplerian period ${\cal P}$, i.e. ${\rm log_{10}}\Delta m = (-2.74\pm0.05){\rm log_{10}}{\cal P}+{\rm log_{10}}(0.20\pm 0.01)$, which is independent of EOS models. Applications of our results to GW170817 is discussed.


title: Misner Gravitational Charges and Variable String Strengths

abstract: As shown recently (arXiv:1903.08668), consistent thermodynamics of the Lorentzian Taub-NUT solutions with Misner strings present can be formulated provided a new pair of conjugate quantities (related to the NUT parameter) $\psi-N$ is introduced. In (arXiv:1903.08668) this pair was calculated from the Euclidean action but no geometrical interpretation for the new quantities was provided. In this paper we propose that the potential $\psi$ should be identified with the surface gravity of the Misner string and the conjugate Misner charge $N$ can be obtained by a Komar-type integration over the tubes surrounding the string singularities. We show that similar tube contributions also modify the Komar formula for the thermodynamic volume. To render the integrals finite we employ the method of Killing co-potentials. By construction the new charges then satisfy the Smarr relation. Equipped with these geometrical notions, we generalize the first law for the (possibly charged) Taub-NUT spacetimes to account for asymmetric distributions of Misner strings and their potential variable strengths.


title: Towards an Understanding of the Resolution Dependence of Core-Collapse Supernova Simulations

abstract: Using our new state-of-the-art core-collapse supernova (CCSN) code Fornax, we explore the dependence upon spatial resolution of the outcome and character of three-dimensional (3D) supernova simulations. For the same 19-M$_{\odot}$ progenitor star, energy and radial binning, neutrino microphysics, and nuclear equation of state, changing only the number of angular bins in the $\theta$ and $\phi$ directions, we witness that our lowest resolution 3D simulation does not explode. However, when jumping progressively up in resolution by factors of two in each angular direction on our spherical-polar grid, models then explode, and explode slightly more vigorously with increasing resolution. This suggests that there can be a qualitative dependence of the outcome of 3D CCSN simulations upon spatial resolution. The critical aspect of higher spatial resolution is the adequate capturing of the physics of neutrino-driven turbulence, in particular its Reynolds stress. The greater numerical viscosity of lower-resolution simulations results in greater drag on the turbulent eddies that embody turbulent stress, and, hence, in a diminution of their vigor. Turbulent stress not only pushes the temporarily stalled shock further out, but bootstraps a concomitant increase in the deposited neutrino power. Both effects together lie at the core of the resolution dependence we observe.


title: Holographic Entropy Cone with Time Dependence in Two Dimensions

abstract: In holographic duality, if a boundary state has a geometric description that realizes the Ryu-Takayanagi proposal then its entanglement entropies must obey certain inequalities that together define the so-called holographic entropy cone. A large family of such inequalities have been proven under the assumption that the bulk geometry is static, using a method involving contraction maps. By using kinematic space techniques, we show that in two boundary (three bulk) dimensions, all entropy inequalities that can be proven in the static case by contraction maps must also hold in holographic states with time dependence.


title: Strange quark suppression from a simultaneous Monte Carlo analysis of parton distributions and fragmentation functions

abstract: We perform the first simultaneous extraction of unpolarized parton distributions and fragmentation functions from a Monte Carlo analysis of inclusive and semi-inclusive deep-inelastic scattering, Drell-Yan lepton-pair production, and single-inclusive $e^+ e^-$ annihilation data. We use data resampling techniques to thoroughly explore the Bayesian posterior distribution of the extracted functions, and use $k$-means clustering on the parameter samples to identify the configurations that give the best description across all reactions. Inclusion of the semi-inclusive data reveals a strong suppression of the strange quark distribution at parton momentum fractions $x \gtrsim 0.01$, in contrast with the ATLAS observation of enhanced strangeness in $W^\pm$ and $Z$ production at the LHC. Our study reveals significant correlations between the strange quark density and the strange $\to$ kaon fragmentation function needed to simultaneously describe semi-inclusive $K^\pm$ production data from COMPASS and inclusive $K^\pm$ spectra in $e^+ e^-$ annihilation from ALEPH and SLD, as well as between the strange and light antiquark densities in the proton.


title: Combined explanations of $(g-2)_\mu$, $(g-2)_e$ and implications for a large muon EDM

abstract: We consider possible beyond-the-Standard-Model (BSM) effects that can accommodate both the long-standing tension in the anomalous magnetic moment of the muon, $a_\mu=(g-2)_\mu/2$, as well as the emerging $2.5\sigma$ deviation in its electron counterpart, $a_e=(g-2)_e/2$. After performing an EFT analysis, we consider BSM physics realized above the electroweak scale and find that a simultaneous explanation becomes possible in models with chiral enhancement. However, this requires a decoupling of the muon and electron BSM sectors to avoid the strong constraints from $\mu\to e\gamma$. In particular, this decoupling implies that there is no reason to expect the muon electric dipole moment (EDM) $d_\mu$ to be correlated with the electron EDM $d_e$, avoiding the very stringent limits for the latter. While some of the parameter space for $d_\mu$ favored by $a_\mu$ could be tested at the $(g-2)_\mu$ experiments at Fermilab and J-PARC, a dedicated muon EDM experiment at PSI would be able to probe most of this region.


title: Dark Matter Signatures of Supermassive Black Hole Binaries

abstract: A natural consequence of the galaxy formation paradigm is the existence of supermassive black hole (SMBH) binaries. Gravitational perturbations from a far-away SMBH companion can induce high orbital eccentricities on dark matter particles orbiting the primary SMBH via the eccentric Kozai-Lidov mechanism. This process yields an influx of dark matter particles into the primary SMBH ergosphere, where test particles linger for long timescales. This influx results in high self-gravitating densities, forming a dark matter clump that is extremely close to the SMBH. In such a situation, the gravitational wave emission between the dark matter clump and the SMBH is potentially detectable by LISA. If dark matter self-annihilates, the high densities of the clump will result in a unique co-detection of gravitational wave emission and high energy electromagnetic signatures.


title: Para-Hermitian Geometries for Poisson-Lie Symmetric $\sigma$-models

abstract: The doubled target space of the fundamental closed string is identified with its phase space and described by an almost para-Hermitian geometry. We explore this setup in the context of group manifolds which admit a maximally isotropic subgroup. This leads to a formulation of the Poisson-Lie $\sigma$-model and Poisson-Lie T-duality in terms of para-Hermitian geometry. The emphasis is put on so called half-integrable setups where only one of the Lagrangian subspaces of the doubled space has to be integrable. Using the dressing coset construction in Poisson-Lie T-duality, we extend our construction to more general coset spaces. This allows to explicitly obtain a huge class of para-Hermitian geometries. Each of them is automatically equipped which a generalized frame field, required for consistent generalized Scherk-Schwarz reductions. As examples we present integrable $\lambda$- and $\eta$-deformations on the three- and two-sphere.


title: Precise predictions for multi-${\rm TeV}$ and ${\rm PeV}$ energy neutrino scattering rates

abstract: The scattering rate of multi-${\rm TeV}$ and ${\rm PeV}$ energy neutrinos is fast becoming an interesting topic in (astro)particle-physics. This is due to experimental progress at Neutrino Telescopes such as IceCube which have begun to gain sensitivity to the flux of neutrinos in this energy range. In view of this, a precise calculation of the scattering rate of neutrinos upon atoms is presented. The two main components of the calculation are the differential cross-section predictions for neutrino scattering upon an atomic nucleus (such as that which make up water), as well as upon atomic electrons. In the first case, the predictions for neutrino-nucleon cross-sections in charged- and neutral-current scattering are refined by including resonant contributions generated within the photon field of the nucleus, which alter the considered distributions by up to $\approx 4\%$. In the latter case, radiative corrections are provided for all $2\to2$ scattering processes of the form $\bar{\nu}_e e^-\to f\bar{f}^{\prime}$. For antineutrino energies of $E_{\bar{\nu}_e}\approx6{\rm~PeV}$, where these processes become resonantly enhanced (the Glashow resonance) and dominate the total cross-section, these corrections amount to $\approx-10\%$.


title: Finding the remnants of the Milky Way's last neutron star mergers

abstract: The discovery of a binary neutron star merger (NSM) through both its gravitational wave and electromagnetic emission has revealed these events to be key sites of r-process nucleosynthesis. Here, we evaluate the prospects of finding the remnants of Galactic NSMs by detecting the gamma-ray decay lines from their radioactive r-process ejecta. We find that $^{126}$Sn, which has several lines in the energy range 415-695 keV and resides close to the second r-process peak, is the most promising isotope, because of its half-life $t_{1/2}=2.30(14)\times 10^{5}$ yr being comparable to the ages of recent NSMs. Using a Monte Carlo procedure, we predict that multiple remnants are detectable as individual sources by next-generation gamma-ray telescopes which achieve sub-MeV line sensitivities of $\sim 10^{-8}$-$10^{-6}$ $\gamma$ cm$^{-2}$ s$^{-1}$. However, given the unknown locations of the remnants, the most promising search strategy is a systematic survey of the Galactic plane and bulge extending to high Galactic latitudes. Individual known supernova remnants which may be mis-classified NSM remnants could also be targeted, especially those located outside the Galactic plane. Detection of a moderate sample of Galactic NSM remnants would provide important clues to unresolved issues such as the production of actinides in NSMs, properties of merging NS binaries, and even help distinguish them from rare supernovae as current Galactic r-process sources. We also investigate the diffuse flux from longer-lived nuclei (e.g. $^{182}$Hf) that could in principle trace the Galactic spatial distribution of NSMs over longer timescales, but find that the detection of the diffuse flux appears challenging even with next-generation telescopes.


title: Novel implementation of the multipole expansion to quarkonium hadronic transitions

abstract: We compute hadronic transitions between heavy quarkonium states with two, or one, pion/eta particles in the final state. We use the multipole expansion but not the twist expansion. The latter cannot be justified for the energy release of hadronic transitions between heavy quarkonium states with different principal quantum numbers. Instead, we use a counting based on the dimension of the interpolating field of the hybrid. This alternative counting allows us to still use chiral low-energy theorems to compute the pion production by local gluonic operators. We explore the phenomenological impact of this counting. Remarkably enough, for the two-pion transitions, we obtain the same predictions for the normalized differential decay rate as those obtained assuming the twist expansion. We implement this computational scheme using the hadronic representation of the effective theory potential NRQCD. We assume that the inverse Bohr radius of the heavy quarkonium is much larger than $\Lambda_{\rm QCD}$ but do not impose any constraint on the relative size of $\Lambda_{\rm QCD}$ and the typical kinetic energy of the bound state.


title: Argyres-Douglas theories and Liouville Irregular States

abstract: We study irregular states of rank-two and three in Liouville theory, based on an ansatz proposed by D. Gaiotto and J. Teschner. Using these irregular states, we evaluate asymptotic expansions of irregular conformal blocks corresponding to the partition functions of $(A_1,A_3)$ and $(A_1,D_4)$ Argyres-Douglas theories for general $\Omega$-background parameters. In the limit of vanishing Liouville charge, our result reproduces strong coupling expansions of the partition functions recently obtained via the Painlev\'e/gauge correspondence. This suggests that the irregular conformal block for one irregular singularity of rank 3 on sphere is also related to Painlev\'e II. We also find that our partition functions are invariant under the action of the Weyl group of flavor symmetries once four and two-dimensional parameters are correctly identified. We finally propose a generalization of this parameter identification to general irregular states of integer rank.


title: Prospects of Cosmic Superstring Detection through Microlensing of Extragalactic Point-Like Sources

abstract: The existence of cosmic superstrings may be probed by astronomical time domain surveys. When crossing the line of sight to point-like sources, strings produce a distinctive microlensing signature. We consider two avenues to hunt for a relic population of superstring loops: frequent monitoring of (1) stars in Andromeda, lensed by loops in the haloes of the Milky-Way and Andromeda and (2) supernovae at cosmological distances, lensed by loops in the intergalactic medium. We assess the potential of such experiments to detect and/or constrain strings with a range of tensions, $10^{-15} \lesssim G \mu/c^2 \lesssim 10^{-6}$. The practical sensitivity is tied to cadence of observations which we explore in detail. We forecast that high-cadence monitoring of $\sim 10^5$ stars on the far side of Andromeda over a year-long period will detect microlensing events if $G\mu/c^2 \sim 10^{-13}$, while $\sim 10^6$ stars will detect events if $10^{-13.5} <G\mu/c^2 < 10^{-11.5}$; the upper and lower bounds of the accessible tension range continue to expand as the number of stars rises. We also analyse the ability to reject models in the absence of fluctuations. While challenging, these studies are within reach of forthcoming time-domain surveys. Supernova observations can hypothetically constrain models with $10^{-12} < G\mu/c^2 < 10^{-6}$ without any optimisation of the survey cadence. However, the event rate forecast suggests it will be difficult to reject models of interest. As a demonstration, we use observations from the Pantheon Type Ia supernova cosmology data-set to place modest constraints on the number density of cosmic superstrings in a poorly tested region of the parameter space.


title: $N_\infty$-operads and associahedra

abstract: We provide a new combinatorial approach to studying the collection of N-infinity-operads in G-equivariant homotopy theory for G a finite cyclic group. In particular, we show that for G the cyclic group of order p^n the natural order on the collection of N-infinity-operads stands in bijection with the poset structure of the (n+1)-associahedron. We further provide a lower bound for the number of possible N-infinity-operads for any finite cyclic group G.


title: Towards the renormalization group flow of Horava gravity in $(3+1)$ dimensions

abstract: We compute the renormalization group running of the Newton constant and the parameter $\lambda$ in $(3+1)$-dimensional projectable Horava gravity. We use the background field method expanding around configurations with flat spatial metric, but non-vanishing shift. This allows us to reduce the number of interaction vertices and thereby drastically simplify the calculations. The gauge invariant $\beta$-function of $\lambda$ has two families of zeros, attractive in the infrared and ultraviolet respectively. They are candidates for the fixed points of the full renormalization group flow of the theory, once the $\beta$-functions for the rest of the couplings are added.


title: Comparaison between Coulomb and Hulth\`en potentials within Bohr Hamiltonian for $\gamma$-rigid nuclei in the presence of minimal length

abstract: In this work we solve the Schr\"odinger equation for Bohr Hamiltonian with Coulomb and Hulth\'en potentials within the formalism of minimal length in order to obtain analytical expressions for the energy eigenvalues and eigenfunctions by means of asymptotic iteration method. The obtained formulas of the energy spectrum and wave functions, are used to calculate excitation energies and transition rates of $\gamma$-rigid nuclei and compared with the experimental data at the shape phase critical point X(3) in nuclei.


title: On the continuous dual of the sequence space $bv$

abstract: Imaninezhad and Miri introduced the sequence space $ d_{\infty} $ in order to characterize the continuous dual of the sequence space $ bv. $ We show by a counterexample that this claimed characterization is false.


title: Toward a Complete Understanding of the Magellanic Stream Formation

abstract: The Magellanic Clouds have lost most of their gas during their passage by the Milky Way, a property that has never been successfully modeled. Here we use accurate and mesh-free hydrodynamic simulations to reproduce the Magellanic Stream and the Magellanic Clouds in the frame of a 'ram-pressure plus collision' model. This model reproduces many of the observed properties of the HI Stream including most of its density profile along its length and its dual filamentary structure. Besides this, ram-pressure combined with Kelvin-Helmholtz instabilities extracts amounts of ionized and HI gas consistent with those observed. The modeled scenario also reproduces the Magellanic Bridge, including the offset between young and old stars, and the collision between the Clouds, which is responsible of the very elongated morphology of the Small Magellanic Cloud along the line of sight. This model has solved most of the mysteries linked to the formation of the Magellanic Stream. The Leading Arm is not reproduced in the current model because it requires an alternative origin.


title: Titan's Dynamic Love Number Implies Stably-Stratified Ocean

abstract: The dynamic quadrupole Love number of Titan measured by \Cassini is $k_\mathrm{2,obs}=0.616\pm 0.067$, strongly indicating a global subsurface ocean. However, the theoretical Love number due to equilibrium tides is at most $k_\mathrm{2,eq}^\mathrm{max}\approx 0.48$ in the absence of an ice shell on top of the ocean. In reality, there is an outer ice shell of thickness $ 100\,\mathrm{km}$, reducing the equilibrium-tide Love number to $k_\mathrm{2,eq}\approx 0.42$. Therefore, other types of tidal response, like dynamic tides, may be also present in Titan. We propose that the ocean is stably stratified. As a result, there exist standing ocean waves (gravity modes) with eigen-frequencies close to the tidal frequency. Such a gravity mode (g-mode) is resonantly excited. It bends the outer ice shell radially and thus enhances the dynamic Love number by $k_\mathrm{2,g}$. In order for $k_\mathrm{2,g}$ to account for the discrepancy between $k_\mathrm{2,eq}$ and $k_\mathrm{2,obs}$, the Brunt-Vaisala frequency in the ocean is required to be $3.3\times 10^{-4}\,\mathrm{rad\, s^{-1}}$. It is compatible with the volatile-rich model for Titan that was proposed to explain the methane-rich atmosphere. The three components of the tidal potential with azimuthal degrees, $m=-2,0,2$, correspond to the three components of the quadrupole Love number, $k_\mathrm{2,-2}$, $k_\mathrm{2,0}$ and $k_\mathrm{2,2}$. They can excite retrograde, axisymmetric and prograde g-modes equally in the absence of rotation. However, Coriolis force induced by Titan's rotation breaks the symmetry among these modes. Most likely, only one of the Love-number components is significantly enhanced by a g-mode, while the other two are still attributed to equilibrium tides. This prediction is testable by observation. If confirmed, the smaller components of the Love number can be used to constrain the thickness of the outer ice shell.


title: Single site factors of Gibbs measures

abstract: It has been an open problem to identify classes of Gibbs measures less regular then H\"older continuous on the full shift which are closed under factor maps. In this article we show that in fact all of the classical uniqueness regimes (Bowen, Walters, and H\"older) from thermodynamic formalism are closed under factor maps between full shifts. In fact we show more generally that the classical uniqueness regimes are closed under factors between shifts of finite type provided the factor map satisfies a suitable mixing in fibers condition.


title: Entropy Production in Open Systems: The Predominant Role of Intra-Environment Correlations

abstract: We show that the entropy production in small open systems coupled to environments made of extended baths is predominantly caused by the displacement of the environment from equilibrium rather than, as often assumed, the mutual information between the system and the environment. The latter contribution is strongly bounded from above by the Araki-Lieb inequality, and therefore is not time-extensive, in contrast to the entropy production itself. We confirm our results with exact numerical calculations of the system-environment dynamics.


title: Real K\"ahler Submanifolds in Codimension $6$

abstract: We show that a real K\"ahler submanifold in codimension $6$ is essentially a holomorphic submanifold of another real K\"ahler submanifold in lower codimension if the second fundamental form is not sufficiently degenerated. We also give a shorter proof of this result when the real K\"ahler submanifold is minimal, using recent results about isometric rigidity.


title: Think Globally, Act Locally: A Deep Neural Network Approach to High-Dimensional Time Series Forecasting

abstract: Forecasting high-dimensional time series plays a crucial role in many applications such as demand forecasting and financial predictions. Modern datasets can have millions of correlated time-series that evolve together, i.e they are extremely high dimensional (one dimension for each individual time-series). There is a need for exploiting global patterns and coupling them with local calibration for better prediction. However, most recent deep learning approaches in the literature are one-dimensional, i.e, even though they are trained on the whole dataset, during prediction, the future forecast for a single dimension mainly depends on past values from the same dimension. In this paper, we seek to correct this deficiency and propose DeepGLO, a deep forecasting model which thinks globally and acts locally. In particular, DeepGLO is a hybrid model that combines a global matrix factorization model regularized by a temporal convolution network, along with another temporal network that can capture local properties of each time-series and associated covariates. Our model can be trained effectively on high-dimensional but diverse time series, where different time series can have vastly different scales, without a priori normalization or rescaling. Empirical results demonstrate that DeepGLO can outperform state-of-the-art approaches; for example, we see more than 25% improvement in WAPE over other methods on a public dataset that contains more than 100K-dimensional time series.


title: Topological pumping of quantum correlations

abstract: Topological pumping and duality transformations are paradigmatic concepts in condensed matter and statistical mechanics. In this paper, we extend the concept of topological pumping of particles to topological pumping of quantum correlations. We propose a scheme to find pumping protocols for highly-correlated states by mapping them to uncorrelated ones. We show that one way to achieve this is to use dualities, because they are non-local transformations that preserve the topological properties of the system. By using them, we demonstrate that topological pumping of kinks and cluster-like excitations can be realized. We find that the entanglement of these highly-correlated excitations is strongly modified during the pumping process and the interactions enhance the robustness against disorder. Our work paves the way to explore topological pumping beyond the notion of particles and opens a new avenue to investigate the relation between correlations and topology.

categories: quant-ph cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el


title: Maximum A Posteriori Probability (MAP) Joint Fine Frequency Offset and Channel Estimation for MIMO Systems with Channels of Arbitrary Correlation

abstract: Channel and frequency offset estimation is a classic topic with a large body of prior work using mainly maximum likelihood (ML) approach together with Cram\'er-Rao Lower bounds (CRLB) analysis. We provide the maximum a posteriori (MAP) estimation solution which is particularly useful for for tracking where previous estimation can be used as prior knowledge. Unlike the ML cases, the corresponding Bayesian Cram\'er-Rao Lower bound (BCRLB) shows clear relation with parameters and a low complexity algorithm achieves the BCRLB in almost all SNR range. We allow the time invariant channel within a packet to have arbitrary correlation and mean. The estimation is based on pilot/training signals. An unexpected result is that the joint MAP estimation is equivalent to an individual MAP estimation of the frequency offset first, again different from the ML results. We provide insight on the pilot/training signal design based on the BCRLB. Unlike past algorithms that trade performance and/or complexity for the accommodation of time varying channels, the MAP solution provides a different route for dealing with time variation. Within a short enough (segment of) packet where the channel and CFO are approximately time invariant, the low complexity algorithm can be employed. Similar to belief propagation, the estimation of the previous (segment of) packet can serve as the prior knowledge for the next (segment of) packet.


title: Wearable Sensor Data Based Human Activity Recognition using Machine Learning: A new approach

abstract: Recent years have witnessed the rapid development of human activity recognition (HAR) based on wearable sensor data. One can find many practical applications in this area, especially in the field of health care. Many machine learning algorithms such as Decision Trees, Support Vector Machine, Naive Bayes, K-Nearest Neighbor, and Multilayer Perceptron are successfully used in HAR. Although these methods are fast and easy for implementation, they still have some limitations due to poor performance in a number of situations. In this paper, we propose a novel method based on the ensemble learning to boost the performance of these machine learning methods for HAR.


title: Electromechanical Brillouin scattering in integrated optomechanical waveguides

abstract: In the well-known stimulated Brillouin scattering (SBS) process, spontaneous acoustic phonons in materials are stimulated by laser light and scatter the latter into a Stokes sideband. SBS becomes more pronounced in optical fibers and has been harnessed to amplify optical signals and even achieve lasing. Exploitation of SBS has recently surged on integrated photonics platforms as simultaneous confinement of photons and phonons in waveguides leads to drastically enhanced interaction. Instead of being optically stimulated, coherent phonons can also be electromechanically excited with very high efficiency as has been exploited in radiofrequency acoustic filters. Here, we demonstrate electromechanically excited Brillouin scattering in integrated optomechanical waveguides made of piezoelectric material aluminum nitride (AlN). Acoustic phonons of 16 GHz in frequency are excited with nanofabricated electromechanical transducers to scatter counter-propagating photons in the waveguide into a single anti-Stokes sideband. We show that phase-matching conditions of Brillouin scattering can be tuned by varying both the optical wavelength and the acoustic frequency to realize tunable single-sideband modulation. Combining Brillouin scattering photonics with nanoelectromechanical systems, our approach provides an efficient interface between microwave and optical photons that will be important for microwave photonics and potentially quantum transduction.


title: Long rainbow arithmetic progressions

abstract: Define $T_k$ as the minimal $t\in \mathbb{N}$ for which there is a rainbow arithmetic progression of length $k$ in every equinumerous $t$-coloring of $[tn]$ for all $n\in \mathbb{N}$. Jungi\'{c}, Licht (Fox), Mahdian, Nesetril and Radoici\'{c} proved that $\lfloor{\frac{k^2}{4}\rfloor}\le T_k$. We almost close the gap between the upper and lower bounds by proving that $T_k \le k^2e^{(\ln\ln k)^2(1+o(1))}$. Conlon, Fox and Sudakov have independently shown a stronger statement that $T_k=O(k^2\log k)$.


title: Convolutional Neural Networks Utilizing Multifunctional Spin-Hall MTJ Neurons

abstract: We propose a new network architecture for standard spin-Hall magnetic tunnel junction-based spintronic neurons that allows them to compute multiple critical convolutional neural network functionalities simultaneously and in parallel, saving space and time. An approximation to the Rectified Linear Unit transfer function and the local pooling function are computed simultaneously with the convolution operation itself. A proof-of-concept simulation is performed on the MNIST dataset, achieving up to 98% accuracy at a cost of less than 1 nJ for all convolution, activation and pooling operations combined. The simulations are remarkably robust to thermal noise, performing well even with very small magnetic layers.


title: When Deep Learning Met Code Search

abstract: There have been multiple recent proposals on using deep neural networks for code search using natural language. Common across these proposals is the idea of $\mathit{embedding}$ code and natural language queries, into real vectors and then using vector distance to approximate semantic correlation between code and the query. Multiple approaches exist for learning these embeddings, including $\mathit{unsupervised}$ techniques, which rely only on a corpus of code examples, and $\mathit{supervised}$ techniques, which use an $\mathit{aligned}$ corpus of paired code and natural language descriptions. The goal of this supervision is to produce embeddings that are more similar for a query and the corresponding desired code snippet. Clearly, there are choices in whether to use supervised techniques at all, and if one does, what sort of network and training to use for supervision. This paper is the first to evaluate these choices systematically. To this end, we assembled implementations of state-of-the-art techniques to run on a common platform, training and evaluation corpora. To explore the design space in network complexity, we also introduced a new design point that is a $\mathit{minimal}$ supervision extension to an existing unsupervised technique. Our evaluation shows that: 1. adding supervision to an existing unsupervised technique can improve performance, though not necessarily by much; 2. simple networks for supervision can be more effective that more sophisticated sequence-based networks for code search; 3. while it is common to use docstrings to carry out supervision, there is a sizeable gap between the effectiveness of docstrings and a more query-appropriate supervision corpus. The evaluation dataset is now available at arXiv:1908.09804


title: Non-Asymptotic Gap-Dependent Regret Bounds for Tabular MDPs

abstract: This paper establishes that optimistic algorithms attain gap-dependent and non-asymptotic logarithmic regret for episodic MDPs. In contrast to prior work, our bounds do not suffer a dependence on diameter-like quantities or ergodicity, and smoothly interpolate between the gap dependent logarithmic-regret, and the $\widetilde{\mathcal{O}}(\sqrt{HSAT})$-minimax rate. The key technique in our analysis is a novel "clipped" regret decomposition which applies to a broad family of recent optimistic algorithms for episodic MDPs.


title: Exploiting the violation of Lorentz symmetry for the planar Hall effect

abstract: The low energy Dirac and Weyl spectra are allowed to violate the Lorentz symmetry and thereby have a tilted energy dispersion. The tilt in the energy dispersion induces a Hall voltage in the plane spanned by the electric field and the tilt velocity. In the presence of a magnetic field the planar Hall conductivity and resistivity show Shubnikov de Haas oscillations. The oscillations in the planar Hall effect can become a fingerprint to spot the anomalous transport in Dirac and Weyl semimetals.

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


title: Enhancing second harmonic generation using dipolar-parity modes in non-planar plasmonic nanocavities

abstract: There is an active demand to develop efficient nanoscale nonlinear sources for applications in photonic circuitry, quantum optics and biosensing. To this end, plasmonic systems have been utilized to boost the nonlinear signal generation, however high efficiencies of frequency conversion for realistic applications remain a challenge. Metal-insulator-metal (MIM) nanocavities are good candidates for strongly concentrating the fields at the nanoscale to enhance the optical nonlinearities, however they typically suffer from the requirement to have a quadrupolar resonance at the emission wavelength. Here, we introduce nonplanar MIM nanocavities with a nonlinear spacer that can strongly enhance the second harmonic generation (SHG) despite of having fundamental and emission modes of the same parity. Our experimental and numerical results indicate that the enhancement is due to the non-planar design of the cavities and the bulk nonlinearity of the spacer layer.


title: On the Linear Speedup Analysis of Communication Efficient Momentum SGD for Distributed Non-Convex Optimization

abstract: Recent developments on large-scale distributed machine learning applications, e.g., deep neural networks, benefit enormously from the advances in distributed non-convex optimization techniques, e.g., distributed Stochastic Gradient Descent (SGD). A series of recent works study the linear speedup property of distributed SGD variants with reduced communication. The linear speedup property enable us to scale out the computing capability by adding more computing nodes into our system. The reduced communication complexity is desirable since communication overhead is often the performance bottleneck in distributed systems. Recently, momentum methods are more and more widely adopted in training machine learning models and can often converge faster and generalize better. For example, many practitioners use distributed SGD with momentum to train deep neural networks with big data. However, it remains unclear whether any distributed momentum SGD possesses the same linear speedup property as distributed SGD and has reduced communication complexity. This paper fills the gap by considering a distributed communication efficient momentum SGD method and proving its linear speedup property.


title: Beta Survival Models

abstract: This article analyzes the problem of estimating the time until an event occurs, also known as survival modeling. We observe through substantial experiments on large real-world datasets and use-cases that populations are largely heterogeneous. Sub-populations have different mean and variance in their survival rates requiring flexible models that capture heterogeneity. We leverage a classical extension of the logistic function into the survival setting to characterize unobserved heterogeneity using the beta distribution. This yields insights into the geometry of the problem as well as efficient estimation methods for linear, tree and neural network models that adjust the beta distribution based on observed covariates. We also show that the additional information captured by the beta distribution leads to interesting ranking implications as we determine who is most-at-risk. We show theoretically that the ranking is variable as we forecast forward in time and prove that pairwise comparisons of survival remain transitive. Empirical results using large-scale datasets across two use-cases (online conversions and retention modeling), demonstrate the competitiveness of the method. The simplicity of the method and its ability to capture skew in the data makes it a viable alternative to standard techniques particularly when we are interested in the time to event and when the underlying probabilities are heterogeneous.