Calculates cross section for a proton-Hydrogen impact ionization interaction
from the n
th state at energy eb
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | eb | Relative collision energy [keV/amu] |
||
integer, | intent(in) | :: | n | Initial atomic energy level/state |
Cross Section []
function p_ioniz_n(eb,n) result(sigma)
!+Calculates cross section for a proton-Hydrogen impact ionization interaction
!+from the `n`th state at energy `eb`
!+
!+###Equation
!+ $$H^+ + H(n) \rightarrow H^+ + H^+ + e$$
!+###References
!+* Eq. 40 and Table 8 in Ref. 2 for \(n=1\) [[atomic_tables(module)]]
!+* Eq. 5 and Table 1 in Ref. 3 for \(n \geq 2\) [[atomic_tables(module)]]
real(Float64), intent(in) :: eb
!+ Relative collision energy [keV/amu]
integer, intent(in) :: n
!+ Initial atomic energy level/state
real(Float64) :: sigma
!+ Cross Section [\(cm^2\)]
select case (n)
case (0)
stop
case (1)
sigma = p_ioniz_1_janev(eb)
case (2)
sigma = p_ioniz_2_omullane(eb)
case (3)
sigma = p_ioniz_3_omullane(eb)
case (4)
sigma = p_ioniz_4_omullane(eb)
case DEFAULT
sigma = p_ioniz_5_omullane(eb)*(n/5.d0)**4
end select
end function p_ioniz_n