title: Option market (in)efficiency and implied volatility dynamics after return jumps

abstract: In informationally efficient financial markets, option prices and this implied volatility should immediately be adjusted to new information that arrives along with a jump in underlying's return, whereas gradual changes in implied volatility would indicate market inefficiency. Using minute-by-minute data on S&P 500 index options, we provide evidence regarding delayed and gradual movements in implied volatility after the arrival of return jumps. These movements are directed and persistent, especially in the case of negative return jumps. Our results are significant when the implied volatilities are extracted from at-the-money options and out-of-the-money puts, while the implied volatility obtained from out-of-the-money calls converges to its new level immediately rather than gradually. Thus, our analysis reveals that the implied volatility smile is adjusted to jumps in underlying's return asymmetrically. Finally, it would be possible to have statistical arbitrage in zero-transaction-cost option markets, but under actual option price spreads, our results do not imply abnormal option returns.


title: HALO7D II: The Halo Velocity Ellipsoid and Velocity Anisotropy with Distant Main Sequence Stars

abstract: The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use of a Bayesian mixture modeling approach for creating the stationary reference frame of distant galaxies. Using the 3D kinematic HALO7D sample, we estimate the parameters of the halo velocity ellipsoid, $\langle v_{\phi} \rangle, \sigma_r, \sigma_{\phi}, \sigma_{\theta}$, and the velocity anisotropy $\beta$. Using the full HALO7D sample, we find $\beta=0.63 \pm 0.05$ at $\langle r \rangle =24$ kpc. We also estimate the ellipsoid parameters for our sample split into three apparent magnitude bins; the posterior medians for these estimates of $\beta$, while consistent with one another, increase as a function of mean sample distance. Finally, we estimate $\beta$ in each of the individual HALO7D fields. We find that the velocity anisotropy $\beta$ can vary from field to field, which suggests that the halo is not phase mixed at $\langle r \rangle =24$ kpc. We explore the $\beta$ variation across the skies of two stellar halos from the \textit{Latte} suite of FIRE-2 simulations, finding that both simulated galaxies show $\beta$ variation over a similar range to the variation observed across the four HALO7D fields. The accretion histories of the two simulated galaxies result in different $\beta$ variation patterns; spatially mapping $\beta$ is thus a way forward in characterizing the accretion history of the Galaxy.


title: Fast, High-Quality Dual-Arm Rearrangement in Synchronous, Monotone Tabletop Setups

abstract: Rearranging objects on a planar surface arises in a variety of robotic applications, such as product packaging. Using two arms can improve efficiency but introduces new computational challenges. This paper studies the structure of dual-arm rearrangement for synchronous, monotone tabletop setups and develops an optimal mixed integer model. It then describes an efficient and scalable algorithm, which first minimizes the cost of object transfers and then of moves between objects. This is motivated by the fact that, asymptotically, object transfers dominate the cost of solutions. Moreover, a lazy strategy minimizes the number of motion planning calls and results in significant speedups. Theoretical arguments support the benefits of using two arms and indicate that synchronous execution, in which the two arms perform together either transfers or moves, introduces only a small overhead. Experiments support these points and show that the scalable method can quickly compute solutions close to the optimal for the considered setup.


title: $K$-theory of valuation rings

abstract: We prove several results showing that the algebraic $K$-theory of valuation rings behave as though such rings were regular Noetherian, in particular an analogue of the Geisser--Levine theorem. We also give some new proofs of known results concerning cdh descent of algebraic $K$-theory.


title: A Proper version of Synthesis-based Sparse Audio Declipper

abstract: Methods based on sparse representation have found great use in the recovery of audio signals degraded by clipping. The state of the art in declipping has been achieved by the SPADE algorithm by Kiti\'c et. al. (LVA/ICA2015). Our recent study (LVA/ICA2018) has shown that although the original S-SPADE can be improved such that it converges significantly faster than the A-SPADE, the restoration quality is significantly worse. In the present paper, we propose a new version of S-SPADE. Experiments show that the novel version of S-SPADE outperforms its old version in terms of restoration quality, and that it is comparable with the A-SPADE while being even slightly faster than A-SPADE.


title: Bounds on the first Betti number - an approach via Schatten norm estimates on semigroup differences

abstract: We derive new estimates for the first Betti number of compact Riemannian manifolds. Our approach relies on the Birman-Schwinger principle and Schatten norm estimates for semigroup differences. In contrast to previous works we do not require any a priori ultracontractivity estimates and we provide bounds which explicitly depend on suitable integral norms of the Ricci tensor.


title: HST/COS observations of the newly discovered obscuring outflow in NGC 3783

abstract: To understand the nature of transient obscuring outflows in active galactic nuclei, we observed the Seyfert 1 galaxy NGC 3783 on two occasions in December 2016 triggered by Swift monitoring indicating strong soft X-ray absorption in November. We obtained ultraviolet spectra using COS on HST and optical spectra using FEROS on the MPG/ESO 2.2-m telescope that were simultaneous with X-ray spectra from XMM-Newton and NuSTAR. We find new components of broad, blue-shifted absorption associated with Ly$\alpha$, \ion{N}{v}, \ion{Si}{iv}, and \ion{C}{iv} in our COS spectra. The absorption extends from velocities near zero in the rest-frame of the host galaxy to $-6200$ $\rm km~s^{-1}$. These features appear for the first time in NGC 3783 at the same time as the heavy soft X-ray absorption seen in the XMM-Newton X-ray spectra. The X-ray absorption has a column density of $\sim 10^{23}~\rm cm^{-2}$, and it partially covers the X-ray continuum source. The X-ray absorption becomes more transparent in the second observation, as does the UV absorption. Combining the X-ray column densities with the UV spectral observations yields an ionization parameter for the obscuring gas of log $\xi =1.84^{+0.4}_{-0.2}$ $\rm erg~cm~s^{-1}$. Despite the high intensity of the UV continuum in NGC 3783, F(1470 \AA)=$8 \times 10^{-14}~\rm erg~cm^{-2}~s^{-1}~\AA^{-1}$, the well known narrow UV absorption lines are deeper than in earlier observations in unobscured states, and low ionization species such as \ion{C}{iii} appear, indicating that the narrow-line gas is more distant from the nucleus and is being shadowed by the gas producing the obscuration. Despite the high continuum flux levels in our observations of NGC 3783, moderate velocities in the UV broad line profiles have substantially diminished. We suggest that a collapse of the broad line region has led to the outburst and triggered the obscuring event.


title: The determinants of academic career advancement: evidence from Italy

abstract: In this work we investigate the determinants of professors' career advancement in Italian universities. From the analyses, it emerges that the fundamental determinant of an academic candidate's success is not scientific merit, but rather the number of years that the candidate has belonged to the same university as the selection committee president. Where applicants have participated in research work with the president, their probability of success also increases significantly. The factors of the years of service and occurrence of joint research for the other commission members also have an effect, however of lesser weight. The specific phenomenon of nepotism, although it exists, seems less important. The scientific quality of the commission members has negligible effect on the expected outcome of the competition, and even more so the geographic location of the university calling for the competition.


title: Vortex modes supported by spin-orbit coupling in a laser with saturable absorption

abstract: We introduce a system of two component two-dimensional (2D) complex Ginzburg-Landau equations (CGLEs) with spin-orbit-coupling (SOC) describing a wide-aperture microcavity laser with saturable gain and absorption. We report families of two-component self-trapped dissipative laser solitons in this system. The SOC terms are represented by the second-order differential operators, which sets the difference, $|\Delta S|=2$, between the vorticities of the two components. We have found stable solitons of two types: vortex-antivortex (VAV) and semi-vortex (SV) bound states, featuring vorticities $\left( -1,+1\right) $ and $\left( 0,2\right) $, respectively. In previous works, 2D localized states of these types were found only in models including a trapping potential, while we are dealing with the self-trapping effect in the latteraly unconfined (free-space) model. The SV states are stable in a narrow interval of values of the gain coefficients. The stability interval is broader for VAV states, and it may be expanded by making SOC stronger (although the system without SOC features a stability interval too). We have found three branches of stationary solutions of both VAV and SV types, two unstable and one stable. The latter one is an attractor, as the unstable states spontaneously transform into the stable one, while retaining vorticities of their components. Unlike previously known 2D localized states, maintained by the combination of the trapping potential and SOC, in the present system the VAV and SV complexes are stable in the absence of diffusion. In contrast with the bright solitons in conservative models, chemical potentials of the dissipative solitons reported here are positive.


title: Quality-of-Service in Multihop Wireless Networks: Diffusion Approximation

abstract: We consider a multihop wireless system. There are multiple source-destination pairs. The data from a source may have to pass through multiple nodes. We obtain a channel scheduling policy which can guarantee end-to-end mean delay for the different traffic streams. We show the stability of the network for this policy by convergence to a fluid limit. It is intractable to obtain the stationary distribution of this network. Thus we also provide a diffusion approximation for this scheme under heavy traffic. We show that the stationary distribution of the scaled process of the network converges to that of the Brownian limit. This theoretically justifies the performance of the system. We provide simulations to verify our claims.


title: A Comparative Measurement Study of Deep Learning as a Service Framework

abstract: Big data powered Deep Learning (DL) and its applications have blossomed in recent years, fueled by three technological trends: a large amount of digitized data openly accessible, a growing number of DL software frameworks in open source and commercial markets, and a selection of affordable parallel computing hardware devices. However, no single DL framework, to date, dominates in terms of performance and accuracy even for baseline classification tasks on standard datasets, making the selection of a DL framework an overwhelming task. This paper takes a holistic approach to conduct empirical comparison and analysis of four representative DL frameworks with three unique contributions. First, given a selection of CPU-GPU configurations, we show that for a specific DL framework, different configurations of its hyper-parameters may have a significant impact on both performance and accuracy of DL applications. Second, to the best of our knowledge, this study is the first to identify the opportunities for improving the training time performance and the accuracy of DL frameworks by configuring parallel computing libraries and tuning individual and multiple hyper-parameters. Third, we also conduct a comparative measurement study on the resource consumption patterns of four DL frameworks and their performance and accuracy implications, including CPU and memory usage, and their correlations to varying settings of hyper-parameters under different configuration combinations of hardware, parallel computing libraries. We argue that this measurement study provides in-depth empirical comparison and analysis of four representative DL frameworks, and offers practical guidance for service providers to deploying and delivering DL as a Service (DLaaS) and for application developers and DLaaS consumers to select the right DL frameworks for the right DL workloads.


