title: The impact of the metallicity and star formation rate on the time-dependent galaxy-wide stellar initial mass function

abstract: The stellar initial mass function (IMF) is commonly assumed to be an invariant probability density distribution function of initial stellar masses being represented by the canonical IMF. As a consequence the galaxy-wide IMF (gwIMF), defined as the sum of the IMFs of all star forming regions, should also be invariant. Recent observational and theoretical results challenge the hypothesis that the gwIMF is invariant. In order to study the possible reasons for this variation we use the IMF determined in resolved star clusters and apply the IGIMF-theory to calculate a grid of gwIMF models for metallicities, -3<[Fe/H]<1, and galaxy-wide star formation rates, $10^{-5}$<SFR<$10^{5}\,\mathrm{M_{\odot}/yr}$. For a galaxy with metallicy [Fe/H]$<0$ and SFR$\,> 1\,M_\odot$/yr, which is a common condition in the early Universe, we find that the gwIMF is top-heavy (more massive stars), when compared to the canonical IMF. For a SFR $< 1\,\mathrm{M_{\odot}/yr}$ the gwIMF becomes top-light regardless of the metallicity. For metallicities $\mathrm{[Fe/H]} > 0$ the gwIMF can become bottom-heavy regardless of the SFR. The IGIMF models predict that massive elliptical galaxies should have formed with a gwIMF that is top-heavy within the first few hundred Myr of the galaxy's life and that it evolves into a bottom-heavy gwIMF in the metal-enriched galactic center. We study the SFR$-$H$\alpha$ relation, its dependency on metallicity and the SFR, the correction factors to the Kennicutt SFR$_{\rm K}-$H$\alpha$ relation, and provide new fitting functions Late-type dwarf galaxies show significantly higher SFRs with respect to Kennicutt SFRs, while star forming massive galaxies have significantly lower SFRs than hitherto thought. This has implications for the gas-consumption time scales and for the main sequence of galaxies. The Leo P and ultra-faint dwarf galaxies are discussed explicitly. [abridged]


title: Magic-angle semimetals

abstract: Breakthroughs in two-dimensional van der Waals heterostructures have revealed that twisting creates a moir\'e pattern that quenches the kinetic energy of electrons, allowing for exotic many-body states. We show that cold-atomic, trapped ion, and metamaterial systems can emulate the effects of a twist in many models from one to three dimensions. Further, we demonstrate at larger angles (and argue at smaller angles) that by considering incommensurate effects, the magic-angle effect becomes a single-particle quantum phase transition (including in a model for twisted bilayer graphene in the chiral limit). We call these models "magic-angle semimetals." Each contains nodes in the band structure and an incommensurate modulation. At magic-angle criticality, we report a nonanalytic density of states, flat bands, multifractal wave functions that Anderson delocalize in momentum space, and an essentially divergent effective interaction scale. As a particular example, we discuss how to observe this effect in an ultracold Fermi gas.


title: Merging black hole binaries with the SEVN code

abstract: Studying the formation and evolution of black hole binaries (BHBs) is essential for the interpretation of current and forthcoming gravitational wave (GW) detections. We investigate the statistics of BHBs that form from isolated binaries, by means of a new version of the SEVN population-synthesis code. SEVN integrates stellar evolution by interpolation over a grid of stellar evolution tracks. We upgraded SEVN to include binary stellar evolution processes and we used it to evolve a sample of $1.5\times{}10^8$ binary systems, with metallicity in the range $\left[10^{-4};4\times 10^{-2}\right]$. From our simulations, we find that the mass distribution of black holes (BHs) in double compact-object binaries is remarkably similar to the one obtained considering only single stellar evolution. The maximum BH mass we obtain is $\sim 30$, $45$ and $55\, \mathrm{M}_\odot$ at metallicity $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. A few massive single BHs may also form ($\lesssim 0.1\%$ of the total number of BHs), with mass up to $\sim 65$, $90$ and $145\, \mathrm{M}_\odot$ at $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. These BHs fall in the mass gap predicted from pair-instability supernovae. We also show that the most massive BHBs are unlikely to merge within a Hubble time. In our simulations, merging BHs like GW151226 and GW170608, form at all metallicities, the high-mass systems (like GW150914, GW170814 and GW170104) originate from metal poor ($Z\lesssim{}6\times 10^{-3}$) progenitors, whereas GW170729-like systems are hard to form, even at $Z = 10^{-4}$. The BHB merger rate in the local Universe obtained from our simulations is $\sim 90 \mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$, consistent with the rate inferred from LIGO-Virgo data.


title: Invisible widths of heavy mesons

abstract: We revisit calculations of invisible widths of heavy mesons in the standard model, which serve as benchmarks for the studies of production of light, long-lived neutral particles in heavy meson decays. We challenge the common assumption that in the standard model these widths are dominated by meson decays into a two-neutrino final state and prove that they are dominated by decays into four-neutrino final states. We show that current estimates of the invisible widths of heavy mesons in the standard model underestimate the effect by orders of magnitude. We examine currently available experimental data on invisible widths and place constraints on the properties of dark photons. We also comment on the invisible widths of the kaons.


title: Star Cluster Formation from Turbulent Clumps. II. Gradual Star Cluster Formation

abstract: We investigate the dynamical evolution of star clusters during their formation, assuming that they are born from a turbulent starless clump of a given mass that is embedded within a parent self-gravitating molecular cloud characterized by a particular mass surface density. In contrast to the standard practice of most $N$-body studies, we do not assume that all stars are formed at once. Rather, we explore the effects of different star formation rates on the global structure and evolution of young embedded star clusters, also considering various primordial binary fractions and mass segregation levels. Our fiducial clumps studied in this paper have initial masses of $M_{\rm cl} = 3000\,M_\odot$, are embedded in ambient cloud environments of $\Sigma_{\rm cloud} = 0.1$ and 1 g cm$^{-2}$, and gradually form stars with an overall efficiency of 50% until the gas is exhausted. We investigate star formation efficiencies per free-fall time in the range $\epsilon_{\rm ff}=0.01$ to 1, and also compare to the instantaneous case ($\epsilon_{\rm ff}=\infty$) of Paper I. We show that most of the interesting dynamical processes that determine the future of the cluster, happen during the early formation phase. In particular, the ejected stellar population is sensitive to the duration of star cluster formation: for example, clusters with longer formation times produce more runaway stars, since these clusters remain in a dense state for longer, thus favouring occurrence of dynamical ejections. We also show that the presence of radial age gradients in star clusters depends sensitively on the star formation efficiency per free fall time, with observed values being matched best by our slowest forming clusters with $\epsilon_{\rm ff}\lesssim0.03$.


title: Large-Scale Poloidal Magnetic Field Dynamo Leads to Powerful Jets in GRMHD Simulations of Black Hole Accretion with Toroidal Field

abstract: Accreting black holes (BHs) launch relativistic collimated jets, across many decades in luminosity and mass, suggesting the jet launching mechanism is universal, robust and scale-free. Theoretical models and general relativistic magnetohydrodynamic (GRMHD) simulations indicate that the key jet-making ingredient is large-scale poloidal magnetic flux. However, its origin is uncertain, and it is unknown if it can be generated in situ or dragged inward from the ambient medium. Here, we use the GPU-accelerated GRMHD code H-AMR to study global 3D BH accretion at unusually high resolutions more typical of local shearing box simulations. We demonstrate that turbulence in a radially-extended accretion disc can generate large-scale poloidal magnetic flux in situ, even when starting from a purely toroidal magnetic field. The flux accumulates around the BH till it becomes dynamically-important, leads to a magnetically arrested disc (MAD), and launches relativistic jets that are more powerful than the accretion flow. The jet power exceeds that of previous GRMHD toroidal field simulations by a factor of 10,000. The jets do not show significant kink or pinch instabilities, accelerate to $\gamma \sim 10$ over 3 decades in distance, and follow a collimation profile similar to the observed M87 jet.


title: Pyaneti: a fast and powerful software suite for multi-planet radial velocity and transit fitting

abstract: Transiting exoplanet parameter estimation from time-series photometry and Doppler spectroscopy is fundamental to study planets' internal structures and compositions. Here we present the code pyaneti, a powerful and user-friendly software suite to perform multi-planet radial velocity and transit data fitting. The code uses a Bayesian approach combined with an MCMC sampling to estimate the parameters of planetary systems. We combine the numerical efficiency of FORTRAN, the versatility of PYTHON, and the parallelization of OpenMP to make pyaneti a fast and easy to use code. The package is freely available at https://github.com/oscaribv/pyaneti.


title: Constraining Sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. III. Candidates from Continued Radial Velocity Tests

abstract: Quasars whose broad emission lines show temporal, bulk radial velocity (RV) shifts have been proposed as candidate sub-parsec (sub-pc), binary supermassive black holes (BSBHs). We identified a sample of 16 BSBH candidates based on two-epoch spectroscopy among 52 quasars with significant RV shifts over a few rest-frame years. The candidates showed consistent velocity shifts independently measured from two broad lines (H$\beta$ and H$\alpha$ or Mg${\rm \,II}$) without significant changes in the broad-line profiles. Here in the third paper of the series, we present further third- and fourth-epoch spectroscopy for 12 of the 16 candidates for continued RV tests, spanning $\sim$5--15 yr in the quasars' rest frames. Cross-correlation analysis of the broad H$\beta$ calibrated against [O${\rm\,III}]\,\lambda 5007$ suggests that 5 of the 12 quasars remain valid as BSBH candidates. They show broad H$\beta$ RV curves that are consistent with binary orbital motion without significant changes in the broad line profiles. Their broad H$\alpha$ (or Mg${\rm \,II}$) lines display RV shifts that are either consistent with or smaller than those seen in broad H$\beta$. The RV shifts can be explained by a $\sim$0.05--0.1 pc BSBH with an orbital period of $\sim$40--130 yr, assuming a mass ratio of 0.5--2 and a circular orbit. However, the parameters are not well constrained given the few epochs that sample only a small portion of the hypothesized binary orbital cycle. The apparent occurrence rate of sub-pc BSBHs is $\lesssim$13$\pm$5% among all SDSS quasars, with no significant difference in the subsets with and without single-epoch broad line velocity offsets. Dedicated long-term spectroscopic monitoring is still needed to further confirm or reject these BSBH candidates.


