title: Resilience of the topological phases to frustration

abstract: Recently it was highlighted that one-dimensional antiferromagnetic spin models with frustrated boundary conditions, i.e. periodic boundary conditions in a ring with an odd number of elements, may show very peculiar behavior. Indeed the presence of frustrated boundary conditions can destroy the local magnetic orders presented by the models when different boundary conditions are taken into account and induce novel phase transitions. Motivated by these results, we analyze the effects of the introduction of frustrated boundary conditions on several models supporting (symmetry protected) topological orders, and compare our results with the ones obtained with different boundary conditions. None of the topological order phases analyzed are altered by this change. This observation leads naturally to the conjecture that topological phases of one-dimensional systems are in general not affected by topological frustration.

categories: cond-mat.stat-mech cond-mat.str-el hep-th quant-ph


title: Rigidification of connective comodules

abstract: Let $\mathbb{k}$ be a commutative ring with global dimension zero. We show that we can rigidify homotopy coherent comodules in connective modules over the Eilenberg-Mac Lane spectrum of $\mathbb{k}$. That is, the $\infty$-category of homotopy coherent comodules is represented by a model category of strict comodules in non-negative chain complexes over $\mathbb{k}$. These comodules are over a coalgebra that is strictly coassociative and simply connected. The rigidification result allows us to derive the notion of cotensor product of comodules and endows the $\infty$-category of comodules with a symmetric monoidal structure via the two-sided cobar resolution.


title: Onset of spontaneous scalarization in spinning Gauss-Bonnet black holes

abstract: It has recently been proved numerically that spinning black holes in Einstein-scalar theories which are characterized by a non-minimal negative coupling of the scalar field to the Gauss-Bonnet invariant of the curved spacetime may develop exponentially growing instabilities. Intriguingly, it has been demonstrated that this tachyonic instability, which marks the onset of the spontaneous scalarization phenomenon in the Einstein-Gauss-Bonnet-scalar theory, characterizes spinning black holes whose dimensionless angular momentum parameter ${\bar a}\equiv a/M$ is larger than some critical value ${\bar a}_{\text{crit}}\simeq0.505$. In the present paper we prove, using {\it analytical} techniques, that the critical rotation parameter which marks the boundary between bald Kerr black holes and hairy (scalarized) spinning black holes in the Einstein-Gauss-Bonnet-scalar theory is given by the exact dimensionless relation ${\bar a}_{\text{crit}}={1\over 2}$.


title: Spontaneous oscillations and negative-conductance transitions in microfluidic networks

abstract: The tendency for flows in microfluidic systems to behave linearly poses a challenge for designing integrated flow control schemes to carry out complex fluid processing tasks. This hindrance has led to the use of numerous external control devices to manipulate flows, thereby thwarting the potential scalability and portability of lab-on-a-chip technology. Here, we devise a microfluidic network exhibiting nonlinear flow dynamics that enable new mechanisms for on-chip flow control. This network is shown to exhibit oscillatory output patterns, bistable flow states, hysteresis, signal amplification, and negative-conductance transitions, all without reliance on dedicated external control hardware, movable parts, flexible components, or oscillatory inputs. These dynamics arise from nonlinear fluid inertia effects in laminar flows that we amplify and harness through the design of the network geometry. We suggest that these results, which are supported by fluid dynamical simulations and theoretical modeling, have the potential to inspire development of new built-in control capabilities, such as on-chip timing and synchronized flow patterns.


title: Phase diagram and compact stars in a holographic QCD model

abstract: A holographic model is used to investigate the thermodynamics and the phase diagram of a heavy quarks system. From such a model we obtain an equation of state and explore its applicability in astrophysical conditions. For this objective, we work in the context of the Einstein-Maxwell-Dilaton (EMD) holographic model for quantum chromodynamics (QCD). At first, we show the existence of a critical point where the first-order transitions line ends, later on, we calculated an analytic expression for the equation of state. Additionally, with the aim of investigating the global properties of compact stars, such as the total gravitational mass and radius, the equation of state is used to solve the Tolman-Oppenheimer-Volkov (TOV) equations for stellar structure. The numerical results show that our equation of state is able to reproduce the expected behavior of hybrid stars. Our main conclusion is that, by using an equation of state emerging in the framework of the EMD holographic model for QCD, it is possible to obtain quark matter properties and that it is also possible to extend the procedure to astrophysical applications.


title: Missing [CII] emission from early galaxies

abstract: ALMA observations have revealed that [CII] 158$\mu$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z>6 galaxies observed in both [CII] and [OIII]. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect [CII] emission in the whole sample, with an extent systematically larger than the [CII] emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40% of the extended [CII] component might be missed at an angular resolution of 0.8$^{\prime\prime}$, implying that $L_{\rm [CII]}$ is underestimated by a factor $\approx2$ in data at low (<7) signal-to-noise ratio . By combining these results, we conclude that $L_{\rm [CII]}$ of z>6 galaxies lies, on average, slightly below the local $L_{\rm [CII]}-SFR$ relation ($\Delta^{z=6-9}=-0.07\pm0.3$), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm [OIII]}/L_{\rm [CII]}<10$, i.e. more in line with current hydrodynamical simulations.


title: GJ 273: On the formation, dynamical evolution and habitability of a planetary system hosted by an M dwarf at 3.75 parsec

abstract: Planets orbiting low-mass stars such as M dwarfs are now considered a cornerstone in the search for life-harbouring planets. GJ273 is a planetary system orbiting an M dwarf only 3.75 pc away, composed of two confirmed planets, GJ273b and GJ273c, and two promising candidates, GJ273d and GJ273e. Planet GJ273b resides in the habitable zone. Currently, due to a lack of observed planetary transits, only the minimum masses of the planets are known.Despite being an interesting system, the GJ273 planetary system is still poorly studied. We aim at precisely determine the physical parameters of the individual planets, in particular to break the mass--inclination degeneracy to accurately determine the mass of the planets. Moreover, we present thorough characterisation of planet GJ273b in terms of its potential habitability. We explored the planetary formation and hydration phases of GJ273 during the first 100 Myr. Then, we analysed the stability of the system. We also searched for regions which may harbour minor bodies such as an asteroid belt and Kuiper belt analogues. We found that the four-planet configuration of the system allows us to break the mass--inclination degeneracy with the following masses: $2.89\leq M_{\mathrm{b}}\leq3.03\,\mathrm{M}_\oplus$, $1.18\leq M_{\mathrm{c}}\leq1.24\,\mathrm{M}_\oplus$, $10.80\leq M_{\mathrm{d}}\leq11.35\,\mathrm{M}_\oplus$, and $9.30\leq M_{\mathrm{e}}\leq9.70\,\mathrm{M}_\oplus$. That is an Earth-mass planet, a super-Earth and two mini-Neptunes. Moreover, GJ273b is found to be an efficient water captor and GJ273c likely a dry planet. Several stable regions are predicted where minor bodies might reside. We comprehensively discuss the habitability of GJ273b.


title: First evidence of inertial modes in $\gamma$ Doradus stars: The core rotation revealed

abstract: Gamma Doradus stars present an incredibly rich pulsation spectra, with gravito-inertial modes, in some cases supplemented with delta Scuti-like pressure modes and in numerous cases with Rossby modes. The present paper aims at showing that, in addition to these modes established in the radiative envelope, pure inertial modes, trapped in the convective core, can be detected in Kepler observations of gamma Doradus stars, thanks to their resonance with the gravito-inertial modes. We start by using a simplified model of perturbations in a full sphere of uniform density. Under these conditions, the spectrum of pure inertial modes is known from analytical solutions of the so-called Poincare equation. We then compute coupling factors which help select the pure inertial modes which interact best with the surrounding dipolar gravito-inertial modes. Using complete calculations of gravito-inertial modes in realistic models of gamma Doradus stars, we are able to show that the pure inertial/gravito-inertial resonances appear as dips in the gravito-inertial mode period spacing series at spin parameters close to those predicted by the simple model. We find the first evidence of such dips in the Kepler gamma Doradus star KIC5608334. Finally, using complete calculations in isolated convective cores, we find that the spin parameters of the pure inertial/gravito-inertial resonances are also sensitive to the density stratification of the convective core. In conclusion, we have discovered that certain dips in gravito-inertial mode period spacings observed in some Kepler stars are in fact the signatures of resonances with pure-inertial modes that are trapped in the convective core. This holds the promise to finally access the central conditions , i.e. rotation and density stratification, of intermediate-mass stars on the main sequence.


title: Non-adiabatic molecular quantum dynamics with quantum computers

abstract: The theoretical investigation of non-adiabatic processes is hampered by the complexity of the coupled electron-nuclear dynamics beyond the Born-Oppenheimer approximation. Classically, the simulation of such reactions is limited by the unfavourable scaling of the computational resources as a function of the system size. While quantum computing exhibits proven quantum advantage for the simulation of real-time dynamics, the study of quantum algorithms for the description of non-adiabatic phenomena is still unexplored. In this work, we propose a quantum algorithm for the simulation of fast non-adiabatic chemical processes together with an initialization scheme for quantum hardware calculations. In particular, we introduce a first-quantization method for the time evolution of a wavepacket on two coupled harmonic potential energy surfaces (Marcus model). In our approach, the computational resources scale polynomially in the system dimensions, opening up new avenues for the study of photophysical processes that are classically intractable.


title: ALMA Observations of Multiple-CO and C Lines Toward the Active Galactic Nucleus of NGC 7469: X-Ray-dominated Region Caught in the Act

abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $\alpha_{\rm CO} = 4.1$ and $\alpha_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments.


title: Coalgebras in the Dwyer-Kan localization of a model category

