title: The Outer Rim Simulation: A Path to Many-Core Supercomputers

abstract: We describe the Outer Rim cosmological simulation, one of the largest high-resolution N-body simulations performed to date, aimed at promoting science to be carried out with large-scale structure surveys. The simulation covers a volume of (4.225Gpc)^3 and evolves more than one trillion particles. It was executed on Mira, a BlueGene/Q system at the Argonne Leadership Computing Facility. We discuss some of the computational challenges posed by a system like Mira, a many-core supercomputer, and how the simulation code, HACC, has been designed to overcome these challenges. We have carried out a large range of analyses on the simulation data and we report on the results as well as the data products that have been generated. The full data set generated by the simulation totals more than 5PB of data, making data curation and data handling a large challenge in of itself. The simulation results have been used to generate synthetic catalogs for large-scale structure surveys, including DESI and eBOSS, as well as CMB experiments. A detailed catalog for the LSST DESC data challenges has been created as well. We publicly release some of the Outer Rim halo catalogs, downsampled particle information, and lightcone data.


title: Neutrino masses, leptogenesis and dark matter

abstract: Despite the lack of evidence of new physics at colliders, neutrino masses, dark matter and matter-antimatter asymmetry of the universe require an extension of the Standard Model. After discussing some new concepts and tools in the description of seesaw neutrino models, such as motion in flavour lepton space and the introduction of the bridging matrix that nicely complements the orthogonal matrix, I discuss recent developments in the connections between neutrino data and scenarios of leptogenesis within some well motivated extensions of the Standard Model. Finally, I briefly review a simple unified picture of neutrino masses, leptogenesis and dark matter based on an extension of the seesaw Lagrangian where a non-renormalizable effective operator coupling right-handed neutrinos to the Higgs is introduced. This operator enhances the right-handed neutrino mixing between a coupled heavy right-handed neutrino and a decoupled one playing the role of dark matter. Interference between the two coupled RH neutrinos generate usual $C\!P$ violation so that one can also achieve successful leptogenesis. Although the dark matter right-handed neutrino escapes direct and collider searches, its decays produce a detectable contribution to the very high energy neutrino flux now discovered by IceCube, so that the model is predictive and can be tested at neutrino telescopes.


title: Transit-period search from single-event space-based data: the role of wide-field surveys

abstract: We investigate the optimization of dataset weighting in searching for the orbital period of transiting planets when high-precision space-based data with a single transit event are combined with (relatively) low-precision ground-based (wide-field) data. The optimization stems from the lack of multiple events in the high-precision data and the likely presence of such events in the low-precision data. With noise minimization, we combined two types of frequency spectra: i) spectra that use two fixed transit parameters (moment of the center of the transit and duration of the event) derived from the space data alone; ii) spectra that result from the traditional weighted box signal search with optimized transit parameters for each trial period. We used many mock signals to test the detection power of the method. Marginal or no detections in the ground-based data may lead to secure detections in the combined data with the above weighting. Depending on the coverage and quality of the ground-based data, transit depths of ~0.05% and periods up to ~100days are accessible by the suggested optimum combination of the data.


title: On the statistical distributions of substance moving through the nodes of a channel of network

abstract: We discuss a model of motion of substance through the nodes of a channel of a network. The channel can be modeled by a chain of urns where each urn can exchange substance with the neighboring urns. In addition the urns can exchange substance with the network nodes and the new point is that we include in the model the possibility for exchange of substance among the urns (nodes) and the environment of the network. We consider stationary regime of motion of substance through a finite channel (stationary regime of exchange of substance along the chain of urns) and obtain a class of statistical distributions of substance in the nodes of the channel. Our attention is focused on this class of distributions and we show that for the case of finite channel the obtained class of distributions contains as particular cases truncated versions of the families of distributions of Katz, Ord, Kemp, etc. The theory for the case of infinite chain of urns is presented in the Appendix.


title: The EPR paradox and quantum entanglement at sub-nucleonic scales

abstract: In 1935, in a paper entitled "Can quantum-mechanical description of reality be considered complete?", Einstein, Podolsky, and Rosen (EPR) formulated an apparent paradox of quantum theory. They considered two quantum systems that were initially allowed to interact, and were then later separated. A measurement of a physical observable performed on one system then had to have an immediate effect on the conjugate observable in the other system - even if the systems were causally disconnected! The authors viewed this as a clear indication of the inconsistency of quantum mechanics. In the parton model of the nucleon formulated by Bjorken, Feynman, and Gribov, the partons (quarks and gluons) are viewed by an external hard probe as independent. The standard argument is that, inside the nucleon boosted to an infinite-momentum frame, the parton probed by a virtual photon with virtuality Q is causally disconnected from the rest of the nucleon during the hard interaction. Yet, the parton and the rest of the nucleon have to form a colour-singlet state due to colour confinement, and so have to be in strongly correlated quantum states - we thus encounter the EPR paradox at the sub-nucleonic scale. In this paper, we propose a resolution of this paradox based on the quantum entanglement of partons. We devise an experimental test of entanglement, and carry it out using data on proton-proton collisions from the Large Hadron Collider (LHC). Our results provide a strong direct indication of quantum entanglement at sub-nucleonic scales.


title: The equilibrium landscape of the Heisenberg spin chain

abstract: We characterise the equilibrium landscape, the entire manifold of local equilibrium states, of an interacting integrable quantum model. Focusing on the isotropic Heisenberg spin chain, we describe in full generality two complementary frameworks for addressing equilibrium ensembles: the functional integral Thermodynamic Bethe Ansatz approach, and the lattice regularisation transfer matrix approach. We demonstrate the equivalence between the two, and in doing so clarify several subtle features of generic equilibrium states. In particular we explain the breakdown of the canonical Y-system, which reflects a hidden structure in the parametrisation of equilibrium ensembles.


title: The Magnetic Properties of Heating Events on High-Temperature Active Region Loops

abstract: Understanding the relationship between the magnetic field and coronal heating is one of the central problems of solar physics. However, studies of the magnetic properties of impulsively heated loops have been rare. We present results from a study of 34 evolving coronal loops observed in the Fe XVIII line component of AIA/SDO 94 A filter images from three active regions with different magnetic conditions. We show that the peak intensity per unit cross-section of the loops depends on their individual magnetic and geometric properties. The intensity scales proportionally to the average field strength along the loop ($B_{avg}$) and inversely with the loop length ($L$) for a combined dependence of $(B_{avg}/L)^{0.52\pm0.13}$. These loop properties are inferred from magnetic extrapolations of the photospheric HMI/SDO line-of-sight and vector magnetic field in three approximations: potential and two Non Linear Force-Free (NLFF) methods. Through hydrodynamic modeling (EBTEL model) we show that this behavior is compatible with impulsively heated loops with a volumetric heating rate that scales as $\epsilon_H\sim B_{avg}^{0.3\pm0.2}/L^{0.2\pm^{0.2}_{0.1}}$.


title: Blueprint for deterministic all-optical switching of magnetization

abstract: We resolve a significant controversy about how to understand and engineer single-shot all-optical switching of magnetization in ferrimagnets using femto- or picosecond-long heat pulses. By realistically modelling a generic ferrimagnet as two coupled macrospins, we comprehensively show that the net magnetization can be reversed via different pathways, using a heat pulse with duration spanning all relevant timescales within the non-adiabatic limit. This conceptual understanding is fully validated by experiments studying the material and optical limits at which the switching process in GdFeCo alloys loses its reliability. Our interpretation and results constitute a blueprint for understanding how deterministic all-optical switching can be achieved in alternative ferrimagnets using short thermal pulses.


title: Energy self-extraction of a Kerr black hole through its frame-dragged force-free magnetosphere

