BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

F. X. Timmes

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-ph.IMastro-phastro-ph.COnucl-exastro-ph.EPastro-ph.GAgr-qcphysics.ed-phphysics.flu-dyn

Catalog footprint

What is connected

39works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

39 published item(s)

preprint2022arXiv

Localised thermonuclear bursts from accreting magnetic white dwarfs

Nova explosions are caused by global thermonuclear runaways triggered in the surface layers of accreting white dwarfs. It has been predicted that localised thermonuclear bursts on white dwarfs can also take place, similar to Type I X-ray bursts observed in accreting neutron stars. Unexplained rapid bursts from the binary system TV Columbae, in which mass is accreted onto a moderately-strong magnetised white dwarf from a low-mass companion, have been observed on several occasions in the past $\approx40$ years. During these bursts the optical/UV luminosity increases by a factor of $>3$ in less than an hour and fades over $\approx10$ hours. Fast outflows have been observed in UV spectral lines, with velocities $>3500$ km s$^{-1}$, comparable to the escape velocity from the white dwarf surface. Here we report on optical bursts observed in TV Columbae as well as in two additional accreting systems, EI Ursae Majoris and ASASSN-19bh. The bursts have a total energy $\approx~10^{-6}$ those of classical nova explosions ("micronovae"), and bear a strong resemblance to Type I X-ray bursts. We exclude accretion or stellar magnetic reconnection events as their origin and suggest thermonuclear runaway events in magnetically-confined accretion columns as a viable explanation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On Trapped Modes In Variable White Dwarfs As Probes Of The $^{12}$C($α, γ$)$^{16}$O Reaction Rate

We seek signatures of the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function in the pulsation periods of carbon-oxygen white dwarf models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function. We find an average spread in relative period shifts of $ΔP/P \simeq \pm$ 2\% for the identified trapped g-modes over the $\pm$ 3$σ$ uncertainty in the $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function -- across the effective temperature range of observed DAV and DBV white dwarfs and for different white dwarf masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed white dwarfs are typically given to 6-7 significant figures of precision. This suggests that an astrophysical constraint on the $^{12}$C$(α,γ)^{16}$O reaction rate could, in principle, be extractable from the period spectrum of observed variable white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Observing intermediate-mass black holes and the upper--stellar-mass gap with LIGO and Virgo

Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap $\sim 60-130~M_\odot$. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing (O4) run, we perform Bayesian analysis on quasi-circular non-precessing, spinning IMBH binaries (IMBHBs) with total masses $50\mbox{--} 500\, M_\odot$, mass ratios 1.25, 4, and 10, and dimensionless spins up to 0.95, and estimate the precision with which the source-frame parameters can be measured. We find that, at $2σ$, the mass of the heavier component of IMBHBs can be constrained with an uncertainty of $\sim 10-40\%$ at a signal-to-noise ratio of $20$. Focusing on the stellar-mass gap with new tabulations of the $^{12}\text{C}(α, γ)^{16} \text{O}$ reaction rate and its uncertanties, we evolve massive helium core stars using \MESA\, to establish the lower and upper edge of the mass gap as $\simeq$\,59$^{+34}_{-13}$\,$M_{\odot}$ and $\simeq$\,139$^{+30}_{-14}$\,$M_{\odot}$ respectively, where the error bars give the mass range that follows from the $\pm 3σ$ uncertainty in the $^{12}\text{C}(α, γ) ^{16} \text{O}$ nuclear reaction rate. We find that high resolution of the tabulated reaction rate and fine temporal resolution are necessary to resolve the peak of the BH mass spectrum. We then study IMBHBs with components lying in the mass gap and show that the O4 run will be able to robustly identify most such systems. Finally, we re-analyse GW190521 with a state-of-the-art aligned-spin waveform model, finding that the primary mass lies in the mass gap with 90\% credibility.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Skye: A Differentiable Equation of State

Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully-ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, non-linear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novel thermodynamic extrapolation scheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in the MESA stellar evolution software instrument by computing white dwarf cooling curves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Laminar Flame Speeds in Degenerate Oxygen-Neon Mixtures