abstract: We show that weak monoidal Quillen equivalences induce equivalences of symmetric monoidal $\infty$-categories with respect to the Dwyer-Kan localization of the symmetric monoidal model categories. The result will induce a Dold-Kan correspondence of coalgebras in $\infty$-categories. Moreover, it shows that Shipley's zig-zag of Quillen equivalences lifts to an explicit symmetric monoidal equivalence of $\infty$-categories for the stable Dold-Kan correspondence. We study homotopy coherent coalgebras associated to a monoidal monoidal category. We show examples when these coalgebras cannot be rigidified. That is, their $\infty$-categories are not equivalent to the Dwyer-Kan localizations of strict coalgebras in the usual monoidal model categories of spectra and of connective discrete $R$-modules.


title: The coalgebraic enrichment of algebras in higher categories

abstract: We prove that given $\mathcal{C}$ a presentably symmetric monoidal $\infty$-category, and any essentially small $\infty$-operad $\mathcal{O}$, the $\infty$-category of $\mathcal{O}$-algebras in $\mathcal{C}$ is enriched, tensored and cotensored over the presentably symmetric monoidal $\infty$-category of $\mathcal{O}$-coalgebras in $\mathcal{C}$. We provide a higher categorical analogue of the universal measuring coalgebra. For categories in the usual sense, the result was proved by Hyland, L\'{o}pez Franco, and Vasilakopoulou.


title: Betelgeuse fainter in the sub-millimetre too: an analysis of JCMT and APEX monitoring during the recent optical minimum

abstract: Betelgeuse is the nearest Red Supergiant star and it underwent an unusually deep minimum at optical wavelengths during its most recent pulsation cycle. We present submillimetre observations taken by the James Clerk Maxwell Telescope and Atacama Pathfinder Experiment over a time span of 13 years including the optical dimming. We find that Betelgeuse has also dimmed by \sim20\% at these longer wavelengths during this optical minimum. Using radiative-transfer models, we show that this is likely due to changes in the photosphere (luminosity) of the star as opposed to the surrounding dust as was previously suggested in the literature.


title: Fast Correlated-Photon Imaging Enhanced by Deep Learning

abstract: Correlated photon pairs, carrying strong quantum correlations, have been harnessed to bring quantum advantages to various fields from biological imaging to range finding. Such inherent non-classical properties support extracting more valid signals to build photon-limited images even in low flux-level, where the shot noise becomes dominant as light source decreases to single-photon level. Optimization by numerical reconstruction algorithms is possible but require thousands of photon-sparse frames, thus unavailable in real time. Here, we present an experimental fast correlated-photon imaging enhanced by deep learning, showing an intelligent computational strategy to discover deeper structure in big data. Convolutional neural network is found being able to efficiently solve image inverse problems associated with strong shot noise and background noise (electronic noise, scattered light). Our results fill the key gap in incompatibility between imaging speed and image quality by pushing low-light imaging technique to the regime of real-time and single-photon level, opening up an avenue to deep leaning-enhanced quantum imaging for real-life applications.

categories: quant-ph eess.IV physics.app-ph physics.optics


title: Gamma-ray-emitting narrow-line Seyfert 1 galaxies: the Swift view

abstract: We report the analysis of all Swift observations available up to 2019 April of $\gamma$-ray-emitting narrow-line Seyfert 1 galaxies (NLSy1). The distribution of X-ray luminosities (and fluxes) indicates that the jet radiation significantly contributes to their X-ray emission, with Doppler boosting making values higher than other radio-loud NLSy1. The 0.3-10 keV photon indices are on average harder with respect to radio-quiet and radio-loud NLSy1, confirming a dominant jet contribution in X-rays. However, the lower variability amplitude with respect to blazars and the softening of the spectrum in some periods suggests that also the corona radiation contributes to the X-ray emission. In optical and ultraviolet (UV) significant flux changes have been observed on daily, weekly, and monthly time-scale, providing a clear indication of the significant contribution of the jet radiation in this part of spectrum. A strong correlation between X-ray, UV, and optical emission and simultaneous flux variations have been observed in 1H 0323+342, SBS 0846+513, PMN J0948+0022 as expected in case the jet radiation is the dominant mechanism. Correlated multiband variability favours the jet-dominated scenario also in FBQS J1644+2619 and PKS 2004-447. The summed X-ray Telescope spectra of 1H 0323+342, SBS 0846+513, PMN J0948+0022, and FBQS J1644+2619 are well fitted by a broken power law with a break around 2 keV. The spectrum above 2 keV is dominated by the non-thermal emission from a beamed relativistic jet, as suggested by the hard photon index. A Seyfert-like feature like the soft X-ray excess has been observed below 2 keV, making these $\gamma$-ray-emitting NLSy1 different from typical blazars.


title: Nonlinear two-photon Rabi-Hubbard model: superradiance and photon/photon-pair Bose-Einstein condensate

abstract: We study the ground state phase diagram of a nonlinear two-photon Rabi-Hubbard (RH) model in one dimension using quantum Monte Carlo (QMC) simulations and density matrix renormalization group (DMRG) calculations. Our model includes a nonlinear photon-photon interaction term. Absent this term, the RH model has only one phase, the normal disordered phase, and suffers from spectral collapse at larger values of the photon-qubit interaction or inter-cavity photon hopping. The photon-photon interaction, no matter how small, stabilizes the system which now exhibits {\it two} quantum phase transitions: Normal phase to {\it photon pair} superfluid (PSF) transition and PSF to single particle superfluid (SPSF). The discrete $Z_4$ symmetry of the Hamiltonian spontaneously breaks in two stages: First it breaks partially as the system enters the PSF and then completely breaks when the system finally enters the SPSF phase. We show detailed numerical results supporting this, and map out the ground state phase diagram.


title: A double integral of dlog forms which is not polylogarithmic

abstract: Feynman integrals are central to all calculations in perturbative Quantum Field Theory. They often give rise to iterated integrals of dlog-forms with algebraic arguments, which in many cases can be evaluated in terms of multiple polylogarithms. This has led to certain folklore beliefs in the community stating that all such integrals evaluate to polylogarithms. Here we discuss a concrete example of a double iterated integral of two dlog-forms that evaluates to a period of a cusp form. The motivic versions of these integrals are shown to be algebraically independent from all multiple polylogarithms evaluated at algebraic arguments. From a mathematical perspective, we study a mixed elliptic Hodge structure arising from a simple geometric configuration in $\mathbb{P}^2$, consisting of a modular plane elliptic curve and a set of lines which meet it at torsion points, which may provide an interesting worked example from the point of view of periods, extensions of motives, and L-functions.


title: Gravitational Collapse in Cubic Horndeski Theories

abstract: We study spherically symmetric gravitational collapse in cubic Horndeski theories of gravity. By varying the coupling constants and the initial amplitude of the scalar field, we determine the region in the space of couplings and amplitudes for which it is possible to construct global solutions to the Horndeski theories. Furthermore, we identify the regime of validity of effective field theory as the sub-region for which a certain weak field condition remains small at all times. We evolve the initial data using the CCZ4 formulation of the Einstein equations and horizon penetrating coordinates without assuming spherical symmetry.


title: Accelerated variational algorithms for digital quantum simulation of many-body ground states

abstract: One of the key applications for the emerging quantum simulators is to emulate the ground state of many-body systems, as it is of great interest in various fields from condensed matter physics to material science. Traditionally, in an analog sense, adiabatic evolution has been proposed to slowly evolve a simple Hamiltonian, initialized in its ground state, to the Hamiltonian of interest such that the final state becomes the desired ground state. Recently, variational methods have also been proposed and realized in quantum simulators for emulating the ground state of many-body systems. Here, we first provide a quantitative comparison between the adiabatic and variational methods with respect to required quantum resources on digital quantum simulators, namely the depth of the circuit and the number of two-qubit quantum gates. Our results show that the variational methods are less demanding with respect to these resources. However, they need to be hybridized with a classical optimization which can converge slowly. Therefore, as the second result of the paper, we provide two different approaches for speeding the convergence of the classical optimizer by taking a good initial guess for the parameters of the variational circuit. We show that these approaches are applicable to a wide range of Hamiltonian and provide significant improvement in the optimization procedure.


title: Empirical completeness assessment of the Gaia DR2, Pan-STARRS 1 and ASAS-SN-II RR Lyrae catalogues

abstract: RR Lyrae stars are an important and widely used tracer of the most ancient populations of our Galaxy, mainly due to their standard candle nature. The availability of large scale surveys of variable stars is allowing us to trace the structure of our entire Galaxy, even in previously inaccessible areas like the Galactic disc. In this work we aim to provide an empirical assessment of the completeness of the three largest RR Lyrae catalogues available: Gaia DR2, PanSTARRS-1 and ASAS-SN-II. Using a joint probabilistic analysis of the three surveys we compute 2D and 3D completeness maps in each survey's full magnitude range. At the bright end (G<13) we find ASAS-SN-II and Gaia are near 100% complete in RRab at high latitude (|b|>20deg); ASAS-SN-II has the best completeness at low latitude for RRab and at all latitudes for RRc. At the faint end (G>13), Gaia DR2 is the most complete catalogue for both RR Lyrae types, at any latitude, with median completeness rates of 95% (RRab) and >85% (RRc) outside the ecliptic plane (|\beta|>25deg). We confirm a high and uniform completeness of PanSTARRS-1 RR Lyrae at 91% (RRab) and 82% (RRc) down to G~18, and provide the first estimate of its completeness at low galactic latitude (|b|<20deg) at an estimated median 65% (RRab) and 50-60% (RRc). Our results are publicly available as 2D and 3D completeness maps, and as functions to evaluate each survey's completeness versus distance or per line-of sight.


