jaxrts.ipd.ipd_stewart_pyatt

jaxrts.ipd.ipd_stewart_pyatt(Zi: float, ne: Quantity, ni: Quantity, Te: Quantity, Ti: Quantity, Zbar: float | None = None, arb_deg: bool = False) → Quantity[source]

Stewart-Pyatt IPD model, using the formulation which can be found, e.g., in [Calisti et al., 2015] or [Röpke et al., 2019].

Note

The Stewart-Pyatt value is always below both the Debye and ion sphere results.

Parameters:

  • Zi – The charge state of the atom (note that this is the state before the ionization).

  • ne – Electron density. Units of 1/[length]**3.

  • ni – Ion density. Units of 1/[length]**3.

  • Te – The electron temperature.

  • Ti – The ion temperature.

  • Zbar (float, optional) – The average ionization of the plasma. If not given, Zi is assumed to be the average ionization of the plasma.

  • arb_deg (bool, default false) – If True the Debye screening length is evaluated using inverse_screening_length_e(), which includes solving the Fermi integral, rather than the classical value.

Returns:

Quantity – The ipd shift in units of electronvolt.