title: The Impact of Renewable Energy Forecasts on Intraday Electricity Prices

abstract: In this paper we study the impact of errors in wind and solar power forecasts on intraday electricity prices. We develop a novel econometric model which is based on day-ahead wholesale auction curves data and errors in wind and solar power forecasts. The model shifts day-ahead supply curves to calculate intraday prices. We apply our model to the German EPEX SPOT SE data. Our model outperforms both linear and non-linear benchmarks. Our study allows us to conclude that errors in renewable energy forecasts exert a non-linear impact on intraday prices. We demonstrate that additional wind and solar power capacities induce non-linear changes in the intraday price volatility. Finally, we comment on economical and policy implications of our findings.


title: Gravitational wave energy budget in strongly supercooled phase transitions

abstract: We derive efficiency factors for the production of gravitational waves through bubble collisions and plasma-related sources in strong phase transitions, and find the conditions under which the bubble collisions can contribute significantly to the signal. We use lattice simulations to clarify the dependence of the colliding bubbles on their initial state. We illustrate our findings in two examples, the Standard Model with an extra $|H|^6$ interaction and a classically scale-invariant $U(1)_{\rm B-L}$ extension of the Standard Model. The contribution to the GW spectrum from bubble collisions is found to be negligible in the $|H|^6$ model, whereas it can play an important role in parts of the parameter space in the scale-invariant $U(1)_{\rm B-L}$ model. In both cases the sound-wave period is much shorter than a Hubble time, suggesting a significant amplification of the turbulence-sourced signal. We find, however, that the peak of the plasma-sourced spectrum is still produced by sound waves with the slower-falling turbulence contribution becoming important off-peak.


title: Constraints on a Gravitational Higgs Mechanism

abstract: We show that it is impossible to improve the high-energy behavior of the tree-level four-point amplitude of a massive spin-2 particle by including the exchange of any number of scalars and vectors in four spacetime dimensions. This constrains possible weakly coupled ultraviolet extensions of massive gravity, ruling out gravitational analogues of the Higgs mechanism based on particles with spins less than two. Any tree-level ultraviolet extension that is Lorentz invariant and unitary must involve additional massive particles with spins greater than or equal to two, as in Kaluza-Klein theories and string theory.


title: Jet grooming through reinforcement learning

abstract: We introduce a novel implementation of a reinforcement learning (RL) algorithm which is designed to find an optimal jet grooming strategy, a critical tool for collider experiments. The RL agent is trained with a reward function constructed to optimize the resulting jet properties, using both signal and background samples in a simultaneous multi-level training. We show that the grooming algorithm derived from the deep RL agent can match state-of-the-art techniques used at the Large Hadron Collider, resulting in improved mass resolution for boosted objects. Given a suitable reward function, the agent learns how to train a policy which optimally removes soft wide-angle radiation, allowing for a modular grooming technique that can be applied in a wide range of contexts. These results are accessible through the corresponding GroomRL framework.


title: Probing IGM accretion onto faint Ly{\alpha} emitters at z~2.8

abstract: (abridged) Observing the signature of accretion from the intergalactic medium (IGM) onto galaxies at z~3 requires the detection of faint (L<<L*) galaxies embedded in a filamentary matrix of low-density, metal-poor gas coherent over hundreds of kpc. We study the gaseous environment of three Lyman$\alpha$ emitters (LAEs) at z=2.7-2.8, found to be aligned in projection with a background QSO over ~250 kpc along the slit of a long-slit spectrum. The lack of detection of the LAEs in deep continuum images and the low inferred Ly$\alpha$ luminosities show the LAEs to be intrinsically faint, low-mass galaxies (L<0.1 L*, M_star< 0.1 M*). An echelle spectrum of the QSO reveals strong Ly-alpha absorption within $\pm200$ km/s from the LAEs. Our absorption line analysis leads to HI column densities in the range of log N(HI) =16-18. Associated absorption from ionic metal species CIV and SiIV constrains the gas metallicities to ~0.01 solar if the gas is optically thin, and possibly as low as ~0.001 solar if the gas is optically thick, assuming photoionization equilibrium. While the inferred metallicities are at least a factor of ten lower than expected metallicities in the interstellar medium (ISM) of these LAEs, they are consistent with the observed chemical enrichment level in the IGM at the same epoch. Total metal abundances and kinematic arguments suggest that these faint galaxies have not been able to affect the properties of their surrounding gas. The projected spatial alignment of the LAEs, together with the kinematic quiescence and correspondence between the LAEs and absorbing gas in velocity space suggests that these observations probe a possible filamentary structure. Taken together with the blue-dominant Ly$\alpha$ emission line profile of one of the objects, the evidence suggests that the absorbing gas is part of an accretion stream of low-metallicity gas in the IGM.


title: On the ratio of lapses in bimetric relativity

abstract: The two lapse functions in the Hassan-Rosen bimetric theory are not independent. Without knowing the relation between them, one cannot evolve the equations in the 3+1 formalism. This work computes the ratio of lapses for the spherically symmetric case, which is a prerequisite for numerical bimetric relativity.


title: On tensionless string field theory in AdS$_3$