title: WD1032+011, an inflated brown dwarf in an old eclipsing binary with a white dwarf

abstract: We present the discovery of only the third brown dwarf known to eclipse a non-accreting white dwarf. Gaia parallax information and multi-colour photometry confirm that the white dwarf is cool (9950$\pm$150K) and has a low mass (0.45$\pm$0.05~MSun), and spectra and lightcurves suggest the brown dwarf has a mass of 0.067 $\pm$0.006 MSun (70 MJup) and a spectral type of L5 $\pm$1. The kinematics of the system show that the binary is likely to be a member of the thick disk and therefore at least 5 Gyr old. The high cadence lightcurves show that the brown dwarf is inflated, making it the first brown dwarf in an eclipsing white dwarf-brown dwarf binary to be so.


title: QED driven QAOA for network-flow optimization

abstract: We present a general framework for modifying quantum approximate optimization algorithms (QAOA) to solve constrained network flow problems. By exploiting an analogy between flow constraints and Gauss's law for electromagnetism, we design lattice quantum electrodynamics (QED) inspired mixing Hamiltonians that preserve flow constraints throughout the QAOA process. This results in an exponential reduction in the size of the configuration space that needs to be explored, which we show through numerical simulations, yields higher quality approximate solutions compared to the original QAOA routine. We outline a specific implementation for edge-disjoint path (EDP) problems related to traffic congestion minimization, numerically analyze the effect of initial state choice, and explore trade-offs between circuit complexity and qubit resources via a particle-vortex duality mapping. Comparing the effect of initial states reveals that starting with an ergodic (unbiased) superposition of solutions yields better performance than beginning with the mixer ground-state, suggesting a departure from the "short-cut to adiabaticity" mechanism often used to motivate QAOA.


title: New production channels for light dark matter in hadronic showers

abstract: Hadronic showers transfer a relevant amount of their energy to electromagnetic subshowers. We show that the generation of "secondary" dark photons in these sub-showers is significant and typically dominates the production at low dark photon masses. The resulting dark photons are however substantially less energetic than the ones originating from mesons decay. We illustrate this point both semi-analytically and through Monte Carlo simulations. Existing limits on vector-mediator scenarios for light dark matter are updated with the inclusion of the new production processes.


title: Correlations between Triaxial Shapes and Formation History of Dark Matter Haloes

abstract: The shape of dark matter haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in dark matter haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the peak-centroid offset. In particular, the correlation between halo ellipticity and halo concentration is nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multi-dimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-size halos. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters.


title: An averaged cluster approach to including chemical short range order in KKR-CPA

abstract: The single-site Korringa-Kohn-Rostoker Coherent Potential Approximation (KKR-CPA) ignores short range ordering present in disordered metallic systems. In this paper, we establish a new technique to fix this shortcoming by embedding an averaged cluster that displays chemical short range order (SRO). The degree of SRO can be tuned by externally defined order parameters. This averaged cluster can be embedded in the single site CPA medium, or a self-consistently obtained effective medium that contains SRO information. The validity of this method is demonstrated by applying it to two alloy systems - the CuZn body centered cubic (BCC) solid solution, and AlCrTiV, a four-element BCC high entropy alloy. A comparison between the non-self-consistent and self-consistent modes is also provided for the two above mentioned systems. We make the code available on the internet. Planned extensions to this work are discussed.


title: Non-bipartite k-common graphs

abstract: A graph H is k-common if the number of monochromatic copies of H in a k-edge-coloring of K_n is asymptotically minimized by a random coloring. For every k, we construct a connected non-bipartite k-common graph. This resolves a problem raised by Jagger, Stovicek and Thomason [Combinatorica 16 (1996), 123-141]. We also show that a graph H is k-common for every k if and only if H is Sidorenko and that H is locally k-common for every k if and only if H is locally Sidorenko.


title: The Gaia-ESO Survey: Membership probabilities for stars in 32 open clusters from 3D kinematics

abstract: The Gaia-ESO Survey (GES) observed many open clusters as part of its programme to spectroscopically characterise the various Milky Way populations. GES spectroscopy and Gaia astrometry from its second data release are used here to assign membership probabilities to targets towards 32 open clusters with ages from 1-3800 Myr, based on maximum likelihood modelling of the 3D kinematics of the cluster and field populations. From a parent catalogue of 14398 individual targets, 5033 stars with uniformly determined 3D velocities, $T_{\rm eff}$, $\log g$ and chemistry are assigned cluster membership with probability $>0.9$, and with an average probability of 0.991. The robustness of the membership probabilities is demonstrated using independent membership criteria (lithium and parallax) in two of the youngest clusters. The addition of radial velocities improves membership discrimination over proper motion selection alone, especially in more distant clusters. The kinematically-selected nature of the membership lists, independent of photometry and chemistry, makes the catalogue a valuable resource for testing stellar evolutionary models and investigating the time evolution of various parameters.


title: Non-equilibrium dynamics of the anyonic Tonks-Girardeau gas at finite temperature

abstract: We derive an exact description of the non-equilibrium dynamics at finite temperature for the anyonic Tonks-Girardeau gas extending the results of Atas et al. [Phys. Rev. A 95, 043622 (2017)] to the case of arbitrary statistics. The one-particle reduced density matrix is expressed as the Fredholm minor of an integral operator with the kernel being the one-particle Green's function of free fermions at finite temperature and the statistics parameter determining the constant in front of the integral operator. We show that the numerical evaluation of this representation using Nystr\"{o}m's method significantly outperforms the other approaches present in the literature when there are no analytical expressions for the overlaps of the wave-functions. We illustrate the distinctive features and novel phenomena present in the dynamics of anyonic systems in two experimentally relevant scenarios: the quantum Newton's cradle setting and the breathing oscillations initiated by a sudden change of the trap frequency.


title: Dynamical Evolution of Gravitational Leptogenesis

abstract: Radiatively-induced gravitational leptogenesis is a potential mechanism to explain the observed matter-antimatter asymmetry of the universe. Gravitational tidal effects at the quantum loop level modify the dynamics of the leptons in curved spacetime and may be encoded in a low-energy effective action Seff. It has been shown in previous work how in a high-scale BSM theory the CP odd curvature-induced interactions in Seff modify the dispersion relations of leptons and antileptons differently in an expanding universe, giving rise to an effective chemical potential and a non-vanishing equilibrium lepton-antilepton asymmetry. In this paper, the CP even curvature interactions are shown to break lepton number current conservation and modify the evolution of the lepton number density as the universe expands. These effects are implemented in a generalised Boltzmann equation and used to trace the dynamical evolution of the lepton number density in different cosmological scenarios. The theory predicts a potentially significant gravitationally-induced lepton-antilepton asymmetry at very early times in the evolution of the universe.


title: Building an irreversible Carnot-like heat engine with an overdamped harmonic oscillator

abstract: We analyse non-equilibrium Carnot-like cycles built with a colloidal particle in a harmonic trap, which is immersed in a fluid that acts as a heat bath. Our analysis is carried out in the overdamped regime. The cycle comprises four branches: two isothermal processes and two \textit{locally} adiabatic ones. In the latter, both the temperature of the bath and the stiffness of the harmonic trap vary in time, but in such a way that the average heat vanishes for all times. All branches are swept at a finite rate and, therefore, the corresponding processes are irreversible, not quasi-static. Specifically, we are interested in optimising the heat engine to deliver the maximum power and characterising the corresponding values of the physical parameters. The efficiency at maximum power is shown to be very close to the Curzon-Ahlborn bound over the whole range of the ratio of temperatures of the two thermal baths, pointing to the near optimality of the proposed protocol.


title: Digit Stability Inference for Iterative Methods Using Redundant Number Representation

abstract: In our recent work on iterative computation in hardware, we showed that arbitrary-precision solvers can perform more favorably than their traditional arithmetic equivalents when the latter's precisions are either under- or over-budgeted for the solution of the problem at hand. Significant proportions of these performance improvements stem from the ability to infer the existence of identical most-significant digits between iterations. This technique uses properties of algorithms operating on redundantly represented numbers to allow the generation of those digits to be skipped, increasing efficiency. It is unable, however, to guarantee that digits will stabilize, i.e., never change in any future iteration. In this article, we address this shortcoming, using interval and forward error analyses to prove that digits of high significance will become stable when computing the approximants of systems of linear equations using stationary iterative methods. We formalize the relationship between matrix conditioning and the rate of growth in most-significant digit stability, using this information to converge to our desired results more quickly. Versus our previous work, an exemplary hardware realization of this new technique achieves an up-to 2.2x speedup in the solution of a set of variously conditioned systems using the Jacobi method.


title: Response by the Montreal AI Ethics Institute to the European Commission's Whitepaper on AI

abstract: In February 2020, the European Commission (EC) published a white paper entitled, On Artificial Intelligence - A European approach to excellence and trust. This paper outlines the EC's policy options for the promotion and adoption of artificial intelligence (AI) in the European Union. The Montreal AI Ethics Institute (MAIEI) reviewed this paper and published a response addressing the EC's plans to build an "ecosystem of excellence" and an "ecosystem of trust," as well as the safety and liability implications of AI, the internet of things (IoT), and robotics. MAIEI provides 15 recommendations in relation to the sections outlined above, including: 1) focus efforts on the research and innovation community, member states, and the private sector; 2) create alignment between trading partners' policies and EU policies; 3) analyze the gaps in the ecosystem between theoretical frameworks and approaches to building trustworthy AI; 4) focus on coordination and policy alignment; 5) focus on mechanisms that promote private and secure sharing of data; 6) create a network of AI research excellence centres to strengthen the research and innovation community; 7) promote knowledge transfer and develop AI expertise through Digital Innovation Hubs; 8) add nuance to the discussion regarding the opacity of AI systems; 9) create a process for individuals to appeal an AI system's decision or output; 10) implement new rules and strengthen existing regulations; 11) ban the use of facial recognition technology; 12) hold all AI systems to similar standards and compulsory requirements; 13) ensure biometric identification systems fulfill the purpose for which they are implemented; 14) implement a voluntary labelling system for systems that are not considered high-risk; 15) appoint individuals to the oversight process who understand AI systems well and are able to communicate potential risks.


