Showing 1–50 of 159 results for author: Thielemann, F

Peculiarities of the chemical enrichment of metal-poor Stars in the Milky Way Galaxy

Authors: T. Mishenina, M. Pignatari, I. Usenko, C. Soubiran, F. -K. Thielemann, A. Yu. Kniazev, S. A. Korotin, T. Gorbaneva

Abstract: The oldest stars in the Milky Way are metal-poor with [Fe/H] < -- 1.0, displaying peculiar elemental abundances compared to solar values. The relative variations in the chemical compositions among stars is also increasing with decreasing stellar metallicity, allowing for the pure signature of unique nucleosynthesis processes to be revealed. In this work, we report the atmospheric parameters, main… ▽ More The oldest stars in the Milky Way are metal-poor with [Fe/H] < -- 1.0, displaying peculiar elemental abundances compared to solar values. The relative variations in the chemical compositions among stars is also increasing with decreasing stellar metallicity, allowing for the pure signature of unique nucleosynthesis processes to be revealed. In this work, we report the atmospheric parameters, main dynamic properties, and the abundances of four metal-poor stars: HE 1523--0901, HD 6268, HD 121135, and HD 195636 (--1.5 > [Fe/H] >--3.0). The abundances were derived from spectra obtained with the HRS echelle spectrograph at the SALT, using both LTE and NLTE approaches. Based on their kinematical properties, we show that HE 1523--0901 and HD 195636 are halo stars with typical high velocities. In particular, HD 121135 displays a peculiar kinematical behaviour, making it unclear whether it is a halo or an accreted star. Furthermore, HD 6268 is possibly a rare prototype of very metal-poor thick disk stars. The abundances derived for our stars are compared with theoretical stellar models and with other stars with similar metallicity values from the literature. HD 121135 is Al-poor and Sc-poor, compared to stars observed in the same metallicity range (--1.62 > [Fe/H] >--1.12). The most metal-poor stars in our sample, HE 1523 -- 0901, HD 6268, and HD 195636, exhibit anomalies that are better explained by supernova models from fast-rotating stellar progenitors for elements up to the Fe group. Compared to other stars in the same metallicity range, their common biggest anomaly is represented by the low Sc abundances. If we consider the elements beyond Zn, HE 1523--0901 can be classified as an r-II star, HD 6268 as an r-I candidate, and HD 195636 and HD 121135 exhibiting a borderline r-process enrichment between limited-r and r-I star. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Comments: 27 pages, 15 figures, accepted in A&A 23.04.2024

Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes

Authors: Domingo García-Senz, Rubén M. Cabezón, Moritz Reichert, Axel S. Lechuga, José A. Escartín, Athanasios Psaltis, Almudena Arcones, Friedrich-Karl Thielemann

Abstract: We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free… ▽ More We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed with the extended network WinNet (> 2000 isotopes) in scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer to WinNet than typical reduced alpha networks. △ Less

Submitted 11 July, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: 19 pages, 21 Figures and 5 Tables. Accepted for publication in A&A

Nuclear Reactions in Evolving Stars

Authors: Friedrich-Karl Thielemann, Thomas Rauscher

Abstract: This chapter will go through the important nuclear reactions in stellar evolution and explosions, passing through the individual stellar burning stages and also explosive burning conditions. To follow the changes in the composition of nuclear abundances requires the knowledge of the relevant nuclear reaction rates. For light nuclei (entering in early stellar burning stages) the resonance density i… ▽ More This chapter will go through the important nuclear reactions in stellar evolution and explosions, passing through the individual stellar burning stages and also explosive burning conditions. To follow the changes in the composition of nuclear abundances requires the knowledge of the relevant nuclear reaction rates. For light nuclei (entering in early stellar burning stages) the resonance density is generally quite low and the reactions are determined by individual resonances, which are best obtained from experiments. For intermediate mass and heavy nuclei the level density is typically sufficient to apply statistical model approaches. For this reason, while we discuss all burning stages and explosive burning, focusing on the reactions of importance, we will for light nuclei refer to the chapters by M. Wiescher, deBoer & Reifarth (Experimental Nuclear Astrophysics) and P. Descouvement (Theoretical Studies of Low-Energy Nuclear Reactions), which display many examples, experimental methods utilized, and theoretical approaches how to predict nuclear reaction rates for light nuclei. For nuclei with sufficiently high level densities we discuss statistical model methods used in present predictions of nuclear reaction cross sections and thermonuclear rates across the nuclear chart, including also the application to nuclei far from stability and fission modes. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: invited contribution on "Nuclear Reactions in Evolving Stars" for the Handbook of Nuclear Physics, 55 pages, 19 figures,

Journal ref: Handbook of Nuclear Physics 2023

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.

Enrichment of the Galactic disc with neutron-capture elements: Gd, Dy, and Th

Authors: T. Mishenina, M. Pignatari, T. Gorbaneva, B. Cote, A. Yague Lopez, F. -K. Thielemann, C. Soubiran

Abstract: The study of the origin of heavy elements is one of the main goals of nuclear astrophysics. In this paper, we present new observational data for the heavy $r$-process elements gadolinium (Gd, Z=64), dysprosium (Dy, Z=66) and thorium (Th, Z=90) in a sample of 276 Galactic disc stars ( --1.0$<$[Fe/H]$<$+0.3). The stellar spectra have a high resolution of 42,000 and 75,000, and the signal-to-noise ra… ▽ More The study of the origin of heavy elements is one of the main goals of nuclear astrophysics. In this paper, we present new observational data for the heavy $r$-process elements gadolinium (Gd, Z=64), dysprosium (Dy, Z=66) and thorium (Th, Z=90) in a sample of 276 Galactic disc stars ( --1.0$<$[Fe/H]$<$+0.3). The stellar spectra have a high resolution of 42,000 and 75,000, and the signal-to-noise ratio higher than 100. The LTE abundances of Gd, Dy and Th have been determined by comparing the observed and synthetic spectra for three Gd lines (149 stars), four Dy lines (152 stars) and the Th line at 4019.13 A (170 stars). For about 70% of the stars in our sample Gd and Dy are measured for the first time, and Th for 95% of the stars. Typical errors vary from 0.07 to 0.16 dex. This paper provides the first extended set of Th observations in the Milky Way disc. Together with europium (Eu, Z = 63) data from our previous studies, we have compared these new observations with nucleosynthesis predictions and Galactic Chemical Evolution simulations. We confirm that [Gd/Fe] and [Dy/Fe] show the same behavior of Eu. We study with GCE simulations the evolution of [Th/Fe] in comparison with [Eu/Fe], showing that unlike Eu either the Th production is metallicity dependent in case of a unique source of the r-process in the Galaxy, or the frequency of the Th-rich r-process source is decreasing with the increasing of [Fe/H]. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 18 pages, 21 figures, accepted in MNRAS, stac2361

Magnetorotational supernovae: A nucleosynthetic analysis of sophisticated 3D models

Authors: Moritz Reichert, Martin Obergaulinger, Miguel-Angel Aloy, Michael Gabler, Almudena Arcones, Friedrich-Karl Thielemann

Abstract: Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements:… ▽ More Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements: i) a prompt ejection of matter right after core bounce, ii) neutron-rich matter that is ejected at late times due to a reconfiguration of the protoneutronstar shape, iii) small amount of mass ejected with high entropies in the center of the jet. We investigate total ejecta yields, including the ones of unstable nuclei such as $^{26}$Al, $^{44}$Ti, $^{56}$Ni, and $^{60}$Fe. The obtained $^{56}$Ni masses vary between $0.01 - 1\,\mathrm{M_\odot}$. The latter maximum is compatible with hypernova observations. Furthermore, all of our models synthesize Zn masses in agreement with observations of old metal-poor stars. We calculate simplified light curves to investigate whether our models can be candidates for superluminous supernovae. The peak luminosities obtained from taking into account only nuclear heating reach up to a few $\sim 10^{43} \,\mathrm{erg\,s^{-1}}$. Under certain conditions, we find a significant impact of the $^{66}$Ni decay chain that can raise the peak luminosity up to $\sim 38\%$ compared to models including only the $^{56}$Ni decay chain. This work reinforces the theoretical evidence on the critical role of magnetorotational supernovae to understand the occurrence of hypernovae, superluminous supernovae, and the synthesis of heavy elements. △ Less

Submitted 16 November, 2022; v1 submitted 23 June, 2022; originally announced June 2022.

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

$r$-process nucleosynthesis from compact binary mergers

Authors: Albino Perego, Friedrich-Karl Thielemann, Gabriele Cescutti

Abstract: The merger of two neutron stars or of a neutron star and a black hole often result in the ejection of a few percents of a solar mass of matter expanding at high speed in space. Being matter coming from the violent disruption of a neutron star, these ejecta are initially very dense, hot and extremely rich in neutrons. The few available protons form heavy nuclei ("seeds") that absorb the more abunda… ▽ More The merger of two neutron stars or of a neutron star and a black hole often result in the ejection of a few percents of a solar mass of matter expanding at high speed in space. Being matter coming from the violent disruption of a neutron star, these ejecta are initially very dense, hot and extremely rich in neutrons. The few available protons form heavy nuclei ("seeds") that absorb the more abundant free neutrons, increasing their size. The neutron density is so high that a substantial number of neutron captures occur before the resulting unstable nuclei can decay toward more stable configurations, converting neutrons into protons. Depending mostly on the initial neutron richness, this mechanism leads to the formation of up to half of the heavy elements that we observe in nature and it is called rapid neutron capture process ("$r$-process"). The prediction of the precise composition of the ejecta requires a detailed knowledge of the properties of very exotic nuclei, that have never been produced in a laboratory. Despite having long been a speculative scenario, nowadays several observational evidences point to compact binary mergers as one of the major sites where heavy elements are formed in the Universe. The most striking one was the detection of a kilonova following the merger of a neutron star binary: the light emitted by this astronomical transient is indeed powered by the radioactive decay of freshly synthesized neutron-rich nuclei and testifies the actual nature of compact binary mergers as cosmic forges. △ Less

