title: CO Multi-line Imaging of Nearby Galaxies (COMING). III. Dynamical effect on molecular gas density and star formation in the barred spiral galaxy NGC 4303

abstract: We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the volume densities of molecular gas $n(\rm{H}_2)$ in the eight different regions within the galactic disk with CO data combined with archival far-ultraviolet and 24 $\mu$m data. We confirmed that SFE in the bar region is lower by 39% than that in the spiral arms. Moreover, velocity-alignment stacking analysis was performed for the spectra in the individual regions. The integrated intensity ratios of $^{12}$CO to $^{13}$CO ($R_{12/13}$) range from 10 to 17 as the results of stacking. Fixing a kinetic temperature of molecular gas, $n(\rm{H}_2)$ was derived from $R_{12/13}$ via non-local thermodynamic equilibrium (non-LTE) analysis. The density $n(\rm{H}_2)$ in the bar is lower by 31-37% than that in the arms and there is a rather tight positive correlation between SFEs and $n(\rm{H}_2)$, with a correlation coefficient of $\sim 0.8$. Furthermore, we found a dependence of $n(\rm{H}_2)$ on the velocity dispersion of inter-molecular clouds ($\Delta V/ \sin i$). Specifically, $n(\rm{H}_2)$ increases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i < 100$ km s$^{-1}$. On the other hand, $n(\rm{H}_2)$ decreases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i > 100$ km s$^{-1}$. These relations indicate that the variations of SFE could be caused by the volume densities of molecular gas, and the volume densities could be governed by the dynamical influence such as cloud-cloud collisions, shear and enhanced inner-cloud turbulence.


title: Dynamic properties of superconductors: Anderson-Bogoliubov mode and Berry phase in BCS and BEC regimes

abstract: We analyze the evolution of the dynamics of a neutral s-wave superconductor between BCS and BEC regimes. We consider 2d case, when BCS-BEC crossover occurs already at weak coupling as a function of the ratio of the two scales -- the Fermi energy $E_F$ and the bound state energy for two fermions in a vacuum, $E_0$. BCS and BEC limits correspond to $E_F \gg E_0$ and $E_F \ll E_0$, respectively. The chemical potential $\mu = E_F-E_0$ changes the sign between the two regimes. We develop an approach to derive the leading terms in the expansion of the effective action in the spatial and time derivative of the slowly varying superconducting order parameter $\Delta (r, \tau)$, and express the action in terms of derivatives of the phase $\phi (r,\tau)$ of $\Delta (r, \tau) = \Delta e^{i\phi (r, \tau)}$. In the long wavelength limit the second gradients $(\nabla \phi)^2$ and ${\dot \phi}^2$, describe Bogoliubov-Anderson mode, which, as we find, does not change through BCS-BEC crossover. The effective action also contains first order gradient term, $i \pi A {\dot \phi}$, which is meaningful only if $\phi$ contains topological defects, such as vortices. We apply our approach to evaluate coefficient $A$ for a moving vortex. We find that $A_{vort}$ has two contributions, $A_{vort}= A_{vort,1} - A_{vort,2}$. One comes from the states away from the vortex core and has $A_{vort,1} = n/2$, where $n$ is the fermion density. The other comes from electronic states inside the vortex core and has $A_{vort,2} = -n_0/2$, where $n_0$ is the fermion density at the vortex core. We interpret this term as the reaction force of normal electrons to the vortex displacement in the limit when the spacing between energy levels is set to zero. The difference $(n-n_0)/2$ changes through the BEC-BCS crossover as $n_0$ nearly compensates $n$ in the BCS regime, but vanishes in the BEC regime.


title: Rediscovery of the Sixth Star Cluster in the Fornax Dwarf Spheroidal Galaxy

abstract: Since first noticed by Shapley in 1939, a faint object coincident with the Fornax dwarf spheroidal has long been discussed as a possible sixth globular cluster system. However, debate has continued over whether this overdensity is a statistical artifact or a blended galaxy group. In this Letter we demonstrate, using deep DECam imaging data, that this object is well resolved into stars and is a bona fide star cluster. The stellar overdensity of this cluster is statistically significant at the level of ~ 6 - 6.7 sigma in several different photometric catalogs including Gaia. Therefore, it is highly unlikely to be caused by random fluctuation. We show that Fornax 6 is a star cluster with a peculiarly low surface brightness and irregular shape, which may indicate a strong tidal influence from its host galaxy. The Hess diagram of Fornax 6 is largely consistent with that of Fornax field stars, but it appears to be slightly bluer. However, it is still likely more metal-rich than most of the globular clusters in the system. Faint clusters like Fornax 6 that orbit and potentially get disrupted in the centers of dwarf galaxies can prove crucial for constraining the dark matter distribution in Milky Way satellites.


title: Scattering mechanisms and electrical transport near an Ising nematic quantum critical point

abstract: Electrical transport properties near an electronic Ising-nematic quantum critical point in two dimensions are of both theoretical and experimental interest. In this work, we derive a kinetic equation valid in a broad regime near the quantum critical point using the memory matrix approach. The formalism is applied to study the effect of the critical fluctuations on the dc resistivity through different scattering mechanisms, including umklapp, impurity scattering, and electron-hole scattering in a compensated metal. We show that electrical transport in the quantum critical regime exhibits a rich behavior that depends sensitively on the scattering mechanism and the band structure. In the case of a single large Fermi surface, the resistivity due to umklapp scattering crosses over from $\rho\sim T^2$ at low temperature to sublinear at high temperature. The crossover temperature scales as $q_0^3$, where $q_0$ is the minimal wavevector for umklapp scattering. Impurity scattering leads to $\rho-\rho_0\sim T^\alpha$ ($\rho_0$ being the residual resistivity), where $\alpha$ is either larger than 2 if there is only a single Fermi sheet present, or 4/3 in the case of multiple Fermi sheets. Finally, in a perfectly compensated metal with an equal density of electrons and holes, the low temperature behavior depends strongly on the structure of "cold spots" on the Fermi surface, where the coupling between the quasiparticles and order parameter fluctuations vanishes by symmetry. In particular, for a system where cold spots are present on some (but not all) Fermi sheets, $\rho\sim T^{5/3}$. At higher temperatures there is a broad crossover regime where $\rho$ either saturates or $\rho\sim T$, depending on microscopic details. We discuss these results in the context of recent quantum Monte Carlo simulations of a metallic Ising nematic critical point, and experiments in certain iron-based superconductors.


title: The consequences of planetary migration on the minor bodies of the early Solar System

abstract: Pebble accretion is an efficient mechanism able to build up the core of the giant planets within the lifetime of the protoplanetary disc gas-phase. The core grows via this process until the protoplanet reaches its pebble isolation mass and starts to accrete gas. During the growth, the protoplanet undergoes a rapid, large-scale, inward migration due to the interactions with the gaseous protoplanetary disc. In our work, we investigate how this early migration would have affected the minor body populations in our solar system. In particular, we focus on the Jupiter Trojans and the Hildas asteroids. We found that a massive and eccentric Hilda group is captured during the migration from a region between 5 and 8 au and subsequently depleted during the late instability of the giant planets. Our simulations also show that inward migration of the giant planets always produces a Jupiter Trojans' leading swarm more populated than the trailing one, with a ratio comparable to the current observed Trojan asymmetry ratio. The in situ formation of Jupiter, on the other hand, produces symmetric leading/trailing swarms. The reason for the asymmetry is the relative drift between the migrating planet and the particles in the coorbital resonance. The capture happens during the growth of Jupiter's core and Trojan asteroids are afterwards carried along during the giant planet's migration to their final orbits. The asymmetry and eccentricity of the captured Trojans correspond well to observations, but their inclinations are near zero and their total mass is 3-4 orders of magnitude higher than the current population. Future modelling will be needed to understand whether the dynamical evolution of the Trojans over billions of years will raise the inclinations and deplete the masses to observed values.


title: Coordinate descent full configuration interaction

abstract: We develop an efficient algorithm, coordinate descent FCI (CDFCI), for the electronic structure ground state calculation in the configuration interaction framework. CDFCI solves an unconstrained non-convex optimization problem, which is a reformulation of the full configuration interaction eigenvalue problem, via an adaptive coordinate descent method with a deterministic compression strategy. CDFCI captures and updates appreciative determinants with different frequencies proportional to their importance. We show that CDFCI produces accurate variational energy for both static and dynamic correlation by benchmarking the binding curve of nitrogen dimer in the cc-pVDZ basis with $10^{-3}$ mHa accuracy. We also demonstrate the efficiency and accuracy of CDFCI for strongly correlated chromium dimer in the Ahlrichs VDZ basis and produces state-of-the-art variational energy.


title: Effective Theory for a Heavy Scalar Coupled to the SM via Vector-Like Quarks

abstract: We illustrate the application of the recently developed SCET$_{\rm BSM}$ framework in the context of a specific model, in which the Standard Model (SM) is supplemented by a heavy scalar $S$ and three generations of heavy, vector-like quarks $\Psi$. We construct the appropriate effective field theory for two-body decays of $S$ into SM particles. We explicitly compute the Wilson coefficients of the SCET$_{\rm BSM}$ operators appearing at leading and next-to-leading order (NLO) in an expansion in powers of $v/M_S$, as well as for a subset of operators arising at NNLO, retaining the full dependence on the ratio $M_S/M_\Psi$. For the phenomenologically most relevant decay channels of the heavy scalar, we study the impact of resummation effects of Sudakov logarithms on the decay rates.


title: Removing instabilities in the hierarchical equations of motion: exact and approximate projection approaches