title: Effective immunity and second waves: a dynamic causal modelling study

abstract: This technical report addresses a pressing issue in the trajectory of the coronavirus outbreak; namely, the rate at which effective immunity is lost following the first wave of the pandemic. This is a crucial epidemiological parameter that speaks to both the consequences of relaxing lockdown and the propensity for a second wave of infections. Using a dynamic causal model of reported cases and deaths from multiple countries, we evaluated the evidence models of progressively longer periods of immunity. The results speak to an effective population immunity of about three months that, under the model, defers any second wave for approximately six months in most countries. This may have implications for the window of opportunity for tracking and tracing, as well as for developing vaccination programmes, and other therapeutic interventions.


abstract: We present Wasserstein Embedding for Graph Learning (WEGL), a novel and fast framework for embedding entire graphs in a vector space, in which various machine learning models are applicable for graph-level prediction tasks. We leverage new insights on defining similarity between graphs as a function of the similarity between their node embedding distributions. Specifically, we use the Wasserstein distance to measure the dissimilarity between node embeddings of different graphs. Unlike prior work, we avoid pairwise calculation of distances between graphs and reduce the computational complexity from quadratic to linear in the number of graphs. WEGL calculates Monge maps from a reference distribution to each node embedding and, based on these maps, creates a fixed-sized vector representation of the graph. We evaluate our new graph embedding approach on various benchmark graph-property prediction tasks, showing state-of-the-art classification performance while having superior computational efficiency. The code is available at https://github.com/navid-naderi/WEGL.


title: Logarithmic Voronoi cells

abstract: We study Voronoi cells in the statistical setting by considering preimages of the maximum likelihood estimator that tessellate an open probability simplex. In general, logarithmic Voronoi cells are convex sets. However, for certain algebraic models, namely finite models, models with ML degree 1, linear models, and log-linear (or toric) models, we show that logarithmic Voronoi cells are polytopes. As a corollary, the algebraic moment map has polytopes for both its fibres and its image, when restricted to the simplex. We also compute non-polytopal logarithmic Voronoi cells using numerical algebraic geometry. Finally, we determine logarithmic Voronoi polytopes for the finite model consisting of all empirical distributions of a fixed sample size. These polytopes are dual to the logarithmic root polytopes of Lie type A, and we characterize their faces.


title: The Role of Verb Semantics in Hungarian Verb-Object Order

abstract: Hungarian is often referred to as a discourse-configurational language, since the structural position of constituents is determined by their logical function (topic or comment) rather than their grammatical function (e.g., subject or object). We build on work by Koml\'osy (1989) and argue that in addition to discourse context, the lexical semantics of the verb also plays a significant role in determining Hungarian word order. In order to investigate the role of lexical semantics in determining Hungarian word order, we conduct a large-scale, data-driven analysis on the ordering of 380 transitive verbs and their objects, as observed in hundreds of thousands of examples extracted from the Hungarian Gigaword Corpus. We test the effect of lexical semantics on the ordering of verbs and their objects by grouping verbs into 11 semantic classes. In addition to the semantic class of the verb, we also include two control features related to information structure, object definiteness and object NP weight, chosen to allow a comparison of their effect size to that of verb semantics. Our results suggest that all three features have a significant effect on verb-object ordering in Hungarian and among these features, the semantic class of the verb has the largest effect. Specifically, we find that stative verbs, such as fed "cover", jelent "mean" and \"ovez "surround", tend to be OV-preferring (with the exception of psych verbs which are strongly VO-preferring) and non-stative verbs, such as b\'ir\'al "judge", cs\"okkent "reduce" and cs\'okol "kiss", verbs tend to be VO-preferring. These findings support our hypothesis that lexical semantic factors influence word order in Hungarian.


title: The dependency of Type Ia Supernova parameters on host galaxy morphology for the Pantheon cosmological sample

abstract: Type Ia Supernovae (SNe Ia) are widely used to measure distances in the Universe. Despite the recent progress achieved in SN Ia standardisation, the Hubble diagram still shows some remaining intrinsic dispersion. The remaining scatter in supernova luminosity could be due to the environmental effects that are not yet accounted for by the current standardisation methods. In this work, we compare the local and global colour $(U-V)$, the local star formation rate, and the host stellar mass to the host galaxy morphology. The observed trends suggest that the host galaxy morphology is a good parameter to characterize the SN Ia environment. Therefore, we study the influence of host galaxy morphology on light-curve parameters of SNe Ia for the Pantheon cosmological supernova sample. We determine the Hubble morphological type of host galaxies for a sub-sample of 330 SNe Ia. We confirm that the SALT2 stretch parameter $x_1$ depends on the host morphology with the $p$-value $\sim10^{-14}$. The supernovae with lower stretch value are hosted mainly by elliptical and lenticular galaxies. No correlation for the SALT2 colour parameter $c$ is found. We also examine Hubble diagram residuals for supernovae hosted by the "Early-type" and "Late-type" morphological groups of galaxies. The analysis reveals that the mean distance modulus residual in early-type galaxies is smaller than the one in late-type galaxies, which means that early-type galaxies contain brighter supernovae. However, we do not observe any difference in the residual dispersion for these two morphological groups. The obtained results are in the line with other analyses showing environmental dependence of SN Ia light-curve parameters and luminosity. We confirm the importance of including a host galaxy parameter into the standardisation procedure of SNe Ia for further cosmological studies.


title: Preserving spectral properties of structured matrices under structured perturbations

abstract: This paper is devoted to the study of preservation of eigenvalues, Jordan structure and complementary invariant subspaces of structured matrices under structured perturbations. Perturbations and structure-preserving perturbations are determined such that a perturbed matrix reproduces a given subspace as an invariant subspace and preserves a pair of complementary invariant subspaces of the unperturbed matrix. These results are further utilized to obtain structure-preserving perturbations which modify certain eigenvalues of a given structured matrix and reproduce a set of desired eigenvalues while keeping the Jordan chains unchanged. Moreover, a no spillover structured perturbation of a structured matrix is obtained whose rank is equal to the number of eigenvalues (including multiplicities) which are modified, and in addition, preserves the rest of the eigenvalues and the corresponding Jordan chains which need not be known. The specific structured matrices considered in this paper form Jordan and Lie algebra corresponding to an orthosymmetric scalar product.


title: Splittings of global Mackey functors and regularity of equivariant Euler classes

abstract: We establish natural splittings for the values of global Mackey functors at orthogonal, unitary and symplectic groups. In particular, the restriction homomorphisms between the orthogonal, unitary and symplectic groups of adjacent dimensions are naturally split epimorphisms. The interest in the splitting comes from equivariant stable homotopy theory. The equivariant stable homotopy groups of every global spectrum form a global Mackey functor, so the splittings imply that certain long exact homotopy group sequences separate into short exact sequences. For the real and complex global Thom spectra $\mathbf{MO}$ and $\mathbf{MU}$, the splittings imply the regularity of various Euler classes related to the tautological representations of $O(n)$ and $U(n)$.


title: SAMBA: Safe Model-Based & Active Reinforcement Learning

abstract: In this paper, we propose SAMBA, a novel framework for safe reinforcement learning that combines aspects from probabilistic modelling, information theory, and statistics. Our method builds upon PILCO to enable active exploration using novel(semi-)metrics for out-of-sample Gaussian process evaluation optimised through a multi-objective problem that supports conditional-value-at-risk constraints. We evaluate our algorithm on a variety of safe dynamical system benchmarks involving both low and high-dimensional state representations. Our results show orders of magnitude reductions in samples and violations compared to state-of-the-art methods. Lastly, we provide intuition as to the effectiveness of the framework by a detailed analysis of our active metrics and safety constraints.


title: A Study of Compositional Generalization in Neural Models

abstract: Compositional and relational learning is a hallmark of human intelligence, but one which presents challenges for neural models. One difficulty in the development of such models is the lack of benchmarks with clear compositional and relational task structure on which to systematically evaluate them. In this paper, we introduce an environment called ConceptWorld, which enables the generation of images from compositional and relational concepts, defined using a logical domain specific language. We use it to generate images for a variety of compositional structures: 2x2 squares, pentominoes, sequences, scenes involving these objects, and other more complex concepts. We perform experiments to test the ability of standard neural architectures to generalize on relations with compositional arguments as the compositional depth of those arguments increases and under substitution. We compare standard neural networks such as MLP, CNN and ResNet, as well as state-of-the-art relational networks including WReN and PrediNet in a multi-class image classification setting. For simple problems, all models generalize well to close concepts but struggle with longer compositional chains. For more complex tests involving substitutivity, all models struggle, even with short chains. In highlighting these difficulties and providing an environment for further experimentation, we hope to encourage the development of models which are able to generalize effectively in compositional, relational domains.


title: Off-policy Bandits with Deficient Support