abstract: It is shown that when only the condition} $0<\Omega_{\rm F}<\Omega_{\rm H}$ is satisfied, the Kerr black hole frame-drags its surrounding force-free magnetosphere with the field-line-angular-velocity (FLAV) $\Omega_{\rm F}$, where $\Omega_{\rm H}$ is the horizon angular-velocity. Then, the zero-angular-momentum-observers (ZAMOs) circulating with the frame-dragging-angular-velocity $\omega$ will see that the `null surface' S$_{\rm N}$ where $\omega_{\rm N}=\Omega_{\rm F}$ always exists. They will see that the outer domain $D_{\rm (out)}$ outside S$_{\rm N}$ is prograde-rotating with $\Omega_{{\rm F} \omega}>0$, whereas the inner domain $D_{\rm (in)}$ inside is retrograde-rotating with $\Omega_{{\rm F} \omega}<0$, where $\Omega_{{\rm F} \omega}=\Omega_{{\rm F}} - \omega$ denotes the ZAMO-FLAV. `This surface' S$_{\rm N}$ must be the magneto-centrifugal divider of the force-free magnetosphere, with a kind of plasma-shed on it. Subsequently, the force-free and freezing-in conditions break down on S$_{\rm N}$, thereby allowing the particle-current sources to be set up on S$_{\rm N}$. `This surface' also is the ZAM-surface S$_{\rm ZAMD}$, on which no flow of angular momentum nor electric current can cross. Because the electric field ${\bf E}_{\rm p}$ reverses sign on S$_{\rm N}$, the Poynting flux reverses direction from outward to inward on S$_{\rm N}$. An electromagnetic self-extraction of energy will be possible only through the frame-dragged magnetosphere, with the inner domain $D_{\rm (in)}$ nested between the horizon and `this surface' S$_{\rm N}$, in order to comply with the first and second laws of thermodynamics.


title: PowerNet: Neural Power Demand Forecasting in Smart Grid

abstract: Power demand forecasting is a critical task for achieving efficiency and reliability in power grid operation. Accurate forecasting allows grid operators to better maintain the balance of supply and demand as well as to optimize operational cost for generation and transmission. This article proposes a novel neural network architecture PowerNet, which can incorporate multiple heterogeneous features, such as historical energy consumption data, weather data, and calendar information, for the power demand forecasting task. Compared to two recent works based on Gradient Boosting Tree (GBT) and Support Vector Regression (SVR), PowerNet demonstrates a decrease of 33.3% and 14.3% in forecasting error, respectively. We further provide empirical results the two operational considerations that are crucial when using PowerNet in practice, i.e., how far in the future the model can forecast with a decent accuracy and how often we should re-train the forecasting model to retain its modeling capability. Finally, we briefly discuss a multilayer anomaly detection approach based on PowerNet.


title: Forecasting Interacting Vacuum-Energy Models using Gravitational Waves

abstract: The physics of the dark sector has remained one of the controversial areas of modern cosmology at present and hence it naturally attracts massive attention to the scientific community. With the developments of the astronomical data, the physics of the dark sector is becoming much more transparent than it was some twenty years back. The detection of gravitational waves (GWs) has now opened a cluster of possibilities in the cosmological regime. Being motivated by the detection of GWs and its possible impact on the physics of dark matter and dark energy, in this work we focus on the interacting dark energy models. Assuming the simplest possibility in which the vacuum energy with equation-of-state $w_x =-1$ is allowed to interact with the pressureless dark matter, we have extracted the constraints of the cosmological parameters. We find that the addition of the GWs data to the CMB measurements significantly improves up to a factor 4 of the parameters space, and up to a factor 2 for the full combination of current cosmological datasets, namely CMB+BAO+Pantheon+RSD+R16+CC+WL. The most affected parameters by the inclusion of the GWs are $\Omega_ch^2$, $\theta_{MC}$, $\xi$, and the derived parameters $\Omega_{m0}$, $\sigma_8$ and $H_0$.


title: Shear-induced instabilities of flows through submerged vegetation

abstract: We consider the instabilities of flows through a submerged canopy, and show how the full governing equations of the fluid-structure interactions can be reduced to a compact framework that captures many key features of vegetative flow. By modelling the canopy as a collection of homogeneous elastic beams, we predict the steady configuration of the plants in response to a unidirectional flow. This treatment couples the beam equations in the canopy to the fluid momentum equations. Our linear stability analysis suggests new insights into the development of instabilities at the surface of the vegetative region. In particular, we show that shear at the top of the canopy is a dominant factor in determining the onset of instabilities known as monami. Based on numerical and asymptotic analysis of the generalised eigenvalue problem, the system is shown to be stable if the canopy is sufficiently sparse or if the plants are sufficiently flexible.


title: Generalized Cline's Formula for G-Drazin Inverse

abstract: Let $R$ be an associative ring with an identity and suppose that $a,b,c,d \in R$ satisfy $bdb = bac,dbd = acd$. If $ac$ has generalized Drazin ( respectively, pseudo Drazin, Drazin) inverse, we prove that $bd$ has generalized Drazin (respectively, pseudo Drazin, Drazin) inverse. In particular, as applications, we obtain new common spectral properties of bounded linear operators over Banach spaces.


title: Deep learning enabled superfast and accurate M^2 evaluation for fiber beams

abstract: We introduce deep learning technique to predict the beam propagation factor M^2 of the laser beams emitting from few-mode fiber for the first time, to the best of our knowledge. The deep convolutional neural network (CNN) is trained with paired data of simulated near-field beam patterns and their calculated M^2 value, aiming at learning a fast and accurate mapping from the former to the latter. The trained deep CNN can then be utilized to evaluate M^2 of the fiber beams from single beam patterns. The results of simulated testing samples have shown that our scheme can achieve an averaged prediction error smaller than 2% even when up to 10 eigenmodes are involved in the fiber. The error becomes slightly larger when heavy noises are added into the input beam patterns but still smaller than 2.5%, which further proves the accuracy and robustness of our method. Furthermore, the M^2 estimation takes only about 5 ms for a prepared beam pattern with one forward pass, which can be adopted for real-time M^2 determination with only one supporting Charge-Coupled Device (CCD). The experimental results further prove the feasibility of our scheme. Moreover, the method we proposed can be confidently extended to other kinds of beams provided that adequate training samples are accessible. Deep learning paves the way to superfast and accurate M^2 evaluation with very low experimental efforts.


title: Quantization of massive spinor fields around global monopoles

abstract: The renormalized quantum stress energy tensor $\left<T_{\mu}^{\nu}\right>_{ren}$ for a massive spinor field around global monopoles is constructed within the framework of Schwinger-DeWitt approximation, valid whenever the Compton length of the quantum field is much less than the characteristic radius of the curvature of the background geometry. The results obtained shows that the quantum massive spinor field in the global monopole spacetime violates all the pointwise energy conditions.


title: Transmission of natural scene images through a multimode fibre

abstract: The optical transport of images through a multimode fibre remains an outstanding challenge with applications ranging from optical communications to neuro-imaging. State of the art approaches either involve measurement and control of the full complex field transmitted through the fibre or, more recently, training of artificial neural networks that however, are typically limited to image classes belong to the same class as the training data set. Here we implement a method that statistically reconstructs the inverse transformation matrix for the fibre. We demonstrate imaging at high frame rates, high resolutions and in full colour of natural scenes, thus demonstrating general-purpose imaging capability. Real-time imaging over long fibre lengths opens alternative routes to exploitation for example for secure communication systems, novel remote imaging devices, quantum state control processing and endoscopy.


title: The quasiclassical theory for interfaces with spin--orbit coupling

abstract: In recent years, substantial progress has been made regarding the application of quasiclassical theory on superconducting hybrid structures. This theoretical framework is reliant on a proper set of boundary conditions in order to describe multilayered systems. With the advent of the field of superconducting spintronics, systems which combine heavy metal layers, in which there is large spin--orbit coupling, with ferromagnets have received a great deal of attention, due to their potential for generating long range triplet superconductivity. In contrast to interfaces of strongly spin polarized materials, which are well understood, a quasiclassical theory for interfaces in systems where there is significant spin--orbit coupling does not yet exist. After reviewing the quasiclassical theory for interfaces, we here solve this problem by deriving a new set of boundary conditions which take spin--orbit coupling explicitly into account. We then go on to apply these boundary conditions to a superconductor-ferromagnet (SF) bilayer and an SFS Josephson weak link, demonstrating the emergence of long range triplet superconductivity in these systems.


title: Disruption of giant molecular clouds and formation of bound star clusters under the influence of momentum stellar feedback

abstract: Energetic feedback from star clusters plays a pivotal role in shaping the dynamical evolution of giant molecular clouds (GMCs). To study the effects of stellar feedback on the star formation efficiency of the clouds and the dynamical response of embedded star clusters, we perform a suite of isolated GMC simulations with star formation and momentum feedback subgrid models using the moving-mesh hydrodynamics code \textsc{Arepo}. The properties of our simulated GMCs span a wide range of initial mass, radius, and velocity configurations. We find that the ratio of the final stellar mass to the total cloud mass, $\epsilon_{\rm int}$, scales strongly with the initial cloud surface density and momentum feedback strength. This correlation is explained by an analytic model that considers force balancing between gravity and momentum feedback. For all simulated GMCs, the stellar density profiles are systematically steeper than that of the gas at the epochs of the peaks of star formation, suggesting a centrally concentrated stellar distribution. We also find that star clusters are always in a sub-virial state with a virial parameter $\sim0.6$ prior to gas expulsion. Both the sub-virial dynamical state and steeper stellar density profiles prevent clusters from dispersal during the gas removal phase of their evolution. The final cluster bound fraction is a continuously increasing function of $\epsilon_{\rm int}$. GMCs with star formation efficiency smaller than 0.5 are still able to form clusters with large bound fractions.


title: The ExoEarth Yield Landscape for Future Direct Imaging Space Telescopes

abstract: The expected yield of potentially Earth-like planets is a useful metric for designing future exoplanet-imaging missions. Recent yield studies of direct-imaging missions have focused primarily on yield methods and trade studies using "toy" models of missions. Here we increase the fidelity of these calculations substantially, adopting more realistic exoplanet demographics as input, an improved target list, and a realistic distribution of exozodi levels. Most importantly, we define standardized inputs for instrument simulations, use these standards to directly compare the performance of realistic instrument designs, include the sensitivity of coronagraph contrast to stellar diameter, and adopt engineering-based throughputs and detector parameters. We apply these new high-fidelity yield models to study several critical design trades: monolithic vs segmented primary mirrors, on-axis vs off-axis secondary mirrors, and coronagraphs vs starshades. We show that as long as the gap size between segments is sufficiently small, there is no difference in yield for coronagraph-based missions with monolithic off-axis telescopes and segmented off-axis telescopes, assuming that the requisite engineering constraints imposed by the coronagraph can be met in both scenarios. We show that there is currently a factor of ~2 yield penalty for coronagraph-based missions with on-axis telescopes compared to off-axis telescopes, and note that there is room for improvement in coronagraph designs for on-axis telescopes. We also reproduce previous results in higher fidelity showing that the yields of coronagraph-based missions continue to increase with aperture size while the yields of starshade-based missions turnover at large apertures if refueling is not possible. Finally, we provide absolute yield numbers with uncertainties that include all major sources of astrophysical noise to guide future mission design.


title: Level Repulsion and Dynamics in the Finite One-Dimensional Anderson Model

abstract: This work shows that dynamical features typical of full random matrices can be observed also in the simple finite one-dimensional (1D) noninteracting Anderson model with nearest neighbor couplings. In the thermodynamic limit, all eigenstates of this model are exponentially localized in configuration space for any infinitesimal onsite disorder strength $W$. But this is not the case when the model is finite and the localization length is larger than the system size $L$, which is a picture that can be experimentally investigated. We analyze the degree of energy-level repulsion, the structure of the eigenstates, and the time evolution of the finite 1D Anderson model as a function of the parameter $\xi \propto (W^2 L)^{-1}$. As $\xi$ increases, all energy-level statistics typical of random matrix theory are observed. The statistics are reflected in the corresponding eigenstates and also in the dynamics. We show that the probability in time to find a particle initially placed on the first site of an open chain decays as fast as in full random matrices and much faster that when the particle is initially placed far from the edges. We also see that at long times, the presence of energy-level repulsion manifests in the form of the correlation hole. In addition, our results demonstrate that the hole is not exclusive to random matrix statistics, but emerges also for $W=0$, when it is in fact deeper.


title: Collider Signals of the Mirror Twin Higgs through the Hypercharge Portal

abstract: We consider the collider signals arising from kinetic mixing between the hypercharge gauge boson of the Standard Model and its twin counterpart in the Mirror Twin Higgs model, in the framework in which the twin photon is massive. Through the mixing, the Standard Model fermions acquire charges under the mirror photon and the mirror Z boson. We determine the current experimental bounds on this scenario, and show that the mixing can be large enough to discover both the twin photon and the twin Z at the LHC, or at a future 100 TeV hadron collider, with dilepton resonances being a particularly conspicuous signal. We show that, in simple models, measuring the masses of both the mirror photon and mirror Z, along with the corresponding event rates in the dilepton channel, overdetermines the system, and can be used to test these theories.


title: Building the Evryscope: Hardware Design and Performance

abstract: The Evryscope is a telescope array designed to open a new parameter space in optical astronomy, detecting short timescale events across extremely large sky areas simultaneously. The system consists of a 780 MPix 22-camera array with an 8150 sq. deg. field of view, 13" per pixel sampling, and the ability to detect objects down to Mg=16 in each 2 minute dark-sky exposure. The Evryscope, covering 18,400 sq.deg. with hours of high-cadence exposure time each night, is designed to find the rare events that require all-sky monitoring, including transiting exoplanets around exotic stars like white dwarfs and hot subdwarfs, stellar activity of all types within our galaxy, nearby supernovae, and other transient events such as gamma ray bursts and gravitational-wave electromagnetic counterparts. The system averages 5000 images per night with ~300,000 sources per image, and to date has taken over 3.0M images, totaling 250TB of raw data. The resulting light curve database has light curves for 9.3M targets, averaging 32,600 epochs per target through 2018. This paper summarizes the hardware and performance of the Evryscope, including the lessons learned during telescope design, electronics design, a procedure for the precision polar alignment of mounts for Evryscope-like systems, robotic control and operations, and safety and performance-optimization systems. We measure the on-sky performance of the Evryscope, discuss its data-analysis pipelines, and present some example variable star and eclipsing binary discoveries from the telescope. We also discuss new discoveries of very rare objects including 2 hot subdwarf eclipsing binaries with late M-dwarf secondaries (HW Vir systems), 2 white dwarf / hot subdwarf short-period binaries, and 4 hot subdwarf reflection binaries. We conclude with the status of our transit surveys, M-dwarf flare survey, and transient detection.


title: Galaxy properties as revealed by MaNGA. I. Constraints on Initial Mass Function and M$_{*}$/L gradients in ellipticals

abstract: We estimate ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distance. This paper is the first of a series and is based on a sample of pure elliptical galaxies at z $\le$ 0.08. We show that the properties of the inner regions of Es with the largest luminosity (L$_r$) and central velocity dispersion ($\sigma_0$) are consistent with those associated with the commonly used Salpeter IMF, whereas a Kroupa-like IMF is a better description at $\sim$ 0.8R/Re (assuming [Ti/Fe] variations are limited). For these galaxies the stellar mass-to-light ratio decreases at most by a factor of 2 from the central regions to Re. In contrast, for lower L$_r$ and $\sigma_0$ galaxies, the IMF is shallower and M$_{*}$/L$_r$ in the central regions is similar to the outskirts. Although a factor of 2 is smaller than previous reports based on a handful of galaxies, it is still large enough to matter for dynamical mass estimates. Accounting self-consistently for these gradients when estimating both M$_{*}$ and M$_{dyn}$ brings the two into good agreement: gradients reduce M$_{dyn}$ by $\sim$ 0.2 dex while only slightly increasing the M$_{*}$ inferred using a Kroupa IMF. This is a different resolution of the M$_{*}$-M$_{dyn}$ discrepancy than has been followed in the recent literature where M$_{*}$ of massive galaxies is increased by adopting a Salpeter IMF while leaving Mdyn unchanged. A companion paper discusses how stellar population differences are even more pronounced if one separates slow from fast rotators.


title: Supernovae Sparked By Dark Matter in White Dwarfs

abstract: It was recently demonstrated that asymmetric dark matter can ignite supernovae by collecting and collapsing inside lone sub-Chandrasekhar mass white dwarfs, and that this may be the cause of Type Ia supernovae. A ball of asymmetric dark matter accumulated inside a white dwarf and collapsing under its own weight, sheds enough gravitational potential energy through scattering with nuclei, to spark the fusion reactions that precede a Type Ia supernova explosion. In this article we elaborate on this mechanism and use it to place new bounds on interactions between nucleons and asymmetric dark matter for masses $m_{X} = 10^{6}-10^{16}$ GeV. Interestingly, we find that for dark matter more massive than $10^{11}$ GeV, Type Ia supernova ignition can proceed through the Hawking evaporation of a small black hole formed by the collapsed dark matter. We also identify how a cold white dwarf's Coulomb crystal structure substantially suppresses dark matter-nuclear scattering at low momentum transfers, which is crucial for calculating the time it takes dark matter to form a black hole. Higgs and vector portal dark matter models that ignite Type Ia supernovae are explored.


title: Massive star formation via torus accretion: the effect of photoionization feedback

abstract: The formation of massive stars is a long standing problem. Although a number of theories of massive star formation exist, ideas appear to converge to a disk-mediated accretion scenario. Here we present radiative hydrodynamic simulations of a star accreting mass via a disk embedded in a torus. We use a Monte Carlo based radiation hydrodynamics code to investigate the impact that ionizing radiation has on the torus. Ionized regions in the torus midplane are found to be either gravitationally trapped or in pressure driven expansion depending on whether or not the size of the ionized region exceeds a critical radius. Trapped Hii regions in the torus plane allow accretion to progress, while expanding Hii regions disrupt the accretion torus preventing the central star from aggregating more mass, thereby setting the star's final mass. We obtain constraints for the luminosities and torus densities that lead to both scenarios.


title: The Multi-Messenger Matrix: the Future of Neutron Star Merger Constraints on the Nuclear Equation of State

abstract: The electromagnetic (EM) signal of a binary neutron star (BNS) merger depends sensitively on the total binary mass, $M_{\rm tot}$, relative to various threshold masses set by the neutron star (NS) equation-of-state (EOS), parameterized through the NS maximum mass, $M_{\rm TOV}$, and characteristic radius, $R_{1.6}$. EM observations of a BNS merger detected through its gravitational wave (GW) emission, which are of sufficient quality to ascertain the identity of the merger remnant, can therefore constrain the values of $M_{\rm TOV}$ and $R_{1.6}$, given the tight connection between $M_{\rm tot}$ and the well-measured chirp mass. We elucidate the present and future landscape of EOS constraints from BNS mergers, introducing the "Multi-Messenger Matrix", a mapping between GW and EM measurables that defines the ranges of event chirp masses which provide the most leverage on constraining the EOS. By simulating a population of BNS mergers drawn from the Galactic double NS mass distribution we show that $\sim 10$ joint detections can constrain $M_{\rm TOV}$ and $R_{1.6}$ to several percent level where systematic uncertainties may become significant. Current EOS constraints imply that most mergers will produce supramassive or hypermassive remnants, a smaller minority (possibly zero) will undergo prompt-collapse, while at most only a few percent of events will form indefinitely stable NSs. In support of the envisioned program, we advocate in favor of LIGO/Virgo releasing chirp mass estimates as early as possible to the scientific community, enabling observational resources to be allocated in the most efficient way to maximize the scientific gain from multi-messenger discoveries.


title: Galaxy properties as revealed by MaNGA II. Differences in stellar populations of slow and fast rotator ellipticals and dependence on environment

abstract: We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $\alpha$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $\sigma_0$. This explains why when one combines SRs and FRs, objects which are small for their $L_r$ and $\sigma_0$ tend to be younger. Their SP gradients are also different. 2) Ignoring the FR/SR dichotomy leads one to conclude that compact galaxies are older than their larger counterparts of the same mass, even though almost the opposite is true for FRs and SRs individually. 3) SRs with $\sigma_0\le 250$ km s$^{-1}$ are remarkably homogeneous within $\sim R_e$: they are old, $\alpha$-enhanced and only slightly super-solar in metallicity. These SRs show no gradients in age and $M_*/L_r$, negative gradients in metallicity, and slightly positive gradients in [$\alpha$/Fe] (the latter are model dependent). SRs with $\sigma_0\ge 250$ km $s^{-1}$ are slightly younger and more metal rich, contradicting previous work suggesting that age increases with $\sigma_0$. They also show larger $M_*/L_r$ gradients. 4) Self-consistently accounting for $M_*/L$ gradients yields $M_{\rm dyn}\approx M_*$ because gradients reduce $M_{\rm dyn}$ by $\sim 0.2$ dex while only slightly increasing the $M_*$ inferred using a Kroupa (not Salpeter) IMF. 5) The FR population all but disappears above $M_*\ge 3\times 10^{11}M_\odot$; this is the same scale at which the size-mass correlation and other scaling relations change. Our results support the finding that this is an important mass scale which correlates with the environment and above which mergers matter.


title: Quantum Theory of Triboelectricity

abstract: We propose a microphysical theory of the triboelectric effect by which mechanical rubbing separates charges across the interface between two materials. Surface electrons are treated as an open system coupled to two baths, corresponding to the bulks. Extending Zel'dovich's theory of bosonic superradiance, we show that motion-induced population inversion can generate an electromotive force. We argue that this is consistent with the basic phenomenology of triboelectrification and triboluminescence as irreversible processes, and we suggest how to carry out more precise experimental tests.

categories: cond-mat.mes-hall cond-mat.stat-mech hep-th quant-ph


title: Real-virtual corrections to Higgs boson pair production at NNLO: three closed top quark loops

abstract: We compute the real-radiation corrections to Higgs boson pair production at next-to-next-to-leading order in QCD, in an expansion for large top quark mass. We concentrate on the radiative corrections to the interference contribution from the next-to-leading order one-particle reducible and the leading order amplitudes. This is a well defined and gauge invariant subset of the full real-virtual corrections to the inclusive cross section. We obtain analytic results for all phase-space master integrals both as an expansion around the threshold and in an exact manner in terms of Goncharov polylogarithms.


title: Towards overcoming the entanglement barrier when simulating long-time evolution

abstract: Quantum many-body systems out of equilibrium pose some of the most intriguing questions in physics. Unfortunately, numerically keeping track of time evolution of states under Hamiltonian dynamics constitutes a severe challenge for all known methods. Prominently, tensor network methods are marred by an entanglement blowup, which allows to simulate systems following global quenches only to constant time. In this work, we present a scheme that allows to significantly extend the simulation time for interacting fermionic or equivalent spin systems. In the past when keeping track of evolution in one-dimensional real space, the subspace parametrised by real-space matrix product states satisfying area laws - often dubbed the "physical corner" of Hilbert space - was chosen as variational set. In contrast, if the manifold containing both tensor network states and fermionic mode transformations is chosen, significantly longer times can be achieved. We argue and our results suggest that in many cases it is genuine correlations between modes that is the actual limiting factor: The system at hand is for intermediate times contained in the "physical corner", but a different one than what is commonly assumed.


title: On the fixed-parameter tractability of the maximum connectivity improvement problem

abstract: In the Maximum Connectivity Improvement (MCI) problem, we are given a directed graph $G=(V,E)$ and an integer $B$ and we are asked to find $B$ new edges to be added to $G$ in order to maximize the number of connected pairs of vertices in the resulting graph. The MCI problem has been studied from the approximation point of view. In this paper, we approach it from the parameterized complexity perspective in the case of directed acyclic graphs. We show several hardness and algorithmic results with respect to different natural parameters. Our main result is that the problem is $W[2]$-hard for parameter $B$ and it is FPT for parameters $|V| - B$ and $\nu$, the matching number of $G$. We further characterize the MCI problem with respect to other complementary parameters.


title: Oscillating solutions for nonlinear equations involving the Pucci's extremal operators

abstract: This paper deals with the following nonlinear equations \[ \mathcal{M}_{\lambda,\Lambda}^\pm(D^2 u)+g(u)=0 \qquad \hbox{ in }\mathbb{R}^N, \] where $\mathcal{M}_{\lambda,\Lambda}^\pm$ are the Pucci's extremal operators, for $N \ge 1$ and under the assumption $g'(0)>0$. We show the existence of oscillating solutions, namely with an unbounded sequence of zeros. Moreover these solutions are periodic, if $N=1$, while they are radial symmetric and decay to zero at infinity with their derivatives, if $N\ge 2$.


title: Accelerating two projection methods via perturbations with application to Intensity-Modulated Radiation Therapy

abstract: Constrained convex optimization problems arise naturally in many real-world applications. One strategy to solve them in an approximate way is to translate them into a sequence of convex feasibility problems via the recently developed level set scheme and then solve each feasibility problem using projection methods. However, if the problem is ill-conditioned, projection methods often show zigzagging behavior and therefore converge slowly. To address this issue, we exploit the bounded perturbation resilience of the projection methods and introduce two new perturbations which avoid zigzagging behavior. The first perturbation is in the spirit of $k$-step methods and uses gradient information from previous iterates. The second uses the approach of surrogate constraint methods combined with relaxed, averaged projections. We apply two different projection methods in the unperturbed version, as well as the two perturbed versions, to linear feasibility problems along with nonlinear optimization problems arising from intensity-modulated radiation therapy (IMRT) treatment planning. We demonstrate that for all the considered problems the perturbations can significantly accelerate the convergence of the projection methods and hence the overall procedure of the level set scheme. For the IMRT optimization problems the perturbed projection methods found an approximate solution up to 4 times faster than the unperturbed methods while at the same time achieving objective function values which were 0.5 to 5.1% lower.


title: Counterexamples to hyperkahler Kirwan surjectivity

abstract: Suppose that M is a complete hyperkahler manifold with a compact Lie group K acting via hyperkahler isometries and with hyperkahler moment map $(\mu_{\mathbb{C}}, \mu_{\mathbb{R}}): M\rightarrow \mathfrak{k}^*\otimes\operatorname{Im}(\mathbb{H})$. It is a long-standing problem to determine when the hyperkahler Kirwan map $H^*_K(M,\mathbb{Q}) \longrightarrow H^*(M//K, \mathbb{Q})$ is surjective. We show that for each $n\geq 2$, the natural $U(n)$-action on $M = T^*(SL_n\times\mathbb{C}^n)$ admits a hyperkahler quotient for which the hyperkahler Kirwan map fails to be surjective. As a tool, we establish a ``Kahler $=$ GIT quotient'' assertion for products of cotangent bundles of reductive groups, equipped with the Kronheimer metric, and representations.


title: Knowing When to Stop: Evaluation and Verification of Conformity to Output-size Specifications

abstract: Models such as Sequence-to-Sequence and Image-to-Sequence are widely used in real world applications. While the ability of these neural architectures to produce variable-length outputs makes them extremely effective for problems like Machine Translation and Image Captioning, it also leaves them vulnerable to failures of the form where the model produces outputs of undesirable length. This behavior can have severe consequences such as usage of increased computation and induce faults in downstream modules that expect outputs of a certain length. Motivated by the need to have a better understanding of the failures of these models, this paper proposes and studies the novel output-size modulation problem and makes two key technical contributions. First, to evaluate model robustness, we develop an easy-to-compute differentiable proxy objective that can be used with gradient-based algorithms to find output-lengthening inputs. Second and more importantly, we develop a verification approach that can formally verify whether a network always produces outputs within a certain length. Experimental results on Machine Translation and Image Captioning show that our output-lengthening approach can produce outputs that are 50 times longer than the input, while our verification approach can, given a model and input domain, prove that the output length is below a certain size.


title: Classifying integrable spin-1/2 chains with nearest neighbour interactions

abstract: We classify all fundamental integrable spin chains with two-dimensional local Hilbert space which have regular R-matrices of difference form. This means that the R-matrix underlying the integrable structures is of the form R(u,v)=R(u-v) and reduces to the permutation operator at some particular point. We find a total of 14 independent solutions, 8 of which correspond to well-known eight or lower vertex models. The remaining 6 models appear to be new and some have peculiar properties such as not being diagonalizable or being nilpotent. Furthermore, for even R-matrices, we find a bijection between solutions of the Yang-Baxter equation and the graded Yang-Baxter equation which extends our results to the graded two-dimensional case.

categories: math-ph cond-mat.stat-mech hep-th math.MP nlin.SI


title: Resolving the degeneracy in top quark Yukawa coupling with Higgs pair production

abstract: The top quark Yukawa coupling ($y_t$) can be modified by two dimension-six operators $\mathcal{O}_H$ and $\mathcal{O}_y$ with the corresponding Wilson coefficients $c_H$ and $c_y$, whose individual contribution cannot be distinguished by measuring $y_t$ alone. However, such a degeneracy can be resolved with Higgs boson pair production. In this work we explore the potential of resolving the degeneracy of the unknown Wilson coefficients $c_H$ and $c_y$ at the 14~TeV LHC and the 100~TeV hadron colliders. Combining the information of the single Higgs production, $t\bar{t}h$ associated production and Higgs pair production, the individual contribution of $c_H$ and $c_y$ to $y_t$ can be separated. Regardless of the value of $c_H$, the Higgs pair production can give a strong constraint on $c_y$ at the 100 TeV hadron collider. We further show that it is possible to differentiate various $c_y$ and $c_t$ values predicted in several benchmark models.


title: Relevant features for Gender Classification in NIR Periocular Images

abstract: Most gender classifications methods from NIR images have used iris information. Recent work has explored the use of the whole periocular iris region which has surprisingly achieve better results. This suggests the most relevant information for gender classification is not located in the iris as expected. In this work, we analyze and demonstrate the location of the most relevant features that describe gender in periocular NIR images and evaluate its influence its classification. Experiments show that the periocular region contains more gender information than the iris region. We extracted several features (intensity, texture, and shape) and classified them according to its relevance using the XgBoost algorithm. Support Vector Machine and nine ensemble classifiers were used for testing gender accuracy when using the most relevant features. The best classification results were obtained when 4,000 features located on the periocular region were used (89.22\%). Additional experiments with the full periocular iris images versus the iris-Occluded images were performed. The gender classification rates obtained were 84.35\% and 85.75\% respectively. We also contribute to the state of the art with a new database (UNAB-Gender). From results, we suggest focussing only on the surrounding area of the iris. This allows us to realize a faster classification of gender from NIR periocular images.


title: Analog Weights in ReRAM DNN Accelerators

abstract: Artificial neural networks have become ubiquitous in modern life, which has triggered the emergence of a new class of application specific integrated circuits for their acceleration. ReRAM-based accelerators have gained significant traction due to their ability to leverage in-memory computations. In a crossbar structure, they can perform multiply-and-accumulate operations more efficiently than standard CMOS logic. By virtue of being resistive switches, ReRAM switches can only reliably store one of two states. This is a severe limitation on the range of values in a computational kernel. This paper presents a novel scheme in alleviating the single-bit-per-device restriction by exploiting frequency dependence of v-i plane hysteresis, and assigning kernel information not only to the device conductance but also partially distributing it to the frequency of a time-varying input. We show this approach reduces average power consumption for a single crossbar convolution by up to a factor of x16 for an unsigned 8-bit input image, where each convolutional process consumes a worst-case of 1.1mW, and reduces area by a factor of x8, without reducing accuracy to the level of binarized neural networks. This presents a massive saving in computing cost when there are many simultaneous in-situ multiply-and-accumulate processes occurring across different crossbars.


title: Universality of the time constant for $2D$ critical first-passage percolation

abstract: We consider first-passage percolation (FPP) on the triangular lattice with vertex weights $(t_v)$ whose common distribution function $F$ satisfies $F(0)=1/2$. This is known as the critical case of FPP because large (critical) zero-weight clusters allow travel between distant points in time which is sublinear in the distance. Denoting by $T(0,\partial B(n))$ the first-passage time from $0$ to $\{x : \|x\|_\infty = n\}$, we show existence of the "time constant'' and find its exact value to be \[ \lim_{n \to \infty} \frac{T(0,\partial B(n))}{\log n} = \frac{I}{2\sqrt{3}\pi} \text{ almost surely}, \] where $I = \inf\{x > 0 : F(x) > 1/2\}$ and $F$ is any critical distribution for $t_v$. This result shows that the time constant is universal and depends only on the value of $I$. Furthermore, we find the exact value of the limiting normalized variance, which is also only a function of $I$, under the optimal moment condition on $F$. The proof method also shows an analogous universality on other two-dimensional lattices, assuming the time constant exists.


title: Mass rigidity for hyperbolic manifolds

abstract: We prove the rigidity of positive mass theorem for asymptotically hyperbolic manifolds. Namely, if the mass equality holds, then the manifold is isometric to hyperbolic space. The result was previously proven for spin manifolds or under special asymptotics.


title: Parameterized algorithms for Partial vertex covers in bipartite graphs

abstract: In the weighted partial vertex cover problem (WPVC), we are given a graph $G=(V,E)$, cost function $c:V\rightarrow N$, profit function $p:E\rightarrow N$, and positive integers $R$ and $L$. The goal is to check whether there is a subset $V'\subseteq V$ of cost at most $R$, such that the total profit of edges covered by $V'$ is at least $L$. In this paper we study the fixed-parameter tractability of WPVC in bipartite graphs (WPVCB). By extending the methods of Amini et al., we show that WPVCB is FPT with respect to $R$ if $c\equiv 1$. On the negative side, it is $W[1]$-hard for arbitrary $c$, even when $p\equiv 1$. In particular, WPVCB is $W[1]$-hard parameterized by $R$. We complement this negative result by proving that for bounded-degree graphs WPVC is FPT with respect to $R$. The same result holds for the case of WPVCB when we allow to take only one fractional vertex. Additionally, we show that WPVC is FPT with respect to $L$. Finally, we discuss a variant of PVCB in which the edges covered are constrained to include a matching of prescribed size and derive a paramterized algorithm for the same.


title: RevealNet: Seeing Behind Objects in RGB-D Scans

abstract: During 3D reconstruction, it is often the case that people cannot scan each individual object from all views, resulting in missing geometry in the captured scan. This missing geometry can be fundamentally limiting for many applications, e.g., a robot needs to know the unseen geometry to perform a precise grasp on an object. Thus, we introduce the task of semantic instance completion: from an incomplete RGB-D scan of a scene, we aim to detect the individual object instances and infer their complete object geometry. This will open up new possibilities for interactions with objects in a scene, for instance for virtual or robotic agents. We tackle this problem by introducing RevealNet, a new data-driven approach that jointly detects object instances and predicts their complete geometry. This enables a semantically meaningful decomposition of a scanned scene into individual, complete 3D objects, including hidden and unobserved object parts. RevealNet is an end-to-end 3D neural network architecture that leverages joint color and geometry feature learning. The fully-convolutional nature of our 3D network enables efficient inference of semantic instance completion for 3D scans at scale of large indoor environments in a single forward pass. We show that predicting complete object geometry improves both 3D detection and instance segmentation performance. We evaluate on both real and synthetic scan benchmark data for the new task, where we outperform state-of-the-art approaches by over 15 in mAP@0.5 on ScanNet, and over 18 in mAP@0.5 on SUNCG.


title: Hadronic and hadron-like physics of Dark Matter

abstract: The problems of simple elementary weakly interacting massive particles (WIMPs) appeal to extend the physical basis for nonbaryonic dark matter. Such extension involves more sophisticated dark matter candidates from physics beyond the Standard Model (BSM) of elementary particles. We discuss several models of dark matter, predicting new colored, hyper-colored or techni-colored particles and their accelerator and non-accelerator probes. The nontrivial properties of the proposed dark matter candidates can shed new light on the dark matter physics. They provide interesting solutions for the puzzles of direct and indirect dark matter search.


title: Knot reversal acts non-trivially on the concordance group of topologically slice knots

abstract: We construct an infinite family of topologically slice knots that are not smoothly concordant to their reverses. More precisely, if T denotes the concordance group of topologically slice knots and R is the involution of T induced by string reversal, then T/Fix(R) contains an infinitely generated free subgroup. The result remains true modulo the subgroup of T generated by knots with trivial Alexander polynomial.


title: A family of MCF solutions for the Heisenberg Group

abstract: The aim of this paper is to investigate the mean curvature flow soliton solutions on the Heisenberg group $\mathcal{H}$ when the initial data is a ruled surface by straight lines. We give a family of those solutions which are generated by $\mathfrak{Iso}_{0}(\mathcal{H})$ (the isometries of $\mathcal{H}$ for which the origin is a fix point). We conclude that the function which describe the motion of these surfaces under MCF, is always a linear affine function. As an application we proof that the Grim Reaper solution evolves from a ruled surface in $\mathcal{H}$. We also provide other examples.


title: CVR: A Continuously Variable Rate LDPC Decoder Using Parity Check Extension for Minimum Latency

abstract: This brief presents a novel IEEE 802.16e (WiMAX) based decoder that performs close to the 5G code but without the expensive hardware re-development cost. The design uses an extension of the existing WiMAX parity check code to reduce the processing latency and power consumption while keeping the decoder throughput at maximum. It achieves similar Frame Error Rate (FER) compared to 5G (0.1dB off), and most notably the error curves trend down like 5G instead flooring. At FER=10^-3 there is 0.1 dB gain in the FER code performance compared to WiMAX. An implementation of the design is a modified version of the existing fully-parallel WiMAX decoder that supports multi-rate codeword size and reduces latency by 33%. Additionally, for SNR greater than 3dB, decoding only the shorter code reduces the power consumption by 36%.


title: A Distributed Method for Fitting Laplacian Regularized Stratified Models

abstract: Stratified models are models that depend in an arbitrary way on a set of selected categorical features, and depend linearly on the other features. In a basic and traditional formulation a separate model is fit for each value of the categorical feature, using only the data that has the specific categorical value. To this formulation we add Laplacian regularization, which encourages the model parameters for neighboring categorical values to be similar. Laplacian regularization allows us to specify one or more weighted graphs on the stratification feature values. For example, stratifying over the days of the week, we can specify that the Sunday model parameter should be close to the Saturday and Monday model parameters. The regularization improves the performance of the model over the traditional stratified model, since the model for each value of the categorical `borrows strength' from its neighbors. In particular, it produces a model even for categorical values that did not appear in the training data set. We propose an efficient distributed method for fitting stratified models, based on the alternating direction method of multipliers (ADMM). When the fitting loss functions are convex, the stratified model fitting problem is convex, and our method computes the global minimizer of the loss plus regularization; in other cases it computes a local minimizer. The method is very efficient, and naturally scales to large data sets or numbers of stratified feature values. We illustrate our method with a variety of examples.


