p_ioniz_5_omullane Function

  1. atomic_tables.f90
  2. atomic_tables
  3. p_ioniz_5_omullane

public function p_ioniz_5_omullane(eb) result(sigma)

Calculates cross section for a proton-Hydrogen impact ionization interaction from the state at energy eb

Equation

References

Arguments

TypeIntentOptionalAttributesName
real(kind=Float64), intent(in) :: eb

Relative collision energy [keV/amu]

Return Value real(kind=Float64)

Cross Section []


Called by

proc~~p_ioniz_5_omullane~~CalledByGraph proc~p_ioniz_5_omullane p_ioniz_5_omullane proc~p_ioniz_n p_ioniz_n proc~p_ioniz_n->proc~p_ioniz_5_omullane proc~p_ioniz p_ioniz proc~p_ioniz->proc~p_ioniz_n proc~write_bb_h_h write_bb_H_H proc~write_bb_h_h->proc~p_ioniz program~generate_tables generate_tables program~generate_tables->proc~write_bb_h_h

Contents

Source Code

p_ioniz_5_omullane

Source Code

function p_ioniz_5_omullane(eb) result(sigma)
    !+Calculates cross section for a proton-Hydrogen impact ionization interaction
    !+from the \(n=5\) state at energy `eb`
    !+
    !+###Equation
    !+ $$H^+ + H(5) \rightarrow H^+ + H^+ + e$$
    !+###References
    !+* Eq. 5 and Table 1 in Ref. 3 [[atomic_tables(module)]]
    real(Float64), intent(in) :: eb
        !+ Relative collision energy [keV/amu]
    real(Float64)             :: sigma
        !+ Cross Section [\(cm^2\)]
    real(Float64), dimension(8), parameter :: b = [1.1297d-2, 1.8685d0,  &
                                                   1.5038d-2, 1.1195d-1, &
                                                   1.0538d0,  8.6096d2,  &
                                                   8.9939d1,  1.9249d4 ]
        !+ Fitting Parameters from Table 1 in Ref. 3

    real(Float64), parameter :: n2 = 25.d0
    real(Float64) :: Ehat
    real(Float64) :: p1, p2, p3

    Ehat = eb*n2
    p1 = b(1)*(n2)**2.0
    p2 = Ehat**b(2) * exp(-b(3)*Ehat) / (1.d0 + b(4)*Ehat**b(5))
    p3 = (b(6)* exp(-b(7)/Ehat) *log(1.d0  +b(8)*Ehat) ) /Ehat
    sigma = 1.0d-16 * p1 * (p2 + p3)

end function p_ioniz_5_omullane