*
Cococubed.com


Variable White Dwarfs

Home

Astronomy research
  Software Infrastructure:
     MESA
     FLASH-X
     STARLIB
     MESA-Web
     starkiller-astro
     My instruments
  White dwarf pulsations:
     12C(α,γ) & overshooting
     Probe of 12C(α,γ)16O
     Impact of 22Ne
     Impact of ν cooling
     Variable white dwarfs
     MC reaction rates
     Micronovae
     Novae
  White dwarf supernova:
     Stable nickel production
     Remnant metallicities
     Colliding white dwarfs
     Merging white dwarfs
     Ignition conditions
     Metallicity effects
     Central density effects
     Detonation density
     Tracer particle burning
     Subsonic burning fronts
     Supersonic fronts
     W7 profiles
  Massive stars:
     Pop III with HST/JWST
     Rotating progenitors
     3D evolution to collapse
     MC reaction rates
     Pre-SN variations
  Massive star supernova:
     Yields of radionuclides
     26Al & 60Fe
     44Ti, 60Co & 56Ni
     SN 1987A light curve
     Constraints on Ni/Fe
     An r-process
     Effects of 12C +12C
  Neutron Stars and Black Holes:
     Black Hole spectrum
     Mass Gap with LVK
     Compact object IMF
     He burn neutron stars
  Neutrino Emission:
     Neutrino emission from stars
     Identifying the Pre-SN
     Neutrino HR diagram
     Pre-SN Beta Processes
     Pre-SN neutrinos
  Stars:
     Hypatia catalog
     SAGB stars
     Nugrid Yields I
     He shell convection
     BBFH at 40 years
     γ-rays within 100 Mpc
     Iron Pseudocarbynes
  Pre-Solar Grains:
     C-rich presolar grains
     SiC Type U/C grains
     Grains from massive stars
     Placing the Sun
     SiC Presolar grains
  Chemical Evolution:
     Radionuclides in 2020s
     Zone models H to Zn
     Mixing ejecta
  Thermodynamics, Opacities & Networks
     Radiative Opacity
     Skye EOS
     Helm EOS
     Five EOSs
     Equations of State
     12C(α,γ)16O Rate
     Proton-rich NSE
     Reaction networks
     Bayesian reaction rates
  Verification Problems:
     Validating an astro code
     Su-Olson
     Cog8
     Mader
     RMTV
     Sedov
     Noh
Software Instruments
AAS Journals
   2024 AAS YouTube
   2024 AAS Peer Review Workshops

2024 ASU Energy in Everyday Life
2024 MESA Classroom
Outreach and Education Materials

Other Stuff:
   Bicycle Adventures
   Illustrations
   Presentations



Contact: F.X.Timmes
my one page vitae,
full vitae,
research statement, and
teaching statement.
Seismic Signatures of the 12C(α,γ)16O Reaction Rate in White Dwarf Models with Overshooting (2023)

In this article, we consider the combined effects that overshooting and the 12C(α,γ)16O reaction rate have on variable white dwarf stellar models. We find that carbon-oxygen white dwarf models continue to yield pulsation signatures of the current experimental 12C(α,γ)16O reaction rate probability distribution function when overshooting is included in the evolution. These signatures hold because the resonating mantle region, encompassing $\simeq$0.2 M$_{\odot}$ in a typical $\simeq$0.6M$_{\odot}$ white dwarf model, still undergoes radiative helium burning during the evolution to a white dwarf. Our specific models show two potential low-order adiabatic g-modes, $g_2$ and $g_6$, that signalize the 12C(α,γ)16O reaction rate probability distribution function. Both g-mode signatures induce average relative period shifts of $\Delta P/P = 0.44 \%$ and $\Delta P/P = 1.33\%$ for $g_2$ and $g_6$ respectively. We find that $g_6$ is a trapped mode, and the $g_2$ period signature is inversely proportional to the 12C(α,γ)16O reaction rate. The $g_6$ period signature generally separates the slower and faster reaction rates, and has a maximum relative period shift of $\Delta P/P = 3.45\%$. We conclude that low-order g-mode periods from carbon-oxygen white dwarfs may still serve as viable probes for the 12C(α,γ)16O reaction rate probability distribution function when overshooting is included in the evolution.

image image
Fig 2. - Composition profiles after the no-overshoot (NOV, left) and overshoot (OV, right) models cool to $T_{\rm eff}$ = 10,000 K. Region boundaries are indicated by vertical blue lines.Green curves represent positive $\sigma_i$ \COrate\ reaction rates, grey curves represent negative $\sigma_i$ \COrate\ reaction rates. For both positive and negative $\sigma_i$, the shading grows fainter the further $\sigma$ is from the standard rate ($\sigma=0$; black curve).

image image
Fig. 4 - Kippenhahn diagrams for the NOV (left) and OV (right) $\sigma=0.0$ models. The x-axis is the respective stellar model's age, the left y-axis is the mass coordinate, the right y-axis is the central mass fraction of the isotopes. Bright green areas represent convection, blue shaded regions depict nuclear burning (see colorbar), white areas represent radiation, and yellow-gold areas represent overshooting (right figure). The light green area shows the hydrogen envelope. The solid pink curve is the central $^{4}$He mass fraction, the solid dark blue curve is the central $^{16}$O mass fraction, and the dark yellow curve is the central $^{12}$C mass fraction. The dashed line shows core helium depletion. The evolution was terminated when $\log(L/L_{\odot})>3.0$ for all stellar models, and the figures are plotted until that point. Annotated is the radiative R2 region's edges and widths. An interactive figure is provided in the online version!