The collapse of degenerate oxygen-neon cores (i.e., electron-capture supernovae or accretion-induced collapse) proceeds through a phase in which a deflagration wave ("flame") forms at or near the center and propagates through the star. In models, the assumed speed of this flame influences whether this process leads to an explosion or to the formation of a neutron star. We calculate the laminar flame speeds in degenerate oxygen-neon mixtures with compositions motivated by detailed stellar evolution models. These mixtures include trace amounts of carbon and have a lower electron fraction than those considered in previous work. We find that trace carbon has little effect on the flame speeds, but that material with electron fraction $Y_e \approx 0.48-0.49$ has laminar flame speeds that are $\approx 2$ times faster than those at $Y_e = 0.5$. We provide tabulated flame speeds and a corresponding fitting function so that the impact of this difference can be assessed via full star hydrodynamical simulations of the collapse process.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Tracer Method for Computing Type Ia Supernova Yields: Burning Model Calibration, Reconstruction of Thickened Flames, and Verification for Planar Detonations

We refine our previously introduced parameterized model for explosive carbon-oxygen fusion during thermonuclear supernovae (SN Ia) by adding corrections to post-processing of recorded Lagrangian fluid element histories to obtain more accurate isotopic yields. Deflagration and detonation products are verified for propagation in a uniform density medium. A new method is introduced for reconstructing the temperature-density history within the artificially thick model deflagration front. We obtain better than 5\% consistency between the electron capture computed by the burning model and yields from post-processing. For detonations, we compare to a benchmark calculation of the structure of driven steady-state planar detonations performed with a large nuclear reaction network and error-controlled integration. We verify that, for steady-state planar detonations down to a density of 5x10^6 g/cc, our post processing matches the major abundances in the benchmark solution typically to better than 10% for times greater than 0.01 s after the shock front passage. As a test case to demonstrate the method, presented here with post-processing for the first time, we perform a two dimensional simulation of a SN Ia in the Chandrasekhar-mass deflagration-detonation transition (DDT) scenario. We find that reconstruction of deflagration tracks leads to slightly more complete silicon burning than without reconstruction. The resulting abundance structure of the ejecta is consistent with inferences from spectroscopic studies of observed SNe Ia. We confirm the absence of a central region of stable Fe-group material for the multi-dimensional DDT scenario. Detailed isotopic yields are tabulated and only change modestly when using deflagration reconstruction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convective Properties of Rotating Two-Dimensional Core-Collapse Supernova Progenitors

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On Measuring the Metallicity of a Type Ia Supernova's Progenitor

In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,Å and a Fe feature around 5200\,Å\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On Variations Of Pre-Supernova Model Properties

We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Properties of Carbon-Oxygen White Dwarfs From Monte Carlo Stellar Models

