title: No-cloning implies unalterability of the past

abstract: A common way of stating the non-cloning theorem -- one of distinguishing characteristics of quantum theory -- is that one cannot make a copy of an arbitrary unknown quantum state. Even though this theorem is an important part of the ongoing discussion of the nature of a quantum state, the role of the theorem in the logical-algebraic approach to quantum theory has not yet been systematically studied. According to the standard point of view (which is in line with the logical tradition), quantum cloning amounts to two classical rules of inference, namely, monotonicity and idempotency of entailment. One can conclude then that the whole of quantum theory should be described through a logic wherein these rules do not hold, which is linear logic. However, in accordance with a supervaluational semantics (that allows one to retain all the theorems of classical logic while admitting `truth-value gaps'), quantum cloning necessitates the permanent loss of the truth values of experimental quantum propositions which violates the unalterability of the past. The present paper demonstrates this.


title: Extension of the HCOOH and CO2 solid-state reaction network during the CO freeze-out stage: inclusion of H2CO

abstract: Formic acid (HCOOH) and carbon dioxide (CO2) are simple species that have been detected in the interstellar medium. The solid-state formation pathways of these species under experimental conditions relevant to prestellar cores are primarily based off of weak infrared transitions of the HOCO complex and usually pertain to the H2O-rich ice phase, and therefore more experimental data are desired. In this article, we present a new and additional solid-state reaction pathway that can form HCOOH and CO2 ice at 10 K 'non-energetically' in the laboratory under conditions related to the "heavy" CO freeze-out stage in dense interstellar clouds, i.e., by the hydrogenation of an H2CO:O2 ice mixture. This pathway is used to piece together the HCOOH and CO2 formation routes when H2CO or CO reacts with H and OH radicals. Temperature programmed desorption - quadrupole mass spectrometry (TPD-QMS) is used to confirm the formation and pathways of newly synthesized ice species as well as to provide information on relative molecular abundances. Reflection absorption infrared spectroscopy (RAIRS) is additionally employed to characterize reaction products and determine relative molecular abundances. We find that for the conditions investigated in conjunction with theoretical results from the literature, H+HOCO and HCO+OH lead to the formation of HCOOH ice in our experiments. Which reaction is more dominant can be determined if the H+HOCO branching ratio is more constrained by computational simulations, as the HCOOH:CO2 abundance ratio is experimentally measured to be around 1.8:1. H+HOCO is more likely than OH+CO (without HOCO formation) to form CO2. Isotope experiments presented here further validate that H+HOCO is the dominant route for HCOOH ice formation in a CO-rich CO:O2 ice mixture that is hydrogenated. These data will help in the search and positive identification of HCOOH ice in prestellar cores.


title: A New Hierarchical Software Architecture Towards Safety-Critical Aspects of a Drone System

abstract: In this paper, a new hierarchical software architecture is proposed to improve the safety and reliability of a safety-critical drone system from the perspective of its source code. The proposed architecture uses formal verification methods to ensure that the implementation of each module satisfies its expected design specification, so that it prevents a drone from crashing due to unexpected software failures. This study builds on top of a formally verified operating system kernel, certified kit operating system (CertiKOS). Since device drivers are considered the most important parts affecting the safety of the drone system, we focus mainly on verifying bus drivers such as the serial peripheral interface and the inter-integrated circuit drivers in a drone system using a rigorous formal verification method. Experiments have been carried out to demonstrate the improvement in reliability in case of device anomalies.


title: Application of the amended Coriolis flowmeter "bubble theory" to sound propagation and attenuation in aerosols and hydrosols

abstract: The existing viscous and incompressible theory of isothermal sound propagation and attenuation in suspensions considers solid particles which are infinitely viscous. We extend the theory by applying the amended Coriolis flowmeter "bubble theory". Here, the drag force is a function of both the fluid and particle Stokes numbers and the particle-to-fluid ratio of the dynamic viscosity [V.Galindo and G.Gerbeth, A note on the force of an accelerating spherical drop at low-Reynolds number, Phys. Fluids A Vol. 5, 3290-3292 (1993)]. Aerosol and hydrosol examples are presented and differences between the original and extended theories are discussed.


title: Osculating Versus Intersecting Circles in Space-Based Microlens Parallax Degeneracies

abstract: I investigate the origin of arc degeneracies in satellite microlens parallax pi_E measurements with only late time data, e.g., t > t0 + t_E as seen from the satellite. I show that these are due to partial overlap of a series of osculating, exactly circular, degeneracies in the pi_E plane, each from a single measurement. In events with somewhat earlier data, these long arcs break up into two arclets, or (with even earlier data) two points, because these earlier measurements give rise to intersecting rather than osculating circles. The two arclets (or points) then constitute one pair of degeneracies in the well-known four-fold degeneracy of space-based microlens parallax. Using this framework of intersecting circles, I show that next-generation microlens satellite experiments could yield good pi_E determinations with only about five measurements per event, i.e., about 30 observations per day to monitor 1500 events per year. This could plausibly be done with a small (hence cheap, in the spirit of Gould & Yee 2012) satellite telescope, e.g., 20 cm.


title: Sherman's inequality and its converse for strongly convex functions with applications to generalized f-divergences

abstract: Considering the weighted concept of majorization, Sherman obtained generalization of majorization inequality for convex functions known as Sherman's inequality. We extend Sherman's result to the class of n-strongly convex functions using extended idea of convexity to the class of strongly convex functions. We also obtaine upper bound for Sherman's inequality, so called the converse Sherman inequality, and as easy consequences we get Jensen's as well as majorization inequality and their conversions for strongly convex functions. Obtained results are stronger versions for analogous results for convex functions. As applications, we introduced a generalized concept of f-divergence and derived some reverse relations for such concept.


title: Element-specific visualization of dynamic magnetic coupling in a Co/Py bilayer microstructure

abstract: We present the element-specific and time resolved visualization of uniform ferromagnetic resonance excitations of a Permalloy (Py) disk - Cobalt (Co) stripe bilayer microstructure. The transverse high frequency component of the resonantly excited magnetization is sampled in the ps regime by a combination of ferromagnetic resonance (FMR) and scanning transmission X-ray microscopy (STXM- FMR) recording snapshots of the local magnetization precession of Py and Co with nanometer spatial resolution. The approach allows us to individually image the resonant dynamic response of each element, and we find that angular momentum is transferred from the Py disk to the Co stripe and vice versa at their respective resonances. The integral (cavity) FMR spectrum of our sample shows an unexpected additional third resonance. This resonance is observed in the STXM-FMR experiments as well and our microscopic findings suggest that it is governed by magnetic exchange between Py and Co, showing for the Co stripe a difference in relative phase of the magnetization due to stray field influence.


title: Input Modeling and Uncertainty Quantification for Improving Volatile Residential Load Forecasting

abstract: Load forecasting has long been recognized as an important building block for all utility operational planning efforts. Over the recent years, it has become ever more challenging to make accurate forecasts due to the proliferation of distributed energy resources, despite the abundance of existing load forecasting methods. In this paper, we identify one drawback suffered by most load forecasting methods: neglect to thoroughly address the impact of input errors on load forecasts. As a potential solution, we propose to incorporate input modeling and uncertainty quantification to improve load forecasting performance via a two-stage approach. The proposed two-stage approach has the following merits. (1) It provides input modeling and quantifies the impact of input errors, rather than neglecting or mitigating the impact, a prevalent practice of existing methods. (2) It propagates the impact of input errors into the ultimate point and interval predictions for the target customer's load to improve predictive performance. (3) A variance-based global sensitivity analysis method is further proposed for input-space dimensionality reduction in both stages to enhance the computational efficiency. Numerical experiments show that the proposed two-stage approach outperforms competing load forecasting methods in terms of both point predictive accuracy and coverage ability of the predictive intervals.


title: Richly Activated Graph Convolutional Network for Action Recognition with Incomplete Skeletons

abstract: Current methods for skeleton-based human action recognition usually work with completely observed skeletons. However, in real scenarios, it is prone to capture incomplete and noisy skeletons, which will deteriorate the performance of traditional models. To enhance the robustness of action recognition models to incomplete skeletons, we propose a multi-stream graph convolutional network (GCN) for exploring sufficient discriminative features distributed over all skeleton joints. Here, each stream of the network is only responsible for learning features from currently unactivated joints, which are distinguished by the class activation maps (CAM) obtained by preceding streams, so that the activated joints of the proposed method are obviously more than traditional methods. Thus, the proposed method is termed richly activated GCN (RA-GCN), where the richly discovered features will improve the robustness of the model. Compared to the state-of-the-art methods, the RA-GCN achieves comparable performance on the NTU RGB+D dataset. Moreover, on a synthetic occlusion dataset, the performance deterioration can be alleviated by the RA-GCN significantly.


title: The onset of $\Lambda\Lambda$ hypernuclear binding

abstract: Binding energies of light, $A\leq 6$, $\Lambda\Lambda$ hypernuclei are calculated using the stochastic variational method in a pionless effective field theory (EFT) approach at leading order with the purpose of assessing critically the onset of binding in the strangeness S=-2 hadronic sector. The EFT input in this sector consists of (i) a $\Lambda\Lambda$ contact term constrained by the $\Lambda\Lambda$ scattering length $a_{\Lambda\Lambda}$, using a range of values compatible with $\Lambda\Lambda$ correlations observed in relativistic heavy ion collisions, and (ii) a $\Lambda\Lambda N$ contact term constrained by the only available $A\leq 6$ $\Lambda\Lambda$ hypernucler binding energy datum of $^{6}_{\Lambda\Lambda}$He. The recently debated neutral three-body and four-body systems $^{3}_{\Lambda\Lambda}$n and $^{4}_{\Lambda\Lambda}$n are found unbound by a wide margin. A relatively large value of $|a_{\Lambda\Lambda}| \gtrsim 1.5$ fm is needed to bind $^{4}_{\Lambda\Lambda}$H, thereby questioning its particle stability. In contrast, the particle stability of the $A=5$ $\Lambda\Lambda$ hypernuclear isodoublet $^{5}_{\Lambda\Lambda} $H--$^{5}_{\Lambda\Lambda}$He is robust, with $\Lambda$ separation energy of order 1 MeV.


title: Dimensional crossover and topological nature of the thin films of a three-dimensional topological insulator by band gap engineering

abstract: Identification and control of topological phases in topological thin films offer great opportunity for fundamental research and the fabrication of topology-based devices. Here, combining molecular beam epitaxy, angle-resolved photoemission spectroscopy and ab-initio calculations, we investigate the electronic structure evolution in (Bi1-xInx)2Se3 films with thickness from 2 to 13 quintuple layers. We identify several phases with their characteristic topological nature and evolution between them, i.e., dimensional crossover from a three-dimensional topological insulator with gapless surface state to its two-dimensional counterpart with gapped surface state, and topological phase transition from topological insulator to a normal semiconductor with increasing In concentration x. Furthermore, by introducing In alloying as an external knob of band gap engineering, we experimentally demonstrated the trivial topological nature of Bi2Se3 thin films (below 6 quintuple layers) as two-dimensional gapped systems, in consistent with our theoretical calculations. Our results provide not only a comprehensive phase diagram of (Bi1-xInx)2Se3 and a route to control its phase evolution, but also a practical way to experimentally determine the topological properties of a gapped compound by topological phase transition and band gap engineering.


title: Towards Comparing Programming Paradigms

abstract: Rapid technological progress in computer sciences finds solutions and at the same time creates ever more complex requirements. Due to an evolving complexity todays programming languages provide powerful frameworks which offer standard solutions for recurring tasks to assist the programmer and to avoid the re-invention of the wheel with so-called out-of-the-box-features. In this paper, we propose a way of comparing different programming paradigms on a theoretical, technical and practical level. Furthermore, the paper presents the results of an initial comparison of two representative programming approaches, both in the closed SAP environment.


title: Global attractors and their upper semicontinuity for a structural damped wave equation with supercritical nonlinearity on $\mathbb{R}^{N}$

abstract: The paper investigates the existence of global attractors and their upper semicontinuity for a structural damped wave equation on $\mathbb{R}^{N}: u_{tt}-\Delta u+(-\Delta)^\alpha u_{t}+u_{t}+u+g(u)=f(x)$, where $\alpha\in (1/2, 1)$ is called a dissipative index. We propose a new method based on the harmonic analysis technique and the commutator estimate to exploit the dissipative effect of the structural damping $(-\Delta)^\alpha u_{t}$ and to overcome the essential difficulty: "both the unbounded domain $\mathbb{R}^N$ and the supercritical nonlinearity cause that the Sobolev embedding loses its compactness"; Meanwhile we show that there exists a supercritical index $p_\alpha\equiv\frac{N+4\alpha}{N-4\alpha}$ depending on $\alpha$ such that when the growth exponent $p$ of the nonlinearity $g(u)$ is up to the supercritical range: $1\leqslant p<p_\alpha$: (i) the IVP of the equation is well-posed and its solution is of additionally global smoothness when $t>0$; (ii) the related solution semigroup possesses a global attractor $\mathcal{A}_\alpha$ in natural energy space for each $\alpha\in (1/2, 1)$; (iii) the family of global attractors $\{\mathcal{A}_\alpha\}_{\alpha\in (1/2, 1) }$ is upper semicontinuous at each point $\alpha_0\in (1/2, 1)$.


title: On the harmonic extension approach to fractional powers in Banach spaces

abstract: We show that fractional powers of general sectorial operators on Banach spaces can be obtained by the harmonic extension approach. Moreover, for the corresponding second order ordinary differential equation with incomplete data describing the harmonic extension we prove existence and uniqueness of a bounded solution (i.e. of the harmonic extension).


title: Functional Form of the Superconducting Critical Temperature from Machine Learning

abstract: Predicting the critical temperature $T_c$ of new superconductors is a notoriously difficult task, even for electron-phonon paired superconductors for which the theory is relatively well understood. Early attempts by McMillan and Allen and Dynes to improve on the weak-coupling BCS formula led to closed-form approximate relations between $T_c$ and various measures of the phonon spectrum and the electron-phonon interaction appearing in Eliashberg theory. Here we propose that these approaches can be improved with the use of machine learning algorithms. As an initial test, we train a model for identifying low-dimensional descriptors using the $T_c < 10$ K data tested by Allen and Dynes, and show that a simple analytical expression thus obtained improves upon the Allen-Dynes fit. Furthermore, the prediction for the recently discovered high $T_c$ material H$_3$S at high pressure is quite reasonable. Interestingly, $T_c$'s for more recently discovered superconducting systems with a more two-dimensional electron-phonon coupling, which do not follow Allen and Dynes' expression, also do not follow our analytic expression. Thus, this machine learning approach appears to be a powerful method for highlighting the need for a new descriptor beyond those used by Allen and Dynes to describe their set of isotropic electron-phonon coupled superconductors. We argue that this machine learning method, and its implied need for a descriptor characterizing Fermi surface properties, represents a promising new approach to superconductor materials discovery which may eventually replace the serendipitous discovery paradigm begun by Kamerlingh Onnes.


title: Improved Beckner-Sobolev inequalities on K\"ahler manifolds

abstract: We prove new Beckner-Sobolev type inequalities on compact K\"{a}hler manifolds with positive Ricci curvature. As an application, we obtain a diameter upper bound that improves the Bonnet-Myers bound.


title: Identifying collaborators in large codebases

abstract: The way developers collaborate inside and particularly across teams often escapes management's attention, despite a formal organization with designated teams being defined. Observability of the actual, organically formed engineering structure provides decision makers invaluable additional tools to manage their talent pool. To identify existing inter and intra-team interactions - and suggest relevant opportunities for suitable collaborations - this paper studies contributors' commit activity, usage of programming languages, and code identifier topics by embedding and clustering them. We evaluate our findings collaborating with the GitLab organization, analyzing 117 of their open source projects. We show that we are able to restore their engineering organization in broad strokes, and also reveal hidden coding collaborations as well as justify in-house technical decisions.


title: Time-Energy Tradeoffs for Evacuation by Two Robots in the Wireless Model

abstract: Two robots stand at the origin of the infinite line and are tasked with searching collaboratively for an exit at an unknown location on the line. They can travel at maximum speed $b$ and can change speed or direction at any time. The two robots can communicate with each other at any distance and at any time. The task is completed when the last robot arrives at the exit and evacuates. We study time-energy tradeoffs for the above evacuation problem. The evacuation time is the time it takes the last robot to reach the exit. The energy it takes for a robot to travel a distance $x$ at speed $s$ is measured as $xs^2$. The total and makespan evacuation energies are respectively the sum and maximum of the energy consumption of the two robots while executing the evacuation algorithm. Assuming that the maximum speed is $b$, and the evacuation time is at most $cd$, where $d$ is the distance of the exit from the origin, we study the problem of minimizing the total energy consumption of the robots. We prove that the problem is solvable only for $bc \geq 3$. For the case $bc=3$, we give an optimal algorithm, and give upper bounds on the energy for the case $bc>3$. We also consider the problem of minimizing the evacuation time when the available energy is bounded by $\Delta$. Surprisingly, when $\Delta$ is a constant, independent of the distance $d$ of the exit from the origin, we prove that evacuation is possible in time $O(d^{3/2}\log d)$, and this is optimal up to a logarithmic factor. When $\Delta$ is linear in $d$, we give upper bounds on the evacuation time.


title: Harvesting Information from Captions for Weakly Supervised Semantic Segmentation

abstract: Since acquiring pixel-wise annotations for training convolutional neural networks for semantic image segmentation is time-consuming, weakly supervised approaches that only require class tags have been proposed. In this work, we propose another form of supervision, namely image captions as they can be found on the Internet. These captions have two advantages. They do not require additional curation as it is the case for the clean class tags used by current weakly supervised approaches and they provide textual context for the classes present in an image. To leverage such textual context, we deploy a multi-modal network that learns a joint embedding of the visual representation of the image and the textual representation of the caption. The network estimates text activation maps (TAMs) for class names as well as compound concepts, i.e. combinations of nouns and their attributes. The TAMs of compound concepts describing classes of interest substantially improve the quality of the estimated class activation maps which are then used to train a network for semantic segmentation. We evaluate our method on the COCO dataset where it achieves state of the art results for weakly supervised image segmentation.


title: Four Lectures on Closed String Field Theory

abstract: The following notes derive from review lectures on closed string field theory given at the Galileo Galilei Institute for Theoretical Physics in March 2019.


title: Boundary control of partial differential equations using frequency domain optimization techniques

abstract: We present a frequency domain based $H_\infty$-control strategy to solve boundary control problems for systems governed by parabolic or hyperbolic partial differential equation, where controllers are constrained to be physically implementable and of simple structure suited for practical applications. The efficiency of our technique is demonstrated by controlling a reaction-diffusion equation with input delay, and a wave equation with boundary anti-damping.


title: Generical behavior of flows strongly monotone with respect to high-rank cones

abstract: We consider a $C^{1,\alpha}$ smooth flow in $\mathbb{R}^n$ which is "strongly monotone" with respect to a cone $C$ of rank $k$, a closed set that contains a linear subspace of dimension $k$ and no linear subspaces of higher dimension. We prove that orbits with initial data from an open and dense subset of the phase space are either pseudo-ordered or convergent to equilibria. This covers the celebrated Hirsch's Generic Convergence Theorem in the case $k=1$, yields a generic Poincar\'{e}-Bendixson Theorem for the case $k=2$, and holds true with arbitrary dimension $k$. Our approach involves the ergodic argument using the $k$-exponential separation and the associated $k$-Lyapunov exponent (that reduces to the first Lyapunov exponent if $k=1$).


title: Controlling secondary flow in Taylor-Couette turbulence through spanwise-varying roughness

abstract: Highly turbulent Taylor-Couette flow with spanwise-varying roughness is investigated experimentally and numerically (direct numerical simulations (DNS) with an immersed boundary method (IBM)) to determine the effects of the spacing and axial width $s$ of the spanwise varying roughness on the total drag and {on} the flow structures. We apply sandgrain roughness, in the form of alternating {rough and smooth} bands to the inner cylinder. Numerically, the Taylor number is $\mathcal{O}(10^9)$ and the roughness width is varied between $0.47\leq \tilde{s}=s/d \leq 1.23$, where $d$ is the gap width. Experimentally, we explore $\text{Ta}=\mathcal{O}(10^{12})$ and $0.61\leq \tilde s \leq 3.74$. For both approaches the radius ratio is fixed at $\eta=r_i/r_o = 0.716$, with $r_i$ and $r_o$ the radius of the inner and outer cylinder respectively. We present how the global transport properties and the local flow structures depend on the boundary conditions set by the roughness spacing $\tilde{s}$. Both numerically and experimentally, we find a maximum in the angular momentum transport as function of $\tilde s$. This can be atributed to the re-arrangement of the large-scale structures triggered by the presence of the rough stripes, leading to correspondingly large-scale turbulent vortices.


title: Phase behavior of flexible and semiflexible polymers in solvents of varying quality

abstract: The interplay of nematic order and phase separation in solutions of semiflexible polymers in solvents of variable quality is investigated by density functional theory (DFT) and molecular dynamics (MD) simulations. We studied coarse-grained models, with a bond-angle potential to control chain stiffness, for chain lengths comparable to the persistence length of the chains. We varied both the density of the monomeric units and the effective temperature that controls the quality of the implicit solvent. For very stiff chains only a single transition from an isotropic fluid to a nematic is found, with a phase diagram of "swan-neck" topology. For less stiff chains, however, also unmixing between isotropic fluids of different concentration, ending in a critical point, occurs for temperatures above a triple point. The associated critical behavior is examined in the MD simulations and found compatible with Ising universality. Apart from this critical behavior, DFT calculations agree qualitatively with the MD simulations.


title: Cosmological singularity as an informational seed for Everything

abstract: It is shown how to place some amount of matter into the cosmological singularity and to encode its state. A free and massless scalar field is considered as a prototype of matter. Two different but coherent approaches to this issue are presented. The expression for the scalar particles' spectral energy density, which is initially encoded at the singularity, is deduced. An informational aspect of the problem is discussed.


title: Almost Unsupervised Text to Speech and Automatic Speech Recognition

abstract: Text to speech (TTS) and automatic speech recognition (ASR) are two dual tasks in speech processing and both achieve impressive performance thanks to the recent advance in deep learning and large amount of aligned speech and text data. However, the lack of aligned data poses a major practical problem for TTS and ASR on low-resource languages. In this paper, by leveraging the dual nature of the two tasks, we propose an almost unsupervised learning method that only leverages few hundreds of paired data and extra unpaired data for TTS and ASR. Our method consists of the following components: (1) a denoising auto-encoder, which reconstructs speech and text sequences respectively to develop the capability of language modeling both in speech and text domain; (2) dual transformation, where the TTS model transforms the text $y$ into speech $\hat{x}$, and the ASR model leverages the transformed pair $(\hat{x},y)$ for training, and vice versa, to boost the accuracy of the two tasks; (3) bidirectional sequence modeling, which addresses error propagation especially in the long speech and text sequence when training with few paired data; (4) a unified model structure, which combines all the above components for TTS and ASR based on Transformer model. Our method achieves 99.84% in terms of word level intelligible rate and 2.68 MOS for TTS, and 11.7% PER for ASR on LJSpeech dataset, by leveraging only 200 paired speech and text data (about 20 minutes audio), together with extra unpaired speech and text data.


title: Serrodyne frequency translation using time-modulated metasurfaces

abstract: Temporally modulated metamaterials have attracted significant attention recently due to their non-reciprocal and frequency converting properties. Here, a transparent, time-modulated metasurface that functions as a serrodyne frequency translator, is reported at X-band frequencies. With a simple biasing architecture, the metasurface provides electrically-tunable transmission phase that covers $360^\circ$. A sawtooth waveform is used to modulate the metasurface, allowing Doppler-like frequency translation. Two such metasurfaces can be cascaded together to achieve magnetless devices that perform either phase or amplitude non-reciprocity.


title: Distributions with Decay and Restriction Problems

abstract: In this paper we introduce a new type of restriction problem, called the \textit{restriction problem with moments}. We show that the surface area measure of the sphere satisfies the $L^p$-$L^2$ restriction problem with moments if $1 \leq p < \frac{2(d+2)}{d+3}$ and that the Frostman measure constructed by Salem satisfies the $L^p$-$L^2$ restriction problem with moments if $1 \leq p < \frac{2(2-2\alpha+\beta)}{4(1-\alpha)+\beta}$ for certain values of $\alpha$ and $\beta$. The main tool to obtain these new type of restriction phenomenon is the notion of distributions with decay in connection with classes of global $L^q$ ultradifferentiable functions. We develop the notion of distributions with decay and use it to define global wavefront sets of classes of function spaces, including $L^p$-Sobolev spaces on \mathbb{R}^d$ as well as global $L^q$-Denjoy Carleman functions. We also introduce the corresponding notion of microglobal regularity. We prove a characterization of distributions (in a given function space) with decay in terms of microglobal regularity in every direction of their Fourier transforms.


title: Ultra-compact photodetection in atomically thin MoSe$_2$

abstract: Excitons in atomically-thin semiconductors interact very strongly with electromagnetic radiation and are necessarily close to a surface. Here, we exploit the deep-subwavelength confinement of surface plasmon polaritons (SPPs) at the edge of a metal-insulator-metal plasmonic waveguide and their proximity of 2D excitons in an adjacent atomically thin semiconductor to build an ultra-compact photodetector. When subject to far-field excitation we show that excitons are created throughout the dielectric gap region of our waveguide and converted to free carriers primarily at the anode of our device. In the near-field regime, strongly confined SPPs are launched, routed and detected in a 20nm narrow region at the interface between the waveguide and the monolayer semiconductor. This leads to an ultra-compact active detector region of only ~0.03$\mu m ^2$ that absorbs 86% of the propagating energy in the SPP. Due to the electromagnetic character of the SPPs, the spectral response is essentially identical to the far-field regime, exhibiting strong resonances close to the exciton energies. While most of our experiments are performed on monolayer thick MoSe$_2$, the photocurrent-per-layer increases super linearly in multilayer devices due to the suppression of radiative exciton recombination. These results demonstrate an integrated device for nanoscale routing and detection of light with the potential for on-chip integration at technologically relevant, few-nanometer length scales.


title: Dark sector evolution in Horndeski models

abstract: We use the Equation of State (EoS) approach to study the evolution of the dark sector in Horndeski models, the most general scalar-tensor theories with second order equations of motion. By including the effects of the dark sector into our code EoS\_class, we demonstrate the numerical stability of the formalism and excellent agreement with results from other publicly available codes for a range of parameters describing the evolution of the function characterising the perturbations for Horndeski models, $\alpha_{\rm x}$, with ${\rm x}=\{{\rm K}, {\rm B}, {\rm M}, {\rm T}\}$. After demonstrating that on sub-horizon scales ($k\gtrsim 10^{-3}~{\rm Mpc}^{-1}$ at $z=0$) velocity perturbations in both the matter and the dark sector are typically subdominant with respect to density perturbations in the equation of state for perturbations, we find an attractor solution for the dark sector gauge-invariant density perturbation $\Delta_{\rm ds}$ by neglecting its time derivatives in the equation describing its time evolution, as commonly done in the well-known quasi-static approximation. Using this result, we provide simplified expressions for the equation-of-state functions: the dark sector entropy perturbations $w_{\rm ds}\Gamma_{\rm ds}$ and anisotropic stress $w_{\rm ds}\Pi_{\rm ds}$. From this we derive a growth factor-like equation for both matter and dark sector and are able to capture the relevant physics for several observables with great accuracy. We finally present new analytical expressions for the well-known modified gravity phenomenological functions $\mu$, $\eta$ and $\Sigma$ for a generic Horndeski model as functions of $\alpha_{\rm x}$. We show that on small scales they reproduce expressions presented in previous works, but on large scales, we find differences with respect to other works.


title: Non-equilibrium probability flux of a thermally driven micromachine

abstract: We discuss the non-equilibrium statistical mechanics of a thermally driven micromachine consisting of three spheres and two harmonic springs [Y. Hosaka et al., J. Phys. Soc. Jpn. 86, 113801 (2017)]. We obtain the non-equilibrium steady state probability distribution function of such a micromachine and calculate its probability flux in the corresponding configuration space. The resulting probability flux can be expressed in terms of a frequency matrix that is used to distinguish between a non-equilibrium steady state and a thermal equilibrium state satisfying detailed balance. The frequency matrix is shown to be proportional to the temperature difference between the spheres. We obtain a linear relation between the eigenvalue of the frequency matrix and the average velocity of a thermally driven micromachine that can undergo a directed motion in a viscous fluid. This relation is consistent with the scallop theorem for a deterministic three-sphere microswimmer.


title: Cutting plane oracles for non-smooth trust-regions

abstract: We prove global convergence of a bundle trust region algorithm for non-smooth non-convex optimization, where cutting planes are generated by oracles respecting four basic rules. The benefit is that convergence theory applies to a large variety of methods encountered in practice. This includes in particular the method of downshifted tangents, for which previously no convergence result in the trust region framework was known. We also show that certain splitting techniques can be seen as special cases of bundle trust region techniques.


title: Fermionic origin of dark energy in the inflationary universe from Unified Spinor Fields

abstract: In this work we explore the boundary conditions in the Einstein-Hilbert action, by considering a displacement from the Riemannian manifold to an extended one. The latter is characterized by including spinor fields into the quantum geometric description of a noncommutative spacetime. These fields are defined on the background spacetime, emerging from the expectation value of the quantum structure of spacetime generated by matrices that comply with a Clifford algebra. We demonstrate that spinor fields are candidate to describe all known interactions in physics, with gravitation included. In this framework we demonstrate that the cosmological constant $\Lambda$, is originated exclusively by massive fermion fields that would be the primordial components of dark energy, during the inflationary expansion of an universe that describes a de Sitter expansion.


title: Extensions of the Schwarzschild solution into regions of non-zero energy density and pressure

abstract: We present solutions of the Einstein equations that extend the static Schwarzschild solution in empty space into regions of non-zero energy density $\rho$ and radial pressure $p= w \rho$, where $w$ is a constant equation of state parameter. For simplicity we focus mainly on solutions with constant $\rho$. For $w=0$ we find solutions both with and without a singularity at the origin. Possible applications to galaxies are considered, where we find enhanced velocity rotation curves towards the edge of a galaxy. We propose that our explicit non-singular solution with $w=-1$ describes the interior of a black hole, which is a form of vacuum energy, and verify that its entropy is consistent with the Bekenstein-Hawking entropy. We propose that this idea can perhaps be applied to dark energy, if one views the latter as arising from black holes as pockets of vacuum energy. We estimate the density of the resulting dark energy to be $\rho_\Lambda \approx 10^{-30} g/{\rm cm}^3$, which is close to the measured value for the observable universe. We also present solutions with non-constant $\rho \propto 1/r^2$.


title: Nonclassical primordial gravitational waves from the initial entangled state

abstract: The nonclassicality of primordial gravitational waves (PGWs) is characterized in terms of sub-Poissonian graviton statistics. The sub-Poissonian statistics are realized when quantum states are squeezed coherent states. In the presence of matter fields, the Universe experiences the squeezed coherent state during inflation. The condition to realize the sub-Poissonian graviton statistics is translated into the frequency range of gravitational waves. If the initial state is the Bunch-Davies vacuum, there is another necessary condition between phases of squeezing and coherent parameters. Here, we extend the initial state to entangled states. We consider $\alpha$-vacua as the initial entangled state that are more general de Sitter invariant vacua than the Bunch-Davies vacuum. We find that, unlike the Bunch-Davies vacuum, PGWs generated in the initial entangled state become sub-Poissonian without requiring the condition between the phases.


title: Re-entrant phase separation in nematically aligning active polar particles

abstract: We present a numerical study of the phase behavior of repulsively interacting active polar particles that align their active velocities nematically. The amplitude of the active velocity, and the noise in its orientational alignment control the active nature of the system. At high values of orientational noise, the structural fluid undergoes a continuous nematic-isotropic transition in active orientation. This transition is well separated from an active phase separation, characterized by the formation of high density hexatic clusters, observed at lower noise strengths. With increasing activity, the system undergoes a re-entrant fluid-phase separation-fluid transition. The phase coexistence at low activity can be understood in terms of motility induced phase separation. In contrast, the re-melting of hexatic clusters, and the collective motion at low orientational noise are dominated by flocking behavior. At high activity, sliding and jamming of polar sub-clusters, formation of grain boundaries, lane formation, and subsequent fragmentation of the polar patches mediate remelting.


title: An Algebraic Approach to Fast Estimation of the Threshold Voltage of Junctionless Double Gate MOSFETs Using the Gram Schmidt Method

abstract: The effect of decreasing Drain-Induced Barrier Lowering (DIBL) is one of the non-desirable short-channel effects in the MOSFETs family, which causes the threshold voltage of the transistor to be reduced by increasing the voltage of the drain. This effect makes it impossible for circuit designers to consider VT as a constant value, and hence, it is necessary to calculate VT as a function of the drain voltage. Therefore, to consider the effect of DIBL in the design of integrated circuits, a large computational burden is imposed on the system, which slows down the simulation process in circuit-level simulators, particularly when a large number of transistors are to be simulated. Accordingly, in this paper, a multiple input single output (MISO) Nonlinear Autoregressive (N-AR) model using the Gram-Schmidt orthogonalization approach is proposed, that calculates the threshold voltage of the new generation of MOSFETs, i.e., Junctionless Double-Gate MOSFETs (JL-DG-MOSFETs), with a high precision and a significant speed-up in the computational procedure of the model. It is shown that, on average, the proposed numerical method is 313 times faster than the state-of-the-art analytical model. The calculated percentage of normalized mean square error between the proposed model and analytical one is 0.435% on average, showing that the proposed approach can be a fast and accurate candidate for replacing the analytical modeling.


title: How is Gaze Influenced by Image Transformations? Dataset and Model

abstract: Data size is the bottleneck for developing deep saliency models, because collecting eye-movement data is very time consuming and expensive. Most of current studies on human attention and saliency modeling have used high quality stereotype stimuli. In real world, however, captured images undergo various types of transformations. Can we use these transformations to augment existing saliency datasets? Here, we first create a novel saliency dataset including fixations of 10 observers over 1900 images degraded by 19 types of transformations. Second, by analyzing eye movements, we find that observers look at different locations over transformed versus original images. Third, we utilize the new data over transformed images, called data augmentation transformation (DAT), to train deep saliency models. We find that label preserving DATs with negligible impact on human gaze boost saliency prediction, whereas some other DATs that severely impact human gaze degrade the performance. These label preserving valid augmentation transformations provide a solution to enlarge existing saliency datasets. Finally, we introduce a novel saliency model based on generative adversarial network (dubbed GazeGAN). A modified UNet is proposed as the generator of the GazeGAN, which combines classic skip connections with a novel center-surround connection (CSC), in order to leverage multi level features. We also propose a histogram loss based on Alternative Chi Square Distance (ACS HistLoss) to refine the saliency map in terms of luminance distribution. Extensive experiments and comparisons over 3 datasets indicate that GazeGAN achieves the best performance in terms of popular saliency evaluation metrics, and is more robust to various perturbations. Our code and data are available at: https://github.com/CZHQuality/Sal-CFS-GAN.


title: A mathematical approach to quantum field theory

abstract: We develop a mathematical theory of quantization of multidimensional variational principles, and compare it with traditional constructions of quantum field theory. We conjecture that mathematical realization of quantum field theory axioms, in general, does not exist.


title: A duality of fields

abstract: It is shown that there exists a duality among fields. If a field is dual to another field, the solution of the field can be obtained from the dual field by the duality transformation. We give a general result on the dual fields. Different fields may have different numbers of dual fields, e.g., the free field and the $\phi^{4}$-field are self-dual, the $\phi^{n}$-field has one dual field, a field with an $n$-term polynomial potential has $n+1$ dual fields, and a field with a nonpolynomial potential may have infinite number of dual fields. All fields which are dual to each other form a duality family. This implies that the field can be classified in the sense of duality, or, the duality family defines a duality class. Based on the duality relation, we can construct a high-efficiency approach for seeking the solution of field equations: solving one field in the duality family, all solutions of other fields in the family are obtained immediately by the duality transformation. As examples, we consider some $\phi^{n}$-fields, general polynomial-potential fields, and the sine-Gordon field.


title: Cosmological Constraints on Chiral Tensor Particles

abstract: We discuss an extended model with chiral tensor particles in the Universe. Their direct influence on the Universe dynamics and their characteristic interactions in the hot Universe plasma, considered in previous publications, are briefly reviewed. A short discussion on the contemporary cosmological bounds on effective number of the relativistic degrees of freedom is provided. Cosmological constraints on the tensor particles interactions strength are obtained, corresponding to different cosmological bounds on the relativistic degrees of freedom and for different assumptions about right handed neutrinos.


title: Solar System and Atomic Stronger Bounds on Exotic Dyonic Matter and Non-associative Quantum Mechanics

abstract: We derive stronger bounds on the magnetic monopole charge of the proton, and hence on non-associative quantum mechanics from measurements of the magnetic fields of the Earth, the Moon and Mars. Limits are several orders of magnitude stronger than the ones inferred from the Hydrogen atomic levels, assuming that electrons do not retain a magnetic monopole charge. Conversely, we show how to estimate bounds on the magnetic charge of a proton correctly, when the magnetic monopole charge of the electron is taken into account.


title: The Effect of Anisotropy on Holographic Entanglement Entropy and Mutual Information

abstract: We study the effect of anisotropy on holographic entanglement entropy (HEE) and holographic mutual information (MI) in the Q-lattice model, by exploring the HEE and MI for infinite strips along arbitrary directions. We find that the lattice always enhances the HEE. The MI, however, is enhanced by lattice for large sub-regions; while for small sub-regions, the MI is suppressed by the lattice. We also discuss how these phenomena result from the deformation of geometry caused by Q-lattices.


title: Occupancy Estimation Using Low-Cost Wi-Fi Sniffers

abstract: Real-time measurements on the occupancy status of indoor and outdoor spaces can be exploited in many scenarios (HVAC and lighting system control, building energy optimization, allocation and reservation of spaces, etc.). Traditional systems for occupancy estimation rely on environmental sensors (CO2, temperature, humidity) or video cameras. In this paper, we depart from such traditional approaches and propose a novel occupancy estimation system which is based on the capture of Wi-Fi management packets from users' devices. The system, implemented on a low-cost ESP8266 microcontroller, leverages a supervised learning model to adapt to different spaces and transmits occupancy information through the MQTT protocol to a web-based dashboard. Experimental results demonstrate the validity of the proposed solution in four different indoor university spaces.


title: Predicting Dynamic Modulus of Asphalt Mixture Using Data Obtained from Indirect Tension Mode of Testing

abstract: Understanding stress-strain behavior of asphalt pavement under repetitive traffic loading is of critical importance to predict pavement performance and service life. For viscoelastic materials, the stress-strain relationship can be represented by the dynamic modulus. The dynamic modulus test in indirect tension mode can be used to measure the modulus of each specific layer of asphalt pavements using representative samples. Dynamic modulus is a function of material properties, loading, and environmental conditions. Developing predictive models for dynamic modulus is efficient and cost effective. This article focuses on developing an accurate Finite Element (FE) model using mixture elastic modulus and asphalt binder properties to predict dynamic modulus of asphalt mix in indirect tension mode. An Artificial Neural Network (ANN) is used to back-calculate the elastic modulus of asphalt mixtures. The developed FE model was verified against experimental results of field cores from nine different pavement sections from five districts in the State of Minnesota. It is demonstrated that the ANN modeling is a powerful tool to back-calculate the elastic modulus and FE model is capable of accurately predicting dynamic modulus.


title: Modern theory of hydraulic fracture modeling with using explicit and implicit schemes

abstract: The paper presents novel results, obtained on the basis of the modified theory of hydraulic fractures (HF). The theory underlines significance of the speed equation. When applied to numerical simulation of HF, the theory reveals three distinct issues: (i) modeling the central part of a HF; (ii) modeling the near-front zone; and (iii) tracing changes in the shape of a fracture contour. Modeling the central part leads to a stiff system of ODE in time, what strongly complicates its integration. For explicit schemes, it requires small time steps to meet the CFL condition. For implicit schemes, it requires proper preconditioners. The gains and flaws of the two strategies are discussed. It is noted that a rough spatial mesh may be used in the central part. Modeling the near-front zone reveals the vital role of the speed equation for HF modeling by any method. Its asymptotic analysis has resulted in the fundamental concept of the universal asymptotic umbrella. For the near-front zone, it is also established that a notable part of the zone adjacent to the front propagates virtually as a simple wave. This implies that the CFL condition of stability for this zone is much less restrictive than for the central part of the fracture. Of essence is also that the wave-like propagation of the near front zone makes preferable upwind schemes of time stepping. On whole, the analysis implies that explicit time stepping may be complementing, competitive and even superior over implicit integration. Tracing changes in the front shape appears merely in 3D problems when the contour changes its form. The problem, being essentially geometrical, it may be solved separately by various methods, including fast marching, level set and the simplest string/marker methods. In 1D cases, it does not arise at all. Quantitative estimations and numerical examples illustrate the theoretical conclusions.


title: Statistical Analysis and Modeling of the Geometry and Topology of Plant Roots

abstract: The root is an important organ of a plant since it is responsible for water and nutrient uptake. Analyzing and modelling variabilities in the geometry and topology of roots can help in assessing the plant's health, understanding its growth patterns, and modeling relations between plant species and between plants and their environment. In this article, we develop a framework for the statistical analysis and modeling of the geometry and topology of plant roots. We represent root structures as points in a tree-shape space equipped with a metric that quantifies geometric and topological differences between pairs of roots. We then use these building blocks to compute geodesics, i.e., optimal deformations under the metric between root structures, and to perform statistical analysis on root populations. We demonstrate the utility of the proposed framework through an application to a dataset of wheat roots grown in different environmental conditions. We also show that the framework can be used in various applications including classification and regression.


title: A Multiparametric Quon Algebra

abstract: The quon algebra is an approach to particle statistics introduced by Greenberg in order to provide a theory in which the Pauli exclusion principle and Bose statistics are violated by a small amount. We generalize these models by introducing a deformation of the quon algebra generated by a collection of operators $\mathtt{a}_i$, $i \in \mathbb{N}^*$ the set of positive integers, on an infinite dimensional module satisfying the $q_{i,j}$-mutator relations $\mathtt{a}_i \mathtt{a}_j^{\dag} - q_{i,j}\, \mathtt{a}_j^{\dag} \mathtt{a}_i = \delta_{i,j}$. The realizability of our model is proved by means of the Aguiar-Mahajan bilinear form on the chambers of hyperplane arrangements. We show that, for suitable values of $q_{i,j}$, the module generated by the particle states obtained by applying combinations of $\mathtt{a}_i$'s and $\mathtt{a}_i^{\dag}$'s to a vacuum state $|0\rangle$ is an indefinite Hilbert module. Furthermore, we refind the extended Zagier's conjecture established independently by Meljanac et al. and by Duchamp et al.


title: Experimental investigation of nonlocal advantage of quantum coherence

abstract: Nonlocal advantage of quantum coherence (NAQC) based on coherence complementarity relations is generally viewed as a stronger nonclassical correlation than Bell nonlocality. An arbitrary two-qubit state with NAQC must be an entangled state, which demonstrates that the criterion of NAQC can also be regarded as an entanglement witness. In this paper, we experimentally investigate the NAQC for Bell-diagonal states with high fidelity in an optics-based platform. We perform local measurements on a subsystem in three mutually unbiased bases and reconstruct the density matrices of the measured states by quantum state tomography process. By analyzing characteristic of the $l_1$ norm, relative entropy and skew information of coherence with parameters of quantum states, NAQC for the quantum states is accurately captured, and it shows that our experimental results are well compatible with the theoretical predictions. It is worth mentioning that quantum states with NAQC would have higher entanglement, and thus NAQC could be expected to be a kind of useful physical resource for quantum information processing.


title: Yang-Baxter random fields and stochastic vertex models

abstract: Bijectivization refines the Yang-Baxter equation into a pair of local Markov moves which randomly update the configuration of the vertex model. Employing this approach, we introduce new Yang-Baxter random fields of Young diagrams based on spin $q$-Whittaker and spin Hall-Littlewood symmetric functions. We match certain scalar Markovian marginals of these fields with (1) the stochastic six vertex model; (2) the stochastic higher spin six vertex model; and (3) a new vertex model with pushing which generalizes the $q$-Hahn PushTASEP introduced recently by Corwin-Matveev-Petrov (arXiv:1811.06475). Our matchings include models with two-sided stationary initial data, and we obtain Fredholm determinantal expressions for the $q$-Laplace transforms of the height functions of all these models. Moreover, we also discover difference operators acting diagonally on spin $q$-Whittaker or (stable) spin Hall-Littlewood symmetric functions.

categories: math.PR cond-mat.stat-mech math-ph math.CO math.MP math.QA


title: Characteristics of extensive air showers around the energy threshold for ground-particle-based gamma-ray observatories

abstract: Very high energy gamma-ray astronomy based on the measurement of air shower particles at ground-level has only recently been established as a viable approach, complementing the well established air Cherenkov technique. This approach requires high (mountain) altitudes and very high surface coverage particle detectors. While in general the properties of air showers are well established for many decades, the extreme situation of ground-level detection of very small showers from low energy primaries has not yet been well characterised for the purposes of gamma-ray astronomy. Here we attempt such a characterisation, with the aim of supporting the optimisation of next-generation gamma-ray observatories based on this technique. We address all of the key ground level observables and provide parameterisations for use in detector optimisation for shower energies around 1 TeV. We emphasise two primary aspects: the need for large area detectors to effectively measure low-energy showers, and the importance of muon identification for the purpose of background rejection.


title: Learning to Regress 3D Face Shape and Expression from an Image without 3D Supervision

abstract: The estimation of 3D face shape from a single image must be robust to variations in lighting, head pose, expression, facial hair, makeup, and occlusions. Robustness requires a large training set of in-the-wild images, which by construction, lack ground truth 3D shape. To train a network without any 2D-to-3D supervision, we present RingNet, which learns to compute 3D face shape from a single image. Our key observation is that an individual's face shape is constant across images, regardless of expression, pose, lighting, etc. RingNet leverages multiple images of a person and automatically detected 2D face features. It uses a novel loss that encourages the face shape to be similar when the identity is the same and different for different people. We achieve invariance to expression by representing the face using the FLAME model. Once trained, our method takes a single image and outputs the parameters of FLAME, which can be readily animated. Additionally we create a new database of faces `not quite in-the-wild' (NoW) with 3D head scans and high-resolution images of the subjects in a wide variety of conditions. We evaluate publicly available methods and find that RingNet is more accurate than methods that use 3D supervision. The dataset, model, and results are available for research purposes at http://ringnet.is.tuebingen.mpg.de.


title: A Search for High-Energy Counterparts to Fast Radio Bursts

abstract: We report on a search for high-energy counterparts to fast radio bursts (FRBs) with the Fermi Gamma-ray Burst Monitor (GBM), Fermi Large Area Telescope (LAT), and the Neil Gehrels Swift Observatory Burst Alert Telescope (BAT). We find no significant associations for any of the 23 FRBs in our sample, but report upper limits to the high-energy fluence for each on timescales of 0.1, 1, 10, and 100 s. We report lower limits on the ratio of the radio to high-energy fluence, $\frac{f_{r}}{f_{\gamma}}$, for timescales of 0.1 and 100 s. We discuss the implications of our non-detections on various proposed progenitor models for FRBs, including analogs of giant pulses from the Crab pulsar and hyperflares from magnetars. This work demonstrates the utility of analyses of high-energy data for FRBs in tracking down the nature of these elusive sources.


title: Outstanding strength, optical characteristics and thermal conductivity of graphene-like BC$_3$ and BC$_6$N semiconductors

abstract: Carbon based two-dimensional (2D) materials with honeycomb lattices, like graphene, polyaniline carbon-nitride (C$_3$N) and boron-carbide (BC$_3$) exhibit exceptional physical properties. On this basis, we propose two novel graphene-like materials with BC$_6$N stoichiometry. We conducted first-principles calculations to explore the stability, mechanical response, electronic, optical and thermal transport characteristics of graphene-like BC$_3$ and BC$_6$N monolayers. The absence of imaginary frequencies in the phonon dispersions confirm dynamical stability of BC$_3$ and BC$_6$N monolayers. Our first principles results reveal that BC$_3$ and BC$_6$N present high elastic moduli of 256 and 305 N/m, and tensile strengths of 29.0 and 33.4 N/m, with room temperature lattice thermal conductivities of 410 and 1710 W/m.K, respectively. Notably, the thermal conductivity of BC$_6$N is one of the highest among all 2D materials. According to electronic structure calculations, monolayers of BC$_3$ and BC$_6$N are indirect and direct bandgap semiconductors, respectively. The optical analysis illustrate that the first absorption peaks along the in-plane polarization for single-layer BC$_3$ and BC$_6$N occur in the visible range of the electromagnetic spectrum. Our results reveal outstandingly high mechanical properties and thermal conductivity along with attractive electronic and optical features of BC$_3$ and BC$_6$N nanosheets and present them as promising candidates to design novel nanodevices.


title: Dealing with Label Scarcity in Computational Pathology: A Use Case in Prostate Cancer Classification

abstract: Large amounts of unlabelled data are commonplace for many applications in computational pathology, whereas labelled data is often expensive, both in time and cost, to acquire. We investigate the performance of unsupervised and supervised deep learning methods when few labelled data are available. Three methods are compared: clustering autoencoder latent vectors (unsupervised), a single layer classifier combined with a pre-trained autoencoder (semi-supervised), and a supervised CNN. We apply these methods on hematoxylin and eosin (H&E) stained prostatectomy images to classify tumour versus non-tumour tissue. Results show that semi-/unsupervised methods have an advantage over supervised learning when few labels are available. Additionally, we show that incorporating immunohistochemistry (IHC) stained data provides an increase in performance over only using H&E.


title: Adaptive Sensor Placement for Continuous Spaces

abstract: We consider the problem of adaptively placing sensors along an interval to detect stochastically-generated events. We present a new formulation of the problem as a continuum-armed bandit problem with feedback in the form of partial observations of realisations of an inhomogeneous Poisson process. We design a solution method by combining Thompson sampling with nonparametric inference via increasingly granular Bayesian histograms and derive an $\tilde{O}(T^{2/3})$ bound on the Bayesian regret in $T$ rounds. This is coupled with the design of an efficent optimisation approach to select actions in polynomial time. In simulations we demonstrate our approach to have substantially lower and less variable regret than competitor algorithms.


title: Quantifying the Plasmonic Character of Optical Excitations in a Molecular J-Aggregate

abstract: The definition of plasmon at the microscopic scale is far from being understood. Yet, it is very important to recognize plasmonic features in optical excitations, as they can inspire new applications and trigger new discoveries by analogy with the rich phenomenology of metal nanoparticle plasmons. Recently, the concepts of plasmonicity index and the generalized plasmonicity index (GPI) have been devised as computational tools to quantify the plasmonic nature of optical excitations. The question may arise whether any strong absorption band, possibly with some sort of collective character in its microscopic origin, shares the status of plasmon. Here we demonstrate that this is not always the case, by considering a well-known class of systems represented by J-aggregates molecular crystals, characterized by the intense J band of absorption. By means of first-principles simulations, based on a many-body perturbation theory formalism, we investigate the optical properties of a J-aggregate made of push-pull organic dyes. We show that the effect of aggregation is to lower the GPI associated with the J-band with respect to the isolated dye one, which corresponds to a nonplasmonic character of the electronic excitations. In order to rationalize our finding, we then propose a simplified one-dimensional theoretical model of the J-aggregate. A useful microscopic picture of what discriminates a collective molecular crystal excitation from a plasmon is eventually obtained.


title: Ricci Reheating

abstract: We present a model for viable gravitational reheating involving a scalar field directly coupled to the Ricci curvature scalar. Crucial to the model is a period of kination after inflation, which causes the Ricci scalar to change sign thus inducing a tachyonic effective mass $m^{2} \propto -H^2$ for the scalar field. The resulting tachyonic growth of the scalar field provides the energy for reheating, allowing for temperatures high enough for thermal leptogenesis. Additionally, the required period of kination necessarily leads to a blue-tilted primordial gravitational wave spectrum with the potential to be detected by future experiments. We find that for reheating temperatures $T_{\rm RH} \lesssim 1$ GeV, the possibility exists for the Higgs field to play the role of the scalar field.


title: Sample Paths Estimates for Stochastic Fast-Slow Systems driven by Fractional Brownian Motion

abstract: We analyze the effect of additive fractional noise with Hurst parameter $H > \frac{1}{2}$ on fast-slow systems. Our strategy is based on sample paths estimates, similar to the approach by Berglund and Gentz in the Brownian motion case. Yet, the setting of fractional Brownian motion does not allow us to use the martingale methods from fast-slow systems with Brownian motion. We thoroughly investigate the case where the deterministic system permits a uniformly hyperbolic stable slow manifold. In this setting, we provide a neighborhood, tailored to the fast-slow structure of the system, that contains the process with high probability. We prove this assertion by providing exponential error estimates on the probability that the system leaves this neighborhood. We also illustrate our results in an example arising in climate modeling, where time-correlated noise processes have become of greater relevance recently.


title: Fast TLB Simulation for RISC-V Systems

abstract: Address translation and protection play important roles in today's processors, supporting multiprocessing and enforcing security. Historically, the design of the address translation mechanisms has been closely tied to the instruction set. In contrast, RISC-V defines its privileged specification in a way that permits a variety of designs. An important part of the design space is the organisation of Translation Lookaside Buffers (TLBs). This paper presents our recent work on simulating TLB behaviours in multi-core RISC-V systems. Our TLB simulation framework allows rapid, flexible and versatile prototyping of various hardware TLB design choices, and enables validation, profiling and benchmarking of software running on RISC-V systems. We show how this framework can be integrated with the dynamic binary translated emulator QEMU to perform online simulation. When simulating complicated multi-level shared TLB designs, the framework runs at around 400 million instructions per second (MIPS) when simulating an 8-core system. The performance overhead compared to unmodified QEMU is only 18% when the benchmark's L1 TLB miss rate is 1%. We also demonstrate how this tool can be used to explore the instruction-set level design space. We test a shared last-level TLB design that is not currently permitted by the RISC-V's privileged specification. We then propose an extension to RISC-V's virtual memory system design based on these experimental results.


title: A chemical kinetics code for modelling exoplanetary atmospheres

abstract: Chemical compositions of exoplanets can provide key insights into their physical processes, and formation and evolutionary histories. Atmospheric spectroscopy provides a direct avenue to probe exoplanetary compositions. However, whether obtained in transit or thermal emission, spectroscopic observations probe limited pressure windows of planetary atmospheres and are directly sensitive to only a limited set of spectroscopically active species. It is therefore critical to have chemical models that can relate retrieved atmospheric compositions to an atmosphere's bulk physical and chemical state. To this end we present LEVI, a new chemical kinetics code for modelling exoplanetary atmospheres. LEVI calculates the gas phase hydrogen, oxygen, carbon, and nitrogen chemistry in planetary atmospheres. Here we focus on hot gas giants. Applying LEVI, we investigate how variations in bulk C/O and N/O affects the observable atmospheric chemistry in hot Jupiters. For typical hot Jupiters we demonstrate the strong sensitivity of molecular detections to the atmospheric C/O. Molecular detections are conversely less sensitive to the atmospheric N/O ratio, although highly super-solar N/O can decrease the C/O required for HCN and NH3 detection. Using a new P-T profile for HD 209458b without a thermal inversion, we evaluate recently reported detection's of CO, H2O and HCN in its day-side atmosphere. We find that our models are consistent with the detected species, albeit with a narrow compositional window around C/O $\sim$ 1. A C/O $\gtrsim$ 0.9 (1.6 times solar) was required to meet the minimum reported value for HCN, while a C/O $\lesssim$ 1 (1.8 times solar) was required to fit the nominal H2O abundance.


title: The Balian-Low theorem for locally compact abelian groups and vector bundles

abstract: Let $\Lambda$ be a lattice in a second countable, locally compact abelian group $G$ with annihilator $\Lambda^{\perp} \subseteq \widehat{G}$. We investigate the validity of the following statement: For every $\eta$ in the Feichtinger algebra $S_0(G)$, the Gabor system $\{ M_{\tau} T_{\lambda} \eta \}_{\lambda \in \Lambda, \tau \in \Lambda^{\perp}}$ is not a frame for $L^2(G)$. When $G = \mathbb{R}$ and $\Lambda = \alpha \mathbb{Z}$, this statement is a variant of the Balian-Low theorem. Extending a result of R. Balan, we show that whether the statement generalizes to $(G,\Lambda)$ is equivalent to the nontriviality of a certain vector bundle over the compact space $(G/\Lambda) \times (\widehat{G}/\Lambda^{\perp})$. We prove this equivalence using a connection between Gabor frames and Heisenberg modules. More specifically, we show that the Zak transform can be viewed as an isomorphism of certain Hilbert $C^*$-modules. As an application, we prove a new Balian-Low theorem for the group $\mathbb{R} \times \mathbb{Q}_p$, where $\mathbb{Q}_p$ denotes the $p$-adic numbers.


title: Convergence of Heuristic Parameter Choice Rules for Convex Tikhonov Regularisation

abstract: We investigate the convergence theory of several known as well as new heuristic parameter choice rules for convex Tikhonov regularisation. The success of such methods is dependent on whether certain restrictions on the noise are satisfied. In the linear theory, such conditions are well understood and hold for typically irregular noise. In this paper, we extend the convergence analysis of heuristic rules using noise restrictions to the convex setting and prove convergence of the aforementioned methods therewith. The convergence theory is exemplified for the case of an ill-posed problem with a diagonal forward operator in $\ell^q$ spaces. Numerical examples also provide further insight.


title: Stochastic precedence and minima among dependent variables

abstract: The notion of stochastic precedence between two random variables emerges as a relevant concept in several fields of applied probability. When one consider a vector of random variables $X_1,...,X_n$, this notion has a preeminent role in the analysis of minima of the type $\min_{j \in A} X_j$ for $A \subset \{1, \ldots n\}$. In such an analysis, however, several apparently controversial aspects can arise (among which phenomena of "non-transitivity"). Here we concentrate attention on vectors of non-negative random variables with absolutely continuous joint distributions, in which a case the set of the multivariate conditional hazard rate (m.c.h.r.) functions can be employed as a convenient method to describe different aspects of stochastic dependence. In terms of the m.c.h.r. functions, we first obtain convenient formulas for the probability distributions of the variables $\min_{j \in A} X_j$ and for the probability of events $\{X_i=\min_{j \in A} X_j\}$. Then we detail several aspects of the notion of stochastic precedence. On these bases, we explain some controversial behavior of such variables and give sufficient conditions under which paradoxical aspects can be excluded. On the purpose of stimulating active interest of readers, we present several comments and pertinent examples.


title: Virtual classes of $\mathbb{G}_\text{m}$-gerbes

abstract: We show that a perfect obstruction theory for a $\mathbb{G}_\text{m}$-gerbe determines a semi-perfect obstruction theory for its base, which is perfect if the gerbe is quasi-compact and affine-pointed. These results streamline the construction of a semi-perfect obstruction theory for the base, and allow us to relate virtual classes of the gerbe and its base.


title: Towards Interlingua Neural Machine Translation

abstract: Common intermediate language representation in neural machine translation can be used to extend bilingual to multilingual systems by incremental training. In this paper, we propose a new architecture based on introducing an interlingual loss as an additional training objective. By adding and forcing this interlingual loss, we are able to train multiple encoders and decoders for each language, sharing a common intermediate representation. Translation results on the low-resourced tasks (Turkish-English and Kazakh-English tasks, from the popular Workshop on Machine Translation benchmark) show the following BLEU improvements up to 2.8. However, results on a larger dataset (Russian-English and Kazakh-English, from the same baselines) show BLEU loses if the same amount. While our system is only providing improvements for the low-resourced tasks in terms of translation quality, our system is capable of quickly deploying new language pairs without retraining the rest of the system, which may be a game-changer in some situations (i.e. in a disaster crisis where international help is required towards a small region or to develop some translation system for a client). Precisely, what is most relevant from our architecture is that it is capable of: (1) reducing the number of production systems, with respect to the number of languages, from quadratic to linear (2) incrementally adding a new language in the system without retraining languages previously there and (3) allowing for translations from the new language to all the others present in the system


title: Space and Time Complexity for Infinite Time Turing Machines

abstract: We consider notions of space complexity for Infinite Time Turing Machines (ITTMs) that were introduced by B. L\"owe and studied further by J. Winter. We answer several open questions about these notions, among them whether low space complexity implies low time complexity (it does not) and whether one of the equalities P=PSPACE, P$_{+}=$PSPACE$_{+}$ and P$_{++}=$PSPACE$_{++}$ holds for ITTMs (all three are false). We also show various separation results between space complexity classes for ITTMs. This considerably expands our earlier observations on the topic in section 7.2.2 of \cite{Ca2}, which appear here as Lemma $6$ up to Corollary $9$.


title: Towards the Efficient Local Tailored Coupled Cluster Approximation and the Peculiar Case of Oxo-Mn(Salen)

abstract: We introduce a new implementation of the coupled cluster method tailored by matrix product states wave functions (DMRG-TCCSD), which employs the local pair natural orbital approach (LPNO). By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of canonical correlation energy for the triplet and singlet state of TME respectively. In case of oxo-Mn(Salen), we found out that the amount of retrieved canonical correlation energy depends on the size of the active space (CAS) - we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-$\zeta$ basis set amounting to 1178 basis functions. Moreover, we examined dependance of the ground state of oxo-Mn(Salen) on CAS composition. We found out that the inclusion of 4d$_{xy}$ orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation the ground state was found to be triplet regardless of the size of the basis set or composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.


title: A complex analysis approach to Atangana-Baleanu fractional calculus

abstract: The standard definition for the Atangana-Baleanu fractional derivative involves an integral transform with a Mittag-Leffler function in the kernel. We show that this integral can be rewritten as a complex contour integral which can be used to provide an analytic continuation of the definition to complex orders of differentiation. We discuss the implications and consequences of this extension, including a more natural formula for the Atangana-Baleanu fractional integral and for iterated Atangana-Baleanu fractional differintegrals.


title: Constraining a black hole companion for M87* through imaging by the Event Horizon Telescope

abstract: The Event Horizon Telescope (EHT), a global very long baseline interferometric array observing at a wavelength of 1.3 mm, detected the first image of the M87 supermassive black hole (SMBH). M87 is a giant elliptical galaxy at the center of Virgo cluster, which is expected to have formed through merging of cluster galaxies. Consequently M87* hosted mergers of black holes through dynamical friction and could have one or multiple binary companions with a low mass ratio at large separations. We show that a long-term monitoring of the M87 SMBH image over $\sim$1 year with absolute positional accuracy of 1$\approx\mu as$ could detect such binary companions and exclude a large parameter space in semi major axis ($a_0$) and mass ratio ($q$), which is currently not constrained. Moreover, the presence of the accretion disk around M87* excludes a binary companion with $a_0\approx$ of order a mili parsec, as otherwise the accretion disk would have been tidally disrupted.


title: Stability of Linear Structural Equation Models of Causal Inference

abstract: We consider the numerical stability of the parameter recovery problem in Linear Structural Equation Model ($\LSEM$) of causal inference. A long line of work starting from Wright (1920) has focused on understanding which sub-classes of $\LSEM$ allow for efficient parameter recovery. Despite decades of study, this question is not yet fully resolved. The goal of this paper is complementary to this line of work; we want to understand the stability of the recovery problem in the cases when efficient recovery is possible. Numerical stability of Pearl's notion of causality was first studied in Schulman and Srivastava (2016) using the concept of condition number where they provide ill-conditioned examples. In this work, we provide a condition number analysis for the $\LSEM$. First we prove that under a sufficient condition, for a certain sub-class of $\LSEM$ that are \emph{bow-free} (Brito and Pearl (2002)), the parameter recovery is stable. We further prove that \emph{randomly} chosen input parameters for this family satisfy the condition with a substantial probability. Hence for this family, on a large subset of parameter space, recovery is numerically stable. Next we construct an example of $\LSEM$ on four vertices with \emph{unbounded} condition number. We then corroborate our theoretical findings via simulations as well as real-world experiments for a sociology application. Finally, we provide a general heuristic for estimating the condition number of any $\LSEM$ instance.


title: Tracing cultural diachronic semantic shifts in Russian using word embeddings: test sets and baselines

abstract: The paper introduces manually annotated test sets for the task of tracing diachronic (temporal) semantic shifts in Russian. The two test sets are complementary in that the first one covers comparatively strong semantic changes occurring to nouns and adjectives from pre-Soviet to Soviet times, while the second one covers comparatively subtle socially and culturally determined shifts occurring in years from 2000 to 2014. Additionally, the second test set offers more granular classification of shifts degree, but is limited to only adjectives. The introduction of the test sets allowed us to evaluate several well-established algorithms of semantic shifts detection (posing this as a classification problem), most of which have never been tested on Russian material. All of these algorithms use distributional word embedding models trained on the corresponding in-domain corpora. The resulting scores provide solid comparison baselines for future studies tackling similar tasks. We publish the datasets, code and the trained models in order to facilitate further research in automatically detecting temporal semantic shifts for Russian words, with time periods of different granularities.


title: Pushing the Limits of the Coronagraphic Occulters on HST/STIS

abstract: The Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) contains the only currently operating coronagraph in space that is not trained on the Sun. In an era of extreme--adaptive-optics--fed coronagraphs, and with the possibility of future space-based coronagraphs, we re-evaluate the contrast performance of the STIS CCD camera. The 50CORON aperture consists of a series of occulting wedges and bars, including the recently commissioned BAR5 occulter. We discuss the latest procedures in obtaining high contrast imaging of circumstellar disks and faint point sources with STIS. For the first time, we develop a noise model for the coronagraph, including systematic noise due to speckles, which can be used to predict the performance of future coronagraphic observations. Further, we present results from a recent calibration program that demonstrates better than $10^{-6}$ point-source contrast at 0.6", ranging to $3\times10^{-5}$ point-source contrast at 0.25". These results are obtained by a combination of sub-pixel grid dithers, multiple spacecraft orientations, and post-processing techniques. Some of these same techniques will be employed by future space-based coronagraphic missions. We discuss the unique aspects of STIS coronagraphy relative to ground-based adaptive-optics--fed coronagraphs.