title: Radiation of scalar modes and the classical double copy

abstract: The double copy procedure relates gauge and gravity theories through color-kinematics replacements and holds for both scattering amplitudes and in classical contexts. Moreover, it has been shown that there is a web of theories whose scattering amplitudes are related through operations that exchange color and kinematic factors. In this paper, we generalize and extend this procedure by showing that the classical perturbative double copy of pions corresponds to special Galileons. We consider point-particles coupled to the relevant scalar fields, and find the leading and next to leading order radiation amplitudes. By considering couplings motivated by those that would arise from extracting the longitudinal modes of the gauge and gravity theories, we are able to map the non-linear sigma model radiation to that of the special Galileon. We also construct the single copy by mapping the bi-adjoint scalar radiation to the non-linear sigma model radiation through generalized color-kinematics replacements.


title: Black-hole evaporation, cosmic censorship, and a quantum lower bound on the Bekenstein-Hawking temperature

abstract: The semi-classical Hawking evaporation process of Reissner-Nordstr\"om black holes is analyzed. It is shown that this quantum mechanism may turn a near-extremal black-hole spacetime with $T_{\text{BH}}\lesssim \hbar^2/G^2M^3$ into an horizonless naked singularity, thus violating the Penrose cosmic censorship conjecture. It is therefore conjectured that, within the framework of a self-consistent quantum theory of gravity, the Bekenstein-Hawking temperature should be bounded from below by the simple relation $T_{\text{BH}}\gtrsim \hbar^2/G^2M^3$.


title: Discrepancies between extinction events and boundary equilibria in reaction networks

abstract: Reaction networks are mathematical models of interacting chemical species that are primarily used in biochemistry. There are two modeling regimes that are typically used, one of which is deterministic and one that is stochastic. In particular, the deterministic model consists of an autonomous system of differential equations, whereas the stochastic system is a continuous time Markov chain. Connections between the two modeling regimes have been studied since the seminal paper by Kurtz (1972), where the deterministic model is shown to be a limit of a properly rescaled stochastic model over compact time intervals. Further, more recent studies have connected the long-term behaviors of the two models when the reaction network satisfies certain graphical properties, such as weak reversibility and a deficiency of zero. These connections have led some to conjecture a link between the long-term behavior of the two models exists, in some sense. In particular, one is tempted to believe that positive recurrence of all states for the stochastic model implies the existence of positive equilibria in the deterministic setting, and that boundary equilibria of the deterministic model imply the occurrence of an extinction event in the stochastic setting. We prove in this paper that these implications do not hold in general, even if restricting the analysis to networks that are bimolecular and that conserve the total mass. In particular, we disprove the implications in the special case of models that have absolute concentration robustness, thus answering in the negative a conjecture stated in the literature in 2014.


title: Probing the environment of high-z quasars using the proximity effect in projected quasar pairs

abstract: We have used spectra of 181 projected quasar pairs at separations $\le 1.5$ arcmins from the Sloan Digital Sky-Survey Data Release 12 in the redshift range of 2.5 to 3.5 to probe the proximity regions of the foreground quasars. We study the proximity effect both in the longitudinal as well as in the transverse directions, by carrying out a comparison of the \lya absorption lines originating from the vicinity of quasars to those originating from the general inter-galactic medium at the same redshift. We found an enhancement in the transmitted flux within 4 Mpc to the quasar in the longitudinal direction. However, the trend is found to be reversed in the transverse direction. In the longitudinal direction, we derived an excess overdensity profile showing an excess up to $r\leq 5$ Mpc after correcting for the quasar's ionization, taking into account the effect of low spectral resolution. This excess overdensity profile matches with the average overdensity profile in the transverse direction without applying any correction for the effect of the quasar's ionization. Among various possible interpretations, we found that the anisotropic obscuration of the quasar's ionization seems to be the most probable explanation. This is also supported by the fact that all of our foreground quasars happen to be Type-I AGNs. Finally, we constrain the average quasar's illumination along the transverse direction as compared to that along the longitudinal direction to be $\leq$27\% (3$\sigma$ confidence level).


title: Nodal points of Weyl semimetals survive the presence of moderate disorder

abstract: In this work we address the physics of individual three dimensional Weyl nodes subject to a moderate concentration of disorder. Previous analysis indicates the presence of a quantum phase transition below which disorder becomes irrelevant and the integrity of sharp nodal points of vanishing spectral density is preserved in this system. This statement appears to be at variance with the inevitable presence of statistically rare fluctuations which cannot be considered as weak and must have strong influence on the system's spectrum, no matter how small the average concentration. We here reconcile the two pictures by demonstrating that rare fluctuation potentials in the Weyl system generate a peculiar type of resonances which carry spectral density in any neighborhood of zero energy, but never at zero. In this way, the vanishing of the DoS for weak disorder survives the inclusion of rare events. We demonstrate this feature by considering three different models of disorder, each emphasizing specific aspects of the problem: a simplistic box potential model, a model with Gaussian distributed disorder, and one with a finite number of $s$-wave scatterers. Our analysis also explains why the protection of the nodal DoS may be difficult to see in simulations of finite size lattices.


title: Static Equilibria of Charged Particles Around Charged Black Holes: Chaos Bound and Its Violations

abstract: We study the static equilibrium of a charged massive particle around a charged black hole, balanced by the Lorentz force. For a given black hole, the equilibrium surface is determined by the charge/mass ratio of the particle. By investigating a large class of charged black holes, we find that the equilibria can be stable, marginal or unstable. We focus on the unstable equilibria which signal chaotic motions and we obtain the corresponding Lyapunov exponents $\lambda$. We find that although $\lambda$ approaches universally the horizon surface gravity $\kappa$ when the equilibria are close to the horizon, the proposed chaotic motion bound $\lambda<\kappa$ is satisfied only by some specific black holes including the RN and RN-AdS black holes. The bound can be violated by a large number of black holes including the RN-dS black holes or black holes in Einstein-Maxwell-Dilaton, Einstein-Born-Infeld and Einstein-Gauss-Bonnet-Maxwell gravities. We find that unstable equilibria can even exist in extremal black holes, implying that the ratio $\lambda/\kappa$ can be arbitrarily large for sufficiently small $\kappa$. Our investigation does suggest a universal bound for sufficiently large $\kappa$, namely $\lambda/\kappa <{\cal C}$ for some order-one constant ${\cal C}$.


title: Europa's Optical Aurora

abstract: Auroral emissions provide opportunities to study the tenuous atmospheres of Solar System satellites, revealing the presence and abundance of molecular and atomic species as well as their spatial and temporal variability. Far-UV aurorae have been used for decades to study the atmospheres of the galilean satellites. Here we present the first detection of Europa's visible-wavelength atomic oxygen aurora at 6300/6364 \AA{} arising from the metastable O$(^1$D) state, observed with the Keck I and Hubble Space Telescopes while Europa was in eclipse by Jupiter on six occasions in February-April 2018. The disk-integrated O($^1$D) brightness varies from $<$500 R up to more than 2 kR between dates, a factor of 15 higher than the OI 1356 \AA{} brightness on average. The ratio of emission at 6300/5577 \AA{} is diagnostic of parent molecule; the 5577 \AA{} emission was not detected in our dataset, which favors O$_2$ as the dominant atmospheric constituent and rules out an O/O$_2$ mixing ratio above 0.35. For an O$_2$ atmosphere and typical plasma conditions at Europa's orbit, the measured surface brightness range corresponds to column densities of 1-9$\times$10$^{14}$ cm$^{-2}$.


title: Global Convergence of Stochastic Gradient Hamiltonian Monte Carlo for Non-Convex Stochastic Optimization: Non-Asymptotic Performance Bounds and Momentum-Based Acceleration

abstract: Stochastic gradient Hamiltonian Monte Carlo (SGHMC) is a variant of stochastic gradient with momentum where a controlled and properly scaled Gaussian noise is added to the stochastic gradients to steer the iterates towards a global minimum. Many works reported its empirical success in practice for solving stochastic non-convex optimization problems, in particular it has been observed to outperform overdamped Langevin Monte Carlo-based methods such as stochastic gradient Langevin dynamics (SGLD) in many applications. Although asymptotic global convergence properties of SGHMC are well known, its finite-time performance is not well-understood. In this work, we study two variants of SGHMC based on two alternative discretizations of the underdamped Langevin diffusion. We provide finite-time performance bounds for the global convergence of both SGHMC variants for solving stochastic non-convex optimization problems with explicit constants. Our results lead to non-asymptotic guarantees for both population and empirical risk minimization problems. For a fixed target accuracy level, on a class of non-convex problems, we obtain complexity bounds for SGHMC that can be tighter than those for SGLD. These results show that acceleration with momentum is possible in the context of global non-convex optimization.


title: Polarization Symmetry Breaking in Nitrogen under High Pressure

abstract: An order parameter with broken polarization symmetry is proposed for the pressure-induced molecular dissociation transitions in nitrogen. The experimental dissociative transition pressure dependence of the calculated order parameter is well described by a power law of the form P=kq2. Our results unveil the important role of symmetry breaking in nitrogen under high pressure and further demonstrate a path from dense molecular fluid to fully polymeric phases by breaking local polarization reversal symmetry.


title: Elasticity and Stability of Shape Changing Structures

abstract: As we enter the age of designer matter - where objects can morph and change shape on command - what tools do we need to create shape-shifting structures? At the heart of an elastic deformation is the combination of dilation and distortion, or stretching and bending. The competition between the latter can cause elastic instabilities, and over the last fifteen years these instabilities have provided a multitude of ways to prescribe and control shape change. Buckling, wrinkling, folding, creasing, snapping have become mechanisms that when harmoniously combined enable mechanical metamaterials, self-folding origami, ultralight and ultrathin kirigami, and structures that appear to grow from one shape to another. In this review, I aim to connect the fundamentals of elastic instabilities to the advanced functionality currently found within mechanical metamaterials.


title: A Two-Step Learning Method For Detecting Landmarks on Faces From Different Domains

