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

We present spectroscopy of the ultra-faint Milky Way satellites Eridanus III (Eri III) and DELVE 1. We identify eight member stars in each satellite and place non-constraining upper limits on their velocity and metallicity dispersions. The brightest star in each object is very metal-poor, at [Fe/H] = -3.1 for Eri III and [Fe/H] = -2.8 for DELVE 1. Both of these stars exhibit large overabundances of carbon and very low abundances of the neutron-capture elements Ba and Sr, and we classify them as CEMP-no stars. Because their metallicities are well below those of the Milky Way globular cluster population, and because no CEMP-no stars have been identified in globular clusters, these chemical abundances could suggest that Eri III and DELVE 1 are dwarf galaxies. On the other hand, the two systems have half-light radii of 8 pc and 6 pc, respectively, which is more compact than any known ultra-faint dwarfs. We conclude that Eri III and DELVE 1 are either the smallest dwarf galaxies yet discovered, or they are representatives of a new class of star clusters that underwent chemical evolution distinct from that of ordinary globular clusters. In the latter scenario, such objects are likely the most primordial star clusters surviving today. These possibilities can be distinguished by future measurements of carbon and/or iron abundances for larger samples of stars or improved stellar kinematics for the two systems., Comment: 21 pages, 7 figures, 6 tables. The full versions of Tables 3 and 6 are included in the source files of this submission. Accepted for publication in ApJ

We present new Galactic reddening maps of the high Galactic latitude sky using DESI imaging and spectroscopy. We directly measure the reddening of 2.6 million stars by comparing the observed stellar colors in $g-r$ and $r-z$ from DESI imaging with the synthetic colors derived from DESI spectra from the first two years of the survey. The reddening in the two colors is on average consistent with the \cite{fitzpatrick_correcting_1999} extinction curve with $R_\mathrm{V}=3.1$. We find that our reddening maps differ significantly from the commonly used \cite{schlegel_maps_1998} (SFD) reddening map (by up to 80 mmag in $E(B-V)$), and we attribute most of this difference to systematic errors in the SFD map. To validate the reddening map, we select a galaxy sample with extinction correction based on our reddening map, and this yields significantly better uniformity than the SFD extinction correction. Finally, we discuss the potential systematic errors in the DESI reddening measurements, including the photometric calibration errors that are the limiting factor on our accuracy. The $E(g-r)$ and $E(g-r)$ maps presented in this work, and for convenience their corresponding $E(B-V)$ maps with SFD calibration, are publicly available., Comment: Submitted to the Open Journal of Astrophysics. Associated data files: https://data.desi.lbl.gov/public/papers/mws/desi_dust/y2/v1/maps/

We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins., Comment: Submitted to ApJ, 23 pages, 13 figures 4 tables

We present the stellar value-added catalogue based on the Dark Energy Spectroscopic Instrument (DESI) Early Data Release. The catalogue contains radial velocity and stellar parameter measurements for $\simeq$ 400,000 unique stars observed during commissioning and survey validation by DESI. These observations were made under conditions similar to the Milky Way Survey (MWS) currently carried out by DESI but also include multiple specially targeted fields, such as those containing well-studied dwarf galaxies and stellar streams. The majority of observed stars have $16

Hypervelocity stars (HVSs) are stars which have been ejected from the Galactic Centre (GC) at velocities of up to a few thousand km/s. They are tracers of the Galactic potential and can be used to infer properties of the GC, such as the initial-mass function and assembly history. HVSs are rare, however, with only about a dozen promising candidates discovered so far. In this work we make use of a novel, highly efficient method to identify new HVS candidates in Gaia. This method uses the nearly radial trajectories of HVSs to infer their distances and velocities based on their position and Gaia proper motion alone. Through comparison of inferred distances with Gaia parallaxes and photometry we identified 600 HVS candidates with G<20 including the previously discovered S5-HVS1, out of which we obtained ground-based follow-up observations for 196 stars. As we found no new HVSs based on their radial velocity, we used detailed HVS ejection simulations to significantly improve previous HVS ejection rate constraints. In particular, the ejection rate of HVSs more massive than 1 $\mathrm{M_\odot}$ cannot be higher than $10^{-5}$ yr$^{-1}$ at $2\sigma$ significance. Additionally, we predict that there are 5-45 unbound HVSs in the complete Gaia catalogue ($1\sigma$ interval), most of which will be main-sequence stars of a few M$_\odot$ at heliocentric distances of tens to hundreds of kpc. By comparing our results to literature HVS candidates, we find an indication of either a time-dependent ejection rate of HVSs or a non-GC origin of previously identified HVS candidates., Comment: 17 pages, 16 figures, accepted for publication in 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

We use Legacy Survey photometric data to probe the stellar halo in multiple directions of the sky using a probabilistic methodology to identify Blue Horizontal Branch (BHB) stars. The measured average radial density profile follows a double power law in the range $ 5 < r_{gc}/{\rm kpc} < 120$, with a density break at $r_{gc}\approx20$ kpc. This description, however, falls short, depending on the chosen line-of-sight, with some regions showing no signature of a break in the profile and a wide range of density slopes, e.g. outer slope $-5.5 \lesssim \alpha_{out} \lesssim -4$, pointing towards a highly anisotropic stellar halo. This explains in part the wide range of density profiles reported in the literature owing to different tracers and sky coverage. Using our detailed 3-D stellar halo density map, we quantify the shape of the Pisces overdensity associated with the transient wake response of the Galaxy's (dark) halo to the Large Magellanic Cloud (LMC). Measured in the LMC's coordinate system, Pisces stands above the background, is 60 degrees long and 25 degrees wide and aligned with the LMC's orbit. This would correspond to a wake width of $\sim 32$ kpc at $\sim 70$ kpc. We do not find a statistically significant signature of the collective response in density as previously reported in the literature measured with K giant stars, despite our larger numbers. We release the catalogue constructed in this study with 95,446 possible BHB stars and their BHB probability., Comment: Accepted for publication in A&A. The catalogue can be found here: https://zenodo.org/records/12155488

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!

