Search for: All records

In recent years, the study of the Milky Way has significantly advanced due to extensive spectroscopic surveys of its stars, complemented by astroseismic and astrometric data. However, it remains disjoint from recent advancements in understanding the physics of the Galactic interstellar medium (ISM). This paper introduces a new model for the chemical evolution of the Milky Way that can be constrained on stellar data, because it combines a state-of-the-art ISM model with a Milky Way stellar disc model. Utilizing a data set of red clump stars from APOGEE, known for their precise ages and metallicities, we concentrate on the last 6 billion years – a period marked by Milky Way’s secular evolution. We examine the oxygen abundance in the low-$\alpha$ disc stars relative to their ages and birth radii, validating or constraining critical ISM parameters that remain largely unexplored in extragalactic observations. The models that successfully reproduce the radius–metallicity distribution and the age–metallicity distribution of stars without violating existing ISM observations indicate a need for modest differential oxygen enrichment in Galactic outflows, meaning that the oxygen abundance of outflows is higher than the local ISM abundance, irrespective of outflow mass loading. The models also suggest somewhat elevated ISM gas velocity dispersion levels over the past 6 billion years compared to galaxies of similar mass. The extra turbulence necessary could result from energy from gas accretion onto the Galaxy, supernovae clustering in the ISM, or increased star formation efficiency per freefall time. This work provides a novel approach to constraining the Galactic ISM and outflows, leveraging the detailed insights available from contemporary Milky Way surveys.

 

We demonstrate the potential of the state-of-the-art OpenAI GPT-4 large language model to engage in meaningful interactions with Astronomy papers using in-context prompting. To optimize for efficiency, we employ a distillation technique that effectively reduces the size of the original input paper by 50%, while maintaining the paragraph structure and overall semantic integrity. We then explore the model’s responses using a multi-document context (ten distilled documents). Our findings indicate that GPT-4 excels in the multi-document domain, providing detailed answers contextualized within the framework of related research findings. Our results showcase the potential of large language models for the astronomical community, offering a promising avenue for further exploration, particularly the possibility of utilizing the models for hypothesis generation.

 

We present the first detailed chemical-abundance analysis of stars from the dwarf-galaxy stellar stream Wukong/LMS-1 covering a wide metallicity range ($-3.5 \lt \rm [Fe/H] \lesssim -1.3$). We find abundance patterns that are effectively indistinguishable from the bulk of Indus and Jhelum, a pair of smaller stellar streams proposed to be dynamically associated with Wukong/LMS-1. We confirmed a carbon-enhanced metal-poor star ($\rm [C/Fe] \gt +0.7$ and $\rm [Fe/H] \sim -2.9$) in Wukong/LMS-1 with strong enhancements in Sr, Y, and Zr, which is peculiar given its solar-level [Ba/Fe]. Wukong/LMS-1 stars have high abundances of α elements up to $\rm [Fe/H] \gtrsim -2$, which is expected for relatively massive dwarfs. Towards the high-metallicity end, Wukong/LMS-1 becomes α-poor, revealing that it probably experienced fairly standard chemical evolution. We identified a pair of N- and Na-rich stars in Wukong/LMS-1, reminiscent of multiple stellar populations in globular clusters. This indicates that this dwarf galaxy contained at least one globular cluster that was completely disrupted in addition to two intact ones previously known to be associated with Wukong/LMS-1, which is possibly connected to similar evidence found in Indus. From these ≥3 globular clusters, we estimate the total mass of Wukong/LMS-1 to be ${\approx }10^{10} \, \mathrm{M}_\odot$, representing ∼1 per cent of the present-day Milky Way. Finally, the [Eu/Mg] ratio in Wukong/LMS-1 continuously increases with metallicity, making this the first example of a dwarf galaxy where the production of r-process elements is clearly dominated by delayed sources, presumably neutron-star mergers.

 

