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_3_5_janev(eb, q) result(sigma)
!+Calculates the excitation cross section for a neutral Hydrogen atom transitioning from
!+the \(n=3\) state to the \(m=5\) state due to a collision an ion with charge `q` at energy `eb`
!+
!+###Equation
!+$$ A^{q+} + H(3) \rightarrow A^{q+} + H(5), q \gt 3 $$
!+
!+###References
!+* Page 142 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 = [190.59d0, 11.096d0, &
2.9098d0, 0.073307d0, &
0.54177d0, -1.2894d0 ]
real(Float64) :: Etil, qf, xsi
qf = real(q)
etil = eb/qf
xsi = 2.0**(0.397*(1 - sqrt(2.0/qf)))
sigma = qf*1.e-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_3_5_janev