title: A numerical approach for solving fractional optimal control problems using modified hat functions

abstract: We introduce a numerical method, based on modified hat functions, for solving a class of fractional optimal control problems. In our scheme, the control and the fractional derivative of the state function are considered as linear combinations of the modified hat functions. The fractional derivative is considered in the Caputo sense while the Riemann-Liouville integral operator is used to give approximations for the state function and some of its derivatives. To this aim, we use the fractional order integration operational matrix of the modified hat functions and some properties of the Caputo derivative and Riemann-Liouville integral operators. Using results of the considered basis functions, solving the fractional optimal control problem is reduced to the solution of a system of nonlinear algebraic equations. An error bound is proved for the approximate optimal value of the performance index obtained by the proposed method. The method is then generalized for solving a class of fractional optimal control problems with inequality constraints. The most important advantages of our method are easy implementation, simple operations, and elimination of numerical integration. Some illustrative examples are considered to demonstrate the effectiveness and accuracy of the proposed technique.


title: Violation of the mean path length invariance property

abstract: The invariance property of the mean path length is an astonishing law of Nature governing the motion of particles inside a disordered material. Whatever the strength of the disorder, the property states that the mean path length is exclusively determined by the ratio between the volume and the surface. Till now, the property has been reported as universal and valid in any kind of disordered medium and also beyond diffusion conditions. Nevertheless, we found out that the property fails in anomalous transport and in other kinds of random walk. By means of Monte Carlo simulations of light transport, we show that, in these cases, the invariance property loses its validity and the mean path length becomes dependent on the diffusive characteristics of the medium. The critical issue of such a violation lies in the breaking of isotropy and homogeneity of the radiance in the whole volume. These results are valid for all natural or artificial phenomena where random walkers, whatever their nature, are able to experiment anomalous transport.


