title: Interpretable Deep Learning for Two-Prong Jet Classification with Jet Spectra

abstract: Classification of jets with deep learning has gained significant attention in recent times. However, the performance of deep neural networks is often achieved at the cost of interpretability. Here we propose an interpretable network trained on the jet spectrum $S_{2}(R)$ which is a two-point correlation function of the jet constituents. The spectrum can be derived from a functional Taylor series of an arbitrary jet classifier function of energy flows. An interpretable network can be obtained by truncating the series. The intermediate feature of the network is an infrared and collinear safe C-correlator which allows us to estimate the importance of a $S_{2}(R)$ deposit at an angular scale R in the classification. The performance of the architecture is comparable to that of a convolutional neural network (CNN) trained on jet images, although the number of inputs and complexity of architecture is significantly simpler than the CNN classifier. We consider two examples: one is the classification of two-prong jets which differ in color charge of the mother particle, and the other is a comparison between Pythia 8 and Herwig 7 generated jets.


title: Short-baseline neutrino oscillations with 3+1 non-unitary mixing

abstract: We consider a scenario with unitary mixing of the three light standard neutrinos and a non-unitary mixing contribution of a heavier massive neutrino that can generate short-baseline neutrino oscillations. We show that this scenario predicts constant flavor-changing probabilities at short-baseline distances. Therefore, it cannot explain the spectral distortions observed in the LSND and MiniBooNE appearance experiments. On the other hand, the survival probabilities oscillate as functions of $L/E$ and could explain oscillations in short-baseline disappearance experiments. We also derive the bounds on the mixing parameters from the existing short-baseline neutrino oscillation data.


title: A violin sonata for reconnection

abstract: The process of magnetic reconnection when studied in Nature or when modeled in 3D simulations differs in one key way from the standard 2D paradigmatic cartoon: it is accompanied by much fluctuations in the electromagnetic fields and plasma properties. We developed a new diagnostics, the topographical fluctuations analysis (TFA) to study the spectrum of fluctuations in the various regions around a reconnection site. We find that fluctuations belong to two very different regimes. The first regime is better known, it develops in the reconnection outflows and is characterized by a strong link between plasma and electromagnetic fluctuations leading to momentum and energy exchanges via anomalous viscosity and resistivity. But there is a second, new, regime: it develops in the inflow and in the region around the separatrix surfaces, including the reconnection diffusion region itself. In this new regime the plasma remains laminar but the electromagnetic fields fluctuates strongly. We present an analogy with the smooth continuous motion of the bow of a violin producing the vibrations of the strings to emit music.


title: Discrimination through optimization: How Facebook's ad delivery can lead to skewed outcomes

abstract: The enormous financial success of online advertising platforms is partially due to the precise targeting features they offer. Although researchers and journalists have found many ways that advertisers can target---or exclude---particular groups of users seeing their ads, comparatively little attention has been paid to the implications of the platform's ad delivery process, comprised of the platform's choices about which users see which ads. It has been hypothesized that this process can "skew" ad delivery in ways that the advertisers do not intend, making some users less likely than others to see particular ads based on their demographic characteristics. In this paper, we demonstrate that such skewed delivery occurs on Facebook, due to market and financial optimization effects as well as the platform's own predictions about the "relevance" of ads to different groups of users. We find that both the advertiser's budget and the content of the ad each significantly contribute to the skew of Facebook's ad delivery. Critically, we observe significant skew in delivery along gender and racial lines for "real" ads for employment and housing opportunities despite neutral targeting parameters. Our results demonstrate previously unknown mechanisms that can lead to potentially discriminatory ad delivery, even when advertisers set their targeting parameters to be highly inclusive. This underscores the need for policymakers and platforms to carefully consider the role of the ad delivery optimization run by ad platforms themselves---and not just the targeting choices of advertisers---in preventing discrimination in digital advertising.


title: GEDI: Gammachirp Envelope Distortion Index for Predicting Intelligibility of Enhanced Speech

abstract: In this study, we propose a new concept, the gammachirp envelope distortion index (GEDI), based on the signal-to-distortion ratio in the auditory envelope, SDRenv to predict the intelligibility of speech enhanced by nonlinear algorithms. The objective of GEDI is to calculate the distortion between enhanced and clean-speech representations in the domain of a temporal envelope extracted by the gammachirp auditory filterbank and modulation filterbank. We also extend GEDI with multi-resolution analysis (mr-GEDI) to predict the speech intelligibility of sounds under non-stationary noise conditions. We evaluate GEDI in terms of speech intelligibility predictions of speech sounds enhanced by a classic spectral subtraction and a Wiener filtering method. The predictions are compared with human results for various signal-to-noise ratio conditions with additive pink and babble noises. The results showed that mr-GEDI predicted the intelligibility curves better than short-time objective intelligibility (STOI) measure, extended-STOI (ESTOI) measure, and hearing-aid speech perception index (HASPI) under pink-noise conditions, and better than HASPI under babble-noise conditions. The mr-GEDI method does not present an overestimation tendency and is considered a more conservative approach than STOI and ESTOI. Therefore, the evaluation with mr-GEDI may provide additional information in the development of speech enhancement algorithms.


title: Strong plasmon-molecule coupling at the nanoscale revealed by first-principles modeling

abstract: Strong light-matter interactions in both the single-emitter and collective strong coupling regimes attract significant attention due to emerging quantum and nonlinear optics applications, as well as opportunities for modifying material-related properties. Further exploration of these phenomena requires an appropriate theoretical methodology, which is demanding since polaritons are at the intersection between quantum optics, solid state physics and quantum chemistry. Fortunately, however, nanoscale polaritons can be realized in small plasmon-molecule systems, which in principle allows treating them using ab initio methods, although this has not been demonstrated to date. Here, we show that time-dependent density-functional theory (TDDFT) calculations can access the physics of nanoscale plasmon-molecule hybrids and predict vacuum Rabi splitting in a system comprising a few-hundred-atom aluminum nanoparticle interacting with one or several benzene molecules. We show that the cavity quantum electrodynamics approach holds down to resonators on the order of a few cubic nanometers, yielding a single-molecule coupling strength exceeding 200 meV due to a massive vacuum field value of 4.5 V/nm. In a broader perspective, our approach enables parameter-free in-depth studies of polaritonic systems, including ground state, chemical and thermodynamic modifications of the molecules in the strong-coupling regime, which may find important use in emerging applications such as cavity enhanced catalysis.


title: The Medical Deconfounder: Assessing Treatment Effects with Electronic Health Records

abstract: The treatment effects of medications play a key role in guiding medical prescriptions. They are usually assessed with randomized controlled trials (RCTs), which are expensive. Recently, large-scale electronic health records (EHRs) have become available, opening up new opportunities for more cost-effective assessments. However, assessing a treatment effect from EHRs is challenging: it is biased by unobserved confounders, unmeasured variables that affect both patients' medical prescription and their outcome, e.g. the patients' social economic status. To adjust for unobserved confounders, we develop the medical deconfounder, a machine learning algorithm that unbiasedly estimates treatment effects from EHRs. The medical deconfounder first constructs a substitute confounder by modeling which medications were prescribed to each patient; this substitute confounder is guaranteed to capture all multi-medication confounders, observed or unobserved (arXiv:1805.06826). It then uses this substitute confounder to adjust for the confounding bias in the analysis. We validate the medical deconfounder on two simulated and two real medical data sets. Compared to classical approaches, the medical deconfounder produces closer-to-truth treatment effect estimates; it also identifies effective medications that are more consistent with the findings in the medical literature.


title: 75 Languages, 1 Model: Parsing Universal Dependencies Universally

abstract: We present UDify, a multilingual multi-task model capable of accurately predicting universal part-of-speech, morphological features, lemmas, and dependency trees simultaneously for all 124 Universal Dependencies treebanks across 75 languages. By leveraging a multilingual BERT self-attention model pretrained on 104 languages, we found that fine-tuning it on all datasets concatenated together with simple softmax classifiers for each UD task can result in state-of-the-art UPOS, UFeats, Lemmas, UAS, and LAS scores, without requiring any recurrent or language-specific components. We evaluate UDify for multilingual learning, showing that low-resource languages benefit the most from cross-linguistic annotations. We also evaluate for zero-shot learning, with results suggesting that multilingual training provides strong UD predictions even for languages that neither UDify nor BERT have ever been trained on. Code for UDify is available at https://github.com/hyperparticle/udify.


title: Agent-based Vs Agent-less Sandbox for Dynamic Behavioral Analysis

abstract: Malicious software is detected and classified by either static analysis or dynamic analysis. In static analysis, malware samples are reverse engineered and analyzed so that signatures of malware can be constructed. These techniques can be easily thwarted through polymorphic, metamorphic malware, obfuscation and packing techniques, whereas in dynamic analysis malware samples are executed in a controlled environment using the sandboxing technique, in order to model the behavior of malware. In this paper, we have analyzed Petya, Spyeye, VolatileCedar, PAFISH etc. through Agent-based and Agentless dynamic sandbox systems in order to investigate and benchmark their efficiency in advanced malware detection.


title: The Landscape of R Packages for Automated Exploratory Data Analysis

abstract: The increasing availability of large but noisy data sets with a large number of heterogeneous variables leads to the increasing interest in the automation of common tasks for data analysis. The most time-consuming part of this process is the Exploratory Data Analysis, crucial for better domain understanding, data cleaning, data validation, and feature engineering. There is a growing number of libraries that attempt to automate some of the typical Exploratory Data Analysis tasks to make the search for new insights easier and faster. In this paper, we present a systematic review of existing tools for Automated Exploratory Data Analysis (autoEDA). We explore the features of twelve popular R packages to identify the parts of analysis that can be effectively automated with the current tools and to point out new directions for further autoEDA development.