abstract: We report on progress in formulating a field theory of tensionless strings in $AdS_3$, starting from the dual large-$N$ symmetric orbifold CFT. We propose a set of field equations which are gauge invariant under the higher spin algebra of the theory, the `Higher Spin Square'. The massless higher spin sector is captured by a Chern-Simons gauge field, while the matter sector is described by unfolded equations similar to those appearing in Vasiliev theory. Our equations incorporate the full perturbative spectrum of the theory, including states coming from the twisted sectors, and capture some of the interactions fixed by gauge invariance. We also discuss the spectrum of the bulk theory and explain how linearization around $AdS_3$ gives rise to the expected set of decoupled wave equations. Our results can be generalized to describe bulk duals of other large-$N$ symmetric orbifolds.


title: Search for high-energy neutrinos from binary neutron star mergers

abstract: To search for transient astrophysical neutrino sources, IceCube's optical and X-ray follow-up program is triggered by two or more neutrino candidates arriving from a similar direction within 100s. However, the rate of such neutrino multiplets was found to be consistent with the expected background of chance coincidences, such that the data does not provide indications for the existence of short-lived transient neutrino sources. Upper limits on the neutrino flux of transient source populations are presented in Aartsen et al. (2019) and we show here how these limits apply to the predicted neutrino emission from binary neutron star mergers.


title: The GREATS H$\beta$+[OIII] Luminosity Function and Galaxy Properties at $\mathbf{z\sim8}$: Walking the Way of JWST

abstract: The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H$\beta$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) survey and from HST, we perform spectral energy distribution (SED) fitting, using synthetic SEDs from a large grid of photoionization models. The deep Spitzer/IRAC data combined with our models exploring a large parameter space enables to constrain the [OIII]+H$\beta$ fluxes and equivalent widths for our sample, as well as the average physical properties of $z\sim8$ galaxies, such as the ionizing photon production efficiency with $\log(\xi_\mathrm{ion}/\mathrm{erg}^{-1}\hspace{1mm}\mathrm{Hz})\geq25.77$. We find a relatively tight correlation between the [OIII]+H$\beta$ and UV luminosity, which we use to derive for the first time the [OIII]+H$\beta$ luminosity function (LF) at $z\sim8$. The $z\sim8$ [OIII]+H$\beta$ LF is higher at all luminosities compared to lower redshift, as opposed to the UV LF, due to an increase of the [OIII]+H$\beta$ luminosity at a given UV luminosity from $z\sim3$ to $z\sim8$. Finally, using the [OIII]+H$\beta$ LF, we make predictions for JWST/NIRSpec number counts of $z\sim8$ galaxies. We find that the current wide-area extragalactic legacy fields are too shallow to use JWST at maximal efficiency for $z\sim8$ spectroscopy even at 1hr depth and JWST pre-imaging to $\gtrsim30$ mag will be required.


title: Differentiable Programming Tensor Networks

abstract: Differentiable programming is a fresh programming paradigm which composes parameterized algorithmic components and trains them using automatic differentiation (AD). The concept emerges from deep learning but is not only limited to training neural networks. We present theory and practice of programming tensor network algorithms in a fully differentiable way. By formulating the tensor network algorithm as a computation graph, one can compute higher order derivatives of the program accurately and efficiently using AD. We present essential techniques to differentiate through the tensor networks contractions, including stable AD for tensor decomposition and efficient backpropagation through fixed point iterations. As a demonstration, we compute the specific heat of the Ising model directly by taking the second order derivative of the free energy obtained in the tensor renormalization group calculation. Next, we perform gradient based variational optimization of infinite projected entangled pair states for quantum antiferromagnetic Heisenberg model and obtain start-of-the-art variational energy and magnetization with moderate efforts. Differentiable programming removes laborious human efforts in deriving and implementing analytical gradients for tensor network programs, which opens the door to more innovations in tensor network algorithms and applications.


title: Cosmic distance determination from photometric redshift samples using BAO peaks only

abstract: The galaxy distributions along the line-of-sight are significantly contaminated by the uncertainty on redshift measurements obtained through multiband photometry, which makes it difficult to get cosmic distance information measured from baryon acoustic oscillations, or growth functions probed by redshift distortions. We investigate the propagation of the uncertainties into large scale clustering by exploiting all known estimators, and propose the wedge approach as a promising analysis tool to extract cosmic distance information still remaining in the photometric galaxy samples. We test our method using simulated galaxy maps with photometric uncertainties of $\sigma_{0} =\left(0.01, 0.02, 0.03\right)$. The measured anisotropy correlation function $\xi$ is binned into the radial direction of $s$ and the angular direction of $\mu$, and the variations of $\xi(s,\mu)$ with perpendicular and radial cosmic distance measures of $D_A$ and $H^{-1}$ are theoretically estimated by an improved RSD model. Although the radial cosmic distance $H^{-1}$ is unable to be probed from any of the three photometric galaxy samples, the perpendicular component of $D_A$ is verified to be accurately measured even after the full marginalisation of $H^{-1}$. We measure $D_A$ with approximately 6% precision which is nearly equivalent to what we can expect from spectroscopic DR12 CMASS galaxy samples.


title: Quench, thermalization and residual entropy across a non-Fermi liquid to Fermi liquid transition

abstract: We study the thermalization, after sudden and slow quenches, of an interacting model having a quantum phase transition from a Sachdev-Ye-Kitaev (SYK) non-Fermi liquid (NFL) to a Fermi liquid (FL). The model has SYK fermions coupled to non-interacting lead fermions and can be realized in a graphene flake connected to external leads. After a sudden quench to the NFL, a thermal state is reached rapidly via collapse-revival oscillations of the quasiparticle residue of the lead fermions. In contrast, the quench to the FL, across the NFL-FL transition, leads to multiple prethermal regimes and much slower thermalization. In the slow quench performed over a time $\tau$, we find that the excitation energy generated has a remarkable intermediate-$\tau$ non-analytic power-law dependence, $\tau^{-\eta}$ with $\eta<1$, which seemingly masks the dynamical manifestation of the initial residual entropy of the SYK fermions. The power-law scaling is expected to eventually break down for $\tau\to\infty$, signaling a violation of adiabaticity, due to the residual entropy present in the SYK fermions.


title: Anti-Turing Machine

abstract: The invention of CPU-centric computing paradigm was incredible breakthrough of computer science that revolutionized our everyday life dramatically. However, the CPU- centric paradigm is based on the Turing machine concept and, as a result, expensive and power-hungry data transferring between the memory and CPU core is inevitable operation. Anti-Turing machine paradigm can be based on two fundamental principles: (1) data-centric computing, and (2) decentralized computing. Anti-Turing machine is able to execute a special type of programs. The commands of such program have to be addressed to the 2D or 3D persistent memory space is able to process data in-place. This program should not define the position or structure of data but it has to define the goal of data processing activity. Generally speaking, it needs to consider the whole memory space like the data transformation space. But the data placement, particular algorithm implementation, and strategy of algorithm execution are out of scope of the program.


title: Carrollian Physics at the Black Hole Horizon

abstract: We show that the geometry of a black hole horizon can be described as a Carrollian geometry emerging from an ultra-relativistic limit where the near-horizon radial coordinate plays the role of a virtual velocity of light tending to zero. We prove that the laws governing the dynamics of a black hole horizon, the null Raychaudhuri and Damour equations, are Carrollian conservation laws obtained by taking the ultra-relativistic limit of the conservation of an energy-momentum tensor; we also discuss their physical interpretation. We show that the vector fields preserving the Carrollian geometry of the horizon, dubbed Carrollian Killing vectors, include BMS-like supertranslations and superrotations and that they have non-trivial associated conserved charges on the horizon. In particular, we build a generalization of the angular momentum to the case of non-stationary black holes. Finally, we discuss the relation of these conserved quantities to the infinite tower of charges of the covariant phase space formalism.


title: Entanglement theory in distributed quantum information processing

abstract: Distributed quantum information processing is a promising platform for scaling up quantum information processing, where small- and intermediate-scale quantum devices are connected by a network of quantum channels for communicating quantum information, so as to cooperate in achieving larger-scale information processing. In such distributed settings, entangled states shared among the multiple devices serve as a resource for achieving nonlocal information processing tasks by local operations and classical communication (LOCC), where transformations of multipartite entangled states play central roles. This thesis analyzes properties of quantum entanglement in these small- and intermediate-scale settings and multipartite settings. The first part of this thesis investigates a communication task, quantum state merging, on the small and intermediate scales. The second part of this thesis analyzes multipartite entanglement in distributed quantum information processing. These analyses clarify fundamental limitations and potential applications of distributed quantum information processing to characterize properties of quantum entanglement in the small- and intermediate-scale settings and multipartite settings, providing a paradigm for investigating multipartite entanglement in distributed quantum information processing over networks beyond the state convertibility under LOCC.


title: How to Optimally Constrain Galaxy Assembly Bias: Supplement Projected Correlation Functions with Count-in-cells Statistics

abstract: Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $\Delta \Sigma(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+\Delta \Sigma(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $\Delta \Sigma(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection.


title: A generalization of the Minkowski distance and a new definition of the ellipse

abstract: In this paper, we generalize the Minkowski distance by defining a new distance function in n-dimensional space, and we show that this function determines also a metric family as the Minkowski distance. Then, we consider three special cases of this family, which generalize the taxicab, Euclidean and maximum metrics respectively, and finally we determine circles of them with their some properties in the real plane. While we determine some properties of circles of the generalized Minkowski distance, we also discover a new definition for the ellipse.


title: A Hybrid Approach to Persistent Coverage in Stochastic Environments

abstract: This paper considers the persistent coverage of a 2-D manifold that has been embedded in 3-D space. The manifold is subject to continual impact by intruders which travel at constant velocities along arbitrarily oriented straight-line trajectories. The trajectories of intruders are estimated online with an extended Kalman filter and their predicted impact points contribute normally distributed decay terms to the coverage level. A formal hybrid control strategy is presented that allows for power-constrained 3-D free-flyer agents to persistently monitor the domain, track and intercept intruders, and periodically deploy from and return to a single charging station on the manifold. Guarantees on intruder interception with respect to agent power lifespans are formally proven. The efficacy of the algorithm is demonstrated through simulation.


title: Eccentric binary black hole mergers in globular clusters hosting intermediate-mass black holes

abstract: Globular clusters (GCs) may harbour intermediate-mass black holes (IMBHs) at their centres. In these dynamically active environments stellar-mass black holes (SBHs) sink to the center soon after formation, due to dynamical friction and start interacting among themselves and with the central IMBH. Likely, some of the SBHs will form bound systems with the IMBH. A fraction of those will be triple systems composed of binary SBHs and the IMBH acting as a third distant perturber. If the SBH binary orbit is sufficiently inclined it can develop Lidov-Kozai (LK) oscillations, which can drive the system to high eccentricities and eventually to a merger due to gravitational wave (GW) emission on short timescales. In this work, we focus on the dynamics of the IMBH-SBH-SBH triples and illustrate that these systems can be possible sources of GWs. A distinctive signature of this scenario is that a considerable fraction of these mergers are highly eccentric when entering the LIGO band ($10$ Hz). Assuming that $\sim 20\%$ of GCs host IMBHs and a GC density in the range $n_{_{\rm GC}}=0.32$-$2.31\,\mathrm{Mpc}^{-3}$, we have estimated a rate $\Gamma=0.06$-$0.46\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ of these events. This suggests that dynamically-driven binary SBH mergers in this scenario could contribute to the merger events observed by LIGO/VIRGO. Full $N$-body simulations of GCs harbouring IMBHs are highly desirable to give a more precise constrain on this scenario.


title: Spectral instability of the peaked periodic wave in the reduced Ostrovsky equation

abstract: We show that the peaked periodic traveling wave of the reduced Ostrovsky equations with quadratic and cubic nonlinearity is spectrally unstable in the space of square integrable periodic functions with zero mean and the same period. The main novelty is that we discover a new instability phenomenon: the instability of the peaked periodic waves is induced by spectrum of a linearized operator which completely covers a closed vertical strip of the complex plane.

categories: math.AP math-ph math.DS math.MP nlin.PS nlin.SI


title: Simple non-perturbative resummation schemes beyond mean-field II: thermodynamics of scalar $\phi^4$ theory in 1+1 dimensions at arbitrary coupling

abstract: Recently, non-perturbative approximate solutions were presented that go beyond the well-known mean-field resummation. In this work, these non-perturbative approximations are used to calculate finite temperature equilibrium properties for scalar $\phi^4$ theory in two dimensions such as the pressure, entropy density and speed of sound. Unlike traditional approaches, it is found that results are well-behaved for arbitrary temperature/coupling strength, are independent of the choice of the renormalization scale $\bar\mu^2$, and are apparently converging as the resummation level is increased. Results also suggest the presence of a possible analytic cross-over from the high-temperature to the low-temperature regime based on the change in the thermal entropy density.


title: Capsule Networks with Max-Min Normalization

abstract: Capsule Networks (CapsNet) use the Softmax function to convert the logits of the routing coefficients into a set of normalized values that signify the assignment probabilities between capsules in adjacent layers. We show that the use of Softmax prevents capsule layers from forming optimal couplings between lower and higher-level capsules. Softmax constrains the dynamic range of the routing coefficients and leads to probabilities that remain mostly uniform after several routing iterations. Instead, we propose the use of Max-Min normalization. Max-Min performs a scale-invariant normalization of the logits that allows each lower-level capsule to take on an independent value, constrained only by the bounds of normalization. Max-Min provides consistent improvement in test accuracy across five datasets and allows more routing iterations without a decrease in network performance. A single CapsNet trained using Max-Min achieves an improved test error of 0.20% on the MNIST dataset. With a simple 3-model majority vote, we achieve a test error of 0.17% on MNIST.


title: Yu-Shiba-Rusinov states in the charge-density modulated superconductor NbSe2

abstract: NbSe$_2$ is a remarkable superconductor in which charge-density order coexists with pairing correlations at low temperatures. Here, we study the interplay of magnetic adatoms and their Yu-Shiba-Rusinov (YSR) bound states with the charge density order. Exploiting the incommensurate nature of the charge-density wave (CDW), our measurements provide a thorough picture of how the CDW affects both the energies and the wavefunctions of the YSR states. Key features of the dependence of the YSR states on adsorption site relative to the CDW are explained by model calculations. Several properties make NbSe$_2$ a promising substrate for realizing topological nanostructures. Our results will be important in designing such systems.


title: Imaging through nonlinear metalens using second harmonic generation

abstract: The abrupt phase change of light at metasurfaces provides high flexibility in wave manipulation without the need of accumulation of propagating phase through dispersive materials. In the linear optical regime, one important application field of metasurfaces is imaging by planar metalenses, which enables device miniaturization and aberration correction compared to conventional optical microlens systems. With the incorporation of nonlinear responses into passive metasurfaces, optical functionalities of metalenses are anticipated to be further enriched, leading to completely new applications areas. Here, we demonstrate imaging with nonlinear metalenses that combine the function of an ultrathin planar lens with simultaneous frequency conversion. With such nonlinear metalenses, we experimentally demonstrate imaging of objects with near infrared light while the image appears in the second harmonic signal of visible frequency range. Furthermore, the functionality of these nonlinear metalenses can be modified by switching the handedness of the circularly polarized fundamental wave, leading to either real or virtual nonlinear image formation. Nonlinear metalenses do not only enable infrared light imaging through a visible detector but also have the ability to modulate nonlinear optical responses through an ultrathin metasurface device while the fundamental wave remains unaffected, which offers the capability of nonlinear information processing with novel optoelectronic devices.


title: Near-Infrared Scanning Cavity Ringdown for Optical Loss Characterization of Supermirrors

abstract: A cavity ringdown system for probing the spatial variation of optical loss across high-reflectivity mirrors is described. This system is employed to examine substrate-transferred crystalline supermirrors and quantify the effect of manufacturing process imperfections. Excellent agreement is observed between the ringdown-generated spatial measurements and differential interference contrast microscopy images. A 2 mm diameter ringdown scan in the center of a crystalline supermirror reveals highly uniform coating properties with excess loss variations below 1 ppm.


title: Kelvin-Helmholtz Instability in Self-Gravitating Streams

abstract: Self-gravitating gaseous filaments exist on many astrophysical scales, from sub-pc filaments in the interstellar medium to Mpc scale streams feeding galaxies from the cosmic web. These filaments are often subject to Kelvin-Helmotz Instability (KHI) due to shearing against a confining background medium. We study the nonlinear evolution of KHI in pressure-confined self-gravitating gas streams initially in hydrostatic equilibrium, using analytic models and hydrodynamic simulations, not including radiative cooling. We derive a critical line-mass as a function of the stream Mach number and density contrast with respect to the background, $\mu_{cr}(M_b,\delta_c)\le 1$, where $\mu=1$ is normalized to the maximal line mass for which initial hydrostatic equilibrium is possible. For $\mu<\mu_{cr}$, KHI dominates the stream evolution. A turbulent shear layer expands into the background and leads to stream deceleration at a similar rate to the non-gravitating case. However, with gravity, penetration of the shear layer into the stream is halted at roughly half the initial stream radius by stabilizing buoyancy forces, significantly delaying total stream disruption. Streams with $\mu_{cr}<\mu\le 1$ fragment and form round, long-lived clumps by gravitational instability (GI), with typical separations roughly 8 times the stream radius, similar to the case without KHI. When KHI is still somewhat effective, these clumps are below the spherical Jeans mass and are partially confined by external pressure, but they approach the Jeans mass as $\mu\rightarrow 1$ and GI dominates. We discuss potential applications of our results to filaments in the ISM and dense streams feeding galaxies at high redshift.


title: Exceeding the Sauter-Schwinger limit of pair production with a quantum gas

abstract: We quantum-simulated particle-antiparticle pair production with a bosonic quantum gas in an optical lattice by emulating the requisite 1d Dirac equation and uniform electric field. We emulated field strengths far in excess of Sauter-Schwinger's limit for pair production in quantum electrodynamics, and therefore readily produced particles from "the Dirac vacuum" in quantitative agreement with theory. The observed process is equivalently described by Landau-Zener tunneling familiar in the atomic physics context.


title: Data Augmentation for Bayesian Deep Learning

abstract: Deep Learning (DL) methods have emerged as one of the most powerful tools for functional approximation and prediction. While the representation properties of DL have been well studied, uncertainty quantification remains challenging and largely unexplored. Data augmentation techniques are a natural approach to provide uncertainty quantification and to incorporate stochastic Monte Carlo search into stochastic gradient descent (SGD) methods. The purpose of our paper is to show that training DL architectures with data augmentation leads to efficiency gains. We use the theory of scale mixtures of normals to derive data augmentation strategies for deep learning. This allows variants of the expectation-maximization and MCMC algorithms to be brought to bear on these high dimensional nonlinear deep learning models. To demonstrate our methodology, we develop data augmentation algorithms for a variety of commonly used activation functions: logit, ReLU, leaky ReLU and SVM. Our methodology is compared to traditional stochastic gradient descent with back-propagation. Our optimization procedure leads to a version of iteratively re-weighted least squares and can be implemented at scale with accelerated linear algebra methods providing substantial improvement in speed. We illustrate our methodology on a number of standard datasets. Finally, we conclude with directions for future research.


title: Physics-Aware Neural Networks for Distribution System State Estimation

abstract: The distribution system state estimation problem seeks to determine the network state from available measurements. Widely used Gauss-Newton approaches are very sensitive to the initialization and often not suitable for real-time estimation. Learning approaches are very promising for real-time estimation, as they shift the computational burden to an offline training stage. Prior machine learning approaches to power system state estimation have been electrical model-agnostic, in that they did not exploit the topology and physical laws governing the power grid to design the architecture of the learning model. In this paper, we propose a novel learning model that utilizes the structure of the power grid. The proposed neural network architecture reduces the number of coefficients needed to parameterize the mapping from the measurements to the network state by exploiting the separability of the estimation problem. This prevents overfitting and reduces the complexity of the training stage. We also propose a greedy algorithm for phasor measuring units placement that aims at minimizing the complexity of the neural network required for realizing the state estimation mapping. Simulation results show superior performance of the proposed method over the Gauss-Newton approach.


title: Vector-scalar mixing to all orders, for an arbitrary gauge model in the generic linear gauge

abstract: I give explicit fromulae for full propagators of vector and scalar fields in a generic spin-1 gauge model quantized in an arbitrary linear covariant gauge. The propagators, expressed in terms of all-order one-particle-irreducible correlation functions, have a remarkably simple form because of constraints originating from Slavnov-Taylor identities of Becchi-Rouet-Stora symmetry. I also determine the behavior of the propagators in the neighborhood of the poles, and give a simple prescription for the coefficients that generalize (to the case with an arbitrary vector-scalar mixing) the standard $\sqrt{\mathcal{Z}}$ factors of Lehmann, Symanzik and Zimmermann. So obtained generalized $\sqrt{\mathcal{Z}}$ factors, are indispensable to the correct extraction of physical amplitudes from the amputated correlation functions in the presence of mixing. The standard $R_\xi$ guauges form a particularly important subclass of gauges considered in this paper. While the tree-level vector-scalar mixing is, by construction, absent in $R_\xi$ gauges, it unavoidably reappears at higher orders. Therefore the prescription for the generalized $\sqrt{\mathcal{Z}}$ factors given in this paper is directly relevant for the extraction of amplitudes in $R_\xi$ gauges.


title: Axonal Conduction Velocity Impacts Neuronal Network Oscillations

abstract: Increasing experimental evidence suggests that axonal action potential conduction velocity is a highly adaptive parameter in the adult central nervous system. Yet, the effects of this newfound plasticity on global brain dynamics is poorly understood. In this work, we analyzed oscillations in biologically plausible neuronal networks with different conduction velocity distributions. Changes of 1-2 (ms) in network mean signal transmission time resulted in substantial network oscillation frequency changes ranging in 0-120 (Hz). Our results suggest that changes in axonal conduction velocity may significantly affect both the frequency and synchrony of brain rhythms, which have well established connections to learning, memory, and other cognitive processes.


title: Large-$N_c$ and Large-$N_F$ Limits of SU($N_c$) Gauge Theories with Fermions in Different Representations

abstract: We present calculations of certain limits of scheme-independent series expansions for the anomalous dimensions of gauge-invariant fermion bilinear operators and for the derivative of the beta function at an infrared fixed point in SU($N_c$) gauge theories with fermions transforming according to two different representations. We first study a theory with $N_f$ fermions in the fundamental representation and $N_{f'}$ fermions in the adjoint or symmetric or antisymmetric rank-2 tensor representation, in the limit $N_c \to \infty$, $N_f \to \infty$ with $N_f/N_c$ fixed and finite. We then study the $N_c \to \infty$ limit of a theory with fermions in the adjoint and rank-2 symmetric or antisymmetric tensor representations.


title: Compliance Shaping for Control of Strength Amplification Exoskeletons with Elastic Cuffs

abstract: Exoskeletons which amplify the strength of their operators can enable heavy-duty manipulation of unknown objects. However, this type of behavior is difficult to accomplish; it requires the exoskeleton to sense and amplify the operator's interaction forces while remaining stable. But, the goals of amplification and robust stability when connected to the operator fundamentally conflict. As a solution, we introduce a design with a spring in series with the force sensitive cuff. This allows us to design an exoskeleton compliance behavior which is nominally passive, even with high amplification ratios. In practice, time delay and discrete time filters prevent our strategy from actually achieving passivity, but the designed compliance still makes the exoskeleton more robust to spring-like human behaviors. Our exoskeleton is actuated by a series elastic actuator (SEA), which introduces another spring into the system. We show that shaping the cuff compliance for the exoskeleton can be made into approximately the same problem as shaping the spring compliance of an SEA. We therefore introduce a feedback controller and gain tuning method which takes advantage of an existing compliance shaping technique for SEAs. We call our strategy the "double compliance shaping" method. With large amplification ratios, this controller tends to amplify nonlinear transmission friction effects, so we additionally propose a "transmission disturbance observer" to mitigate this drawback. Our methods are validated on a single-degree-of-freedom elbow exoskeleton.


title: Broadband Phonon Scattering in PbTe-based Materials Driven Near the Ferroelectric Phase Transition by Strain or Alloying

abstract: The major obstacle in the design of materials with low lattice thermal conductivity is the difficulty in efficiently scattering phonons across the entire frequency spectrum. Using first principles calculations, we show that driving PbTe materials to the brink of the ferroelectric phase transition could be a powerful strategy to solve this problem. We illustrate this concept by applying tensile [001] strain to PbTe and its alloys with another rock-salt IV-VI material, PbSe; and by alloying PbTe with a rhombohedral IV-VI material, GeTe. This induces extremely soft optical modes at the zone center, which increase anharmonic acoustic-optical coupling and decrease phonon lifetimes at all frequencies. We predict that PbTe, Pb(Se,Te) and (Pb,Ge)Te alloys driven close to the phase transition in the described manner will have considerably lower lattice thermal conductivity than that of PbTe (by a factor of 2-3). The proposed concept may open new opportunities for the development of more efficient thermoelectric materials.


title: Violation of lepton number in 3 units

abstract: The number of leptons may or may not be a conserved quantity. The Standard Model predicts that it is (in perturbative processes), but there is the well known possibility that new physics violates lepton number in one or two units. The first case ($\Delta L=1$) is associated to proton decay into mesons plus a lepton or an anti-lepton, while the second one ($\Delta L=2$) is usually associated to Majorana neutrino masses and neutrinoless double beta decay. It is also conceivable that leptons can only be created or destroyed in groups of three ($\Delta L=3$). Colliders and proton decay experiments can explore this possibility.


title: Stochastic phase-cohesiveness of discrete-time Kuramoto oscillators in a frequency-dependent tree network

abstract: This paper presents the notion of stochastic phase-cohesiveness based on the concept of recurrent Markov chains and studies the conditions under which a discrete-time stochastic Kuramoto model is phase-cohesive. It is assumed that the exogenous frequencies of the oscillators are combined with random variables representing uncertainties. A bidirectional tree network is considered such that each oscillator is coupled to its neighbors with a coupling law which depends on its own noisy exogenous frequency. In addition, an undirected tree network is studied. For both cases, a sufficient condition for the common coupling strength and a necessary condition for the sampling-period are derived such that the stochastic phase-cohesiveness is achieved. The analysis is performed within the stochastic systems framework and validated by means of numerical simulations.


title: Supernova Remnants in High Definition

abstract: Supernova remnants (SNRs) offer the means to study SN explosions, dynamics, and shocks at sub-parsec scales. X-ray observations probe the hot metals synthesized in the explosion and the TeV electrons accelerated by the shocks, and thus they are key to test recent, high-fidelity three-dimensional SN simulations. In this white paper, we discuss the major advances possible with X-ray spectro-imaging at arcsecond scales, with a few eV spectral resolution and a large effective area. These capabilities would revolutionize SN science, offering a three-dimensional view of metals synthesized in explosions and enabling population studies of SNRs in Local Group galaxies. Moreover, this future X-ray mission could detect faint, narrow synchrotron filaments and shock precursors that will constrain the diffusive shock acceleration process.


title: CONICOPF: Conic relaxations for AC optimal power flow computations

abstract: Computational speed and global optimality are key needs for practical algorithms for the optimal power flow problem. Two convex relaxations offer a favorable trade-off between the standard second-order cone and the standard semidefinite relaxations for large-scale meshed networks in terms of optimality gap and computation time: the tight-and-cheap relaxation (TCR) and the quadratic convex relaxation (QCR). We compare these relaxations on 60 PGLib-OPF test cases with up to 1,354 buses under three operating conditions and show that TCR dominates QCR on all 20 typical (TYP) test cases, on 18 out of 20 active power increase (API), and 12 out of 20 small angle difference (SAD). Selected state-of-the-art conic relaxations are implemented in the new MATLAB-based package CONICOPF available on GitHub.


title: Identification and Estimation of a Partially Linear Regression Model using Network Data

abstract: I study a regression model in which one covariate is an unknown function of a latent driver of link formation in a network. Rather than specify and fit a parametric network formation model, I introduce a new method based on matching pairs of agents with similar columns of the squared adjacency matrix, the ijth entry of which contains the number of other agents linked to both agents i and j. The intuition behind this approach is that for a large class of network formation models the columns of the squared adjacency matrix characterize all of the identifiable information about individual linking behavior. In this paper, I describe the model, formalize this intuition, and provide consistent estimators for the parameters of the regression model. Auerbach (2021) considers inference and an application to network peer effects.


title: Boundedness of a class of discretized reaction-diffusion systems

abstract: Although the spatially continuous version of the reaction-diffusion equation has been well studied, in some instances a spatially-discretized representation provides a more realistic approximation of biological processes. Indeed, mathematically the discretized and continuous systems can lead to different predictions of biological dynamics. It is well known in the continuous case that the incorporation of diffusion can cause diffusion-driven blow-up with respect to the $L^{\infty}$ norm. However, this does not imply diffusion-driven blow-up will occur in the discretized version of the system. For example, in a continuous reaction-diffusion system with Dirichlet boundary conditions and nonnegative solutions, diffusion-driven blow up occurs even when the total species concentration is non-increasing. For systems that instead have homogeneous Neumann boundary conditions, it is currently unknown whether this deviation between the continuous and discretized system can occur. Therefore, it is worth examining the discretized system independently of the continuous system. Since no criteria exist for the boundedness of the discretized system, the focus of this paper is to determine sufficient conditions to guarantee the system with diffusion remains bounded for all time. We consider reaction-diffusion systems on a 1D domain with homogeneous Neumann boundary conditions and non-negative initial data and solutions. We define a Lyapunov-like function and show that its existence guarantees that the discretized reaction-diffusion system is bounded. These results are considered in the context of three example systems for which Lyapunov-like functions can and cannot be found.


title: Single dot photoluminescence excitation spectroscopy in the telecommunication spectral range

abstract: Single dot photoluminescence excitation spectroscopy provides an insight into energy structure of individual quantum dots, energy transfer processes within and between the dots and their surroundings. The access to single dot energy structure is vital for further development of telecom-based quantum emitters, like single photon sources or entangled pair of photons. However, application of single dot photoluminescence excitation spectroscopy is limited mainly to dots emitting below 1 $\mu$m, while nanostructures optically active in the telecommunication windows of 1.3 and 1.55 $\mu$m are of particular interest, as they correspond to the desirable wavelengths in nanophotonic applications. This report presents an approach to photoluminescence excitation spectroscopy covering this application-relevant spectral range on single dot level. Experimental details are discussed, including issues related to the tunable excitation source and its spectral filtering, and illustrated with examples of photoluminescence excitation spectroscopy results from single quantum dots emitting in both the 1.3 and 1.55 $\mu$m spectral ranges.


title: Polynomial chaos expansions for dependent random variables

abstract: Polynomial chaos expansions (PCE) are well-suited to quantifying uncertainty in models parameterized by independent random variables. The assumption of independence leads to simple strategies for evaluating PCE coefficients. In contrast, the application of PCE to models of dependent variables is much more challenging. Three approaches can be used. The first approach uses mapping methods where measure transformations, such as the Nataf and Rosenblatt transformation, can be used to map dependent random variables to independent ones; however we show that this can significantly degrade performance since the Jacobian of the map must be approximated. A second strategy is the class of dominating support methods which build PCE using independent random variables whose distributional support dominates the support of the true dependent joint density; we provide evidence that this approach appears to produce approximations with suboptimal accuracy. A third approach, the novel method proposed here, uses Gram-Schmidt orthogonalization (GSO) to numerically compute orthonormal polynomials for the dependent random variables. This approach has been used successfully when solving differential equations using the intrusive stochastic Galerkin method, and in this paper we use GSO to build PCE using a non-intrusive stochastic collocation method. The stochastic collocation method treats the model as a black box and builds approximations of model output from a set of samples. Building PCE from samples can introduce ill-conditioning which does not plague stochastic Galerkin methods. To mitigate this ill-conditioning we generate weighted Leja sequences, which are nested sample sets, to build accurate polynomial interpolants. We show that our proposed approach produces PCE which are orders of magnitude more accurate than PCE constructed using mapping or dominating support methods.


title: Modern Asset Theory: A Framework for Successful Active Management

abstract: Active management is a term that has many meanings and we have found the defining characteristics needed for success as an "active manager" elusive within the literature. In this paper we offer a set of criteria that defines an active manager and his success. In order to facilitate this, we introduce several definitions, which lead to a logically coherent evaluation framework. We expand on the definitions of these six key concepts: Introduce a specific concept of Ruin, Assets, Risk, Discount rate, Margin of safety, and Optimization. Through these definitions a strong defense of active management emerges. Furthermore, to the extent one chooses to limit the definitions we offer, our framework reduces to that of Modern Portfolio Theory. Overall, we have aimed to construct a robust expansion of a framework for active management, one that has been found wanting in the current literature.


title: Extreme Starlight Polarization in a Region with Highly Polarized Dust Emission

abstract: Galactic dust emission is polarized at unexpectedly high levels, as revealed by Planck. The origin of the observed $\simeq 20\%$ polarization fractions can be identified by characterizing the properties of optical starlight polarization in a region with maximally polarized dust emission. We measure the R-band linear polarization of 22 stars in a region with a submillimeter polarization fraction of $\simeq 20%$. A subset of 6 stars is also measured in the B, V and I bands to investigate the wavelength dependence of polarization. We find that starlight is polarized at correspondingly high levels. Through multiband polarimetry we find that the high polarization fractions are unlikely to arise from unusual dust properties, such as enhanced grain alignment. Instead, a favorable magnetic field geometry is the most likely explanation, and is supported by observational probes of the magnetic field morphology. The observed starlight polarization exceeds the classical upper limit of $\left[p_V/E\left(B-V\right)\right]_{\rm max} = 9$%mag$^{-1}$ and is at least as high as 13%mag$^{-1}$ that was inferred from a joint analysis of Planck data, starlight polarization and reddening measurements. Thus, we confirm that the intrinsic polarizing ability of dust grains at optical wavelengths has long been underestimated.


title: Free-electron effects on optical absorption of hybrid perovskite CH$_3$NH$_3$PbI$_3$ from first principles

abstract: Hybrid organic-inorganic perovskites, such as methyl-ammonium lead tri-iodide (MAPbI$_3$), are interesting candidates for efficient absorber materials in next-generation solar cells, partly due to an unusual combination of low exciton binding energy and strong optical absorption. Excitonic effects in this material have been subject to debate both for experiment and theory, indicating a need for better understanding of the screening mechanisms that act upon the electron-hole interaction. Here we use cutting-edge first-principles theoretical spectroscopy, based on density-functional and many-body perturbation theory, to study atomic geometries, electronic structure, and optical properties of three MAPbI$_3$ polymorphs and find good agreement with earlier results and experiment. We then study the influence of free electrons on the electron-hole interaction and show that this explains consistently smaller exciton binding energies, compared to those in the material without free electrons. Interestingly, we also find that the absorption line shape strongly resembles that of the spectrum without free electrons up to high free electron concentrations. We explain this unexpected behavior by formation of Mahan excitons that dominate the absorption edge, making it robust against free-electron induced changes observed in other semiconductors.


title: Supermassive Black Hole Feedback

abstract: Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes.


title: Can reverberation-measured quasars be used for cosmology?

abstract: Quasars have been proposed as a new class of standard candles analogous to Supernovae, as their large redshift range and high luminosities make them excellent candidates. Reverberation mapping (RM) method enables to estimate the distance to the source from the time delay measurement of the emission lines with respect to the continuum, since the time delay depends on the absolute luminosity of the source. This radius-luminosity (RL) relation exhibits a low scatter and it was thus proposed for cosmological purposes. However, in recent years the increase in the studied sample, and in particular the inclusion of highly accreting QSO has increased the dispersion in the RL relation, with many objects showing time delays shorter than the expected. Using H$\beta$ RM measurements for 117 sources with $0.002<z<0.9$ and $41.5<\mathrm{log}~L_{5100}<45.9$, we find a correction for the time delay based on the dimensionless accretion rate considering the virial factor anti-correlated with the FWHM of H$\beta$. This correction decreases the scattering in the accretion parameters compared with the typical value used, which is directly reflected by suppressing the RL relation dispersion. We also confirm the anti-correlation between the excess of variability and the accretion parameters. With this correction we are able to build the Hubble diagram and estimate the cosmological constants $\Omega_m$ and $\Omega_\Lambda$, which are consistent with the standard cosmological model at 2$\sigma$ confidence level. Therefore, RM results can be used to constrain cosmological models in the future.


title: Symbolic Regression Methods for Reinforcement Learning

abstract: Reinforcement learning algorithms can solve dynamic decision-making and optimal control problems. With continuous-valued state and input variables, reinforcement learning algorithms must rely on function approximators to represent the value function and policy mappings. Commonly used numerical approximators, such as neural networks or basis function expansions, have two main drawbacks: they are black-box models offering little insight into the mappings learned, and they require extensive trial and error tuning of their hyper-parameters. In this paper, we propose a new approach to constructing smooth value functions in the form of analytic expressions by using symbolic regression. We introduce three off-line methods for finding value functions based on a state-transition model: symbolic value iteration, symbolic policy iteration, and a direct solution of the Bellman equation. The methods are illustrated on four nonlinear control problems: velocity control under friction, one-link and two-link pendulum swing-up, and magnetic manipulation. The results show that the value functions yield well-performing policies and are compact, mathematically tractable, and easy to plug into other algorithms. This makes them potentially suitable for further analysis of the closed-loop system. A comparison with an alternative approach using neural networks shows that our method outperforms the neural network-based one.


title: Distributed estimation and control of node centrality in undirected asymmetric networks

abstract: Measures of node centrality that describe the importance of a node within a network are crucial for understanding the behavior of social networks and graphs. In this paper, we address the problems of distributed estimation and control of node centrality in undirected graphs with asymmetric weight values. In particular, we focus our attention on $\alpha$-centrality, which can be seen as a generalization of eigenvector centrality. In this setting, we first consider a distributed protocol where agents compute their $\alpha$-centrality, focusing on the convergence properties of the method; then, we combine the estimation method with a consensus algorithm to achieve a consensus value weighted by the influence of each node in the network. Finally, we formulate an $\alpha$-centrality control problem which is naturally decoupled and, thus, suitable for a distributed setting and we apply this formulation to protect the most valuable nodes in a network against a targeted attack, by making every node in the network equally important in terms of {\alpha}-centrality. Simulations results are provided to corroborate the theoretical findings.


title: Evanescent Gravitational Waves

abstract: We describe the properties of evanescent gravitational waves (EGWs)---wave solutions of Einstein equations which decay exponentially in some direction while propagating in another. Evanescent waves are well-known in acoustics and optics and have recently received much attention due to their extraordinary properties such as their transverse spin and spin-momentum locking. We show that EGWs possess similarly remarkable properties, carrying transverse spin angular momenta and driving freely falling test masses along in elliptical trajectories. Hence, test masses on a plane transverse to the direction of propagation exhibit correlated vector and scalar-like deformation---correlations which can be used to distinguish it from modified gravity. We demonstrate that EGWs are present and dominant in the vicinity of sub-wavelength sources such as orbiting binaries.


title: Rapid and unambiguous determination of the open time for waterborne paint films with Laser Speckle Imaging

abstract: In drying coatings, the time span between the deposition of the wet film and the formation of an uniform solid harbors a procession of complex phenomena each of which carves the final properties of the dry paint film. Attaining a better control of the drying process requires the availability of a tool that allows monitoring the attending dynamics. The experimental access to the dynamics underlying the drying process is challenging as paint films are often opaque and the phenomena accompanying the drying process usually populate a multitude of time and length scales. To circumvent these obstacles we deploy Laser Speckle Imaging (LSI) as a technique that allow us to quantitatively probe the internal dynamics of drying waterborne paints and unambiguously determine a crucial handling parameter such as the open time. We develop a set of scaling relations which accurately predict the behavior of the experimentally determined open time as a function of a parameter set that governs the drying process. Additionally, we harness the big temporal dynamic range of LSI to monitor phenomena that inhabit the later stages of the drying process. Remarkably, we not only capture the film formation time, the time point in which phase inversion takes place, but also microscopically detect the deformation of polymeric particles that precedes it.


title: Astro2020: Hot Drivers of Stellar Feedback from 10 to 10,000 pc

abstract: Stellar feedback -- stars regulating further star formation through the injection of energy and momentum into the interstellar medium -- operates through a complex set of processes that originate in star clusters but shape entire galaxies. A mature theory of stellar feedback is essential to a complete theory of star and galaxy formation, but the energy and momentum injected by hot gas into its surroundings remains unclear. With a next-generation X-ray observatory, breakthrough progress can be made through precision measurements of the temperature, density, velocity, and abundances of hot gas on scales of star clusters to galactic superwinds.


title: Astro2020: How does dust escape from galaxies?

abstract: Whenever gas is blown out of a galaxy, chances are that it contains some cosmic dust. This dust is an important part of the metals budget for the circumgalactic and intergalactic media, and traces the outflow and stripping history of the galaxy. The dust is also interesting in its own right, as dust plays an essential role in many astrophysical processes. We have only begun to learn about circumgalactic dust, and in particular we do not know how (and when) it escapes its host galaxy. Here we describe the prospects for measuring the dust mass and properties around many individual galaxies, which will form the basis for understanding how the dust got there.


title: Does History Repeat Itself? Periodic Time Cosmology

abstract: It has been suggested that the cosmic history might repeat in cycles, with an infinite series of similar aeons in the past and the future. Here, we instead propose that the cosmic history repeats itself exactly, constructing a universe on a periodic temporal history, which we call Periodic Time Cosmology. In particular, the primordial power spectrum, convolved with the transfer function throughout the cosmic history, would form the next aeon's primordial power spectrum. By matching the big bang to the infinite future using a conformal rescaling (a la Penrose), we uniquely determine the primordial power spectrum, in terms of the transfer function up to two free parameters. While nearly scale invariant with a red tilt on large scales, using Planck and Baryonic Acoustic Oscillation observations, we find the minimal model is disfavoured compared to a power-law power spectrum at $5.1\sigma$. However, extensions of $\Lambda$CDM cosmic history change the large scale transfer function and can provide better relative fits to the data. For example, the best fit seven parameter model for our Periodic Time Cosmology, with $w=-1.024$ for dark energy equation of state, is only disfavoured relative to a power-law power spectrum (with the same number of parameters) at $1.8\sigma$ level. Therefore, consistency between cosmic history and initial conditions provides a viable description of cosmological observations in the context of Periodic Time Cosmology.


title: On the generalized distributive set of a finite nearfield

abstract: For any nearfield $(R,+, \circ)$, denote by $D(R)$ the set of all distributive elements of $R$. Let $R$ be a finite Dickson nearfield that arises from Dickson pair $(q,n)$. For a given pair $(\alpha, \beta) \in R^2$ we study the generalized distributive set $ D(\alpha, \beta)$ where $"\circ"$ is the multiplication of the Dickson nearfield. We find that $ D(\alpha, \beta)$ is not in general a subfield of the finite field $\mathbb{F}_{q^n}$. In contrast to the situation for $D(R)$, we also find that $D(\alpha, \beta)$ is not in general a subnearfield of $R$. We obtain sufficient conditions on $\alpha, \beta$ for $ D(\alpha, \beta)$ to be a subfield of $\mathbb{F}_{q^n}$ and derive an algorithm that tests if $D(\alpha, \beta)$ is a subfield of $\mathbb{F}_{q^n}$ or not. We also study the notions of $R$-dimension, $R$-basis, seed sets and seed number of $R$-subgroups of the Beidleman near-vector spaces $R^m$ where $m$ is a positive integer. Finally we determine the maximal $R$-dimension of $gen(v_1,v_2)$ for $v_1,v_2 \in R^m$, where $gen(v_1,v_2)$ is the smallest $R$-subgroup containing the vectors $v_1$ and $v_2$.


title: Algebras of Continuous Fourier Multipliers on Variable Lebesgue Spaces

abstract: We show that several definitions of algebras of continuous Fourier multipliers on variable Lebesgue spaces over the real line are equivalent under some natural assumptions on variable exponents. Some of our results are new even in the case of standard Lebesgue spaces and give answers on two questions about algebras of continuous Fourier multipliers on Lebesgue spaces over the real line posed by H. Mascarenhas, P. Santos and M. Seidel.


title: Unruh effect for fermions from the Zubarev density operator

abstract: Using the Zubarev quantum-statistical density operator, we calculated the corrections to the energy-momentum tensor of a massless fermion gas associated with acceleration. It is shown that when fourth-order corrections are taken into account, the energy-momentum tensor in the laboratory frame is equal to zero at a proper temperature measured by a comoving observer equal to Unruh temperature. Consequently, the Minkowski vacuum is visible to the accelerated observer as a medium filled with a heat bath of particles with the Unruh temperature, which is the essence of the Unruh effect.


title: CUR Decompositions, Approximations, and Perturbations

abstract: This article discusses a useful tool in dimensionality reduction and low-rank matrix approximation called the CUR decomposition. Various viewpoints of this method in the literature are synergized and are compared and contrasted; included in this is a new characterization of exact CUR decompositions. A novel perturbation analysis is performed on CUR approximations of noisy versions of low-rank matrices, which compares them with the putative CUR decomposition of the underlying low-rank part. Additionally, we give new column and row sampling results which allow one to conclude that a CUR decomposition of a low-rank matrix is attained with high probability. We then illustrate the stability of these sampling methods under the perturbations studied before, and provide numerical illustrations of the methods and bounds discussed.


title: Spin-mechanics with levitating ferromagnetic particles

abstract: We propose and demonstrate first steps towards schemes where the librational mode of levitating ferromagnets is strongly coupled to the electronic spin of Nitrogen-Vacancy (NV) centers in diamond. Experimentally, we levitate ferromagnets in a Paul trap and employ magnetic fields to attain oscillation frequencies in the hundreds of kHz range with Q factors close to $10^4$. These librational frequencies largely exceed the decoherence rate of NV centers in typical CVD grown diamonds offering prospects for sideband resolved operation. We also prepare and levitate composite diamond-ferromagnet particles and demonstrate both coherent spin control of the NV centers and read-out of the particle libration using the NV spin. Our results will find applications in ultra-sensitive gyroscopy and bring levitating objects a step closer to spin-mechanical experiments at the quantum level.


title: Modified exchange interaction between magnetic impurities in spin-orbit coupled quantum wires

abstract: Indirect exchange interaction between magnetic impurities in one dimensional systems is a matter of long discussions since Kittel has established that in the asymptotic limit it decays as the inverse of distance x between the impurities. In this work we investigate this problem in a quantum wire with Rashba spin-orbit coupling (SOC). By employing a second order perturbation theory we find that one additional oscillatory term appears in each of the RKKY, the Dzaloshinkii-Moryia and the Ising couplings. Remarkably, these extra terms resulting from the spin precession of the conduction electrons induced by the SOC do not decay as in the usual RKKY interaction. We show that these extra oscillations arise from the finite momenta band splitting induced by the spin-orbit coupling that modifies the spin-flip scatterings occurring at the Fermi energy. Our findings open up a new perspective in the long-distance magnetic interactions in solid state systems.


title: Quality of Experience from Cache Hierarchies: Keep your low-bitrate close, and high-bitrate closer

abstract: Recent studies into streaming media delivery suggest that performance gains from cache hierarchies such as Information-Centric Networks (ICNs) may be negated by Dynamic Adaptive Streaming (DAS), the de facto method for retrieving multimedia content. The bitrate adaptation mechanisms that drive video streaming clash with caching hierarchies in ways that affect users' Quality of Experience (QoE). Cache performance also diminishes as consumers dynamically select content encoded at different bitrates. In this paper we use the evidence to draw a novel insight: in a cache hierarchy for adaptive streaming content, bitrates should be prioritized over or alongside popularity and hit rates. We build on this insight to propose RippleCache as a family of cache placement schemes that safeguard high-bitrate content at the edge and push low-bitrate content into the network core. Doing so reduces contention of cache resources, as well as congestion in the network. To validate RippleCache claims we construct two separate implementations. We design RippleClassic as a benchmark solution that optimizes content placement by maximizing a measure for cache hierarchies shown to have high correlation with QoE. In addition, our lighter-weight RippleFinder is then re-designed with distributed execution for application in large-scale systems. RippleCache performance gains are reinforced by evaluations in NS-3 against state-of-the-art baseline approaches, using standard measures of QoE as defined by the DASH Industry Forum. Measurements show that RippleClassic and RippleFinder deliver content that suffers less oscillation and rebuffering, as well as the highest levels of video quality, indicating overall improvements to QoE.


title: Deterministic Field-free Switching of a Perpendicularly Magnetized Ferromagnetic Layer via the Joint Effects of Dzyaloshinskii-Moriya Interaction and Field-like Spin-orbit Torque: An Appraisal

abstract: Field-free switching of perpendicularly magnetized ferromagnetic layer by spin orbit torque (SOT) from the spin Hall effect (SHE) is of great interest in the applications of magnetic memory devices. In this paper, we investigate the deterministic SOT switching through the joint effects of Dzyaloshinskii-Moriya Interaction (DMI) and field-like torque (FLT) by micromagnetic simulations. We confirmed that within a certain range of DMI values and charge current densities, it is possible to deterministically switch the magnetization without the assistance external magnetic field. We show that the FLT could play an adverse role in blocking and slowing down the magnetization switching under certain cases of DMI and charge current-driven field-free switching. However, in other cases, FLT can assist DMI on the deterministic field-free SOT switching. In addition, it is found that FLT can effectively expand the current density window for a deterministic field-free SOT switching in the presence of DMI.


title: The development of nations conditions the disease space

abstract: Using the economic complexity methodology on data for disease prevalence in 195 countries during the period of 1990-2016, we propose two new metrics for quantifying the relatedness between diseases, or the `disease space' of countries. With these metrics, we analyze the geography of diseases and empirically investigate the effect of economic development on the health complexity of countries. We show that a higher income per capita increases the complexity of countries' diseases. Furthermore, we build a disease-level index that links a disease to the average level of GDP per capita of the countries that have prevalent cases of the disease. With this index, we highlight the link between economic development and the complexity of diseases and illustrate, at the disease-level, how increases in income per capita are associated with more complex diseases


