title: Parametric Systems: Verification and Synthesis

abstract: In this paper we study possibilities of using hierarchical reasoning, symbol elimination and model generation for the verification of parametric systems, where the parameters can be constants or functions. Our goal is to automatically provide guarantees that such systems satisfy certain safety or invariance conditions. We analyze the possibility of automatically generating such guarantees in the form of constraints on parameters. We illustrate our methods on several examples


title: Rational hyperbolic discounting

abstract: How much should you receive in a week to be indifferent to \$ 100 in six months? Note that the indifference requires a rule to ensure the similarity between early and late payments. Assuming that rational individuals have low accuracy, then the following rule is valid: if the amounts to be paid are much less than the personal wealth, then the $q$-exponential discounting guarantees indifference in several periods. Thus, the discounting can be interpolated between hyperbolic and exponential functions due to the low accuracy to distinguish time averages when the payments have low impact on personal wealth. Therefore, there are physical conditions that allow the hyperbolic discounting regardless psycho-behavioral assumption.


title: Accelerating Stochastic Quantum Chemistry

abstract: The convergence of full configuration interaction quantum Monte Carlo (FCIQMC) is accelerated using a quasi-Newton propagation (QN) which can also be applied to coupled cluster Monte Carlo (CCMC). The computational scaling of this optimised propagation is O(1), keeping the additional computational cost to a bare minimum. Its effects are investigated deterministically and stochastically on a model system, the uniform electron gas, with Hilbert space size up to $10^{40}$ and shown to accelerate convergence of the instantaneous projected energy by over an order of magnitude in the FCIQMC test case. Its capabilities are then demonstrated with FCIQMC on an archetypical quantum chemistry problem, the chromium dimer, in an all-electron basis set with Hilbert space size of about $10^{22}$ yielding highly accurate FCI energies.


title: The Minimal Abstract Robust Subdifferential

abstract: This paper introduces new subdifferential concepts together with their main properties and place with respect to other subdifferential notions.


title: Burglary in London: Insights from Statistical Heterogeneous Spatial Point Processes

abstract: To obtain operational insights regarding the crime of burglary in London we consider the estimation of effects of covariates on the intensity of spatial point patterns. By taking into account localised properties of criminal behaviour, we propose a spatial extension to model-based clustering methods from the mixture modelling literature. The proposed Bayesian model is a finite mixture of Poisson generalised linear models such that each location is probabilistically assigned to one of the clusters. Each cluster is characterised by the regression coefficients which we subsequently use to interpret the localised effects of the covariates. Using a blocking structure of the study region, our approach allows specifying spatial dependence between nearby locations. We estimate the proposed model using Markov Chain Monte Carlo methods and provide a Python implementation.


title: Diabatic Hamiltonian Construction in van der Waals heterostructure complexes

