Function adas.qlpr
def qlpr(iz, nlower, nupper, zimp, mass, energy, excitation=True)
-
Evaluate Lodge-Percival-Richards ion impact excitation cross section.
Parameters
iz
:int
- recombining charge of target
nlower
:int
- lower n for transition
nupper
:int
- upper n for transition
zimp
:int
- charge of projectile
mass
:float
- projectile mass (amu)
energy
:float, array
- set of energies (eV/amu)
excitation
:bool
- if True return the excitation cross section if False return the de-excitation cross section
Returns
oa
:float, array
- set of cross sections (cm**2)
Notes
Calls a fortran based shared object file - not pure python.
References
The Lodge-Percival-Richards ion impact excitation cross section is described in J. Phys. B. (1976) v9, p239 with a scaling to an arbitrary projectile charge following recommendations from Reinnhold, Olsen & Fritsch, Phys. Rev. A (1990) v41, p4837.
Version History
- Martin O'Mullane, 23-05-2022
- First version
Example
Calculate the excitation cross section of Ne10+ 3-4 transition from a deuteron impact.
>>> import adas as adas >>> import numpy as np >>> energy=np.geomspace(100.0, 1e5, 4) >>> adas.qlpr(10, 3, 4, 1, 2.0, energy) array([1.03632792e-19, 7.27556948e-19, 4.08154637e-18, 5.09102351e-18])