K. Setoodehnia

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 3He could be produced in situ in solarcomposition envelopes accreted with slow rates (M < ˙ 10−10 M yr−1) by cold (TWD < 107 K) CO WDs, and that convection was triggered by 3He burning before the nova outburst in that case. In addition, we now find that the interplay between the 3He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g. M˙ = 10−10 M yr−1, by the 1.15 M ONe WD with a relatively high initial central temperature, e.g. TWD = 15 × 106 K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. We present detailed nucleosynthesis calculations based on the post-processing technique, and demonstrate in which way much simpler single zone T and ρ trajectories extracted from the multi-zone stellar evolution simulations can be used, in lieu of full multi-zone simulations, to analyse the sensitivity of nova abundance predictions on nuclear reaction rate uncertainties. Trajectories for both CO and ONe nova models for different central temperatures and accretion rates are provided. We compare our nova simulations with observations of novae and pre-solar grains believed to originate in novae.

Constraining nova observables: direct measurements of resonance strengths in 33S(p,\gamma)34Cl

The 33S(p,\gamma)34Cl reaction is important for constraining predictions of certain isotopic abundances in oxygen-neon novae. Models currently predict as much as 150 times the solar abundance of 33S in oxygen-neon nova ejecta. This overproduction factor may, however, vary by orders of magnitude due to uncertainties in the 33S(p,\gamma)34Cl reaction rate at nova peak temperatures. Depending on this rate, 33S could potentially be used as a diagnostic tool for classifying certain types of presolar grains. Better knowledge of the 33S(p,\gamma)34Cl rate would also aid in interpreting nova observations over the S-Ca mass region and contribute to the firm establishment of the maximum endpoint of nova nucleosynthesis. Additionally, the total S elemental abundance which is affected by this reaction has been proposed as a thermometer to study the peak temperatures of novae. Previously, the 33S(p,\gamma)34Cl reaction rate had only been studied directly down to resonance energies of 432 keV. However, for nova peak temperatures of 0.2-0.4 GK there are 7 known states in 34Cl both below the 432 keV resonance and within the Gamow window that could play a dominant role. Direct measurements of the resonance strengths of these states were performed using the DRAGON recoil separator at TRIUMF. Additionally two new states within this energy region are reported. Several hydrodynamic simulations have been performed, using all available experimental information for the 33S(p,\gamma)34Cl rate, to explore the impact of the remaining uncertainty in this rate on nucleosynthesis in nova explosions. These calculations give a range of ~ 20-150 for the expected 33S overproduction factor, and a range of ~ 100-450 for the 32S/33S ratio expected in ONe novae.

Nuclear structure of 30 S and its implications for nucleosynthesis in classical novae

The uncertainty in the 29P(p,gamma)30S reaction rate over the temperature range of 0.1 - 1.3 GK was previously determined to span ~4 orders of magnitude due to the uncertain location of two previously unobserved 3+ and 2+ resonances in the 4.7 - 4.8 MeV excitation region in 30S. Therefore, the abundances of silicon isotopes synthesized in novae, which are relevant for the identification of presolar grains of putative nova origin, were uncertain by a factor of 3. To investigate the level structure of 30S above the proton threshold (4394.9(7) keV), a charged-particle spectroscopy and an in-beam gamma-ray spectroscopy experiments were performed. Differential cross sections of the 32S(p,t)30S reaction were measured at 34.5 MeV. Distorted wave Born approximation calculations were performed to constrain the spin-parity assignments of the observed levels. An energy level scheme was deduced from gamma-gamma coincidence measurements using the 28Si(3He,n-gamma)30S reaction. Spin-parity assignments based on measurements of gamma-ray angular distributions and gamma-gamma directional correlation from oriented nuclei were made for most of the observed levels of 30S. As a result, the resonance energies corresponding to the excited states in 4.5 MeV - 6 MeV region, including the two astrophysically important states predicted previously, are measured with significantly better precision than before. The uncertainty in the rate of the 29P(p,gamma)30S reaction is substantially reduced over the temperature range of interest. Finally, the influence of this rate on the abundance ratios of silicon isotopes synthesized in novae are obtained via 1D hydrodynamic nova simulations.

