Monica Valluri

We present the catalog of RR Lyrae stars (RRLs) observed in the first year of operations of the Dark Energy Spectroscopic Instrument (DESI) survey. This catalog contains 6240 RRLs out to ∼120 kpc from the Galactic center and over 12,000 individual epochs with homogeneously derived stellar atmospheric parameters. We introduce a novel methodology to model the cyclical variation of the spectroscopic properties of RRLs from single-epoch measurements. We employ this method to infer the systemic velocities and mean temperatures of fundamental- and first-overtone-mode RRLs in our sample (without distinguishing between individual spectral lines). For fundamental-mode pulsators, we obtain radial velocity curves with amplitudes of ∼30–80 km s ⁻¹ and effective temperature curves with 300–1000 K variations, whereas for first-overtone pulsators these amplitudes are ∼20 km s ⁻¹ and ∼600 K, respectively. We use our sample to study the metallicity distribution of the halo and its dependence on Galactocentric distance ( R _GC ). Using a radius-dependent mixture model, we split the data into chemodynamically distinct components and find that our inner halo sample ( R _GC ≲ 50 kpc) is predominantly composed of stars with [Fe/H] ∼ −1.5 and largely radial orbits (with an anisotropy parameter β ∼ 0.94), which we associate with the Gaia–Sausage–Enceladus merger. Stars in the halo field exhibit a broader and more metal-poor [Fe/H] distribution with more circular orbits ( β ∼ 0.39). The metallicity gradient of the metal-rich and metal-poor components is found to be 0.005 and 0.010 dex kpc ⁻¹ , respectively. Our catalog highlights DESI’s tantalizing potential for studying the Milky Way and the pulsation properties of RRLs in the era of large spectroscopic surveys.

Dark matter halos exhibit slow tumbling known as “figure rotation.” Figure rotation occurs in roughly 80% of isolated ΛCDM halos and is relatively unaffected by baryonic feedback, making it an attractive and relatively unexplored test for ΛCDM. However, galaxy mergers are capable of exerting strong torques on a host halo’s quadrupole, warranting an investigation of whether such interactions are capable of influencing figure rotation. Using isolated N- body simulations and six TNG50 halos containing Milky Way analog galaxies, we investigate figure rotation in the context of massive mergers. For isolated galaxies we find that massive mergers elicit rotation of the host halo in both a rapid transient phase during the satellite infall with pattern speeds between 10° and 60° Gyr ⁻¹ , and postmerger as a steady tumbling with pattern speeds between 2° and 40° Gyr ⁻¹ lasting for at least 6–8 Gyr. Both transient and long-lived tumbling depend on the trajectory of the infalling satellite. We find that figure rotation can be induced by radial mergers with no orbital angular momentum, suggesting that tidal torques from the satellite are a driving mechanism. Meanwhile, quadrupole tumbling in TNG50 halos is modified in both pattern speed (3/6 halos) and rotation axis orientation (6/6 halos) during major mergers. Quadrupole tumbling remains stable after the last major merger in five of our six halos. Our results suggest that, rather than tracing the early tidal shear, figure rotation at present day may instead retain a dynamical memory of a halo’s last major merger.

Stellar streams are the result of a host galaxy’s gravitational potential tidally disrupting satellite dwarf galaxies and globular clusters (GCs), causing them to grow leading and trailing tidal tails. The C-19 stellar stream is an extremely metal-poor stellar population, showing chemical abundance patterns characteristic of a GC. However, its large velocity dispersion is difficult to reconcile with a conventional, purely baryonic, disrupting-GC progenitor. Current techniques for stream characterization are primarily applied to Gaia DR3, relying heavily on proper-motion measurements. Using the Dark Energy Spectroscopic Instrument (DESI), which provides radial velocities and metallicities for over 10 million stars and reaches significantly fainter magnitudes than comparable surveys, we employ a mixture model approach to jointly characterize stream populations in proper motions, radial velocities, and metallicities against a Milky Way halo background. By applying this framework to the C-19 stellar stream, we identify a total of 47 spectroscopically confirmed member stars, of which 41 are newly identified and only 6 were previously reported in the literature. In this work, we measure a velocity dispersion of 7 . 8 − 1.3 + 1.5 km s ⁻¹ and a mean metallicity of [Fe/H] = − 3.3 6 − 0.10 + 0.12 . We further identify a novel “spur” feature within the stream. We conclude that our measurements are in line with previous works identifying C-19 as a “hot,” metal-poor stream. In forthcoming work, we will apply this approach to many more streams in the DESI footprint, enabling population-level comparisons with predictions from simulations.