abstract: A diabatization method is developed for the approximated description of the photoinduced charge separation/transfer processes in the van der Waals (vdW) heterostructure complex, which is based on the wavefunction projection approach using a plane wave basis set in the framework of the single-particle picture. We build the diabatic Hamiltonian for the description of the interlayer photoinduced hole-transfer process of the two-dimensional vdW MoS2/WS2 heterostructure complexes. The diabatic Hamiltonian gives the energies of the localized valence band states (located at MoS2 and valence band states (located at WS2), as well as the couplings between them. The wavefunction projection method provides a practical and reasonable approach to construct the diabatic model in the description of photoinduced charge transfer processes in the vdW heterostructure complexes.


title: Chiral Magnonic Edge States in Ferromagnetic Skyrmion Crystals Controlled by Magnetic Fields

abstract: Achieving control over magnon spin currents in insulating magnets - where dissipation due to Joule heating is highly suppressed - is an active area of research that could lead to energy-efficient spintronics applications. However, magnon spin currents supported by conventional systems with uniform magnetic order have proven hard to control. An alternative approach that relies on topologically protected magnonic edge states of spatially periodic magnetic textures has recently emerged. A prime example of such textures is the ferromagnetic skyrmion crystal which hosts chiral edge states providing a platform for magnon spin currents. Here, we show, for the first time, an external magnetic field can drive a topological phase transition in the spin wave spectrum of a ferromagnetic skyrmion crystal. The topological phase transition is signaled by the closing of a low-energy bulk magnon gap at a critical field. In the topological phase, below the critical field, two topologically protected chiral magnonic edge states lie within this gap, but they unravel in the trivial phase, above the critical field. Remarkably, the topological phase transition involves an inversion of two magnon bands that at the $\Gamma$ point correspond to the breathing and anticlockwise modes of the skyrmions in the crystal. Our findings suggest that an external magnetic field could be used as a knob to switch on and off magnon spin currents carried by topologically protected chiral magnonic edge states.


title: Syntactic Interpolation for Tense Logics and Bi-Intuitionistic Logic via Nested Sequents

abstract: We provide a direct method for proving Craig interpolation for a range of modal and intuitionistic logics, including those containing a "converse" modality. We demonstrate this method for classical tense logic, its extensions with path axioms, and for bi-intuitionistic logic. These logics do not have straightforward formalisations in the traditional Gentzen-style sequent calculus, but have all been shown to have cut-free nested sequent calculi. The proof of the interpolation theorem uses these calculi and is purely syntactic, without resorting to embeddings, semantic arguments, or interpreted connectives external to the underlying logical language. A novel feature of our proof includes an orthogonality condition for defining duality between interpolants.


title: Charge-carrier thermalization in bulk and monolayer CdTe from first principles

abstract: While cadmium telluride (CdTe) thin films are being used in solar cell prototyping for decades, the recent advent of two-dimensional (2D) materials challenges the fundamental limit for thickness of conventional CdTe layers. Here, we report our theoretical predictions on photocarrier dynamics in an ultimately thin (about 1 nm) CdTe slab. It corresponds to a layer that is just a single unit cell thick, when the bulk parent crystal in the zinc blende phase is cleaved along the [110] facet. Using an \textit{ab-initio} method based on density functional theory (DFT) and the Boltzmann equation in the relaxation time approximation (RTA), we determine the thermalization time for charge carriers excited to a certain energy for instance through laser irradiation. Our calculations include contributions arising from all phonon branches in the first Brillouin zone (BZ), thus capturing all relevant inter- and intraband carrier transitions due to electron-phonon scattering. We find that the photocarrier thermalization time is strongly reduced, by one order of magnitude for holes and by three orders of magnitude for electrons, once the CdTe crystal is thinned down from the bulk to a monolayer. Most surprisingly, the electron thermalization time becomes independent of the electron excess energy up to about 0.5~eV, when counted from the conduction band minimum (CBM). We relate this peculiar behavior to the degenerate and nearly parabolic lowest conduction band that yields a constant density of states (DOS) in the 2D limit. Our findings may be useful for designing novel CdTe-based optoelectronic devices, which employ nonequilibrium photoexcited carriers to improve the performance.


title: Working principle and demonstrator of microwave-multiplexing for the HOLMES experiment microcalorimeters

abstract: The determination of the neutrino mass is an open issue in modern particle physics and astrophysics. The direct mass measurement is the only theory-unrelated experimental tool capable to probe such quantity. The HOLMES experiment aims to measure the end-point energy of the electron capture (EC) decay of $^{163}$Ho with a statistical sensitivity on the neutrino mass as low as $\sim 1$ eV/c$^2$. In order to acquire the large needed statistics, by keeping the pile-up contribution as low as possible, 1024 transition edge sensors (TESs) with high energy and time resolutions will be employed. Microcalorimeter and bolometer arrays based on transition edge sensor with thousands of pixels are under development for several space-based and ground-based applications, including astrophysics, nuclear and particle physics, and materials science. The common necessary challenge is to develop pratical multiplexing techniques in order to simplify the cryogenics and readout systems. Despite the various multiplexing variants which are being developed have been successful, new approaches are needed to enable scaling to larger pixel counts and faster sensors, as requested for HOLMES, reducing also the cost and complexity of readout. A very novel technique that meets all of these requirements is based on superconducting microwave resonators coupled to radio-frequency Superconducting Quantum Interference Devices, in which the the changes in the TES input current is tranduced to a change in phase of a microwave signal. In this work we introduce the basics of this technique, the design and development of the first two-channel read out system and its performances with the first TES detectors specifically designed for HOLMES. In the last part we explain how to extend this approach scaling to 1024 pixels.


title: Computing period integrals of rigid double octic Calabi-Yau threefolds with Picard-Fuchs operator

abstract: We present a method for numerical computation of period integrals of a rigid Calabi-Yau threefold using Picard-Fuchs operator of a one-parameter smoothing. Our method gives a possibility of computing the lattice of period integrals of a rigid double octic without any explicit knowledge of its geometric properties, employing only simple facts from the theory of Fuchsian equations and computations in MAPLE with a library for differential equations. As a surprising consequence we also get approximations of additional integrals related to a singular (nodal) model of considered Calabi-Yau threefold.


title: Measuring productivity dispersion: a parametric approach using the L\'{e}vy alpha-stable distribution

abstract: It is well-known that value added per worker is extremely heterogeneous among firms, but relatively little has been done to characterize this heterogeneity more precisely. Here we show that the distribution of value-added per worker exhibits heavy tails, a very large support, and consistently features a proportion of negative values, which prevents log transformation. We propose to model the distribution of value added per worker using the four parameter L\'evy stable distribution, a natural candidate deriving from the Generalised Central Limit Theorem, and we show that it is a better fit than key alternatives. Fitting a distribution allows us to capture dispersion through the tail exponent and scale parameters separately. We show that these parametric measures of dispersion are at least as useful as interquantile ratios, through case studies on the evolution of dispersion in recent years and the correlation between dispersion and intangible capital intensity.


title: Recalibrating the Cosmic Star Formation History

abstract: The calibrations linking observed luminosities to the star formation rate depend on the assumed stellar population synthesis model, initial mass function, star formation and metal enrichment history, and whether reprocessing by dust and gas is included. Consequently the shape and normalisation of the inferred cosmic star formation history is sensitive to these assumptions. Using v2.2.1 of the Binary Population and Spectral Synthesis (\bpass) model we determine a new set of calibration coefficients for the ultraviolet, thermal-infrared, and, hydrogen recombination lines. These ultraviolet and thermal infrared coefficients are 0.15-0.2 dex higher than those widely utilised in the literature while the H$\alpha$ coefficient is $\sim 0.35$ dex larger. These differences arise in part due to the inclusion binary evolution pathways but predominantly reflect an extension in the IMF to 300 $M_{\odot}$ and a change in the choice of reference metallicity. We use these new coefficients to recalibrate the cosmic star formation history, and find improved agreement between the integrated cosmic star formation history and the in-situ measured stellar mass density as a function of redshift. However, these coefficients produce new tension between star formation rate densities inferred from the ultraviolet and thermal-infrared and those from H$\alpha$.


title: Non-Uniform Time-Step Deep Q-Network for Carrier-Sense Multiple Access in Heterogeneous Wireless Networks

abstract: This paper investigates a new class of carrier-sense multiple access (CSMA) protocols that employ deep reinforcement learning (DRL) techniques, referred to as carrier-sense deep-reinforcement learning multiple access (CS-DLMA). The goal of CS-DLMA is to enable efficient and equitable spectrum sharing among a group of co-located heterogeneous wireless networks. Existing CSMA protocols, such as the medium access control (MAC) of WiFi, are designed for a homogeneous network in which all nodes adopt the same protocol. Such protocols suffer from severe performance degradation in a heterogeneous environment where there are nodes adopting other MAC protocols. CS-DLMA aims to circumvent this problem by making use of DRL. In particular, this paper adopts alpha-fairness as the general objective of CS-DLMA. With alpha-fairness, CS-DLMA can achieve a range of different objectives when coexisting with other MACs by changing the value of alpha. A salient feature of CS-DLMA is that it can achieve these objectives without knowing the coexisting MACs through a learning process based on DRL. The underpinning DRL technique in CS-DLMA is deep Q-network (DQN). However, the conventional DQN algorithms are not suitable for CS-DLMA due to their uniform time-step assumption. In CSMA protocols, time steps are non-uniform in that the time duration required for carrier sensing is smaller than the duration of data transmission. This paper introduces a non-uniform time-step formulation of DQN to address this issue. Our simulation results show that CS-DLMA can achieve the general alpha-fairness objective when coexisting with TDMA, ALOHA, and WiFi protocols by adjusting its own transmission strategy. Interestingly, we also find that CS-DLMA is more Pareto efficient than other CSMA protocols when coexisting with WiFi.


title: Charge and statistics of lattice quasiholes from density measurements: a Tree Tensor Network study

abstract: We numerically investigate the properties of the quasihole excitations above the bosonic fractional Chern insulator state at filling $\nu = 1/2$, in the specific case of the Harper-Hofstadter Hamiltonian with hard-core interactions. For this purpose we employ a Tree Tensor Network technique, which allows us to study systems with up to $N=18$ particles on a $16 \times 16$ lattice and experiencing an additional harmonic confinement. First, we observe the quantization of the quasihole charge at fractional values and its robustness against the shape and strength of the impurity potentials used to create and localize such excitations. Then, we numerically characterize quasihole anyonic statistics by applying a discretized version of the relation connecting the statistics of quasiholes in the lowest Landau level to the depletions they create in the density profile [Macaluso et al., arXiv:1903.03011]. Our results give a direct proof of the anyonic statistics for quasiholes of fractional Chern insulators, starting from a realistic Hamiltonian. Moreover, they provide strong indications that this property can be experimentally probed through local density measurements, making our scheme readily applicable in state-of-the-art experiments with ultracold atoms and superconducting qubits.


title: Correlated Gaussians and low-discrepancy sequences

abstract: Within the Correlated Gaussian Method the parameters of the Gaussian basis functions are often chosen stochastically using pseudo-random sequences. We show that alternative low-discrepancy sequences, also known as quasi-random sequences, provide bases of better quality.


title: Phosphine as a Biosignature Gas in Exoplanet Atmospheres

abstract: A long-term goal of exoplanet studies is the identification and detection of biosignature gases. Beyond the most discussed biosignature gas O$_2$, only a handful of gases have been considered in detail. Here we evaluate phosphine (PH$_3$). On Earth, PH$_3$ is associated with anaerobic ecosystems, and as such is a potential biosignature gas on anoxic exoplanets. We simulate CO$_2-$ and H$_2-$dominated habitable terrestrial planet atmospheres. We find that PH$_3$ can accumulate to detectable concentrations on planets with surface production fluxes of 10$^{10}$-10$^{14}$ cm$^{-2}$ s$^{-1}$ (corresponding to surface concentrations of 10s of ppb to 100s of ppm), depending on atmospheric composition and UV flux. While high, such surface fluxes are comparable to the global terrestrial production rate of CH$_4$ (10$^{11}$ cm$^{-2}$ s$^{-1}$) and below the maximum local terrestrial PH$_3$ production rate (10$^{14}$ cm$^{-2}$ s$^{-1}$). As with other gases, PH$_3$ can more readily accumulate on low-UV planets, e.g. planets orbiting quiet M-dwarfs or with a photochemical UV shield. If detected, PH$_3$ is a promising biosignature gas, as it has no known abiotic false positives on terrestrial planets that could generate the high fluxes required for detection. PH$_3$ also has 3 strong spectral features such that in any atmosphere scenario 1 of the 3 will be unique compared to other dominant spectroscopic molecules. PH$_3$'s weakness as a biosignature gas is its high reactivity, requiring high outgassing for detectability. We calculate that 10s of hours of JWST time are required for a potential detection of PH$_3$. Yet because PH$_3$ is spectrally active in the same wavelength regions as other atmospherically important molecules (e.g., H$_2$O and CH$_4$), searches for PH$_3$ can be carried out at no additional observational cost to searches for other molecules relevant to exoplanet habitability.


title: Identifying Exclusive Displaced Hadronic Signatures in the Forward Region of the LHC

abstract: The LHCb detector provides accurate vertex reconstruction and hadronic particle identification, which make the experiment an ideal place to look for light long-lived particles (LLP) decaying into Standard Model (SM) hadrons. In contrast with the typical search strategy relying on energetic jets and a high multiplicity of tracks from the LLP decay, LHCb can identify LLPs in exclusive, specific hadronic final states. To illustrate the idea, we study the sensitivity of LHCb to an exotic Higgs decay $h\to SS$, followed by the displaced decay of GeV-scale scalars into charged kaons $S\to K^+K^-$. We show that the reconstruction of kaon vertices in narrow invariant mass windows can efficiently eliminate the combinatorial backgrounds from $B$-meson decays. While the same signal is extremely difficult to probe in the existing displaced jet searches at ATLAS/CMS, the LHCb search we propose can probe the branching ratio BR$(h\to SS)$ down to $0.1\%$ ($0.02\%$) level with $15$ ($300$) fb$^{-1}$ of data. We also apply this projected bound to two scenarios with Higgs portal couplings, where the scalar mediator $S$ either couples to a) the SM quarks only, or b) to both quarks and leptons in the minimal flavor violation paradigm. In both scenarios we compare the reach of our proposed search with the expected constraints from ATLAS and CMS on the invisible Higgs width and with the constraints from rare B-decays studies at LHCb. We find that for 1 GeV $< m_S < $ 2 GeV and $0.5~{\rm mm} \lesssim c \tau \lesssim 10$ mm our proposed search will be competitive with the ATLAS and CMS projections, while at the same time providing crucial information of the hadronic interactions of $S$, which can not be obtained from the {\it indirect} measurement of the Higgs invisible width.


title: The D_8-tower of weak Jacobi forms and applications

abstract: We construct a tower of arithmetic generators of the bigraded polynomial ring J_{*,*}^{w, O}(D_n) of weak Jacobi modular forms invariant with respect to the full orthogonal group O(D_n) of the root lattice D_n for 2\le n\le 8. This tower corresponds to the tower of strongly reflective modular forms on the orthogonal groups of signature (2,n) which determine the Lorentzian Kac-Moody algebras related to the BCOV (Bershadsky-Cecotti-Ooguri-Vafa)-analytic torsions. We prove that the main three generators of index one of the graded ring satisfy a special system of modular differential equations. We found also a general modular differential equation of the generator of weight 0 and index 1 which generates the automorphic discriminant of the moduli space of Enriques surfaces.


title: Theory of the tertiary instability and the Dimits shift from reduced drift-wave models

abstract: Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity $U(x)$ with respect to the radial coordinate $x$. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $\gamma_{\rm TI}$ is derived within the modified Hasegawa--Wakatani model. We show that $\gamma_{\rm TI}$ equals the primary-instability growth rate plus a term that depends on the local $U''$; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry--Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear $|U'|$ is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.


title: Assessment of the Systematic Uncertainties in the Cosmological Analysis of the SDSS Supernovae Photometric Sample

abstract: Improvement in the precision of measurements of cosmological parameters with Type Ia Supernovae (SNIa) is expected to come from large photometrically identified (photometric) SN samples. Here we re-analyse the SDSS photometric SN sample, with roughly 700 high-quality, likely but unconfirmed SNIa light-curves, to develop new analysis tools aimed at evaluating systematic uncertainties on the dark energy equation-of-state parameter $w$. Since we require a spectroscopically measured host galaxy redshift for each SN, we determine the associated selection efficiency of host galaxies in order to simulate bias corrections. We determine that the mis-association rate of host galaxies is $0.6\%$; ignoring this effect in simulated bias corrections leads to a $w$-bias of $\Delta w = +0.0007$, where $w$ is evaluated from SNIa and priors from measurements of baryon acoustic oscillations and the cosmic microwave background. We assess the uncertainty in our modeling of the host galaxy selection efficiency and find the associated $w$ uncertainty to be $-0.0072$. Finally, we explore new core collapse (CC) models in simulated training samples and find that adjusting the CC luminosity distribution to be in agreement with previous Pan-STARRS analyses yields a better match to the SDSS data. The impact of ignoring this adjustment is $\Delta w = -0.0109$; the impact of replacing the new CC models with those used by Pan-STARRS is $\Delta w = -0.0028$. These systematic uncertainties are subdominant to the statistical constraints from the SDSS sample, but must be considered in future photometric analyses of large SN samples such as those from DES, LSST and WFIRST.


