Calculates the total charge exchange cross section for a Neutral Hydrogen atom
in the state colliding with a fully stripped Carbon ion at energy eb
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | eb | Relative collision energy [keV/amu] |
Cross Section []
function C6_cx_2_adas(eb) result(sigma)
!+ Calculates the total charge exchange cross section for a Neutral Hydrogen atom
!+in the \(n=2\) state colliding with a fully stripped Carbon ion at energy `eb`
!+
!+###Equation
!+$$ C^{6+} + H(2) \rightarrow C^{5+} + H^+ $$
!+
!+###References
!+* Ref. 4 [[atomic_tables(module)]]
real(Float64), intent(in) :: eb
!+ Relative collision energy [keV/amu]
real(Float64) :: sigma
!+ Cross Section [\(cm^2\)]
real(Float64), dimension(11), parameter :: A = [9.151879441d5, -2.134573133d6, &
2.223792624d6, -1.362648703d6, &
5.438401343d5, -1.477110500d5, &
2.764972254d4, -3.522105245d3, &
2.921934171d2, -1.425552507d1, &
3.106007048d-1 ]
real(Float64) :: e, l, p, m
e = max(eb*1.d3,1.5d3)*2.0**2
if(eb.le.700) then
l = log10(e)
p = A(1) + A(2)*l + A(3)*l**2 + A(4)*l**3 + &
A(5)*l**4 + A(6)*l**5 + A(7)*l**6 + &
A(8)*l**7 + A(9)*l**8 + A(10)*l**9 + A(11)*l**10
else
l = log10(700*1.d3*2.0**2)
p = A(1) + A(2)*l + A(3)*l**2 + A(4)*l**3 + &
A(5)*l**4 + A(6)*l**5 + A(7)*l**6 + &
A(8)*l**7 + A(9)*l**8 + A(10)*l**9 + A(11)*l**10
m = A(2) + 2*A(3)*l + 3*A(4)*l**2 + &
4*A(5)*l**3 + 5*A(6)*l**4 + 6*A(7)*l**5 + &
7*A(8)*l**6 + 8*A(9)*l**7 + 9*A(10)*l**8 + 10*A(11)*l**9
p = p + m*(log10(e) - l)
endif
sigma = 10.d0**p
end function C6_cx_2_adas