title: Self-intersection of the Fallback Stream in Tidal Disruption Events

abstract: We propose a semi-analytical model for the self-intersection of the fallback stream in tidal disruption events (TDEs). When the initial periapsis is less than about 15 gravitational radii, a large fraction of the shocked gas is unbound in the form of a collision-induced outflow (CIO). This is because large apsidal precession causes the stream to self-intersect near the local escape speed at radius much below the apocenter. The rest of the fallback gas is left in more tightly bound orbits and quickly joins the accretion flow. We propose that the CIO is responsible for reprocessing the hard emission from the accretion flow into the optical band. This picture naturally explains the large photospheric radius (or low blackbody temperature) and typical widths of the H and/or He emission lines seen in optical TDEs. We predict the CIO-reprocessed spectrum in the infrared to be L_\nu ~ \nu^{~0.5}, shallower than a blackbody. The partial sky coverage of the CIO also provides a unification of the diverse X-ray behaviors of optical TDEs. According to this picture, optical surveys filter out a large fraction of TDEs with low-mass blackholes due to lack of a reprocessing layer, and the volumetric rate of optical TDEs is nearly flat wrt. the blackhole mass for M < 10^7 solar masses. This filtering causes the optical TDE rate to be lower than the total rate by a factor of ~10 or more. When the CIO is decelerated by the ambient medium, radio emission at the level of that in ASASSN-14li may be produced, but the timescales and peak luminosities can be highly diverse. Finally, our method paves the way for global simulations of the disk formation process by injecting gas at the intersection point according to the prescribed velocity and density profiles.