title: Haldane and Dimer phases in a frustrated spin chain: an exact groundstate and associated topological phase transition

abstract: A Heisenberg spin-$s$ chain with alternating ferromagnetic ($-J_1^F<0$) and antiferromagnetic ($J_1^A>0$) nearest-neighbor (NN) interactions, exhibits the Dimer and spin-$2s$ Haldane phases in the limits $J_1^F/J_1^A \rightarrow 0$ and $J_1^F/J_1^A \rightarrow \infty$ respectively. These two phases are understood to be topologically equivalent. Induction of the frustration through the next nearest-neighbor ferromagnetic interaction ($-J_2^F<0$) produces a very rich quantum phase diagram. With frustration, the whole phase diagram is divided into a ferromagnetic (FM) and a nonmagnetic (NM) phase. For $s=1/2$, the full NM phase is seen to be of Haldane-Dimer type, but for $s>1/2$, a spiral phase comes between the FM and the Haldane-Dimer phases. The study of a suitably defined string-order parameter and spin-gap at the phase boundary indicates that the Haldane-Dimer and spiral phases have different topological characters. We also find that, along the $J_2^F=\frac 12 J_1^F$ line in the NM phase, an NN dimer state is the {\it exact} groundstate, provided $J_1^A>J_C=\kappa J_1^F$ where $\kappa \le s + h$ for applied magnetic field $h$. Without magnetic field, the position of $J_C$ is on the FM-NM phase boundary when $s=1/2$, but for $s>1/2$, the location of $J_C$ is on the phase separation line between the Haldane-Dimer and spiral phases.


title: ewN2HDECAY - A program for the Calculation of Electroweak One-Loop Corrections to Higgs Decays in the Next-to-Minimal Two-Higgs-Doublet Model Including State-of-the-Art QCD Corrections