abstract: Learning effective contextual-bandit policies from past actions of a deployed system is highly desirable in many settings (e.g. voice assistants, recommendation, search), since it enables the reuse of large amounts of log data. State-of-the-art methods for such off-policy learning, however, are based on inverse propensity score (IPS) weighting. A key theoretical requirement of IPS weighting is that the policy that logged the data has "full support", which typically translates into requiring non-zero probability for any action in any context. Unfortunately, many real-world systems produce support deficient data, especially when the action space is large, and we show how existing methods can fail catastrophically. To overcome this gap between theory and applications, we identify three approaches that provide various guarantees for IPS-based learning despite the inherent limitations of support-deficient data: restricting the action space, reward extrapolation, and restricting the policy space. We systematically analyze the statistical and computational properties of these three approaches, and we empirically evaluate their effectiveness. In addition to providing the first systematic analysis of support-deficiency in contextual-bandit learning, we conclude with recommendations that provide practical guidance.


title: Goodness-of-Fit Test for Mismatched Self-Exciting Processes

abstract: Recently there have been many research efforts in developing generative models for self-exciting point processes, partly due to their broad applicability for real-world applications. However, rarely can we quantify how well the generative model captures the nature or ground-truth since it is usually unknown. The challenge typically lies in the fact that the generative models typically provide, at most, good approximations to the ground-truth (e.g., through the rich representative power of neural networks), but they cannot be precisely the ground-truth. We thus cannot use the classic goodness-of-fit (GOF) test framework to evaluate their performance. In this paper, we develop a GOF test for generative models of self-exciting processes by making a new connection to this problem with the classical statistical theory of Quasi-maximum-likelihood estimator (QMLE). We present a non-parametric self-normalizing statistic for the GOF test: the Generalized Score (GS) statistics, and explicitly capture the model misspecification when establishing the asymptotic distribution of the GS statistic. Numerical simulation and real-data experiments validate our theory and demonstrate the proposed GS test's good performance.


title: The Bethe-Salpeter Equation Formalism: From Physics to Chemistry

abstract: The Bethe-Salpeter equation (BSE) formalism is steadily asserting itself as a new efficient and accurate tool in the ensemble of computational methods available to chemists in order to predict optical excitations in molecular systems. In particular, the combination of the so-called $GW$ approximation, giving access to reliable ionization energies and electron affinities, and the BSE formalism, able to model UV/Vis spectra, has shown to provide accurate singlet excitation energies with a typical error of $0.1$--$0.3$ eV. With a similar computational cost as time-dependent density-functional theory (TD-DFT), BSE is able to provide an accuracy on par with the most accurate global and range-separated hybrid functionals without the unsettling choice of the exchange-correlation functional, resolving further known issues (\textit{e.g.}, charge-transfer excitations). In this \textit{Perspective} article, we provide a historical overview of BSE, with a particular focus on its condensed-matter roots. We also propose a critical review of its strengths and weaknesses in different chemical situations.

categories: physics.chem-ph cond-mat.mtrl-sci cond-mat.str-el physics.comp-ph


title: Real-time 3D Nanoscale Coherent Imaging via Physics-aware Deep Learning

abstract: Phase retrieval, the problem of recovering lost phase information from measured intensity alone, is an inverse problem that is widely faced in various imaging modalities ranging from astronomy to nanoscale imaging. The current process of phase recovery is iterative in nature. As a result, the image formation is time-consuming and computationally expensive, precluding real-time imaging. Here, we use 3D nanoscale X-ray imaging as a representative example to develop a deep learning model to address this phase retrieval problem. We introduce 3D-CDI-NN, a deep convolutional neural network and differential programming framework trained to predict 3D structure and strain solely from input 3D X-ray coherent scattering data. Our networks are designed to be "physics-aware" in multiple aspects; in that the physics of x-ray scattering process is explicitly enforced in the training of the network, and the training data are drawn from atomistic simulations that are representative of the physics of the material. We further refine the neural network prediction through a physics-based optimization procedure to enable maximum accuracy at lowest computational cost. 3D-CDI-NN can invert a 3D coherent diffraction pattern to real-space structure and strain hundreds of times faster than traditional iterative phase retrieval methods, with negligible loss in accuracy. Our integrated machine learning and differential programming solution to the phase retrieval problem is broadly applicable across inverse problems in other application areas.

categories: eess.IV cond-mat.mtrl-sci cs.LG physics.app-ph


title: On the Complexity of Solving Generic Over-determined Bilinear Systems

abstract: In this paper, we study the complexity of solving generic over-determined bilinear systems over a finite field $\mathbb{F}$. Given a generic bilinear sequence $B \in \mathbb{F}[\mathbf{x},\mathbf{y}]$, with respect to a partition of variables $\mathbf{x}$, $\mathbf{y}$, we show that, the solutions of the system $B= \mathbf{0}$ can be efficiently found on the $\mathbb{F}[\mathbf{y}]$-module generated by $B$. Following this observation, we propose three variations of Gr\"obner basis algorithms, that only involve multiplication by monomials in they-variables, namely, $\mathbf{y}$-XL, based on the XL algorithm, $\mathbf{y}$-MLX, based on the mutant XL algorithm, and $\mathbf{y}$-HXL, basedon a hybrid approach. We define notions of regularity for over-determined bilinear systems,that capture the idea of genericity, and we develop the necessary theoretical tools to estimate the complexity of the algorithms for such sequences. We also present extensive experimental results, testing our conjecture, verifying our results, and comparing the complexity of the various methods.


title: Novel Heavy-Quark Physics Phenomena

abstract: We review the current understanding of heavy quark parton distributions in nucleons and their impact on deep inelastic scattering, collider physics, and other processes at high energies. The determination of the heavy-quark parton distribution functions is particularly significant for the analysis of hard processes at LHC energies, including the forward rapidity high $x_\mathrm{F}$ domain. The contribution of "intrinsic" heavy quarks, which are multiply connected to the valence quarks of nucleons, is reviewed within non-perturbative physics which provides new information on the fundamental structure of hadrons in QCD. A new prediction for the non-perturbative intrinsic charm-anticharm asymmetry of the proton eigenstate has recently been obtained from a QCD lattice gauge theory calculation of the proton's $G_\mathrm{E}^p(Q^2)$ form factor. This form factor only arises from non-valence quarks and anti-quarks if they have different contributions in the proton's eigenstate. This result, together with the exclusive and inclusive connection and analytic constraints on the form of hadronic structure functions from Light-Front Holographic QCD (LFHQCD) predicts a significant non-perturbative $c(x,Q) - \bar{c}(x,Q)$ asymmetry in the proton structure function at high $x$, consistent with the dynamics predicted by intrinsic charm models. Recent ATLAS data on the associated production of prompt photons and charm-quark jets in $pp$ collisions at $\sqrt{s} = 8$ TeV has provided new constraints on non-perturbative intrinsic charm and tests of the LGTH predictions. We also focus on other experimental observables which have high sensitivity to the intrinsic heavy contributions to PDFs.


title: The transport-structural correspondence across the nematic phase transition probed by elasto-x-ray diffraction

abstract: Electronic nematicity in iron pnictide materials is coupled to both the lattice and the conducting electrons, which allows both structural and transport observables to probe nematic fluctuations and the order parameter. Here we combine simultaneous transport and x-ray diffraction measurements with in-situ tunable strain (elasto-XRD) to measure the temperature dependence of the shear modulus and elastoresistivity above the nematic transition and the spontaneous orthorhombicity and resistivity anisotropy below the nematic transition, all within a single sample of $Ba(Fe_{0.96}Co_{0.04})_{2} As_{2}$. The ratio of transport to structural quantities is nearly temperature-independent over a 74 K range and agrees between the ordered and disordered phases. These results show that elasto-XRD is a powerful technique to probe the nemato-elastic and nemato-transport couplings, which have important implications to the nearby superconductivity. It also enables the measurement in the large strain limit, where the breakdown of mean field description reveals the intertwined nature of nematicity.


title: Algebraic and combinatorial expansion in random simplicial complexes

abstract: In this paper we consider the expansion properties and the spectrum of the combinatorial Laplace operator of a $d$-dimensional Linial-Meshulam random simplicial complex, above the cohomological connectivity threshold. We consider the spectral gap of the Laplace operator and the Cheeger constant as this was introduced by Parzanchevski, Rosenthal and Tessler ($Combinatorica$ 36, 2016). We show that with high probability the spectral gap of the random simplicial complex as well as the Cheeger constant are both concentrated around the minimum co-degree of among all $d-1$-faces. Furthermore, we consider a generalisation of a random walk on such a complex and show that the associated conductance is with high probability bounded away from 0.


title: Real-Time Regression with Dividing Local Gaussian Processes

abstract: The increased demand for online prediction and the growing availability of large data sets drives the need for computationally efficient models. While exact Gaussian process regression shows various favorable theoretical properties (uncertainty estimate, unlimited expressive power), the poor scaling with respect to the training set size prohibits its application in big data regimes in real-time. Therefore, this paper proposes dividing local Gaussian processes, which are a novel, computationally efficient modeling approach based on Gaussian process regression. Due to an iterative, data-driven division of the input space, they achieve a sublinear computational complexity in the total number of training points in practice, while providing excellent predictive distributions. A numerical evaluation on real-world data sets shows their advantages over other state-of-the-art methods in terms of accuracy as well as prediction and update speed.


title: Agent Modelling under Partial Observability for Deep Reinforcement Learning

abstract: Modelling the behaviours of other agents is essential for understanding how agents interact and making effective decisions. Existing methods for agent modelling commonly assume knowledge of the local observations and chosen actions of the modelled agents during execution. To eliminate this assumption, we extract representations from the local information of the controlled agent using encoder-decoder architectures. Using the observations and actions of the modelled agents during training, our models learn to extract representations about the modelled agents conditioned only on the local observations of the controlled agent. The representations are used to augment the controlled agent's decision policy which is trained via deep reinforcement learning; thus, during execution, the policy does not require access to other agents' information. We provide a comprehensive evaluation and ablations studies in cooperative, competitive and mixed multi-agent environments, showing that our method achieves higher returns than baseline methods which do not use the learned representations.