abstract: The detection of fiducial points on faces has significantly been favored by the rapid progress in the field of machine learning, in particular in the convolution networks. However, the accuracy of most of the detectors strongly depends on an enormous amount of annotated data. In this work, we present a domain adaptation approach based on a two-step learning to detect fiducial points on human and animal faces. We evaluate our method on three different datasets composed of different animal faces (cats, dogs, and horses). The experiments show that our method performs better than state of the art and can use few annotated data to leverage the detection of landmarks reducing the demand for large volume of annotated data.


title: The Three Hundred project: a large catalogue of theoretically modelled galaxy clusters for cosmological and astrophysical applications

abstract: We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.


title: Global solutions for a chemotaxis hyperbolic-parabolic system on networks with nonhomogeneous boundary conditions

abstract: In this paper we study a semilinear hyperbolic-parabolic system as a model for some chemotaxis phenomena evolving on networks; we consider transmission conditions at the inner nodes which preserve the fluxes and non- homogeneous boundary conditions having in mind phenomena with inflow of cells and food providing at the network exits. We give some conditions on the boundary data which ensure the existence of stationary solutions and we prove that these ones are asymptotic profiles for a class of global solutions.


title: Visual-Quality-Driven Learning for Underwater Vision Enhancement

abstract: The image processing community has witnessed remarkable advances in enhancing and restoring images. Nevertheless, restoring the visual quality of underwater images remains a great challenge. End-to-end frameworks might fail to enhance the visual quality of underwater images since in several scenarios it is not feasible to provide the ground truth of the scene radiance. In this work, we propose a CNN-based approach that does not require ground truth data since it uses a set of image quality metrics to guide the restoration learning process. The experiments showed that our method improved the visual quality of underwater images preserving their edges and also performed well considering the UCIQE metric.


title: Quark jets scattering from a gluon field: from saturation to high $p_t$

abstract: We continue our studies of possible generalization of the Color Glass Condensate (CGC) effective theory of high energy QCD to include the high $p_t$ (or equivalently large $x$) QCD dynamics as proposed in [JJM-elastic]. Here we consider scattering of a quark from both the small and large $x$ gluon degrees of freedom in a proton or nucleus target and derive the full scattering amplitude by including the interactions between the small and large $x$ gluons of the target. We thus generalize the standard eikonal approximation for parton scattering which can now be deflected by a large angle (and therefore have large $p_t$) and also lose a significant fraction of its longitudinal momentum (unlike the eikonal approximation). The corresponding production cross section can thus serve as the starting point toward derivation of a general evolution equation that would contain DGLAP evolution equation at large $Q^2$ and the JIMWLK evolution equation at small $x$. This amplitude can also be used to construct the quark Feynman propagator which is the first ingredient needed to generalize the Color Glass Condensate (CGC) effective theory of high energy QCD to include the high $p_t$ dynamics. We outline how it can be used to compute observables in the large $x$ (high $p_t$) kinematic region where the standard Color Glass Condensate formalism breaks down.


title: Iterative phase retrieval in coherent diffractive imaging: practical issues

abstract: In this work, issues in phase retrieval in the coherent diffractive imaging (CDI) technique, from discussion on parameters for setting up a CDI experiment to evaluation of the goodness of the final reconstruction, are discussed. The distribution of objects under study by CDI often cannot be cross-validated by another imaging technique. It is therefore important to make sure that the developed CDI procedure delivers an artifact-free object reconstruction. Critical issues that can lead to artifacts are presented and recipes on how to avoid them are provided.

categories: physics.optics eess.IV physics.comp-ph physics.data-an


title: The Main Decomposition of Finite-Dimensional Protori

abstract: A protorus is a compact connected abelian group. We use a result on finite rank torsion-free abelian groups and Pontryagin Duality to considerably generalize a well-known factorization of a finite-dimensional protorus into a product of a torus and a torus-free complementary factor. We also classify by types the solenoids of Hewitt and Ross.


title: On $p$-adic valuations of colored $p$-ary partitions

abstract: Let $m\in\N_{\geq 2}$ and for given $k\in\N_{+}$ consider the sequence $(A_{m,k}(n))_{n\in\N}$ defined by the power series expansion $$ \prod_{n=0}^{\infty}\frac{1}{\left(1-x^{m^{n}}\right)^{k}}=\sum_{n=0}^{\infty}A_{m,k}(n)x^{n}. $$ The number $A_{m,k}(n)$ counts the number of representations of $n$ as sums of powers of $m$, where each summand has one among $k$ colors. In this note we prove that for each $p\in\mathbb{P}_{\geq 3}$ and $s\in\N_{+}$, the $p$-adic valuation of the number $A_{p,(p-1)(p^s-1)}(n)$ is equal to 1 for $n\geq p^s$. We also obtain some results concerning the behaviour of the sequence $(\nu_{p}(A_{p,(p-1)(up^s-1)}(n)))_{n\in\N}$ for fixed $u\in\{2,\ldots,p-1\}$ and $p\geq 3$. Our results generalize the earlier findings obtained for $p=2$ by Gawron, Miska and the first author.


title: Occlusion-Aware Risk Assessment for Autonomous Driving in Urban Environments

abstract: Navigating safely in urban environments remains a challenging problem for autonomous vehicles. Occlusion and limited sensor range can pose significant challenges to safely navigate among pedestrians and other vehicles in the environment. Enabling vehicles to quantify the risk posed by unseen regions allows them to anticipate future possibilities, resulting in increased safety and ride comfort. This paper proposes an algorithm that takes advantage of the known road layouts to forecast, quantify, and aggregate risk associated with occlusions and limited sensor range. This allows us to make predictions of risk induced by unobserved vehicles even in heavily occluded urban environments. The risk can then be used either by a low-level planning algorithm to generate better trajectories, or by a high-level one to plan a better route. The proposed algorithm is evaluated on intersection layouts from real-world map data with up to five other vehicles in the scene, and verified to reduce collision rates by 4.8x comparing to a baseline method while improving driving comfort.


title: Magnetic Nanoparticles in a Nematic Channel: A One-Dimensional Study

abstract: We study a ferromagnetic suspension or a suspension of magnetic nanoparticles in an anisotropic nematic medium, in three different one-dimensional variational settings, ordered in terms of increasing complexity. The three models are featured by a nematic energy, a magnetic energy and a magneto-nematic coupling energy and the experimentally observed patterns are modelled as local or global energy minimizers. We numerically observe polydomains with distinct states of magnetization for weak to moderate magneto-nematic coupling in our models. We demonstrate that these polydomains are stabilised by lowering the temperature (as in Mertelj et al., 2013) and that the polydomain structures lose stability as the magneto-nematic coupling increases. Some exact solutions for prototypical situations are also obtained.


title: Coronal hard X-ray sources revisited

abstract: This paper reports on the re-analysis of solar flares in which the hard X-rays (HXRs) come predominantly from the corona rather than from the more usual chromospheric footpoints. All of the 26 previously analyzed event time intervals, over 13 flares, are re-examined for consistency with a flare model in which electrons are accelerated near the top of a magnetic loop that has a sufficiently high density to stop most of the electrons by Coulomb collisions before they can reach the footpoints. Of particular importance in the previous analysis was the finding that the length of the coronal HXR source increased with energy in the 20 - 30 keV range. However, after allowing for the possibility that footpoint emission at the higher energies affects the inferred length of the coronal HXR source, and using analysis techniques that suppress the possible influence of such footpoint emission, we conclude that there is no longer evidence that the length of the HXR coronal sources increase with increasing energy. In fact, for the 6 flares and 12 time intervals that satisfied our selection criteria, the loop lengths decreased on average by 1.0 +/- 0.2 arcsec between 20 and 30 keV, with a standard deviation of 3.5 arcsec. We find strong evidence that the peak of the coronal HXR source increases in altitude with increasing energy. For the thermal component of the emission, this is consistent with the standard CHSKP flare model in which magnetic reconnection in a coronal current sheet results in new hot loops being formed at progressively higher altitudes. The explanation for the nonthermal emission is not so clear.


title: Efficient Global Optimization using Deep Gaussian Processes

abstract: Efficient Global Optimization (EGO) is widely used for the optimization of computationally expensive black-box functions. It uses a surrogate modeling technique based on Gaussian Processes (Kriging). However, due to the use of a stationary covariance, Kriging is not well suited for approximating non stationary functions. This paper explores the integration of Deep Gaussian processes (DGP) in EGO framework to deal with the non-stationary issues and investigates the induced challenges and opportunities. Numerical experimentations are performed on analytical problems to highlight the different aspects of DGP and EGO.


title: Association and Simpson conversion in $2 \times 2 \times 2$ contingency tables

