Showing 1–50 of 160 results for author: Hirschi, R

Grids of stellar models with rotation VIII: Models from 1.7 to 500 $M_\odot$ at metallicity $Z = 10^{-5}$

Authors: Yves Sibony, Kendall G. Shepherd, Norhasliza Yusof, Raphael Hirschi, Caitlan Chambers, Sophie Tsiatsiou, Devesh Nandal, Luca Sciarini, Facundo D. Moyano, Jérôme Bétrisey, Gaël Buldgen, Cyril Georgy, Sylvia Ekström, Patrick Eggenberger, Georges Meynet

Abstract: Grids of stellar evolution models with rotation using the Geneva stellar evolution code (Genec) have been published for a wide range of metallicities. We introduce the last remaining grid of Genec models, with a metallicity of $Z=10^{-5}$. We study the impact of this extremely metal-poor initial composition on various aspects of stellar evolution, and compare it to the results from previous grids… ▽ More Grids of stellar evolution models with rotation using the Geneva stellar evolution code (Genec) have been published for a wide range of metallicities. We introduce the last remaining grid of Genec models, with a metallicity of $Z=10^{-5}$. We study the impact of this extremely metal-poor initial composition on various aspects of stellar evolution, and compare it to the results from previous grids at other metallicities. We provide electronic tables that can be used to interpolate between stellar evolution tracks and for population synthesis. Using the same physics as in the previous papers of this series, we computed a grid of stellar evolution models with Genec spanning masses between 1.7 and 500 $M_\odot$, with and without rotation, at a metallicity of $Z=10^{-5}$. Due to the extremely low metallicity of the models, mass-loss processes are negligible for all except the most massive stars. For most properties (such as evolutionary tracks in the Hertzsprung-Russell diagram, lifetimes, and final fates), the present models fit neatly between those previously computed at surrounding metallicities. However, specific to this metallicity is the very large production of primary nitrogen in moderately rotating stars, which is linked to the interplay between the hydrogen- and helium-burning regions. The stars in the present grid are interesting candidates as sources of nitrogen-enrichment in the early Universe. Indeed, they may have formed very early on from material previously enriched by the massive short-lived Population III stars, and as such constitute a very important piece in the puzzle that is the history of the Universe. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted in A&A. 20 pages, 14 figures, 3 tables

3D simulations of convective shell Neon-burning in a massive star

Authors: C. Georgy, F. Rizzuti, R. Hirschi, V. Varma, W. D. Arnett, C. Meakin, M. Mocak, A. StJ. Murphy, T. Rauscher

Abstract: The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulen… ▽ More The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same timescale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 17 pages, 20 figures, accepted for publication in MNRAS

Rapidly rotating Population III stellar models as a source of primary nitrogen

Authors: Sophie Tsiatsiou, Yves Sibony, Devesh Nandal, Luca Sciarini, Yutaka Hirai, Sylvia Ekstrom, Eoin Farrell, Laura Murphy, Arthur Choplin, Raphael Hirschi, Cristina Chiappini, Boyuan Liu, Volker Bromm, Jose Groh, Georges Meynet

Abstract: The first stars might have been fast rotators. This would have important consequences for their radiative, mechanical and chemical feedback. We discuss the impact of fast initial rotation on the evolution of massive Population III models and on their nitrogen and oxygen stellar yields. We explore the evolution of Population III stars with initial masses in the range of 9Msol < Mini < 120Msol start… ▽ More The first stars might have been fast rotators. This would have important consequences for their radiative, mechanical and chemical feedback. We discuss the impact of fast initial rotation on the evolution of massive Population III models and on their nitrogen and oxygen stellar yields. We explore the evolution of Population III stars with initial masses in the range of 9Msol < Mini < 120Msol starting with an initial rotation on the Zero Age Main Sequence equal to 70% of the critical one. We find that with the physics of rotation considered here, our rapidly-rotating Population III stellar models do not follow a homogeneous evolution. They lose very little mass in case mechanical winds are switched on when the surface rotation becomes equal or larger than the critical velocity. Impact on the ionising flux appears modest when compared to moderately-rotating models. Fast rotation favours, in models with initial masses above ~20Msol, the appearance of a very extended intermediate convective zone around the H-burning shell during the core He-burning phase. This shell has important consequences on the sizes of the He- and CO-cores and thus impacts the final fate of stars. Moreover, it has a strong impact on nucleosynthesis boosting the production of primary 14N. Fast initial rotation impacts significantly the chemical feedback of Population III stars. Observations of extremely metal-poor stars and/or starbursting regions are essential to provide constraints on the properties of the first stars. △ Less

Submitted 16 May, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: Accepted in A&A. Pages 17. Figures 15. Tables 2

Very Massive Star Models: I. Impact of Rotation and Metallicity and Comparisons with Observations

Authors: S. Martinet, G. Meynet, S. Ekström, C. Georgy, R. Hirschi

Abstract: In addition to being spectacular objects, Very Massive Stars (VMS) are suspected to have a tremendous impact on their environment and on the whole cosmic evolution. The nucleosynthesis both during their advanced stages and their final explosion may contribute greatly to the overall enrichment of the Universe. Their resulting supernovae are candidates for the most superluminous events and their ext… ▽ More In addition to being spectacular objects, Very Massive Stars (VMS) are suspected to have a tremendous impact on their environment and on the whole cosmic evolution. The nucleosynthesis both during their advanced stages and their final explosion may contribute greatly to the overall enrichment of the Universe. Their resulting supernovae are candidates for the most superluminous events and their extreme conditions also lead to very important radiative and mechanical feedback effects, from local to cosmic scale. We explore the impact of rotation and metallicity on the evolution of very massive stars across cosmic times. With the recent implementation of an equation of state in the GENEC stellar evolution code, appropriate for describing the conditions in the central regions of very massive stars in the advanced phases, we present new results on VMS evolution from Population III to solar metallicity. Low metallicity VMS models are highly sensitive to rotation, while the evolution of higher metallicity models is dominated by mass loss effects. The mass loss affects strongly their surface velocity evolution, breaking quickly at high metallicity while reaching the critical velocity for low metallicity models. The comparison to observed VMS in the LMC shows that the mass loss prescriptions used for these models are compatible with observed mass loss rates. In our framework for modelling rotation, our models of VMS need a high initial velocity to reproduce the observed surface velocities. The surface enrichment of these VMS is difficult to explain with only one initial composition, and could suggest multiple populations in the R136 cluster. At a metallicity typical of R136, only our non- or slowly rotating VMS models may produce Pair Instability supernovae. The most massive black holes that can be formed are less massive than about 60 M$_\odot$. △ Less

Submitted 4 September, 2023; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: 13 pages, 11 figures

Journal ref: A&A 679, A137 (2023)

Stellar Wind Yields of Very Massive Stars

Authors: Erin R. Higgins, Jorick S. Vink, Raphael Hirschi, Alison M. Laird, Gautham N. Sabhahit

Abstract: The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMS, M>100M) 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 com… ▽ More The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMS, M>100M) 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 - 500M, 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 anti-correlations, such as C-N and Na-O, in globular clusters. We find that for VMS 95% of the total wind yields is produced on the main sequence, while only ~5% 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, 200M stars eject 100 times more of each heavy element in their winds than 50M stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50M stars. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS. 14 pages, 10 figures

3D stellar evolution: hydrodynamic simulations of a complete burning phase in a massive star

Authors: F. Rizzuti, R. Hirschi, W. D. Arnett, C. Georgy, C. Meakin, A. StJ. Murphy, T. Rauscher, V. Varma

Abstract: 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… ▽ More 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$_\odot$ 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. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 12 pages, 15 figures. Accepted for publication in MNRAS

The nuclear reaction network WinNet

Authors: M. Reichert, C. Winteler, O. Korobkin, A. Arcones, J. Bliss, M. Eichler, U. Frischknecht, C. Fröhlich, R. Hirschi, M. Jacobi, J. Kuske, G. Martínez-Pinedo, D. Martin, D. Mocelj, T. Rauscher, F. -K. Thielemann

Abstract: We present the state-of-the-art single-zone nuclear reaction network WinNet that is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we desc… ▽ More We present the state-of-the-art single-zone nuclear reaction network WinNet that is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication, WinNet is publicly available and open source at GitHub and Zenodo. △ Less

Submitted 31 August, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Spin-down and reduced mass loss in early-type stars with large-scale magnetic fields

Authors: Z. Keszthelyi, A. de Koter, Y. Götberg, G. Meynet, S. A. Brands, V. Petit, M. Carrington, A. David-Uraz, S. T. Geen, C. Georgy, R. Hirschi, J. Puls, K. J. Ramalatswa, M. E. Shultz, A. ud-Doula