title: Recurrent Embedding Aggregation Network for Video Face Recognition

abstract: Recurrent networks have been successful in analyzing temporal data and have been widely used for video analysis. However, for video face recognition, where the base CNNs trained on large-scale data already provide discriminative features, using Long Short-Term Memory (LSTM), a popular recurrent network, for feature learning could lead to overfitting and degrade the performance instead. We propose a Recurrent Embedding Aggregation Network (REAN) for set to set face recognition. Compared with LSTM, REAN is robust against overfitting because it only learns how to aggregate the pre-trained embeddings rather than learning representations from scratch. Compared with quality-aware aggregation methods, REAN can take advantage of the context information to circumvent the noise introduced by redundant video frames. Empirical results on three public domain video face recognition datasets, IJB-S, YTF, and PaSC show that the proposed REAN significantly outperforms naive CNN-LSTM structure and quality-aware aggregation methods.


title: Quantum noise limited nanoparticle detection with exposed-core fiber

abstract: Label-free biosensors are important tools for clinical diagnostics and for studying biology at the single molecule level. The development of optical label-free sensors has allowed extreme sensitivity, but can expose the biological sample to photodamage. Moreover, the fragility and complexity of these sensors can be prohibitive to applications. To overcome these problems, we develop a quantum noise limited exposed-core fiber sensor providing robust platform for label-free biosensing with a natural path toward microfluidic integration. We demonstrate the detection of single nanoparticles down to 25 nm in radius with optical intensities beneath known biophysical damage thresholds.