We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding time-scales of τin = 1.01 ± 0.13 Gyr for GSE and $\tau _\text{in} = 3.08^{+3.19}_{-1.16}$ Gyr for Wukong/LMS-1. GSE formed stars for $\tau _\text{tot} = 5.40^{+0.32}_{-0.31}$ Gyr, sustaining star formation for ∼1.5–2 Gyr after its first infall into the Milky Way ∼10 Gyr ago. Our fit suggests that star formation lasted for $\tau _\text{tot} = 3.36^{+0.55}_{-0.47}$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass M⋆ predicted by simulations and semi-analytic models of galaxy formation. Our inferred values of the outflow mass-loading factor reasonably match $\eta \propto M_\star ^{-1/3}$ as predicted by galactic wind models. Our fitting method is based only on Poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes (20 ≤ N ≤ 2000) and measurement uncertainties (0.01 ≤ σ[α/Fe], σ[Fe/H] ≤ 0.5; $0.02 \le \sigma _{\log _{10}(\text{age})} \le 1$). Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space.

 

Carbon-enhanced metal-poor (CEMP) stars comprise almost a third of stars with [Fe/H] < −2, although their origins are still poorly understood. It is highly likely that one sub-class (CEMP-s stars) is tied to mass-transfer events in binary stars, while another sub-class (CEMP-no stars) are enriched by the nucleosynthetic yields of the first generations of stars. Previous studies of CEMP stars have primarily concentrated on the Galactic halo, but more recently they have also been detected in the thick disc and bulge components of the Milky Way. Gaia DR3 has provided an unprecedented sample of over 200 million low-resolution (R ≈ 50) spectra from the BP and RP photometers. Training on the CEMP catalogue from the SDSS/SEGUE database, we use XGBoost to identify the largest all-sky sample of CEMP candidate stars to date. In total, we find 58 872 CEMP star candidates, with an estimated contamination rate of 12 per cent. When comparing to literature high-resolution catalogues, we positively identify 60–68 per cent of the CEMP stars in the data, validating our results and indicating a high completeness rate. Our final catalogue of CEMP candidates spans from the inner to outer Milky Way, with distances as close as r ∼ 0.8 kpc from the Galactic centre, and as far as r > 30 kpc. Future higher resolution spectroscopic follow-up of these candidates will provide validations of their classification and enable investigations of the frequency of CEMP-s and CEMP-no stars throughout the Galaxy, to further constrain the nature of their progenitors.

 

Abstract Modern Galactic surveys have revealed an ancient merger that dominates the stellar halo of our galaxy (Gaia–Sausage–Enceladus, GSE). Using chemical abundances and kinematics from the H3 Survey, we identify 5559 halo stars from this merger in the radial range r Gal = 6–60kpc. We forward model the full selection function of H3 to infer the density profile of this accreted component of the stellar halo. We consider a general ellipsoid with principal axes allowed to rotate with respect to the galactocentric axes, coupled with a multiply broken power law. The best-fit model is a triaxial ellipsoid (axes ratios 10:8:7) tilted 25° above the Galactic plane toward the Sun and a doubly broken power law with breaking radii at 12 kpc and 28 kpc. The doubly broken power law resolves a long-standing dichotomy in literature values of the halo breaking radius, being at either ∼15 kpc or ∼30 kpc assuming a singly broken power law. N -body simulations suggest that the breaking radii are connected to apocenter pile-ups of stellar orbits, and so the observed double-break provides new insight into the initial conditions and evolution of the GSE merger. Furthermore, the tilt and triaxiality of the stellar halo could imply that a fraction of the underlying dark matter halo is also tilted and triaxial. This has important implications for dynamical mass modeling of the galaxy as well as direct dark matter detection experiments. 

