Source author record

F. Herwig

F. Herwig appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

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

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

preprint2019arXiv

Experimentally well-constrained masses of $^{27}$P and $^{27}$S: Implications for studies of explosive binary systems

The mass of $^{27}$P was predicted to impact the X-ray burst (XRB) model predictions of burst light curves and the composition of the burst ashes. To address the uncertainties and inconsistencies in the reported $^{27}$P masses in literature, a wealth of information has been extracted from the $β$-decay spectroscopy of the drip-line nucleus $^{27}$S. We determine the most precise mass excess of $^{27}$P to date to be $-659(9)$~keV, which is 63~keV (2.3$σ$) higher than the AME2016 recommended value of $-722(26)$~keV. The experimentally unknown mass excess of $^{27}$S was estimated to be 17030(400)~keV in AME2016, and we constrain this mass to be 17678(77)~keV based on the measured $β$-delayed two-proton energy. In the temperature region of $(0.06-0.3)$~GK, the $^{26}$Si$(p,γ)^{27}$P reaction rate determined in this work is significantly lower than the rate recommended in the reaction rate libraries, up to two orders of magnitude around 0.1~GK. The impact of these newly determined masses and well-constrained rate on the modeling of the explosive astrophysical scenarios has been explored by hydrodynamic nova and post-processing XRB models. No substantial change was found in the nova contribution to the synthesis of galactic $^{26}$Al or in the XRB energy generation rate, but we found that the calculated abundances of $^{26}$Al and $^{26}$Si at the last stage of XRB are increased by a factor of 2.4. We also conclude that $^{27}$S is not a significant waiting point in the rapid proton capture process.

preprint2016arXiv

NuGrid stellar data set I. Stellar yields from H to Bi for stars with metallicities Z = 0.02 and Z = 0.01

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

preprint2015arXiv

Carbon-rich presolar grains from massive stars. Subsolar 12C/13C and 14N/15N ratios and the mystery of 15N

Carbon-rich grains with isotopic anomalies compared to the Sun are found in primitive meteorites. They were made by stars, and carry the original stellar nucleosynthesis signature. Silicon carbide grains of Type X and C, and low-density graphites condensed in the ejecta of core-collapse supernovae. We present a new set of models for the explosive He shell and compare them with the grains showing 12C/13C and 14N/15N ratios lower than solar. In the stellar progenitor H was ingested into the He shell and not fully destroyed before the explosion. Different explosion energies and H concentrations are considered. If the SN shock hits the He-shell region with some H still present, the models can reproduce the C and N isotopic signatures in C-rich grains. Hot-CNO cycle isotopic signatures are obtained, including a large production of 13C and 15N. The short-lived radionuclides 22Na and 26Al are increased by orders of magnitude. The production of radiogenic 22Ne from the decay of 22Na in the He shell might solve the puzzle of the Ne-E(L) component in low-density graphite grains. This scenario is attractive for the SiC grains of type AB with 14N/15N ratios lower than solar, and provides an alternative solution for SiC grains originally classified as nova grains. Finally, this process may contribute to the production of 14N and 15N in the Galaxy, helping to produce the 14N/15N ratio in the solar system.

preprint2015arXiv

H ingestion into He-burning convection zones in super-AGB stellar models as a potential site for intermediate neutron-density nucleosynthesis

We investigate the evolution of super-AGB thermal pulse (TP) stars for a range of metallicities (Z) and explore the effect of convective boundary mixing (CBM). With decreasing metallicity and evolution along the TP phase, the He-shell flash and the third dredge-up (TDU) occur closer together in time. After some time (depending upon the CBM parameterisation), efficient TDU begins while the pulse-driven convection zone (PDCZ) is still present, causing a convective exchange of material between the PDCZ and the convective envelope. This results in the ingestion of protons into the convective He-burning pulse. Even small amounts of CBM encourage the interaction of the convection zones leading to transport of protons from the convective envelope into the He layer. H-burning luminosities exceed $10^9$ (in some cases $10^{10}$) $\mathrm{L}_\odot$. We also calculate models of dredge-out in the most massive super-AGB stars and show that the dredge-out phenomenon is another likely site of convective-reactive H-$^{12}$C combustion. We discuss the substantial uncertainties of stellar evolution models under these conditions. Nevertheless, the simulations suggest that in the convective-reactive H-combustion regime of H ingestion the star may encounter conditions for the intermediate neutron capture process (i process). We speculate that some CEMP-s/r stars could originate in i-process conditions in the H-ingestion phases of low-Z SAGB stars. This scenario would however suggest a very low electron-capture supernova rate from super-AGB stars. We also simulate potential outbursts triggered by such H-ingestion events, present their light curves and briefly discuss their transient properties.