Abstract: Magnetism can greatly impact the evolution of stars. In some stars with OBA spectral types there is direct evidence via the Zeeman effect for stable, large-scale magnetospheres, which lead to the spin-down of the stellar surface and reduced mass loss. So far, a comprehensive grid of stellar structure and evolution models accounting for these effects was lacking. For this reason, we computed and st… ▽ More Magnetism can greatly impact the evolution of stars. In some stars with OBA spectral types there is direct evidence via the Zeeman effect for stable, large-scale magnetospheres, which lead to the spin-down of the stellar surface and reduced mass loss. So far, a comprehensive grid of stellar structure and evolution models accounting for these effects was lacking. For this reason, we computed and studied models with two magnetic braking and two chemical mixing schemes in three metallicity environments with the MESA software instrument. We find notable differences between the subgrids, which affects the model predictions and thus the detailed characterisation of stars. We are able to quantify the impact of magnetic fields in terms of preventing quasi-chemically homogeneous evolution and producing slowly-rotating, nitrogen-enriched ("Group 2") stars. Our model grid is fully open access and open source. △ Less

Submitted 13 November, 2022; originally announced November 2022.

Comments: to appear in "Winds of Stars and Exoplanets" Proceedings IAUS 370, 2022, eds.: A.A. Vidotto, L. Fossati, J.S. Vink

The effects of surface fossil magnetic fields on massive star evolution: IV. Grids of models at Solar, LMC, and SMC metallicities

Authors: Z. Keszthelyi, A. de Koter, Y. Götberg, G. Meynet, S. A. Brands, V. Petit, M. Carrington, A. David-Uraz, S. T. Geen, C. Georgy, R. Hirschi, J. Puls, K. J. Ramalatswa, M. E. Shultz, A. ud-Doula

Abstract: Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to… ▽ More Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3-60 M$_\odot$), surface equatorial magnetic field strength (0-50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly-rotating, nitrogen-enriched ("Group 2") stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of i) revisiting the derived stellar parameters of known magnetic stars, and ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes. △ Less

Submitted 28 October, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS. A full reproduction package is available on Zenodo at https://doi.org/10.5281/zenodo.7069766

Realistic 3D hydrodynamics simulations find significant turbulent entrainment in massive stars

Authors: F. Rizzuti, R. Hirschi, C. Georgy, W. D. Arnett, C. Meakin, A. StJ. Murphy

Abstract: Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional 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 multi-dimensional processes over short timescal… ▽ More Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional 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 multi-dimensional processes over short timescales, 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. △ Less

Submitted 16 August, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 7 pages, 6 figures, accepted for publication in MNRAS

The impact of $^{17}$O$+α$ reaction rate uncertainties on the s-process in rotating massive stars

Authors: J. Frost-Schenk, P. Adsley, A. M. Laird, R. Longland, C. Angus, C. Barton, A. Choplin, C. Aa. Diget, R. Hirschi, C. Marshall, F. Portillo Chaves, K. Setoodehnia

Abstract: Massive stars are crucial to galactic chemical evolution for elements heavier than iron. Their contribution at early times in the evolution of the Universe, however, is unclear due to poorly constrained nuclear reaction rates. The competing $^{17}$O($α,γ$)$^{21}$Ne and $^{17}$O($α,n$)$^{20}$Ne reactions strongly impact weak s-process yields from rotating massive stars at low metallicities. Abundan… ▽ More Massive stars are crucial to galactic chemical evolution for elements heavier than iron. Their contribution at early times in the evolution of the Universe, however, is unclear due to poorly constrained nuclear reaction rates. The competing $^{17}$O($α,γ$)$^{21}$Ne and $^{17}$O($α,n$)$^{20}$Ne reactions strongly impact weak s-process yields from rotating massive stars at low metallicities. Abundant $^{16}$O absorbs neutrons, removing flux from the s-process, and producing $^{17}$O. The $^{17}$O($α,n$)$^{20}$Ne reaction releases neutrons, allowing continued s-process nucleosynthesis, if the $^{17}$O($α,γ$)$^{21}$Ne reaction is sufficiently weak. While published rates are available, they are based on limited indirect experimental data for the relevant temperatures and, more importantly, no uncertainties are provided. The available nuclear physics has been evaluated, and combined with data from a new study of astrophysically relevant $^{21}$Ne states using the $^{20}$Ne($d,p$)$^{21}$Ne reaction. Constraints are placed on the ratio of the ($α,n$)/($α,γ$) reaction rates with uncertainties on the rates provided for the first time. The new rates favour the ($α,n$) reaction and suggest that the weak s-process in rotating low-metallicity stars is likely to continue up to barium and, within the computed uncertainties, even to lead. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: The MNRAS way of doing a sideways table forces a new page but apparently this is a known behaviour

Very massive star winds as sources of the short-lived radioactive isotope $^{26}$Al

Authors: S. Martinet, G. Meynet, D. Nandal, S. Ekström, C. Georgy, L. Haemmerlé, R. Hirschi, N. Yusof, M. Gounelle, V. Dwarkadas

Abstract: The $^{26}$Al short-lived radioactive nuclide is the source of the observed galactic diffuse $γ$-ray emission at 1.8 MeV. While different sources of $^{26}$Al have been explored, such as AGB stars, massive stars winds, and supernovae, the contribution of very massive stars has never been studied. We study the stellar wind contribution of very massive stars, i.e stars with initial masses between 15… ▽ More The $^{26}$Al short-lived radioactive nuclide is the source of the observed galactic diffuse $γ$-ray emission at 1.8 MeV. While different sources of $^{26}$Al have been explored, such as AGB stars, massive stars winds, and supernovae, the contribution of very massive stars has never been studied. We study the stellar wind contribution of very massive stars, i.e stars with initial masses between 150 and 300 M$_\odot$, to the enrichment in $^{26}$Al of the galactic interstellar medium. We discuss the production of $^{26}$Al by studying rotating and non-rotating very massive stellar models with initial masses between 150 and 300 M$_\odot$ for metallicities Z=0.006, 0.014, and 0.020. We confront this result to a simple Milky Way model taking into account both the metallicity and the star formation rate gradients. We obtain that very massive stars in the Z=0.006-0.020 metallicity range might be very significant contributors to the $^{26}$Al enrichment of the interstellar medium. Typically, the contribution of the winds of massive stars to the total quantity of $^{26}$Al in the Galaxy increases by 150\% when very massive stars are considered. Very massive stars, despite their rarity, might be important contributors to $^{26}$Al and overall very important actors for nucleosynthesis in the Galaxy. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 10 pages, 6 figures

Journal ref: A&A 664, A181 (2022)

The p-process in exploding rotating massive stars

Authors: A. Choplin, S. Goriely, R. Hirschi, N. Tominaga, G. Meynet

Abstract: The p-process nucleosynthesis can explain proton-rich isotopes that are heavier than iron, which are observed in the Solar System, but discrepancies still persist and important questions concerning the astrophysical site(s) of the p-process remain unanswered. We investigate how the p-process operates in exploding rotating massive stars that have experienced an enhanced s-process nucleosynthesis du… ▽ More The p-process nucleosynthesis can explain proton-rich isotopes that are heavier than iron, which are observed in the Solar System, but discrepancies still persist and important questions concerning the astrophysical site(s) of the p-process remain unanswered. We investigate how the p-process operates in exploding rotating massive stars that have experienced an enhanced s-process nucleosynthesis during their life through rotational mixing. We computed 25 $M_{\odot}$ stellar models at a metallicity of $Z=10^{-3}$ with different initial rotation velocities and rates for the uncertain $^{17}$O($α$,$γ$)$^{21}$Ne reaction. The nucleosynthesis calculation, followed with a network of 737 isotopes, was coupled to stellar evolution, and the p-process nucleosynthesis was calculated in post-processing during both the final evolutionary stages and spherical explosions of various energies. In our models, the p-nuclides are mainly synthesized during the explosion, but not much during the ultimate hydrostatic burning stages. The p-process yields mostly depend on the initial number of trans-iron seeds, which in turn depend on the initial rotation. We found that the impact of rotation on the p-process is comparable to the impact of rotation on the s-process. From no to fast rotation, the s-process yields of nuclides with mass number $A<140$ increase by $3-4$ dex, and so do the p-process yields. Fast rotation with a lower $^{17}$O($α,γ$) rate significantly produces s- and p-nuclides with $A\geq140$. Our results suggest that the contribution of core-collapse supernovae from massive stars to the solar (and Galactic) p-nuclei has been underestimated in the past, and more specifically, that the contribution from massive stars with sub-solar metallicities may even dominate. A more detailed study including stellar models with a wide range of masses and metallicities remains to be performed. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 12 pages, 13 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 661, A86 (2022)

Grids of stellar models with rotation VI: Models from 0.8 to 120 $M_\odot$ at a metallicity Z = 0.006

Authors: Patrick Eggenberger, Sylvia Ekström, Cyril Georgy, Sébastien Martinet, Camilla Pezzotti, Devesh Nandal, Georges Meynet, Gaël Buldgen, Sébastien Salmon, Lionel Haemmerlé, André Maeder, Raphael Hirschi, Norhasliza Yusof, José Groh, Eoin Farrell, Laura Murphy, Arthur Choplin