Abstract Some studies of stars’ multielement abundance distributions suggest at least 5–7 significant dimensions, but others show that many elemental abundances can be predicted to high accuracy from [Fe/H] and [Mg/Fe] (or [Fe/H] and age) alone. We show that both propositions can be, and are, simultaneously true. We adopt a machine-learning technique known as normalizing flow to reconstruct the probability distribution of Milky Way disk stars in the space of 15 elemental abundances measured by APOGEE. Conditioning on T eff and log g minimizes the differential systematics. After further conditioning on [Fe/H] and [Mg/Fe], the residual scatter for most abundances is σ [ X /H] ≲ 0.02 dex, consistent with APOGEE’s reported statistical uncertainties of ∼0.01–0.015 dex and intrinsic scatter of 0.01–0.02 dex. Despite the small scatter, residual abundances display clear correlations between elements, which we show are too large to be explained by measurement uncertainties or by the finite sampling noise. We must condition on at least seven elements to reduce the correlations to a level consistent with the observational uncertainties. Our results demonstrate that cross-element correlations are a much more sensitive probe of a hidden structure than dispersion, and they can be measured precisely in a large sample even if the star-by-star measurement noise is comparable to the intrinsic scatter. We conclude that many elements have an independent story to tell, even for the mundane disk stars and elements produced by the core-collapse and Type Ia supernovae. The only way to learn these lessons is to measure the abundances directly, and not merely infer them. 

Recent observations of the stellar halo have uncovered the debris of an ancient merger, Gaia–Sausage–Enceladus (GSE), estimated to have occurred ≳8 Gyr ago. Follow-up studies have associated GSE with a large-scale tilt in the stellar halo that links two well-known stellar overdensities in diagonally opposing octants of the Galaxy (the Hercules–Aquila Cloud and Virgo Overdensity; HAC and VOD). In this paper, we study the plausibility of such unmixed merger debris persisting over several gigayears in the Galactic halo. We employ the simulated stellar halo from Naidu et al., which reproduces several key properties of the merger remnant, including the large-scale tilt. By integrating the orbits of these simulated stellar halo particles, we show that adoption of a spherical halo potential results in rapid phase mixing of the asymmetry. However, adopting a tilted halo potential preserves the initial asymmetry in the stellar halo for many gigayears. The asymmetry is preserved even when a realistic growing disk is added to the potential. These results suggest that HAC and VOD are long-lived structures that are associated with GSE and that the dark matter halo of the Galaxy is tilted with respect to the disk and aligned in the direction of HAC–VOD. Such halo–disk misalignment is common in modern cosmological simulations. Lastly, we study the relationship between the local and global stellar halo in light of a tilted global halo comprised of highly radial orbits. We find that the local halo offers a dynamically biased view of the global halo due to its displacement from the Galactic center.

 

Abstract Several lines of evidence suggest that the Milky Way underwent a major merger at z ∼ 2 with the Gaia-Sausage-Enceladus (GSE) galaxy. Here we use H3 Survey data to argue that GSE entered the Galaxy on a retrograde orbit based on a population of highly retrograde stars with chemistry similar to the largely radial GSE debris. We present the first tailored N -body simulations of the merger. From a grid of ≈500 simulations we find that a GSE with M ⋆ = 5 × 10 8 M ⊙ , M DM = 2 × 10 11 M ⊙ best matches the H3 data. This simulation shows that the retrograde stars are stripped from GSE’s outer disk early in the merger. Despite being selected purely on angular momenta and radial distributions, this simulation reproduces and explains the following phenomena: (i) the triaxial shape of the inner halo, whose major axis is at ≈35° to the plane and connects GSE’s apocenters; (ii) the Hercules-Aquila Cloud and the Virgo Overdensity, which arise due to apocenter pileup; and (iii) the 2 Gyr lag between the quenching of GSE and the truncation of the age distribution of the in situ halo, which tracks the lag between the first and final GSE pericenters. We make the following predictions: (i) the inner halo has a “double-break” density profile with breaks at both ≈15–18 kpc and 30 kpc, coincident with the GSE apocenters; and (ii) the outer halo has retrograde streams awaiting discovery at >30 kpc that contain ≈10% of GSE’s stars. The retrograde (radial) GSE debris originates from its outer (inner) disk—exploiting this trend, we reconstruct the stellar metallicity gradient of GSE (−0.04 ± 0.01 dex r 50 − 1 ). These simulations imply that GSE delivered ≈20% of the Milky Way’s present-day dark matter and ≈50% of its stellar halo.