jaxrts.ionization.solve_BU

jaxrts.ionization.solve_BU(element_list: list[~jaxrts.elements.Element], T_e: ~pint.registry.Quantity, n_e: ~pint.registry.Quantity, ion_number_densities: ~pint.registry.Quantity, continuum_lowering: list[~pint.registry.Quantity] | None = None, chem_pot_ideal: ~pint.registry.Quantity = <Quantity(0, 'electron_volt')>) → tuple[Quantity, Quantity, Quantity][source]

Solve the Bethe-Uhlenbeck equation using Planck-Larkin partition sums.

Thin wrapper around solve_ionization() using bu_balance_term().

Parameters:

element_list (list[Element])
T_e (Quantity) – Temperature.
n_e (Quantity) – Free-electron density (initial estimate; caller should iterate for self-consistency, e.g. via calculate_mean_free_charge_saha()).
ion_number_densities (Quantity)
continuum_lowering (list[Quantity] or None) – Shift of the binding energy. Default None corresponds to no shift.
chem_pot_ideal (Quantity) – Electron chemical potential \(\mu\) (default 0 eV).

Returns:

ionised_number_densities, n_e, Z_mean – See solve_ionization().