Search

We construct time-evolving gravitational potential models for a Milky Way-mass galaxy from the FIRE-2 suite of cosmological-baryonic simulations using basis function expansions. These models capture the angular variation with spherical harmonics for the halo and azimuthal harmonics for the disk, and the radial or meridional plane variation with splines. We fit low-order expansions (4 angular/harmonic terms) to the galaxy's potential for each snapshot, spaced roughly 25 Myr apart, over the last 4 Gyr of its evolution, then extract the forces at discrete times and interpolate them between adjacent snapshots for forward orbit integration. Our method reconstructs the forces felt by simulation particles with high fidelity, with 95% of both stars and dark matter, outside of self-gravitating subhalos, exhibiting errors $\leq4\%$ in both the disk and the halo. Imposing symmetry on the model systematically increases these errors, particularly for disk particles, which show greater sensitivity to imposed symmetries. The majority of orbits recovered using the models exhibit positional errors $\leq10\%$ for 2-3 orbital periods, with higher errors for orbits that spend more time near the galactic center. Approximate integrals of motion are retrieved with high accuracy even with a larger potential sampling interval of 200 Myr. After 4 Gyr of integration, 43% and 70% of orbits have total energy and angular momentum errors within 10%, respectively. Consequently, there is higher reliability in orbital shape parameters such as pericenters and apocenters, with errors $\sim10\%$ even after multiple orbital periods. These techniques have diverse applications, including studying satellite disruption in cosmological contexts., Comment: 28 pages, 13 figures, and 4 tables. Submitted to APJ

A $^6$Li-doped plastic scintillator with pulse-shape discrimination capabilities, commercially identified as EJ-299-50, has been developed and produced at the kilogram-scale. A total of 44 kg-scale bars of dimensions 5.5 cm $\times$ 5.5 cm $\times$ 50 cm of this material have been characterized. Optical properties like light output and effective attenuation length have been found to be comparable to $^6$Li-doped liquid scintillators. The scintillator EJ-299-50 shows good neutron detection capabilities with an effective efficiency for capture on $^6$Li of approximately 85%. Stability tests performed on two formulation variations showed no intrinsic degradation in the material or optical properties during several months of observations.

In the cold dark matter paradigm, our Galaxy is predicted to contain >10000 dark matter subhaloes in the $10^5-10^8M_\odot$ range which should be completely devoid of stars. Stellar streams are sensitive to the presence of these subhaloes, which can create small-scale features in streams if they pass closely enough. Modelling these encounters can therefore, potentially recover the subhalo's properties. In this work, we demonstrate this for streams generated in numerical simulations, modelled on eccentric orbits in a realistic Milky Way potential, which includes the Large Magellanic Cloud and the subhalo itself. We focus on a mock model of the ATLAS-Aliqa Uma stream and inject a $10^7 M_\odot$ subhalo, creating a similar discontinuous morphology to current observations. We then explore how well subhalo properties are recovered using mock stream observations, consisting of no observational errors, as well as assuming realistic observational setups. These setups include present day style observations, and what will be possible with 4MOST and Gaia DR5 in the future. We show that we can recover all parameters describing the impact even with uncertainties matching existing data, including subhalo positions, velocities, mass and scale radius. Modelling the subhalo on an orbit instead of assuming an impulse approximation, we greatly reduce the degeneracy between subhalo mass and velocity seen in previous works. However, we find a slight bias in the subhalo mass (~0.1 dex). This demonstrates that we should be able to reliably extract the properties of subhaloes with stellar streams in the near future., Comment: 13 pages, 14 figures. Submitted to MMRAS. Comments welcome!

We examine the gamma-ray signal from dark matter (DM) annihilation from analogues of the Sagittarius (Sgr) dwarf spheroidal galaxy in the Auriga cosmological simulations. For velocity-dependent annihilation cross sections, we compute emissions from simulated Sgr subhalos and from the Milky Way (MW) foreground. In addition to the annihilation signals from DM particles bound to Sgr, we consider for the first time the annihilation of DM particles bound to the MW that overlap spatially with Sgr. For p-wave models this contribution can enhance the signal by over an order of magnitude, while for d-wave models the enhancement can be over two orders of magnitude. For Sommerfeld models, the corresponding emission decreases by up to nearly an order of magnitude. For Sommerfeld and s-wave models, the Sgr source can be visible above the MW foreground emission, while for p and d-wave models, the signal towards Sgr is most likely dominated by foreground MW emission. We interpret our results within the context of the observed gamma-ray emission from Sgr, finding that the templates from simulations likely have spatial morphology that is too extended to explain the point-like emission that is observed., Comment: 16 pages, 5 figures, 2 tables, to be submitted to JCAP

