title: Constraints on the magnetic field within a stratified outer core

abstract: Mounting evidence from both seismology and experiments on core composition suggests the existence of a layer of stably stratified fluid at the top of Earth's outer core. In this work we examine the structure of the geomagnetic field within such a layer, building on the important but little known work of Malkus (1979). We assume (i) an idealised magnetostrophic spherical model of the geodynamo neglecting inertia, viscosity and the solid inner core, and (ii) a strongly stratified layer of constant depth immediately below the outer boundary within which there is no spherically radial flow. Due to the restricted dynamics, Malkus showed that the geomagnetic field must obey a certain condition which is a more restrictive version of the condition of Taylor (1963). The nonlinear nature of these constraints makes finding a magnetic field that obeys them, here termed a Malkus state, a challenging task. Nevertheless, such Malkus states when constrained further by geomagnetic observations have the potential to probe the interior of the core. By focusing on a particular class of magnetic fields for which the Malkus constraints are linear, we describe a constructive method that turns any purely-poloidal field into an exact Malkus state by adding a suitable toroidal field. We consider poloidal fields following a prescribed smooth profile within the core that match observation-derived models of the magnetic field in either epoch 2015 or the 10000-yr time averaged field. Multiple possible solutions for the toroidal field exist, hence we determine the Malkus state of miumum toroidal energy and we find that it has a strong azimuthal toroidal field, larger than the observed poloidal component at the core-mantle boundary. For the 2015 field for a layer of depth 300 km, we estimate a root mean squared azimuthal toroidal field of 3 mT with a pointwise maximum of 8 mT occurring at a depth of about 70 km.


title: Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae

abstract: We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.


title: Distinct reduction of Knight shift in superconducting state of Sr$_2$RuO$_4$ under uniaxial strain

abstract: Shortly after the discovery of superconductivity in Sr$_2$RuO$_4$, spin-triplet pairing was proposed and further corroborated by a constant Knight shift ($K$) across the transition temperature (T$_c$). However, a recent experiment observed a drop in $K$ at T$_c$ which becomes larger under uniaxial strain, ruling out several spin-triplet scenarios. Here we show that inter-orbital spin-triplet pairing features a d-vector that rotates when uniaxial strain is applied, leading to a larger drop in the spin polarization perpendicular to the strain direction, distinct from spin-singlet pairing. We propose that anisotropic spin polarization under strain will ultimately differentiate triplet vs.~singlet pairing.


title: Electrons in the Young Solar Wind: First Results from the Parker Solar Probe

abstract: The Solar Wind Electrons Alphas and Protons experiment on the Parker Solar Probe (PSP) mission measures the three-dimensional electron velocity distribution function. We derive the parameters of the core, halo, and strahl populations utilizing a combination of fitting to model distributions and numerical integration for $\sim 100,000$ electron distributions measured near the Sun on the first two PSP orbits, which reached heliocentric distances as small as $\sim 0.17$ AU. As expected, the electron core density and temperature increase with decreasing heliocentric distance, while the ratio of electron thermal pressure to magnetic pressure ($\beta_e$) decreases. These quantities have radial scaling consistent with previous observations farther from the Sun, with superposed variations associated with different solar wind streams. The density in the strahl also increases; however, the density of the halo plateaus and even decreases at perihelion, leading to a large strahl/halo ratio near the Sun. As at greater heliocentric distances, the core has a sunward drift relative to the proton frame, which balances the current carried by the strahl, satisfying the zero-current condition necessary to maintain quasi-neutrality. Many characteristics of the electron distributions near perihelion have trends with solar wind flow speed, $\beta_e$, and/or collisional age. Near the Sun, some trends not clearly seen at 1 AU become apparent, including anti-correlations between wind speed and both electron temperature and heat flux. These trends help us understand the mechanisms that shape the solar wind electron distributions at an early stage of their evolution.


title: Assessing the best edit in perturbation-based iterative refinement algorithms to compute the median string

abstract: Strings are a natural representation of biological data such as DNA, RNA and protein sequences. The problem of finding a string that summarizes a set of sequences has direct application in relative compression algorithms for genome and proteome analysis, where reference sequences need to be chosen. Median strings have been used as representatives of a set of strings in different domains. However, several formulations of those problems are NP-Complete. Alternatively, heuristic approaches that iteratively refine an initial coarse solution by applying edit operations have been proposed. Recently, we investigated the selection of the optimal edit operations to speed up convergence without spoiling the quality of the approximated median string. We propose a novel algorithm that outperforms state of the art heuristic approximations to the median string in terms of convergence speed by estimating the effect of a perturbation in the minimization of the expressions that define the median strings. We present corpus of comparative experiments to validate these results.


title: Instanton Corrections and Emergent Strings

abstract: We study limits of infinite distance in the moduli space of 4d $\mathcal{N} = 2$ string compactifications, in which instanton effects dominate. We first consider trajectories in the hypermultiplet moduli space of type IIB Calabi-Yau compactifications. We observe a correspondence between towers of D-brane instantons and D-brane 4d strings, such that the lighter the string the more relevant the instanton effects are. The dominant instantons modify the classical trajectory such that the lightest D-brane string becomes tensionless even faster, while the other strings are prevented to go below the fundamental string tension. This lightest string is dual to a fundamental type IIB string and realises the Emergent String Conjecture. We also consider the vector multiplet moduli space of type I string theory on $K3 \times T^2$ , where quantum corrections can also become significant. Naively, we only find trajectories that correspond to decompactification limits, in apparent contradiction with the picture obtained in some dual setup.


title: Strain-Engineered High Responsivity MoTe2 Photodetector for Silicon Photonic Integrated Circuits

abstract: In integrated photonics, specific wavelengths are preferred such as 1550 nm due to low-loss transmission and the availability of optical gain in this spectral region. For chip-based photodetectors, layered two-dimensional (2D) materials bear scientific and technologically-relevant properties leading to strong light-matter-interaction devices due to effects such as reduced coulomb screening or excitonic states. However, no efficient photodetector in the telecommunication C-band using 2D materials has been realized yet. Here, we demonstrate a MoTe2-based photodetector featuring strong photoresponse (responsivity = 0.5 A/W) operating at 1550nm on silicon photonic waveguide enabled by engineering the strain (4%) inside the photo-absorbing transition-metal-dichalcogenide film. We show that an induced tensile strain of ~4% reduces the bandgap of MoTe2 by about 0.2 eV by microscopically measuring the work-function across the device. Unlike Graphene-based photodetectors relying on a gapless band structure, this semiconductor-2D material detector shows a ~100X improved dark current enabling an efficient noise-equivalent power of just 90 pW/Hz^0.5. Such strain-engineered integrated photodetector provides new opportunities for integrated optoelectronic systems.


title: OE-CAM: A Hybrid Opto-Electronic Content Addressable Memory

abstract: A content addressable memory (CAM) is a type of memory that implements a parallel search engine at its core. A CAM takes as an input a value and outputs the address where this value is stored in case of a match. CAMs are used in a wide range of applications including networking, cashing, neuromorphic associative memories, multimedia, and data analytics. Here, we introduce a novel opto-electronic CAM (OE-CAM) utilizing the integrated silicon photonic platform. In our approach, we explore the performance of an experimental OE-CAM and verify the efficiency of the device at 25Gbit/s while maintaining the bit integrity under noise conditions. We show that OE-CAM enables a) two orders of magnitude faster search functionality resulting in b) a five orders of magnitude lower power-delay-product compared to CAMs implementations based on other emerging technologies. This remarkable performance potential is achieved by utilizing i) a high parallelism of wavelength-division-multiplexing in the optical domain, combined with ii) 10s of GHz-fast opto-electronic components, packaged in iii) integrated photonics for 10-100s ps-short communication delays. We further verify the upper optical input power limit of this OE-CAM to be given by parasitic nonlinearities inside the silicon waveguides, and the minimal detectable optical power at the back-end photoreceivers responsivity given channel noise. Such energy-efficient and short-delay OE-CAMs could become a key component of functional photonic-augmented ASICS, co-processors, or smart sensors.


title: Artificial Synapse with Mnemonic Functionality using GSST-based Photonic Integrated Memory

abstract: Machine-learning tasks performed by neural networks demonstrated useful capabilities for producing reliable, and repeatable intelligent decisions. Integrated photonics, leveraging both component miniaturization and the wave-nature of the signals, can potentially outperform electronics architectures when performing inference tasks. However, the missing photon-photon force challenges non-volatile photonic device-functionality required for efficient neural networks. Here we present a novel concept and optimization of multi-level discrete-state non-volatile photonic memory based on an ultra-compact (<4um) hybrid phase change material GSST-silicon Mach Zehnder modulator, with low insertion losses (3dB), to serve as node in a photonic neural network. An optimized electro-thermal switching mechanism, induced by Joule heating through tungsten contacts, is engineered. This operation allows to change the phase of the GSST film thus providing weight updating functionality to the network. We show that a 5 V pulse-train (<1 us, 20 pulses) applied to a serpentine contact produces crystallization and a single pulse of longer duration (2 us) amorphization, used to set the analog synaptic weights of a neuron. Emulating an opportunely trained 100x100 fully connected multilayered perceptron neural network with this weighting functionality embedded as photonic memory, shows up to 93% inference accuracy and robustness towards noise when performing predictions of unseen data


title: Reinforcement learning for bandwidth estimation and congestion control in real-time communications