title: Self-propelled motion of a rigid body inside a density dependent incompressible fluid

abstract: This paper is devoted to the existence of a weak solution to a system describing a self-propelled motion of a rigid body in a viscous fluid in the whole $\mathbb{R}^3$. The fluid is modelled by the incompressible nonhomogeneous Navier-Stokes system with a nonnegative density. The motion of the rigid body is described by the balance of linear and angular momentum. We consider the case where slip is allowed at the fluid-solid interface through Navier condition and prove the global existence of a weak solution.


title: A one-loop exact quantization of Chern-Simons theory

abstract: We examine Chern-Simons theory as a deformation of a 3-dimensional BF theory that is partially holomorphic and partially topological. In particular, we introduce a novel gauge that leads naturally to a one-loop exact quantization of this BF theory and Chern-Simons theory. This approach illuminates several important features of Chern-Simons theory, notably the bulk-boundary correspondence of Chern-Simons theory with chiral WZW theory. In addition to rigorously constructing the theory, we also explain how it applies to a large class of closely related 3-dimensional theories and some of the consequences for factorization algebras of observables.


title: R-SQAIR: Relational Sequential Attend, Infer, Repeat

abstract: Traditional sequential multi-object attention models rely on a recurrent mechanism to infer object relations. We propose a relational extension (R-SQAIR) of one such attention model (SQAIR) by endowing it with a module with strong relational inductive bias that computes in parallel pairwise interactions between inferred objects. Two recently proposed relational modules are studied on tasks of unsupervised learning from videos. We demonstrate gains over sequential relational mechanisms, also in terms of combinatorial generalization.


title: Efficient random number generation techniques for CMOS SPAD array based devices

abstract: This work presents new techniques to produce true random bits by exploiting single photon time of arrival. Two FPGA-based QRNG devices are presented: Randy which uses one discrete SPAD and LinoSPAD which uses a CMOS SPAD array, along with a time-to-digital converter (TDC). Post-processing procedures are explained in order to extract randomness taking care of SPAD and TDC non-idealities. These procedures are based on the application of Peres [Y. Peres, Ann. Statist.20, 590 (1992)] and Zhou-Bruk [H. Zhou and J. Bruck, arXiv:1209.0726 (2012)] extraction algorithms. Achieved generation rates are 1.8 Mbit/s for Randy device and 310 Mbit/s for LinoSPAD device. Randy QRNG also features a real time procedure which was used for the realization of fundamental tests of physics.


title: Predicting dynamical system evolution with residual neural networks

abstract: Forecasting time series and time-dependent data is a common problem in many applications. One typical example is solving ordinary differential equation (ODE) systems $\dot{x}=F(x)$. Oftentimes the right hand side function $F(x)$ is not known explicitly and the ODE system is described by solution samples taken at some time points. Hence, ODE solvers cannot be used. In this paper, a data-driven approach to learning the evolution of dynamical systems is considered. We show how by training neural networks with ResNet-like architecture on the solution samples, models can be developed to predict the ODE system solution further in time. By evaluating the proposed approaches on three test ODE systems, we demonstrate that the neural network models are able to reproduce the main dynamics of the systems qualitatively well. Moreover, the predicted solution remains stable for much longer times than for other currently known models.


title: The Heavy Quark Expansion for Inclusive Semileptonic Charm Decays Revisited

abstract: The Heavy Quark Expansion (HQE) has become an extremely powerful tool in flavor physics. For charm decays, where the expansion parameters $\alpha_s(m_c)$ and $\Lambda_{\rm QCD}/m_c$ are bigger than for bottom decays, it remains to be seen if the HQE can be applied with similar success. Nevertheless, to make optimal use of the plethora of data already available and coming in the near future, a better understanding of HQE for charm decays is crucial. This paper discusses in detail how the HQE for charm decays is set up, what is the role of four-quark (weak annihilation) operators and how this compares to the well understood bottom decays. Subtleties concerning radiative corrections and the charm mass scheme are briefly discussed. An experimental study of the relevant HQE hadronic matrix elements will then show if the HQE expansion for charm converges well enough. Besides serving as an important cross check for inclusive $B$ decays, in the end, this study might open the road for inclusive $|V_{cs}|$ and $|V_{cd}|$ extractions.


title: Branching of twins in shape memory alloys revisited

abstract: We study the branching of twins appearing in shape memory alloys at the interface between austenite and martensite. In the framework of three-dimensional non-linear elasticity theory, we propose an explicit, low-energy construction of the branched microstructure, generally applicable to any shape memory material without restrictions on the symmetry class of martensite or on the geometric parameters of the interface. We show that the suggested construction follows the expected energy scaling law, i.e., that (for the surface energy of the twins being sufficiently small) the branching leads to energy reduction. Furthermore, the construction can be modified to capture different features of experimentally observed microstructures without violating this scaling law. By using a numerical procedure, we demonstrate that the proposed construction is able to predict realistically the twin width and the number of branching generations in a Cu-Al-Ni single crystal.


title: Solving A Class of Mean-Field LQG Problems

abstract: In this work, we study a class of mean-field linear quadratic Gaussian (LQG) problems. Under suitable conditions, explicit solutions of the distribution-dependent optimal control problems are obtained. Riccati systems are derived by directly solving the associated master equations. Some extensions on controls with partial observations are also considered.


title: Spin-reversal energy barriers of 305 K for Fe$^{2+}$ $d^{6}$ ions with linear ligand coordination

abstract: A remarkably large magnetic anisotropy energy of 305 K is computed by quantum chemistry methods for divalent Fe$^{2+}$ $d^6$ substitutes at Li-ion sites with $D_{6h}$ point-group symmetry within the solid-state matrix of Li$_3$N. This is similar to values calculated by the same approach and confirmed experimentally for linearly coordinated monovalent Fe$^{1+}$ $d^7$ species, among the largest so far in the research area of single-molecule magnets. Our ab initio results therefore mark a new exciting exploration path in the search for superior single-molecule magnets, rooted in the $d_{xy}^{1.5} d_{x^2-y^2}^{1.5} d_{z^2}^1 d_{yz}^1 d_{zx}^1$ configuration of $d^6$ transition-metal ions with linear or quasilinear nearest-neighbor coordination. This $d^6$ axial anisotropy may be kept robust even for symmetries lower than $D_{6h}$, provided the ligand and farther-neighbor environment is engineered such that the $d_{xy}^{1.5} d_{x^2-y^2}^{1.5} d_{z^2}^1 d_{yz}^1 d_{zx}^1$ - $d_{xy}^{1} d_{x^2-y^2}^{1} d_{z^2}^2 d_{yz}^1 d_{zx}^1$ splitting remains large enough.


title: Formation of Primordial Black Holes from Warm Inflation

abstract: Primordial Black Holes (PBHs) serve as a unique probe to the physics of the early Universe, particularly inflation. In light of this, we study the formation of PBHs by the collapse of overdense perturbations generated during a model of warm inflation. For our model, we find that the primordial curvature power spectrum is red-tilted (spectral index $n_s<1$) at the large scale (small $k$) and is consistent with the $n_s-r$ values allowed from the CMB observations. Along with that, it has a blue-tilt ($n_s>1$) for the small PBH scales (large $k$), with a sufficiently large amplitude of the primordial curvature power spectrum required to form PBHs. These features originate because of the inflaton's coupling with the other fields during warm inflation. We discuss the role of the inflaton dissipation to the enhancement in the primordial power spectrum at the PBH scales. We find that for some parameter range of our warm inflation model, PBHs with mass $\sim 10^3$ g can be formed with significant abundance. Such tiny mass PBHs have a short lifetime $\sim 10^{-19}$ s and would have evaporated into Hawking radiation in the early Universe. Further in this study, we discuss the evaporation constraints on the initial mass fraction of the generated PBHs and the possibility of Planck mass PBH relics to constitute the dark matter.


title: Kondo physics of magnetic adatoms on metallic surfaces when the onset of the surface conduction density of states crosses the Fermi level

abstract: We study the role of the onset of Shockley states, $D_s$, belonging to (111) surfaces of Cu, Ag and Au in the Kondo effect when a magnetic impurity is deposited on them. When $D_s$ approaches to the Fermi level, $E_F$, thing that can be done by compressing (stretching) the metallic sample, we found that most of the thermodynamic and dynamic properties of the impurity are affected in a non trivial way. We model the system by a generic Anderson impurity model and solve it by using the numerical renormalization group, NRG, technique. In particular, the impurity contribution to magnetic susceptibility and entropy as a function of temperature exhibit negative values and goes to zero slowly in a logarithmic shape. Furthermore, we found a suppression of the spectral density weight at the Fermi level when $D_s\sim E_F$ even in the Kondo regime. As a consequence, the conductance through the impurity is strongly reduced by near $25\%$ of the unitary value $2e^2/h$. Finally, we analyze these features in realistic systems like Co on Ag(111) reported in the literature.


title: Defence Against the Modern Arts: the Curse of Statistics "Score-based likelihood ratios"

abstract: For several decades, legal and scientific scholars have argued that conclusions from forensic examinations should be supported by statistical data and reported within a probabilistic framework. Multiple models have been proposed to quantify the probative value of forensic evidence. Unfortunately, several of these models rely on ad-hoc strategies that are not scientifically sound. The opacity of the technical jargon used to present these models and their results, and the complexity of the techniques involved make it very difficult for the untrained user to separate the wheat from the chaff. This series of papers is intended to help forensic scientists and lawyers recognise limitations and issues in tools proposed to interpret the results of forensic examinations. This paper focuses on tools that have been proposed to leverage the use of similarity scores to assess the probative value of forensic findings. We call this family of tools "score-based likelihood ratios". In this paper, we present the fundamental concepts on which these tools are built, we describe some specific members of this family of tools, and we explore their convergence to the Bayes factor through an intuitive geometrical approach and through simulations. Finally, we discuss their validation and their potential usefulness as a decision-making tool in forensic science.


