title: Monotonicity properties arising in a simple model of Wolbachia invasion for wild mosquito populations

abstract: In this paper, we propose a simplified bidimensional Wolbachia infestation model in a population of Aedes aegypti mosquitoes, preserving the main features associated with the biology of this species that can be found in higher-dimensional models. Namely, our model represents the maternal transmission of the Wolbachia symbiont, expresses the reproductive phenotype of cytoplasmic incompatibility, accounts for different fecundities and mortalities of infected and wild insects, and exhibits the bistable nature leading to the so-called principle of competitive exclusion. Since Wolbachia-based biocontrol is now accepted as an ecologically friendly and potentially cost-effective method for prevention and control of dengue and other arboviral infections, it is essential to have reduced models with the main biological characteristics of Aedes aegypti in the presence of Wolbachia-carriers because such models help to simplify the mathematical analysis for determining appropriate biocontrol strategies. Using tools borrowed from monotone dynamical system theory, in the proposed model, we prove the existence of an invariant threshold manifold that allows us to provide practical recommendations for performing single and periodic releases of Wolbachia-carrying mosquitoes, seeking the eventual elimination of wild insects that are capable of transmitting infections to humans. We illustrate these findings with numerical simulations using parameter values corresponding to the wMelPop strain of Wolbachia that is considered the best virus blocker but induces fitness loss in its carriers.


title: Quickest Detection of COVID-19 Pandemic Onset

abstract: This paper develops an easily-implementable version of Page's CUSUM quickest-detection test, designed to work in certain composite hypothesis scenarios with time-varying data statistics. The decision statistic can be cast in a recursive form and is particularly suited for on-line analysis. By back-testing our approach on publicly-available COVID-19 data we find reliable early warning of infection flare-ups, in fact sufficiently early that the tool may be of use to decision-makers on the timing of restrictive measures that may in the future need to be taken.


title: On Rastall gravity formulation as a $f(R,\mathcal{L}_m)$ and a $f(R,T)$ theory

abstract: Rastall introduced a stress-energy tensor whose divergence is proportional to the gradient of the Ricci scalar. This proposal leads to a change in the form of the field equations of General Relativity, but it preserves the number of degrees of freedom. Rastall's field equations can be either interpreted as GR with a redefined SET, or it can imply different physical consequences inside the matter sector. We investigate limits under which the Rastall field equations can be directly derived from an action, in particular from two $f(R)$-gravity extensions: $f(R,\mathcal L_m)$ and $f(R,T)$. We show that there are similarities between these theories, but the Rastall SET cannot be fully recovered from them, apart from certain particular cases here discussed. It is remarkable that a simple, covariant and invertible redefinition of the SET, as the one proposed by Rastall, is hard to be directly implemented in the action.


title: Effects of the Background Turbulence on the Relaxation of Ion Temperature Anisotropy in Space Plasmas

abstract: Turbulence in space plasmas usually exhibits two regimes separated by a spectral break that divides the so called inertial and kinetic ranges. Large scale magnetic fluctuations are dominated by non-linear MHD wave-wave interactions following a -5/3 or -2 slope power-law spectrum. After the break, at scales in which kinetic effects take place, the magnetic spectrum follows a steeper power-law $k^{-\alpha}$ shape given by a spectral index $\alpha > 5/3$. Despite its ubiquitousness, the possible effects of a turbulent background spectrum in the quasilinear relaxation of solar wind temperatures are usually not considered. In this work, a quasilinear kinetic theory is used to study the evolution of the proton temperatures in an initially turbulent collisionless plasma composed by cold electrons and bi-Maxwellian protons, in which electromagnetic waves propagate along a background magnetic field. Four wave spectrum shapes are compared with different levels of wave intensity. We show that a sufficient turbulent magnetic power can drive stable protons to transverse heating, resulting in an increase in the temperature anisotropy and the reduction of the parallel proton beta. Thus, stable proton velocity distribution can evolve in such a way as to develop kinetic instabilities. This may explain why the constituents of the solar wind can be observed far from thermodynamic equilibrium and near the instability thresholds.


title: Synthetic Image Rendering Solves Annotation Problem in Deep Learning Nanoparticle Segmentation

abstract: Nanoparticles occur in various environments as a consequence of man-made processes, which raises concerns about their impact on the environment and human health. To allow for proper risk assessment, a precise and statistically relevant analysis of particle characteristics (such as e.g. size, shape and composition) is required that would greatly benefit from automated image analysis procedures. While deep learning shows impressive results in object detection tasks, its applicability is limited by the amount of representative, experimentally collected and manually annotated training data. Here, we present an elegant, flexible and versatile method to bypass this costly and tedious data acquisition process. We show that using a rendering software allows to generate realistic, synthetic training data to train a state-of-the art deep neural network. Using this approach, we derive a segmentation accuracy that is comparable to man-made annotations for toxicologically relevant metal-oxide nanoparticle ensembles which we chose as examples. Our study paves the way towards the use of deep learning for automated, high-throughput particle detection in a variety of imaging techniques such as microscopies and spectroscopies, for a wide variety of studies and applications, including the detection of plastic micro- and nanoparticles.

categories: cs.LG cond-mat.mtrl-sci cs.CV eess.IV physics.app-ph


title: Mechanical, optoelectronic and thermoelectric properties of half-Heusler p-type semiconductor BaAgP: A DFT investigation

abstract: We have explored the mechanical, electronic, optical and thermoelectric properties of p-type half-Heusler compound BaAgP for the first time using density functional theory based calculations. The mechanical and dynamical stability of this compound is confirmed by studying the Born stability criteria and phonon dispersion curve, respectively. It is soft, ductile and elastically anisotropic. The atomic bonding along a-axis is stronger than that along c-axis. The calculated electronic structure reveals that the studied compound is an indirect band gap semiconductor. The analysis of charge density distribution map and Mulliken population reveals that the bonding in BaAgP is a mixture of covalent and ionic. The optical features confirm that BaAgP is optically anisotropic. The high absorption coefficient and low reflectivity in the visible to ultraviolet region make this compound a possible candidate for solar cell and optoelectronic device applications. The thermoelectric properties have been evaluated by solving the Boltzmann semi-classical transport equations. The calculated power factor at 1000K along a-axis is 35.2 micro-W/cmK2 (with tau=10-14 s) which is ~3.5 times larger than that of SnSe, a promising layered thermoelectric materials. The thermoelectric figure of merit, ZT of BaAgP is 0.44 which is small due to high thermal conductivity. So the reduction of thermal conductivity is essential to enhance thermoelectric performance of BaAgP in device applications.


title: Digital-Analog Quantum Simulations Using The Cross-Resonance Effect

abstract: Digital-analog quantum computation aims to reduce the currently infeasible resource requirements needed for near-term quantum information processing by replacing sequences of one- and two-qubit gates with a unitary transformation generated by the systems' underlying Hamiltonian. Inspired by this paradigm, we consider superconducting architectures and extend the cross-resonance effect, up to first order in perturbation theory, from a two-qubit interaction to an analog Hamiltonian acting on 1D chains and 2D square lattices which, in an appropriate reference frame, results in a purely two-local Hamiltonian. By augmenting the analog Hamiltonian dynamics with single-qubit gates we show how one may generate a larger variety of distinct analog Hamiltonians. We then synthesize unitary sequences, in which we toggle between the various analog Hamiltonians as needed, simulating the dynamics of Ising, $XY$, and Heisenberg spin models. Our dynamics simulations are Trotter error-free for the Ising and $XY$ models in 1D. We also show that the Trotter errors for 2D $XY$ and 1D Heisenberg chains are reduced, with respect to a digital decomposition, by a constant factor. In order to realize these important near-term speedups, we discuss the practical considerations needed to accurately characterize and calibrate our analog Hamiltonians for use in quantum simulations. We conclude with a discussion of how the Hamiltonian toggling techniques could be extended to derive new analog Hamiltonians which may be of use in more complex digital-analog quantum simulations for various models of interacting spins.


title: Planning Folding Motion with Simulation in the Loop Using Laser Forming Origami and Thermal Behaviors as an Example

abstract: Designing a robot or structure that can fold itself into a target shape is a process that involves challenges originated from multiple sources. For example, the designer of rigid self-folding robots must consider foldability from geometric and kinematic aspects to avoid self-intersection and undesired deformations. Recent works have shown success in estimating foldability of a design using robot motion planners. However, many foldable structures are actuated using physically coupled reactions (i.e., folding originated from thermal, chemical, or electromagnetic loads). Therefore, a reliable foldability analysis must consider additional constraints that resulted from these critical phenomena. This work investigates the idea of efficiently incorporating computationally expensive physics simulation within the folding motion planner to provide a better estimation of the foldability. In this paper, we will use laser forming origami as an example to demonstrate the benefits of considering the properties beyond geometry. We show that the design produced by the proposed method can be folded more efficiently.


title: Understanding the Distributional Aspects of Microcredit Expansions

abstract: Various poverty reduction strategies are being implemented in the pursuit of eliminating extreme poverty. One such strategy is increased access to microcredit in poor areas around the world. Microcredit, typically defined as the supply of small loans to underserved entrepreneurs that originally aimed at displacing expensive local money-lenders, has been both praised and criticized as a development tool (Banerjee et al., 2015b). This paper presents an analysis of heterogeneous impacts from increased access to microcredit using data from three randomised trials. In the spirit of recognising that in general the impact of a policy intervention varies conditional on an unknown set of factors, particular, we investigate whether heterogeneity presents itself as groups of winners and losers, and whether such subgroups share characteristics across RCTs. We find no evidence of impacts, neither average nor distributional, from increased access to microcredit on consumption levels. In contrast, the lack of average effects on profits seems to mask heterogeneous impacts. The findings are, however, not robust to the specific machine learning algorithm applied. Switching from the better performing Elastic Net to the worse performing Random Forest leads to a sharp increase in the variance of the estimates. In this context, methods to evaluate the relative performing machine learning algorithm developed by Chernozhukov et al. (2019) provide a disciplined way for the analyst to counter the uncertainty as to which algorithm to deploy.