abstract: We present in this paper our new program package ewN2HDECAY for the calculation of the partial decay widths and branching ratios of the Higgs bosons of the Next-to-Minimal 2-Higgs Doublet Model (N2HDM). The N2HDM is based on a general CP-conserving 2HDM which is extended by a real scalar singlet field. The program computes the complete electroweak one-loop corrections to all non-loop-induced two-body on-shell Higgs boson decays in the N2HDM and combines them with the state-of-the-art QCD corrections that are already implemented in the existing program N2HDECAY. Most of the independent input parameters of the electroweak sector of the N2HDM are renormalized in an on-shell scheme. The soft-$\mathbb{Z}_2$-breaking squared mass scale $m_{12}^2$ and the vacuum expectation value $v_S$ of the $SU(2)_L$ singlet field, however, are renormalized with $\overline{\text{MS}}$ conditions, while for the four scalar mixing angles $\alpha _i$ ($i=1,2,3$) and $\beta$ of the N2HDM, several different renormalization schemes are applied. By giving out the leading-order and the loop-corrected partial decay widths separately from the branching ratios, the program ewN2HDECAY not only allows for phenomenological analyses of the N2HDM at highest precision, it can also be used for a study of the impact of the electroweak corrections and the remaining theoretical uncertainty due to missing higher-order corrections based on a change of the renormalization scheme. The input parameters are then consistently calculated with a parameter conversion routine when switching from one renormalization scheme to the other. The latest version of the program ewN2HDECAY can be downloaded from the URL \href{https://github.com/marcel-krause/ewN2HDECAY}{https://github.com/marcel-krause/ewN2HDECAY}.


title: Target-Tailored Source-Transformation for Scene Graph Generation

abstract: Scene graph generation aims to provide a semantic and structural description of an image, denoting the objects (with nodes) and their relationships (with edges). The best performing works to date are based on exploiting the context surrounding objects or relations,e.g., by passing information among objects. In these approaches, to transform the representation of source objects is a critical process for extracting information for the use by target objects. In this work, we argue that a source object should give what tar-get object needs and give different objects different information rather than contributing common information to all targets. To achieve this goal, we propose a Target-TailoredSource-Transformation (TTST) method to efficiently propagate information among object proposals and relations. Particularly, for a source object proposal which will contribute information to other target objects, we transform the source object feature to the target object feature domain by simultaneously taking both the source and target into account. We further explore more powerful representations by integrating language prior with the visual context in the transformation for the scene graph generation. By doing so the target object is able to extract target-specific information from the source object and source relation accordingly to refine its representation. Our framework is validated on the Visual Genome bench-mark and demonstrated its state-of-the-art performance for the scene graph generation. The experimental results show that the performance of object detection and visual relation-ship detection are promoted mutually by our method.


title: Impact Dynamics of Moons Within a Planetary Potential

abstract: Current lunar origin scenarios suggest that Earth's Moon may have resulted from the merger of two (or more) smaller moonlets. Dynamical studies of multiple moons find that these satellite systems are not stable, resulting in moonlet collision or loss of one or more of the moonlets. We perform Smoothed Particle Hydrodynamic (SPH) impact simulations of two orbiting moonlets inside the planetary gravitational potential and find that the classical outcome of two bodies impacting in free space is altered as erosive mass loss is more significant with decreasing distance to the planet. Depending on the conditions of accretion, each moonlet could have a distinct isotopic signature, therefore, we assess the initial mixing during their merger, in order to estimate whether future measurements of surface variations could distinguish between lunar origin scenarios (single vs. multiple moonlets). We find that for comparable-size impacting bodies in the accretionary regime, surface mixing is efficient, but in the hit-and-run regime, only a small amount of material is transferred between the bodies. However, sequences of hit-and-run impacts are expected, which will enhance the surface mixing. Overall, our results show that large scale heterogeneities can arise only from the merger of drastically different component masses. Surfaces of moons resulting from the merger of comparable-sized components have little material heterogeneities, and such impacts are preferred, as the relatively massive impactor generates more melt, extending the lunar magma ocean phase.


title: Stratified Whitehead's theorem and knot invariants

abstract: By considering homotopies that preserve the stratification, one obtains a natural notion of homotopy for stratified spaces. In this short note, we introduce invariants of stratified homotopy, the stratified homotopy groups. We show that they satisify a stratified version of Whitehead's theorem. As an example, we introduce a complete knot invariant defined via the stratified homotopy groups ----- En consid\'erant des homotopies pr\'eservant la stratification, on obtient une notion naturelle d'homotopie pour les espaces stratifi\'es. Dans cette note, on pr\'esente des invariants d'homotopie stratifi\'ee, les groupes d'homotopie stratifi\'es. On montre que ces groupes d'homotopie stratifi\'es v\'erifient un analogue stratifi\'e au th\'eor\`eme de Whitehead. Comme illustration, on pr\'esente un invariant de noeud complet d\'efini \`a partir des groupes d'homotopie stratifi\'es.


title: Data-driven discovery of coordinates and governing equations

abstract: The discovery of governing equations from scientific data has the potential to transform data-rich fields that lack well-characterized quantitative descriptions. Advances in sparse regression are currently enabling the tractable identification of both the structure and parameters of a nonlinear dynamical system from data. The resulting models have the fewest terms necessary to describe the dynamics, balancing model complexity with descriptive ability, and thus promoting interpretability and generalizability. This provides an algorithmic approach to Occam's razor for model discovery. However, this approach fundamentally relies on an effective coordinate system in which the dynamics have a simple representation. In this work, we design a custom autoencoder to discover a coordinate transformation into a reduced space where the dynamics may be sparsely represented. Thus, we simultaneously learn the governing equations and the associated coordinate system. We demonstrate this approach on several example high-dimensional dynamical systems with low-dimensional behavior. The resulting modeling framework combines the strengths of deep neural networks for flexible representation and sparse identification of nonlinear dynamics (SINDy) for parsimonious models. It is the first method of its kind to place the discovery of coordinates and models on an equal footing.


title: Efficient LBM on GPUs for dense moving objects using immersed boundary condition

abstract: There exists an increasing interest for using immersed boundary methods (IBMs) (Peskin 2000) to model moving objects in computational fluid dynamics. Indeed, this approach is particularly efficient, because the fluid mesh does not require to be body-fitted or to adjust dynamically to the motion of the body. Frequently, IBMs are implemented in combination with the lattice Boltzmann methods (LBM) (Kr\"uger 2016). They fit elegantly into the framework of this method, and yield impressive parallel performances. It has also become quite common to accelerate LBM simulations with the use of Graphics Processing Units (GPUs) (T\"olke 2010), as the underlying algorithm adjusts naturally to the architecture of such platforms. It is not uncommon that speedups of an order of magnitude, or more, at equal financial cost or energy consumption are observed, as compared to classical CPUs. IBM algorithms are however more difficult to adapt to GPUs, because their complex memory access pattern conflicts with a GPU's strategy of broadcasting data to a large number of GPU cores in single memory accesses. In the existing literature, GPU implementations of LBM-IBM codes are therefore restricted to situations in which the immersed surfaces are very small compared to the total number of fluid cells (Valero-Lara 2014), as is often the case in exterior flow simulations around an obstacle. This assumption is however not valid in many other cases of interest. We propose a new method for the implementation of a LBM-IBM on GPUs in the CUDA language, which allows to handle a substantially larger immersed surfaces with acceptable performance than previous implementations.


title: Microscopyc description of the Photothermal increase of temperature for metallic nanoparticles excited with short-laser pulses

abstract: The pulsed photothermal phenomenon due to optically absorbed energy, result from non-radiative decay mechanisms, which in nature, these imply the temporal change in the local free energy and thus a temporary change in the local temperature. At the nanoscale, this is a prediction broadly described in terms of macroscopic variables. Here we introduce a formalism based on the Jarzynski's physical statistics description, for interpreting the equilibrium free energy difference between two configurations of a metallic nanoparticle. In this way, within a finite-time span, we describe the temporal increase of the local free energy and thus of the local temperature arising from temporarily bringing the sample far from thermal equilibrium. The result is an expression for which one can get the photothermally induced change of temperature for a metallic nanoparticle. For practical purposes, we limited the study to Au nanoparticles. The study is made as function of the particle size and the optical properties for wavelengths spanning the optical range. The assessment indicates that, for nanoparticles with radii shorter than 40 nm, the temperature change is strongly dependent on particle size and on the illumination wavelength. While, for near 40 nm particle radii, the current description and the known formalism, based on macroscopic variables, predict the very same temperature change. At the closing we discuss additional possible thermodynamic consequences associated to the scale considerations.


title: Inclined turbulent thermal convection in liquid sodium

abstract: Inclined turbulent thermal convection by large Rayleigh numbers in extremely small-Prandtl-number fluids is studied based on results of both, measurements and high-resolution numerical simulations. The Prandtl number $Pr\approx0.0093$ considered in the experiments and the Large-Eddy Simulations (LES) and $Pr=0.0094$ considered in the Direct Numerical Simulations (DNS) correspond to liquid sodium, which is used in the experiments. Also similar are the studied Rayleigh numbers, which are, respectively, $Ra=1.67\times10^7$ in the DNS, $Ra=1.5\times10^7$ in the LES and $Ra=1.42\times10^7$ in the measurements. The working convection cell is a cylinder with equal height and diameter, where one circular surface is heated and another one is cooled. The cylinder axis is inclined with respect to the vertical and the inclination angle varies from $\beta=0^\circ$, which corresponds to a Rayleigh-B\'enard configuration (RBC), to $\beta=90^\circ$, as in a vertical convection (VC) setup. The turbulent heat and momentum transport as well as time-averaged and instantaneous flow structures and their evolution in time are studied in detail, for different inclination angles, and are illustrated also by supplementary videos, obtained from the DNS and experimental data. To investigate the scaling relations of the mean heat and momentum transport in the limiting cases of RBC and VC configurations, additional measurements are conducted for about one decade of the Rayleigh numbers around $Ra=10^7$ and $Pr\approx0.009$. With respect to the turbulent heat transport in inclined thermal convection by low $Pr$, a similarity of the global flow characteristics for the same value of $RaPr$ is proposed and analysed, based on the above simulations and measurements and on complementary DNS for $Ra=1.67\times10^6$, $Pr=0.094$ and $Ra=10^9$, $Pr=1$.


title: Multidimensional Capacitive Sensing for Robot-Assisted Dressing and Bathing

abstract: Robotic assistance presents an opportunity to benefit the lives of many people with physical disabilities, yet accurately sensing the human body and tracking human motion remain difficult for robots. We present a multidimensional capacitive sensing technique that estimates the local pose of a human limb in real time. A key benefit of this sensing method is that it can sense the limb through opaque materials, including fabrics and wet cloth. Our method uses a multielectrode capacitive sensor mounted to a robot's end effector. A neural network model estimates the position of the closest point on a person's limb and the orientation of the limb's central axis relative to the sensor's frame of reference. These pose estimates enable the robot to move its end effector with respect to the limb using feedback control. We demonstrate that a PR2 robot can use this approach with a custom six electrode capacitive sensor to assist with two activities of daily living-dressing and bathing. The robot pulled the sleeve of a hospital gown onto able-bodied participants' right arms, while tracking human motion. When assisting with bathing, the robot moved a soft wet washcloth to follow the contours of able-bodied participants' limbs, cleaning their surfaces. Overall, we found that multidimensional capacitive sensing presents a promising approach for robots to sense and track the human body during assistive tasks that require physical human-robot interaction.


title: Measurement of the Gamma Ray Background in the Davis Cavern at the Sanford Underground Research Facility

abstract: Deep underground environments are ideal for low background searches due to the attenuation of cosmic rays by passage through the earth. However, they are affected by backgrounds from $\gamma$-rays emitted by $^{40}$K and the $^{238}$U and $^{232}$Th decay chains in the surrounding rock. The LUX-ZEPLIN (LZ) experiment will search for dark matter particle interactions with a liquid xenon TPC located within the Davis campus at the Sanford Underground Research Facility, Lead, South Dakota, at the 4,850-foot level. In order to characterise the cavern background, in-situ $\gamma$-ray measurements were taken with a sodium iodide detector in various locations and with lead shielding. The integral count rates (0--3300~keV) varied from 596~Hz to 1355~Hz for unshielded measurements, corresponding to a total flux in the cavern of $1.9\pm0.4$~$\gamma~$cm$^{-2}$s$^{-1}$. The resulting activity in the walls of the cavern can be characterised as $220\pm60$~Bq/kg of $^{40}$K, $29\pm15$~Bq/kg of $^{238}$U, and $13\pm3$~Bq/kg of $^{232}$Th.


title: Point Cloud Oversegmentation with Graph-Structured Deep Metric Learning

abstract: We propose a new supervized learning framework for oversegmenting 3D point clouds into superpoints. We cast this problem as learning deep embeddings of the local geometry and radiometry of 3D points, such that the border of objects presents high contrasts. The embeddings are computed using a lightweight neural network operating on the points' local neighborhood. Finally, we formulate point cloud oversegmentation as a graph partition problem with respect to the learned embeddings. This new approach allows us to set a new state-of-the-art in point cloud oversegmentation by a significant margin, on a dense indoor dataset (S3DIS) and a sparse outdoor one (vKITTI). Our best solution requires over five times fewer superpoints to reach similar performance than previously published methods on S3DIS. Furthermore, we show that our framework can be used to improve superpoint-based semantic segmentation algorithms, setting a new state-of-the-art for this task as well.


title: On Lie superalgebroids of Poisson type

abstract: The purpose of this work is to study Lie superalgebroid structures on the space of superdifferential $1$-forms over the supermanifolds whose superfunctions are the differential forms on its underlying manifold. These superalgbroids are constructed from graded Poisson structures defined on the latter superalgebra.


title: CMB Constraints on the Stochastic Gravitational-Wave Background at Mpc scales

abstract: We present robust constraints on the stochastic gravitational waves (GWs) at Mpc scales from the cosmic microwave background (CMB) data. CMB constraints on GWs are usually characterized as the tensor-to-scalar ratio, assuming specifically a power-law form of the primordial spectrum, and are obtained from the angular spectra of CMB. Here, we relax the assumption of the power-law form, and consider to what extent one can constrain a monochromatic GW at shorter wavelengths. Previously, such a constraint has been derived at the wavelengths larger than the resolution scale of the CMB measurements, typically above $100$Mpc (below $10^{-16}$Hz in frequency). However, GWs whose wavelength is much shorter than $100$Mpc can imprint a small but non-negligible signal on CMB anisotropies at observed angular scales, $\ell<1000$. Here, using the CMB temperature, polarization, and lensing data set, we obtain the best constraints to date at $10^{-16}-10^{-14}$Hz of the GWs produced before the time of decoupling, which are tighter than those derived from the astrometric measurements and upper bounds on extra radiations. In the future, the constraints on GWs at Mpc scales will be further improved by several orders of magnitude with the precision $B$-mode measurement on large scales, $\ell<100$.


title: Evidence for the Cross-correlation between Cosmic Microwave Background Polarization Lensing from POLARBEAR and Cosmic Shear from Subaru Hyper Suprime-Cam

abstract: We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.


title: Collective spinon spin wave in a magnetized U(1) spin liquid

abstract: We study the transverse dynamical spin susceptibility of the two dimensional U(1) spinon Fermi surface spin liquid in a small applied Zeeman field. We show that both short-range interactions, present in a generic Fermi liquid, as well as gauge fluctuations, characteristic of the U(1) spin liquid, qualitatively change the result based on the frequently assumed non-interacting spinon approximation. Short-range interaction leads to a new collective mode: a "spinon spin wave" which splits off from the two-spinon continuum at small momentum and disperses downward. Gauge fluctuations renormalize the susceptibility, providing non-zero power law weight in the region outside the spinon continuum and giving the spin wave a finite lifetime, which scales as momentum squared. We also study the effect of Dzyaloshinskii-Moriya anisotropy on the zero momentum susceptibility, which is measured in electron spin resonance (ESR), and obtain a resonance linewidth linear in temperature and varying as $B^{2/3}$ with magnetic field $B$ at low temperature. Our results form the basis for a theory of inelastic neutrons scattering, ESR, and resonant inelastic x-ray scattering (RIXS) studies of this quantum spin liquid state.


title: Unified Superradiant phase transitions

abstract: We prove, by means of a unified treatment, that the superradiant phase transitions of Dicke and classical oscillator limits of simple light-matter models are indeed of the same type. We show that the mean-field approximation is exact in both cases, and compute the structure and location of the transitions in parameter space. We extend this study to a fuller range of models, paying special attention to symmetry considerations. We uncover general features of the phase structure in the space of parameters of these models.


title: An Improved Upper Bound for the Ring Loading Problem

abstract: The Ring Loading Problem emerged in the 1990s to model an important special case of telecommunication networks (SONET rings) which gained attention from practitioners and theorists alike. Given an undirected cycle on $n$ nodes together with non-negative demands between any pair of nodes, the Ring Loading Problem asks for an unsplittable routing of the demands such that the maximum cumulated demand on any edge is minimized. Let $L$ be the value of such a solution. In the relaxed version of the problem, each demand can be split into two parts where the first part is routed clockwise while the second part is routed counter-clockwise. Denote with $L^*$ the maximum load of a minimum split routing solution. In a landmark paper, Schrijver, Seymour and Winkler [SSW98] showed that $L \leq L^* + 1.5D$, where $D$ is the maximum demand value. They also found (implicitly) an instance of the Ring Loading Problem with $L = L^* + 1.01D$. Recently, Skutella [Sku16] improved these bounds by showing that $L \leq L^* + \frac{19}{14}D$, and there exists an instance with $L = L^* + 1.1D$. We contribute to this line of research by showing that $L \leq L^* + 1.3D$. We also take a first step towards lower and upper bounds for small instances.


title: Astro2020 Science White Paper: Cold Gas Outflows, Feedback, and the Shaping of Galaxies

abstract: There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present observational and modeling challenges. Here we outline several key open questions, briefly discussing the requirements of the observations necessary to make progress, and the relevance of several existing and planned facilities. It is clear that galaxy outflows, and particularly cold outflows, will remain a topic of active research for the next decade and beyond.


title: Reversible Pebbling Game for Quantum Memory Management

abstract: Quantum memory management is becoming a pressing problem, especially given the recent research effort to develop new and more complex quantum algorithms. The only existing automatic method for quantum states clean-up relies on the availability of many extra resources. In this work, we propose an automatic tool for quantum memory management. We show how this problem exactly matches the reversible pebbling game. Based on that, we develop a SAT-based algorithm that returns a valid clean-up strategy, taking the limitations of the quantum hardware into account. The developed tool empowers the designer with the flexibility required to explore the trade-off between memory resources and number of operations. We present three show-cases to prove the validity of our approach. First, we apply the algorithm to straight-line programs, widely used in cryptographic applications. Second, we perform a comparison with the existing approach, showing an average improvement of 52.77%. Finally, we show the advantage of using the tool when synthesizing a quantum circuit on a constrained near-term quantum device.


title: Group-wise classification approach to improve Android malicious apps detection accuracy

abstract: In the fast-growing smart devices, Android is the most popular OS, and due to its attractive features, mobility, ease of use, these devices hold sensitive information such as personal data, browsing history, shopping history, financial details, etc. Therefore, any security gap in these devices means that the information stored or accessing the smart devices are at high risk of being breached by the malware. These malware are continuously growing and are also used for military espionage, disrupting the industry, power grids, etc. To detect these malware, traditional signature matching techniques are widely used. However, such strategies are not capable to detect the advanced Android malicious apps because malware developer uses several obfuscation techniques. Hence, researchers are continuously addressing the security issues in the Android based smart devices. Therefore, in this paper using Drebin benchmark malware dataset we experimentally demonstrate how to improve the detection accuracy by analyzing the apps after grouping the collected data based on the permissions and achieved 97.15% overall average accuracy. Our results outperform the accuracy obtained without grouping data (79.27%, 2017), Arp, et al. (94%, 2014), Annamalai et al. (84.29%, 2016), Bahman Rashidi et al. (82%, 2017)) and Ali Feizollah, et al. (95.5%, 2017). The analysis also shows that among the groups, Microphone group detection accuracy is least while Calendar group apps are detected with the highest accuracy, and with the highest accuracy, and for the best performance, one shall take 80-100 features.


title: Projective geometry of Wachspress coordinates

abstract: We show that there is a unique hypersurface of minimal degree passing through the non-faces of a polytope which is defined by a simple hyperplane arrangement. This generalizes the construction of the adjoint curve of a polygon by Wachspress in 1975. The defining polynomial of our adjoint hypersurface is the adjoint polynomial introduced by Warren in 1996. This is a key ingredient for the definition of Wachspress coordinates, which are barycentric coordinates on an arbitrary convex polytope. The adjoint polynomial also appears both in algebraic statistics, when studying the moments of uniform probability distributions on polytopes, and in intersection theory, when computing Segre classes of monomial schemes. We describe the Wachspress map, the rational map defined by the Wachspress coordinates, and the Wachspress variety, the image of this map. The inverse of the Wachspress map is the projection from the linear span of the image of the adjoint hypersurface. To relate adjoints of polytopes to classical adjoints of divisors in algebraic geometry, we study irreducible hypersurfaces that have the same degree and multiplicity along the non-faces of a polytope as its defining hyperplane arrangement. We list all finitely many combinatorial types of polytopes in dimensions two and three for which such irreducible hypersurfaces exist. In the case of polygons, the general such curves< are elliptic. In the three-dimensional case, the general such surfaces are either K3 or elliptic.


title: A Recoverable Mutex Algorithm with Sub-logarithmic RMR on Both CC and DSM

abstract: In light of recent advances in non-volatile main memory technology, Golab and Ramaraju reformulated the traditional mutex problem into the novel {\em Recoverable Mutual Exclusion} (RME) problem. In the best known solution for RME, due to Golab and Hendler from PODC 2017, a process incurs at most $O(\frac{\log n}{\log \log n})$ remote memory references (RMRs) per passage, where a passage is an interval from when a process enters the Try section to when it subsequently returns to Remainder. Their algorithm, however, guarantees this bound only for cache-coherent (CC) multiprocessors, leaving open the question of whether a similar bound is possible for distributed shared memory (DSM) multiprocessors. We answer this question affirmatively by designing an algorithm that satisfies the same complexity bound as Golab and Hendler's for both CC and DSM multiprocessors. Our algorithm has some additional advantages over Golab and Hendler's: (i) its Exit section is wait-free, (ii) it uses only the Fetch-and-Store instruction, and (iii) on a CC machine our algorithm needs each process to have a cache of only $O(1)$ words, while their algorithm needs $O(n)$ words.


title: The Kramers problem for SDEs driven by small, accelerated L\'evy noise with exponentially light jumps

abstract: We establish Freidlin-Wentzell results for a nonlinear ordinary differential equation starting close to the stable state $0$, say, subject to a perturbation by a stochastic integral which is driven by an $\varepsilon$-small and $(1/\varepsilon)$-accelerated L\'evy process with exponentially light jumps. For this purpose we derive a large deviations principle for the stochastically perturbed system using the weak convergence approach developed by Budhiraja, Dupuis, Maroulas and collaborators in recent years. In the sequel we solve the associated asymptotic first escape problem from the bounded neighborhood of $0$ in the limit as $\varepsilon \rightarrow 0$ which is also known as the Kramers problem in the literature.


title: On the origin of very massive stars around NGC 3603

abstract: The formation of the most massive stars in the Universe remains an unsolved problem. Are they able to form in relative isolation in a manner similar to the formation of solar-type stars, or do they necessarily require a clustered environment? In order to shed light on this important question, we study the origin of two very massive stars (VMS): the O2.5If*/WN6 star RFS7 ($\sim$100 $M_{\odot}$), and the O3.5If* star RFS8 ($\sim$70 $M_{\odot}$), found within $\approx$ 53 and 58 pc respectively from the Galactic massive young cluster NGC 3603, using Gaia data. RFS7 is found to exhibit motions resembling a runaway star from NGC 3603. This is now the most massive runaway star candidate known in the Milky Way. Although RFS8 also appears to move away from the cluster core, it has proper-motion values that appear inconsistent with being a runaway from NGC 3603 at the $3\sigma$ level (but with substantial uncertainties due to distance and age). Furthermore, no evidence for a bow-shock or a cluster was found surrounding RFS8 from available near-infrared photometry. In summary, whilst RFS7 is likely a runaway star from NGC 3603, making it the first VMS runaway in the Milky Way, RFS8 is an extremely young ($\sim$2 Myr) VMS, which might also be a runaway, but this would need to be established from future spectroscopic and astrometric observations, as well as precise distances. If RFS8 were still not meeting the criteria for being a runaway from NGC 3603 from such future data, this would have important ramifications for current theories of massive star formation, as well as the way the stellar initial mass function (IMF) is sampled.


title: One Solution to the Mass Budget Problem for Planet Formation: Optically Thick Disks with Dust Scattering

abstract: ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks.


title: Uniform limit of discrete convex functions

abstract: We consider mesh functions which are discrete convex in the sense that their central second order directional derivatives are positive. Analogous to the case of a uniformly bounded sequence of convex functions, we prove that the uniform limit on compact subsets of discrete convex mesh functions which are uniformly bounded is a continuous convex function. Furthermore, if the discrete convex mesh functions interpolate boundary data of a continuous convex function and their Monge-Ampere masses are uniformly bounded, the limit function satisfies the boundary condition strongly. The domain of the solution needs not be uniformly convex. The result is applied to the convergence of some numerical methods for the Monge-Ampere equation.


title: Primordial black holes from the QCD epoch: Linking dark matter, baryogenesis and anthropic selection

abstract: If primordial black holes (PBHs) formed at the quark-hadron epoch, their mass must be close to the Chandrasekhar limit, this also being the characteristic mass of stars. If they provide the dark matter (DM), the collapse fraction must be of order the cosmological baryon-to-photon ratio $\sim 10^{-9}$, which suggests a scenario in which a baryon asymmetry is produced efficiently in the outgoing shock around each PBH and then propagates to the rest of the Universe. We suggest that the temperature increase in the shock provides the ingredients for hot spot electroweak baryogenesis. This also explains why baryons and DM have comparable densities, the precise ratio depending on the size of the PBH relative to the cosmological horizon at formation. The observed value of the collapse fraction and baryon asymmetry depends on the amplitude of the curvature fluctuations which generate the PBHs and may be explained by an anthropic selection effect associated with the existence of galaxies. We propose a scenario in which the quantum fluctuations of a light stochastic spectator field during inflation generate large curvature fluctuations in some regions, with the stochasticity of this field providing the basis for the required selection. Finally, we identify several observational predictions of our scenario that should be testable within the next few years. In particular, the PBH mass function could extend to sufficiently high masses to explain the black hole coalescences observed by LIGO/Virgo.


title: Normal Approximation for Stochastic Gradient Descent via Non-Asymptotic Rates of Martingale CLT

abstract: We provide non-asymptotic convergence rates of the Polyak-Ruppert averaged stochastic gradient descent (SGD) to a normal random vector for a class of twice-differentiable test functions. A crucial intermediate step is proving a non-asymptotic martingale central limit theorem (CLT), i.e., establishing the rates of convergence of a multivariate martingale difference sequence to a normal random vector, which might be of independent interest. We obtain the explicit rates for the multivariate martingale CLT using a combination of Stein's method and Lindeberg's argument, which is then used in conjunction with a non-asymptotic analysis of averaged SGD proposed in [PJ92]. Our results have potentially interesting consequences for computing confidence intervals for parameter estimation with SGD and constructing hypothesis tests with SGD that are valid in a non-asymptotic sense.


title: Private Information Retrieval from Heterogeneous Uncoded Storage Constrained Databases with Reduced Sub-Messages

abstract: We propose capacity-achieving schemes for private information retrieval (PIR) from uncoded databases (DBs) with both homogeneous and heterogeneous storage constraints. In the PIR setting, a user queries a set of DBs to privately download a message, where privacy implies that no one DB can infer which message the user desires. In general, a PIR scheme is comprised of storage placement and delivery designs. Previous works have derived the capacity, or infimum download cost, of PIR with uncoded storage placement and also sufficient conditions of a storage placement design to meet capacity. However, the currently proposed storage placement designs require splitting each message into an exponential number of sub-messages with respect to the number of DBs. In this work, when DBs have the same storage constraint, we propose two simple storage placement designs that satisfy the capacity conditions. Then, for more general heterogeneous storage constraints, we translate the storage placement design process into a "filling problem". We design an iterative algorithm to solve the filling problem where, in each iteration, messages are partitioned into sub-messages and stored at subsets of DBs. All of our proposed storage placement designs require a number of sub-messages per message at most equal to the number of DBs.


title: Large deviations in models of growing clusters with symmetry-breaking transitions

abstract: We analyse large deviations of the magnetisation in two models of growing clusters. The models have symmetry-breaking transitions, so the typical magnetisation of a growing cluster may be either positive or negative, with equal probability. For large clusters, the magnetisation obeys a large deviation principle. We show that the corresponding rate function is zero for values of the magnetisation that are intermediate between the two steady state values, which means that fluctuations with these values of the magnetisation are much less unlikely than previously thought. We show that their probabilities decay as power laws in the cluster size, instead of the exponential scaling that would be expected from the large deviation principle. We discuss how this observation is related to dynamical phase coexistence phenomena. We also comment on the typical size of magnetisation fluctuations.


title: Local reactive boundary scheme for lattice Boltzmann method

abstract: In this paper, a boundary scheme is proposed for the two-dimensional five-velocity (D2Q5) lattice Boltzmann method with heterogeneous surface reaction, in which the unknown distribution function is determined locally based on the kinetic flux of the incident particles. Compared with previous boundary schemes, the proposed scheme has a clear physical picture that reflects the consumption and production in the reaction. Furthermore, the scheme only involves local information of boundary nodes such that it can be easily applied to complex geometric structures. In order to validate the accuracy of the scheme, some benchmark tests, including the convection-diffusion problems in straight and inclined channels are conducted. Numerical results are in excellent agreement with the analytical solutions, and the convergence tests demonstrate that second-order spatial accuracy is achieved for straight walls, and the order of accuracy is between 1.5 and 2.0 for general inclined walls. Finally, we simulated the density driving flow with dissolution reactions in a two-dimensional cylindrical array, and the results agree well with those in previous studies


title: A Truncation Error Estimation Scheme for the Finite Volume Method on Unstructured Meshes

abstract: This work is an attempt to develop an approximate scheme for estimating the volume-based truncation errors in the finite volume analysis of laminar flows. The volume-based truncation error is the net flow error across the faces of a control volume. Unfortunately, truncation error is not a natural outcome of the finite volume solution and needs to be estimated separately. Previous works in the literature estimate truncation error using either higher order interpolation schemes, higher order discretization schemes, or neglected terms in the discretization scheme. The first two approaches become complicated on general unstructured meshes and the third approach provides inaccurate results. This work proposes a truncation error estimation scheme, which is based on the third approach, but provides more accurate results compared to the existing results in the literature. The potential application of such a truncation error estimation scheme is in mesh adaptation.


title: Towards sodium combustion modelling with liquid water

abstract: Solid and liquid sodium combustion with liquid water occurs through a thin gas layer where exothermic reactions happen with sodium and water vapors. It thus involves multiple interfaces separating liquid and gas in the presence of surface tension, phase transition and surface reactions. The gas phase reaction involves compressible effects resulting in possible shock wave appearance in both gas and liquid phases. To understand and predict the complexity of sodium combustion with water a diffuse interface flow model is built. This formulation enables flow resolution in multidimension in the presence of complex motion, such as for example Leidenfrost-type thermo-chemical flow. More precisely sodium drop autonomous motion on the liquid surface is computed. Various modelling and numerical issues are present and addressed in the present contribution. In the author's knowledge, the first computed results of such type of combustion phenomenon in multidimensions are presented in this paper thanks to the diffuse interface approach. Explosion phenomenon is addressed as well and is reproduced at least qualitatively thanks to extra ingredients such as turbulent mixing of sodium and water vapors in the gas film and delayed ignition. Shock wave emission from the thermo-chemical Leidenfrost-type flow is observed as reported in related experiments.


title: Automatic Differentiation using Operator Overloading (ADOO) for implicit resolution of hyperbolic single phase and two-phase flow models

abstract: Implicit time integration schemes are widely used in computational fluid dynamics numerical codes to speed-up computations. Indeed, implicit schemes usually allow for less stringent time-step stability constraints than their explicit counterpart. The derivation of an implicit scheme is however a challenging and time-consuming task, increasing substantially with the model equations complexity since this method usually requires a fairly accurate evaluation of the spatial scheme's matrix Jacobian. This article presents a flexible method to overcome the difficulties associated to the computation of the derivatives, based on the forward mode of automatic differentiation using operator overloading. Flexibility and simplicity of the method are illustrated through implicit resolution of various flow models of increasing complexity such as the compressible Euler equations, a two-phase flow model in full equilibrium (Le Martelot, et al., 2014) and a symmetric variant (Saurel , et al., 2003) of the two-phase flow model of (Baer \& Nunziato, 1986) dealing with mixtures in total disequilibrium.


title: Supervised parallel-in-time algorithm for long-time Lagrangian simulations of stochastic dynamics: Application to hydrodynamics

abstract: Lagrangian particle methods based on detailed atomic and molecular models are powerful computational tools for studying the dynamics of microscale and nanoscale systems. However, the maximum time step is limited by the smallest oscillation period of the fastest atomic motion, rendering long-time simulations very expensive. To resolve this bottleneck, we propose a supervised parallel-in-time algorithm for stochastic dynamics (SPASD) to accelerate long-time Lagrangian particle simulations. Our method is inspired by bottom-up coarse-graining projections that yield mean-field hydrodynamic behavior in the continuum limit. Here as an example, we use the dissipative particle dynamics (DPD) as the Lagrangian particle simulator that is supervised by its macroscopic counterpart, i.e., the Navier-Stokes simulator. The low-dimensional macroscopic system (here, the Navier-Stokes solver) serves as a predictor to supervise the high-dimensional Lagrangian simulator, in a predictor-corrector type algorithm. The results of the Lagrangian simulation then correct the mean-field prediction and provide the proper microscopic details (e.g., consistent fluctuations, correlations, etc.). The unique feature that sets SPASD apart from other multiscale methods is the use of a low-fidelity macroscopic model as a predictor. The macro-model can be approximate and even inconsistent with the microscale description, but SPASD anticipates the deviation and corrects it internally to recover the true dynamics.


title: A coupled finite-volume solver for numerical simulation of electrically-driven flows

abstract: The accuracy and stability of implicit CFD codes are frequently impaired by the decoupling between variables, which can ultimately lead to numerical divergence. Coupled solvers, which solve all the governing equations simultaneously, have the potential to fix this problem. In this work, we report the implementation of coupled solvers for transient and steady-state electrically-driven flow simulations in the finite-volumes framework. The numerical method, developed in OpenFOAM, is generic for Newtonian and viscoelastic fluids and is formulated for the Poisson-Nernst-Planck and Poisson-Boltzmann models. The resulting coupled systems of equations are solved efficiently with PETSc library. The performance of the coupled solvers is assessed in two test cases: induced-charge electroosmosis of a Newtonian fluid around a cylinder; electroosmotic flow of a PTT viscoelastic fluid in a contraction/expansion microchannel. The coupled solvers are more accurate in transient simulations and allow the use of larger time-steps without numerical divergence. For steady-state simulations, the coupled solvers converge in fewer iterations than segregated solvers. Although coupled solvers are much slower in a per time-step basis, the overall speedup factor obtained in this study reached a maximum value of ~100, where the highest factors have been obtained with semi-coupled solvers, which drop some coupling terms between equations. While further research is needed to improve the efficiency of the matrix solving stage, coupled solvers are already superior to segregated solvers in a number of cases.


title: Actinide-rich and Actinide-poor $r$-Process Enhanced Metal-Poor Stars do not Require Separate $r$-Process Progenitors

abstract: The astrophysical production site of the heaviest elements in the universe remains a mystery. Incorporating heavy element signatures of metal-poor, $r$-process enhanced stars into theoretical studies of $r$-process production can offer crucial constraints on the origin of heavy elements. In this study, we introduce and apply the "Actinide-Dilution with Matching" model to a variety of stellar groups ranging from actinide-deficient to actinide-enhanced to empirically characterize $r$-process ejecta mass as a function of electron fraction. We find that actinide-boost stars do not indicate the need for a unique and separate $r$-process progenitor. Rather, small variations of neutron richness within the same type of $r$-process event can account for all observed levels of actinide enhancements. The very low-$Y_e$, fission-cycling ejecta of an $r$-process event need only constitute 10-30% of the total ejecta mass to accommodate most actinide abundances of metal-poor stars. We find that our empirical $Y_e$ distributions of ejecta are similar to those inferred from studies of GW170817 mass ejecta ratios, which is consistent with neutron-star mergers being a source of the heavy elements in metal-poor, $r$-process enhanced stars.


title: Exact Correlation Functions for Dual-Unitary Lattice Models in 1+1 Dimensions

abstract: We consider a class of quantum lattice models in $1+1$ dimensions represented as local quantum circuits that enjoy a particular "dual-unitarity" property. In essence, this property ensures that both the evolution "in time" and that "in space" are given in terms of unitary transfer matrices. We show that for this class of circuits, generically non-integrable, one can compute explicitly all dynamical correlations of local observables. Our result is exact, non-pertubative, and holds for any dimension $d$ of the local Hilbert space. In the minimal case of qubits ($d = 2$) we also present a classification of all dual-unitary circuits which allows us to single out a number of distinct classes for the behaviour of the dynamical correlations. We find "non-interacting" classes, where all correlations are preserved, the ergodic and mixing one, where all correlations decay, and, interestingly, also classes that are are both interacting and non-ergodic.

categories: cond-mat.stat-mech hep-th nlin.CD quant-ph


title: CAN-NER: Convolutional Attention Network for Chinese Named Entity Recognition

abstract: Named entity recognition (NER) in Chinese is essential but difficult because of the lack of natural delimiters. Therefore, Chinese Word Segmentation (CWS) is usually considered as the first step for Chinese NER. However, models based on word-level embeddings and lexicon features often suffer from segmentation errors and out-of-vocabulary (OOV) words. In this paper, we investigate a Convolutional Attention Network called CAN for Chinese NER, which consists of a character-based convolutional neural network (CNN) with local-attention layer and a gated recurrent unit (GRU) with global self-attention layer to capture the information from adjacent characters and sentence contexts. Also, compared to other models, not depending on any external resources like lexicons and employing small size of char embeddings make our model more practical. Extensive experimental results show that our approach outperforms state-of-the-art methods without word embedding and external lexicon resources on different domain datasets including Weibo, MSRA and Chinese Resume NER dataset.


title: Unsupervised Latent Tree Induction with Deep Inside-Outside Recursive Autoencoders

abstract: We introduce deep inside-outside recursive autoencoders (DIORA), a fully-unsupervised method for discovering syntax that simultaneously learns representations for constituents within the induced tree. Our approach predicts each word in an input sentence conditioned on the rest of the sentence and uses inside-outside dynamic programming to consider all possible binary trees over the sentence. At test time the CKY algorithm extracts the highest scoring parse. DIORA achieves a new state-of-the-art F1 in unsupervised binary constituency parsing (unlabeled) in two benchmark datasets, WSJ and MultiNLI.


title: Liouville type theorems and regularity of solutions to degenerate or singular problems part I: even solutions

abstract: We consider a class of equations in divergence form with a singular/degenerate weight $$-\mathrm{div}(|y|^a A(x,y)\nabla u)=|y|^a f(x,y)\; \quad\textrm{or} \ \textrm{div}(|y|^aF(x,y))\;.$$ Under suitable regularity assumptions for the matrix $A$ and $f$ (resp. $F$) we prove H\"older continuity of solutions which are even in $y\in\mathbb{R}$, and possibly of their derivatives up to order two or more (Schauder estimates). In addition, we show stability of the $C^{0,\alpha}$ and $C^{1,\alpha}$ a priori bounds for approximating problems in the form $$-\mathrm{div}((\varepsilon^2+y^2)^{a/2} A(x,y)\nabla u)=(\varepsilon^2+y^2)^{a/2} f(x,y)\; \quad\textrm{or} \ \textrm{div}((\varepsilon^2+y^2)^{a/2}F(x,y))$$ as $\varepsilon\to 0$. Finally, we derive $C^{0,\alpha}$ and $C^{1,\alpha}$ bounds for inhomogenous Neumann boundary problems as well. Our method is based upon blow-up and appropriate Liouville type theorems.


title: HopSkipJumpAttack: A Query-Efficient Decision-Based Attack

abstract: The goal of a decision-based adversarial attack on a trained model is to generate adversarial examples based solely on observing output labels returned by the targeted model. We develop HopSkipJumpAttack, a family of algorithms based on a novel estimate of the gradient direction using binary information at the decision boundary. The proposed family includes both untargeted and targeted attacks optimized for $\ell_2$ and $\ell_\infty$ similarity metrics respectively. Theoretical analysis is provided for the proposed algorithms and the gradient direction estimate. Experiments show HopSkipJumpAttack requires significantly fewer model queries than Boundary Attack. It also achieves competitive performance in attacking several widely-used defense mechanisms. (HopSkipJumpAttack was named Boundary Attack++ in a previous version of the preprint.)


title: Fitting stochastic predator-prey models using both population density and kill rate data

abstract: Most mechanistic predator-prey modelling has involved either parameterization from process rate data or inverse modelling. Here, we take a median road: we aim at identifying the potential benefits of combining datasets, when both population growth and predation processes are viewed as stochastic. We fit a discrete-time, stochastic predator-prey model of the Leslie type to simulated time series of densities and kill rate data. Our model has both environmental stochasticity in the growth rates and interaction stochasticity, i.e., a stochastic functional response. We examine what the kill rate data brings to the quality of the estimates, and whether estimation is possible (for various time series lengths) solely with time series of population counts or biomass data. Both Bayesian and frequentist estimation are performed, providing multiple ways to check model identifiability. The Fisher Information Matrix suggests that models with and without kill rate data are all identifiable, although correlations remain between parameters that belong to the same functional form. However, our results show that if the attractor is a fixed point in the absence of stochasticity, identifying parameters in practice requires kill rate data as a complement to the time series of population densities, due to the relatively flat likelihood. Only noisy limit cycle attractors can be identified directly from population count data (as in inverse modelling), although even in this case, adding kill rate data - including in small amounts - can make the estimates much more precise. Overall, we show that under process stochasticity in interaction rates, interaction data might be essential to obtain identifiable dynamical models for multiple species. These results may extend to other biotic interactions than predation, for which similar models combining interaction rates and population counts could be developed.


title: Purveyors of fine halos: Re-assessing globular cluster contributions to the Milky Way halo build-up with SDSS-IV

abstract: There is ample evidence in the Milky Way for globular cluster (GC) disruption. Hence one may expect that also part of the Galactic halo field stars may once have formed in GCs. We quantify the fraction of halo stars donated by GCs by searching for stars that bear the unique chemical fingerprints typical for a subset of GC stars often dubbed `second-generation stars'. These are stars showing light element abundance anomalies such as a pronounced CN-band strength accompanied by weak CH-bands. Based on this indicator, past studies have placed the fraction of halo stars with a GC origin between a few to up to 50%. Using low-resolution spectra from the most recent data release of the latest extension of the Sloan Digital Sky Survey (SDSS-IV), we were able to identify 118 metal-poor ($-1.8\le$[Fe/H]$\le -1.3$) CN-strong stars in a sample of 4470 halo giant stars out to 50 kpc. This results in an observed fraction of these stars of 2.6$\pm$0.2%. Using an updated formalism to account for the fraction of stars lost early on in the GCs' evolution we estimate the fraction of the halo that stems from disrupted clusters to be 11$\pm$1%. This number represents the case that stars lost from GCs were entirely from the first generation and is thus merely an upper limit. Our conclusions are sensitive to our assumptions of the mass lost early on from the first generation formed in the GCs, the ratio of first-to-second generation stars, and other GC parameters. We carefully test the influence of varying these parameters on the final result and find that, under realistic scenarios, the above fraction depends on the main assumptions at less than 10%. We further recover a flat trend in this fraction with Galactocentric radius, with a marginal indication of a rise beyond 30 kpc that could reflect the ex-situ origin of the outer halo. (abridged)


title: Learning Shared Encoding Representation for End-to-End Speech Recognition Models

abstract: In this work, we learn a shared encoding representation for a multi-task neural network model optimized with connectionist temporal classification (CTC) and conventional framewise cross-entropy training criteria. Our experiments show that the multi-task training not only tackles the complexity of optimizing CTC models such as acoustic-to-word but also results in significant improvement compared to the plain-task training with an optimal setup. Furthermore, we propose to use the encoding representation learned by the multi-task network to initialize the encoder of attention-based models. Thereby, we train a deep attention-based end-to-end model with 10 long short-term memory (LSTM) layers of encoder which produces 12.2\% and 22.6\% word-error-rate on Switchboard and CallHome subsets of the Hub5 2000 evaluation.


title: TLS 1.3 for engineers: An exploration of the TLS 1.3 specification and OpenJDK's Java implementation

abstract: The Internet delivered in excess of forty terabytes per second in 2017 (Cisco, 2018), and over half of today's Internet traffic is encrypted (Sandvine, 2018); enabling trade worth trillions of dollars (Statista, 2017). Yet, the underlying encryption technology is only understood by a select few. This manuscript broadens understanding by exploring TLS, an encryption technology used to protect application layer communication (including HTTP, FTP and SMTP traffic), and by examining Oracle's Java implementation. We focus on the most recent TLS release, namely, version 1.3, which is defined by RFC 8446.


title: Ionization driven intrinsic absorption line variability of BAL quasars in the Stripe 82 region

abstract: We investigate the connection between the intrinsic C IV absorption line variability and the continuum flux changes of broad absorption line (BAL) quasars using a sample of 78 sources in the Stripe 82 region. The absorption trough variability parameters are measured using the archival multi-epoch spectroscopic data from the Sloan Digital Sky Survey (SDSS), and the continuum flux variability parameters are estimated from the photometric light curves obtained by the SDSS and the Catalina Real-Time Survey (CRTS) surveys. We find evidence for weak correlations (\rho ~ 0.3) between the intrinsic C IV absorption line variability and the quasar continuum variability for the final sample of 78 BAL quasars. The correlation strengths improve (\rho ~0.5) for the "high-SNR" sample sources that have higher spectral signal-to-noise ratio. Using two sub-sets of the high-SNR sample differing on the absorption trough depth, we find that the shallow trough sub-set shows an even stronger correlation (\rho ~ 0.6), whereas the deep trough sub-set does not show any correlation between the absorption line variability and the continuum variability. These results point to the important role of saturation effects in the correlation between the absorption line variability and the continuum variability of BAL quasars. Considering other effects that can also smear the correlation, we conclude that the actual correlation between the absorption line and continuum variability is even stronger.


title: Two Peculiar Classes of Solvable Systems Featuring 2 Dependent Variables Evolving in Discrete-Time via 2 Nonlinearly-Coupled First-Order Recursion Relations

abstract: In this paper we identify certain peculiar systems of 2 discrete-time evolution equations,x~n = F^(n)(x1; x2) , n = 1, 2 , which are algebraically solvable. Here l is the "discrete-time" independent variable taking integer values (l = 0, 1, 2,...), xn = xn (l) are 2 dependent variables, and x~n = xn (l + 1) are the corresponding 2 updated variables. In a previous paper the 2 functions F^(n)(x1, x2), n = 1, 2 were defined as follows: F^(n)(x1; x2) = P2 (xn, xn+1), n = 1,2 mod[2]; with P2 (x1; x2) a specific second-degree homogeneous polynomials in the 2 (indistinguishable!) dependent variables x1 (`) and x2 (l). In the present paper we further clarify some aspects of that model and we present its extension to the case when F(n)(x1, x2) = Q^(n)_k(x1, x2), n = 1, 2 mod[2], with Q(n)_k(x1; x2) a specific homogeneous function of arbitrary (integer ) degree k (hence a polynomial of degree k when k > 0) in the 2 dependent variables x1 (l) and x2 (l).


