EOS project update

Ionization of the envelope

Equation of State

  • MESA-OPAL-SCVH eos covers region -11.0 ⇐ logRho ⇐ 3.64; 2.4 ⇐ logT ⇐ 7.8
  • I carved out the region -11.0 ⇐ logRho ⇐ 0; 2.4 ⇐ logT ⇐ 7.8
  • Then extended the EOS to lower temperature using ideal gas, with gamma = 5/3 and mu = 2.21 (0.69 + 0.29 * 4 + 0.02 * 18)

Result

  • All EOS tables are ready

Next Step

Perhaps future project?

  • If we want to run simulation with a lower ambient, we might need an EOS at lower density (-18 ⇐ logRho ⇐ -11)
  • To extend to lower density region (-18 ⇐ logRho ⇐ -11), we need to use PTEH eos, which is also part of MESA-eos.
  • However, the code of PTEH is different from OPAL-SCVH, so some work are needed to get the additional outputs (non-trivial)
  • Current python code I'm using to get the current EOS can be extended to implement the additional region (not too difficult)

Comments

No comments.