title: Best Match Graphs and Reconciliation of Gene Trees with Species Trees

abstract: A wide variety of problems in computational biology, most notably the assessment of orthology, are solved with the help of reciprocal best matches. Using an evolutionary definition of best matches that captures the intuition behind the concept we clarify rigorously the relationships between reciprocal best matches, orthology, and evolutionary events under the assumption of duplication/loss scenarios. We show that the orthology graph is a subgraph of the reciprocal best match graph (RBMG). We furthermore give conditions under which an RBMG that is a cograph identifies the correct orthlogy relation. Using computer simulations we find that most false positive orthology assignments can be identified as so-called good quartets -- and thus corrected -- in the absence of horizontal transfer. Horizontal transfer, however, may introduce also false-negative orthology assignments.


title: Blue-Green Emission from Epitaxial Yet Cation-Disordered ZnGeN$_{2-x}$O$_x$

abstract: ZnGeN$_2$ offers a low-cost alternative to InGaN with the potential for bandgap tuning to span the green gap using cation site ordering. The addition of oxygen on the anion site creates an additional degree of electronic tunability. Here, we investigate the structure and optoelectronic properties of an epitaxial ZnGeN$_{2-x}$O$_{x}$ thin film library grown by combinatorial co-sputtering on c-Al$_2$O$_3$. Samples exhibit X-ray diffraction patterns and X-ray pole figures characteristic of a wurtzite (cation-disordered) structure with the expected 6-fold in-plane symmetry. Transmission electron microscopy reveals a semi-coherent interface with periodic dislocations that relieve strain from the large lattice mismatch, and confirms the in-plane and out-of-plane crystallographic orientation. Room-temperature photoluminescence exhibits peaks between 2.4 and 2.8 eV which are consistent with a sharp absorption onset observed by UV-vis spectroscopy. These results demonstrate low-cost synthesis of optically active yet cation disordered ZnGeN$_{2-x}$O$_{x}$, indicating a path toward application as a blue-green emitter.