Submitted 17 December, 2021; v1 submitted 19 September, 2021; originally announced September 2021.

Comments: 56 pages; 14 figures, Chapter of "Handbook of Gravitational Wave Astronomy", published by Springer

Correlations of r-Process Elements in Very Metal-Poor Stars as Clues to their Nucleosynthesis Sites

Authors: Khalil Farouqi, Friedrich-Karl Thielemann, Stephan Rosswog, Karl-Ludwig Kratz

Abstract: Various nucleosynthesis studies have pointed out that the r-process elements in very metal-poor (VMP) halo stars might have different origins. By means of familiar concepts from statistics (correlations, cluster analysis, rank tests of elemental abundances), we look for causally correlated elemental abundance patterns and attempt to link them to astrophysical events. Some of these events produce t… ▽ More Various nucleosynthesis studies have pointed out that the r-process elements in very metal-poor (VMP) halo stars might have different origins. By means of familiar concepts from statistics (correlations, cluster analysis, rank tests of elemental abundances), we look for causally correlated elemental abundance patterns and attempt to link them to astrophysical events. Some of these events produce the r-process elements jointly with iron, while others do not have any significant iron contribution. In the early stage of our Galaxy, at least three r-process nucleosynthesis sites have been active. The first two produce and eject iron and the majority of the lighter r-process elements. We assign them to two different types of core-collapse events, not identical to regular core-collapse supernovae (CCSNe), which produce only light trans-Fe elements. The third category is characterized by a strong r-process and responsible for the major fraction of the heavy main r-process elements without a significant co-production of Fe. It does not appear to be connected to CCSNe, in fact the Fe found in the related r-process enriched stars must come from previously occurring CCSNe. The existence of actinide boost stars indicates a further division among strong r-process sites. We assign these two strong r-process sites to neutron star mergers without fast black hole formation and to events where the ejecta are dominated by black hole accretion disk outflows. Indications from the lowest-metallicity stars hint at a connection to massive single stars (collapsars) forming black holes in the early Galaxy. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: 40 pages, 25 figures, 10 tables

Journal ref: A&A 663, A70 (2022)

Stellar Mass Black Hole Formation and Multi-messenger Signals from Three Dimensional Rotating Core-Collapse Supernova Simulations

Authors: Kuo-Chuan Pan, Matthias Liebendörfer, Sean Couch, Friedrich-Karl Thielemann

Abstract: We present self-consistent, 3D core-collapse supernova simulations of a 40 Msun progenitor model using the isotropic diffusion source approximation for neutrino transport and an effective general relativistic potential up to $\sim0.9$~s~postbounce. We consider three different rotational speeds with initial angular velocities of $Ω_0=0$,~0.5, and~1~rad~s$^{-1}$ and investigate the impact of rotatio… ▽ More We present self-consistent, 3D core-collapse supernova simulations of a 40 Msun progenitor model using the isotropic diffusion source approximation for neutrino transport and an effective general relativistic potential up to $\sim0.9$~s~postbounce. We consider three different rotational speeds with initial angular velocities of $Ω_0=0$,~0.5, and~1~rad~s$^{-1}$ and investigate the impact of rotation on shock dynamics, black hole formation, and gravitational wave signals. The rapidly-rotating model undergoes an early explosion at $\sim 250$~ms postbounce and shows signs of the low $T/|W|$ instability. We do not find black hole formation in this model within $\sim 460$~ms postbounce. In contrast, we find black hole formation at 776~ms~postbounce and 936~ms~postbounce for the non-rotating and slowly-rotating models, respectively. The slowly-rotating model explodes at $\sim 650$~ms postbounce, and the subsequent fallback accretion onto the proto-neutron star (PNS) results in BH formation. In addition, the standing~accretion~shock~instability induces rotation of the proto-neutron star in the model that started with a non-rotating progenitor. Assuming conservation of specific angular momentum during black hole formation, this corresponds to a black~hole spin parameter of $a=J/M=0.046$. However, if no explosion sets in, all the angular momentum will eventually be accreted by the BH, resulting in a non-spinning BH. The successful explosion of the slowly-rotating model drastically slows down the accretion onto the PNS, allowing continued cooling and contraction that results in an extremely high gravitational-wave frequency ($f\sim3000$~Hz) at black~hole formation, while the non-rotating model generates gravitational wave signals similar to our corresponding 2D simulations. △ Less

Submitted 20 April, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: 17 pages, 11 figures. Accepted for publication in ApJ

Heavy Elements Nucleosynthesis On Accreting White Dwarfs: building seeds for the p-process

Authors: Umberto Battino, Marco Pignatari, Claudia Travaglio, Claudia Lederer-Woods, Pavel Denissenkov, Falk Herwig, Friedrich-Karl Thielemann, Thomas Rauscher

Abstract: The origin of the proton-rich trans-iron isotopes in the solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor's rapidly accreting white dwarf phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf models with initial masses >… ▽ More The origin of the proton-rich trans-iron isotopes in the solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor's rapidly accreting white dwarf phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf models with initial masses >~ 0.85Msun using the stellar code MESA. The near-surface layers of the 1, 1.26, 1.32 and 1.38Msun models are most representative of the regions in which the bulk of the p nuclei are produced during SNIa explosions, and for these models we also calculate the neutron-capture nucleosynthesis in the external layers. Contrary to previous rapidly-accreting white dwarf models at lower mass, we find that the H-shell ashes are the main site of n-capture nucleosynthesis. We find high neutron densities up to several 10^15 cm^-3 in the most massive WDs. Through the recurrence of the H-shell ashes these intermediate neutron densities can be sustained effectively for a long time leading to high neutron exposures with a strong production up to Pb. Both the neutron density and the neutron exposure increase with increasing the mass of the accreting WD. Finally, the SNIa nucleosynthesis is calculated using the obtained abundances as seeds. We obtain solar to super-solar abundances for p-nuclei with A>96. Our models show that SNIa are a viable p-process production site. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS. 49 pages, 17 figures, 5 tables

Core-collapse Supernova Explosions Driven by the Hadron-quark Phase Transition as a Rare $r$-process Site

Authors: Tobias Fischer, Meng-Ru Wu, Benjamin Wehmeyer, Niels-Uwe F. Bastian, Gabriel Martínez-Pinedo, Friedrich-Karl Thielemann

Abstract: Supernova explosions of massive stars are one of the primary sites for the production of the elements in the universe. Up to now, stars with zero-age main-sequence masses in the range of 35--50~$M_\odot$ had mostly been representing the failed supernova explosion branch. In contrast, it has been demonstrated recently that the appearance of exotic phases of hot and dense matter, associated with a s… ▽ More Supernova explosions of massive stars are one of the primary sites for the production of the elements in the universe. Up to now, stars with zero-age main-sequence masses in the range of 35--50~$M_\odot$ had mostly been representing the failed supernova explosion branch. In contrast, it has been demonstrated recently that the appearance of exotic phases of hot and dense matter, associated with a sufficiently strong phase transition from nuclear matter to the quark-gluon plasma at high baryon density, can trigger supernova explosions of such massive supergiant. Here, we present the first results obtained from an extensive nucleosynthesis analysis for material being ejected from the surface of the newly born proto-neutron star of such supernova explosions. These ejecta contain an early neutron-rich component and a late-time high-entropy neutrino-driven wind. The nucleosynthesis robustly overcomes the production of nuclei associated with the second $r$-process peak, at nuclear mass number $A\simeq 130$, and proceeds beyond the formation of the third peak ($A\simeq 195$) to the actinides. These yields may account for metal-poor star observations concerning $r$-process elements such as strontium and europium in the Galaxy at low metalicity, while the actinide yields suggests that this source may be a candidate contributing to the abundances of radioactive $^{244}$Pu measured in deep-sea sediments on Earth. △ Less

Submitted 7 May, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: 17 pages, 7 figures, The Astrophysical Journal, Volume 894, Issue 1, id.9

Journal ref: Astrophys. J. 894 (2020) 9

PUSHing core-collapse supernovae to explosions in spherical symmetry IV: Explodability, remnant properties and nucleosynthesis yields of low metallicity stars

Authors: Kevin Ebinger, Sanjana Curtis, Somdutta Ghosh, Carla Fröhlich, Matthias Hempel, Albino Perego, Matthias Liebendörfer, Friedrich-Karl Thielemann

Abstract: In this fourth paper of the series, we use the parametrized, spherically symmetric explosion method PUSH to perform a systematic study of two sets of non-rotating stellar progenitor models. Our study includes pre-explosion models with metallicities Z=0 and Z=Z$_{\odot}\times 10^{-4}$ and covers a progenitor mass range from 11 up to 75 M$_\odot$. We present and discuss the explosion properties of a… ▽ More In this fourth paper of the series, we use the parametrized, spherically symmetric explosion method PUSH to perform a systematic study of two sets of non-rotating stellar progenitor models. Our study includes pre-explosion models with metallicities Z=0 and Z=Z$_{\odot}\times 10^{-4}$ and covers a progenitor mass range from 11 up to 75 M$_\odot$. We present and discuss the explosion properties of all models and predict remnant (neutron star or black hole) mass distributions within this approach. We also perform systematic nucleosynthesis studies and predict detailed isotopic yields as function of the progenitor mass and metallicity. We present a comparison of our nucleosynthesis results with observationally derived $^{56}$Ni ejecta from normal core-collapse supernovae and with iron-group abundances for metal-poor star HD~84937. Overall, our results for explosion energies, remnant mass distribution, $^{56}$Ni mass, and iron group yields are consistent with observations of normal CCSNe. We find that stellar progenitors at low and zero metallicity are more prone to BH formation than those at solar metallicity, which allows for the formation of BHs in the mass range observed by LIGO/VIRGO. △ Less