title: Inverse-chirp signals and spontaneous scalarisation with self-interacting potentials in stellar collapse

abstract: We study how the gravitational wave signal from stellar collapse in scalar-tensor gravity varies under the influence of scalar self-interaction. To this end, we extract the gravitational radiation from numerical simulations of stellar collapse for a range of potentials with higher-order terms in addition to the quadratic mass term. Our study includes collapse to neutron stars and black holes and we find the strong inverse-chirp signals obtained for the purely quadratic potential to be exceptionally robust under changes in the potential at higher orders; quartic and sextic terms in the potential lead to noticeable differences in the wave signal only if their contribution is amplified, implying a relative fine-tuning to within 5 or more orders of magnitude between the mass and self-interaction parameters.


title: Effective Aspects of Bernoulli Randomness

abstract: In this paper, we study Bernoulli random sequences, i.e., sequences that are Martin-L\"of random with respect to a Bernoulli measure $\mu_p$ for some $p\in[0,1]$, where we allow for the possibility that $p$ is noncomputable. We focus in particular on the case in which the underlying Bernoulli parameter $p$ is proper (that is, Martin-L\"of random with respect to some computable measure). We show for every Bernoulli parameter $p$, if there is a sequence that is both proper and Martin-L\"of random with respect to $\mu_p$, then $p$ itself must be proper, and explore further consequences of this result. We also study the Turing degrees of Bernoulli random sequences, showing, for instance, that the Turing degrees containing a Bernoulli random sequence do not coincide with the Turing degrees containing a Martin-L\"of random sequence. Lastly, we consider several possible approaches to characterizing blind Bernoulli randomness, where the corresponding Martin-L\"of tests do not have access to the Bernoulli parameter $p$, and show that these fail to characterize blind Bernoulli randomness.