title: Solvable Dynamical Systems in the Plane with Polynomial Interactions

abstract: In this paper we report a few examples of algebraically solvable dynamical systems characterized by 2 coupled Ordinary Differential Equations which read as follows: x_n = P(n) (x1, x2) , n = 1, 2 , with P(n) (x1, x2) specific polynomials of relatively low degree in the 2 dependent variables x1 = x1 (t) and x2 = x2 (t) . These findings are obtained via a new twist of a recent technique to identify dynamical systems solvable by algebraic operations, themselves explicitly identified as corresponding to the time evolutions of the zeros of polynomials the coefficients of which evolve according to algebraically solvable (systems of) evolution equations.


title: Some Algebraically Solvable Two-Dimensional Dynamical Systems with Polynomial Interactions

abstract: We tersely review a recently introduced technique to identify systems of two nonlinearly-coupled Ordinary Di{\S}erential Equations (ODEs) solvable by algebraic operations; and we report some specifc examples of this kind, namely systems of 2 first-order ODEs with polynomial right-hand sides, x_ n= P(n)(x1, x2) , n = 1, 2 , satisfied by the 2 (possibly complex ) dependent variables xn = xn (t). Here P(n)(x1, x2) indicates some specific polynomial. These examples are analogous, but different, from those previously reported.