title: Combinatorial Inscribability Obstructions for Higher-Dimensional Polytopes

abstract: For $3$-dimensional convex polytopes, inscribability is a classical property that is relatively well-understood due to its relation with Delaunay subdivisions of the plane and hyperbolic geometry. In particular, inscribability can be tested in polynomial time, and for every $f$-vector of $3$-polytopes, there exists an inscribable polytope with that $f$-vector. For higher-dimensional polytopes, much less is known. Of course, for any inscribable polytope, all of its lower-dimensional faces need to be inscribable, but this condition does not appear to be very strong. We observe non-trivial new obstructions to the inscribability of polytopes that arise when imposing that a certain inscribable face be inscribed. Using this obstruction, we show that the duals of the $4$-dimensional cyclic polytopes with at least $8$ vertices---all of whose faces are inscribable---are not inscribable. This result is optimal in the following sense: We prove that the duals of the cyclic $4$-polytopes with up to $7$ vertices are, in fact, inscribable. Moreover, we interpret this obstruction combinatorially as a forbidden subposet of the face lattice of a polytope, show that $d$-dimensional cyclic polytopes with at least $d+4$ vertices are not circumscribable, and that no dual of a neighborly $4$-polytope with $8$ vertices, that is, no polytope with $f$-vector $(20,40,28,8)$, is inscribable.


title: Semi-Automatic Crowdsourcing Tool for Online Food Image Collection and Annotation

abstract: Assessing dietary intake accurately remains an open and challenging research problem. In recent years, image-based approaches have been developed to automatically estimate food intake by capturing eat occasions with mobile devices and wearable cameras. To build a reliable machine-learning models that can automatically map pixels to calories, successful image-based systems need large collections of food images with high quality groundtruth labels to improve the learned models. In this paper, we introduce a semi-automatic system for online food image collection and annotation. Our system consists of a web crawler, an automatic food detection method and a web-based crowdsoucing tool. The web crawler is used to download large sets of online food images based on the given food labels. Since not all retrieved images contain foods, we introduce an automatic food detection method to remove irrelevant images. We designed a web-based crowdsourcing tool to assist the crowd or human annotators to locate and label all the foods in the images. The proposed semi-automatic online food image collection system can be used to build large food image datasets with groundtruth labels efficiently from scratch.


title: As You Are, So Shall You Move Your Head: A System-Level Analysis between Head Movements and Corresponding Traits and Emotions

abstract: Identifying physical traits and emotions based on system-sensed physical activities is a challenging problem in the realm of human-computer interaction. Our work contributes in this context by investigating an underlying connection between head movements and corresponding traits and emotions. To do so, we utilize a head movement measuring device called eSense, which gives acceleration and rotation of a head. Here, first, we conduct a thorough study over head movement data collected from 46 persons using eSense while inducing five different emotional states over them in isolation. Our analysis reveals several new head movement based findings, which in turn, leads us to a novel unified solution for identifying different human traits and emotions through exploiting machine learning techniques over head movement data. Our analysis confirms that the proposed solution can result in high accuracy over the collected data. Accordingly, we develop an integrated unified solution for real-time emotion and trait identification using head movement data leveraging outcomes of our analysis.


title: Intracavity-squeezed optomechanical cooling

abstract: Quantum ground-state cooling of macroscopic mechanical resonators is of essential importance to both fundamental physics and applied science. Conventional method of laser cooling is limited by the quantum backaction, which requires mechanical sideband resolved in order to cool to ground state. This work presents an idea to break the quantum backaction limit by engineering intracavity optical squeezing. It gives rise to quantum interference for all the dissipation channels, and under certain circumstances can totally remove the influence of the cavity dissipation and the resultant quantum backaction, with much lower cooling limit irrespective of the sideband resolution. We show that our scheme enables ground-state cooling in the highly unresolved sideband limit and it also works beyond the weak coupling regime, which provides the opportunity for quantum manipulation of macroscopic mechanical systems.


title: Decoupling Hierarchical Recurrent Neural Networks With Locally Computable Losses

abstract: Learning long-term dependencies is a key long-standing challenge of recurrent neural networks (RNNs). Hierarchical recurrent neural networks (HRNNs) have been considered a promising approach as long-term dependencies are resolved through shortcuts up and down the hierarchy. Yet, the memory requirements of Truncated Backpropagation Through Time (TBPTT) still prevent training them on very long sequences. In this paper, we empirically show that in (deep) HRNNs, propagating gradients back from higher to lower levels can be replaced by locally computable losses, without harming the learning capability of the network, over a wide range of tasks. This decoupling by local losses reduces the memory requirements of training by a factor exponential in the depth of the hierarchy in comparison to standard TBPTT.


title: Strongly Convex Optimization for Joint Fractal Feature Estimation and Texture Segmentation

abstract: The present work investigates the segmentation of textures by formulating it as a strongly convex optimization problem, aiming to favor piecewise constancy of fractal features (local variance and local regularity) widely used to model real-world textures in numerous applications very different in nature. Two objective functions combining these two features are compared, referred to as joint and coupled, promoting either independent or co-localized changes in local variance and regularity. To solve the resulting convex nonsmooth optimization problems, because the processing of large size images and databases are targeted, two categories of proximal algorithms (dual forward-backward and primal-dual), are devised and compared. An in-depth study of the objective functions, notably of their strong convexity, memory and computational costs, permits to propose significantly accelerated algorithms. A class of synthetic models of piecewise fractal texture is constructed and studied. They enable, by means of large-scale Monte-Carlo simulations, to quantify the benefits in texture segmentation of combining local regularity and local variance (as opposed to regularity only) while using strong-convexity accelerated primal-dual algorithms. Achieved results also permit to discuss the gains/costs in imposing co-localizations of changes in local regularity and local variance in the problem formulation. Finally, the potential of the proposed approaches is illustrated on real-world textures taken from a publicly available and documented database.


title: Magnetization governed magnetoresistance anisotropy in topological semimetal CeBi

abstract: Magnetic topological semimetals, the latest member of topological quantum materials, are attracting extensive attention as they may lead to topologically-driven spintronics. Currently, magnetotransport investigations on these materials are focused on anomalous Hall effect. Here, we report on the magnetoresistance anisotropy of topological semimetal CeBi, which has tunable magnetic structures arising from localized Ce 4f electrons and exhibits both negative and positive magnetoresistances, depending on the temperature. We found that the angle dependence of the negative magnetoresistance, regardless of its large variation with the magnitude of the magnetic field and with temperature, is solely dictated by the field-induced magnetization that is orientated along a primary crystalline axis and flops under the influence of a rotating magnetic field. The results reveal the strong interaction between conduction electrons and magnetization in CeBi. They also indicate that magnetoresistance anisotropy can be used to uncover the magnetic behavior and the correlation between transport phenomena and magnetism in magnetic topological semimetals.


title: Vector bundles of non-negative curvature over cohomogeneity one manifolds

abstract: We provide several results on the existence of metrics of non-negative sectional curvature on vector bundles over certain cohomogeneity one manifolds and homogeneous spaces up to suitable stabilization. Beside explicit constructions of the metrics, this is achieved by identifying equivariant structures upon these vector bundles via a comparison of their equivariant and non-equivariant K-theory. For this, in particular, we transcribe equivariant K-theory to equivariant rational cohomology and investigate surjectivity properties of induced maps in the Borel fibration via rational homotopy theory.


title: Extending powers of pseudo-Anosovs

abstract: Biringer, Johnson, and Minsky showed that a pseudo-Anosov map on a boundary component of an irreducible 3-manifold has a power that partially extends to the interior if and only if the (un)stable laminations of $f$ is an $\mathbb{R}$-projective limit of meridians. We prove that the power required for a pseudo-Anosov map to partially extend is not universally bounded. We construct a family of pseudo-Anosov maps $f_i$ for all $i=1,2...$ on a boundary component of a family of irreducible 3-manifolds $M_i$ such that $f_i^i$ partially extends to the interior of $M_i$ but $f_i^j$ does not for $j<i$.


title: Efficient and Adaptive Kernelization for Nonlinear Max-margin Multi-view Learning

abstract: Existing multi-view learning methods based on kernel function either require the user to select and tune a single predefined kernel or have to compute and store many Gram matrices to perform multiple kernel learning. Apart from the huge consumption of manpower, computation and memory resources, most of these models seek point estimation of their parameters, and are prone to overfitting to small training data. This paper presents an adaptive kernel nonlinear max-margin multi-view learning model under the Bayesian framework. Specifically, we regularize the posterior of an efficient multi-view latent variable model by explicitly mapping the latent representations extracted from multiple data views to a random Fourier feature space where max-margin classification constraints are imposed. Assuming these random features are drawn from Dirichlet process Gaussian mixtures, we can adaptively learn shift-invariant kernels from data according to Bochners theorem. For inference, we employ the data augmentation idea for hinge loss, and design an efficient gradient-based MCMC sampler in the augmented space. Having no need to compute the Gram matrix, our algorithm scales linearly with the size of training set. Extensive experiments on real-world datasets demonstrate that our method has superior performance.


title: The moment problem on curves with bumps

abstract: The power moments of a positive measure on the real line or the circle are characterized by the non-negativity of an infinite matrix, Hankel, respectively Toeplitz, attached to the data. Except some fortunate configurations, in higher dimensions there are no non-negativity criteria for the power moments of a measure to be supported by a prescribed closed set. We combine two well studied fortunate situations, specifically a class of curves in two dimensions classified by Scheiderer and Plaumann, and compact, basic semi-algebraic sets, with the aim at enlarging the realm of geometric shapes on which the power moment problem is accessible and solvable by non-negativity certificates.


title: Comet C/2017 S3 (PanSTARRS): Outbursts and Disintegration