title: Seismic Facies Analysis: A Deep Domain Adaptation Approach

abstract: Deep neural networks (DNNs) can learn accurately from large quantities of labeled input data, but often fail to do so when labelled data are scarce. DNNs sometimes fail to generalize ontest data sampled from different input distributions. Unsupervised Deep Domain Adaptation (DDA)techniques have been proven useful when no labels are available, and when distribution shifts are observed in the target domain (TD). In the present study, experiments are performed on seismic images of the F3 block 3D dataset from offshore Netherlands (source domain; SD) and Penobscot 3D survey data from Canada (target domain; TD). Three geological classes from SD and TD that have similar reflection patterns are considered. A deep neural network architecture named EarthAdaptNet (EAN) is proposed to semantically segment the seismic images when few classes have data scarcity, and we use a transposed residual unit to replace the traditional dilated convolution in the decoder block. The EAN achieved a pixel-level accuracy >84% and an accuracy of ~70% for the minority classes, showing improved performance compared to existing architectures. In addition, we introduce the CORAL (Correlation Alignment) method to the EAN to create an unsupervised deep domain adaptation network (EAN-DDA) for the classification of seismic reflections from F3 and Penobscot, to demonstrate possible approaches when labelled data are unavailable. Maximum class accuracy achieved was ~99% for class 2 of Penobscot, with an overall accuracy>50%. Taken together, the EAN-DDA has the potential to classify target domain seismic facies classes with high accuracy.

categories: physics.geo-ph cs.AI cs.CV cs.LG eess.IV


title: Comparison between the Thomas-Fermi and Hartree-Fock-Bogoliubov Methods in the Inner Crust of a Neutron Star: The Role of Pairing Correlations

abstract: We investigated the role of a pairing correlation in the chemical composition of the inner crust of a neutron star with the extended Thomas-Fermi method, using the Strutinsky integral correction. We compare our results with the fully self-consistent Hartree-Fock-Bogoliubov approach, showing that the resulting discrepancy, apart from the very low density region, is compatible with the typical accuracy we can achieve with standard mean-field methods.


title: Intrinsic Image Decomposition using Paradigms

abstract: Intrinsic image decomposition is the classical task of mapping image to albedo. The WHDR dataset allows methods to be evaluated by comparing predictions to human judgements ("lighter", "same as", "darker"). The best modern intrinsic image methods learn a map from image to albedo using rendered models and human judgements. This is convenient for practical methods, but cannot explain how a visual agent without geometric, surface and illumination models and a renderer could learn to recover intrinsic images. This paper describes a method that learns intrinsic image decomposition without seeing WHDR annotations, rendered data, or ground truth data. The method relies on paradigms - fake albedos and fake shading fields - together with a novel smoothing procedure that ensures good behavior at short scales on real images. Long scale error is controlled by averaging. Our method achieves WHDR scores competitive with those of strong recent methods allowed to see training WHDR annotations, rendered data, and ground truth data. Because our method is unsupervised, we can compute estimates of the test/train variance of WHDR scores; these are quite large, and it is unsafe to rely small differences in reported WHDR.


title: Optimal Quantum Thermometry with Coarse-grained Measurements

abstract: Precise thermometry for quantum systems is important to the development of new technology, and understanding the ultimate limits to precision presents a fundamental challenge. It is well known that optimal thermometry requires projective measurements of the total energy of the sample. However, this is infeasible in even moderately-sized systems, where realistic energy measurements will necessarily involve some coarse graining. Here, we explore the precision limits for temperature estimation when only coarse-grained measurements are available. Utilizing tools from signal processing, we derive the structure of optimal coarse-grained measurements and find that good temperature estimates can generally be attained even with a small number of outcomes. We apply our results to many-body systems and nonequilibrium thermometry. For the former, we focus on interacting spin lattices, both at and away from criticality, and find that the Fisher-information scaling with system size is unchanged after coarse-graining. For the latter, we consider a probe of given dimension interacting with the sample, followed by a measurement of the probe. We derive an upper bound on arbitrary, nonequilibrium strategies for such probe-based thermometry and illustrate it for thermometry on a Bose-Einstein condensate using an atomic quantum-dot probe.


title: A novel approach to the localization and the estimate of radioactivity in contaminated waste packages via imaging techniques

abstract: Dismantling nuclear power plants entails the production of a large amount of contaminated (or potentially contaminated) material whose disposal is of crucial importance. Most of the end products have to be stored in special repositories but it can happen that some of them are slightly contaminated or not contaminated at all, making it possible to free release them. One possible approach to free release measurements uses Large Clearance Monitors, chambers surrounded by plastic scintillation detectors that allow a decision about the clearance of waste packages up to 1000 kg. Due to the composite nature of the detectors in a Large Clearance Monitor, it is easy to imagine that one can apply 3D imaging algorithms to localize radioactive sources inside a waste package. In this work we will show how a special algorithm that maximizes the conditional informational entropy allows decisions about the clearance of portions of the sample.


title: A modified bond model for describing isotropic linear elastic material behaviour with the particle method

abstract: Particle based methods such as the Discrete Element Method and the Lattice Spring Method may be used for describing the behaviour of isotropic linear elastic materials. However, the common bond models employed to describe the interaction between particles restrict the range of Poisson's ratio that can be represented. In this paper, to overcome the restriction, a modified bond model that includes the coupling of shear strain energy of neighbouring bonds is proposed. The coupling is described by a multi-bond term that enables the model to distinguish between shear deformations and rigid-body rotations. The positive definiteness of the strain energy function of the modified bond model is verified. To validate the model, uniaxial tension, pure shear, pure bending and cantilever bending tests are performed. Comparison of the particle displacements with continuum mechanics solution demonstrates the ability of the model to describe the behaviour of isotropic linear elastic material for values of Poisson's ratio in the range $0 \leq \nu < 0.5$.


title: Mean field limit of Ensemble Square Root Filters -- discrete and continuous time

abstract: Consider the class of Ensemble Square Root filtering algorithms for the numerical approximation of the posterior distribution of nonlinear Markovian signals partially observed with linear observations corrupted with independent measurement noise. We analyze the asymptotic behavior of these algorithms in the large ensemble limit both in discrete and continuous time. We identify limiting mean-field processes on the level of the ensemble members, prove corresponding propagation of chaos results and derive associated convergence rates in terms of the ensemble size. In continuous time we also identify the stochastic partial differential equation driving the distribution of the mean-field process and perform a comparison with the Kushner-Stratonovich equation.


title: Quantifying Ionospheric Effects on Global 21-cm Observations

abstract: We modelled the two major layer of Earth's ionosphere, the F-layer and the D-layer, by a simplified spatial model with temporal variance to study the chromatic ionospheric effects on global 21-cm observations. From the analyses, we found that the magnitude of the ionospheric disruptions due to ionospheric refraction and absorption can be greater than the expected global 21-cm signal, and the variation of its magnitude can differ, depending on the ionospheric conditions. Within the parameter space adopted in the model, the shape of the global 21-cm signal is distorted after propagating through the ionosphere, while its amplitude is weakened. It is observed that the ionospheric effects do not cancel out over time, and thus should be accounted for in the foreground calibration at each timestep to account for the chromaticity introduced by the ionosphere.


title: COVID-19 and Mental Health/Substance Use Disorders on Reddit: A Longitudinal Study

abstract: COVID-19 pandemic has adversely and disproportionately impacted people suffering from mental health issues and substance use problems. This has been exacerbated by social isolation during the pandemic and the social stigma associated with mental health and substance use disorders, making people reluctant to share their struggles and seek help. Due to the anonymity and privacy they provide, social media emerged as a convenient medium for people to share their experiences about their day to day struggles. Reddit is a well-recognized social media platform that provides focused and structured forums called subreddits, that users subscribe to and discuss their experiences with others. Temporal assessment of the topical correlation between social media postings about mental health/substance use and postings about Coronavirus is crucial to better understand public sentiment on the pandemic and its evolving impact, especially related to vulnerable populations. In this study, we conduct a longitudinal topical analysis of postings between subreddits r/depression, r/Anxiety, r/SuicideWatch, and r/Coronavirus, and postings between subreddits r/opiates, r/OpiatesRecovery, r/addiction, and r/Coronavirus from January 2020 - October 2020. Our results show a high topical correlation between postings in r/depression and r/Coronavirus in September 2020. Further, the topical correlation between postings on substance use disorders and Coronavirus fluctuates, showing the highest correlation in August 2020. By monitoring these trends from platforms such as Reddit, epidemiologists, and mental health professionals can gain insights into the challenges faced by communities for targeted interventions.


title: The speed of random walk on Galton-Watson trees with vanishing conductances

abstract: In this paper we consider random walks on Galton-Watson trees with random conductances. On these trees, the distance of the walker to the root satisfies a law of large numbers with limit the effective velocity, or speed of the walk. We study the regularity of the speed as a function of the distribution of conductances, in particular when the distribution of conductances converges to a non-elliptic limit.


title: Rethinking Weight Decay For Efficient Neural Network Pruning

abstract: Introduced in the late 1980s for generalization purposes, pruning has now become a staple for compressing deep neural networks. Despite many innovations in recent decades, pruning approaches still face core issues that hinder their performance or scalability. Drawing inspiration from early work in the field, and especially the use of weight decay to achieve sparsity, we introduce Selective Weight Decay (SWD), which carries out efficient, continuous pruning throughout training. Our approach, theoretically grounded on Lagrangian smoothing, is versatile and can be applied to multiple tasks, networks, and pruning structures. We show that SWD compares favorably to state-of-the-art approaches, in terms of performance-to-parameters ratio, on the CIFAR-10, Cora, and ImageNet ILSVRC2012 datasets.