title: Quantum Double Models coupled to matter fields: a detailed review for a dualization procedure

abstract: In this paper, we investigate how it is possible to define a new class of lattice gauge models based on a dualization procedure of a previous generalization of the Kitaev Quantum Double Models. In the case of this previous generalization that will be used as a basis, it was defined by adding new qudits (which were denoted as matter fields in reference to some works) to the lattice vertices with the intention of, for instance, interpreting its models as Kitaev Quantum Double Models coupled with Potts ones. Now, with regard to the generalization that we investigate here, which we want to define as the dual of this previous one, these new qudits were added to the lattice faces. And as the coupling between gauge and matter qudits of the previous generalization was performed by a gauge group action, we show that the dual behaviour of these two generalizations was achieved by coupling these same qudits in the second one through a gauge group co-action homomorphism. One of the most striking dual aspects of these two generalizations is that, in both, part of the quasiparticles that were inherited from the Kitaev Quantum Double Models become confined when these action and co-action are non-trivial. But the big news here is that, in addition to the group homomorphism (that defines this gauge group co-action) allows us to classify all the different models of this second generalization, this same group homomorphism also suggests that all these models can be interpreted as two-dimensional restrictions of the $ 2 $-lattice gauge theories.

categories: quant-ph cond-mat.str-el hep-th math-ph math.MP


