Source author record

Falk Herwig

Falk 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

47works
10topics
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

47 published item(s)

preprint2026arXiv

Propagating Uncertainties from Nuclear Physics to Gamma-rays in Core Collapse Supernovae

Nuclear yields are powerful probes of supernova explosions, their engines and their progenitors. In addition, as we improve our understanding of these explosions, we can use nuclear yields to probe dense matter and neutrino physics, both of which play a critical role in the central supernova engine. Especially with upcoming gamma-ray detectors that can directly detect radioactive isotopes out to increasing distances from gamma-rays emitted during their decay, nuclear yields have the potential to provide some of the most direct probes of supernova engines and stellar burning. To utilize these probes, we must understand and limit the uncertainties in their production. Uncertainties in the nuclear physics can be minimized by combining both laboratory experiments and nuclear theory. Similarly, astrophysical uncertainties caused by simplified explosion trajectories can be minimized by higher-fidelity stellar-evolution and supernova-engine models. This paper reviews the physics and astrophysics uncertainties in modeling nucleosynthetic yields, identifying the key areas of study needed to maximize the potential of supernova yields as probes of astrophysical transients and dense-matter physics.

preprint2026arXiv

The NuGrid AGB Evolution and Nucleosynthesis Data Set

Asymptotic Giant Branch (AGB) stars play a key role in the chemical evolution of galaxies. These stars are the fundamental stellar site for the production of light elements such as C, N and F, and half of the elements heavier than Fe via the slow neutron capture process (s-process). Hence, detailed computational models of AGB stars' evolution and nucleosynthesis are essential for galactic chemical evolution. In this work, we discuss the progress in updating the NuGrid data set of AGB stellar models and abundance yields. All stellar models have been computed using the MESA stellar evolution code, coupled with the post-processing mppnp code to calculate the full nucleosynthesis. The final data set will include the initial masses Mini/Msun = 1, 1.65, 2, 3, 4, 5, 6 and 7 for initial metallicities Z = 0.0001, 0.001, 0.006, 0.01, 0.02 and 0.03. Observed s-process abundances on the surfaces of evolved stars as well as the typical light elements in the composition of H-deficient post-AGB stars are reproduced. A key short-term goal is to complete and expand the AGB stars data set for the full metallicity range. Chemical yield tables are provided for the available models.

preprint2021arXiv

3D1D hydro-nucleosynthesis simulations. I. Advective-reactive post-processing method and its application to H ingestion into He-shell flash convection in rapidly accreting white dwarfs

We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the $^{12}$C(p,$γ$)$^{13}$N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the $X_{\mathrm{H}}$ profiles between the up- and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the up- and downstream. For example, $^{89}$Kr and $^{90}$Kr with half-lives of 3.18 min and 32.3 s differ by a factor 2-10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050.

preprint2020arXiv

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

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.

preprint2019arXiv

Chromium Nucleosynthesis and Silicon-Carbon Shell Mergers in Massive Stars

We analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] $\approx 0$ reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 $M_\odot$ model at metallicity $Z=0.01$, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including $^{51}$V and $^{52}$Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 $M_\odot$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 $M_\odot$ model at $Z=0.01$ in the yields interpolation scheme of our GCE model for stars in between 15 and 25 $M_\odot$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] $\approx$ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si-C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the pre-supernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si-shell could launch a convective-reactive global oscillation, if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.

preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

preprint2019arXiv

The impact of (n,$γ$) reaction rate uncertainties on the predicted abundances of i-process elements with $32\leq Z\leq 48$ in the metal-poor star HD94028

Several anomalous elemental abundance ratios have been observed in the metal-poor star HD94028. We assume that its high [As/Ge] ratio is a product of a weak intermediate (i) neutron-capture process. Given that observational errors are usually smaller than predicted nuclear physics uncertainties, we have first set up a benchmark one-zone i-process nucleosynthesis simulation results of which provide the best fit to the observed abundances. We have then performed Monte Carlo simulations in which 113 relevant (n,$γ$) reaction rates of unstable species were randomly varied within Hauser-Feshbach model uncertainty ranges for each reaction to estimate the impact on the predicted stellar abundances. One of the interesting results of these simulations is a double-peaked distribution of the As abundance, which is caused by the variation of the $^{75}$Ga (n,$γ$) cross section. This variation strongly anti-correlates with the predicted As abundance, confirming the necessity for improved theoretical or experimental bounds on this cross section. The $^{66}$Ni (n,$γ$) reaction is found to behave as a major bottleneck for the i-process nucleosynthesis. Our analysis finds the Pearson product-moment correlation coefficient $r_\mathrm{P} > 0.2$ for all of the i-process elements with $32 \leq Z \leq 42$, with significant changes in their predicted abundances showing up when the rate of this reaction is reduced to its theoretically constrained lower bound. Our results are applicable to any other stellar nucleosynthesis site with the similar i-process conditions, such as Sakurai's object (V4334 Sagittarii) or rapidly-accreting white dwarfs.

preprint2018arXiv

H-He Shell Interactions and Nucleosynthesis in Massive Population III Stars

We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code.

preprint2016arXiv

Idealised hydrodynamic simulations of turbulent oxygen-burning shell convection in 4π geometry