Abstract: Context: Grids of stellar models, computed with the same physical ingredients, allow one to study the impact of a given physics on a broad range of initial conditions and are a key ingredient for modeling the evolution of galaxies. Aims: We present a grid of single star models for masses between 0.8 and 120 $M_\odot$, with and without rotation for a mass fraction of heavy element Z=0.006, represen… ▽ More Context: Grids of stellar models, computed with the same physical ingredients, allow one to study the impact of a given physics on a broad range of initial conditions and are a key ingredient for modeling the evolution of galaxies. Aims: We present a grid of single star models for masses between 0.8 and 120 $M_\odot$, with and without rotation for a mass fraction of heavy element Z=0.006, representative of the Large Magellanic Cloud (LMC). Methods: We used the Geneva stellar evolution code. The evolution was computed until the end of the central carbon-burning phase, the early asymptotic giant branch phase, or the core helium-flash for massive, intermediate, and low mass stars, respectively. Results: The outputs of the present stellar models are well framed by the outputs of the two grids obtained by our group for metallicities above and below the one considered here. The models of the present work provide a good fit to the nitrogen surface enrichments observed during the main sequence for stars in the LMC with initial masses around 15 $M_\odot$. They also reproduce the slope of the luminosity function of red supergiants of the LMC well, which is a feature that is sensitive to the time-averaged mass loss rate over the red supergiant phase. The most massive black hole that can be formed from the present models at Z=0.006 is around 55 $M_\odot$. No model in the range of mass considered will enter into the pair-instability supernova regime, while the minimal mass to enter the region of pair pulsation instability is around 60 $M_\odot$ for the rotating models and 85 $M_\odot$ for the nonrotating ones. Conclusions: The present models are of particular interest for comparisons with observations in the LMC and also in the outer regions of the Milky Way. We provide public access to numerical tables that can be used for computing interpolated tracks and for population synthesis studies. △ Less

Submitted 31 January, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 11 pages, 11 figures, published in A&A

Journal ref: 2021, A&A 652, A137

Grids of stellar models with rotation VII: Models from 0.8 to 300 M$_\odot$ at super-solar metallicity (Z = 0.020)

Authors: Norhasliza Yusof, Raphael Hirschi, Patrick Eggenberger, Sylvia Ekström, Cyril Georgy, Yves Sibony, Paul A. Crowther, Georges Meynet, Hasan Abu Kassim, Wan Aishah Wan Harun, André Maeder, Jose H. Groh, Eoin Farrell, Laura Murphy

Abstract: We present a grid of stellar models at super-solar metallicity (Z = 0.020) extending the previous grids of Geneva models at solar and sub-solar metallicities. A metallicity of Z = 0.020 was chosen to match that of the inner Galactic disk. A modest increase of 43% (=0.02/0.014) in metallicity compared to solar models means that the models evolve similarly to solar models but with slightly larger ma… ▽ More We present a grid of stellar models at super-solar metallicity (Z = 0.020) extending the previous grids of Geneva models at solar and sub-solar metallicities. A metallicity of Z = 0.020 was chosen to match that of the inner Galactic disk. A modest increase of 43% (=0.02/0.014) in metallicity compared to solar models means that the models evolve similarly to solar models but with slightly larger mass loss. Mass loss limits the final total masses of the super-solar models to 35 M$_\odot$ even for stars with initial masses much larger than 100 M$_\odot$. Mass loss is strong enough in stars above 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) to remove the entire hydrogen-rich envelope. Our models thus predict SNII below 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) and SNIb (possibly SNIc) above that. We computed both isochrones and synthetic clusters to compare our super-solar models to the Westerlund 1 (Wd1) massive young cluster. A synthetic cluster combining rotating and non-rotating models with an age spread between log10 (age/yr) = 6.7 and 7.0 is able to reproduce qualitatively the observed populations of WR, RSG and YSG stars in Wd1, in particular their simultaneous presence at log10(L/L$_\odot$) = 5-5.5. The quantitative agreement is imperfect and we discuss the likely causes: synthetic cluster parameters, binary interactions, mass loss and their related uncertainties. In particular, mass loss in the cool part of the HRD plays a key role. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: 15 pages, 12 figures. Accepted for publication in MNRAS

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Authors: R. Andrassy, J. Higl, H. Mao, M. Mocák, D. G. Vlaykov, W. D. Arnett, I. Baraffe, S. W. Campbell, T. Constantino, P. V. F. Edelmann, T. Goffrey, T. Guillet, F. Herwig, R. Hirschi, L. Horst, G. Leidi, C. Meakin, J. Pratt, F. Rizzuti, F. K. Roepke, P. Woodward

Abstract: Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o… ▽ More Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study. △ Less

Submitted 26 January, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 20 pages, 19 figures, new Fig. B.2, new clarifications, electronic materials updated, accepted for publication in A&A

Journal ref: A&A 659, A193 (2022)

UVES analysis of red giants in the bulge globular cluster NGC 6522

Authors: B. Barbuy, E. Cantelli, L. Muniz, S. O. Souza, C. Chiappini, R. Hirschi, G. Cescutti, M. Pignatari, S. Ortolani, L. Kerber, F. F. S. Maia, E. Bica, E. Depagne

Abstract: NGC 6522 is a moderately metal-poor bulge globular cluster ([Fe/H]$\sim$$-$1.0), and it is a well-studied representative among a number of moderately metal-poor blue horizontal branch clusters located in the bulge. The NGC 6522 abundance pattern can give hints on the earliest chemical enrichment in the central Galaxy. The aim of this study is to derive abundances of the light elements C and N; alp… ▽ More NGC 6522 is a moderately metal-poor bulge globular cluster ([Fe/H]$\sim$$-$1.0), and it is a well-studied representative among a number of moderately metal-poor blue horizontal branch clusters located in the bulge. The NGC 6522 abundance pattern can give hints on the earliest chemical enrichment in the central Galaxy. The aim of this study is to derive abundances of the light elements C and N; alpha elements O, Mg, Si, Ca, and Ti; odd-Z elements Na and Al; neutron-capture elements Y, Zr, Ba, La, and Nd; and the r-process element Eu. We verify if there are first- and second-generation stars: we find clear evidence of Na-Al, Na-N, and Mg-Al correlations, while we cannot identify the Na-O anti-correlation from our data. High-resolution spectra of six red giants in the bulge globular cluster NGC 6522 were obtained at the 8m VLT UT2-Kueyen telescope in FLAMES+UVES configuration. In light of Gaia data, it turned out that two of them are non-members, but these were also analysed. Spectroscopic parameters were derived through the excitation and ionisation equilibrium of FeI and FeII lines from UVES spectra. The abundances were obtained with spectrum synthesis. The present analysis combined with previous UVES results gives a mean radial velocity of vrhel = -15.62+-7.7 km.s-1 and a metallicity of [Fe/H] = -1.05+-0.20 for NGC 6522. Mean abundances of alpha elements for the present four member stars are enhanced with [O/Fe]=+0.38, [Mg/Fe]=+0.28, [Si/Fe]=+0.19, and [Ca/Fe]=+0.13, together with the iron-peak element [Ti/Fe]=+0.13, and the r-process element [Eu/Fe]=+0.40.The neutron-capture elements Y, Zr, Ba, and La show enhancements in the +0.08 < [Y/Fe] < +0.90, 0.11 < [Zr/Fe] < +0.50, 0.00 < [Ba/Fe] < +0.63, 0.00 < [La/Fe] < +0.45, and -0.10 < [Nd/Fe] < +0.70 ranges. We also discuss the spread in heavy-element abundances. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: Accepted for publication in Astronomy & Astrophysics, 27 pages, 18 figures

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

Authors: L. Horst, R. Hirschi, P. V. F. Edelmann, R. Andrassy, F. K. Roepke

Abstract: Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers. Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional si… ▽ More Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers. Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M$_\odot$ star model. The results obtained with the low-Mach AUSM$^{+}$-up solver were compared to results when using its non low-Mach variant AUSM$_\mathrm{B}^{+}$-up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from $180 \times 90^2$ to $810 \times 540^2$ cells and the nuclear energy release was boosted by factors of $3 \times 10^3$, $1 \times 10^4$, and $3 \times 10^4$ to study the dependence of the results on these parameters. Results. The boosted energy input results in convection at Mach numbers in the range of $10^{-2}$ to $10^{-3}$. Standard mixing-length theory (MLT) predicts convective velocities of about $1.6 \times 10^{-4}$ if no boosting is applied. Simulations with AUSM$^{+}$-up show a Kolmogorov-like inertial range in the kinetic energy spectrum that extends further toward smaller scales compared with its non low-Mach variant. The kinetic energy dissipation of the AUSM$^{+}$-up solver already converges at a lower resolution compared to AUSM$^{+}_{\mathrm{B}}$ -up. The extracted entrainment rates at the boundaries of the convection zone are well represented by the bulk Richardson entrainment law and the corresponding fitting parameters are in agreement with published results for carbon shell burning. △ Less

Submitted 5 July, 2021; originally announced July 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 653, A55 (2021)

Evolution of Wolf-Rayet stars as black hole progenitors

Authors: Erin R. Higgins, Andreas A. C. Sander, Jorick S. Vink, Raphael Hirschi