title: Discovering Common Change-Point Patterns in Functional Connectivity Across Subjects

abstract: This paper studies change-points in human brain functional connectivity (FC) and seeks patterns that are common across multiple subjects under identical external stimulus. FC relates to the similarity of fMRI responses across different brain regions when the brain is simply resting or performing a task. While the dynamic nature of FC is well accepted, this paper develops a formal statistical test for finding {\it change-points} in times series associated with FC. It represents short-term connectivity by a symmetric positive-definite matrix, and uses a Riemannian metric on this space to develop a graphical method for detecting change-points in a time series of such matrices. It also provides a graphical representation of estimated FC for stationary subintervals in between the detected change-points. Furthermore, it uses a temporal alignment of the test statistic, viewed as a real-valued function over time, to remove inter-subject variability and to discover common change-point patterns across subjects. This method is illustrated using data from Human Connectome Project (HCP) database for multiple subjects and tasks.


title: New type II Cepheids from VVV data towards the Galactic center

abstract: The Galactic center (GC) is the densest region of the Milky Way. Variability surveys towards the GC potentially provide the largest number of variable stars per square degree within the Galaxy. However, high stellar density is also a drawback due to blending. Moreover, the GC is affected by extreme reddening, therefore near infrared observations are needed. We plan to detect new variable stars towards the GC, focusing on type II Cepheids (T2Cs) which have the advantage of being brighter than RR Lyrae stars. We perform parallel Lomb-Scargle and Generalized Lomb-Scargle periodogram analysis of the $K_s$-band time series of the VISTA variables in the Via Lactea survey, to detect periodicities. We employ statistical parameters to clean our sample. We take account of periods, light amplitudes, distances, and proper motions to provide a classification of the candidate variables. We detected 1,019 periodic variable stars, of which 164 are T2Cs, 210 are Miras and 3 are classical Cepheids. We also found the first anomalous Cepheid in this region. We compare their photometric properties with overlapping catalogs and discuss their properties on the color-magnitude and Bailey diagrams. We present the most extensive catalog of T2Cs in the GC region to date. Offsets in E($J-K_s$) and in the reddening law cause very large ($\sim$1-2 kpc) uncertainties on distances in this region. We provide a catalog which will be the starting point for future spectroscopic surveys in the innermost regions of the Galaxy.


title: Polyhedral Properties of the Induced Cluster Subgraphs

abstract: A cluster graph is a graph whose every connected component is a complete graph. Given a simple undirected graph $G$, a subset of vertices inducing a cluster graph is called an independent union of cliques (IUC), and the IUC polytope associated with $G$ is defined as the convex hull of the incidence vectors of all IUCs in the graph. The {\sc Maximum IUC} problem, which is to find a maximum-cardinality IUC in a graph, finds applications in network-based data analysis. In this paper, we derive several families of facet-defining valid inequalities for the IUC polytope. We also give a complete description of this polytope for some special classes of graphs. We establish computational complexity of the separation problem for most of the considered families of valid inequalities and explore the effectiveness of employing the corresponding cutting planes in an integer (linear) programming framework for the {\sc Maximum IUC} problem through computational experiments.


title: StartupBR: Higher Education's Influence on Social Networks and Entrepreneurship in Brazil

abstract: Developing and middle-income countries increasingly empha-size higher education and entrepreneurship in their long-term develop-ment strategy. Our work focuses on the influence of higher education institutions (HEIs) on startup ecosystems in Brazil, an emerging economy. First, we describe regional variability in entrepreneurial network characteristics. Then we examine the influence of elite HEIs in economic hubs on entrepreneur networks. Second, we investigate the influence ofthe academic trajectories of startup founders, including their courses of study and HEIs of origin, on the fundraising capacity of startups. Given the growing capability of social media databases such as Crunchbase and LinkedIn to provide startup and individual-level data, we draw on computational methods to mine data for social network analysis. We find that HEI quality and the maturity of the ecosystem influence startup success. Our network analysis illustrates that elite HEIs have powerful influences on local entrepreneur ecosystems. Surprisingly, while the most nationally prestigious HEIs in the South and Southeast have the longest geographical reach, their network influence still remains local.


title: Flare-productive active regions

abstract: Strong solar flares and coronal mass ejections, here defined not only as the bursts of electromagnetic radiation but as the entire process in which magnetic energy is released through magnetic reconnection and plasma instability, emanate from active regions (ARs) in which high magnetic non-potentiality resides in a wide variety of forms. This review focuses on the formation and evolution of flare-productive ARs from both observational and theoretical points of view. Starting from a general introduction of the genesis of ARs and solar flares, we give an overview of the key observational features during the long-term evolution in the pre-flare state, the rapid changes in the magnetic field associated with the flare occurrence, and the physical mechanisms behind these phenomena. Our picture of flare-productive ARs is summarized as follows: subject to the turbulent convection, the rising magnetic flux in the interior deforms into a complex structure and gains high non-potentiality; as the flux appears on the surface, an AR with large free magnetic energy and helicity is built, which is represented by delta-sunspots, sheared polarity inversion lines, magnetic flux ropes, etc; the flare occurs when sufficient magnetic energy has accumulated, and the drastic coronal evolution affects magnetic fields even in the photosphere. We show that the improvement of observational instruments and modeling capabilities has significantly advanced our understanding in the last decades. Finally, we discuss the outstanding issues and future perspective and further broaden our scope to the possible applications of our knowledge to space-weather forecasting, extreme events in history, and corresponding stellar activities.


title: Arbitrary-time error suppression for Markovian adiabatic quantum computing using stabilizer subspace codes

abstract: Adiabatic quantum computing (AQC) can be protected against thermal excitations via an encoding into error detecting codes, supplemented with an energy penalty formed from a sum of commuting Hamiltonian terms. Earlier work showed that it is possible to suppress the initial thermally induced excitation out of the encoded ground state, in the case of local Markovian environments, by using an energy penalty strength that grows only logarithmically in the system size, at a fixed temperature. The question of whether this result applies beyond the initial time was left open. Here we answer this in the affirmative. We show that thermal excitations out of the encoded ground state can be suppressed at arbitrary times under the additional assumption that the total evolution time is polynomial in the system size. Thus, computational problems that can be solved efficiently using AQC in a closed system setting, can still be solved efficiently subject to coupling to a thermal environment. Our construction uses stabilizer subspace codes, which require at least $4$-local interactions to achieve this result.


title: The Origin of [CII] 158um Emission toward the HII Region Complex S235

abstract: Although the 2P3/2-2P1/2 transition of [CII] at 158um is known to be an excellent tracer of active star formation, we still do not have a complete understanding of where within star formation regions the emission originates. Here, we use SOFIA upGREAT observations of [CII] emission toward the HII region complex Sh2-235 (S235) to better understand in detail the origin of [CII] emission. We complement these data with a fully-sampled Green Bank Telescope radio recombination line map tracing the ionized hydrogen gas. About half of the total [CII] emission associated with S235 is spatially coincident with ionized hydrogen gas, although spectroscopic analysis shows little evidence that this emission is coming from the ionized hydrogen volume. Velocity-integrated [CII] intensity is strongly correlated with WISE 12um intensity across the entire complex, indicating that both trace ultra-violet radiation fields. The 22um and radio continuum intensities are only correlated with [CII] intensity in the ionized hydrogen portion of the S235 region and the correlations between the [CII] and molecular gas tracers are poor across the region. We find similar results for emission averaged over a sample of external galaxies, although the strength of the correlations is weaker. Therefore, although many tracers are correlated with the strength of [CII] emission, only WISE 12um emission is correlated on small-scales of the individual HII region S235 and also has a decent correlation at the scale of entire galaxies. Future studies of a larger sample of Galactic HII regions would help to determine whether these results are truly representative.


title: From Trigroups To Leibniz 3-Algebras

abstract: In this paper, we study the category of trigroups as a generalization of the notion of digroup [4] and analyze their relationship with 3-racks [1] and Leibniz 3-algebras [6]. Trigroups are essentially associative trioids in which there are bar-units and bar-inverses. We prove that 3-racks can be constructed by conjugating trigroups. We also prove that trigroups equipped with a smooth manifold structure produce Leibniz 3-algebras via their associated Lie 3-racks.


title: A Perturbative Approach to the Tunneling Phenomena

abstract: The double-well potential is a good example, where we can compute the splitting in the bound state energy of the system due to the tunneling effect with various methods, namely WKB or instanton calculations. All these methods are non-perturbative and there is a common belief that it is difficult to find the splitting in the energy due to the barrier penetration from a perturbative analysis. However, we will illustrate by explicit examples containing singular potentials (e.g., Dirac delta potentials supported by points and curves and their relativistic extensions)that it is possible to find the splitting in the bound state energies by developing some kind of perturbation method.


title: On finding the smallest happy numbers of any heights