title: HERBS I: Metallicity and alpha enhancement along the Galactic bulge minor axis

abstract: To better understand the origin and evolution of the Milky Way bulge, we have conducted a survey of bulge red giant branch and clump stars using the HERMES spectrograph on the Anglo-Australian Telescope. We targeted ARGOS survey stars with pre-determined bulge memberships, covering the full metallicity distribution function. The spectra have signal-to-noise ratios comparable to, and were analysed using the same methods as the GALAH survey. In this work, we present the survey design, stellar parameters, distribution of metallicity and alpha-element abundances along the minor bulge axis at latitudes $b$ = $-10^{\circ}, -7.5^{\circ}$ and $-5^{\circ}$. Our analysis of ARGOS stars indicates that the centroids of ARGOS metallicity components should be located $\approx$0.09 dex closer together. The vertical distribution of $\alpha$-element abundances is consistent with the varying contributions of the different metallicity components. Closer to the plane, alpha abundance ratios are lower as the metal-rich population dominates. At higher latitudes, the alpha abundance ratios increase as the number of metal-poor stars increases. However, we find that the trend of alpha-enrichment with respect to metallicity is independent of latitude. Comparison of our results with those of GALAH DR2 revealed that for [Fe/H] $\approx -0.8$, the bulge shares the same abundance trend as the high-$\alpha$ disk population. However, the metal-poor bulge population ([Fe/H] $\lesssim -0.8$) show enhanced alpha abundance ratios compared to the disk/halo. These observations point to fairly rapid chemical evolution in the bulge, and that the metal-poor bulge population does not share the same similarity with the disk as the more metal-rich populations.