title: Extension of the Hertz model for accounting to surface tension in nano-indentation tests of soft materials

abstract: The contact between a spherical indenter and a solid is considered. A numerical finite element model (F. E. M) to taking into account the surface tension of the solid is presented and assessed. It is shown that for nano-indentation of soft materials, the surface tension of the solid influences significantly the reaction force due to indentation. The validity of the classical Hertz model is defined. In very good approximation, the force vs. indentation depth curve can be fitted by a power law function $F=a^\delta b$ where $F$ denotes the force acting on the indentor, $d$ the indentation depth, $a$ and $b\in\ ]1,1.5]$ are constants depending on the materials and the size of the indentor.


title: Indentation, elasticity and surface tension

abstract: The classical models of Hertz, Sneddon and Boussinesq provide solutions for problems of indentation of a semi-infinite elastic massif by a sphere, a sphere or a cone and a flat punch. Although these models have been widely tested, it appears that at small scales and for flexible materials, surface tension can contribute to considerably to the mechanical response to indentation. The scales are typically those of the less than one micron for an elastomer and less than one millimetre for a gel. The exploitation of certain experimental results of microscopy or nanoindentation remain approximate due to the absence of models incorporating the effect of surface tension.


title: Violation of CHSH inequality and marginal laws in mixed sequential measurements with order effects

