Search

The Large Magellanic Cloud (LMC) is a Milky Way (MW) satellite that is massive enough to gravitationally attract the MW disc and inner halo, causing significant motion of the inner MW with respect to the outer halo. In this work, we probe this interaction by constructing a sample of 9,866 blue horizontal branch (BHB) stars with radial velocities from the DESI spectroscopic survey out to 120 kpc from the Galactic centre. This is the largest spectroscopic set of BHB stars in the literature to date, and it contains four times more stars with Galactocentric distances beyond 50 kpc than previous BHB catalogues. Using the DESI BHB sample combined with SDSS BHBs, we measure the bulk radial velocity of stars in the outer halo and observe that the velocity in the Southern Galactic hemisphere is different by 3.7$\sigma$ from the North. Modelling the projected velocity field shows that its dipole component is directed at a point 22 degrees away from the LMC along its orbit, which we interpret as the travel direction of the inner MW. The velocity field includes a monopole term that is -24 km/s, which we refer to as compression velocity. This velocity is significantly larger than predicted by the current models of the MW and LMC interaction. This work uses DESI data from its first two years of observations, but we expect that with upcoming DESI data releases, the sample of BHB stars will increase and our ability to measure the MW-LMC interaction will improve significantly., Comment: 22 pages, 19 figures. Submitted to MNRAS

In a hierarchically formed Universe, galaxies accrete smaller systems that tidally disrupt as they evolve in the host's potential. We present a complete catalogue of disrupting galaxies accreted onto Milky Way-mass haloes from the Auriga suite of cosmological magnetohydrodynamic zoom-in simulations. We classify accretion events as intact satellites, stellar streams, or phase-mixed systems based on automated criteria calibrated to a visually classified sample, and match accretions to their counterparts in haloes re-simulated at higher resolution. Most satellites with a bound progenitor at the present day have lost substantial amounts of stellar mass -- 67 per cent have $f_\text{bound} < 0.97$ (our threshold to no longer be considered intact), while 53 per cent satisfy a more stringent $f_\text{bound} < 0.8$. Streams typically outnumber intact systems, contribute a smaller fraction of overall accreted stars, and are substantial contributors at intermediate distances from the host centre ($\sim$0.1 to $\sim$0.7$R_\text{200m}$, or $\sim$35 to $\sim$250 kpc for the Milky Way). We also identify accretion events that disrupt to form streams around massive intact satellites instead of the main host. Streams are more likely than intact or phase-mixed systems to have experienced preprocessing, suggesting this mechanism is important for setting disruption rates around Milky Way-mass haloes. All of these results are preserved across different simulation resolutions, though we do find some hints that satellites disrupt more readily at lower resolution. The Auriga haloes suggest that disrupting satellites surrounding Milky Way-mass galaxies are the norm and that a wealth of tidal features waits to be uncovered in upcoming surveys., Comment: 16+4 pages, 13+2 figures, 1+1 tables. Submitted to MNRAS, comments are welcome

Galaxies like the Milky Way are surrounded by complex populations of satellites at all stages of tidal disruption. In this paper, we present a dynamical study of the disrupting satellite galaxies in the Auriga simulations that are orbiting 28 distinct Milky Way-mass hosts across three resolutions. We find that the satellite galaxy populations are highly disrupted. The majority of satellites that remain fully intact at present day were accreted recently without experiencing more than one pericentre ($n_{\rm peri} \lesssim 1$) and have large apocentres ($r_{\rm apo} \gtrsim 200$ kpc) and pericentres ($r_{\rm peri} \gtrsim 50$ kpc). The remaining satellites have experienced significant tidal disruption and, given full knowledge of the system, would be classified as stellar streams. We find stellar streams in Auriga across the range of pericentres and apocentres of the known Milky Way dwarf galaxy streams and, interestingly, overlapping significantly with the Milky Way intact satellite population. We find no significant change in satellite orbital distributions across resolution. However, we do see substantial halo-to-halo variance of $(r_\text{peri}, r_\text{apo})$ distributions across host galaxies, as well as a dependence of satellite orbits on host halo mass - systems disrupt at larger pericentres and apocentres in more massive hosts. Our results suggest that either cosmological simulations (including, but not limited to, Auriga) are disrupting satellites far too readily, or that the Milky Way's satellites are more disrupted than current imaging surveys have revealed. Future observing facilities and careful mock observations of these systems will be key to revealing the nature of this apparent discrepancy., Comment: 16 pages, 10 figures, 1 table, submitted to MNRAS