abstract: We study a generalisation of Simpson reversal (also known as Simpson's paradox or the Yule-Simpson effect) to $2 \times 2 \times 2$ contingency tables and characterise the cases for which it can and cannot occur with two combinatorial-geometric lemmas. We also present a conjecture based on some computational experiments on the expected likelihood of such events.


title: Role of electron inertia and electron/ion finite Larmor radius effects in low-beta, magneto-Rayleigh-Taylor instability

abstract: The magneto-Rayleigh-Taylor (MRT) instability has been investigated in great detail in previous work using magnetohydrodynamic and kinetic models for low-beta plasmas. The work presented here extends previous studies of this instability to regimes where finite-Larmor-Radius (FLR) effects may be important. Comparisons of the MRT instability are made using a 5-moment and a 10-moment two-fluid model, the two fluids being ions and electrons. The 5-moment model includes Hall stabilization whereas the 10-moment model includes Hall and FLR stabilization. Results are presented for these two models using different electron mass to understand the role of electron inertia in the late-time nonlinear evolution of the MRT instability. For the 5-moment model, the late-time nonlinear MRT evolution does not significantly depend on the electron inertia. However, when FLR stabilization is important, the 10-moment results show that a lower ion-to-electron mass ratio (i.e. larger electron inertia) under-predicts the energy in high-wavenumber modes due to larger FLR stabilization.


title: Au Nanocluster Growth on Graphene Supported on Ni(111)

abstract: Low energy alkali ion scattering is used to investigate the deposition of Au onto a single layer of graphene grown onto Ni(111) by chemical vapor deposition. The yield of 3.0 keV Na$^+$ singly scattered from Au as a function of coverage indicates that it grows in a Volmer-Weber mode forming nanoclusters that increase in size with the amount of deposition. The neutralization probability of the scattered Na$^+$ is high for the smallest clusters and decreases as they increase in size. This is presumably caused by the cluster edge atoms being positively charged combined with the fact that the ratio of edge to center atoms decreases with size, which is similar to the behavior of Au nanoclusters on oxide substrates. In addition, oxygen is intercalated under the graphene film to decouple it from the substrate, but no changes in the growth mode or neutralization probability are observed.


title: Connection Coefficients for Higher-order Bernoulli and Euler Polynomials: A Random Walk Approach

abstract: We consider the use of random walks as an approach to obtain connection coefficients for higher-order Bernoulli and Euler polynomials. In particular, we consider the cases of a $1$-dimensional linear reflected Brownian motion and of a $3$-dimensional Bessel process. Considering the successive hitting times of two, three, and four fixed levels by these random walks yields non-trivial identities that involve higher-order Bernoulli and Euler polynomials.


title: Electrically Tunable Low Density Superconductivity in a Monolayer Topological Insulator

abstract: The capability to switch electrically between superconducting and insulating states of matter represents a novel paradigm in the state-of-the-art engineering of correlated electronic systems. An exciting possibility is to turn on superconductivity in a topologically non-trivial insulator, which provides a route to search for non-Abelian topological states. However, existing demonstrations of superconductor-insulator switches have involved only topologically trivial systems, and even those are rare due to the stringent requirement to tune the carrier density over a wide range. Here we report reversible, in-situ electrostatic on off switching of superconductivity in a recently established quantum spin Hall insulator, namely monolayer tungsten ditelluride (WTe2). Fabricated into a van der Waals field effect transistor, the monolayer's ground state can be continuously gate-tuned from the topological insulating to the superconducting state, with critical temperatures Tc up to ~ 1 Kelvin. The critical density for the onset of superconductivity is estimated to be ~ 5 x 10^12 cm^-2, among the lowest for two-dimensional (2D) superconductors. Our results establish monolayer WTe2 as a material platform for engineering novel superconducting nanodevices and topological phases of matter.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el cond-mat.supr-con


title: 3d TQFTs from Argyres-Douglas theories

abstract: We construct a new class of three-dimensional topological quantum field theories (3d TQFTs) by considering generalized Argyres-Douglas theories on $S^1 \times M_3$ with a non-trivial holonomy of a discrete global symmetry along the $S^1$. For the minimal choice of the holonomy, the resulting 3d TQFTs are non-unitary and semisimple, thus distinguishing themselves from theories of Chern-Simons and Rozansky-Witten types respectively. Changing the holonomy performs a Galois transformation on the TQFT, which can sometimes give rise to more familiar unitary theories such as the $(G_2)_1$ and $(F_4)_1$ Chern-Simons theories. Our construction is based on an intriguing relation between topologically twisted partition functions, wild Hitchin characters, and chiral algebras which, when combined together, relate Coulomb branch and Higgs branch data of the same 4d $\mathcal{N}=2$ theory. We test our proposal by applying localization techniques to the conjectural $\mathcal{N}=1$ UV Lagrangian descriptions of the $(A_1,A_2)$, $(A_1,A_3)$ and $(A_1,D_3)$ theories.

categories: hep-th math-ph math.AG math.MP math.QA math.RT


title: Non-equilibrium dynamics of isostatic spring networks

abstract: Marginally stable systems exhibit rich critical mechanical behavior. Such isostatic assemblies can be driven out of equilibrium by internal activity, but it remains unclear how the isostatic and critical nature of such systems affects their non-equilibrium dynamics. Here, we investigate the influence of the isostatic threshold on the non-equilibrium dynamics of active diluted spring networks. In our model, heterogeneously distributed active noise sources drive the system into a non-equilibrium steady state. We quantify the non-equilibrium dynamics of nearest-neighbor network nodes by the characteristic cycling frequency $\omega$---a measure of the circulation of the associated phase space currents. The distribution of these nearest-neighbor cycling frequencies exhibits critical scaling, which we describe using a mean-field approach. Overall, our work provides a theoretical approach to elucidate the role of marginality in active disordered systems.


title: Jump to better conclusions: SCAN both left and right

abstract: Lake and Baroni (2018) recently introduced the SCAN data set, which consists of simple commands paired with action sequences and is intended to test the strong generalization abilities of recurrent sequence-to-sequence models. Their initial experiments suggested that such models may fail because they lack the ability to extract systematic rules. Here, we take a closer look at SCAN and show that it does not always capture the kind of generalization that it was designed for. To mitigate this we propose a complementary dataset, which requires mapping actions back to the original commands, called NACS. We show that models that do well on SCAN do not necessarily do well on NACS, and that NACS exhibits properties more closely aligned with realistic use-cases for sequence-to-sequence models.


title: Emergence of micro-frequency comb via limit cycles in dissipatively coupled condensates

abstract: Self-sustained oscillations, limit cycles, are a fundamental phenomenon unique to nonlinear dynamic systems of high-dimensional phase space. They enable understanding of a wide range of cyclic processes in natural, social and engineering systems. Here we show that limit cycles form in coupled polariton cavities following the breaking of Josephson coupling, leading to frequency-comb emission. The limit cycles and destruction of Josephson coupling both appear due to interplay between strong polariton-polariton interaction and a dissipative contribution to the cavity coupling. The resulting nonlinear dynamics of the condensates is characterized by asymmetric population distribution and nontrivial average phase difference between the two condensates, and by time-periodic modulation of their amplitudes and phases. The latter is manifested by coherent emission of new equidistant frequency components. The emission spectrum resembles that of a micro-frequency comb, but originates from a fundamentally different mechanism than that of existing frequency combs. It allows non-resonant excitation with a power input much below the conventional semiconductor laser threshold. The comb line spacing is determined by the interaction and coupling strengths, and is adjustable up to multi-terahertz frequency. The work establishes coupled polariton cavities as an experimental platform for rich nonlinear dynamic phenomena.


title: Deep Spectral Correspondence for Matching Disparate Image Pairs

abstract: A novel, non-learning-based, saliency-aware, shape-cognizant correspondence determination technique is proposed for matching image pairs that are significantly disparate in nature. Images in the real world often exhibit high degrees of variation in scale, orientation, viewpoint, illumination and affine projection parameters, and are often accompanied by the presence of textureless regions and complete or partial occlusion of scene objects. The above conditions confound most correspondence determination techniques by rendering impractical the use of global contour-based descriptors or local pixel-level features for establishing correspondence. The proposed deep spectral correspondence (DSC) determination scheme harnesses the representational power of local feature descriptors to derive a complex high-level global shape representation for matching disparate images. The proposed scheme reasons about correspondence between disparate images using high-level global shape cues derived from low-level local feature descriptors. Consequently, the proposed scheme enjoys the best of both worlds, i.e., a high degree of invariance to affine parameters such as scale, orientation, viewpoint, illumination afforded by the global shape cues and robustness to occlusion provided by the low-level feature descriptors. While the shape-based component within the proposed scheme infers what to look for, an additional saliency-based component dictates where to look at thereby tackling the noisy correspondences arising from the presence of textureless regions and complex backgrounds. In the proposed scheme, a joint image graph is constructed using distances computed between interest points in the appearance (i.e., image) space. Eigenspectral decomposition of the joint image graph allows for reasoning about shape similarity to be performed jointly, in the appearance space and eigenspace.


title: On Boltzmann vs. Gibbs and the Equilibrium in Statistical Mechanics

abstract: In a recent article, Werndl and Frigg discuss the relationship between the Boltzmannian and Gibbsian framework of statistical mechanics, addressing in particular the question when equilibrium values calculated in both frameworks agree. In this paper, I address conceptual confusions that could arise from their discussion, concerning in particular the authors' use of "Boltzmann equilibrium". I also clarify the status of the Khinchin condition for the equivalence of Boltzmannian and Gibbsian, and show that it follows under the assumptions proposed by Werndl and Frigg from standard arguments in probability theory.


title: The closed loop between opinion formation and personalised recommendations

abstract: In online platforms, recommender systems are responsible for directing users to relevant contents. In order to enhance the users' engagement, recommender systems adapt their output to the reactions of the users, who are in turn affected by the recommended contents. In this work, we study a tractable analytical model of a user that interacts with an online news aggregator, with the purpose of making explicit the feedback loop between the evolution of the user's opinion and the personalised recommendation of contents. More specifically, we assume that the user is endowed with a scalar opinion about a certain issue and seeks news about it on a news aggregator: this opinion is influenced by all received news, which are characterized by a binary position on the issue at hand. The user is affected by a confirmation bias, that is, a preference for news that confirm her current opinion. The news aggregator recommends items with the goal of maximizing the number of user's clicks (as a measure of her engagement): in order to fulfil its goal, the recommender has to compromise between exploring the user's preferences and exploiting what it has learned so far. After defining suitable metrics for the effectiveness of the recommender systems (such as the click-through rate) and for its impact on the opinion, we perform both extensive numerical simulations and a mathematical analysis of the model. We find that personalised recommendations markedly affect the evolution of opinions and favor the emergence of more extreme ones: the intensity of these effects is inherently related to the effectiveness of the recommender. We also show that by tuning the amount of randomness in the recommendation algorithm, one can seek a balance between the effectiveness of the recommendation system and its impact on the opinions.


title: Computational Arithmetic of Modular Forms (Course notes)

abstract: These course notes are about computing modular forms and some of their arithmetic properties. Their aim is to explain and prove the modular symbols algorithm in as elementary and as explicit terms as possible, and to enable the devoted student to implement it over any ring (such that a sufficient linear algebra theory is available in the chosen computer algebra system). The chosen approach is based on group cohomology and along the way the needed tools from homological algebra are provided.


title: Arrow(s) of Time without a Past Hypothesis

abstract: The paper discusses recent proposals by Carroll and Chen, as well as Barbour, Koslowski, and Mercati to explain the (thermodynamic) arrow of time without a Past Hypothesis, i.e., the assumption of a special (low-entropy) initial state of the universe. After discussing the role of the Past Hypothesis and the controversy about its status, we explain why Carroll's model - which establishes an arrow of time as typical - can ground sensible predictions and retrodictions without assuming something akin to a Past Hypothesis. We then propose a definition of a Boltzmann entropy for a classical $N$-particle system with gravity, suggesting that a Newtonian gravitating universe might provide a relevant example of Carroll's entropy model. This invites comparison with the work of Barbour, Koslowski, and Mercati that identifies typical arrows of time in a relational formulation of classical gravity on shape space. We clarify the difference between this gravitational arrow in terms of shape complexity and the entropic arrow in absolute spacetime and work out the key advantages of the relationalist theory. We end by pointing out why the entropy concept relies on absolute scales and is thus not relational.


title: Characterization of queer supercrystals

abstract: We provide a characterization of the crystal bases for the quantum queer superalgebra recently introduced by Grantcharov et al.. This characterization is a combination of local queer axioms generalizing Stembridge's local axioms for crystal bases for simply-laced root systems, which were recently introduced by Assaf and Oguz, with further axioms and a new graph $G$ characterizing the relations of the type $A$ components of the queer crystal. We provide a counterexample to Assaf's and Oguz' conjecture that the local queer axioms uniquely characterize the queer supercrystal. We obtain a combinatorial description of the graph $G$ on the type $A$ components by providing explicit combinatorial rules for the odd queer operators on certain highest weight elements.


title: Reversing the asymmetry in data exfiltration

abstract: Preventing data exfiltration from computer systems typically depends on perimeter defences, but these are becoming increasingly fragile. Instead we suggest an approach in which each at-risk document is supplemented by many fake versions of itself. An attacker must either exfiltrate all of them; or try to discover which is the real one while operating within the penetrated system, and both are difficult. Creating and maintaining many fakes is relatively inexpensive, so the advantage that typically accrues to an attacker now lies with the defender. We show that algorithmically generated fake documents are reasonably difficult to detect using algorithmic analytics.


title: Semantic WordRank: Generating Finer Single-Document Summarizations

abstract: We present Semantic WordRank (SWR), an unsupervised method for generating an extractive summary of a single document. Built on a weighted word graph with semantic and co-occurrence edges, SWR scores sentences using an article-structure-biased PageRank algorithm with a Softplus function adjustment, and promotes topic diversity using spectral subtopic clustering under the Word-Movers-Distance metric. We evaluate SWR on the DUC-02 and SummBank datasets and show that SWR produces better summaries than the state-of-the-art algorithms over DUC-02 under common ROUGE measures. We then show that, under the same measures over SummBank, SWR outperforms each of the three human annotators (aka. judges) and compares favorably with the combined performance of all judges.


title: Analysis of variable-step/non-autonomous artificial compression methods

abstract: A standard artificial compression (AC) method for incompressible flow is $$ \frac{u_{n+1}^{\varepsilon }-u_{n}^{\varepsilon }}{k}+u_{n+1}^{\varepsilon }\cdot \nabla u_{n+1}^{\varepsilon }+{\frac{1}{2}}u_{n+1}^{\varepsilon }\nabla \cdot u_{n+1}^{\varepsilon }+\nabla p_{n+1}^{\varepsilon }-\nu \Delta u_{n+1}^{\varepsilon }=f\text{ ,} \\ \varepsilon \frac{p_{n+1}^{\varepsilon }-p_{n}^{\varepsilon }}{k} +\nabla \cdot u_{n+1}^{\varepsilon }=0 $$ for, typically, $\varepsilon =k$ (timestep). It is fast, efficient and stable with accuracy $O(\varepsilon +k)$. For adaptive (and thus variable) timestep $k_{n}$ (and thus $\varepsilon =\varepsilon _{n}$) its long time stability is unknown. For variable $k,\varepsilon $ this report shows how to adapt a standard AC method to recover a provably stable method. For the associated continuum AC model, we prove convergence of the $\varepsilon =\varepsilon (t)$\ artificial compression model to a weak solution of the incompressible Navier-Stokes equations as $\varepsilon =\varepsilon (t)\rightarrow 0$. The analysis is based on space-time Strichartz estimates for a non-autonomous acoustic equation. Variable $\varepsilon ,k$ numerical tests in $2d$ and $3d$ are given for the new AC method.


title: Feynman's solution of the quintessential problem in solid state physics

abstract: Two of the most influential ideas developed by Richard Feynman are the Feynman diagram technique and his variational approach. The former provides a powerful tool to construct a systematic expansion for a generic interacting system, while the latter allows optimization of a perturbation theory using a variational principle. Here we show that combining a variational approach with a new diagrammatic quantum Monte Carlo method, both based on the Feynman's original ideas, results in a powerful and accurate solver to the generic solid state problem, in which a macroscopic number of electrons interact by the long range Coulomb repulsion. We apply the solver to the quintessential problem of solid state, the uniform electron gas (UEG), which is at the heart of the density functional theory (DFT) success in describing real materials, yet it has not been adequately solved for over 90 years. While some wave-function properties, like the ground state energy, have been very accurately calculated by the diffusion Monte Carlo method (DMC), the static and dynamic response functions, which are directly accessed by the experiment, and are crucial for development of non-local and non-adiabatic density functionals, remain poorly understood. Our method allows us to calculate the momentum-frequency resolved spin response functions for the first time, and to improve on the precision of the charge response function. The accuracy of both response functions is sufficiently high, so as to uncover previously missed fine structure in these responses. This method can be straightforwardly applied to a large number of moderately interacting electron systems in the thermodynamic limit, including realistic models of metallic and semiconducting solids.


title: Conversion of any finite bandwidth optical field into a shape invariant beam

abstract: We discuss a method to transform any optical field, with finite frequency bandwidth, into a shape invariant beam with transverse scaling, dependent on the propagation distance. The method consists in modulating the field with a quadratic phase of appropriate curvature radius. As a particular application, we employ the method to extend the existence region of a finite non-diffracting field to an unbounded domain.


title: The major merger origin of the Andromeda II kinematics

abstract: Prolate rotation (i.e. rotation around the long axis) has been reported for two Local Group dwarf galaxies: Andromeda II, a dwarf spheroidal satellite of M31, and Phoenix, a transition type dwarf galaxy. The prolate rotation may be an exceptional indicator of a past major merger between dwarf galaxies. We showed that this type of rotation cannot be obtained in the tidal stirring scenario, in which the satellite is transformed from disky to spheroidal by tidal forces of the host galaxy. However, we successfully reproduced the observed Andromeda II kinematics in controlled, self-consistent simulations of mergers between equal-mass disky dwarf galaxies on a radial or close-to-radial orbit. In simulations including gas dynamics, star formation and ram pressure stripping, we are able to reproduce more of the observed properties of Andromeda II: the unusual rotation, the bimodal star formation history and the spatial distribution of the two stellar populations, as well as the lack of gas. We support this scenario by demonstrating the merger origin of prolate rotation in the cosmological context for sufficiently resolved galaxies in the Illustris large-scale cosmological hydrodynamical simulation.


title: Inversion and magnetic quantum oscillations in symmetric periodic Anderson model

abstract: We study the symmetric periodic Anderson model of the conduction electrons hybridized with the localized correlated electrons on square lattice. Using the canonical representation of electrons by Kumar, we do a self-consistent theory of its effective charge and spin dynamics, which produces an insulating ground state that undergoes continuous transition from the Kondo singlet to N\'eel antiferromagnetic phase with decreasing hybridization, and uncovers two inversion transitions for the charge quasiparticles. With suitably inverted quasiparticle bands for moderate to weaker effective Kondo couplings, this effective charge dynamics in the magnetic field coupled to electronic motion produces magnetic quantum oscillations with frequency corresponding to the half Brillouin zone.


title: Vertical dispersion of Lagrangian tracers in fully developed stably stratified turbulence

abstract: We study the effect of different forcing functions and of the local gradient Richardson number Rig on the vertical mixing of Lagrangian tracers in stably stratified turbulence under the Boussinesq approximation, and present a wave and continuous-time random walk model for single- and two- particle vertical dispersion. The model consists of a random superposition of linear waves with their amplitude based on the observed Lagrangian spectrum of vertical velocity, and a random walk process to capture overturning that depends on the statistics of Rig among other Eulerian quantities. The model is in good agreement with direct numerical simulations of stratified turbulence, where single- and two-particle dispersion differs from the homogeneous and isotropic case. Moreover, the model gives insight into the mixture of linear and non-linear physics in the problem, as well as on the different processes responsible for vertical turbulent dispersion.


title: A Conceptual Approach to Complex Model Management with Generalized Modelling Patterns and Evolutionary Identification

abstract: Complex systems' modeling and simulation are powerful ways to investigate a multitude of natural phenomena providing extended knowledge on their structure and behavior. However, enhanced modeling and simulation require integration of various data and knowledge sources, models of various kinds (data-driven models, numerical models, simulation models, etc.), intelligent components in one composite solution. Growing complexity of such composite model leads to the need of specific approaches for management of such model. This need extends where the model itself becomes a complex system. One of the important aspects of complex model management is dealing with the uncertainty of various kinds (context, parametric, structural, input/output) to control the model. In the situation where a system being modeled, or modeling requirements change over time, specific methods and tools are needed to make modeling and application procedures (meta-modeling operations) in an automatic manner. To support automatic building and management of complex models we propose a general evolutionary computation approach which enables managing of complexity and uncertainty of various kinds. The approach is based on an evolutionary investigation of model phase space to identify the best model's structure and parameters. Examples of different areas (healthcare, hydrometeorology, social network analysis) were elaborated with the proposed approach and solutions.


title: Magnetically induced termination of giant planet formation

abstract: Here a physical model for terminating giant planet formation is outlined and compared to other methods of late-stage giant planet formation. As has been pointed out before, gas accreting into a gap and onto the planet will encounter the planetary dynamo-generated magnetic field. The planetary magnetic field produces an effective cross section through which gas is accreted. Gas outside this cross section is recycled into the protoplanetary disk, hence only a fraction of mass that is accreted into the gap remains bound to the planet. This cross section inversely scales with the planetary mass, which naturally leads to stalled planetary growth late in the formation process. We show that this method naturally leads to Jupiter-mass planets and does not invoke any artificial truncation of gas accretion, as has been done in some previous population synthesis models. The mass accretion rate depends on the radius of the growing planet after the gap has opened, and we show that so-called hot-start planets tend to become more massive than cold-start planets. When this result is combined with population synthesis models, it might show observable signatures of cold-start versus hot-start planets in the exoplanet population.


title: A series of series topologies on $\mathbb{N}$

abstract: Each series $\sum_{n=1}^\infty a_n$ of real positive terms gives rise to a topology on $\mathbb{N} = \{1,2,3,...\}$ by declaring a proper subset $A\subseteq \mathbb{N}$ to be closed if $\sum_{n\in A} a_n < \infty$. We explore the relationship between analytic properties of the series and topological properties on $\mathbb{N}$. In particular, we show that, up to homeomorphism, $|\mathbb{R}|$-many topologies are generated. We also find an uncountable family of examples $\{\mathbb{N}_\alpha\}_{\alpha \in [0,1]}$ with the property that for any $\alpha < \beta$, there is a continuous bijection $\mathbb{N}_\beta\rightarrow \mathbb{N}_\alpha$, but the only continuous functions $\mathbb{N}_\alpha\rightarrow \mathbb{N}_\beta$ are constant.


title: Are object detection assessment criteria ready for maritime computer vision?

abstract: Maritime vessels equipped with visible and infrared cameras can complement other conventional sensors for object detection. However, application of computer vision techniques in maritime domain received attention only recently. The maritime environment offers its own unique requirements and challenges. Assessment of the quality of detections is a fundamental need in computer vision. However, the conventional assessment metrics suitable for usual object detection are deficient in the maritime setting. Thus, a large body of related work in computer vision appears inapplicable to the maritime setting at the first sight. We discuss the problem of defining assessment metrics suitable for maritime computer vision. We consider new bottom edge proximity metrics as assessment metrics for maritime computer vision. These metrics indicate that existing computer vision approaches are indeed promising for maritime computer vision and can play a foundational role in the emerging field of maritime computer vision.


title: Search for a time-varying electron antineutrino signal at Daya Bay

abstract: A search for a time-varying $\bar{\nu}_{e}$ signal was performed with 621 days of data acquired by the Daya Bay Reactor Neutrino Experiment over 704 calendar days. The time spectrum of the measured $\overline{\nu}_e$ flux normalized to its prediction was analyzed with a Lomb-Scargle periodogram, which yielded no significant signal for periods ranging from 2 hours to nearly 2 years. The normalized time spectrum was also fit for a sidereal modulation under the Standard Model extension (SME) framework to search for Lorentz and CPT violation (LV-CPTV). Limits were obtained for all six flavor pairs $\bar{e}\bar{\mu}$, $\bar{e}\bar{\tau}$, $\bar{\mu}\bar{\tau}$, $\bar{e}\bar{e},\bar{\mu}\bar{\mu}$ and $\bar{\tau}\bar{\tau}$ by fitting them one at a time, constituting the first experimental constraints on the latter three. Daya Bay's high statistics and unique layout of multiple directions from three pairs of reactors to three experimental halls allowed the simultaneous constraint of individual SME LV-CPTV coefficients without assuming others contribute negligibly, a first for a neutrino experiment.


title: A general method for transforming non-physical configurations in BPS states

abstract: In this work, we apply the so-called BPS method in order to obtain topological defects for a complex scalar field Lagrangian introduced by Trullinger and Subbaswamy. The BPS approach led us to compute new analytical solutions for this model.In our investigation, we found analytical configurations which satisfy the BPS first-order differential equations but do not obey the equations of motion of the model. Such defects were named non-physical ones. In order to recover the physical meaning of these defects, we proposed a procedure which can transform them into BPS states of new scalar field models. The new models here founded were applied in the context of hybrid cosmological scenarios, where we derived cosmological parameters compatible with the observed Universe. Such a methodology opens a new window to connect different two scalar fields systems and can be implemented in several distinct applications such as Bloch Branes, Lorentz and Symmetry Breaking Scenarios, Q-Balls, Oscillons, Cosmological Contexts, and Condensed Matter Systems.


title: Compact QC-LDPC Block and SC-LDPC Convolutional Codes for Low-Latency Communications

abstract: Low decoding latency and complexity are two important requirements of channel codes used in many applications, like machine-to-machine communications. In this paper, we show how these requirements can be fulfilled by using some special quasi-cyclic low-density parity-check block codes and spatially coupled low-density parity-check convolutional codes that we denote as compact. They are defined by parity-check matrices designed according to a recent approach based on sequentially multiplied columns. This method allows obtaining codes with girth up to 12. Many numerical examples of practical codes are provided.


title: Creating Fair Models of Atherosclerotic Cardiovascular Disease Risk

abstract: Guidelines for the management of atherosclerotic cardiovascular disease (ASCVD) recommend the use of risk stratification models to identify patients most likely to benefit from cholesterol-lowering and other therapies. These models have differential performance across race and gender groups with inconsistent behavior across studies, potentially resulting in an inequitable distribution of beneficial therapy. In this work, we leverage adversarial learning and a large observational cohort extracted from electronic health records (EHRs) to develop a "fair" ASCVD risk prediction model with reduced variability in error rates across groups. We empirically demonstrate that our approach is capable of aligning the distribution of risk predictions conditioned on the outcome across several groups simultaneously for models built from high-dimensional EHR data. We also discuss the relevance of these results in the context of the empirical trade-off between fairness and model performance.


title: The Variable Pulsar Wind Nebula of PSR J1809-1917

abstract: PSR J1809-1917 is a young ($\tau=51$ kyr) energetic ($\dot{E}=1.8\times10^{36}$ erg s$^{-1}$) radio pulsar powering a pulsar wind nebula (PWN). We report on the results of three Chandra X-ray Observatory observations which show that the PWN consists of a small ($\sim 20''$) bright compact nebula (CN) and faint extended emission seen up to $2'$ from the pulsar. The CN is elongated in the northeast-southwest direction and exhibits morphological and flux variability on a timescale of a few months. We also find evidence of small arcsecond-scale jets extending from the pulsar along the same direction, and exhibiting a hard power-law (PL) spectrum with photon index $\Gamma_{\rm jet}=1.2\pm0.1$. The more extended emission and CN share the same symmetry axis, which is also aligned with the direction toward the TeV $\gamma$-ray source HESS J1809--193, supporting their association. The spectrum of the extended nebula (EN) fits an absorbed PL with about the same slope as that of the CN, $\Gamma_{\rm CN}\approx\Gamma_{\rm EN}=1.55\pm0.09$; no spectral changes across the EN's 2 pc extent are seen. The total PWN 0.5-8 keV luminosity is $L_{\rm PWN}\approx 9\times10^{32}$ erg s$^{-1}$, about half of which is due to the EN.


title: Extragalactic Searches for Dark Matter Annihilation

abstract: We are at the dawn of a data-driven era in astrophysics and cosmology. A large number of ongoing and forthcoming experiments combined with an increasingly open approach to data availability offer great potential in unlocking some of the deepest mysteries of the Universe. Among these is understanding the nature of dark matter (DM)---one of the major unsolved problems in particle physics. Characterizing DM through its astrophysical signatures will require a robust understanding of its distribution in the sky and the use of novel statistical methods. The first part of this thesis describes the implementation of a novel statistical technique which leverages the "clumpiness" of photons originating from point sources (PSs) to derive the properties of PS populations hidden in astrophysical datasets. This is applied to data from the Fermi satellite at high latitudes ($|b| > 30$\deg) to characterize the contribution of PSs of extragalactic origin. We find that the majority of extragalactic gamma-ray emission can be ascribed to unresolved PSs having properties consistent with known sources such as active galactic nuclei. This leaves considerably less room for significant dark matter contribution. The second part of this thesis poses the question: "what is the best way to look for annihilating dark matter in extragalactic sources?" and attempts to answer it by constructing a pipeline to robustly map out the distribution of dark matter outside the Milky Way using galaxy group catalogs. This framework is then applied to Fermi data and existing group catalogs to search for annihilating dark matter in extragalactic galaxies and clusters.


title: Hausdorff dimension of Furstenberg-type sets associated to families of affine subspaces

abstract: We show that if $B \subset \mathbb{R}^n$ and $E \subset A(n,k)$ is a nonempty collection of $k$-dimensional affine subspaces of $\mathbb{R}^n$ such that every $P \in E$ intersects $B$ in a set of Hausdorff dimension at least $\alpha$ with $k-1 < \alpha \leq k$, then $\dim B \geq \alpha +\dim E/(k+1)$, where $\dim$ denotes the Hausdorff dimension. This estimate generalizes the well known Furstenberg-type estimate that every $\alpha$-Furstenberg set in the plane has Hausdorff dimension at least $\alpha + 1/2$. More generally, we prove that if $B$ and $E$ are as above with $0 < \alpha \leq k$, then $\dim B \geq \alpha +(\dim E-(k-\lceil \alpha \rceil)(n-k))/(\lceil \alpha \rceil+1)$. We also show that this bound is sharp for some parameters. As a consequence, we prove that for any $1 \leq k<n$, the union of any nonempty $s$-Hausdorff dimensional family of $k$-dimensional affine subspaces of $\mathbb{R}^n$ has Hausdorff dimension at least $k+\frac{s}{k+1}$.


title: Lifetime renormalization of weakly anharmonic superconducting qubits: I. Role of number non-conserving terms

abstract: The dynamics of a weakly anharmonic superconducting qubit in a complex electromagnetic environment is generally well-described by an effective multimode Kerr Hamiltonian at sufficiently weak excitation. This Hamiltonian can be embedded in a master equation with losses determined by the details of the electromagnetic environment. Recent experiments indicate, however, that when a superconducting circuit is driven with microwave signals the observed relaxation rates appear to be substantially different from expectations based on the electromagnetic environment of the qubit alone. This issue is a limiting factor in the optimization of superconducting qubit readout schemes. We claim here that an effective master equation with drive-power dependent parameters is an efficient approach to model such quantum dynamics. In this sequence of papers, we derive effective master equations, whose parameters exhibit nonlinear dependence on the excitation level of the circuit as well as the electromagnetic environment of the qubit. We show that the number non-conserving terms in the qubit nonlinearity generally lead to a renormalization of dissipative parameters of the effective master equation, while the number conserving terms give rise to a renormalization of the system frequencies. Here, in Part I, we consider the excitation-relaxation dynamics of a transmon qubit that is prepared in a certain initial state, but is not driven otherwise. A unitary transformation technique is introduced to study the renormalization of i) qubit relaxation due to coupling to a generic bath and ii) Purcell decay. Analytic expressions are provided for the dependence of the nonlinear dissipative terms on the details of the electromagnetic environment of the qubit. The perturbation technique based on unitary transformations developed here is generalized to the continuously driven case in Part II.