Abstract: Evolved Wolf-Rayet stars form a key aspect of massive star evolution, and their strong outflows determine their final fates. In this study, we calculate grids of stellar models for a wide range of initial masses at five metallicities (ranging from solar down to just 2% solar). We compare a recent hydrodynamically-consistent wind prescription with two earlier frequently-used wind recipes in stellar… ▽ More Evolved Wolf-Rayet stars form a key aspect of massive star evolution, and their strong outflows determine their final fates. In this study, we calculate grids of stellar models for a wide range of initial masses at five metallicities (ranging from solar down to just 2% solar). We compare a recent hydrodynamically-consistent wind prescription with two earlier frequently-used wind recipes in stellar evolution and population synthesis modelling, and we present the ranges of maximum final masses at core He-exhaustion for each wind prescription and metallicity Z. Our model grids reveal qualitative differences in mass-loss behaviour of the wind prescriptions in terms of "convergence". Using the prescription from Nugis & Lamers the maximum stellar black hole is found to converge to a value of 20-30Msun, independent of host metallicity, however when utilising the new physically-motivated prescription from Sander & Vink there is no convergence to a maximum black hole mass value. The final mass is simply larger for larger initial He-star mass, which implies that the upper black hole limit for He-stars below the pair-instability gap is set by prior evolution with mass loss, or the pair instability itself. Quantitatively, we find the critical Z for pair-instability (Z_PI) to be as high as 50% Zsolar, corresponding to the host metallicity of the LMC. Moreover, while the Nugis & Lamers prescription would not predict any black holes above the approx 130Msun pair-instability limit, with Sander & Vink winds included, we demonstrate a potential channel for very massive helium stars to form such massive black holes at ~2% Zsolar or below. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS. 17 pages, 16 figures

Convective core entrainment in 1D main sequence stellar models

Authors: L. J. A. Scott, R. Hirschi, C. Georgy, W. D. Arnett, C. Meakin, E. A. Kaiser, S. Ekström, N. Yusof

Abstract: 3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M$_{\odot}$ at solar metallicity ($Z=0.014$) and compare them to previous generations of mod… ▽ More 3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M$_{\odot}$ at solar metallicity ($Z=0.014$) and compare them to previous generations of models and observations on the main sequence. The boundary penetrability, quantified by the bulk Richardson number, $Ri_{\mathrm{B}}$, varies with mass and to a smaller extent with time. The variation of $Ri_{\mathrm{B}}$ with mass is due to the mass dependence of typical convective velocities in the core and hence the luminosity of the star. The chemical gradient above the convective core dominates the variation of $Ri_{\mathrm{B}}$ with time. An entrainment law method can therefore explain the apparent mass dependence of convective boundary mixing through $Ri_{\mathrm{B}}$. New models including entrainment can better reproduce the mass dependence of the main sequence width using entrainment law parameters $A \sim 2 \times 10^{-4}$ and $n=1$. We compare these empirically constrained values to the results of 3D hydrodynamics simulations and discuss implications. △ Less

Submitted 11 March, 2021; v1 submitted 10 March, 2021; originally announced March 2021.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

Convective core sizes in rotating massive stars: I. Constraints from solar metallicity OB field stars

Authors: S. Martinet, G. Meynet, S. Ekström, S. Simón-Díaz, G. Holgado, N. Castro, C. Georgy, P. Eggenberger, G. Buldgen, S. Salmon, R. Hirschi, J. Groh, E. Farrell, L. Murphy

Abstract: Spectroscopic studies of Galactic O and B stars show that many stars with masses above 8 M$_{\odot}$ are observed in the HR diagram just beyond the Main-Sequence (MS) band predicted by stellar models computed with a moderate overshooting. This may be an indication that the convective core sizes in stars in the upper part of the HR diagram are larger than predicted by these models. Combining stella… ▽ More Spectroscopic studies of Galactic O and B stars show that many stars with masses above 8 M$_{\odot}$ are observed in the HR diagram just beyond the Main-Sequence (MS) band predicted by stellar models computed with a moderate overshooting. This may be an indication that the convective core sizes in stars in the upper part of the HR diagram are larger than predicted by these models. Combining stellar evolution models and spectroscopic parameters derived for a large sample of Galactic O and B stars, including brand new information about their projected rotational velocities, we reexamine the question of the convective core size in MS massive stars. We confirm that for stars more massive than about 8 M$_{\odot}$, the convective core size at the end of the MS phase increases more rapidly with the mass than in models computed with a constant step overshoot chosen to reproduce the main sequence width in the low mass range (around 2 M$_{\odot}$). This conclusion is valid for both the cases of non-rotating models and rotating models either with a moderate or a strong angular momentum transport. The increase of the convective core mass with the mass obtained from the TAMS position is, however, larger than the one deduced from the surface velocity drop for masses above about 15 M$_{\odot}$. Although observations available at the moment cannot decide what is the best choice between the core sizes given by the TAMS and the velocity drop, we discuss different methods to get out of this dilemma. At the moment, comparisons with eclipsing binaries seem to favor the solution given by the velocity drop. While we confirm the need for larger convective cores at higher masses, we find tensions in-between different methods for stars more massive than 15 M$_{\odot}$. The use of single-aged stellar populations (non-interacting binaries or stellar clusters) would be a great asset to resolve this tension. △ Less

Submitted 5 March, 2021; originally announced March 2021.

Comments: 13 pages, 8 figures

Journal ref: A&A 648, A126 (2021)

Constraints on stellar rotation from the evolution of Sr and Ba in the Galactic halo

Authors: F. Rizzuti, G. Cescutti, F. Matteucci, A. Chieffi, R. Hirschi, M. Limongi, A. Saro

Abstract: 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 t… ▽ More 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. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: 13 pages, 7 figures

Grids of stellar models with rotation: V. Models from 1.7 to 120 Msun at zero metallicity

Authors: Laura J. Murphy, Jose H. Groh, Sylvia Ekstrom, Georges Meynet, Camila Pezzotti, Cyril Georgy, Arthur Choplin, Patrick Eggenberger, Eoin Farrell, Lionel Haemmerle, Raphael Hirschi, Andre Maeder, Sebastien Martinet

Abstract: Understanding the nature of the first stars is key to understanding the early universe. With new facilities such as JWST we may soon have the first observations of the earliest stellar populations, but to understand these observations we require detailed theoretical models. Here we compute a grid of stellar evolution models using the Geneva code with the aim to improve our understanding of the evo… ▽ More Understanding the nature of the first stars is key to understanding the early universe. With new facilities such as JWST we may soon have the first observations of the earliest stellar populations, but to understand these observations we require detailed theoretical models. Here we compute a grid of stellar evolution models using the Geneva code with the aim to improve our understanding of the evolution of zero-metallicity stars, with particular interest in how rotation affects surface properties, interior structure, and metal enrichment. We produce a range of models of initial masses (Mini) from 1.7 Msun to 120 Msun, focusing on massive models of 9 Msun < Mini < 120 Msun. Our grid includes models with and without rotation, with rotating models having an initial velocity of 40% of the critical velocity. We find that rotation strongly impacts the evolution of the first stars, mainly through increased core size and stronger H-burning shells during core He-burning. Without radiative mass loss, angular momentum builds at the surface in rotating models, thus models of initial masses Mini > 60 Msun reach critical rotation on the main sequence and experience mass loss. We find that rotational mixing strongly affects metal enrichment, but does not always increase metal production as we see at higher metallicities. This is because rotation leads to an earlier CNO boost to the H shell during He-burning, which may hinder metal enrichment depending on initial mass and rotational velocity. Electronic tables of this new grid of Population III models are publicly available. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 20 pages, 15 figures. Accepted for publication in MNRAS

Is GW190521 the merger of black holes from the first stellar generations?

Authors: Eoin J. Farrell, Jose H. Groh, Raphael Hirschi, Laura Murphy, Etienne Kaiser, Sylvia Ekström, Cyril Georgy, Georges Meynet

Abstract: GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair-instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime and produce a black hole with a mass in the mass gap by fallback. We presen… ▽ More GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair-instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the Geneva and MESA stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0-0.0004 have three properties which favour higher BH masses as compared to higher metallicity models. These are (i) lower mass-loss rates during the post-MS phase, (ii) a more compact star disfavouring binary interaction and (iii) possible H-He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70-75 Msun. Uncertainties related to convective mixing, mass loss, H-He shell interactions and pair-instability pulsations may increase this limit to ~85 Msun. △ Less

Submitted 6 January, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 6 pages, 3 figures, accepted for publication in MNRAS Letters

Relative Importance of Convective Uncertainties in Massive Stars

Authors: Etienne A. Kaiser, Raphael Hirschi, W. David Arnett, Cyril Georgy, Laura J. A. Scott, Andrea Cristini