We investigate properties of carbon-oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M$_{\odot}$ models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95\% confidence interval to be $ΔM_{\rm 1TP}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $Δ$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $Δ\log(T_{\rm c}/{\rm K}) \approx$ 0.013 for the central temperature, $Δ\log(ρ_{\rm c}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $ΔY_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $ΔX_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $ΔX_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $ΔX_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($α,γ)^{16}\rm{O}$, triple-$α$, and $^{14}$N($p,γ)^{15}\rm{O}$ reaction rates dominate these variations. We also consider a grid of 1 to 6 M$_{\odot}$ models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial-final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Importance of Computation in Astronomy Education

Computational skills are required across all astronomy disciplines. Many students enter degree programs without sufficient skills to solve computational problems in their core classes or contribute immediately to research. We recommend advocacy for computational literacy, familiarity with fundamental software carpentry skills, and mastery of basic numerical methods by the completion of an undergraduate degree in Astronomy. We recommend the AAS Education Task Force advocate for a significant increase in computational literacy. We encourage the AAS to modestly fund efforts aimed at providing Open Education Resources (OER) that will significantly impact computational literacy in astronomy education.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraints on explosive silicon burning in core-collapse supernovae from measured Ni/Fe ratios

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $η\approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling Astrophysical Explosions with Sustained Exascale Computing

Our understanding of stars and their fates is based on coupling observations to theoretical models. Unlike laboratory physicists, we cannot perform experiments on stars, but rather must patiently take what nature allows us to observe. Simulation offers a means of virtual experimentation, enabling a detailed understanding of the most violent ongoing explosions in the Universe---the deaths of stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On Carbon Burning in Super Asymptotic Giant Branch Stars

We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In general terms, we find these properties of carbon burning in SAGB models are not a strong function of the initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically derive an approximate ignition density, $ρ_{ign} \approx 2.1 \times 10^6$ g cm$^{-3}$, to predict the location of first carbon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of $ΔM_{\rm ZAMS}$/$Δf_{\rm{ov}}\approx$ 1.6 $M_{\odot}$. For zero overshoot, $f_{\rm{ov}}$=0.0, our models in the ZAMS mass range $\approx$ 8.9 to 11 $M_{\odot}$ show off-center carbon ignition. For canonical amounts of overshooting, $f_{\rm{ov}}$=0.016, the off-center carbon ignition range shifts to $\approx$ 7.2 to 8.8 $M_{\odot}$. Only systems with $f_{\rm{ov}}$ $\geq 0.01$ and ZAMS mass $\approx$ 7.2-8.0 $M_{\odot}$ show carbon burning is quenched a significant distance from the center. These results suggest a careful assessment of overshoot modeling approximations on claims that carbon burning quenches an appreciable distance from the center of the carbon core.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Three Dimensional Evolution to Core Collapse of a Massive Star

We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasekhar) mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion (turbulent fluctuations) generated by 3D convection is substantial at the point of collapse. We examine the impact of such physically-realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage. We conclude that non-spherical progenitor structure should not be ignored, and has a significant and favorable impact on the likelihood for neutrino-driven explosions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Silicon Group Elements Ejected by Supernovae Type Ia

There is compelling evidence that the peak brightness of a Type Ia supernova is affected by the electron fraction Ye at the time of the explosion. The electron fraction is set by the aboriginal composition of the white dwarf and the reactions that occur during the pre explosive convective burning. To date, determining the makeup of the white dwarf progenitor has relied on indirect proxies, such as the average metallicity of the host stellar population. In this paper, we present analytical calculations supporting the idea that the electron fraction of the progenitor systematically influences the nucleosynthesis of silicon group ejecta in Type Ia supernovae. In particular, we suggest the abundances generated in quasi nuclear statistical equilibrium are preserved during the subsequent freezeout. This allows one to potential recovery of Ye at explosion from the abundances recovered from an observed spectra. We show that measurement of 28Si, 32S, 40Ca, and 54Fe abundances can be used to construct Ye in the silicon rich regions of the supernovae. If these four abundances are determined exactly, they are sufficient to recover Ye to 6 percent. This is because these isotopes dominate the composition of silicon-rich material and iron rich material in quasi nuclear statistical equilibrium. Analytical analysis shows that the 28Si abundance is insensitive to Ye, the 32S abundance has a nearly linear trend with Ye, and the 40Ca abundance has a nearly quadratic trend with Ye. We verify these trends with post-processing of 1D models and show that these trends are reflected in model synthetic spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical Methods for Thermonuclear Reaction Rates and Nucleosynthesis Simulations

Rigorous statistical methods for estimating thermonuclear reaction rates and nucleosynthesis are becoming increasingly established in nuclear astrophysics. The main challenge being faced is that experimental reaction rates are highly complex quantities derived from a multitude of different measured nuclear parameters (e.g., astrophysical S-factors, resonance energies and strengths, particle and gamma-ray partial widths). We discuss the application of the Monte Carlo method to two distinct, but related, questions. First, given a set of measured nuclear parameters, how can one best estimate the resulting thermonuclear reaction rates and associated uncertainties? Second, given a set of appropriate reaction rates, how can one best estimate the abundances from nucleosynthesis (i.e., reaction network) calculations? The techniques described here provide probability density functions that can be used to derive statistically meaningful reaction rates and final abundances for any desired coverage probability. Examples are given for applications to s-process neutron sources, core-collapse supernovae, classical novae, and big bang nucleosynthesis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stellar Abundances in the Solar Neighborhood: The Hypatia Catalog