abstract: The Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha camera on the SOlar and Heliospheric Observer (SOHO) satellite observed the hydrogen coma of comet C/2017 S3 (PanSTARRS) for the last month of its activity from 2018 July 4 to August 4 and what appears to have been its final disintegration just 11 days before its perihelion on August 15. The hydrogen coma indicated water production had a small outburst on July 8 at a heliocentric distance of 1.1AU and then a much larger one on July 20 at 0.8 AU. Over the following two weeks the water production dropped by more than a factor of ten after which it was no longer detectable. The behavior is reminiscent of comet C/1999 S4 (LINEAR) in 2000, which had a few small outbursts on its inbound orbit and a major outburst at a heliocentric distance of about 0.8 AU, which was close to its perihelion, followed by its complete disintegration that was documented by several sets of observations including SWAN. C/2017 S3 (PanSTARRS) however had a much larger water production rate than C/ 1999 S4 (LINEAR). Here we estimate the size of the nucleus of C/2017 S3 just before its final outburst and apparent disintegration was estimated using the total amount of water produced during its last weeks for a range of values of the refractory/ice ratio in the nucleus. We also determine the size distribution of the disintegrating particles as the comet faded.


title: Adversarial Colorization Of Icons Based On Structure And Color Conditions

abstract: We present a system to help designers create icons that are widely used in banners, signboards, billboards, homepages, and mobile apps. Designers are tasked with drawing contours, whereas our system colorizes contours in different styles. This goal is achieved by training a dual conditional generative adversarial network (GAN) on our collected icon dataset. One condition requires the generated image and the drawn contour to possess a similar contour, while the other anticipates the image and the referenced icon to be similar in color style. Accordingly, the generator takes a contour image and a man-made icon image to colorize the contour, and then the discriminators determine whether the result fulfills the two conditions. The trained network is able to colorize icons demanded by designers and greatly reduces their workload. For the evaluation, we compared our dual conditional GAN to several state-of-the-art techniques. Experiment results demonstrate that our network is over the previous networks. Finally, we will provide the source code, icon dataset, and trained network for public use.


title: On the Price of Independence for Vertex Cover, Feedback Vertex Set and Odd Cycle Transversal

abstract: Let $vc(G)$, $fvs(G)$ and $oct(G)$, respectively, denote the size of a minimum vertex cover, minimum feedback vertex set and minimum odd cycle transversal in a graph $G$. One can ask, when looking for these sets in a graph, how much bigger might they be if we require that they are independent; that is, what is the price of independence? If $G$ has a vertex cover, feedback vertex set or odd cycle transversal that is an independent set, then we let $ivc(G)$, $ifvs(G)$ or $ioct(G)$, respectively, denote the minimum size of such a set. Similar to a recent study on the price of connectivity (Hartinger et al. EuJC 2016), we investigate for which graphs $H$ the values of $ivc(G)$, $ifvs(G)$ and $ioct(G)$ are bounded in terms of $vc(G)$, $fvs(G)$ and $oct(G)$, respectively, when the graph $G$ belongs to the class of $H$-free graphs. We find complete classifications for vertex cover and feedback vertex set and an almost complete classification for odd cycle transversal (subject to three non-equivalent open cases). We also investigate for which graphs $H$ the values of $ivc(G)$, $ifvs(G)$ and $ioct(G)$ are equal to $vc(G)$, $fvs(G)$ and $oct(G)$, respectively, when the graph $G$ belongs to the class of $H$-free graphs. We find a complete classification for vertex cover and almost complete classifications for feedback vertex set (subject to one open case) and odd cycle transversal (subject to three open cases).


title: Quasar and galaxy classification in Gaia Data Release 2

abstract: We construct a supervised classifier based on Gaussian Mixture Models to probabilistically classify objects in Gaia data release 2 (GDR2) using only photometric and astrometric data in that release. The model is trained empirically to classify objects into three classes -- star, quasar, galaxy -- for G<=14.5 mag down to the Gaia magnitude limit of G=21.0 mag. Galaxies and quasars are identified for the training set by a cross-match to objects with spectroscopic classifications from the Sloan Digital Sky Survey. Stars are defined directly from GDR2. When allowing for the expectation that quasars are 500 times rarer than stars, and galaxies 7500 times rarer than stars (the class imbalance problem), samples classified with a threshold probability of 0.5 are predicted to have purities of 0.43 for quasars and 0.28 for galaxies, and completenesses of 0.58 and 0.72 respectively. The purities can be increased up to 0.60 by adopting a higher threshold. Not accounting for this expected low frequency of extragalactic objects (the class prior) would give both erroneously optimistic performance predictions and severely impure samples. Applying our model to all 1.20 billion objects in GDR2 with the required features, we classify 2.3 million objects as quasars and 0.37 million objects as galaxies (with individual probabilities above 0.5). The small number of galaxies is due to the strong bias of the satellite detection algorithm and on-ground data selection against extended objects. We infer the true number of quasars and galaxies -- as these classes are defined by our training set -- to be 690,000 and 110,000 respectively (+/- 50%). The aim of this work is to see how well extragalactic objects can be classified using only GDR2 data. Better classifications should be possible with the low resolution spectroscopy (BP/RP) planned for GDR3.


title: Vector-borne epidemics driven by human mobility

abstract: Vector-borne epidemics are the result of the combination of different factors such as the crossed contagions between humans and vectors, their demographic distribution and human mobility among others. The current availability of information about the former ingredients demands their incorporation to current mathematical models for vector-borne disease transmission. Here, relying on metapopulation dynamics, we propose a framework whose results are in fair agreement with those obtained from mechanistic simulations. This framework allows us to derive an expression of the epidemic threshold capturing with high accuracy the conditions leading to the onset of epidemics. Driven by these insights, we obtain a prevalence indicator to rank the patches according to the risk of being affected by a vector-borne disease. We illustrate the utility of this epidemic risk indicator by reproducing the spatial distribution Dengue cases reported in the city of Santiago de Cali (Colombia) from 2015 to 2016.


title: Direct measurements of the properties of Thick-GEM reflective photocathodes

abstract: In the context of the development of novel Thick GEM based detectors of single photons, the high resolution optical system, nicknamed Leopard, providing a detailed surface scanning of the Thick GEM electron multipliers, has been used for a set of systematic measurements of key Thick GEM properties. These results are reported and discussed. They confirm by direct observation Thick GEM properties previously inferred by indirect measurements and answer to relevant questions related to the use of Thick GEMs as photocathode substrates in novel gaseous photon detectors.


title: The chaotic wind of WR 40 as probed by BRITE

abstract: Among Wolf-Rayet stars, those of subtype WN8 are the intrinsically most variable. We have explored the long-term photometric variability of the brightest known WN8 star, WR 40, through four contiguous months of time-resolved, single-passband optical photometry with the BRIght Target Explorer (BRITE) nanosatellite mission. The Fourier transform of the observed light-curve reveals that the strong light variability exhibited by WR 40 is dominated by many randomly-triggered, transient, low-frequency signals. We establish a model in which the whole wind consists of stochastic clumps following an outflow visibility promptly rising to peak brightness upon clump emergence from the optically thick pseudo-photosphere in the wind, followed by a gradual decay according to the right-half of a Gaussian. Free electrons in each clump scatter continuum light from the star. We explore a scenario where the clump size follows a power-law distribution, and another one with an ensemble of clumps of constant size. Both scenarios yield simulated light curves morphologically resembling the observed light curve remarkably well, indicating that one cannot uniquely constrain the details of clump size distribution with only a photometric light curve. Nevertheless, independent evidence favours a negative-index power law, as seen in many other astrophysical turbulent media.


title: Improved Material Decomposition with a Two-step Regularization for spectral CT

abstract: One of the advantages of spectral computed tomography (CT) is it can achieve accurate material components using the material decomposition methods. The image-based material decomposition is a common method to obtain specific material components, and it can be divided into two steps: image reconstruction and post material decomposition. To obtain accurate material maps, the image reconstruction method mainly focuses on improving image quality by incorporating regularization priors. Very recently, the regularization priors are introduced into the post material decomposition procedure in the iterative image-based methods. Since the regularization priors can be incorporated into image reconstruction and post image-domain material decomposition, the performance of regularization by combining these two cases is still an open problem. To realize this goal, the material accuracy from those steps are first analyzed and compared. Then, to further improve the accuracy of decomposition materials, a two-step regularization based method is developed by incorporating priors into image reconstruction and post material decomposition. Both numerical simulation and preclinical mouse experiments are performed to demonstrate the advantages of the two-step regularization based method in improving material accuracy.


title: Structures of solid hydrogen at 300K

abstract: We present results predicting experimentally measurable structural quantities from molecular dynamics studies of hydrogen. In doing this, we propose a paradigm shift for experimentalists -- that the predictions from such calculations should be seen as the most likely hypotheses. Specifically, the experimental results should be aiming to distinguish between the candidate low-energy structures, rather than aiming to solve the simplest structure consistent with the data. We show that the room temperature X-ray diffraction patterns for hydrogen phases I, III, IV and V are very similar, with only small peaks denoting symmetry-breaking from the hcp Phase I. Because they incorporate atomic displacements the XRD patterns implied by molecular dynamics calculations are very different from those arising from the static minimum enthalpy structures found by structure searching. Simulations also show that within Phase I the molecular becomes increasingly confined to the basal plane and suggest the possibility of an unusual critical point terminating the Phase I-III boundary line.


title: Finite element approximation of Lyapunov equations related to parabolic stochastic PDEs

abstract: A numerical analysis for the fully discrete approximation of an operator Lyapunov equation related to linear SPDEs (stochastic partial differential equations) driven by multiplicative noise is considered. The discretization of the Lyapunov equation in space is given by finite elements and in time by a semiimplicit Euler scheme. The main result is the derivation of the rate of convergence in operator norm. Moreover, it is shown that the solution of the equation provides a representation of a quadratic and path dependent functional of the SPDE solution. This fact yields a deterministic numerical method to compute such functionals. As a secondary result, weak error rates are established for a fully discrete finite element approximation of the SPDE with respect to this functional. This is obtained as a consequence of the approximation analysis of the Lyapunov equation. It is the first weak convergence analysis for fully discrete finite element approximations of SPDEs driven by multiplicative noise that obtains double the strong rate of convergence, especially for path dependent functionals and smooth spatial noise. Numerical experiments illustrate the results empirically and it is demonstrated that the deterministic method has advantages over Monte Carlo sampling in a stability context.


