Search for: All records

Recent observational studies have uncovered a small number of very metal-poor (VMP) stars with cold kinematics in the Galactic disc and bulge. However, their origins remain enigmatic. We select a total of 138 Milky Way (MW) analogues from the TNG50 cosmological simulation based on their z = 0 properties: discy morphology, stellar mass, and local environment. In order to make more predictive statements for the MW, we further limit the spatial volume coverage of stellar populations in galaxies to that targeted by the upcoming 4MOST high-resolution survey of the Galactic disc and bulge. We find that across all galaxies, ∼20 per cent of VMP ([Fe/H] < −2) stars belong to the disc, with some analogues reaching 30 per cent. About 50 ± 10 per cent of the VMP disc stars are, on average, older than 12.5 Gyr and ∼70 ± 10 per cent come from accreted satellites. A large fraction of the VMP stars belong to the halo (∼70) and have a median age of 12 Gyr. Our results with the TNG50 cosmological simulation confirm earlier findings with simulations of fewer individual galaxies, and suggest that the stellar disc of the MW is very likely to host significant amounts of very- and extremely-metal-poor stars that, although mostly of ex situ origin, can also form in situ, reinforcing the idea of the existence of a primordial Galactic disc.

The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMSs, M>100 $\rm {\rm M}_{\odot }$) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50–500 $\rm {\rm M}_{\odot }$, including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne–Na and Mg–Al cycles. VMS with enhanced winds eject 5–10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anticorrelations, such as C–N and Na–O, in globular clusters. We find that for VMS 95 per cent of the total wind yields is produced on the main sequence, while only ∼ 5 per cent is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf–Rayet winds had been suggested to be responsible for galactic 26Al enrichment. Finally, 200 $\rm {\rm M}_{\odot }$ stars eject 100 times more of each heavy element in their winds than 50 $\rm {\rm M}_{\odot }$ stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50 $\rm {\rm M}_{\odot }$ stars.

Theoretical physical-chemical models for the formation of planetary systems depend on data quality for the Sun’s composition, that of stars in the solar neighbourhood, and of the estimated ’pristine’ compositions for stellar systems. The effective scatter and the observational uncertainties of elements within a few hundred parsecs from the Sun, even for the most abundant metals like carbon, oxygen and silicon, are still controversial. Here we analyse the stellar production and the chemical evolution of key elements that underpin the formation of rocky (C, O, Mg, Si) and gas/ice giant planets (C, N, O, S). We calculate 198 galactic chemical evolution (GCE) models of the solar neighbourhood to analyse the impact of different sets of stellar yields, of the upper mass limit for massive stars contributing to GCE (Mup) and of supernovae from massive-star progenitors which do not eject the bulk of the iron-peak elements (faint supernovae). Even considering the GCE variation produced via different sets of stellar yields, the observed dispersion of elements reported for stars in the Milky Way (MW) disc is not reproduced. Among others, the observed range of super-solar [Mg/Si] ratios, sub-solar [S/N], and the dispersion of up to 0.5 dex for [S/Si] challenge our models. The impact of varying Mup depends on the adopted supernova yields. Thus, observations do not provide a constraint on the Mup parametrization. When including the impact of faint supernova models in GCE calculations, elemental ratios vary by up to 0.1–0.2 dex in the MW disc; this modification better reproduces observations.

The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer simulations have not been able to fully address whether helium spectral features would form. This is because helium can not be sufficiently excited to form spectral features by thermal processes, but can be excited by collisions with non-thermal electrons, which most studies have neglected. We carry out a full non-local thermodynamic equilibrium radiative transfer simulation for an instance of a double detonation explosion model, and include a non-thermal treatment of fast electrons. We find a clear He i λ10830 feature which is strongest in the first few days after explosion and becomes weaker with time. Initially this feature is blended with the Mg ii λ10927 feature but over time separates to form a secondary feature to the blue wing of the Mg ii λ10927 feature. We compare our simulation to observations of iPTF13ebh, which showed a similar feature to the blue wing of the Mg ii λ10927 feature, previously identified as C i. Our simulation shows a good match to the evolution of this feature and we identify it as high velocity He i λ10830. This suggests that He i λ10830 could be a signature of the double detonation scenario.

We study the effects of the time-variable properties of thermonuclear X-ray bursts on modelling their millisecond-period burst oscillations. We apply the pulse profile modelling technique that is being used in the analysis of rotation-powered millisecond pulsars by the Neutron Star Interior Composition Explorer to infer masses, radii, and geometric parameters of neutron stars. By simulating and analysing a large set of models, we show that overlooking burst time-scale variability in temperatures and sizes of the hot emitting regions can result in substantial bias in the inferred mass and radius. To adequately infer neutron star properties, it is essential to develop a model for the time-variable properties or invest a substantial amount of computational time in segmenting the data into non-varying pieces. We discuss prospects for constraints from proposed future X-ray telescopes.

We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]),α-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low-αdisk, high-αdisk, and total population independently. The low-αdisk exhibits a negative radial metallicity gradient of −0.06 ± 0.001 dex kpc⁻¹, which flattens with distance from the midplane. The high-αdisk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]–[Fe/H] plane and in the [Mg/Fe]–age relation persist across the entire disk. The age estimates have typical uncertainties of ∼0.15 in log(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy.

Many of the short-lived radioactive nuclei that were present in the early solar system can be produced in massive stars. In the first paper in this series, we focused on the production of²⁶Al in massive binaries. In our second paper, we considered rotating single stars; two more short-lived radioactive nuclei,³⁶Cl and⁴¹Ca; and the comparison to the early solar system data. In this work, we update our previous conclusions by further considering the impact of binary interactions. We used the MESA stellar evolution code with an extended nuclear network to compute massive (10–80M_⊙), binary stars at various initial periods and solar metallicity (Z= 0.014), up to the onset of core collapse. The early solar system abundances of²⁶Al and⁴¹Ca can be matched self-consistently by models with initial masses ≥25M_⊙, while models with initial primary masses ≥35M_⊙can also match³⁶Cl. Almost none of the models provide positive net yields for¹⁹F, while for²²Ne the net yields are positive from 30M_⊙and higher. This leads to an increase by a factor of approximately 4 in the amount of²²Ne produced by a stellar population of binary stars, relative to single stars. In addition, besides the impact on the stellar yields, our 10M_⊙primary star undergoing Case A mass transfer ends its life as a white dwarf instead of as a core-collapse supernova. This demonstrates that binary interactions can also strongly impact the evolution of stars close to the supernova boundary.

Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multidimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M⊙ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution; thus, revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions, and compact remnants.

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.

We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the ppmstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a $1.2\, {\rm M}_{\odot }$ upper RGB star. We perform simulations for different grid resolutions (7683, 15363, and 28803), a range of driving luminosities, and two different stratifications (corresponding to the bump luminosity and the tip of the RGB). Our RGB tip simulations can be directly performed at the nominal luminosity, circumventing the need for extrapolations to lower luminosities. A rich, continuous spectrum of IGWs is observed, with a significant amount of total power contained at high wavenumbers. By following the time evolution of a passive dye in the stable layers, we find that IGW mixing in our simulations is weaker than predicted by a simple analytical prescription based on shear mixing and not efficient enough to explain the missing RGB extra mixing. However, we may be underestimating the efficiency of IGW mixing given that our simulations include a limited portion of the convective envelope. Quadrupling its radial extent compared to our fiducial set-up increases convective velocities by up to a factor 2 and IGW velocities by up to a factor 4. We also report the formation of a $\sim 0.2\, H_P$ penetration zone and evidence that IGWs are excited by plumes that overshoot into the stable layers.