title: Enforcing Statistical Constraints in Generative Adversarial Networks for Modeling Chaotic Dynamical Systems

abstract: Simulating complex physical systems often involves solving partial differential equations (PDEs) with some closures due to the presence of multi-scale physics that cannot be fully resolved. Therefore, reliable and accurate closure models for unresolved physics remains an important requirement for many computational physics problems, e.g., turbulence simulation. Recently, several researchers have adopted generative adversarial networks (GANs), a novel paradigm of training machine learning models, to generate solutions of PDEs-governed complex systems without having to numerically solve these PDEs. However, GANs are known to be difficult in training and likely to converge to local minima, where the generated samples do not capture the true statistics of the training data. In this work, we present a statistical constrained generative adversarial network by enforcing constraints of covariance from the training data, which results in an improved machine-learning-based emulator to capture the statistics of the training data generated by solving fully resolved PDEs. We show that such a statistical regularization leads to better performance compared to standard GANs, measured by (1) the constrained model's ability to more faithfully emulate certain physical properties of the system and (2) the significantly reduced (by up to 80%) training time to reach the solution. We exemplify this approach on the Rayleigh-Benard convection, a turbulent flow system that is an idealized model of the Earth's atmosphere. With the growth of high-fidelity simulation databases of physical systems, this work suggests great potential for being an alternative to the explicit modeling of closures or parameterizations for unresolved physics, which are known to be a major source of uncertainty in simulating multi-scale physical systems, e.g., turbulence or Earth's climate.