The Early Data Release (EDR) of the Dark Energy Spectroscopic Instrument (DESI) comprises spectroscopy obtained from 2020 December 14 to 2021 June 10. White dwarfs were targeted by DESI both as calibration sources and as science targets and were selected based on Gaia photometry and astrometry. Here we present the DESI EDR white dwarf catalogue, which includes 2706 spectroscopically confirmed white dwarfs of which approximately 1630 (roughly 60 per cent) have been spectroscopically observed for the first time, as well as 66 white dwarf binary systems. We provide spectral classifications for all white dwarfs, and discuss their distribution within the Gaia Hertzsprung-Russell diagram. We provide atmospheric parameters derived from spectroscopic and photometric fits for white dwarfs with pure hydrogen or helium photospheres, a mixture of those two, and white dwarfs displaying carbon features in their spectra. We also discuss the less abundant systems in the sample, such as those with magnetic fields, and cataclysmic variables. The DESI EDR white dwarf sample is significantly less biased than the sample observed by the Sloan Digital Sky Survey, which is skewed to bluer and therefore hotter white dwarfs, making DESI more complete and suitable for performing statistical studies of white dwarfs., Comment: 19 pages, 17 figures, not including appendix. Submitting to MNRAS. Comments and suggestions welcome

Blue Horizontal Branch (BHB) stars, excellent distant tracers for probing the Milky Way's halo density profile, are distinguished in the $(g-r)_0$ vs $(i-z)_0$ color space from another class of stars, blue straggler stars (BSs). We develop a Bayesian mixture model to classify BHB stars using high-precision photometry data from the Dark Energy Survey Data Release 2 (DES DR2). We select $\sim2100$ highly-probable BHBs based on their $griz$ photometry and the associated uncertainties, and use these stars to map the stellar halo over the Galactocentric radial range $20 \lesssim R \lesssim 70$ kpc. After excluding known stellar overdensities, we find that the number density $n_\star$ of BHBs can be represented by a power law density profile $n_\star \propto R^{-\alpha}$ with an index of $\alpha=4.28_{-0.12}^{+0.13}$, consistent with existing literature values. In addition, we examine the impact of systematic errors and the spatial inhomogeneity on the fitted density profile. Our work demonstrates the effectiveness of high-precision $griz$ photometry in selecting BHB stars. The upcoming photometric survey from the Rubin Observatory, expected to reach depths 2-3 magnitudes greater than DES during its 10-year mission, will enable us to investigate the density profile of the Milky Way's halo out to the virial radius, unravelling the complex processes of formation and evolution in our Galaxy., Comment: 22 pages, 18 figures. Submitted AAS Journal. Comments Welcome Code used in this work can be found at: https://github.com/sazabi4/Yu2024_BHB/

Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5\%) observed in our 300S stars is less than the 50\% of stars found enriched in Milky Way GCs of comparable present-day mass ($\sim10^{4.5}$\msun). We calculate the mass of 300S's progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S's progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of \textit{s}-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted., Comment: 15 pages, 7 tables, 8 figures. Accepted for publication in MNRAS

We present spectroscopic data for 16369 stellar targets within and/or toward 38 dwarf spheroidal galaxies and faint star clusters within the Milky Way halo environment. All spectra come from observations with the multi-object, fiber-fed echelle spectrographs M2FS at the Magellan/Clay telescope or Hectochelle at the MMT, reaching a typical limiting magnitude G < 21. Data products include processed spectra from all observations and catalogs listing estimates -- derived from template model fitting -- of line-of-sight velocity (median uncertainty 1.1 km/s) effective temperature (234 K), (base10 logarithm of) surface gravity (0.52 dex in cgs units), [Fe/H] (0.38 dex) and [Mg/Fe] (0.24 dex) abundance ratios. The sample contains multi-epoch measurements for 3720 sources, with up to 15 epochs per source, enabling studies of intrinsic spectroscopic variability. The sample contains 6078 likely red giant stars (based on surface gravity), and 4494 likely members (based on line-of-sight velocity and Gaia-measured proper motion) of the target systems. The number of member stars per individual target system ranges from a few, for the faintest systems, to ~ 850 for the most luminous. For most systems, our new samples extend over wider fields than have previously been observed; of the likely members in our samples, 823 lie beyond twice the projected halflight radius of their host system, and 42 lie beyond 5 Rhalf., Comment: Published September 2023 in The Astrophysical Journal Supplement Series, processed spectra and catalogs publicly available at the Zenodo database doi:10.5281/zenodo.7837922