abstract: We model a typical Bell-test experimental situation by considering that Alice and Bob perform incompatible measurements in a sequential way, with mixed orders of execution. After emphasizing that order effects will generally produce a violation of the marginal laws, we derive an upper limit for the observed correlations. More precisely, when Alice's and Bob's measurements are compatible, the marginal laws are obeyed and Tsirelson's bound limits the quantum correlations in the Bell-CHSH inequality to $2\sqrt{2}$. On the other hand, when Alice and Bob perform incompatible mixed sequential measurements, the marginal laws are typically violated and the upper limit for the correlations is pushed up to $2\sqrt{3}$. Considering that significant violations of the marginal laws (also called no-signaling conditions) have been observed in the data of numerous Bell-test experiments, the present analysis provides a possible mechanism for their appearance, when the protocols are such that Alice's and Bob's measurements can be assumed to be performed in a mixed sequential way. We however emphasize that this does not imply that a communication with superluminal effective speed would be possible.


title: Observation of a thermoelectric Hall plateau in the extreme quantum limit

abstract: The thermoelectric Hall effect is the generation of a transverse heat current upon applying an electric field in the presence of a magnetic field. Here we demonstrate that the thermoelectric Hall conductivity $\alpha_{xy}$ in the three-dimensional Dirac semimetal ZrTe$_5$ acquires a robust plateau in the extreme quantum limit of magnetic field. The plateau value is independent of the field strength, disorder strength, carrier concentration, or carrier sign. We explain this plateau theoretically and show that it is a unique signature of three-dimensional Dirac or Weyl electrons in the extreme quantum limit. We further find that other thermoelectric coefficients, such as the thermopower and Nernst coefficient, are greatly enhanced over their zero-field values even at relatively low fields.


title: Probing 3D Structure with a Large MUSE Mosaic: Extending the Mass Model of Frontier Field Abell 370

abstract: We present an updated strong-lensing analysis of the massive cluster Abell 370 (A370), continuing the work first presented in Lagattuta et al. (2017). In this new analysis, we take advantage of the deeper imaging data from the Hubble Space Telescope (HST) Frontier Fields program, as well as a large spectroscopic mosaic obtained with the Multi-Unit Spectroscopic Explorer (MUSE). Thanks to the extended coverage of this mosaic, we probe the full 3D distribution of galaxies in the field, giving us a unique picture of the extended structure of the cluster and its surroundings. Our final catalog contains 584 redshifts, representing the largest spectroscopic catalog of A370 to date. Constructing the model, we measure a total mass distribution that is quantitatively similar to our previous work -- though to ensure a low rms error in the model fit, we invoke a significantly large external shear term. Using the redshift catalog, we search for other bound groups of galaxies, which may give rise to a more physical interpretation of this shear. We identify three structures in narrow redshift ranges along the line of sight, highlighting possible infalling substructures into the main cluster halo. We also discover additional substructure candidates in low-resolution imaging at larger projected radii. More spectroscopic coverage of these regions (pushing close to the A370 virial radius) and more extended, high-resolution imaging will be required to investigate this possibility, further advancing the analysis of these interesting developments.


title: Stars that Move Together Were Born Together

abstract: It is challenging to reliably identify stars that were born together outside of actively star-forming regions and bound stellar systems. However, co-natal stars should be present throughout the Galaxy, and their demographics can shed light on the clustered nature of star formation and the dynamical state of the disk. In previous work we presented a set of simulations of the Galactic disk that followed the clustered formation and dynamical evolution of 4 billion individual stars over the last 5 Gyr. The simulations predict that a high fraction of co-moving stars with physical and 3D velocity separation of $\Delta r < 20$ pc and $\Delta v < 1.5$ km s$^{-1}$ are co-natal. In this \textit{Letter}, we use \textit{Gaia} DR2 and LAMOST DR4 data to identify and study co-moving pairs. We find that the distribution of relative velocities and separations of pairs in the data is in good agreement with the predictions from the simulation. We identify 111 co-moving pairs in the Solar neighborhood with reliable astrometric and spectroscopic measurements. These pairs show a strong preference for having similar metallicities when compared to random field pairs. We therefore conclude that these pairs were very likely born together. The simulations predict that co-natal pairs originate preferentially from high-mass and relatively young ($< 1$ Gyr) star clusters. \textit{Gaia} will eventually deliver well-determined metallicities for the brightest stars, enabling the identification of thousands of co-natal pairs due to disrupting star clusters in the solar neighborhood.