title: On the dynamics of non-rigid asteroid rotation

abstract: We have presented in this communication a new solving procedure for the dynamics of non-rigid asteroid rotation, considering the final spin state of rotation for a small celestial body (asteroid). The last condition means the ultimate absence of the applied external torques (including short-term effect from torques during collisions, long-term YORP effect, etc.). Fundamental law of angular momentum conservation has been used for the aforementioned solving procedure. The system of Euler equations for dynamics of non-rigid asteroid rotation has been explored with regard to the existence of an analytic way of presentation of the approximated solution. Despite of various perturbations (such as collisions, YORP effect) which destabilize the rotation of asteroid via deviating from the current spin state, the inelastic (mainly, tidal) dissipation reduces kinetic energy of asteroid. So, evolution of the spinning asteroid should be resulting by the rotation about maximal-inertia axis with the proper spin state corresponding to minimal energy with a fixed angular momentum. Basing on the aforesaid assumption (component K_1 is supposed to be fluctuating near the given appropriate constant of the fixed angular momentum), we have obtained that 2-nd component K_2 is the solution of appropriate Riccati ordinary differential equation of 1-st order, whereas component K_3 should be determined via expression for K_2.


title: Design and Experimental Validation of Tube-based MPC for Timed-constrained Robot Planning

abstract: This paper deals with the design and experimental validation of a state-of-the art tube-based Model Predictive Control (MPC) for achieving time-constrained tasks. Given the uncertain nonlinear dynamics of the robot as well as a high-level task written in Metric Interval Temporal Logic (MITL), the goal is to design a feedback control law that guarantees the satisfaction of the task. The workspace is divided into Regions of Interest (RoI) and contains also unsafe regions (obstacles) that the robot should not visit. The feedback control law consists of two terms: a control input which is the outcome of a Finite Horizon Optimal Control (FHOCP); and a state feedback law that guarantees that the nominal trajectories are bounded within a tube centered along the nominal trajectories. The aforementioned control law guarantees that the robot is safely navigated through the RoI within certain time bounds. The proposed framework can handle the rich expressiveness of MITL and is experimentally tested with a Nexus mobile robot in our lab facilities. The experimental results show that the proposed framework is promising for solving real-life robotic as well as industrial problems.


