Calculates the excitation cross section for a neutral Hydrogen atom transitioning from
the state to the state due to a collision an ion with charge q
at energy eb
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | eb | Collision energy [keV] |
||
integer, | intent(in) | :: | q | Ion charge |
Cross Section []
function Aq_excit_1_5_janev(eb, q) result(sigma)
!+Calculates the excitation cross section for a neutral Hydrogen atom transitioning from
!+the \(n=1\) state to the \(m=5\) state due to a collision an ion with charge `q` at energy `eb`
!+
!+###Equation
!+$$ A^{q+} + H(1) \rightarrow A^{q+} + H(5), q \gt 4 $$
!+
!+###References
!+* Page 136 in Ref. 5 [[atomic_tables(module)]]
!+
real(Float64), intent(in) :: eb
!+ Collision energy [keV]
integer, intent(in) :: q
!+ Ion charge
real(Float64) :: sigma
!+ Cross Section [\(cm^2\)]
real(Float64), dimension(6), parameter :: A = [0.56565d0, 67.333d0, &
55.290d0, 2.1595d-4, &
0.081624d0, -2.1971d0 ]
real(Float64) :: Etil, xsi, qf
qf = real(q)
etil = eb/qf
xsi = 2.**(0.5238*(1 - sqrt(2.0/qf)))
sigma = qf*1.d-16*xsi*A(1)*(exp(-A(2)/etil)*log(1 + A(3)*etil)/etil &
+ A(4)*exp(-A(5)*etil)/etil**A(6))
end function Aq_excit_1_5_janev