The Gaia mission has provided us full astrometric solutions for over 1.5B sources. However, only the brightest 34M of those have radial velocity measurements. As a proof of concept, this paper aims to close that gap, by obtaining radial velocity estimates from the low-resolution BP/RP spectra that Gaia now provides. These spectra are currently published for about 220M sources, with this number increasing to the full $\sim 2$B Gaia sources with Gaia Data Release 4. To obtain the radial velocity measurements, we fit Gaia BP/RP spectra with models based on a grid of synthetic spectra, with which we obtain the posterior probability on the radial velocity for each object. Our measured velocities show systematic biases that depend mainly on colours and magnitudes of stars. We correct for these effects by using external catalogues of radial velocity measurements. We present in this work a catalogue of about 6.4M sources with our most reliable radial velocity measurements and uncertainties $<300$ km s$^{-1}$ obtained from the BP/RP spectra. About 23% of these have no previous radial velocity measurement in Gaia RVS. Furthermore, we provide an extended catalogue containing all 125M sources for which we were able to obtain radial velocity measurements. The latter catalogue, however, also contains a fraction of measurements for which the reported radial velocities and uncertainties are inaccurate. Although typical uncertainties in the catalogue are significantly higher compared to those obtained with precision spectroscopy instruments, the number of potential sources for which this method can be applied is orders of magnitude higher than any previous radial velocity catalogue. Further development of the analysis could therefore prove extremely valuable in our understanding of Galactic dynamics., Comment: 15 pages, 18 figures, accepted in A&A

We present stellar parameters and chemical abundances of 30 elements for five stars located at large radii (3.5-10.7 times the half-light radius) in the Sextans dwarf spheroidal galaxy. We selected these stars using proper motions, radial velocities, and metallicities, and we confirm them as metal-poor members of Sextans with -3.34 <= [Fe/H] <= -2.64 using high-resolution optical spectra collected with the Magellan Inamori Kyocera Echelle spectrograph. Four of the five stars exhibit normal abundances of C (-0.34 <= [C/Fe] <= +0.36), mild enhancement of the alpha elements Mg, Si, Ca, and Ti ([alpha/Fe] = +0.12 +/- 0.03), and unremarkable abundances of Na, Al, K, Sc, V, Cr, Mn, Co, Ni, and Zn. We identify three chemical signatures previously unknown among stars in Sextans. One star exhibits large overabundances ([X/Fe] > +1.2) of C, N, O, Na, Mg, Si, and K, and large deficiencies of heavy elements ([Sr/Fe] = -2.37 +/- 0.25, [Ba/Fe] = -1.45 +/- 0.20, [Eu/Fe] < +0.05), establishing it as a member of the class of carbon-enhanced metal-poor stars with no enhancement of neutron-capture elements. Three stars exhibit moderate enhancements of Eu (+0.17 <= [Eu/Fe] <= +0.70), and the abundance ratios among 12 neutron-capture elements are indicative of r-process nucleosynthesis. Another star is highly enhanced in Sr relative to heavier elements ([Sr/Ba] = +1.21 +/- 0.25). These chemical signatures can all be attributed to massive, low-metallicity stars or their end states. Our results, the first for stars at large radius in Sextans, demonstrate that these stars were formed in chemically inhomogeneous regions, such as those found in ultra-faint dwarf galaxies., Comment: Accepted for publication in the Astrophysical Journal (22 pages, 6 figures, 1 machine-readable table)

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 observations from the Dark Energy Spectroscopic Instrument (DESI) will significantly increase the numbers of known extremely metal-poor stars by a factor of ~ 10, improving the sample statistics to study the early chemical evolution of the Milky Way and the nature of the first stars. In this paper we report high signal-to-noise follow-up observations of 9 metal-poor stars identified during the DESI commissioning with the Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy (OSIRIS) instrument on the 10.4m Gran Telescopio Canarias (GTC). The analysis of the data using a well-vetted methodology confirms the quality of the DESI spectra and the performance of the pipelines developed for the data reduction and analysis of DESI data., Comment: 14 pages, 4 figures, accepted for publication in ApJ, data available from https://doi.org/10.5281/zenodo.8363303

We present Magellan/M2FS spectroscopy of four recently discovered Milky Way star clusters (Gran 3/Patchick~125, Gran 4, Garro 01, LP 866) and two newly discovered open clusters (Gaia 9, Gaia 10) at low Galactic latitudes. We measure line-of-sight velocities and stellar parameters ([Fe/H], $\log{g}$, $T_{\rm eff}$, [Mg/Fe]) from high resolution spectroscopy centered on the Mg triplet and identify 20-80 members per star cluster. We determine the kinematics and chemical properties of each cluster and measure the systemic proper motion and orbital properties by utilizing Gaia astrometry. We find Gran 3 to be an old, metal-poor (mean metallicity of [Fe/H]=-1.84) globular cluster located in the Galactic bulge on a retrograde orbit. Gran 4 is an old, metal-poor ([Fe/H]=-1.84) globular cluster with a halo-like orbit that happens to be passing through the Galactic plane. The orbital properties of Gran 4 are consistent with the proposed LMS-1/Wukong and/or Helmi streams merger events. Garro 01 is metal-rich ([Fe/H]=-0.30) and on a near circular orbit in the outer disk but its classification as an open cluster or globular cluster is ambiguous. . Gaia 9 and Gaia 10 are among the most distant known open clusters at $R_{GC}\sim 18, 21.2~kpc$ and most metal-poor with [Fe/H]~-0.50,-0.46 for Gaia 9 and Gaia 10, respectively. LP 866 is a nearby, metal-rich open cluster ([Fe/H]$=+0.1$). The discovery and confirmation of multiple star clusters in the Galactic plane shows the power of {\it Gaia} astrometry and the star cluster census remains incomplete., Comment: 20 pages, 15 figures, accepted to MNRAS, associated data products available at https://doi.org/10.5281/zenodo.7809128

