Search for: All records

Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, pulsar timing, and eclipsing binary stars in the solar neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density and infer the amount of dark matter (DM) by accounting for the mass in stars, gas, and dust. Leveraging FIRE-2 simulations of Milky Way–mass galaxies we compare vertical acceleration profiles between cold DM (CDM) and self-interacting DM (SIDM) with a constant cross section of 1 cm²g⁻¹across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos manifest a more oblate shape and consistently exhibit higher local stellar and DM densities and steeper vertical acceleration gradients, up to 10%–30% steeper near the solar neighborhood. However, similar magnitude changes can arise from azimuthal variations in the baryonic components at a fixed radius and external influences like mergers, making it difficult to distinguish between CDM and SIDM using acceleration measurements in a single galaxy.

 

We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (f_diss= 0.75). In this paper, we examine the properties of dwarf galaxies withM_*∼ 10⁵–10⁹M_⊙in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm²g⁻¹. On the contrary, models with (σ/m) = 10 cm²g⁻¹produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm²g⁻¹are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm²g⁻¹,f_diss= 0.75 are in tension in massive dwarfs (M_halo∼ 10¹¹M_⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures.

 

The third data release (DR3) of Gaia has provided a fivefold increase in the number of radial velocity measurements of stars, as well as a stark improvement in parallax and proper motion measurements. To help with studies that seek to test models and interpret Gaia DR3, we present nine Gaia synthetic surveys, based on three solar positions in three Milky Way-mass galaxies of theLattesuite of theFire-2 cosmological simulations. These synthetic surveys match the selection function, radial velocity measurements, and photometry of Gaia DR3, adapting the code baseAnanke, previously used to match the Gaia DR2 release by Sanderson et al. The synthetic surveys are publicly available and can be found athttp://ananke.hub.yt/. Similarly to the previous release ofAnanke, these surveys are based on cosmological simulations and thus are able to model nonequilibrium dynamical effects, making them a useful tool in testing and interpreting Gaia DR3.

 

In this paper, we construct the circular velocity curve of the Milky Way out to ∼30 kpc, providing an updated model of the dark matter density profile. We derive precise parallaxes for 120 309 stars with a data-driven model, using APOGEE DR17 spectra combined with GaiaDR3, 2MASS, and WISE photometry. At outer galactic radii up to 30 kpc, we find a significantly faster decline in the circular velocity curve compared to the inner parts. This decline is better fit with a cored Einasto profile with a slope parameter $0.91^{+0.04}_{-0.05}$ than a generalized Navarro–Frenk–White (NFW) profile. The virial mass of the best-fitting dark matter halo profile is only $1.81^{+0.06}_{-0.05}\times 10^{11}$ M⊙, significantly lower than what a generalized NFW profile delivers. We present a study of the potential systematics, affecting mainly large radii. Such a low mass for the Galaxy is driven by the functional forms tested, given that it probes beyond our measurements. It is found to be in tension with mass measurements from globular clusters, dwarf satellites, and streams. Our best-fitting profile also lowers the expected dark matter annihilation signal flux from the galactic centre by more than an order of magnitude, compared to an NFW profile-fit. In future work, we will explore profiles with more flexible functional forms to more fully leverage the circular velocity curve and observationally constrain the properties of the Milky Way’s dark matter halo.

 

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

 

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm hdbscan on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the hdbscan member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

 

Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e. proto-galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: (i) can have a stellar mass range between 1 × 108 < M⋆ < 2 × 1010 [M⊙], a virial mass range between 3 × 1010 < M⋆ < 6 × 1011 [M⊙], and be as young as 8 ≲ Age ≲ 12.8 [Gyr] (1 ≲ z ≲ 6); (ii) are pre-dominantly centrally concentrated, with $\sim 50~{{\ \rm per\ cent}}$ of the stars contained within 5–10 kpc; (iii) on average show weak but systematic net rotation in the plane of the host’s disc at z = 0 (i.e. 0.25 ≲ 〈κ/κdisc〉 ≲ 0.8); (iv) present [α/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host’s old disc; and (v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that $\sim 60~{{\ \rm per\ cent}}$ of the proto-Milky Way galaxies are comprised by 1 dominant system (1/5 ≲M⋆/M⋆, proto-MilkyWay≲ 4/5) and 4–5 lower mass systems (M⋆/M⋆, proto-MilkyWay≲ 1/10); the other $\sim 40~{{\ \rm per\ cent}}$ are comprised by 2 dominant systems and 3–4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems.

 

Nyx is a nearby, prograde, and high-eccentricity stellar stream physically contained in the thick disk, but its origin is unknown. Nyx could be the remnant of a disrupted dwarf galaxy, in which case the associated dark matter substructure could affect terrestrial dark matter direct-detection experiments. Alternatively, Nyx could be a signature of the Milky Way’s disk formation and evolution. To determine the origin of Nyx, we obtained high-resolution spectroscopy of 34 Nyx stars using Keck/HIRES and Magellan/MIKE. A differential chemical abundance analysis shows that most Nyx stars reside in a metal-rich ([Fe/H] > −1) high-αcomponent that is chemically indistinguishable from the thick disk. This rules out the originally suggested scenario that Nyx is the remnant of a single massive dwarf galaxy merger. However, we also identify 5 substantially more metal-poor stars ([Fe/H] ∼ −2.0) whose chemical abundances are similar to those of the metal-weak thick disk. It remains unclear how stars that are chemically identical to the thick disk can be on such prograde, high-eccentricity orbits. We suggest two most likely scenarios: that Nyx is the result of an early minor dwarf galaxy merger, or that it is a record of the early spin-up of the Milky Way disk—although neither perfectly reproduces the chemodynamic observations. The most likely formation scenarios suggest that future spectroscopic surveys should find Nyx-like structures outside of the solar neighborhood.

 

Abstract In the Λ-Cold Dark Matter model of the universe, galaxies form in part through accreting satellite systems. Previous works have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way–like galaxies in the Latte suite of FIRE-2 cosmological simulations. The resolution of these simulations allows a comparison of contributions from satellites above M * ≳ 10 × 7 M ⊙ , enabling the analysis of observable properties for disrupted satellites in a fully self-consistent and cosmological context. Our results show that the time of accretion and the stellar mass of an accreted satellite are fundamental parameters that in partnership dictate the resulting spatial distribution, orbital energy, and [ α /Fe]-[Fe/H] compositions of the stellar debris of such mergers at present day. These parameters also govern the resulting dynamical state of an accreted galaxy at z = 0, leading to the expectation that the inner regions of the stellar halo ( R GC ≲ 30 kpc) should contain fully phase-mixed debris from both lower- and higher-mass satellites. In addition, we find that a significant fraction of the lower-mass satellites accreted at early times deposit debris in the outer halo ( R GC > 50 kpc) that are not fully phased-mixed, indicating that they could be identified in kinematic surveys. Our results suggest that, as future surveys become increasingly able to map the outer halo of our Galaxy, they may reveal the remnants of long-dead dwarf galaxies whose counterparts are too faint to be seen in situ in higher-redshift surveys.