pisn
This test case evolves an initialy 200 \({\rm M}_\odot\) star from ZAMS untill it undergoes a pair instability supernovae (PISN)
Physical checks
This tracks the central value of (Gamma1 - 4/3) and the point when the pressure weighted integral of Gamma1-4/3 (integral_gamma1) first drops below 0. This value should be 0 or negative otherwise we do not have a PISN but instead PPISN (Renzo et al (2020))
Inlists
This test case has seven parts.
Part 1 (
inlist_make_late_pre_zams
) This creates a 200 \({\rm M}_\odot\) Z=1.6*10-3 modelPart 2 (
inlist_to_zams_header
) This evolves the model to just before the MS startsPart 3 (
inlist_to_end_core_he_burn
) This takes the model up to the end of core helium burningPart 4 (
inlist_remove_envelope_header
) This removes whats left of the hydrogren envelopePart 5 (
inlist_to_end_core_c_burn
) This takes the model up to the end of carbon burningPart 6 (
inlist_convert
) This prepares the model for a PISN, by switching from cell-faced to cell-centered hydrodynamics, turning off rotation, and turning on the AMR mesh controls.Part 7 (
inlist_finish
) This evoles the model through the PISN by taking it through carbon burning, explosive oxygen ignition, and stops once the star reaches a positive total energy.
Last-Updated: 23Jun2021 (MESA 21a860) by rjfarmer