The Hercules ultra-faint dwarf galaxy (UFD) has long been hypothesized to be tidally disrupting, yet no conclusive evidence has been found for tidal disruption owing partly to difficulties in identifying Hercules member stars. In this work, we present a homogeneous re-analysis of new and existing observations of Hercules, including the detection of a new potential member star located $\sim$1 $^{\circ}$ ($\sim1.7$ kpc) west of the center of the system. In addition to measuring the line-of-sight velocity gradient, we compare predictions from dynamical models of stream formation to these observations. We report an updated velocity dispersion measurement based on 28 stars, $1.9^{+0.6}_{-0.6}$ km s$^{\rm -1}$, which is significantly lower than previous measurements. We find that the line-of-sight velocity gradient is $1.8^{+1.8}_{-1.8}$ km s$^{\rm -1}$ kpc$^{\rm -1}$ along the major axis of Hercules, consistent with zero within 1 $\sigma$. Our dynamical models of stream formation, on the other hand, can reproduce the morphology of the Hercules UFD, specifically the misalignment between the elongation and the orbital motion direction. Additionally, these dynamical models indicate that any radial velocity gradient from tidal disruption would be too small, $0.00^{+0.97}_{-0.91}$ km s$^{\rm -1}$ kpc$^{\rm -1}$, to be detectable with current sample sizes. Combined with our analysis of the tidal radius evolution of the system as a function of its orbital phase, we argue that it is likely that Hercules is indeed currently undergoing tidal disruption in its extended stellar halo with a line-of-sight velocity gradient too small to be detected with current observational datasets., Comment: 21 pages, 11 figures. Accepted for publication in ApJ

The third data release (DR3) of Gaia has provided a five-fold increase in the number of radial velocity measurements of stars, as well as a stark improvement in parallax and proper motion measurements. To help with studies that seek to test models and interpret Gaia DR3, we present nine Gaia synthetic surveys, based on three solar positions in three Milky-Way-mass galaxies of the Latte suite of the Fire-2 cosmological simulations. These synthetic surveys match the selection function, radial velocity measurements, and photometry of Gaia DR3, adapting the code base Ananke, previously used to match the Gaia DR2 release in Sanderson et al. 2020. The synthetic surveys are publicly available and can be found at http://ananke.hub.yt/. Similarly to the previous release of Ananke, these surveys are based on cosmological simulations and thus able to model non-equilibrium dynamical effects, making them a useful tool in testing and interpreting Gaia DR3., Comment: 17 pages, 7 tables, 6 figures

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe., Comment: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

We present a 6-D map of the Orphan-Chenab (OC) stream by combining the data from Southern Stellar Stream Spectroscopic Survey ($S^5$) and {\it Gaia}. We reconstruct the proper motion, radial velocity, distance, on-sky track and stellar density along the stream with spline models. The stream has a total luminosity of $M_V=-8.2$ and metallicity of $\mathrm{[Fe/H]}=-1.9$, similar to classical Milky Way (MW) satellites like Draco. The stream shows drastic changes in its physical width varying from 200 pc to 1 kpc, but a constant line of sight velocity dispersion of 5 km/ss. Despite the large apparent variation in the stellar number density along the stream, the flow rate of stars along the stream is remarkably constant. We model the 6-D stream track by a Lagrange-point stripping method with a flexible MW potential in the presence of a moving extended Large Magellanic Cloud (LMC). This allows us to constrain the mass profile of the MW within the distance range 15.6 < r < 55.5 kpc, with the best measured enclosed mass of $(2.85\pm 0.1)\times 10^{11}\,M_\odot$ within 32.4 kpc. Our stream measurements are highly sensitive to the LMC mass profile with the most precise measurement of its enclosed mass made at 32.8 kpc, $(7.02\pm 0.9)\times10^{10}\, {\rm M}_\odot$. We also detect that the LMC dark matter halo extends to at least 53 kpc. The fitting of the OC stream allows us to constrain the past LMC trajectory and the degree of dynamical friction it experienced. We demonstrate that the stars in the OC stream show large energy and angular momentum spreads caused by LMC perturbation., Comment: accepted to MNRAS; data released with the paper is available on Zenodo https://doi.org/10.5281/zenodo.7222654