This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4π$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The $1536^3$ simulation exhibits an entrainment rate at the upper CB of $1.33\times10^{-6}~M_\odot~\mathrm{s}^{-1}$. The $768^3$ simulation with the same heating rate agrees within 17 per cent. The entrainment rate at the upper convective boundary is found to scale linearly with the driving luminosity and with the cube of the shear velocity at the upper boundary, while the radial RMS fluid velocity scales with the cube root of the driving luminosity, as expected. The mixing is analysed in a 1D diffusion framework, resulting in a simple model for CB mixing. The analysis confirms previous findings that limiting the MLT mixing length to the distance to the CB in 1D simulations better represents the spherically-averaged radial velocity profiles from the 3D simulations and provides an improved determination of the reference diffusion coefficient $D_0$ for the exponential diffusion CB mixing model in 1D. From the 3D simulation data we adopt as the convective boundary the location of the maximum gradient in the horizontal velocity component which has $2σ$ spatial fluctuations of $\approx0.17 H_P$ . The exponentially decaying diffusion CB mixing model with $f = 0.03$ reproduces the spherically-averaged 3D abundance profiles.

preprint2016arXiv

Mass and Metallicity Requirement in Stellar Models for Galactic Chemical Evolution Applications

We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We ran several simulations where we sampled subsets of stellar models to explore the impact of different grid resolutions. Stellar yields from low- and intermediate-mass stars and from Type Ia supernovae have been included in our simulations, but with a fixed grid resolution. We compared our results with the stellar abundances observed in the Milky Way for O, Na, Mg, Si, Ca, Ti, and Mn. Our results suggest that the range of metallicity considered is more important than the number of metallicities within that range, which only affects our numerical predictions by about 0.1 dex. We found that our predictions at [Fe/H] < -2 are very sensitive to the metallicity range and the mass sampling used for the lowest metallicity included in the set of yields. Variations between results can be as high as 0.8 dex, for any remnant mass prescription. At higher [Fe/H], we found that the required number of masses depends on the element of interest and on the remnant mass prescription. With a monotonic remnant mass prescription where every model explodes as a core-collapse supernova, the mass resolution induces variations of 0.2 dex on average. But with a remnant mass prescription that includes islands of non-explodability, the mass resolution can cause variations of about 0.2 to 0.7 dex depending on the choice of metallicity range. With such a prescription, explosive or non-explosive models can be missed if not enough masses are selected, resulting in over- or under-estimations of the mass ejected by massive stars.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2016arXiv

The Diverse Origins of Neutron-Capture Elements in the Metal-Poor Star HD 94028: Possible Detection of Products of i-process Nucleosynthesis

We present a detailed analysis of the composition and nucleosynthetic origins of the heavy elements in the metal-poor ([Fe/H]=-1.62+/-0.09) star HD94028. Previous studies revealed that this star is mildly enhanced in elements produced by the slow neutron-capture process (s-process; e.g., [Pb/Fe]=+0.79+/-0.32) and rapid neutron-capture process (r-process; e.g., [Eu/Fe]=+0.22+/-0.12), including unusually large molybdenum ([Mo/Fe]=+0.97+/-0.16) and ruthenium ([Ru/Fe]=+0.69+/-0.17) enhancements. However, this star is not enhanced in carbon ([C/Fe]=-0.06+/-0.19). We analyze an archival near-ultraviolet spectrum of HD94028, collected using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope, and other archival optical spectra collected from ground-based telescopes. We report abundances or upper limits derived from 64 species of 56 elements. We compare these observations with s-process yields from low-metallicity AGB evolution and nucleosynthesis models. No combination of s- and r-process patterns can adequately reproduce the observed abundances, including the super-solar [As/Ge] ratio (+0.99+/-0.23) and the enhanced [Mo/Fe] and [Ru/Fe] ratios. We can fit these features when including an additional contribution from the intermediate neutron-capture process (i process), which perhaps operated by the ingestion of H in He-burning convective regions in massive stars, super-AGB stars, or low-mass AGB stars. Currently, only the i process appears capable of consistently producing the super-solar [As/Ge] ratios and ratios among neighboring heavy elements found in HD94028. Other metal-poor stars also show enhanced [As/Ge] ratios, hinting that operation of the i process may have been common in the early Galaxy.

preprint2016arXiv

The effect of a wider initial separation on common envelope binary interaction simulations

We present hydrodynamic simulations of the common envelope binary interaction between a giant star and a compact companion carried out with the adaptive mesh refinement code ENZO and the smooth particle hydrodynamics code PHANTOM. These simulations mimic the parameters of one of the simulations by Passy et al., but assess the impact of a larger, more realistic initial orbital separation on the simulation outcome. We conclude that for both codes the post-common envelope separation is somewhat larger and the amount of unbound mass slightly greater when the initial separation is wide enough that the giant does not yet overflow or just overflows its Roche lobe. PHANTOM has been adapted to the common envelope problem here for the first time and a full comparison with ENZO is presented, including an investigation of convergence as well as energy and angular momentum conservation. We also set our simulations in the context of past simulations. This comparison reveals that it is the expansion of the giant before rapid in-spiral and not spinning up of the star that causes a larger final separation. We also suggest that the large range in unbound mass for different simulations is difficult to explain and may have something to do with simulations that are not fully converged.

preprint2016arXiv

Type Ia Supernova Explosions from Hybrid Carbon-Oxygen-Neon White Dwarf Progenitors

Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C-O core inside an O-Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration to detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced with the MESA stellar evolution code and include the effects of the Urca process, and we map the progenitor model to the FLASH grid. We performed a suite of DDT simulations over a range of ignition conditions consistent with the progenitor's thermal and convective structure assuming multiple ignition points. To compare the results from these hybrid WD stars to previous results from C-O white dwarfs, we construct a set of C-O WD models with similar properties and similarly simulate a suite of explosions. We find that despite significant variability within each suite, trends distinguishing the explosions are apparent in their $^{56}$Ni yields and the kinetic properties of the ejecta. We comment on the feasibility of these explosions as the source of some classes of observed subluminous events.