abstract: Bandwidth estimation and congestion control for real-time communications (i.e., audio and video conferencing) remains a difficult problem, despite many years of research. Achieving high quality of experience (QoE) for end users requires continual updates due to changing network architectures and technologies. In this paper, we apply reinforcement learning for the first time to the problem of real-time communications (RTC), where we seek to optimize user-perceived quality. We present initial proof-of-concept results, where we learn an agent to control sending rate in an RTC system, evaluating using both network simulation and real Internet video calls. We discuss the challenges we observed, particularly in designing realistic reward functions that reflect QoE, and in bridging the gap between the training environment and real-world networks.


title: Channel Estimation, Interference Cancellation, and Symbol Detection for Communications on Overlapping Channels

abstract: In this paper, we propose the joint interference cancellation, fast fading channel estimation, and data symbol detection for a general interference setting where the interfering source and the interfered receiver are unsynchronized and occupy overlapping channels of different bandwidths. The interference must be canceled before the channel estimation and data symbol detection of the desired communication are performed. To this end, we have to estimate the Effective Interference Coefficients (EICs) and then the desired fast fading channel coefficients. We construct a two-phase framework where the EICs and desired channel coefficients are estimated using the joint maximum likelihood-maximum a posteriori probability (JML-MAP) criteria in the first phase; and the MAP based data symbol detection is performed in the second phase. Based on this two-phase framework, we also propose an iterative algorithm for interference cancellation, channel estimation and data detection. We analyze the channel estimation error, residual interference, symbol error rate (SER) achieved by the proposed framework. We then discuss how to optimize the pilot density to achieve the maximum throughput. Via numerical studies, we show that our design can effectively mitigate the interference for a wide range of SNR values, our proposed channel estimation and symbol detection design can achieve better performances compared to the existing method. Moreover, we demonstrate the improved performance of the iterative algorithm with respect to the non-iterative counterpart.


title: Magnetic Soft Charges, Dual Supertranslations and 't Hooft Line Dressings

abstract: We construct the Faddeev-Kulish asymptotic states in a quantum field theory of electric and magnetic charges. We find that there are two kind of dressings: apart from the well known (electric) Wilson line dressing, there is a magnetic counterpart which can be written as a 't Hooft line operator. The 't Hooft line dressings are charged under the magnetic large gauge transformation (LGT), but are neutral under electric LGT. This is in contrast to the Faddeev-Kulish dressings of electrons, which can be written as a Wilson line operator and are charged under electric LGT but neutral under magnetic LGT. With these dressings and the corresponding construction of the coherent states, the infrared finiteness of the theory of electric and magnetic charges is guaranteed. Even in the absence of magnetic monopoles, the electric and magnetic soft modes exhibit the electromagnetic duality of vacuum Maxwell theory. Using only the asymptotic form of three-point interactions in a field theory of electric and magnetic charges, we show that the leading magnetic dressings, like the leading electric ones, are exact in the field theory of electric and magnetic charges, in accordance with a conjecture of Strominger. We then extend the construction to perturbative quantum gravity in asymptotically flat spacetime, and construct gravitational 't Hooft line dressings that are charged under dual supertranslations. The duality in the quantum theory between the electric and magnetic soft charges and their dressings is thus made manifest.


title: Intrinsic Topological Transforms via the Distance Kernel Embedding

abstract: Topological transforms are parametrized families of topological invariants, which, by analogy with transforms in signal processing, are much more discriminative than single measurements. The first two topological transforms to be defined were the Persistent Homology Transform and Euler Characteristic Transform, both of which apply to shapes embedded in Euclidean space. The contribution of this paper is to define topological transforms that depend only on the intrinsic geometry of a shape, and hence are invariant to the choice of embedding. To that end, given an abstract metric measure space, we define an integral operator whose eigenfunctions are used to compute sublevel set persistent homology. We demonstrate that this operator, which we call the distance kernel operator, enjoys desirable stability properties, and that its spectrum and eigenfunctions concisely encode the large-scale geometry of our metric measure space. We then define a number of topological transforms using the eigenfunctions of this operator, and observe that these transforms inherit many of the stability and injectivity properties of the distance kernel operator.


title: An Estimate of the Surface Pollution of the Arctic Sea Ice

abstract: The Arctic sea ice represents an important energy reservoir for the climate of the northern hemisphere. The shrinking of the polar ice in the past decades decreases the stored energy and raises serious concerns about future climate changes.[1-4] Model calculations of the present authors [5,6] suggest that half of the global warming during the past fifty years is directly related to the retreat of the sea ice, while the cause is not well understood, e.g. the role of surface pollution [7-10]. We have analysed the reported annual melting and freezing data of the northern sea ice in the years 1979 to 2018 [11] to gain some insight. Two features can be deduced from our simple model: (i) recent results [12,13] are confirmed that approximately 60 % of the loss of sea ice stems from energy transport to the arctic region. (ii) We find evidence that the remaining part of the ice retreat originates from an increasing surface absorption of solar radiation, obviously due to the rising surface pollution of the sea ice. While the phenomenon was previously considered by several authors in a qualitative way, our analysis contributes semi-quantitative information on the situation. We estimate that the relevant fall-out of light absorbing aerosols onto the sea ice increased by 17 +/- 5 % during the past fifty years. A deposition of additional 3 +/- 1 % of solar radiation in the melting region results that accounts for the ice retreat. Recalling the important role of the ice loss for the terrestrial climate,[3,5,9] the precipitation of air pollution in the Arctic seems to be an important factor for the global warming.


title: Ab initio models of atomic nuclei: challenges and new ideas

abstract: This review presents some of the challenges in constructing models of atomic nuclei starting from theoretical descriptions of the strong interaction between nucleons. The focus is on statistical computing and methods for analyzing the link between bulk properties of atomic nuclei, such as radii and binding energies, and the underlying microscopic description of the nuclear interaction. The importance of careful model calibration and uncertainty quantification of theoretical predictions is highlighted.


title: Simulated Direct Imaging Detection of Water Vapor For Exo-Earths

abstract: Habitable planets are often defined as terrestrial worlds capable of maintaining surface liquid water. As a result, atmospheric water vapor can be a critical indicator of habitability. Thus, habitability-themed exoplanet investigations emphasize detection of water vapor signatures for their targets. A variety of concept missions for exoplanet direct imaging in reflected light have seen recent study, including the HabEx and LUVOIR concepts. Here, it is important to understand how direct imaging in reflected light --- coupled with moderate-resolution spectroscopy --- could be used to detect various water vapor amounts in Earth-like exoplanetary atmospheres. To investigate water vapor detection for terrestrial exoplanets, we generated reflectance spectra over a grid of water vapor column masses and used an instrument model to explore requisite integration times for spectral feature detection at either visible or near-infrared wavelengths. Lower-resolution near-infrared spectroscopy is generally optimal for detecting water vapor in the atmospheres of Earth-like exoplanets when using direct imaging in reflected light. This holds true for dry or cold terrestrial planets, whose atmospheres would contain relatively little water vapor. Atmospheres richer in water vapor, such as planets undergoing a moist or runaway greenhouse, could have water vapor efficiently detected at visible wavelengths. Understanding details like an exoplanet's size, temperature, and location relative to the habitable zone can aid in determining appropriate wavelength ranges for atmospheric characterization. Overall, water vapor detection for Earth-like exoplanets is quite feasible for future direct imaging missions.


title: The East Asian Observatory SCUBA--2 survey of the COSMOS field: unveiling 1147 bright sub-millimeter sources across 2.6 square degrees