We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to > 110 kpc) and surviving only at larger apocenters (> 40 kpc) than those observed in the Milky Way. We find that the population of high-stellar mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as satellites in DES-like observations, since their tidal tails are too low-surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low surface brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations., Comment: 24 pages, 13 figures, 3 tables, submitted to ApJ

It has recently been shown that the Large Magellanic Cloud (LMC) has a substantial effect on the Milky Way's stellar halo and stellar streams. Here, we explore how deformations of the Milky Way and LMC's dark matter haloes affect stellar streams, and whether these effects are observable. In particular, we focus on the Orphan-Chenab (OC) stream which passes particularly close to the LMC, and spans a large portion of the Milky Way's halo. We represent the Milky Way--LMC system using basis function expansions that capture their evolution in an $N$-body simulation. We present the properties of this system, such as the evolution of the densities and force fields of each galaxy. The OC stream is evolved in this time-dependent, deforming potential, and we investigate the effects of the various moments of the Milky Way and the LMC. We find that the simulated OC stream is strongly influenced by the deformations of both the Milky Way and the LMC, and that this effect is much larger than current observational errors. In particular, the Milky Way dipole has the biggest impact on the stream, followed by the evolution of the LMC's monopole, and the LMC's quadrupole. Detecting these effects would confirm a key prediction of collisionless, cold dark matter, and would be a powerful test of alternative dark matter and alternative gravity models., Comment: 17 pages, 10 figures. Accepted for publication in MNRAS

The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics., Comment: White paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). 88 pages, 9 figures. Comments welcome

Machine Learning is becoming a pivotal aspect of many systems today, offering newfound performance on classification and prediction tasks, but this rapid integration also comes with new unforeseen vulnerabilities. To harden these systems the ever-growing field of Adversarial Machine Learning has proposed new attack and defense mechanisms. However, a great asymmetry exists as these defensive methods can only provide security to certain models and lack scalability, computational efficiency, and practicality due to overly restrictive constraints. Moreover, newly introduced attacks can easily bypass defensive strategies by making subtle alterations. In this paper, we study an alternate approach inspired by honeypots to detect adversaries. Our approach yields learned models with an embedded watermark. When an adversary initiates an interaction with our model, attacks are encouraged to add this predetermined watermark stimulating detection of adversarial examples. We show that HoneyModels can reveal 69.5% of adversaries attempting to attack a Neural Network while preserving the original functionality of the model. HoneyModels offer an alternate direction to secure Machine Learning that slightly affects the accuracy while encouraging the creation of watermarked adversarial samples detectable by the HoneyModel but indistinguishable from others for the adversary., Comment: Published in: MILCOM 2021 - 2021 IEEE Military Communications Conference (MILCOM)

We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ${\sim}10-50$ kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass-metallicity relation, the luminosities of the progenitors of the DG streams range between Carina and Ursa Major I ($-9.5\lesssim M_V\lesssim-5.5$). Four of the six GC streams have mean metallicities of [Fe/H]$< -2$, more metal-poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the 300S and Jet GC streams are the only streams on retrograde orbits in our dozen stream sample. Finally, we compare the orbital properties of the streams with known DGs and GCs in the MW, finding several possible associations. Some streams appear to have been accreted with the recently discovered Gaia-Enceladus-Sausage system, and others suggest that GCs were formed in and accreted together with the progenitors of DG streams whose stellar masses are similar to Draco to Carina ($\sim10^5-10^6M_\odot$)., Comment: 24 pages, 8 Figures, 2 Tables. Accepted for publication in ApJ