preprint2016arXiv

Uncertainties in Galactic Chemical Evolution Models

We use a simple one-zone galactic chemical evolution model to quantify the uncertainties generated by the input parameters in numerical predictions, for a galaxy with properties similar to those of the Milky Way. We compiled several studies from the literature to gather the current constraints for our simulations regarding the typical value and uncertainty of seven basic parameters, which are: the lower and upper mass limits of the stellar initial mass function (IMF), the slope of the high-mass end of the stellar IMF, the slope of the delay-time distribution function of Type Ia supernovae (SNe Ia), the number of SNeIa per Msun formed, the total stellar mass formed, and the final mass of gas. We derived a probability distribution function to express the range of likely values for every parameter, which were then included in a Monte Carlo code to run several hundred simulations with randomly selected input parameters. This approach enables us to analyze the predicted chemical evolution of 16 elements in a statistical way by identifying the most probable solutions along with their 68% and 95% confidence levels. Our results show that the overall uncertainties are shaped by several input parameters that individually contribute at different metallicities, and thus at different galactic ages. The level of uncertainty then depends on the metallicity and is different from one element to another. Among the seven input parameters considered in this work, the slope of the IMF and the number of SNe Ia are currently the two main sources of uncertainty.

preprint2016arXiv

White Paper on Nuclear Astrophysics

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

preprint2015arXiv

The primordial and evolutionary abundance variations in globular-cluster stars: a problem with two unknowns

We demonstrate that among the potential sources of the primordial abundance variations of the proton-capture elements in globular-cluster stars proposed so far, such as the hot-bottom burning in massive AGB stars and H burning in the convective cores of supermassive and fast-rotating massive MS stars, only the supermassive MS stars with M > 10,000 Msun can explain all the observed abundance correlations without any fine-tuning of model parameters. We use our assumed chemical composition for the pristine gas in M13 (NGC6205) and its mixtures with 50% and 90% of the material partially processed in H burning in the 60,000 Msun MS model star as the initial compositions for the normal, intermediate and extreme populations of low-mass stars in this globular cluster, as suggested by its O-Na anti-correlation. We evolve these stars from the zero-age MS to the RGB tip with the thermohaline and parametric prescriptions for the RGB extra mixing. We find that the 3He-driven thermohaline convection cannot explain the evolutionary decline of [C/Fe] in M13 RGB stars, which, on the other hand, is well reproduced with the universal values for the mixing depth and rate calibrated using the observed decrease of [C/Fe] with MV in the globular cluster NGC5466 that does not have the primordial abundance variations.

preprint2014arXiv

Barium Isotopic Composition of Mainstream Silicon Carbides from Murchison: Constraints for s-Process Nucleosynthesis in AGB Stars

We present barium, carbon, and silicon isotopic compositions of 38 acid-cleaned presolar SiC grains from Murchison. Comparison with previous data shows that acid washing is highly effective in removing barium contamination. Strong depletions in $δ$($^{138}$Ba/$^{136}$Ba) values are found, down to $-$400 permil, which can only be modeled with a flatter $^{13}$C profile within the $^{13}$C pocket than is normally used. The dependence of $δ$($^{138}$Ba/$^{136}$Ba) predictions on the distribution of $^{13}$C within the pocket in AGB models allows us to probe the $^{13}$C profile within the $^{13}$C pocket and the pocket mass in asymptotic giant branch (AGB) stars. In addition, we provide constraints on the $^{22}$Ne$(α,n)^{25}$Mg rate in the stellar temperature regime relevant to AGB stars, based on $δ$($^{134}$Ba/$^{136}$Ba) values of mainstream grains. We found two nominally mainstream grains with strongly negative $δ$($^{134}$Ba/$^{136}$Ba) values that cannot be explained by any of the current AGB model calculations. Instead, such negative values are consistent with the intermediate neutron capture process ($i$-process), which is activated by the Very Late Thermal Pulse (VLTP) during the post-AGB phase and characterized by a neutron density much higher than the $s$-process. These two grains may have condensed around post-AGB stars. Finally, we report abundances of two $p$-process isotopes, $^{130}$Ba and $^{132}$Ba, in single SiC grains. These isotopes are destroyed in the $s$-process in AGB stars. By comparing their abundances with respect to that of $^{135}$Ba, we conclude that there is no measurable decay of $^{135}$Cs ($t_{1/2}$= 2.3 Ma) to $^{135}$Ba in individual SiC grains, indicating condensation of barium, but not cesium into SiC grains before $^{135}$Cs decayed.

preprint2014arXiv

Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

preprint2014arXiv

Global non-spherical oscillations in 3-D $4π$ simulations of the H-ingestion flash

We performed 3-D simulations of proton-rich material entrainment into \czw-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurai's object. Observations of the transient nature and anomalous abundance features are available to validate our method and assumptions, with the aim to apply them to very low metallicity stars in the future. We include nuclear energy feedback from H burning and cover the full $4π$ geometry of the shell. Runs on $768^3$ and $1536^3$ grids agree well with each other and have been followed for $1500\mathrm{min}$ and $1200\mathrm{min}$. After a $850\mathrm{min}$ long quiescent entrainment phase the simulations enter into a global non-spherical oscillation that is launched and sustained by individual ignition events of H-rich fluid pockets. Fast circumferential flows collide at the antipode and cause the formation and localized ignition of the next H-overabundant pocket. The cycle repeats for more than a dozen times while its amplitude decreases. During the global oscillation the entrainment rate increases temporarily by a factor $\approx 100$. Entrained entropy quenches convective motions in the upper layer until the burning of entrained H establishes a separate convection zone. The lower-resolution run hints at the possibility that another global oscillation, perhaps even more violent will follow. The location of the H-burning convection zone agrees with a 1-D model in which the mixing efficiency is calibrated to reproduce the light curve. The simulations have been performed at the NSF Blue Waters supercomputer at NCSA.