Using 64,000 halo K giants from the Dark Energy Spectroscopic Instrument second Data Release (DR2), we decompose the Milky Way (MW) stellar halo between 3 and 160 kpc into metal-rich (MR) and metal-poor (MP) components via a Gaussian mixture model. The two populations are nearly equal in number but chemically and kinematically distinct: MR stars occupy highly radial orbits with velocity anisotropy of β ≈ 0.94 and metallicity dispersion σ _[Fe/H] ≈ 0.17 dex, without obvious dependence on distance, and are mainly contributed by Gaia-Sausage/Enceladus (GSE) debris. The MR component dominates the inner 30 kpc and reemerges beyond 50 kpc, implying GSE debris can extend to ∼70–80 kpc. MP stars exhibit a weaker radial bias of β ≈ 0.46, decreasing to −0.5 beyond 80 kpc, and with a larger metallicity dispersion of σ _[Fe/H] ≈ 0.46 dex, showing signatures of multiple minor mergers. Both components exhibit net prograde rotation at ∼10–30 kpc with a stronger azimuthal signal in the MP population. The nonequilibrium motions of the outer halo (>50 kpc) are quantified with a dipole-plus-contraction velocity field. We find that the outer halo is simultaneously contracting ( μ _compr = −19 km s ⁻¹ , distance-independent) and subject to reflex motions ( μ _dipole increases from −19 to −44 km s ⁻¹ with radius), reflecting the perturbation from the Large Magellanic Cloud (LMC). We also confirm a linear dependence of mean polar velocity for the outer stellar halo on μ _dipole , a direct consequence of the LMC and MW interaction. Our results provide a quantitative distance-resolved description of the MW’s last major accretion event and its ongoing response to the first infall of the LMC.

The Milky Way Backup Program (MWBP), a survey currently underway with the Dark Energy Spectroscopic Instrument (DESI) on the Nicholas U. Mayall 4 m Telescope, works at the margins of the DESI Main surveys to obtain spectra of millions of additional stars from the Gaia catalog. Efficiently utilizing times between ∼12° and 18° twilight and poor weather conditions, the MWBP extends the range of stellar sources studied to both brighter magnitudes and lower Galactic latitude and declination than the stars studied in DESI’s Main Milky Way Survey. While the MWBP prioritizes candidate giant stars selected from the Gaia catalog (using color and parallax criteria), it also includes an unbiased sample of bright stars (i.e., 11.2 ≲ G < 16 mag) as well as fainter sources (to G ≲ 19 mag). As of 2025 March 1, the survey had obtained spectra of ∼7 million stars, approximately 1.2 million of which are included in the DESI Data Release 1. The DESI spectra cover the wavelength range from 3600 to 9800 Å at a resolution λ /Δ λ varying from 2000 to 5000. The full survey, when completed, will cover an area of more than 21,000 deg ² and include approximately 10 million Gaia sources, roughly equal to the number of stellar spectra obtained through the DESI Main Survey, while only utilizing ≈9% of all DESI observing time. This paper provides an overview of the MWBP, describing the target selection, observing strategy, and an introduction to the resulting data.

We investigate the role of hierarchical assembly in the angular momentum (AM) evolution of galaxies using a sample of 471 Milky Way-mass galaxies from the TNG-50 simulation. While galaxy orientation is often attributed to tidal torques and the cooling of gas within halos, we demonstrate that galaxy reorientation (tilting) is a common consequence of satellite accretion. Specifically, 80 ± 2 % of galaxies show alignment between their present-day AM and the orbital AM of their most massive (dominant) merger progenitor. This reorientation typically results in changes of around 50% in the galaxies’ specific AM, with the most significant shifts occurring in galaxies that were initially highly misaligned. We find only a weak influence from the second most massive merger, and negligible impacts from surviving satellites. We show that accreted stellar halos encode the history of this reorientation. Driven by the same accretion event, the main bodies of galaxies and their stellar halos tend to co-align, with 81 ± 2 % of TNG-50 stellar halos showing prograde rotation relative to the galaxy. This signature will be detectable through major-axis kinematics with 30-meter class telescopes for Milky Way mass galaxies, offering a valuable observational test of this picture. While halo rotation directly constrains the specific AM of mergers within the last ∼ 7 ,Gyr, this kinematic ‘memory’ is largely erased for older and more radial events. Consequently, the Milky Way itself appears to be a notable exception to the general merger-driven trend: TNG-50 analogs with early, radial, and low angular momentum dominant mergers affect present-day disk orientation minimally. The current MW disk orientation may instead reflect the accumulated influences of gas accretion or dark matter torques.