Abstract: In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary th… ▽ More In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses $15$, $20$ and $25\,\rm{M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main-sequence with more CBM, which is more in agreement with observations. Furthermore during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the $^{12}$C to $^{16}$O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70\%$ for the core masses and the total mass of the star. △ Less

Submitted 17 July, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 24 pages, 13 figures, 3 tables, accepted for publication in MNRAS. Stellar tracks and profiles are available at http://www.astro.keele.ac.uk/shyne/datasets/convective-boundary-mixing-uncertainties

Re-evaluation of the $^{22}$Ne($α,γ$)$^{26}$Mg and $^{22}$Ne($α,n$)$^{25}$Mg reaction rates

Authors: Philip Adsley, Umberto Battino, Andreas Best, Antonio Caciolli, Alessandra Guglielmetti, Gianluca Imbriani, Heshani Jayatissa, Marco La Cognata, Livio Lamia, Eliana Masha, Cristian Massimi, Sara Palmerini, Ashley Tattersall, Raphael Hirschi

Abstract: The competing $^{22}$Ne($α,γ$)$^{26}$Mg and $^{22}$Ne($α,n$)$^{25}$Mg reactions control the production of neutrons for the weak $s$-process in massive and AGB stars. In both systems, the ratio between the corresponding reaction rates strongly impacts the total neutron budget and strongly influences the final nucleosynthesis. The $^{22}$Ne($α,γ$)$^{26}$Mg and $^{22}$Ne($α,n$)$^{25}$Mg reaction rate… ▽ More The competing $^{22}$Ne($α,γ$)$^{26}$Mg and $^{22}$Ne($α,n$)$^{25}$Mg reactions control the production of neutrons for the weak $s$-process in massive and AGB stars. In both systems, the ratio between the corresponding reaction rates strongly impacts the total neutron budget and strongly influences the final nucleosynthesis. The $^{22}$Ne($α,γ$)$^{26}$Mg and $^{22}$Ne($α,n$)$^{25}$Mg reaction rates was re-evaluated by using newly available information on $^{26}$Mg given by various recent experimental studies. Evaluations of The evaluated $^{22}$Ne($α,γ$)$^{26}$Mg reaction rate remains substantially similar to that of Longland {\it et al.} but, including recent results from Texas A\&M, the $^{22}$Ne($α,n$)$^{25}$Mg reaction rate is lower at a range of astrophysically important temperatures. Stellar models computed with NEWTON and MESA predict decreased production of the weak branch $s$-process due to the decreased efficiency of $^{22}$Ne as a neutron source. Using the new reaction rates in the MESA model results in $^{96}$Zr/$^{94}$Zr and $^{135}$Ba/$^{136}$Ba ratios in much better agreement with the measured ratios from presolar SiC grains. △ Less

Submitted 21 March, 2021; v1 submitted 29 May, 2020; originally announced May 2020.

Comments: 23 pages, 15 figures, updated

Nucleosynthesis in early rotating massive stars and chemical composition of CEMP stars

Authors: Arthur Choplin, Raphael Hirschi

Abstract: The first massive stars triggered the onset of chemical evolution by releasing the first metals (elements heavier than helium) in the Universe. The nature of these stars and how the early chemical enrichment took place is still largely unknown. Rotational-induced mixing in the stellar interior can impact the nucleosynthesis during the stellar life of massive stars and lead to stellar ejecta having… ▽ More The first massive stars triggered the onset of chemical evolution by releasing the first metals (elements heavier than helium) in the Universe. The nature of these stars and how the early chemical enrichment took place is still largely unknown. Rotational-induced mixing in the stellar interior can impact the nucleosynthesis during the stellar life of massive stars and lead to stellar ejecta having various chemical compositions. We present low and zero-metallicity 20, 25 and 40 $M_{\odot}$ stellar models with various initial rotation rates and assumptions for the nuclear reactions rates. With increasing initial rotation, the yields of light (from $\sim$ C to Al) and trans-iron elements are boosted. The trans-iron elements (especially elements heavier than Ba) are significantly affected by the nuclear reaction uncertainties. The chemical composition of the observed CEMP (carbon-enhanced metal-poor) stars CS29528-028 and HE0336+0113 are consistent with the chemical composition of the material ejected by a fast rotating 40~$M_{\odot}$ model. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: Proceeding of the Nuclear Physics in Astrophysics IX conference, to appear in the Journal of Physics: Conference Series

The Formation of a 70 Msun Black Hole at High Metallicity

Authors: K. Belczynski, R. Hirschi, E. A. Kaiser, Jifeng Liu, J. Casares, Youjun Lu, R. O'Shaughnessy, A. Heger, S. Justham, R. Soria

Abstract: A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if s… ▽ More A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios. △ Less

Submitted 11 January, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Comments: ApJ: after 2nd review

The contribution from rotating massive stars to the enrichment in Sr and Ba of the Milky Way

Authors: Federico Rizzuti, Gabriele Cescutti, Francesca Matteucci, Alessandro Chieffi, Raphael Hirschi, Marco Limongi

Abstract: Most neutron capture elements have a double production by r- and s-processes, but the question of production sites is complex and still open. Recent studies show that including stellar rotation can have a deep impact on nucleosynthesis. We studied the evolution of Sr and Ba in the Milky Way. A chemical evolution model was employed to reproduce the Galactic enrichment. We tested two different nucle… ▽ More Most neutron capture elements have a double production by r- and s-processes, but the question of production sites is complex and still open. Recent studies show that including stellar rotation can have a deep impact on nucleosynthesis. We studied the evolution of Sr and Ba in the Milky Way. A chemical evolution model was employed to reproduce the Galactic enrichment. We tested two different nucleosynthesis prescriptions for s-process in massive stars, adopted from the Geneva group and the Rome group. Rotation was taken into account, studying the effects of stars without rotation or rotating with different velocities. We also tested different production sites for the r-process: magneto rotational driven supernovae and neutron star mergers. The evolution of the abundances of Sr and Ba is well reproduced. The comparison with the most recent observations shows that stellar rotation is a good assumption, but excessive velocities result in overproduction of these elements. In particular, the predicted evolution of the [Sr/Ba] ratio at low metallicity does not explain the data at best if rotation is not included. Adopting different rotational velocities for different stellar mass and metallicity better explains the observed trends. Despite the differences between the two sets of adopted stellar models, both show a better agreement with the data assuming an increase of rotational velocity toward low metallicity. Assuming different r-process sources does not alter this conclusion. △ Less

Submitted 10 September, 2019; originally announced September 2019.

Comments: 13 pages, 11 figures

Impact of Uncertainties in Astrophysical Reaction Rates on Nucleosynthesis in the $νp$ Process

Authors: T. Rauscher, N. Nishimura, G. Cescutti, R. Hirschi, A. St. J. Murphy, C. Fröhlich

Abstract: The $νp$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a larg… ▽ More The $νp$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a large reaction network. As the detailed astrophysical conditions of the $νp$ process still are unknown, a parameter study was performed, with 23 trajectories covering a large range of entropies and $Y_\mathrm{e}$. The resulting abundance uncertainties are given for each trajectory. The $νp$ process has been speculated to contribute to the light $p$ nuclides but it was not possible so far to reproduce the solar isotope ratios. It is found that it is possible to reproduce the solar $^{92}$Mo/$^{94}$Mo abundance ratio within nuclear uncertainties, even within a single trajectory. The solar values of the abundances in the Kr-Sr region relative to the Mo region, however, cannot be achieved within a single trajectory. They may still be obtained from a weighted superposition of different trajectories, though, depending on the actual conditions in the production site. For a stronger constraint of the required conditions, it would be necessary to reduce the uncertainties in the 3$α$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK. △ Less

Submitted 12 September, 2019; v1 submitted 7 September, 2019; originally announced September 2019.

Comments: 6 pages, 1 figure, 2 tables; to appear in JPS Conference Proceedings (OMEG15, 2019); summarizes and expands on arXiv:1907.13129 with an additional discussion of isotopic ratios of p-nuclides; v2: typos fixed

Journal ref: JPS Conf. Proc. 31 (2020) 011026

The s process in rotating low-mass AGB stars. Nucleosynthesis calculations in models matching asteroseismic constraints

Authors: J. W. den Hartogh, R. Hirschi, M. Lugaro, C. L. Doherty, U. Battino, F. Herwig, M. Pignatari, P. Eggenberger

Abstract: Aims. We investigate the s-process during the AGB phase of stellar models whose cores are enforced to rotate at rates consistent with asteroseismology observations of their progenitors and successors. Methods. We calculated new 2M$_{\odot}$, Z=0.01 models, rotating at 0, 125, and 250 km/s at the start of main sequence. An artificial, additional viscosity was added to enhance the transport of angul… ▽ More Aims. We investigate the s-process during the AGB phase of stellar models whose cores are enforced to rotate at rates consistent with asteroseismology observations of their progenitors and successors. Methods. We calculated new 2M$_{\odot}$, Z=0.01 models, rotating at 0, 125, and 250 km/s at the start of main sequence. An artificial, additional viscosity was added to enhance the transport of angular momentum in order to reduce the core rotation rates to be in agreement with asteroseismology observations. We compared rotation rates of our models with observed rotation rates during the MS up to the end of core He burning, and the white dwarf phase. Results. We present nucleosynthesis calculations for these rotating AGB models that were enforced to match the asteroseismic constraints on rotation rates of MS, RGB, He-burning, and WD stars. In particular, we calculated one model that matches the upper limit of observed rotation rates of core He-burning stars and we also included a model that rotates one order of magnitude faster than the upper limit of the observations. The s-process production in both of these models is comparable to that of non-rotating models. Conclusions. Slowing down the core rotation rate in stars to match the above mentioned asteroseismic constraints reduces the rotationally induced mixing processes to the point that they have no effect on the s-process nucleosynthesis. This result is independent of the initial rotation rate of the stellar evolution model. However, there are uncertainties remaining in the treatment of rotation in stellar evolution, which need to be reduced in order to confirm our conclusions, including the physical nature of our approach to reduce the core rotation rates of our models, and magnetic processes. △ Less

Submitted 24 August, 2019; originally announced August 2019.

Comments: accepted for publication in A&A

Journal ref: A&A 629, A123 (2019)

Uncertainties in $ν$p-process nucleosynthesis from Monte Carlo variation of reaction rates

Authors: N. Nishimura, T. Rauscher, R. Hirschi, G. Cescutti, A. St. J. Murphy, C. Fröhlich

Abstract: It has been suggested that a $ν$p process can occur when hot, dense, and proton-rich matter is expanding within a strong flux of anti-neutrinos. In such an environment, proton-rich nuclides can be produced in sequences of proton captures and (n,p) reactions, where the free neutrons are created in situ by $\overlineν_\mathrm{e}+\mathrm{p} \rightarrow \mathrm{n}+\mathrm{e}^+$ reactions. The detailed… ▽ More It has been suggested that a $ν$p process can occur when hot, dense, and proton-rich matter is expanding within a strong flux of anti-neutrinos. In such an environment, proton-rich nuclides can be produced in sequences of proton captures and (n,p) reactions, where the free neutrons are created in situ by $\overlineν_\mathrm{e}+\mathrm{p} \rightarrow \mathrm{n}+\mathrm{e}^+$ reactions. The detailed hydrodynamic evolution determines where the nucleosynthesis path turns off from N = Z line and how far up the nuclear chart it runs. In this work, the uncertainties on the final isotopic abundances stemming from uncertainties in the nuclear reaction rates were investigated in a large-scale Monte Carlo approach, simultaneously varying ten thousand reactions. A large range of model conditions was investigated because a definitive astrophysical site for the $ν$p process has not yet been identified. The present parameter study provides, for each model, identification of the key nuclear reactions dominating the uncertainty for a given nuclide abundance. As all rates appearing in the $ν$p process involve unstable nuclei, and thus only theoretical rates are available, the final abundance uncertainties are larger than those for nucleosynthesis processes closer to stability. Nevertheless, most uncertainties remain below a factor of three in trajectories with robust nucleosynthesis. More extreme conditions allow production of heavier nuclides but show larger uncertainties because of the accumulation of the uncertainties in many rates and because the termination of nucleosynthesis is not at equilibrium conditions. It is also found that the solar ratio of the abundances of ${}^{92}$Mo and ${}^{94}$Mo could be reproduced within uncertainties. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 19 pages, 12 figures, 8 tables, accepted for publication in MNRAS

Journal ref: MNRAS 489 1379 (2019)

Impact of Uncertainties in Nuclear Reaction Cross Sections on p-Nucleosynthesis in Thermonuclear Supernovae

Authors: T. Rauscher, N. Nishimura, G. Cescutti, R. Hirschi, A. St. J. Murphy, C. Travaglio

Abstract: The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucle… ▽ More The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucleosynthesis. Temperature-dependent rate uncertainties combining realistic experimental and theoretical uncertainties are used. From detailed statistical analyses uncertainties in the final abundances are derived as probability density distributions. Furthermore, based on rate and abundance correlations an automated procedure identifies the most important reactions in complex flow patterns from superposition of many zones or tracers. The method so far was already applied to a number of nucleosynthesis processes. Here we focus on the production of p-nuclei in white dwarfs exploding as thermonuclear (type Ia) supernovae. We find generally small uncertainties in the final abundances despite of the dominance of theoretical nuclear uncertainties. A separate analysis of low- and high-density regions indicates that the total uncertainties are dominated by the high-density regions. △ Less

Submitted 16 August, 2019; v1 submitted 22 July, 2019; originally announced July 2019.

Comments: 6 pages, 3 figures; invited talk at 13th Int. Conf. on Nucleus-Nucleus Collisions, Saitama, Japan, 2018. To be published in JPS Conf. Proceedinngs; summarizes results from arXiv:1807.11475; v2: updated reference and receipt date; as accepted for publication

NuGrid stellar data set - III. Updated low-mass AGB models and s-process nucleosynthesis with metallicities Z=0.01, Z=0.02 and Z=0.03

Authors: Umberto Battino, Ashley Tattersall, Claudia Lederer-Woods, Falk Herwig, Pavel Denissenkov, Raphael Hirschi, Reto Trappitsch, Jacqueline W. den Hartogh, Marco Pignatari

Abstract: The production of the neutron-capture isotopes beyond iron that we observe today in the solar system is the result of the combined contribution of the r-process, the s- process and possibly the i-process. Low-mass AGB (2 < M/Msun < 3) and massive (M >10 Msun ) stars have been identified as the sites of the s-process. In this work we consider the evolution and nucleosynthesis of low-mass AGB stars.… ▽ More The production of the neutron-capture isotopes beyond iron that we observe today in the solar system is the result of the combined contribution of the r-process, the s- process and possibly the i-process. Low-mass AGB (2 < M/Msun < 3) and massive (M >10 Msun ) stars have been identified as the sites of the s-process. In this work we consider the evolution and nucleosynthesis of low-mass AGB stars. We provide an update of the NuGrid Set models, adopting the same general physics assumptions but using an updated convective-boundary mixing model accounting for the contribution from internal gravity waves. The combined data set includes the initial masses Mzams/Msun = 2, 3 for Z = 0.03, 0.02, 0.01. These models are computed with the MESA stellar code and the evolution is followed up to the end of the AGB phase. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. The convective boundary mixing model leads to the formation of a 13C-pocket three times wider compared to the one obtained in the previous set of models, bringing the simulation results now in closer agreement with observations. We also discuss the potential impact of other processes inducing mixing, like rotation, adopting parametric models compatible with theory and observations. Complete yield data tables, derived data products and online analytic data access are provided. △ Less

Submitted 5 June, 2019; originally announced June 2019.

Comments: Submitted to MNRAS

Grids of stellar models with rotation IV. Models from 1.7 to 120 Msun at a metallicity Z = 0.0004

Authors: Jose H. Groh, Sylvia Ekström, Cyril Georgy, Georges Meynet, Arthur Choplin, Patrick Eggenberger, Raphael Hirschi, Andre Maeder, Laura J. Murphy, Ioana Boian, Eoin J. Farrell

Abstract: The effects of rotation on stellar evolution are particularly important at low metallicity, when mass loss by stellar winds diminishes and the surface enrichment due to rotational mixing becomes relatively more pronounced than at high metallicities. Here we investigate the impact of rotation and metallicity on stellar evolution. Using a similar physics as in our previous large grids of models at Z… ▽ More The effects of rotation on stellar evolution are particularly important at low metallicity, when mass loss by stellar winds diminishes and the surface enrichment due to rotational mixing becomes relatively more pronounced than at high metallicities. Here we investigate the impact of rotation and metallicity on stellar evolution. Using a similar physics as in our previous large grids of models at Z=0.002 and Z=0.014, we compute stellar evolution models with the Geneva code for rotating and nonrotating stars with initial masses (Mini) between 1.7 and 120 Msun and Z=0.0004 (1/35 solar). This is comparable to the metallicities of the most metal poor galaxies observed so far, such as I Zw 18. Concerning massive stars, both rotating and nonrotating models spend most of their core-helium burning phase with an effective temperature higher than 8000 K. Stars become red supergiants only at the end of their lifetimes, and few RSGs are expected. Our models predict very few to no classical Wolf-Rayet stars as a results of weak stellar winds at low metallicity. The most massive stars end their lifetimes as luminous blue supergiants or luminous blue variables, a feature that is not predicted by models with higher metallicities. Interestingly, due to the behavior of the intermediate convective zone, the mass domain of stars producing pair-instability supernovae is smaller at Z=0.0004 than at Z=0.002. We find that during the main sequence phase, the ratio between nitrogen and carbon abundances (N/C) remains unchanged for nonrotating models. However, N/C increases by factors of 10-20 in rotating models at the end of the MS. Cepheids coming from stars with Mini > 4-6 Msun are beyond the core helium burning phase and spend little time in the instability strip. Since they would evolve towards cooler effective temperatures, these Cepheids should show an increase of the pulsation period as a function of age. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: 19 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 627, A24 (2019)

Constraining transport of angular momentum in stars: Combining asteroseismic observations of core helium burning stars and white dwarfs

Authors: J. W. den Hartogh, P. Eggenberger, R. Hirschi

Abstract: Transport of angular momentum has been a challenging topic within the stellar evolution community, even more since the recent asteroseismic surveys. All published studies on rotation using asteroseismic observations show a discrepancy between the observed and calculated rotation rates, indicating there is an undetermined process of angular momentum transport active in these stars. We aim to constr… ▽ More Transport of angular momentum has been a challenging topic within the stellar evolution community, even more since the recent asteroseismic surveys. All published studies on rotation using asteroseismic observations show a discrepancy between the observed and calculated rotation rates, indicating there is an undetermined process of angular momentum transport active in these stars. We aim to constrain the efficiency of this process by investigating rotation rates of 2.5 M$_{\odot}$ stars. First, we investigated whether the Tayler-Spruit dynamo could be responsible for the extra transport of angular momentum for stars with an initial mass of 2.5 M$_{\odot}$. Then, by computing rotating models including a constant additional artificial viscosity, we determined the efficiency of the missing process of angular momentum transport by comparing the models to the asteroseismic observations of core helium burning stars. Parameter studies were performed to investigate the effect of the stellar evolution code used, initial mass, and evolutionary stage. We evolved our models into the white dwarf phase, and provide a comparison to white dwarf rotation rates. The Tayler-Spruit dynamo is unable to provide enough transport of angular momentum to reach the observed values of the core helium burning stars investigated in this paper. We find that a value for the additional artificial viscosity $ν_{\rm{add}}$ around 10$^7$ cm$^2$ s$^{-1}$ provides enough transport of angular momentum. However, the rotational period of these models is too high in the white dwarf phase to match the white dwarf observations. From this comparison we infer that the efficiency of the missing process must decrease during the core helium burning phase. When excluding the $ν_{\rm{add}}$ during core helium burning phase, we can match the rotational periods of both the core helium burning stars and white dwarfs. △ Less

Submitted 12 February, 2019; originally announced February 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 622, A187 (2019)

Catching Element Formation In The Act

Authors: Chris L. Fryer, Frank Timmes, Aimee L. Hungerford, Aaron Couture, Fred Adams, Wako Aoki, Almudena Arcones, David Arnett, Katie Auchettl, Melina Avila, Carles Badenes, Eddie Baron, Andreas Bauswein, John Beacom, Jeff Blackmon, Stephane Blondin, Peter Bloser, Steve Boggs, Alan Boss, Terri Brandt, Eduardo Bravo, Ed Brown, Peter Brown, Steve Bruenn. Carl Budtz-Jorgensen, Eric Burns , et al. (194 additional authors not shown)

Abstract: Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray… ▽ More Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 14 pages including 3 figures

Report number: LA-UR-18-29748

Dependence of Convective Boundary Mixing on Boundary Properties and Turbulence Strength

Authors: A. Cristini, R. Hirschi, C. Meakin, D. Arnett, C. Georgy, I. Walkington

Abstract: Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with th… ▽ More Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the PROMPI code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier-Stokes framework. Both the vertical RMS velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number. Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

3D Simulations and MLT: II. Onsager's Ideal Turbulence

Authors: W. David Arnett, Raphael Hirschi, Simon W. Campbell, Miroslav Mocák, Cyril Georgy, Casey Meakin, Andrea Cristini, Laura J. A. Scott, Etienne A. Kaiser, Maxime Viallet

Abstract: We simulate stellar convection at high Reynolds number (Re$\lesssim$7000) with causal time stepping but no explicit viscosity. We use the 3D Euler equations with shock capturing (Colella & Woodward 1984). Anomalous dissipation of turbulent kinetic energy occurs as an emergent feature of advection ("Onsager damping"), caused by the moderate shocks which terminate the turbulent kinetic energy spectr… ▽ More We simulate stellar convection at high Reynolds number (Re$\lesssim$7000) with causal time stepping but no explicit viscosity. We use the 3D Euler equations with shock capturing (Colella & Woodward 1984). Anomalous dissipation of turbulent kinetic energy occurs as an emergent feature of advection ("Onsager damping"), caused by the moderate shocks which terminate the turbulent kinetic energy spectrum; see also (Perry 2021). In strongly stratified stellar convection the asymptotic limit for the global damping length of turbulent kinetic energy is $\ell_d \sim \langle u^3 \rangle /\langle ε\rangle$. This "dissipative anomaly" (Onsager 1949) fixes the value of the "mixing length parameter", $α= \ell_{\rm MLT}/H_P =\overline{\langleΓ_1\rangle}$, which is $\sim\, 5/3$ for complete ionization. The estimate is numerically robust, agrees to within 10% with estimates from stellar evolution with constant $α$. For weak stratification $\ell_d$ shrinks to the depth of a thin convective region. Our flows are filamentary, produce surfaces of separation at boundary layers, resolve the energy-containing eddies, and develop a turbulent cascade down to the grid scale which agrees with the $4096^3$ direct numerical simulation of Kaneda (2003). The cascade converges quickly, and satisfies a power-law velocity spectrum similar to Kolmogorov (1941). Our flows exhibit intermittency, anisotropy, and interactions between coherent structures, features missing from K41 theory. We derive a dissipation rate from Reynolds stresses which agrees with (i) our flows, (ii) experiment (Warhaft 2002), and (iii) high Re simulations of the Navier-Stokes equations (Iyer, et al. 2018). △ Less

Submitted 23 April, 2021; v1 submitted 10 October, 2018; originally announced October 2018.

Comments: 31 pages, 6 figures. The time evolution and the fly-through movies may be found at \url{http://www.astro.keele.ac.uk/shyne/321D/convection-and-convective-boundary-mixing/visualisations}

3D Simulations and MLT: I. Renzini's Critique

Authors: W. David Arnett, Casey Meakin, Raphael Hirschi, Andrea Cristini, Cyril Georgy, Simon Campbell, Laura Scott, Etienne Kaiser

Abstract: Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES)… ▽ More Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES) automatically connect large scale turbulence to a Kolmogorov cascade below the grid scale, allowing turbulent boundary layers to remove singularities that appear in the theory. Interactions between coherent structures give multi-modal behavior, driving intermittency and fluctuations. Reynolds averaging (RA) allows us to abstract the essential features of this dynamical behavior of boundaries which are appropriate to stellar evolution, and consider how they relate static boundary conditions (Richardson, Schwarzschild or Ledoux). We clarify several questions concerning when and why MLT works, and does not work, using both analytical theory and 3D high resolution numerical simulations. The composition gradients and boundary layer structure which are produced by our simulations suggest a self-consistent approach to boundary layers, removing the need for ad hoc procedures for 'convective overshooting' and `semi-convection'. In a companion paper we quantify the adequacy of our numerical resolution, determine of the length scale of dissipation (the `mixing length') without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and extend MLT to deal with strong stratification. △ Less

