Aq_excit_2_5_janev Function

public function Aq_excit_2_5_janev(eb, q) result(sigma)

Calculates the excitation cross section for a neutral Hydrogen atom transitioning from the state to the state due to a collision an ion with charge q at energy eb

Equation

References

Arguments

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

Collision energy [keV]

integer, intent(in) :: q

Ion charge

Return Value real(kind=Float64)

Cross Section []


Called by

proc~~aq_excit_2_5_janev~~CalledByGraph proc~aq_excit_2_5_janev Aq_excit_2_5_janev proc~aq_excit_2_10_janev Aq_excit_2_10_janev proc~aq_excit_2_10_janev->proc~aq_excit_2_5_janev proc~aq_excit_2_9_janev Aq_excit_2_9_janev proc~aq_excit_2_9_janev->proc~aq_excit_2_5_janev proc~aq_excit_2_8_janev Aq_excit_2_8_janev proc~aq_excit_2_8_janev->proc~aq_excit_2_5_janev proc~aq_excit_2_7_janev Aq_excit_2_7_janev proc~aq_excit_2_7_janev->proc~aq_excit_2_5_janev proc~aq_excit_2_6_janev Aq_excit_2_6_janev proc~aq_excit_2_6_janev->proc~aq_excit_2_5_janev proc~aq_excit_2_janev Aq_excit_2_janev proc~aq_excit_2_janev->proc~aq_excit_2_5_janev proc~aq_excit_2_janev->proc~aq_excit_2_10_janev proc~aq_excit_2_janev->proc~aq_excit_2_9_janev proc~aq_excit_2_janev->proc~aq_excit_2_8_janev proc~aq_excit_2_janev->proc~aq_excit_2_7_janev proc~aq_excit_2_janev->proc~aq_excit_2_6_janev proc~aq_excit_n Aq_excit_n proc~aq_excit_n->proc~aq_excit_2_janev proc~aq_excit_n_m Aq_excit_n_m proc~aq_excit_n_m->proc~aq_excit_n proc~aq_excit Aq_excit proc~aq_excit->proc~aq_excit_n proc~write_bb_h_aq write_bb_H_Aq proc~write_bb_h_aq->proc~aq_excit program~generate_tables generate_tables program~generate_tables->proc~write_bb_h_aq

Contents

Source Code


Source Code

function Aq_excit_2_5_janev(eb, q) result(sigma)
    !+Calculates the excitation cross section for a neutral Hydrogen atom transitioning from
    !+the \(n=2\) state to the \(m=5\) state due to a collision an ion with charge `q` at energy `eb`
    !+
    !+###Equation
    !+$$ A^{q+} + H(2) \rightarrow A^{q+} + H(5), q \gt 3 $$
    !+
    !+###References
    !+* Page 138 in Ref. 5 [[atomic_tables(module)]]
    !+
    real(Float64), intent(in) :: eb
        !+ Collision energy [keV]
    integer, intent(in)       :: q
        !+ Ion charge
    real(Float64)             :: sigma
        !+ Cross Section [\(cm^2\)]

    real(Float64), dimension(6), parameter :: A = [18.264d0, 18.973d0,   &
                                                   2.9056d0, 0.013701d0, &
                                                   0.31711d0, -1.4775d0 ]
    real(Float64) :: etil, qf, xsi

    qf = real(q)
    etil = eb/qf
    xsi = 2.0**(0.5238*(1 - sqrt(2.0/qf)))
    sigma = qf*1.e-16*xsi*A(1)*(exp(-A(2)/etil)*log(1 + A(3)*etil)/etil &
          + A(4)*exp(-A(5)*etil)/etil**A(6))

end function Aq_excit_2_5_janev