make_co_wd¶
This test builds a carbon-oxygen white dwarf starting from a specified pre main-sequence mass. It does this by running through a series of inlists for different stages of evolution. By default, the first two steps (pre main-sequence and main-sequence through core helium burning) are skipped so that the test will run in a reasonable amount of time, and the test produces a 0.6 \({\rm M}_\odot\) white dwarf from a roughly 3 \({\rm M}_\odot\) progenitor.
Starting from the beginning to produce a white dwarf descended from a
progenitor of different mass or metallicity can be accomplished by
setting the environment variable MESA_RUN_OPTIONAL=t
and editing
the values assigned to initial_mass
or initial_z
in
inlist_common
.
Larger values for initial_mass
can be used to produce more massive
carbon-oxygen white dwarfs, but note that this test is not designed to
produce a realistic initial-to-final-mass relation (IFMR). Instead,
this test is designed to robustly produce usable white dwarf models as
soon as evolution has yielded the desired core structure, and so this
procedure truncates the AGB evolution just before reaching the
thermally-pulsing AGB (TP-AGB) stage, which can be difficult and/or
resource-intensive to model. Skipping the TP-AGB stage means that the
procedure in this test case will generally require a somewhat larger
progenitor mass to produce a given white dwarf mass than what would be
expected from a realistic IFMR.
The five steps in this procedure are:
1. inlist_zams
(optional)¶
This step starts from the pre main-sequence to produce a ZAMS model
with the initial mass and metallicity specified in inlist_common
.
2. inlist_to_end_he_core_burn
(optional)¶
This step evolves from ZAMS through the end of core helium burning.
3. inlist_co_core
¶
This step evolves to near the end of shell helium burning, producing the final C/O core profile in the interior. This step terminates when the He layer mass falls below 0.04 \({\rm M}_\odot\), avoiding the TP-AGB phase that would begin when the He layer becomes smaller. After the AGB envelope is removed in the next step, residual burning in the He layer will reduce its final mass to around 0.01 \({\rm M}_\odot\).
4. inlist_remove_env
¶
This step removes the AGB envelope using the star_relax_to_star_cut
method in src/run_star_extras.f90
. The location of the cut is
specified using x_ctrl(1)
to leave \(10^{-3}\) \({\rm M}_\odot\) of the
hydrogen envelope. After the relaxation ends, residual burning will
reduce the final hydrogen envelope mass to around
\(10^{-4}\) \({\rm M}_\odot\) (for the default 0.6 \({\rm M}_\odot\) white dwarf).
5. inlist_settle
¶
This step turns on diffusion in the young proto-WD model to allow the model to settle into a stratified envelope structure with a pure hydrogen atmosphere. This final step ends when the white dwarf has cooled down to reach a luminosity of 1 \({\rm L}_\odot\).