title: On moments and strong local H\"older regularity of solutions of stochastic differential equations and of their spatial derivative processes

abstract: Spatial differentiability of solutions of stochastic differential equations (SDEs) is a classical question in stochastic analysis. The case of coefficients with globally Lipschitz continuous derivatives is well understood in the literature. Counterexamples with smooth and bounded coefficients demonstrate that the non-globally Lipschitz case is more subtle. In this article we establish conditions, including a suitable local monotonicity property, which provide existence of continuously differentiable solutions of SDEs, moment estimates and strong local H\"older regularity.


title: Explaining Reinforcement Learning to Mere Mortals: An Empirical Study

abstract: We present a user study to investigate the impact of explanations on non-experts' understanding of reinforcement learning (RL) agents. We investigate both a common RL visualization, saliency maps (the focus of attention), and a more recent explanation type, reward-decomposition bars (predictions of future types of rewards). We designed a 124 participant, four-treatment experiment to compare participants' mental models of an RL agent in a simple Real-Time Strategy (RTS) game. Our results show that the combination of both saliency and reward bars were needed to achieve a statistically significant improvement in mental model score over the control. In addition, our qualitative analysis of the data reveals a number of effects for further study.


title: An Interaction Framework for Studying Co-Creative AI

abstract: Machine learning has been applied to a number of creative, design-oriented tasks. However, it remains unclear how to best empower human users with these machine learning approaches, particularly those users without technical expertise. In this paper we propose a general framework for turn-based interaction between human users and AI agents designed to support human creativity, called {co-creative systems}. The framework can be used to better understand the space of possible designs of co-creative systems and reveal future research directions. We demonstrate how to apply this framework in conjunction with a pair of recent human subject studies, comparing between the four human-AI systems employed in these studies and generating hypotheses towards future studies.