Submitted 10 October, 2018; originally announced October 2018.

The s-process nucleosynthesis in low mass stars: impact of the uncertainties in the nuclear physics determined by Monte Carlo variations

Authors: G. Cescutti, R. Hirschi, N. Nishimura, T. Rauscher, J. den Hartogh, A. St. J. Murphy, S. Cristallo

Abstract: We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to eva… ▽ More We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to evaluate the reactions that impact more significantly the final abundances. We found that beta-decay rate uncertainties affect typically nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are approximately 50%. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Comments: 3 pages, 1 figure; to appear in Springer Proceedings of Physics as proceedings of the XV Int. Conf. "Nuclei in the Cosmos", L'Aquila, Italy, June 2018. This conference contribution presents results published in MNRAS 478 (2018) 4101

Relative Importance of Convective Uncertainties

Authors: Etienne A. Kaiser, Raphael Hirschi, W. David Arnett, Andrea Cristini, Cyril Georgy, Laura J. A. Scott

Abstract: Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of… ▽ More Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of extra mixing, and their impact on the early evolution of massive stars. With increasing convective boundary mixing (CBM), we find a convergence of the two different boundary locations, a decreasing blue to red super giant ratio and a reduced importance of semiconvection. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 3 pages, 1 figure; to appear in Springer Proceedings in Physics (Proc. of Itl. Conf. "Nuclei in the Cosmos XV", LNGS Assergi, Italy, June 2018)