preprint2014arXiv

Hybrid C-O-Ne White Dwarfs as Progenitors of Diverse SNe Ia

When carbon is ignited off-center in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the center. Whether the C flame will actually be able to reach the center 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 center. 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-center 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 center. 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

Hydrodynamic simulations of H entrainment at the top of He-shell flash convection

We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4π$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurai's object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 \pm 1.48 \times 10^{-13}M_\odot \mathrm{/s}$.

preprint2014arXiv

MESA and NuGrid simulations of classical novae: CO and ONe nova nucleosynthesis

Classical novae are the result of thermonuclear flashes of hydrogen accreted by CO or ONe white dwarfs, leading eventually to the dynamic ejection of the surface layers. These are observationally known to be enriched in heavy elements, such as C, O and Ne that must originate in layers below the H-flash convection zone. Building on our previous work, we now present stellar evolution simulations of ONe novae and provide a comprehensive comparison of our models with published ones. Some of our models include exponential convective boundary mixing to account for the observed enrichment of the nova ejecta even when accreted material has a solar abundance distribution. Our models produce maximum temperature evolution profiles and nucleosynthesis yields in good agreement with models that generate enriched ejecta by assuming that the accreted material was pre-mixed. We confirm for ONe novae the result we reported previously, i.e.\ we found that $^3$He could be produced {\it in situ} in solar-composition envelopes accreted with slow rates ($\dot{M} < 10^{-10}\,M_\odot/\mbox{yr}$) by cold ($T_{\rm WD} < 10^7$ K) CO WDs, and that convection was triggered by $^3$He burning before the nova outburst in that case. In addition, we now find that the interplay between the $^3$He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g.\ $\dot{M} = 10^{-10}\,M_\odot/\mbox{yr}$, by the $1.15\,M_\odot$ ONe WD with a relatively high initial central temperature, e.g.\ $T_{\rm WD} = 15\times 10^6$ K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. (Abridged)

preprint2014arXiv

The dependence of the evolution of SN type Ia progenitors on the C burning rate uncertainty and parameters of convective boundary mixing

Evolution of a supernova type Ia progenitor requires formation of a CO white dwarf, which implies a dependence on the C burning rate (CBR). It can also be affected by the recently identified possibility of C flame quenching by convective boundary mixing. We present first results of our study of the combined effect of these two potential sources of uncertainty on the SN Ia progenitor evolution. We consider the possibility that the CBR is higher than its currently recommended value by as much as a factor of 1000 if unidentified resonances are important, or that it is significantly lower because of the hindrance effect. For stellar models that assume the Schwarzschild boundary for convection, the maximum initial mass for the formation of CO WDs increases from M_i ~ 5.5 Msun for the CBR factor of 1000 to M_i > 7.0 Msun for the CBR factor of 0.01. For C-flame quenching models, hybrid C-O-Ne WDs form for a range of initial mass of Delta M_i ~ 1 Msun, which increases a fraction of stars that form WDs capable of igniting C in a thermonuclear runaway. The most extreme case is found for the CBR factor of 0.1 that is supported by the hindrance model. This nuclear physics assumption, combined with C flame quenching, leads to the formation of a hybrid C-O-Ne WD with a mass of 1.3 Msun. Such WDs do not need to accrete much mass to reach the Chandrasekhar limit.

preprint2013arXiv

Advanced burning stages and fate of 8-10 Mo stars