title: PARyOpt: A software for Parallel Asynchronous Remote Bayesian Optimization

abstract: PARyOpt is a python based implementation of the Bayesian optimization routine designed for remote and asynchronous function evaluations. Bayesian optimization is especially attractive for computational optimization due to its low cost function footprint as well as the ability to account for uncertainties in data. A key challenge to efficiently deploy any optimization strategy on distributed computing systems is the synchronization step, where data from multiple function calls is assimilated to identify the next campaign of function calls. Bayesian optimization provides an elegant approach to overcome this issue via asynchronous updates. We formulate, develop and implement a parallel, asynchronous variant of Bayesian optimization. The framework is robust and resilient to external failures. We show how such asynchronous evaluations help reduce the total optimization wall clock time for a suite of test problems. Additionally, we show how the software design of the framework allows easy extension to response surface reconstruction (Kriging), providing a high performance software for autonomous exploration. The software is available on PyPI, with examples and documentation.


title: Prediction and estimation consistency of sparse multi-class penalized optimal scoring

abstract: Sparse linear discriminant analysis via penalized optimal scoring is a successful tool for classification in high-dimensional settings. While the variable selection consistency of sparse optimal scoring has been established, the corresponding prediction and estimation consistency results have been lacking. We bridge this gap by providing probabilistic bounds on out-of-sample prediction error and estimation error of multi-class penalized optimal scoring allowing for diverging number of classes.