title: Hear "No Evil", See "Kenansville": Efficient and Transferable Black-Box Attacks on Speech Recognition and Voice Identification Systems

abstract: Automatic speech recognition and voice identification systems are being deployed in a wide array of applications, from providing control mechanisms to devices lacking traditional interfaces, to the automatic transcription of conversations and authentication of users. Many of these applications have significant security and privacy considerations. We develop attacks that force mistranscription and misidentification in state of the art systems, with minimal impact on human comprehension. Processing pipelines for modern systems are comprised of signal preprocessing and feature extraction steps, whose output is fed to a machine-learned model. Prior work has focused on the models, using white-box knowledge to tailor model-specific attacks. We focus on the pipeline stages before the models, which (unlike the models) are quite similar across systems. As such, our attacks are black-box and transferable, and demonstrably achieve mistranscription and misidentification rates as high as 100% by modifying only a few frames of audio. We perform a study via Amazon Mechanical Turk demonstrating that there is no statistically significant difference between human perception of regular and perturbed audio. Our findings suggest that models may learn aspects of speech that are generally not perceived by human subjects, but that are crucial for model accuracy. We also find that certain English language phonemes (in particular, vowels) are significantly more susceptible to our attack. We show that the attacks are effective when mounted over cellular networks, where signals are subject to degradation due to transcoding, jitter, and packet loss.


title: Contextualising and Aligning Security Metrics and Business Objectives: a GQM-based Methodology

abstract: Pre-defined security metrics suffer from the problem of contextualisation, i.e. a lack of adaptability to particular organisational contexts - domain, technical infrastructure, stakeholders, business process, etc. Adapting metrics to an organisational context is essential (1) for the metrics to align with business requirements (2) for decision makers to maintain relevant security goals based on measurements from the field. In this paper we propose SYMBIOSIS, a methodology that defines a goal elicitation and refinement process mapping business objectives to security measurement goals via the use of systematic templates that capture relevant context elements (business goals, purpose, stakeholders, system scope). The novel contribution of SYMBIOSIS is the well-defined process, which enforces that (1) metrics align with business objectives via a top-down derivation that refines top-level business objectives to a manageable granularity (2) the impact of metrics on business objectives is explicitly traced via a bottom-up feedback mechanism, allowing an incremental approach where feedback from metrics influences business goals, and vice-versa. In this paper, we discuss the findings from applying SYMBIOSIS to three case studies of known security incidents. Our analysis shows how the aforementioned pitfalls of security metrics development processes affected the outcome of these high-profile security incidents and how SYMBIOSIS addresses such issues.


title: The generalized Holmgren problem for elliptic equation with several singular coefficients

abstract: Recently found all the fundamental solutions of a multidimensional singular elliptic equation are expressed in terms of the well-known Lauricella hypergeometric function in many variables. In this paper, we find a unique solution of the generalized Holmgren problem for an elliptic equation with several singular coefficients in explicit form. When finding a solution, we use decomposition formulas and some adjacent relations for the Lauricella hypergeometric function in many variables.


title: Automating dynamic consent decisions for the processing of social media data in health research

abstract: Social media have become a rich source of data, particularly in health research. Yet, the use of such data raises significant ethical questions about the need for the informed consent of those being studied. Consent mechanisms, if even obtained, are typically broad and inflexible, or place a significant burden on the participant. Machine learning algorithms show much promise for facilitating a 'middle ground' approach: using trained models to predict and automate granular consent decisions. Such techniques, however, raise a myriad of follow-on ethical and technical considerations. In this paper, we present an exploratory user study (n = 67) in which we find that we can predict the appropriate flow of health-related social media data with reasonable accuracy, while minimising undesired data leaks. We then attempt to deconstruct the findings of this study, identifying and discussing a number of real-world implications if such a technique were put into practice.


title: Backpropagation Algorithms and Reservoir Computing in Recurrent Neural Networks for the Forecasting of Complex Spatiotemporal Dynamics

abstract: We examine the efficiency of Recurrent Neural Networks in forecasting the spatiotemporal dynamics of high dimensional and reduced order complex systems using Reservoir Computing (RC) and Backpropagation through time (BPTT) for gated network architectures. We highlight advantages and limitations of each method and discuss their implementation for parallel computing architectures. We quantify the relative prediction accuracy of these algorithms for the longterm forecasting of chaotic systems using as benchmarks the Lorenz-96 and the Kuramoto-Sivashinsky (KS) equations. We find that, when the full state dynamics are available for training, RC outperforms BPTT approaches in terms of predictive performance and in capturing of the long-term statistics, while at the same time requiring much less training time. However, in the case of reduced order data, large scale RC models can be unstable and more likely than the BPTT algorithms to diverge. In contrast, RNNs trained via BPTT show superior forecasting abilities and capture well the dynamics of reduced order systems. Furthermore, the present study quantifies for the first time the Lyapunov Spectrum of the KS equation with BPTT, achieving similar accuracy as RC. This study establishes that RNNs are a potent computational framework for the learning and forecasting of complex spatiotemporal systems.


title: On the Measurement of Sachs Form Factors in Processes without and with Proton Spin Flip

abstract: The physical meaning of the decomposition of the Rosenbluth formula into two terms containing only squares of Sachs form factors has been established. A new method has been proposed for their independent measurement in the $e \vec{p} \to e \vec{p}$ elastic process when the initial proton at rest is fully polarized along the direction of motion of the final proton.


title: Improving Gradient Estimation in Evolutionary Strategies With Past Descent Directions

abstract: Evolutionary Strategies (ES) are known to be an effective black-box optimization technique for deep neural networks when the true gradients cannot be computed, such as in Reinforcement Learning. We continue a recent line of research that uses surrogate gradients to improve the gradient estimation of ES. We propose a novel method to optimally incorporate surrogate gradient information. Our approach, unlike previous work, needs no information about the quality of the surrogate gradients and is always guaranteed to find a descent direction that is better than the surrogate gradient. This allows to iteratively use the previous gradient estimate as surrogate gradient for the current search point. We theoretically prove that this yields fast convergence to the true gradient for linear functions and show under simplifying assumptions that it significantly improves gradient estimates for general functions. Finally, we evaluate our approach empirically on MNIST and reinforcement learning tasks and show that it considerably improves the gradient estimation of ES at no extra computational cost.


title: The density of complex zeros of random sums

abstract: Let $\{\eta_{j}\}_{j = 0}^{N}$ be a sequence of independent, identically distributed random complex Gaussian variables, and let $\{f_{j} (z)\}_{j = 0}^{N}$ be a sequence of given analytic functions that are real-valued on the real number line. We prove an exact formula for the expected density of the distribution of complex zeros of the random equation $\sum_{j = 0}^{N} \eta_{j} f_{j} (z) = \mathbf{K}$, where $\mathbf{K} \in \mathds{C}$. The method of proof employs a formula for the expected absolute value of quadratic forms of Gaussian random variables. We then obtain the limiting behaviour of the density function as $N$ tends to infinity and provide numerical computations for the density function and empirical distributions for random sums with certain functions $f_{j} (z)$. Finally, we study the case when the $f_{j} (z)$ are polynomials orthogonal on the real line and the unit circle.


title: Fast and Bayes-consistent nearest neighbors

abstract: Research on nearest-neighbor methods tends to focus somewhat dichotomously either on the statistical or the computational aspects -- either on, say, Bayes consistency and rates of convergence or on techniques for speeding up the proximity search. This paper aims at bridging these realms: to reap the advantages of fast evaluation time while maintaining Bayes consistency, and further without sacrificing too much in the risk decay rate. We combine the locality-sensitive hashing (LSH) technique with a novel missing-mass argument to obtain a fast and Bayes-consistent classifier. Our algorithm's prediction runtime compares favorably against state of the art approximate NN methods, while maintaining Bayes-consistency and attaining rates comparable to minimax. On samples of size $n$ in $\R^d$, our pre-processing phase has runtime $O(d n \log n)$, while the evaluation phase has runtime $O(d\log n)$ per query point.


title: A Test for Shared Patterns in Cross-modal Brain Activation Analysis

abstract: Determining the extent to which different cognitive modalities (understood here as the set of cognitive processes underlying the elaboration of a stimulus by the brain) rely on overlapping neural representations is a fundamental issue in cognitive neuroscience. In the last decade, the identification of shared activity patterns has been mostly framed as a supervised learning problem. For instance, a classifier is trained to discriminate categories (e.g. faces vs. houses) in modality I (e.g. perception) and tested on the same categories in modality II (e.g. imagery). This type of analysis is often referred to as cross-modal decoding. In this paper we take a different approach and instead formulate the problem of assessing shared patterns across modalities within the framework of statistical hypothesis testing. We propose both an appropriate test statistic and a scheme based on permutation testing to compute the significance of this test while making only minimal distributional assumption. We denote this test cross-modal permutation test (CMPT). We also provide empirical evidence on synthetic datasets that our approach has greater statistical power than the cross-modal decoding method while maintaining low Type I errors (rejecting a true null hypothesis). We compare both approaches on an fMRI dataset with three different cognitive modalities (perception, imagery, visual search). Finally, we show how CMPT can be combined with Searchlight analysis to explore spatial distribution of shared activity patterns.


title: Counting independent dominating sets in linear polymers

abstract: A vertex subset $W\subseteq V$ of the graph $G = (V,E)$ is an independent dominating set if every vertex in $V\setminus W$ is adjacent to at least one vertex in $W$ and the vertices of $W$ are pairwise non-adjacent. We enumerate independent dominating sets in several classes of graphs made by a linear or cyclic concatenation of basic building blocks. Explicit recurrences are derived for the number of independent dominating sets of these kind of graphs. Generating functions for the number of independent dominating sets of triangular and squares cacti chain are also computed.


title: The Natural Selection of Conservative Science