abstract: This paper focuses on finding the smallest happy number for each height in any numerical base. Using the properties of the height, we deduce a recursive relationship between the smallest happy number and the height where the initial height is function of the numerical base. With the usage of the recursive relationship, we build an algorithm that exploits the properties of the height in order to find all of those smallest happy numbers with unknown height. However, with the modular arithmetic, we conclude on an equation that calculates the smallest happy numbers based on known heights for binary and ternary bases.


title: Simplified model and analysis of a pair of coupled thermo-optical MEMS oscillators

abstract: Motivated by the dynamics of micro-scale oscillators with thermo-optical feedback, a simplified third order model, capturing the key features of these oscillators is developed, where each oscillator consists of a displacement variable coupled to a temperature variable. Further, the dynamics of a pair of such oscillators coupled via a linear spring is analyzed. The analytical procedures used are the variational equation method and the two-variable expansion method. It is shown that the analytical results are in agreement with the results of numerical integration. The bifurcation structure of the system is revealed through a bifurcation diagram.


title: Study of $\eta$-meson Decays at KLOE/KLOE-2

abstract: Working at a $\phi$-factory, KLOE/KLOE-2 has collected a very large and clean sample of $\eta$ decays, which allow to perform a variety of precision measurements. This report presents the recent results on the $\eta$-meson decays from KLOE/KLOE-2 experiments, including the search of the $C$-violating decay $\eta\rightarrow\gamma\gamma\gamma$, the determination of the $CP$ asymmetry in $\eta\rightarrow\pi^+\pi^-e^+e^-$, the precision measurement of the $\eta\rightarrow\pi^+\pi^-\pi^0$ Dalitz plot distribution, the search for the $CP$-violating decay $\eta\rightarrow\pi^+\pi^-$ and the precision measurement of ${\cal O}(p^6)$ dominated decay $\eta\rightarrow\gamma\gamma\pi^0$.


title: The nature of the infrared counterpart and of the optical nebula associated with the Central Compact Object in Vela Jr.\

abstract: The X-ray source CXO\, J085201.4$-$461753in the few kyr-old Vela Jr. supernova remnant (G266.2-1.2) belongs to the peculiar class of isolated neutron stars dubbed "Central Compact Objects" (CCOs). At variance with the other CCOs, which are only detected in the X-rays, \cxo\ has been possibly detected also at other wavelengths. In particular, a candidate near-infrared counterpart (H = 21.6 $\pm$0.1) was detected by the Very Large Telescope (VLT) in addition to a 6\arcsec-wide nebula detected in H$\alpha$, interpreted as a velocity-driven bow-shock. Here, we present new near-infrared and optical VLT observations of the candidate counterpart to \cxo\ and its nebula to confirm the proposed associations. Moreover, we used archival Chandra observations to measure for the first time the CXO\, J085201.4$-$461753 proper motion. The comparison between the two sets of near-infrared VLT images, taken 10.56 year apart, do not show evidence of proper motion for the candidate counterpart to CXO\, J085201.4$-$461753, expected from its 4\arcmin\ offset from the SNR geometrical centre, with a $3\sigma$ upper limit of $\sim 10$ mas yr$^{-1}$. This is much smaller than the expected proper motion of $\sim$50--100 mas yr$^{-1}$, which, in turn, is below the $3\sigma$ upper limit of $\sim$ 300 mas yr$^{-1}$ that we obtained with Chandra. At the same time, VLT spectroscopy resolved the nebula optical emission, which is dominated by the NII doublet at 6548 and 6584 \AA\ and not by H$\alpha$. To summarise, we conclude that the proposed near-infrared counterpart is likely unassociated with \cxo\ and that the nebula is not a velocity-driven bow-shock.


title: Repeating Fast Radio Bursts from Magnetars with Low Magnetospheric Twist

abstract: We analyze the statistics of pulse arrival times in fast radio burst (FRB) 121102 and demonstrate that they are remarkably similar to statistics of magnetar high-energy short bursts. Motivated by this correspondence, we propose that repeating FRBs are generated during short bursts in the closed field line zone of magnetar magnetospheres via a pulsar-like emission mechanism. Crustal slippage events dislocate field line foot points, initiating intense particle acceleration and pair production, giving rise to coherent radio emission similar to that generated near pulsar polar caps. We argue that the energetics of FRB 121102 can be readily accounted for if the efficiency of the conversion of Poynting flux into coherent radio emission is $\sim10^{-4}-10^{-2}$, values consistent with empirical efficiencies of radio emission in pulsars and radio-loud magnetars. Such a mechanism could operate only in magnetars with preexisting low twist of the magnetosphere, so that the charge density in the closed zone is initially insufficient to screen the electric field provoked by the wiggling of magnetic field lines and is low enough to let $\sim 1$ GHz radio emission escape the magnetosphere, which can explain the absence of FRBs from known magnetars. The pair cascades crowd the closed flux tubes with plasma, screening the accelerating electric field, thus limiting the radio pulse duration to $\sim1$ ms. Within the framework of our model, the current dataset of the polarization angle variation in FRB 121102 suggests a magnetic obliquity $\alpha\lesssim40^\circ$ and viewing angle $\zeta$ with respect to the spin axis $\alpha<\zeta<180^\circ-\alpha$.


title: Ultrafast Polarization of an Electron Beam in an Intense Bi-chromatic Laser Field

abstract: Here, we demonstrate the radiative polarization of high-energy electron beams in collisions with ultrashort pulsed bi-chromatic laser fields. Employing a Boltzmann kinetic approach for the electron distribution allows us to simulate the beam polarization over a wide range of parameters and determine the optimum conditions for maximum radiative polarization. Those results are contrasted with a Monte-Carlo algorithm where photon emission and associated spin effects are treated fully quantum mechanically using spin-dependent photon emission rates. The latter method includes realistic focusing laser fields, which allows us to simulate a near-term experimentally feasible scenario of a 8 GeV electron beam scattering from a 1 PW laser pulse and provide a measurement that would verify the ultrafast radiative polarization in high-intensity laser pulses that we predict. Aspects of spin dependent radiation reaction are also discussed, with spin polarization leading to a measurable (5%) splitting of the energies of spin-up and spin-down electrons.


title: Phonon Lifetimes and Thermal Conductivity of the Molecular Crystal ${\alpha}$-RDX

abstract: The heat transfer properties of the organic molecular crystal ${\alpha}$-RDX were studied using three phonon-based thermal conductivity models. It was found that the widely used Peierls-Boltzmann model for thermal transport in crystalline materials breaks down for ${\alpha}$-RDX. We show this breakdown is due to a large degree of anharmonicity that leads to a dominance of diffusive-like carriers. Despite being developed for disordered systems, the Allen-Feldman theory for thermal conductivity actually gives the best description of thermal transport. This is likely because diffusive carriers contribute to over 95% of the thermal conductivity in ${\alpha}$-RDX. The dominance of diffusive carriers is larger than previously observed in other fully ordered crystalline systems. These results indicate than van-der Waals bonded organic crystalline solids conduct heat in a manner more akin to amorphous materials than simple atomic crystals.


title: Causal relationship between eWOM topics and profit of rural tourism at Japanese Roadside Stations "MICHINOEKI"

abstract: Affected by urbanization, centralization and the decrease of overall population, Japan has been making efforts to revitalize the rural areas across the country. One particular effort is to increase tourism to these rural areas via regional branding, using local farm products as tourist attractions across Japan. Particularly, a program subsidized by the government called Michinoeki, which stands for 'roadside station', was created 20 years ago and it strives to provide a safe and comfortable space for cultural interaction between road travelers and the local community, as well as offering refreshment, and relevant information to travelers. However, despite its importance in the revitalization of the Japanese economy, studies with newer technologies and methodologies are lacking. Using sales data from establishments in the Kyushu area of Japan, we used Support Vector to classify content from Twitter into relevant topics and studied their causal relationship to the sales for each establishment using LiNGAM, a linear non-gaussian acyclic model built for causal structure analysis, to perform an improved market analysis considering more than just correlation. Under the hypotheses stated by the LiNGAM model, we discovered a positive causal relationship between the number of tweets mentioning those establishments, specially mentioning deserts, a need for better access and traf^ic options, and a potentially untapped customer base in motorcycle biker groups.


title: Community Detection and Growth Potential Prediction from Patent Citation Networks

abstract: The scoring of patents is useful for technology management analysis. Therefore, a necessity of developing citation network clustering and prediction of future citations for practical patent scoring arises. In this paper, we propose a community detection method using the Node2vec. And in order to analyze growth potential we compare three ''time series analysis methods'', the Long Short-Term Memory (LSTM), ARIMA model, and Hawkes Process. The results of our experiments, we could find common technical points from those clusters by Node2vec. Furthermore, we found that the prediction accuracy of the ARIMA model was higher than that of other models.


title: Quantum Frequency Conversion of a Quantum Dot Single-Photon Source on a Nanophotonic Chip

abstract: Single self-assembled InAs/GaAs quantum dots are promising bright sources of indistinguishable photons for quantum information science. However, their distribution in emission wavelength, due to inhomogeneous broadening inherent to their growth, has limited the ability to create multiple identical sources. Quantum frequency conversion can overcome this issue, particularly if implemented using scalable chip-integrated technologies. Here, we report the first demonstration of quantum frequency conversion of a quantum dot single-photon source on a silicon nanophotonic chip. Single photons from a quantum dot in a micropillar cavity are shifted in wavelength with an on-chip conversion efficiency of approximately 12 %, limited by the linewidth of the quantum dot photons. The intensity autocorrelation function g(2)(tau) for the frequency-converted light is antibunched with g(2)(0) = 0.290 +/- 0.030, compared to the before-conversion value g(2)(0) = 0.080 +/- 0.003. We demonstrate the suitability of our frequency conversion interface as a resource for quantum dot sources by characterizing its effectiveness across a wide span of input wavelengths (840 nm to 980 nm), and its ability to achieve tunable wavelength shifts difficult to obtain by other approaches.


title: Truly Optimal Euclidean Spanners