abstract: The hierarchical equations of motion (HEOM) provide a numerically exact approach for computing the reduced dynamics of a quantum system linearly coupled to a bath. We have found that HEOM contains temperature-dependent instabilities that grow exponentially in time. In the case of continuous-bath models, these instabilities may be delayed to later times by increasing the hierarchy dimension; however, for systems coupled to discrete, non-dispersive modes, increasing the hierarchy dimension does little to alleviate the problem. We show that these instabilities can also be removed completely at a potentially much lower cost via projection onto the space of stable eigenmodes; furthermore, we find that for discrete-bath models at zero temperature, the remaining projected dynamics computed with few hierarchy levels are essentially identical to the exact dynamics that otherwise might require an intractably large number of hierarchy levels for convergence. Recognizing that computation of the eigenmodes might be prohibitive, e.g. for large or strongly-coupled models, we present a Prony filtration algorithm that may be useful as an alternative for accomplishing this projection when diagonalization is too costly. We present results demonstrating the efficacy of HEOM projected via diagonalization and Prony filtration. We also discuss issues associated with the nonnormality of HEOM.


title: Spectra of networks containing short loops

abstract: The spectrum of the adjacency matrix plays several important roles in the mathematical theory of networks and in network data analysis, for example in percolation theory, community detection, centrality measures, and the theory of dynamical systems on networks. A number of methods have been developed for the analytic computation of network spectra, but they typically assume that networks are locally tree-like, meaning that the local neighborhood of any node takes the form of a tree, free of short loops. Empirically observed networks, by contrast, often have many short loops. Here we develop an approach for calculating the spectra of networks with short loops using a message passing method. We give example applications to some previously studied classes of networks.


title: Categorical Saito theory, I: A comparison result

abstract: In this paper, we present an explicit cyclic minimal $A_\infty$ model for the category of matrix factorizations $\MF(W)$ of an isolated hypersurface singularity. The key observation is to use Kontsevich's deformation quantization technique. Pushing this idea further, we use the Tsygan formality map to obtain a comparison theorem that the categorical Variation of Semi-infinite Hodge Structure of $\MF(W)$ is isomorphic to Saito's original geometric construction in primitive form theory. An immediate corollary of this comparison result is that the analogue of Caldararu's conjecture holds for the category $\MF(W)$.


title: Gravitational Waves from Compact Dark Objects in Neutron Stars

abstract: Dark matter could be composed of compact dark objects (CDOs). We find that the oscillation of CDOs inside neutron stars can be a detectable source of gravitational waves (GWs). The GW strain amplitude depends on the mass of the CDO, and its frequency is typically in the range 3-5 kHz as determined by the central density of the star. In the best cases, LIGO may be sensitive to CDO masses greater than or of order $10^{-8}$ solar masses.


title: On the information content of the difference from hamiltonian evolution

abstract: A dissipative version of hamiltonian mechanics is proposed via a principle of minimal information content of the deviation from hamiltonian evolution. We show that we can cover viscosity, plasticity, damage and unilateral contact. This article continues arXiv:1807.10480


title: Topological Phase Transition in Superconductors with Mirror Symmetry

abstract: We provide analytical and numerical evidence that the attractive two-dimensional Kitaev model on a lattice with mirror symmetry demonstrates unusual 'intrinsic' phase at half filling. This phase emerges in the phase diagram at the boundary separating two topological superconductors with opposite Chern numbers and exists due to condensation of non-zero momentum Cooper pairs. Unlike Fulde-Ferrell-Larkin-Ovchinnikov superconductivity, the Cooper pairs momenta are lying along two lines in the Brillouin zone meaning simultaneous condensation of a continuum of Cooper pairs.


title: Towards Jointly Optimal Placement and Delivery: To Code or Not to Code in Wireless Caching Networks

abstract: Coded caching techniques have received significant attention lately due to their provable gains in reducing the cost of data delivery in wireless networks. These gains, however, have only been demonstrated under the assumption of a free placement phase. This unrealistic assumption poses a significant limitation, especially in cases where aggressive placement strategies can lead to a significant transmission cost that may even be higher than the corresponding cost of the delivery phase. In this paper, we relax this assumption and propose a general caching framework that captures the transmission cost of the two phases, and hence, results in minimizing the overall rate of the caching network. We model the dynamic nature of the network through a cost structure that allows for varying the network architecture and cost per transmission, across the placement and delivery phases. We start with the scenario where the individual users have no limit on the available caching memory and characterize the jointly optimal solution as a function of the different parameters in our cost structure. Then, we characterize the effect of memory constraints on the optimal solution in certain special cases. Interestingly, our results identify regions where the uncoded caching scheme outperforms its coded counterpart. Further, coded caching is shown to offer performance gains only when the network architecture during the placement phase is different from that during the delivery phase.


title: Contrastive Variational Autoencoder Enhances Salient Features

abstract: Variational autoencoders are powerful algorithms for identifying dominant latent structure in a single dataset. In many applications, however, we are interested in modeling latent structure and variation that are enriched in a target dataset compared to some background---e.g. enriched in patients compared to the general population. Contrastive learning is a principled framework to capture such enriched variation between the target and background, but state-of-the-art contrastive methods are limited to linear models. In this paper, we introduce the contrastive variational autoencoder (cVAE), which combines the benefits of contrastive learning with the power of deep generative models. The cVAE is designed to identify and enhance salient latent features. The cVAE is trained on two related but unpaired datasets, one of which has minimal contribution from the salient latent features. The cVAE explicitly models latent features that are shared between the datasets, as well as those that are enriched in one dataset relative to the other, which allows the algorithm to isolate and enhance the salient latent features. The algorithm is straightforward to implement, has a similar run-time to the standard VAE, and is robust to noise and dataset purity. We conduct experiments across diverse types of data, including gene expression and facial images, showing that the cVAE effectively uncovers latent structure that is salient in a particular analysis.


title: Dynamics of Hydrogen-oxygen-argon cellular detonations with a constant mean lateral strain rate

abstract: The present work revisits the problem of modelling the real gaseous detonation dynamics at the macro-scale by simple steady one-dimensional (1D) models. Experiments of detonations propagating in channels with exponentially expanding cross-sections were conducted in the H2/O2/Ar reactive system. Steady detonation waves were obtained at the macro-scale, with cellular structures characterized by reactive transverse waves. For all the mixtures studied, the dependence of the mean detonation speed was found to be in excellent agreement with first principles predictions of quasi-1D detonation dynamics with lateral strain rate predicted from detailed chemical kinetic models. This excellent agreement departs from the earlier experiments of Radulescu and Borzou (2018) in more unstable detonations. The excellent agreement is likely due to the much longer reaction zone lengths of argon diluted hydrogen-oxygen detonations at low pressures, as compared with the characteristic induction zone lengths. While the cellular instability modifies the detonation induction zone, the detonation dynamics at the macro-scale are arguably controlled by its hydrodynamic thickness. Near the limit, minor discrepancy is observed, with the experimental detonations typically continuing to propagate to slightly higher lateral strain rates and higher velocity deficits.


title: On a special kind of integral

abstract: In the world of mathematical analysis, many counterintuitive answers arise from the manipulation of seemingly unrelated concepts, ideas, or functions. For example, Euler showed that $e^{i\pi} + 1 = 0$, whereas Gauss proved that the area underneath $y = e^{-x^2}$ spanning the whole real axis is $ \sqrt{\pi} $. In this paper, we will determine the closed-form solution of the improper integral \[ I_n = \int_{0}^{\infty} \frac{\ln{x}}{x^n+1} dx, \ \forall n \in \mathbb{R} \text{, with}\ n > 1. \] Determining closed-form solutions of improper integrals have real implications not only in easing the solving of similar, yet more difficult integrals, but also in speeding up numerical approximations of the answer by making them more efficient. Following our calculations, we derived the formula \[ I_n = \int_{0}^{\infty} \frac{\ln{x}}{x^n+1} dx = -\frac{\pi^2}{n^2}\cot{\frac{\pi}{n}}\csc{\frac{\pi}{n}} = -\frac{d}{dn} \Bigg[ \Gamma\Big(1-\frac{1}{n}\Big) \Gamma\Big(\frac{1}{n}\Big) \Bigg]. \]


title: Progressively Growing Generative Adversarial Networks for High Resolution Semantic Segmentation of Satellite Images

abstract: Machine learning has proven to be useful in classification and segmentation of images. In this paper, we evaluate a training methodology for pixel-wise segmentation on high resolution satellite images using progressive growing of generative adversarial networks. We apply our model to segmenting building rooftops and compare these results to conventional methods for rooftop segmentation. We present our findings using the SpaceNet version 2 dataset. Progressive GAN training achieved a test accuracy of 93% compared to 89% for traditional GAN training.


title: Ultra-low damping in lift-off structured yttrium iron garnet thin films

abstract: We show that using maskless photolithography and the lift-off technique, patterned yttrium iron garnet thin films possessing ultra-low Gilbert damping can be accomplished. The films of 70 nm thickness were grown on (001)-oriented gadolinium gallium garnet by means of pulsed laser deposition, and they exhibit high crystalline quality, low surface roughness, and the effective magnetization of 127 emu/cm3. The Gilbert damping parameter is as low as 5x10-4. The obtained structures have well-defined sharp edges which along with good structural and magnetic film properties pave a path in the fabrication of high-quality magnonic circuits and oxide-based spintronic devices.


title: Quantifying the Loss of Information from Binning List-Mode Data

abstract: List-mode data is increasingly being uesd in SPECT and PET imaging, among other imaging modalities. However, there are still many imaging designs that effectively bin list-mode data before image reconstruction or other estimation tasks are performed. Intuitively, the binning operation should result in a loss of information. In this work we show that this is true for Fisher information and provide a computational method for quantifying the information loss. In the end we find that the information loss depends on three factors. The first factor is related to the smoothness of the mean data function for the list-mode data. The second factor is the actual object being imaged. Finally, the third factor is the binning scheme in relation to the other two factors.


title: Inequalities and Approximations for Fisher Information in the Presence of Nuisance Parameters

abstract: Many imaging systems are used to estimate a vector of parameters associated with the object being imaged. In many cases there are other parameters in the model for the imaging data that are not of interest for the task at hand. We refer to these as nuisance parameters and use them to form the components of the nuisance parameter vector. If we have a prior probability distribution function (PDF) for the nuisance parameter vector, then we may mariginalize over the nuisance parameters to produce a conditional PDF for the data that only depends on the parameters of interest. We will examine this approach to develop inequalities and approximations for the FIM when the data is affected by nuisance parameters.