We compile spectroscopic abundance data from 84 literature sources for 50 elements across 3058 stars in the solar neighborhood, within 150 pc of the Sun, to produce the Hypatia Catalog. We evaluate the variability of the spread in abundance measurements reported for the same star by different surveys. We also explore the likely association of the star within the Galactic disk, the corresponding observation and abundance determination methods for all catalogs in Hypatia, the influence of specific catalogs on the overall abundance trends, and the effect of normalizing all abundances to the same solar scale. The resulting large number of stellar abundance determinations in the Hypatia Catalog are analyzed only for thin-disk stars with observations that are consistent between literature sources. As a result of our large dataset, we find that the stars in the solar neighborhood may be reveal an asymmetric abundance distribution, such that a [Fe/H]-rich group near to the mid-plane is deficient in Mg, Si, S, Ca, Sc II, Cr II, and Ni as compared to stars further from the plane. The Hypatia Catalog has a wide number of applications, including exoplanet hosts, thick and thin disk stars, or stars with different kinematic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

STARLIB: A Next-Generation Reaction-Rate Library for Nuclear Astrophysics

STARLIB is a next-generation, all-purpose nuclear reaction-rate library. For the first time, this library provides the rate probability density at all temperature grid points for convenient implementation in models of stellar phenomena. The recommended rate and its associated uncertainties are also included. Currently, uncertainties are absent from all other rate libraries, and, although estimates have been attempted in previous evaluations and compilations, these are generally not based on rigorous statistical definitions. A common standard for deriving uncertainties is clearly warranted. STARLIB represents a first step in addressing this deficiency by providing a tabular, up-to-date database that supplies not only the rate and its uncertainty but also its distribution. Because a majority of rates are lognormally distributed, this allows the construction of rate probability densities from the columns of STARLIB. This structure is based on a recently suggested Monte Carlo method to calculate reaction rates, where uncertainties are rigorously defined. In STARLIB, experimental rates are supplemented with: (i) theoretical TALYS rates for reactions for which no experimental input is available, and (ii) laboratory and theoretical weak rates. STARLIB includes all types of reactions of astrophysical interest to Z = 83, such as (p,g), (p,a), (a,n), and corresponding reverse rates. Strong rates account for thermal target excitations. Here, we summarize our Monte Carlo formalism, introduce the library, compare methods of correcting rates for stellar environments, and discuss how to implement our library in Monte Carlo nucleosynthesis studies. We also present a method for accessing STARLIB on the Internet and outline updated Monte Carlo-based rates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mixing of Clumpy Supernova Ejecta into Molecular Clouds

Several lines of evidence, from isotopic analyses of meteorites to studies of the Sun's elemental and isotopic composition, indicate that the solar system was contaminated early in its evolution by ejecta from a nearby supernova (SN). Previous models have invoked SN material being injected into an extant protoplanetary disk, or isotropically expanding ejecta sweeping over a distant (>10 pc) cloud core, simultaneously enriching it and triggering its collapse. Here we consider a new astrophysical setting: the injection of clumpy SN ejecta, as observed in the Cas A SN remnant, into the molecular gas at the periphery of an HII region created by the SN's progenitor star. To track these interactions we have conducted a suite of high-resolution (1500^3 effective) 3D simulations that follow the evolution of individual clumps as they move into molecular gas. Even at these high resolutions, our simulations do not quite achieve numerical convergence, due to the challenge of properly resolving the small-scale mixing of ejecta and molecular gas, although they do allow some robust conclusions to be drawn. Isotropically exploding ejecta do not penetrate into the molecular cloud, but, if cooling is properly accounted for, clumpy ejecta penetrate to distances ~10^18 cm and mix effectively with star-forming molecular gas. The ~2 M_\odot high-metallicity ejecta from a core-collapse SN is likely to mix with ~2 \times 10^4 M_\odot of molecular gas. Thus all stars forming late (~5 Myr) in the evolution of an HII region may be contaminated by SN ejecta at a level ~10^-4. This level of contamination is consistent with the abundances of short-lived radionuclides and possibly some stable isotopic shifts in the early solar system, and is potentially consistent with the observed variability in stellar elemental abundances. SN contamination of forming planetary systems may be a common, universal process.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Theoretical Studies of Accretion of Matter onto White Dwarfs and the Single Degenerate Scenario for Supernovae of Type Ia