The stellar mass range 8<M/Mo<12 corresponds to the most massive AGB stars and the most numerous massive stars. It is host to a variety of supernova progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB star to massive star and find that a propagating neon-oxygen burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core collapse supernova (FeCCSN) progenitors. Of the models that ignite neon burning off-center, the 9.5Mo model would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8Mo model, however, fails to reach the center as the URCA process and an extended (0.6 Mo) region of low Ye (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses. We also present models of an 8.75Mo super-AGB star through its entire thermal pulse phase until electron captures on 20Ne begin at its center and of a 12Mo star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-supernova structure. Finally, we compare our results to the observed neutron star mass distribution.

preprint2013arXiv

Evidence for Radiogenic Sulfur-32 in Type AB Presolar Silicon Carbide Grains?

We report C, Si, and S isotope measurements on 34 presolar silicon carbide grains of Type AB, characterized by 12C/13C < 10. Nitrogen, Mg-Al-, and Ca-Ti-isotopic compositions were measured on a subset of these grains. Three grains show large 32S excesses, a signature that has been previously observed for grains from supernovae (SNe). Enrichments in 32S may be due to contributions from the Si/S zone and the result of S molecule chemistry in still unmixed SN ejecta or due to incorporation of radioactive 32Si from C-rich explosive He shell ejecta. However, a SN origin remains unlikely for the three AB grains considered here, because of missing evidence for 44Ti, relatively low 26Al/27Al ratios (a few times 10-3), and radiogenic 32S along with low 12C/13C ratios. Instead, we show that born-again asymptotic giant branch (AGB) stars that have undergone a very-late thermal pulse (VLTP), known to have low 12C/13C ratios and enhanced abundances of the light s-process elements, can produce 32Si, which makes such stars attractive sources for AB grains with 32S excesses. This lends support to the proposal that at least some AB grains originate from born-again AGB stars, although uncertainties in the born-again AGB star models and possible variations of initial S-isotopic compositions in the parent stars of AB grains make it difficult to draw a definitive conclusion.

preprint2013arXiv

Relics of ancient post-AGB stars in a primitive meteorite

Graphite is one of the many presolar circumstellar condensate species found in primitive meteorites. While the isotopic compositions of low-density graphite grains indicate an origin in core-collapse supernovae, some high-density grains have extreme isotopic anomalies in C, Ca and Ti, which cannot be explained by envelope predictions of asymptotic giant branch (AGB) stars or theoretical supernova models. The Ca and Ti isotopic anomalies, however, match the predictions of He-shell abundances in AGB stars. In this study, we show that the C, Ca, and Ti isotopic anomalies are consistent with nucleosynthesis predictions of the H-ingestion phase during a very late thermal pulse (VLTP) event in post-AGB stars. The low $^{12}$C/$^{13}$C isotopic ratios in these grains are a result of abundant $^{12}$C efficiently capturing the protons that are being ingested during the VLTP. Very high neutron densities of $\sim 10^{15}$ cm$^{-3}$, typical of the $i$-process, are achieved during this phase in post-AGB stars. The large $^{42,43,44}$Ca excesses in some graphite grains are indicative of neutron capture nucleosynthesis during VLTP. The comparison of VLTP nucleosynthesis calculations to the graphite data also indicate that apparent anomalies in the Ti isotopic ratios are due to large contributions from $^{46,48}$Ca, which cannot be resolved from the isobars $^{46,48}$Ti during the measurements. We conclude that presolar graphite grains with moderate to extreme Ca and Ti isotopic anomalies originate in post-AGB stars that suffer a very late thermal pulse.

preprint2013arXiv

Reproducing the observed abundances in RCB and HdC stars with post-double degenerate merger models - constraints on merger and post-merger simulations and physics processes

The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars that are most likely the result of He-CO WD mergers. They display extremely low oxygen isotopic ratios, 16O/18O ~ 1 - 10, 12C/13C>=100, and enhancements up to 2.6dex in F and in s-process elements from Zn to La, compared to solar. These abundances provide stringent constraints on the physical processes during and after the double-degenerate merger. As shown before O-isotopic ratios observed in RCB stars cannot result from the dynamic double-degenerate merger phase, and we investigate now the role of the long-term 1D spherical post-merger evolution and nucleosynthesis based on realistic hydrodynamic merger progenitor models. We adopt a model for extra envelope mixing to represent processes driven by rotation originating in the dynamical merger. Comprehensive nucleosynthesis post-processing simulations for these stellar evolution models reproduce, for the first time, the full range of the observed abundances for almost all the elements measured in RCB stars: 16O/18O ratios between 9 and 15, C-isotopic ratios above 100, and ~1.4 - 2.35dex F enhancements, along with enrichments in s-process elements. The nucleosynthesis processes in our models constrain the length and temperature in the dynamic merger shell-of-fire feature as well as the envelope mixing in the post-merger phase. s-process elements originate either in the shell-of-fire merger feature or during the post-merger evolution, but the contribution from the AGB progenitors is negligible. The post-merger envelope mixing must eventually cease ~ 10^6yr after the dynamic merger phase, before the star enters the RCB phase.

preprint2013arXiv

The C-flame Quenching by Convective Boundary Mixing in Super-AGB Stars and the Formation of Hybrid C/O/Ne White Dwarfs and SN Progenitors

After off-center C ignition in the cores of super-AGB stars the C flame propagates all the way down to the center, trailing behind it the C-shell convective zone, and thus building a degenerate ONe core. This standard picture is obtained in stellar evolution simulations if the bottom C-shell convection boundary is assumed to be a discontinuity associated with a strict interpretation of the Schwarzschild condition for convective instability. However, this boundary is prone to additional mixing processes, such as thermohaline convection and convective boundary mixing. Using hydrodynamic simulations we show that, contrary to previous results, thermohaline mixing is too inefficient to interfere with the C-flame propagation. However, even a small amount of convective boundary mixing removes the physical conditions required for the C-flame propagation all the way to the center. This result holds even if we allow for some turbulent heat transport in the CBM region. As a result, super AGB stars build in their interiors hybrid C-O-Ne degenerate cores composed of a relatively large CO core (M_CO ~ 0.2 M_sun) surrounded by a thick ONe zone (M_ONe ~ 0.85 M_sun) with another thin CO layer above. If exposed by mass loss, these cores will become hybrid C-O-Ne white dwarfs. Otherwise, the ignition of C-rich material in the central core, surrounded by the thick ONe zone, may trigger a thermonuclear supernova explosion. The quenching of the C-flame may have implications for the ignition mechanism of SN Ia in the double-degenerate merger scenario.

preprint2012arXiv

Do R Coronae Borealis Stars Form from Double White Dwarf Mergers?

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

preprint2012arXiv

From the CMD of Omega Centauri and (super-)AGB stellar models to a Galactic plane passage gas purging chemical evolution scenario

[Abbreviated] We have investigated the color-magnitude diagram of Omega Centauri and find that the blue main sequence (bMS) can be reproduced only by models that have a of helium abundance in the range Y=0.35-$0.40. To explain the faint subgiant branch of the reddest stars ("MS-a/RG-a" sequence), isochrones for the observed metallicity ([Fe/H]\approx0.7) appear to require both a high age (~13Gyr) and enhanced CNO abundances ([CNO/Fe]\approx0.9$). Y~0.35 must also be assumed in order to counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively blue turnoff of this stellar population. This suggest a short chemical evolution period of time (<1Gyr) for Omega Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium abundances, along with [N/Fe]~2 and substantial O depletions if uncertainties in the treatment of convection are fully taken into account. These abundance features distinguish the bMS stars from the dominant [Fe/H] $\approx1.7$ population. The most massive super-AGB stellar models (M_zams>=6.8M_sun, M_He,core>=1.245M_sun) predict too large N-enhancements, which limits their role in contributing to the extreme populations. We show quantitatively that highly He- and N-enriched AGB ejecta have particularly efficient cooling properties. Based on these results and on the reconstruction of the orbit of Omega Cen with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the chemical evolution of this cluster. Our model addresses the formation and properties of the bMS population (including their central location in the cluster). We follow our model descriptively through four passage events, which could explain not only some key properties of the bMS, but also of the MS-a/RGB-a and the s-enriched stars.

preprint2012arXiv

MESA Models of Classical Nova Outbursts: The Multicycle Evolution and Effects of Convective Boundary Mixing

Novae are cataclysmic variables driven by accretion of H-rich material onto a white-dwarf (WD) star from its low-mass main-sequence binary companion. New time-domain observational capabilities, such as the Palomar Transient Factory and Pan-STARRS, have revealed a diversity of their behaviour that should be theoretically addressed. Nova outbursts depend sensitively on nuclear physics data, and more readily available nova simulations are needed in order to effectively prioritize experimental effort in nuclear astrophysics. In this paper we use the MESA stellar evolution code to construct multicycle nova evolution sequences with CO WD cores. We explore a range of WD masses and accretion rates as well as the effect of different cooling times before the onset of accretion. In addition, we study the dependence on the elemental abundance distribution of accreted material and convective boundary mixing at the core-envelope interface. Models with such convective boundary mixing display an enrichment of the accreted envelope with C and O from the underlying white dwarf that is commensurate with observations. We compare our results with the previous work and investigate a new scenario for novae with the 3He-triggered convection.

preprint2012arXiv

Nova Framework: A New Tool For Modeling of Nova Outbursts and Nucleosynthesis

Classical novae are the results of surface thermonuclear explosions of H-rich material accreted by white dwarfs (WDs) from their low-mass main-sequence or red-giant binary companions. Chemical composition analysis of their ejecta shows that nova outbursts occur on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs, and that there is cross-boundary mixing between the accreted envelope and underlying WD. We have combined the stellar evolution code MESA and post-processing nucleosynthesis tools of NuGrid into a framework that allows to produce up-to-date models of nova outbursts and compute detailed nucleosynthesis in novae occurring on CO and ONe WDs. The convective boundary mixing (CBM) in our 1D numerical simulations is implemented using a diffusion coefficient that is exponentially decreasing with a distance below the bottom of the convective envelope. This MESA CBM prescription is based on the findings in 3D hydrodynamic simulations that the velocity field, and along with it the mixing expressed in terms of a diffusion coefficient, decays exponentially in the stable layer adjacent to a convective boundary. The framework can also use the commonly adopted 1D nova model in which the CBM is mimicked by assuming that the accreted envelope has been pre-enriched with WD's material. In this preliminary report, we present the most interesting new results related to CO and ONe nova outbursts that have been obtained with the Nova Framework.

preprint2012arXiv

The Response of Giant Stars To Dynamical-Timescale Mass Loss

We study the response of giant stars to mass loss. One-dimensional simulations of red and asymptotic giant branch stars with mass loss rates from $10^{-3}$ up to a few \msun/yr show in no case any significant radius increase. The largest radius increase of 0.2% was found in the case with the lowest mass loss rate. For dynamical-timescale mass loss rates, that may be encountered during a common envelope phase, the evolution is not adiabatic. The superadiabatic outer layer of the giant's envelope has a local thermal timescale comparable to the dynamical timescale. Therefore, this layer has enough time to readjust thermally. Moreover, the giant star is driven out of hydrostatic equilibrium and evolves dynamically. In these cases no increase of the stellar radius with respect to its initial value is found. If the mass loss rate is high enough, the superadiabaticity of the outer layer is lost progressively and a radiative zone forms due to a combination of thermal and dynamical readjustment. Conditions for unstable mass transfer based on adiabatic mass loss models that predict a significant radius increase, may need to be re-evaluated.

preprint2011arXiv

[O/Fe] Estimates for Carbon-Enhanced Metal-Poor Stars from Near-IR Spectroscopy

We report on oxygen abundances determined from medium-resolution near-IR spectroscopy for a sample of 57 carbon-enhanced metal-poor (CEMP) stars selected from the Hamburg/ESO survey. The majority of our program stars exhibit oxygen-to-iron ratios in the range +0.5 < [O/Fe]< +2.0. The [O/Fe] values for this sample are statistically compared to available high-resolution estimates for known CEMP stars, as well as to high-resolution estimates for a set of carbon-normal metal-poor stars. Carbon, nitrogen, and oxygen abundance patterns for a sub-sample of these stars are compared to yield predictions for very metal-poor asymptotic giant-branch abundances in the recent literature. We find that the majority of our sample exhibit patterns that are consistent with previously studied CEMP stars having s-process-element enhancements, and thus have very likely been polluted by carbon- and oxygen-enhanced material transferred from a metal-poor asymptotic giant-branch companion.

preprint2011arXiv

Convective-reactive proton-C12 combustion in Sakurai's object (V4334 Sagittarii) and implications for the evolution and yields from the first generations of stars

Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(α, n)O16 reaction.

preprint2011arXiv

Simulating the Common Envelope Phase of a Red Giant Using SPH and Uniform Grid Codes

We use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 \rsun, are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.

preprint2010arXiv

Fluorine and Sodium in C-rich low-metallicity stars

We present the N, O, F and Na abundance and 12C/13C isotopic ratio measurements or upper limits for a sample of 10 C-rich, metal-poor giant stars, eight enhanced in s-process (CEMP-s) elements and two poor in n-capture elements (CEMP-no). The abundances are derived from IR, K-band, high-resolution CRIRES@VLT spectra obtained. The metallicity of our sample ranges from [Fe/H]=-3.4 to -1.3. F abundance could be measured only in two CEMP-s stars. With [F/Fe]=0.64, one is mildly F-overabundant, while the other is F-rich, at [F/Fe]=1.44. For the remaining eight objects, including both CEMP-no in our sample, only upper limits on F abundance could be placed. Our measurements and upper limits show that there is a spread in [F/C+N] ratio in CEMP-s stars as predicted by theory. Predictions from nucleosynthetic models for low-mass, low-metallicity Asymptotic Giant Branch stars, account for the derived F abundances, while the upper limits on F content derived for most of the stars are lower than the predicted values. The measured Na content is accounted for by AGB models in the 1.25 to 1.7 Msun range, confirming that the stars responsible for the peculiar abundance pattern observed in CEMP-s stars are low-mass, low-metallicity AGB stars, in agreement with the most accepted astrophysical scenario. We conclude that the mechanism of F production in current state-of-the-art low-metallicity low-mass AGB models needs further scrutiny and that F measurements in a larger number of metal-poor stars are needed to better constraint the models.

preprint2010arXiv

Modules for Experiments in Stellar Astrophysics (MESA)

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very-low mass to massive stars, including advanced evolutionary phases. MESA star solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (http://mesa.sourceforge.net/).

preprint2010arXiv

On the alpha formalism for the common envelope interaction

The α-formalism is a common way to parametrize the common envelope interaction between a giant star and a more compact companion. The α parameter describes the fraction of orbital energy released by the companion that is available to eject the giant star's envelope. By using new, detailed stellar evolutionary calculations we derive a user-friendly prescription for the λ parameter and an improved approximation for the envelope binding energy, thus revising the α equation. We then determine α both from simulations and observations in a self consistent manner. By using our own stellar structure models as well as population considerations to reconstruct the primary's parameters at the time of the common envelope interaction, we gain a deeper understanding of the uncertainties. We find that systems with very low values of q (the ratio of the companion's mass to the mass of the primary at the time of the common envelope interaction) have higher values of α. A fit to the data suggests that lower mass companions are left at comparable or larger orbital separations to more massive companions. We conjecture that lower mass companions take longer than a stellar dynamical time to spiral in to the giant's core, and that this is key to allowing the giant to use its own thermal energy to help unbind its envelope. As a result, although systems with light companions might not have enough orbital energy to unbind the common envelope, they might stimulate a stellar reaction that results in the common envelope ejection.

preprint2010arXiv

The effect of 12C + 12C rate uncertainties on the weak s-process component

The contribution by massive stars (M > 9 solar masses) to the weak s-process component of the solar system abundances is primarily due to the 22Ne neutron source, which is activated near the end of helium-core burning. The residual 22Ne left over from helium-core burning is then reignited during carbon burning, initiating further s-processing that modifies the isotopic distribution. This modification is sensitive to the stellar structure and the carbon burning reaction rate. Recent work on the 12C + 12C reaction suggests that resonances located within the Gamow peak may exist, causing a strong increase in the astrophysical S-factor and consequently the reaction rate. To investigate the effect of an increased rate, 25 solar mass stellar models with three different carbon burning rates, at solar metallicity, were generated using the Geneva Stellar Evolution Code (GENEC) with nucleosynthesis post-processing calculated using the NuGrid Multi-zone Post-Processing Network code (MPPNP). The strongest rate caused carbon burning to occur in a large convective core rather than a radiative one. The presence of this large convective core leads to an overlap with the subsequent convective carbon-shell, significantly altering the initial composition of the carbon-shell. In addition, an enhanced rate causes carbon-shell burning episodes to ignite earlier in the evolution of the star, igniting the 22Ne source at lower temperatures and reducing the neutron density.

preprint2007arXiv

Very Large Excesses of 18O in Hydrogen-Deficient Carbon and R Coronae Borealis Stars: Evidence for White Dwarf Mergers

We have found that at least seven hydrogen-deficient carbon (HdC) and R Coronae Borealis (RCB) stars, have 16O/18O ratios close to and in some cases less than unity, values that are orders of magnitude lower than measured in other stars (the Solar value is 500). Greatly enhanced 18O is evident in every HdC and RCB we have measured that is cool enough to have detectable CO bands. The three HdC stars measured have 16O/18O < 1, lower values than any of the RCB stars. These discoveries are important clues in determining the evolutionary pathways of HdC and RCB stars, for which two models have been proposed: the double degenerate (white dwarf (WD) merger), and the final helium-shell flash (FF). No overproduction of 18O is expected in the FF scenario. We have quantitatively explored the idea that HdC and RCB stars originate in the mergers of CO- and He-WDs. The merger process is estimated to take only a few days, with accretion rates of 150 Msun/ yr producing temperatures at the base of the accreted envelope of 1.2 - 1.9 x 10^8 K. Analysis of a simplified one-zone calculation shows that nucleosynthesis in the dynamically accreting material may provide a suitable environment for a significant production of 18O, leading to very low values of 16O/18O, similar to those observed. We also find qualitative agreement with observed values of 12C/13C and with the CNO elemental ratios. H-admixture during the accretion process from the small H-rich C/O WD envelope may play an important role in producing the observed abundances. Overall our analysis shows that WD mergers may very well be the progenitors of O18-rich RCB and HdC stars, and that more detailed simulations and modeling are justified.

preprint2006arXiv

A Search for Nitrogen-Enhanced Metal-Poor Stars

Theoretical models of very metal-poor intermediate-mass Asymptotic Giant Branch (AGB) stars predict a large overabundance of primary nitrogen. The very metal-poor, carbon-enhanced, s-process-rich stars, which are thought to be the polluted companions of now-extinct AGB stars, provide direct tests of the predictions of these models. Recent studies of the carbon and nitrogen abundances in metal-poor stars have focused on the most carbon-rich stars, leading to a potential selection bias against stars that have been polluted by AGB stars that produced large amounts of nitrogen, and hence have small [C/N] ratios. We call these stars Nitrogen-Enhanced Metal-Poor (NEMP) stars, and define them as having [N/Fe] > +0.5 and [C/N] < -0.5. In this paper, we report on the [C/N] abundances of a sample of 21 carbon-enhanced stars, all but three of which have [C/Fe] < +2.0. If NEMP stars were made as easily as Carbon-Enhanced Metal-Poor (CEMP) stars, then we expected to find between two and seven NEMP stars. Instead, we found no NEMP stars in our sample. Therefore, this observational bias is not an important contributor to the apparent dearth of N-rich stars. Our [C/N] values are in the same range as values reported previously in the literature (-0.5 to +2.0), and all stars are in disagreement with the predicted [C/N] ratios for both low-mass and high-mass AGB stars. We suggest that the decrease in [C/N] from the low-mass AGB models is due to enhanced extra-mixing, while the lack of NEMP stars may be caused by unfavorable mass ratios in binaries or the difficulty of mass transfer in binary systems with large mass ratios.

preprint2005arXiv

An Extremely Large Excess of 18O in the Hydrogen-Deficient Carbon Star, HD 137613

We report the discovery of a uniquely large excess of 18O in the hydrogen-deficient carbon (HdC) star, HD 137613, based on a spectrum of the first overtone bands of CO at 2.3-2.4 micron in which three strong absorption bands of 12C18O are clearly present. Bands of 12C16O also are present but no bands of 13C16O or 12C17O are seen. We estimate an isotopic ratio 16O/18O ~1. The Solar value of this ratio is ~500. Neither He-core burning nor He-shell flash burning can produce the isotopic ratios of oxygen and carbon observed in HD 137613. However, a remarkable similarity exists between the observed abundances and those found in the outer layers of the broad He-shell of early-AGB stars, soon after the end of He-core burning. It is not known how the outer envelope down to the He-shell could be lost but some mechanism of enhanced mass loss must be involved. HD 137613 may be a post-early-AGB star with the outer layers of the former He-burning shell as its photosphere. The unusual elemental abundances of the HdC stars resemble those of the R Coronae Borealis (RCB) stars, but HdC stars do not produce clouds of dust that produce declines in brightness. None of the other RCB or HdC stars observed shows significant 18O.

preprint2004arXiv

Evolution and Yields of Extremely Metal Poor Intermediate Mass Stars

Intermediate mass stellar evolution tracks from the main sequence to the tip of the AGB for five initial masses (2 to 6Msun) and metallicity Z=0.0001 have been computed. The detailed 1D structure and evolution models include exponential overshooting, mass loss and a detailed nucleosynthesis network with updated nuclear reaction rates. The network includes a two-particle heavy neutron sink for approximating neutron density in the He-shell flash. It is shown how the neutron-capture nucleosynthesis is important in models of very low metallicity for the formation of light neutron-heavy species, like sodium or the heavy neon and magnesium isotopes. The models have high resolution, as required for modeling the third dredge-up. All sequences have been followed from the pre-main sequence to the end of the AGB when all envelope mass is lost. Detailed structural and chemical model properties as well as yields are presented. This set of stellar models is based on standard assumptions and updated input physics. It can be confronted with observations of extremely-metal poor stars and may be used to assess the role of AGB stars in the origin of abundance anomalies of some Globular Cluster members of correspondingly low metallicity.

preprint2002arXiv

Past and future evolution of Sakurai's object

We present a stellar evolution model sequence of the past and future evolution of the post-AGB born again star Sakurai's object (V4334 Sgr). In order to match the short evolutionary time scale we have to assume that the convective ingestion of hydrogen-rich envelope material into the He-flash convection zone proceeds slower than predicted by the mixing length theory. For the future we predict a swift second evolution through the central star region before a second born-again evolution occurs.

preprint2001arXiv

Stellar evolution and nucleosynthesis of Post-AGB Stars

I discuss recent new models of post-Asymptotic Giant Branch stellar evolution. These models aim to clarify the evolutionary origin and status of a variety of hydrogen-deficient post-AGB stars such as central stars of planetary nebulae of Wolf-Rayet spectral type, PG1159 stars or Sakurai's object. Starting with AGB models with overshoot such stars can evolve through one of four distinct channels. Each of these channels has typical abundance patterns depending on the relative timing of the departure from the AGB and the occurrence of the last thermal pulse. I discuss the responsible mechanisms and observational counterparts.