Calculates cross section for a proton-Hydrogen charge exchange interaction
from the state to the state at energy Erel
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | Erel | Relative collision energy [keV/amu] |
Cross Section []
function p_cx_1_1_adas(Erel) result(sigma)
!+Calculates cross section for a proton-Hydrogen charge exchange interaction
!+from the \(n=1\) state to the \(m=1\) state at energy `Erel`
!+
!+###Equation
!+ $$H^+ + H(1) \rightarrow H(1) + H^+$$
!+###References
!+* Ref. 4 [[atomic_tables(module)]]
real(Float64), intent(in) :: Erel
!+ Relative collision energy [keV/amu]
real(Float64) :: sigma
!+ Cross Section [\(cm^2\)]
real(Float64), dimension(7), parameter :: a = [-3.496092687d2, 4.724931484d2, &
-2.720493064d2, 8.158564625d1, &
-1.339790721d1, 1.138706949d0, &
-3.914774156d-2 ]
real(Float64) :: e, ee, fac, l, p
e = Erel*1.d3
if(e.ge.1.d3) then
ee = max(e,1.0)
fac = 1.d0
else
ee = 1.0d3
fac = Erel**(-0.2)
endif
l = log10(ee)
p = a(1) + a(2)*l + a(3)*l**2.0 + a(4)*l**3.0 + &
a(5)*l**4.0 + a(6)*l**5.0 + a(7)*l**6.0
sigma = fac*(10.d0**p)
end function p_cx_1_1_adas