title: Quantifying the impact of network structure on speed and accuracy in collective decision-making

abstract: Found in varied contexts from neurons to ants to fish, binary decision-making is one of the simplest forms of collective computation. In this process, information collected by individuals about an uncertain environment is accumulated to guide behavior at the aggregate scale. We study binary decision-making dynamics in networks responding to inputs with small signal-to-noise ratios, looking for quantitative measures of collectivity that control decision-making performance. We find that decision accuracy is controlled largely by three factors: the leading eigenvalue of the network adjacency matrix, the corresponding eigenvector's participation ratio, and distance from the corresponding symmetry-breaking bifurcation. This allows us to predict how decision-making performance scales in large networks based on their spectral properties. Specifically, we explore the effects of localization caused by the hierarchical assortative structure of a "rich club" topology. This gives insight into the tradeoffs involved in the higher-order structure found in living networks performing collective computations.


title: Barrier Functions in Cascaded Controller: Safe Quadrotor Control

abstract: Safe control for inherently unstable systems such as quadrotors is crucial. Imposing multiple dynamic constraints simultaneously on the states for safety regulation can be a challenging problem. In this paper, we propose a quadratic programming (QP) based approach on a cascaded control architecture for quadrotors to enforce safety. Safety regions are constructed using control barrier functions (CBF) while explicitly considering the nonlinear underactuated dynamics of the quadrotor. The safety regions constructed using CBFs establish a non-conservative forward invariant safe region for quadrotor navigation. Barriers imposed across the cascaded architecture allows independent safety regulation in quadrotor's altitude and lateral domains. Despite barriers appearing in a cascaded fashion, we show preservation of safety for quadrotor motion in SE(3). We demonstrate the feasibility of our method on a quadrotor in simulation with static and dynamic constraints enforced on position and velocity spaces simultaneously.