title: Characterization of spin wave propagation in (111) YIG thin films with large anisotropy

abstract: We report on long-range spin wave (SW) propagation in nanometer-thick yttrium iron garnet (YIG) film with an ultralow Gilbert damping. The knowledge of a wavenumber value $|\vec{k}|$ is essential for designing SW devices. Although determining the wavenumber $|\vec{k}|$ in experiments like Brillouin light scattering spectroscopy is straightforward, quantifying the wavenumber in all-electrical experiments has not been widely commented upon so far. We analyze magnetostatic spin wave (SW) propagation in YIG films in order to determine the SW wavenumber $|\vec{k}|$ excited by the coplanar waveguide. We show that it is crucial to consider the influence of magnetic anisotropy fields present in YIG thin films for precise determination of SW wavenumber. With the proposed methods we find that experimentally derived values of $|\vec{k}|$ are in perfect agreement with that obtained from electromagnetic simulation only if anisotropy fields are included.


title: Super-Eddington accretion discs with advection and outflows around magnetized neutron stars

abstract: We present a model for a super-Eddington accretion disc around a magnetized neutron star taking into account advection of heat and the mass loss by the wind. The model is semi-analytical and predicts radial profiles of all basic physical characteristics of the accretion disc. The magnetospheric radius is found as an eigenvalue of the problem. When the inner disc is in radiation-pressure-dominated regime but does not reach its local Eddington limit, advection is mild, and the radius of the magnetosphere depends weakly on the accretion rate. Once approaching the local Eddington limit, the disc becomes advection-dominated, and the scaling for the magnetospheric radius with the mass accretion rate is similar to the classical Alfven relation. Allowing for the mass loss in a wind leads to an increase of the magnetospheric radius. Our model may be applied to a large variety of magnetized neutron stars accreting close to or above their Eddington limits: ultra-luminous X-ray pulsars, Be/X-ray binaries in outbursts, and other systems. In the context of our model we discuss the observational properties of NGC 5907~X-1, the brightest ultra-luminous pulsar known so far, and NGC 300~ULX-1 which is apparently a Be/X-ray binary experiencing a very bright super-Eddington outburst.


title: Salus: Fine-Grained GPU Sharing Primitives for Deep Learning Applications

abstract: GPU computing is becoming increasingly more popular with the proliferation of deep learning (DL) applications. However, unlike traditional resources such as CPU or the network, modern GPUs do not natively support fine-grained sharing primitives. Consequently, implementing common policies such as time sharing and preemption are expensive. Worse, when a DL application cannot completely use a GPU's resources, the GPU cannot be efficiently shared between multiple applications, leading to GPU underutilization. We present Salus to enable two GPU sharing primitives: fast job switching and memory sharing, in order to achieve fine-grained GPU sharing among multiple DL applications. Salus implements an efficient, consolidated execution service that exposes the GPU to different DL applications, and enforces fine-grained sharing by performing iteration scheduling and addressing associated memory management issues. We show that these primitives can then be used to implement flexible sharing policies such as fairness, prioritization, and packing for various use cases. Our integration of Salus with TensorFlow and evaluation on popular DL jobs show that Salus can improve the average completion time of DL training jobs by $3.19\times$, GPU utilization for hyper-parameter tuning by $2.38\times$, and GPU utilization of DL inference applications by $42\times$ over not sharing the GPU and $7\times$ over NVIDIA MPS with small overhead.


title: Ultrafast X-Ray Induced Changes of the Electronic and Magnetic Response of Solids Due to Valence Electron Redistribution

abstract: We report a novel mechanism, consisting of redistribution of valence electrons near the Fermi level, during interactions of intense femtosecond X-ray pulses with a Co/Pd multilayer. The changes in Co 3d valence shell occupation were directly revealed by fluence-dependent changes of the Co L$_3$ X-ray absorption and magnetic circular dichroism spectra near the excitation threshold. The valence shell redistribution arises from inelastic scattering of high energy Auger electrons and photoelectrons that lead to transient holes below and electrons above the Fermi level on the femtosecond time scale. The valence electron reshuffling effect scales with the energy deposited by X-rays and within 17 fs extends to valence states within 2 eV of the Fermi level. As a consequence the sample demagnetizes by more than twenty percent due to magnon generation.


title: Spatial segregation of dust grains in transition disks. SPHERE observations of 2MASS J16083070-3828268 and RXJ1852.3-3700

abstract: Context. The mechanisms governing the opening of cavities in transition disks are not fully understood. Several processes have been proposed but their occurrence rate is still unknown. Aims. We present spatially resolved observations of two transition disks and aim at constraining their vertical and radial structure using multiwavelength observations. Methods. We have obtained near-IR scattered light observations with VLT/SPHERE of the transition disks J1608 and J1852. We complement our datasets with ALMA observations and with unresolved photometric observations covering a wide range of wavelengths. We performed radiative transfer modeling to analyze the morphology of the disks and compare the results with a sample of 20 other transition disks observed with both SPHERE and ALMA. Results. The scattered light image of J1608 reveals a very inclined disk, with two bright lobes and a large cavity. J1852 shows an inner ring extending beyond the coronagraphic radius up to 15au, a gap and a second ring at 42au. Our radiative transfer model of J1608 indicates that the millimeter-sized grains are less extended vertically and radially than the micron-sized grains, indicating advanced settling and radial drift. We find good agreement with the observations of J1852 with a similar model, but due to the low inclination of the system, the model remains partly degenerate. The analysis of 22 transition disks shows that, in general, the cavities observed in scattered light are smaller than the ones detected at millimeter wavelengths. Conclusions. The analysis of a sample of transition disks indicates that the small grains can flow inward of the region where millimeter grains are trapped. While 15 out of the 22 cavities in our sample could be explained by a planet of less than 13 Jupiter masses, the others either require the presence of a more massive companion or of several low-mass planets.


title: Global labor flow network reveals the hierarchical organization and dynamics of geo-industrial clusters in the world economy

abstract: Groups of firms often achieve a competitive advantage through the formation of geo-industrial clusters. Although many exemplary clusters, such as Hollywood or Silicon Valley, have been frequently studied, systematic approaches to identify and analyze the hierarchical structure of the geo-industrial clusters at the global scale are rare. In this work, we use LinkedIn's employment histories of more than 500 million users over 25 years to construct a labor flow network of over 4 million firms across the world and apply a recursive network community detection algorithm to reveal the hierarchical structure of geo-industrial clusters. We show that the resulting geo-industrial clusters exhibit a stronger association between the influx of educated-workers and financial performance, compared to existing aggregation units. Furthermore, our additional analysis of the skill sets of educated-workers supplements the relationship between the labor flow of educated-workers and productivity growth. We argue that geo-industrial clusters defined by labor flow provide better insights into the growth and the decline of the economy than other common economic units.


title: A Spatial Filtering Approach to Biological Patterning

abstract: Interactions between neighboring cells are essential for generating or refining patterns in a number of biological systems. We propose a discrete filtering approach to predict how networks of cells modulate spatially varying input signals to produce more complicated or precise output signals. The interconnections between cells determine the set of spatial modes that are amplified or suppressed based on the coupling and internal dynamics of each cell, analogously to the way a traditional digital filter modifies the frequency components of a discrete signal. We apply the framework to two systems in developmental biology: the Notch-Delta interaction that shapes \textit{Drosophila} wing veins and the Sox9/Bmp/Wnt network responsible for digit formation in vertebrate limbs. The latter case study demonstrates that Turing-like patterns may occur even in the absence of instabilities. Results also indicate that developmental biological systems may be inherently robust to both correlated and uncorrelated noise sources. Our work shows that a spatial frequency-based interpretation simplifies the process of predicting patterning in living organisms when both environmental influences and intercellular interactions are present.


title: Gauge Equivariant Convolutional Networks and the Icosahedral CNN

abstract: The principle of equivariance to symmetry transformations enables a theoretically grounded approach to neural network architecture design. Equivariant networks have shown excellent performance and data efficiency on vision and medical imaging problems that exhibit symmetries. Here we show how this principle can be extended beyond global symmetries to local gauge transformations. This enables the development of a very general class of convolutional neural networks on manifolds that depend only on the intrinsic geometry, and which includes many popular methods from equivariant and geometric deep learning. We implement gauge equivariant CNNs for signals defined on the surface of the icosahedron, which provides a reasonable approximation of the sphere. By choosing to work with this very regular manifold, we are able to implement the gauge equivariant convolution using a single conv2d call, making it a highly scalable and practical alternative to Spherical CNNs. Using this method, we demonstrate substantial improvements over previous methods on the task of segmenting omnidirectional images and global climate patterns.


title: The impact of the crust equation of state on the analysis of GW170817

abstract: The detection of GW170817, the first neutron star-neutron star merger observed by Advanced LIGO and Virgo, and its following analyses represent the first contributions of gravitational wave (GW) data to understanding dense matter. Parameterizing the high density section of the equation of state (EOS) of both neutron stars through spectral decomposition, and imposing a lower limit on the maximum mass value, led to an estimate of the stars' radii of $R_1 = 11.9_{- 1.4}^{+ 1.4}$ km and $R_2 = 11.9_{- 1.4}^{+ 1.4}$ km. These values do not, however, take into account any uncertainty owed to the choice of the crust low-density EOS, which was fixed to reproduce the SLy EOS model. We here re-analyze GW170817 data and establish that different crust models do not strongly impact the mass or tidal deformability of a neutron star: it is impossible to distinguish between low-density models with GW analysis. However, the crust does have an effect on the inferred radius. We predict the systematic error due to this effect using neutron star structure equations, and compare the prediction to results from full parameter estimation runs. For GW170817, this systematic error affects the radius estimate by 0.3 km, approximately $3\%$ of the NS radii.


title: How to make the Earth orbit the Sun in 1614