image
Fig. 6 - \textit{Top to bottom:} Mass fractions of $^{12}$C and $^{16}$O; $B$ profile; normalized weight function profile for the g$_6$ mode; normalized weight function profile for the g$_2$mode. The left column shows the NOV results and the right column shows the OV results. Both figures are for $\sigma=0.0$, $T_{\rm eff}$=10,000 K. The R1-R4 region boundaries are indicated by dashed vertical lines. An interactive figure is provided in the online version. Its functionality compares the NOV and OV diagrams, as structured in this figure, for any given $\sigma_i$ 12C(α,γ)16O reaction rate. An interactive figure is provided in the online version!

image
Fig. 8 - Adiabatic g$_2$ (left) and $g_6$ (right) mode signatures for the NOV and OV sets, at $T_{\rm eff}$=11,500 K (bright green) and $T_{\rm eff}$=10,000 K (blue) respectively.





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

In this article, we seek signatures of the current experimental $^{12}$C$(\alpha,\gamma)^{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 white dwarf models and the current experimental $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate probability distribution function. We find an average spread in relative period shifts of $\Delta P/P \simeq \pm$ 2% for the identified trapped g-modes over the $\pm$3$\sigma$ uncertainty in the $^{12}$C$(\alpha,\gamma)^{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$(\alpha,\gamma)^{16}$O reaction rate could, in principle, be extractable from the period spectrum of observed variable white dwarfs.

image
Left: $^{12}$C($\alpha,\gamma$) reaction rate ratios, $\sigma_i/\sigma_0$, as a function of temperature. For our models, $\sigma_i$ spans -3.0 to 3.0 in 0.5 step increments, with $\sigma_0$ being the current nominal rate. Negative $\sigma_i$ are gray curves and positive $\sigma_i$ are green curves. The $\pm$1,2,3 $\sigma_i$ curves are labeled. The blue band show the range of temperatures encountered during core and shell He burning. Right: Mass fraction profiles of the evolutionary DAV models resulting from the $^{12}$C($\alpha,\gamma$) reaction rate uncertainties $\sigma_i$ after each model has cooled to T$_{\rm eff}$ = 11,500 K. The nominal $\sigma$=0 reaction rate is the black curve, negative $\sigma_i$ are gray curves and positive $\sigma_i$ are green curves. Solid curves are for $^{12}$C and $^{16}$O, dashed curves are for $^{1}$H and $^{4}$He. The trace isotopes $^{14}$N, $^{20}$Ne, $^{22}$Ne, and $^{56}$Fe are labeled. Key regions and transitions are also labeled.

image
Integer multiples $q$ of the radial wavelength $\lambda_r$ profiles versus radius for the $g_5$, and $g_{10}$ modes, for $\sigma_0$. The trapped $g_5$ mode is shown by the dark blue solid line, and the dotted curves depict the trapped $g_{10}$ mode. Solid black segments depict the width of the regions R1, R2, R3, and R4, as defined by distance between labeled chemical transitions.

image
Top panels: Relative period differences from $\sigma$=0 for the $g_5, g_6$, and $g_{10}$ modes as the evolutionary DAVs cool. The range of T$_{\rm eff}$ in observed DAV WDs is marked, with the vertical dashed black line selecting the T$_{\rm eff}$ = 11,500 K midpoint. Bottom panel: Relative period differences for $g_5$ versus the $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction rate uncertainties $\sigma_i$ at T$_{\rm eff}$ =11,500 K. Scatter points are the raw data values and the curve is a polynomial fit.

image
Relative period differences for model sequences of varying WD masses, shell masses, and classes. The wd_builder DAVs were measured at T$_{\rm eff}$ = 11,500K and the wd_builder DBVs were measured at T$_{\rm eff}$ = 25,000 K. Labeled are the $\pm$3$\sigma$ uncertainties on the experimental $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate in steps of 0.5$\sigma$, the trapped g-mode that most distinctly probes the radiatively-formed, carbon-rich region R2, and the mean relative period difference for each sequence.






image






On The Impact Of 22Ne On The Pulsation Periods Of Carbon-Oxygen White Dwarfs With Helium Dominated Atmospheres (2021)

In this article, we explore changes in the low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_{\odot}$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of 22Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6-7 significant figures of precision. Usually white dwarf models without 22Ne are fit to the observed periods and other properties. The root-mean-square residuals to the ≈ 150-400 s low-order g-mode periods are typically in the range of $\sigma_{\rm rms}$ $\lesssim$ 0.3 s, for a fit precision of $\sigma_{\rm rms}/ P$ $\lesssim$ 0.3%. We find average relative period shifts of $\Delta P/P$ ≈ $\pm$ 0.5% for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of $\Delta P/P$ depending on the specific g-mode, abundance of 22Ne, effective temperature, and mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when 22Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of 22Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including 22Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model fitting processes.


image
DB WD from a 2.1 M$_{\odot}$, Z=0.02, ZAMS model
image
element diffusion in the 0.560 M$_{\odot}$ WD
image
propagation diagram
image
approximation to Brunt-Väisälä frequency
image
mass-radius with markers of key transitions
image
period evolution
image
period changes in the 0.560 M$_{\odot}$ WD with zero or supersolar 22Ne
image
what causes the period changes
image
period changes in the 0.526 and 0.729 M$_{\odot}$ WD with zero and supersolar 22Ne






The Impact of White Dwarf Luminosity Profiles on Oscillation Frequencies (2018)

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$.

In this article we explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021.

This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction.


image



white dwarf structure
image
propagation diagram
image
mode frequency differences
image
weight function shifts
image
white dwarf cooling