Uncertainties in the production of p nuclides in SN Ia determined by Monte Carlo variations

Authors: T. Rauscher, N. Nishimura, R. Hirschi, A. St. J. Murphy, G. Cescutti, C. Travaglio

Abstract: Several thousand tracers from a 2D model of a thermonuclear supernova were used in a Monte Carlo post-processing approach to determine p-nuclide abundance uncertainties originating from nuclear physics uncertainties in the reaction rates. Several thousand tracers from a 2D model of a thermonuclear supernova were used in a Monte Carlo post-processing approach to determine p-nuclide abundance uncertainties originating from nuclear physics uncertainties in the reaction rates. △ Less

Submitted 28 July, 2018; originally announced July 2018.

Comments: 3 pages, 1 figure; to appear in Springer Proceedings of Physics as proceedings of the XV Int. Conf. "Nuclei in the Cosmos", L'Aquila, Italy, June 2018. This conference contribution presents results published in MNRAS 474 (2018) 3133

Non standard s-process in massive rotating stars. Yields of $10-150$ $M_{\odot}$ models at $Z=10^{-3}$

Authors: Arthur Choplin, Raphael Hirschi, Georges Meynet, Sylvia Ekstroöm, Cristina Chiappini, Alison Laird

Abstract: Recent studies show that rotation significantly affects the s-process in massive stars. We provide tables of yields for non-rotating and rotating massive stars between 10 and 150 $M_{\odot}$ at $Z=10^{-3}$ ([Fe/H] $=-1.8$). Tables for different mass cuts are provided. The complete s-process is followed during the whole evolution with a network of 737 isotopes, from Hydrogen to Polonium. A grid of… ▽ More Recent studies show that rotation significantly affects the s-process in massive stars. We provide tables of yields for non-rotating and rotating massive stars between 10 and 150 $M_{\odot}$ at $Z=10^{-3}$ ([Fe/H] $=-1.8$). Tables for different mass cuts are provided. The complete s-process is followed during the whole evolution with a network of 737 isotopes, from Hydrogen to Polonium. A grid of stellar models with initial masses of 10, 15, 20, 25, 40, 60, 85, 120 and 150 $M_{\odot}$ and with an initial rotation rate of both 0 or 40$~\%$ of the critical velocity was computed. Three extra models were computed in order to investigate the effect of faster rotation (70$~\%$ of the critical velocity) and of a lower $^{17}$O($α,γ$) reaction rate. At the considered metallicity, rotation has a strong impact on the production of s-elements for initial masses between 20 and 60 $M_{\odot}$. In this range, the first s-process peak is boosted by $2-3$ dex if rotation is included. Above 60 $M_{\odot}$, s-element yields of rotating and non-rotating models are similar. Increasing the initial rotation from 40$~\%$ to 70$~\%$ of the critical velocity enhances the production of $40 \lesssim Z \lesssim 60$ elements by $\sim 0.5-1$ dex. Adopting a reasonably lower $^{17}$O($α,γ$) rate in the fast rotating model (70$~\%$ of the critical velocity) boosts again the yields of s-elements with $55 \lesssim Z \lesssim 82$ by about 1 dex. In particular, a modest amount of Pb is produced. Together with s-elements, some light elements (particularly fluorine) are strongly overproduced in rotating models. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: 16 pages, 14 figures, 4 tables. Accepted for publication in A&A

