title: Digital Twin Enabled Smart Control Engineering as an Industrial AI: A New Framework and A Case Study

abstract: In the way towards Industry 4.0, the complexity of the industrial systems increases due to the presence of multiple agents, Cyber-Physical Systems, distributed sensing, and big data introducing unknown dynamics that affect the production goals of the manufacturing processes. Thus, Digital Twin is a breaking technology corresponding to the capacity of developing a virtual representation of any complex system in order to perform design, analysis, and behavior prediction tasks that enhance the understanding of these systems through new enabling capabilities like real-time analytics, parallel sensing, or Smart Control Engineering. In this paper, a novel framework is proposed for the design and implementation of Digital Twin applications to the development of Smart Control Engineering. The steps of this framework involve system documentation, multidomain simulation, behavioral matching, and real-time monitoring. This framework is applied to develop the Digital Twin for a real-time vision feedback infrared temperature uniformity control. The obtained results show that Digital Twin is a fundamental part of the transformation into Industry 4.0.


title: GPU-Accelerated Drug Discovery with Docking on the Summit Supercomputer: Porting, Optimization, and Application to COVID-19 Research

abstract: Protein-ligand docking is an in silico tool used to screen potential drug compounds for their ability to bind to a given protein receptor within a drug-discovery campaign. Experimental drug screening is expensive and time consuming, and it is desirable to carry out large scale docking calculations in a high-throughput manner to narrow the experimental search space. Few of the existing computational docking tools were designed with high performance computing in mind. Therefore, optimizations to maximize use of high-performance computational resources available at leadership-class computing facilities enables these facilities to be leveraged for drug discovery. Here we present the porting, optimization, and validation of the AutoDock-GPU program for the Summit supercomputer, and its application to initial compound screening efforts to target proteins of the SARS-CoV-2 virus responsible for the current COVID-19 pandemic.


title: Model independent determination of the spin of the $\Omega^{-}$ and its polarization alignment in $\psi(3686)\rightarrow\Omega^{-}\bar{\Omega}^{+}$

abstract: We present an analysis of the process $\psi(3686) \to \Omega^- \bar{\Omega}^+$ ($\Omega^-\to K^-\Lambda$, $\bar{\Omega}^+\to K^+\bar{\Lambda}$, $\Lambda\to p\pi^-$, $\bar{\Lambda}\to \bar{p}\pi^+$) based on a data set of $448\times 10^6$ $\psi(3686)$ decays collected with the BESIII detector at the BEPCII electron-positron collider. The helicity amplitudes for the process $\psi(3686) \to \Omega^- \bar{\Omega}^+$ and the decay parameters of the subsequent decay $\Omega^-\to K^-\Lambda$ $(\bar{\Omega}^+\to K^+\bar{\Lambda})$ are measured for the first time by a fit to the angular distribution of the complete decay chain. The branching fraction of $\psi(3686) \to \Omega^- \bar{\Omega}^+$ is measured to be $(5.82\pm 0.12\pm 0.24)\times 10^{-5}$, with an improved precision compared to previous measurements.


title: DRIVE: A Digital Network Oracle for Cooperative Intelligent Transportation Systems

abstract: In a world where Artificial Intelligence revolutionizes inference, prediction and decision-making tasks, Digital Twins emerge as game-changing tools. A case in point is the development and optimization of Cooperative Intelligent Transportation Systems (C-ITSs): a confluence of cyber-physical digital infrastructure and (semi)automated mobility. Herein we introduce Digital Twin for self-dRiving Intelligent VEhicles (DRIVE). The developed framework tackles shortcomings of traditional vehicular and network simulators. It provides a flexible, modular, and scalable implementation to ensure large-scale, city-wide experimentation with a moderate computational cost. The defining feature of our Digital Twin is a unique architecture allowing for submission of sequential queries, to which the Digital Twin provides instantaneous responses with the "state of the world", and hence is an Oracle. With such bidirectional interaction with external intelligent agents and realistic mobility traces, DRIVE provides the environment for development, training and optimization of Machine Learning based C-ITS solutions.


title: Fast Bayesian Estimation of Spatial Count Data Models

abstract: Spatial count data models are used to explain and predict the frequency of phenomena such as traffic accidents in geographically distinct entities such as census tracts or road segments. These models are typically estimated using Bayesian Markov chain Monte Carlo (MCMC) simulation methods, which, however, are computationally expensive and do not scale well to large datasets. Variational Bayes (VB), a method from machine learning, addresses the shortcomings of MCMC by casting Bayesian estimation as an optimisation problem instead of a simulation problem. Considering all these advantages of VB, a VB method is derived for posterior inference in negative binomial models with unobserved parameter heterogeneity and spatial dependence. P\'olya-Gamma augmentation is used to deal with the non-conjugacy of the negative binomial likelihood and an integrated non-factorised specification of the variational distribution is adopted to capture posterior dependencies. The benefits of the proposed approach are demonstrated in a Monte Carlo study and an empirical application on estimating youth pedestrian injury counts in census tracts of New York City. The VB approach is around 45 to 50 times faster than MCMC on a regular eight-core processor in a simulation and an empirical study, while offering similar estimation and predictive accuracy. Conditional on the availability of computational resources, the embarrassingly parallel architecture of the proposed VB method can be exploited to further accelerate its estimation by up to 20 times.


title: A Dynamic Choice Model with Heterogeneous Decision Rules: Application in Estimating the User Cost of Rail Crowding

abstract: Crowding valuation of subway riders is an important input to various supply-side decisions of transit operators. The crowding cost perceived by a transit rider is generally estimated by capturing the trade-off that the rider makes between crowding and travel time while choosing a route. However, existing studies rely on static compensatory choice models and fail to account for inertia and the learning behaviour of riders. To address these challenges, we propose a new dynamic latent class model (DLCM) which (i) assigns riders to latent compensatory and inertia/habit classes based on different decision rules, (ii) enables transitions between these classes over time, and (iii) adopts instance-based learning theory to account for the learning behaviour of riders. We use the expectation-maximisation algorithm to estimate DLCM, and the most probable sequence of latent classes for each rider is retrieved using the Viterbi algorithm. The proposed DLCM can be applied in any choice context to capture the dynamics of decision rules used by a decision-maker. We demonstrate its practical advantages in estimating the crowding valuation of an Asian metro's riders. To calibrate the model, we recover the daily route preferences and in-vehicle crowding experiences of regular metro riders using a two-month-long smart card and vehicle location data. The results indicate that the average rider follows the compensatory rule on only 25.5% of route choice occasions. DLCM estimates also show an increase of 47% in metro riders' valuation of travel time under extremely crowded conditions relative to that under uncrowded conditions.


title: Resurgence in the Bi-Yang-Baxter Model

abstract: We study the integrable bi-Yang-Baxter deformation of the $SU(2)$ principal chiral model (PCM) and its finite action uniton solutions. Under an adiabatic compactification on an $S^1$, we obtain a quantum mechanics with an elliptic Lam\'e-like potential. We perform a perturbative calculation of the ground state energy in this quantum mechanics to large orders obtaining an asymptotic series. Using the Borel-Pad\'e technique, we determine the expected locations of branch cuts in the Borel plane of the perturbative series and show that they match the values of the uniton actions. Therefore, we can match the non-perturbative contributions to the energy with the uniton solutions which fractionate upon adiabatic compactification. An off-shoot of the WKB analysis, is to identify the quadratic differential of this deformed PCM with that of an $\mathcal{N}=2$ Seiberg-Witten theory. This can be done either as an $N_f=4$ $SU(2)$ theory or as an elliptic quiver $SU(2)\times SU(2)$ theory. The mass parameters of the gauge theory are given by the deformation parameters of the PCM.


title: On the spectral type of rank one flows and Banach problem with calculus of generalized Riesz products on the real line