abstract: We present sensitive 850$\mu$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $\sigma_{850\mu{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $\sigma_{850\mu{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$\sigma$ and >4.3$\sigma$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$\mu$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.


title: Detection of large-scale noisy multi-periodic patterns with discrete double Fourier transform

abstract: In many processes, the variations in underlying characteristics can be approximated by noisy multi-periodic patterns. If large-scale patterns are superimposed by a noise with long-range correlations, the detection of multi-periodic patterns becomes especially challenging. To solve this problem, we developed a discrete double Fourier transform (DDFT). DDFT is based on the equidistance property of harmonics generated by multi-periodic patterns in the discrete Fourier transform (DFT) spectra. As the large-scale patterns generate long enough equidistant series, they can be detected by the iteration of the primary DFT. DDFT is defined as Fourier transform of intensity spectral harmonics or of their functions. It comprises widely used cepstrum transform as a particular case. We present also the relevant analytical criteria for the assessment of statistical significance of peak harmonics in DDFT spectra in the presence of noise. DDFT technique was tested by extensive numerical simulations. The practical applications of DDFT technique are illustrated by the analysis of variations in solar wind speed related to solar rotation and by the study of large-scale multi-periodic patterns in DNA sequences. The latter application can be considered as generic example for the general spectral analysis of symbolic sequences. The results are compared with those obtained by the cepstrum transform. The mutual combination of DFT and DDFT provides an efficient technique to search for noisy large-scale multi-periodic patterns.


title: Estimating Large Mixed-Frequency Bayesian VAR Models

abstract: We discuss the issue of estimating large-scale vector autoregressive (VAR) models with stochastic volatility in real-time situations where data are sampled at different frequencies. In the case of a large VAR with stochastic volatility, the mixed-frequency data warrant an additional step in the already computationally challenging Markov Chain Monte Carlo algorithm used to sample from the posterior distribution of the parameters. We suggest the use of a factor stochastic volatility model to capture a time-varying error covariance structure. Because the factor stochastic volatility model renders the equations of the VAR conditionally independent, settling for this particular stochastic volatility model comes with major computational benefits. First, we are able to improve upon the mixed-frequency simulation smoothing step by leveraging a univariate and adaptive filtering algorithm. Second, the regression parameters can be sampled equation-by-equation in parallel. These computational features of the model alleviate the computational burden and make it possible to move the mixed-frequency VAR to the high-dimensional regime. We illustrate the model by an application to US data using our mixed-frequency VAR with 20, 34 and 119 variables.


title: Collective dynamics of pedestrians in a non-panic evacuation scenario

abstract: We present a study of pedestrian motion along a corridor in a non-panic regime (e.g., schools, hospitals or airports). Such situations have been discussed so far within the Social Force Model (SFM). We suggest to enrich this model by interactions based on the velocity of the particles and some randomness, both of which we introduce using the ideas of the Vicsek Model (VM). This new model allows to introduce fluctuations for a given average speed and geometry, and considering that the alignment interactions are modulated by an external control parameter (the noise $\eta$) allows to introduce phase transitions between ordered and disordered states. We have compared simulations of pedestrian motion along a corridor using (a) the VM with two boundary conditions (periodic and bouncing back) and with or without desired direction of motion, (b) the SFM, and (c) the new model SFM+VM. The study of steady-state configurations in the VM with confined geometry shows the expected bands perpendicular to the motion direction, while in the SFM and SFM+VM particles order in stripes of a given width $w$ along the direction of motion. The results in the SFM+VM case show that $w(t)\simeq t^\alpha$ has a diffusive-like behavior at low noise $\eta$ (dynamic exponent $\alpha \approx 1/2$), while it is sub-diffusive at high values of external noise ($\alpha \approx 1/4$). We observe the order-disorder transition in the VM with both boundary conditions, but the application of a desired direction condition inhibits the existence of disorder as expected. For the SFM+VM case we find a susceptibility maximum which increases with system size as a function of noise strength indicative of a order-disorder transition in the whole range of densities and speeds studied. From our results we conclude that the new SFM+VM model is a well-suited model to describe non-panic evacuation with diverse degrees of disorder.


title: Large-Scale Semi-Supervised Learning via Graph Structure Learning over High-Dense Points

abstract: We focus on developing a novel scalable graph-based semi-supervised learning (SSL) method for a small number of labeled data and a large amount of unlabeled data. Due to the lack of labeled data and the availability of large-scale unlabeled data, existing SSL methods usually encounter either suboptimal performance because of an improper graph or the high computational complexity of the large-scale optimization problem. In this paper, we propose to address both challenging problems by constructing a proper graph for graph-based SSL methods. Different from existing approaches, we simultaneously learn a small set of vertexes to characterize the high-dense regions of the input data and a graph to depict the relationships among these vertexes. A novel approach is then proposed to construct the graph of the input data from the learned graph of a small number of vertexes with some preferred properties. Without explicitly calculating the constructed graph of inputs, two transductive graph-based SSL approaches are presented with the computational complexity in linear with the number of input data. Extensive experiments on synthetic data and real datasets of varied sizes demonstrate that the proposed method is not only scalable for large-scale data, but also achieve good classification performance, especially for extremely small number of labels.


title: Matter Mixing in Aspherical Core-collapse Supernovae: Three-dimensional Simulations with Single Star and Binary Merger Progenitor Models for SN 1987A

abstract: We perform three-dimensional hydrodynamic simulations of aspherical core-collapse supernovae focusing on the matter mixing in SN 1987A. The impacts of four progenitor (pre-supernova) models and parameterized aspherical explosions are investigated. The four pre-supernova models include a blue supergiant (BSG) model based on a slow merger scenario developed recently for the progenitor of SN 1987A (Urushibata et al. 2018). The others are a BSG model based on a single star evolution and two red supergiant (RSG) models. Among the investigated explosion (simulation) models, a model with the binary merger progenitor model and with an asymmetric bipolar-like explosion, which invokes a jetlike explosion, best reproduces constraints on the mass of high velocity $^{56}$Ni, as inferred from the observed [Fe II] line profiles. The advantage of the binary merger progenitor model for the matter mixing is the flat and less extended $\rho \,r^3$ profile of the C+O core and the helium layer, which may be characterized by the small helium core mass. From the best explosion model, the direction of the bipolar explosion axis (the strongest explosion direction), the neutron star (NS) kick velocity, and its direction are predicted. Other related implications and future prospects are also given.


title: Streaming Classification of Variable Stars

abstract: In the last years, automatic classification of variable stars has received substantial attention. Using machine learning techniques for this task has proven to be quite useful. Typically, machine learning classifiers used for this task require to have a fixed training set, and the training process is performed offline. Upcoming surveys such as the Large Synoptic Survey Telescope (LSST) will generate new observations daily, where an automatic classification system able to create alerts online will be mandatory. A system with those characteristics must be able to update itself incrementally. Unfortunately, after training, most machine learning classifiers do not support the inclusion of new observations in light curves, they need to re-train from scratch. Naively re-training from scratch is not an option in streaming settings, mainly because of the expensive pre-processing routines required to obtain a vector representation of light curves (features) each time we include new observations. In this work, we propose a streaming probabilistic classification model; it uses a set of newly designed features that work incrementally. With this model, we can have a machine learning classifier that updates itself in real time with new observations. To test our approach, we simulate a streaming scenario with light curves from CoRot, OGLE and MACHO catalogs. Results show that our model achieves high classification performance, staying an order of magnitude faster than traditional classification approaches.

categories: astro-ph.IM astro-ph.GA astro-ph.SR cs.LG


title: Dynamical Coulomb Blockade as a Local Probe for Quantum Transport

abstract: Quantum fluctuations are imprinted with valuable information about transport processes. Experimental access to this information is possible, but challenging. We introduce the dynamical Coulomb blockade (DCB) as a local probe for fluctuations in a scanning tunneling microscope (STM) and show that it provides information about the conduction channels. In agreement with theoretical predictions, we find that the DCB disappears in a single-channel junction with increasing transmission following the Fano factor, analogous to what happens with shot noise. Furthermore we demonstrate local differences in the DCB expected from changes in the conduction channel configuration. Our experimental results are complemented by ab initio transport calculations that elucidate the microscopic nature of the conduction channels in our atomic-scale contacts. We conclude that probing the DCB by STM provides a technique complementary to shot noise measurements for locally resolving quantum transport characteristics.


title: Max-compound Cox processes. III

abstract: Extreme values are considered in samples with random size that has a mixed Poisson distribution being generated by a doubly stochastic Poisson process. We prove some inequalities providing bounds on the rate of convergence in limit theorems for the distributions of max-compound Cox processes.


title: The role of environment in galaxy evolution in the SERVS Survey I: density maps and cluster candidates

abstract: We use photometric redshifts derived from new $u$-band through 4.5$\mu$m Spitzer IRAC photometry in the 4.8\,deg$^2$ of the XMM-LSS field to construct surface density maps in the redshift range 0.1-1.5. Our density maps show evidence for large-scale structure in the form of filaments spanning several tens of Mpc. Using these maps, we identify 339 overdensities that our simulated lightcone analysis suggests are likely associated with dark matter haloes with masses, $M_{\rm halo}$, log($M_{\rm halo}/M_{\odot})>$13.7. From this list of overdensities we recover 43 of 70 known X-ray detected and spectroscopically confirmed clusters. The missing X-ray clusters are largely at lower redshifts and lower masses than our target log($M_{\rm halo}/M_{\odot})>$13.7. The bulk of the overdensities are compact, but a quarter show extended morphologies which include likely projection effects, clusters embedded in apparent filaments as well as at least one potential cluster merger (at $z\sim1.28$). The strongest overdensity in our highest redshift slice (at $z\sim1.5$) shows a compact red galaxy core potentially implying a massive evolved cluster.


title: Radon numbers grow linearly

abstract: Define the $k$-th Radon number $r_k$ of a convexity space as the smallest number (if it exists) for which any set of $r_k$ points can be partitioned into $k$ parts whose convex hulls intersect. Combining the recent abstract fractional Helly theorem of Holmsen and Lee with earlier methods of Bukh, we prove that $r_k$ grows linearly, i.e., $r_k\le c(r_2)\cdot k$.


title: The Ryu-Takayanagi Formula from Quantum Error Correction: An Algebraic Treatment of the Boundary CFT

abstract: It was recently shown by Harlow that any quantum error correcting code, satisfying the same complementary recovery properties as AdS/CFT, will obey a version of the Ryu-Takayanagi formula. In his most general result, Harlow allowed the bulk algebras to have nontrivial center, which was necessary for the "area operator" in this Ryu-Takayanagi formula to be nontrivial. However, the boundary Hilbert space was still assumed to factorise into Hilbert spaces associated with complementary boundary regions. We extend this work to include more general boundary theories, such as gauge theories, where the subalgebras associated with boundary regions may also have nontrivial center. We show the equivalence of a set of four conditions for a bulk algebra to be reconstructable from a boundary algebra, and then show that complementary recovery implies that the algebraic boundary entropy obeys a Ryu-Takayanagi formula. In contrast, we show that the distillable boundary entropy does not obey any such formula. If an additional "log dim R" term is added to the algebraic entropy, it will still obey a Ryu-Takayanagi formula, with a different area operator. However, since the "log dim R" term is a sum over local boundary contributions, we argue that it can only be related to the regularisation of the area at the bulk cut-off.


title: Learning from Interventions using Hierarchical Policies for Safe Learning

abstract: Learning from Demonstrations (LfD) via Behavior Cloning (BC) works well on multiple complex tasks. However, a limitation of the typical LfD approach is that it requires expert demonstrations for all scenarios, including those in which the algorithm is already well-trained. The recently proposed Learning from Interventions (LfI) overcomes this limitation by using an expert overseer. The expert overseer only intervenes when it suspects that an unsafe action is about to be taken. Although LfI significantly improves over LfD, the state-of-the-art LfI fails to account for delay caused by the expert's reaction time and only learns short-term behavior. We address these limitations by 1) interpolating the expert's interventions back in time, and 2) by splitting the policy into two hierarchical levels, one that generates sub-goals for the future and another that generates actions to reach those desired sub-goals. This sub-goal prediction forces the algorithm to learn long-term behavior while also being robust to the expert's reaction time. Our experiments show that LfI using sub-goals in a hierarchical policy framework trains faster and achieves better asymptotic performance than typical LfD.


title: A bi-integrated model for coupling lot-sizing and cutting-stock problems

abstract: In this paper, a framework that addresses the core of the papermaking process is proposed, starting from the production of jumbos and ending with the paper sheets used in daily life. The first phase of the process is modelled according to a lot-sizing problem, where the quantities of jumbos are determined in order to meet the demand of the entire chain. The second phase follows a one-dimensional cutting-stock formulation, where these jumbos are cut into smaller reels of predetermined lengths. Some of these are intended to fulfil a portfolio of orders, while others are used as raw material for the third phase of the process, when the reels are cut into sheets with specific dimensions and demands, following a two-dimensional cutting-stock problem. The model is called the Bi-Integrated Model, since it is composed of two integrated models. The heuristic method developed uses the Simplex Method with column generation for the two cutting-stock phases and applies the Relax-and-Fix technique to obtain the rounded-integer solution. Computational experiments comparing the solutions of the Bi-Integrated Model to other strategies of modelling the production process indicate average cost gains reaching 26.63%. Additional analyses of the model behaviour under several situations resulted in remarkable findings.


title: Electric field fluctuations in the two-dimensional Coulomb fluid

abstract: The structure factor for electric field correlations in the two dimensional Coulomb fluid is simulated and compared to theories of the dielectric function. Singular changes in the structure factor occur at the BKT insulator to conductor transition, as well as at a higher temperature correlation transition between a poor electrolyte and perturbed Debye-H\"uckel fluid. Structure factors are found to differ in the canonical and grand canonical ensembles, with the poor electrolyte showing full ensemble inequivalence. We identify mechanisms of `underscreening' and `pinch point' scattering that are relevant to experiments on ionic liquids and artificial spin ice respectively.


title: Magnetic connectivity of the ecliptic plane within 0.5 AU : PFSS modeling of the first PSP encounter

abstract: We compare magnetic field measurements taken by the FIELDS instrument on Parker Solar Probe (PSP) during its first solar encounter to predictions obtained by Potential Field Source Surface (PFSS) modeling. Ballistic propagation is used to connect the spacecraft to the source surface. Despite the simplicity of the model, our results show striking agreement with PSPs first observations of the heliospheric magnetic field from 0.5 AU (107.5 Rs) down to 0.16 AU (35.7 Rs). Further, we show the robustness of the agreement is improved both by allowing the photospheric input to the model to vary in time, and by advecting the field from PSP down to the PFSS model domain using in situ PSP/SWEAP measurements of the solar wind speed instead of assuming it to be constant with longitude and latitude. We also explore the source surface height parameter (RSS) to the PFSS model finding that an extraordinarily low source surface height (1.3-1.5Rs) predicts observed small scale polarity inversions which are otherwise washed out with regular modeling parameters. Finally, we extract field line traces from these models. By overlaying these on EUV images we observe magnetic connectivity to various equatorial and mid-latitude coronal holes indicating plausible magnetic footpoints and offering context for future discussions of sources of the solar wind measured by PSP.


title: Recursive Network Estimation From Binary-Valued Observation Data

abstract: This paper studies the problem of recursively estimating the weighted adjacency matrix of a network out of a temporal sequence of binary-valued observations. The observation sequence is generated from nonlinear networked dynamics in which agents exchange and display binary outputs. Sufficient conditions are given to ensure stability of the observation sequence and identifiability of the system parameters. It is shown that stability and identifiability can be guaranteed under the assumption of independent standard Gaussian disturbances. Via a maximum likelihood approach, the estimation problem is transformed into an optimization problem, and it is verified that its solution is the true parameter vector under the independent standard Gaussian assumption. A recursive algorithm for the estimation problem is then proposed based on stochastic approximation techniques. Its strong consistency is established and convergence rate analyzed. Finally, numerical simulations are conducted to illustrate the results and to show that the proposed algorithm is insensitive to small unmodeled factors.


title: Microlensing and Photon Bunching: The impact of decoherence

abstract: Gravitational microlensing within the Galaxy offers the prospect of probing the details of distant stellar sources, as well as revealing the distribution of compact (and potentially non-luminous) masses along the line-of-sight. Recently, it has been suggested that additional constraints on the lensing properties can be determined through the measurement of the time delay between images through the correlation of the bunching of photon arrival times; an application of the Hanbury-Brown Twiss effect. In this paper, we revisit this analysis, examining the impact of decoherence of the radiation from a spatially extended source along the multiple paths to an observer. The result is that, for physically reasonable situations, such decoherence completely erases any correlation that could otherwise be used to measure the gravitational lensing time delay. Indeed, the divergent light paths traverse extremely long effective baselines at the lens plane, corresponding to extremes of angular resolving power well beyond those attainable with any terrestrial technologies; the drawback being that few conceivable celestial objects would be sufficiently compact with high enough surface brightness to yield usable signals.


title: Semiregular and strongly irregular boundary points for $p$-harmonic functions on unbounded sets in metric spaces

abstract: The trichotomy between regular, semiregular, and strongly irregular boundary points for $p$-harmonic functions is obtained for unbounded open sets in complete metric spaces with a doubling measure supporting a $p$-Poincar\'e inequality, $1<p<\infty$. We show that these are local properties. We also deduce several characterizations of semiregular points and strongly irregular points. In particular, semiregular points are characterized by means of capacity, $p$-harmonic measures, removability, and semibarriers.


title: Physics-Informed Machine Learning with Conditional Karhunen-Lo\`eve Expansions

abstract: We present a new physics-informed machine learning approach for the inversion of PDE models with heterogeneous parameters. In our approach, the space-dependent partially-observed parameters and states are approximated via Karhunen-Lo\`eve expansions (KLEs). Each of these KLEs is then conditioned on their corresponding measurements, resulting in low-dimensional models of the parameters and states that resolve observed data. Finally, the coefficients of the KLEs are estimated by minimizing the norm of the residual of the PDE model evaluated at a finite set of points in the computational domain, ensuring that the reconstructed parameters and states are consistent with both the observations and the PDE model to an arbitrary level of accuracy. In our approach, KLEs are constructed using the eigendecomposition of covariance models of spatial variability. For the model parameters, we employ a parameterized covariance model calibrated on parameter observations; for the model states, the covariance is estimated from a number of forward simulations of the PDE model corresponding to realizations of the parameters drawn from their KLE. We apply the proposed approach to identifying heterogeneous log-diffusion coefficients in diffusion equations from spatially sparse measurements of the log-diffusion coefficient and the solution of the diffusion equation. We find that the proposed approach compares favorably against state-of-the-art point estimates such as maximum a posteriori estimation and physics-informed neural networks.


title: Let's Get Dirty: GAN Based Data Augmentation for Camera Lens Soiling Detection in Autonomous Driving

abstract: Wide-angle fisheye cameras are commonly used in automated driving for parking and low-speed navigation tasks. Four of such cameras form a surround-view system that provides a complete and detailed view of the vehicle. These cameras are directly exposed to harsh environmental settings and can get soiled very easily by mud, dust, water, frost. Soiling on the camera lens can severely degrade the visual perception algorithms, and a camera cleaning system triggered by a soiling detection algorithm is increasingly being deployed. While adverse weather conditions, such as rain, are getting attention recently, there is only limited work on general soiling. The main reason is the difficulty in collecting a diverse dataset as it is a relatively rare event. We propose a novel GAN based algorithm for generating unseen patterns of soiled images. Additionally, the proposed method automatically provides the corresponding soiling masks eliminating the manual annotation cost. Augmentation of the generated soiled images for training improves the accuracy of soiling detection tasks significantly by 18% demonstrating its usefulness. The manually annotated soiling dataset and the generated augmentation dataset will be made public. We demonstrate the generalization of our fisheye trained GAN model on the Cityscapes dataset. We provide an empirical evaluation of the degradation of the semantic segmentation algorithm with the soiled data.


title: A Verified Optimizer for Quantum Circuits

abstract: We present VOQC, the first fully verified optimizer for quantum circuits, written using the Coq proof assistant. Quantum circuits are expressed as programs in a simple, low-level language called SQIR, a simple quantum intermediate representation, which is deeply embedded in Coq. Optimizations and other transformations are expressed as Coq functions, which are proved correct with respect to a semantics of SQIR programs. SQIR uses a semantics of matrices of complex numbers, which is the standard for quantum computation, but treats matrices symbolically in order to reason about programs that use an arbitrary number of quantum bits. SQIR's careful design and our provided automation make it possible to write and verify a broad range of optimizations in VOQC, including full-circuit transformations from cutting-edge optimizers.


title: Quality Selection in Two-Sided Markets: A Constrained Price Discrimination Approach

abstract: Online platforms collect rich information about participants and then share some of this information back with them to improve market outcomes. In this paper we study the following information disclosure problem in two-sided markets: If a platform wants to maximize revenue, which sellers should the platform allow to participate, and how much of its available information about participating sellers' quality should the platform share with buyers? We study this information disclosure problem in the context of two distinct two-sided market models: one in which the platform chooses prices and the sellers choose quantities (similar to ride-sharing), and one in which the sellers choose prices (similar to e-commerce). Our main results provide conditions under which simple information structures commonly observed in practice, such as banning certain sellers from the platform while not distinguishing between participating sellers, maximize the platform's revenue. The platform's information disclosure problem naturally transforms into a constrained price discrimination problem where the constraints are determined by the equilibrium outcomes of the specific two-sided market model being studied. We analyze this constrained price discrimination problem to obtain our structural results.


title: Multiple Anchor Learning for Visual Object Detection

abstract: Classification and localization are two pillars of visual object detectors. However, in CNN-based detectors, these two modules are usually optimized under a fixed set of candidate (or anchor) bounding boxes. This configuration significantly limits the possibility to jointly optimize classification and localization. In this paper, we propose a Multiple Instance Learning (MIL) approach that selects anchors and jointly optimizes the two modules of a CNN-based object detector. Our approach, referred to as Multiple Anchor Learning (MAL), constructs anchor bags and selects the most representative anchors from each bag. Such an iterative selection process is potentially NP-hard to optimize. To address this issue, we solve MAL by repetitively depressing the confidence of selected anchors by perturbing their corresponding features. In an adversarial selection-depression manner, MAL not only pursues optimal solutions but also fully leverages multiple anchors/features to learn a detection model. Experiments show that MAL improves the baseline RetinaNet with significant margins on the commonly used MS-COCO object detection benchmark and achieves new state-of-the-art detection performance compared with recent methods.


title: Heavy Quarkonium in Extreme Conditions

abstract: In this report we review recent progress achieved in the understanding of heavy quarkonium under extreme conditions from a theory perspective. Its focus lies both on quarkonium properties in thermal equilibrium, as well as recent developments towards a genuine real-time description, valid also out-of-equilibrium. We will give an overview of the theory tools developed and deployed over the last decade, including effective field theories, lattice field theory simulations, modern methods for spectral reconstructions and the the open-quantum systems paradigm. The report will discuss in detail the concept of quarkonium melting, providing the reader with a contemporary perspective. In order to judge where future progress is needed we will also discuss recent results from experiments and phenomenological modeling of quarkonium in relativistic heavy-ion collisions.


title: Deep Model Compression Via Two-Stage Deep Reinforcement Learning

abstract: Besides accuracy, the model size of convolutional neural networks (CNN) models is another important factor considering limited hardware resources in practical applications. For example, employing deep neural networks on mobile systems requires the design of accurate yet fast CNN for low latency in classification and object detection. To fulfill the need, we aim at obtaining CNN models with both high testing accuracy and small size to address resource constraints in many embedded devices. In particular, this paper focuses on proposing a generic reinforcement learning-based model compression approach in a two-stage compression pipeline: pruning and quantization. The first stage of compression, i.e., pruning, is achieved via exploiting deep reinforcement learning (DRL) to co-learn the accuracy and the FLOPs updated after layer-wise channel pruning and element-wise variational pruning via information dropout. The second stage, i.e., quantization, is achieved via a similar DRL approach but focuses on obtaining the optimal bits representation for individual layers. We further conduct experimental results on CIFAR-10 and ImageNet datasets. For the CIFAR-10 dataset, the proposed method can reduce the size of VGGNet by 9x from 20.04MB to 2.2MB with a slight accuracy increase. For the ImageNet dataset, the proposed method can reduce the size of VGG-16 by 33x from 138MB to 4.14MB with no accuracy loss.


title: WISPR Imaging of a Pristine CME

abstract: The Wide-field Imager for Solar Probe (WISPR) on board the Parker Solar Probe (PSP) observed a CME on 2018 November 01, the first day of the initial PSP encounter. The speed of the CME, approximately 200-300 km s$^{-1}$ in the WISPR field of view, is typical of slow, streamer blowout CMEs. This event was also observed by the LASCO coronagraphs. WISPR and LASCO view remarkably similar structures that enable useful cross-comparison between the two data sets as well as stereoscopic imaging of the CME. Analysis is extended to lower heights by linking the white-light observations to EUV data from AIA, which reveal a structure that erupts more than a full day earlier before the CME finally gathers enough velocity to propagate outward. This EUV feature appears as a brightness enhancement in cooler temperatures such as 171 \AA, but as a cavity in nominal coronal temperatures such as 193 \AA. By comparing this circular, dark feature in 193 \AA \ to the dark, white-light cavity at the center of the eruption in WISPR and LASCO, it can be seen that this is one coherent structure that exists prior to the eruption in the low corona before entering the heliosphere and likely corresponds to the core of the magnetic flux rope. It is also believed that the relative weakness of the event contributed to the clarity of the flux rope in WISPR, as the CME did not experience impulsive forces or strong interaction with external structures that can lead to more complex structural evolution.


title: Compositional Temporal Visual Grounding of Natural Language Event Descriptions

abstract: Temporal grounding entails establishing a correspondence between natural language event descriptions and their visual depictions. Compositional modeling becomes central: we first ground atomic descriptions "girl eating an apple," "batter hitting the ball" to short video segments, and then establish the temporal relationships between the segments. This compositional structure enables models to recognize a wider variety of events not seen during training through recognizing their atomic sub-events. Explicit temporal modeling accounts for a wide variety of temporal relationships that can be expressed in language: e.g., in the description "girl stands up from the table after eating an apple" the visual ordering of the events is reversed, with first "eating an apple" followed by "standing up from the table." We leverage these observations to develop a unified deep architecture, CTG-Net, to perform temporal grounding of natural language event descriptions to videos. We demonstrate that our system outperforms prior state-of-the-art methods on the DiDeMo, Tempo-TL, and Tempo-HL temporal grounding datasets.


title: Metrizability of holonomy invariant projective deformation of sprays

abstract: In this paper, we consider projective deformation of the geodesic system of Finsler spaces by holonomy invariant functions: Starting by a Finsler spray $S$ and a holonomy invariant function $P$, we investigate the metrizability property of the projective deformation $\widetilde{S}=S-2\lambda P C$. We prove that for any holonomy invariant nontrivial function $P$ and for almost every value $\lambda\in R$, such deformation is not Finsler metrizable. We identify the cases where such deformation can lead to a metrizable spray: in these cases, the holonomy invariant function is necessarily one of the principal curvatures of the geodesic structure.


title: A Survey of Game Theoretic Approaches for Adversarial Machine Learning in Cybersecurity Tasks

abstract: Machine learning techniques are currently used extensively for automating various cybersecurity tasks. Most of these techniques utilize supervised learning algorithms that rely on training the algorithm to classify incoming data into different categories, using data encountered in the relevant domain. A critical vulnerability of these algorithms is that they are susceptible to adversarial attacks where a malicious entity called an adversary deliberately alters the training data to misguide the learning algorithm into making classification errors. Adversarial attacks could render the learning algorithm unsuitable to use and leave critical systems vulnerable to cybersecurity attacks. Our paper provides a detailed survey of the state-of-the-art techniques that are used to make a machine learning algorithm robust against adversarial attacks using the computational framework of game theory. We also discuss open problems and challenges and possible directions for further research that would make deep machine learning-based systems more robust and reliable for cybersecurity tasks.


title: Extending the Morphological Hit-or-Miss Transform to Deep Neural Networks

abstract: While most deep learning architectures are built on convolution, alternative foundations like morphology are being explored for purposes like interpretability and its connection to the analysis and processing of geometric structures. The morphological hit-or-miss operation has the advantage that it takes into account both foreground and background information when evaluating target shape in an image. Herein, we identify limitations in existing hit-or-miss neural definitions and we formulate an optimization problem to learn the transform relative to deeper architectures. To this end, we model the semantically important condition that the intersection of the hit and miss structuring elements (SEs) should be empty and we present a way to express Don't Care (DNC), which is important for denoting regions of an SE that are not relevant to detecting a target pattern. Our analysis shows that convolution, in fact, acts like a hit-miss transform through semantic interpretation of its filter differences. On these premises, we introduce an extension that outperforms conventional convolution on benchmark data. Quantitative experiments are provided on synthetic and benchmark data, showing that the direct encoding hit-or-miss transform provides better interpretability on learned shapes consistent with objects whereas our morphologically inspired generalized convolution yields higher classification accuracy. Last, qualitative hit and miss filter visualizations are provided relative to single morphological layer.


title: The effect of task and training on intermediate representations in convolutional neural networks revealed with modified RV similarity analysis

abstract: Centered Kernel Alignment (CKA) was recently proposed as a similarity metric for comparing activation patterns in deep networks. Here we experiment with the modified RV-coefficient (RV2), which has very similar properties as CKA while being less sensitive to dataset size. We compare the representations of networks that received varying amounts of training on different layers: a standard trained network (all parameters updated at every step), a freeze trained network (layers gradually frozen during training), random networks (only some layers trained), and a completely untrained network. We found that RV2 was able to recover expected similarity patterns and provide interpretable similarity matrices that suggested hypotheses about how representations are affected by different training recipes. We propose that the superior performance achieved by freeze training can be attributed to representational differences in the penultimate layer. Our comparisons of random networks suggest that the inputs and targets serve as anchors on the representations in the lowest and highest layers.


title: A Catalog of Emission-Line Galaxies from the Faint Infrared Grism Survey: Studying Environmental Influence on Star Formation

abstract: We present a catalog of 208 $0.3 < z < 2.1$ Emission Line Galaxies (ELG) selected from 1D slitless spectroscopy obtained using Hubble's WFC3 G102 grism, as part of the Faint Infrared Grism Survey (FIGS). We identify ELG candidates by searching for significant peaks in all continuum-subtracted G102 spectra, and, where possible, confirm candidates by identifying consistent emission lines in other available spectra or with published spectroscopic redshifts. We provide derived emission line fluxes and errors, redshifts, and equivalent widths (EW) for H$\alpha$ $\lambda6563$, [OIII]$\lambda\lambda4959,5007$, and [OII]$\lambda\lambda3727$ emission lines, for emission line galaxies down to AB(F105W) $ > 28$ and $> 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ line flux. We use the resulting line catalog to investigate a possible relationship between line emission and a galaxy's environment. We use 7th-nearest-neighbor distances to investigate the typical surroundings of ELGs compared to non-ELGs, and we find that [OIII] emitters are preferentially found at intermediate galaxy densities near galaxy groups. We characterize these ELGs in terms of the galaxy specific star formation rate (SSFR) versus stellar mass, and find no significant influence of environment on that relation. We calculate star formation rates (SFR), and find no dependence of SFR on local galaxy surface density for $0.3 < z < 0.8$ H$\alpha$ emitters and for $0.8<z<1.3$ [OIII] emitters. We find similar rates of close-pair interaction between ELGs and non-ELGs. For galaxy surface densities $\Sigma \leq 30$ Mpc$^{-2}$, we find no consistent effect of environment on star formation.


title: Correspondent Banking Networks: Theory and Experiment

abstract: We employ the mathematical programming approach in conjunction with the graph theory to study the structure of correspondent banking networks. Optimizing the network requires decisions to be made to onboard, terminate or restrict the bank relationships to optimize the size and overall risk of the network. This study provides theoretical foundation to detect the components, the removal of which does not affect some key properties of the network such as connectivity and diameter. We find that the correspondent banking networks have a feature we call k-accessibility, which helps to drastically reduce the computational burden required for finding the above mentioned components. We prove a number of fundamental theorems related to k-accessible directed graphs, which should be also applicable beyond the particular problem of financial networks. The theoretical findings are verified through the data from a large international bank.


title: Evaluation Metrics for Item Recommendation under Sampling

abstract: The task of item recommendation requires ranking a large catalogue of items given a context. Item recommendation algorithms are evaluated using ranking metrics that depend on the positions of relevant items. To speed up the computation of metrics, recent work often uses sampled metrics where only a smaller set of random items and the relevant items are ranked. This paper investigates sampled metrics in more detail and shows that sampled metrics are inconsistent with their exact version. Sampled metrics do not persist relative statements, e.g., 'algorithm A is better than B', not even in expectation. Moreover the smaller the sampling size, the less difference between metrics, and for very small sampling size, all metrics collapse to the AUC metric.


title: A New Likely Redback Millisecond Pulsar Binary with a Massive Neutron Star: 4FGL J2333.1--5527

abstract: We present the discovery of a likely new redback millisecond pulsar binary associated with the \emph{Fermi} $\gamma$-ray source 4FGL J2333.1--5527. Using optical photometric and spectroscopic observations from the SOAR telescope, we identify a low-mass, main sequence-like companion in a 6.9-hr, highly inclined orbit around a suspected massive neutron star primary. Archival XMM-Newton X-ray observations show this system has a hard power-law spectrum $\Gamma = 1.6\pm0.3$ and $L_X \sim 5 \times 10^{31}$ erg s$^{-1}$, consistent with redback millisecond pulsar binaries. Our data suggest that for secondary masses typical of redbacks, the mass of the neutron star is likely well in excess of $\sim1.4\,M_{\odot}$, but future timing of the radio pulsar is necessary to bolster this tentative conclusion. This work shows that a bevy of nearby compact binaries still await discovery, and that unusually massive neutron stars continue to be common in redbacks.


title: Gaussian graphical models with toric vanishing ideals

abstract: Gaussian graphical models are semi-algebraic subsets of the cone of positive definite covariance matrices. They are widely used throughout natural sciences, computational biology and many other fields. Computing the vanishing ideal of the model gives us an implicit description of the model. In this paper, we resolve two conjectures of Sturmfels and Uhler from \cite{BS n CU}. In particular, we characterize those graphs for which the vanishing ideal of the Gaussian graphical model is generated in degree $1$ and $2$. These turn out to be the Gaussian graphical models whose ideals are toric ideals, and the resulting graphs are the $1$-clique sums of complete graphs.


title: On the asymptotic behavior of the $q$-analog of Kostant's partition function

abstract: Kostant's partition function counts the number of distinct ways to express a weight of a classical Lie algebra $\mathfrak{g}$ as a sum of positive roots of $\mathfrak{g}$. We refer to each of these expressions as decompositions of a weight. Our main result considers an infinite family of weights, irrespective of Lie type, for which we establish a closed formula for the $q$-analog of Kostant's partition function and then prove that the (normalized) distribution of the number of positive roots in the decomposition of any of these weights converges to a Gaussian distribution as the rank of the Lie algebra goes to infinity. We also extend these results to the highest root of the classical Lie algebras and we end our analysis with some directions for future research.


title: Masur-Veech volumes and intersection theory: the principal strata of quadratic differentials

abstract: We describe a conjectural formula via intersection numbers for the Masur-Veech volumes of strata of quadratic differentials with prescribed zero orders, and we prove the formula for the case when the zero orders are odd. For the principal strata of quadratic differentials with simple zeros, the formula reduces to compute the top Segre class of the quadratic Hodge bundle, which can be further simplified to certain linear Hodge integrals. An appendix proves that the intersection of this class with $\psi$-classes can be computed by Eynard-Orantin topological recursion. As applications, we analyze numerical properties of Masur-Veech volumes, area Siegel-Veech constants and sums of Lyapunov exponents of the principal strata for fixed genus and varying number of zeros, which settles the corresponding conjectures due to Grivaux-Hubert, Fougeron, and elaborated in [the7]. We also describe conjectural formulas for area Siegel-Veech constants and sums of Lyapunov exponents for arbitrary affine invariant submanifolds, and verify them for the principal strata.


title: Signatures of Clouds in Hot Jupiter Atmospheres: Modeled High Resolution Emission Spectra from 3D General Circulation Models

abstract: Observations of scattered light and thermal emission from hot Jupiter exoplanets have suggested the presence of inhomogeneous aerosols in their atmospheres. 3D general circulation models (GCMs) that attempt to model the effects of aerosols have been developed to understand the physical processes that underlie their dynamical structures. In this work, we investigate how different approaches to aerosol modeling in GCMs of hot Jupiters affect high-resolution thermal emission spectra throughout the duration of the planet's orbit. Using results from a GCM with temperature-dependent cloud formation, we calculate spectra of a representative hot Jupiter with different assumptions regarding the vertical extent and thickness of clouds. We then compare these spectra to models in which clouds are absent or simply post-processed (i.e., added subsequently to the completed clear model). We show that the temperature-dependent treatment of clouds in the GCM produces high-resolution emission spectra that are markedly different from the clear and post-processed cases -- both in the continuum flux levels and line profiles -- and that increasing the vertical extent and thickness of clouds leads to bigger changes in these features. We evaluate the net Doppler shifts of the spectra induced by global winds and the planet's rotation and show that they are strongly phase-dependent, especially for models with thicker and more extended clouds. This work further demonstrates the importance of radiative feedback in cloudy atmospheric models of hot Jupiters, as this can have a significant impact on interpreting spectroscopic observations of exoplanet atmospheres.


title: A Computational Approach to Classifying Low Rank Modular Categories

abstract: This paper introduces a computational approach to classifying low rank modular categories up to their modular data. The modular data of a modular category is a pair of matrices, $(S,T)$. Virtually all the numerical information of the category is contained within or derived from the modular data. The modular data satisfy a variety of criteria that Bruillard, Ng, Rowell, and Wang call the admissibility criteria. Of note is the Galois group of the $S$ matrix is an abelian group that acts faithfully on the columns of the eigenvalue matrix, $s = (\frac{S_{ij}}{S_{0j}})$. This gives an injection from Gal$(\mathbb{Q}(S),\mathbb{Q}) \to $ Sym$_r$, where $r$ is the rank of the category. Our approach begins by listing all the possible abelian subgroups of Sym$_6$ and building all the possible modular data for each group. We run each set of modular data through a series of Gr\"obner basis calculations until we either find a contradiction or solve for the modular data.


title: A Unified Switching System Perspective and O.D.E. Analysis of Q-Learning Algorithms

abstract: In this paper, we introduce a unified framework for analyzing a large family of Q-learning algorithms, based on switching system perspectives and ODE-based stochastic approximation. We show that the nonlinear ODE models associated with these Q-learning algorithms can be formulated as switched linear systems, and analyze their asymptotic stability by leveraging existing switching system theories. Our approach provides the first O.D.E. analysis of the asymptotic convergence of various Q-learning algorithms, including asynchronous Q-learning and averaging Q-learning. We also extend the approach to analyze Q-learning with linear function approximation and derive a new sufficient condition for its convergence.


title: Global-scale Observations and Modeling of Far-Ultraviolet Airglow During Twilight

abstract: The NASA Global-scale Observations of the Limb and Disk (GOLD) ultraviolet imaging spectrograph performs observations of upper-atmosphere airglow from the sunlit disk and limb of the Earth, in order to infer quantities such as the composition and temperature of the thermosphere. To interpret the measurements, the observational and solar illumination geometry must be considered. We use forward models of upper atmosphere density and composition, photoelectron impact, airglow emissions, radiative transfer, and line-of-sight integration, to describe the expected observations, and here test those calculations against observations near the terminator, and near the limb. On the night side of the terminator, broad regions of faint airglow are seen, particularly near the winter solstice. These are caused by photoelectrons that were transported along field lines from magnetically conjugate areas in the other hemisphere, where those areas are still illuminated. We perform model calculations to demonstrate that this process is the source of the emission, and obtain good agreement with its morphology and intensity. In some regions, the observed emissions are not as intense as the model simulations. Some of the reductions in electron flux are explained by changes in magnetic field strength; in other cases, particularly at high magnetic latitude, the cause is unknown, but must occur along extended field lines as they reach into the plasma sheet.


title: Multivariate Rational Approximation

abstract: We present two approaches for computing rational approximations to multivariate functions, motivated by their effectiveness as surrogate models for high-energy physics (HEP) applications. Our first approach builds on the Stieltjes process to efficiently and robustly compute the coefficients of the rational approximation. Our second approach is based on an optimization formulation that allows us to include structural constraints on the rational approximation, resulting in a semi-infinite optimization problem that we solve using an outer approximation approach. We present results for synthetic and real-life HEP data, and we compare the approximation quality of our approaches with that of traditional polynomial approximations.


title: The impact of automation and connectivity on traffic flow and CO2 emissions. A detailed microsimulation study

abstract: Many of the anticipated advantages of connected and automated vehicles or automated vehicles without connectivity (CAVs and AVs respectively) on congestion and energy consumption are questionable. Some studies provide quantitative answers to the above questions through microsimulation but they systematically ignore the realistic simulation of vehicle dynamics, driver behaviour or instantaneous emissions estimates, mostly due to the overall increased complexity of the transport systems and the need for low computational demand on large-scale simulations. However, recent studies question the capability of common car-following models to produce realistic vehicle dynamics or driving behaviour, which directly impacts emissions estimations as well. This work presents a microsimulation study that contributes on the topic, using a scenario-based approach to give insights regarding the impact of CAVs and AVs on the evolution of emissions over a highway network. The motivation here is to answer whether the different driving behaviours produce significant differences in emissions during rush hours, and how significant is the impact of detailed vehicle dynamics simulation and instantaneous emissions in the outcome.


title: Geometry of martensite needles in shape memory alloys

abstract: We study the geometry of needle-shaped domains in shape-memory alloys. Needle-shaped domains are ubiquitously found in martensites around macroscopic interfaces between regions which are laminated in different directions, or close to macroscopic austenite/twinned-martensite interfaces. Their geometry results from the interplay of the local nonconvexity of the effective energy density with long-range (linear) interactions mediated by the elastic strain field, and is up to now poorly understood. We present a two-dimensional shape optimization model based on finite elasticity and discuss its numerical solution. Our results indicate that the tapering profile of the needles can be understood within finite elasticity, but not with linearized elasticity. The resulting tapering and bending reproduce the main features of experimental observations on NiAl.


title: Optical phased arrays for LIDAR: beam steering via tunable plasmonic metasurfaces

abstract: Controlling the phase and amplitude of light emitted by the elements (i.e., pixels) of an optical phased array is of paramount importance to realizing dynamic beam steering for LIDAR applications. In this paper, we propose a plasmonic pixel composed of a metallic nanoantenna covered by a thin oxide layer, and a conductive oxide, e.g., ITO, for use in a reflectarray metasurface. By considering voltage biasing of the nanoantenna via metallic connectors, and exploiting the carrier refraction effect in the metal-oxide-semiconductor capacitor in the accumulation and depletion regions, our simulations predict control of the reflection coefficient phase over a range $>330^{\circ}$ with a nearly constant magnitude. We discuss the physical mechanism underlying the optical response, the effect of the connectors, and propose strategies to maximize the magnitude of the reflection coefficient and to achieve dual-band operation. The suitability of our plasmonic pixel design for beam steering in LIDAR is demonstrated via 3D-FDTD simulations.


title: Enhancing Stratospheric Weather Analyses and Forecasts by Deploying Sensors from a Weather Balloon

abstract: The ability to analyze and forecast stratospheric weather conditions is fundamental to addressing climate change. However, our capacity to collect data in the stratosphere is limited by sparsely deployed weather balloons. We propose a framework to collect stratospheric data by releasing a contrail of tiny sensor devices as a weather balloon ascends. The key machine learning challenges are determining when and how to deploy a finite collection of sensors to produce a useful data set. We decide when to release sensors by modeling the deviation of a forecast from actual stratospheric conditions as a Gaussian process. We then implement a novel hardware system that is capable of optimally releasing sensors from a rising weather balloon. We show that this data engineering framework is effective through real weather balloon flights, as well as simulations.


title: On coefficients of Poincar\'e series and single-valued periods of modular forms

abstract: We prove that the field generated by the Fourier coefficients of weakly holomorphic Poincar\'e series of a given level $\Gamma_0(N)$ and weight $k\ge 2$ coincides with the field generated by the single-valued periods of a certain motive attached to $\Gamma_0(N)$. This clarifies the arithmetic nature of such Fourier coefficients and generalises previous formulas of Brown and Acres--Broadhurst giving explicit series expansions for the single-valued periods of some modular forms. Our proof is based on Bringmann--Ono's construction of harmonic lifts of Poincar\'e series.


title: Smoothing the gap between NP and ER

abstract: We study algorithmic problems that belong to the complexity class of the existential theory of the reals (ER). A problem is ER-complete if it is as hard as the problem ETR and if it can be written as an ETR formula. Traditionally, these problems are studied in the real RAM, a model of computation that assumes that the storage and comparison of real-valued numbers can be done in constant space and time, with infinite precision. The complexity class ER is often called a real RAM analogue of NP, since the problem ETR can be viewed as the real-valued variant of SAT. In this paper we prove a real RAM analogue to the Cook-Levin theorem which shows that ER membership is equivalent to having a verification algorithm that runs in polynomial-time on a real RAM. This gives an easy proof of ER-membership, as verification algorithms on a real RAM are much more versatile than ETR-formulas. We use this result to construct a framework to study ER-complete problems under smoothed analysis. We show that for a wide class of ER-complete problems, its witness can be represented with logarithmic input-precision by using smoothed analysis on its real RAM verification algorithm. This shows in a formal way that the boundary between NP and ER (formed by inputs whose solution witness needs high input-precision) consists of contrived input. We apply our framework to well-studied ER-complete recognition problems which have the exponential bit phenomenon such as the recognition of realizable order types or the Steinitz problem in fixed dimension.


title: Angular Visual Hardness

abstract: Recent convolutional neural networks (CNNs) have led to impressive performance but often suffer from poor calibration. They tend to be overconfident, with the model confidence not always reflecting the underlying true ambiguity and hardness. In this paper, we propose angular visual hardness (AVH), a score given by the normalized angular distance between the sample feature embedding and the target classifier to measure sample hardness. We validate this score with an in-depth and extensive scientific study, and observe that CNN models with the highest accuracy also have the best AVH scores. This agrees with an earlier finding that state-of-art models improve on the classification of harder examples. We observe that the training dynamics of AVH is vastly different compared to the training loss. Specifically, AVH quickly reaches a plateau for all samples even though the training loss keeps improving. This suggests the need for designing better loss functions that can target harder examples more effectively. We also find that AVH has a statistically significant correlation with human visual hardness. Finally, we demonstrate the benefit of AVH to a variety of applications such as self-training for domain adaptation and domain generalization.


title: Natural Alpha Embeddings

abstract: Learning an embedding for a large collection of items is a popular approach to overcome the computational limitations associated to one-hot encodings. The aim of item embedding is to learn a low dimensional space for the representations, able to capture with its geometry relevant features or relationships for the data at hand. This can be achieved for example by exploiting adjacencies among items in large sets of unlabelled data. In this paper we interpret in an Information Geometric framework the item embeddings obtained from conditional models. By exploiting the $\alpha$-geometry of the exponential family, first introduced by Amari, we introduce a family of natural $\alpha$-embeddings represented by vectors in the tangent space of the probability simplex, which includes as a special case standard approaches available in the literature. A typical example is given by word embeddings, commonly used in natural language processing, such as Word2Vec and GloVe. In our analysis, we show how the $\alpha$-deformation parameter can impact on standard evaluation tasks.


title: A high-order discontinuous Galerkin method for nonlinear sound waves

abstract: We propose a high-order discontinuous Galerkin scheme for nonlinear acoustic waves on polytopic meshes. To model sound propagation with and without losses, we use Westervelt's nonlinear wave equation with and without strong damping. Challenges in the numerical analysis lie in handling the nonlinearity in the model, which involves the derivatives in time of the acoustic velocity potential, and in preventing the equation from degenerating. We rely in our approach on the Banach fixed-point theorem combined with a stability and convergence analysis of a linear wave equation with a variable coefficient in front of the second time derivative. By doing so, we derive an a priori error estimate for Westervelt's equation in a suitable energy norm for the polynomial degree $p \geq 2$. Numerical experiments carried out in two-dimensional settings illustrate the theoretical convergence results. In addition, we demonstrate efficiency of the method in a three-dimensional domain with varying medium parameters, where we use the discontinuous Galerkin approach in a hybrid way.


title: Bound states and the potential parameter spectrum

abstract: In this article, we answer the following question: If the wave equation possesses bound states but it is exactly solvable for only a single non-zero energy, can we find all bound state solutions (energy spectrum and associated wavefunctions)? To answer this question, we use the "tridiagonal representation approach" to solve the wave equation at the given energy by expanding the wavefunction in a series of energy-dependent square integrable basis functions in configuration space. The expansion coefficients satisfy a three-term recursion relation, which is solved in terms of orthogonal polynomials. Depending on the selected energy we show that one of the potential parameters must assume a value from within a discrete set called the "potential parameter spectrum" (PPS). This discrete set is obtained from the spectrum of the above polynomials and can be either a finite or an infinite discrete set. Inverting the relation between the energy and the PPS gives the bound state energy spectrum. Therefore, the answer to the above question is affirmative.


title: Sub-linear RACE Sketches for Approximate Kernel Density Estimation on Streaming Data

abstract: Kernel density estimation is a simple and effective method that lies at the heart of many important machine learning applications. Unfortunately, kernel methods scale poorly for large, high dimensional datasets. Approximate kernel density estimation has a prohibitively high memory and computation cost, especially in the streaming setting. Recent sampling algorithms for high dimensional densities can reduce the computation cost but cannot operate online, while streaming algorithms cannot handle high dimensional datasets due to the curse of dimensionality. We propose RACE, an efficient sketching algorithm for kernel density estimation on high-dimensional streaming data. RACE compresses a set of N high dimensional vectors into a small array of integer counters. This array is sufficient to estimate the kernel density for a large class of kernels. Our sketch is practical to implement and comes with strong theoretical guarantees. We evaluate our method on real-world high-dimensional datasets and show that our sketch achieves 10x better compression compared to competing methods.


title: Magnetotransport and internodal tunnelling in Weyl semimetals

abstract: Internodal dynamics of quasiparticles in Weyl semimetals manifest themselves in hydrodynamic, transport and thermodynamic phenomena and are essential for potential valleytronic applications of these systems. In an external magnetic field, coherent quasiparticle tunnelling between the nodes modifies the quasiparticle dispersion and, in particular, opens gaps in the dispersion of quasiparticles at the zeroth Landau level. We study magnetotransport in a Weyl semimetal taking into account mechanisms of quasiparticle scattering both affected by such gaps and independent of them. We compute the longitudal resistivity of a disordered Weyl semimetal with two nodes in a strong magnetic field microscopically and demonstrate that in a broad range of magnetic fields it has a strong angular dependence $\rho(\eta)\propto C_1+C_2 \cos^2\eta$, where $\eta$ is the angle between the field and the separation between the nodes in momentum space. The first term is determined by the coherent internodal tunnelling and is important only at angles $\eta$ close to $\pi/2$. This contribution depends exponentially on the magnetic field, $\propto \exp\left(-B_0/B\right)$. The second term is weakly dependent on the magnetic field for realistic concentrations of the impurities in a broad interval of fields.

categories: cond-mat.mes-hall cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.str-el


title: Gobi: WebAssembly as a Practical Path to Library Sandboxing

abstract: Software based fault isolation (SFI) is a powerful approach to reduce the impact of security vulnerabilities in large C/C++ applications like Firefox and Apache. Unfortunately, practical SFI tools have not been broadly available. Developing SFI toolchains are a significant engineering challenge. Only in recent years have browser vendors invested in building production quality SFI tools like Native Client (NaCl) to sandbox code. Further, without committed support, these tools are not viable, e.g. NaCl has been discontinued, orphaning projects that relied on it. WebAssembly (Wasm) offers a promising solution---it can support high performance sandboxing and has been embraced by all major browser vendors---thus seems to have a viable future. However, Wasm presently only offers a solution for sandboxing mobile code. Providing SFI for native application, such as C/C++ libraries requires additional steps. To reconcile the different worlds of Wasm on the browser and native platforms, we present Gobi. Gobi is a system of compiler changes and runtime support that can sandbox normal C/C++ libraries with Wasm---allowing them to be compiled and linked into native applications. Gobi has been tested on libjpeg, libpng, and zlib. Based on our experience developing Gobi, we conclude with a call to arms to the Wasm community and SFI research community to make Wasm based module sandboxing a first class use case and describe how this can significantly benefit both communities. Addendum: This short paper was originally written in January of 2019. Since then, the implementation and design of Gobi has evolved substantially as some of the issues raised in this paper have been addressed by the Wasm community. Nevertheless, several challenges still remain. We have thus left the paper largely intact and only provide a brief update on the state of Wasm tooling as of November 2019 in the last section.


title: The onset of electron-only reconnection

abstract: Motivated by recent observations of "electron-only" magnetic reconnection, without an ion-scale sheet or ion outflows, in both the Earth's magnetosheath and in numerical simulations, we study the formation and reconnection of electron-scale current sheets at low plasma beta. We first show that ideal sheets collapse to thicknesses much smaller than the ion scales, by deriving an appropriate analogue of the Chapman-Kendall collapse solution. Second, we show that, in practice, reconnection onset happens in these collapsing sheets once they reach a critical aspect ratio, because the tearing instability then becomes faster than their collapse timescale. We show that this can happen for sheet thicknesses larger than the ion scale or at only a few times the electron scale, depending on plasma parameters and the aspect ratio of the collapsing structure, thereby unifying the usual picture of ion-coupled reconnection and the new regime of electron-only reconnection. %and thus the usual picture of ion-coupled reconnection and the new regime of electron-only reconnection can be unified. We derive relationships between plasma beta, ion-to-electron temperature ratio, the aspect ratio, electron outflow velocity, and the final thickness of the sheets, and thus determine under what circumstances electron-scale sheets form and reconnect.


title: A new algorithm to classify chiral polytopes with a given automorphism group

abstract: We present a new algorithm to compute all the chiral polytopes that have a given group $G$ as full automorphism group. This algorithm uses a new set of generators that characterize the group, all of them except one being involutions. It permits to compute all chiral polytopes of groups that were previously unreachable by other known algorithms.


title: Simplified Action Decoder for Deep Multi-Agent Reinforcement Learning

abstract: In recent years we have seen fast progress on a number of benchmark problems in AI, with modern methods achieving near or super human performance in Go, Poker and Dota. One common aspect of all of these challenges is that they are by design adversarial or, technically speaking, zero-sum. In contrast to these settings, success in the real world commonly requires humans to collaborate and communicate with others, in settings that are, at least partially, cooperative. In the last year, the card game Hanabi has been established as a new benchmark environment for AI to fill this gap. In particular, Hanabi is interesting to humans since it is entirely focused on theory of mind, i.e., the ability to effectively reason over the intentions, beliefs and point of view of other agents when observing their actions. Learning to be informative when observed by others is an interesting challenge for Reinforcement Learning (RL): Fundamentally, RL requires agents to explore in order to discover good policies. However, when done naively, this randomness will inherently make their actions less informative to others during training. We present a new deep multi-agent RL method, the Simplified Action Decoder (SAD), which resolves this contradiction exploiting the centralized training phase. During training SAD allows other agents to not only observe the (exploratory) action chosen, but agents instead also observe the greedy action of their team mates. By combining this simple intuition with best practices for multi-agent learning, SAD establishes a new SOTA for learning methods for 2-5 players on the self-play part of the Hanabi challenge. Our ablations show the contributions of SAD compared with the best practice components. All of our code and trained agents are available at https://github.com/facebookresearch/Hanabi_SAD.


title: Explicit form of the effective evolution equation for the randomly forced Schr\"{o}dinger equation with quadratic nonlinearity

abstract: An effective equation describes a weakly nonlinear wave field evolution governed by nonlinear dispersive PDEs \emph{via} the set of its resonances in an arbitrary big but finite domain in the Fourier space. We consider the Schr\"{o}dinger equation with quadratic nonlinearity including small external random forcing/dissipation. An effective equation is deduced explicitly for each case of monomial quadratic nonlinearities $ u^2, \, \bar{u}u, \, \bar{u}^2$ and the sets of resonance clusters are studied. In particular, we demonstrate that the nonlinearity $\bar{u}^2$ generates no 3-wave resonances and its effective equation is degenerate while in two other cases the sets of resonances are not empty. Possible implications for wave turbulence theory are briefly discussed.


title: Hierarchical Indian Buffet Neural Networks for Bayesian Continual Learning

abstract: We place an Indian Buffet process (IBP) prior over the structure of a Bayesian Neural Network (BNN), thus allowing the complexity of the BNN to increase and decrease automatically. We further extend this model such that the prior on the structure of each hidden layer is shared globally across all layers, using a Hierarchical-IBP (H-IBP). We apply this model to the problem of resource allocation in Continual Learning (CL) where new tasks occur and the network requires extra resources. Our model uses online variational inference with reparameterisation of the Bernoulli and Beta distributions, which constitute the IBP and H-IBP priors. As we automatically learn the number of weights in each layer of the BNN, overfitting and underfitting problems are largely overcome. We show empirically that our approach offers a competitive edge over existing methods in CL.