title: Zero Queueing for Multi-Server Jobs

abstract: Cloud computing today is dominated by multi-server jobs. These are jobs that request multiple servers simultaneously and hold onto all of these servers for the duration of the job. Multi-server jobs add a lot of complexity to the traditional one-job-per-server model: an arrival might not "fit" into the available servers and might have to queue, blocking later arrivals and leaving servers idle. From a queueing perspective, almost nothing is understood about multi-server job queueing systems; even understanding the exact stability region is a very hard problem. In this paper, we investigate a multi-server job queueing model under scaling regimes where the number of servers in the system grows. Specifically, we consider a system with multiple classes of jobs, where jobs from different classes can request different numbers of servers and have different service time distributions, and jobs are served in first-come-first-served order. The multi-server job model opens up new scaling regimes where both the number of servers that a job needs and the system load scale with the total number of servers. Within these scaling regimes, we derive the first results on stability, queueing probability, and the transient analysis of the number of jobs in the system for each class. In particular we derive sufficient conditions for zero queueing. Our analysis introduces a novel way of extracting information from the Lyapunov drift, which can be applicable to a broader scope of problems in queueing systems.


title: Scale estimation and data-driven tuning constant selection for M-quantile regression

abstract: M-quantile regression is a general form of quantile-like regression which usually utilises the Huber influence function and corresponding tuning constant. Estimation requires a nuisance scale parameter to ensure the M-quantile estimates are scale invariant, with several scale estimators having previously been proposed. In this paper we assess these scale estimators and evaluate their suitability, as well as proposing a new scale estimator based on the method of moments. Further, we present two approaches for estimating data-driven tuning constant selection for M-quantile regression. The tuning constants are obtained by i) minimising the estimated asymptotic variance of the regression parameters and ii) utilising an inverse M-quantile function to reduce the effect of outlying observations. We investigate whether data-driven tuning constants, as opposed to the usual fixed constant, for instance, at c=1.345, can improve the efficiency of the estimators of M-quantile regression parameters. The performance of the data-driven tuning constant is investigated in different scenarios using model-based simulations. Finally, we illustrate the proposed methods using a European Union Statistics on Income and Living Conditions data set.


title: General Solutions with Hidden Variables to the Abel Differential Equation

abstract: Using appropriate transformation, by the coefficient decomposition method, general solutions with hidden variables (parameters) to the Abel differential equation are obtained. In the process of solving, a set of entangled function pairs is discovered.


title: Delay Constrained Buffer-Aided Relay Selection in the Internet of Things with Decision-Assisted Reinforcement Learning

abstract: This paper investigates the reinforcement learning for the relay selection in the delay-constrained buffer-aided networks. The buffer-aided relay selection significantly improves the outage performance but often at the price of higher latency. On the other hand, modern communication systems such as the Internet of Things often have strict requirement on the latency. It is thus necessary to find relay selection policies to achieve good throughput performance in the buffer-aided relay network while stratifying the delay constraint. With the buffers employed at the relays and delay constraints imposed on the data transmission, obtaining the best relay selection becomes a complicated high-dimensional problem, making it hard for the reinforcement learning to converge. In this paper, we propose the novel decision-assisted deep reinforcement learning to improve the convergence. This is achieved by exploring the a-priori information from the buffer-aided relay system. The proposed approaches can achieve high throughput subject to delay constraints. Extensive simulation results are provided to verify the proposed algorithms.


title: The search for extratidal star candidates around Galactic globular clusters NGC 2808, NGC 6266, and NGC 6397 with {\textit{Gaia}} DR2 astrometry

abstract: Context. Extratidal stars are stellar bodies that end up outside the tidal radius of a cluster as a result of internal processes or external forces acting upon it. The presence and spatial distribution of these stars can give us insights into the past evolution of a cluster inside our Galaxy. Aims. Previous works suggest that globular clusters, when explored in detail, show evidence of extratidal stars. We aim to search for possible extratidal stars in the Galactic globular clusters NGC 6397, NGC 2808, and NGC 6266 using the photometry and proper motion measurements from Gaia DR2 database (Gaia Collaboration et al. 2018). Results. Finally, 120, 126, and 107 extratidal candidate stars were found lying outside the tidal radius of the globular clusters NGC 6397, NGC 2808, and NGC 6266, respectively. 70%, 25.4%, and 72.9% of the extratidal stars found are located outside the Jacobi radius of NGC 6397, NGC 2808, and NGC 6266, respectively. The spatial distribution of the extratidal stars belonging to NGC 6397 appears S-like, extending along the curved leading and trailing arms. NGC 2808 has an overdensity of stars in the trailing part of the cluster and NGC 6266 seems to have overdensities of extratidal stars in its eastern and northern sides. Conclusions. Proper motions and color-magnitude diagrams can be used to identify extratidal candidate stars around GCs. Nonetheless, depending on how different the kinematics and stellar populations of a cluster are compared to the Milky Way field, the fraction of contamination can be larger. All three clusters are found to have extratidal stars outside their tidal radii. For NGC 6397 and NGC 2808, these stars may be the result of a combined effect of the disc shocks and tidal disruptions. For NGC 6266, the distribution of extratidal stars is symmetrical around it, most likely indicating that the cluster has an extended stellar envelope.


title: A small actively-controlled high-resolution spectrograph based on off-the-shelf components

abstract: We present the design and testing of a prototype in-plane echelle spectrograph based on an actively controlled fibre-fed double-pass design. This system aims to be small and efficient with the minimum number of optical surfaces - currently a collimator/camera lens, cross-dispersing prism, grating and a reflector to send light to the detector. It is built from catalogue optical components and has dimensions of approximately 20x30 cm. It works in the optical regime with a resolution of >70,000. The spectrograph is fed by a bifurcated fibre with one fibre to a telescope and the other used to provide simultaneous Thorium Argon light illumination for wavelength calibration. The positions of the arc lines on the detector are processed in real time and commercial auto-guiding software is used to treat the positions of the arc lines as guide stars. The guiding software sends any required adjustments to mechanical piezo-electric actuators which move the mirror sending light to the camera removing any drift in the position of the arc lines. The current configuration using an sCMOS detector provides a precision of 3.5 milli-pixels equivalent to 4 m/s in a standard laboratory environment.


title: Multi-Scale Speaker Diarization With Neural Affinity Score Fusion

abstract: Identifying the identity of the speaker of short segments in human dialogue has been considered one of the most challenging problems in speech signal processing. Speaker representations of short speech segments tend to be unreliable, resulting in poor fidelity of speaker representations in tasks requiring speaker recognition. In this paper, we propose an unconventional method that tackles the trade-off between temporal resolution and the quality of the speaker representations. To find a set of weights that balance the scores from multiple temporal scales of segments, a neural affinity score fusion model is presented. Using the CALLHOME dataset, we show that our proposed multi-scale segmentation and integration approach can achieve a state-of-the-art diarization performance.


title: Space Lower Bounds for Graph Stream Problems

abstract: This work concerns with proving space lower bounds for graph problems in the streaming model. It is known that computing the length of shortest path between two nodes in the streaming model requires $\Omega(n)$ space, where $n$ is the number of nodes in the graph. We study the problem of finding the depth of a given node in a rooted tree in the streaming model. For this problem we prove a tight single pass space lower bound and a multipass space lower bound. As this is a special case of computing shortest paths on graphs, the above lower bounds also apply to the shortest path problem in the streaming model. The results are obtained by using known communication complexity lower bounds or by constructing hard instances for the problem. Additionally, we apply the techniques used in proving the above lower bound results to prove space lower bounds (single and multipass) for other graph problems like finding min $s-t$ cut, detecting negative weight cycle and finding whether two nodes lie in the same strongly connected component.


title: Computation capacities of a broad class of signaling networks are higher than their communication capacities

abstract: Due to structural and functional abnormalities or genetic variations and mutations, there may be dysfunctional molecules within an intracellular signaling network that do not allow the network to correctly regulate its output molecules, such as transcription factors. This disruption in signaling interrupts normal cellular functions and may eventually develop some pathological conditions. In this paper, computation capacity of signaling networks is introduced as a fundamental limit on signaling capability and performance of such networks. The computation capacity measures the maximum number of computable inputs, that is, the maximum number of input values for which the correct functional output values can be recovered from the erroneous network outputs, when the network contains some dysfunctional molecules. This contrasts with the conventional communication capacity that measures instead the maximum number of input values that can be correctly distinguished based on the erroneous network outputs. The computation capacity is higher than the communication capacity, if the network response function is not a one-to-one function of the input signals. By explicitly incorporating the effect of signaling errors that result in the network dysfunction, the computation capacity provides more information about the network and its malfunction. Two examples of signaling networks are studied here, one regulating caspase3 and another regulating NFkB, for which computation and communication capacities are analyzed. Higher computation capacities are observed for both networks. One biological implication of this finding is that signaling networks may have more capacity than that specified by the conventional communication capacity metric. The effect of feedback is also studied. In summary, this paper reports findings on a new fundamental feature of the signaling capability of cell signaling networks.


title: On barren plateaus and cost function locality in variational quantum algorithms

abstract: Variational quantum algorithms rely on gradient based optimization to iteratively minimize a cost function evaluated by measuring output(s) of a quantum processor. A barren plateau is the phenomenon of exponentially vanishing gradients in sufficiently expressive parametrized quantum circuits. It has been established that the onset of a barren plateau regime depends on the cost function, although the particular behavior has been demonstrated only for certain classes of cost functions. Here we derive a lower bound on the variance of the gradient, which depends mainly on the width of the circuit causal cone of each term in the Pauli decomposition of the cost function. Our result further clarifies the conditions under which barren plateaus can occur.


title: Gas compression and stellar feedback in the tidally interacting and ram-pressure stripped Virgo spiral galaxy NGC 4654

