title: Generating Master Faces for Use in Performing Wolf Attacks on Face Recognition Systems

abstract: Due to its convenience, biometric authentication, especial face authentication, has become increasingly mainstream and thus is now a prime target for attackers. Presentation attacks and face morphing are typical types of attack. Previous research has shown that finger-vein- and fingerprint-based authentication methods are susceptible to wolf attacks, in which a wolf sample matches many enrolled user templates. In this work, we demonstrated that wolf (generic) faces, which we call "master faces," can also compromise face recognition systems and that the master face concept can be generalized in some cases. Motivated by recent similar work in the fingerprint domain, we generated high-quality master faces by using the state-of-the-art face generator StyleGAN in a process called latent variable evolution. Experiments demonstrated that even attackers with limited resources using only pre-trained models available on the Internet can initiate master face attacks. The results, in addition to demonstrating performance from the attacker's point of view, can also be used to clarify and improve the performance of face recognition systems and harden face authentication systems.


title: Continuity equations for entanglement

abstract: We introduce a complex purity density and its associated current for pure states of continuous variable systems. The scheme is constructed by analogy with the notions of probability density and probability current. Taking advantage of the formal continuity equations obtained this way we can introduce an entanglement subdynamics. We suggest the use of the dimensionality of this subdynamics as a potential measure of the complexity of entanglement. The scheme also provides insights into the relation between the global and local aspects of quantum correlations.


title: Optical computation of a spin glass dynamics with tunable complexity

abstract: Spin Glasses (SG) are paradigmatic models for physical, computer science, biological and social systems. The problem of studying the dynamics for SG models is NP hard, i.e., no algorithm solves it in polynomial time. Here we implement the optical simulation of a SG, exploiting the N segments of a wavefront shaping device to play the role of the spin variables, combining the interference at downstream of a scattering material to implement the random couplings between the spins (the J ij matrix) and measuring the light intensity on a number P of targets to retrieve the energy of the system. By implementing a plain Metropolis algorithm, we are able to simulate the spin model dynamics, while the degree of complexity of the potential energy landscape and the region of phase diagram explored is user-defined acting on the ratio the P/N = \alpha. We study experimentally, numerically and analytically this peculiar system displaying a paramagnetic, a ferromagnetic and a SG phase, and we demonstrate that the transition temperature T g to the glassy phase from the paramagnetic phase grows with \alpha. With respect to standard in silico approach, in the optical SG interaction terms are realized simultaneously when the independent light rays interferes at the target screen, enabling inherently parallel measurements of the energy, rather than computations scaling with N as in purely in silico simulations.


title: Elliptic problems with unknowns on the boundary and irregular boundary data

abstract: We consider an elliptic problem with unknowns on the boundary of the domain of the elliptic equation and suppose that the right-hand side of this equation is square integrable and that the boundary data are arbitrary (specifically, irregular) distributions. We investigate local (up to the boundary) properties of generalized solutions to the problem in Hilbert distribution spaces that belong to the refined Sobolev scale. These spaces are parametrized with a real number and a function that varies slowly at infinity. The function parameter refines the number order of the space. We prove theorems on local regularity and a local a priori estimate of generalized solutions to the problem under investigation. These theorems are new for Sobolev spaces as well.


title: Causal Inference with Deep Causal Graphs

abstract: Parametric causal modelling techniques rarely provide functionality for counterfactual estimation, often at the expense of modelling complexity. Since causal estimations depend on the family of functions used to model the data, simplistic models could entail imprecise characterizations of the generative mechanism, and, consequently, unreliable results. This limits their applicability to real-life datasets, with non-linear relationships and high interaction between variables. We propose Deep Causal Graphs, an abstract specification of the required functionality for a neural network to model causal distributions, and provide a model that satisfies this contract: Normalizing Causal Flows. We demonstrate its expressive power in modelling complex interactions and showcase applications of the method to machine learning explainability and fairness, using true causal counterfactuals.


title: DreamCoder: Growing generalizable, interpretable knowledge with wake-sleep Bayesian program learning

abstract: Expert problem-solving is driven by powerful languages for thinking about problems and their solutions. Acquiring expertise means learning these languages -- systems of concepts, alongside the skills to use them. We present DreamCoder, a system that learns to solve problems by writing programs. It builds expertise by creating programming languages for expressing domain concepts, together with neural networks to guide the search for programs within these languages. A ``wake-sleep'' learning algorithm alternately extends the language with new symbolic abstractions and trains the neural network on imagined and replayed problems. DreamCoder solves both classic inductive programming tasks and creative tasks such as drawing pictures and building scenes. It rediscovers the basics of modern functional programming, vector algebra and classical physics, including Newton's and Coulomb's laws. Concepts are built compositionally from those learned earlier, yielding multi-layered symbolic representations that are interpretable and transferrable to new tasks, while still growing scalably and flexibly with experience.


title: Asymptotic regularity and attractors for slightly compressible Brinkman-Forcheimer equations

abstract: Slightly compressible Brinkman-Forchheimer equations in a bounded 3D domain with Dirichlet boundary conditions are considered. These equations model fluids motion in porous media. The dissipativity of these equations in higher order energy spaces is obtained and regularity and smoothing properties of the solutions are studied. In addition, the existence of a global and an exponential attractors for these equations in a natural phase space is verified.


title: Pixel Invisibility: Detecting Objects Invisible in Color Images

abstract: Despite recent success of object detectors using deep neural networks, their deployment on safety-critical applications such as self-driving cars remains questionable. This is partly due to the absence of reliable estimation for detectors' failure under operational conditions such as night, fog, dusk, dawn and glare. Such unquantifiable failures could lead to safety violations. In order to solve this problem, we created an algorithm that predicts a pixel-level invisibility map for color images that does not require manual labeling - that computes the probability that a pixel/region contains objects that are invisible in color domain, during various lighting conditions such as day, night and fog. We propose a novel use of cross modal knowledge distillation from color to infra-red domain using weakly-aligned image pairs from the day and construct indicators for the pixel-level invisibility based on the distances of their intermediate-level features. Quantitative experiments show the great performance of our pixel-level invisibility mask and also the effectiveness of distilled mid-level features on object detection in infra-red imagery.


title: Equivalence of weak and viscosity solutions in fractional non-homogeneous problems

abstract: We establish the equivalence between the notions of weak and viscosity solutions for non-homogeneous equations whose main operator is the fractional p-Laplacian and the lower order term depends on $x$, $u$ and $D_s^p u$, being the last one a type of fractional derivative.


title: Plasmon-Driven Hot Electron Transfer at Atomically Sharp Metal-Semiconductor Nanojunctions

abstract: Recent advances in guiding and localizing light at the nanoscale exposed the enormous potential of ultra-scaled plasmonic devices. In this context, the decay of surface plasmons to hot carriers triggers a variety of applications in boosting the efficiency of energy-harvesting, photo-catalysis and photo-detection. However, a detailed understanding of plasmonic hot carrier generation and particularly the transfer at metal-semiconductor interfaces is still elusive. In this paper, we introduce a monolithic metal-semiconductor (Al-Ge) heterostructure device, providing a platform to examine surface plasmon decay and hot electron transfer at an atomically sharp Schottky nanojunction. The gated metal-semiconductor heterojunction device features electrostatic control of the Schottky barrier height at the Al-Ge interface, enabling hot electron filtering. The ability of momentum matching and to control the energy distribution of plasmon-driven hot electron injection is demonstrated by controlling the interband electron transfer in Ge leading to negative differential resistance.


title: COALA: Co-Aligned Autoencoders for Learning Semantically Enriched Audio Representations

abstract: Audio representation learning based on deep neural networks (DNNs) emerged as an alternative approach to hand-crafted features. For achieving high performance, DNNs often need a large amount of annotated data which can be difficult and costly to obtain. In this paper, we propose a method for learning audio representations, aligning the learned latent representations of audio and associated tags. Aligning is done by maximizing the agreement of the latent representations of audio and tags, using a contrastive loss. The result is an audio embedding model which reflects acoustic and semantic characteristics of sounds. We evaluate the quality of our embedding model, measuring its performance as a feature extractor on three different tasks (namely, sound event recognition, and music genre and musical instrument classification), and investigate what type of characteristics the model captures. Our results are promising, sometimes in par with the state-of-the-art in the considered tasks and the embeddings produced with our method are well correlated with some acoustic descriptors.


title: Catplayinginthesnow: Impact of Prior Segmentation on a Model of Visually Grounded Speech

abstract: The language acquisition literature shows that children do not build their lexicon by segmenting the spoken input into phonemes and then building up words from them, but rather adopt a top-down approach and start by segmenting word-like units and then break them down into smaller units. This suggests that the ideal way of learning a language is by starting from full semantic units. In this paper, we investigate if this is also the case for a neural model of Visually Grounded Speech trained on a speech-image retrieval task. We evaluated how well such a network is able to learn a reliable speech-to-image mapping when provided with phone, syllable, or word boundary information. We present a simple way to introduce such information into an RNN-based model and investigate which type of boundary is the most efficient. We also explore at which level of the network's architecture such information should be introduced so as to maximise its performances. Finally, we show that using multiple boundary types at once in a hierarchical structure, by which low-level segments are used to recompose high-level segments, is beneficial and yields better results than using low-level or high-level segments in isolation.


title: Low-energy collisions between carbon atoms and oxygen molecules in a magnetic trap

abstract: Trapping of atoms and molecules in electrostatic, magnetic and optical traps has enabled studying atomic and molecular interactions on a timescale of many seconds, allowing observations of ultra-cold collisions and reactions. Here we report the first magnetic deceleration and trapping of neutral carbon atoms in a static magnetic trap. When co-trapping the carbon atoms with oxygen molecules in a superconducting trap, the carbon signal decays in a non-exponential manner, consistent with losses resulting from atom-molecule collisions. Our findings pave the way to studying both elastic and inelastic collisions of species that cannot be laser cooled, and specifically may facilitate the observation of reactions at low temperatures, such as C + O2 --> CO + O, which is important in interstellar chemistry.