title: Photoionization cross sections of the 5S${}_{1/2}$ and 5P${}_{3/2}$ states of Rb in simultaneous magneto-optical trapping of Rb and Hg

abstract: We report the measurement of the photoionization cross sections of the 5S${}_{1/2}$ and 5P${}_{3/2}$ states of ${}^{87}$Rb in a two-species Hg and Rb magneto-optical trap (MOT) by the cooling laser for Hg. The photoionization cross sections of Rb in the 5S${}_{1/2}$ and 5P${}_{3/2}$ states at 253.7~nm are determined to be $1^{+4.3}_{-1}\times10^{-20}~\text{cm}^2$ and $4.63(30)\times 10^{-18}\text{cm}^2$, respectively. To measure the 5S${}_{1/2}$ and 5P${}_{3/2}$ states fractions in the MOT we detected photoionization rate of the 5P${}_{3/2}$ state by an additional 401.5~nm laser. The photoionization cross section of Rb in the 5P${}_{3/2}$ state at 401.5~nm is determined to be $\text{1.18(10)}\times10^{-17}~\text{cm}^2$.


title: Adapting a plant tissue model to animal development: introducing cell sliding into VirtualLeaf

abstract: Cell-based, mathematical modeling of collective cell behavior has become a prominent tool in developmental biology. Cell-based models represent individual cells as single particles or as sets of interconnected particles, and predict the collective cell behavior that follows from a set of interaction rules. In particular, vertex-based models are a popular tool for studying the mechanics of confluent, epithelial cell layers. They represent the junctions between three (or sometimes more) cells in confluent tissues as point particles, connected using structural elements that represent the cell boundaries. A disadvantage of these models is that cell-cell interfaces are represented as straight lines. This is a suitable simplification for epithelial tissues, where the interfaces are typically under tension, but this simplification may not be appropriate for mesenchymal tissues or tissues that are under compression, such that the cell-cell boundaries can buckle. In this paper we introduce a variant of VMs in which this and two other limitations of VMs have been resolved. The new model can also be seen as on off-the-lattice generalization of the Cellular Potts Model. It is an extension of the open-source package VirtualLeaf, which was initially developed to simulate plant tissue morphogenesis where cells do not move relative to one another. The present extension of VirtualLeaf introduces a new rule for cell-cell shear or sliding, from which T1 and T2 transitions emerge naturally, allowing application of VirtualLeaf to problems of animal development. We show that the updated VirtualLeaf yields different results than the traditional vertex-based models for differential-adhesion-driven cell sorting and for the neighborhood topology of soft cellular networks.


title: Strictification tensor product of 2-categories

abstract: Given 2-categories $\mathcal{C}$ and $\mathcal{D}$, let $\textrm{Lax}(\mathcal{C},\mathcal{D})$ denote the 2-category of lax functors, lax natural transformations and modifications, and $[\mathcal{C},\mathcal{D}]_\mathrm{lnt}$ its full sub-2-category of (strict) 2-functors. We give two isomorphic constructions of a 2-category $\mathcal{C}\boxtimes\mathcal{D}$ satisfying $\textrm{Lax}(\mathcal{C},\textrm{Lax}(\mathcal{D},\mathcal{E})) \cong [\mathcal{C}\boxtimes \mathcal{D},\mathcal{E}]_\mathrm{lnt}$, hence generalising the case of the free distributive law $1\boxtimes 1$. We also discuss dual constructions.


title: Lower central series, surface braid groups, surjections and permutations

abstract: Generalising previous results on classical braid groups by Artin and Lin, we determine the values of m, n $\in$ N for which there exists a surjection between the n-and m-string braid groups of an orientable surface without boundary. This result is essentially based on specific properties of their lower central series, and the proof is completely combinatorial. We provide similar but partial results in the case of orientable surfaces with boundary components and of non-orientable surfaces without boundary. We give also several results about the classification of different representations of surface braid groups in symmetric groups.


title: $^{7}$Li evolution in the thin and thick disks of the Milky Way

abstract: Recent detection of the isotope $^7$Be (that decays into $^7$Li) in the outburst of classical novae confirms the suggestion made in the 70s that novae could make $^7$Li. We reconsidered the role of novae as producers of Li by means of a detailed model of the chemical evolution of the Milky Way. We showed that novae could be the Galactic lithium source accounting for the observed increase of Li abundances in the thin disk with stellar metallicity. The best agreement with the upper envelope of the observed Li abundances is obtained for a delay time of ~ 1 Gyr for the nova production and for an effective Li yield of $1.8 (\pm 0.6) \times 10^{-5}$ $M_{\odot}$ over the whole nova life. Lithium in halo stars is depleted by ~ 0.35 dex, assuming the pristine abundance from standard big bang nucleosynthesis. We elaborate a model that matches the pristine stellar abundances, assuming that all stars are depleted by 0.35 dex. In this case, the delay remains the same, but the Li yields are doubled compared to the previous case. This model also has the merit to match the Li abundance in meteorites and young TTauri stars. The thick disk model, adopting the parameters derived for the thin disk, can also explain the absence of increase of Li abundances in its stars. The thick disk is old, but formed and evolved in a time shorter than that required by novae to contribute significantly to Li. Therefore, no Li enhancement is expected in thick disk stars. We show that the almost constant Li abundance in the thick disk results from the compensation of stellar astration by spallation processes.


title: On a class of spaces of skew-symmetric forms related to Hamiltonian systems of conservation laws

abstract: It was shown in \cite{FPV} that the classification of $n$-component systems of conservation laws possessing a third-order Hamiltonian structure reduces to the following algebraic problem: classify $n$-planes $H$ in $\wedge^2(V_{n+2})$ such that the induced map $Sym^2H\longrightarrow \wedge^4V_{n+2}$ has 1-dimensional kernel generated by a non-degenerate quadratic form on $H^*$. This problem is trivial for $n=2, 3$ and apparently wild for $n\geq 5$. In this paper we address the most interesting borderline case $n=4$. We prove that the variety $\mathcal{V}$ parametrizing those 4-planes $H$ is an irreducible 38-dimensional $PGL(V_6)$-invariant subvariety of the Grassmannian $G(4, \wedge^2V_6)$. With every $H\in\mathcal{V}$ we associate a {\it characteristic} cubic surface $S_H\subset \mathbb{P} H$, the locus of rank 4 two-forms in $ H$. We demonstrate that the induced characteristic map $\sigma: \mathcal{V} / PGL(V_6) \dashrightarrow \mathcal{M}_c,$ where $\mathcal{M}_c$ denotes the moduli space of cubic surfaces in $\mathbb{P}^3$, is dominant, hence generically finite. A complete classification of 4-planes $H\in\mathcal{V}$ with the reducible characteristic surface $S_H$ is given.


title: Dreaming neural networks: forgetting spurious memories and reinforcing pure ones

abstract: The standard Hopfield model for associative neural networks accounts for biological Hebbian learning and acts as the harmonic oscillator for pattern recognition, however its maximal storage capacity is $\alpha \sim 0.14$, far from the theoretical bound for symmetric networks, i.e. $\alpha =1$. Inspired by sleeping and dreaming mechanisms in mammal brains, we propose an extension of this model displaying the standard on-line (awake) learning mechanism (that allows the storage of external information in terms of patterns) and an off-line (sleep) unlearning$\&$consolidating mechanism (that allows spurious-pattern removal and pure-pattern reinforcement): this obtained daily prescription is able to saturate the theoretical bound $\alpha=1$, remaining also extremely robust against thermal noise. Both neural and synaptic features are analyzed both analytically and numerically. In particular, beyond obtaining a phase diagram for neural dynamics, we focus on synaptic plasticity and we give explicit prescriptions on the temporal evolution of the synaptic matrix. We analytically prove that our algorithm makes the Hebbian kernel converge with high probability to the projection matrix built over the pure stored patterns. Furthermore, we obtain a sharp and explicit estimate for the "sleep rate" in order to ensure such a convergence. Finally, we run extensive numerical simulations (mainly Monte Carlo sampling) to check the approximations underlying the analytical investigations (e.g., we developed the whole theory at the so called replica-symmetric level, as standard in the Amit-Gutfreund-Sompolinsky reference framework) and possible finite-size effects, finding overall full agreement with the theory.


title: The Primordial Black Hole Dark Matter - LISA Serendipity

abstract: There has recently been renewed interest in the possibility that the dark matter in the universe consists of primordial black holes (PBHs). Current observational constraints leave only a few PBH mass ranges for this possibility. One of them is around $10^{-12} M_\odot$. If PBHs with this mass are formed due to an enhanced scalar-perturbation amplitude, their formation is inevitably accompanied by the generation of gravitational waves (GWs) with frequency peaked in the mHz range, precisely around the maximum sensitivity of the LISA mission. We show that, if these primordial black holes are the dark matter, LISA will be able to detect the associated GW power spectrum. Although the GW source signal is intrinsically non-Gaussian, the signal measured by LISA is a sum of the signal from a large number of independent sources suppressing the non-Gaussianity at detection to an unobservable level. We also discuss the effect of the GW propagation in the perturbed universe. PBH dark matter generically leads to a detectable, purely isotropic, Gaussian and unpolarised GW signal, a prediction that is testable with LISA.


title: An Automated Singularity-Capturing Scheme for Fractional Differential Equations