title: Undecidability of $\mathbb Q^{(2)}$

abstract: It is shown that the compositum $ \mathbb Q^{(2)}$ of all degree 2 extensions of $\mathbb Q$ has undecidable theory.


title: Higher Derivative Corrections to Lower Order RG Flow Equations

abstract: We show that the RG flow equation for the cosmological constant (CC) receives contributions (in addition to those coming from the CC the Einstein-Hilbert term and $R^{2}$ and $R_{\mu\nu}^{2}$ terms) only from terms with just two powers of curvature, but having also powers of the covariant derivative, in the Wilsonian effective action. In pure gravity our argument implies that just considering $f(R)$ theories will miss this effect which arises from terms such as $"R"\square^{n}"R",\,n=0,1,2,\ldots$. We expect similar contributions for the flow equation of the Einstein-Hilbert term as well. Finally we argue that the perturbative ghosts coming from curvature squared terms in the action are in fact spurious since they are at the cutoff scale and can be removed by (cutoff dependent) field redefinitions.


title: The twisted cohomological equation over the partially hyperbolic flow

abstract: Let $\mathbb{G}$ be a higher-rank connected semisimple Lie group with finite center and without compact factors. In any unitary representation $(\pi, \mathcal{H})$ of $\mathbb{G}$ without non-trivial $\mathbb{G}$-fixed vectors, we study the twisted cohomological equation $(X+m)f=g$, where $m\in\mathbb{R}$ and $X$ is in a $\mathbb{R}$-split Cartan subalgebra of $\text{Lie}(\mathbb{G})$. We characterize the obstructions to solving the cohomological equation, construct smooth solutions of the cohomological equation and obtain tame Sobolev estimates for $f$. We also study common solution to (the infinitesimal version of) the twisted cocycle equation $(X+m)g_1=(\mathfrak{v}+m_1)g_2$, where $\mathfrak{v}$ is nilpotent or in a $\mathbb{R}$-split Cartan subalgebra, $m,m_1\in\mathbb{R}$. This is the first paper studying general twisted equations. Compared to former papers, a new technique in representation theory is developed by Mackey theory and Mellin transform.