preprint2015arXiv

The i-process and CEMP-r/s stars

We investigate whether the anomalous elemental abundance patterns in some of the C-enhanced metal-poor-s+r (CEMP-r/s) stars are consistent with predictions of nucleosynthesis yields from the i-process, a neutron-capture regime at neutron densities intermediate between those typical for the slow (s) and rapid (r) processes. Conditions necessary for the i-process are expected to be met at multiple stellar sites, such as the He-core and He-shell flashes in low-metallicity low-mass stars, super-AGB and post-AGB stars, as well as low-metallicity massive stars. We have found that single-exposure one-zone simulations of the i-process reproduce the abundance patterns in some of the CEMP-r/s stars much better than the model that assumes a superposition of yields from s- and r-process sources. Our previous study of nuclear data uncertainties relevant to the i-process revealed that they could have a significant impact on the i-process yields obtained in our idealized one-zone calculations, leading, for example, to ~0.7dex uncertainty in our predicted [Ba/La] ratio. Recent 3D hydrodynamic simulations of convection driven by a He-shell flash in post-AGB Sakurai's object have discovered a new mode of non-radial instabilities: the Global Oscillation of Shell H-ingestion. This has demonstrated that spherically symmetric stellar evolution simulations cannot be used to accurately model physical conditions for the i-process.

preprint2014arXiv

Hybrid C-O-Ne white dwarfs as progenitors of type Ia supernovae: dependence on Urca process and mixing assumptions

When carbon is ignited off-centre in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the centre. Whether the C flame will actually be able to reach the centre depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the centre. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations, the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-centre C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the centre. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

preprint2014arXiv

New insights on Ba over-abundance in open clusters. Evidence for the intermediate neutron-capture process at play?

Recently an increasing number of studies were devoted to measure the abundances of neutron-capture elements heavier than iron in stars belonging to Galactic Open Clusters (OCs). OCs span a sizeable range in metallicity -0.6<[Fe/H]<+0.4), and they show abundances of light elements similar to disk stars of the same age. A different pattern is observed for heavy elements. A large scatter is observed for Ba, with most OCs showing [Ba/Fe] and [Ba/La] overabundant with respect to the Sun. The origin of this overabundance is not clearly understood. With the goal of providing new observational insights we determined radial velocities, atmospheric parameters and chemical composition of 27 giant stars members of five OCs: Cr 110, Cr 261, NGC 2477, NGC 2506 and NGC 5822. We used high-resolution spectra obtained with the UVES spectrograph at ESO Paranal. We perform a detailed spectroscopic analysis of these stars to measure the abundance of up to 22 elements per star. We study the dependence of element abundance on metallicity and age with unprecedented detail, complementing our analysis with data culled from the literature. We confirm the trend of Ba overabundance in OCs, and show its large dispersion for clusters younger than ~4 Gyr. Finally, the implications of our results for stellar nucleosynthesis are discussed. We show in this work that the Ba enrichment compared to other neutron-capture elements in OCs cannot be explained by the contributions from the slow neutron-capture process and the rapid neutron-capture process. Instead, we argue that this anomalous signature can be explained by assuming an additional contribution by the intermediate neutron-capture process.

preprint2014arXiv

The Very Fast Evolution of the VLTP Object V4334 Sgr