The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterising remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way's accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution $\Lambda$CDM magneto-hydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum., Comment: Accepted for publication in MNRAS, 31 pages, 27 figues, 7 tables. The mock catalogues are available at https://data.desi.lbl.gov/public/papers/mws/auridesi/v1

In the conventional approach to decomposing a rotation curve into a set of contributions from mass model components, the measurements of the rotation curve at different radii are taken to be independent. It is clear, however, that radial correlations are present in such data, for instance (but not only) because the orbital speed depends on the mass distribution at all (or, minimally, inner) radii. We adopt a very simple parametric form for a covariance matrix and constrain its parameters using Gaussian process regression. Applied to the rotation curve of the Milky Way, this suggests the presence of correlations between neighbouring rotation curve points with amplitudes $<10\,\mathrm{km}\,\mathrm{s}^{-1}$ over length scales of $1.5$-$2.5\,\mathrm{kpc}$ regardless of the assumed dark halo component. We show that accounting for such covariance can result in a $\sim 50$ per cent lower total mass estimate for the Milky Way than when it is neglected, and that the statistical uncertainty associated with the covariance is comparable to or exceeds the total systematic uncertainty budget. Our findings motivate including more detailed treatment of rotation curve covariance in future analyses., Comment: MNRAS Letters, accepted. 8 pages, 5 figures. Interactive html version of all figures are available as ancillary files (see right panel)

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 three orders of magnitude. For Sommerfeld and s-wave models, the corresponding emission does not significantly change. For the Sommerfeld model, the Sgr source can be visible above the MW foreground emission, while for s, 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. We find that, given the background emission estimated from this region, the templates from simulations likely have spatial morphology that is too extended to explain the point-like emission that is observed., Comment: 18 pages, 6 figures, 2 tables; added figures 4 and 6; JCAP accepted version

We present the new ARTEMIS Emulator suite of high resolution (baryon mass of $2.23 \times 10^{4}$ $h^{-1}$M$_{\odot}$) zoom-in simulations of Milky Way mass systems. Here, three haloes from the original ARTEMIS sample have been rerun multiple times, systematically varying parameters for the stellar feedback model, the density threshold for star formation, the reionisation redshift and the assumed warm dark matter (WDM) particle mass (assuming a thermal relic). From these simulations emulators are trained for a wide range of statistics that allow for fast predictions at combinations of parameters not originally sampled, running in $\sim 1$ms (a factor of $\sim 10^{11}$ faster than the simulations). In this paper we explore the dependence of the central haloes' stellar mass on the varied parameters, finding the stellar feedback parameters to be the most important. When constraining the parameters to match the present-day stellar mass halo mass relation inferred from abundance matching we find that there is a strong degeneracy in the stellar feedback parameters, corresponding to a freedom in formation time of the stellar component for a fixed halo assembly history. We additionally explore the dependence of the satellite stellar mass function, where it is found that variations in stellar feedback, the reionisation redshift and the WDM mass all have a significant effect. The presented emulators are a powerful tool which allows for fundamentally new ways of analysing and interpreting cosmological hydrodynamic simulations. Crucially, allowing their free (subgrid) parameters to be varied and marginalised, leading to more robust constraints and predictions., Comment: 18 pages, 9 figures, accepted to MNRAS