We use six years of data from the Dark Energy Survey to perform a detailed photometric characterization of the Phoenix stellar stream, a 15-degree long, thin, dynamically cold, low-metallicity stellar system in the southern hemisphere. We use natural splines, a non-parametric modeling technique, to simultaneously fit the stream track, width, and linear density. This updated stream model allows us to improve measurements of the heliocentric distance ($17.4 \pm 0.1\,{\rm (stat.)} \pm 0.8\,{\rm (sys.)}$ kpc) and distance gradient ($-0.009 \pm 0.006$ kpc deg$^{-1}$) of Phoenix, which corresponds to a small change of $0.13 \pm 0.09$ kpc in heliocentric distance along the length of the stream. We measure linear intensity variations on degree scales, as well as deviations in the stream track on $\sim 2$-degree scales, suggesting that the stream may have been disturbed during its formation and/or evolution. We recover three peaks and one gap in linear intensity along with fluctuations in the stream track. Compared to other thin streams, the Phoenix stream shows more fluctuations and, consequently, the study of Phoenix offers a unique perspective on gravitational perturbations of stellar streams. We discuss possible sources of perturbations to Phoenix including baryonic structures in the Galaxy and dark matter subhalos., Comment: 18 Pages, 7 Figures, Matches published version

Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (${S}^5$), Gaia EDR3, and the Dark Energy Survey (DES), to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from 14 to $19 \times 10^{10}\ \mathrm{M}_{\odot}$, probed over a range of closest approach times and distances. With the most constraining stream (Orphan-Chenab), we measure an LMC mass of $18.8^{+ 3.5}_{- 4.0} \times 10^{10}\ \mathrm{M}_{\odot}$, probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC's influence on structures in the Milky Way. Using this sample of streams, we find that the LMC's effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on the impulse approximation and we show that the LMC's effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick., Comment: 23 pages, 10 figures, submitted to ApJ

We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia EDR3 proper motions, detects a clear velocity gradient in Antlia 2 and a tentative velocity gradient in Crater 2. The velocity gradient magnitudes and directions are consistent with particle stream simulations of tidal disruption. Furthermore, the orbit and kinematics of Antlia 2 require a model that includes the reflex motion of the Milky Way induced by the Large Magellanic Cloud. We also find that Antlia 2's metallicity was previously overestimated, so it lies on the empirical luminosity-metallicity relation and is likely only now experiencing substantial stellar mass loss. Current dynamical models of Antlia 2 require it to have lost over 90% of its stars to tides, in tension with the low stellar mass loss implied by the updated metallicity. Overall, the new kinematic measurements support a tidal disruption scenario for the origin of these large and extended dwarf spheroidal galaxies., Comment: 23 pages (+appendix), 11 figures, accepted to ApJ. Updated to accepted version (1 new figure, minor changes to text) Video of Antlia 2 tidal disruption here: https://www.youtube.com/watch?v=wvt-Q4kRq90

The recently discovered Indus stellar stream exhibits a diverse chemical signature compared to what is found for most other streams due to the abundances of two outlier stars, Indus$\_$0 and Indus$\_$13. Indus$\_$13, exhibits an extreme enhancement in rapid neutron-capture ($r$-)process elements with $\mathrm{[Eu/Fe]} = +1.81$. It thus provides direct evidence of the accreted nature of $r$-process enhanced stars. In this paper we present a detailed chemical analysis of the neutron-capture elements in Indus$\_$13, revealing the star to be slightly actinide poor. The other outlier, Indus$\_0$, displays a globular cluster-like signature with high N, Na, and Al abundances, while the rest of the Indus stars show abundances compatible with a dwarf galaxy origin. Hence, Indus$\_0$ provides the first chemical evidence of a fully disrupted dwarf containing a globular cluster. We use the chemical signature of the Indus stars to discuss the nature of the stream progenitor which was likely a chemically evolved system, with a mass somewhere in the range from Ursa Minor to Fornax., Comment: 14 pages, 4 figures, and 5 tables, accepted for publication in ApJ

We perform a detailed photometric and astrometric analysis of stars in the Jet stream using data from the first data release of the DECam Local Volume Exploration Survey (DELVE) DR1 and \emph{Gaia} EDR3. We discover that the stream extends over $\sim 29 ^{\circ}$ on the sky (increasing the known length by $18^{\circ}$), which is comparable to the kinematically cold Phoenix, ATLAS, and GD-1 streams. Using blue horizontal branch stars, we resolve a distance gradient along the Jet stream of 0.2 kpc/deg, with distances ranging from $D_\odot \sim 27-34$ kpc. We use natural splines to simultaneously fit the stream track, width, and intensity to quantitatively characterize density variations in the Jet stream, including a large gap, and identify substructure off the main track of the stream. Furthermore, we report the first measurement of the proper motion of the Jet stream and find that it is well-aligned with the stream track suggesting the stream has likely not been significantly perturbed perpendicular to the line of sight. Finally, we fit the stream with a dynamical model and find that the stream is on a retrograde orbit, and is well fit by a gravitational potential including the Milky Way and Large Magellanic Cloud. These results indicate the Jet stream is an excellent candidate for future studies with deeper photometry, astrometry, and spectroscopy to study the potential of the Milky Way and probe perturbations from baryonic and dark matter substructure., Comment: 19 pages, 10 figures, Accepted in ApJ