abstract: NGC 4654 is a Virgo galaxy seen almost face-on, which undergoes nearly edge-on gas ram pressure stripping and a fly-by gravitational interaction with another massive galaxy, NGC 4639. NGC 4654 shows a strongly compressed gas region near the outer edge of the optical disk, with HI surface densities (HSDR), exceeding the canonical value of 10-15 Msun/pc2. New IRAM 30m HERA CO(2-1) data of NGC 4654 are used to study the physical conditions of the ISM. The CO-to-H$_2$ conversion factor was estimated and found to be one to two times the Galactic value with significant decrease in the ratio between the molecular fraction and the total ISM pressure in the HSDR, self-gravitating gas, a Toomre parameter below $Q=1$ and star-formation efficiency 1.5-2 times higher. Analytical models were used to reproduce radial profiles of the SFR and the atomic and molecular surface densities. A Toomre parameter of $\rm Q \sim 0.8$ combined with an increase in the velocity dispersion of 5 km/s are necessary conditions to simultaneously reproduce the gas surface densities and the SFR. A dynamical model was used to reproduce the gas distribution of NGC 4654. The comparison between the velocity dispersion given by the moment 2 map and the intrinsic 3D velocity dispersion from the model were used to discriminate between regions of broader linewidths caused by a real increase in the velocity dispersion and those caused by an unresolved velocity gradient only. We found that the 5 km/s increase in the intrinsic velocity dispersion is compatible with observations. During a period of gas compression through external interactions, the gas surface density is enhanced, leading to an increased SFR and stellar feedback. Under the influence of stellar feedback, the gas density increases only moderately. The stellar feedback acts as a regulator of star-formation, increasing the turbulent velocity within the region.


title: The BDX-MINI detector for Light Dark Matter search at JLab

abstract: This paper describes the design and performance of a compact detector, BDX-MINI, that incorporates all features of a concept that optimized the detection of light dark matter produced by electrons in a beam dump. It represents a reduced version of the future BDX experiment expected to run at JLAB. BDX-MINI was exposed to penetrating particles produced by a 2.176 GeV electron beam incident on the beam dump of Hall A at Jefferson Lab. The detector consists of 30.5 kg of PbWO4 crystals with sufficient material following the beam dump to eliminate all known particles except neutrinos. The crystals are read out using silicon photomultipliers. Completely surrounding the detector are a passive layer of tungsten and two active scintillator veto systems, which are also read out using silicon photomultipliers. The design was validated and the performance of the robust detector was shown to be stable during a six month period during which the detector was operated with minimal access.


title: New insights into the use of Ultra Long Period Cepheids as cosmological standard candles

abstract: Ultra Long Period Cepheids (ULPs) are pulsating variable stars with a period longer than 80d and have been hypothesized to be the extension of the Classical Cepheids (CCs) at higher masses and luminosities. If confirmed as standard candles, their intrinsic luminosities, 1 to 3 mag brighter than typical CCs, would allow to reach the Hubble flow and, in turn, to determine the Hubble constant, H_0, in one step, avoiding the uncertainties associated with the calibration of primary and secondary indicators. To investigate the accuracy of ULPs as cosmological standard candles, we first collect all the ULPs known in the literature. The resulting sample includes 63 objects with a very large metallicity spread with 12 + log([O/H]) ranging from 7.2 to 9.2 dex. The analysis of their properties in the VI period-Wesenheit plane and in the color-magnitude diagram (CMD) supports the hypothesis that the ULPs are the extension of CCs at longer periods, higher masses and luminosities, even if, additional accurate and homogeneous data and a devoted theoretical scenario are needed to get firm conclusions. Finally, the three M31 ULPs, 8-0326, 8-1498 and H42, are investigated in more detail. For 8-1498 and H42, we cannot confirm their nature as ULPs, due to the inconsistency between their position in the CMD and the measured periods. For 8-0326, the light curve model fitting technique applied to the available time-series data allows us to constrain its intrinsic stellar parameters, distance and reddening.


title: Ambiguity through the lens of measure theory

abstract: In this paper, we establish a strong link between the ambiguity for finite words of a B\"uchi automaton and the ambiguity for infinite words of the same automaton. This link is based on measure theory. More precisely, we show that such an automaton is unambiguous, in the sense that no finite word labels two runs with the same starting state and the same ending state if and only if for each state, the set of infinite sequences labelling two runs starting from that state has measure zero. The measure used to define these negligible sets, that is sets of measure zero, can be any measure computed by a weighted automaton which is compatible with the B\"uchi automaton. This latter condition is very natural: the measure must put weight on cylinders [w] where w is the label of some run in the B\"uchi automaton.


title: On chiral responses to geometric torsion

abstract: We show that geometric torsion does not lead to new chiral dissipationless transport effects. Instead apparent response to torsion can be viewed as a manifestation of the chiral vortical effect.


title: Atom Cloud Detection Using a Deep Neural Network

abstract: We use a deep neural network to detect and place region-of-interest boxes around ultracold atom clouds in absorption and fluorescence images---with the ability to identify and bound multiple clouds within a single image. The neural network also outputs segmentation masks that identify the size, shape and orientation of each cloud from which we extract the clouds' Gaussian parameters. This allows 2D Gaussian fits to be reliably seeded thereby enabling fully automatic image processing.


title: Operando optical tracking of single-particle ion dynamics and phase transitions in battery electrodes

abstract: Key to advancing lithium-ion battery technology, and in particular fast charging capabilities, is our ability to follow and understand the dynamic processes occurring in operating materials under realistic conditions, in real time, and on the nano- to meso-scale. Currently, operando imaging of lithium-ion dynamics requires sophisticated synchrotron X-ray or electron microscopy techniques, which do not lend themselves to high-throughput material screening. This limits rapid and rational materials improvements. Here we introduce a simple lab-based, optical interferometric scattering microscope to resolve nanoscopic lithium-ion dynamics in battery materials and apply it to follow the repeated cycling of the archetypical cathode material Li$_\textit{x}$CoO$_2$. The method allows us to visualise directly the insulator-metal, solid solution and lithium ordering phase transitions in this material. We determine rates of lithium insertion and removal at the single-particle level and identify different mechanisms that occur on charge vs. discharge. Finally, we capture the dynamic formation of domain boundaries between different crystal orientations associated with the monoclinic lattice distortion at around Li$_{0.5}$CoO$_2$. The high throughput nature of our methodology allows many particles to be sampled across the entire electrode and, moving forward, will enable exploration of the role of dislocations, morphologies and cycling rate on battery degradation. The generality of our imaging concept means that it can be applied to study any battery electrode, and more broadly, systems where the transport of ions is associated with electronic or structural changes, including nanoionic films, ionic conducting polymers, photocatalytic materials and memristors.


title: Improving RNN-T ASR Accuracy Using Context Audio

abstract: We present a training scheme for streaming automatic speech recognition (ASR) based on recurrent neural network transducers (RNN-T) which allows the encoder network to learn to exploit context audio from a stream, using segmented or partially labeled sequences of the stream during training. We show that the use of context audio during training and inference can lead to word error rate reductions of more than 6% in a realistic production setting for a voice assistant ASR system. We investigate the effect of the proposed training approach on acoustically challenging data containing background speech and present data points which indicate that this approach helps the network learn both speaker and environment adaptation. To gain further insight into the ability of a long short-term memory (LSTM) based ASR encoder to exploit long-term context, we also visualize RNN-T loss gradients with respect to the input.


title: A sharp $L^{10}$ decoupling for the twisted cubic

abstract: We prove a sharp $l^{10}(L^{10})$ decoupling for the moment curve in $\mathbb{R}^3$. The proof involves a two-step decoupling combined with new incidence estimates for planks, tubes and plates.


title: Qubit-excitation-based adaptive variational quantum eigensolver

abstract: Molecular simulations with the variational quantum eigensolver (VQE) are a promising application for emerging noisy intermediate-scale quantum computers. Constructing accurate molecular ans\"atze that are easy to optimize and implemented by shallow quantum circuits is crucial for the successful implementation of such simulations. Ans\"atze are, generally, constructed as series of fermionic-excitation evolutions. Instead, we demonstrate the usefulness of constructing ans\"atze with "qubit-excitation evolutions", which, contrary to fermionic excitation evolutions, obey "qubit commutation relations". We show that qubit excitation evolutions, despite the lack of some of the physical features of fermionic excitation evolutions, accurately construct ans\"atze, while requiring asymptotically fewer gates. Utilizing qubit excitation evolutions, we introduce the qubit-excitation-based adaptive (QEB-ADAPT)-VQE protocol. The QEB-ADAPT-VQE is a modification of the ADAPT-VQE that performs molecular simulations using a problem-tailored ansatz, grown iteratively by appending evolutions of qubit excitation operators. By performing classical numerical simulations for small molecules, we benchmark the QEB-ADAPT-VQE, and compare it against the original fermionic-ADAPT-VQE and the qubit-ADAPT-VQE. In terms of circuit efficiency and convergence speed, we demonstrate that the QEB-ADAPT-VQE outperforms the qubit-ADAPT-VQE, which to our knowledge was the previous most circuit-efficient scalable VQE protocol for molecular simulations.


title: Performance vs. Spectral Properties For Single-Sideband Continuous Phase Modulation

abstract: This study revokes the performance of continuous phase modulation (CPM) able to generate a single-sideband (SSB) spectrum directly. This signal is analyzed in terms of modulation indices, pulse lengths, and pulse widths, all of which affect error probability, bandwidth, SSB property, and receiver complexity. The error probability performance is based on an approximation of the minimum Euclidean distance. A numerical power spectral density calculation for this particular SSB modulation in terms of modulation index is presented. Reasonable tradeoffs in designing modulation schemes have been proposed using multi-objective optimization to ensure sizable improvements in bit error rate (BER), spectral efficiencies, and complexity without losing the property of being a SSB signal. Performance comparisons are made with known CPM schemes, e.g., Gaussian Minimum Shift Keying (GMSK) and Raised Cosine based CPM (RC)