Submitted 16 February, 2020; v1 submitted 20 October, 2019; originally announced October 2019.

Comments: 22 pages, 16 figures, 6 tables. Accepted to ApJ

Using failed supernovae to constrain the Galactic r-process element production

Authors: B. Wehmeyer, C. Frohlich, B. Côté, M. Pignatari, F. -K. Thielemann

Abstract: Rapid neutron capture process (r-process) elements have been detected in a large fraction of metal-poor halo stars, with abundances relative to iron (Fe) that vary by over two orders of magnitude. This scatter is reduced to less than a factor of 3 in younger Galactic disc stars. The large scatter of r-process elements in the early Galaxy suggests that the r-process is made by rare events, like com… ▽ More Rapid neutron capture process (r-process) elements have been detected in a large fraction of metal-poor halo stars, with abundances relative to iron (Fe) that vary by over two orders of magnitude. This scatter is reduced to less than a factor of 3 in younger Galactic disc stars. The large scatter of r-process elements in the early Galaxy suggests that the r-process is made by rare events, like compact binary mergers and rare sub-classes of supernovae. Although being rare, neutron star mergers alone have difficulties to explain the observed enhancement of r-process elements in the lowest metallicity stars compared to Fe. The supernovae producing the two neutron stars already provide a substantial Fe abundance where the r-process ejecta from the merger would be injected. In this work we investigate another complementary scenario, where the r-process occurs in neutron star-black hole mergers in addition to neutron star mergers. Neutron star-black hole mergers would eject similar amounts of r-process matter as neutron star mergers, but only the neutron star progenitor would have produced Fe. Furthermore, a reduced efficiency of Fe production from single stars significantly alters the age-metallicity relation, which shifts the onset of r-process production to lower metallicities. We use the high-resolution [(20 pc)3/cell] inhomogeneous chemical evolution tool `ICE' to study the outcomes of these effects. In our simulations, an adequate combination of neutron star mergers and neutron star-black hole mergers qualitatively reproduces the observed r-process abundances in the Galaxy. △ Less

Submitted 15 August, 2019; originally announced August 2019.

Journal ref: 2019 Monthly Notices of the Royal Astronomical Society, Volume 487, Issue 2, p.1745-1753

Enrichment of the Galactic disc with neutron-capture elements: Mo and Ru

Authors: Tamara Mishenina, Marco Pignatari, Tatiana Gorbaneva, Claudia Travaglio, Benoit Côté, Friedrich-Karl Thielemann, Caroline Soubiran

Abstract: We present new observational data for the heavy elements molybdenum (Mo, Z = 42) and ruthenium (Ru, Z = 44) in F-, G-, and K-stars belonging to different substructures of the Milky Way. The range of metallicity covered is --1.0 $<$ [Fe/H] $<$ +0.3. The spectra of Galactic disc stars have a high resolution of 42,000 and 75,000 and signal-to-noise ratio better than 100. Mo and Ru abundances were der… ▽ More We present new observational data for the heavy elements molybdenum (Mo, Z = 42) and ruthenium (Ru, Z = 44) in F-, G-, and K-stars belonging to different substructures of the Milky Way. The range of metallicity covered is --1.0 $<$ [Fe/H] $<$ +0.3. The spectra of Galactic disc stars have a high resolution of 42,000 and 75,000 and signal-to-noise ratio better than 100. Mo and Ru abundances were derived by comparing the observed and synthetic spectra in the region of Mo I lines at 5506, 5533 Å~ for 209 stars and Ru I lines at 4080, 4584, 4757 Å~ for 162 stars using the LTE approach. For all the stars, the Mo and Ru abundance determinations are obtained for the first time with an average error of 0.14 dex. This is the first extended sample of stellar observations for Mo and Ru in the Milky Way disc, and together with earlier observations in halo stars it is pivotal in providing a complete picture of the evolution of Mo and Ru across cosmic timescales. The Mo and Ru abundances were compared with those of the neutron-capture elements (Sr, Y, Zr, Ba, Sm, Eu). The complex nucleosynthesis history of Mo and Ru is compared with different Galactic Chemical Evolution (GCE) simulations. In general, present theoretical GCE simulations show underproduction of Mo and Ru at all metallicities compared to observations. This highlights a significant contribution of nucleosynthesis processes not yet considered in our simulations. A number of possible scenarios are discussed. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: 23 pages, 12 figures, 7 tables, accepted in MNRAS

Ultra-heavy cosmic-ray science--Are r-process nuclei in the cosmic rays produced in supernovae or binary neutron star mergers?

Authors: W. R. Binns, M. H. Israel, B. F. Rauch, A. C. Cummings, A. J. Davis, A. W. Labrador, R. A. Leske, R. A Mewaldt, E. C. Stone, M. E. Wiedenbeck, T. J. Brandt, E. R. Christian, J. T. Link, J. W. Mitchell, G. A. de Nolfo, T. T. von Rosenvinge, K. Sakai, M. Sasaki, C. J. Waddington, H. T. Janka, A. L. Melott, G. M. Mason, E-S. Seo, J. H. Adams, F-K. Thielemann , et al. (3 additional authors not shown)

Abstract: The recent detection of 60Fe in the cosmic rays provides conclusive evidence that there is a recently synthesized component (few MY) in the GCRs (Binns et al. 2016). In addition, these nuclei must have been synthesized and accelerated in supernovae near the solar system, probably in the Sco-Cen OB association subgroups, which are about 100 pc distant from the Sun. Recent theoretical work on the pr… ▽ More The recent detection of 60Fe in the cosmic rays provides conclusive evidence that there is a recently synthesized component (few MY) in the GCRs (Binns et al. 2016). In addition, these nuclei must have been synthesized and accelerated in supernovae near the solar system, probably in the Sco-Cen OB association subgroups, which are about 100 pc distant from the Sun. Recent theoretical work on the production of r-process nuclei appears to indicate that it is difficult for SNe to produce the solar system abundances relative to iron of r-process elements with high atomic number (Z), including the actinides (Th, U, Np, Pu, and Cm). Instead, it is believed by many that the heaviest r-process nuclei, or perhaps even all r-process nuclei, are produced in binary neutron star mergers. Since we now know that there is at least a component of the GCRs that has been recently synthesized and accelerated, models of r-process production by SNe and BNSM can be tested by measuring the relative abundances of these ultra-heavy r-process nuclei, and especially the actinides, since they are radioactive and provide clocks that give the time interval from nucleosynthesis to detection at Earth. Since BNSM are believed to be much less frequent in our galaxy than SNe (roughly 1000 times less frequent, the ratios of the actinides, each with their own half-life, will enable a clear determination of whether the heaviest r-process nuclei are synthesized in SNe or in BNSM. In addition, the r-process nuclei for the charge range from 34 to 82 can be used to constrain models of r-process production in BNSM and SNe. Thus, GCRs become a multi-messenger component in the study of BNSM and SNe. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: Astro2020 Science White Paper

Near-Earth Supernova Explosions: Evidence, Implications, and Opportunities

Authors: Brian D. Fields, John R. Ellis, Walter R. Binns, Dieter Breitschwerdt, Georgia A. de Nolfo, Roland Diehl, Vikram V. Dwarkadas, Adrienne Ertel, Thomas Faestermann, Jenny Feige, Caroline Fitoussi, Priscilla Frisch, David Graham, Brian Haley, Alexander Heger, Wolfgang Hillebrandt, Martin H. Israel, Thomas Janka, Michael Kachelriess, Gunther Korschinek, Marco Limongi, Maria Lugaro, Franciole Marinho, Adrian Melott, Richard A. Mewaldt , et al. (14 additional authors not shown)

Abstract: There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact o… ▽ More There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures. Astro2020 Science White Paper submitted to the 2020 Decadal Survey on Astronomy and Astrophysics

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

Enrichment of the Galactic disc with neutron capture elements: Sr

Authors: T. Mishenina, M. Pignatari, T. Gorbaneva, S. Bisterzo, C. Travaglio, F. K. Thielemann, C. Soubiran

Abstract: The enrichment history of heavy neutron-capture elements in the Milky Way disc provides fundamental information about the chemical evolution of our Galaxy and about the stellar sources that made those elements. In this work we give new observational data for Sr, the element at the first neutron-shell closure beyond iron, N=50, based on the analysis of the high resolution spectra of 276 Galactic di… ▽ More The enrichment history of heavy neutron-capture elements in the Milky Way disc provides fundamental information about the chemical evolution of our Galaxy and about the stellar sources that made those elements. In this work we give new observational data for Sr, the element at the first neutron-shell closure beyond iron, N=50, based on the analysis of the high resolution spectra of 276 Galactic disc stars. The Sr abundance was derived by comparing the observed and synthetic spectra in the region of the SrI 4607 A line, making use of the LTE approximation. NLTE corrections lead to an increase of the abundance estimates obtained under LTE, but for these lines they are minor near solar metallicity. The average correction that we find is 0.151 dex. The star that is mostly affected is HD 6582, with a 0.244 dex correction. The behavior of the Sr abundance as a function of metallicity is discussed within a stellar nucleosynthesis context, in comparison with the abundance of the heavy neutron-capture elements Ba (Z=56) and Eu (Z=63). The comparison of the observational data with the current GCE models confirm that the s-process contributions from Asymptotic Giant Branch stars and from massive stars are the main sources of Sr in the Galactic disc and in the Sun, while different nucleosynthesis sources can explain the high [Sr/Ba] and [Sr/Eu] ratios observed in the early Galaxy. △ Less

Submitted 25 January, 2019; originally announced January 2019.

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

Origin of the Heaviest Elements: the Rapid Neutron-Capture Process

Authors: John J. Cowan, Christopher Sneden, James E. Lawler, Ani Aprahamian, Michael Wiescher, Karlheinz Langanke, Gabriel Martínez-Pinedo, Friedrich-Karl Thielemann

Abstract: The production of about half of the heavy elements found in nature is assigned to a specific astrophysical nucleosynthesis process: the rapid neutron capture process (r-process). Although this idea has been postulated more than six decades ago, the full understanding faces two types of uncertainties/open questions: (a) The nucleosynthesis path in the nuclear chart runs close to the neutron-drip li… ▽ More The production of about half of the heavy elements found in nature is assigned to a specific astrophysical nucleosynthesis process: the rapid neutron capture process (r-process). Although this idea has been postulated more than six decades ago, the full understanding faces two types of uncertainties/open questions: (a) The nucleosynthesis path in the nuclear chart runs close to the neutron-drip line, where presently only limited experimental information is available, and one has to rely strongly on theoretical predictions for nuclear properties. (b) While for many years the occurrence of the r-process has been associated with supernovae, more recent studies have cast substantial doubts on this environment. Alternative scenarios include the mergers of neutron stars, neutron-star black hole mergers, but possibly also rare classes of supernovae as well as hypernovae/collapsars with polar jet ejecta and also accretion disk outflows related to the collapse of fast rotating massive stars with high magnetic fields. Stellar r-process abundance observations, have provided insights into, and constraints on the frequency of and conditions in the responsible stellar production sites. One of them, neutron star mergers, was just identified and related to the Gravitational Wave event GW170817. High resolution observations, increasingly more precise due to improved experimental atomic data, have been particularly important in defining the heavy element abundance patterns of the old halo stars, and thus determining the extent, and nature, of the earliest nucleosynthesis in our Galaxy. Combining new results and important breakthroughs in the related nuclear, atomic and astronomical fields of science, this review attempts to provide an answer to the question "How Were the Elements from Iron to Uranium Made?" (Abridged) △ Less

Submitted 1 February, 2021; v1 submitted 5 January, 2019; originally announced January 2019.

Comments: 96 pages, 42 figures, 2 tables, author copy manuscript published in Reviews of Modern Physics

Journal ref: Rev. Mod. Phys. 93, 15002 (2021)

Observatory science with eXTP

Authors: Jean J. M. in 't Zand, Enrico Bozzo, Jinlu Qu, Xiang-Dong Li, Lorenzo Amati, Yang Chen, Immacolata Donnarumma, Victor Doroshenko, Stephen A. Drake, Margarita Hernanz, Peter A. Jenke, Thomas J. Maccarone, Simin Mahmoodifar, Domitilla de Martino, Alessandra De Rosa, Elena M. Rossi, Antonia Rowlinson, Gloria Sala, Giulia Stratta, Thomas M. Tauris, Joern Wilms, Xuefeng Wu, Ping Zhou, Iván Agudo, Diego Altamirano , et al. (159 additional authors not shown)

Abstract: In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to stu… ▽ More In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

Explosive Nucleosynthesis: What we learned and what we still do not understand

Authors: Friedrich-Karl Thielemann

Abstract: This review touches on historical aspects, going back to the early days of nuclear astrophysics, initiated by B$^2$FH and Cameron, discusses (i) the required nuclear input from reaction rates and decay properties up to the nuclear equation of state, continues (ii) with the tools to perform nucleosynthesis calculations and (iii) early parametrized nucleosynthesis studies, before (iv) reliable stell… ▽ More This review touches on historical aspects, going back to the early days of nuclear astrophysics, initiated by B$^2$FH and Cameron, discusses (i) the required nuclear input from reaction rates and decay properties up to the nuclear equation of state, continues (ii) with the tools to perform nucleosynthesis calculations and (iii) early parametrized nucleosynthesis studies, before (iv) reliable stellar models became available for the late stages of stellar evolution. It passes then through (v) explosive environments from core-collapse supernovae to explosive events in binary systems (including type Ia supernovae and compact binary mergers), and finally (vi) discusses the role of all these nucleosynthesis production sites in the evolution of galaxies. The focus is put on the comparison of early ideas and present, very recent, understanding. △ Less

Submitted 19 September, 2018; originally announced September 2018.

Comments: 11 pages, to appear in Springer Proceedings in Physics (Proc. of Intl. Conf. "Nuclei in the Cosmos XV", LNGS Assergi, Italy, June 2018)

Women Scientists Who Made Nuclear Astrophysics

Authors: C. V. Hampton, M. Lugaro, P. Papakonstantinou, P. G. Isar, B. Nordström, N. Özkan, M. Aliotta, A. Ćiprijanović, S. Curtis, M. Di Criscienzo, J. den Hartogh, A. S. Font, A. Kankainen, C. Kobayashi, C. Lederer-Woods, E. Niemczura, T. Rauscher, A. Spyrou, S. Van Eck, M. Yavahchova, W. Chantereau, S. E. de Mink, E. Kaiser, F. -K. Thielemann, C. Travaglio , et al. (2 additional authors not shown)

Abstract: Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models… ▽ More Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models for women scientists; instead, it often renders these women invisible. In response to this situation, we present a selection of twelve outstanding women who helped to develop nuclear astrophysics. △ Less

Submitted 25 August, 2018; originally announced September 2018.

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

Core-collapse supernovae in the hall of mirrors. A three-dimensional code-comparison project

Authors: Rubén M. Cabezón, Kuo-Chuan Pan, Matthias Liebendörfer, Takami Kuroda, Kevin Ebinger, Oliver Heinimann, Friedrich-Karl Thielemann, Albino Perego

Abstract: Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotatio… ▽ More Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotation on the simulation of CCSNe in 3D. We use four different hydrodynamics codes in this work (ELEPHANT, FLASH, fGR1, and SPHYNX) in combination with two different neutrino treatments, the isotropic diffusion source approximation (IDSA) and two-moment M1, and three different gravity treatments: Newtonian, 1D General Relativity (GR) correction, and full GR). Additional parameters discussed in this study are the inclusion of neutrino-electron scattering via a parametrized deleptonization (PD) and the influence of rotation. The four codes compared in this work include Eulerian and fully Lagrangian (smoothed particle hydrodynamics) codes for the first time. They show agreement in the overall evolution of the collapse phase and early post-bounce within the range of 10% (20% in some cases). The comparison of the different neutrino treatments highlights the need to further investigate the antineutrino luminosities in IDSA, which tend to be relatively high. We also demonstrate the requirement for a more detailed heavy-lepton neutrino leakage. When comparing with a full GR code, including an M1 transport method, we confirm the influence of neutrino-electron scattering during the collapse phase, which is adequately captured by the PD scheme. Also, the effective GR potential reproduces the overall dynamic evolution correctly in all Newtonian codes. Additionally, we verify that rotation aids the shock expansion and estimate the overall angular momentum losses for each code in rotating scenarios. △ Less

Submitted 18 December, 2018; v1 submitted 24 June, 2018; originally announced June 2018.

Comments: 20 pages, 11 figures, 4 Tables. Published in A&A

Journal ref: A&A 619, A118 (2018)

PUSHing Core-Collapse Supernovae to Explosions in Spherical Symmetry III: Nucleosynthesis Yields

Authors: Sanjana Curtis, Kevin Ebinger, Carla Fröhlich, Matthias Hempel, Albino Perego, Matthias Liebendörfer, Friedrich-Karl Thielemann

Abstract: In a previously presented proof-of-principle study, we established a parametrized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae for a wide range of pre-explosion models. The method is based on the neutrino-driven mechanism and follows collapse, bounce and explosion. There are two crucial aspects of our model for nucleosynthesis predictio… ▽ More In a previously presented proof-of-principle study, we established a parametrized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae for a wide range of pre-explosion models. The method is based on the neutrino-driven mechanism and follows collapse, bounce and explosion. There are two crucial aspects of our model for nucleosynthesis predictions. First, the mass cut and explosion energy emerge simultaneously from the simulation (determining, for each stellar model, the amount of Fe-group ejecta). Second, the interactions between neutrinos and matter are included consistently (setting the electron fraction of the innermost ejecta). In the present paper, we use the successful explosion models from Ebinger et al. (2018) which include two sets of pre-explosion models at solar metallicity, with combined masses between 10.8 and 120 M$_{\odot}$. We perform systematic nucleosynthesis studies and predict detailed isotopic yields. The resulting $^{56}$Ni ejecta are in overall agreement with observationally derived values from normal core-collapse supernovae. The Fe-group yields are also in agreement with derived abundances for metal-poor star HD84937. We also present a comparison of our results with observational trends in alpha element to iron ratios. △ Less

Submitted 10 December, 2018; v1 submitted 1 May, 2018; originally announced May 2018.

Comments: 24 pages, 10 figures, 8 tables (Accepted, ApJ)

PUSHing Core-Collapse Supernovae to Explosions in Spherical Symmetry II: Explodability and Global Properties

Authors: Kevin Ebinger, Sanjana Curtis, Carla Fröhlich, Matthias Hempel, Albino Perego, Matthias Liebendörfer, Friedrich-Karl Thielemann

Abstract: In a previously presented proof-of-principle study we established a parametrized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae. The present paper goes beyond a specific application that is able to reproduce observational properties of SN1987A and performs a systematic study of the explosion properties for an extensive set of non-rotating… ▽ More In a previously presented proof-of-principle study we established a parametrized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae. The present paper goes beyond a specific application that is able to reproduce observational properties of SN1987A and performs a systematic study of the explosion properties for an extensive set of non-rotating, solar metallicity stellar progenitor models in the mass range from 10.8 to 120 M$_\odot$.This includes the transition from neutron stars to black holes as the final result of the collapse of massive stars, and the relation of the latter to supernovae and faint/failed supernovae. The present paper provides the basis for extended nucleosynthesis predictions in a forthcoming paper to be employed in galactic evolution models. △ Less

Submitted 16 February, 2020; v1 submitted 9 April, 2018; originally announced April 2018.

Comments: 24 pages, 18 figures, 5 tables. Accepted to ApJ

Mass Measurements of Neutron-Deficient Y, Zr, and Nb Isotopes and Their Impact on $rp$ and $νp$ Nucleosynthesis Processes

Authors: Y. M. Xing, K. A. Li, Y. H. Zhang, X. H. Zhou, M. Wang, Yu. A. Litvinov, K. Blaum, S. Wanajo, S. Kubono, G. Martínez-Pinedo, A. Sieverding, R. J. Chen, P. Shuai, C. Y. Fu, X. L. Yan, W. J. Huang, X. Xu, X. D. Tang, H. S. Xu, T. Bao, X. C. Chen, B. S. Gao, J. J. He, Y. H. Lam, H. F. Li , et al. (26 additional authors not shown)

Abstract: Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than… ▽ More Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $α$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $νp$-process simulations. △ Less

Submitted 6 April, 2018; originally announced April 2018.

A full general relativistic neutrino radiation-hydrodynamics simulation of a collapsing very massive star and the formation of a black hole

Authors: Takami Kuroda, Kei Kotake, Tomoya Takiwaki, Friedrich-Karl Thielemann

Abstract: We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leadin… ▽ More We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leading to black hole (BH) formation. Our results show that the BH formation occurs at a post-bounce time of ~300 ms for the 70 Msun star. This is significantly earlier than those in the literature where lower mass progenitors were employed. At a few ~10 ms before BH formation, we find that the stalled bounce shock is revived by intense neutrino heating from the very hot PNS, which is aided by violent convection behind the shock. In the context of 3D-GR core-collapse modeling with multi-energy neutrino transport, our numerical results present the first evidence to validate a fallback BH formation scenario of the 70 Msun star. △ Less

Submitted 4 April, 2018; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS letter

Neutron Star Mergers and Nucleosynthesis of Heavy Elements

Authors: F. -K. Thielemann, M. Eichler, I. V. Panov, B. Wehmeyer

Abstract: The existence of neutron star mergers has been supported since the discovery of the binary pulsar and the observation of its orbital energy loss, consistent with General Relativity. They are considered nucleosynthesis sites of the rapid neutron-capture process (r-process), which is responsible for creating approximately half of all heavy elements beyond Fe and is the only source of elements beyond… ▽ More The existence of neutron star mergers has been supported since the discovery of the binary pulsar and the observation of its orbital energy loss, consistent with General Relativity. They are considered nucleosynthesis sites of the rapid neutron-capture process (r-process), which is responsible for creating approximately half of all heavy elements beyond Fe and is the only source of elements beyond Pb and Bi. Detailed nucleosynthesis calculations based on the decompression of neutron star matter are consistent with solar r-process abundances of heavy nuclei. Neutron star mergers have also been identified with short-duration γ-ray bursts via their IR afterglow. The high neutron densities in ejected matter permit a violent r-process, leading to fission cycling of the heaviest nuclei in regions far from (nuclear) stability. Uncertainties in several nuclear properties affect the abundance distributions. The modeling of astrophysical events also depends on the hydrodynamic treatment, the occurrence of a neutrino wind after the merger and before the possible emergence of a black hole, and the properties of black hole accretion disks. We discuss the effect of nuclear and modeling uncertainties and conclude that binary compact mergers are probably a (or the) dominant site of the production of r-process nuclei in our Galaxy. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: preprint version of a review to appear in Ann. Rev. Nucl. Part. Sci. 67 (2017), literature review up to February 2017

Equation of State Dependent Dynamics and Multimessenger Signals from Stellar-mass Black Hole Formation

Authors: Kuo-Chuan Pan, Matthias Liebendörfer, Sean M. Couch, Friedrich-Karl Thielemann

Abstract: We investigate axisymmetric black hole~(BH) formation and its gravitational wave (GW) and neutrino signals with self-consistent core-collapse supernova simulations of a non-rotating $40~M_\odot$ progenitor star using the isotropic diffusion source approximation for the neutrino transport and a modified gravitational potential for general relativistic effects. We consider four different neutron sta… ▽ More We investigate axisymmetric black hole~(BH) formation and its gravitational wave (GW) and neutrino signals with self-consistent core-collapse supernova simulations of a non-rotating $40~M_\odot$ progenitor star using the isotropic diffusion source approximation for the neutrino transport and a modified gravitational potential for general relativistic effects. We consider four different neutron star (NS) equations of state~(EoS): LS220, SFHo, BHB$Λφ$ and DD2, and study the impact of the EoS on BH formation dynamics and GW emission. We find that the BH formation time is sensitive to the EoS from 460 to $>1300$~ms and is delayed in multiple dimensions for $\sim~100-250$~ms due to the finite entropy effects. Depending on the EoS, our simulations show the possibility that shock revival can occur along with the collapse of the proto-neutron star~(PNS) to a BH. The gravitational waveforms contain four major features that are similar to previous studies but show extreme values: (1)~a low frequency signal ($\sim~300-500$~Hz) from core-bounce and prompt convection, (2)~a strong signal from the PNS g-mode oscillation among other features, (3)~a high frequency signal from the PNS inner-core convection, and (4)~signals from the standing accretion shock instability and convection. The peak frequency at the onset of BH formation reaches to $\sim~2.3$~kHz. The characteristic amplitude of a 10~kpc object at peak frequency is detectable but close to the noise threshold of the Advanced~LIGO and KAGRA, suggesting that the next generation gravitational wave detector will need to improve the sensitivity at the kHz domain to better observe stellar-mass BH formation from core-collapse supernovae or failed supernovae. △ Less

Submitted 9 April, 2018; v1 submitted 4 October, 2017; originally announced October 2017.

Comments: 8 figures. Accepted by ApJ

Nucleosynthesis in 2D Core-Collapse Supernovae of 11.2 and 17.0 M$_{\odot}$ Progenitors: Implications for Mo and Ru Production

Authors: M. Eichler, K. Nakamura, T. Takiwaki, T. Kuroda, K. Kotake, M. Hempel, R. Cabezón, M. Liebendörfer, F-K. Thielemann

Abstract: Core-collapse supernovae are the first polluters of heavy elements in the galactic history. As such, it is important to study the nuclear compositions of their ejecta, and understand their dependence on the progenitor structure (e.g., mass, compactness, metallicity). Here, we present a detailed nucleosynthesis study based on two long-term, two-dimensional core-collapse supernova simulations of a 1… ▽ More Core-collapse supernovae are the first polluters of heavy elements in the galactic history. As such, it is important to study the nuclear compositions of their ejecta, and understand their dependence on the progenitor structure (e.g., mass, compactness, metallicity). Here, we present a detailed nucleosynthesis study based on two long-term, two-dimensional core-collapse supernova simulations of a 11.2 M$_{\odot}$ and a 17.0 M$_{\odot}$ star. We find that in both models nuclei well beyond the iron group (up to $Z \approx 44$) can be produced, and discuss in detail also the nucleosynthesis of the p-nuclei $^{92,94}$Mo and $^{96,98}$Ru. While we observe the production of $^{92}$Mo and $^{94}$Mo in slightly neutron-rich conditions in both simulations, $^{96,98}$Ru can only be produced efficiently via the $ν$p-process. Furthermore, the production of Ru in the $ν$p-process heavily depends on the presence of very proton-rich material in the ejecta. This disentanglement of production mechanisms has interesting consequences when comparing to the abundance ratios between these isotopes in the solar system and in presolar grains. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: 48 pages, 19 figures, accepted for publication in: J. Phys. G: Nucl. Part. Phys

Observing the metal-poor solar neighbourhood: a comparison of galactic chemical evolution predictions

Authors: T. Mishenina, M. Pignatari, B. Cot'e, F. -K. Thielemann, C. Soubiran, N. Basak, T. Gorbaneva, S. A. Korotin, V. V. Kovtyukh, B. Wehmeyer, S. Bisterzo, C. Travaglio, B. K. Gibson, C. Jordan, A. Paul, C. Ritter, F. Herwig

Abstract: Atmospheric parameters and chemical compositions for ten stars with metallicities in the region of -2.2< [Fe/H] <-0.6 were precisely determined using high resolution, high signal to noise, spectra. For each star the abundances, for 14 to 27 elements, were derived using both LTE and NLTE approaches. In particular, differences by assuming LTE or NLTE are about 0.10 dex; depending on [Fe/H], Teff, gr… ▽ More Atmospheric parameters and chemical compositions for ten stars with metallicities in the region of -2.2< [Fe/H] <-0.6 were precisely determined using high resolution, high signal to noise, spectra. For each star the abundances, for 14 to 27 elements, were derived using both LTE and NLTE approaches. In particular, differences by assuming LTE or NLTE are about 0.10 dex; depending on [Fe/H], Teff, gravity and element lines used in the analysis. We find that the O abundance has the largest error, ranging from 0.10 and 0.2 dex. The best measured elements are Cr, Fe, and Mn; with errors etween 0.03 and 0.11 dex. The stars in our sample were included in previous different observational work. We provide a consistent data analysis. The data dispersion introduced in the literature by different techniques and assumptions used by the different authors is within the observational errors, excepting for HD103095. We compare these results with stellar observations from different data sets and a number of theoretical galactic chemical evolution (GCE) simulations. We find a large scatter in the GCE results, used to study the origin of the elements. Within this scatter as found in previous GCE simulations, we cannot reproduce the evolution of the elemental ratios [Sc/Fe], [Ti/Fe], and [V/Fe] at different metallicities. The stellar yields from core collapse supernovae (CCSN) are likely primarily responsible for this discrepancy. Possible solutions and open problems are discussed. △ Less

Submitted 10 May, 2017; originally announced May 2017.

Comments: 25 pages, 7 figures

Multi-dimensional Core-Collapse Supernova Simulations with Neutrino Transport

Authors: Kuo-Chuan Pan, Matthias Liebendörfer, Matthias Hempel, Friedrich-Karl Thielemann

Abstract: We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitati… ▽ More We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitation of GPUs, but also sophisticated approximations to capture the essential microphysics. We demonstrate that the IDSA is an elegant and efficient neutrino radiation transfer scheme, which is portable to multiple hydrodynamics codes and fast enough to investigate long-term evolutions in two and three dimensions. Simulations with a 40 solar mass progenitor are presented in both FLASH (1D and 2D) and ELEPHANT (3D) as an extreme test condition. It is found that the black hole formation time is delayed in multiple dimensions and we argue that the strong standing accretion shock instability before black hole formation will lead to strong gravitational waves. △ Less

Submitted 23 January, 2017; originally announced January 2017.

Comments: 3 pages, proceedings for Nuclei in the Cosmos XIV, Niigata, Japan (2016)

PUSHing Core-Collapse Supernovae to Explosions in Spherical Symmetry: Nucleosynthesis Yields

Authors: Sanjana Sinha, Carla Fröhlich, Kevin Ebinger, Albino Perego, Matthias Hempel, Marius Eichler, Matthias Liebendörfer, Friedrich-Karl Thielemann

Abstract: Core-collapse supernovae (CCSNe) are the extremely energetic deaths of massive stars. They play a vital role in the synthesis and dissemination of many heavy elements in the universe. In the past, CCSN nucleosynthesis calculations have relied on artificial explosion methods that do not adequately capture the physics of the innermost layers of the star. The PUSH method, calibrated against SN1987A,… ▽ More Core-collapse supernovae (CCSNe) are the extremely energetic deaths of massive stars. They play a vital role in the synthesis and dissemination of many heavy elements in the universe. In the past, CCSN nucleosynthesis calculations have relied on artificial explosion methods that do not adequately capture the physics of the innermost layers of the star. The PUSH method, calibrated against SN1987A, utilizes the energy of heavy-flavor neutrinos emitted by the proto-neutron star (PNS) to trigger parametrized explosions. This makes it possible to follow the consistent evolution of the PNS and to ensure a more accurate treatment of the electron fraction of the ejecta. Here, we present the Iron group nucleosynthesis results for core-collapse supernovae, exploded with PUSH, for two different progenitor series. Comparisons of the calculated yields to observational metal-poor star data are also presented. Nucleosynthesis yields will be calculated for all elements and over a wide range of progenitor masses. These yields can be immensely useful for models of galactic chemical evolution. △ Less

Submitted 18 January, 2017; originally announced January 2017.

Comments: 3 pages, 3 figures, poster presentation to appear in the proceedings of the 14th International Symposium on Nuclei in the Cosmos (NIC-XIV), Ed. S. Kubono, JPS (Japan Physical Society)

The intermediate r-process in core-collapse supernovae driven by the magneto-rotational instability

Authors: Nobuya Nishimura, Hidetomo Sawai, Tomoya Takiwaki, Shoichi Yamada, Friedrich-Karl Thielemann

Abstract: We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability. A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the magneto-rotational instability. Neutrino-heating dominated explosions, enhanced by magne… ▽ More We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability. A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the magneto-rotational instability. Neutrino-heating dominated explosions, enhanced by magnetic fields, showed mildly neutron-rich ejecta producing nuclei up to $A \sim 130$ (i.e. the weak r-process), while explosion models with stronger magnetic fields reproduce a solar-like r-process pattern. More commonly seen abundance patterns in our models are in between the weak and regular r-process, producing lighter and intermediate mass nuclei. These {\it intermediate r-processes} exhibit a variety of abundance distributions, compatible with several abundance patterns in r-process-enhanced metal-poor stars. The amount of Eu ejecta $\sim 10^{-5} M_\odot$ in magnetically-driven jets agrees with predicted values in the chemical evolution of early galaxies. In contrast, neutrino-heating dominated explosions have a significant amount of Fe ($^{56}{\rm Ni}$) and Zn, comparable to regular supernovae and hypernovae, respectively. These results indicate magneto-rotational supernovae can produce a wide range of heavy nuclei from iron-group to r-process elements, depending on the explosion dynamics. △ Less

Submitted 2 February, 2017; v1 submitted 7 November, 2016; originally announced November 2016.

Comments: 8 pages, 5 figures, ApJL, accepted; v2: text and some plots are revised based on referee's comments, but the main conclusions remain the same

Journal ref: The Astrophysical Journal Letters, 836:L21, 2017

Explosion Dynamics of Parametrized Spherically Symmetric Core-Collapse Supernova Simulations

Authors: Kevin Ebinger, Sanjana Sinha, Carla Fröhlich, Albino Perego, Matthias Hempel, Marius Eichler, Jordi Casanova, Matthias Liebendörfer, Friedrich-Karl Thielemann

Abstract: We report on a method, PUSH, for triggering core-collapse supernova (CCSN) explosions of massive stars in spherical symmetry. This method provides a framework to study many important aspects of core collapse supernovae: the effects of the shock passage through the star, explosive supernova nucleosynthesis and the progenitor-remnant connection. Here we give an overview of the method, compare the re… ▽ More We report on a method, PUSH, for triggering core-collapse supernova (CCSN) explosions of massive stars in spherical symmetry. This method provides a framework to study many important aspects of core collapse supernovae: the effects of the shock passage through the star, explosive supernova nucleosynthesis and the progenitor-remnant connection. Here we give an overview of the method, compare the results to multi-dimensional simulations and investigate the effects of the progenitor and the equation of state on black hole formation. △ Less

Submitted 12 October, 2016; originally announced October 2016.

Comments: Proceedings for Nuclei in the Cosmos XIV, Niigata, Japan (2016)

Towards generating a new supernova equation of state: A systematic analysis of cold hybrid stars

Authors: Oliver Heinimann, Matthias Hempel, Friedrich-Karl Thielemann

Abstract: The hadron-quark phase transition in core-collapse supernovae (CCSNe) has the potential to trigger explosions in otherwise nonexploding models. However, those hybrid supernova equations of state (EOS) shown to trigger an explosion do not support the observational 2 M$_\odot$ neutron star maximum mass constraint. In this work, we analyze cold hybrid stars by the means of a systematic parameter scan… ▽ More The hadron-quark phase transition in core-collapse supernovae (CCSNe) has the potential to trigger explosions in otherwise nonexploding models. However, those hybrid supernova equations of state (EOS) shown to trigger an explosion do not support the observational 2 M$_\odot$ neutron star maximum mass constraint. In this work, we analyze cold hybrid stars by the means of a systematic parameter scan for the phase transition properties, with the aim to develop a new hybrid supernova EOS. The hadronic phase is described with the state-of-the-art supernova EOS HS(DD2), and quark matter by an EOS with a constant speed of sound (CSS) of $c_{QM}^2=1/3$. We find promising cases which meet the 2 M$_\odot$ criterion and are interesting for CCSN explosions. We show that the very simple CSS EOS is transferable into the well-known thermodynamic bag model, important for future application in CCSN simulations. In the second part, the occurrence of reconfinement and multiple phase transitions is discussed. In the last part, the influence of hyperons in our parameter scan is studied. Including hyperons no change in the general behavior is found, except for overall lower maximum masses. In both cases (with and without hyperons) we find that quark matter with $c_{QM}^2=1/3$ can increase the maximum mass only if reconfinement is suppressed or if quark matter is absolutely stable. △ Less

Submitted 22 November, 2016; v1 submitted 31 August, 2016; originally announced August 2016.

Comments: 14 pages, 11 figures, v2: matches published version

Journal ref: Phys. Rev. D 94, 103008 (2016)

Application of a Theory and Simulation based Convective Boundary Mixing model for AGB Star Evolution and Nucleosynthesis

Authors: U. Battino, M. Pignatari, C. Ritter, F. Herwig, P. Denisenkov, J. W. Den Hartogh, R. Trappitsch, R. Hirschi, B. Freytag, F. Thielemann, B. Paxton

Abstract: The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence… ▽ More The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up where the 13C pocket for the s process in AGB stars forms. In this work we apply a CBM model motivated by simulations and theory to models with initial mass $M = 2$ and $M = 3M_\odot$, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundance of 12C and 16O are increased by CBM at the bottom of pulse-driven convection zone. This mixing is affecting the $^{22}Ne(α,n)^{25}Mg$ activation and the s-process effciency in the 13C-pocket. In our model CBM at the bottom of the convective envelope during the third dredgeup represents gravity wave mixing. We take further into account that hydrodynamic simulations indicate a declining mixing efficiency already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the 13C-pocket with a mass of $\approx 10^{-4}M_\odot$. The final $s$-process abundances are characterized by 0.36 < [s=Fe] < 0.78 and the heavy-to-light s-process ratio is 0.23 < [hs=ls] < 0.45. Finally, we compare our results with stellar observations, pre-solar grain measurements and previous work. △ Less

Submitted 19 May, 2016; originally announced May 2016.

Comments: Submitted to ApJ on 11-24-2015. Accepted on 5-17-2016 (Manuscript #: ApJ101257)

Hot third family of compact stars and the possibility of core-collapse supernova explosions

Authors: Matthias Hempel, Oliver Heinimann, Andrey Yudin, Igor Iosilevskiy, Matthias Liebendörfer, Friedrich-Karl Thielemann

Abstract: A phase transition to quark matter can lead to interesting phenomenological consequences in core-collapse supernovae, e.g., triggering an explosion in spherically symmetric models. However, until now this explosion mechanism was only shown to be working for equations of state that are in contradiction with recent pulsar mass measurements. Here, we identify that this explosion mechanism is related… ▽ More A phase transition to quark matter can lead to interesting phenomenological consequences in core-collapse supernovae, e.g., triggering an explosion in spherically symmetric models. However, until now this explosion mechanism was only shown to be working for equations of state that are in contradiction with recent pulsar mass measurements. Here, we identify that this explosion mechanism is related to the existence of a third family of compact stars. For the equations of state investigated, the third family is only pronounced in the hot, early stages of the protocompact star and absent or negligibly small at zero temperature and thus represents a novel kind of third family. This interesting behavior is a result of unusual thermal properties induced by the phase transition, e.g., characterized by a decrease of temperature with increasing density for isentropes, and can be related to a negative slope of the phase transition line in the temperature-pressure phase diagram. △ Less

Submitted 11 November, 2016; v1 submitted 20 November, 2015; originally announced November 2015.

Comments: 6 pages, 4 figures, v2: matches published version

Journal ref: Phys. Rev. D 94, 103001 (2016)

$s$-process production in rotating massive stars at solar and low metallicities

Authors: Urs Frischknecht, Raphael Hirschi, Marco Pignatari, André Maeder, George Meynet, Cristina Chiappini, Friedrich-Karl Thielemann, Thomas Rauscher, Cyril Georgy, Sylvia Ekström

Abstract: Rotation was shown to have a strong impact on the structure and light element nucleosynthesis in massive stars. In particular, models including rotation can reproduce the primary nitrogen observed in halo extremely metal-poor (EMP) stars. Additional exploratory models showed that rotation may enhance $s$-process production at low metallicity. Here we present a large grid of massive star models i… ▽ More Rotation was shown to have a strong impact on the structure and light element nucleosynthesis in massive stars. In particular, models including rotation can reproduce the primary nitrogen observed in halo extremely metal-poor (EMP) stars. Additional exploratory models showed that rotation may enhance $s$-process production at low metallicity. Here we present a large grid of massive star models including rotation and a full $s$-process network to study the impact of rotation on the weak $s$-process. We explore the possibility of producing significant amounts of elements beyond the strontium peak, which is where the weak $s$-process usually stops. We used the Geneva stellar evolution code coupled to an enlarged reaction network with 737 nuclear species up to bismuth to calculate $15-40\,\text{M}_\odot$ models at four metallicities ($Z = 0.014,10^{-3}$, $10^{-5}$, and $10^{-7}$) from the main sequence up to the end of oxygen burning. We confirm that rotation-induced mixing between the convective H-shell and He-core enables an important production of primary $^{14}$N and $^{22}$Ne and $s$-process at low metallicity. At low metallicity, even though the production is still limited by the initial number of iron seeds, rotation enhances the $s$-process production, even for isotopes heavier than strontium, by increasing the neutron to seed ratio. The increase in this ratio is a direct consequence of the primary production of $^{22}$Ne. Despite nuclear uncertainties affecting the $s$-process production and stellar uncertainties affecting the rotation-induced mixing, our results show a robust production of $s$ process at low metallicity when rotation is taken into account. Considering models with a distribution of initial rotation rates enables to reproduce the observed large range of the [Sr/Ba] ratios in (carbon-enhanced and normal) EMP stars. △ Less

Submitted 18 November, 2015; originally announced November 2015.

Comments: 26 pages, 15 figures, 8 tables, accepted for publication in MNRAS

Journal ref: MNRAS 456 (2016) 1803-1825

Nucleosynthesis in the Ejecta of Neutron Star Mergers

Authors: Dirk Martin, Albino Perego, Almudena Arcones, Oleg Korobkin, Friedrich-Karl Thielemann

Abstract: Heavy elements like gold, platinum or uranium are produced in the r-process, which needs neutron-rich and explosive environments. Neutron star mergers are a promising candidate for an r-process site. They exhibit three different channels for matter ejection fulfilling these conditions: dynamic ejecta due to tidal torques, neutrino-driven winds and evaporating matter from the accretion disk. We pre… ▽ More Heavy elements like gold, platinum or uranium are produced in the r-process, which needs neutron-rich and explosive environments. Neutron star mergers are a promising candidate for an r-process site. They exhibit three different channels for matter ejection fulfilling these conditions: dynamic ejecta due to tidal torques, neutrino-driven winds and evaporating matter from the accretion disk. We present a first study of the integrated nucleosynthesis for a neutrino-driven wind from a neutron star merger with a hyper-massive neutron star. Trajectories from a recent hydrodynamical simulation are divided into four different angle regions and post-processed with a reaction network. We find that the electron fraction varies around $Y_e \approx 0.1 - 0.4$, but its distribution differs for every angle of ejection. Hence, the wind ejecta do not undergo a robust r-process, but rather possess distinct nucleosynthesis yields depending on the angle range. Compared to the dynamic ejecta, a smaller amount of neutron-rich matter gets unbound, but the production of lighter heavy elements with $A \lesssim 130$ in the neutrino-driven wind can complement the strong r-process of the dynamic ejecta. △ Less

Submitted 25 September, 2015; originally announced September 2015.

Comments: 6 pages, 4 figures, accepted for publication on PoS

Mn abundances in the stars of the Galactic disc with metallicities -1.0 < [Fe/H] < 0.3

Authors: T. Mishenina, T. Gorbaneva, M. Pignatari, F. -K. Thielemann, S. A. Korotin

Abstract: In this work we present and discuss the observations of the Mn abundances for 247 FGK dwarfs, located in the Galactic disc with metallicity -1<Fe/H]<+0.3. The observed stars belong to the substructures of the Galaxy thick and thin discs, and to the Hercules stream. The observations were conducted using the 1.93 m telescope at Observatoire de Haute-Provence (OHP, France) equipped with the echelle t… ▽ More In this work we present and discuss the observations of the Mn abundances for 247 FGK dwarfs, located in the Galactic disc with metallicity -1<Fe/H]<+0.3. The observed stars belong to the substructures of the Galaxy thick and thin discs, and to the Hercules stream. The observations were conducted using the 1.93 m telescope at Observatoire de Haute-Provence (OHP, France) equipped with the echelle type spectrographs ELODIE and SOPHIE. The abundances were derived under the LTE approximation, with an average error for the [Mn/Fe] ratio of 0.10 dex. For most of the stars in the sample Mn abundances are not available in the literature. We obtain an evolution of [Mn/Fe] ratio with the metallicity [Fe/H] consistent with previous data compilations. In particular, within the metallicity range covered by our stellar sample the [Mn/Fe] ratio is increasing with the increase of metallicity. This due to the contribution to the Galactic chemical evolution of Mn and Fe from thermonuclear supernovae. We confirm the baseline scenario where most of the Mn in the Galactic disc and in the Sun is made by thermonuclear supernovae. In particular, the effective contribution from core-collapse supernovae to the Mn in the Solar system is about 10-20%. However, present uncertainties affecting the production of Mn and Fe in thermonuclear supernovae are limiting the constraining power of the observed [Mn/Fe] trend in the Galactic discs on, e.g., the frequency of different thermonuclear supernovae populations. The different production of these two elements in different types of thermonuclear supernovae needs to be disentangled by the dependence of their relative production on the metallicity of the supernova progenitor. △ Less

Submitted 22 September, 2015; v1 submitted 17 September, 2015; originally announced September 2015.

Comments: 17 pages, 5 figures, revised version, accepted for publication in MNRAS

Type Ia Supernovae: Can Coriolis force break the symmetry of the gravitational confined detonation explosion mechanism?

Authors: D. García-Senz, R. M. Cabezón, I. Domínguez, F. K. Thielemann

Abstract: Nowadays the number of models aimed at explaining the Type Ia supernova phenomenon is high and discriminating between them is a must-do. In this work we explore the influence of rotation in the evolution of the nuclear flame which drives the explosion in the so called gravitational confined detonation models. Assuming that the flame starts in a point-like region slightly above the center of the wh… ▽ More Nowadays the number of models aimed at explaining the Type Ia supernova phenomenon is high and discriminating between them is a must-do. In this work we explore the influence of rotation in the evolution of the nuclear flame which drives the explosion in the so called gravitational confined detonation models. Assuming that the flame starts in a point-like region slightly above the center of the white dwarf (WD) and adding a moderate amount of angular velocity to the star we follow the evolution of the deflagration using a smoothed particle hydrodynamics code. We find that the results are very dependent on the angle between the rotational axis and the line connecting the initial bubble of burned material with the center of the white dwarf at the moment of the ignition. The impact of rotation is larger for angles close to 90° because the Coriolis force on a floating element of fluid is maximum, and its principal effect is to break the symmetry of the deflagration. Such symmetry breaking weakens the convergence of the nuclear flame at the antipodes of the initial ignition volume, changing the environmental conditions around the convergence region with respect to non-rotating models. These changes seem to disfavor the emergence of a detonation in the compressed volume at the antipodes, thus compromising the viability of the so called gravitational confined detonation mechanism. △ Less

Submitted 21 July, 2015; originally announced July 2015.

Comments: 14 pages, 12 figures, submitted to ApJ

Neutrino-driven winds in the aftermath of a neutron star merger: nucleosynthesis and electromagnetic transients

Authors: Dirk Martin, Albino Perego, Almudena Arcones, Friedrich-Karl Thielemann, Oleg Korobkin, Stephan Rosswog

Abstract: We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynth… ▽ More We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the polar angle. Matter in excess of up to $9 \cdot 10^{-3} M_\odot$ becomes unbound until $\sim 200~{\rm ms}$. Due to electron fractions of $Y_{\rm e} \approx 0.2 - 0.4$ mainly nuclei with mass numbers $A < 130$ are synthesized, complementing the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue band after about $4~{\rm h}$. Furthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta lead to a second peak in the infrared after $3-4~{\rm d}$. △ Less

Submitted 24 September, 2015; v1 submitted 16 June, 2015; originally announced June 2015.

Comments: 15 pages, 18 figures, 2 tables, accepted by ApJ

Report number: INT-PUB-15-021

Journal ref: ApJ 813 (2015) 2

Two-Dimensional Core-Collapse Supernova Simulations with the Isotropic Diffusion Source Approximation for Neutrino Transport

Authors: Kuo-Chuan Pan, Matthias Liebendörfer, Matthias Hempel, Friedrich-Karl Thielemann

Abstract: The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation-hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn (1985) and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped particle and streaming… ▽ More The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation-hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn (1985) and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped particle and streaming particle components. Heavy neutrinos are described by a leakage scheme. Unlike the "ray-by-ray" approach in some other multi-dimensional supernova models, we use cylindrical coordinates and solve the trapped particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We provide an IDSA verification by performing 1D and 2D simulations with 15 and 20 $M_\odot$ progenitors from Woosley et al.~(2007) and discuss the difference of our IDSA results with those existing in the literature. Additionally, we perform Newtonian 1D and 2D simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 $M_\odot$ progenitors from Woosley et al.~(2002) with the HS(DD2) equation of state. General relativistic effects are neglected. We obtain robust explosions with diagnostic energies $E_{\rm dig} \gtrsim 0.1- 0.5$~B for all considered 2D models within approximately $100-300$ milliseconds after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affect the postbounce evolution, e.g.~the ignorance of neutrino-electron scattering during collapse will lead to a stronger explosion. △ Less

Submitted 30 November, 2015; v1 submitted 11 May, 2015; originally announced May 2015.

Comments: 23 pages. Accepted for publication in ApJ

A primordial r-process?

Authors: T. Rauscher, J. H. Applegate, J. J. Cowan, F. -K. Thielemann, M. Wiescher

Abstract: Baryon density inhomogeneities in the early universe can give rise to a floor of heavy elements (up to $A\approx 270$) produced in a primordial r-process with fission cycling. A parameter study with variation of the global baryon to photon ratio $η$, and under inclusion of neutron diffusion effects was performed. New results concerning the dependence of the results on nuclear physics parameters ar… ▽ More Baryon density inhomogeneities in the early universe can give rise to a floor of heavy elements (up to $A\approx 270$) produced in a primordial r-process with fission cycling. A parameter study with variation of the global baryon to photon ratio $η$, and under inclusion of neutron diffusion effects was performed. New results concerning the dependence of the results on nuclear physics parameters are presented. △ Less

Submitted 17 April, 2015; originally announced April 2015.

Comments: 6 pages, 1 table, 2 figures; talk at the European Workshop on Heavy Element Nucleosynthesis, Budapest 1994; appeared in the proceedings, uploaded for easy access. in Proc. European Workshop on Heavy Element Nucleosynthesis, Budapest 1994, eds. E. Somorjai und Z. Fülöp (Institute of Nuclear Research of the Hungarian Academy of Sciences, Debrecen 1994), p. 121

Primordial heavy element production

Authors: T. Rauscher, F. -K. Thielemann

Abstract: A number of possible mechanisms have been suggested to generate density inhomogeneities in the early Universe which could survive until the onset of primordial nucleosynthesis and generate neutron-rich regions. In this work we are not concerned with how the inhomogeneities were generated but we want to focus on the effect of such inhomogeneities on primordial nucleosynthesis. One of the proposed s… ▽ More A number of possible mechanisms have been suggested to generate density inhomogeneities in the early Universe which could survive until the onset of primordial nucleosynthesis and generate neutron-rich regions. In this work we are not concerned with how the inhomogeneities were generated but we want to focus on the effect of such inhomogeneities on primordial nucleosynthesis. One of the proposed signatures of inhomogeneity, the synthesis of very heavy elements by neutron capture in a primordial r-process, was analyzed for varying baryon to photon ratios $η$ and fluctuation length scales $L$. △ Less

Submitted 17 April, 2015; originally announced April 2015.

Comments: 4 pages, 1 figure; talk at the ESO/EIPC Workshop on "Light Element Abundances", Marciana Marina, Isola d'Elba 21-26 May 1994; appeared in the workshop proceedings, uploaded for easy access. in Proceedings of the ESO/EIPC Workshop "The Light Element Abundances", ed. P. Crane (Springer, Berlin Heidelberg 1995), p. 31

Galactic evolution of rapid neutron capture process abundances: the inhomogeneous approach

Authors: B. Wehmeyer, M. Pignatari, F. -K. Thielemann

Abstract: For the origin of heavy r-process elements, different sources have been proposed, e.g., core-collapse supernovae or neutron star mergers. Old metal-poor stars carry the signature of the astrophysical source(s). Among the elements dominantly made by the r-process, europium (Eu) is relatively easy to observe. In this work we simulate the evolution of europium in our galaxy with the inhomogeneous che… ▽ More For the origin of heavy r-process elements, different sources have been proposed, e.g., core-collapse supernovae or neutron star mergers. Old metal-poor stars carry the signature of the astrophysical source(s). Among the elements dominantly made by the r-process, europium (Eu) is relatively easy to observe. In this work we simulate the evolution of europium in our galaxy with the inhomogeneous chemical evolution model 'ICE', and compare our results with spectroscopic observations. We test the most important parameters affecting the chemical evolution of Eu: (a) for neutron star mergers the coalescence time scale of the merger ($t_{\mathrm{coal}}$) and the probability to experience a neutron star merger event after two supernova explosions occurred and formed a double neutron star system ($P_{\mathrm{NSM}}$) and (b) for the sub-class of magneto-rotationally driven supernovae ("Jet-SNe"), their occurrence rate compared to standard supernovae ($P_{\mathrm{Jet-SN}}$). We find that the observed [Eu/Fe] pattern in the galaxy can be reproduced by a combination of neutron star mergers and magneto-rotationally driven supernovae as r-process sources. While neutron star mergers alone seem to set in at too high metallicities, Jet-SNe provide a cure for this deficiency at low metallicities. Furthermore, we confirm that local inhomogeneities can explain the observed large spread in the europium abundances at low metallicities. We also predict the evolution of [O/Fe] to test whether the spread in $α$-elements for inhomogeneous models agrees with observations and whether this provides constraints on supernova explosion models and their nucleosynthesis. △ Less

Submitted 16 September, 2015; v1 submitted 30 January, 2015; originally announced January 2015.

Comments: 11 Figures, Submitted to MNRAS

Journal ref: MNRAS, 452, 1970W (2015)

The r-process nucleosynthesis in the various jet-like explosions of magnetorotational core-collapse supernovae

Authors: Nobuya Nishimura, Tomoya Takiwaki, Friedrich-Karl Thielemann

Abstract: The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., pr… ▽ More The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena. △ Less

Submitted 31 August, 2015; v1 submitted 26 January, 2015; originally announced January 2015.

Comments: 25 pages, 17 figures, 6 tables, ApJ in press, numerical data are available at http://www.astro.keele.ac.uk/~nobuya/mrsn and https://github.com/nnobuya/mrsn

Journal ref: ApJ, 810, 109 (2015)

Pushing 1D CCSNe to explosions: model and SN 1987A

Authors: A. Perego, M. Hempel, C. Fröhlich, K. Ebinger, M. Eichler, J. Casanova, M. Liebendoerfer, F. -K. Thielemann

Abstract: We report on a method, PUSH, for triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH such that the observables of SN1987A are reproduced. Our simulations are based on the general relativistic hydrodynamics code AGILE combined with the detailed neutrino transport scheme IDSA for electron neutrinos and ALS f… ▽ More We report on a method, PUSH, for triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH such that the observables of SN1987A are reproduced. Our simulations are based on the general relativistic hydrodynamics code AGILE combined with the detailed neutrino transport scheme IDSA for electron neutrinos and ALS for the muon and tau neutrinos. To trigger explosions in the otherwise non-exploding simulations, we rely on the neutrino-driven mechanism. The PUSH method locally increases the energy deposition in the gain region through energy deposition by the heavy neutrino flavors. Our setup allows us to model the explosion for several seconds after core bounce. We explore the progenitor range 18-21M$_{\odot}$. Our studies reveal a distinction between high compactness (HC) and low compactness (LC) progenitor models, where LC models tend to explore earlier, with a lower explosion energy, and with a lower remnant mass. HC models are needed to obtain explosion energies around 1 Bethe, as observed for SN1987A. However, all the models with sufficiently high explosion energy overproduce $^{56}$Ni. We conclude that fallback is needed to reproduce the observed nucleosynthesis yields. The nucleosynthesis yields of $^{57-58}$Ni depend sensitively on the electron fraction and on the location of the mass cut with respect to the initial shell structure of the progenitor star. We identify a progenitor and a suitable set of PUSH parameters that fit the explosion properties of SN1987A when assuming 0.1M$_{\odot}$ of fallback. We predict a neutron star with a gravitational mass of 1.50M$_{\odot}$. We find correlations between explosion properties and the compactness of the progenitor model in the explored progenitors. However, a more complete analysis will require the exploration of a larger set of progenitors with PUSH. △ Less

Submitted 22 May, 2015; v1 submitted 12 January, 2015; originally announced January 2015.

Comments: revised version as accepted by ApJ (results unchanged, text modified for clarification, a few references added); 26 pages, 20 figures

Journal ref: Astrophys. J. 806, 275 (2015)