Calculates total cross section for proton-Hydrogen charge exchange interactions from the n
state at energy Erel
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | Erel | Relative collision energy [keV/amu] |
||
integer, | intent(in) | :: | n | Initial atomic energy level/state |
Cross Section []
function p_cx_janev(Erel,n) result(sigma)
!+Calculates total cross section for proton-Hydrogen charge exchange interactions from the `n` state at energy `Erel`
!+
!+###Equation
!+ $$H^+ + H(n) \rightarrow H(m) + H^+$$
!+###References
!+* Eq. 44 and Table 9 in Ref. 2 [[atomic_tables(module)]]
real(Float64), intent(in) :: Erel
!+ Relative collision energy [keV/amu]
integer, intent(in) :: n
!+ Initial atomic energy level/state
real(Float64) :: sigma
!+ Cross Section [\(cm^2\)]
integer :: i
i = min(n,4)
select case (i)
case (0)
stop
case (1)
sigma = p_cx_1_janev(Erel)
case (2)
sigma = p_cx_2_janev(Erel)
case (3)
sigma = p_cx_3_janev(Erel)
case DEFAULT
sigma = p_cx_n_janev(Erel, n)
end select
end function p_cx_janev