title: A Liouville theorem on asymptotically Calabi spaces

abstract: In this paper, we will study harmonic functions on the complete and incomplete spaces with nonnegative Ricci curvature which exhibit inhomogeneous collapsing behaviors at infinity. The main result states that any nonconstant harmonic function on such spaces yields a definite exponential growth rate which depends explicitly on the geometric data at infinity.


title: Network Diffusions via Neural Mean-Field Dynamics

abstract: We propose a novel learning framework based on neural mean-field dynamics for inference and estimation problems of diffusion on networks. Our new framework is derived from the Mori-Zwanzig formalism to obtain an exact evolution of the node infection probabilities, which renders a delay differential equation with memory integral approximated by learnable time convolution operators, resulting in a highly structured and interpretable RNN. Directly using cascade data, our framework can jointly learn the structure of the diffusion network and the evolution of infection probabilities, which are cornerstone to important downstream applications such as influence maximization. Connections between parameter learning and optimal control are also established. Empirical study shows that our approach is versatile and robust to variations of the underlying diffusion network models, and significantly outperform existing approaches in accuracy and efficiency on both synthetic and real-world data.


title: Noise2Inpaint: Learning Referenceless Denoising by Inpainting Unrolling

abstract: Deep learning based image denoising methods have been recently popular due to their improved performance. Traditionally, these methods are trained in a supervised manner, requiring a set of noisy input and clean target image pairs. More recently, self-supervised approaches have been proposed to learn denoising from only noisy images. These methods assume that noise across pixels is statistically independent, and the underlying image pixels show spatial correlations across neighborhoods. These methods rely on a masking approach that divides the image pixels into two disjoint sets, where one is used as input to the network while the other is used to define the loss. However, these previous self-supervised approaches rely on a purely data-driven regularization neural network without explicitly taking the masking model into account. In this work, building on these self-supervised approaches, we introduce Noise2Inpaint (N2I), a training approach that recasts the denoising problem into a regularized image inpainting framework. This allows us to use an objective function, which can incorporate different statistical properties of the noise as needed. We use algorithm unrolling to unroll an iterative optimization for solving this objective function and train the unrolled network end-to-end. The training paradigm follows the masking approach from previous works, splitting the pixels into two disjoint sets. Importantly, one of these is now used to impose data fidelity in the unrolled network, while the other still defines the loss. We demonstrate that N2I performs successful denoising on real-world datasets, while better preserving details compared to its purely data-driven counterpart Noise2Self.

categories: eess.IV cs.CV cs.LG eess.SP physics.med-ph


title: Iterative trajectory reweighting for estimation of equilibrium and non-equilibrium observables

abstract: We present two algorithms by which a set of short, unbiased trajectories can be iteratively reweighted to obtain various observables. The first algorithm estimates the stationary (steady state) distribution of a system by iteratively reweighting the trajectories based on the average probability in each state. The algorithm applies to equilibrium or non-equilibrium steady states, exploiting the `left' stationarity of the distribution under dynamics -- i.e., in a discrete setting, when the column vector of probabilities is multiplied by the transition matrix expressed as a left stochastic matrix. The second procedure relies on the `right' stationarity of the committor (splitting probability) expressed as a row vector. The algorithms are unbiased, do not rely on computing transition matrices, and make no Markov assumption about discretized states. Here, we apply the procedures to a one-dimensional double-well potential, and to a 208$\mu$s atomistic Trp-cage folding trajectory from D.E. Shaw Research.


title: On the Divergence of First Order Resonance Widths at Low Eccentricities

abstract: Orbital resonances play an important role in the dynamics of planetary systems. Classical theoretical analyses found in textbooks report that libration widths of first order mean motion resonances diverge for nearly circular orbits. Here we examine the nature of this divergence with a non-perturbative analysis of a few first order resonances interior to a Jupiter-mass planet. We show that a first order resonance has two branches, the pericentric and the apocentric resonance zone. As the eccentricity approaches zero, the centers of these zones diverge away from the nominal resonance location but their widths shrink. We also report a novel finding of "bridges" between adjacent first order resonances: at low eccentricities, the apocentric libration zone of a first order resonance smoothly connects with the pericentric libration zone of the neighboring first order resonance. These bridges may facilitate resonant migration across large radial distances in planetary systems, entirely in the low eccentricity regime.


abstract: In this paper, we consider Rota-Baxter operators on involutive associative algebras. We define cohomology for Rota-Baxter operators on involutive algebras that governs the formal deformation of the operator. This cohomology can be seen as the Hochschild cohomology of a certain involutive associative algebra with coefficients in a suitable involutive bimodule. We also relate this cohomology with the cohomology of involutive dendriform algebras. Finally, we show that the standard Fard-Guo construction of the functor from the category of dendriform algebras to the category of Rota-Baxter algebras restricts to the involutive case.


title: Interpretable multimodal fusion networks reveal mechanisms of brain cognition

abstract: Multimodal fusion benefits disease diagnosis by providing a more comprehensive perspective. Developing algorithms is challenging due to data heterogeneity and the complex within- and between-modality associations. Deep-network-based data-fusion models have been developed to capture the complex associations and the performance in diagnosis has been improved accordingly. Moving beyond diagnosis prediction, evaluation of disease mechanisms is critically important for biomedical research. Deep-network-based data-fusion models, however, are difficult to interpret, bringing about difficulties for studying biological mechanisms. In this work, we develop an interpretable multimodal fusion model, namely gCAM-CCL, which can perform automated diagnosis and result interpretation simultaneously. The gCAM-CCL model can generate interpretable activation maps, which quantify pixel-level contributions of the input features. This is achieved by combining intermediate feature maps using gradient-based weights. Moreover, the estimated activation maps are class-specific, and the captured cross-data associations are interest/label related, which further facilitates class-specific analysis and biological mechanism analysis. We validate the gCAM-CCL model on a brain imaging-genetic study, and show gCAM-CCL's performed well for both classification and mechanism analysis. Mechanism analysis suggests that during task-fMRI scans, several object recognition related regions of interests (ROIs) are first activated and then several downstream encoding ROIs get involved. Results also suggest that the higher cognition performing group may have stronger neurotransmission signaling while the lower cognition performing group may have problem in brain/neuron development, resulting from genetic variations.


title: Consistent Recalibration Models and Deep Calibration

abstract: Consistent Recalibration models (CRC) have been introduced to capture in necessary generality the dynamic features of term structures of derivatives' prices. Several approaches have been suggested to tackle this problem, but all of them, including CRC models, suffered from numerical intractabilities mainly due to the presence of complicated drift terms or consistency conditions. We overcome this problem by machine learning techniques, which allow to store the crucial drift term's information in neural network type functions. This yields first time dynamic term structure models which can be efficiently simulated.


title: Crystallographic orientation errors in mechanical exfoliation

abstract: We evaluate the effect of mechanical exfoliation of van der Waals materials on crystallographic orientations of the resulting flakes. Flakes originating from a single crystal of graphite, whose orientation is confirmed using STM, are studied using facet orientations and electron back-scatter diffraction (EBSD). While facets exhibit a wide distribution of angles after a single round of exfoliation ($ \sigma \approx 5^o $), EBSD shows that the true crystallographic orientations are more narrowly distributed ($ \sigma \approx 1.5^o $), and facets have an approximately error from the true orientation. Furthermore, we find that the majority of graphite fractures are along armchair lines, and that the cleavage process results in an increase of the zigzag lines portion. Our results place values on the rotation caused by a single round of the exfoliation process, and suggest that when a 1-2 degree precision is necessary, the orientation of a flake can be gauged by the orientation of the macroscopic single crystal from which it was exfoliated.


title: Blind Numerology Identification for Mixed Numerologies

abstract: 5G New Radio (NR) introduces new flexibility that different numerologies can be selected to meet the requirements of a wide variety of services. For this new structure, blind numerology identification can increase system efficiency. Therefore, we propose a blind identification method for mixed numerologies. An autocorrelation method is applied in the time domain by correlating the cyclic prefix (CP) signal of the candidate numerology in the received composite signal for numerology type identification. Then, the location of each numerology in the frequency domain is identified by the variance difference in the power spectral density (PSD) of the subbands, on which different numerologies are occupied. The simulation results are obtained under additive white Gaussian noise (AWGN) and frequencyselective channels. The obtained results show that the proposed method has a robust identification accuracy and a satisfactory BER performance as compared to the non-blind identification approach in the conventional mixed-numerology system.


title: Stability of metric measure spaces with integral Ricci curvature bounds

abstract: In this article we study stability and compactness w.r.t. measured Gromov-Hausdorff convergence of smooth metric measure spaces with integral Ricci curvature bounds. More precisely, we prove that a sequence of $n$-dimensional Riemannian manifolds subconverges to a metric measure space that satisfies the curvature-dimension condition $CD(K,n)$ in the sense of Lott-Sturm-Villani provided the $L^p$-norm for $p>\frac{n}{2}$ of the part of the Ricci curvature that lies below $K$ converges to $0$. The results also hold for sequences of general smooth metric measure spaces $(M,g_M, e^{-f}\mbox{vol}_M)$ where Bakry-Emery curvature replaces Ricci curvature. Corollaries are a Brunn-Minkowski-type inequality, a Bonnet-Myers estimate and a statement on finiteness of the fundamental group. Together with a uniform noncollapsing condition the limit even satisfies the Riemannian curvature-dimension condition $RCD(K,N)$. This implies volume and diameter almost rigidity theorems.