abstract: Solutions to fractional models inherently exhibit non-smooth behavior, which significantly deteriorates the accuracy and therefore efficiency of existing numerical methods. We develop a two-stage data-infused computational framework for accurate time-integration of single- and multi-term fractional differential equations. In the first stage, we formulate a self-singularity-capturing scheme, given available/observable data for diminutive time. In this approach, the fractional differential equation provides the necessary knowledge/insight on how the hidden singularity can bridge between the initial and the subsequent short-time solution data. We develop a new self-singularity-capturing finite-difference algorithm for automatic determination of the underlying power-law singularities nearby the initial data, employing gradient descent optimization. In the second stage, we can utilize the multi-singular behavior of solution in a variety of numerical methods, without resorting to making any ad-hoc/uneducated guesses for the solution singularities. Particularly, we employed an implicit finite-difference method, where the captured singularities, in the first stage, are taken into account through some Lubich-like correction terms, leading to an accuracy of order $\mathcal{O}(\Delta t^{3-\alpha})$. Our computational results demonstrate that the developed framework can either fully capture or successfully control the solution error in the time-integration of fractional differential equations, especially in the presence of strong multi-singularities.


title: Machine learning uncovers new cosmological information

abstract: Large cosmological datasets have been probing the properties of our universe and constraining the parameters of dark matter and dark energy with increasing precision. Deep learning techniques have shown potential to be smarter, and to greatly outperform human-designed statistics.


title: A 3D Nonlinear Maxwell's Equations Solver Based On A Hybrid Numerical Method

abstract: In this paper we explore the possibility for solving the 3D Maxwell's equations in the presence of nonlinear and/or inhomogeneous material response. We propose using a hybrid approach which combines a bound- ary integral representation with a domain-based method. This hybrid approach has previously been successfully applied to 1D linear and non- linear transient wave scattering problems. The basic idea of the approach is to propagate the Maxwell's equations inside the scattering objects for- ward in time by using a domain-based method, while a boundary integral representation of the electromagnetic field is used to supply the domain- based method with the required surface values. Thus no grids outside the scattering objects are needed and this greatly reduces the computational cost and complexity.


title: A Preconditioned Multiple Shooting Shadowing Algorithm for the Sensitivity Analysis of Chaotic Systems

abstract: We propose a preconditioner that can accelerate the rate of convergence of the Multiple Shooting Shadowing (MSS) method. This recently proposed method can be used to compute derivatives of time-averaged objectives (also known as sensitivities) to system parameter(s) for chaotic systems. We propose a block diagonal preconditioner, which is based on a partial singular value decomposition of the MSS constraint matrix. The preconditioner can be computed using matrix-vector products only (i.e. it is matrix-free) and is fully parallelised in the time domain. Two chaotic systems are considered, the Lorenz system and the 1D Kuramoto Sivashinsky equation. Combination of the preconditioner with a regularisation method leads to tight bracketing of the eigenvalues to a narrow range. This combination results in a significant reduction in the number of iterations, and renders the convergence rate almost independent of the number of degrees of freedom of the system, and the length of the trajectory that is used to compute the time-averaged objective. This can potentially allow the method to be used for large chaotic systems (such as turbulent flows) and optimal control applications. The singular value decomposition of the constraint matrix can also be used to quantify the effect of the system condition on the accuracy of the sensitivities. In fact, neglecting the singular modes affected by noise, we recover the correct values of sensitivity that match closely with those obtained with finite differences for the Kuramoto Sivashinsky equation in the light turbulent regime.


title: A dynamical view on stellar metallicity gradient diversity across the Hubble sequence with CALIFA

abstract: We analyze radial stellar metallicity and kinematic profiles out to 1Re in 244 CALIFA galaxies ranging from morphological type E to Sd, to study the evolutionary mechanisms of stellar population gradients. We find that linear metallicity gradients exhibit a clear correlation with galaxy morphological type - with early type galaxies showing the steepest gradients. We show that the metallicity gradients simply reflect the local mass-metallicity relation within a galaxy. This suggests that the radial stellar population distribution within a galaxys effective radius is primarily a result of the \emph{in-situ} local star formation history. In this simple picture, the dynamically derived stellar surface mass density gradient directly predicts the metallicity gradient of a galaxy. We show that this correlation and its scatter can be reproduced entirely by using independent empirical galaxy structural and chemical scaling relations. Using Schwarzschild dynamical models, we also explore the link between a galaxys local stellar populations and their orbital structures. We find that galaxies angular momentum and metallicity gradients show no obvious causal link. This suggests that metallicity gradients in the inner disk are not strongly shaped by radial migration, which is confirmed by the lack of correlation between the metallicity gradients and observable probes of radial migration in the galaxies, such as bars and spiral arms. Finally, we find that galaxies with positive metallicity gradients become increasingly common towards low mass and late morphological types - consistent with stellar feedback more efficiently modifying the baryon cycle in the central regions of these galaxies.


title: Testing Primordial Black Holes as Dark Matter through LISA

abstract: The idea that primordial black holes (PBHs) can comprise most of the dark matter of the universe has recently reacquired a lot of momentum. Observational constraints, however, rule out this possibility for most of the PBH masses, with a notable exception around $10^{-12} M_\odot$. These light PBHs may be originated when a sizeable comoving curvature perturbation generated during inflation re-enters the horizon during the radiation phase. During such a stage, it is unavoidable that gravitational waves (GWs) are generated. Since their source is quadratic in the curvature perturbations, these GWs are generated fully non-Gaussian. Their frequency today is about the mHz, which is exactly the range where the LISA mission has the maximum of its sensitivity. This is certainly an impressive coincidence. We show that this scenario of PBHs as dark matter can be tested by LISA by measuring the GW two-point correlator. On the other hand, we show that the short observation time (as compared to the age of the universe) and propagation effects of the GWs across the perturbed universe from the production point to the LISA detector suppress the bispectrum to an unobservable level. This suppression is completely general and not specific to our model.