title: Weak Electric-Field Detection with Sub-1 Hz Resolution at Radio Frequencies Using A Rydberg Atom-Based Mixer

abstract: Rydberg atoms have been used for measuring radio-frequency (RF) electric (E)-fields due to their strong dipole moments over the frequency range of 500 MHz-1 THz. For this, electromagnetically induced transparency (EIT) within the Autler-Townes (AT) regime is used such that the detected E-field is proportional to AT splitting. However, for weak E-fields AT peak separation becomes unresolvable thus limiting the minimum detectable E-field. Here, we demonstrate using the Rydberg atoms as an RF mixer for weak E-field detection well below the AT regime with frequency discrimination better than 1 Hz resolution. Two E-fields incident on a vapor cell filled with cesium atoms are used. One E-field at 19.626000 GHz drives the 34D_(5/2)->5P_(3/2) Rydberg transition and acts as a local oscillator (LO) and a second signal E-field (Sig) of interest is at 19.626090 GHz. In the presence of the LO, the Rydberg atoms naturally down convert the Sig field to a 90 kHz intermediate frequency (IF) signal. This IF signal manifests as an oscillation in the probe laser intensity through the Rydberg vapor and is easily detected with a photodiode and lock-in amplifier. In the configuration used here, E-field strength down to ? 46 mV/m +/-2 mV/m were detected. Furthermore, neighboring fields 0.1 Hz away and equal in strength to Sig could be discriminated without any leakage into the lock-in signal. For signals 1 Hz away and as high as +60 dB above Sig, leakage into the lock-in signal could be kept below -3 dB.


title: SLING: Using Dynamic Analysis to Infer Program Invariants in Separation Logic

abstract: We introduce a new dynamic analysis technique to discover invariants in separation logic for heap-manipulating programs. First, we use a debugger to obtain rich program execution traces at locations of interest on sample inputs. These traces consist of heap and stack information of variables that point to dynamically allocated data structures. Next, we iteratively analyze separate memory regions related to each pointer variable and search for a formula over predefined heap predicates in separation logic to model these regions. Finally, we combine the computed formulae into an invariant that describes the shape of explored memory regions. We present SLING, a tool that implements these ideas to automatically generate invariants in separation logic at arbitrary locations in C programs, e.g., program pre and postconditions and loop invariants. Preliminary results on existing benchmarks show that SLING can efficiently generate correct and useful invariants for programs that manipulate a wide variety of complex data structures.


title: Graph Temporal Logic Inference for Classification and Identification

abstract: Inferring spatial-temporal properties from data is important for many complex systems, such as additive manufacturing systems, swarm robotic systems and biological networks. Such systems can often be modeled as a labeled graph where labels on the nodes and edges represent relevant measurements such as temperatures and distances. We introduce graph temporal logic (GTL) which can express properties such as "whenever a node's label is above 10, for the next 3 time units there are always at least two neighboring nodes with an edge label of at most 2 where the node labels are above 5". This paper is a first attempt to infer spatial (graph) temporal logic formulas from data for classification and identification. For classification, we infer a GTL formula that classifies two sets of graph temporal trajectories with minimal misclassification rate. For identification, we infer a GTL formula that is informative and is satisfied by the graph temporal trajectories in the dataset with high probability. The informativeness of a GTL formula is measured by the information gain with respect to given prior knowledge represented by a prior probability distribution. We implement the proposed approach to classify the graph patterns of tensile specimens built from selective laser sintering (SLS) process with varying strengths, and to identify informative spatial-temporal patterns from experimental data of the SLS cooldown process and simulation data of a swarm of robots.


title: Computing 3 point correlation function randoms counts without the randoms catalogue

abstract: As we move towards future galaxy surveys, the three-point statistics will be increasingly leveraged to enhance the constraining power of the data on cosmological parameters. An essential part of the three-point function estimation is performing triplet counts of synthetic data points in random catalogues. Since triplet counting algorithms scale at best as $\mathcal{O}(N^2\log N)$ with the number of particles and the random catalogues are typically at least 50 times denser than the data; this tends to be by far the most time-consuming part of the measurements. Here we present a simple method of computing the necessary triplet counts involving uniform random distributions through simple one-dimensional integrals. The method speeds up the computation of the three-point function by orders of magnitude, eliminating the need for random catalogues, with the simultaneous pair and triplet counting of the data points alone being sufficient.


title: Monte Carlo study on the detection of classical order by disorder in real antiferromagnetic Ising pyrochlores

abstract: We use Monte Carlo simulations to evaluate the feasibility of detecting thermal order by disorder in real antiferromagnetic Ising pyrochlores, frustrated by a magnetic field applied in the $[110]$ direction. Building on an ideal system with only nearest-neighbour exchange interactions and a perfectly oriented field, we consider the effects of dipolar interactions and field misalignment. Our approach is special in that it relies more in the possibility to switch on and off the entropic drive towards order than in the absence of (or immunity to) a particular perturbation. It can then be applied, in principle, to other uncontrolled interactions expected to be naturally present in real {magnetic} materials. We establish the conditions under which entropic effects can be discerned from an interaction drive towards order, show how to use neutron scattering as a means to unveil this mechanism, and discuss possible materials where to test these ideas.


title: Instance and Output Optimal Parallel Algorithms for Acyclic Joins