V4334 Sgr (Sakurai's object) is an enigmatic evolved star that underwent a very late thermal pulse a few years before its discovery in 1996. It ejected a new, hydrogen-deficient nebula in the process. Emission lines from the newly ejected gas were first discovered in 1998 (He I 1083 nm) and 2001 (optical). We have monitored the optical emission spectrum since. From 2001 through 2007 the optical spectrum showed an exponential decline in flux, consistent with a shock that occurred around 1998 and started cooling soon after that. In this paper we show that since 2008 the line fluxes have been continuously rising again. Our preliminary interpretation is that this emission comes from a region close to the central star, and is excited by a second shock. This shock may have been induced by an increase in the stellar mass loss and wind velocity associated with a rise in the stellar temperature.

preprint2013arXiv

Production of carbon-rich presolar grains from massive stars

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($α$,$γ$)$^{16}$O reaction relative to the rest of the $α$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

preprint2013arXiv

Silicon carbide grains of type C provide evidence for the production of the unstable isotope $^{32}$Si in supernovae

Carbon-rich grains are observed to condense in the ejecta of recent core-collapse supernovae, within a year after the explosion. Silicon carbide grains of type X are C-rich grains with isotpic signatures of explosive supernova nucleosynthesis have been found in primitive meteorites. Much rarer silicon carbide grains of type C are a special sub-group of SiC grains from supernovae. They show peculiar abundance signatures for Si and S, isotopically heavy Si and isotopically light S, which appear to to be in disagreement with model predictions. We propose that C grains are formed mostly from C-rich stellar material exposed to lower SN shock temperatures than the more common type X grains. In this scenario, extreme $^{32}$S enrichments observed in C grains may be explained by the presence of short-lived $^{32}$Si ($τ$$_{1/2}$ = 153 years) in the ejecta, produced by neutron capture processes starting from the stable Si isotopes. No mixing from deeper Si-rich material and/or fractionation of Si from S due to molecular chemistry is needed to explain the $^{32}$S enrichments. The abundance of $^{32}$Si in the grains can provide constraints on the neutron density reached during the supernova explosion in the C-rich He shell material. The impact of the large uncertainty of the neutron capture cross sections in the $^{32}$Si region is discussed.

preprint2013arXiv

Systematic and correlated nuclear uncertainties in the i-process at the neutron shell closure N = 82

Nuclear astrophysics simulations aiming to study the origin of the elements in stars require a multitude of nuclear physics input. Both systematic model dependent and statistically correlated uncertainties need to be considered. An application where realistic uncertainty assessments are especially important is the intermediate neutron capture process or i process: a neutron capture regime with neutron densities intermediate between the slow and rapid processes. Accordingly, the main network flux proceeds on the neutron-rich unstable isotopes up to 4-5 species off the valley of stability. The i process has been clearly identified to be active in post-AGB stars during the Very Late Thermal Pulse H-ingestion event, and a recent work infers about its important role in early generations of stars. Here we demonstrate the effect of propagating systematic nuclear uncertainties from different theoretical models to final abundances for a region around the 2$^\mathrm{nd}$ peak at $A-Z=80$ for elemental ratio predictions involving Ba, La and Eu in i-process conditions. These elements are used to distinguish different n-capture contributions observed in low-metallicity stars. For the simple 1-zone model adopted here, predictions vary as much as a factor of 22 in on possible observational plane ([La/Eu] vs. [Ba/La]). To consider statistically correlated uncertainties, we similarly perform a full nuclear physics uncertainty study within a given Hauser-Feshbach model and demonstrate the role of correlations on the final stellar abundance uncertainties. We show that in i-process conditions the main result of neglecting correlations is to underestimate the impact of nuclear uncertainties on the final nucleosynthesis yields by as much as two orders of magnitude. (continued in document)

preprint2012arXiv

The 12C + 12C reaction and the impact on nucleosynthesis in massive stars

Despite much effort in the past decades, the C-burning reaction rate is uncertain by several orders of magnitude, and the relative strength between the different channels 12C(12C,alpha)20Ne, 12C(12C,p)23Na and 12C(12C,n)23Mg is poorly determined. Additionally, in C-burning conditions a high 12C+12C rate may lead to lower central C-burning temperatures and to 13C(alpha,n)16O emerging as a more dominant neutron source than 22Ne(alpha,n)25Mg, increasing significantly the s-process production. This is due to the rapid decrease of the 13N(gamma,p)12C with decreasing temperature, causing the 13C production via 13N(beta+)13C. Presented here is the impact of the 12C+12C reaction uncertainties on the s-process and on explosive p-process nucleosynthesis in massive stars, including also fast rotating massive stars at low metallicity. Using various 12C+12C rates, in particular an upper and lower rate limit of ~ 50000 higher and ~ 20 lower than the standard rate at 5*10^8 K, five 25 Msun stellar models are calculated. The enhanced s-process signature due to 13C(alpha,n)16O activation is considered, taking into account the impact of the uncertainty of all three C-burning reaction branches. Consequently, we show that the p-process abundances have an average production factor increased up to about a factor of 8 compared to the standard case, efficiently producing the elusive Mo and Ru proton-rich isotopes. We also show that an s-process being driven by 13C(alpha,n)16O is a secondary process, even though the abundance of 13C does not depend on the initial metal content. Finally, implications for the Sr-peak elements inventory in the Solar System and at low metallicity are discussed.

preprint2012arXiv

The effect of 12C + 12C rate uncertainties on the evolution and nucleosynthesis of massive stars

[Shortened] The 12C + 12C fusion reaction has been the subject of considerable experimental efforts to constrain uncertainties at temperatures relevant for stellar nucleosynthesis. In order to investigate the effect of an enhanced carbon burning rate on massive star structure and nucleosynthesis, new stellar evolution models and their yields are presented exploring the impact of three different 12C + 12C reaction rates. Non-rotating stellar models were generated using the Geneva Stellar Evolution Code and were later post-processed with the NuGrid Multi-zone Post-Processing Network tool. The enhanced rate causes core carbon burning to be ignited more promptly and at lower temperature. This reduces the neutrino losses, which increases the core carbon burning lifetime. An increased carbon burning rate also increases the upper initial mass limit for which a star exhibits a convective carbon core. Carbon shell burning is also affected, with fewer convective-shell episodes and convection zones that tend to be larger in mass. Consequently, the chance of an overlap between the ashes of carbon core burning and the following carbon shell convection zones is increased, which can cause a portion of the ashes of carbon core burning to be included in the carbon shell. Therefore, during the supernova explosion, the ejecta will be enriched by s-process nuclides synthesized from the carbon core s process. The yields were used to estimate the weak s-process component in order to compare with the solar system abundance distribution. The enhanced rate models were found to produce a significant proportion of Kr, Sr, Y, Zr, Mo, Ru, Pd and Cd in the weak component, which is primarily the signature of the carbon-core s process. Consequently, it is shown that the production of isotopes in the Kr-Sr region can be used to constrain the 12C + 12C rate using the current branching ratio for a- and p-exit channels.

preprint2006arXiv

Detection of a Far-Infrared Bow-Shock Nebula Around R Hya: the First MIRIAD Results

We present the first results of the MIRIAD (MIPS [Multiband Imaging Photometer for Spitzer] Infra-Red Imaging of AGB [asymptotic giant branch] Dustshells) project using the Spitzer Space Telescope. The primary aim of the project is to probe the material distribution in the extended circumstellar envelopes (CSE) of evolved stars and recover the fossil record of their mass loss history. Hence, we must map the whole of the CSEs plus the surrounding sky for background subtraction, while avoiding the central star that is brighter than the detector saturation limit. With our unique mapping strategy, we have achieved better than one MJy/sr sensitivity in three hours of integration and successfully detected a faint (< 5 MJy/sr), extended (~400 arcsec) far-infrared nebula around the AGB star R Hya. Based on the parabolic structure of the nebula, the direction of the space motion of the star with respect to the nebula shape, and the presence of extended H alpha emission co-spatial to the nebula, we suggest that the detected far-IR nebula is due to a bow shock at the interface of the interstellar medium and the AGB wind of this moving star. This is the first detection of the stellar-wind bow-shock interaction for an AGB star and exemplifies the potential of Spitzer as a tool to examine the detailed structure of extended far-IR nebulae around bright central sources. \