title: Quantitative differentiation of protein aggregates from other subvisible particles in viscous mixtures through holographic characterization

abstract: We demonstrate the use of holographic video microscopy to detect individual subvisible particles dispersed in biopharmaceutical formulations and to differentiate them based on material characteristics measured from their holograms. The result of holographic analysis is a precise and accurate measurement of the concentrations and size distributions of multiple classes of subvisible contaminants dispersed in the same product simultaneously. We demonstrate this analytical technique through measurements on model systems consisting of human IgG aggregates in the presence of common contaminants such as silicone oil emulsion droplets and fatty acids. Holographic video microscopy also clearly identifies metal particles and air bubbles. Being able to differentiate and characterize the individual components of such heterogeneous dispersions provides a basis for tracking other factors that influence the stability of protein formulations including handling and degradation of surfactant and other excipients.


title: Double-valued strong-coupling corrections to Bardeen-Cooper-Schrieffer ratios

abstract: Experimental discovery of near-room-temperature (NRT) superconductivity in highly-compressed H$_3$S, LaH$_1$$_0$ and YH$_6$ restores fundamental interest to electron-phonon pairing mechanism in superconductors. One of prerequisites of phonon-mediated NRT superconductivity in highly-compressed hydrides is strong electron-phonon interaction, which can be quantified by dimensionless ratios of Bardeen-Cooper-Schrieffer (BCS) theory vs $k{_B}T{_c}$/(${\hslash}$${\omega}{_l}{_n}$), where Tc is the critical temperature and ${\omega}{_l}{_n}$ is the logarithmic phonon frequency (Mitrovic et al. 1984 Phys. Rev. B 29 184). However, all known strong-coupling correction functions for BCS ratios are applicable for $k{_B}T{_c}$/(${\hslash}$${\omega}{_l}{_n}$)<0.20, which is not high enough $k{_B}T{_c}$/(${\hslash}$${\omega}{_l}{_n}$) range for NRT superconductors, because the latter exhibit variable values of 0.13 < $k{_B}T{_c}$/(${\hslash}$${\omega}{_l}{_n}$) < 0.32. In this paper, we reanalyze full experimental dataset (including data for highly-compressed H3S) and find that strong-coupling correction functions for the gap-to-critical-temperature ratio and for the specific-heat-jump ratio are double-valued nearly-linear functions of $k{_B}T{_c}$/(${\hslash}$${\omega}{_l}{_n}$).


title: On Lipschitz Regularization of Convolutional Layers using Toeplitz Matrix Theory

abstract: This paper tackles the problem of Lipschitz regularization of Convolutional Neural Networks. Lipschitz regularity is now established as a key property of modern deep learning with implications in training stability, generalization, robustness against adversarial examples, etc. However, computing the exact value of the Lipschitz constant of a neural network is known to be NP-hard. Recent attempts from the literature introduce upper bounds to approximate this constant that are either efficient but loose or accurate but computationally expensive. In this work, by leveraging the theory of Toeplitz matrices, we introduce a new upper bound for convolutional layers that is both tight and easy to compute. Based on this result we devise an algorithm to train Lipschitz regularized Convolutional Neural Networks.


title: Exciton-exciton interaction in transition metal dichalcogenide monolayers and van der Waals heterostructures

abstract: Due to a strong Coulomb interaction, excitons dominate the excitation kinetics in 2D materials. While Coulomb-scattering between electrons has been well studied, the interaction of excitons is more challenging and remains to be explored. As neutral composite bosons consisting of electrons and holes, excitons show a non-trivial scattering dynamics. Here, we study on microscopic footing exciton-exciton interaction in transition-metal dichalcogenides and related van der Waals heterostructures. We demonstrate that the crucial criterion for efficient scattering is a large electron/hole mass asymmetry giving rise to internal charge inhomogeneities of excitons and emphasizing their cobosonic substructure. Furthermore, both exchange and direct exciton-exciton interactions are boosted by enhanced exciton Bohr radii. We also predict an unexpected temperature dependence that is usually associated to phonon-driven scattering and we reveal an orders of magnitude stronger interaction of interlayer excitons due to their permanent dipole moment. The developed approach can be generalized to arbitrary material systems and will help to study strongly correlated exciton systems, such as moire super lattices.


title: Systematic error analysis in the absolute hydrogen gas jet polarimeter at RHIC

abstract: The Polarized Atomic Hydrogen Gas Jet Target polarimeter (HJET) is operated at the Relativistic Heavy Ion Collider (RHIC) since 2004 to measure the absolute polarization of each colliding proton beam. Polarimeter detectors and data acquisition were upgraded in 2015 to increase the solid angle, energy range, and to improve the energy and time resolution. These upgrades along with an improved beam intensity and polarization allowed us to greatly reduce the statistical and systematic errors for the proton polarization measurements in RHIC Runs 15 ($E_\text{beam}\!=\!100~\text{GeV}$) and 17 ($255~\text{GeV}$). For a typical 8 hour RHIC store, the measured proton beam average polarization was about $P_\text{beam}\!\sim\!55\pm2.0_\text{stat}\pm0.3_\text{syst}\%$. The elastic $\mathit{pp}$ analyzing power, $A_\text{N}(t)\!\sim\!0.04$, was determined with a precision of about $|\delta A_\text{N}(t)|\!\sim\!0.0002$ in the momentum transfer squared range $0.001\!<\!-t\!<\!0.020\,\text{GeV}^2$. In this paper we present a detailed systematic error analysis of the polarization measurements at HJET. Perspectives of using the HJET based absolute polarimeters in the future Electron Ion Collider (EIC) will be also discussed.


title: A question of Bukh on sums of dilates

abstract: We answer in the affirmative a question of Bukh on the cardinality of the dilate sum $A + 2 \cdot A$.


title: A novel *R-based perspective on solving ordinary differential equations

abstract: The real numbers, it is taught at universities, correspond to our idea of a continuum, although the hyperreal numbers are located ``in between'' the real numbers. The number $x + dx$, where $dx$ should be an infinitesimal number and $x$ real, is infinitesimally close to $x$ but ``infinitely'' far away from all other real numbers. Analogously: If $f'(x_0)$ and $f(x_0)$ are given for a differentiable function $f:\mathbb{R}\rightarrow\mathbb{R}$ at $x_0\in\mathbb{R}$, we can not determine $f(x)$ at {\em any} point $x\in \mathbb{R}$ different from $x_0$. These points seem to be ``infinitely'' far away. That is one conceptual problem of solving differential equations in numerical mathematics. In this article, we will present a numerical algorithm to solve very simple initial value problems. However, the change of paradigm is, that we will not ``leave'' the point $x_0$. Solving ordinary differential equations is like searching for ``recipes'' $f$. Instead of trying to find these recipes for values $x\in\mathbb{R}$, we will learn them from special relations in the ``monad'' of $x_0$.


title: Gauge-invariant quadratic approximation of quasi-local mass and its relation with Hamiltonian for gravitational field

abstract: Gauge invariant, Hamiltonian formulation of field dynamics within a compact region $\Sigma$ with boundary $\partial \Sigma$ is given for the gravitational field linearized over a Kottler metric. The boundary conditions which make the system autonomous are discussed and the corresponding Hamiltonian functional $\mathcal{H}_\text{Invariant}$ is calculated. It is shown that, under specific boundary conditions, the quasi-local Geroch--Hawking mass $\mathcal{H}_\text{Hawking}$ reduces to $\mathcal{H}_\text{Invariant}$ in the weak field approximation. This observation is a quasi-local version of the classical Brill--Deser result.


title: Nearly Linear Row Sampling Algorithm for Quantile Regression

abstract: We give a row sampling algorithm for the quantile loss function with sample complexity nearly linear in the dimensionality of the data, improving upon the previous best algorithm whose sampling complexity has at least cubic dependence on the dimensionality. Based upon our row sampling algorithm, we give the fastest known algorithm for quantile regression and a graph sparsification algorithm for balanced directed graphs. Our main technical contribution is to show that Lewis weights sampling, which has been used in row sampling algorithms for $\ell_p$ norms, can also be applied in row sampling algorithms for a variety of loss functions. We complement our theoretical results by experiments to demonstrate the practicality of our approach.


title: Cosmological Mass of the Photon and Dark Energy as its Bose-Einstein Condensate in de Sitter Space