We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Freeze-out yields of radioactivities in core-collapse supernovae

We explore the nucleosynthesis trends from two mechanisms during freeze-out expansions in core-collapse supernovae. The first mechanism is related to the convection and instabilities within homogeneous stellar progenitor matter that is accreted through the supernova shock. The second mechanism is related to the impact of the supersonic wind termination shock (reverse shock) within the tumultuous inner regions of the ejecta above the proto-neutron star. Our results suggest that isotopes in the mass range 12<=A<=122 that are produced during the freeze-out expansions may be classified in two families. The isotopes of the first family manifest a common mass fraction evolutionary profile, whose specific shape per isotope depends on the characteristic transition between two equilibrium states (equilibrium state transition) during each type of freeze-out expansion. The first family includes the majority of isotopes in this mass range. The second family is limited to magic nuclei and isotopes in their locality, which do not sustain any transition, become nuclear flow hubs, and dominate the final composition. We use exponential and power-law adiabatic profiles to identify dynamic large-scale and small-scale equilibrium patterns among nuclear reactions. (truncated abstract) In addition, we introduce non-monotonic parameterized profiles to probe the impact of the reverse shock and multi-dimensional explosion asymmetries on nucleosynthesis. Cases are shown in which the non-monotonic profiles favor the production of radioactivities. Non-monotonic freeze-out profiles involve longer non-equilibrium nucleosynthesis intervals compared with the exponential and power-law profiles, resulting in mass fraction trends and yield distributions that may not be achieved by the monotonic freeze-out profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Remnants of Binary White Dwarf Mergers

We carry out a comprehensive smooth particle hydrodynamics simulation survey of double-degenerate white dwarf binary mergers of varying mass combinations in order to establish correspondence between initial conditions and remnant configurations. We find that all but one of our simulation remnants share general properties such as a cold, degenerate core surrounded by a hot disk, while our least massive pair of stars forms only a hot disk. We characterize our remnant configurations by the core mass, the rotational velocity of the core, and the half-mass radius of the disk. We also find that some of our simulations with very massive constituent stars exhibit helium detonations on the surface of the primary star before complete disruption of the secondary. However, these helium detonations are insufficiently energetic to ignite carbon, and so do not lead to prompt carbon detonations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

56Ni Production in Double Degenerate White Dwarf Collisions

We present a comprehensive study of white dwarf collisions as an avenue for creating type Ia supernovae. Using a smooth particle hydrodynamics code with a 13-isotope, α-chain nuclear network, we examine the resulting 56Ni yield as a function of total mass, mass ratio, and impact parameter. We show that several combinations of white dwarf masses and impact parameters are able to produce sufficient quantities of 56Ni to be observable at cosmological distances. We find the 56Ni production in double-degenerate white dwarf collisions ranges from sub-luminous to the super-luminous, depending on the parameters of the collision. For all mass pairs, collisions with small impact parameters have the highest likelihood of detonating, but 56Ni production is insensitive to this parameter in high-mass combinations, which significantly increases their likelihood of detection. We also find that the 56Ni dependence on total mass and mass ratio is not linear, with larger mass primaries producing disproportionately more 56Ni than their lower mass secondary counterparts, and symmetric pairs of masses producing more 56Ni than asymmetric pairs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evaluating Systematic Dependencies of Type Ia Supernovae: The Influence of Deflagration to Detonation Density

