[mesa-users] EoS & partial ionization

[Ken S. <ke...@as...>]

Hi all.  I've got a question about how the MESA EoS handles partial
ionization.  In particular, are the changes in free electron fraction
taken into account when the EoS returns, e.g., chiT and chiRho?  As a
test, here are the results for two calls to the EoS for X_12C = X_16O
= 0.5 and rho = 2e-8 g/cm^3:

T = 6.5e4 K
P = 9.022e4 dyne/cm^2
free_e = 0.3463
chiT = 2.499
mu = 1.745

T = 6.4e4 K
P = 8.577e4 dyne/cm^2
free_e = 0.3364
chiT = 2.482
mu = 1.745

A crude calculation from these two calls gives dlnP/dlnT = 3.26, which
is significantly different from the returned value of chiT = 2.5.
Also, the mean molecular weight that is returned seems to assume full
ionization; it should be 2.39 for the first call and 2.44 for the
second.


Thanks,
Ken