title: Existence and uniqueness of solutions to the time-dependent Kohn-Sham equations coupled with classical nuclear dynamics

abstract: We prove existence and uniqueness of solutions to the initial-value problem associated with a class of time-dependent Kohn-Sham equations coupled with Newtonian nuclear dynamics. We consider a pure power exchange term within a generalisation of the Local Density Approximation (LDA), identifying a range of exponents for the existence and uniqueness of $H^2$ solutions to the Kohn-Sham equations.


title: Search for trilepton resonances from chargino and neutralino pair production in $\sqrt{s}$ = 13 TeV $pp$ collisions with the ATLAS detector

abstract: A search is performed for the electroweak pair production of charginos and associated production of a chargino and neutralino, each of which decays through an $R$-parity-violating coupling into a lepton and a $W$, $Z$, or Higgs boson. The trilepton invariant-mass spectrum is constructed from events with three or more leptons, targeting chargino decays that include an electron or muon and a leptonically decaying $Z$ boson. The analyzed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collision data produced by the Large Hadron Collider at a center-of-mass energy of $\sqrt{s}$ = 13 TeV and collected by the ATLAS experiment between 2015 and 2018. The data are found to be consistent with predictions from the Standard Model. The results are interpreted as limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model. Limits are also set on the production of charginos and neutralinos for a Minimal Supersymmetric Standard Model with an approximate $B$-$L$ symmetry. Charginos and neutralinos with masses between 100 GeV and 1100 GeV are excluded depending on the assumed decay branching fractions into a lepton (electron, muon, or $\tau$-lepton) plus a boson ($W$, $Z$, or Higgs).


title: On Intersection Graph of Dihedral Group

abstract: Let $G$ be a finite group. The intersection graph of $G$ is a graph whose vertex set is the set of all proper non-trivial subgroups of $G$ and two distinct vertices $H$ and $K$ are adjacent if and only if $H\cap K \neq \{e\}$, where $e$ is the identity of the group $G$. In this paper, we investigate some properties and exploring some topological indices such as Wiener, Hyper-Wiener, first and second Zagreb, Schultz, Gutman and eccentric connectivity indices of the intersection graph of $D_{2n}$ for $n=p^2$, $p$ is prime. We also find the metric dimension and the resolving polynomial of the intersection graph of $D_{2p^2}$.


title: Two-Step Meta-Learning for Time-Series Forecasting Ensemble

abstract: Amounts of historical data collected increase and business intelligence applicability with automatic forecasting of time series are in high demand. While no single time series modeling method is universal to all types of dynamics, forecasting using an ensemble of several methods is often seen as a compromise. Instead of fixing ensemble diversity and size, we propose to predict these aspects adaptively using meta-learning. Meta-learning here considers two separate random forest regression models, built on 390 time-series features, to rank 22 univariate forecasting methods and recommend ensemble size. The forecasting ensemble is consequently formed from methods ranked as the best, and forecasts are pooled using either simple or weighted average (with a weight corresponding to reciprocal rank). The proposed approach was tested on 12561 micro-economic time-series (expanded to 38633 for various forecasting horizons) of M4 competition where meta-learning outperformed Theta and Comb benchmarks by relative forecasting errors for all data types and horizons. Best overall results were achieved by weighted pooling with a symmetric mean absolute percentage error of 9.21% versus 11.05% obtained using the Theta method.


title: A Function Based on Chebyshev Polynomials as an Alternative to the Sinc Function in FIR Filter Design

abstract: The sinc function is often used as the basis for the design of discrete linear-phase FIR filters. However the Fourier transform of the truncated sinc function exhibits ripple in the pass band due to the Gibbs phenomenon. This paper introduces an alternative function based on Chebyshev polynomials whose Fourier transform decreases monotonically in the pass band. Furthermore this function features an intrinsic window function with an adjustable parameter influencing the Fourier transform in the transition and stop bands.


title: The evolution of the gas fraction of quiescent galaxies modeled as a consequence of their creation rate

abstract: We discuss the evolution of the interstellar medium of quiescent galaxies, currently emerging from recent analyses, with the help of a simple model based on well-established empirical relations such as the stellar mass functions and the main sequence of star formation. This model is meant to describe observed quantities without making specific assumptions on the nature of quenching processes, but relying on their observable consequences. We find that the high gas fractions seen or suggested at high redshift in quiescent galaxies, and their apparent mild evolution at early times, can be mostly attributed to a progenitor effect where recently quenched galaxies with ~10% gas fractions dominate the quiescent galaxy population until z~1. In the same context, the much lower gas and dust fractions measured in local early-type galaxies are interpreted as the product of the steady depletion of their interstellar medium on a ~2 Gyr timescale, coupled with a higher fraction of more gas-exhaustive events.


title: Bubble deformation by a turbulent flow

abstract: We investigate the modes of deformation of an initially spherical bubble immersed in a homogeneous and isotropic turbulent background flow. We perform direct numerical simulations of the two-phase incompressible Navier-Stokes equations, considering a low-density bubble in the high density turbulent flow at various Weber number (the ratio of turbulent and surface tension forces) using the air-water density ratio. We discuss a theoretical framework for the bubble deformation in a turbulent flow using a spherical harmonic decomposition. We propose, for each mode of bubble deformation, a forcing term given by the statistics of velocity and pressure fluctuations, evaluated on a sphere of the same radius. This approach formally relates the bubble deformation and the background turbulent velocity fluctuations, in the limit of small deformations. The growth of the total surface deformation and of each individual mode is computed from the direct numerical simulations using an appropriate Voronoi decomposition of the bubble surface. We show that two successive temporal regimes occur: the first regime corresponds to deformations driven only by inertial forces, with the interface deformation growing linearly in time, in agreement with the model predictions, whereas the second regime results from a balance between inertial forces and surface tension. The transition time between the two regimes is given by the period of the first Rayleigh mode of bubble oscillation. We discuss how our approach can be used to relate the bubble lifetime to the turbulence statistics and eventually show that at high Weber number, bubble lifetime can be deduced from the statistics of turbulent fluctuations at the bubble scale.


title: Graph Signal Recovery Using Restricted Boltzmann Machines

abstract: We propose a model-agnostic pipeline to recover graph signals from an expert system by exploiting the content addressable memory property of restricted Boltzmann machine and the representational ability of a neural network. The proposed pipeline requires the deep neural network that is trained on a downward machine learning task with clean data, data which is free from any form of corruption or incompletion. We show that denoising the representations learned by the deep neural networks is usually more effective than denoising the data itself. Although this pipeline can deal with noise in any dataset, it is particularly effective for graph-structured datasets.


title: AZP: Automatic Specialization for Zero Values in Gaming Applications

abstract: Recent research has shown that dynamic zeros in shader programs of gaming applications can be effectively leveraged with a profile-guided, code-versioning transform. This transform duplicates code, specializes one path assuming certain key program operands, called versioning variables, are zero, and leaves the other path unspecialized. Dynamically, depending on the versioning variable's value, either the specialized fast path or the default slow path will execute. Prior work applied this transform manually and showed promising gains on gaming applications. In this paper, we present AZP, an automatic compiler approach to perform the above code-versioning transform. Our framework automatically determines which versioning variables or combinations of them are profitable, and determines the code region to duplicate and specialize (called the versioning scope). AZP takes operand zero value probabilities as input and it then uses classical techniques such as constant folding and dead-code elimination to determine the most profitable versioning variables and their versioning scopes. This information is then used to affect the final transform in a straightforward manner. We demonstrate that AZP is able to achieve an average speedup of 16.4% for targeted shader programs, amounting to an average frame-rate speedup of 3.5% across a collection of modern gaming applications on an NVIDIA GeForce RTX 2080 GPU GPU.


title: Mass transfer to freely suspended particles at high P\'eclet number

abstract: In a theoretical analysis, we generalise well known asymptotic results to obtain expressions for the rate of transfer of material from the surface of an arbitrary, rigid particle suspended in an open pathline flow at large P\'eclet number, $\textrm{Pe}$. The flow may be steady or periodic in time. We apply this result to numerically evaluate expressions for the surface flux to a freely suspended, axisymmetric ellipsoid (spheroid) in Stokes flow driven by a steady linear shear. We complement these analytical predictions with numerical simulations conducted over a range of $\textrm{Pe} = 10^1 - 10^4$ and confirm good agreement at large P\'eclet number. Our results allow us to examine the influence of particle shape upon the surface flux for a broad class of flows. When the background flow is irrotational, the surface flux is steady and is prescribed by three parameters only: the P\'eclet number, the particle aspect ratio and the strain topology. We observe that slender prolate spheroids tend to experience a higher surface flux compared to oblate spheroids with equivalent surface area. For rotational flows, particles may begin to spin or tumble, which may suppress or augment the convective transfer due to a realignment of the particle with respect to the strain field.


title: On the infinite Borwein product raised to a positive real power

abstract: In this paper, we study properties of the coefficients appearing in the $q$-series expansion of $\prod_{n\ge 1}[(1-q^n)/(1-q^{pn})]^\delta$, the infinite Borwein product for an arbitrary prime $p$, raised to an arbitrary positive real power $\delta$. We use the Hardy--Ramanujan--Rademacher circle method to give an asymptotic formula for the coefficients. For $p=3$ we give an estimate of their growth which enables us to partially confirm an earlier conjecture of the first author concerning an observed sign pattern of the coefficients when the exponent $\delta$ is within a specified range of positive real numbers. We further establish some vanishing and divisibility properties of the coefficients of the cube of the infinite Borwein product. We conclude with an Appendix presenting several new conjectures on precise sign patterns of infinite products raised to a real power which are similar to the conjecture we made in the $p=3$ case.


title: Mapping Gaia parallax systematic errors over the sky with faint Milky Way stars

