Exploring the fate of primordial discs in Milky Way-sized galaxies with the GigaEris simulation
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The chemical abundances of a galaxy’s metal-poor stellar population can be used to investigate the earliest stages of its formation and chemical evolution. The Magellanic Clouds are the most massive of the Milky Way’s satellite galaxies and are thought to have evolved in isolation until their recent accretion by the Milky Way. Unlike the Milky Way’s less massive satellites, little is known about the Magellanic Clouds’ metal-poor stars. We have used the mid-infrared metal-poor star selection of Schlaufman & Casey and archival data to target nine LMC and four SMC giants for high-resolution Magellan/MIKE spectroscopy. These nine LMC giants with -2.4 ≲ [Fe/H] ≲ -1.5 and four SMC giants with -2.6 ≲ [Fe/H] ≲ -2.0 are the most metal-poor stars in the Magellanic Clouds yet subject to a comprehensive abundance analysis. While we find that at constant metallicity these stars are similar to Milky Way stars in their α, light, and iron-peak elemental abundances, both the LMC and SMC are enhanced relative to the Milky Way in the r-process element europium. These abundance offsets are highly significant, equivalent to 3.9σ for the LMC, 2.7σ for the SMC, and 5.0σ for the complete Magellanic Cloud sample. We propose that the r-process enhancement of the Magellanic Clouds’ metal-poor stellar population is a result of the Magellanic Clouds’ isolated chemical evolution and long history of accretion from the cosmic web combined with r-process nucleosynthesis on a timescale longer than the core-collapse supernova timescale but shorter than or comparable to the thermonuclear (i.e., Type Ia) supernova timescale.
It has been suggested that a class of chemically peculiar metal-poor stars called iron-rich metal-poor (IRMP) stars formed from molecular cores with metal contents dominated by thermonuclear supernova nucleosynthesis. If this interpretation is accurate, then IRMP stars should be more common in environments where thermonuclear supernovae were important contributors to chemical evolution. Conversely, IRMP stars should be less common in environments where thermonuclear supernovae were not important contributors to chemical evolution. At constant [Fe/H] ≲ -1, the Milky Way's satellite classical dwarf spheroidal (dSph) galaxies and the Magellanic Clouds have lower [α/Fe] than the Milky Way field and globular cluster populations. This difference is thought to demonstrate the importance of thermonuclear supernova nucleosynthesis for the chemical evolution of the Milky Way's satellite classical dSph galaxies and the Magellanic Clouds. We use data from the Sloan Digital Sky Survey Apache Point Observatory Galactic Evolution Experiment and Gaia to infer the occurrence of IRMP stars in the Milky Way's satellite classical dSph galaxies η dSph and the Magellanic Clouds η Mag , as well as in the Milky Way field η MWF and globular cluster populations η MWGC . In order of decreasing occurrence, we find η dSph = $0.07$$^{+0.02}_{-0.02}$, η Mag = $0.037$$^{+0.007}_{-0.006}$, η MWF = $0.0013$$^{+0.0006}_{-0.0005}$, and a 1σ upper limit η MWGC < 0.00057. These occurrences support the inference that IRMP stars formed in environments dominated by thermonuclear supernova nucleosynthesis and that the time lag between the formation of the first and second stellar generations in globular clusters was longer than the thermonuclear supernova delay time.
We present interstellar gas-phase abundances of zinc and chromium in the Milky Way for comparison with the abundances of these elements in damped Ly-alpha systems at redshifts near z approximately 2. We make use of both IUE (International Ultraviolet Explorer) and GHRS (Goddard High Resolution Spectrograph) archival information in deriving the Milky Way values and find that the average difference in (Zn/H) and (Cr/Zn) between the Milky Way and damped Ly-alpha systems at large N(H1) is about 0.9 dex. The Zn/H values we find for the Milky Way are approximately -0.2 to -0.7 dex below solar and are due in large part to depletion of Zn onto dust grains. The amount of this gas-phase depletion is variable and depends upon the H1 column density and molecular fraction of the gas. Milky Way sight lines having low molecular fractions have values of (Zn/H) within approximately 0.2 dex of the solar value. Combined with larger values of (Cr/Zn) found for the damped Ly-alpha systems, this result confirms that measures of (Zn/H) should accurately reflect the lower gas-phase abundances of Zn found in these systems.
APOGEE is a high-resolution (R ~ 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all of the main regions of the Milky Way; in particular, by operating in the H band, APOGEE is uniquely able to probe the dust-hidden inner regions of the Milky Way that are best accessed from the Southern Hemisphere. In this paper we present the targeting strategy of APOGEE-2S, with special attention to documenting modifications to the original, previously published plan. The motivation for these changes is explained as well as an assessment of their effectiveness in achieving their intended scientific objective. In anticipation of this being the last paper detailing APOGEE targeting, we present an accounting of all such information complete through the end of the APOGEE-2S project; this includes several main survey programs dedicated to exploration of major stellar populations and regions of the Milky Way, as well as a full list of programs contributing to the APOGEE database through allocations of observing time by the Chilean National Time Allocation Committee and the Carnegie Institution for Science. This work was presented along with a companion article, Beaton et al. (2021), presenting the final target selection strategy adopted for APOGEE-2 in the Northern Hemisphere.
The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. Methods. A method known as “transfer learning” is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the F IRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia /RAVE data set, which improves sensitivity to properties of the real Milky Way. Results. The result is a catalog that identifies ~767 000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.
Recent observations of the Milky Way and galaxies at high redshifts suggest that galaxy discs were already in place soon after the Big Bang. While the gas infall history of the Milky Way in the inner disc has long been assumed to be characterised by a short accretion timescale, this has not been directly constrained using observations. Using data for the inner regions of the Milky Way recently produced by APOGEE and Gaia and of unprecedented quantity and quality, we aim to derive strong constraints on the infall history of the inner (< 6 kpc) Galaxy (with a focus on stars between 4 and 6 kpc, which we show is an appropriate proxy for the entire inner disc). We implemented gas infall into a chemical evolution model of the Galaxy disc, and used a Schmidt–Kennicutt law to connect the infall to the star formation. We explore a number of models, and two different formulations of the infall law. In one formulation, the infall is non-parametric, and in the other the infall has an explicitly exponential form. We fit the model parameters to the time–[Si/Fe] distribution of solar vicinity stars, and the metallicity and [Si/Fe] distribution function of stars with a galactocentric radius of between 4 and 6 kpc from APOGEE. Our results point to a fast, early gas accretion, and an upper limit on the accretion timescale of around 2 Gyr in the inner disc of the Milky Way. This suggests that at least half the baryons were in place within 2–3 Gyr of the Big Bang, and that half the stars of the inner disc formed within the first 5 Gyr, during the thick disc formation phase. This implies that the stellar mass of the inner disc is dominated by the thick disc, supporting our previous work, and that the gas accretion onto the inner disc was rapid and early.
Context. Satellite galaxies experience multiple physical processes when interacting with their host halos, often leading to the quenching of star formation. In the Local Group, satellite quenching has been shown to be highly efficient, affecting nearly all satellites except the most massive ones. While recent surveys study Milky Way-analogs to assess how representative our Local Group is, the dominant physical mechanisms behind satellite quenching in Milky Way-mass halos remain under debate. Aims. We analyze satellite quenching within the same Milky Way-mass halo simulated using various widely used astrophysical codes, each using different hydrodynamic methods and implementing different supernovae feedback recipes. The goal is to determine whether quenched fractions, quenching timescales, and the dominant quenching mechanisms are consistent across codes or if they show sensitivity to the specific hydrodynamic method and supernovae feedback physics employed. Methods. We used a subset of high-resolution cosmological zoom-in simulations of a Milky Way-mass halo from the multiple-code AGORA CosmoRun suite. Our analysis focuses on comparing satellite quenching across the different models and against observational data. We also analyzed the dominant mechanisms driving satellite quenching in each model. Results. We find that the quenched fraction is consistent with the latest SAGA Survey results within its 1σ host-to-host scatter across all the models. Regarding quenching timescales, all the models reproduce the trend observed in the ELVES survey, Local Group observations, and previous simulations: The less massive the satellite, the shorter its quenching timescale. All of our models converge on the dominant quenching mechanisms: Strangulation halts cold gas accretion in all satellites, while ram pressure stripping is the predominant mechanism for gas removal, and it is particularly effective in satellites with M * <10 8 M ⊙ . Nevertheless, the efficiency of the stripping mechanisms differs among the codes, showing a strong sensitivity to the different supernovae feedback implementations and/or hydrodynamic methods employed.
The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.
ABSTRACT We present the first high-resolution chemical abundances of seven stars in the recently discovered high-energy stream Typhon. Typhon stars have apocentres r ≳ 100 kpc, making this the first detailed chemical picture of the Milky Way’s very distant stellar halo. Though the sample size is limited, we find that Typhon’s chemical abundances are more like a dwarf galaxy than a globular cluster, showing a metallicity dispersion and no presence of multiple stellar populations. Typhon stars display enhanced α-element abundances and increasing r-process abundances with increasing metallicity. The high-α abundances suggest a short star formation duration for Typhon, but this is at odds with expectations for the distant Milky Way halo and the presence of delayed r-process enrichment. If the progenitor of Typhon is indeed a new dwarf galaxy, possible scenarios explaining this apparent contradiction include a dynamical interaction that increases Typhon’s orbital energy, a burst of enhanced late-time star formation that raises [α/Fe], and/or group pre-processing by another dwarf galaxy before infall into the Milky Way. Alternatively, Typhon could be the high-energy tail of a more massive disrupted dwarf galaxy that lost energy through dynamical friction. We cannot clearly identify a known low-energy progenitor of Typhon in the Milky Way, but 70 per cent of high-apocentre stars in cosmological simulations are from high-energy tails of large dwarf galaxies. Typhon’s surprising combination of kinematics and chemistry thus underscores the need to fully characterize the dynamical history and detailed abundances of known substructures before identifying the origin of new substructures.
Within the ΛCDM cosmology, dark matter haloes are composed of both a smooth component and a population of smaller gravitationally bound subhaloes. These components are often treated as a single halo when properties, such as density profiles, are extracted from simulations. Recent work has shown that density profiles change substantially when subhalo mass is excluded. In this paper, we expand on this result by analysing three specific host halo properties – concentration (cNFW), spin (λB), and shape (c/a) – when calculated only from the smooth component of the halo. This analysis is performed on both Milky Way-mass haloes and cluster-mass haloes in high-resolution zoom-in N-body simulations. We find that when subhaloes are excluded, the median value of (1) cNFW is enhanced by $\approx 30\pm 11$ and $\approx 77\pm 8.1~{{\ \rm per\ cent}}$ for Milky Way-mass ($10^{12.1}\, \text{M}_\odot$) and cluster-mass ($10^{14.8}\, \text{M}_\odot$) haloes, respectively, (2) λB is reduced for Milky Way-mass by $\approx 11\pm 4.9~{{\ \rm per\ cent}}$ and cluster-mass haloes by $\approx 27\pm 3.5~{{\ \rm per\ cent}}$. Additionally, with the removal of subhaloes, cluster-mass haloes tend to become more spherical as the ratio of minor-to-major axis, c/a, increases by $\approx 11\pm 3.6~{{\ \rm per\ cent}}$, whereas Milky Way-mass haloes remain approximately the same shape with c/a changed by $\approx 1.0\pm 5.8~{{\ \rm per\ cent}}$. Fractional changes of each of these properties depend primarily on the amount of mass in subhaloes and, to a lesser extent, mass accretion history. Our findings demonstrate that the properties of the smooth components of dark matter haloes are biased relative to the total halo mass.
What role the Large Magellanic Cloud (LMC), a dwarf irregular galaxy, plays in understanding infrared luminous galaxies is discussed. There are two main reasons the LMC may prove helpful. One, the LMC is only 55 kpc away, very nearby compared to much rarer high luminosity systems. Second, the environment in the LMC is distinctly different than in the Milky Way, at least those parts of the Milky Way interior to the sun, where most of the studies of massive star formation were concentrated. The LMC is an interacting system with a large amount of neutral hydrogen that is pushed around by the galaxy's encounter with the Milky Way. Perhaps a good understanding of star formation process in the LMC will provide guidance in the study of the infrared luminous galaxies. Two questions which will be addressed are: how is star formation in the LMC similar to the Milky Way Galaxy, and how is it different?
The sensitivity of direct detection experiments depends on the phase-space distribution of dark matter near the Sun, which can be modeled theoretically using cosmological hydrodynamical simulations of Milky Way–like galaxies. However, capturing the halo-to-halo variation in the local dark matter speeds—a necessary step for quantifying the astrophysical uncertainties that feed into experimental results—requires a sufficiently large sample of simulated galaxies, which has been a challenge. In this Letter, we quantify this variation with nearly 100 Milky Way–like galaxies from the tng50 simulation, the largest sample to date at this resolution. Moreover, we introduce a novel phase-space scaling procedure that endows every system with a reference frame that accurately reproduces the local standard-of-rest speed of our Galaxy, providing a principled way of extrapolating the simulation results to real-world data. The ensemble of predicted speed distributions is well characterized by the standard halo model, a Maxwell-Boltzmann distribution truncated at the escape speed, though the individual distributions can deviate from it, especially at high speeds. The dark matter–nucleon cross section limits placed by these speed distributions vary by ∼ 60% about the median. This places the 1𝜎 astrophysical uncertainty at or below the level of the systematic uncertainty of current ton-scale detectors, even down to the energy threshold. The predicted uncertainty remains unchanged when subselecting on those TNG 50 galaxies with merger histories similar to the Milky Way. Tabulated speed distributions, as well as Maxwell-Boltzmann fits, are provided for use in computing direct detection bounds or projecting sensitivities.
The bulge is the oldest component of the Milky Way. Since numerous simulations of Milky Way formation have predicted that the oldest stars at a given metallicity are found on tightly bound orbits, the Galaxy’s oldest stars are likely metal-poor stars in the inner bulge with small apocenters (i.e., R {sub apo} ≲ 4 kpc). In the past, stars with these properties have been impossible to find due to extreme reddening and extinction along the line of sight to the inner bulge. We have used the mid-infrared metal-poor star selection of Schlaufman & Casey (2014) on Spitzer/Galactic Legacy Infrared Mid-Plane Survey Extraordinaire data to overcome these problems and target candidate inner bulge metal-poor giants for moderate-resolution spectroscopy with Anglo-Australian Telescope/AAOmega. We used those data to select three confirmed metal-poor giants ([Fe/H] = −3.15, −2.56, −2.03) for follow-up high-resolution Magellan/Magellan Inamori Kyocera Echelle spectroscopy. A comprehensive orbit analysis using Gaia DR2 astrometry and our measured radial velocities confirms that these stars are tightly bound inner bulge stars. We determine the elemental abundances of each star and find high titanium and iron-peak abundances relative to iron in our most metal-poor star. We propose that the distinct abundance signature we detect is a product of nucleosynthesis in the Chandrasekhar-mass thermonuclear supernova of a CO white dwarf accreting from a helium star with a delay time of about 10 Myr. Even though chemical evolution is expected to occur quickly in the bulge, the intense star formation in the core of the nascent Milky Way was apparently able to produce at least one Chandrasekhar-mass thermonuclear supernova progenitor before chemical evolution advanced beyond [Fe/H]∼−3.