[Mesa-users] PPISN evolution with large nuclear network: timestep shrinks in late burning stages

[翟朔]

Dear all,



I am currently simulating the evolution of pulsational pair-instability supernova (PPISN) progenitors using MESA version r11701. 

My setup is based on the same inlists and run_star_extras.f90 modifications as used in Farmer et al. (2020). 

I employ a large nuclear reaction network (280 isotopes) to resolve detailed nucleosynthesis during advanced burning stages.

However, I’ve encountered a critical issue: during late evolution (log Tc > 9.4, log ρc ~ 8), the timestep drops sharply (e.g., log dt ≈ -21), 

causing the model to either freeze or progress extremely slowly. 

This does not occur when using a smaller reaction network under otherwise identical conditions.


In addition to standard timestep and mesh controls, I have tried tuning the following parameters:
- Pextra_factor
- opacity_max
- mesh_delta_coeff
- max_dq
- min_D_mix
Despite these adjustments, the timestep collapse persists when using the large reaction network.

To facilitate feedback, I’ve uploaded a compressed test folder to Google Drive. It includes:
- the custom nuclear reaction network,
- initial metallicity settings,
- relevant inlists,


- history.data and log.txt from the problematic run.

You can download the test setup here:  

https://drive.google.com/file/d/15fuWWYKaEFllDBMEoFTPU1WhEM3IsNg2/view?usp=sharing



If further output (e.g., selected profiles) would be helpful, I’m happy to share.  
Any suggestions on how to stabilize timestep behavior in this context would be greatly appreciated.





Thanks and regards,

Shuo Zhai,

The National Astronomical Observatories of the Chinese Academy of Sciences
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20250623/3f9d2e3f/attachment.htm>