jaxrts.saha.calculate_mean_free_charge_saha

jaxrts.saha.calculate_mean_free_charge_saha(plasma_state, use_ipd: bool = False, use_chem_pot: bool = False, use_distribution: bool = False, exclude_non_negative_energies: bool = True)[source]

Calculates the mean charge of each ion in a plasma using the Saha-Boltzmann equation.

Parameters:

plasma_state (PlasmaState) – The plasma state object.
use_ipd (bool) – If true, the ipd correction of the plasma state is used to reduce the continuum. Note: the IPD can very much depend on the ionization state. this could result in some circular dependency.
use_chem_pot (bool) – If true, the chemical potential of the plasma state is used, instead of the non-degenerate limiting case.
use_distribution – If true, the ipd is evaluated for each ion species separately.
exclude_non_negative_energies (bool, default = True) – If true, bound states for which the ionization energy is pushed into the continuum are removed from the calculation and do not appear with their Boltzmann factors in the Saha equations.

Returns:

jnp.ndarray – The charge state distribution in fractions.
jnp.ndarray – An array containing the mean charge of each ion in the plasma.