abstract: Euclidean spanners are important geometric structures, having found numerous applications over the years. Cornerstone results in this area from the late 80s and early 90s state that for any $d$-dimensional $n$-point Euclidean space, there exists a $(1+\epsilon)$-spanner with $nO(\epsilon^{-d+1})$ edges and lightness $O(\epsilon^{-2d})$. Surprisingly, the fundamental question of whether or not these dependencies on $\epsilon$ and $d$ for small $d$ can be improved has remained elusive, even for $d = 2$. This question naturally arises in any application of Euclidean spanners where precision is a necessity. The state-of-the-art bounds $nO(\epsilon^{-d+1})$ and $O(\epsilon^{-2d})$ on the size and lightness of spanners are realized by the {\em greedy} spanner. In 2016, Filtser and Solomon proved that, in low dimensional spaces, the greedy spanner is near-optimal. The question of whether the greedy spanner is truly optimal remained open to date. The contribution of this paper is two-fold. We resolve these longstanding questions by nailing down the exact dependencies on $\epsilon$ and $d$ and showing that the greedy spanner is truly optimal. Specifically, for any $d= O(1), \epsilon = \Omega({n}^{-\frac{1}{d-1}})$: - We show that any $(1+\epsilon)$-spanner must have $n \Omega(\epsilon^{-d+1})$ edges, implying that the greedy (and other) spanners achieve the optimal size. - We show that any $(1+\epsilon)$-spanner must have lightness $\Omega(\epsilon^{-d})$, and then improve the upper bound on the lightness of the greedy spanner from $O(\epsilon^{-2d})$ to $O(\epsilon^{-d})$. We then complement our negative result for the size of spanners with a rather counterintuitive positive result: Steiner points lead to a quadratic improvement in the size of spanners! Our bound for the size of Steiner spanners is tight as well (up to lower-order terms).


title: Dynamic Mini-batch SGD for Elastic Distributed Training: Learning in the Limbo of Resources

abstract: With an increasing demand for training powers for deep learning algorithms and the rapid growth of computation resources in data centers, it is desirable to dynamically schedule different distributed deep learning tasks to maximize resource utilization and reduce cost. In this process, different tasks may receive varying numbers of machines at different time, a setting we call elastic distributed training. Despite the recent successes in large mini-batch distributed training, these methods are rarely tested in elastic distributed training environments and suffer degraded performance in our experiments, when we adjust the learning rate linearly immediately with respect to the batch size. One difficulty we observe is that the noise in the stochastic momentum estimation is accumulated over time and will have delayed effects when the batch size changes. We therefore propose to smoothly adjust the learning rate over time to alleviate the influence of the noisy momentum estimation. Our experiments on image classification, object detection and semantic segmentation have demonstrated that our proposed Dynamic SGD method achieves stabilized performance when varying the number of GPUs from 8 to 128. We also provide theoretical understanding on the optimality of linear learning rate scheduling and the effects of stochastic momentum.


title: Microscopic laws vs. Macroscopic laws: Perspectives from kinetic theory and hydrodynamics

abstract: Reductionism is a prevalent viewpoint in science according to which all physical phenomena can be understood from fundamental laws of physics. Anderson [Science, 177, 393 (1972)], Laughlin and Pines [PNAS, 97, 28 (2000)], and others have countered this viewpoint and argued in favour hierarchical structure of the universe and laws. In this paper we advance the latter perspective by showing that some of the complex flow properties derived using hydrodynamic equations (macroscopic laws) are very difficult, if not impossible, to describe in microscopic framework---kinetic theory. These properties include Kolmogorov's theory of turbulence, turbulence dissipation and diffusion, and dynamic pressure. We also provide several other examples of hierarchical description.

categories: physics.flu-dyn cond-mat.stat-mech nlin.CD physics.hist-ph quant-ph


title: Field-theoretical approach to the Casimir-like interaction in a one-dimensional Bose gas

abstract: We study the fluctuation-induced interaction between two impurities in a weakly-interacting one-dimensional Bose gas using the field theoretical approach. At separations between impurities shorter and of the order of the healing length of the system, the induced interaction has a classical origin and behaves exponentially. At separations longer than the healing length, the interaction is of a quantum origin and scales as the third power of the inverse distance. Finite temperature destroys the quasi-long-range order of the Bose gas and, accordingly, the induced interaction becomes exponentially suppressed beyond the thermal length. We obtain analytical expressions for the induced interaction at zero and finite temperature that are valid at arbitrary distances. We discuss experimental realizations as well as possible formation of bound states of two impurities, known as bipolarons.


title: The rarefied elliptic Bailey lemma and the Yang-Baxter equation

abstract: An elliptic Bailey lemma is formulated on the basis of the univariate rarefied elliptic beta integral. It leads to a generalized operator star-triangle relation and a new solution of the Yang-Baxter equation written as an integral operator with a rarefied elliptic hypergeometric kernel.


title: One-step deposition of nano-to-micron-scalable, high-quality digital image correlation patterns for high-strain in-situ multi-microscopy testing

abstract: Digital Image Correlation (DIC) is of vital importance in the field of experimental mechanics, yet, producing suitable DIC patterns for demanding in-situ mechanical tests remains challenging, especially for ultra-fine patterns, despite the large number of patterning techniques in the literature. Therefore, we propose a simple, flexible, one-step technique (only requiring a conventional deposition machine) to obtain scalable, high-quality, robust DIC patterns, suitable for a range of microscopic techniques, by deposition of a low melting temperature solder alloy in so-called 'island growth' mode, without elevating the substrate temperature. Proof of principle is shown by (near-)room-temperature deposition of InSn patterns, yielding highly dense, homogeneous DIC patterns over large areas with a feature size that can be tuned from as small as 10nm to 2um and with control over the feature shape and density by changing the deposition parameters. Pattern optimization, in terms of feature size, density, and contrast, is demonstrated for imaging with atomic force microscopy, scanning electron microscopy (SEM), optical microscopy and profilometry. Moreover, the performance of the InSn DIC patterns and their robustness to large deformations is validated in two challenging case studies of in-situ micro-mechanical testing: (i) self-adaptive isogeometric digital height correlation of optical surface height profiles of a coarse, bimodal InSn pattern providing microscopic 3D deformation fields (illustrated for delamination of aluminum interconnects on a polyimide substrate) and (ii) DIC on SEM images of a much finer InSn pattern allowing quantification of high strains near fracture locations (illustrated for rupture of a Fe foil). As such, the high controllability, performance and scalability of the DIC patterns offers a promising step towards more routine DIC-based in-situ micro-mechanical testing.


title: Quasi-Periodicity of Supermassive Binary Black Hole Accretion Approaching Merger

abstract: In this paper we continue the first ever study of magnetized mini-disks coupled to circumbinary accretion in a supermassive binary black hole (SMBBH) approaching merger reported in Bowen et al. 2018. We extend this simulation from 3 to 12 binary orbital periods. We find that relativistic SMBBH accretion acts as a resonant cavity, where quasi-periodic oscillations tied to the the frequency at which the black hole's orbital phase matches a non-linear $m=1$ density feature, or ``lump'', in the circumbinary accretion disk permeate the system. The rate of mass accretion onto each of the mini-disks around the black holes is modulated at the beat frequency between the binary frequency and the lump's mean orbital frequency, i.e., $\Omega_{\rm beat} = \Omega_{\rm bin} - \bar{\Omega}_{\rm lump}$, while the total mass accretion rate of this equal-mass binary is modulated at two different frequencies, $\gtrsim \bar{\Omega}_{\rm lump}$ and $\approx 2 \Omega_{\rm beat}$. The instantaneous rotation rate of the lump itself is also modulated at two frequencies close to the modulation frequencies of the total accretion rate, $\bar{\Omega}_{\rm lump}$ and $2 \Omega_{\rm beat}$. Because of the compact nature of the mini-disks in SMBBHs approaching merger, the inflow times within the mini-disks are comparable to the period on which their mass-supply varies, so that their masses---and the accretion rates they supply to their black holes---are strongly modulated at the same frequency. In essence, the azimuthal symmetry of the circumbinary disk is broken by the dynamics of orbits near a binary, and this $m=1$ asymmetry then drives quasi-periodic variation throughout the system, including both accretion and disk-feeding. In SMBBHs approaching merger, such time variability could introduce distinctive, increasingly rapid, fluctuations in their electromagnetic emission.


title: Measurement of the Spin Tune Using the Coherent Spin Motion of Polarized Proton in a Storage Ring

abstract: This paper reports the first spin tune measurement at high energies (24 GeV and 255 GeV) with a driven coherent spin motion. To maintain polarization in a polarized proton collider, it is important to know the spin tune of the polarized proton beam, which is defined as the number of full spin precessions per revolution. A nine-magnet spin flipper has demonstrated high spin-flip efficiency in the presence of two Siberian snakes [1]. The spin flipper drives a spin resonance with a given frequency (or tune) and strength. When the drive tune is close to the spin tune, the proton spin direction is not vertical anymore, but precesses around the vertical direction. By measuring the precession frequency of the horizontal component the spin tune can be precisely measured. A driven coherent spin motion and fast turn-by-turn polarization measurement are keys to the measurement. The vertical spin direction is restored after turning the spin flipper off and the polarization value is not affected by the measurement. The fact that this manipulation preserves the polarization makes it possible to measure the spin tune during operation of a high energy accelerator.


title: Traveling concentration pulses of bacteria in a generalized Keller-Segel model

abstract: We formulate the Smoluchowski equation for a run-and-tumble particle. It includes the mean tumble rate in a chemical field, for which we derive a Markovian response theory. Using a multipole expansion and a reaction-diffusion equation for the chemoattractant field, we derive a polarization extended model, which also includes the recently discovered angle bias. In the adiabatic limit we recover generalized Keller-Segel equations with diffusion and chemotactic coefficients that depend on the microscopic swimming parameters. By requiring the tumble rate to be positive, our model possesses an upper bound of the chemotactic drift velocity, which is no longer singular as in the original Keller-Segel equations. Using the Keller-Segel model, we present an extensive study of traveling bacterial concentration pulses demonstrating how speed, width, and height of the pulse depend on the microscopic parameters. Most importantly, we discover a maximum number of bacteria that the pulse can sustain - the maximum carrying capacity. Finally, we obtain a remarkably good match to experimental results on the traveling bacterial pulse. It does not require a second, signaling chemical field nor a singular chemotactic drift velocity.