Journal ref: A&A 618, A133 (2018)

Uncertainties in s-process nucleosynthesis in low mass stars determined from Monte Carlo variations

Authors: G. Cescutti, R. Hirschi, N. Nishimura, J. W. den Hartogh, T. Rauscher, A. St. J. Murphy, S. Cristallo

Abstract: The main s-process taking place in low mass stars produces about half of the elements heavier than iron. It is therefore very important to determine the importance and impact of nuclear physics uncertainties on this process. We have performed extensive nuclear reaction network calculations using individual and temperature-dependent uncertainties for reactions involving elements heavier than iron,… ▽ More The main s-process taking place in low mass stars produces about half of the elements heavier than iron. It is therefore very important to determine the importance and impact of nuclear physics uncertainties on this process. We have performed extensive nuclear reaction network calculations using individual and temperature-dependent uncertainties for reactions involving elements heavier than iron, within a Monte Carlo framework. Using this technique, we determined the uncertainty in the main s-process abundance predictions due to nuclear uncertainties link to weak interactions and neutron captures on elements heavier than iron. We also identified the key nuclear reactions dominating these uncertainties. We found that $β$-decay rate uncertainties affect only a few nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are in general small (less than 50%). Three key reactions, nevertheless, stand out because they significantly affect the uncertainties of a large number of nuclides. These are $^{56}$Fe(n,$γ$), $^{64}$Ni(n,$γ$), and $^{138}$Ba(n,$γ$). We discuss the prospect of reducing uncertainties in the key reactions identified in this study with future experiments. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 28 pages, 17 figures, 8 tables, accepted for publication in MNRAS

Journal ref: Mon. Not. Roy. Astron. Soc. 478 (2018) 4101

OGLE14-073 - a promising pair-instability supernova candidate

Authors: Alexandra Kozyreva, Markus Kromer, Ulrich M. Noebauer, Raphael Hirschi

Abstract: The recently discovered bright type II supernova OGLE14-073 evolved very slowly. The light curve rose to maximum for 90 days from discovery and then declined at a rate compatible with the radioactive decay of 56Co. In this study, we show that a pair-instability supernova is a plausible mechanism for this event. We calculate explosion models and light curves with the radiation hydrodynamics code ST… ▽ More The recently discovered bright type II supernova OGLE14-073 evolved very slowly. The light curve rose to maximum for 90 days from discovery and then declined at a rate compatible with the radioactive decay of 56Co. In this study, we show that a pair-instability supernova is a plausible mechanism for this event. We calculate explosion models and light curves with the radiation hydrodynamics code STELLA starting from two M(ZAMS) = 150 solar masses, Z=0.001 progenitors. We obtain satisfactory fits to OGLE14-073 broadband light curves by including additional 56Ni in the centre of the models and mixing hydrogen down into the inner layers of the ejecta to a radial mass coordinate of 10 solar masses. The extra 56Ni required points to a slightly more massive progenitor star. The mixing of hydrogen could be due to large scale mixing during the explosion. We also present synthetic spectra for our models simulated with the Monte Carlo radiative transfer code ARTIS. The synthetic spectra reproduce the main features of the observed spectra of OGLE14-073. We conclude that OGLE14-073 is one of the most promising candidates for a pair-instability explosion. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 10 pages, 8 figures. Accepted for publication in MNRAS

Progenitors of Core-Collapse Supernovae

Authors: Raphael Hirschi, David Arnett, Andrea Cristini, Cyril Georgy, Casey Meakin, Ian Walkington

Abstract: Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss som… ▽ More Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss some of the major uncertainties in the modelling of massive stars, with a particular emphasis on the treatment of convection in 1D stellar evolution codes. Finally, we present new 3D hydrodynamic simulations of convection in carbon burning and list key points to take from 3D hydrodynamic studies for the development of new prescriptions for convective boundary mixing in 1D stellar evolution codes. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 10 pages, 5 figures, "SN 1987A, 30 years later", Proceedings IAU Symposium No. 331, 2017; A. Marcowith, G. Dubner, A. Ray, A. Bykov, & M. Renaud, eds

Sensitivity to neutron captures and beta-decays of the enhanced s-process in rotating massive stars at low metallicities

Authors: N. Nishimura, R. Hirschi, T. Rauscher

Abstract: The s-process in massive stars, producing nuclei up to $A\approx 90$, has a different behaviour at low metallicity if stellar rotation is significant. This enhanced s-process is distinct from the s-process in massive stars around solar metallicity, and details of the nucleosynthesis are poorly known. We investigated nuclear physics uncertainties in the enhanced s-process in metal-poor stars within… ▽ More The s-process in massive stars, producing nuclei up to $A\approx 90$, has a different behaviour at low metallicity if stellar rotation is significant. This enhanced s-process is distinct from the s-process in massive stars around solar metallicity, and details of the nucleosynthesis are poorly known. We investigated nuclear physics uncertainties in the enhanced s-process in metal-poor stars within a Monte-Carlo framework. We applied temperature-dependent uncertainties of reaction rates, distinguishing contributions from the ground state and from excited states. We found that the final abundance of several isotopes shows uncertainties larger than a factor of 2, mostly due to the neutron capture uncertainties. A few nuclei around branching points are affected by uncertainties in the $β$-decay. △ Less

Submitted 15 February, 2018; originally announced February 2018.

Comments: 3 pages, 3 figures, published in the Proceedings of "Nuclear Physics in Astrophysics Conference (NPA VII)"; see arXiv:1701.00489, for the completed results

Journal ref: 2018 J. Phys.: Conf. Ser. 940 012051

Impacts of nuclear-physics uncertainties in the s-process determined by Monte-Carlo variations

Authors: N. Nishimura, G. Cescutti, R. Hirschi, T. Rauscher, J. den Hartogh, A. St. J. Murphy

Abstract: The s-process, a production mechanism based on slow-neutron capture during stellar evolution, is the origin of about half the elements heavier than iron. Abundance predictions for s-process nucleosynthesis depend strongly on the relevant neutron-capture and $β$-decay rates, as well as on the details of the stellar model being considered. Here, we have used a Monte-Carlo approach to evaluate the nu… ▽ More The s-process, a production mechanism based on slow-neutron capture during stellar evolution, is the origin of about half the elements heavier than iron. Abundance predictions for s-process nucleosynthesis depend strongly on the relevant neutron-capture and $β$-decay rates, as well as on the details of the stellar model being considered. Here, we have used a Monte-Carlo approach to evaluate the nuclear uncertainty in s-process nucleosynthesis. We considered the helium burning of massive stars for the weak s-process and low-mass asymptotic-giant-branch stars for the main s-process. Our calculations include a realistic and general prescription for the temperature dependent uncertainty for the reaction cross sections. We find that the adopted uncertainty for (${\rm n},γ$) rates, tens of per cent on average, effects the production of s-process nuclei along the line of $β$-stability, and that the uncertainties in $β$-decay from excited state contributions, has the strongest impact on branching points. △ Less

Submitted 15 February, 2018; originally announced February 2018.

Comments: 6 pages, 4 figures, 2 tables, the Proceedings of "the 2017 Symposium on Nuclear Data"; a supplementary article of arXiv:1701.00489

Journal ref: JAEA-Conf 2018-001

The First 3D Simulations of Carbon Burning in a Massive Star

Authors: Andrea Cristini, Casey Meakin, Raphael Hirschi, David Arnett, Cyril Georgy, Maxime Viallet

Abstract: We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs… ▽ More We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs via turbulent entrainment as in the case of oxygen burning. The expansion of the boundary into the surrounding stable region and the entrainment rate are smaller at the bottom boundary because it is stiffer than the upper boundary. The results of this and similar studies call for improved CBM prescriptions in 1D stellar evolution models. △ Less

Submitted 19 December, 2017; originally announced December 2017.

Comments: 5 pages, 3 figures. Published in IAUS 329 on 28/07/17

Journal ref: Proceedings of the International Astronomical Union, 12(S329), 237-241 (2017)