Calculates an array of cross sections for proton-Hydrogen impact ionization interactions
from the n = 1..n_max
state at energy eb
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | eb | Relative collision energy [keV/amu] |
||
integer, | intent(in) | :: | n_max | Number of initial atomic energy level/state |
Array of cross sections where the index refers to the n
'th state []
function p_ioniz(eb,n_max) result(sigma)
!+Calculates an array of cross sections for proton-Hydrogen impact ionization interactions
!+from the n = 1..`n_max` state at energy `eb`
!+
!+###Equation
!+ $$H^+ + H(n=1..n_{max}) \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_max
!+ Number of initial atomic energy level/state
real(Float64), dimension(n_max) :: sigma
!+ Array of cross sections where the index refers to the `n`'th state [\(cm^2\)]
integer :: i
do i=1,n_max
sigma(i) = p_ioniz_n(eb,i)
enddo
end function p_ioniz