abstract: In 1614 Johann Georg Locher, a student of the Jesuit astronomer Christoph Scheiner, proposed a physical mechanism to explain how the Earth could orbit the sun. An orbit, Locher said, is a perpetual fall. He proposed this despite the fact that he rejected the Copernican system, citing problems with falling bodies and the sizes of stars under that system. In 1651 and again in 1680, Jesuit writers Giovanni Battista Riccioli and Athanasius Kircher, respectively, considered and rejected outright Locher's idea of an orbit as a perpetual fall. Thus this important concept of an orbit was proposed, considered, and rejected well before Isaac Newton would use an entirely different physics to make the idea that an orbit is a perpetual fall the common way of envisioning and explaining orbits.


title: Security-Aware Synthesis Using Delayed-Action Games

abstract: Stochastic multiplayer games (SMGs) have gained attention in the field of strategy synthesis for multi-agent reactive systems. However, standard SMGs are limited to modeling systems where all agents have full knowledge of the state of the game. In this paper, we introduce delayed-action games (DAGs) formalism that simulates hidden-information games (HIGs) as SMGs, where hidden information is captured by delaying a player's actions. The elimination of private variables enables the usage of SMG off-the-shelf model checkers to implement HIGs. Furthermore, we demonstrate how a DAG can be decomposed into subgames that can be independently explored, utilizing parallel computation to reduce the model checking time, while alleviating the state space explosion problem that SMGs are notorious for. In addition, we propose a DAG-based framework for strategy synthesis and analysis. Finally, we demonstrate applicability of the DAG-based synthesis framework on a case study of a human-on-the-loop unmanned-aerial vehicle system under stealthy attacks, where the proposed framework is used to formally model, analyze and synthesize security-aware strategies for the system.


title: The number of ergodic measures for transitive subshifts under the regular bispecial condition

