Calculates cross section for a proton-Hydrogen impact ionization interaction
from the state at energy eb
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | eb | Relative collision energy [keV/amu] |
Cross Section []
function p_ioniz_3_omullane(eb) result(sigma)
!+Calculates cross section for a proton-Hydrogen impact ionization interaction
!+from the \(n=3\) state at energy `eb`
!+
!+###Equation
!+ $$H^+ + H(3) \rightarrow H^+ + H^+ + e$$
!+###References
!+* Eq. 5 and Table 1 in Ref. 3 [[atomic_tables(module)]]
real(Float64), intent(in) :: eb
!+ Relative collision energy [keV/amu]
real(Float64) :: sigma
!+ Cross Section [\(cm^2\)]
real(Float64), dimension(8), parameter :: b = [1.1076d-2, 1.6197d0, &
6.7154d-3, 5.1188d-3, &
1.8549d0, 2.3696d2, &
7.8286d1, 1.0926d23 ]
!+ Fitting Parameters from Table 1 in Ref. 3
real(Float64), parameter :: n2 = 9.d0
real(Float64) :: Ehat
real(Float64) :: p1, p2, p3
Ehat = eb*n2
p1 = b(1)*(n2)**2.0
p2 = Ehat**b(2) * exp(-b(3)*Ehat) / (1.d0 + b(4)*Ehat**b(5))
p3 = (b(6)* exp(-b(7)/Ehat) *log(1.d0 +b(8)*Ehat) ) /Ehat
sigma = 1.0d-16 * p1 * (p2 + p3)
end function p_ioniz_3_omullane