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 model - Part 2 (
inlist_to_zams_header
) This evolves the model to just before the MS starts - Part 3 (
inlist_to_end_core_he_burn
) This takes the model up to the end of core helium burning - Part 4 (
inlist_remove_envelope_header
) This removes whats left of the hydrogren envelope - Part 5 (
inlist_to_end_core_c_burn
) This takes the model up to the end of carbon burning - Part 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