We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretical framework for a well-controlled statistical study of two-dimensional simulations of thermonuclear supernovae with randomized initial conditions that can, with a particular choice of transition density, produce a similar average and range of Ni-56 masses to those inferred from observations. Within this framework, we utilize a more realistic "simmered" white dwarf progenitor model with a flame model and energetics scheme to calculate the amount of Ni-56 and NSE material synthesized for a suite of simulated explosions in which the transition density is varied in the range 1-3x10^7 g/cc. We find a quadratic dependence of the NSE yield on the log of the transition density, which is determined by the competition between plume rise and stellar expansion. By considering the effect of metallicity on the transition density, we find the NSE yield decreases by 0.055 +/- 0.004 solar masses for a 1 solar metallicity increase evaluated about solar metallicity. For the same change in metallicity, this result translates to a 0.067 +/- 0.004 solar mass decrease in the Ni-56 yield, slightly stronger than that due to the variation in electron fraction from the initial composition. Observations testing the dependence of the yield on metallicity remain somewhat ambiguous, but the dependence we find is comparable to that inferred from some studies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On variations of the brightness of type Ia supernovae with the age of the host stellar population

Recent observational studies of type Ia supernovae (SNeIa) suggest correlations between the peak brightness of an event and the age of the progenitor stellar population. This trend likely follows from properties of the progenitor white dwarf (WD), such as central density, that follow from properties of the host stellar population. We present a statistically well-controlled, systematic study utilizing a suite of multi-dimensional SNeIa simulations investigating the influence of central density of the progenitor WD on the production of Fe-group material, particularly radioactive Ni-56, which powers the light curve. We find that on average, as the progenitor's central density increases, production of Fe-group material does not change but production of Ni-56 decreases. We attribute this result to a higher rate of neutronization at higher density. The central density of the progenitor is determined by the mass of the WD and the cooling time prior to the onset of mass transfer from the companion, as well as the subsequent accretion heating and neutrino losses. The dependence of this density on cooling time, combined with the result of our central density study, offers an explanation for the observed age-luminosity correlation: a longer cooling time raises the central density at ignition thereby producing less Ni-56 and thus a dimmer event. While our ensemble of results demonstrates a significant trend, we find considerable variation between realizations, indicating the necessity for averaging over an ensemble of simulations to demonstrate a statistically significant result.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Evaluating Systematic Dependencies of Type Ia Supernovae: The Influence of Progenitor Ne22 Content on Dynamics

We present a theoretical framework for formal study of systematic effects in Supernovae Type Ia (SN Ia) that utilizes 2-d simulations to implement a form of the deflagration-detonation transition (DDT) explosion scenario. The framework is developed from a randomized initial condition that leads to a sample of simulated SN Ia whose Ni56 masses have a similar average and range to those observed, and have many other modestly realistic features such as the velocity extent of intermediate mass elements. The intended purpose is to enable statistically well-defined studies of both physical and theoretical parameters of the SN Ia explosion simulation. We present here a thorough description of the outcome of the SN Ia explosions produced by our current simulations. A first application of this framework is utilized to study the dependence of the SN Ia on the Ne22 content, which is known to be directly influenced by the progenitor stellar population's metallicity. Our study is very specifically tailored to measure how the Ne22 content influences the competition between the rise of plumes of burned material and the expansion of the star before these plumes reach DDT conditions. This competition controls the amount of material in nuclear statistical equilibrium (NSE) and therefore Ni56 produced by setting the density at which nucleosynthesis takes place during the detonation phase of the explosion. Although the outcome following from any particular ignition condition can change dramatically with Ne22 content, with a sample of 20 ignition conditions we find that the systematic change in the expansion of the star prior to detonation is not large enough to compete with the dependence on initial neutron excess discussed by Timmes, Brown & Truran (2003). (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On Type Ia Supernovae From The Collisions of Two White Dwarfs