abstract: Massively parallel join algorithms have received much attention in recent years, while most prior work has focused on worst-optimal algorithms. However, the worst-case optimality of these join algorithms relies on hard instances having very large output sizes, which rarely appear in practice. A stronger notion of optimality is {\em output-optimal}, which requires an algorithm to be optimal within the class of all instances sharing the same input and output size. An even stronger optimality is {\em instance-optimal}, i.e., the algorithm is optimal on every single instance, but this may not always be achievable. In the traditional RAM model of computation, the classical Yannakakis algorithm is instance-optimal on any acyclic join. But in the massively parallel computation (MPC) model, the situation becomes much more complicated. We first show that for the class of r-hierarchical joins, instance-optimality can still be achieved in the MPC model. Then, we give a new MPC algorithm for an arbitrary acyclic join with load $O ({\IN \over p} + {\sqrt{\IN \cdot \OUT} \over p})$, where $\IN,\OUT$ are the input and output sizes of the join, and $p$ is the number of servers in the MPC model. This improves the MPC version of the Yannakakis algorithm by an $O (\sqrt{\OUT \over \IN} )$ factor. Furthermore, we show that this is output-optimal when $\OUT = O(p \cdot \IN)$, for every acyclic but non-r-hierarchical join. Finally, we give the first output-sensitive lower bound for the triangle join in the MPC model, showing that it is inherently more difficult than acyclic joins.


title: Observing exoplanets in the near-infrared from a high altitude balloon platform

abstract: Although there exists a large sample of known exoplanets, little spectroscopic data exists that can be used to study their global atmospheric properties. This deficiency can be addressed by performing phase-resolved spectroscopy -- continuous spectroscopic observations of a planet's entire orbit about its host star -- of transiting exoplanets. Planets with characteristics suitable for atmospheric characterization have orbits of several days, thus phase curve observations are highly resource intensive, especially for shared use facilities. In this work, we show that an infrared spectrograph operating from a high altitude balloon platform can perform phase-resolved spectroscopy of hot Jupiter-type exoplanets with performance comparable to a space-based telescope. Using the EXoplanet Climate Infrared TElescope (EXCITE) experiment as an example, we quantify the impact of the most important systematic effects that we expect to encounter from a balloon platform. We show an instrument like EXCITE will have the stability and sensitivity to significantly advance our understanding of exoplanet atmospheres. Such an instrument will both complement and serve as a critical bridge between current and future space-based near infrared spectroscopic instruments.


title: Investigation of High Energy Behaviour of HERA Data

abstract: We analyse the high precision HERA $F_2$ data in the low-$x$ regions, $x<0,01$ and the very-low-$x$, $x<0.001$, regions using $\lambda$-Fits. $\lambda$ is a measure of the rate of rise of $F_2$ defined by $F_2 \propto (1/x)^{\lambda}$. We show that $\lambda$ determined in these two regions, at various $Q^2$ values, is systematically smaller in the very-low-$x$ region as compared to the low-$x$ region. We discuss some possible physical interpretations of this effect.


title: Desorption Kinetics and Binding Energies of Small Hydrocarbons

abstract: Small hydrocarbons are an important organic reservoir in protostellar and protoplanetary environments. Constraints on desorption temperatures and binding energies of such hydrocarbons are needed for accurate predictions of where these molecules exist in the ice vs. gas-phase during the different stages of star and planet formation. Through a series of temperature programmed desorption (TPD) experiments, we constrain the binding energies of 2 and 3-carbon hydrocarbons (C$_{2}$H$_{2}$ - acetylene, C$_{2}$H$_{4}$ - ethylene, C$_{2}$H$_{6}$ - ethane, C$_{3}$H$_{4}$ - propyne, C$_{3}$H$_{6}$ - propene, and C$_{3}$H$_{8}$ - propane) to 2200-4200 K in the case of pure amorphous ices, to 2400-4400 K on compact amorphous H$_{2}$O, and to 2800-4700 K on porous amorphous H$_{2}$O. The 3-carbon hydrocarbon binding energies are always larger than the 2-carbon hydrocarbon binding energies. Within the 2- and 3-carbon hydrocarbon families, the alkynes (i.e., least-saturated) hydrocarbons exhibit the largest binding energies, while the alkane and alkene binding energies are comparable. Binding energies are $\sim$5-20% higher on water ice substrates compared to pure ices, which is a small increase compared to what has been measured for other volatile molecules such as CO and N$_{2}$. Thus in the case of hydrocarbons, H$_{2}$O has a less pronounced effect on sublimation front locations (i.e., snowlines) in protoplanetary disks.


title: Automorphisms of Kronrod-Reeb graphs of Morse functions on 2-sphere

abstract: Let $M$ be a compact two-dimensional manifold and, $f \in C^{\infty}(M,\mathbb{R})$ be a Morse function, and $\Gamma_f$ be its Kronrod-Reeb graph. Denote by $\mathcal{O}_{f}=\{f \circ h \mid h \in \mathcal{D}\}$ the orbit of $f$ with respect to the natural right action of the group of diffeomorphisms $\mathcal{D}$ on $C^{\infty}(M,\mathbb{R})$, and by $\mathcal{S}(f)=\{h\in\mathcal{D} \mid f \circ h = f\}$ the corresponding stabilizer of this function. It is easy to show that each $h\in\mathcal{S}(f)$ induces a homeomorphism of $\Gamma_f$. Let also $\mathcal{D}_{\mathrm{id}}(M)$ be the identity path component of $\mathcal{D}(M)$, $\mathcal{S}'(f)= \mathcal{S}(f) \cap \mathcal{D}_{\mathrm{id}}(M)$ be group of diffeomorphisms of $M$ preserving $f$ and isotopic to identity map, and $G_f$ be the group of homeomorphisms of the graph $\Gamma_f$ induced by diffeomorphisms belonging to $\mathcal{S}'(f)$. This group is one of the key ingredients for calculating the homotopy type of the orbit $\mathcal{O}_{f}$. Recently the authors described the structure of groups $G_f$ for Morse functions on all orientable surfaces distinct from $2$-torus $T^2$ and $2$-sphere $S^2$. The present paper is devoted to the case $M=S^{2}$. In this situation $\Gamma_f$ is always a tree, and therefore all elements of the group $G_f$ have a common fixed subtree $\mathrm{Fix}(G_f)$, which may even consist of a unique vertex. Our main result calculates the groups $G_f$ for all Morse functions $f:S^{2}\to\mathbb{R}$ whose fixed subtree $\mathrm{Fix}(G_f)$ consists of more than one point.


title: Pre-trained Language Model Representations for Language Generation

abstract: Pre-trained language model representations have been successful in a wide range of language understanding tasks. In this paper, we examine different strategies to integrate pre-trained representations into sequence to sequence models and apply it to neural machine translation and abstractive summarization. We find that pre-trained representations are most effective when added to the encoder network which slows inference by only 14%. Our experiments in machine translation show gains of up to 5.3 BLEU in a simulated resource-poor setup. While returns diminish with more labeled data, we still observe improvements when millions of sentence-pairs are available. Finally, on abstractive summarization we achieve a new state of the art on the full text version of CNN/DailyMail.


title: On extensions of partial isometries

abstract: In this paper we define a notion of S-extension for a metric space and study minimality and coherence of S-extensions. We show that every S-extension can be identified with an algebraic object. We use this algebraic representation to give a complete characterization of all finite minimal S-extensions of a given finite metric space and a complete characterization of all minimal coherent S-extensions. We also define a notion of ultraextensive metric spaces and show that every countable metric space can be extended to a countable ultraextensive metric space. %As an application, we show that every countable subset of the Urysohn metric space can be extended to a countable dense ultraextensive subset of the Urysohn space. We also show that the isometry group of an infinite ultraextensive metric space has a dense locally finite subgroup, generalizing several previously known results. We also study compact ultrametric spaces and show that every compact ultrametric space can be extended to a compact ultraextensive ultrametric space.


title: Tunable magnetism of a hexagonal Anderson droplet on the triangular lattice

abstract: Motivated by recent progress on quantum engineered Kondo lattices, we numerically investigated the local magnetic properties of a hexagonal Anderson droplet consisting of multiple rings of magnetic atoms periodically arrayed on a triangular lattice. We demonstrated the tunability of the magnetic properties via their evolution with the droplet geometry for two types of systems with distinct local orbital occupancy profile. We found that the local susceptibility of the droplet center of some types of droplets can be remarkably enhanced in contrast to the conventionally rapid decrease due to spin correlations of surrounding droplet rings. The tunability of the magnetic properties is attributed to the charge redistribution with varying the droplet geometry enforced by the confined lattice with open boundary. Our simulations complements the exploration of the novel artificial tunability of engineered lattice systems.


title: Large homogeneous subgraphs in bipartite graphs with forbidden induced subgraphs

abstract: For a bipartite graph G, let h(G) be the largest t such that either G or the bipartite complement of G contain K_{t,t}. For a class F of graphs, let h(F)= min {h(G): G\in F}. We say that a bipartite graph H is strongly acyclic if neither H nor its bipartite complement contain a cycle. By Forb(n, H) we denote a set of bipartite graphs with parts of sizes n each, that do not contain H as an induced bipartite subgraph respecting the sides. One can easily show that h(Forb(n,H))= O(n^{1-s}) for a positive s if H is not strongly acyclic. Here, we prove that h(Forb(n, H)) is linear in n for all strongly acyclic graphs except for four graphs.


title: Optical Bloch modeling of few-cycle laser-induced electron dynamics in dielectrics

abstract: We develop a new model of laser-matter interaction based on Optical Bloch Equations, which includes photo-ionization, impact ionization, and various relaxation processes typical of dielectric materials. This approach is able to describe the temporal evolution of the electron dynamics in the conduction band driven by few-cycle laser pulses of any wavelength. Moreover, the nonlinear polarization response of both centrosymmetric and non-centrosymmetric materials can be described while ensuring the proper selection rules for the harmonics emission.