abstract: Social epistemologists have argued that high risk, high reward science has an important role to play in scientific communities. Recently, though, it has also been argued that various scientific fields seem to be trending towards conservatism -- the increasing production of what Kuhn (1970) would have called `normal science'. This paper will explore a possible explanation for this sort of trend: that the process by which scientific research groups form, grow, and dissolve might be inherently hostile to high risk science. In particular, I employ a paradigm developed by Smaldino and McElreath (2016) that treats a scientific community as a population undergoing selection. As will become clear, perhaps counter-intuitively this sort of process in some ways promotes high risk, high reward science. But, as I will point out, high risk high reward science is, in general, the sort of thing that is hard to repeat. While more conservative scientists will be able to train students capable of continuing their successful projects, and so create thriving lineages, successful risky science may not be the sort of thing one can easily pass on. In such cases, the structure of scientific communities selects against high risk, high rewards projects. More generally, this paper makes clear that there are at least two processes to consider in thinking about how incentives shape scientific communities -- the process by which individual scientists make choices about their careers and research, and the selective process governing the formation of new research groups.


title: A Geometric Model of Opinion Polarization

abstract: We introduce a simple, geometric model of opinion polarization. It is a model of political persuasion, as well as marketing and advertising, utilizing social values. It focuses on the interplay between different topics and persuasion efforts. We demonstrate that societal opinion polarization often arises as an unintended byproduct of influencers attempting to promote a product or idea. We discuss a number of mechanisms for the emergence of polarization involving one or more influencers, sending messages strategically, heuristically, or randomly. We also examine some computational aspects of choosing the most effective means of influencing agents, and the effects of those strategic considerations on polarization.


title: Maximum conditional entropy Hamiltonian Monte Carlo sampler

abstract: The performance of Hamiltonian Monte Carlo (HMC) sampler depends critically on some algorithm parameters such as the total integration time and the numerical integration stepsize. The parameter tuning is particularly challenging when the mass matrix of the HMC sampler is adapted. We propose in this work a Kolmogorov-Sinai entropy (KSE) based design criterion to optimize these algorithm parameters, which can avoid some potential issues in the often used jumping-distance based measures. For near-Gaussian distributions, we are able to derive the optimal algorithm parameters with respect to the KSE criterion analytically. As a byproduct the KSE criterion also provides a theoretical justification for the need to adapt the mass matrix in HMC sampler. Based on the results, we propose an adaptive HMC algorithm, and we then demonstrate the performance of the proposed algorithm with numerical examples.


title: exBERT: A Visual Analysis Tool to Explore Learned Representations in Transformers Models

abstract: Large language models can produce powerful contextual representations that lead to improvements across many NLP tasks. Since these models are typically guided by a sequence of learned self attention mechanisms and may comprise undesired inductive biases, it is paramount to be able to explore what the attention has learned. While static analyses of these models lead to targeted insights, interactive tools are more dynamic and can help humans better gain an intuition for the model-internal reasoning process. We present exBERT, an interactive tool named after the popular BERT language model, that provides insights into the meaning of the contextual representations by matching a human-specified input to similar contexts in a large annotated dataset. By aggregating the annotations of the matching similar contexts, exBERT helps intuitively explain what each attention-head has learned.


title: Keeping an Eye on DBI: Power-counting for small-$c_s$ Cosmology

abstract: Inflationary mechanisms for generating primordial fluctuations ultimately compute them as the leading contributions in a derivative expansion, with corrections controlled by powers of derivatives like the Hubble scale over Planck mass: $H/M_p$. At face value this derivative expansion breaks down for models with a small sound speed, $c_s$, to the extent that $c_s \ll 1$ is obtained by having higher-derivative interactions like $\mathfrak{L}_{\rm eff} \sim (\partial \Phi)^4$ compete with lower-derivative propagation. This concern arises more generally for models whose lagrangian is given as a function $P(X)$ for $X = -\partial_\mu \Phi \partial^\mu \Phi$ --- including in particular DBI models for which $P(X) \propto \sqrt{1-kX}$ --- since these keep all orders in $\partial \Phi$ while dropping $\partial^n \Phi$ for $n > 1$. We here find a sensible power-counting scheme for DBI models that gives a controlled expansion in powers of three types of small parameters: $H/M_p$, slow-roll parameters (possibly) and $c_s \ll 1$. We do not find a similar expansion framework for generic small-$c_s$ or $P(X)$ models. Our power-counting result quantifies the theoretical error for any prediction (such as for inflationary correlation functions) by fixing the leading power of these small parameters that is dropped when not computing all graphs (such as by restricting to the classical approximation); a prerequisite for meaningful comparisons with observations. The new power-counting regime arises because small $c_s$ alters the kinematics of free fluctuations in a way that changes how interactions scale at low energies, in particular allowing $1-c_s$ to be larger than derivative-measuring quantities like $(H/M_p)^2$.


title: On the Automorphism Group of a Polynomial Differential Ring in Two Variables

abstract: We consider differential rings of the form (K[x; y];D), where K is an algebraically closed field of characteristic zero and D : K[x; y] \to K[x; y] is a K-derivation. We study the Automorphism Group of such a ring and give criteria for deciding whether that group is an algebraic group. In most cases, from that study we deduce a primary classification of this type of differential ring up to conjugation with a polynomial automorphism.


title: Dynamical coherence of partially hyperbolic diffeomorphisms on nilmanifolds isotopic to Anosov

abstract: The purpose of this article is to obtain dynamically coherence of partially hyperbolic diffeomorphisms in certain classes of Anosov diffeomorphisms on nilmanifolds, extending a result due to T. Fisher, R. Potrie and M. Sambarino [FPS] on the torus.


title: Augmented Hard Example Mining for Generalizable Person Re-Identification

abstract: Although the performance of person re-identification (Re-ID) has been much improved by using sophisticated training methods and large-scale labelled datasets, many existing methods make the impractical assumption that information of a target domain can be utilized during training. In practice, a Re-ID system often starts running as soon as it is deployed, hence training with data from a target domain is unrealistic. To make Re-ID systems more practical, methods have been proposed that achieve high performance without information of a target domain. However, they need cumbersome tuning for training and unusual operations for testing. In this paper, we propose augmented hard example mining, which can be easily integrated to a common Re-ID training process and can utilize sophisticated models without any network modification. The method discovers hard examples on the basis of classification probabilities, and to make the examples harder, various types of augmentation are applied to the examples. Among those examples, excessively augmented ones are eliminated by a classification based selection process. Extensive analysis shows that our method successfully selects effective examples and achieves state-of-the-art performance on publicly available benchmark datasets.


title: Scale invariant solids

abstract: Scale invariance (SI) can in principle be realized in the elastic response of solid materials. There are two basic options: that SI is a manifest symmetry or that it is spontaneously broken. The manifest case corresponds physically to the existence of a non-trivial infrared fixed point with phonons among its degrees of freedom. We use simple bottom-up AdS/CFT constructions to model this case. We characterize the types of possible elastic response and discuss how the sound speeds can be realistic, that is, sufficiently small compared to the speed of light. We also study the spontaneously broken case using Effective Field Theory (EFT) methods. We present a new one-parameter family of nontrivial EFTs that includes the previously known `conformal solid' as a particular case as well as others which display small sound speeds. We also point out that an emergent Lorentz invariance at low energies could affect by order-one factors the relation between sound speeds and elastic moduli.


title: NGTS and WASP photometric recovery of a single-transit candidate from TESS

abstract: The Transiting Exoplanet Survey Satellite (\tess) produces a large number of single-transit event candidates, since the mission monitors most stars for only $\sim$27\,days. Such candidates correspond to long-period planets or eclipsing binaries. Using the \tess\ Sector 1 full-frame images, we identified a 7750\,ppm single-transit event with a duration of 7\,hours around the moderately evolved F-dwarf star \tic\ (Tmag=10.23, \teff=6280$\pm{85}$\,K). Using archival WASP photometry we constrained the true orbital period to one of three possible values. We detected a subsequent transit-event with NGTS, which revealed the orbital period to be 38.20\,d. Radial velocity measurements from the CORALIE Spectrograph show the secondary object has a mass of $M_2$= $0.148\pm{0.003}$\,M$_{\odot}$, indicating this system is an F-M eclipsing binary. The radius of the M-dwarf companion is $R_2$ = $0.171\pm{0.003}$\,R$_{\odot}$, making this one of the most well characterised stars in this mass regime. We find that its radius is 2.3-$\sigma$ lower than expected from stellar evolution models.


title: Shape Constrained Network for Eye Segmentation in the Wild

abstract: Semantic segmentation of eyes has long been a vital pre-processing step in many biometric applications. Majority of the works focus only on high resolution eye images, while little has been done to segment the eyes from low quality images in the wild. However, this is a particularly interesting and meaningful topic, as eyes play a crucial role in conveying the emotional state and mental well-being of a person. In this work, we take two steps toward solving this problem: (1) We collect and annotate a challenging eye segmentation dataset containing 8882 eye patches from 4461 facial images of different resolutions, illumination conditions and head poses; (2) We develop a novel eye segmentation method, Shape Constrained Network (SCN), that incorporates shape prior into the segmentation network training procedure. Specifically, we learn the shape prior from our dataset using VAE-GAN, and leverage the pre-trained encoder and discriminator to regularise the training of SegNet. To improve the accuracy and quality of predicted masks, we replace the loss of SegNet with three new losses: Intersection-over-Union (IoU) loss, shape discriminator loss and shape embedding loss. Extensive experiments shows that our method outperforms state-of-the-art segmentation and landmark detection methods in terms of mean IoU (mIoU) accuracy and the quality of segmentation masks. The eye segmentation database is available at https://www.dropbox.com/s/yvveouvxsvti08x/Eye_Segmentation_Database.zip?dl=0.


title: The Pentagon as a Substructure Lattice of Models of Peano Arithmetic

