Function adas.qh
def qh(epro, ttar, amd, amr, ethra, alfa, ea, oa)
-
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