title: Reconfigurable Hardware Implementation of the Successive Overrelaxation Method

abstract: In this chapter, we study the feasibility of implementing SOR in reconfigurable hardware. We use Handel-C, a higher level design tool, to code our design, which is analyzed, synthesized, and placed and routed using the FPGAs proprietary software (DK Design Suite, Xilinx ISE 8.1i, and Quartus II 5.1). We target Virtex II Pro, Altera Stratix, and Spartan3L, which is embedded in the RC10 FPGA-based system from Celoxica. We report our timing results when targeting Virtex II Pro and compare them to software version results written in C++ and running on a general purpose processor (GPP).


title: Bit-Swap: Recursive Bits-Back Coding for Lossless Compression with Hierarchical Latent Variables

abstract: The bits-back argument suggests that latent variable models can be turned into lossless compression schemes. Translating the bits-back argument into efficient and practical lossless compression schemes for general latent variable models, however, is still an open problem. Bits-Back with Asymmetric Numeral Systems (BB-ANS), recently proposed by Townsend et al. (2019), makes bits-back coding practically feasible for latent variable models with one latent layer, but it is inefficient for hierarchical latent variable models. In this paper we propose Bit-Swap, a new compression scheme that generalizes BB-ANS and achieves strictly better compression rates for hierarchical latent variable models with Markov chain structure. Through experiments we verify that Bit-Swap results in lossless compression rates that are empirically superior to existing techniques. Our implementation is available at https://github.com/fhkingma/bitswap.

