jaxrts.ionization.solve_saha

jaxrts.ionization.solve_saha(element_list: list[Element], T_e: Quantity, ion_number_densities: Quantity, continuum_lowering: list[Quantity] | None = None, exclude_non_negative_energies: bool = True) → tuple[Quantity, Quantity, Quantity][source]

Solve the classic (non-degenerate) Saha equation.

Thin wrapper around solve_ionization() using saha_balance_term() to build the balance terms.

Parameters:

  • element_list (list[Element]) – List of elements which should be considered.

  • T_e (Quantity) – Temperature.

  • ion_number_densities (Quantity) – Total number density per species. Has to have the same length as element_list.

  • continuum_lowering (list[Quantity] or None) – Per-species IPD shift (default: zero for all levels). Has to have the same length as element_list.

  • exclude_non_negative_energies (bool) – If True (default), pressure-ionised levels (shifted \(E_b \leq 0\)) are excluded from the balance equations.

Returns:

ionised_number_densities, n_e, Z_mean – See solve_ionization().

Examples using jaxrts.ionization.solve_saha

Saha equation for hot jupiter materials

Saha equation for hot jupiter materials