Calculates the total charge exchange cross section for a Neutral Hydrogen atom
in the state colliding with a fully stripped Boron ion at energy eb
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | eb | Relative collision energy [keV/amu] |
Cross Section []
function B5_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 Boron ion at energy `eb`
!+
!+###Equation
!+$$ B^{5+} + H(2) \rightarrow B^{4+} + 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(10), parameter :: A = [6.603246818d1, -3.072575676d2, &
5.030801019d2, -4.585636345d2, &
2.568666393d2, -9.185150382d1, &
2.100012584d1, -2.964174788d0, &
2.346396110d-1, -7.943766873d-3]
real(Float64) :: e, l, p
e = max(eb*1.d3,10.0)
l = log10(e)
if(e.le.1.d5) then
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
sigma = 10.d0**p
else
sigma = 0.d0
endif
end function B5_cx_2_adas