title: A Unified Batch Online Learning Framework for Click Prediction

abstract: We present a unified framework for Batch Online Learning (OL) for Click Prediction in Search Advertisement. Machine Learning models once deployed, show non-trivial accuracy and calibration degradation over time due to model staleness. It is therefore necessary to regularly update models, and do so automatically. This paper presents two paradigms of Batch Online Learning, one which incrementally updates the model parameters via an early stopping mechanism, and another which does so through a proximal regularization. We argue how both these schemes naturally trade-off between old and new data. We then theoretically and empirically show that these two seemingly different schemes are closely related. Through extensive experiments, we demonstrate the utility of of our OL framework; how the two OL schemes relate to each other and how they trade-off between the new and historical data. We then compare batch OL to full model retrains, and show how online learning is more robust to data issues. We also demonstrate the long term impact of Online Learning, the role of the initial Models in OL, the impact of delays in the update, and finally conclude with some implementation details and challenges in deploying a real world online learning system in production. While this paper mostly focuses on application of click prediction for search advertisement, we hope that the lessons learned here can be carried over to other problem domains.


title: Boundary value problems in Lipschitz domains for equations with lower order coefficients

abstract: We use the method of layer potentials to study the $R_2$ Regularity problem and the $D_2$ Dirichlet problem for second order elliptic equations of the form $\mathcal{L}u=0$, with lower order coefficients, in bounded Lipschitz domains. For $R_2$ we establish existence and uniqueness assuming that $\mathcal{L}$ is of the form $\mathcal{L}u=-\text{div}(A\nabla u+bu)+c\nabla u+du$, where the matrix $A$ is uniformly elliptic and H\"older continuous, $b$ is H\"older continuous, and $c,d$ belong to Lebesgue classes and they satisfy either the condition $d\geq\text{div}b$, or $d\geq\text{div}c$ in the sense of distributions. In particular, $A$ is not assumed to be symmetric, and there is no smallness assumption on the norms of the lower order coefficients. We also show existence and uniqueness for $D_2$ for the adjoint equations $\mathcal{L}^tu=0$.


title: The Simple Chromatic Number of $(m,n)$-Mixed Graphs

abstract: An $(m,n)$-mixed graph generalizes the notions of oriented graphs and edge-coloured graphs to a graph object with $m$ arc types and $n$ edge types. A simple colouring of such a graph is a non-trivial homomorphism to a reflexive target. We find that simple chromatic number of complete $(m,n)$-mixed graphs can be found in polynomial time. For planar graphs and $k$-trees ($k \geq 3$) we find that allowing the target to be reflexive does not lower the chromatic number of the respective family of $(m,n)$-mixed graphs. This implies that the search for universal targets for such families may be restricted to simple cliques.


title: Improved Basic Block Reordering

abstract: Basic block reordering is an important step for profile-guided binary optimization. The state-of-the-art goal for basic block reordering is to maximize the number of fall-through branches. However, we demonstrate that such orderings may impose suboptimal performance on instruction and I-TLB caches. We propose a new algorithm that relies on a model combining the effects of fall-through and caching behavior. As details of modern processor caching is quite complex and often unknown, we show how to use machine learning in selecting parameters that best trade off different caching effects to maximize binary performance. An extensive evaluation on a variety of applications, including Facebook production workloads, the open-source compilers Clang and GCC, and SPEC CPU benchmarks, indicate that the new method outperforms existing block reordering techniques, improving the resulting performance of applications with large code size. We have open sourced the code of the new algorithm as a part of a post-link binary optimization tool, BOLT.


title: A scalable method to find the shortest path in a graph with circuits of memristors

abstract: Finding the shortest path in a graph has applications to a wide range of optimization problems. However, algorithmic methods scale with the size of the graph in terms of time and energy. We propose a method to solve the shortest path problem using circuits of nanodevices called memristors and validate it on graphs of different sizes and topologies. It is both valid for an experimentally derived memristor model and robust to device variability. The time and energy of the computation scale with the length of the shortest path rather than with the size of the graph, making this method particularly attractive for solving large graphs with small path lengths.


title: The First Glitch in a Central Compact Object Pulsar: 1E 1207.4-5209

abstract: Since its discovery as a pulsar in 2000, the central compact object (CCO) 1E 1207.4-5209 in the supernova remnant PKS 1209-51/52 had been a stable 0.424 s rotator with an extremely small spin-down rate and weak (Bs ~ 9E10 G) surface dipole magnetic field. In 2016 we observed a glitch from 1E 1207.4-5209 of at least Delta f/f = (2.8+/-0.4)E-9, which is typical in size for the general pulsar population. However, glitch activity is closely correlated with spin-down rate fdot, and pulsars with fdot as small as that of 1E 1207.4-5209 are never seen to glitch. Unlike in glitches of ordinary pulsars, there may have been a large increase in fdot as well. The thermal X-ray spectrum of 1E 1207.4-5209, with its unique cyclotron absorption lines that measure the surface magnetic field strength, did not show any measurable change after the glitch, which rules out a major disruption in the dipole field as a cause or result of the glitch. A leading theory of the origin and evolution of CCOs, involving prompt burial of the magnetic field by fall-back of supernova ejecta, might hold the explanation for the glitch.


title: Supervised machine learning of ultracold atoms with speckle disorder

abstract: We analyze how accurately supervised machine learning techniques can predict the lowest energy levels of one-dimensional noninteracting ultracold atoms subject to the correlated disorder due to an optical speckle field. Deep neural networks with different numbers of hidden layers and neurons per layer are trained on large sets of instances of the speckle field, whose energy levels have been preventively determined via a high-order finite difference technique. The Fourier components of the speckle field are used as feature vector to represent the speckle-field instances. A comprehensive analysis of the details that determine the possible success of supervised machine learning tasks, namely the depth and the width of the neural network, the size of the training set, and the magnitude of the regularization parameter, is presented. It is found that ground state energies of previously unseen instances can be predicted with essentially arbitrary accuracy. First and second excited state energies can be predicted too, albeit with slightly lower accuracy and using more layers of hidden neurons.