abstract: Parallaxes measured by the Gaia mission have huge significance for astronomy, but parallaxes in Gaia DR2 are known to have systematic errors that depend on the source position and other quantities. We use the abundant information in faint Milky Way stars, along with the GOG simulation of the Gaia catalog, to probe the spatial dependence of Gaia DR2 parallax systematic errors in an empirical way. The parallax signal, concentrated in thick disk turnoff stars with magnitude G ~ 17, is sufficient to construct maps of the parallax systematic error over the majority of the sky. These maps show a locally regular "waffle pattern" on ~1 degree scales following Gaia scan directions, stronger linear "scar" features, and coherent variations on larger scales. The parallax bias maps also retain traces of astrophysical effects such as dust clouds. The waffle pattern, known from earlier maps of the Magellanic Clouds, extends over the entire sky; its local rms amplitude averages 15 microarcsec and varies by about a factor of two. The strength of this pattern increases by a factor ~6 from magnitude G = 13 to G = 20. Correlations with parallaxes of quasars and of stars with independent distance estimates support our bias estimates. Using similar methods, we map systematic errors in the proper motion and examine the relationship with the parallax systematics. We provide a code package to access and query our bias maps. Similar tests on the general stellar population should be useful in quantifying systematic errors in future Gaia releases.


title: Hund's coupling-assisted ferromagnetic percolation transition in a multiorbital flat band

abstract: By connecting Hund's physics with flat band physics, we establish an exact result for studying ferromagnetism in a multiorbital system. We consider a two-layer model consisting of a $p_x$, $p_y$-orbital honeycomb lattice layer and an $f$-orbital triangular lattice layer with sites aligned with the centers of the honeycomb plaquettes. The system features a flat band that admits a percolation representation for an appropriate chemical potential difference between the two layers. In this representation, the ground state space is spanned by maximum-spin clusters of localized single-particle states, and averaging over the ground states yields a correlated percolation problem with weights due to the spin degeneracy of the clusters. A paramagnetic-ferromagnetic transition occurs as the band approaches half filling and the ground states become dominated by states with a large maximum-spin cluster, as shown by Monte Carlo simulation.


title: New Spectral Analysis Results Within The Scope Of Extended Matter Research in The AR Lacertae Active Binary System

abstract: Within the scope of extended matter research, we present new spectral analysis results of an active binary system AR Lac. The low and high resolution spectra of this system were taken during the period 2013-2016. The evaluation of low dispersion spectra together with the B, V, Rc, Ic and WISE photometric data showed that AR Lac has an excess radiation in W2 band. In addition, the spectral energy distribution and the minima depth ratios of the light curves of this active binary system were studied to examine the flux contributions of the components of the system depending on wavelengths and on orbital phase. Furthermore, high resolution spectral analysis showed evidence of prominence-like structures and a possible extended matter around the cooler component of AR Lac binary system.


title: Scattering of Gravitons and Spinning Massive States from Compact Numerators

abstract: We provide a new efficient diagrammatic tool, in the context of the scattering equations, for computation of covariant $D$-dimensional tree-level $n$-point amplitudes with pairs of spinning massive particles using compact exponential numerators. We discuss how this framework allows non-integer spin extensions of recurrence relations for amplitudes developed for integer spin. Our results facilitate the on-going program for generating observables in classical general relativity from on-shell tree amplitudes through the Kawai-Lewellen-Tye relations and generalized unitarity.


title: Axisymmetric models for neutron star merger remnants with realistic thermal and rotational profiles

abstract: Merging neutron stars are expected to produce hot, metastable remnants in rapid differential rotation, which subsequently cool and evolve into rigidly rotating neutron stars or collapse to black holes. Studying this metastable phase and its further evolution is essential for the prediction and interpretation of the electromagnetic, neutrino, and gravitational signals from such a merger. In this work, we model binary neutron star merger remnants and propose new rotation and thermal laws that describe post-merger remnants. Our framework is capable to reproduce quasi-equilibrium configurations for generic equations of state, rotation and temperature profiles, including nonbarotropic ones. We demonstrate that our results are in agreement with numerical relativity simulations concerning bulk remnant properties like the mass, angular momentum, and the formation of a massive accretion disk. Because of the low computational cost for our axisymmetric code compared to full 3+1-dimensional simulations, we can perform an extensive exploration of the binary neutron star remnant parameter space studying several hundred thousand configurations for different equation of states.


title: Evidence for chromium hydride in the atmosphere of Hot Jupiter WASP-31b

abstract: Context. The characterisation of exoplanet atmospheres has shown a wide diversity in compositions. Hot Jupiters have the appropriate temperatures to host metallic compounds, which should be detectable through transmission spectroscopy. Aims. We aim to detect exotic species in transmission spectra of hot Jupiters, specifically WASP-31b, by testing a variety of chemical species to explain the spectrum. Methods. We conduct a re-analysis of publicly available transmission data of WASP-31b using the Bayesian retrieval framework TauREx II. We retrieve various combinations of the opacities of 25 species to determine the minimum set needed to fit the observed spectrum. Results. We report evidence for the spectroscopic signatures of chromium hydride (CrH), H2O and K in WASP-31b. Compared to a flat model without any signatures, a CrH-only model is preferred with a statistical significance of ~3.9 sigma. A model consisting of both CrH+H2O is found with ~2.6 sigma and ~3 sigma confidence over a CrH-only and H2O-only model respectively. Weak evidence for the addition of K is found at ~2.2 sigma over the H2O+CrH model, although the fidelity of the data point associated with this signature was questioned in earlier studies. Finally, the inclusion of collision-induced absorption and a Rayleigh scattering slope is found with ~3.5 sigma confidence over a flat model. This analysis presents the first evidence for signatures of CrH in a hot Jupiter atmosphere. At a retrieved temperature of ~1481 K, the atmosphere of WASP-31b is hot enough to host gaseous Cr-bearing species. The retrieved abundances agree well with predictions from thermal equilibrium chemistry and with the abundance found in an L-type brown dwarf atmosphere. However, additional retrievals using VLT FORS2 data lead to a non-detection of CrH. Future coverage with JWST has the potential to confirm the detection and discover other CrH features.


title: A look at the Hubble speed from first principles

abstract: We introduce a novel way of measuring $H_0$ from a combination of independent geometrical datasets, with no need of calibration nor of the choice of a cosmological model. We build on the {\it distance duality relation} which sets the ratio of the luminosity and angular diameter distance to a fixed scaling in redshift for any metric theory of gravity with standard photon propagation and constitutes a founding block of any theory describing our Universe. Our method provides the unprecedented possibility of determining $H_0$ from first principles, unleashing the measurement of this fundamental constant from calibration and assumption of a cosmological model. We find $H_0=69.5 \pm 1.7$ km/s/Mpc at $68\%$ C.L. showing that the Hubble constant can be constrained at percent level with minimal assumptions.


title: Influence of shear rate and surface chemistry on thrombus formation in micro-crevice

abstract: Thromboembolic complications remain a central issue in management of patients on mechanical circulatory support. Despite the best practices employed in design and manufacturing of modern ventricular assist devices, complexity and modular nature of these systems often introduces internal steps and crevices in the flow path which can serve as nidus for thrombus formation. Thrombotic potential is influenced by multiple factors including the characteristics of the flow and surface chemistry of the biomaterial. This study explored these elements in the setting of blood flow over a micro-crevice using a multi-constituent numerical model of thrombosis. The simulations reproduced the platelet deposition patterns observed experimentally and elucidated the role of flow, shear rate, and surface chemistry in shaping the deposition. The results offer insights for design and operation of blood-contacting devices.


title: Resolving the Dust-to-Metals Ratio and CO-to-H$_2$ Conversion Factor in the Nearby Universe

abstract: We investigate the relationship between the dust-to-metals ratio (D/M) and the local interstellar medium environment at ~2 kpc resolution in five nearby galaxies: IC342, M31, M33, M101, and NGC628. A modified blackbody model with a broken power-law emissivity is used to model the dust emission from 100 to 500 um observed by Herschel. We utilize the metallicity gradient derived from auroral line measurements in HII regions whenever possible. Both archival and new CO rotational line and HI 21 cm maps are adopted to calculate gas surface density, including new wide field CO and HI maps for IC342 from IRAM and the VLA, respectively. We experiment with several prescriptions of CO-to-H$_2$ conversion factor, and compare the resulting D/M-metallicity and D/M-density correlations, both of which are expected to be non-negative from depletion studies. The D/M is sensitive to the choice of the conversion factor. The conversion factor prescriptions based on metallicity only yield too much molecular gas in the center of IC342 to obtain the expected correlations. Among the prescriptions tested, the one that yields the expected correlations depends on both metallicity and surface density. The 1-$\sigma$ range of the derived D/M spans 0.40-0.58. Compared to chemical evolution models, our measurements suggest that the dust growth time scale is much shorter than the dust destruction time scale. The measured D/M is consistent with D/M in galaxy-integrated studies derived from infrared dust emission. Meanwhile, the measured D/M is systematically higher than the D/M derived from absorption, which likely indicates a systematic offset between the two methods.


title: MRAC-RL: A Framework for On-Line Policy Adaptation Under Parametric Model Uncertainty

abstract: Reinforcement learning (RL) algorithms have been successfully used to develop control policies for dynamical systems. For many such systems, these policies are trained in a simulated environment. Due to discrepancies between the simulated model and the true system dynamics, RL trained policies often fail to generalize and adapt appropriately when deployed in the real-world environment. Current research in bridging this sim-to-real gap has largely focused on improvements in simulation design and on the development of improved and specialized RL algorithms for robust control policy generation. In this paper we apply principles from adaptive control and system identification to develop the model-reference adaptive control & reinforcement learning (MRAC-RL) framework. We propose a set of novel MRAC algorithms applicable to a broad range of linear and nonlinear systems, and derive the associated control laws. The MRAC-RL framework utilizes an inner-loop adaptive controller that allows a simulation-trained outer-loop policy to adapt and operate effectively in a test environment, even when parametric model uncertainty exists. We demonstrate that the MRAC-RL approach improves upon state-of-the-art RL algorithms in developing control policies that can be applied to systems with modeling errors.