title: Non-Hermitian topology of spontaneous magnon decay

abstract: Spontaneous magnon decay is a generic feature of the magnetic excitations of anisotropic magnets and isotropic magnets with non-collinear order. In this paper, we argue that the effect of interactions on one-magnon states can, under many circumstances, be treated in terms of an effective, energy independent, non-Hermitian Hamiltonian for the magnons. In the vicinity of Dirac or Weyl touching points, we show that the spectral function has a characteristic anisotropy arising from topologically protected exceptional points or lines in the non-Hermitian spectrum. Such features can, in principle, be detected using inelastic neutron scattering or other spectroscopic probes. We illustrate this physics through a concrete example: a honeycomb ferromagnet with Dzyaloshinskii-Moriya exchange. We perform interacting spin wave calculations of the structure factor and spectral function of this model, showing good agreement with results from a simple effective non-Hermitian model for the splitting of the Dirac point. Finally, we argue that the zoo of known topological protected magnon band structures may serve as a nearly ideal platform for realizing and exploring non-Hermitian physics in solid-state systems.


title: Envelope Convection, Surface Magnetism and Spots in A- and late B-type Stars

abstract: Weak magnetic fields have recently been detected in a number of A-type stars, including Vega and Sirius. At the same time, space photometry observations of A- and late B-type stars from Kepler and TESS have highlighted the existence of rotational modulation of surface features akin to stellar spots. Here we explore the possibility that surface magnetic spots might be caused by the presence of small envelope convective layers at or just below the stellar surface, caused by recombination of H and He. Using 1D stellar evolution calculations and assuming an equipartition dynamo, we make simple estimates of field strength at the photosphere. For most models the largest effects are caused by a convective layer driven by second helium ionization. While it is difficult to predict the geometry of the magnetic field, we conclude that the majority of intermediate-mass stars should have dynamo-generated magnetic fields of order a few gauss at the surface. These magnetic fields can appear at the surface as bright spots, and cause photometric variability via rotational modulation, which could also be wide-spread in A-stars. The amplitude of surface magnetic fields and their associated photometric variability is expected to decrease with increasing stellar mass and surface temperature, so that magnetic spots and their observational effects should be much harder to detect in late B-type stars.


title: Almost linear Haldane pseudopotentials and emergent conformal block wave-functions in a Landau level

abstract: We consider a two dimensional (2D) model of particles interacting in a Landau level. We work in a finite disk geometry and take the particles to interact with a linearly decreasing two-body Haldane pseudo-potential. We show that the ground state subspace of this model is spanned by the wave-functions that can be written as polynomial conformal blocks (of an arbitrary conformal field theory) consistent with the filling fraction (scaling dimension). To remove degeneracies, we then add a quadratic perturbation to the Hamiltonian and show that; 1. Conformal blocks constructed using the Moore-Read construction (e.g. Laughlin, Pfaffian, and Read-Rezayi states) remain exact eigenstates of this model in the thermodynamic limit and 2. By tuning an externally imposed single-body $-L_z^2$ potential we can enforce Moore-Read conformal blocks to become exact ground states of this model in the thermodynamic limit. We cannot rule out the possibility of residual degeneracies in this limit. This model has no filling dependence and is comprised only from two-body long-range interactions and external single-body potentials. Our results provide insight into how conformal block wave-functions can emerge in a Landau level.


title: A planetesimal orbiting within the debris disc around a white dwarf star

abstract: Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength.


title: Environment-Induced Rabi Oscillations in the Optomechanical Boson-Boson Model

abstract: We analyze the strong-coupling dynamics of a driven harmonic oscillator whose energy is modulated by a continuum of other bosonic modes. This type of system-bath interaction appears, for example, in optomechanical or equivalent circuit QED setups, where the frequency of a confined photonic mode depends linearly on a fluctuating boundary. Compared to the canonical spin-boson model, where coupling to bath modes only leads to decoherence, the role of the environment in such systems is more complex, since it also provides the only source of nonlinearity. We show that even for an unstructured bath, these environment-induced nonlinearities can dominate over decoherence processes resulting in Rabi oscillations and the formation of highly non-classical states. These findings provide important insights into the non-Markovian dynamics of higher-dimensional open quantum systems and for realizing few-photon optical nonlinearities through strong interactions with a bath.


title: Low-Depth Optical Neural Networks

abstract: Optical neural network (ONN) is emerging as an attractive proposal for machine-learning applications, enabling high-speed computation with low-energy consumption. However, there are several challenges in applying ONN for industrial applications, including the realization of activation functions and maintaining stability. In particular, the stability of ONNs decrease with the circuit depth, limiting the scalability of the ONNs for practical uses. Here we demonstrate how to compress the circuit depth of ONN to scale only logarithmically, leading to an exponential gain in terms of noise robustness. Our low-depth (LD) ONN is based on an architecture, called Optical CompuTing Of dot-Product UnitS (OCTOPUS), which can also be applied individually as a linear perceptron for solving classification problems. Using the standard data set of Letter Recognition, we present numerical evidence showing that LD-ONN can exhibit a significant gain in noise robustness, compared with a previous ONN proposal based on singular-value decomposition [Nature Photonics 11, 441 (2017)].


title: Small-time expansion of the Fokker-Planck kernel for space and time dependent diffusion and drift coefficients

abstract: We study the general solution of the Fokker-Planck equation in d dimensions with arbitrary space and time dependent diffusion matrix and drift term. We show how to construct the solution, for arbitrary initial distributions, as an asymptotic expansion for small time. This generalizes the well-known asymptotic expansion of the heat-kernel for the Laplace operator on a general Riemannian manifold. We explicitly work out the general solution to leading and next-to-leading order in this small-time expansion, as well as to next-to-next-to-leading order for vanishing drift. We illustrate our results on a several examples.


title: The impact of magnetic fields on cold streams feeding galaxies

abstract: High redshift, massive halos are observed to have sustained, high star formation rates, which require that the amount of cold gas in the halo is continuously replenished. The cooling time scale for the hot virialized halo gas is too long to provide the source of cold gas. Supersonic, cold streams have been invoked as a mechanism for feeding massive halos at high redshift and deliver the cold gas required for continued star formation at the rates observed. This mechanism for replenishing the cold gas reservoir is motivated by some cosmological simulations. However, the cold streams are likely to be subject to the supersonic version of the Kelvin-Helmholtz instability (KHI), which eventually leads to stream disruption. Cosmological simulations have yet to obtain the spatial resolution required for understanding the detailed stability properties of cold streams. In this paper, we consider instead an idealized model of magnetized cold streams that we spatially resolve. Using linear theory we show how magnetic fields with dynamically important field strengths do not inhibit the KHI but rather enhance its growth rate. We perform nonlinear simulations of magnetized stream disruption and find that magnetic fields can nevertheless increase stream survival times by suppressing the mixing rate of cold gas with the circumgalactic medium. We find that magnetic fields can allow streams to survive $\sim 2-8$ times longer and, consequently, that streams $\sim 2-8$ times thinner can reach the central galaxy if the magnetic field strength is $\sim 0.3-0.8 \mu$G.


title: Multi-oscillating Boson Stars

abstract: We propose that stable boson stars generically fall within an infinite-parameter family of solutions that oscillate on any number of non-commensurate frequencies. We numerically construct two-frequency solutions and explore their parameter space. These solutions merge with the standard boson star family in the limit where the non-dominating frequencies are turned off. We find that for a fixed energy, these two-frequency solutions can differ considerably in size from standard boson stars.


title: Blue Lurkers: Hidden Blue Stragglers on the M67 Main Sequence Identified from their Kepler/K2 Rotation Periods

abstract: At an age of 4 Gyr, typical solar-type stars in M67 have rotation rates of 20-30 days. Using K2 Campaign 5 and 16 light curves and the spectral archive of the WIYN Open Cluster Study, we identify eleven three-dimensional kinematic members of M67 with anomalously fast rotation periods of 2-8 days, implying ages of less than 1 Gyr. We hypothesize that these anomalously fast rotators have been spun up by mass transfer, mergers, or stellar collisions during dynamical encounters within the last Gyr, and thus represent lower-luminosity counterparts to the blue straggler stars. These 11 candidate post-interaction stellar systems have much in common with the blue stragglers including a high binary fraction (73%), a number of long-period, low-eccentricity binary systems, and in at least one case a UV excess consistent with the presence of a hot white dwarf companion. The identification of these 11 systems provides the first picture of the low-luminosity end of the blue straggler distribution, providing new constraints for detailed binary evolution models and cluster population studies. This result also clearly demonstrates the need to properly account for the impact of binaries on stellar evolution, as significant numbers of post-interaction binaries likely exist on cluster main sequences and in the field. These stars are not always easy to identify, but make up ~10% of the spectroscopic binary population among the solar-type stars in M67.


title: Entropy, Entanglement and Swampland Bounds in DS/dS

abstract: We calculate the entanglement entropy of the de-Sitter (dS) static patch in the context of the DS/dS correspondence. Interestingly, we find that there exists a one parameter family of bulk minimal surfaces that all have the same area. Two of them have appeared earlier in the literature. All of them correctly calculate the dS entropy. One surface yields the entanglement between the two different CFTs that provide the holographic dual of the bulk DS geometry. The second surface describes the entanglement across the horizon in the boundary static patch. The other surfaces describe a mixture of these two concepts. We also show that in the presence of extra matter fields the former entanglement entropy always exceeds the dS entropy. We interpret this result in the context of entropy bounds in de Sitter space and the swampland program.