Function `adas.r8yip`

def r8yip(xi, delta)

Evaluates the impact parameter cross-section second bessel integral.

Parameters

xi : float: z|1/kn - 1/kn'|/a0
delta : float: Rc|kn - kn'|

Returns

y : float: Burgess Y value.

Notes

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

This is a support function for ECIP evaluations. It evaluates the Burgess Y-integral, equations C14-C15 in the reference. The inputs are defined (equation C9) in terms of the initial and final electron wave numbers (kn, kn') with a0 being the Bohr radius. Rc is the distance of closest approach.

References

A Burgess and H P Summers, Mon. Not. T. Astron. Soc., 174, p345 (1976)

Version History

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

Example

Arguments are of typical magnitude.

>>> import adas as adas
>>> adas.r8yip(0.0, 0.226)
1.5439344203999998