Calculates cross section for a proton-Hydrogen impact excitation transition from
the state to 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_excit_1_5_janev(eb) result(sigma)
!+Calculates cross section for a proton-Hydrogen impact excitation transition from
!+the \(n=1\) state to the \(m=5\) state at energy `eb`
!+
!+###Equation
!+$$ H^+ + H(1) \rightarrow H^+ + H(5) $$
!+
!+###References
!+* Eq. 30 and Table 5 in Ref. 2 [[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.2449d0, 3.0826d-4, &
0.31063d0, -2.4161d0, &
0.024664d0, -6.3726d0, &
32.291d0, 0.21176d0 ]
!+ Fitting parameters from Table 5 in Ref. 2
sigma = 1.d-16 * b(1) * (b(2)*exp(-b(3)*eb)/ &
(eb**b(4)+b(5)*eb**b(6)) + &
exp(-b(7)/eb)*log(1.+b(8)*eb)/eb)
end function p_excit_1_5_janev