As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($\Delta t\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$\sigma$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own., Comment: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

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

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science., Comment: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

We use Gaia EDR3 data to identify stars associated with six classical dwarf spheroidals (Draco, Ursa Minor, Sextans, Sculptor, Fornax, Carina) at their outermost radii, beyond their nominal King stellar limiting radius. For all of the dSphs examined, we find radial velocity matches with stars residing beyond the King limiting radius and with $> 50\%$ astrometric probability (four in Draco, two in Ursa Minor, eight in Sextans, two in Sculptor, twelve in Fornax, and five in Carina), indicating that these stars are associated with their respective dwarf spheroidals (dSphs) at high probability. We compare the positions of our candidate "extra-tidal" stars with the orbital tracks of the galaxies, and identify stars, both with and without radial velocity matches, that are consistent with lying along the orbital track of the satellites. However, given the small number of candidate stars, we cannot make any conclusive statements about the significance of these spatially correlated stars. Cross matching with publicly available catalogs of RR Lyrae, we find one RR Lyrae candidate with $> 50\%$ astrometric probability outside the limiting radius in each of Sculptor and Fornax, two such candidates in Draco, nine in Ursa Minor, seven in Sextans, and zero in Carina. Follow-up spectra on all of our candidates, including possible metallicity information, will help confirm association with their respective dSphs, and could represent evidence for extended stellar halos or tidal debris around these classical dSphs., Comment: 19 pages, 9 figures, 6 tables, accepted for publication in MNRAS

There is increasing evidence that a substantial fraction of Milky Way satellite galaxies align in a rotationally-supported plane of satellites, a rare configuration in cosmological simulations of galaxy formation. It has been suggested that other Milky Way substructures (namely young halo globular clusters and stellar/gaseous streams) similarly tend to align with this plane, accordingly dubbed the Vast Polar Structure (VPOS). Using systemic proper motions inferred from Gaia data, we find that globular cluster orbital poles are not clustered in the VPOS direction, though the population with the highest VPOS membership fraction is the young halo clusters (~30%). We additionally provide a current census of stellar streams, including new streams discovered using the Dark Energy Survey and Gaia datasets, and find that stellar stream normals are also not clustered in the direction of the VPOS normal. We also find that, based on orbit modeling, there is a likely association between NGC 3201 and the Gj\"{o}ll stellar stream and that, based on its orbital pole, NGC 4147 is likely not a Sagittarius globular cluster. That the Milky Way's accreted globular clusters and streams do not align in the same planar configuration as its satellites suggests that the plane of satellites is either a particularly stable orbital configuration or a population of recently accreted satellites. Neither of these explanations is particularly likely in light of other recent studies, leaving the plane of satellites problem as one of the more consequential open problems in galaxy formation and cosmology., Comment: 10+8 pages, 4 figures, 4+5 tables. Accepted for publication in MNRAS

We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $\beta$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $\beta\sim -2$ at $r\sim20$ kpc to $\beta\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $\beta(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $\beta(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system., Comment: 16 pages, 8 figures. Key results are summarized in Figures 6 and 7. Accepted to MNRAS

We use the Fermi-LAT gamma-ray observatory to search for gamma-ray emission from four nearby, debris disk-hosting main sequence stars: $\tau$ Ceti, $\epsilon$ Eridani, Fomalhaut, and Vega. For three stars ($\tau$ Ceti, Fomalhaut, and Vega), we establish upper limits that are consistent with theoretical expectations. For $\epsilon$ Eridani, we find a possible spatially coincident source with a soft energy spectrum of $dN/dE \sim E^{-3.6}$. However, at this stage we are unable to rule out that this emission is due to a more extended feature in the diffuse background. In the interpretation that the emission is due to $\epsilon$ Eridani, the $> 100$ MeV gamma-ray luminosity is $\sim 10^{27}$ erg/s $\simeq 3\times 10^{-7}$ L$_\odot$, which is $\sim 10^{10}$ times the gamma-ray luminosity from the disk of the quiet Sun. We find $\lesssim 2 \sigma$ evidence of source variability over a $\sim 7$ year timescale. In the interpretation that the gamma-ray emission from $\epsilon$ Eridani itself, we consider two possible models: 1) cosmic-ray collisions with solid bodies in the debris disk which extends out $\sim$60 AU from the host star, and 2) emission from the stellar activity. For the former model, assuming a total disk mass consistent with infrared measurements, we find that the size distribution of bodies is steeper than expected for a collisional cascade state. If confirmed as being associated with $\epsilon$ Eridani, this would be the first indication of gamma-ray emission from the vicinity of a main sequence star other than the Sun., Comment: 13 pages, 6 figures. Key results are summarized in Table 2 and Figure 5. Accepted to ApJ

The recent discovery of gravitational waves (GWs) from stellar-mass binary black holes (BBHs) provided direct evidence of the existence of these systems. BBH lenses would have gravitational microlensing signatures that are distinct from single-lens signals. We apply Bayesian statistics to examine the distinguishability of BBH microlensing events from single-lens events under ideal observing conditions, using the photometric capabilities of the Korean Microlensing Telescope Network. Given one year of observations, a source star at the Galactic Centre, a GW150914-like BBH lens (total mass 65M$_\odot$, mass ratio 0.8) at half that distance, and an impact parameter of 0.4 Einstein radii, we find that binarity is detectable for BBHs with separations down to 0.0250 Einstein radii, which is nearly 3.5 times greater than the maximum separation for which such BBHs would merge within the age of the Universe. Microlensing searches are thus sensitive to more widely separated BBHs than GW searches, perhaps allowing the discovery of BBH populations produced in different channels of binary formation., Comment: 5 pages, 2 figures