Very metal-poor stars ($\rm[Fe/H] < -2$) in the Milky Way are fossil records of early chemical evolution and the assembly and structure of the Galaxy. However, they are rare and hard to find. Gaia DR3 has provided over 200 million low-resolution ($R \approx 50$) XP spectra, which provides an opportunity to greatly increase the number of candidate metal-poor stars. In this work, we utilise the \texttt{XGBoost} classification algorithm to identify $\sim$200,000 very metal-poor star candidates. Compared to past work, we increase the candidate metal-poor sample by about an order of magnitude, with comparable or better purity than past studies. Firstly, we develop three classifiers for bright stars ($BP$ $<$ 16). They are Classifier-T (for Turn-off stars), Classifier-GC (for Giant stars with high completeness), and Classifier-GP (for Giant stars with high purity) with expected purity of 52\%/45\%/76\% and completeness of 32\%/93\%/66\% respectively. These three classifiers obtained a total of 11,000/111,000/44,000 bright metal-poor candidates. We apply model-T and model-GP on faint stars ($BP$ $>$ 16) and obtain 38,000/41,000 additional metal-poor candidates with purity 29\%/52\%, respectively. We make our metal-poor star catalogs publicly available, for further exploration of the metal-poor Milky Way., Comment: Zenodo Link To Our Tables: https://doi.org/10.5281/zenodo.8360958

Accreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the `stellar haloes' of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass-metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L > L_host/100. We find that the MW stellar halo has N~1-3 massive progenitors (L > 10^8 L_Sun) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of two. Future datasets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass-scale for both the MW and its satellites., Comment: 13 pages, 11 figures. MNRAS in press

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator., Comment: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

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 describe the spectroscopic data processing pipeline of the Dark Energy Spectroscopic Instrument (DESI), which is conducting a redshift survey of about 40 million galaxies and quasars using a purpose-built instrument on the 4-m Mayall Telescope at Kitt Peak National Observatory. The main goal of DESI is to measure with unprecedented precision the expansion history of the Universe with the Baryon Acoustic Oscillation technique and the growth rate of structure with Redshift Space Distortions. Ten spectrographs with three cameras each disperse the light from 5000 fibers onto 30 CCDs, covering the near UV to near infrared (3600 to 9800 Angstrom) with a spectral resolution ranging from 2000 to 5000. The DESI data pipeline generates wavelength- and flux-calibrated spectra of all the targets, along with spectroscopic classifications and redshift measurements. Fully processed data from each night are typically available to the DESI collaboration the following morning. We give details about the pipeline's algorithms, and provide performance results on the stability of the optics, the quality of the sky background subtraction, and the precision and accuracy of the instrumental calibration. This pipeline has been used to process the DESI Survey Validation data set, and has exceeded the project's requirements for redshift performance, with high efficiency and a purity greater than 99 percent for all target classes., Comment: AJ, revised version, 55 pages, 55 figures, 4 tables

Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-structures are predicted to emerge following the dissolution of a satellite, when its tidal debris is given time to wind up, thin out and fold. Additionally, the observed energy and angular momentum (E, L_z) distribution appears more prograde at high energies, possibly revealing the original orbital angular momentum of the in-falling galaxy. The energy distribution of the debris is strongly asymmetric with a peak at low E -- which, we surmise, may be evidence of the dwarf's rapid sinking -- and riddled with wrinkles and bumps. If these small-scale energy inhomogeneities have been seeded during or immediately after the interaction with the Milky Way, and are not due to the spatial restriction of our study, then making use of the (v_r,r) chevrons to constrain the time of the merger becomes cumbersome. Nonetheless, we demonstrate that similar phase-space and (E,L_z) sub-structures are present in numerical simulations of galaxy interactions, both in bespoke N-body runs and in cosmological hydrodynamical zoom-in suites. The remnant traces of the progenitor's disruption and the signatures of the on-going phase-mixing discovered here will not only help to constrain the properties of our Galaxy's most important interaction, but also can be used as a novel tool to map out the Milky Way's current gravitational potential and its perturbations., Comment: 16 pages, 14 figures, submitted to MNRAS. Comments are welcome!

In 2021 May, the Dark Energy Spectroscopic Instrument (DESI) began a 5 yr survey of approximately 50 million total extragalactic and Galactic targets. The primary DESI dark-time targets are emission line galaxies (ELGs), luminous red galaxies (LRGs) and quasars (QSOs). In bright time, DESI will focus on two surveys known as the Bright Galaxy Survey (BGS) and the Milky Way Survey (MWS). DESI also observes a selection of "secondary" targets for bespoke science goals. This paper gives an overview of the publicly available pipeline (desitarget) used to process targets for DESI observations. Highlights include details of the different DESI survey targeting phases, the targeting ID (TARGETID) used to define unique targets, the bitmasks used to indicate a particular type of target, the data model and structure of DESI targeting files, and examples of how to access and use the desitarget code base. This paper will also describe "supporting" DESI target classes, such as standard stars, sky locations, and random catalogs that mimic the angular selection function of DESI targets. The DESI target selection pipeline is complex and sizable; this paper attempts to summarize the most salient information required to understand and work with DESI targeting data., Comment: AJ, accepted, 27 pages, 4 figures, 10 tables, one of a suite of 8 papers detailing targeting for DESI. Minor textual updates to better match the final, accepted version. Also added two missing co-authors

We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes |b|>20 degrees, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also include several high-completeness samples of rare stellar types, including white dwarfs, low-mass stars within 100pc of the Sun, and horizontal branch stars. We summarize the potential of DESI to advance understanding of Galactic structure and stellar evolution. We introduce the final definitions of the main MWS target classes and estimate the number of stars in each class that will be observed. We describe our pipelines for deriving radial velocities, atmospheric parameters, and chemical abundances. We use ~500,000 spectra of unique stellar targets from the DESI Survey Validation program (SV) to demonstrate that our pipelines can measure radial velocities to ~1 km/s and [Fe/H] accurate to ~0.2 dex for typical stars in our main sample. We find the stellar parameter distributions from ~100 sq. deg of SV observations with >90% completeness on our main sample are in good agreement with expectations from mock catalogs and previous surveys., Comment: Accepted for publication in ApJ, 44 pages, 25 figures, 4 tables, one of a suite of 8 papers detailing targeting for DESI; v2 added links to data shown in figures, added citations to other DESI papers, corrected author list and minor typos; v3 fixed minor errors in Fig. 6 and clarified associated text; v4 updated to include minor changes in response to review

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

Context. The presence of [$\alpha$/Fe]-[Fe/H] bi-modality in the Milky Way disc has animated the Galactic archaeology community since more than two decades. Aims. Our goal is to investigate the chemical, temporal, and kinematical structure of the Galactic discs using abundances, kinematics, and ages derived self-consistently with the new Bayesian framework SAPP. Methods. We employ the public Gaia-ESO spectra, as well as Gaia EDR3 astrometry and photometry. Stellar parameters and chemical abundances are determined for 13 426 stars using NLTE models of synthetic spectra. Ages are derived for a sub-sample of 2 898 stars, including subgiants and main-sequence stars. The sample probes a large range of Galactocentric radii, $\sim$ 3 to 12 kpc, and extends out of the disc plane to $\pm$ 2 kpc. Results. Our new data confirm the known bi-modality in the [Fe/H] - [$\alpha$/Fe] space, which is often viewed as the manifestation of the chemical thin and thick discs. The over-densities significantly overlap in metallicity, age, and kinematics, and none of these is a sufficient criterion for distinguishing between the two disc populations. Different from previous studies, we find that the $\alpha$-poor disc population has a very extended [Fe/H] distribution and contains $\sim$ 20$\%$ old stars with ages of up to $\sim$ 11 Gyr. Conclusions. Our results suggest that the Galactic thin disc was in place early, at look-back times corresponding to redshifts z $\sim$ 2 or more. At ages $\sim$ 9 to 11 Gyr, the two disc structures shared a period of co-evolution. Our data can be understood within the clumpy disc formation scenario that does not require a pre-existing thick disc to initiate a formation of the thin disc. We anticipate that a similar evolution can be realised in cosmological simulations of galaxy formation., Comment: 14 pages, 10 figures, re-submitted to A&A

We present the first global analysis of the impact of the April 2022 cuts to the future pensions of members of the Universities Superannuation Scheme. For the 196,000 active members, if Consumer Price Inflation (CPI) remains at its historic average of 2.5%, the distribution of the range of cuts peaks between 30%-35%. This peak increases to 40%-45% cuts if CPI averages 3.0%. The global loss across current USS scheme members, in today's money, is calculated to be 16-18 billion GBP, with most of the 71,000 staff under the age of 40 losing between 100k-200k GBP each, for CPI averaging 2.5%-3.0%. A repeated claim made during the formal consultation by the body representing university management (Universities UK) that those earning under 40k GBP would receive a "headline" cut of 12% to their future pension is shown to be a serious underestimate for realistic CPI projections., Comment: 29 pp, 7 figures

We present the results of fitting a flexible stellar stream density model to a collection of thirteen streams around the Milky Way, using photometric data from DES, DECaLS, and Pan-STARRS. We construct density maps for each stream and characterise their tracks on the sky, width, and distance modulus curves along the length of each stream. We use these measurements to compute lengths and total luminosities of streams and identify substructures. Several streams show prominent substructures, such as stream broadening, gaps, large deviations of stream tracks and sharp changes in stream densities. Examining the group of streams as a population, as expected we find that streams with globular cluster progenitors are typically narrower than those with dwarf galaxy progenitors, with streams around 100 pc wide showing overlap between the two populations. We also note the average luminosity of globular cluster streams is significantly lower than the typical luminosity of intact globular clusters. The likely explanation is that observed globular cluster streams preferentially come from lower luminosity and lower density clusters. The stream measurements done in a uniform manner presented here will be helpful for more detailed stream studies such as identifying candidate stream members for spectroscopic follow up and stellar stream dynamical modeling., Comment: 25 pages, 40 figures, 4 tables. Accepted for publication in MNRAS

Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA's Gaia mission (Gaia EDR3) and radial velocity information from the SDSS-IV APOGEE survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 13 stars) of the OC stream in the $\alpha$ (O, Mg, Ca, Si, Ti, S), odd-Z (Al, K, V), Fe-peak (Fe, Ni, Mn, Co, Cr) and neutron capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a mono-metallic population and have a median metallicity of --1.92~dex with a dispersion of 0.28 dex. Our results also indicate that the $\alpha$-elements are depleted compared to the known Milky Way populations and that its [Mg/Al] abundance ratio is not consistent with second generation stars from globular clusters. The detailed chemical pattern of these stars indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ~10$^6$ M$_\odot$., Comment: 13 Pages; 4 Figures. Submitted to AAS Journals, comments welcome

Gaia Data Release 3 has provided the astronomical community with the largest stellar spectroscopic survey to date ($>$ 220 million sources). The low resolution (R$\sim$50) blue photometer (BP) and red photometer (RP) spectra will allow for the estimation of stellar atmospheric parameters such as effective temperature, surface gravity and metallicity. We create mock Gaia BP/RP spectra and use Fisher information matrices to probe the resolution limit of stellar parameter measurements using BP/RP spectra. The best-case scenario uncertainties that this analysis provides are then used to produce a mock-observed stellar population in order to evaluate the false positive rate (FPR) of identifying extremely metal-poor (EMP) stars. We conclude that the community will be able to confidently identify metal-poor stars at magnitudes brighter than $G = 16$ using BP/RP spectra. At fainter magnitudes true detections will start to be overwhelmed by false positives. When adopting the commonly-used $G < 14$ limit for metal-poor star searches, we find a FPR for the low-metallicity regimes [Fe/H] $<$ -2, -2.5 and -3 of just 14$\%$, 33$\%$ and 56$\%$ respectively, offering the potential for significant improvements on previous targeting campaigns. Additionally, we explore the chemical sensitivity obtainable directly from BP/RP spectra for Carbon and $\alpha$-elements. We find an absolute Carbon abundance uncertainty of $\sigma_{A(C)} < 1$ dex for Carbon-enriched metal-poor (CEMP) stars, indicating the potential to identify a CEMP stellar population for follow-up confirmation with higher resolution spectroscopy. Finally, we find that large uncertainties in $\alpha$-element abundance measurements using BP/RP spectra means that efficiently obtaining these abundances will be challenging., Comment: 12 pages, 9 figures, 1 table. Accepted for publication in MNRAS

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

Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science., Comment: Contributed white paper to Snowmass 2021, CF03; minor revisions

We combine the Dark Energy Camera Legacy Survey (DECaLS) DR8 photometry with Gaia photometry to study the 3-D structure of Bootes I, Draco, Ursa Minor, Sextans and Sculptor dwarf galaxies using blue horizontal branch (BHB) stars as distance indicators. We construct a new colour-absolute magnitude of BHB stars that we use to measure the distance gradients within the body of the dwarf galaxies. We detect a statistically significant non-zero gradient only in Sextans and Sculptor. Through modeling of the gradient and 2-D density of the systems by triaxial Plummer models we find that the distance gradients in both dwarf galaxies are inconsistent with prolate shape, but compatible with oblate or triaxial shapes. In order to explain the observed gradients, oblate models of Sextans and Sculptor need to have a significant intrinsic ellipticity larger than $0.47$ for Sextans and $0.46$ for Sculptor. The flattened oblate shape may imply a significant anisotropy in velocity distribution in order to be consistent with the lack of significant velocity gradients in these systems., Comment: 21 pages, 7 figures

We present a RR Lyrae (RRL) catalogue based on the combination of the third data release of the Zwicky Transient Facility (ZTF DR3) and \textit{Gaia} EDR3. We use a multi-step classification pipeline relying on the Fourier decomposition fitting to the multi-band ZTF light curves and random forest classification. The resulting catalogue contains 71,755 RRLs with period and light curve parameter measurements and has completeness of 0.92 and purity of 0.92 with respect to the SOS \textit{Gaia} DR2 RRLs. The catalogue covers the Northern sky with declination $\geq -28^\circ$, its completeness is $\gtrsim 0.8$ for heliocentric distance $\leq 80$~kpc, and the most distant RRL at 132~kpc. Compared with several other RRL catalogues covering the Northern sky, our catalogue has more RRLs around the Galactic halo and is more complete at low Galactic latitude areas. Analysing the spatial distribution of RRL in the catalogue reveals the previously known major over-densities of the Galactic halo, such as the Virgo over-density and the Hercules-Aquila Cloud, with some evidence of an association between the two. We also analyse the Oosterhoff fraction differences throughout the halo, comparing it with the density distribution, finding increasing Oosterhoff I fraction at the elliptical radii between 16 and 32 kpc and some evidence of different Oosterhoff fractions across various halo substructures.

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 measure rotational broadening in spectra taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to characterise the relationship between stellar multiplicity and rotation. We create a sample of 2786 giants and 24 496 dwarfs with stellar parameters and multiple radial velocities from the APOGEE pipeline, projected rotation speeds \vsini\ determined from our own pipeline, and distances, masses, and ages measured by Sanders \& Das. We use the statistical distribution of the maximum shift in the radial velocities, \drvm, as a proxy for the close binary fraction to explore the interplay between stellar evolution, rotation, and multiplicity. Assuming that the minimum orbital period allowed is the critical period for Roche Lobe overflow and rotational synchronization, we calculate theoretical upper limits on expected \vsini\ and \drvm\ values. These expectations agree with the positive correlation between the maximum \drvm\ and \vsini\ values observed in our sample as a function of \logg. We find that the fast rotators in our sample have a high occurrence of short-period ($\log(P/\text{d})\lesssim 4$) companions. We also find that old, rapidly-rotating main sequence stars have larger completeness-corrected close binary fractions than their younger peers. Furthermore, rapidly-rotating stars with large \drvm\ consistently show differences of 1-10 Gyr between the predicted gyrochronological and measured isochronal ages. These results point towards a link between rapid rotation and close binarity through tidal interactions. We conclude that stellar rotation is strongly correlated with stellar multiplicity in the field, and caution should be taken in the application of gyrochronology relations to cool stars., Comment: 12 pages, 9 figures; accepted by MNRAS

The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low mass globular cluster. In this work we use Magellan/MIKE high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($\sigma(\mathrm{[Fe~II/H]}) = 0.04^{+0.11}_{-0.03}$) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr II/Fe] spanning $\sim$1 dex, while [Ba II/Fe] shows nearly no intrinsic spread ($\sigma(\mathrm{[Ba~II/H]}) = {0.03}^{+0.10}_{-0.02}$). This abundance signature is best interpreted as slow neutron capture element production from a massive fast-rotating metal-poor star ($15-20 \mathrm{M}_\odot$, $v_\mathrm{ini}/v_\mathrm{crit} = 0.4$, $[\mathrm{Fe/H}] = -3.8$). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Neutron capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant ($A(\mathrm{Li}) = 3.1 \pm 0.1$). At $[\mathrm{Fe/H}] = -2.93$ it is among the most metal-poor lithium-rich giants known., Comment: Accepted to ApJ 2021-07-09

We present new MMT/Hectochelle spectroscopic measurements for 257 stars observed along the line of sight to the ultra-faint dwarf galaxy Triangulum II. Combining with results from previous Keck/DEIMOS spectroscopy, we obtain a sample that includes 16 likely members of Triangulum II, with up to 10 independent redshift measurements per star. To this multi-epoch kinematic data set we apply methodology that we develop in order to infer binary orbital parameters from sparsely sampled radial velocity curves with as few as two epochs. For a previously-identified (spatially unresolved) binary system in Tri~II, we infer an orbital solution with period $296.0_{-3.3}^{+3.8} \rm~ days$ , semi-major axis $1.12^{+0.41}_{-0.24}\rm~AU$, and a systemic velocity $ -380.0 \pm 1.7 \rm~km ~s^{-1}$ that we then use in the analysis of Tri~II's internal kinematics. Despite this improvement in the modeling of binary star systems, the current data remain insufficient to resolve the velocity dispersion of Triangulum II. We instead find a 95% confidence upper limit of $\sigma_{v} \lesssim 3.4 \rm ~km~s^{-1}$.

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 use a geometric method to derive (two-dimensional) separation functions amongst pairs of objects within populations of specified position function $\mathrm{d} N/\mathrm{d} \vec{R}$. We present analytic solutions for separation functions corresponding to a uniform surface density within a circular field, a Plummer sphere (viewed in projection), and the mixture thereof -- including contributions from binary objects within both sub-populations. These results enable inferences about binary object populations via direct modeling of object position and pair separation data, without resorting to standard estimators of the two-point correlation function. Analyzing mock data sets designed to mimic known dwarf spheroidal galaxies, we demonstrate the ability to recover input properties including the number of wide binary star systems and, in cases where the number of resolved binary pairs is assumed to be $\gtrsim$ a few hundred, characteristic features (e.g., steepening and/or truncation) of their separation function. Combined with forthcoming observational capabilities, this methodology opens a window onto the formation and/or survival of wide binary populations in dwarf galaxies, and offers a novel probe of dark matter substructure on the smallest galactic scales., Comment: 20 pages, 15 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

We report the serendipitous discovery of a late-type giant star that exhibited a smooth, eclipse-like drop in flux to a depth of 97 per cent. Minimum flux occurred in April 2012 and the total event duration was a few hundred days. Light curves in V, I and K$_s$ from the Optical Gravitational Lensing Experiment and VISTA Variables in the Via Lactea surveys show a remarkably achromatic event. During 17 years of observational coverage of this source only one such event was detected. The physical properties of the giant star itself appear somewhat unusual, which may ultimately provide a clue towards the nature of the system. By modelling the event as an occultation by an object that is elliptical in projection with uniform transparency, we place limits on its physical size and velocity. We find that the occultation is unlikely to be due to a chance alignment with a foreground object. We consider a number of possible candidates for the occulter, which must be optically thick and possess a radius or thickness in excess of 0.25 au. None are completely satisfactory matches to all the data. The duration, depth and relative achromaticity of the dip mark this out as an exceptionally unusual event, whose secret has still not been fully revealed. We find two further candidates in the VVV survey and we suggest that these systems, and two previously known examples, may point to a broad class of long period eclipsing binaries wherein a giant star is occulted by a circumsecondary disc., Comment: 18 pages, 11 figures. Accepted for publication in the MNRAS main journal

We present a search for "hyper-compact" star clusters in the Milky Way using a combination of Gaia and the Dark Energy Camera Legacy Survey (DECaLS). Such putative clusters, with sizes of ~1 pc and containing 500-5000 stars, are expected to remain bound to intermediate-mass black holes (Mbh~10^3-10^5 M-sun) that may be accreted into the Milky Way halo within dwarf satellites. Using the semi-analytic model SatGen we find an expected ~100 wandering intermediate-mass black holes with if every infalling satellite hosts a black hole. We do not find any such clusters in our search. Our upper limits rule out 100% occupancy, but do not put stringent constraints on the occupation fraction. Of course, we need stronger constraints on the properties of the putative star clusters, including their assumed sizes as well as the fraction of stars that would be compact remnants., Comment: 20 pages, 8 figures, under review at ApJ

The Fornax dwarf spheroidal galaxy has an anomalous number of globular clusters, five, for its stellar mass. There is a longstanding debate about a potential sixth globular cluster (Fornax~6) that has recently been `rediscovered' in DECam imaging. We present new Magellan/M2FS spectroscopy of the Fornax~6 cluster and Fornax dSph. Combined with literature data we identify $\sim15-17$ members of the Fornax~6 cluster that this overdensity is indeed a star cluster and associated with the Fornax dSph. The cluster is significantly more metal-rich (mean metallicity of $\overline{\rm [Fe/H]}=-0.71\pm0.05$) than the other five Fornax globular clusters ($-2.5<[Fe/H]<-1.4$) and more metal-rich than the bulk of Fornax. We measure a velocity dispersion of $5.6_{-1.6}^{+2.0}\,{\rm km \, s^{-1}}$ corresponding to anomalously high mass-to-light of 15$<$M/L$<$258 at 90\% confidence when calculated assuming equilibrium. Two stars inflate this dispersion and may be either Fornax field stars or as yet unresolved binary stars. Alternatively the Fornax~6 cluster may be undergoing tidal disruption. Based on its metal-rich nature, the Fornax 6 cluster is likely younger than the other Fornax clusters, with an estimated age of $\sim2$ Gyr when compared to stellar isochrones. The chemodynamics and star formation history of Fornax shows imprints of major events such as infall into the Milky Way, multiple pericenter passages, star formation bursts, and/or potential mergers or interactions. Any of these events may have triggered the formation of the Fornax~6 cluster., Comment: 17 pages, 9 figures. Accepted to ApJ. Data catalogs included. Comments welcome