Using the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we examine the six-dimensional space of the anticenter region of the Milky Way stellar disk (150° < Galactic longitude < 220°) using 61,883 main-sequence turnoff stars. We focus on two well-known stellar overdensities in the anticenter: the Monoceros Ring (MRi) and Anticenter Stream (ACS). We find that the MRi overdensity has kinematic signatures consistent with a tidally induced spiral arm, a type of dynamic spiral arm created by an interaction with a satellite galaxy, most likely the Sagittarius dwarf spheroidal galaxy (Sgr). We use the kinematics of the MRi to calculate the two most recent passage times of Sgr, finding 0.25 ± 0.09 Gyr and 1.10 ± 0.23 Gyr from the present day. We validate that the ACS is kinematically decoupled from the MRi because they are moving in opposite radial and vertical directions. We find that the kinematics associated with the ACS extends beyond our defined overdensity. The features we see in the ACS region are likely part of a broader distribution of stars with the same kinematic signature as detected in other places, like the vertical wave in the outer disk and phase spiral.

In 2021 May the Dark Energy Spectroscopic Instrument (DESI) collaboration began a 5 yr spectroscopic redshift survey to produce a detailed map of the evolving three-dimensional structure of the Universe between z = 0 and z ≈ 4. DESI’s principal scientific objectives are to place precise constraints on the equation of state of dark energy, the gravitationally driven growth of large-scale structure, and the sum of the neutrino masses, and to explore the observational signatures of primordial inflation. We present DESI DR1, which consists of all data acquired during the first 13 months of the DESI main survey, as well as a uniform reprocessing of the DESI Survey Validation data, which were previously made public in the DESI Early Data Release. The DR1 main survey includes high-confidence redshifts for 18.7M objects, of which 13.1M are spectroscopically classified as galaxies, 1.6M as quasars, and 4M as stars, making DR1 the largest sample of extragalactic redshifts ever assembled. We summarize the DR1 observations, the spectroscopic data-reduction pipeline and data products, large-scale structure catalogs, value-added catalogs, and describe how to access and interact with the data. In addition to fulfilling its core cosmological objectives with unprecedented precision, we expect DR1 to enable a wide range of transformational astrophysical studies and discoveries.

A recent study of the compact elliptical galaxy NGC 4486B using JWST Near-Infrared Spectrograph integral-field-unit-mode kinematics confirmed a supermassive black hole (SMBH) of mass M BH = 3.6 ± 0.7 × 1 0 8 (∼8% of the stellar mass). In addition to its double nucleus, the nuclear kinematics show pronounced asymmetries: a velocity dispersion peak displaced by 6 pc from the galaxy center and a ∼16 km s ⁻¹ offset in the mean stellar line-of-sight velocity near the SMBH. We examine the origin of the 12 pc double nucleus and show that the observations favor an SMBH surrounded by an eccentric nuclear disk (END). One proposed formation pathway for ENDs involves a gravitational-wave (GW) recoil of the SMBH following a binary merger. Our orbit superposition models contain ∼50% retrograde stars at the edge of the nuclear region, in agreement with END formation simulations. We infer a premerger mass ratio q > 0.15 and a recoil kick of ∼340 km s ⁻¹ . Our N -body simulations show that with such a kick, the SMBH returns to the center within ∼30 Myr. Its flat central core is also consistent with earlier “binary black hole scouring.” We test two alternative mechanisms—buoyancy-driven oscillations and a premerger SMBH binary—but neither reproduces the observed offsets, favoring the GW kick scenario. Our N -body simulations also show that a prograde SMBH binary in a rotating host can become trapped in a corotation resonance, delaying coalescence. Although NGC 4486B is an old galaxy near the Virgo Cluster center, its SMBH appears to have merged only recently, making its nucleus a rare nearby laboratory for postmerger SMBH dynamics.

Using K giants from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we measure the shape, orientation, radial profile, and density anisotropies of the Milky Way (MW) stellar halo over 8 kpc < r _GC < 200 kpc. We identify a triaxial stellar halo (axis ratio 10:8:7), 43° tilted from the disk, showing two break radii at ∼16 and ∼76 kpc, likely associated with Gaia-Sausage/Enceladus and the Large Magellanic Cloud (LMC), respectively. The inner stellar halo (oblate and aligned with the disk, whereas the outer stellar halo becomes prolate and perpendicular to the disk, consistent with the vast polar structure of MW satellites. The twisted halo may arise from the disk−halo angular momentum shift triggered by the infall of a massive satellite. The anisotropic density distribution of the stellar halo is also measured, with successful reidentification of the Hercules-Aquila Cloud North/South (HAC-N/S) overdensity and the Virgo overdensity (VOD). Break radii are found at 15 and 30 kpc for VOD and HAC-N/S, respectively. We identify the LMC transient density wake with a break radius at 60 kpc in the Pisces overdensity region. We also find new observational evidence of the LMC collective density wake, by showing a break radius at ∼100 kpc in the northern Galactic cap with a clear density peak at 90 kpc. In the end, we found that more metal-poor halo stars are more radially extended. Our results provide important clues to the assembly and evolution of the MW stellar halo under the standard cosmic structure formation framework.