title: Production of photonic universal quantum gates enhanced by machine learning

abstract: We introduce photonic architectures for universal quantum computation. The first step is to produce a resource state which is a superposition of the first four Fock states with a probability $\geq 10^{-2}$, an increase by a factor of $10^4$ over standard sequential photon-subtraction techniques. The resource state is produced with near-perfect fidelity from a quantum gadget that uses displaced squeezed vacuum states, interferometers and photon-number resolving detectors. The parameters of this gadget are trained using machine learning algorithms for variational circuits. We discuss in detail various aspects of the non-Gaussian state preparation resulting from the numerical experiments. We then propose a notion of resource farms where these gadgets are stacked in parallel, to increase the success probability further. We find a trade-off between the success probability of the farm, the error tolerance, and the number of gadgets. Using the resource states in conventional gate teleportation techniques we can then implement weak tuneable cubic phase gates. The numerical tools that have been developed could potentially be useful for other applications in photonics as well.


title: Loosely-bound objects produced in nuclear collisions at the LHC

abstract: Loosely-bound objects such as light nuclei are copiously produced in proton-proton and nuclear collisions at the Large Hadron Collider (LHC), despite the fact that typical energy scales in such collisions exceed the binding energy of the objects by orders of magnitude. In this review we summarise the experimental observations, put them into context of previous studies at lower energies, and discuss the underlying physics. Most of the data discussed here were taken by the ALICE Collaboration during LHC Run1, which started in 2009 and ended in 2013. Specifically we focus on the production of (anti-)nuclei and (anti-)hypernuclei. Also included are searches for exotic objects like the H-dibaryon, a possible $uuddss$ hexaquark state, or also a possible bound state of a $\Lambda$ hyperon and a neutron. Furthermore, the study of hyperon-nucleon and hyperon-hyperon interactions through measurements of correlations are briefly discussed, especially in connection with the possible existence of loosely-bound states composed of these baryons. In addition, some results in the strange and charmed hadron sector are presented, to show the capabilities for future measurements on loosely-bound objects in this direction. Finally, perspectives are given for measurements in the currently ongoing Run2 period of the LHC and in the future LHC Run3.


title: Automatic Program Synthesis of Long Programs with a Learned Garbage Collector

abstract: We consider the problem of generating automatic code given sample input-output pairs. We train a neural network to map from the current state and the outputs to the program's next statement. The neural network optimizes multiple tasks concurrently: the next operation out of a set of high level commands, the operands of the next statement, and which variables can be dropped from memory. Using our method we are able to create programs that are more than twice as long as existing state-of-the-art solutions, while improving the success rate for comparable lengths, and cutting the run-time by two orders of magnitude. Our code, including an implementation of various literature baselines, is publicly available at https://github.com/amitz25/PCCoder


title: SAFE: A Neural Survival Analysis Model for Fraud Early Detection

abstract: Many online platforms have deployed anti-fraud systems to detect and prevent fraudulent activities. However, there is usually a gap between the time that a user commits a fraudulent action and the time that the user is suspended by the platform. How to detect fraudsters in time is a challenging problem. Most of the existing approaches adopt classifiers to predict fraudsters given their activity sequences along time. The main drawback of classification models is that the prediction results between consecutive timestamps are often inconsistent. In this paper, we propose a survival analysis based fraud early detection model, SAFE, which maps dynamic user activities to survival probabilities that are guaranteed to be monotonically decreasing along time. SAFE adopts recurrent neural network (RNN) to handle user activity sequences and directly outputs hazard values at each timestamp, and then, survival probability derived from hazard values is deployed to achieve consistent predictions. Because we only observe the user suspended time instead of the fraudulent activity time in the training data, we revise the loss function of the regular survival model to achieve fraud early detection. Experimental results on two real world datasets demonstrate that SAFE outperforms both the survival analysis model and recurrent neural network model alone as well as state-of-the-art fraud early detection approaches.


title: Fair lending needs explainable models for responsible recommendation

abstract: The financial services industry has unique explainability and fairness challenges arising from compliance and ethical considerations in credit decisioning. These challenges complicate the use of model machine learning and artificial intelligence methods in business decision processes.


title: Thermal Switching of Lasing Regimes in Heavily Doped Er3+ Fiber Lasers

abstract: A pulsed regime of short-cavity, heavily erbium-doped fiber lasers is of high interest for its possible applications in telecommunications and sensorics. Here, we demonstrate these lasers in two configurations, distributed feedback laser and compare it with a classic Fabry-Perot type laser. We have managed to create lasers that function stably with cavities as small as 50 mm. Pulse properties such as amplitude, frequency and duration, are in a good agreement with our theoretical analysis, which takes into account spontaneous emission. We report the observation of the thermal switching effect, which consists of the pulsing regime changing to CW upon cooling the laser cavities down to the liq-uid nitrogen temperature. We theoretically show that this effect may be explained by weakening of the up-conversion process responsible for the pulsed regime. The slowing of the up-conversion processes is due to the energy mismatch in this process, which is overcome by interaction with phonons. At low temperatures, the number of phonons decreas-es and pulsing switches off.


title: Zero-Shot Cross-lingual Classification Using Multilingual Neural Machine Translation

abstract: Transferring representations from large supervised tasks to downstream tasks has shown promising results in AI fields such as Computer Vision and Natural Language Processing (NLP). In parallel, the recent progress in Machine Translation (MT) has enabled one to train multilingual Neural MT (NMT) systems that can translate between multiple languages and are also capable of performing zero-shot translation. However, little attention has been paid to leveraging representations learned by a multilingual NMT system to enable zero-shot multilinguality in other NLP tasks. In this paper, we demonstrate a simple framework, a multilingual Encoder-Classifier, for cross-lingual transfer learning by reusing the encoder from a multilingual NMT system and stitching it with a task-specific classifier component. Our proposed model achieves significant improvements in the English setup on three benchmark tasks - Amazon Reviews, SST and SNLI. Further, our system can perform classification in a new language for which no classification data was seen during training, showing that zero-shot classification is possible and remarkably competitive. In order to understand the underlying factors contributing to this finding, we conducted a series of analyses on the effect of the shared vocabulary, the training data type for NMT, classifier complexity, encoder representation power, and model generalization on zero-shot performance. Our results provide strong evidence that the representations learned from multilingual NMT systems are widely applicable across languages and tasks.


title: Low-Temperature and Dynamic Magnetism of Highly Frustrated 5d2 Li4MgOsO6 Polymorphs in Comparison with 5d3 Li3Mg2OsO6

abstract: Geometric magnetic frustration (GMF) has attracted substantial interest due to the exotic physics and rich phase diagrams revealed by the cancellation of normally-dominant magnetic interactions, giving impetus for the search for novel frustrated systems, most often based on antiferromagnetic correlations between magnetic ions decorating triangular or tetrahedral lattices. We report here low-temperature magnetic susceptibility and muon spin relaxation results on Li4MgOsO6 and Li3Mg2OsO6, members of the A5BO6 "rock salt ordered" family of frustrated materials. In Li3Mg2OsO6 we find spin freezing below 12K. In Li4MgOsO6, which can crystallize into either orthorhombic Fddd or monoclinic C2/m crystal symmetries depending on synthesis conditions, we find magnetism consistent with glassy-like behavior dominating below 2K, with partial ordering and evidence for dynamics at somewhat higher temperatures.


title: Complex electric double layers in charged topological colloids

abstract: Charged surfaces in contact with liquids containing ions are accompanied in equilibrium by an electric double layer consisting of a layer of electric charge on the surface that is screened by a diffuse ion cloud in the bulk fluid. This screening cloud determines not only the interactions between charged colloidal particles or polyelectrolytes and their self-assembly into ordered structures, but it is also pivotal in understanding energy storage devices, such as electrochemical cells and supercapacitors. However, little is known to what spatial complexity the electric double layers can be designed. Here, we show that electric double layers of non-trivial topology and geometry -including tori, multi-tori and knots- can be realised in charged topological colloidal particles, using numerical modelling within a mean-field Poisson-Boltzmann theory. We show that the complexity of double layers -including geometry and topology- can be tuned by changing the Debye screening length of the medium, or by changing the shape and topology of the (colloidal) particle. More generally, this work is an attempt to introduce concepts of topology in the field of charged colloids, which could lead to novel exciting material design paradigms.

categories: cond-mat.soft cond-mat.stat-mech physics.bio-ph physics.chem-ph


title: Global and short-range entanglement properties in excited, many-body localized spin chains

abstract: We explore the use of short-range entanglement measures, such as concurrence and negativity, and global entanglement measures such as geometric entanglement, as indicators of many-body localization (MBL) in the spectra of disordered spin systems. From the perspective of entanglement monogamy, the two types of entanglement behave oppositely in the thermalized and MBL phases. In a recent work, the concurrence of subsystems, a measure of local entanglement, was used in a study of many-body localization in a one-dimensional spin-$1/2$ system (Bera and Lakshminarayan, 2016). We show numerically that the negativity displays notably similar behavior for this system, with the advantage that it can also be extended to systems of higher local dimension. We then demonstrate this extension in practice by using it to predict the existence of an MBL phase in a disordered a spin-1 system. In terms of global entanglement, the geometric entanglement of both spin-$1/2$ and spin-1 systems is also shown to behave as a complementary indicator of the MBL phenomenon.


title: Hyperplane Sections of Determinantal Varieties over Finite Fields and Linear Codes

abstract: We determine the number of ${\mathbb{F}}_q$-rational points of hyperplane sections of classical determinantal varieties defined by the vanishing of minors of a fixed size of a generic matrix, and identify sections giving the maximum number of ${\mathbb{F}}_q$-rational points. Further we consider similar questions for sections by linear subvarieties of a fixed codimension in the ambient projective space. This is closely related to the study of linear codes associated to determinantal varieties, and the determination of their weight distribution, minimum distance and generalized Hamming weights. The previously known results about these are generalized and expanded significantly