title: Long Short Term Memory Networks for Bandwidth Forecasting in Mobile Broadband Networks under Mobility

abstract: Bandwidth forecasting in Mobile Broadband (MBB) networks is a challenging task, particularly when coupled with a degree of mobility. In this work, we introduce HINDSIGHT++, an open-source R-based framework for bandwidth forecasting experimentation in MBB networks with Long Short Term Memory (LSTM) networks. We instrument HINDSIGHT++ following an Automated Machine Learning (AutoML) paradigm to first, alleviate the burden of data preprocessing, and second, enhance performance related aspects. We primarily focus on bandwidth forecasting for Fifth Generation (5G) networks. In particular, we leverage 5Gophers, the first open-source attempt to measure network performance on operational 5G networks in the US. We further explore the LSTM performance boundaries on Fourth Generation (4G) commercial settings using NYU-METS, an open-source dataset comprising of hundreds of bandwidth traces spanning different mobility scenarios. Our study aims to investigate the impact of hyperparameter optimization on achieving state-of-the-art performance and beyond. Results highlight its significance under 5G scenarios showing an average Mean Absolute Error (MAE) decrease of near 30% when compared to prior state-of-the-art values. Due to its universal design, we argue that HINDSIGHT++ can serve as a handy software tool for a multitude of applications in other scientific fields.


title: Free Mode Removal and Mode Decoupling for Simulating General Superconducting Quantum Circuits

abstract: Superconducting quantum circuits is one of the leading candidates for a universal quantum computer. Designing novel qubit and multiqubit superconducting circuits requires the ability to simulate and analyze the properties of a general circuit. In particular, going outside the transmon approach, we cannot make assumptions on anharmonicity, thus precluding blackbox quantization approaches and necessitating the formal circuit quantization approach. We consider and solve two issues involved in simulating general superconducting circuits. One of the issues is the handling of free modes in the circuit, that is, circuit modes with no potential term in the Hamiltonian. Another issue is circuit size, namely the challenge of simulating strongly coupled multimode circuits. The main mathematical tool we use to address these issues is the linear canonical transformation in the setting of quantum mechanics. We address the first issue by giving a provably correct algorithm for removing free modes by performing a linear canonical transformation to completely decouple the free modes from other circuit modes. We address the second by giving a series of different linear canonical transformations to reduce intermode couplings, thereby reducing the problem to the weakly coupled case and greatly mitigating the overhead for classical simulation. We benchmark our decoupling methods by applying them to the circuit of two inductively coupled fluxonium qubits, obtaining several orders of magnitude reduction in the size of the Hilbert space that needs to be simulated.


title: Gravitational Production of a Conformal Dark Sector

abstract: Dark sectors with purely gravitational couplings to the Standard Model are unavoidably populated from the SM plasma by graviton exchange, and naturally provide DM candidates. We examine the production in the relativistic regime where the dark sector is approximately scale invariant, providing general analytical formulas that depend solely on the central charge of the dark sector. We then assess the relevance of interactions that can lead to a variety of phenomena including thermalisation, nonperturbative mass gaps, out-of-equilibrium phase transitions and cannibalism in the dark sector. As an illustrative example we consider the dark glueball scenario in this light and show it to be a viable DM candidate due to the suppression of gravitational production. We go on to extend these results to strongly coupled CFTs and their holographic duals at large-N with the dark dilaton as the DM candidate.


title: Exploring Simple Siamese Representation Learning

abstract: Siamese networks have become a common structure in various recent models for unsupervised visual representation learning. These models maximize the similarity between two augmentations of one image, subject to certain conditions for avoiding collapsing solutions. In this paper, we report surprising empirical results that simple Siamese networks can learn meaningful representations even using none of the following: (i) negative sample pairs, (ii) large batches, (iii) momentum encoders. Our experiments show that collapsing solutions do exist for the loss and structure, but a stop-gradient operation plays an essential role in preventing collapsing. We provide a hypothesis on the implication of stop-gradient, and further show proof-of-concept experiments verifying it. Our "SimSiam" method achieves competitive results on ImageNet and downstream tasks. We hope this simple baseline will motivate people to rethink the roles of Siamese architectures for unsupervised representation learning. Code will be made available.


title: Extending rationally connected fibrations from ample subvarieties

abstract: Using deformation theory of rational curves, we prove a conjecture of Sommese on the extendability of morphisms from ample subvarieties when the morphism is a smooth (or mildly singular) fibration with rationally connected fibers. We apply this result in the context of Fano fibrations and prove a classification theorem for projective bundle and quadric fibration structures on ample subvarieties.


title: Learn to Bind and Grow Neural Structures

abstract: Task-incremental learning involves the challenging problem of learning new tasks continually, without forgetting past knowledge. Many approaches address the problem by expanding the structure of a shared neural network as tasks arrive, but struggle to grow optimally, without losing past knowledge. We present a new framework, Learn to Bind and Grow, which learns a neural architecture for a new task incrementally, either by binding with layers of a similar task or by expanding layers which are more likely to conflict between tasks. Central to our approach is a novel, interpretable, parameterization of the shared, multi-task architecture space, which then enables computing globally optimal architectures using Bayesian optimization. Experiments on continual learning benchmarks show that our framework performs comparably with earlier expansion based approaches and is able to flexibly compute multiple optimal solutions with performance-size trade-offs.


title: Vector computation

abstract: Quantum physical resources are directional quantities that can be formalized by unit vectors or the associated orthogonal projection operators. When compared to classical computational states which are elements of (power) sets vector computations offer (dis)advantages.


title: Scalable Local Timestepping on Octree Grids

abstract: Numerical solutions of hyperbolic partial differential equations(PDEs) are ubiquitous in science and engineering. Method of lines is a popular approach to discretize PDEs defined in spacetime, where space and time are discretized independently. When using explicit timesteppers on adaptive grids, the use of a global timestep-size dictated by the finest grid-spacing leads to inefficiencies in the coarser regions. Even though adaptive space discretizations are widely used in computational sciences, temporal adaptivity is less common due to its sophisticated nature. In this paper, we present highly scalable algorithms to enable local timestepping (LTS) for explicit timestepping schemes on fully adaptive octrees. We demonstrate the accuracy of our methods as well as the scalability of our framework across 16K cores in TACC's Frontera. We also present a speed up estimation model for LTS, which predicts the speedup compared to global timestepping (GTS) with an average of 0.1 relative error.


title: Efficient computation of baryon interpolating fields in Lattice QCD

abstract: In this work we present an efficient construction of baryon interpolating fields for lattice QCD computations of two and three point functions. These are essential building blocks of computations of nucleon parton distribution functions (PDFs), generalized parton distribution functions (GPDs) and transverse momentum dependent distributions functions (TMDs). Lattice QCD computations of these quantities can provide additional input to assist with the global fits on experimental data for determining TMDs, GPDs and PDFs.


title: Design Approach for Additive Manufacturing in Spare Part Supply Chains

abstract: In the current industrial revolution, additive manufacturing (AM) embodies a promising technology that can enhance the effectiveness, adaptability, and competitiveness of supply chains (SCs). Moreover, it facilitates the development of distributed SCs, thereby enhancing product availability, inventory levels, and lead time. However, the wide adoption of AM in industrial SCs creates various challenges, leading to new difficulties for SC design. In this context, this paper proposes a new design approach to AM SCs using optimization methods. More specifically, the proposed approach, comprising the p-median and mixed-integer linear programming models, considers the decision of deploying productive resources (3D printers) in specific locations of generic spare part SCs. The approach was evaluated in a real-world use case of an elevator maintenance service provider. The obtained results demonstrated the promising capabilities of the proposed design approach in managing the challenges arising from the forthcoming widespread use of 3D printers in manufacturing SCs.


title: Korn inequalities for incompatible tensor fields in three space dimensions with conformally invariant dislocation energy

abstract: Let $\Omega \subset \mathbb{R}^3$ be an open and bounded set with Lipschitz boundary and outward unit normal $\nu$. For $1<p<\infty$ we establish an improved version of the generalized $L^p$-Korn inequality for incompatible tensor fields $P$ in the new Banach space $$ W^{1,\,p,\,r}_0(\operatorname{dev}\operatorname{sym}\operatorname{Curl}; \Omega,\mathbb R^{3\times3}) = \{ P \in L^p(\Omega,\mathbb R^{3\times3}) \mid \operatorname{dev} \operatorname{sym} \operatorname{Curl} P \in L^r(\Omega,\mathbb R^{3\times3}),\ \operatorname{dev} \operatorname{sym} (P \times \nu) = 0 \text{ on $\partial \Omega$}\} $$ where $$ r \in [1, \infty), \qquad \frac1r \le \frac1p + \frac13, \qquad r >1 \quad \text{if $p = \frac32$.}$$ Specifically, there exists a constant $c=c(p,\Omega,r)>0$ such that the inequality \[ \|P \|_{L^p}\leq c\,\left(\|\operatorname{sym} P \|_{L^p} + \|\operatorname{dev}\operatorname{sym} \operatorname{Curl} P \|_{L^{r}}\right) \] holds for all tensor fields $P\in W^{1,\,p, \, r}_0(\operatorname{dev}\operatorname{sym}\operatorname{Curl})$. Here, $\operatorname{dev} X := X -\frac13 \operatorname{tr}(X)\,\mathbb{1}$ denotes the deviatoric (trace-free) part of a $3 \times 3$ matrix $X$ and the boundary condition is understood in a suitable weak sense.


title: S2D2: Small-scale Significant substructure DBSCAN Detection I. NESTs detection in 2D star-forming regions