We use $Gaia$ eDR3 data and legacy spectroscopic surveys to map the Milky Way disc substructure towards the Galactic Anticenter at heliocentric distances $d\geq10\,\rm{kpc}$. We report the discovery of multiple previously undetected new filaments embedded in the outer disc in highly extincted regions. Stars in these over-densities have distance gradients expected for disc material and move on disc-like orbits with $v_{\phi}\sim170-230\,\rm{km\,s^{-1}}$, showing small spreads in energy. Such a morphology argues against a quiescently growing Galactic thin disc. Some of these structures are interpreted as excited outer disc material, kicked up by satellite impacts and currently undergoing phase-mixing ("feathers"). Due to the long timescale in the outer disc regions, these structures can stay coherent in configuration space over several Gyrs. We nevertheless note that some of these structures could also be folds in the perturbed disc seen in projection from the Sun's location. A full 6D phase-space characterization and age dating of these structure should help distinguish between the two possible morphologies., Comment: 5 pages, 4 figures, submitted to MNRAS

Gaia DR2 published positions, parallaxes and proper motions for an unprecedented 1,331,909,727 sources, revolutionising the field of Galactic dynamics. We complement this data with the Astrometry Spread Function (ASF), the expected uncertainty in the measured positions, proper motions and parallax for a non-accelerating point source. The ASF is a Gaussian function for which we construct the 5D astrometric covariance matrix as a function of position on the sky and apparent magnitude using the Gaia DR2 scanning law and demonstrate excellent agreement with the observed data. This can be used to answer the question `What astrometric covariance would Gaia have published if my star was a non-accelerating point source?'. The ASF will enable characterisation of binary systems, exoplanet orbits, astrometric microlensing events and extended sources which add an excess astrometric noise to the expected astrometry uncertainty. By using the ASF to estimate the unit weight error (UWE) of Gaia DR2 sources, we demonstrate that the ASF indeed provides a direct probe of the excess source noise. We use the ASF to estimate the contribution to the selection function of the Gaia astrometric sample from a cut on astrometric_sigma5d_max showing high completeness for $G<20$ dropping to $<1\%$ in underscanned regions of the sky for $G=21$. We have added an ASF module to the Python package SCANNINGLAW (https://github.com/gaiaverse/scanninglaw) through which users can access the ASF., Comment: 17 pages, 18 figures, published in MNRAS

We perform consistent reductions and measurements for three ultra-faint dwarf galaxies (UFDs): Bo\"otes I, Leo IV and Leo V. Using the public archival data from the GIRAFFE spectrograph on the Very Large Telescope (VLT), we locate new members and provide refined measurements of physical parameters for these dwarf galaxies. We identify nine new Leo IV members and four new Leo V members, and perform a comparative analysis of previously discovered members. Additionally, we identify one new binary star in both Leo IV and Leo V. After removing binary stars, we recalculate the velocity dispersions of Bo\"otes I and Leo IV to be 5.1$^{+0.7}_{-0.8}$ and 3.4$^{+1.3}_{-0.9}$ km s$^{-1}$, respectively; We do not resolve the Leo V velocity dispersion. We identify a weak velocity gradient in Leo V that is $\sim$4$\times$ smaller than the previously calculated gradient and that has a corresponding position angle which differs from the literature value by $\sim$120 deg. Combining the VLT data with previous literature, we re-analyze the Bo\"otes I metallicity distribution function and find that a model including infall of pristine gas while Bo\"otes I was forming stars best fits the data. Our analysis of Leo IV, Leo V and other UFDs will enhance our understanding of these enigmatic stellar populations and contribute to future dark matter studies. This is the first in a series of papers examining thirteen UDFs observed with VLT/GIRAFFE between 2009 and 2017. Similar analyses of the remaining ten UFDs will be presented in forthcoming papers., Comment: 29 pages, 18 figures, 8 tables. Accepted for publication on ApJ. Code and data can be found at https://github.com/sazabi4/vlt_dwarf/