title: Strong regularization by Brownian noise propagating through a weak H{\"o}rmander structure

abstract: We establish strong uniqueness for a class of degenerate SDEs of weak H{\"o}rmander type under suitable H{\"o}lder regularity conditions for the associated drift term. Our approach relies on the Zvonkin transform which requires to exhibit good smoothing properties of the underlying parabolic PDE with rough, here H{\"o}lder, drift coefficients and source term. Such regularizing effects are established through a perturbation technique (forward parametrix approach) which also heavily relies on appropriate duality properties on Besov spaces. For the method employed, we exhibit some sharp thresholds on the H{\"o}lder exponents for the strong uniqueness to hold.


title: Suzuki functor at the critical level

abstract: In this paper we define and study a critical-level generalization of the Suzuki functor, relating the affine general linear Lie algebra to the rational Cherednik algebra of type A. Our main result states that this functor induces a surjective algebra homomorphism from the centre of the completed universal enveloping algebra at the critical level to the centre of the rational Cherednik algebra at t=0. We use this homomorphism to obtain several results about the functor. We compute it on Verma modules, Weyl modules, and their restricted versions. We describe the maps between endomorphism rings induced by the functor and deduce that every simple module over the rational Cherednik algebra lies in its image. Our homomorphism between the two centres gives rise to a closed embedding of the Calogero-Moser space into the space of opers on the punctured disc. We give a partial geometric description of this embedding.


title: Sharp Schauder Estimates for some Degenerate Kolmogorov Equations

abstract: We provide here some sharp Schauder estimates for degenerate PDEs of Kolmogorov type when the coefficients lie in some suitable anisotropic H{\"o}lder spaces and the first order term is non-linear and unbounded. We proceed through a perturbative approach based on forward parametrix expansions. Due to the low regularizing properties of the degenerate variables, for the procedure to work, we heavily exploit duality results between appropriate Besov spaces. Our method can be seen as constructive and provides, even in the non-degenerate case, an alternative approach to Schauder estimates.


title: Leveraging Gaussian Process and Voting-Empowered Many-Objective Evaluation for Fault Identification

abstract: Using piezoelectric impedance/admittance sensing for structural health monitoring is promising, owing to the simplicity in circuitry design as well as the high-frequency interrogation capability. The actual identification of fault location and severity using impedance/admittance measurements, nevertheless, remains to be an extremely challenging task. A first-principle based structural model using finite element discretization requires high dimensionality to characterize the high-frequency response. As such, direct inversion using the sensitivity matrix usually yields an under-determined problem. Alternatively, the identification problem may be cast into an optimization framework in which fault parameters are identified through repeated forward finite element analysis which however is oftentimes computationally prohibitive. This paper presents an efficient data-assisted optimization approach for fault identification without using finite element model iteratively. We formulate a many-objective optimization problem to identify fault parameters, where response surfaces of impedance measurements are constructed through Gaussian process-based calibration. To balance between solution diversity and convergence, an -dominance enabled many-objective simulated annealing algorithm is established. As multiple solutions are expected, a voting score calculation procedure is developed to further identify those solutions that yield better implications regarding structural health condition. The effectiveness of the proposed approach is demonstrated by systematic numerical and experimental case studies.


title: Supersolid behaviour of a dipolar Bose-Einstein condensate confined in a tube

abstract: Motivated by a recent experiment [L.Chomaz et al., Nature Physics 14, 442 (2018)], we perform numerical simulations of a dipolar Bose-Einstein Condensate (BEC) in a tubular confinement at T=0 within Density Functional Theory, where the beyond-mean-field correction to the ground state energy is included in the Local Density Approximation. We study the excitation spectrum of the system by solving the corresponding Bogoliubov-de Gennes equations. The calculated spectrum shows a roton minimum, and the roton gap decreases by reducing the effective scattering length. As the roton gap disappears, the system spontaneously develops in its ground-state a periodic, linear structure formed by denser clusters of atomic dipoles immersed in a dilute superfluid background. This structure shows the hallmarks of a supersolid system, i.e. (i) a finite non-classical translational inertia along the tube axis and (ii) the appearance, besides the phonon mode, of the Nambu-Goldstone gapless mode corresponding to phase fluctuations, and related to the spontaneous breaking of the gauge symmetry. A further decrease in the scattering length eventually leads to the formation of a periodic linear array of self-bound droplets.


title: A priori estimates for elliptic equations with reaction terms involving the function and its gradient

abstract: We study local and global properties of solutions of --$\Delta$u = u p + M ||u| q in a domain $\Omega$ of R N , in the range min{p, q} > 1 and M $\in$ R. We prove a priori estimates and existence or non-existence of ground states.


title: Optimal Output Regulation for Square, Over-Actuated and Under-Actuated Linear Systems

abstract: This paper considers two different problems in trajectory tracking control for linear systems. First, if the control is not unique which is most input energy efficient. Second, if exact tracking is infeasible which control performs most accurately. These are typical challenges for over-actuated systems and for under-actuated systems, respectively. We formulate both goals as optimal output regulation problems. Then we contribute two new sets of regulator equations to output regulation theory that provide the desired solutions. A thorough study indicates solvability and uniqueness under weak assumptions. E.g., we can always determine the solution of the classical regulator equations that is most input energy efficient. This is of great value if there are infinitely many solutions. We derive our results by a linear quadratic tracking approach and establish a useful link to output regulation theory.


title: Global null-controllability and nonnegative-controllability of slightly superlinear heat equations

abstract: We consider the semilinear heat equation posed on a smooth bounded domain $\Omega$ of $\mathbb{R}^{N}$ with Dirichlet or Neumann boundary conditions. The control input is a source term localized in some arbitrary nonempty open subset $\omega$ of $\Omega$. The goal of this paper is to prove the uniform large time global null-controllability for semilinearities $f(s) = \pm |s| \log^{\alpha}(2+|s|)$ where $\alpha \in [3/2,2)$ which is the case left open by Enrique Fernandez-Cara and Enrique Zuazua in 2000. It is worth mentioning that the free solution (without control) can blow-up. First, we establish the small-time global nonnegative-controllability (respectively nonpositive-controllability) of the system, i.e., one can steer any initial data to a nonnegative (respectively nonpositive) state in arbitrary time. In particular, one can act locally thanks to the control term in order to prevent the blow-up from happening. The proof relies on precise observability estimates for the linear heat equation with a bounded potential $a(t,x)$. More precisely, we show that observability holds with a sharp constant of the order $\exp\left(C |a|\_{\infty}^{1/2}\right)$ for nonnegative initial data. This inequality comes from a new $L^1$ Carleman estimate. A Kakutani's fixed point argument enables to go back to the semilinear heat equation. Secondly, the uniform large time null-controllability result comes from three ingredients: the global nonnegative-controllability, a comparison principle between the free solution and the solution to the underlying ordinary differential equation which provides the convergence of the free solution toward $0$ in $L^{\infty}(\Omega)$-norm, and the local null-controllability of the semilinear heat equation.


title: Approximate Bayesian Computation via Population Monte Carlo and Classification

abstract: Approximate Bayesian computation (ABC) methods can be used to sample from posterior distributions when the likelihood function is unavailable or intractable, as is often the case in biological systems. ABC methods suffer from inefficient particle proposals in high dimensions, and subjectivity in the choice of summary statistics, discrepancy measure, and error tolerance. Sequential Monte Carlo (SMC) methods have been combined with ABC to improve the efficiency of particle proposals, but suffer from subjectivity and require many simulations from the likelihood function. Likelihood-Free Inference by Ratio Estimation (LFIRE) leverages classification to estimate the posterior density directly but does not explore the parameter space efficiently. This work proposes a classification approach that approximates population Monte Carlo (PMC), where model class probabilities from classification are used to update particle weights. This approach, called Classification-PMC, blends adaptive proposals and classification, efficiently producing samples from the posterior without subjectivity. We show through a simulation study that Classification-PMC outperforms two state-of-the-art methods: ratio estimation and SMC ABC when it is computationally difficult to simulate from the likelihood.


title: Indirect Probe of Electroweak-Interacting Particles with Mono-Lepton Signatures at Hadron Colliders

abstract: Electroweak-interacting massive particles (EWIMPs) exist in a broad class of new physics models beyond the Standard Model. Searching for such particles is one of most primary goal at the LHC and future colliders. The EWIMP generally affects the LHC signatures through quantum corrections even without direct productions. By measuring the Standard Model processes precisely, we can indirectly probe the EWIMPs. In this paper, we study the current constraint and future prospect of the EWIMPs by using the precision measurements of mono-lepton production from the charged Drell-Yan processes at hadron colliders. We found the mono-lepton signature can be a better probe than dilepton signature from the neutral Drell-Yan processes.


title: High resolution ultra-sensitive trace gas detection by use of cavity-position-modulated sub-Doppler NICE-OHMS - Application to detection of acetylene in human breath

abstract: A sensitive high resolution sub-Doppler detecting spectrometer, based on noise-immune cavity-enhanced optical heterodyne molecular spectrometry (NICE-OHMS), for trace gas detection of species whose transitions have severe spectral overlap with abundant concomitant species is presented. It is designed around a NICE-OHMS instrumentation utilizing balanced detection that provides shot-noise limited Doppler-broadened detection. By synchronous dithering the positions of the two cavity mirrors, the effect of residual etalons between the cavity and other surface in the system could be reduced. An Allan deviation of the absorption coefficient of coefficient of $2.2 \times 10^{-13}$ $\text{cm}^{-1}$ at 60 s, which, for the targeted transition in $\text{C}_{2}\text{H}_{2}$, corresponds to a $3\sigma$ detection sensitivity of 130 ppt, is demonstrated. It is shown that despite significant spectral interference from $\text{CO}_{2}$ at the targeted transition, which precludes Db detection of $\text{C}_{2}\text{H}_{2}$, acetylene could be detected in exhaled breath of healthy smokers.


title: Equilibrium to Einstein: Entanglement, Thermodynamics, and Gravity

abstract: Here we develop the connection between thermodynamics, entanglement, and gravity. By attributing thermodynamics to timeslices of a causal diamond, we show that the Clausius relation $T\Delta S_{\text{rev}}=Q$, where $\Delta S_{\text{rev}}$ is the reversible entropy change, gives rise to the non-linear gravitational equations of motion for a wide class of diffeomorphism invariant theories. We then compare the Clausius relation to the first law of causal diamond mechanics (FLCD), a geometric identity and necessary ingredient in deriving Jacobson's entanglement equilibrium proposal -- the entanglement entropy of a spherical region with a fixed volume is maximal in vacuum. Specifically we show that the condition of fixed volume can be understood as subtracting the irreversible contribution to the thermodynamic entropy. This provides a "reversible thermodynamic process" interpretation of the FLCD, and that the condition of entanglement equilibrium may be regarded as equilibrium thermodynamics for which the Clausius relation holds. Finally, we extend the entanglement equilibrium proposal to the timelike stretched horizons of future lightcones, providing an entanglement interpretation of stretched lightcone thermodynamics.


title: Magneto-optical properties of Au upon the injection of hot spin-polarized electrons across Fe/Au(001) interfaces

abstract: We demonstrate a novel method for the excitation of sizable magneto-optical effects in Au by means of the laser-induced injection of hot spin-polarized electrons in Au/Fe/MgO(001) heterostructures. It is based on the energy- and spin-dependent electron transmittance of Fe/Au interface which acts as a spin filter for non-thermalized electrons optically excited in Fe. We show that after crossing the interface, majority electrons propagate through the Au layer with the velocity on the order of 1 nm/fs (close to the Fermi velocity) and the decay length on the order of 100 nm. Featuring ultrafast functionality and requiring no strong external magnetic fields, spin injection results in a distinct magneto-optical response of Au. We develop a formalism based on the phase of the transient complex MOKE response and demonstrate its robustness in a plethora of experimental and theoretical MOKE studies on Au, including our ab initio calculations. Our work introduces a flexible tool to manipulate magneto-optical properties of metals on the femtosecond timescale that holds high potential for active magneto-photonics, plasmonics, and spintronics.


title: Dark Matter Searches at Colliders

abstract: Colliders, among the most successful tools of particle physics, have revealed much about matter. This review describes how colliders contribute to the search for particle dark matter, focusing on the highest-energy collider currently in operation, the Large Hadron Collider (LHC) at CERN. In the absence of hints about the character of interactions between dark matter and standard matter, this review emphasizes what could be observed in the near future, presents the main experimental challenges, and discusses how collider searches fit into the broader field of dark matter searches. Finally, it highlights a few areas to watch for the future LHC program.


title: On the long time behavior of a tumor growth model

abstract: We consider the problem of the long time dynamics for a diffuse interface model for tumor growth. The model describes the growth of a tumor surrounded by host tissues in the presence of a nutrient and consists in a Cahn-Hilliard-type equation for the tumor phase coupled with a reaction-diffusion equation for the nutrient concentration. We prove that, under physically motivated assumptions on parameters and data, the corresponding initial-boundary value problem generates a dissipative dynamical system that admits the global attractor in a proper phase space.


title: Eco-PANDA: A Computationally Economic, Geometrically Converging, Dual Optimization Method on Time-Varying Undirected Graphs

abstract: In this paper we consider distributed convex optimization over time-varying undirected graphs. We propose a linearized version of primarily averaged network dual ascent (PANDA) while requiring less computational costs. The proposed method, economic primarily averaged network dual ascent (Eco-PANDA), provably converges at R-linear rate to the optimal point given that the agents' objective functions are strongly convex and have Lipschitz continuous gradients. Therefore, the method is competitive, in terms of type of rate, with both DIGing and PANDA. The proposed method halves the communication costs of methods like DIGing while still converging R-linearly and having the same per iterate complexity.


title: Few-shot 3D Multi-modal Medical Image Segmentation using Generative Adversarial Learning

abstract: We address the problem of segmenting 3D multi-modal medical images in scenarios where very few labeled examples are available for training. Leveraging the recent success of adversarial learning for semi-supervised segmentation, we propose a novel method based on Generative Adversarial Networks (GANs) to train a segmentation model with both labeled and unlabeled images. The proposed method prevents over-fitting by learning to discriminate between true and fake patches obtained by a generator network. Our work extends current adversarial learning approaches, which focus on 2D single-modality images, to the more challenging context of 3D volumes of multiple modalities. The proposed method is evaluated on the problem of segmenting brain MRI from the iSEG-2017 and MRBrainS 2013 datasets. Significant performance improvement is reported, compared to state-of-art segmentation networks trained in a fully-supervised manner. In addition, our work presents a comprehensive analysis of different GAN architectures for semi-supervised segmentation, showing recent techniques like feature matching to yield a higher performance than conventional adversarial training approaches. Our code is publicly available at https://github.com/arnab39/FewShot_GAN-Unet3D


title: The exact discontinuity of a partial wave along the left-hand cut and the exact $N/D$ method in non-relativistic scattering

abstract: We first deduce the analytical continuation in the complex planes of the initial and final three-momenta of the Lippmann-Schwinger equation in coupled or uncoupled partial-wave amplitudes. This result allows us to deduce a master equation whose solution is the exact discontinuity of the on-shell partial-wave amplitudes along the left-hand cut. This equation is always a linear non-singular integral equation whose solution is fixed exclusively by the knowledge of the potential, applicable to either regular or singular potentials. The capability of calculating exactly this discontinuity allows one to settle the exact $N/D $ method in two-body non-relativistic scattering for coupled and uncoupled waves. We exemplify this new advance in scattering theory by explicitly checking the agreement between the Lippmann-Schwinger equation with the corresponding solutions of the exact $N/D$ method for some examples that involve regular and singular potentials, either attractive or repulsive.

categories: hep-ph cond-mat.quant-gas hep-th nucl-th quant-ph


title: Seismicity in sheared granular matter

abstract: We report on experiments investigating the dynamics of a slider that is pulled by a spring across a granular medium consisting of a vertical layer of photo-elastic disks. The motion proceeds through a sequence of discrete events, analogous to seismic shocks, in which elastic energy stored in the spring is rapidly released. We measure the statistics of several properties of the individual events: the energy loss in the spring, the duration of the movement, and the temporal profile of the slider motion. We also study certain conditional probabilities and the statistics of mainshock-aftershock sequences. At low driving rates, we observe crackling with Omori-Utsu, Bath, and waiting time laws similar to those observed in seismic dynamics. At higher driving rates, where the sequence of events shows strong periodicity, we observe scaling laws and asymmetrical event shapes that are clearly distinguishable from those in the crackling regime.


title: MCell-R: A particle-resolution network-free spatial modeling framework

abstract: Spatial heterogeneity can have dramatic effects on the biochemical networks that drive cell regulation and decision-making. For this reason, a number of methods have been developed to model spatial heterogeneity and incorporated into widely used modeling platforms. Unfortunately, the standard approaches for specifying and simulating chemical reaction networks become untenable when dealing with multi-state, multi-component systems that are characterized by combinatorial complexity. To address this issue, we developed MCell-R, a framework that extends the particle-based spatial Monte Carlo simulator, MCell, with the rule-based model specification and simulation capabilities provided by BioNetGen and NFsim. The BioNetGen syntax enables the specification of biomolecules as structured objects whose components can have different internal states that represent such features as covalent modification and conformation and which can bind components of other molecules to form molecular complexes. The network-free simulation algorithm used by NFsim enables efficient simulation of rule-based models even when the size of the network implied by the biochemical rules is too large to enumerate explicitly, which frequently occurs in detailed models of biochemical signaling. The result is a framework that can efficiently simulate systems characterized by combinatorial complexity at the level of spatially-resolved individual molecules over biologically relevant time and length scales.


title: Prediction of Discretization of GMsFEM using Deep Learning

abstract: In this paper, we propose a deep-learning-based approach to a class of multiscale problems. THe Generalized Multiscale Finite Element Method (GMsFEM) has been proven successful as a model reduction technique of flow problems in heterogeneous and high-contrast porous media. The key ingredients of GMsFEM include mutlsicale basis functions and coarse-scale parameters, which are obtained from solving local problems in each coarse neighborhood. Given a fixed medium, these quantities are precomputed by solving local problems in an offline stage, and result in a reduced-order model. However, these quantities have to be re-computed in case of varying media. The objective of our work is to make use of deep learning techniques to mimic the nonlinear relation between the permeability field and the GMsFEM discretizations, and use neural networks to perform fast computation of GMsFEM ingredients repeatedly for a class of media. We provide numerical experiments to investigate the predictive power of neural networks and the usefulness of the resultant multiscale model in solving channelized porous media flow problems.


title: Reply to "Comment on 'Piezoelectricity in planar boron nitride via a geometric phase'"

abstract: In relation to our original paper [M. Droth et al., Phys. Rev. B 94, 075404 (2016)], the Comment by Li et al. [Phys. Rev. B 98, 167403 (2018)] claims to have identified a "mistake in constructing the adiabatic process of the piezoelectricity". More specifically, they write that in our original work "the erroneous usage of the polarization difference formula in Eq. (4) leads to an invalid analytical expression of piezoelectric constant in Eq. (12) of Reference 1". We explain below why these and other minor claims in the Comment are unwarranted, and why we maintain that our result is correct and physically sound.


title: Enabling Factorized Piano Music Modeling and Generation with the MAESTRO Dataset

abstract: Generating musical audio directly with neural networks is notoriously difficult because it requires coherently modeling structure at many different timescales. Fortunately, most music is also highly structured and can be represented as discrete note events played on musical instruments. Herein, we show that by using notes as an intermediate representation, we can train a suite of models capable of transcribing, composing, and synthesizing audio waveforms with coherent musical structure on timescales spanning six orders of magnitude (~0.1 ms to ~100 s), a process we call Wave2Midi2Wave. This large advance in the state of the art is enabled by our release of the new MAESTRO (MIDI and Audio Edited for Synchronous TRacks and Organization) dataset, composed of over 172 hours of virtuosic piano performances captured with fine alignment (~3 ms) between note labels and audio waveforms. The networks and the dataset together present a promising approach toward creating new expressive and interpretable neural models of music.


title: Electron Monte Carlo Simulations of Nanoporous Si Thin Films -- The Influence of Pore-Edge Charges

abstract: Electron transport within nanostructures can be important to varied engineering applications, such as thermoelectrics and nanoelectronics. In theoretical studies, electron Monte Carlo simulations are widely used as an alternative approach to solving the electron Boltzmann transport equation, where the energy-dependent electron scattering, exact structure shape, and detailed electric field distribution can be fully incorporated. In this work, such electron Monte Carlo simulations are employed to predict the electrical conductivity of periodic nanoporous Si films that have been widely studied for thermoelectric applications. The focus is on the influence of pore-edge charges on the electron transport. The results are further compared to our previous modeling [Hao et al., J. Appl. Phys. 121, 094308 (2017)], where the pore-edge electric field has its own scattering rate to be added to the scattering rates of other mechanisms.


title: Efficient Constrained Signal Reconstruction by Randomized Epigraphical Projection

abstract: This paper proposes a randomized optimization framework for constrained signal reconstruction, where the word "constrained" implies that data-fidelity is imposed as a hard constraint instead of adding a data-fidelity term to an objective function to be minimized. Such formulation facilitates the selection of regularization terms and hyperparameters, but due to the non-separability of the data-fidelity constraint, it does not suit block-coordinate-wise randomization as is. To resolve this, we give another expression of the data-fidelity constraint via epigraphs, which enables to design a randomized solver based on a stochastic proximal algorithm with randomized epigraphical projection. Our method is very efficient especially when the problem involves non-structured large matrices. We apply our method to CT image reconstruction, where the advantage of our method over the deterministic counterpart is demonstrated.


title: Chromatic homotopy is algebraic when $p > n^{2}+n+1$

abstract: We show that if $E$ is a $p$-local Landweber exact homology theory of height $n$ and $p > n^2+n+1$, then there exists an equivalence $h \mathcal{S}p_{E} \simeq h\mathcal{D}(E_{*}E)$ between homotopy categories of $E$-local spectra and differential $E_{*}E$-comodules, generalizing Bousfield's and Franke's results to heights $n > 1$.


title: ALMA detects a radial disk wind in DG Tau

abstract: Aims: We aim to use the high spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA) to map the flow pattern of molecular gas near DG Tau and its disk, a young stellar object driving a jet and a molecular outflow. Methods: We use observations from ALMA in the J = 2 - 1 transition of 12CO, 13CO, and C18O to study the Keplerian disk of DG Tau and outflows that may be related to the disk and the jet. Results: We find a new wind component flowing radially at a steep angle (~25 deg from the vertical) above the disk with a velocity of ~ 3.1 km/s. It continues the trend of decreasing velocity for increasing distance from the jet axis ("onion-like velocity structure"). Conclusions: The new component is located close to the protostellar disk surface and may be related to photoevaporative winds.


title: Linear functionals preserving the units of a Riesz space

abstract: Let $L$ be a Riesz space with a strong unit $e>0$.$\mathfrak{\ }$We show that a unital linear functional $H:L\rightarrow \mathbb{R}$ satisfies $% H\left( u\right) \neq 0$ for any strong unit $u\in L$ if and only if $H$ acts like a Riesz homomorphism on every $e$-clean vector subspace of $L.$ We deduce that $L$ is $e$-clean if and only if any unital linear functional $% H:L\rightarrow \mathbb{R}$ such that $H\left( u\right) \neq 0$ for any strong unit $u\in L$ is a Riesz homomorphism.


title: On the End-Vertex Problem of Graph Searches

abstract: End vertices of graph searches can exhibit strong structural properties and are crucial for many graph algorithms. The problem of deciding whether a given vertex of a graph is an end-vertex of a particular search was first introduced by Corneil, K\"ohler and Lanlignel in 2010. There they showed that this problem is in fact NP-complete for LBFS on weakly chordal graphs. A similar result for BFS was obtained by Charbit, Habib and Mamcarz in 2014. Here, we prove that the end-vertex problem is NP-complete for MNS on weakly chordal graphs and for MCS on general graphs. Moreover, building on previous results, we show that this problem is linear for various searches on split and unit interval graphs.


title: Particle physics experiments based on the AWAKE acceleration scheme

abstract: New particle acceleration schemes open up exciting opportunities, potentially providing more compact or higher-energy accelerators. The AWAKE experiment at CERN is currently taking data to establish the method of proton-driven plasma wakefield acceleration. A second phase aims to demonstrate that bunches of about $10^9$ electrons can be accelerated to high energy, preserving emittance and that the process is scalable with length. With this, an electron beam of $O$(50 GeV) could be available for new fixed-target or beam-dump experiments searching for the hidden sector, like dark photons. The rate of electrons on target could be increased by a factor of more than 1000 compared to currently available, leading to a corresponding increase in sensitivity to new physics. Such a beam could also be brought into collision with a high-power laser and thereby probe the completely unmeasured region of strong fields at values of the Schwinger critical field. An ultimate goal is to produce an electron beam of $O$(3 TeV) and collide with an LHC proton beam. This very high energy electron-proton collider would probe a new regime in which the structure of matter is completely unknown.


title: Low dose CT reconstruction assisted by an image manifold prior

abstract: X-ray Computed Tomography (CT) is an important tool in medical imaging to obtain a direct visualization of patient anatomy. However, the x-ray radiation exposure leads to the concern of lifetime cancer risk. Low-dose CT scan can reduce the radiation exposure to patient while the image quality is usually degraded due to the appearance of noise and artifacts. Numerous studies have been conducted to regularize CT image for better image quality. Yet, exploring the underlying manifold where real CT images residing on is still an open problem. In this paper, we propose a fully data-driven manifold learning approach by incorporating the emerging deep-learning technology. An encoder-decoder convolutional neural network has been established to map a CT image to the inherent low-dimensional manifold, as well as to restore the CT image from its corresponding manifold representation. A novel reconstruction algorithm assisted by the leant manifold prior has been developed to achieve high quality low-dose CT reconstruction. In order to demonstrate the effectiveness of the proposed framework, network training, testing, and comprehensive simulation study have been performed using patient abdomen CT images. The trained encoder-decoder CNN is capable of restoring high-quality CT images with average error of ~20 HU. Furthermore, the proposed manifold prior assisted reconstruction scheme achieves high-quality low-dose CT reconstruction, with average reconstruction error of < 30 HU, more than five times and two times lower than that of filtered back projection method and total-variation based iterative reconstruction method, respectively.


title: Semi-simple enlargement of the $\mathfrak{bms}_3$ algebra from a $\mathfrak{so}(2,2)\oplus\mathfrak{so}(2,1)$ Chern-Simons theory

abstract: In this work we present a BMS-like ansatz for a Chern-Simons theory based on the semi-simple enlargement of the Poincar\'e symmetry, also known as AdS-Lorentz algebra. We start by showing that this ansatz is general enough to contain all the relevant stationary solutions of this theory and provides with suitable boundary conditions for the corresponding gauge connection. We find an explicit realization of the asymptotic symmetry at null infinity, which defines a semi-simple enlargement of the $\mathfrak{bms}_3$ algebra and turns out to be isomorphic to three copies of the Virasoro algebra. The flat limit of the theory is discussed at the level of the action, field equations, solutions and asymptotic symmetry.


title: First Results from ABRACADABRA-10 cm: A Search for Sub-$\mu$eV Axion Dark Matter

abstract: The axion is a promising dark matter candidate, which was originally proposed to solve the strong-CP problem in particle physics. To date, the available parameter space for axion and axion-like particle dark matter is relatively unexplored, particularly at masses $m_a\lesssim1\,\mu$eV. ABRACADABRA is a new experimental program to search for axion dark matter over a broad range of masses, $10^{-12}\lesssim m_a\lesssim10^{-6}$ eV. ABRACADABRA-10 cm is a small-scale prototype for a future detector that could be sensitive to the QCD axion. In this Letter, we present the first results from a 1 month search for axions with ABRACADABRA-10 cm. We find no evidence for axion-like cosmic dark matter and set 95% C.L. upper limits on the axion-photon coupling between $g_{a\gamma\gamma}<1.4\times10^{-10}$ GeV$^{-1}$ and $g_{a\gamma\gamma}<3.3\times10^{-9}$ GeV$^{-1}$ over the mass range $3.1\times10^{-10}$ eV - $8.3\times10^{-9}$ eV. These results are competitive with the most stringent astrophysical constraints in this mass range.


title: Reflexive polytopes arising from bipartite graphs with $\gamma$-positivity associated to interior polynomials

abstract: In this paper, we introduce polytopes ${\mathcal B}_G$ arising from root systems $B_n$ and finite graphs $G$, and study their combinatorial and algebraic properties. In particular, it is shown that ${\mathcal B}_G$ is reflexive if and only if $G$ is bipartite. Moreover, in the case, ${\mathcal B}_G$ has a regular unimodular triangulation. This implies that the $h^*$-polynomial of ${\mathcal B}_G$ is palindromic and unimodal when $G$ is bipartite. Furthermore, we discuss stronger properties, namely the $\gamma$-positivity and the real-rootedness of the $h^*$-polynomials. In fact, if $G$ is bipartite, then the $h^*$-polynomial of ${\mathcal B}_G$ is $\gamma$-positive and its $\gamma$-polynomial is given by an interior polynomial (a version of the Tutte polynomial for a hypergraph). The $h^*$-polynomial is real-rooted if and only if the corresponding interior polynomial is real-rooted. From a counterexample to Neggers--Stanley conjecture, we construct a bipartite graph $G$ whose $h^*$-polynomial is not real-rooted but $\gamma$-positive, and coincides with the $h$-polynomial of a flag triangulation of a sphere.


title: H-He Shell Interactions and Nucleosynthesis in Massive Population III Stars

abstract: We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code.


title: Wireless Terahertz System Architectures for Networks Beyond 5G

abstract: The present white paper focuses on the system requirements of TERRANOVA. Initially details the key use cases for the TERRANOVA technology and presents the description of the network architecture. In more detail, the use cases are classified into two categories, namely backhaul & fronthaul and access and small cell backhaul. The first category refers to fibre extender, point-to-point and redundancy applications, whereas the latter is designed to support backup connection for small and medium-sized enterprises (SMEs), internet of things (IoT) dense environments, data centres, indoor wireless access, ad hoc networks, and last mile access. Then, it provides the networks architecture for the TERRANOVA system as well as the network elements that need to be deployed. The use cases are matched to specific technical scenarios, namely outdoor fixed point-to-point (P2P), outdoor/indoor individual point-to-multipoint (P2MP), and outdoor/indoor "quasi"-omnidirection, and the key performance requirements of each scenario are identified. Likewise, we present the breakthrough novel technology concepts, including the joint design of baseband signal processing for the complete optical and wireless link, the development of broadband and spectrally efficient RF-frontends for frequencies >275 GHz, as well as channel modelling, waveforms, antenna array and multiple-access schemes design, which we are going to use in order to satisfy the presented requirements. Next, an overview of the required new functionalities in both physical (PHY) layer and medium access control (MAC) layers in the TERRANOVA system architecture will be given. Finally, the individual enablers of the TERRANOVA system are combined to develop particular candidate architectures for each of the three technical scenarios.


title: Automated Solution of First Order Factorizable Systems of Differential Equations in One Variable

abstract: We present an algorithm which allows to solve analytically linear systems of differential equations which factorize to first order. The solution is given in terms of iterated integrals over an alphabet where its structure is implied by the coefficient matrix of the differential equations. These systems appear in a large variety of higher order calculations in perturbative Quantum Field Theories. We apply this method to calculate the master integrals of the three--loop massive form factors for different currents, as an illustration, and present the results for the vector form factors in detail. Here the solution space emerging is given by the cyclotomic harmonic polylogarithms and their associated special constants. No special basis representation of the master integrals is needed. The algorithm can be applied as well to more general cases factorizing at first order, which are based on more general alphabets, iterated integrals and associated constants.


title: High Open Circuit Voltages in pin-Type Perovskite Solar Cells through Strontium Addition

abstract: The incorporation of even small amounts of strontium (Sr) into lead-based quadruple cation hybrid perovskite solar cells results in a systematic increase of the open circuit voltage (Voc) in pin-type perovskite solar cells. We demonstrate via transient and absolute photoluminescence (PL) experiments how the incorporation of Sr significantly reduces the non-radiative recombination losses in the neat perovskite layer and specifically at the perovskite/C60 interface. We show that Sr segregates at the perovskite surface, where it induces important changes of morphology and energetics. Notably, the Sr-enriched surface exhibits a wider band gap and a more n-type character, accompanied with significantly stronger surface band bending. As a result, we observe a significant increase of the quasi-Fermi level splitting in the neat perovskite by reduced surface recombination and more importantly, a strong reduction of losses attributed to non-radiative recombination at the interface to the C60 electron-transporting layer. The resulting solar cells exhibited a Voc of 1.18 V, which could be further improved to nearly 1.23 V through addition of a thin polymer interlayer, bringing the non-radiative voltage loss to only 110 meV. Our work shows that simply adding a small amount of Sr to the precursor solutions induces a beneficial surface modification in the perovskite, without requiring any post treatment, resulting in high efficiency solar cells with power conversion efficiency (PCE) up to 20.3%. Our results demonstrate very high Voc values and efficiencies in Sr-containing quadruple cation perovskite pin solar cells and highlight the imperative importance of addressing and minimizing the recombination losses at the interface between perovskite and charge transporting layer.


title: Learning Gaussian Processes by Minimizing PAC-Bayesian Generalization Bounds

abstract: Gaussian Processes (GPs) are a generic modelling tool for supervised learning. While they have been successfully applied on large datasets, their use in safety-critical applications is hindered by the lack of good performance guarantees. To this end, we propose a method to learn GPs and their sparse approximations by directly optimizing a PAC-Bayesian bound on their generalization performance, instead of maximizing the marginal likelihood. Besides its theoretical appeal, we find in our evaluation that our learning method is robust and yields significantly better generalization guarantees than other common GP approaches on several regression benchmark datasets.


title: Learning Comment Generation by Leveraging User-Generated Data

abstract: Existing models on open-domain comment generation are difficult to train, and they produce repetitive and uninteresting responses. The problem is due to multiple and contradictory responses from a single article, and by the rigidity of retrieval methods. To solve this problem, we propose a combined approach to retrieval and generation methods. We propose an attentive scorer to retrieve informative and relevant comments by leveraging user-generated data. Then, we use such comments, together with the article, as input for a sequence-to-sequence model with copy mechanism. We show the robustness of our model and how it can alleviate the aforementioned issue by using a large scale comment generation dataset. The result shows that the proposed generative model significantly outperforms strong baseline such as Seq2Seq with attention and Information Retrieval models by around 27 and 30 BLEU-1 points respectively.


title: Concepts of Architecture, Structure and System

abstract: The current ISO standards pertaining to the Concepts of System and Architecture express succinct definitions of these two key terms that lend themselves to practical application and can be understood through elementary mathematical foundations. The current work of the ISO/IEC Working Group 42 is seeking to refine and elaborate the existing standards. This position paper revisits the fundamental concepts underlying both of these key terms and offers an approach to: (i) refine and exemplify the term 'fundamental concepts' in the current ISO definition of Architecture, (ii) exploit existing standards for the term 'concept', and (iii) introduce a new concept, Architectural Structure, that can serve to unify the current terminology at a fundamental level. Precise elementary examples are used in to conceptualise the approach offered.


title: The problem with probabilistic DAG automata for semantic graphs

abstract: Semantic representations in the form of directed acyclic graphs (DAGs) have been introduced in recent years, and to model them, we need probabilistic models of DAGs. One model that has attracted some attention is the DAG automaton, but it has not been studied as a probabilistic model. We show that some DAG automata cannot be made into useful probabilistic models by the nearly universal strategy of assigning weights to transitions. The problem affects single-rooted, multi-rooted, and unbounded-degree variants of DAG automata, and appears to be pervasive. It does not affect planar variants, but these are problematic for other reasons.


title: A catalog of 29 open clusters and associations observed by the Kepler and K2 Missions

abstract: Over the past nine years, the Kepler and K2 Missions have carried out high precision photometric monitoring of more than half a million stars. Among these targets are 29 clusters and associations, with ages from 1 Myr to over 11 Gyr. We have generated a catalog of Kepler/K2 clusters, including basic information about the observations as well as the degree of scientific attention paid to them thus far. We present the catalog and a discussion of current and future cluster science endeavors as a way to motivate the astronomical community to continue mining these exquisite datasets.


title: Superconductivity in the Hubbard model: a hidden-order diagnostics from the Luther-Emery phase on ladders

abstract: Short-range antiferromagnetic correlations are known to open a spin gap in the repulsive Hubbard model on ladders with $M$ legs, when $M$ is even. We show that the spin gap originates from the formation of correlated pairs of electrons with opposite spin, captured by the hidden ordering of a spin-parity operator. Since both spin gap and parity vanish in the two-dimensional limit, we introduce the fractional generalization of spin parity and prove that it remains finite in the thermodynamic limit. Our results are based upon variational wave functions and Monte Carlo calculations: performing a finite size-scaling analysis with growing $M$, we show that the doping region where the parity is finite coincides with the range in which superconductivity is observed in two spatial dimensions. Our observations support the idea that superconductivity emerges out of spin gapped phases on ladders, driven by a spin-pairing mechanism, in which the ordering is conveniently captured by the finiteness of the fractional spin-parity operator.


title: RG transport theory for open quantum systems: Charge fluctuations in multilevel quantum dots in and out of equilibrium

abstract: We present the real-time renormalization group (RTRG) method as a method to describe the stationary state current through generic multi-level quantum dots with a complex setup in nonequilibrium. The employed approach consists of a very rudiment approximation for the RG equations which neglects all vertex corrections while it provides a means to compute the effective dot Liouvillian self-consistently. Being based on a weak-coupling expansion in the tunneling between dot and reservoirs, the RTRG approach turns out to reliably describe charge fluctuations in and out of equilibrium for arbitrary coupling strength, even at zero temperature. We confirm this in the linear response regime with a benchmark against highly-accurate numerically renormalization group data in the exemplary case of three-level quantum dots. For small to intermediate bias voltages and weak Coulomb interactions, we find an excellent agreement between RTRG and functional renormalization group data, which can be expected to be accurate in this regime. As a consequence, we advertise the presented RTRG approach as an efficient and versatile tool to describe charge fluctuations theoretically in quantum dot systems.


title: Computing the endomorphism ring of an ordinary abelian surface over a finite field

abstract: We present a new algorithm for computing the endomorphism ring of an ordinary abelian surface over a finite field which is subexponential and generalizes an algorithm of Bisson and Sutherland for elliptic curves. The correctness of this algorithm only requires the heuristic assumptions required by the algorithm of Biasse and Fieker which computes the class group of an order in a number field in subexponential time. Thus we avoid the multiple heuristic assumptions on isogeny graphs and polarized class groups which were previously required. The output of the algorithm is an ideal in the maximal totally real subfield of the endomorphism algebra, generalizing the elliptic curve case.


title: Toward Creating Subsurface Camera

abstract: In this article, the framework and architecture of Subsurface Camera (SAMERA) is envisioned and described for the first time. A SAMERA is a geophysical sensor network that senses and processes geophysical sensor signals, and computes a 3D subsurface image in-situ in real-time. The basic mechanism is: geophysical waves propagating/reflected/refracted through subsurface enter a network of geophysical sensors, where a 2D or 3D image is computed and recorded; a control software may be connected to this network to allow view of the 2D/3D image and adjustment of settings such as resolution, filter, regularization and other algorithm parameters. System prototypes based on seismic imaging have been designed. SAMERA technology is envisioned as a game changer to transform many subsurface survey and monitoring applications, including oil/gas exploration and production, subsurface infrastructures and homeland security, wastewater and CO2 sequestration, earthquake and volcano hazard monitoring. The system prototypes for seismic imaging have been built. Creating SAMERA requires an interdisciplinary collaboration and transformation of sensor networks, signal processing, distributed computing, and geophysical imaging.


title: Adversarial Risk and Robustness: General Definitions and Implications for the Uniform Distribution

abstract: We study adversarial perturbations when the instances are uniformly distributed over $\{0,1\}^n$. We study both "inherent" bounds that apply to any problem and any classifier for such a problem as well as bounds that apply to specific problems and specific hypothesis classes. As the current literature contains multiple definitions of adversarial risk and robustness, we start by giving a taxonomy for these definitions based on their goals, we identify one of them as the one guaranteeing misclassification by pushing the instances to the error region. We then study some classic algorithms for learning monotone conjunctions and compare their adversarial risk and robustness under different definitions by attacking the hypotheses using instances drawn from the uniform distribution. We observe that sometimes these definitions lead to significantly different bounds. Thus, this study advocates for the use of the error-region definition, even though other definitions, in other contexts, may coincide with the error-region definition. Using the error-region definition of adversarial perturbations, we then study inherent bounds on risk and robustness of any classifier for any classification problem whose instances are uniformly distributed over $\{0,1\}^n$. Using the isoperimetric inequality for the Boolean hypercube, we show that for initial error $0.01$, there always exists an adversarial perturbation that changes $O(\sqrt{n})$ bits of the instances to increase the risk to $0.5$, making classifier's decisions meaningless. Furthermore, by also using the central limit theorem we show that when $n\to \infty$, at most $c \cdot \sqrt{n}$ bits of perturbations, for a universal constant $c< 1.17$, suffice for increasing the risk to $0.5$, and the same $c \cdot \sqrt{n} $ bits of perturbations on average suffice to increase the risk to $1$, hence bounding the robustness by $c \cdot \sqrt{n}$.


title: Kalman Gradient Descent: Adaptive Variance Reduction in Stochastic Optimization

abstract: We introduce Kalman Gradient Descent, a stochastic optimization algorithm that uses Kalman filtering to adaptively reduce gradient variance in stochastic gradient descent by filtering the gradient estimates. We present both a theoretical analysis of convergence in a non-convex setting and experimental results which demonstrate improved performance on a variety of machine learning areas including neural networks and black box variational inference. We also present a distributed version of our algorithm that enables large-dimensional optimization, and we extend our algorithm to SGD with momentum and RMSProp.


title: Phase field modelling of surfactants in multi-phase flow

abstract: A diffuse interface model for surfactants in multi-phase flow with three or more fluids is derived. A system of Cahn-Hilliard equations is coupled with a Navier-Stokes system and an advection-diffusion equation for the surfactant ensuring thermodynamic consistency. By an asymptotic analysis the model can be related to a moving boundary problem in the sharp interface limit, which is derived from first principles. Results from numerical simulations support the theoretical findings. The main novelties are centred around the conditions in the triple junctions where three fluids meet. Specifically the case of local chemical equilibrium with respect to the surfactant is considered, which allows for interfacial surfactant flow through the triple junctions.


title: Abelian Anti-Powers in Infinite Words

abstract: An abelian anti-power of order $k$ (or simply an abelian $k$-anti-power) is a concatenation of $k$ consecutive words of the same length having pairwise distinct Parikh vectors. This definition generalizes to the abelian setting the notion of a $k$-anti-power, as introduced in [G. Fici et al., Anti-powers in infinite words, J. Comb. Theory, Ser. A, 2018], that is a concatenation of $k$ pairwise distinct words of the same length. We aim to study whether a word contains abelian $k$-anti-powers for arbitrarily large $k$. S. Holub proved that all paperfolding words contain abelian powers of every order [Abelian powers in paper-folding words. J. Comb. Theory, Ser. A, 2013]. We show that they also contain abelian anti-powers of every order.


title: Scale transformations in metric-affine geometry

abstract: This article presents an exhaustive classification of metric-affine theories according to their scale symmetries. First it is clarified that there are three relevant definitions of a scale transformation. These correspond to a projective transformation of the connection, a rescaling of the orthonormal frame, and a combination of the two. The most general second order quadratic metric-affine action, including the parity-violating terms, is constructed in each of the three cases. The results can be straightforwardly generalised by including higher derivatives, and implemented in the general metric-affine, teleparallel, and symmetric teleparallel geometries.


title: A Maximum Linear Arrangement Problem on Directed Graphs

abstract: We propose a new arrangement problem on directed graphs, Maximum Directed Linear Arrangement (MaxDLA). This is a directed variant of a similar problem for undirected graphs, in which however one seeks maximum and not minimum; this problem known as the Minimum Linear Arrangement Problem (MinLA) has been much studied in the literature. We establish a number of theorems illustrating the behavior and complexity of MaxDLA. First, we relate MaxDLA to Maximum Directed Cut (MaxDiCut) by proving that every simple digraph $D$ on $n$ vertices satisfies $\frac{n}{2}$$maxDiCut(D) \leq MaxDLA(D) \leq (n-1)MaxDiCut(D)$. Next, we prove that MaxDiCut is NP-Hard for planar digraphs (even with the added restriction of maximum degree 15); it follows from the above bounds that MaxDLA is also NP-Hard for planar digraphs. In contrast, Hadlock (1975) and Dorfman and Orlova (1972) showed that the undirected Maximum Cut problem is solvable in polynomial time on planar graphs. On the positive side, we present a polynomial-time algorithm for solving MaxDLA on orientations of trees with degree bounded by a constant, which translates to a polynomial-time algorithm for solving MinLA on the complements of those trees. This pairs with results by Goldberg and Klipker (1976), Shiloach (1979) and Chung (1984) solving MinLA in polynomial time on trees. Finally, analogues of Harper's famous isoperimetric inequality for the hypercube, in the setting of MaxDLA, are shown for tournaments, orientations of graphs with degree at most two, and transitive acyclic digraphs.


title: Towards Principled Uncertainty Estimation for Deep Neural Networks

abstract: When the cost of misclassifying a sample is high, it is useful to have an accurate estimate of uncertainty in the prediction for that sample. There are also multiple types of uncertainty which are best estimated in different ways, for example, uncertainty that is intrinsic to the training set may be well-handled by a Bayesian approach, while uncertainty introduced by shifts between training and query distributions may be better-addressed by density/support estimation. In this paper, we examine three types of uncertainty: model capacity uncertainty, intrinsic data uncertainty, and open set uncertainty, and review techniques that have been derived to address each one. We then introduce a unified hierarchical model, which combines methods from Bayesian inference, invertible latent density inference, and discriminative classification in a single end-to-end deep neural network topology to yield efficient per-sample uncertainty estimation in a detection context. This approach addresses all three uncertainty types and can readily accommodate prior/base rates for binary detection. We then discuss how to extend this model to a more generic multiclass recognition context.


title: On a cohomological generalization of the Shafarevich conjecture for K3 surfaces

abstract: The Shafarevich conjecture for K3 surfaces asserts the finiteness of isomorphism classes of K3 surfaces over a fixed number field admitting good reduction away from a fixed finite set of finite places. Andr\'{e} proved this conjecture for polarized K3 surfaces of fixed degree, and recently She proved it for polarized K3 surfaces of unspecified degree. In this paper, we prove a certain generalization of their results, which is stated by the unramifiedness of l-adic \'{e}tale cohomology groups for K3 surfaces over finitely generated fields of characteristic 0. As a corollary, we get the original Shafarevich conjecture for K3 surfaces without assuming the extendability of polarization, which is stronger than the results of Andr\'{e} and She. Moreover, as an application, we get the finiteness of twists of K3 surfaces via a finite extension of characteristic 0 fields.


title: Quantum theory of superfluorescence based on two-point correlation functions

abstract: Irradiation of a medium by short intense pulses from x-ray / XUV free electron lasers can result in saturated photoionization of inner electronic shells. As a result an inversion of populations between core levels appears. The resulting fluorescent radiation can be amplified during its propagation through the inverted medium and results in intense, quasi transform-limited radiation bursts. While the optical counterpart of this phenomena, known as superfluorescence, was intensively investigated, a generalized treatment is needed in the x-ray / XUV domain, where the dynamics of pumping and evolution due to fast decay processes play a crucial role. To provide a general theoretical approach, we start from the fundamental, quantized minimal coupling Hamiltonian of light-matter interaction and after a series of approximations arrive at a closed system of equations for the two-point correlation function of atomic coherences and the two-time correlation function of the emitted field. The obtained formalism enables us to investigate collective spontaneous emission in various regimes. It is extended consistently to include incoherent processes that are relevant in the x-ray / XUV domain. These processes are introduced into the formalism by corresponding Lindblad superoperators. The connection to other approaches is discussed and numerical examples related to recent experiments are presented.


title: Three Mechanisms of Weight Decay Regularization

abstract: Weight decay is one of the standard tricks in the neural network toolbox, but the reasons for its regularization effect are poorly understood, and recent results have cast doubt on the traditional interpretation in terms of $L_2$ regularization. Literal weight decay has been shown to outperform $L_2$ regularization for optimizers for which they differ. We empirically investigate weight decay for three optimization algorithms (SGD, Adam, and K-FAC) and a variety of network architectures. We identify three distinct mechanisms by which weight decay exerts a regularization effect, depending on the particular optimization algorithm and architecture: (1) increasing the effective learning rate, (2) approximately regularizing the input-output Jacobian norm, and (3) reducing the effective damping coefficient for second-order optimization. Our results provide insight into how to improve the regularization of neural networks.


title: Assessing Generalization in Deep Reinforcement Learning

abstract: Deep reinforcement learning (RL) has achieved breakthrough results on many tasks, but agents often fail to generalize beyond the environment they were trained in. As a result, deep RL algorithms that promote generalization are receiving increasing attention. However, works in this area use a wide variety of tasks and experimental setups for evaluation. The literature lacks a controlled assessment of the merits of different generalization schemes. Our aim is to catalyze community-wide progress on generalization in deep RL. To this end, we present a benchmark and experimental protocol, and conduct a systematic empirical study. Our framework contains a diverse set of environments, our methodology covers both in-distribution and out-of-distribution generalization, and our evaluation includes deep RL algorithms that specifically tackle generalization. Our key finding is that `vanilla' deep RL algorithms generalize better than specialized schemes that were proposed specifically to tackle generalization.


title: Scaling Gaussian Process Regression with Derivatives

abstract: Gaussian processes (GPs) with derivatives are useful in many applications, including Bayesian optimization, implicit surface reconstruction, and terrain reconstruction. Fitting a GP to function values and derivatives at $n$ points in $d$ dimensions requires linear solves and log determinants with an ${n(d+1) \times n(d+1)}$ positive definite matrix -- leading to prohibitive $\mathcal{O}(n^3d^3)$ computations for standard direct methods. We propose iterative solvers using fast $\mathcal{O}(nd)$ matrix-vector multiplications (MVMs), together with pivoted Cholesky preconditioning that cuts the iterations to convergence by several orders of magnitude, allowing for fast kernel learning and prediction. Our approaches, together with dimensionality reduction, enables Bayesian optimization with derivatives to scale to high-dimensional problems and large evaluation budgets.


title: Anisotropic Correlated Electronic Structure of Colossal Thermopower Marcasite FeSb$_2$

abstract: Iron antimonide (FeSb$_2$) is a mysterious material with peculiar colossal thermopower of about $-45$ mV/K at 10 K. However, a unified microscopic description of this phenomenon is far from being achieved. The understanding of the electronic structure in details is crucial in identifying the microscopic mechanism of FeSb$_2$ thermopower. Combining angle-resolved photoemission spectroscopy (ARPES) and first-principles calculations we find that the spectrum of FeSb$_2$ consists of two bands near the Fermi energy: the nondispersive strongly renormalized $\alpha$-band, and the hole-like $\beta$-band that intersects the first one at $\Gamma$ and Y points of the Brillouin zone. Our study reveals the presence of sizable correlations, predominantly among electrons derived from Fe-3d states, and considerable anisotropy in the electronic structure of FeSb$_2$. These key ingredients are of fundamental importance in the description of colossal thermopower in FeSb$_2$.


title: Study of doubly heavy tetraquarks in Lattice QCD

abstract: We present results of a lattice calculation of tetraquark states with quark contents $q_1q_2\bar{Q}\bar{Q}, \, q_1,q_2 \subset u,d,s,c$ and $Q \equiv b,c$ in both spin zero ($J=0$) and spin one ($J=1$) sectors. These calculations are performed on three dynamical $N_f = 2 + 1 + 1$ highly improved staggered quark ensembles at lattice spacings of about 0.12, 0.09 and 0.06 fm. We use the overlap action for light to charm quarks while a non-relativistic action with non-perturbatively improved coefficients with terms up to $\mathcal{O}(\alpha_s v^4)$ is employed for the bottom quark. While considering two heavy quarks as charm or bottom, we calculate the energy levels of various four-quark configurations with light quark masses ranging from the physical strange quark mass to that of the corresponding physical pion mass. Results for the spin one states show the presence of ground state energy levels which are below their respective thresholds for all the light flavor combinations with both doubly heavy quarks and particularly for the bottom quarks. Further, we identify a trend that the energy splittings, defined as the energy difference between the ground state energy levels and their respective thresholds, increase with decreasing the light quark masses and are maximum at the physical point for all the spin one states. The rate of increase is however dependent on the light quark configuration of the particular spin one state. We also present a study of hadron mass relations involving tetraquarks, baryons and mesons arising in the limit of infinitely heavy quark and find that these relations are more compatible with the heavy quark limit in the bottom sector but deviate substantially in the charm sector. The ground state spectra of the spin zero tetraquark states with various flavor combinations are seen to lie above their respective thresholds.


title: Compressive Sampling Approach for Image Acquisition with Lensless Endoscope

abstract: The lensless endoscope is a promising device designed to image tissues in vivo at the cellular scale. The traditional acquisition setup consists in raster scanning during which the focused light beam from the optical fiber illuminates sequentially each pixel of the field of view (FOV). The calibration step to focus the beam and the sampling scheme both take time. In this preliminary work, we propose a scanning method based on compressive sampling theory. The method does not rely on a focused beam but rather on the random illumination patterns generated by the single-mode fibers. Experiments are performed on synthetic data for different compression rates (from 10 to 100% of the FOV).