abstract: I develop a physical picture of dark energy (DE) based on fundamental principles and constants of quantum mechanics (QM) and general relativity (GR) theories. It derives from a conjecture of non-zero masses for nearly standard-model photons, based on QM localization at a cosmological scale. Dark energy is associated with de Sitter space and that has a fundamentally invariant event horizon. I conceive of DE as a Bose-Einstein condensate (BEC) of cosmologically massive photons and I estimate fundamentally the binding energy per particle originating from an effectively attractive QM potential in that BEC. Since massive photons may stand at rest in a de Sitter universe with flat spatial geometry, I solve the time-independent Schroedinger equation for a non-relativistic attractive spherical-well potential self-confining at the de Sitter horizon. The minimal critical potential depth that binds a particle state at the top of that well, combined with the prototypical condition of dark energy-pressure relation in the standard flat $\Lambda$-CDM model, provides an estimate of the photon mass, $m_g$. That is supported by an independent calculation of the vacuum energy of the BEC in a de Sitter static metric with coordinate-time slicing. These estimates provide compatible accounts of dark energy condensation, bridging a chasm between nuclear and cosmological scales. Most notably, I consider a system of cosmological units, or `$g$-units,' that complements the fundamental system of Planck units in various ways. The geometric mean of Planck and $g$-mass turns out to be remarkably close to current estimates of neutrino masses, suggesting that even masses of the lightest known fermions may be deeply related to both GR and QM fundamental constants $\Lambda$, $G$, $c$ and $h$. I finally relate the fine structure constant to $m_g$ and a Coulomb potential confined within a corresponding cosmological horizon.


title: Dynamical Instantons and Activated Processes in Mean-Field Glass Models

abstract: We focus on the energy landscape of a simple mean-field model of glasses and analyze activated barrier-crossing by combining the Kac-Rice method for high-dimensional Gaussian landscapes with dynamical field theory. In particular, we consider Langevin dynamics at low temperature in the energy landscape of the pure spherical $p$-spin model. We select as initial condition for the dynamics one of the many unstable index-1 saddles in the vicinity of a reference local minimum. We show that the associated dynamical mean-field equations admit two solutions: one corresponds to falling back to the original reference minimum, and the other to reaching a new minimum past the barrier. By varying the saddle we scan and characterize the properties of such minima reachable by activated barrier-crossing. Finally, using time-reversal transformations, we construct the two-point function dynamical instanton of the corresponding activated process.


title: Cold Molecular Gas and Free-Free Emission from Hot, Dust-Obscured Galaxies at $z$~3

abstract: We report on observations of redshifted CO(1-0) line emission and observed-frame $\rm\sim$ 30GHz radio continuum emission from five ultra-luminous, mid-IR selected hot, Dust-Obscured Galaxies (Hot DOGs) at $z\rm\gtrsim$ 3 using the Karl G. Jansky Very Large Array. We detect CO(1-0) line emission in all five Hot DOGs, with one of them at high signal to noise. We analyse FIR-radio spectral energy distributions, including dust, free-free and synchrotron emission for the galaxies. We find that most of the 115 GHz rest-frame continuum is mostly due to synchrotron or free-free emission, with only a potentially small contribution from thermal emission. We see a deficit in the rest-frame 115 GHz continuum emission compared to dusty star-forming galaxies (DSFGs) and sub-millimetre galaxies (SMGs) at high redshift, suggesting that Hot DOGs do not have similar cold gas reserves compared with star-forming galaxies. One target, W2305-0039, is detected in the FIRST 1.4 GHz survey, and is likely to possess compact radio jets. We compare to the FIR-radio correlation, and find that at least half of the Hot DOGs in our sample are radio-quiet with respect to normal galaxies. These findings suggest that Hot DOGs have comparably less cold molecular gas than star-forming galaxies at lower, $z\rm\sim$ 2 redshifts, and are dominated by powerful, yet radio-quiet AGN.


title: Propargylimine in the laboratory and in space: millimetre-wave spectroscopy and first detection in the ISM

abstract: Small imines containing up to three carbon atoms are present in the interstellar medium. As alkynyl compounds are abundant in this medium, propargylimine thus represents a promising candidate for a new interstellar detection. The goal of the present work is to perform a comprehensive laboratory investigation of the rotational spectrum of propargylimine in its ground vibrational state in order to obtain a highly precise set of rest-frequencies and to search it in space. The rotational spectra of $E$ and $Z$ geometrical isomers of propargylimine have been recorded in laboratory in the 83-500 GHz frequency interval. The measurements have been performed using a source-modulation millimetre-wave spectrometer equipped with a pyrolysis system for the production of unstable species. High-level ab initio calculations were performed to assist the analysis and to obtain reliable estimates for an extended set of spectroscopic quantities. The improved spectral data allow us to perform a successful search for this new imine in the quiescent G+0.693-0.027 molecular cloud. Eighteen lines of $Z$-propargylimine have been detected at level $>2.5 \sigma$, resulting in a column density estimate of $N = (0.24\pm 0.02)\times 10^{14}$ cm$^{-2}$. An upper limit was retrieved for the higher-energy $E$ isomer, which was not detected in the data. The fractional abundance (w.r.t. H$_2$) derived for $Z$-propargylimine is $1.8\times 10^{-10}$. We discuss the possible formation routes by comparing the derived abundance with those measured in the source for possible chemical precursors.


title: Assumption-lean inference for generalised linear model parameters

abstract: Inference for the parameters indexing generalised linear models is routinely based on the assumption that the model is correct and a priori specified. This is unsatisfactory because the chosen model is usually the result of a data-adaptive model selection process, which may induce excess uncertainty that is not usually acknowledged. Moreover, the assumptions encoded in the chosen model rarely represent some a priori known, ground truth, making standard inferences prone to bias, but also failing to give a pure reflection of the information that is contained in the data. Inspired by developments on assumption-free inference for so-called projection parameters, we here propose novel nonparametric definitions of main effect estimands and effect modification estimands. These reduce to standard main effect and effect modification parameters in generalised linear models when these models are correctly specified, but have the advantage that they continue to capture respectively the primary (conditional) association between two variables, or the degree to which two variables interact (in a statistical sense) in their effect on outcome, even when these models are misspecified. We achieve an assumption-lean inference for these estimands (and thus for the underlying regression parameters) by deriving their influence curve under the nonparametric model and invoking flexible data-adaptive (e.g., machine learning) procedures.


title: On the Loss Landscape of Adversarial Training: Identifying Challenges and How to Overcome Them

abstract: We analyze the influence of adversarial training on the loss landscape of machine learning models. To this end, we first provide analytical studies of the properties of adversarial loss functions under different adversarial budgets. We then demonstrate that the adversarial loss landscape is less favorable to optimization, due to increased curvature and more scattered gradients. Our conclusions are validated by numerical analyses, which show that training under large adversarial budgets impede the escape from suboptimal random initialization, cause non-vanishing gradients and make the model find sharper minima. Based on these observations, we show that a periodic adversarial scheduling (PAS) strategy can effectively overcome these challenges, yielding better results than vanilla adversarial training while being much less sensitive to the choice of learning rate.


title: Generators in $\mathcal{Z}$-stable C*-algebras of real rank zero

abstract: We show that every separable C*-algebra of real rank zero that tensorially absorbs the Jiang-Su algebra contains a dense set of generators. It follows that in every classifiable, simple, nuclear C*-algebra, a generic element is a generator.


title: Analytic insights on the information content of new observables

abstract: Uncertainty quantification has emerged as a rapidly growing field in nuclear science. Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experiments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of limited value. Also, one often deals with theoretical constructs that involve fundamental quantities that are not accessible to experiment or observation. Particularly relevant in this context is the pressure of pure neutron matter. In this contribution we develop an analytic framework to answer the question of "How can new data reduce uncertainties of current theoretical models?" [P.-G. Reinhard and W. Nazarewicz, Phys. Rev. C81, 051303(R) (2010)]. Simple and insightful expressions are obtained to quantify the impact of one or two new observables on theoretical uncertainties in two critical quantities: the slope of the symmetry energy at saturation density and the pressure of pure neutron matter at twice nuclear matter saturation density.


title: Lerch's $\Phi$ and the Polylogarithm at the Positive Integers

abstract: We review the closed-forms of the partial Fourier sums associated with $HP_k(n)$ and create an asymptotic expression for $HP(n)$ as a way to obtain formulae for the full Fourier series (if $b$ is such that $|b|<1$, we get a surprising pattern, $HP(n) \sim H(n)-\sum_{k\ge 2}(-1)^k\zeta(k)b^{k-1}$). Finally, we use the found Fourier series formulae to obtain the values of the Lerch transcendent function, $\Phi(e^m,k,b)$, and by extension the polylogarithm, $\mathrm{Li}_{k}(e^{m})$, at the positive integers $k$.


title: Submillimetre and far-infrared spectroscopy of monodeuterated amidogen radical (NHD): improved rest-frequencies for astrophysical observations

abstract: Observations of ammonia in interstellar environments have revealed high levels of deuteration, and all its D-containing variants, including ND$_3$, have been detected in cold prestellar cores and around young protostars. The observation of these deuterated isotopologues is very useful to elucidate the chemical and physical processes taking place during the very early stages of star formation, as the abundance of deuterated molecules is highly enhanced in dense and cold gas. Nitrogen hydride radicals are key species lying at the very beginning of the reaction pathway leading to the formation of NH$_3$ and organic molecules of pre-biotic interest, but relatively little information is known about their D-bearing isotopologues. To date, only ND has been detected in the interstellar gas. To aid the identification of further deuterated nitrogen radicals, we have thoroughly re-investigated the rotational spectrum of NHD employing two different instruments: a frequency-modulation submillimetre spectrometer operating in the THz region and a synchrotron-based Fourier Transform infrared spectrometer operating in the 50-240 cm$^{-1}$ wavelength range. NHD was produced in a plasma of NH$_3$ and D$_2$. A wide range of rotational energy levels has been probed thanks to the observation of high $N$ (up to 15) and high $K_a$ (up to 9) transitions. A global analysis including our new data and those already available in the literature has provided a comprehensive set of very accurate spectroscopic parameters. A highly reliable line catalogue has been generated to assist archival data searches and future astronomical observations of NHD at submillimetre and THz regimes.


title: Computing Shortest Paths and Diameter in the Hybrid Network Model

abstract: The $\mathsf{HYBRID}$ model, introduced in [Augustine et al., SODA '20], provides a theoretical foundation for networks that allow multiple communication modes. The model follows the principles of synchronous message passing, whereas nodes are allowed to use \textit{two} fundamentally different communication modes. First, a local mode where nodes may exchange arbitrary information per round over edges of a local communication graph $G$ (akin to the $\mathsf{LOCAL}$ model). Second, a global mode where every node may exchange $O(\log n)$ messages of size $O(\log n)$ bits per round with arbitrary nodes in the network. The $\mathsf{HYBRID}$ model intends to reflect the conditions of many real hybrid networks, where high-bandwidth but inherently local communication is combined with highly flexible global communication with restricted bandwidth. We continue to explore the power and limitations of the $\mathsf{HYBRID}$ model by investigating the complexity of computing shortest paths and diameter of the local communication graph $G$. We improve on the known upper bound for the exact all pairs shortest paths problem (APSP) from [Augustine et al., SODA '20] and provide algorithms to approximate solutions for the $k$ source shortest paths problem ($k$-SSSP). We demonstrate that our results for APSP and $k$-SSP are almost tight (up to poly-logarithmic factors). Furthermore, we give an improved algorithm for the exact single source shortest paths problem for graphs with large diameter. For the problem of approximating the diameter of the local communication network we give the first non-trivial upper bound. This upper bound is complemented by a lower bound for the exact diameter problem.


title: Machine Common Sense

abstract: Machine common sense remains a broad, potentially unbounded problem in artificial intelligence (AI). There is a wide range of strategies that can be employed to make progress on this challenge. This article deals with the aspects of modeling commonsense reasoning focusing on such domain as interpersonal interactions. The basic idea is that there are several types of commonsense reasoning: one is manifested at the logical level of physical actions, the other deals with the understanding of the essence of human-human interactions. Existing approaches, based on formal logic and artificial neural networks, allow for modeling only the first type of common sense. To model the second type, it is vital to understand the motives and rules of human behavior. This model is based on real-life heuristics, i.e., the rules of thumb, developed through knowledge and experience of different generations. Such knowledge base allows for development of an expert system with inference and explanatory mechanisms (commonsense reasoning algorithms and personal models). Algorithms provide tools for a situation analysis, while personal models make it possible to identify personality traits. The system so designed should perform the function of amplified intelligence for interactions, including human-machine.


title: Mukai's Program (reconstructing a K3 surface from a curve) via wall-crossing, II

abstract: Let $C$ be a curve on a K3 surface $X$ with Picard group $\mathbb{Z}.[C]$. Mukai's program seeks to recover $X$ from $C$ by exhibiting it as a Fourier-Mukai partner to a Brill-Noether locus of vector bundles on $C$. We use wall-crossing in the space of Bridgeland stability conditions to prove this for genus $\ge14$. This paper deals with the case $g-1$ prime left over from Paper I.


title: Reanalysis of low-energy electron-proton scattering data and proton radius

abstract: We reanalyze electron-proton scattering data in the momentum transfer region $Q^2 < 1\ GeV^2$, which were used to determine proton electric radius, with main focus on the Mainz experiment (Bernauer et al., 2010). We examine dependence of fit results and $\chi^2$ on the cut-off in $Q^2$ and degree of fitting polynomial, using pseudo-data and compare it with the case of real experimental data. We find that under some assumptions, the data could be consistent with the radius, obtained from muonic hydrogen.


title: Bubbly drag reduction using a hydrophobic inner cylinder in Taylor-Couette turbulence

abstract: In this study we experimentally investigate bubbly drag reduction in a highly turbulent flow of water with dispersed air at $5.0 \times 10^{5} \leq \text{Re} \leq 1.7 \times 10^{6}$ over a non-wetting surface containing micro-scale roughness. To do so, the Taylor-Couette geometry is used, allowing for both accurate global drag and local flow measurements. The inner cylinder - coated with a rough, hydrophobic material - is rotating, whereas the smooth outer cylinder is kept stationary. The crucial control parameter is the air volume fraction $\alpha$ present in the working fluid. For small volume fractions ($\alpha < {4}\,\%$), we observe that the surface roughness from the coating increases the drag. For large volume fractions of air ($\alpha \geq 4\,\%$), the drag decreases compared to the case with both the inner and outer cylinders uncoated, i.e. smooth and hydrophilic, using the same volume fraction of air. This suggests that two competing mechanisms are at place: on the one hand the roughness invokes an extension of the log-layer - resulting in an increase in drag - and on the other hand there is a drag-reducing mechanism of the hydrophobic surface interacting with the bubbly liquid. The balance between these two effects determines whether there is overall drag reduction or drag enhancement. For further increased bubble concentration $\alpha = {6}\,\%$ we find a saturation of the drag reduction effect. Our study gives guidelines for industrial applications of bubbly drag reduction in hydrophobic wall-bounded turbulent flows.


title: Reciprocal Adversarial Learning via Characteristic Functions

abstract: Generative adversarial nets (GANs) have become a preferred tool for tasks involving complicated distributions. To stabilise the training and reduce the mode collapse of GANs, one of their main variants employs the integral probability metric (IPM) as the loss function. This provides extensive IPM-GANs with theoretical support for basically comparing moments in an embedded domain of the \textit{critic}. We generalise this by comparing the distributions rather than their moments via a powerful tool, i.e., the characteristic function (CF), which uniquely and universally comprising all the information about a distribution. For rigour, we first establish the physical meaning of the phase and amplitude in CF, and show that this provides a feasible way of balancing the accuracy and diversity of generation. We then develop an efficient sampling strategy to calculate the CFs. Within this framework, we further prove an equivalence between the embedded and data domains when a reciprocal exists, where we naturally develop the GAN in an auto-encoder structure, in a way of comparing everything in the embedded space (a semantically meaningful manifold). This efficient structure uses only two modules, together with a simple training strategy, to achieve bi-directionally generating clear images, which is referred to as the reciprocal CF GAN (RCF-GAN). Experimental results demonstrate the superior performances of the proposed RCF-GAN in terms of both generation and reconstruction.


title: Infrared wavefront sensing for adaptive optics assisted Galactic Center observations with the VLT interferometer and GRAVITY: operation and results

abstract: This article describes the operation of the near-infrared wavefront sensing based Adaptive Optics (AO) system CIAO. The Coud\'e Infrared Adaptive Optics (CIAO) system is a central auxiliary component of the Very Large Telescope (VLT) interferometer (VLTI). It enables in particular the observations of the Galactic Center (GC) using the GRAVITY instrument. GRAVITY is a highly specialized beam combiner, a device that coherently combines the light of the four 8-m telescopes and finally records interferometric measurements in the K-band on 6 baselines simultaneously. CIAO compensates for phase disturbances caused by atmospheric turbulence, which all four 8 m Unit Telescopes (UT) experience during observation. Each of the four CIAO units generates an almost diffraction-limited image quality at its UT, which ensures that maximum flux of the observed stellar object enters the fibers of the GRAVITY beam combiner. We present CIAO performance data obtained in the first 3 years of operation as a function of weather conditions. We describe how CIAO is configured and used for observations with GRAVITY. In addition, we focus on the outstanding features of the near-infrared sensitive Saphira detector, which is used for the first time on Paranal, and show how it works as a wavefront sensor detector.


title: Does Cascading Schmitt-Trigger Stages Improve the Metastable Behavior?

abstract: Schmitt-Trigger stages are the method of choice for robust discretization of input voltages with excessive transition times or significant noise. However, they may suffer from metastability. Based on the experience that the cascading of flip-flop stages yields a dramatic improvement of their overall metastability hardness, in this paper we elaborate on the question whether the cascading of Schmitt-Trigger stages can obtain a similar gain. We perform a theoretic analysis that is backed up by an existing metastability model for a single Schmitt-Trigger stage and elaborate some claims about the behavior of a Schmitt-Trigger cascade. These claims suggest that the occurrence of metastability is indeed reduced from the first stage to the second which suggests an improvement. On the downside, however, it becomes clear that metastability can still not be completely ruled out, and in some cases the behavior of the cascade may be less beneficial for a given application, e.g. by introducing seemingly acausal transitions. We validate our findings by extensive HSPICE simulations in which we directly cover our most important claims.


title: The Limiting Poisson Law of Massive MIMO Detection with Box Relaxation

abstract: Estimating a binary vector from noisy linear measurements is a prototypical problem for MIMO systems. A popular algorithm, called the box-relaxation decoder, estimates the target signal by solving a least squares problem with convex constraints. This paper shows that the performance of the algorithm, measured by the number of incorrectly-decoded bits, has a limiting Poisson law. This occurs when the sampling ratio and noise variance, two key parameters of the problem, follow certain scalings as the system dimension grows. Moreover, at a well-defined threshold, the probability of perfect recovery is shown to undergo a phase transition that can be characterized by the Gumbel distribution. Numerical simulations corroborate these theoretical predictions, showing that they match the actual performance of the algorithm even in moderate system dimensions.


title: Searching for the Shortest Path to the Point of Voltage Collapse on the Algebraic Manifold

abstract: Voltage instability is one of the main causes of power system blackouts. Emerging technologies such as renewable energy integration, distributed energy resources and demand responses may introduce significant uncertainties in analyzing of system-wide voltage stability. This paper starts with summarizing different known voltage instability mechanisms, and then focuses on a class of voltage instability which is induced by the singular surface of the algebraic manifold. We argue and demonstrate that this class can include both dynamic and static voltage instabilities. To determine the minimum distance to the point of voltage collapse, a new formulation is proposed on the algebraic manifold. This formulation is further converted into an optimal control framework for identifying the path with minimum distance on the manifold. Comprehensive numerical studies are conducted on some manifolds of different power system test cases and demonstrate that the proposed method yields candidates for the local shortest paths to the singular surface on the manifold for both the dynamic model and the static model. Simulations show that the proposed method can identify shorter paths on the manifold than the paths associated with the minimum Euclidean distances. Furthermore, the proposed method always locates the right path ending at the correct singular surface which is responsible for the voltage instability; while the Euclidean distance formulation can mistakenly find solutions on the wrong singular surface. A broad range of potential applications using the proposed method are also discussed.


title: A symmetric function of increasing forests

abstract: For an indifference graph $G$ we define a symmetric function of increasing spanning forests of $G$. We prove that this symmetric function satisfies certain linear relations, which are also satisfied by the chromatic quasisymmetric function and unicellular LLT polynomials. As a consequence we give a combinatorial interpretation of the coefficients of the LLT polynomial in the elementary basis (up to a factor of a power of $(q-1)$), strengthening the description given by Alexandersson and Sulzgruber.


title: Spherical Motion Dynamics: Learning Dynamics of Neural Network with Normalization, Weight Decay, and SGD

abstract: In this work, we comprehensively reveal the learning dynamics of neural network with normalization, weight decay (WD), and SGD (with momentum), named as Spherical Motion Dynamics (SMD). Most related works study SMD by focusing on "effective learning rate" in "equilibrium" condition, where weight norm remains unchanged. However, their discussions on why equilibrium condition can be reached in SMD is either absent or less convincing. Our work investigates SMD by directly exploring the cause of equilibrium condition. Specifically, 1) we introduce the assumptions that can lead to equilibrium condition in SMD, and prove that weight norm can converge at linear rate with given assumptions; 2) we propose "angular update" as a substitute for effective learning rate to measure the evolving of neural network in SMD, and prove angular update can also converge to its theoretical value at linear rate; 3) we verify our assumptions and theoretical results on various computer vision tasks including ImageNet and MSCOCO with standard settings. Experiment results show our theoretical findings agree well with empirical observations.


title: Combinatorial Group Testing and Sparse Recovery Schemes with Near-Optimal Decoding Time

abstract: In the long-studied problem of combinatorial group testing, one is asked to detect a set of $k$ defective items out of a population of size $n$, using $m \ll n$ disjunctive measurements. In the non-adaptive setting, the most widely used combinatorial objects are disjunct and list-disjunct matrices, which define incidence matrices of test schemes. Disjunct matrices allow the identification of the exact set of defectives, whereas list disjunct matrices identify a small superset of the defectives. Apart from the combinatorial guarantees, it is often of key interest to equip measurement designs with efficient decoding algorithms. The most efficient decoders should run in sublinear time in $n$, and ideally near-linear in the number of measurements $m$. In this work, we give several constructions with an optimal number of measurements and near-optimal decoding time for the most fundamental group testing tasks, as well as for central tasks in the compressed sensing and heavy hitters literature. For many of those tasks, the previous measurement-optimal constructions needed time either quadratic in the number of measurements or linear in the universe size. Most of our results are obtained via a clean and novel approach which avoids list-recoverable codes or related complex techniques which were present in almost every state-of-the-art work on efficiently decodable constructions of such objects.


title: An approximation algorithm for joint caching and recommendations in cache networks

abstract: Streaming platforms, like Netflix and YouTube, strive to offer high streaming quality (SQ), in terms of bitrate, delays, etc., to their users. Meanwhile, a significant share of content consumption of these platforms is heavily influenced by recommendations. In this setting, the user's overall experience is a product of both the user's interest in a recommended content, i.e., the recommendation quality (RQ), and the SQ of this content. However, network decisions (like caching) that affect the SQ are usually made without considering the recommender's actions. Likewise, recommendations are chosen independently of the potential delivery quality. In this paper, we define a metric of streaming experience (MoSE) that captures the fundamental tradeoff between the SQ and RQ. We aim to jointly optimize caching and recommendations in a generic network of caches, with the objective of maximizing this metric. This is in line with the recent trend for content providers to simultaneously act as Content Delivery Network owners, implying that the same entity may handle both caching and recommendation decisions. We formulate this joint optimization problem and prove that it can be approximated up to a constant factor. To the best of our knowledge, this is the first polynomial algorithm to achieve a constant approximation ratio for the joint problem. Moreover, our numerical experiments show important performance gains of our algorithm over baseline schemes and existing algorithms in the literature.


title: Starlink Satellite Brightness Before VisorSat

abstract: The mean of 830 visual magnitudes adjusted to a distance of 550 km (the operational altitude) is 4.63 +/-0.02. The data on DarkSat, the low-albedo satellite, indicate that it is fainter than the others by 1.6 magnitudes or 78%. However, there is considerable uncertainty in this value due to the small number of observations. Some satellites were observed to flare by 10,000 times their normal brightness. These statistics can serve as a baseline for assessing the reduced brightness of the VisorSat design for future Starlink satellites.


title: QCD traveling waves phenomenology revisited

abstract: In this paper we review and update the Amaral-Gay Ducati-Betemps-Soyez saturation model, by testing it against the recent H1-ZEUS combined data on deep inelastic scattering, including heavy quarks in the dipole amplitude. We obtain that this model, which is based on traveling wave solutions of the Balitsky-Kovchegov equation and built in the momentum space framework, yields very accurate descriptions of the reduced cross section, $\sigma_{r}(x,y,Q^{2})$, as well as DIS structure functions such as $F_{2}(x,Q^{2})$ and $F_{L}(x,Q^{2})$, all measured at HERA. Additionally, it provides good descriptions of heavy quark structure functions, $F_{2}^{cc}$ and $F_{2}^{bb}$ at small-$x$ and $Q^{2}\lesssim 60$ GeV$^{2}$. We also use the improved model to make predictions for structure functions to be measured in the near future at LHeC.


title: Spectroscopic factors, overlaps, and isospin symmetry from an $R$-matrix point of view

abstract: Background: Spectroscopic factors, overlaps, and isospin symmetry are often used in conjunction with single-particle wave functions for the phenomenological analysis of nuclear structure and reactions. Many differing prescriptions for connecting these quantities to physically relevant asymptotic normalization constants or widths are available in the literature, but their relationship and degree of validity are not always clear. Purpose: This paper derives relationships among the above quantities of interest using well-defined methodology and starting assumptions. Method: $R$-matrix theory is used as the primary tool to interoperate between the quantities of interest to this work. Particular attention is paid to effects arising from beyond the nuclear surface, where isospin symmetry is strongly violated. Results: Relationships among the quantities of interest are derived. Example applications of these methods to mirror levels in nucleon+${}^{12}{\rm C}$, nucleon+${}^{16}{\rm O}$, and nucleon+${}^{26}{\rm Al}$ are presented. A new approach to multi-level mirror symmetry is derived and applied to the first three $2^+$ states of ${}^{18}{\rm O}$ and ${}^{18}{\rm Ne}$. Conclusions: The relationship between the quantities of interest is clarified and certain procedures are recommended. It is found that the asymptotic normalization constant of the second $2^+$ state in ${}^{18}{\rm Ne}$ deduced from the mirror state in ${}^{18}{\rm O}$ is significantly larger than found in previous work. This finding has the effect of increasing the ${}^{17}{\rm F}(p,\gamma){}^{18}{\rm Ne}$ reaction rate in novae.


title: Finding the Loops that Matter

abstract: The Loops that Matter method (Schoenberg et. al, 2019) for understanding model behavior provides metrics showing the contribution of the feedback loops in a model to behavior at each point in time. To provide these metrics, it is necessary find the set of loops on which to compute them. We show in this paper the necessity of including loops that are important at different points in the simulation. These important loops may not be independent of one another and cannot be determined from static analysis of the model structure. We then describe an algorithm that can be used to discover the most important loops in models that are too feedback rich for exhaustive loop discovery. We demonstrate the use of this algorithm in terms of its ability to find the most explanatory loops, and its computational performance for large models. By using this approach, the Loops that Matter method can be applied to models of any size or complexity.


title: Walking in the Shadow: A New Perspective on Descent Directions for Constrained Minimization

abstract: Descent directions such as movement towards Descent directions, including movement towards Frank-Wolfe vertices, away-steps, in-face away-steps and pairwise directions, have been an important design consideration in conditional gradient descent (CGD) variants. In this work, we attempt to demystify the impact of the movement in these directions towards attaining constrained minimizers. The optimal local direction of descent is the directional derivative (i.e., shadow) of the projection of the negative gradient. We show that this direction is the best away-step possible, and the continuous-time dynamics of moving in the shadow is equivalent to the dynamics of projected gradient descent (PGD), although it's non-trivial to discretize. We also show that Frank-Wolfe (FW) vertices correspond to projecting onto the polytope using an "infinite" step in the direction of the negative gradient, thus providing a new perspective on these steps. We combine these insights into a novel Shadow-CG method that uses FW and shadow steps, while enjoying linear convergence, with a rate that depends on the number of breakpoints in its projection curve, rather than the pyramidal width. We provide a linear bound on the number of breakpoints for simple polytopes and present scaling-invariant upper bounds for general polytopes based on the number of facets. We exemplify the benefit of using Shadow-CG computationally for various applications, while raising an open question about tightening the bound on the number of breakpoints for general polytopes.


title: Random graphs with arbitrary clustering and their applications

abstract: The structure of many real networks is not locally tree-like and hence, network analysis fails to characterise their bond percolation properties. In a recent paper [P. Mann, V. A. Smith, J. B. O. Mitchell, and S. Dobson, Percolation in random graphs with higher-order clustering, arXiv e-prints, p. arXiv:2006.06744, June 2020.], we developed analytical solutions to the percolation properties of random networks with homogeneous clustering (clusters whose nodes are degree-equivalent). In this paper, we extend this model to investigate networks that contain clusters whose nodes are not degree-equivalent, including multilayer networks. Through numerical examples we show how this method can be used to investigate the properties of random complex networks with arbitrary clustering, extending the applicability of the configuration model and generating function formulation.


title: Temperature-Induced Lifshitz Transition and Possible Excitonic Instability in ZrSiSe

abstract: The nodal-line semimetals have attracted immense interest due to the unique electronic structures such as the linear dispersion and the vanishing density of states as the Fermi energy approaching the nodes. Here, we report temperature-dependent transport and scanning tunneling microscope (spectroscopy) (STM[S]) measurements on nodal-line semimetal ZrSiSe.Our experimental results and theoretical analyses consistently demonstrate that the temperature induces Lifshitz transitions at 80 and 106 K in ZrSiSe, which results in the transport anomalies at the same temperatures. More strikingly, we observe a V-shaped dip structure around Fermi energy from the STS spectrum at low temperature,which can be attributed to co-effect of the spin-orbit coupling and excitonic instability. Our observations indicate the correlation interaction may play an important role in ZrSiSe, which owns the quasi-two-dimensional electronic structures.


title: Generating Reliable and Efficient Predictions of Human Motion: A Promising Encounter between Physics and Neural Networks

abstract: Generating accurate and efficient predictions for the motion of the humans present in the scene is key to the development of effective motion planning algorithms for robots moving in promiscuous areas, where wrong planning decisions could generate safety hazard or simply make the presence of the robot "socially" unacceptable. Our approach to predict human motion is based on a neural network of a peculiar kind. Contrary to conventional deep neural networks, our network embeds in its structure the popular Social Force Model, a dynamic equation describing the motion in physical terms. This choice allows us to concentrate the learning phase in the aspects, which are really unknown (i.e., the model's parameters) and to keep the structure of the network simple and manageable. As a result, we are able to obtain a good prediction accuracy with a small synthetically generated training set, and the accuracy remains acceptable even when the network is applied in scenarios quite different from those for which it was trained. Finally, the choices of the network are "explainable", as they can be interpreted in physical terms. Comparative and experimental results prove the effectiveness of the proposed approach.


title: Multimode cold-damping optomechanics with delayed feedback

abstract: We investigate the role of time delay in cold-damping optomechanics with multiple mechanical resonances. For instantaneous electronic response, it was recently shown in \textit{Phys. Rev. Lett. \textbf{123}, 203605 (2019)}, that a single feedback loop is sufficient to simultaneously remove thermal noise from many mechanical modes. While the intrinsic delayed response of the electronics can induce single mode and mutual heating between adjacent modes, we propose to counteract such detrimental effects by introducing an additional time delay to the feedback loop. For lossy cavities and broadband feedback, we derive analytical results for the final occupancies of the mechanical modes within the formalism of quantum Langevin equations. For modes that are frequency degenerate collective effects dominate, mimicking behavior similar to Dicke super- and subradiance. These analytical results, corroborated with numerical simulations of both transient and steady state dynamics, allow to find suitable conditions and strategies for efficient single or multimode feedback optomechanics.


title: General Relativistic Mean-field Dynamo Model for Proto-neutron Stars

abstract: Neutron stars, and magnetars in particular, are known to host the strongest magnetic fields in the Universe. The origin of these strong fields is a matter of controversy. In this preliminary work, via numerical simulations, we study, for the first time in non-ideal general relativistic magnetohydrodynamic (GRMHD) regime, the growth of the magnetic field due to the action of the mean-field dynamo due to sub-scale, unresolved turbulence. The dynamo process, combined with the differential rotation of the (proto-)star, is able to produce an exponential growth of any initial magnetic seed field up to the values required to explain the observations. By varying the dynamo coefficient we obtain different growth rates. We find a quasi-linear dependence of the growth rates on the intensity of the dynamo. Furthermore, the time interval in which exponential growth occurs and the growth rates also seems to depend on the initial configuration of the magnetic field.


title: SD-RSIC: Summarization Driven Deep Remote Sensing Image Captioning

abstract: Deep neural networks (DNNs) have been recently found popular for image captioning problems in remote sensing (RS). Existing DNN based approaches rely on the availability of a training set made up of a high number of RS images with their captions. However, captions of training images may contain redundant information (they can be repetitive or semantically similar to each other), resulting in information deficiency while learning a mapping from the image domain to the language domain. To overcome this limitation, in this paper, we present a novel Summarization Driven Remote Sensing Image Captioning (SD-RSIC) approach. The proposed approach consists of three main steps. The first step obtains the standard image captions by jointly exploiting convolutional neural networks (CNNs) with long short-term memory (LSTM) networks. The second step, unlike the existing RS image captioning methods, summarizes the ground-truth captions of each training image into a single caption by exploiting sequence to sequence neural networks and eliminates the redundancy present in the training set. The third step automatically defines the adaptive weights associated to each RS image to combine the standard captions with the summarized captions based on the semantic content of the image. This is achieved by a novel adaptive weighting strategy defined in the context of LSTM networks. Experimental results obtained on the RSCID, UCM-Captions and Sydney-Captions datasets show the effectiveness of the proposed approach compared to the state-of-the-art RS image captioning approaches. The code of the proposed approach is publicly available at https://gitlab.tubit.tu-berlin.de/rsim/SD-RSIC.


title: Calibration of the von Wolffersdorff model using Genetic Algorithms

abstract: This article proposes an optimization framework, based on Genetic Algorithms (GA), to calibrate the constitutive law of von Wolffersdorff. This constitutive law is known as Sand Hypoplasticity (SH), and allows for robust and accurate modeling of the soil behavior but requires a complex calibration involving eight parameters. The proposed optimization can automatically fit these parameters from the results of an oedometric and a triaxial drained compression test, by combining the GA with a numerical solver that integrates the SH in the test conditions. By repeating the same calibration several times, the stochastic nature of the optimizer enables the uncertainty quantification of the calibration parameters and allows studying their relative importance on the model prediction. After validating the numerical solver on the ExCaliber-Laboratory software from the SoilModels' website, the GA calibration is tested on a synthetic dataset to analyze the convergence and the statistics of the results. In particular, a correlation analysis reveals that two couples of the eight model parameters are strongly correlated. Finally, the calibration procedure is tested on the results from von Wolffersdorff, 1996, and Herle & Gudehus, 1999, on the Hochstetten sand. The model parameters identified by the Genetic Algorithm optimization improves the matching with the experimental data and hence lead to a better calibration.


title: An efficient application of Bayesian optimization to an industrial MDO framework for aircraft design

abstract: The multi-level, multi-disciplinary and multi-fidelity optimization framework developed at Bombardier Aviation has shown great results to explore efficient and competitive aircraft configurations. This optimization framework has been developed within the Isight software, the latter offers a set of ready-to-use optimizers. Unfortunately, the computational effort required by the Isight optimizers can be prohibitive with respect to the requirements of an industrial context. In this paper, a constrained Bayesian optimization optimizer, namely the super efficient global optimization with mixture of experts, is used to reduce the optimization computational effort. The obtained results showed significant improvements compared to two of the popular Isight optimizers. The capabilities of the tested constrained Bayesian optimization solver are demonstrated on Bombardier research aircraft configuration study cases.


title: Efficient Ab-Initio Molecular Dynamic Simulations by Offloading Fast Fourier Transformations to FPGAs

abstract: A large share of today's HPC workloads is used for Ab-Initio Molecular Dynamics (AIMD) simulations, where the interatomic forces are computed on-the-fly by means of accurate electronic structure calculations. They are computationally intensive and thus constitute an interesting application class for energy-efficient hardware accelerators such as FPGAs. In this paper, we investigate the potential of offloading 3D Fast Fourier Transformations (FFTs) as a critical routine of plane-wave-based electronic structure calculations to FPGA and in conjunction demonstrate the tolerance of these simulations to lower precision computations.


title: Implementation of the orthodoxy test as a validity check on experimental field emission data

abstract: In field electron emission (FE) studies, it is important to check and analyse the quality and validity of experimental current-voltage data, which is usually plotted in one of a small number of standard forms. These include the so-called Fowler-Nordheim (FN), Millikan-Lauritsen (ML) and Murphy-Good (MG) plots. The Field Emission Orthodoxy Test is a simple quantitative test that aims to check for the reasonableness of the values of the parameter "scaled field" that can be extracted from these plots. This is done in order to establish whether characterization parameters extracted from the plot will be reliable or, alternative, likely to be spurious. This paper summarises the theory behind the orthodoxy test, for each of the plot forms, and confirms that it is easy to apply it to the newly developed MG plot. A simple web tool has been developed that extracts scaled-field values from any of these three plot forms, and tests for lack of field emission orthodoxy.


title: Depth Uncertainty in Neural Networks

abstract: Existing methods for estimating uncertainty in deep learning tend to require multiple forward passes, making them unsuitable for applications where computational resources are limited. To solve this, we perform probabilistic reasoning over the depth of neural networks. Different depths correspond to subnetworks which share weights and whose predictions are combined via marginalisation, yielding model uncertainty. By exploiting the sequential structure of feed-forward networks, we are able to both evaluate our training objective and make predictions with a single forward pass. We validate our approach on real-world regression and image classification tasks. Our approach provides uncertainty calibration, robustness to dataset shift, and accuracies competitive with more computationally expensive baselines.


title: A practical model of twin-beam experiments for sub-shot-noise absorption measurements

abstract: Quantum-intensity-correlated twin beams of light can be used to measure absorption with precision beyond the classical shot-noise limit. The degree to which this can be achieved with a given estimator is defined by the quality of the twin-beam intensity correlations, which is quantified by the noise reduction factor. We derive an analytical model of twin-beam experiments, incorporating experimental parameters such as the relative detection efficiency of the beams, uncorrelated optical noise, and uncorrelated detector noise. We show that for twin beams without excessive noise, measured correlations can be improved by increasing the detection efficiency of each beam, notwithstanding this may unbalance detection efficiency. However, for beams with excess intensity or other experimental noise, one should balance detection efficiency, even at the cost of reducing detection efficiency -- we specifically define these noise conditions and verify our results with statistical simulation. This has application in design and optimization of absorption spectroscopy and imaging experiments.


title: Initial-boundary value problem for a time-fractional subdiffusion equation with an arbitrary elliptic differential operator

abstract: An initial-boundary value problem for a time-fractional subdiffusion equation with an arbitrary order elliptic differential operator is considered. Uniqueness and existence of the classical solution of the posed problem are proved by the classical Fourier method. Sufficient conditions for the initial function and for the right-hand side of the equation are indicated, under which the corresponding Fourier series converge absolutely and uniformly. In the case of an initial-boundary value problem on N-dimensional torus, one can easily see that these conditions are not only sufficient, but also necessary.


title: The permanent tide and the International Height Reference Frame IHRF

abstract: The International Height Reference System (IHRS) adopted by the IAG in 2015 contains two novelties. First, the mean-tide system is adopted for handling the permanent tide. Second, the reference level of height system is defined by the equipotential surface where the geopotential has a conventional value W0=62636853.4 m2s-2. This value was first determined empirically to provide a good approximation to the global mean sea level and then adopted as a reference value by convention. I analyse the tidal aspects of the reference level based on W0. The W0 is by definition independent of the tidal system adopted for the equipotential surface, but for different tidal systems, different functions are involved in the W of the equation W=W0. I find that the empirical determination of the adopted estimate W0 is inconsistent from the viewpoint of tidal systems. However, the adopted estimate and the consistent estimate from the same data round off to the same value. I propose a simplified definition of IHRF geopotential numbers that would make it possible to transform between IHRF and zero-tide geopotential numbers using a simple datum-difference surface. IHRF should adopt a conventional (best) estimate of the permanent tidal potential such as is contained in the IERS Conventions and use it as a basis for other conventional formulas. The tide-free coordinates of the ITRF, and tide-free Global Geopotential Models are central in the modelling of geopotential for the purposes of the IHRF. I present a set of correction formulas that can be used to move to the zero-tide model before, in the middle, or after the processing, and finally to the mean-tide IHRF. To reduce the confusion around the multitude of tidal systems I propose that modelling should primarily done in the zero-tide system, with the mean-tide potential as an add-on.


title: A self-lensing supermassive binary black hole at radio frequencies: the story of Spikey continues

abstract: The quasar J1918+4937 was recently suggested to harbour a milliparsec-separation binary supermassive black hole (SMBH), based upon modeling the narrow spike in its high-cadence Kepler optical light curve. Known binary SMBHs are extremely rare, and the tight constraints on the physical and geometric parameters of this object are unique. The high-resolution radio images of J1918+4937 obtained with very long baseline interferometry (VLBI) indicate a rich one-sided jet structure extending to 80 milliarcseconds. Here we analyse simultaneously-made sensitive 1.7- and 5-GHz archive VLBI images as well as snapshot 8.4/8.7-GHz VLBI images of J1918+4937, and show that the appearance of the wiggled jet is consistent with the binary scenario. We develop a jet structural model that handles eccentric orbits. By applying this model to the measured VLBI component positions, we constrain the inclination of the radio jet, as well as the spin angle of the jet emitter SMBH. We find the jet morphological model is consistent with the optical and radio data, and that the secondary SMBH is most likely the jetted one in the system. Furthermore, the decade-long 15-GHz radio flux density monitoring data available for J1918+4937 are compatible with a gradual overall decrease in the the total flux density caused by a slow secular change of the jet inclination due to the spin-orbit precession. J1918+4937 could be an efficient high-energy neutrino source if the horizon of the secondary SMBH is rapidly rotating.


title: Functional linear regression with truncated signatures

abstract: We place ourselves in a functional regression setting and propose a novel methodology for regressing a real output on vector-valued functional covariates. This methodology is based on the notion of signature, which is a representation of a function as an infinite series of its iterated integrals. The signature depends crucially on a truncation parameter for which an estimator is provided, together with theoretical guarantees. An empirical study on both simulated and real-world datasets shows that the resulting methodology is competitive with traditional functional linear models, in particular when the functional covariates take their values in a high dimensional space.


title: Polariton gap and gap-stripe solitons in Zeeman lattices

abstract: We predict that spatially modulated Zeeman splitting resulting in the formation of Zeeman lattice can be used for creation of localized self-sustained excitations in spinor polariton condensates with dominant repulsive interactions. In such lattices, the phenomenon of TE-TM splitting, playing the role of effective spin-orbit interaction, leads to the emergence of the stripe phase and formation of stable gap-stripe solitons with complex intrinsic structure resulting from the presence of two characteristic spatial scales, one of which is set by the period of Zeeman lattice, while other is set by the momentum in the depth of the Brillouin zone, at which such solitons bifurcate from the linear spectrum. Gap-stripe polariton solitons can be excited by suitable resonant pump.


title: Taxonomy and Practical Evaluation of Primality Testing Algorithms

abstract: Modern cryptography algorithms are commonly used to ensure information security. Prime numbers are needed in many asymmetric cryptography algorithms. For example, RSA algorithm selects two large prime numbers and multiplies to each other to obtain a large composite number whose factorization is very difficult. Producing a prime number is not an easy task as they are not distributed regularly through integers. Primality testing algorithms are used to determine whether a particular number is prime or composite. In this paper, an intensive survey is thoroughly conducted among the several primality testing algorithms showing the pros and cons, the time complexity, and a brief summary of each algorithm. Besides, an implementation of these algorithms is accomplished using Java and Python as programming languages to evaluate the efficiency of both the algorithms and the programming languages.


title: Analytic Gradients and Derivative Couplings for Configuration Interaction with All Single Excitations and One Double Excitation -- En Route to Nonadiabatic Dynamics

abstract: We present analytic gradients and derivative couplings for the simplest possible multireference configuration interaction method, CIS-1D, an electronic structure ansatz that includes all single excitations and one lone double excitation on top of a Hartree-Fock reference state. We show that the resulting equations are numerically stable and require the evaluation of a similar number of integrals as compared to standard CIS theory; one can easily differentiate the required frontier orbitals ($\mathtt{h}$ and $\mathtt{l}$) with minimal cost. The resulting algorithm has been implemented within the Q-Chem electronic structure package and should be immediately useful for understanding photochemistry with $S_{0}$-$S_{1}$ crossings.


title: Modeling Joint Lives within Families

abstract: Family history is usually seen as a significant factor insurance companies look at when applying for a life insurance policy. Where it is used, family history of cardiovascular diseases, death by cancer, or family history of high blood pressure and diabetes could result in higher premiums or no coverage at all. In this article, we use massive (historical) data to study dependencies between life length within families. If joint life contracts (between a husband and a wife) have been long studied in actuarial literature, little is known about child and parents dependencies. We illustrate those dependencies using 19th century family trees in France, and quantify implications in annuities computations. For parents and children, we observe a modest but significant positive association between life lengths. It yields different estimates for remaining life expectancy, present values of annuities, or whole life insurance guarantee, given information about the parents (such as the number of parents alive). A similar but weaker pattern is observed when using information on grandparents.


title: Characterization of SARS-CoV-2 Dynamics in the Host

abstract: While many epidemiological models have being proposed to understand and handle COVID-19, too little has been invested to understand how the virus replicates in the human body and potential antiviral can be used to control the replication cycle. In this work, using a control theoretical approach, validated mathematical models of SARS-CoV-2 in humans are properly characterized. A complete analysis of the main dynamic characteristic is developed based on the reproduction number. The equilibrium regions of the system are fully characterized, and the stability of such a regions, formally established. Mathematical analysis highlights critical conditions to decrease monotonically SARS-CoV-2 in the host, such conditions are relevant to tailor future antiviral treatments. Simulation results show the potential benefits of the aforementioned system characterization.


title: Deep unfolding of the weighted MMSE beamforming algorithm

abstract: Downlink beamforming is a key technology for cellular networks. However, computing the transmit beamformer that maximizes the weighted sum rate subject to a power constraint is an NP-hard problem. As a result, iterative algorithms that converge to a local optimum are used in practice. Among them, the weighted minimum mean square error (WMMSE) algorithm has gained popularity, but its computational complexity and consequent latency has motivated the need for lower-complexity approximations at the expense of performance. Motivated by the recent success of deep unfolding in the trade-off between complexity and performance, we propose the novel application of deep unfolding to the WMMSE algorithm for a MISO downlink channel. The main idea consists of mapping a fixed number of iterations of the WMMSE algorithm into trainable neural network layers, whose architecture reflects the structure of the original algorithm. With respect to traditional end-to-end learning, deep unfolding naturally incorporates expert knowledge, with the benefits of immediate and well-grounded architecture selection, fewer trainable parameters, and better explainability. However, the formulation of the WMMSE algorithm, as described in Shi et al., is not amenable to be unfolded due to a matrix inversion, an eigendecomposition, and a bisection search performed at each iteration. Therefore, we present an alternative formulation that circumvents these operations by resorting to projected gradient descent. By means of simulations, we show that, in most of the settings, the unfolded WMMSE outperforms or performs equally to the WMMSE for a fixed number of iterations, with the advantage of a lower computational load.


title: Quantum-enhanced interferometry with large heralded photon-number states

abstract: Quantum phenomena such as entanglement can improve fundamental limits on the sensitivity of a measurement probe. In optical interferometry, a probe consisting of $N$ entangled photons provides up to a $\sqrt{N}$ enhancement in phase sensitivity compared to a classical probe of the same energy. Here, we employ high-gain parametric down-conversion sources and photon-number-resolving detectors to perform interferometry with heralded quantum probes of sizes up to $N=8$ (i.e. measuring up to 16-photon coincidences). Our probes are created by injecting heralded photon-number states into an interferometer, and in principle provide quantum-enhanced phase sensitivity even in the presence of significant optical loss. Our work paves the way towards quantum-enhanced interferometry using large entangled photonic states.


title: Revisiting the minimum length in the Schwinger-Keldysh formalism

abstract: The existence of a minimum length in quantum gravity is investigated by computing the in-in expectation value of the proper distance in the Schwinger-Keldysh formalism. No minimum geometrical length is found for arbitrary gravitational theories to all orders in perturbation theory. Using non-perturbative techniques, we also show that neither the conformal sector of general relativity nor higher-derivative gravity features a minimum length. A minimum length scale, on the other hand, seems to always be present when one considers in-out amplitudes, from which one could extract the energy scale of scattering processes.


title: Symmetry and Isoperimetry for Riemannian Surfaces

abstract: For a domain $\Omega$ in a geodesically convex surface, we introduce a scattering energy $\mathcal{E}(\Omega)$, which measures the asymmetry of $\Omega$ by quantifying its incompatibility with an isometric circle action. We prove several sharp quantitative isoperimetric inequalities involving $\mathcal{E}(\Omega)$ and characterize the domains with vanishing scattering energy by their convexity and rotational symmetry. We also give a new proof of the sharp Sobolev inequality for Riemannian surfaces which is independent of the isoperimetric inequality.


title: Graded identities with involution for the algebra of upper triangular matrices

abstract: Let $F$ be a field of characteristic zero. We prove that if a group grading on $UT_m(F)$ admits a graded involution then this grading is a coarsening of a $\mathbb{Z}^{\lfloor\frac{m}{2}\rfloor}$-grading on $UT_m(F)$ and the graded involution is equivalent to the reflection or symplectic involution on $UT_m(F)$. A finite basis for the $(\mathbb{Z}^{\lfloor\frac{m}{2}\rfloor},\ast)$-identities is exhibited for the reflection and symplectic involutions and the asymptotic growth of the $(\mathbb{Z}^{\lfloor\frac{m}{2}\rfloor},\ast)$-codimensions is determined. As a consequence we prove that for any $G$-grading on $UT_m(F)$ and any graded involution the $(G,\ast)$-exponent is $m$ if $m$ is even and either $m$ or $m+1$ if $m$ is odd. For the algebra $UT_3(F)$ there are, up to equivalence, two non-trivial gradings that admit a graded involution: the canonical $\mathbb{Z}$-grading and the $\mathbb{Z}_2$-grading induced by $(0,1,0)$. We determine a basis for the $(\mathbb{Z}_2,\ast)$-identities and prove that the exponent is $3$. Hence we conclude that the ordinary $\ast$-exponent for $UT_3(F)$ is $3$.


title: Bayesian Neural Network via Stochastic Gradient Descent

abstract: The goal of bayesian approach used in variational inference is to minimize the KL divergence between variational distribution and unknown posterior distribution. This is done by maximizing the Evidence Lower Bound (ELBO). A neural network is used to parametrize these distributions using Stochastic Gradient Descent. This work extends the work done by others by deriving the variational inference models. We show how SGD can be applied on bayesian neural networks by gradient estimation techniques. For validation, we have tested our model on 5 UCI datasets and the metrics chosen for evaluation are Root Mean Square Error (RMSE) error and negative log likelihood. Our work considerably beats the previous state of the art approaches for regression using bayesian neural networks.


title: Pseudo-Sylvester domains and skew Laurent polynomials over firs

abstract: Building on recent work of Jaikin-Zapirain, we provide a homological criterion for a ring to be a pseudo-Sylvester domain, that is, to admit a division ring of fractions over which all stably full matrices become invertible. We use the criterion to study skew Laurent polynomial rings over free ideal rings (firs). As an application of our methods, we prove that crossed products of division rings with free-by-{infinite cyclic} and surface groups are pseudo-Sylvester domains unconditionally and Sylvester domains if and only if they admit stably free cancellation. This relies on the recent proof of the Farrell--Jones conjecture for normally poly-free groups and extends previous results of Linnell--L\"uck and Jaikin-Zapirain on universal localizations and universal fields of fractions of such crossed products.


title: Experimental demonstration of phase estimation advantage in presence of depolarizing noise by using coherent measurements

abstract: We report an experimental investigation of the role of measurement in quantum metrology when the states of the probes are mixed. In particular, we investigated optimized local measurements and general global projective measurements, involving entangling operations, on noisy Werner states of polarization entangled photons. We demonstrate experimentally that global measurement presents an advantage in parameter estimation with respect to the optimized local strategy. Moreover, the global strategy provides unambiguous information about the parameter of interest even when the amount of noise is not well characterized. This shows that the coherence in quantum operations, such as the Bell-state projection device used in our protocol, can be used to further boost the quantum advantage in metrology and play a fundamental role in the design of future quantum measurement devices.


title: A Machine Learning-Based Migration Strategy for Virtual Network Function Instances

abstract: With the growing demand for data connectivity, network service providers are faced with the task of reducing their capital and operational expenses while simultaneously improving network performance and addressing the increased demand. Although Network Function Virtualization (NFV) has been identified as a promising solution, several challenges must be addressed to ensure its feasibility. In this paper, we address the Virtual Network Function (VNF) migration problem by developing the VNF Neural Network for Instance Migration (VNNIM), a migration strategy for VNF instances. The performance of VNNIM is further improved through the optimization of the learning rate hyperparameter through particle swarm optimization. Results show that the VNNIM is very effective in predicting the post-migration server exhibiting a binary accuracy of 99.07% and a delay difference distribution that is centered around a mean of zero when compared to the optimization model. The greatest advantage of VNNIM, however, is its run-time efficiency highlighted through a run-time analysis.


title: Interaction Networks: Using a Reinforcement Learner to train other Machine Learning algorithms

abstract: The wiring of neurons in the brain is more flexible than the wiring of connections in contemporary artificial neural networks. It is possible that this extra flexibility is important for efficient problem solving and learning. This paper introduces the Interaction Network. Interaction Networks aim to capture some of this extra flexibility. An Interaction Network consists of a collection of conventional neural networks, a set of memory locations, and a DQN or other reinforcement learner. The DQN decides when each of the neural networks is executed, and on what memory locations. In this way, the individual neural networks can be trained on different data, for different tasks. At the same time, the results of the individual networks influence the decision process of the reinforcement learner. This results in a feedback loop that allows the DQN to perform actions that improve its own decision-making. Any existing type of neural network can be reproduced in an Interaction Network in its entirety, with only a constant computational overhead. Interaction Networks can then introduce additional features to improve performance further. These make the algorithm more flexible and general, but at the expense of being harder to train. In this paper, thought experiments are used to explore how the additional abilities of Interaction Networks could be used to improve various existing types of neural networks. Several experiments have been run to prove that the concept is sound. These show that the basic idea works, but they also reveal a number of challenges that do not appear in conventional neural networks, which make Interaction Networks very hard to train. Further research needs to be done to alleviate these issues. A number of promising avenues of research to achieve this are outlined in this paper.


title: Evolution of Group-Theoretic Cryptology Attacks using Hyper-heuristics

abstract: In previous work, we developed a single Evolutionary Algorithm (EA) to solve random instances of the Anshel-Anshel-Goldfeld (AAG) key exchange protocol over polycyclic groups. The EA consisted of six simple heuristics which manipulated strings. The present work extends this by exploring the use of hyper-heuristics in group-theoretic cryptology for the first time. Hyper-heuristics are a way to generate new algorithms from existing algorithm components (in this case the simple heuristics), with the EAs being one example of the type of algorithm which can be generated by our hyper-heuristic framework. We take as a starting point the above EA and allow hyper-heuristics to build on it by making small tweaks to it. This adaptation is through a process of taking the EA and injecting chains of heuristics built from the simple heuristics. We demonstrate we can create novel heuristic chains, which when placed in the EA create algorithms which out-perform the existing EA. The new algorithms solve a markedly greater number of random AAG instances than the EA for harder instances. This suggests the approach could be applied to many of the same kinds of problems, providing a framework for the solution of cryptology problems over groups. The contribution of this paper is thus a framework to automatically build algorithms to attack cryptology problems.


title: Modified Bohmian quantum potential due to the second quantization of Schrodinger equation

abstract: A causal interpretation of the quantum world needs second quantization in order to cover phenomena like creation and annihilation of particles, and this leads to quantum field theory. The causal effects of the second quantization can be described through a new quantum potential. In this article we have dealt with the second quantization of Schrodinger and its effects on the path of a particle. This generalization leads to a modified Schrodinger which affects the particle through a modified quantum potential and a new term in the continuity equation. We have shown that these effects can provide a framework for the explanation of the creation and annihilation phenomena and other effects of quantum field theory on the development of a particle.


title: p-Diamond, Si, GaN and InGaAs TeraFETs

abstract: p-diamond field effect transistors (FETs) featuring large effective mass, long momentum relaxation time and high carrier mobility are a superb candidate for plasmonic terahertz (THz) applications. Previous studies have shown that p-diamond plasmonic THz FETs (TeraFETs) could operate in plasmonic resonant mode at a low frequency window of 200 GHz to ~600 GHz, thus showing promising potential for beyond 5G sub-THz applications. In this work, we explore the advantages of p-diamond transistors over n-diamond, Si, GaN and InGaAs TeraFETs and estimate the minimum mobility required for the resonant plasmons. Our numerical simulation shows that the p-diamond TeraFET has a relatively low minimum resonant mobility, and thus could enable resonant detection. The diamond response characteristics can be adjusted by changing operating temperature. A decrease of temperature from 300 K to 77 K improves the detection performance of TeraFETs. At both room temperature and 77 K, the p-diamond TeraFET presents a high detection sensitivity in a large dynamic range. When the channel length is reduced to 20 nm, the p-diamond TeraFET exhibits the highest DC response among all types of TeraFETs in a large frequency window.


title: Visualizing Quantum Phases And Identifying Quantum Phase Transitions By Nonlinear Dimensionality Reduction

abstract: Identifying quantum phases and phase transitions is key to understand complex phenomena in statistical physics. In this work, we propose an unconventional strategy to access quantum phases and phase transitions by visualization based on the distribution of ground states in Hilbert space. By mapping the quantum states in Hilbert space onto a two-dimensional feature space using an unsupervised machine learning method, distinct phases can be directly specified and quantum phase transitions can be well identified. Our proposal is benchmarked on gapped, critical, and topological phases in several strongly correlated spin systems. As this proposal directly learns quantum phases and phase transitions from the distributions of the quantum states, it does not require priori knowledge of order parameters of physical systems, which thus indicates a perceptual route to identify quantum phases and phase transitions particularly in complex systems by visualization through learning.


title: A sparse equidistribution result for $(\mathrm{SL}(2,\mathbb{R})/\Gamma_0)^n$

abstract: Let $G=\mathrm{SL}(2,\mathbb{R})^n$, let $\Gamma=\Gamma_0^n$, where $\Gamma_0$ is a co-compact lattice in $\mathrm{SL}(2,\mathbb{R})$, let $F(\mathbf{x})$ be a non-singular quadratic form and let $u(x_1,...,x_n)$ denote the unipotent elements in $G$ which generate the standard $n$ dimensional horospherical subgroup, consisting of $2\times 2$ upper triangular unipotent matrices in each co-ordinate. We prove that in absence of any local obstructions for $F$, given any $x_0\in G/\Gamma$, the sparse subset $\{u(\mathbf{x})x_0:\in\mathbb{Z}^n, F(\mathbf{x})=0\}$ equidistributes in $G/\Gamma$ as long as $n\geq 481$, independent of the spectral gap of $\Gamma_0$.


title: Structure and overstability of resistive modes with runaway electrons

abstract: We investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge whose properties depend on the ratio of the Alfv\'en velocity to the runaway electron convection speed. Due to the contribution from this layer, both the tearing mode and kink mode will have a real frequency in addition to a growth rate. The derived dispersion relation has been compared with numerical results using both a simplified eigenvalue calculation and a M3D-C1 linear simulation, and good agreement is found in both cases.