Instant crash for massive planets with Zalmoxis

[lorlo9999]

Problem description

Running simulations with Zalmoxis turned on crashes instantly for more massive planets. It runs fine for a planet (TOI-431) configuration file with mass set to 1.5 Earth masses but crashes instantly if the same planet is run with a mass of 3 Earth masses.

Evidence

I attach here the configuration files and the logs for the normal case ('full') and for the lower mass version ('_lowmass'). The breaking difference happens in this step of the log:

Broken:

[ INFO ] Using Zalmoxis to solve for interior structure
[ INFO ] Total target planet mass (dry mass + volatiles): 1.833404e+25 kg with EOS: core=Seager2007:iron, mantle=WolfBower2018:MgSiO3, ice=(none)
[ INFO ] Volatile mass: 0.0 kg
[ INFO ] Target planet mass (dry mass): 1.833404e+25 kg
[ WARNING ] RTPress100TPa EOS folder not found at /home5/s4220110/FWL_DATA/EOS_material_properties/EOS_RTPress_melt_100TPa. RTPress100TPa:MgSiO3 will fail if used. Run proteus get to download EOS data.
[ WARNING ] RTPress100TPa melt Cp table not found at /home5/s4220110/FWL_DATA/EOS_material_properties/EOS_RTPress_melt_100TPa/heat_capacity_melt.dat. Adiabatic mode will fall back to constant Cp for this EOS.
[ INFO ] Found solution for interior structure with Zalmoxis
[ INFO ] Interior (dry calculated mass) mass: 0.0 kg or approximately 0.00 M_earth
[ INFO ] Interior radius: 6.42e+07 m or 10.14 R_earth
[ INFO ] Core radius: 0.00e+00 or 0.00 R_earth
[ INFO ] Core-mantle boundary mass: 0.00e+00 kg
[ INFO ] Mantle density at the core-mantle boundary: 8.46e+03 kg/m^3
[ INFO ] Core density at the core-mantle boundary: 8.30e+03 kg/m^3
[ INFO ] Pressure at the core-mantle boundary: 3.26e+09 Pa
[ INFO ] Pressure at the center: 3.26e+09 Pa
[ INFO ] Average density: 0.00e+00 kg/m^3
[ INFO ] Core-mantle boundary mass fraction: nan
[ INFO ] Core radius fraction: 0.0000
[ INFO ] Inner mantle radius fraction: 0.0000
[ INFO ] Overall Convergence Status: False with Pressure: False, Density: True, Mass: False

Working version:

[ INFO ] Using Zalmoxis to solve for interior structure
[ INFO ] Total target planet mass (dry mass + volatiles): 8.958e+24 kg with EOS: core=Seager2007:iron, mantle=WolfBower2018:MgSiO3, ice=(none)
[ INFO ] Volatile mass: 0.0 kg
[ INFO ] Target planet mass (dry mass): 8.958e+24 kg
[ WARNING ] RTPress100TPa EOS folder not found at /home5/s4220110/FWL_DATA/EOS_material_properties/EOS_RTPress_melt_100TPa. RTPress100TPa:MgSiO3 will fail if used. Run proteus get to download EOS data.
[ WARNING ] RTPress100TPa melt Cp table not found at /home5/s4220110/FWL_DATA/EOS_material_properties/EOS_RTPress_melt_100TPa/heat_capacity_melt.dat. Adiabatic mode will fall back to constant Cp for this EOS.
[ INFO ] Found solution for interior structure with Zalmoxis
[ INFO ] Interior (dry calculated mass) mass: 8.952018349575473e+24 kg or approximately 1.50 M_earth
[ INFO ] Interior radius: 7.42e+06 m or 1.17 R_earth
[ INFO ] Core radius: 3.75e+06 or 0.59 R_earth
[ INFO ] Core-mantle boundary mass: 2.97e+24 kg
[ INFO ] Mantle density at the core-mantle boundary: 5.72e+03 kg/m^3
[ INFO ] Core density at the core-mantle boundary: 1.25e+04 kg/m^3
[ INFO ] Pressure at the core-mantle boundary: 1.98e+11 Pa
[ INFO ] Pressure at the center: 6.05e+11 Pa
[ INFO ] Average density: 5.22e+03 kg/m^3
[ INFO ] Core-mantle boundary mass fraction: 0.332
[ INFO ] Core radius fraction: 0.5051
[ INFO ] Inner mantle radius fraction: 0.5051
[ INFO ] Overall Convergence Status: True with Pressure: True, Density: True, Mass: True

My computer

Habrok cluster with Python 3.12

Relevant people

@timlichtenberg

toi431-full.log
TOI-431-full.txt

toi431-lowmass.log
TOI-431-lowmass.txt

[timlichtenberg]

Thanks! Will check out asap. In the meantime you can revert back the main version one step and/or run in "self" mode.