categories: cs.LG cs.AI cs.IT math.IT stat.CO stat.ML


title: Photo-induced two-body loss of ultracold molecules

abstract: The lifetime of nonreactive ultracold bialkali gases was conjectured to be limited by sticky collisions amplifying three-body loss. We show that the sticking times were previously overestimated and do not support this hypothesis. We find that electronic excitation of NaK+NaK collision complexes by the trapping laser leads to the experimentally observed two-body loss. We calculate the excitation rate with a quasiclassical, statistical model employing ab initio potentials and transition dipole moments. Using longer laser wavelengths or repulsive box potentials may suppress the losses.


title: Inferring Concise Specifications of APIs

abstract: Modern software relies on libraries and uses them via application programming interfaces (APIs). Correct API usage as well as many software engineering tasks are enabled when APIs have formal specifications. In this work, we analyze the implementation of each method in an API to infer a formal postcondition. Conventional wisdom is that, if one has preconditions, then one can use the strongest postcondition predicate transformer (SP) to infer postconditions. However, SP yields postconditions that are exponentially large, which makes them difficult to use, either by humans or by tools. Our key idea is an algorithm that converts such exponentially large specifications into a form that is more concise and thus more usable. This is done by leveraging the structure of the specifications that result from the use of SP. We applied our technique to infer postconditions for over 2,300 methods in seven popular Java libraries. Our technique was able to infer specifications for 75.7% of these methods, each of which was verified using an Extended Static Checker. We also found that 84.6% of resulting specifications were less than 1/4 page (20 lines) in length. Our technique was able to reduce the length of SMT proofs needed for verifying implementations by 76.7% and reduced prover execution time by 26.7%.