abstract: Wilke proved in 1977 that every countable model ${\mathcal M}$ of Peano Arithmetic has an elementary end extension ${\mathcal N}$ such that the interstructure lattice Lt(${\mathcal N} / {\mathcal M}$) is the pentagon lattice ${\mathbf N}_5$. This theorem implies that every countable nonstandard $\mathcal M$ has an elementary cofinal extension such that Lt(${\mathcal N} / {\mathcal M}) \cong {\mathbf N}_5$. It is proved here that if ${\mathcal M} \prec {\mathcal N}$ and Lt(${\mathcal N} / {\mathcal M}) \cong {\mathbf N}_5$, then ${\mathcal N}$ is either an end or a cofinal extension of ${\mathcal M}$. In contrast, there are ${\mathcal M}^* \prec {\mathcal N}^* \models {\mathsf PA}^*$ such that Lt(${\mathcal N} / {\mathcal M}) \cong {\mathbf N}_5$ and ${\mathcal N}^*$ is neither an end nor a cofinal extension of ${\mathcal M}^*$.


title: EFT of Interacting Spin-2 Fields

abstract: We consider the effective field theory of multiple interacting massive spin-2 fields. We focus on the case where the interactions are chosen so that the cutoff is the highest possible, and highlight two distinct classes of theories. In the first class, the mass eigenstates only interact through potential operators that carry no derivatives in unitary gauge at leading order. In the second class, a specific kinetic mixing between the mass eigenstates is included non-linearly. Performing a decoupling and ADM analysis, we point out the existence of a ghost present at a low scale for the first class of interactions. For the second class of interactions where kinetic mixing is included, we derive the full $\Lambda_3$ decoupling limit and confirm the absence of any ghosts. Nevertheless both formulations can be used to consistently describe an EFT of interacting massive spin-2 fields which, for a suitable technically natural tuning of the EFT, have the same strong coupling scale $\Lambda_3$. We identify the generic form of EFT corrections in each case. By using Galileon Duality transformations for the specific case of two massive spin-2 fields with suitable couplings, the decoupling limit theory is shown to be a bi-Galileon.


title: Black hole spontaneous scalarisation with a positive cosmological constant

abstract: A scalar field non-minimally coupled to certain geometric [or matter] invariants which are sourced by [electro]vacuum black holes (BHs) may spontaneously grow around the latter, due to a tachyonic instability. This process is expected to lead to a new, dynamically preferred, equilibrium state: a scalarised BH. The most studied geometric [matter] source term for such spontaneous BH scalarisation is the Gauss-Bonnet quadratic curvature [Maxwell invariant]. This phenomenon has been mostly analysed for asymptotically flat spacetimes. Here we consider the impact of a positive cosmological constant, which introduces a cosmological horizon. The cosmological constant does not change the local conditions on the scalar coupling for a tachyonic instability of the scalar-free BHs to emerge. But it leaves a significant imprint on the possible new scalarised BHs. It is shown that no scalarised BH solutions exist, under a smoothness assumption, if the scalar field is confined between the BH and cosmological horizons. Admitting the scalar field can extend beyond the cosmological horizon, we construct new scalarised BHs. These are asymptotically de Sitter in the (matter) Einstein-Maxwell-scalar model, with only mild difference with respect to their asymptotically flat counterparts. But in the (geometric) extended-scalar-tensor-Gauss-Bonnet-scalar model, they have necessarily non-standard asymptotics, as the tachyonic instability dominates in the far field. This interpretation is supported by the analysis of a test tachyon on a de Sitter background.


title: Improvements of upper curvature bounds

abstract: We show that any space with a positive upper curvature bound has in a small neighborhood of any point a closely related metric with a negative upper curvature bound.


title: Smoothly bounded domains covering compact manifolds

abstract: We show that if a bounded domain in complex Euclidean space with $\mathcal{C}^{1,1}$ boundary covers a compact manifold, then the domain is biholomorphic to the unit ball.


title: Universal quantum logic in hot silicon qubits

abstract: Quantum computation requires many qubits that can be coherently controlled and coupled to each other. Qubits that are defined using lithographic techniques are often argued to be promising platforms for scalability, since they can be implemented using semiconductor fabrication technology. However, leading solid-state approaches function only at temperatures below 100 mK, where cooling power is extremely limited, and this severely impacts the perspective for practical quantum computation. Recent works on spins in silicon have shown steps towards a platform that can be operated at higher temperatures by demonstrating long spin lifetimes, gate-based spin readout, and coherent single-spin control, but the crucial two-qubit logic gate has been missing. Here we demonstrate that silicon quantum dots can have sufficient thermal robustness to enable the execution of a universal gate set above one Kelvin. We obtain single-qubit control via electron-spin-resonance (ESR) and readout using Pauli spin blockade. We show individual coherent control of two qubits and measure single-qubit fidelities up to 99.3 %. We demonstrate tunability of the exchange interaction between the two spins from 0.5 up to 18 MHz and use this to execute coherent two-qubit controlled rotations (CROT). The demonstration of `hot' and universal quantum logic in a semiconductor platform paves the way for quantum integrated circuits hosting the quantum hardware and their control circuitry all on the same chip, providing a scalable approach towards practical quantum information.


title: Medium-Induced Transverse Momentum Broadening via Forward Dijet Correlations

abstract: Dijet azimuthal angle correlation is arguably one of the most direct probes of the medium-induced broadening effects. The evidence for such broadening, however, is not yet clearly observed within the precision of current mid-rapidity measurements at RHIC and the LHC. We show that the dijet correlation in forward rapidity from the future LHC RUN3, aided by forward detector upgrades, can reveal this broadening thanks to the steeper jet spectra, suppressed vacuum radiations and lower underlying event background, with a sensitivity comparable to that of the future high-luminosity Au+Au run at RHIC. Dijet correlation measurements at the two facilities together can provide powerful constraints on the temperature dependence of medium transport properties.


title: The Emergence of Compositional Languages for Numeric Concepts Through Iterated Learning in Neural Agents

abstract: Since first introduced, computer simulation has been an increasingly important tool in evolutionary linguistics. Recently, with the development of deep learning techniques, research in grounded language learning has also started to focus on facilitating the emergence of compositional languages without pre-defined elementary linguistic knowledge. In this work, we explore the emergence of compositional languages for numeric concepts in multi-agent communication systems. We demonstrate that compositional language for encoding numeric concepts can emerge through iterated learning in populations of deep neural network agents. However, language properties greatly depend on the input representations given to agents. We found that compositional languages only emerge if they require less iterations to be fully learnt than other non-degenerate languages for agents on a given input representation.


title: Ultrafast creation of overlapping Rydberg electrons in an atomic BEC and Mott-insulator lattice

abstract: An array of ultracold atoms in an optical lattice (Mott insulator) excited to a state where single electron wave-functions spatially overlap would represent a new and ideal platform to simulate exotic electronic many-body phenomena in the condensed phase. However, this highly excited non-equilibrium system is expected to be so short-lived that it has eluded observation so far. Here, we demonstrate the first step toward its realization by exciting high-lying electronic (Rydberg) states of the atomic Mott insulator with a coherent ultrashort laser pulse. Beyond a threshold principal quantum number where Rydberg orbitals of neighboring lattice sites overlap with each other, the atoms efficiently undergo spontaneous Penning ionization resulting in a drastic change of ion-counting statistics, sharp increase of avalanche ionization and the formation of an ultracold plasma. These observations signal the actual creation of exotic electronic states with overlapping wave functions, which is further confirmed by a significant difference in ionization dynamics between a Bose-Einstein condensate and a Mott insulator.

categories: physics.atom-ph cond-mat.quant-gas physics.chem-ph quant-ph


title: Lusin and Suslin properties of function spaces

abstract: A topological space is $Suslin$ ($Lusin$) if it is a continuous (and bijective) image of a Polish space. For a Tychonoff space $X$ let $C_p(X)$, $C_k(X)$ and $C_{{\downarrow}F}(X)$ be the space of continuous real-valued functions on $X$, endowed with the topology of pointwise convergence, the compact-open topology, and the Fell hypograph topology, respectively. For a metrizable space $X$ we prove the equivalence of the following statements: (1) $X$ is $\sigma$-compact, (2) $C_p(X)$ is Suslin, (3) $C_k(X)$ is Suslin, (4) $C_{{\downarrow}F}(X)$ is Suslin, (5) $C_p(X)$ is Lusin, (6) $C_k(X)$ is Lusin, (7) $C_{{\downarrow}F}(X)$ is Lusin, (8) $C_p(X)$ is $F_\sigma$-Lusin, (9) $C_k(X)$ is $F_\sigma$-Lusin, (10) $C_{{\downarrow}F}(X)$ is $C_{\delta\sigma}$-Lusin. Also we construct an example of a sequential $\aleph_0$-space $X$ with a unique non-isolated point such that the function spaces $C_p(X)$, $C_k(X)$ and $C_{{\downarrow}F}(X)$ are not Suslin.


title: Topology change of levels sets in Morse theory

abstract: Classical Morse theory proceeds by considering sublevel sets $f^{-1}(-\infty, a]$ of a Morse function $f: M \to R$, where $M$ is a smooth finite-dimensional manifold. In this paper, we study the topology of the level sets $f^{-1}(a)$ and give conditions under which the topology of $f^{-1}(a)$ changes when passing a critical value. We show that for a general class of functions, which includes all exhaustive Morse function, the topology of a regular level $f^{-1}(a)$ always changes when passing a single critical point, unless the index of the critical point is half the dimension of the manifold $M$. When $f$ is a natural Hamiltonian on a cotangent bundle, we obtain more precise results in terms of the topology of the configuration space. (Counter-)examples and applications to celestial mechanics are also discussed.


title: Optimal Approximation of Doubly Stochastic Matrices

abstract: We consider the least-squares approximation of a matrix C in the set of doubly stochastic matrices with the same sparsity pattern as C. Our approach is based on applying the well-known Alternating Direction Method of Multipliers (ADMM) to a reformulation of the original problem. Our resulting algorithm requires an initial Cholesky factorization of a positive definite matrix that has the same sparsity pattern as C + I followed by simple iterations whose complexity is linear in the number of nonzeros in C, thus ensuring excellent scalability and speed. We demonstrate the advantages of our approach in a series of experiments on problems with up to 82 million nonzeros; these include normalizing large scale matrices arising from the 3D structure of the human genome, clustering applications, and the SuiteSparse matrix library. Overall, our experiments illustrate the outstanding scalability of our algorithm; matrices with millions of nonzeros can be approximated in a few seconds on modest desktop computing hardware.