Spins and parities of astrophysically important {sup 30}S states from {sup 28}Si({sup 3}He, n{gamma}){sup 30}S

The structure of proton-unbound {sup 30}S states strongly determines the thermonuclear {sup 29}P(p,{gamma}){sup 30}S reaction rate at temperatures characteristic of explosive hydrogen burning in classical novae and type I x-ray bursts. Specifically, the rate had been previously predicted to be dominated by two low-lying, unobserved, levels in the E{sub x}= 4.7-4.8 MeV region, with spin and parity assignments of 3{sup +} and 2{sup +}. In recent experimental work, two candidate levels were observed with energies of 4.699 MeV and 4.814 MeV, but no experimental information on their spins and parities was obtained. We have performed an in-beam {gamma}-ray spectroscopy study of {sup 30}S with the {sup 28}Si({sup 3}He, n{gamma}){sup 30}S reaction. By constructing the decay schemes of proton-unbound states with a {gamma}-{gamma} coincidence analysis of their decay {gamma} rays, their J{sup {pi}} values were inferred from a comparison to the known decay schemes of the corresponding mirror states in {sup 30}Si. For the two aforementioned states, our results strongly corroborate the spin-parity assignments assumed in recent evaluations of the {sup 29}P(p,{gamma}){sup 30}S reaction rate.

Production of 26Al in stellar hydrogen-burning environments: spectroscopic properties of states in 27Si

Model predictions of the amount of the radioisotope 26Al produced in hydrogen-burning environments require reliable estimates of the thermonuclear rates for the 26gAl(p,{\gamma})27Si and 26mAl(p,{\gamma})27Si reactions. These rates depend upon the spectroscopic properties of states in 27Si within about 1 MeV of the 26gAl+p threshold (Sp = 7463 keV). We have studied the 28Si(3He,{\alpha})27Si reaction at 25 MeV using a high-resolution quadrupole-dipole-dipole-dipole magnetic spectrograph. For the first time with a transfer reaction, we have constrained J{\pi} values for states in 27Si over Ex = 7.0 - 8.1 MeV through angular distribution measurements. Aside from a few important cases, we generally confirm the energies and spin-parity assignments reported in a recent {\gamma}-ray spectroscopy study. The magnitudes of neutron spectroscopic factors determined from shell-model calculations are in reasonable agreement with our experimental values extracted using this reaction.

Proton‐Rich Sulphur and Nucleosynthesis in Classical Novae

The structure of proton-unbound states in {sup 30}{sub S} and {sup 31}{sub S} is important for determining the {sup 29}P(p,{gamma}){sup 30}S and {sup 30}P(p,{gamma}){sup 31}S reaction rates, which influence explosive hydrogen burning in classical novae. The former reaction rate in this temperature regime had been previously predicted to be dominated by two low-lying, unobserved, J{sup {pi}} = 3{sup +} and 2{sup +} resonances in {sup 30}S. To search for evidence for these levels, the structure of {sup 30}{sub S} was studied using the {sup 32}S(p,t){sup 30}S reaction with a magnetic spectrograph. We provide an update on the status of the ongoing analysis and some preliminary results.

Structure of {sup 30}S with {sup 32}S(p,t){sup 30}S and the thermonuclear {sup 29}P(p,{gamma}){sup 30}S reaction rate

The structure of proton unbound

Study of Astrophysically Important Resonant States in 30S by the 32S(p,t)30S Reaction

^{30}\mathrm{S}

Toward precise QEC values for the superallowed 0+→0 + β decays of T=2 nuclides: The masses of Na20, Al24, P28, and Cl32

states is important for determining the

Measurement of the 17 O ( p , γ ) 18 F reaction rate at astrophysically relevant energies

^{29}\mathrm{P}