title: Magnetization dynamics modulated by Dzyaloshinskii-Moriya interaction in the double-interface spin transfer torque magnetic tunnel junction

abstract: Currently double-interface MTJs have been developed for enhancing the thermal stability barrier in small technology node. Dzyaloshinskii-Moriya interaction (DMI) inevitably exists in such devices due to the use of the heavy-metal/ferromagnet structures. Previous studies have demonstrated the detrimental effect of DMI on the conventional single-interface spin transfer torque (STT) MTJs. Here in this work we will prove that the detrimental effect of the DMI could be almost eliminated in the double-interface STT-MTJ. This conclusion is attributed to the suppressing effect of the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction on the DMI. Detailed mechanisms are analyzed based on the theoretical models and micromagnetic simulation results. Our work highlights the importance of appropriately controlling the DMI in two free layers of the double-interface STT-MTJ.


title: Observations of Ultrafast Kelvin Wave Breaking in the Mars Thermosphere

abstract: Ultrafast Kelvin Waves (UFKWs) have been recently discovered at Mars using measurements by the Mars Atmosphere and Evolution (MAVEN) Accelerometer (ACC) instrument and the Mars Reconnaissance Orbiter (MRO) Mars Climate Sound (MCS) instrument [1]. UFKWs are eastward propagating, have a 2-3 sol period with sufficiently long vertical wavelengths to propagate into the thermosphere, and are predicted by classical wave theory to be equatorially trapped. These prior measurements characterized UFKWs at two relatively narrow altitude bands, with one measurement altitude near 80 km and a second near 150 km. Thermospheric density profiles from solar occultation (SO) measurements made by the Extreme Ultraviolet Monitor (EUVM) onboard MAVEN [2] provide a means to characterize multi-sol period waves, including UFKWs, as a function of altitude. In this study, we present first-ever observations of UFKWs as a function of altitude, including observations of rapid wave energy dissipation near 170 km indicative of wave breaking.


