One-dimensional (1D) stellar evolution models are widely used across various astrophysical fields, however they are still dominated by important uncertainties that deeply affect their predictive power. Among those, the merging of independent convective regions is a poorly understood phenomenon predicted by some 1D models but whose occurrence and impact in real stars remain very uncertain. Being an intrinsically multi-D phenomenon, it is challenging to predict the exact behaviour of shell mergers with 1D models. In this work, we conduct a detailed investigation of a multiple shell merging event in a 20 M_⊙ star using 3D hydrodynamic simulations. Making use of the active tracers for composition and the nuclear network included in the 3D model, we study the merging not only from a dynamical standpoint but also considering its nucleosynthesis and energy generation. Our simulations confirm the occurrence of the merging also in 3D, but reveal significant differences from the 1D case. Specifically, we identify entrainment and the erosion of stable regions as the main mechanisms that drive the merging, we predict much faster convective velocities compared to the mixing-length-theory velocities, and observe multiple burning phases within the same merged shell, with important effects for the chemical composition of the star, which presents a strongly asymmetric (dipolar) distribution. We expect that these differences will have important effects on the final structure of massive stars and thus their final collapse dynamics and possible supernova explosion, subsequently affecting the resulting nucleosynthesis and remnant.

Context. Recent observations by JWST have revealed supersolar ¹⁴N abundances in galaxies at very high redshift. On the other hand, these galaxies show subsolar metallicity. The observed N/O ratios are difficult to reproduce in the framework of chemical evolution models for the Milky Way. Aims. Our aim is to reproduce these high N/O ratios with chemical evolution models, assuming different histories of star formation triggering galactic winds, coupled with detailed nucleosynthesis prescriptions for ¹⁴N, ¹²C, ¹⁶O, and ⁵⁶Fe. Methods. We computed several models for small galaxies (10⁹−10¹⁰ M_⊙) with a high star formation efficiency and strong galactic winds. These winds are assumed to be differential, mainly carrying out the products of the explosion of core-collapse supernovae. Results. We find that only models with high star formation rates, a normal initial mass function, and differential galactic winds can reproduce the observed chemical abundances. We also find that, with the same assumptions about star formation and galactic winds, but with a very rapid formation resulting from fast gas infall, we can also reproduce the estimated ages of these objects. We find no necessity to invoke peculiar nucleosynthesis from population III stars, very massive stars, and supermassive stars.

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 multi-dimensional 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.

Recent studies show that the chemical evolution of Sr and Ba in the Galaxy can be explained if different production sites, hosting r- and s-processes, are taken into account. However, the question of unambiguously identifying these sites is still unsolved. Massive stars are shown to play an important role in the production of s-material if rotation is considered. In this work, we study in detail the contribution of rotating massive stars to the production of Sr and Ba, in order to explain their chemical evolution, but also to constrain the rotational behaviour of massive stars. A stochastic chemical evolution model was employed to reproduce the enrichment of the Galactic halo. We developed new methods for model-data comparison which help to objectively compare the stochastic results to the observations. We employed these methods to estimate the value of free parameters which describe the rotation of massive stars, assumed to be dependent on the stellar metallicity. We constrain the parameters using the observations for Sr and Ba. Employing these parameters for rotating massive stars in our stochastic model, we are able to correctly reproduce the chemical evolution of Sr and Ba, but also Y, Zr, and La. The data supports a decrease of both the mean rotational velocities and their dispersion with increasing metallicity. Our results show that a metallicity-dependent rotation is a necessary assumption to explain the s-process in massive stars. Our novel methods of model-data comparison represent a promising tool for future galactic chemical evolution studies.

Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multidimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multidimensional processes over short time-scales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here, we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars, and black holes physics.

Context. The chemical abundances of elements such as barium and the lanthanides are essential to understand the nucleosynthesis of heavy elements in the early Universe as well as the contribution of different neutron capture processes (for example slow versus rapid) at different epochs. Aims. The Chemical Evolution of R-process Elements in Stars (CERES) project aims to provide a homogeneous analysis of a sample of metal-poor stars ([Fe/H]<−1.5) to improve our understanding of the nucleosynthesis of neutron capture elements, in particular the r-process elements, in the early Galaxy. Methods. Our data consist of a sample of high resolution and high signal-to-noise ratio UVES spectra. The chemical abundances were derived through spectrum synthesis, using the same model atmospheres and stellar parameters as derived in the first paper of the CERES series. Results. We measured chemical abundances or upper limits of seven heavy neutron capture elements (Ba, La, Ce, Pr, Nd, Sm, and Eu) for a sample of 52 metal-poor giant stars. We estimated through the mean shift clustering algorithm that at [Ba/H]=−2.4 and [Fe/H]=−2.4 a variation in the trend of [X/Ba], with X=La,Nd,Sm,Eu, versus [Ba/H] occurs. This result suggests that, for [Ba/H]<−2.4, Ba nucleosynthesis in the Milky Way halo is primarily due to the r-process, while for [Ba/H]>−2.4 the effect of the s-process contribution begins to be visible. In our sample, stars with [Ba/Eu] compatible with a Solar System pure r-process value (hereafter, r-pure) do not show any particular trend compared to other stars, suggesting r-pure stars may form in similar environments to stars with less pure r-process enrichments. Conclusions. Homogeneous investigations of high resolution and signal-to-noise ratio spectra are crucial for studying the heavy elements formation, as they provide abundances that can be used to test nucleosynthesis models as well as Galactic chemical evolution models.