We explore collisions between two white dwarfs as a pathway for making Type Ia Supernovae (SNIa). White dwarf number densities in globular clusters allow 10-100 redshift <1 collisions per year, and observations by (Chomiuk et al. 2008) of globular clusters in the nearby S0 galaxy NGC 7457 have detected what is likely to be a SNIa remnant. We carry out simulations of the collision between two 0.6 solar mass white dwarfs at various impact parameters and mass resolutions. For impact parameters less than half the radius of the white dwarf, we find such collisions produce approximately 0.4 solar masses of Ni56, making such events potential candidates for underluminous SNIa or a new class of transients between Novae and SNIa.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Effects of the pep Nuclear Reaction and Other Improvements in the Nuclear Reaction Rate Library on Simulations of the Classical Nova Outburst

We have continued our studies of the Classical Nova outburst by evolving TNRs on 1.25Msun and 1.35Msun WDs (ONeMg composition) under conditions which produce mass ejection and a rapid increase in the emitted light, by examining the effects of changes in the nuclear reaction rates on both the observable features and the nucleosynthesis during the outburst. In order to improve our calculations over previous work, we have incorporated a modern nuclear reaction network into our hydrodynamic computer code. We find that the updates in the nuclear reaction rate libraries change the amount of ejected mass, peak luminosity, and the resulting nucleosynthesis. In addition, as a result of our improvements, we discovered that the pep reaction was not included in our previous studies of CN explosions. Although the energy production from this reaction is not important in the Sun, the densities in WD envelopes can exceed $10^4$ gm cm$^{-3}$ and the presence of this reaction increases the energy generation during the time that the p-p chain is operating. The effect of the increased energy generation is to reduce the evolution time to the peak of the TNR and, thereby, the accreted mass as compared to the evolutionary sequences done without this reaction included. As expected from our previous work, the reduction in accreted mass has important consequences on the characteristics of the resulting TNR and is discussed in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Mixing by Non-linear Gravity Wave Breaking on a White Dwarf Surface

We present the results of a simulation of a wind-driven non-linear gravity wave breaking on the surface of a white dwarf. The ``wind'' consists of H/He from an accreted envelope, and the simulation demonstrates that this breaking wave mechanism can produce a well-mixed layer of H/He with C/O from the white dwarf above the surface. Material from this mixed layer may then be transported throughout the accreted envelope by convection, which would enrich the C/O abundance of the envelope as is expected from observations of novae.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The Neutron Star and Black Hole Initial Mass Function

Using recently calculated models for massive stellar evolution and supernovae coupled to a model for Galactic chemical evolution, neutron star and black hole birth functions (number of neutron stars and black holes as a function of their mass) are determined for the Milky Way Galaxy. For those stars that explode as Type II supernovae, the models give birth functions that are bimodal with peaks at 1.27 and 1.76 M\sun \ and average masses within those peaks of 1.28 and 1.73 M\sun. For those stars that explode as Type Ib there is a narrower spread of remnant masses, the average being 1.32 M\sun, and less evidence for bimodality. These values will be increased, especially in the more massive Type II supernovae, if significant accretion continues during the initial launching of the shock, and the number of heavier neutron stars could be depleted by black hole formation. The principal reason for the dichotomy in remnant masses for Type II is the difference in the presupernova structure of stars above and below 19 M\sun, the mass separating stars that burn carbon convectively from those that produce less carbon and burn radiatively. The Type Ib's and the lower mass group of the Type II's compare favorably with measured neutron star masses, and in particular to the Thorsett et al. (1993) determination of the average neutron star mass in 17 systems; 1.35 $\pm$ 0.27 M\sun. Variations in the exponent of a Salpeter initial mass function are shown not to affect the locations of the two peaks in the distribution function, but do affect their relative amplitudes. Sources of uncertainty, in particular placement of the mass cut and sensitivity to the explosion energy, are discussed, and estimates of the total number of neutron stars and black holes in the Galaxy are given. Accretion induced collapse should give a unique gravitational mass of 1.27

0 citations0 reviewsNo code linkedOpen access