title: Improving strong scaling of the Conjugate Gradient method for solving large linear systems using global reduction pipelining

abstract: This paper presents performance results comparing MPI-based implementations of the popular Conjugate Gradient (CG) method and several of its communication hiding (or 'pipelined') variants. Pipelined CG methods are designed to efficiently solve SPD linear systems on massively parallel distributed memory hardware, and typically display significantly improved strong scaling compared to classic CG. This increase in parallel performance is achieved by overlapping the global reduction phase (MPI_Iallreduce) required to compute the inner products in each iteration by (chiefly local) computational work such as the matrix-vector product as well as other global communication. This work includes a brief introduction to the deep pipelined CG method for readers that may be unfamiliar with the specifics of the method. A brief overview of implementation details provides the practical tools required for implementation of the algorithm. Subsequently, easily reproducible strong scaling results on the US Department of Energy (DoE) NERSC machine 'Cori' (Phase I - Haswell nodes) on up to 1024 nodes with 16 MPI ranks per node are presented using an implementation of p(l)-CG that is available in the open source PETSc library. Observations on the staggering and overlap of the asynchronous, non-blocking global communication phases with communication and computational kernels are drawn from the experiments.


title: Domino successive-deviation ghost imaging

abstract: Traditional ghost imaging acquires images via the correlation of the intensity fluctuations of reference patterns and bucket values, and can even generate positive-negative images by conditionally averaging partial patterns. Here, we propose a domino successive-deviation ghost imaging method, which owns a good image quality comparable to that of differential ghost imaging with real-time fast computation, and a better robustness in practical scenarios where measurement noise and the instability of illuminating source coexists. Furthermore, it happens to generate real-time positive and negative images, giving a new insight into physical essence of positive-negative ghost image phenomenon. Both simulation and experimental results have demonstrated the feasibility of our approach. Therefore, this work complements the theory of ghost imaging and opens a door to practical applications of real-time single-pixel imaging.