title: Observation of tunable single-atom Yu-Shiba-Rusinov states

abstract: The coupling of a spin to an underlying substrate is the basis for a plethora of phenomena. In the case of a metallic substrate, Kondo screening of the adatom magnetic moment can occur. As the substrate turns superconducting, an intriguing situation emerges where the pair breaking due to the adatom spins leads to Yu-Shiba-Rusinov bound states, but also intertwines with Kondo phenomena. Through scanning tunneling spectroscopy, we analyze the interdependence of Kondo screening and superconductivity. Our data obtained on single Fe adatoms on Nb(110) show that the coupling and the resulting YSR states are strongly adsorption site-dependent and reveal a quantum phase transition at a Kondo temperature comparable to the superconducting gap. The experimental signatures are rationalized by combined density functional theory and continuous-time quantum Monte-Carlo calculations to rigorously treat magnetic and hybridization effects on equal footing.


title: Constructing exchangeable pairs by diffusion on manifolds and its application

abstract: We construct a continuous family of exchangeable pairs by perturbing the random variable through diffusion processes on manifold in order to apply Stein method to certain geometric settings. We compare our perturbation by diffusion method with other approaches of building exchangeable pairs and show that our perturbation scheme cooperates with the infinitesimal version of Stein's method harmoniously. More precisely, our exchangeable pairs satisfy a key condition in the infinitesimal Stein's method in general. Based on the exchangeable pairs, we are able to extend the approximate normality of eigenfunctions of Laplacian on compact manifold to eigenfunctions of Witten Laplacian, which is of the form:$\Delta_w = \Delta - \nabla H$. We then apply our abstract theorem to recover a central limit result of linear statistics on sphere. Finally, we prove an an infinitesimal version of Stein's method for exponential distribution and combine it with our continuous family of exchangeable pairs to extend an approximate exponentiality result of $|Tr U|^2$, where $Tr U$ is the trace of the first power of a matrix $U$ sampled from the Haar measure of unitary group, to arbitrary power and its analog for general circular ensemble.


title: Robust Compressed Sensing using Generative Models

abstract: The goal of compressed sensing is to estimate a high dimensional vector from an underdetermined system of noisy linear equations. In analogy to classical compressed sensing, here we assume a generative model as a prior, that is, we assume the vector is represented by a deep generative model $G: \mathbb{R}^k \rightarrow \mathbb{R}^n$. Classical recovery approaches such as empirical risk minimization (ERM) are guaranteed to succeed when the measurement matrix is sub-Gaussian. However, when the measurement matrix and measurements are heavy-tailed or have outliers, recovery may fail dramatically. In this paper we propose an algorithm inspired by the Median-of-Means (MOM). Our algorithm guarantees recovery for heavy-tailed data, even in the presence of outliers. Theoretically, our results show our novel MOM-based algorithm enjoys the same sample complexity guarantees as ERM under sub-Gaussian assumptions. Our experiments validate both aspects of our claims: other algorithms are indeed fragile and fail under heavy-tailed and/or corrupted data, while our approach exhibits the predicted robustness.


title: Selective Question Answering under Domain Shift

abstract: To avoid giving wrong answers, question answering (QA) models need to know when to abstain from answering. Moreover, users often ask questions that diverge from the model's training data, making errors more likely and thus abstention more critical. In this work, we propose the setting of selective question answering under domain shift, in which a QA model is tested on a mixture of in-domain and out-of-domain data, and must answer (i.e., not abstain on) as many questions as possible while maintaining high accuracy. Abstention policies based solely on the model's softmax probabilities fare poorly, since models are overconfident on out-of-domain inputs. Instead, we train a calibrator to identify inputs on which the QA model errs, and abstain when it predicts an error is likely. Crucially, the calibrator benefits from observing the model's behavior on out-of-domain data, even if from a different domain than the test data. We combine this method with a SQuAD-trained QA model and evaluate on mixtures of SQuAD and five other QA datasets. Our method answers 56% of questions while maintaining 80% accuracy; in contrast, directly using the model's probabilities only answers 48% at 80% accuracy.


title: Constraining X-ray reflection in the low-luminosity AGN NGC 3718 using NuSTAR and XMM--Newton

abstract: One distinctive feature of low-luminosity active galactic nuclei (LLAGN) is the relatively weak reflection features they may display in the X-ray spectrum, which can result from the disappearance of the torus with decreasing accretion rates. Some material, however, must surround the active nucleus, i.e., the accretion flow itself and, possibly, a flattened-out or thinned torus. In this work, we study whether reflection is indeed absent or undetectable due to its intrinsically weak features together with the low statistics inherent to LLAGN. Here we focus on NGC 3718 ($L/L_{\rm Edd}\sim10^{-5}$) combining observations from XMM--Newton and the deepest to date NuSTAR (0.5--79 keV) spectrum of a LLAGN, to constrain potential reflectors, and analyze how the fitted coronal parameters depend on the reflection model. We test models representing both an accretion disc (Relxill) and a torus-like (MYTorus and Borus) neutral reflector. From a statistical point of view, reflection is not required, but its inclusion allows to place strong constraints on the geometry and physical features of the surroundings: both neutral reflectors (torus) tested should be Compton thin ($N_H<10^{23.2}$cm$^{-2}$) and preferentially cover a large fraction of the sky. If the reflected light instead arises from an ionized reflector, a highly ionized case is preferred. These models produce an intrinsic power-law spectral index in the range [1.81--1.87], where the torus models result in steeper slopes. The cut-off energy of the power-law emission also changes with the inclusion of reflection models, resulting in constrained values for the disc reflectors and unconstrained values for torus reflectors.


title: Visualization for Histopathology Images using Graph Convolutional Neural Networks

abstract: With the increase in the use of deep learning for computer-aided diagnosis in medical images, the criticism of the black-box nature of the deep learning models is also on the rise. The medical community needs interpretable models for both due diligence and advancing the understanding of disease and treatment mechanisms. In histology, in particular, while there is rich detail available at the cellular level and that of spatial relationships between cells, it is difficult to modify convolutional neural networks to point out the relevant visual features. We adopt an approach to model histology tissue as a graph of nuclei and develop a graph convolutional network framework based on attention mechanism and node occlusion for disease diagnosis. The proposed method highlights the relative contribution of each cell nucleus in the whole-slide image. Our visualization of such networks trained to distinguish between invasive and in-situ breast cancers, and Gleason 3 and 4 prostate cancers generate interpretable visual maps that correspond well with our understanding of the structures that are important to experts for their diagnosis.


title: Medium Amplitude Parallel Superposition (MAPS) Rheology, Part 2: Experimental Protocols and Data Analysis

abstract: An experimental protocol is developed to directly measure the new material functions revealed by medium amplitude parallel superposition (MAPS) rheology. This experimental protocol measures the medium amplitude response of a material to a simple shear deformation composed of three sine waves at different frequencies. Imposing this deformation and measuring the mechanical response reveals a rich data set consisting of up to 19 measurements of the third order complex modulus at distinct three-frequency coordinates. We discuss how the choice of the input frequencies influences the features of the MAPS domain studied by the experiment. A polynomial interpolation method for reducing the bias of measured values from spectral leakage and variance due to noise is discussed, including a derivation of the optimal range of amplitudes for the input signal. This leads to the conclusion that conducting the experiment in a stress-controlled fashion possesses a distinct advantage to the strain-controlled mode. The experimental protocol is demonstrated through measurements of the MAPS response of a model complex fluid: a surfactant solution of wormlike micelles. The resulting data set is indeed large and feature-rich, while still being acquired in a time comparable to similar medium amplitude oscillatory shear (MAOS) experiments. We demonstrate that the data represents measurements of an intrinsic material function by studying its internal consistency, its compatibility with low-frequency predictions for Coleman-Noll simple fluids, and its agreement with data obtained via MAOS amplitude sweeps. Finally, the data is compared to predictions from the corotational Maxwell model to demonstrate the power of MAPS rheology in determining whether a constitutive model is consistent with a material's time-dependent response.


title: Rubidium Isotope Shift Measurement using Noisy Lasers

