jaxrts.ionization.calculate_mean_free_charge_saha

jaxrts.ionization.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]

Mean ionization charge per species via Saha-Boltzmann.

Iterates up to 6 times for self-consistency between the ionization state and the IPD / chemical potential.

Parameters:

  • plasma_state (PlasmaState)

  • use_ipd (bool) – Apply IPD continuum-lowering. Note: the IPD often depends on the ionization state, creating a circular dependency that the iteration resolves.

  • use_chem_pot (bool) – Use the generalised Saha equation with chemical potential.

  • use_distribution (bool) – Evaluate IPD per-species rather than from the mean state.

  • exclude_non_negative_energies (bool) – Exclude pressure-ionised bound states.

Returns:

  • charge_distribution (jnp.ndarray) – The charge state distribution in fractions.

  • Z_mean (jnp.ndarray) – An array containing the mean charge of each ion in the plasma.

See also

jaxrts.ionization.solve_saha

Function used to solve the saha equation

Examples using jaxrts.ionization.calculate_mean_free_charge_saha

Saha Ionization Map

Saha Ionization Map