abstract: If $\mathcal{A}$ is a finite set (alphabet), the shift dynamical system consists of the space $\mathcal{A}^{\mathbb{N}}$ of sequences with entries in $\mathcal{A}$, along with the left shift operator $S$. Closed $S$-invariant subsets are called subshifts and arise naturally as encodings of other systems. In this paper, we study the number of ergodic measures for transitive subshifts under a condition ("regular bispecial condition") on the possible extensions of words in the associated language. Our main result shows that under this condition, the subshift can support at most $\frac{K+1}{2}$ ergodic measures, where $K$ is the limiting value of $p(n+1)-p(n)$, and $p$ is the complexity function of the language. As a consequence, we answer a question of Boshernitzan from `84, providing a combinatorial proof for the bound on the number of ergodic measures for interval exchange transformations.


title: Extreme Tensoring for Low-Memory Preconditioning

abstract: State-of-the-art models are now trained with billions of parameters, reaching hardware limits in terms of memory consumption. This has created a recent demand for memory-efficient optimizers. To this end, we investigate the limits and performance tradeoffs of memory-efficient adaptively preconditioned gradient methods. We propose extreme tensoring for high-dimensional stochastic optimization, showing that an optimizer needs very little memory to benefit from adaptive preconditioning. Our technique applies to arbitrary models (not necessarily with tensor-shaped parameters), and is accompanied by regret and convergence guarantees, which shed light on the tradeoffs between preconditioner quality and expressivity. On a large-scale NLP model, we reduce the optimizer memory overhead by three orders of magnitude, without degrading performance.


title: Testing the Reliability of Interpretable Neural Networks in Geoscience Using the Madden-Julian Oscillation

abstract: We test the reliability of two neural network interpretation techniques, backward optimization and layerwise relevance propagation, within geoscientific applications by applying them to a commonly studied geophysical phenomenon, the Madden-Julian Oscillation. The Madden-Julian Oscillation is a multi-scale pattern within the tropical atmosphere that has been extensively studied over the past decades, which makes it an ideal test case to ensure the interpretability methods can recover the current state of knowledge regarding its spatial structure. The neural networks can, indeed, reproduce the current state of knowledge and can also provide new insights into the seasonality of the Madden-Julian Oscillation and its relationships with atmospheric state variables. The neural network identifies the phase of the Madden-Julian Oscillation twice as accurately as linear regression, which means that nonlinearities used by the neural network are important to the structure of the Madden-Julian Oscillation. Interpretations of the neural network show that it accurately captures the spatial structures of the Madden-Julian Oscillation, suggest that the nonlinearities of the Madden-Julian Oscillation are manifested through the uniqueness of each event, and offer physically meaningful insights into its relationship with atmospheric state variables. We also use the interpretations to identify the seasonality of the Madden-Julian Oscillation, and find that the conventionally defined extended seasons should be shifted later by one month. More generally, this study suggests that neural networks can be reliably interpreted for geoscientific applications and may there


title: Learning Theory and Support Vector Machines - a primer

abstract: The main goal of statistical learning theory is to provide a fundamental framework for the problem of decision making and model construction based on sets of data. Here, we present a brief introduction to the fundamentals of statistical learning theory, in particular the difference between empirical and structural risk minimization, including one of its most prominent implementations, i.e. the Support Vector Machine.


title: Value constrained model-free continuous control

abstract: The naive application of Reinforcement Learning algorithms to continuous control problems -- such as locomotion and manipulation -- often results in policies which rely on high-amplitude, high-frequency control signals, known colloquially as bang-bang control. Although such solutions may indeed maximize task reward, they can be unsuitable for real world systems. Bang-bang control may lead to increased wear and tear or energy consumption, and tends to excite undesired second-order dynamics. To counteract this issue, multi-objective optimization can be used to simultaneously optimize both the reward and some auxiliary cost that discourages undesired (e.g. high-amplitude) control. In principle, such an approach can yield the sought after, smooth, control policies. It can, however, be hard to find the correct trade-off between cost and return that results in the desired behavior. In this paper we propose a new constraint-based reinforcement learning approach that ensures task success while minimizing one or more auxiliary costs (such as control effort). We employ Lagrangian relaxation to learn both (a) the parameters of a control policy that satisfies the desired constraints and (b) the Lagrangian multipliers for the optimization. Moreover, we demonstrate that we can satisfy constraints either in expectation or in a per-step fashion, and can even learn a single policy that is able to dynamically trade-off between return and cost. We demonstrate the efficacy of our approach using a number of continuous control benchmark tasks, a realistic, energy-optimized quadruped locomotion task, as well as a reaching task on a real robot arm.


title: Filaments in VIPERS: galaxy quenching in the infalling regions of groups

abstract: We study the quenching of galaxies in different environments and its evolution in the redshift range $0.43 \leq z \leq 0.89$. For this purpose, we identify galaxies inhabiting filaments, the isotropic infall region of groups, the field, and groups in the VIMOS Public Extragalactic Redshift Survey (VIPERS). We classify galaxies as quenched (passive), through their $NUV-r$ versus $r-K$ colours. We study the fraction of quenched galaxies ($F_{r}$) as a function of stellar mass and environment at two redshift intervals. Our results confirm that stellar mass is the dominant factor determining galaxy quenching over the full redshift range explored. We find compelling evidence of evolution in the quenching of intermediate mass galaxies $(9.3 \leq \log(M_{*}/M_{\odot}) \leq 10.5)$ for all environments. For this mass range, $F_{r}$ is largest for galaxies in groups and smallest for galaxies in the field, while galaxies in filaments and in the isotropic infall regions appear to have intermediate values with the exception of the high redshift bin, where the latter show similar fraction of quenched galaxies as in the field. Galaxies in filaments are systematically more quenched than their counterparts infalling from other directions, in agreement to similar results found at low redshift. The least massive galaxies in our samples do not show evidence of internal or environmental quenching.


title: Modular Inflation at Higher Level $N$

abstract: We introduce the framework of modular inflation with level structure, generalizing the level one theory considered previously to higher levels. We analyze the modular structure of CMB observables in this framework and show that the nontrivial geometry of the target space suffices to ensure the almost holomorphic modularity of the relevant parameters. We further introduce a concrete class of models based on hauptmodul functions that provide generators of the corresponding inflationary potentials at level $N>1$. The phenomenology of this class of models provides targets for ground-based CMB experiments in the immediate future. In the framework of our models we also discuss the status of the quantum gravity conjectures that have been formulated in the context of the swampland.


title: 2D Semiconductor Nonlinear Plasmonic Modulators

abstract: A plasmonic modulator is a device that controls the amplitude or phase of propagating plasmons. In a pure plasmonic modulator, the presence or absence of a pump plasmonic wave controls the amplitude of a probe plasmonic wave through a channel. This control has to be mediated by an interaction between disparate plasmonic waves, typically requiring the integration of a nonlinear material. In this work, we demonstrate the first 2D semiconductor nonlinear plasmonic modulator based on a WSe2 monolayer integrated on top of a lithographically defined metallic waveguide. We utilize the strong coupling between the surface plasmon polaritons, SPPs, and excitons in the WSe2 to give a 73 percent change in transmission through the device. We demonstrate control of the propagating SPPs using both optical and SPP pumps, realizing the first demonstration of a 2D semiconductor nonlinear plasmonic modulator, with a modulation depth of 4.1 percent, and an ultralow switching energy estimated to be 40 aJ.


title: Time-aware Test Case Execution Scheduling for Cyber-Physical Systems

abstract: Testing cyber-physical systems involves the execution of test cases on target-machines equipped with the latest release of a software control system. When testing industrial robots, it is common that the target machines need to share some common resources, e.g., costly hardware devices, and so there is a need to schedule test case execution on the target machines, accounting for these shared resources. With a large number of such tests executed on a regular basis, this scheduling becomes difficult to manage manually. In fact, with manual test execution planning and scheduling, some robots may remain unoccupied for long periods of time and some test cases may not be executed. This paper introduces TC-Sched, a time-aware method for automated test case execution scheduling. TC-Sched uses Constraint Programming to schedule tests to run on multiple machines constrained by the tests' access to shared resources, such as measurement or networking devices. The CP model is written in SICStus Prolog and uses the Cumulatives global constraint. Given a set of test cases, a set of machines, and a set of shared resources, TC-Sched produces an execution schedule where each test is executed once with minimal time between when a source code change is committed and the test results are reported to the developer. Experiments reveal that TC-Sched can schedule 500 test cases over 100 machines in less than 4 minutes for 99.5% of the instances. In addition, TC-Sched largely outperforms simpler methods based on a greedy algorithm and is suitable for deployment on industrial robot testing.


title: The effect of plasma-core induced self-guiding on phase matching of high-order harmonic generation in gases

abstract: In this work we numerically study a self-guiding process in which ionization plays a dominant role and analyze its effect on high-order harmonic generation (HHG) in gases. Although this type of self-guiding --- termed as plasma-core induced self-guiding in previous works --- limits the achievable cut-off by regulating the intensity of the laser beam, it provides favorable conditions for phase matching, which is indispensable for high-flux gas high-harmonic sources. To underline the role of self-guiding in efficient HHG, we investigate the time-dependent phase matching conditions in the guided beam and show how the spatio-temporally constant fundamental intensity contributes to the constructive build-up of the harmonic field in a broad photon-energy range up to the provided cut-off.


title: Crowdsourced PAC Learning under Classification Noise

abstract: In this paper, we analyze PAC learnability from labels produced by crowdsourcing. In our setting, unlabeled examples are drawn from a distribution and labels are crowdsourced from workers who operate under classification noise, each with their own noise parameter. We develop an end-to-end crowdsourced PAC learning algorithm that takes unlabeled data points as input and outputs a trained classifier. Our three-step algorithm incorporates majority voting, pure-exploration bandits, and noisy-PAC learning. We prove several guarantees on the number of tasks labeled by workers for PAC learning in this setting and show that our algorithm improves upon the baseline by reducing the total number of tasks given to workers. We demonstrate the robustness of our algorithm by exploring its application to additional realistic crowdsourcing settings.


title: Derivative-based global sensitivity analysis for models with high-dimensional inputs and functional outputs

abstract: We present a framework for derivative-based global sensitivity analysis (GSA) for models with high-dimensional input parameters and functional outputs. We combine ideas from derivative-based GSA, random field representation via Karhunen--Lo\`{e}ve expansions, and adjoint-based gradient computation to provide a scalable computational framework for computing the proposed derivative-based GSA measures. We illustrate the strategy for a nonlinear ODE model of cholera epidemics and for elliptic PDEs with application examples from geosciences and biotransport.


title: Some properties of coefficients of cyclotomic polynomials

abstract: This paper investigates coefficients of cyclotomic polynomials theoretically and experimentally. We prove the following result. {{\em If $n=p_1\ldots p_k$ where $p_i$ are odd primes and $p_1<p_2<\ldots<p_r<p_1+p_2<p_{r+1}<\ldots<p_t$ with $t\geq 3$ odd, then the numbers $-(r-2),-(r-3),\ldots, r-2, r-1$ are all coefficients of the cyclotomic polynomial $\Phi_{2n}$. Furthermore, if $1+p_r<p_1+p_2$ then $1-r$ is also a coefficient of $\Phi_{2n}$.} In the experimental part, in two instances we present computational evidence for asymptotic symmetry between distribution of positive and negative coefficients, and state the resulting conjectures.}


title: Collective Motion in the Interfacial and Interior Regions of Supported Polymer Films and its Relation to Relaxation

abstract: To understand the role of collective motion in the often large changes in interfacial molecular mobility observed in polymer films, we investigate the extent of collective motion in the interfacial regions of a thin supported polymer film and within the film interior by molecular dynamics simulation. Contrary to commonly stated expectations, we find that the extent of collective motion, as quantified by string-like molecular exchange motion, is similar in magnitude in polymer-air interfacial layer to that in the film interior, and both are generally distinct from the bulk material. This finding is consistent with Adam-Gibbs description of the segmental dynamics within mesoscopic film regions where the extent of collective motion is related to the configurational entropy of the film as whole rather than a locally defined extent of collective motion or configurational entropy.


title: Field-induced double dome and Bose-Einstein condensation in the crossing quantum spin chain system AgVOAsO4

abstract: We present inelastic neutron scattering data on the quantum paramagnet AgVOAsO4 that establish the system is a S=1/2 alternating spin chain compound and provide a direct measurement of the spin gap. We also present experimental evidence for two different types of field-induced magnetic order between mu_0H_c1 = 8.4T and mu_0H_c2 = 48.9T, which may be related to Bose-Einstein condensation (BEC) of triplons. Thermodynamic measurements in magnetic fields up to 60T and temperatures down to 0.1K reveal a H-T phase diagram consisting of a dome encapsulating two ordered phases with maximum ordering temperatures of 3.8K and 5.3K respectively. This complex phase diagram is not expected for a single-Q BEC system and therefore establishes AgVOAsO4 as a promising multi-Q BEC candidate capable of hosting exotic vortex phases.


title: On the selectivity of KcsA potassium channel: asymptotic analysis and computation

abstract: Potassium (K$^+$) channels regulate the flux of K$^+$ ions through cell membranes and plays significant roles in many physiological functions. This work studies the KcsA potassium channel, including the selectivity and current-voltage (IV) relations. A modified Poisson-Nernst-Planck system is employed, which include the size effect by Bikerman model and solvation energy by Born model. The selectivity of KcsA for various ions (K$^+$, Na$^+$, Cl$^-$, Ca$^{2+}$ and Ba$^{2+}$) is studied analytically, and the profiles of concentrations and electric potential are provided. The selectivity is mainly influenced by permanent negative charges in filter of channel and the ion sizes. K$^+$ is always selected compared with Na$^+$ (or Cl$^-$), as smaller ion size of Na$^+$ causes larger solvation energy. There is a transition for selectivity among K$^+$ and divalent ions (Ca$^{2+}$ and Ba$^{2+}$), when negative charge in filter exceeds a critical value determined by ion size. This explains why divalent ions can block the KcsA channel. The profiles and IV relations are studied by analytical, numerical and hybrid methods, and are cross-validated. The results show the selectivity of the channel and also the saturation of IV curve. A simple strategy is given to compute IV relations analytically, as first approximation. The numerical method deals with general structure or parameters, but the limitations and difficulties of pure numerical simulation are also pointed out. The hybrid method provides IV relations most effectively for comparison. The reason for saturation of IV relation is illustrated, and the IV curve shows agreement with the profile and scale of experimental results.


title: Optimal Budget-Feasible Mechanisms for Additive Valuations

abstract: In this paper, we show a tight approximation guarantee for budget-feasible mechanisms with an additive buyer. We propose a new simple randomized mechanism with approximation ratio of $2$, improving the previous best known result of $3$. Our bound is tight with respect to either the optimal offline benchmark, or its fractional relaxation. We also present a simple deterministic mechanism with the tight approximation guarantee of $3$ against the fractional optimum, improving the best known result of $(2+ \sqrt{2})$ for the weaker integral benchmark.


title: The Spin Budget of the Proton at NNLO and Beyond

abstract: We revisit the scale evolution of the quark and gluon spin contributions to the proton spin, $\frac{1}{2}\Delta \Sigma$ and $\Delta G$, using the three-loop results for the spin-dependent evolution kernels available in the literature. We argue that the evolution of the quark spin contribution may actually be extended to four-loop order, and that to all orders a single anomalous dimension governs the evolution of both $\Delta \Sigma$ and $\Delta G$. We present analytical solutions of the evolution equations for $\Delta \Sigma$ and $\Delta G$ and investigate their scale dependence both to large and down to lower "hadronic" scales. We find that the solutions remain perturbatively stable even to low scales, where they come closer to simple quark model expectations. We discuss a curious scenario for the proton spin, in which even the gluon spin contribution is essentially scale independent and has a finite asymptotic value as the scale becomes large. We finally also show that perturbative three-loop evolution leads to a larger spin contribution of strange anti-quarks than of strange quarks.


title: Scattered moonlight observations with X-Shooter: Implications for the aerosol properties at Cerro Paranal and the ESO sky background model

abstract: Estimating the sky background is critical for ground-based astronomical research. In the optical, scattered moonlight dominates the sky background, when the moon is above the horizon. The most uncertain component of a scattered moonlight model is the aerosol scattering. The current sky background model for Cerro Paranal uses an extrapolated aerosol extinction curve. With a set of X-Shooter sky observations, we have tested the current model as well as determined the aerosol extinction from the ultra-violet to near-infrared. To our knowledge, this is the first time that a scattered moonlight model has been used for this purpose. These observations were taken of blank sky, during three different lunar phases, and at six different angular distances from the moon for each night/lunar phase. Overall, the current model does reproduce the observations for average conditions decently well. Using a set of sky background models with varying aerosol distributions to compare with the observations, we found the most likely aerosol extinction curves, phase functions, and volume densities for the three nights of observations and compare them with the current model. While there were some degeneracies in the aerosol scattering properties, the extinction curves tend to flatten towards redder wavelengths and are overall less steep compared to the extrapolated curve used in the current model. Also, the current model had significantly less coarse particles compared to the favored volume densities from the X-Shooter data. For the three nights of sky observations, the aerosol size distributions differed, most likely reflecting the changes in atmospheric conditions and aerosol content. In short, the current sky background model is in fair agreement with the observations, and we have determined better aerosol extinction curves and phase functions for Cerro Paranal.


title: Whether Lyra's Manifold Itself is a Hidden Source of Dark Energy

abstract: In the course of investigation of some interesting cosmic string universes in the five dimensional Lyra manifold it is excitingly found that the geometry itself of Lyra manifold behaves as a new source of dark energy and this energy takes a form similar to that of quintessence in most of the cases, though in one case the dark energy comes out to be that of the cosmological constant type. The behaviour of the universes and their contribution to the process of evolution are examined. Further study of such type of universes will be very helpful in explaining the present accelerated expansion behaviour of the universe.


title: A Tunable Loss Function for Binary Classification

abstract: We present $\alpha$-loss, $\alpha \in [1,\infty]$, a tunable loss function for binary classification that bridges log-loss ($\alpha=1$) and $0$-$1$ loss ($\alpha = \infty$). We prove that $\alpha$-loss has an equivalent margin-based form and is classification-calibrated, two desirable properties for a good surrogate loss function for the ideal yet intractable $0$-$1$ loss. For logistic regression-based classification, we provide an upper bound on the difference between the empirical and expected risk at the empirical risk minimizers for $\alpha$-loss by exploiting its Lipschitzianity along with recent results on the landscape features of empirical risk functions. Finally, we show that $\alpha$-loss with $\alpha = 2$ performs better than log-loss on MNIST for logistic regression.


title: Regularity of extremal solutions of nonlocal elliptic systems

abstract: We examine regularity of the extremal solution of nonlinear nonlocal eigenvalue problem \begin{eqnarray} \left\{ \begin{array}{lcl} \hfill \mathcal L u &=& \lambda F(u,v) \qquad \text{in} \ \ \Omega, \\ \hfill \mathcal L v &=& \gamma G(u,v) \qquad \text{in} \ \ \Omega, \\ \hfill u,v &=&0 \qquad \qquad \text{on} \ \ \mathbb R^n\setminus\Omega , \end{array}\right. \end{eqnarray} with an integro-differential operator, including the fractional Laplacian, of the form \begin{equation*}\label{} \mathcal L(u (x))= \lim_{\epsilon\to 0} \int_{\mathbb R^n\setminus B_\epsilon(x) } [u(x) - u(z)] J(z-x) dz , \end{equation*} when $J$ is a nonnegative measurable even jump kernel. In particular, we consider jump kernels of the form of $J(y)=\frac{a(y/|y|)}{|y|^{n+2s}}$ where $s\in (0,1)$ and $a$ is any nonnegative even measurable function in $L^1(\mathbb {S}^{n-1})$ that satisfies ellipticity assumptions. We first establish stability inequalities for minimal solutions of the above system for a general nonlinearity and a general kernel. Then, we prove regularity of the extremal solution in dimensions $n < 10s$ and $ n<2s+\frac{4s}{p\mp 1}[p+\sqrt{p(p\mp1)}]$ for the Gelfand and Lane-Emden systems when $p>1$ (with positive and negative exponents), respectively. When $s\to 1$, these dimensions are optimal. However, for the case of $s\in(0,1)$ getting the optimal dimension remains as an open problem. Moreover, for general nonlinearities, we consider gradient systems and we establish regularity of the extremal solution in dimensions $n<4s$. As far as we know, this is the first regularity result on the extremal solution of nonlocal system of equations.


title: A Case for Superconducting Accelerators

abstract: As the scaling of conventional CMOS-based technologies slows down, there is growing interest in alternative technologies that can improve performance or energy-efficiency. Superconducting circuits based on Josephson Junction (JJ) is an emerging technology that can provide devices which can be switched with pico-second latencies and consuming two orders of magnitude lower switching energy compared to CMOS. While JJ-based circuits can provide high operating frequency and energy-efficiency, this technology faces three critical challenges: limited device density and lack of area-efficient technology for memory structures, reduced gate fanout compared to CMOS, and new failure modes of Flux-Traps that occurs due to the operating environment. The lack of dense memory technology restricts the use of superconducting technology in the near term to application domains that have high compute intensity but require negligible amount of memory. In this paper, we study the use of superconducting technology to build an accelerator for SHA-256 engines commonly used in Bitcoin mining applications. We show that merely porting existing CMOS-based accelerator to superconducting technology provides 10.6X improvement in energy efficiency. Redesigning the accelerator to suit the unique constraints of superconducting technology (such as low fanout) improves the energy efficiency to 12.2X. We also investigate solutions to make the accelerator tolerant of new fault modes and show how this fault-tolerant design can be leveraged to reduce the operating current, thereby increasing the overall energy-efficiency to 46X compared to CMOS. Our paper also develops a workflow for evaluating area, performance, and power for accelerators built in superconducting technology, and this workflow can help other researchers explore designs using this technology.


title: Positive Lyapunov Exponents and a Large Deviation Theorem for Continuum Anderson Models, Briefly

abstract: In this short note, we prove positivity of the Lyapunov exponent for 1D continuum Anderson models by leveraging some classical tools from inverse spectral theory. The argument is much simpler than the existing proof due to Damanik--Sims--Stolz, and it covers a wider variety of random models. Along the way we note that a Large Deviation Theorem holds uniformly on compacts.


title: Statistical Equilibrium of the Korteweg-de Vries and Benjamin-Ono Unidirectional Soliton Models

abstract: The Korteweg-de Vries and Benjamin-Ono nonlinear wave equations can describe solitary waves, all of which propagate in the same direction and which emerge from collisions with their shapes unchanged. There are technical challenges to giving a description of the thermal equilibrium of a dilute gas of such solitons. These challenges can be overcome by considering a stochastic source of new solitons located at $x=0$. The total intensity and the momentum distributions of the solitons emitted by such sources determine the equilibrium distribution of soliton momenta for $x>0$.


title: Constraints on bosonic dark matter from ultralow-field nuclear magnetic resonance

abstract: The nature of dark matter, the invisible substance making up over $80\%$ of the matter in the Universe, is one of the most fundamental mysteries of modern physics. Ultralight bosons such as axions, axion-like particles or dark photons could make up most of the dark matter. Couplings between such bosons and nuclear spins may enable their direct detection via nuclear magnetic resonance (NMR) spectroscopy: as nuclear spins move through the galactic dark-matter halo, they couple to dark-matter and behave as if they were in an oscillating magnetic field, generating a dark-matter-driven NMR signal. As part of the Cosmic Axion Spin Precession Experiment (CASPEr), an NMR-based dark-matter search, we use ultralow-field NMR to probe the axion-fermion "wind" coupling and dark-photon couplings to nuclear spins. No dark matter signal was detected above background, establishing new experimental bounds for dark-matter bosons with masses ranging from $1.8\times 10^{-16}$ to $7.8\times 10^{-14}$ eV.


title: Program Equivalence for Algebraic Effects via Modalities

abstract: This dissertation is concerned with the study of program equivalence and algebraic effects as they arise in the theory of programming languages. Algebraic effects represent impure behaviour in a functional programming language, such as input and output, exceptions, nondeterminism etc. all treated in a generic way. Program equivalence aims to identify which programs can be considered equal in some sense. This question has been studied for a long time but has only recently been extended to languages with algebraic effects, which are a newer development. Much work remains to be done in order to understand program equivalence in the presence of algebraic effects. In particular, there is no characterisation of contextual equivalence using a logic. We define a logic whose formulas express properties of higher-order programs with algebraic effects. We then investigate three notions of program equivalence for algebraic effects: logical equivalence induced by the aforementioned logic, applicative bisimilarity and contextual equivalence. For the programming language used in this dissertation, we prove that they all coincide. Therefore, the main novel contribution of the dissertation is defining the first logic for algebraic effects whose induced program equivalence coincides with contextual equivalence.


title: Weighted Tensor Completion for Time-Series Causal Inference

abstract: Marginal Structural Models (MSM) are the most popular models for causal inference from time-series observational data. However, they have two main drawbacks: (a) they do not capture subject heterogeneity, and (b) they only consider fixed time intervals and do not scale gracefully with longer intervals. In this work, we propose a new family of MSMs to address these two concerns. We model the potential outcomes as a three-dimensional tensor of low rank, where the three dimensions correspond to the agents, time periods and the set of possible histories. Unlike the traditional MSM, we allow the dimensions of the tensor to increase with the number of agents and time periods. We set up a weighted tensor completion problem as our estimation procedure, and show that the solution to this problem converges to the true model in an appropriate sense. Then we show how to solve the estimation problem, providing conditions under which we can approximately and efficiently solve the estimation problem. Finally we propose an algorithm based on projected gradient descent, which is easy to implement, and evaluate its performance on a simulated dataset.


title: Collisional dynamics of multiple dark solitons in a toroidal Bose-Einstein condensate: Quasiparticle picture

abstract: We study the collisional dynamics of multiple dark solitons in a Bose-Einstein condensate confined by a toroidal trap. We assume a tight enough confinement in the radial direction to prevent possible dissipative effects due to the presence of solitonic vortices. Analytical expressions for the initial order parameters with imprinted phases are utilized to generate different initial arrays of solitons, for which the time-dependent Gross-Pitaevskii equation is numerically solved. Given that the soliton velocity is conserved due to the lack of dissipation, we are able to apply a simple quasiparticle description of the soliton dynamics. In fact, the trajectory equations are written in terms of the velocities and the angular shifts produced at each collision, in analogy to the infinite one-dimensional system. To calculate the angular shifts, we directly extract them from the trajectories given by the Gross-Pitaevskii simulations and, on the other hand, we show that accurate values can be analytically obtained by adapting a formula valid for the infinite one-dimensional system that involves the healing length, which in our inhomogeneous system must be evaluated in terms of the sound velocity along the azimuthal direction. We further show that very good estimates of such a sound velocity can be directly determined by using the ground state density profile and the values of the imprinted phases. We discuss the possible implementation of the system here proposed using the current experimental techniques.


title: A consistent reduction of the two-layer shallow-water equations to an accurate one-layer spreading model

abstract: The gravity-driven spreading of one fluid in contact with another fluid is of key importance to a range of topics. To describe these phenomena, the two-layer shallow-water equations is commonly employed. When one layer is significantly deeper than the other, it is common to approximate the system with the much simpler one-layer shallow water equations. So far, it has been assumed that this approximation is invalid near shocks, and one has applied additional front conditions for the shock speed. In this paper, we prove mathematically that an effective one-layer model can be derived from the two-layer equations that correctly captures the behaviour of shocks and contact discontinuities without any additional closure relations. The proof yields a novel formulation of an effective one-layer shallow water model. The result shows that simplification to an effective one-layer model is well justified mathematically and can be made without additional knowledge of the shock behaviour. The shock speed in the proposed model is consistent with empirical models and identical to the front conditions that have been found theoretically by e.g. von K\'arm\'an and by Benjamin. This suggests that the breakdown of the shallow-water equations in the vicinity of shocks is less severe than previously thought. We further investigate the applicability of the shallow water framework to shocks by studying shocks in one-dimensional lock-exchange/lock-release. We derive expressions for the Froude number that are in good agreement with the widely employed expression by Benjamin. We then solve the equations numerically to illustrate how quickly the proposed model converges to solutions of the full two-layer shallow-water equations. We also compare numerical results using our model with results from dam break experiments. Predictions from the one-layer model are found to be in good agreement with experiments.


title: Seamless Rate Adaptation at Receiver Side for Indoor Visible Light Communications by Using Raptor Codes

abstract: Visible light communication (VLC), which uses LEDs for both illumination and communication purposes, has recently emerged as either an alternative or a support system for radio frequency (RF) communication. VLC channel inside a room is highly variable since it is pretty much depends on optical line of sight (LOS) and field of view (FOV) and it causes the received signal to vary dramatically. Therefore, to achieve spectral efficient communication, using rate adaptation methods are essential. In this study, we propose a rateless codes based rate adaptation method which adjusts transmission rate at receiver side and does not require channel state information (CSI) at the transmitter side. Spectral efficiency of proposed method is analyzed and compared with adaptive modulation and coding (AMC) rate adaptation method widely used in literature. For simulation work, we consider a realistic communication system model by using commercially available LEDs and also fulfilling illumination requirements. To overcome multi-path effects, optical orthogonal frequency-division multiplexing (OFDM) is used. Simulation results indicate that proposed rate adaptive method can achieve seamless rate adaptation and has higher spectral efficiency compared to AMC in presence of significant signal quality variations caused by mobility in a model room.


title: Confidence regions and minimax rates in outlier-robust estimation on the probability simplex

abstract: We consider the problem of estimating the mean of a distribution supported by the $k$-dimensional probability simplex in the setting where an $\varepsilon$ fraction of observations are subject to adversarial corruption. A simple particular example is the problem of estimating the distribution of a discrete random variable. Assuming that the discrete variable takes $k$ values, the unknown parameter $\boldsymbol \theta$ is a $k$-dimensional vector belonging to the probability simplex. We first describe various settings of contamination and discuss the relation between these settings. We then establish minimax rates when the quality of estimation is measured by the total-variation distance, the Hellinger distance, or the $\mathbb L^2$-distance between two probability measures. We also provide confidence regions for the unknown mean that shrink at the minimax rate. Our analysis reveals that the minimax rates associated to these three distances are all different, but they are all attained by the sample average. Furthermore, we show that the latter is adaptive to the possible sparsity of the unknown vector. Some numerical experiments illustrating our theoretical findings are reported.


title: Spin evolution and feedback of supermassive black holes in cosmological simulations

abstract: It is well established that the properties of supermassive black holes and their host galaxies are correlated through scaling relations. While hydrodynamical cosmological simulations have begun to account for the co-evolution of BHs and galaxies, they typically have neglected the BH spin, even though it may play an important role in modulating the growth and feedback of BHs. Here we introduce a new sub-grid model for the BH spin evolution in the moving-mesh code {\small AREPO} in order to improve the physical faithfulness of the BH modelling in galaxy formation simulations. We account for several different channels of spin evolution, in particular gas accretion through a Shakura-Sunyaev $\alpha$-disc, chaotic accretion, and BH mergers. For BH feedback, we extend the IllustrisTNG model, which considers two different BH feedback modes, a thermal quasar mode for high accretion states and a kinetic mode for low Eddington ratios, with a self-consistent accounting of spin-dependent radiative efficiencies and thus feedback strength. We find that BHs with mass $M_{\rm{bh}}\lesssim 10^{8}\, {\rm M}_\odot$ reach high spin values as they typically evolve in the coherent gas accretion regime. On the other hand, BHs with mass $M_{\rm{bh}}\gtrsim 10^{8}\, {\rm M}_\odot$ have lower spins as BH mergers become more frequent, and their accretion discs fragment due to self-gravity, inducing chaotic accretion. We also explore the hypothesis that the transition between the quasar and kinetic feedback modes is mediated by the accretion mode of the BH disc itself, i.e.~the kinetic feedback mode is activated when the disc enters the self-gravity regime. We find excellent agreement between the galaxy and BH populations for this approach and the fiducial TNG model with no spin evolution. Furthermore, our new approach alleviates a tension in the galaxy morphology -- colour relation of the original TNG model.


title: Survivability of planetary systems in young and dense star clusters

abstract: We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced. Four, five, or size planets per star were introduced in orbit around the 500 solar-like stars with a maximum orbital separation of 400au. Our study focuses on the production of free-floating planets. A total of 357 become unbound from a total of 2522 planets in the initial conditions of the simulation. Of these, 281 leave the cluster within the crossing timescale of the star cluster; the others remain bound to the cluster as free-floating intra-cluster planets. Five of these free-floating intra-cluster planets are captured at a later time by another star. The two main mechanisms by which planets are lost from their host star, ejection upon a strong encounter with another star or internal planetary scattering, drive the evaporation independent of planet mass of orbital separation at birth. The effect of small perturbations due to slow changes in the cluster potential are important for the evolution of planetary systems. In addition, the probability of a star to lose a planet is independent of the planet mass and independent of its initial orbital separation. As a consequence, the mass distribution of free-floating planets indistinguishable from the mass distribution of planets bound to their host star.


title: Modified second-order generalized integrators with modified frequency locked loop for fast harmonics estimation of distorted single-phase signals (LONG VERSION)

abstract: This paper proposes modified Second-Order Generalized Integrators (mSOGIs) for a fast estimation of all harmonic components of arbitrarily distorted single-phase signals such as voltages or currents in power systems. The estimation is based on the internal model principle leading to an overall observer system consisting of parallelized mSOGIs. The observer is tuned by pole placement. For a constant fundamental frequency, the observer is capable of estimating all harmonic components with prescribed settling time by choosing the observer poles appropriately. For time-varying fundamental frequencies, the harmonic estimation is combined with a modified Frequency Locked Loop (mFLL) with gain normalization, sign-correct anti-windup and rate limitation. The estimation performances of the proposed parallelized mSOGIs with and without mFLL are illustrated and validated by measurement results. The results are compared to standard approaches such as parallelized standard SOGIs (sSOGIs) and adaptive notch filters (ANFs).


title: Holographic entanglement contour, bit threads, and the entanglement tsunami

abstract: We study the entanglement contour, a quasi-local measure of entanglement, and propose a generic formula for the contour in 1+1d quantum systems. We use this formalism to investigate the real space entanglement structure of various static CFTs as well as local and global quantum quenches. The global quench elucidates the spatial distribution of entanglement entropy in strongly interacting CFTs and clarifies the interpretation of the entanglement tsunami picture. The entanglement tsunami effectively characterizes the non-local growth of entanglement entropy while the contour characterizes the local propagation of entanglement. We generalize the formula for the entanglement contour to arbitrary dimensions and entangling surface geometries using bit threads, and are able to realize a holographic contour for logarithmic negativity and the entanglement of purification by restricting the bulk spacetime to the entanglement wedge. Furthermore, we explore the connections between the entanglement contour, bit threads, and entanglement density in kinematic space.


title: Electron reconstruction and identification in the ATLAS experiment using the 2015 and 2016 LHC proton-proton collision data at $\sqrt{s}$ = 13 TeV

abstract: Algorithms used for the reconstruction and identification of electrons in the central region of the ATLAS detector at the Large Hadron Collider (LHC) are presented in this paper; these algorithms are used in ATLAS physics analyses that involve electrons in the final state and which are based on the 2015 and 2016 proton-proton collision data produced by the LHC at $\sqrt{s}$ = 13 TeV. The performance of the electron reconstruction, identification, isolation, and charge identification algorithms is evaluated in data and in simulated samples using electrons from $Z\rightarrow ee$ and $J/\psi\rightarrow ee$ decays. Typical examples of combinations of electron reconstruction, identification, and isolation operating points used in ATLAS physics analyses are shown.


title: High repetition rate (>= kHz) targets and optics from liquid microjets for the study and application of high intensity laser-plasma interactions

abstract: High intensity laser-plasma interactions produce a wide array of energetic particles and beams with promising applications. Unfortunately, high repetition rate and high average power requirements for many applications are not satisfied by the lasers, optics, targets, and diagnostics currently employed. Here, we address the need for high repetition rate targets and optics through the use of liquids. A novel nozzle assembly is used to generate high-velocity, laminar-flowing liquid microjets which are compatible with a low-vacuum environment, generate little to no debris, and exhibit precise positional and dimensional tolerances. Jets, droplets, submicron thick sheets, and other configurations are characterized with pump-probe shadowgraphy to evaluate their use as targets. To demonstrate a high repetition rate, consumable liquid optical element, we present a plasma mirror created by a submicron thick liquid sheet. This plasma mirror provides etalon-like anti-reflection properties in the low-field of 0.1% and high reflectivity as a plasma, 69%, at a repetition rate of 1 kHz. Practical considerations of fluid compatibility, in-vacuum operation, and estimates of maximum repetition rate, in excess of 10 kHz, are addressed. The targets and optics presented here enable the use of relativistically intense lasers at high average power and make possible many long proposed applications.


title: Enhancing entanglement detection of quantum optical frequency combs via stimulated emission

abstract: We investigate the performance of a certain nonclassicality identifier, expressed via integrated second-order intensity moments of optical fields, in revealing bipartite entanglement of quantum-optical frequency combs (QOFCs), which are generated in both spontaneous and stimulated parametric down-conversion processes. We show that, by utilizing that nonclassicality identifier, one can well identify the entanglement of the QOFC directly from the experimentally measured intensity moments without invoking any state reconstruction techniques or homodyne detection. Moreover, we demonstrate that the stimulated generation of the QOFC improves the entanglement detection of these fields with the nonclassicality identifier. Additionally, we show that the nonclassicality identifier can be expressed in a factorized form of detectors quantum efficiencies and the number of modes, if the QOFC consists of many copies of the same two-mode twin beam. As an example, we apply the nonclassicality identifier to two specific types of QOFC, where: (i) the QOFC consists of many independent two-mode twin beams with non-overlapped spatial frequency modes, and (ii) the QOFC contains entangled spatial frequency modes which are completely overlapped, i.e., each mode is entangled with all the remaining modes in the system. We show that, in both cases, the nonclassicality identifier can reveal bipartite entanglement of the QOFC including noise, and that it becomes even more sensitive for the stimulated processes.


title: Geometry of B\"acklund Transformations I: Generality

abstract: Using Cartan's Method of Equivalence, we prove an upper bound for the generality of generic rank-1 B\"acklund transformations relating two hyperbolic Monge-Amp\`ere systems. In cases when the B\"acklund transformation admits a symmetry group whose orbits have codimension 1, 2, or 3, we obtain classification results and new examples of auto-B\"acklund transformations.


title: Cost Analysis of Nondeterministic Probabilistic Programs

abstract: We consider the problem of expected cost analysis over nondeterministic probabilistic programs, which aims at automated methods for analyzing the resource-usage of such programs. Previous approaches for this problem could only handle nonnegative bounded costs. However, in many scenarios, such as queuing networks or analysis of cryptocurrency protocols, both positive and negative costs are necessary and the costs are unbounded as well. In this work, we present a sound and efficient approach to obtain polynomial bounds on the expected accumulated cost of nondeterministic probabilistic programs. Our approach can handle (a) general positive and negative costs with bounded updates in variables; and (b) nonnegative costs with general updates to variables. We show that several natural examples which could not be handled by previous approaches are captured in our framework. Moreover, our approach leads to an efficient polynomial-time algorithm, while no previous approach for cost analysis of probabilistic programs could guarantee polynomial runtime. Finally, we show the effectiveness of our approach by presenting experimental results on a variety of programs, motivated by real-world applications, for which we efficiently synthesize tight resource-usage bounds.


title: Magnetic Properties and Flow Angle of the Inverse Evershed Flow at Its Downflow Points

abstract: We determined the direction and strength of the photospheric and lower chromospheric magnetic field in the umbra and penumbra of a sunspot from inversions of spectropolarimetric observations of photospheric lines at 617\,nm and 1565\,nm, and the chromospheric \ion{Ca}{ii} IR line at 854\,nm, respectively. We compare the magnetic field vector with the direction of 75 flow channels that harbor the chromospheric inverse Evershed effect (IEF) near their downflow points (DFPs) in the sunspot's penumbra. The azimuth and inclination of the IEF channels to the line of sight (LOS) were derived from spatial maps of the LOS velocity and line-core intensity of the \ion{Ca}{ii} IR line and a thermal inversion of the \ion{Ca}{ii} IR spectra to obtain temperature cubes. We find that the flow direction of the IEF near the DFPs is aligned with the photospheric magnetic field to within about $\pm$\,15\,deg. The IEF flow fibrils make an angle of 30--90\,deg to the local vertical with an average value of about 65\,deg. The average field strength at the DFPs is about 1.3\,kG. Our findings suggest that the IEF in the lower chromosphere is a field-aligned siphon flow, where the larger field strength at the inner footpoints together with the lower temperature in the penumbra causes the necessary gas pressure difference relative to the outer footpoints in the hotter quiet Sun with lower magnetic field strength. The IEF connects to magnetic field lines that are not horizontal like for the regular photospheric Evershed flow, but which continue upwards into the chromosphere indicating an "uncombed" penumbral structure.


title: Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation

abstract: Geometric integrators of the Schr\"{o}dinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to systems whose Hamiltonian is separable into a kinetic and potential terms. Here, we describe several implicit geometric integrators applicable to both separable and non-separable Hamiltonians, and, in particular, to the nonadiabatic molecular Hamiltonian in the adiabatic representation. These integrators combine the dynamic Fourier method with recursive symmetric composition of the trapezoidal rule or implicit midpoint method, which results in an arbitrary order of accuracy in the time step. Moreover, these integrators are exactly unitary, symplectic, symmetric, time-reversible, and stable, and, in contrast to the split-operator algorithm, conserve energy exactly, regardless of the accuracy of the solution. The order of convergence and conservation of geometric properties are proven analytically and demonstrated numerically on a two-surface NaI model in the adiabatic representation. Although each step of the higher order integrators is more costly, these algorithms become the most efficient ones if higher accuracy is desired; a thousand-fold speedup compared to the second-order trapezoidal rule (the Crank-Nicolson method) was observed for wavefunction convergence error of $10^{-10}$. In a companion paper [J. Roulet, S. Choi, and J. Van\'{\i}\v{c}ek (2019)], we discuss analogous, arbitrary-order compositions of the split-operator algorithm and apply both types of geometric integrators to a higher-dimensional system in the diabatic representation.


title: High-order mass-lumped schemes for nonlinear degenerate elliptic equations

abstract: We present and analyse a numerical framework for the approximation of nonlinear degenerate elliptic equations of the Stefan or porous medium types. This framework is based on piecewise constant approximations for the functions, which we show are essentially necessary to obtain convergence and error estimates. Convergence is established without regularity assumption on the solution. A detailed analysis is then performed to understand the design properties that enable a scheme, despite these piecewise constant approximations and the degeneracy of the model, to satisfy high-order error estimates if the solution is piecewise smooth. Numerical tests, based on continuous and discontinuous approximation methods, are provided on a variety of 1D and 2D problems, showing the influence on the convergence rate of the nature of the degeneracy and of the design choices.


title: A Privacy-Preserving Traffic Monitoring Scheme via Vehicular Crowdsourcing

abstract: The explosive growth of vehicle amount has given rise to a series of traffic problems, such as traffic congestion, road safety, and fuel waste. Collecting vehicles' speed information is an effective way to monitor the traffic condition and avoid vehicles being congested, which however may bring threats to vehicles' location and trajectory privacy. Motivated by the fact that traffic monitoring does not need to know each individual vehicle's speed and the average speed would be sufficient, we propose a privacy-preserving traffic monitoring (PPTM) scheme to aggregate vehicles' speeds at different locations. In PPTM, the roadside unit (RSU) collects vehicles' speed information at multiple road segments, and further cooperates with a service provider to calculate the average speed information for every road segment. To preserve vehicles' privacy, both homomorphic Paillier cryptosystem and super-increasing sequence are adopted. A comprehensive security analysis indicates that the proposed PPTM can preserve vehicles' identities, speeds, locations, and trajectories privacy from being disclosed. In addition, extensive simulations are conducted to validate the effectiveness and efficiency of the proposed PPTM scheme.


title: Learning Generative Models of Structured Signals from Their Superposition Using GANs with Application to Denoising and Demixing

abstract: Recently, Generative Adversarial Networks (GANs) have emerged as a popular alternative for modeling complex high dimensional distributions. Most of the existing works implicitly assume that the clean samples from the target distribution are easily available. However, in many applications, this assumption is violated. In this paper, we consider the observation setting when the samples from target distribution are given by the superposition of two structured components and leverage GANs for learning the structure of the components. We propose two novel frameworks: denoising-GAN and demixing-GAN. The denoising-GAN assumes access to clean samples from the second component and try to learn the other distribution, whereas demixing-GAN learns the distribution of the components at the same time. Through extensive numerical experiments, we demonstrate that proposed frameworks can generate clean samples from unknown distributions, and provide competitive performance in tasks such as denoising, demixing, and compressive sensing.


title: Effective field analysis using the full angular spin-orbit torque magnetometry dependence

abstract: Spin-orbit torques promise ultra-efficient magnetization switching used for advanced devices based on emergent quasi-particles such as domain walls and skyrmions. Recently, the spin structure dynamics, materials and systems with tailored spin-orbit torques are being developed. A method, which allows one to detect the acting torques in a given system as a function of the magnetization direction is the torque-magnetometry method based on a higher harmonics analysis of the anomalous Hall-effect. Here we show that the effective fields acting on magnetic domain walls that govern the efficiency of their dynamics require a sophisticated analysis taking into account the full angular dependence of the torques. Using a 1-D model we compared the spin orbit torque efficiencies by depinning measurements and spin torque magnetometry. We show that the effective fields can be accurately determined and we find good agreement. Thus our method allows us now to rapidly screen materials and predict the resulting quasi-particle dynamics.


title: Nonsteady-state diffusion in two-dimensional periodic channels

abstract: The dynamics of a freely diffusing particle in a two-dimensional channel with cross sectional area $A(x)$, can be effectively described by a one-dimensional diffusion equation under the action of a potential of mean force $U(x)=-k_BT\ln [A(x)]$ (where $k_BT$ is the thermal energy) in a system with a spatially-dependent diffusion coefficient $D(x)$. Several attempts to derive expressions relating $D(x)$ to $A(x)$ and its derivatives have been made, which were based on considering stationary flows in periodic channels. Here, we take an alternative approach and consider non-steady state single particle diffusion in an open periodic channel. The approach allows us to express $D(x)$ as a series of terms of increasing powers of $\epsilon$ - a parameter associated with the aspect ratio of the channel. When the expansion is truncated at the leading term, we recover the expression suggested by Zwanzig [J. Phys. Chem. {\bf 96}, 3926 (1992)] for $D(x)$. Furthermore, comparison of the first few terms in our expansion for $D(x)$ with the one proposed by Kalinay and Percus [Phys. Rev. E {\bf 74}, 041203 (2006)] shows that they are consistent with each other. In the limit of long wavelength channels ($\epsilon\ll 1$), the expansion converges rapidly and the leading approximation provides a very accurate description of the two-dimensional dynamics. For short wavelength channels, the expansion does not converge and the validity of the effective one-dimensional description is questionable.