abstract: Data students collect from the typical advanced undergraduate laboratory on Saturated Absorption Spectroscopy (SAS) of rubidium can be used to measure the isotope shift and thus leads to an estimate of the isotopic ground state energy shift. This helps students refine their `picture' of the atomic ground state. We describe theoretically why this laboratory works well with free-running laser diodes, demonstrate it experimentally using these lasers tuned to either principal near-infrared transitions, and show an extension of the laboratory using the modulation transfer spectroscopy method.


title: Tell Me Something I Don't Know: Randomization Strategies for Iterative Data Mining

abstract: There is a wide variety of data mining methods available, and it is generally useful in exploratory data analysis to use many different methods for the same dataset. This, however, leads to the problem of whether the results found by one method are a reflection of the phenomenon shown by the results of another method, or whether the results depict in some sense unrelated properties of the data. For example, using clustering can give indication of a clear cluster structure, and computing correlations between variables can show that there are many significant correlations in the data. However, it can be the case that the correlations are actually determined by the cluster structure. In this paper, we consider the problem of randomizing data so that previously discovered patterns or models are taken into account. The randomization methods can be used in iterative data mining. At each step in the data mining process, the randomization produces random samples from the set of data matrices satisfying the already discovered patterns or models. That is, given a data set and some statistics (e.g., cluster centers or co-occurrence counts) of the data, the randomization methods sample data sets having similar values of the given statistics as the original data set. We use Metropolis sampling based on local swaps to achieve this. We describe experiments on real data that demonstrate the usefulness of our approach. Our results indicate that in many cases, the results of, e.g., clustering actually imply the results of, say, frequent pattern discovery.


title: Quantifying and mapping covalent bond scission during elastomer fracture

abstract: Many new soft but tough rubbery materials have been recently discovered and new applications such as flexible prosthetics, stretchable electrodes or soft robotics continuously emerge. Yet, a credible multi-scale quantitative picture of damage and fracture of these materials has still not emerged, due to our fundamental inability to disentangle the irreversible scission of chemical bonds along the fracture path from dissipation by internal molecular friction. Here, by coupling new fluorogenic mechanochemistry with quantitative confocal microscopy mapping, we uncover how many and where covalent bonds are broken as an elastomer fractures. Our measurements reveal that bond scission near the crack plane can be delocalized over up to hundreds of micrometers and increase by a factor of 100 depending on temperature and stretch rate, pointing to an intricated coupling between strain rate dependent viscous dissipation and strain dependent irreversible network scission. These findings paint an entirely novel picture of fracture in soft materials, where energy dissipated by covalent bond scission accounts for a much larger fraction of the total fracture energy than previously believed. Our results pioneer the sensitive, quantitative and spatially-resolved detection of bond scission to assess material damage in a variety of soft materials and their applications.


title: Quantum wakes in lattice fermions

abstract: The wake following a vessel in water is a signature interference effect of moving bodies, and, as described by Lord Kelvin, is contained within a constant universal angle. However, wakes may accompany different kinds of moving disturbances in other situations and even in lattice systems. Here, we investigate the effect of moving disturbances on a Fermi lattice gas of ultracold atoms and analyze the novel types of wake patterns that may occur. We show how at half-filling, the wake angles are dominated by the ratio of the hopping energy to the velocity of the disturbance and on the angle of motion relative to the lattice direction. Moreover, we study the difference between wakes left behind a moving particle detector versus that of a moving potential or a moving particle extractor. We show that these scenarios exhibit dramatically different behavior at half-filling, with the "measurement wake" following an idealized detector vanishing, though the motion of the detector does still leaves a trace through a "fluctuation wake." Finally, we discuss the experimental requirements to observe our predictions in ultracold fermionic atoms in optical lattices.

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


title: Cosmological singularities and 2-dimensional dilaton gravity

abstract: We study Big-Bang or -Crunch cosmological singularities in 2-dimensional dilaton-gravity-scalar theories, in general obtained by dimensional reduction of higher dimensional theories. The dilaton potential encodes information about the asymptotic data defining the theories, and encompasses various families such as flat space, $AdS$, conformally $AdS$ as arising from nonconformal branes, and more general nonrelativistic theories. We find a kind of universal near singularity behaviour independent of the dilaton potential, giving universal interrelations between the exponents defining the time behaviour near the cosmological singularity. More detailed analysis using a scaling ansatz enables finding various classes of cosmological backgrounds, recovering known examples such as the $AdS$ Kasner singularity as well finding as new ones. We give some comments on the dual field theory from this point of view.


title: Untangling tradeoffs between recurrence and self-attention in neural networks

abstract: Attention and self-attention mechanisms, are now central to state-of-the-art deep learning on sequential tasks. However, most recent progress hinges on heuristic approaches with limited understanding of attention's role in model optimization and computation, and rely on considerable memory and computational resources that scale poorly. In this work, we present a formal analysis of how self-attention affects gradient propagation in recurrent networks, and prove that it mitigates the problem of vanishing gradients when trying to capture long-term dependencies by establishing concrete bounds for gradient norms. Building on these results, we propose a relevancy screening mechanism, inspired by the cognitive process of memory consolidation, that allows for a scalable use of sparse self-attention with recurrence. While providing guarantees to avoid vanishing gradients, we use simple numerical experiments to demonstrate the tradeoffs in performance and computational resources by efficiently balancing attention and recurrence. Based on our results, we propose a concrete direction of research to improve scalability of attentive networks.


title: Continuous quantitative Helly-type results

abstract: Brazitikos' results on quantititative Helly-type theorems (for the volume and for the diameter) rely on the work of Srivastava on sparsification of John's decompositions. We change this technique by a stronger recent result due to Friedland and Youssef. This, together with an appropriate selection in the accuracy of the approximation, allow us to obtain Helly-type versions which are sensitive to the number of convex sets involved.


title: Guiding Optimizations with Meliora: A Deep Walk down Memory Lane

abstract: Performance models can be very useful for understanding the behavior of applications and hence can help guide design and optimization decisions. Unfortunately, performance modeling of nontrivial computations typically requires significant expertise and human effort. Moreover, even when performed by experts, it is necessarily limited in scope, accuracy, or both. However, since models are not typically available, programmers, compilers or autotuners cannot use them easily to guide optimizations and are limited to heuristic-based methods that potentially take a lot of time to perform unnecessary transformations. We believe that streamlining model generation and making it scalable (both in terms of human effort and code size) would enable dramatic improvements in compilation techniques, as well as manual optimization and autotuning. To that end, we are building the Meliora code analysis infrastructure for machine learning-based performance model generation of arbitrary codes based on static analysis of intermediate language representations. We demonstrate good accuracy in matching known codes and show how Meliora can be used to optimize new codes though reusing optimization knowledge, either manually or in conjunction with an autotuner. When autotuning, Meliora eliminates or dramatically reduces the empirical search space, while generally achieving competitive performance.


title: A demographic microsimulation model with an integrated household alignment method

abstract: Many dynamic microsimulation models have shown their ability to reasonably project detailed population and households using non-data based household formation and dissolution rules. Although, those rules allow modellers to simplify changes in the household construction, they typically fall short in replicating household projections or if applied retrospectively the observed household numbers. Consequently, such models with biased estimation for household size and other household related attributes lose their usefulness in applications that are sensitive to household size, such as in travel demand and housing demand modelling. Nonetheless, these demographic microsimulation models with their associated shortcomings have been commonly used to assess various planning policies which can result in misleading judgements. In this paper, we contribute to the literature of population microsimulation by introducing a fully integrated system of models for different life event where a household alignment method adjusts household size distribution to closely align with any given target distribution. Furthermore, some demographic events that are generally difficult to model, such as incorporating immigrant families into a population, can be included. We illustrated an example of the household alignment method and put it to test in a dynamic microsimulation model that we developed using dymiumCore, a general-purpose microsimulation toolkit in R, to show potential improvements and weaknesses of the method. The implementation of this model has been made publicly available on GitHub.


title: SPEED: Secure, PrivatE, and Efficient Deep learning

abstract: We introduce a deep learning framework able to deal with strong privacy constraints. Based on collaborative learning, differential privacy and homomorphic encryption, the proposed approach advances state-of-the-art of private deep learning against a wider range of threats, in particular the honest-but-curious server assumption. We address threats from both the aggregation server, the global model and potentially colluding data holders. Building upon distributed differential privacy and a homomorphic argmax operator, our method is specifically designed to maintain low communication loads and efficiency. The proposed method is supported by carefully crafted theoretical results. We provide differential privacy guarantees from the point of view of any entity having access to the final model, including colluding data holders, as a function of the ratio of data holders who kept their noise secret. This makes our method practical to real-life scenarios where data holders do not trust any third party to process their datasets nor the other data holders. Crucially the computational burden of the approach is maintained reasonable, and, to the best of our knowledge, our framework is the first one to be efficient enough to investigate deep learning applications while addressing such a large scope of threats. To assess the practical usability of our framework, experiments have been carried out on image datasets in a classification context. We present numerical results that show that the learning procedure is both accurate and private.


title: Infrared-detected AGNs in the Local Universe

abstract: Spitzer/IRAC color selection is a promising technique to identify hot accreting nuclei, that is to say AGN, in galaxies. We demonstrate this using a small sample of SAURON galaxies, and explore this further. The goal of this study is to find a simple and efficient way to reveal optically obscured nuclear accretion in (nearby) galaxies. We apply an infrared selection method to the Spitzer Survey of Stellar Structures in Galaxies (S4G) sample of more than 2500 galaxies, together with its extension sample of more than 400 galaxies. We use the Spitzer colors to find galaxies in the S$^{4}$G survey containing a hot core, suggesting the presence of a strong AGN, and study the detection fraction as a function of morphological type. We test this infrared color selection method by examining the radio properties of the galaxies, using the VLA NVSS and FIRST surveys. Using the radio data, we demonstrate that galaxies displaying hot mid-infrared nuclei stand out as being (candidate) active galaxies. When using, instead of Spitzer, colors from the lower spatial resolution WISE mission, we reproduce these results. Hence multi-band infrared imaging represents a useful tool to uncover optically obscured nuclear activity in galaxies.


title: Uniqueness of a three-dimensional stochastic differential equation

abstract: In order to extend the study of uniqueness property of multi-dimensional systems of stochastic differential equations, in this paper, we look at the following three-dimensional system of equations, of which the two-dimensional case was well-studied before: $dX_t=Y_tdt\quad, dY_t=Z_tdt,\quad dZ_t=|X_t|^{\alpha}dB_t$. We proved that if $(X_0,Y_0,Z_0)\neq(0,0,0)$, and $\frac{3}{4}<\alpha<1$, then the system of equations has a unique solution in the strong sense.


title: Further Novel reductions of Kamp\'e de F\'eriet function

abstract: In a recent paper, Rathie and Pogany established thirty two novel and general reductions of two and three variables generalized hypergeometric functions. In this paper we provide twenty four further novel and general reduction formulas. The results are established by the application of Beta and Gamma integral methods to the three identities involving products of generalized hypergeometric functions obtained earlier by Kim and Rathie. As special cases, we mention some interesting results.