abstract: The spatial and dynamical structure of star-forming regions can help provide insights on stellar formation patterns. The amount of data from current and upcoming surveys calls for robust and objective procedures to detect structure, so the results can be statistically analysed and different regions compared. We provide the community with a tool able to detect the small scale significant structure, above random expectation, in star-forming regions, which could be the imprint of the stellar formation process. The tool makes use of the one point correlation function and of nearest neighbour statistics to determine the parameters for the DBSCAN algorithm. The procedure successfully detects significant small scale substructures in heterogeneous regions, fulfilling the goals it was designed for, and providing very reliable structures. The analysis of regions close to complete spatial randomness ($Q \in [0.7,0.87]$) shows that, even when some structure is present and recovered, it is hardly distinguishable from spurious detection in homogeneous regions due to projection effects. Interpretation should thus be done with care. For concentrated regions, we detect a main structure surrounded by smaller ones, corresponding to the core plus some Poisson fluctuations around it. We argue that these structures do not correspond to the small compact regions we are looking for. In some realistic cases, a more complete hierarchical, multi-scale analysis would be needed to capture the complexity of the region. We have developed implementations of our procedure, and a catalogue of the NESTs (Nested Elementary STructures) detected by it in four star-forming regions (Taurus, IC 348, Upper Scorpius, and Carina), which are publicly available to the community. Implementations of the 3D, and up to 6D versions of the procedure including proper movements are in progress, and will be provided as future work.


title: Design of Experiments for Verifying Biomolecular Networks

abstract: There is a growing trend in molecular and synthetic biology of using mechanistic (non machine learning) models to design biomolecular networks. Once designed, these networks need to be validated by experimental results to ensure the theoretical network correctly models the true system. However, these experiments can be expensive and time consuming. We propose a design of experiments approach for validating these networks efficiently. Gaussian processes are used to construct a probabilistic model of the discrepancy between experimental results and the designed response, then a Bayesian optimization strategy used to select the next sample points. We compare different design criteria and develop a stopping criterion based on a metric that quantifies this discrepancy over the whole surface, and its uncertainty. We test our strategy on simulated data from computer models of biochemical processes.


title: Robust smoothed canonical correlation analysis for functional data

abstract: This paper provides robust estimators for the first canonical correlation and directions of random elements on Hilbert separable spaces by using robust association and scale measures combined with basis expansion and/or penalizations as a regularization tool. Under regularity conditions, the resulting estimators are consistent.


title: Deep learning insights into cosmological structure formation

abstract: The evolution of linear initial conditions present in the early universe into extended halos of dark matter at late times can be computed using cosmological simulations. However, a theoretical understanding of this complex process remains elusive; in particular, the role of anisotropic information in the initial conditions in establishing the final mass of dark matter halos remains a long-standing puzzle. Here, we build a deep learning framework to investigate this question. We train a three-dimensional convolutional neural network (CNN) to predict the mass of dark matter halos from the initial conditions, and quantify in full generality the amounts of information in the isotropic and anisotropic aspects of the initial density field about final halo masses. We find that anisotropies add a small, albeit statistically significant amount of information over that contained within spherical averages of the density field about final halo mass. However, the overall scatter in the final mass predictions does not change qualitatively with this additional information, only decreasing from 0.9 dex to 0.7 dex. Given such a small improvement, our results demonstrate that isotropic aspects of the initial density field essentially saturate the relevant information about final halo mass. Therefore, instead of searching for information directly encoded in initial conditions anisotropies, a more promising route to accurate, fast halo mass predictions is to add approximate dynamical information based e.g. on perturbation theory. More broadly, our results indicate that deep learning frameworks can provide a powerful tool for extracting physical insight into cosmological structure formation.


title: Completeness of the Gaia-verse III: using hidden states to infer gaps, detection efficiencies and the scanning law from the DR2 light curves

abstract: The completeness of the Gaia catalogues heavily depends on the status of that space telescope through time. Stars are only published with each of the astrometric, photometric and spectroscopic data products if they are detected a minimum number of times. If there is a gap in scientific operations, a drop in the detection efficiency or Gaia deviates from the commanded scanning law, then stars will miss out on potential detections and thus be less likely to make it into the Gaia catalogues. We lay the groundwork to retrospectively ascertain the status of Gaia throughout the mission from the tens of individual measurements of the billions of stars, by developing novel methodologies to infer both the orientation and angular velocity of Gaia through time and gaps and efficiency drops in the detections. We have applied these methodologies to the Gaia DR2 variable star epoch photometry -- which are the only publicly available Gaia time-series at the present time -- and make the results publicly available. We accompany these results with a new Python package scanninglaw (https://github.com/gaiaverse/scanninglaw) that you can use to easily predict Gaia observation times and detection probabilities for arbitrary locations on the sky.


title: D3-branes wrapped on a spindle

abstract: We construct supersymmetric $AdS_3\times \Sigma$ solutions of minimal gauged supergravity in $D=5$, where $\Sigma$ is a two-dimensional orbifold known as a spindle. Remarkably, these uplift on $S^5$, or more generally on any regular Sasaki-Einstein manifold, to smooth solutions of type IIB supergravity. The solutions are dual to $d=2$, $\mathcal{N}=(0,2)$ SCFTs and we show that the central charge for the gravity solution agrees with a field theory calculation associated with D3-branes wrapped on $\Sigma$. Unlike for smooth $\Sigma$ the superconformal R-symmetry mixes with the $U(1)$ isometry of the spindle.


title: Blandford-Znajek jets in galaxy formation simulations: method and implementation

abstract: Jets launched by active galactic nuclei (AGN) are believed to play a significant role in shaping the properties of galaxies and provide an energetically viable mechanism through which galaxies can become quenched. Here we present a novel AGN feedback model, which we have incorporated into the AREPO code, that evolves the black hole mass and spin as the accretion flow proceeds through a thin $\alpha$-disc which we self-consistently couple to a Blandford-Znajek jet. We apply our model to the central region of a typical radio-loud Seyfert galaxy embedded in a hot circumgalactic medium (CGM). We find that jets launched into high pressure environments thermalise efficiently due to the formation of recollimation shocks and the vigorous instabilities that these shocks excite increase the efficiency of the mixing of CGM and jet material. The beams of more overpressured jets, however, are not as readily disrupted by instabilities so the majority of the momentum flux at the jet base is retained out to the head, where the jet terminates in a reverse shock. All jets entrain a significant amount of cold circumnuclear disc material which, while energetically insignificant, dominates the lobe mass together with the hot, entrained CGM material. The jet power evolves significantly due to effective self-regulation by the black hole, fed by secularly-driven, intermittent mass flows. The direction of jets launched directly into the circumnuclear disc changes considerably due to effective Bardeen-Petterson torquing. Interestingly, these jets obliterate the innermost regions of the disc and drive large-scale, multi-phase, turbulent, bipolar outflows.


title: Detection of the Milky Way reflex motion due to the Large Magellanic Cloud infall

abstract: The Large Magellanic Cloud is the most massive satellite galaxy of the Milky Way, with an estimated mass exceeding a tenth of the mass of the Milky Way. Just past its closest approach of about 50 kpc, and flying by the Milky Way at an astonishing speed of 327 km/s, the Large Magellanic Cloud can affect our Galaxy in a number of ways, including dislodging the Milky Way disc from the Galactic centre-of-mass. Here, we report evidence that the Milky Way disc is moving with respect to stellar tracers in the outer halo ($40<r<120$ kpc) at $v_{\rm travel}=32^{+4}_{-4}$ km/s, in the direction $(\ell,b)_{\rm apex}=(56^{+9}_{-9},-34^{+10}_{-9})$ degrees, which points at an earlier location on the LMC trajectory. The resulting reflex motion is detected in the kinematics of outer halo stars and Milky Way satellite galaxies with accurate distances, proper motions and line-of-sight velocities. Our results indicate that dynamical models of our Galaxy cannot neglect gravitational perturbations induced by the Large Magellanic Cloud infall, nor can observations of the stellar halo be treated in a reference frame that does not correct for disc reflex motion. Future spectroscopic surveys of the stellar halo combined with Gaia astrometry will allow for sophisticated modelling of the Large Magellanic Cloud trajectory across the Milky Way, constraining the dark matter distribution in both galaxies with unprecedented detail.


title: Interstellar and Circumgalactic Properties of an Unseen $z=6.84$ Galaxy: Abundances, Ionization, and Heating in the Earliest Known Quasar Absorber

abstract: We analyze relative abundances and ionization conditions in a strong absorption system at z=6.84, seen in the spectrum of the z=7.54 background quasar ULAS J134208.10+092838.61. Singly ionized C, Si, Fe, Mg, and Al measurements are consistent with a warm neutral medium that is metal-poor but not chemically pristine. Firm non-detections of C IV and Si IV imply that any warm ionized phase of the IGM or CGM has not yet been enriched past the ultra-metal-poor regime (<0.001Z_{solar}), unlike lower redshift DLAs where these lines are nearly ubiquitous. Relative abundances of the heavy elements 794 Myr after the Big Bang resemble those of metal-poor damped Lyman Alpha systems at intermediate redshift and Milky Way halo stars, and show no evidence of enhanced [alpha/Fe], [C/Fe] or other signatures of yields dominated by massive stars. A detection of the CII* fine structure line reveals local sources of excitation from heating, beyond the level of photo-excitation supplied by the CMB. We estimate the total and [CII] cooling rates, balancing against ISM heating sources to develop an heuristic two-phase model of the neutral medium. The implied heating requires a surface density of star formation slightly exceeding that of the Milky Way but not at the level of a strong starburst. For a typical (assumed) NHI=10^{20.6}, an abundance of [Fe/H]=-2.2 matches the columns of species in the neutral phase. To remain undetected in C IV, a warm ionized phase would either need much lower [C/H]<-4.2 over an absorption path of 1 kpc, or else a very small absorption path (a few pc). While still speculative, these results suggest a significant reduction in heavy element enrichment outside of neutral star forming regions of the ISM, as would be expected in early stages of galactic chemical evolution.