Function `adas.r8necip`

def r8necip(iz, xi, zeta, te, alfred=False)

Evaluates the ECIP (exchange classical impact parameter) electron impact ionization rate.

Parameters

iz : int: recombined ion charge.
xi : float: effective ionization potential (Rydbergs).
zeta : float: effective number of equivalent electrons.
te : float: evaluate rate at these temperatures (Kelvin).
alfred : bool: if set to True also return the scaled 3-body recombination coefficient. Default is False.

Returns

rate : float: ionization rate coefficient (cm3/s)
r3b : float: optionally return the 3-body recombination coefficient (cm6/s) if alfred is True.

Notes

Calls a fortran based shared object file - not pure python.

References

The primary references are two internal laboratory reports and a journal article: A Burgess, AERE R-4818 (1964) H P Summers, Appleton Laboratory IM-367 (1974) A Burgess and H P Summers, Mon. Not. T. Astron. Soc., 174, p345 (1976)

Version History

Martin O'Mullane, 02-08-2022
- First version

Example

Calculate the ionization rate from n=3 level of neutral hydrogen. In the example the ionization potential and n=3 energy are in cm-1 but must be converted to Rydbergs for use.

>>> import adas as adas
>>> xi = (109679.0 - 97492.3) / 109737.0
>>> adas.r8necip(0, xi, 1.0, [5e4, 10e4, 100e4])
array([3.13826466e-07, 3.85232294e-07, 3.19024094e-07])