abstract: It is shown that a certain class of Riesz product type measures on $\mathbb{R}$ is realized a spectral type of rank one flows. As a consequence, we will establish that some class of rank one flows has a singular spectrum. Some of the results presented here are even new for the $\mathbb{Z}$-action. Our method is based, on one hand, on the extension of Bourgain-Klemes-Reinhold-Peyri\`ere method, and on the other hand, on the extension of the Central Limit Theorem approach to the real line which gives a new extension of Salem-Zygmund Central Limit Theorem. We extended also a formula for Radon-Nikodym derivative between two generalized Riesz products obtained by el Abdalaoui-Nadkarni and a formula of Mahler measure of the spectral type of rank one flow but in the weak form. We further present an affirmative answer to the flow version of the Banach problem, and we discuss some issues related to flat trigonometric polynomials on the real line in connection with the famous Banach-Rhoklin problem in the spectral theory of dynamical systems.


title: Resistivity saturation in an electron-doped cuprate

abstract: We report the observation of resistivity saturation in lightly doped ($x\sim 0.10)$ as-grown samples of the electron-doped cuprate La$_{2-x}$Ce$_x$CuO$_4$ (LCCO). The saturation occurs at resistivity values roughly consistent with the phenomenological Mott-Ioffe-Regel criterion once the low effective carrier density of these materials is included in the analysis. These results imply that, at least for light doping, the high-temperature metallic phase of these materials is not necessarily strange and may be understood as simply a low-density metal.


title: Is GW170817 a Multimessenger Neutron Star-Primordial Black Hole Merger?

abstract: We investigate the possibility of the gravitational-wave event GW170817 being a light, solar-mass black hole (BH) - neutron star (NS) merger. We explore two exotic scenarios involving primordial black holes (PBH) that could produce such an event, taking into account available observational information on NGC 4993. First, we entertain the possibility of dynamical NS-PBH binary formation where a solar-mass PBH and a NS form a binary through gravitational interaction. We find that while dynamical NS-PBH formation could account for the GW170817 event, the rate is highly dependent on unknown density contrast factors and could potentially be affected by galaxy mergers. We also find that PBH-PBH binaries would likely have a larger merger rate, assuming the density contrast boost factor of an order similar to the NS-PBH case. These exotic merger formations could provide new channels to account for the volumetric rate of compact-object mergers reported by LIGO/Virgo. Secondly, we consider the case where one of the NS's in a binary NS system is imploded by a microscopic PBH. We find that the predicted rate for NS implosion into a BH is very small, at least for the specific environment of NGC 4993. We point out that similar existing (e.g. GW190425 and GW190814) and future observations will shed additional light on these scenarios.

categories: astro-ph.HE astro-ph.CO astro-ph.GA gr-qc hep-ph


title: The age-chemical abundance structure of the Galactic disc II: $\alpha$-dichotomy and thick disc formation

abstract: We extend our previous work on the age-chemical abundance structure of the Galactic outer disc to the inner disc (4 < r < 8 kpc) based on the SDSS/APOGEE survey. Different from the outer disc, the inner disc stars exhibit a clear bimodal distribution in the [Mg/Fe]-[Fe/H] plane. While a number of scenarios have been proposed in the literature, it remains challenging to recover this bimodal distribution with theoretical models. To this end, we present a chemical evolution model embedding a complex multi-phase inner disc formation scenario that matches the observed bimodal [Mg/Fe]-[Fe/H] distribution. In this scenario, the formation of the inner disc is dominated by two main starburst episodes 6 Gyr apart with secular, low-level star formation activity in between. In our model, the first starburst occurs at early cosmic times (t~1 Gyr) and the second one 6 Gyr later at a cosmic time of t~7 Gyr. Both these starburst episodes are associated with gas accretion events in our model, and are quenched rapidly. The first starburst leads to the formation of the high-$\alpha$ sequence, and the second starburst leads to the formation of the metal-poor low-$\alpha$ sequence. The metal-rich low-$\alpha$ stars, instead, form during the secular evolution phase between the two bursts. Our model shows that the $\alpha$-dichotomy originates from the rapid suppression of star formation after the first starburst. The two starburst episodes are likely to be responsible for the formation of the geometric thick disc (z >1 kpc), with the old inner thick disc and the young outer thick disc forming during the first and the second starbursts, respectively.


title: Very wide companion fraction from Gaia DR2: a weak or no enhancement for hot jupiter hosts, and a strong enhancement for contact binaries

abstract: There is an ongoing debate on whether hot jupiter hosts are more likely to be found in wide binaries with separations of $\gtrsim 100$ AU. In this paper, we search for comoving, very wide companions with separations of $10^3-10^4$ AU for hot jupiter hosts and main-sequence contact binaries in Gaia DR2, and compare the very wide companion fractions with their object-by-object-matched field star samples. We find that $11.9\pm 2.5$% of hot jupiter hosts and $14.1\pm 1.0$% of contact binaries have companions at separations of $10^3-10^4$ AU. While the very wide companion fraction of hot jupiter hosts is a factor of $1.9\pm0.5$ larger than their matched field star sample, it is consistent, within $\sim1\sigma$, with that of matched field stars if the matching is only with field stars without close companions (within $\sim50$ AU) as is the case for hot jupiter hosts. The very wide companion fraction of contact binaries is a factor of $3.1\pm0.5$ larger than their matched field star sample, suggesting that the formation and evolution of contact binaries are either tied to or correlated with the presence of wide companions. In contrast, the weak enhancement of very wide companion fraction for hot jupiter hosts implies that the formation of hot jupiters is not as sensitive to those environment properties. Our results also hint that the occurrence rates of dual hot jupiter hosts and dual contact binaries may be higher than the expected values from random pairing of field stars, which may be due to their underlying metallicity and age dependence.


title: Microscopic bath effects on noise spectra in semiconductor quantum dot qubits

abstract: When a system is thermally coupled to only a small part of a larger bath, statistical fluctuations of the temperature (more precisely, the internal energy) of this "sub-bath" around the mean temperature defined by the larger bath can become significant. We show that these temperature fluctuations generally give rise to 1/f-like noise power spectral density from even a single two-level system. We extend these results to a distribution of fluctuators, finding the corresponding modification to the Dutta-Horn relation. Then we consider the specific situation of charge noise in silicon quantum dot qubits and show that recent experimental data [E. J. Connors, et al., Phys. Rev. B 100, 165305 (2019)] can be modeled as arising from as few as two two-level fluctuators, and accounting for sub-bath size improves the quality of the fit.


title: Frustration-induced anomalous transport and strong photon decay in waveguide QED

abstract: We study the propagation of photons in a one-dimensional environment consisting of two non-interacting species of photons frustratingly coupled to a single spin-1/2. The ultrastrong frustrated coupling leads to an extreme mixing of the light and matter degrees of freedom, resulting in the disintegration of the spin and a breakdown of the "dressed-spin", or polaron, description. Using a combination of numerical and analytical methods, we show that the elastic response becomes increasingly weak at the effective spin frequency, showing instead an increasingly strong and broadband response at higher energies. We also show that the photons can decay into multiple photons of smaller energies. The total probability of these inelastic processes can be as large as the total elastic scattering rate, or half of the total scattering rate, which is as large as it can be. The frustrated spin induces strong anisotropic photon-photon interactions that are dominated by inter-species interactions. Our results are relevant to state-of-the-art circuit and cavity quantum electrodynamics experiments.


title: Towards a multi-tracer timeline of star formation in the LMC -- I.\ Deriving the lifetimes of H\,{\sc i} clouds

abstract: The time-scales associated with the various stages of the star formation process remain poorly constrained. This includes the earliest phases of star formation, during which molecular clouds condense out of the atomic interstellar medium. We present the first in a series of papers with the ultimate goal of compiling the first multi-tracer timeline of star formation, through a comprehensive set of evolutionary phases from atomic gas clouds to unembedded young stellar populations. In this paper, we present an empirical determination of the lifetime of atomic clouds using the Uncertainty Principle for Star Formation formalism, based on the de-correlation of H$\alpha$ and H\,{\sc i} emission as a function of spatial scale. We find an atomic gas cloud lifetime of 48$\substack{+13\\-8}$\,Myr. This timescale is consistent with the predicted average atomic cloud lifetime in the LMC (based on galactic dynamics) that is dominated by the gravitational collapse of the mid-plane ISM. We also determine the overlap time-scale for which both H\,{\sc i} and H$\alpha$ emission are present to be very short ($t_{over}<1.7$\,Myr), consistent with zero, indicating that there is a near-to-complete phase change of the gas to a molecular form in an intermediary stage between H\,{\sc i} clouds and H\,{\sc ii} regions. We utilise the time-scales derived in this work to place empirically determined limits on the time-scale of molecular cloud formation. By performing the same analysis with and without the 30 Doradus region included, we find that the most extreme star forming environment in the LMC has little effect on the measured average atomic gas cloud lifetime. By measuring the lifetime of the atomic gas clouds, we place strong constraints on the physics that drives the formation of molecular clouds and establish a solid foundation for the development of a multi-tracer timeline of star formation in the LMC.


title: Sensitivities of KM3NeT on decaying dark matter

abstract: The discovery of high-energy astrophysical neutrinos by IceCube has opened a new window to the Universe. However, the origin of these neutrinos is still a mystery, and some of them could be a result of dark matter interactions such as decay. Next generation gigaton water-Cherenkov neutrino telescope, KM3NeT, is expected to offer significantly improved energy resolution in the cascade channel, and advantageous viewing condition to the Galactic Center; both important for searches of dark matter decay signals. We study the sensitivity of KM3NeT on dark matter decays by performing a mock likelihood analysis for both cascade and track type events, taking into account both angular and energy information. We find that, combining both channels, KM3NeT is expected to produce world leading limits on dark matter decay lifetime in the PeV mass range, and could test some of the dark matter hints in the current IceCube data.


title: Soft de Sitter Effective Theory

abstract: Calculating the quantum evolution of a de Sitter universe on superhorizon scales is notoriously difficult. To address this challenge, we introduce the Soft de Sitter Effective Theory (SdSET). This framework holds for superhorizon modes whose comoving momentum is far below the UV scale, which is set by the inverse comoving horizon. The SdSET is formulated using the same approach that yields the Heavy Quark Effective Theory. The degrees of freedom that capture the long wavelength dynamics are identified with the growing and decaying solutions to the equations of motion. The operator expansion is organized using a power counting scheme, and loops can be regulated while respecting the low energy symmetries. For massive quantum fields in a fixed de Sitter background, power counting implies that all interactions beyond the horizon are irrelevant. Alternatively, if the fields are very light, the leading interactions are at most marginal, and resumming the associated logarithms using (dynamical) renormalization group techniques yields the evolution equation for canonical stochastic inflation. The SdSET is also applicable to models where gravity is dynamical, including inflation. In this case, diffeomorphism invariance ensures that all interactions are irrelevant, trivially implying the all-orders conservation of adiabatic density fluctuations and gravitational waves. We briefly touch on the application to slow-roll eternal inflation by identifying novel relevant operators. This work serves to demystify many aspects of perturbation theory outside the horizon, and has a variety of applications to problems of cosmological interest.


title: Axion and neutrino bounds improved with new calibrations of the tip of the red-giant branch using geometric distance determinations

abstract: The brightness of the tip of the red-giant branch (TRGB) allows one to constrain novel energy losses that would lead to a larger core mass at helium ignition and thus to a brighter TRGB than expected by standard stellar models. The required absolute TRGB calibrations strongly improve with reliable geometric distances that have become available for the galaxy NGC 4258 that hosts a water megamaser and to the Large Magellanic Cloud based on 20 detached eclipsing binaries. Moreover, we revise a previous TRGB calibration in the globular cluster $\omega$ Centauri with a recent kinematical distance determination based on Gaia DR2 data. All of these calibrations have similar uncertainties and they agree with each other and with recent dedicated stellar models. Using NGC 4258 as the cleanest extra-galactic case, we thus find an updated constraint on the axion-electron coupling of $g_{ae}<1.6\times10^{-13}$ and $\mu_\nu<1.5\times10^{-12}\mu_{\rm B}$ (95\% CL) on a possible neutrino dipole moment, whereas $\omega$ Centauri as the best galactic target provides instead $g_{ae}<1.3\times10^{-13}$ and $\mu_\nu<1.2\times10^{-12}\mu_{\rm B}$. The reduced observational errors imply that stellar evolution theory and bolometric corrections begin to dominate the overall uncertainties.


title: The Past History of Galaxy Clusters told by their present neighbors

abstract: Galaxy clusters can play a key role in modern cosmology provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of clusters to their formation history are dominated by halo mass - environment relations. This paper presents a mass-free correlation between the present neighbor distribution of cluster-size halos and the latter mass assembly history. From the Big Multidark simulation, we extract two large samples of random halos with masses ranging from Virgo to Coma cluster sizes. Additionally, to find the main environmental culprit for the formation history of the Virgo cluster, we compare the Virgo-size halos to 200 Virgo-like halos extracted from simulations that resemble the local Universe. The number of neighbors at different cluster-centric distances permits discriminating between clusters with different mass accretion histories. Similarly to Virgo-like halos, clusters with numerous neighbors within a distance of about 2 times their virial radius experience a transition at z~1 between an active period of mass accretion, relative to the mean, and a quiet history. On the contrary, clusters with few neighbors share an opposite trend: from passive to active assembly histories. Additionally, clusters with massive companions within about 4 times their virial radius tend to have recent active merging histories. Therefore, the radial distribution of cluster neighbors provides invaluable insights into the past history of these objects.


title: Precise dark matter relic abundance in decoupled sectors

abstract: Dark matter (DM) as a thermal relic of the primordial plasma is increasingly pressured by direct and indirect searches, while the same production mechanism in a decoupled sector is much less constrained. We extend the standard treatment of the freeze-out process to such scenarios and perform precision calculations of the annihilation cross-section required to match the observed DM abundance. We demonstrate that the difference to the canonical value of this 'thermal cross-section' is generally sizeable, and can reach orders of magnitude. Our results directly impact the interpretation of DM searches in hidden sector scenarios.


title: A Comparison of the Stellar, CO and Dust-Continuum Emission from Three, Star-Forming HUDF Galaxies at $z\sim 2$

abstract: We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk.


title: The Morphology-Density relationship in 1<z<2 clusters

abstract: The morphology-density relationship states that dense cosmic environments such as galaxy clusters have an overabundance of quiescent elliptical galaxies, but it is unclear at which redshift this relationship is first established. We study the morphology of 4 clusters with $1.2<z<1.8$ using HST imaging and the morphology computation code statmorph. By comparing median morphology of cluster galaxies to CANDELS field galaxies using Monte Carlo analysis, we find that 2 out of 4 clusters (at z=1.19 and z=1.75) have an established morphology-density relationship with more than $3\sigma$ significance. $\sim$50% of galaxies in these clusters are bulge-dominated compared to $\sim$30% in the field, and they are significantly more compact. This result is more significant for low-mass galaxies with $\log M/M_\odot \lessapprox 10.5$, showing that low-mass galaxies are affected the most in clusters. We also find an intriguing system of two z $\approx$ 1.45 clusters at a unusually small separation 2D separation of $3'$ and 3D separation of $\approx73$ Mpc that exhibit no morphology-density relationship but have enhanced merger signatures. We conclude that the environmental mechanism responsible for the morphology-density relationship is 1) already active as early as z=1.75, 2) forms compact, bulge-dominated galaxies and 3) affects primarily low-mass galaxies. However, there is a significant degree of intracluster variance that may depend on the larger cosmological environment in which the cluster is embedded.


title: Parsimony in Model Selection: Tools for Assessing Fit Propensity

abstract: Theories can be represented as statistical models for empirical testing. There is a vast literature on model selection and multimodel inference that focuses on how to assess which statistical model, and therefore which theory, best fits the available data. For example, given some data, one can compare models on various information criterion or other fit statistics. However, what these indices fail to capture is the full range of counterfactuals. That is, some models may fit the given data better not because they represent a more correct theory, but simply because these models have more fit propensity - a tendency to fit a wider range of data, even nonsensical data, better. Current approaches fall short in considering the principle of parsimony (Occam's Razor), often equating it with the number of model parameters. Here we offer a toolkit for researchers to better study and understand parsimony through the fit propensity of Structural Equation Models. We provide an R package (ockhamSEM) built on the popular lavaan package. To illustrate the importance of evaluating fit propensity, we use ockhamSEM to investigate the factor structure of the Rosenberg Self-Esteem Scale.


title: Stirred and shaken: dynamical behavior of boson stars and dark matter cores

abstract: Bosonic fields can give rise to self-gravitating structures. These are interesting hypothetical new "dark matter stars" and good descriptions of dark matter haloes if the fields are very light. We study the dynamical response of Newtonian boson stars (NBS) when excited by external matter (stars, planets or black holes) in their vicinities. Our setup can describe the interaction between a massive black hole and the surrounding environment, shortly after the massive body has undergone a "kick", due to the collapse of baryonic matter at the galactic center, or dark matter depletion as a reaction to an inspiralling binary. We perform the first self-consistent calculation of dynamical friction acting on moving bodies in these backgrounds. Binaries close to coalescence "stir" the NBS core, and backreaction affects gravitational waveforms at leading $-6PN$ order with respect to the dominant quadrupolar term; the coefficient is too small to allow detection by next-generation interferometers. We also show that the gravitational collapse to a supermassive black hole at the center of a NBS is accompanied by only a small change in the surrounding core. The NBS eventually gets accreted, but for astrophysical parameters this occurs only after several Hubble times.


title: GeV-scale neutrinos: interactions with mesons and DUNE sensitivity

abstract: The simplest extension of the SM to account for the observed neutrino masses and mixings is the addition of at least two singlet fermions (or right-handed neutrinos). If their masses lie at or below the GeV scale, such new fermions would be produced in meson decays. Similarly, provided they are sufficiently heavy, their decay channels may involve mesons in the final state. Although the couplings between mesons and heavy neutrinos have been computed previously, significant discrepancies can be found in the literature. The aim of this paper is to clarify such discrepancies and provide consistent expressions for all relevant effective operators involving mesons with masses up to 2 GeV. Moreover, the effective Lagrangians obtained for both the Dirac and Majorana scenarios are made publicly available as FeynRules models so that fully differential event distributions can be easily simulated. As an application of our setup, we numerically compute the expected sensitivity of the DUNE near detector to these heavy neutral leptons.


title: The Role of Stellar Feedback in the Chemical Evolution of a Low Mass Dwarf Galaxy

abstract: We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust grains, ionization and associated heating, Lyman-Werner (LW) dissociation of H$_2$, and radiation pressure; and winds from massive main sequence (neglecting their energy input) and asymptotic giant branch (AGB) stars. Radiative transfer is done by ray tracing. We consider the effects each of these processes have on regulating the star formation rate, global properties, multi-phase interstellar medium (ISM), and driving of galactic winds. We follow individual metal species from distinct nucleosynthetic enrichment channels (AGB winds, massive star stellar winds, core collapse and Type Ia supernovae) and pay particular attention to how these feedback processes regulate metal mixing in the ISM, the metal content of outflows, and the stellar abundance patterns in our galaxy. We find that---for a low-metallicity, low-mass dwarf galaxy ---stellar radiation, particularly ionizing radiation and LW radiation, are important sources of stellar feedback whose effects dominate over photoelectric heating and HI radiation pressure. However, feedback is coupled non-linearly, and the inclusion or exclusion of each process produces non-negligible effects. We find strong variations with: the star formation history; the ejection fractions of metals, mass, and energy; and the distribution of elements from different nucleosynthetic sources in both the gas and stars.


title: Phenomenology of Magnetic Black Holes with Electroweak-Symmetric Coronas

abstract: Magnetically charged black holes (MBHs) are interesting solutions of the Standard Model and general relativity. They may possess a "hairy" electroweak-symmetric corona outside the event horizon, which speeds up their Hawking radiation and leads them to become nearly extremal on short timescales. Their masses could range from the Planck scale up to the Earth mass. We study various methods to search for primordially produced MBHs and estimate the upper limits on their abundance. We revisit the Parker bound on magnetic monopoles and show that it can be extended by several orders of magnitude using the large-scale coherent magnetic fields in Andromeda. This sets a mass-independent constraint that MBHs have an abundance less than $6 \times 10^{-3}$ times that of dark matter. MBHs can also be captured in astrophysical systems like the Sun, the Earth, or neutron stars. There, they can become non-extremal either from merging with an oppositely charged MBH or absorbing nucleons. The resulting Hawking radiation can be detected as neutrinos, photons, or heat. High-energy neutrino searches in particular can set a stronger bound than the Parker bound for some MBH masses, down to an abundance $10^{-7}$ of dark matter.


title: Computer-Aided Personalized Education

abstract: The shortage of people trained in STEM fields is becoming acute, and universities and colleges are straining to satisfy this demand. In the case of computer science, for instance, the number of US students taking introductory courses has grown three-fold in the past decade. Recently, massive open online courses (MOOCs) have been promoted as a way to ease this strain. This at best provides access to education. The bigger challenge though is coping with heterogeneous backgrounds of different students, retention, providing feedback, and assessment. Personalized education relying on computational tools can address this challenge. While automated tutoring has been studied at different times in different communities, recent advances in computing and education technology offer exciting opportunities to transform the manner in which students learn. In particular, at least three trends are significant. First, progress in logical reasoning, data analytics, and natural language processing has led to tutoring tools for automatic assessment, personalized instruction including targeted feedback, and adaptive content generation for a variety of subjects. Second, research in the science of learning and human-computer interaction is leading to a better understanding of how different students learn, when and what types of interventions are effective for different instructional goals, and how to measure the success of educational tools. Finally, the recent emergence of online education platforms, both in academia and industry, is leading to new opportunities for the development of a shared infrastructure. This CCC workshop brought together researchers developing educational tools based on technologies such as logical reasoning and machine learning with researchers in education, human-computer interaction, and cognitive psychology.


title: CMB constraints on a physical model of reionization

abstract: We study constraints on allowed reionization histories by comparing predictions of a physical semi-numerical model with secondary temperature and polarization anisotropies of the cosmic microwave background (CMB). Our model has four free parameters characterizing the evolution of ionizing efficiency $\zeta$ and the minimum mass $M_{\mathrm{min}}$ of haloes that can produce ionizing radiation. Comparing the model predictions with the presently available data of the optical depth $\tau$ and kinematic Sunyaev-Zeldovich signal, we find that we can already rule out a significant region of the parameter space. We limit the duration of reionization $\Delta z=1.30^{+0.19}_{-0.60}$ ($\Delta z < 2.9$ at $99\%$ C.L.), one of the tightest constraints on the parameter. The constraints mildly favour $M_{\mathrm{min}} \gtrsim 10^9 \mathrm{M}_{\odot}$ (at $68\%$ C.L.) at $z \sim 8$, thus indicating the presence of reionization feedback. Our analysis provides an upper bound on the secondary $B$-mode amplitude $D_{l=200}^{BB}<18$ nK$^2$ at $99\%$ C.L. We also study how the constraints can be further tightened with upcoming space and ground-based CMB missions. Our study, which relies solely on CMB data, has implications not only for upcoming CMB surveys for detecting primordial gravitational waves but also redshifted 21 cm studies.


title: Massive molecular gas reservoir around the central AGN in the CARLA J1103+3449 cluster at z=1.44

abstract: Passive early-type galaxies dominate cluster cores at z $\lesssim$1.5. At higher redshift, cluster core galaxies are observed to have still on-going star-formation, fuelled by cold molecular gas. We measure the molecular gas reservoir of the central region around the radio-loud AGN in the cluster CARLA J1103+3449 at z=1.44 with NOEMA. The AGN synchrotron emission dominates the continuum emission at 94.48 GHz, and we measure its flux at the AGN position and at the position of two radio jets. Combining our measurements with published results over the range 4.71 GHz-94.5 GHz, we obtain a flat spectral index $\alpha = 0.14 \pm 0.03$ for the AGN core emission, and a steeper index $\alpha = 1.43 \pm 0.04$ and $\alpha = 1.15 \pm 0.04$ at positions close to the western and eastern lobe, respectively. The total spectral index is $\alpha = 0.92 \pm 0.02$ over the range 73.8 MHz-94.5 GHz. We detect two CO(2-1) emission lines, both blue-shifted with respect to the AGN. Their emission corresponds to two regions, ~17 kpc south-east and ~14 kpc south-west of the AGN, not associated with galaxies. In these two regions, we find a total massive molecular gas reservoir of $M_{gas}$ = 3.9 $\pm$ 0.4 $10^{10} M_{\odot}$, which dominates (~ 60%) the central total molecular gas reservoir. These results can be explained by massive cool gas flows in the center of the cluster. The AGN early-type host is not yet quenched; its star formation rate is consistent with being on the main sequence of star-forming galaxies in the field (SFR~30-140 $M_{\odot}$/yr), and the cluster core molecular gas reservoir is expected to feed the AGN and the host star-formation before quiescence. The other cluster confirmed members show star formation rates at ~2 $\sigma$ below the field main sequence at similar redshifts and do not have molecular gas masses larger than galaxies of similar stellar mass in the field.


title: An algebraic groups perspective on Erd\H{o}s-Ko-Rado

abstract: We give a proof of the Erd\H{o}s-Ko-Rado Theorem using the Borel Fixed Point Theorem from algebraic group theory. This perspective gives a strong analogy between the Erd\H{o}s-Ko-Rado Theorem and (generalizations of) the Gerstenhaber Theorem on spaces of nilpotent matrices.


title: Soft displaced leptons at the LHC

abstract: Soft displaced leptons are representative collider signatures of compressed dark sectors with feeble couplings to the standard model. Prime targets are dark matter scenarios where co-scattering or co-annihilation sets the relic abundance upon freeze-out. At the LHC, searches for soft displaced leptons are challenged by a large background from hadron or tau lepton decays. In this article, we present an analysis tailored for displaced leptons with a low transverse momentum threshold at 20 GeV. Using a neural network, we perform a comprehensive analysis of the event kinematics, including a study of the expected detection efficiencies and backgrounds at small momenta. Our results show that weak-scale particles decaying into soft leptons with decay lengths between 1mm and 1m can be probed with LHC Run 2 data. This motivates the need for dedicated triggers that maximize the sensitivity to displaced soft leptons.


title: Central Charges for AdS Black Holes

abstract: Nontrivial diffeomorphisms act on the horizon of a generic 4D black holes and create distinguishing features referred to as soft hair. Amongst these are a left-right pair of Virasoro algebras with associated charges that reproduce the Bekenstein-Hawking entropy for Kerr black holes. In this paper we show that if one adds a negative cosmological constant, there is a similar set of infinitesimal diffeomorphisms that act non-trivially on the horizon. The algebra of these diffeomorphisms gives rise to a central charge. Adding a boundary counterterm, justified to achieve integrability, leads to well-defined central charges with cL = cR. The macroscopic area law for Kerr-AdS black holes follows from the assumption of a Cardy formula governing the black hole microstates.


title: Cryptanalysis of Quantum Secure Direct Communication Protocol with Mutual Authentication Based on Single Photons and Bell States

abstract: Recently, Yan et al. proposed a quantum secure direct communication (QSDC) protocol with authentication using single photons and Einstein-Podolsky-Rosen (EPR) pairs (Yan et al., CMC-Computers, Materials \& Continua, 63(3), 2020). In this work, we show that the QSDC protocol is not secure against intercept-and-resend attack and impersonation attack. An eavesdropper can get the full secret message by applying these attacks. We propose a modification of this protocol, which defeats the above attacks along with all the familiar attacks.


title: Separate Universe calibration of the dependence of halo bias on cosmic web anisotropy

abstract: We use the Separate Universe technique to calibrate the dependence of linear and quadratic halo bias $b_1$ and $b_2$ on the local cosmic web environment of dark matter haloes. We do this by measuring the response of halo abundances at fixed mass and cosmic web tidal anisotropy $\alpha$ to an infinite wavelength initial perturbation. We augment our measurements with an analytical framework developed in earlier work which exploits the near-Lognormal shape of the distribution of $\alpha$ and results in very high precision calibrations. We present convenient fitting functions for the dependence of $b_1$ and $b_2$ on $\alpha$ over a wide range of halo mass for redshifts $0\leq z\leq1$. Our calibration of $b_2(\alpha)$ is the first demonstration to date of the dependence of non-linear bias on the local web environment. Motivated by previous results which showed that $\alpha$ is the primary indicator of halo assembly bias for a number of halo properties beyond halo mass, we then extend our analytical framework to accommodate the dependence of $b_1$ and $b_2$ on any such secondary property which has, or can be monotonically transformed to have, a Gaussian distribution. We demonstrate this technique for the specific case of halo concentration, finding good agreement with previous results. Our calibrations will be useful for a variety of halo model analyses focusing on galaxy assembly bias, as well as analytical forecasts of the potential for using $\alpha$ as a segregating variable in multi-tracer analyses.


title: Two dialects for KZB equations: generating one-loop open-string integrals

abstract: Two different constructions generating the low-energy expansion of genus-one configuration-space integrals appearing in one-loop open-string amplitudes have been put forward in \rcites{Mafra:2019xms, *Mafra:2019ddf, Broedel:2019gba}. We are going to show that both approaches can be traced back to an elliptic system of Knizhnik--Zamolodchikov--Bernard(KZB) type on the twice-punctured torus. We derive an explicit all-multiplicity representation of the elliptic KZB system for a vector of iterated integrals with an extra marked point and explore compatibility conditions for the two sets of algebra generators appearing in the two differential equations.


title: Two-Loop Rational Terms in Yang-Mills Theories

abstract: Scattering amplitudes in $D$ dimensions involve particular terms that originate from the interplay of UV poles with the $D-4$ dimensional parts of loop numerators. Such contributions can be controlled through a finite set of process-independent rational counterterms, which make it possible to compute loop amplitudes with numerical tools that construct the loop numerators in four dimensions. Building on a recent study [1] of the general properties of two-loop rational counterterms, in this paper we investigate their dependence on the choice of renormalisation scheme. We identify a nontrivial form of scheme dependence, which originates from the interplay of mass and field renormalisation with the $D-4$ dimensional parts of loop numerators, and we show that it can be controlled through a new kind of one-loop counterterms. This guarantees that the two-loop rational counterterms for a given renormalisable theory can be derived once and for all in terms of generic renormalisation constants, which can be adapted a posteriori to any scheme. Using this approach, we present the first calculation of the full set of two-loop rational counterterms in Yang-Mills theories. The results are applicable to SU(N) and U(1) gauge theories coupled to $n_{f}$ fermions with arbitrary masses.


title: Towards an Understanding of Residual Networks Using Neural Tangent Hierarchy (NTH)

abstract: Gradient descent yields zero training loss in polynomial time for deep neural networks despite non-convex nature of the objective function. The behavior of network in the infinite width limit trained by gradient descent can be described by the Neural Tangent Kernel (NTK) introduced in \cite{Jacot2018Neural}. In this paper, we study dynamics of the NTK for finite width Deep Residual Network (ResNet) using the neural tangent hierarchy (NTH) proposed in \cite{Huang2019Dynamics}. For a ResNet with smooth and Lipschitz activation function, we reduce the requirement on the layer width $m$ with respect to the number of training samples $n$ from quartic to cubic. Our analysis suggests strongly that the particular skip-connection structure of ResNet is the main reason for its triumph over fully-connected network.


title: Discovery of the optical afterglow and host galaxy of short GRB181123B at $z =1.754$: Implications for Delay Time Distributions

abstract: We present the discovery of the optical afterglow and host galaxy of the {\it Swift} short-duration gamma-ray burst, GRB\,181123B. Observations with Gemini-North starting at $\approx 9.1$~hr after the burst reveal a faint optical afterglow with $i\approx25.1$~mag, at an angular offset of 0.59 $\pm$ 0.16$''$ from its host galaxy. Using $grizYJHK$ observations, we measure a photometric redshift of the host galaxy of $z = 1.77^{+0.30}_{-0.17}$. From a combination of Gemini and Keck spectroscopy of the host galaxy spanning 4500-18000~\AA , we detect a single emission line at 13390~\AA, inferred as H$\beta$ at $z = 1.754 \pm 0.001$ and corroborating the photometric redshift. The host galaxy properties of GRB\,181123B are typical to those of other SGRB hosts, with an inferred stellar mass of $\approx 1.7 \times 10^{10}\,M_{\odot}$, mass-weighted age of $\approx 0.9$~Gyr and optical luminosity of $\approx 0.9L^{*}$. At $z=1.754$, GRB\,181123B is the most distant secure SGRB with an optical afterglow detection, and one of only three at $z>1.5$. Motivated by a growing number of high-$z$ SGRBs, we explore the effects of a missing $z>1.5$ SGRB population among the current {\it Swift} sample on delay time distribution models. We find that log-normal models with mean delay times of $\approx 4-6$~Gyr are consistent with the observed distribution, but can be ruled out to $95\%$ confidence with an additional $\approx1-5$~{\it Swift} SGRBs recovered at $z>1.5$. In contrast, power-law models with $\propto$ $t^{-1}$ are consistent with the redshift distribution and can accommodate up to $\approx30$ SGRBs at these redshifts. Under this model, we predict that $\approx 1/3$ of the current {\it Swift} population of SGRBs is at $z>1$. The future discovery or recovery of existing high-$z$ SGRBs will provide significant discriminating power on their delay time distributions, and thus their formation channels.


title: Identification of CaCuSi$_4$O$_{10}$ (Egyptian blue) in the "Birch. Spring" painting by Robert Falk (1907) using photoluminescence

abstract: We have detected Egyptian blue pigment in the paint layer of the "Birch. Spring" painting by Robert Falk (1907); we have also found this pigment in the paints of the sketch drawn on the canvas back side. This is probably the first discovery of Egyptian blue in a 20th century work of art. We have analyzed a modern commercial Egyptian blue pigment (Kremer) and found it to be suitable as a standard for photoluminescence spectral analysis. The characteristic photoluminescence band of CaCuSi$_4$O$_{10}$ reaches maximum at the wavelength of about 910 nm. The luminescence is efficiently excited by incoherent green or blue light. The study demonstrates that the photoluminescence spectral micro analysis using excitation by incoherent light can be efficiently used for the identification of luminescent pigments in paint layers of artworks.


title: Heisenberg-Kitaev model in a magnetic field: $1/S$ expansion

abstract: The exact solution of Kitaev's spin-$1/2$ honeycomb spin-liquid model has sparked an intense search for Mott insulators hosting bond-dependent Kitaev interactions, of which $\mathrm{Na}_{2}\mathrm{IrO}_{3}$ and $\alpha-\mathrm{RuCl}_{3}$ are prime examples. Subsequently, it has been proposed that also spin-$1$ and spin-$3/2$ analogs of Kitaev interactions may occur in materials with strong spin-orbit coupling. As a minimal model to describe these Kitaev materials, we study the Heisenberg-Kitaev Hamiltonian in a consistent $1/S$ expansion, with $S$ being the spin size. We present a comprehensive study of this model in the presence of an external magnetic field applied along two different directions, [001] and [111], for which an intricate classical phase diagram has been reported. In both settings, we employ spin-wave theory in a number of ordered phases to compute phase boundaries at the next-to-leading order in $1/S$ and show that quantum corrections substantially modify the classical phase diagram. More broadly, our work presents a consistent route to investigate the leading quantum corrections in spin models that break spin-rotational symmetry.


title: Line Positions and Intensities of the {\nu}$_4$ Band of Methyl Iodide Using Mid-Infrared Optical Frequency Comb Fourier Transform Spectroscopy

abstract: We use optical frequency comb Fourier transform spectroscopy to measure high-resolution spectra of iodomethane, CH$_3$I in the C-H stretch region from 2800 to 3160 cm$^{-1}$. The fast-scanning Fourier transform spectrometer with auto-balanced detection is based on a difference frequency generation comb with repetition rate, f$_{rep}$, of 125 MHz. A series of spectra with sample point spacing equal to f$_{rep}$ are measured at different f$_{rep}$ settings and interleaved to yield sampling point spacing of 11 MHz. Iodomethane is introduced into a 76 m long multipass absorption cell by its vapor pressure at room temperature. The measured spectrum contains three main ro-vibrational features: the parallel vibrational overtone and combination bands centered around 2850 cm$^{-1}$, the symmetric stretch ${\nu}_1$ band centered at 2971 cm$^{-1}$, and the asymmetric stretch ${\nu}_4$ band centered at 3060 cm$^{-1}$. The spectra of the ${\nu}_4$ band and the nearby ${\nu}_3$+${\nu}_4$-${\nu}_3$ hot band are simulated using PGOPHER and a new assignment of these bands is presented. The resolved ro-vibrational structures are used in a least square fit together with the microwave data to provide the upper state parameters. We assign 2603 transitions to the ${\nu}_4$ band with standard deviation (observed - calculated) of 0.00034 cm$^{-1}$, and 831 transitions to the ${\nu}_3$+${\nu}_4$-${\nu}_3$ hot band with standard deviation of 0.00084 cm$^{-1}$. The hyperfine splittings due to the ${^{127}}$I nuclear quadrupole moment are observed for transitions with J$\leq$2xK. Finally, intensities of 157 isolated transitions in the ${\nu}_4$ band are reported for the first time using the Voigt line shape as a model in multispectral fitting.


title: ViscNet: Neural network for predicting the fragility index and the temperature-dependency of viscosity

abstract: Viscosity ($\eta$) is one of the most important properties of disordered matter. The temperature-dependence of viscosity is used to adjust process variables for glass-making, from melting to annealing. The aim of this work was to develop a physics-informed machine learning model capable of predicting $\eta(T)$ of oxide liquids. Instead of predicting the viscosity itself, the NN predicts the parameters of the MYEGA viscosity equation: the liquid's fragility index, the glass transition temperature, and the asymptotic viscosity. With these parameters, $\eta$ can be computed at any temperature of interest, with the advantage of good extrapolation capabilities inherent to the MYEGA equation. The dataset was collected from the SciGlass database; only oxide liquids with enough data points in the high and low viscosity regions were selected, resulting in a final dataset with 17,584 data points containing 847 different liquids. About 600 features were engineered from the liquids' chemical composition and 35 of these features were selected using a feature selection protocol. The hyperparameter (HP) tuning of the NN was performed in a set of experiments using both random search and Bayesian strategies, where a total of 700 HP sets were tested. The most successful HP sets were further tested using 10-fold cross-validation, and the one with the lowest average validation loss was selected as the best set. The final trained NN was tested with a test dataset of 85 liquids with different compositions than those used for training and validating the NN. The $R^2$ for the test dataset's prediction was 0.97. This work introduces three advantages: the model can predict viscosity as well as the liquids' glass transition temperature and fragility index; the model is designed and trained with a focus on extrapolation; finally, the model is available as free and open-source software licensed under the GPL3.


title: Machine Learning Classification of Kuiper Belt Populations

abstract: In the outer solar system, the Kuiper Belt contains dynamical sub-populations sculpted by a combination of planet formation and migration and gravitational perturbations from the present-day giant planet configuration. The subdivision of observed Kuiper Belt objects (KBOs) into different dynamical classes is based on their current orbital evolution in numerical integrations of their orbits. Here we demonstrate that machine learning algorithms are a promising tool for reducing both the computational time and human effort required for this classification. Using a Gradient Boosting Classifier, a type of machine learning regression tree classifier trained on features derived from short numerical simulations, we sort observed KBOs into four broad, dynamically distinct populations - classical, resonant, detached, and scattering - with a >97 per cent accuracy for the testing set of 542 securely classified KBOs. Over 80 per cent of these objects have a $>3\sigma$ probability of class membership, indicating that the machine learning method is classifying based on the fundamental dynamical features of each population. We also demonstrate how, by using computational savings over traditional methods, we can quickly derive a distribution of class membership by examining an ensemble of object clones drawn from the observational errors. We find two major reasons for misclassification: inherent ambiguity in the orbit of the object - for instance, an object that is on the edge of resonance - and a lack of representative examples in the training set. This work provides a promising avenue to explore for fast and accurate classification of the thousands of new KBOs expected to be found by surveys in the coming decade.


title: Property G and the $4$--genus

abstract: We say a null-homologous knot $K$ in a $3$--manifold $Y$ has Property G, if the properties about the Thurston norm and fiberedness of the complement of $K$ is preserved under the zero surgery on $K$. In this paper, we will show that, if the smooth $4$--genus of $K\times\{0\}$ (in a certain homology class) in $(Y\times[0,1])\#N\overline{\mathbb CP^2}$, where $Y$ is a rational homology sphere, is smaller than the Seifert genus of $K$, then $K$ has Property G. When the smooth $4$--genus is $0$, $Y$ can be taken to be any closed, oriented $3$--manifold.


title: Learning excursion sets of vector-valued Gaussian random fields for autonomous ocean sampling

abstract: Improving and optimizing oceanographic sampling is a crucial task for marine science and maritime resource management. Faced with limited resources in understanding processes in the water-column, the combination of statistics and autonomous systems provide new opportunities for experimental design. In this work we develop efficient spatial sampling methods for characterizing regions defined by simultaneous exceedances above prescribed thresholds of several responses, with an application focus on mapping coastal ocean phenomena based on temperature and salinity measurements. Specifically, we define a design criterion based on uncertainty in the excursions of vector-valued Gaussian random fields, and derive tractable expressions for the expected integrated Bernoulli variance reduction in such a framework. We demonstrate how this criterion can be used to prioritize sampling efforts at locations that are ambiguous, making exploration more effective. We use simulations to study and compare properties of the considered approaches, followed by results from field deployments with an autonomous underwater vehicle as part of a study mapping the boundary of a river plume. The results demonstrate the potential of combining statistical methods and robotic platforms to effectively inform and execute data-driven environmental sampling.


title: Collisional loss of one-dimensional fermions near a $p$-wave Feshbach resonance

abstract: We study collisional loss of a quasi-one-dimensional (1D) spin-polarized Fermi gas near a $p$-wave Feshbach resonance in ultracold $^6$Li atoms. We measure the location of the $p$-wave resonance in quasi-1D and observe a confinement-induced shift and broadening. We find that the three-body loss coefficient $L_3$ as a function of the quasi-1D confinement has little dependence on confinement strength. We also analyze the atom loss with a two-step cascade three-body loss model in which weakly bound dimers are formed prior to their loss arising from atom-dimer collisions. Our data are consistent with this model. We also find a possible suppression in the rate of dimer relaxation with strong quasi-1D confinement. We discuss the implications of these measurements for observing $p$-wave pairing in quasi-1D.


title: Learning while Respecting Privacy and Robustness to Distributional Uncertainties and Adversarial Data

abstract: Data used to train machine learning models can be adversarial--maliciously constructed by adversaries to fool the model. Challenge also arises by privacy, confidentiality, or due to legal constraints when data are geographically gathered and stored across multiple learners, some of which may hold even an "anonymized" or unreliable dataset. In this context, the distributionally robust optimization framework is considered for training a parametric model, both in centralized and federated learning settings. The objective is to endow the trained model with robustness against adversarially manipulated input data, or, distributional uncertainties, such as mismatches between training and testing data distributions, or among datasets stored at different workers. To this aim, the data distribution is assumed unknown, and lies within a Wasserstein ball centered around the empirical data distribution. This robust learning task entails an infinite-dimensional optimization problem, which is challenging. Leveraging a strong duality result, a surrogate is obtained, for which three stochastic primal-dual algorithms are developed: i) stochastic proximal gradient descent with an $\epsilon$-accurate oracle, which invokes an oracle to solve the convex sub-problems; ii) stochastic proximal gradient descent-ascent, which approximates the solution of the convex sub-problems via a single gradient ascent step; and, iii) a distributionally robust federated learning algorithm, which solves the sub-problems locally at different workers where data are stored. Compared to the empirical risk minimization and federated learning methods, the proposed algorithms offer robustness with little computation overhead. Numerical tests using image datasets showcase the merits of the proposed algorithms under several existing adversarial attacks and distributional uncertainties.


title: Levy flights in steep potential wells: Langevin modeling versus direct response to energy landscapes

abstract: We investigate the non-Langevin relative of the L\'{e}vy-driven Langevin random system, under an assumption that both systems share a common (asymptotic, stationary, steady-state) target pdf. The relaxation to equilibrium in the fractional Langevin-Fokker-Planck scenario results from an impact of confining conservative force fields on the random motion. A non-Langevin alternative has a built-in direct response of jump intensities to energy (potential) landscapes in which the process takes place. We revisit the problem of L\'{e}vy flights in superharmonic potential wells, with a focus on the extremally steep well regime, and address the issue of its (spectral) "closeness" to the L\'{e}vy jump-type process confined in a finite enclosure with impenetrable (in particular reflecting) boundaries. The pertinent random system "in a box/interval" is expected to have a fractional Laplacian with suitable boundary conditions as a legitimate motion generator. The problem is, that in contrast to amply studied Dirichlet boundary problems, a concept of reflecting boundary conditions and the path-wise implementation of the pertinent random process in the vicinity of (or sharply at) reflecting boundaries are not unequivocally settled for L\'{e}vy processes. This ambiguity extends to fractional motion generators, for which nonlocal analogs of Neumann conditions are not associated with path-wise reflection scenarios at the boundary, respecting the impenetrability assumption.

categories: cond-mat.stat-mech math-ph math.MP math.SP quant-ph


title: Emergent classical universes from initial quantum states in a tomographical description

abstract: Quantum and classical physical states are represented in a unified way when they are described by symplectic tomography. Therefore this representation allows us to study directly the necessary conditions for a classical universe to emerge from a quantum state. In the previous study on the de Sitter universe this was done by comparing the classical limit of the quantum tomograms with the tomograms resulting from the classical cosmological equations. In this paper we first review these results and extend them to all the de Sitter models. We show further that these tomograms can be obtained directly from transposing the Wheeler-De Witt equation to the tomographic variables. Subsequently, because the classic limits of the quantum tomograms are identified with their asymptotic expressions, we found the necessary conditions to extend the previous results by taking the tomograms of the WKB solutions of the Wheeler-DeWitt equation with a any potential. Furthermore in the previous works we found that the de Sitter models undergo the quantum-to-classical transition when the cosmological constant decays to its present value, we discuss at the end how far we can extend this result to more general models. In the conclusions, after discussing any improvements and developments of the results of this work, we sketch a phenomenological approach from which to extract information about the initial states of the universe.


title: TripMD: Driving patterns investigation via Motif Analysis

abstract: Processing driving data and investigating driving behavior has been receiving an increasing interest in the last decades, with applications ranging from car insurance pricing to policy making. A common strategy to analyze driving behavior is to study the maneuvers being performance by the driver. In this paper, we propose TripMD, a system that extracts the most relevant driving patterns from sensor recordings (such as acceleration) and provides a visualization that allows for an easy investigation. Additionally, we test our system using the UAH-DriveSet dataset, a publicly available naturalistic driving dataset. We show that (1) our system can extract a rich number of driving patterns from a single driver that are meaningful to understand driving behaviors and (2) our system can be used to identify the driving behavior of an unknown driver from a set of drivers whose behavior we know.


title: Learning to Optimize Power Distribution Grids using Sensitivity-Informed Deep Neural Networks

abstract: Deep learning for distribution grid optimization can be advocated as a promising solution for near-optimal yet timely inverter dispatch. The principle is to train a deep neural network (DNN) to predict the solutions of an optimal power flow (OPF), thus shifting the computational effort from real-time to offline. Nonetheless, before training this DNN, one has to solve a large number of OPFs to create a labeled dataset. Granted the latter step can still be prohibitive in time-critical applications, this work puts forth an original technique for improving the prediction accuracy of DNNs by taking into account the sensitivities of the OPF minimizers with respect to the OPF parameters. By expanding on multiparametric programming, it is shown that although inverter control problems may exhibit dual degeneracy, the required sensitivities do exist in general and can be computed readily using the output of any standard quadratic program (QP) solver. Numerical tests showcase that sensitivity-informed deep learning can enhance prediction accuracy in terms of mean square error (MSE) by 2-3 orders of magnitude at minimal computational overhead. Improvements are more significant in the small-data regime, where a DNN has to learn to optimize using a few examples. Beyond multiparametric QPs, the approach is currently being generalized to parametric (non)-convex optimization problems.


title: Robust and effective eSIF preconditioning for general SPD matrices

abstract: We propose an unconditionally robust and highly effective preconditioner for general symmetric positive definite (SPD) matrices based on structured incomplete factorization (SIF), called enhanced SIF (eSIF) preconditioner. The original SIF strategy proposed recently derives a structured preconditioner by applying block diagonal preprocessing to the matrix and then compressing appropriate scaled off-diagonal blocks. Here, we use an enhanced scaling-and-compression strategy to design the new eSIF preconditioner. Some subtle modifications are made, such as the use of two-sided block triangular preprocessing. A practical multilevel eSIF scheme is then designed. We give rigorous analysis for both the enhanced scaling-and-compression strategy and the multilevel eSIF preconditioner. The new eSIF framework has some significant advantages and overcomes some major limitations of the SIF strategy. (i) With the same tolerance for compressing the off-diagonal blocks, the eSIF preconditioner can approximate the original matrix to a much higher accuracy. (ii) The new preconditioner leads to much more significant reductions of condition numbers due to an accelerated magnification effect for the decay in the singular values of the scaled off-diagonal blocks. (iii) With the new preconditioner, the eigenvalues of the preconditioned matrix are much better clustered around $1$. (iv) The multilevel eSIF preconditioner is further unconditionally robust or is guaranteed to be positive definite without the need of extra stabilization, while the multilevel SIF preconditioner has a strict requirement in order to preserve positive definiteness. Comprehensive numerical tests are used to show the advantages of the eSIF preconditioner in accelerating the convergence of iterative solutions.


title: Detection as Regression: Certified Object Detection by Median Smoothing

abstract: Despite the vulnerability of object detectors to adversarial attacks, very few defenses are known to date. While adversarial training can improve the empirical robustness of image classifiers, a direct extension to object detection is very expensive. This work is motivated by recent progress on certified classification by randomized smoothing. We start by presenting a reduction from object detection to a regression problem. Then, to enable certified regression, where standard mean smoothing fails, we propose median smoothing, which is of independent interest. We obtain the first model-agnostic, training-free, and certified defense for object detection against $\ell_2$-bounded attacks. The code for all experiments in the paper is available at http://github.com/Ping-C/CertifiedObjectDetection .


title: PINT: Probabilistic In-band Network Telemetry

abstract: Commodity network devices support adding in-band telemetry measurements into data packets, enabling a wide range of applications, including network troubleshooting, congestion control, and path tracing. However, including such information on packets adds significant overhead that impacts both flow completion times and application-level performance. We introduce PINT, an in-band telemetry framework that bounds the amount of information added to each packet. PINT encodes the requested data on multiple packets, allowing per-packet overhead limits that can be as low as one bit. We analyze PINT and prove performance bounds, including cases when multiple queries are running simultaneously. PINT is implemented in P4 and can be deployed on network devices. Using real topologies and traffic characteristics, we show that PINT concurrently enables applications such as congestion control, path tracing, and computing tail latencies, using only sixteen bits per packet, with performance comparable to the state of the art.


title: Space-time smoothing models for sub-national measles routine immunization coverage estimation with complex survey data

abstract: Despite substantial advances in global measles vaccination, measles disease burden remains high in many low- and middle-income countries. A key public health strategy for controling measles in such high-burden settings is to conduct supplementary immunization activities (SIAs) in the form of mass vaccination campaigns, in addition to delivering scheduled vaccination through routine immunization (RI) programs. To achieve balanced implementations of RI and SIAs, robust measurement of sub-national RI-specific coverage is crucial. In this paper, we develop a space-time smoothing model for estimating RI-specific coverage of the first dose of measles-containing-vaccines (MCV1) at sub-national level using complex survey data. The application that motivated this work is estimation of the RI-specific MCV1 coverage in Nigeria's 36 states and the Federal Capital Territory. Data come from four Demographic and Health Surveys, three Multiple Indicator Cluster Surveys, and two National Nutrition and Health Surveys conducted in Nigeria between 2003 and 2018. Our method incorporates information from the SIA calendar published by the World Health Organization and accounts for the impact of SIAs on the overall MCV1 coverage, as measured by cross-sectional surveys. The model can be used to analyze data from multiple surveys with different data collection schemes and construct coverage estimates with uncertainty that reflects the various sampling designs. Implementation of our method can be done efficiently using integrated nested Laplace approximation (INLA).


title: Dipolar localization of waves in twisted phononic crystal plates

abstract: The localization of waves in two-dimensional clusters of scatterers arranged in relatively twisted lattices is studied by multiple scattering theory. It is found that, for a given frequency, it is always possible to find localized modes for a discrete set of rotation angles, analogous to the so-called "magic angles" recently found in two-dimensional materials like graphene. Similarly, for small rotations of the lattices, a large number of resonant frequencies is found, whose position strongly depends on the rotation angle. Moreover, for angles close to those that make the two lattices commensurable a single mode appears that can be easily tuned by the rotation angle. Unlike other twisted materials, where the properties of the bilayers are mainly explained in terms of the dispersion relation of the individual lattices, the special angles in these clusters happen because of the formation of dipolar scatterers due to the relative rotation between the two lattices, enhancing therefore their interaction. While the presented results are valid for any type of wave, the specific case of flexural waves in thin elastic plates is numerically studied, and the different modes found are comprehensively explained in terms of the interaction between pairs of scatterers. The analysis presented here shows that these structures are promising candidates for the inverse design of tunable wave-trapping devices for classical waves.


title: FEDReD II : 3D Extinction Map with 2MASS and Gaia DR2 data

abstract: Aims. We aim to map the 3D distribution of the interstellar extinction of the Milky Way disk up to distances larger than those probed with the Gaia parallax alone. Methods. We apply the FEDReD (Field Extinction-Distance Relation Deconvolver) algorithm to the 2MASS near-infrared photometry together with the Gaia DR2 astrometry and photometry. This algorithm uses a Bayesian deconvolution approach, based on an empirical HR-diagram representative of the local thin disk, in order to map the extinction as a function of distance of various fields of view. Results. We analysed more than 5.6 million stars to obtain an extinction map of the entire Galactic disk within $|b|<0.24^{\circ}$. This map provides information up to $5~\mathrm{kpc}$ in the direction of the Galactic centre and at more than $7~\mathrm{kpc}$ in the direction of the anticentre. This map reveals the complete shape of structures known locally, such as the Vela complex or the split of the local arm. Furthermore our extinction map shows many large "clean bubbles" especially one in the Sagittarius -- Carina complex, and four others which define a structure that we nickname the butterfly.


title: Chern Insulators and Topological Flat-bands in Magic-angle Twisted Bilayer Graphene

abstract: Magic-angle twisted bilayer graphene (MA-TBG) exhibits intriguing quantum phase transitions triggered by enhanced electron-electron interactions when its flat-bands are partially filled. However, the phases themselves and their connection to the putative non-trivial topology of the flat bands are largely unexplored. Here we report transport measurements revealing a succession of doping-induced Lifshitz transitions that are accompanied by van Hove singularities (VHS) which facilitate the emergence of correlation-induced gaps and topologically non-trivial sub-bands. In the presence of a magnetic field, well quantized Hall plateaus at filling of 1, 2, 3 carriers per moir\'e-cell reveal the sub-band topology and signal the emergence of Chern insulators with Chern-numbers, ! = !, !, !, respectively. Surprisingly, for magnetic fields exceeding 5T we observe a VHS at a filling of 3.5, suggesting the possibility of a fractional Chern insulator. This VHS is accompanied by a crossover from low-temperature metallic, to high-temperature insulating behavior, characteristic of entropically driven Pomeranchuk-like transitions,


title: On exponential bases and frames with non-linear phase functions and some applications

abstract: In this paper, we study the spectrality and frame-spectrality of exponential systems of the type $E(\Lambda,\varphi) = \{e^{2\pi i \lambda\cdot\varphi(x)}: \lambda\in\Lambda\}$ where the phase function $\varphi$ is a Borel measurable which is not necessarily linear. A complete characterization of pairs $(\Lambda,\varphi)$ for which $E(\Lambda,\varphi)$ is an orthogonal basis or a frame for $L^{2}(\mu)$ is obtained. In particular, we show that the middle-third Cantor measures and the unit disc, each admits an orthogonal basis with a certain non-linear phase. Under a natural regularity condition on the phase functions, when $\mu$ is the Lebesgue measure on $[0,1]$ and $\Lambda= {\mathbb{Z}},$ we show that only the standard phase functions $\varphi(x) = \pm x$ are the only possible functions that give rise to orthonormal bases. Surprisingly, however we prove that there exist a greater degree of flexibility, even for continuously differentiable phase functions in higher dimensions. For instance, we were able to describe a large class of functions $\varphi$ defined on ${\mathbb{R}}^{d}$ such that the system $E(\Lambda,\varphi)$ is an orthonormal basis for $L^{2}[0,1]^{d}$ when $d\geq2.$ Moreover, we discuss how our results apply to the discretization problem of unitary representations of locally compact groups for the construction of orthonormal bases. Finally, we conclude the paper by stating several open problems.


title: Cooperative half-guards in art galleries

abstract: In any simple polygonal art gallery with n walls, we show that it is possible to place floor(n/2)-1 guards whose range of vision is 180 degrees in such a way that every interior point of the gallery can be seen by one of them, and such that the mutual visibility graph formed by the guards is connected. This upper bound is tight, in that there exist galleries which require this number of guards, and equals the known result for guards with full 360 degree range of vision. We also show that for orthogonal art galleries, this result may be improved to floor(n/2)-2 guards with 180 degree range of vision.


title: Skin effect and winding number in disordered non-Hermitian systems

abstract: Unlike their Hermitian counterparts, non-Hermitian (NH) systems may display an exponential sensitivity to boundary conditions and an extensive number of edge-localized states in systems with open boundaries, a phenomena dubbed the "non-Hermitian skin effect." The NH skin effect is one of the primary challenges to defining a topological theory of NH Hamiltonians, as the sensitivity to boundary conditions invalidates the traditional bulk-boundary correspondence. The NH skin effect has recently been connected to the winding number, a topological invariant unique to NH systems. In this paper, we extend the definition of the winding number to disordered NH systems by generalizing established results on disordered Hermitian topological insulators. Our real-space winding number is self-averaging, continuous as a function of the parameters in the problem, and remains quantized even in the presence of strong disorder. We verify that our real-space formula still predicts the NH skin effect, allowing for the possibility of predicting and observing the NH skin effect in strongly disordered NH systems. As an application we apply our results to predict a NH Anderson skin effect where a skin effect is developed as disorder is added to a clean system, and to explain recent results in optical funnels.


title: Unconventional superconductivity in a strongly correlated band-insulator without doping

abstract: We present a novel route for attaining unconventional superconductivity (SC) in a strongly correlated system without doping. In a simple model of a correlated band insulator (BI) at half-filling we demonstrate, based on a generalization of the projected wavefunctions method, that SC emerges when e-e interactions and the bare band-gap are both much larger than the kinetic energy, provided the system has sufficient frustration against the magnetic order. As the interactions are tuned, SC appears sandwiched between the correlated BI followed by a paramagnetic metal on one side, and a ferrimagnetic metal, antiferromagnetic (AF) half-metal, and AF Mott insulator phases on the other side.


title: Infrared fixed points of higher-spin fermions in topological semimetals

abstract: We determine the fate of interacting fermions described by the Hamiltonian $H=\textbf{p}\cdot \textbf{J}$ in three-dimensional topological semimetals with linear band crossing, where $\textbf{p}$ is momentum and $\textbf{J}$ are the spin-$j$ matrices for half-integer pseudospin $j\geq 3/2$. While weak short-range interactions are irrelevant at the crossing point due to the vanishing density of states, weak long-range Coulomb interactions lead to a renormalization of the band structure. Using a self-consistent perturbative renormalization group approach, we show that band crossings of the type $\textbf{p}\cdot \textbf{J}$ are unstable for $j\leq 7/2$. Instead, through an intriguing interplay between cubic crystal symmetry, band topology, and interaction effects, the system is attracted to a variety of infrared fixed points. We also unravel several other properties of higher-spin fermions for general $j$, such as the relation between fermion self-energy and free energy, or the vanishing of the renormalized charge. An $\text{O}(3)$ symmetric fixed point composed of equal chirality Weyl fermions is stable for $j\leq 7/2$ and very likely so for all $j$. We then explore the rich fixed point structure for $j=5/2$ in detail. We find additional attractive fixed points with enhanced $\text{O}(3)$ symmetry that host both emergent Weyl or massless Dirac fermions, and identify a puzzling, infrared stable, anisotropic fixed point without enhanced symmetry in close analogy to the known case of $j=3/2$.


title: When we can trust computers (and when we can't)

abstract: With the relentless rise of computer power, there is a widespread expectation that computers can solve the most pressing problems of science, and even more besides. We explore the limits of computational modelling and conclude that, in the domains of science and engineering that are relatively simple and firmly grounded in theory, these methods are indeed powerful. Even so, the availability of code, data and documentation, along with a range of techniques for validation, verification and uncertainty quantification, are essential for building trust in computer generated findings. When it comes to complex systems in domains of science that are less firmly grounded in theory, notably biology and medicine, to say nothing of the social sciences and humanities, computers can create the illusion of objectivity, not least because the rise of big data and machine learning pose new challenges to reproducibility, while lacking true explanatory power. We also discuss important aspects of the natural world which cannot be solved by digital means. In the long-term, renewed emphasis on analogue methods will be necessary to temper the excessive faith currently placed in digital computation.


title: Meta-active Learning in Probabilistically-Safe Optimization

abstract: Learning to control a safety-critical system with latent dynamics (e.g. for deep brain stimulation) requires taking calculated risks to gain information as efficiently as possible. To address this problem, we present a probabilistically-safe, meta-active learning approach to efficiently learn system dynamics and optimal configurations. We cast this problem as meta-learning an acquisition function, which is represented by a Long-Short Term Memory Network (LSTM) encoding sampling history. This acquisition function is meta-learned offline to learn high quality sampling strategies. We employ a mixed-integer linear program as our policy with the final, linearized layers of our LSTM acquisition function directly encoded into the objective to trade off expected information gain (e.g., improvement in the accuracy of the model of system dynamics) with the likelihood of safe control. We set a new state-of-the-art in active learning for control of a high-dimensional system with altered dynamics (i.e., a damaged aircraft), achieving a 46% increase in information gain and a 20% speedup in computation time over baselines. Furthermore, we demonstrate our system's ability to learn the optimal parameter settings for deep brain stimulation in a rat's brain while avoiding unwanted side effects (i.e., triggering seizures), outperforming prior state-of-the-art approaches with a 58% increase in information gain. Additionally, our algorithm achieves a 97% likelihood of terminating in a safe state while losing only 15% of information gain.


title: Maxwell equations in a curved spacetime: Spin optics approximation

abstract: We study propagation of high-frequency electromagnetic waves in a curved spacetime. We demonstrate how a modification of the standard geometric optics allows one to include the helicity dependent corrections into the equations of motion of circularly polarized beams of radiation. As a result, polarized light rays are still null but not geodesic curves. To achieve these results we construct null frames associated with a set of (non-geodesic) null rays and use these frames for description of the high-frequency wave propagation. We call this approach spin optics approximation. It is completely covariant and it can be used in an arbitrary time-dependent gravitational field.


title: Predictive Analytics for Water Asset Management: Machine Learning and Survival Analysis

abstract: Understanding performance and prioritizing resources for the maintenance of the drinking-water pipe network throughout its life-cycle is a key part of water asset management. Renovation of this vital network is generally hindered by the difficulty or impossibility to gain physical access to the pipes. We study a statistical and machine learning framework for the prediction of water pipe failures. We employ classical and modern classifiers for a short-term prediction and survival analysis to provide a broader perspective and long-term forecast, usually needed for the economic analysis of the renovation. To enrich these models, we introduce new predictors based on water distribution domain knowledge and employ a modern oversampling technique to remedy the high imbalance coming from the few failures observed each year. For our case study, we use a dataset containing the failure records of all pipes within the water distribution network in Barcelona, Spain. The results shed light on the effect of important risk factors, such as pipe geometry, age, material, and soil cover, among others, and can help utility managers conduct more informed predictive maintenance tasks.


title: Electric Vehicle transition in the UK

abstract: This paper provides an overview of the electric vehicle transition in the UK. The spatial disparity in the uptake of BEVs across the different regions is analysed using historic BEV sales. A forecast for future growth in BEVs (ignoring the impact of Covid-19) is performed using an s-curve model. Currently, South East England and Greater London have the highest BEV sales as a percentage of new vehicle sales. The spatial distribution of EV chargers across the different regions is also analysed. The spatial analysis clearly shows the regional disparity in the uptake of EV. South East England has the highest number of public chargers excluding Greater London. However, if we consider the number of EVs in that region, it has the second-lowest ratio of approx. 1 charger per 10 BEV. The lowest ratio being 0.8 in the West Midlands.


title: Transfer Learning for Motor Imagery Based Brain-Computer Interfaces: A Complete Pipeline

abstract: Transfer learning (TL) has been widely used in motor imagery (MI) based brain-computer interfaces (BCIs) to reduce the calibration effort for a new subject, and demonstrated promising performance. While a closed-loop MI-based BCI system, after electroencephalogram (EEG) signal acquisition and temporal filtering, includes spatial filtering, feature engineering, and classification blocks before sending out the control signal to an external device, previous approaches only considered TL in one or two such components. This paper proposes that TL could be considered in all three components (spatial filtering, feature engineering, and classification) of MI-based BCIs. Furthermore, it is also very important to specifically add a data alignment component before spatial filtering to make the data from different subjects more consistent, and hence to facilitate subsequential TL. Offline calibration experiments on two MI datasets verified our proposal. Especially, integrating data alignment and sophisticated TL approaches can significantly improve the classification performance, and hence greatly reduces the calibration effort.


title: On Cokriging, Neural Networks, and Spatial Blind Source Separation for Multivariate Spatial Prediction

abstract: Multivariate measurements taken at irregularly sampled locations are a common form of data, for example in geochemical analysis of soil. In practical considerations predictions of these measurements at unobserved locations are of great interest. For standard multivariate spatial prediction methods it is mandatory to not only model spatial dependencies but also cross-dependencies which makes it a demanding task. Recently, a blind source separation approach for spatial data was suggested. When using this spatial blind source separation method prior the actual spatial prediction, modelling of spatial cross-dependencies is avoided, which in turn simplifies the spatial prediction task significantly. In this paper we investigate the use of spatial blind source separation as a pre-processing tool for spatial prediction and compare it with predictions from Cokriging and neural networks in an extensive simulation study as well as a geochemical dataset.


title: Resonator networks for factoring distributed representations of data structures

abstract: The ability to encode and manipulate data structures with distributed neural representations could qualitatively enhance the capabilities of traditional neural networks by supporting rule-based symbolic reasoning, a central property of cognition. Here we show how this may be accomplished within the framework of Vector Symbolic Architectures (VSA) (Plate, 1991; Gayler, 1998; Kanerva, 1996), whereby data structures are encoded by combining high-dimensional vectors with operations that together form an algebra on the space of distributed representations. In particular, we propose an efficient solution to a hard combinatorial search problem that arises when decoding elements of a VSA data structure: the factorization of products of multiple code vectors. Our proposed algorithm, called a resonator network, is a new type of recurrent neural network that interleaves VSA multiplication operations and pattern completion. We show in two examples -- parsing of a tree-like data structure and parsing of a visual scene -- how the factorization problem arises and how the resonator network can solve it. More broadly, resonator networks open the possibility to apply VSAs to myriad artificial intelligence problems in real-world domains. A companion paper (Kent et al., 2020) presents a rigorous analysis and evaluation of the performance of resonator networks, showing it out-performs alternative approaches.


title: Sharp Analysis of Smoothed Bellman Error Embedding

abstract: The \textit{Smoothed Bellman Error Embedding} algorithm~\citep{dai2018sbeed}, known as SBEED, was proposed as a provably convergent reinforcement learning algorithm with general nonlinear function approximation. It has been successfully implemented with neural networks and achieved strong empirical results. In this work, we study the theoretical behavior of SBEED in batch-mode reinforcement learning. We prove a near-optimal performance guarantee that depends on the representation power of the used function classes and a tight notion of the distribution shift. Our results improve upon prior guarantees for SBEED in ~\citet{dai2018sbeed} in terms of the dependence on the planning horizon and on the sample size. Our analysis builds on the recent work of ~\citet{Xie2020} which studies a related algorithm MSBO, that could be interpreted as a \textit{non-smooth} counterpart of SBEED.


title: Deep Reinforcement Learning and its Neuroscientific Implications

abstract: The emergence of powerful artificial intelligence is defining new research directions in neuroscience. To date, this research has focused largely on deep neural networks trained using supervised learning, in tasks such as image classification. However, there is another area of recent AI work which has so far received less attention from neuroscientists, but which may have profound neuroscientific implications: deep reinforcement learning. Deep RL offers a comprehensive framework for studying the interplay among learning, representation and decision-making, offering to the brain sciences a new set of research tools and a wide range of novel hypotheses. In the present review, we provide a high-level introduction to deep RL, discuss some of its initial applications to neuroscience, and survey its wider implications for research on brain and behavior, concluding with a list of opportunities for next-stage research.


title: Resource-Aware Protocols for Network Cost-Sharing Games

abstract: We study the extent to which decentralized cost-sharing protocols can achieve good price of anarchy (PoA) bounds in network cost-sharing games with $n$ agents. We focus on the model of resource-aware protocols, where the designer has prior access to the network structure and can also increase the total cost of an edge(overcharging), and we study classes of games with concave or convex cost functions. We first consider concave cost functions and our main result is a cost-sharing protocol for symmetric games on directed acyclic graphs that achieves a PoA of $2+\varepsilon$ for some arbitrary small positive $\varepsilon$, which improves to $1+\varepsilon$ for games with at least two players. We also achieve a PoA of 1 for series-parallel graphs and show that no protocol can achieve a PoA better than $\Omega(\sqrt{n})$ for multicast games. We then also consider convex cost functions and prove analogous results for series-parallel networks and multicast games, as well as a lower bound of $\Omega(n)$ for the PoA on directed acyclic graphs without the use of overcharging.