Function adas.qh

Evaluate Maxwell averaged total ionisation, excitation and charge exchange rate coefficients. The incident particle is a monoenergetic beam and the target is a Maxwell distribution. Arbitrary relative speeds are allowed.

Parameters

epro : float

incident particle energy (eV/amu)

ttar : float

Maxwell temperature of target particles (eV)

amd : float

donor mass (amu)

amr : float

receiver mass (amu)

ethra : float

threshold energy for the reaction (eV)

alfa : float

power law slope in energy for extrapolation

ea : float, array

set of energies (eV/amu)

oa : float, array

set of cross sections (cm**2)

Returns

qf : float

rate coefficient (cm3 s-1).

Notes

Pure python translation of the adas310 function, qh.for. This version expects that the energies and cross sections, and possible scalings thereof, be input rather than being read from an adf02 file as in the fortran code.

ea, oa, ethra and alfa are usually stored in adf02 format.

The donor and receiver masses are interchangeable thus allowing either particle to be the incident or target particle.

epro and ttar must be scalars.

Version History

• Martin O'Mullane, 07-04-2021
  • First version

Example

Calculate the rate coefficient for a 2keV/amu proton undergoing charge exchange with a background hydrogen neutral at 5keV.

>>> import adas as adas
>>> import numpy as np
>>> file='/home/adas/adas/adf24/scx#h0/scx#h0_ornl#h1.dat'
>>> fulldata=adas.xxdata_24(file)
>>> te=np.asarray([1.0e3, 1.0e4, 5.0e4])
>>> energy=fulldata['eea'][:,0]
>>> sigma=fulldata['scx'][:,0]
>>> adas.qh(2e3, 5e3, 1.0, 1.0, 13.6, 1, energy, sigma)
9.668911637676699e-08