p_excit_3_6_janev Function

public function p_excit_3_6_janev(eb) result(sigma)

Calculates cross section for a proton-Hydrogen impact excitation transition from the state to the state at energy eb

Equation

References

Arguments

Type IntentOptional AttributesName
real(kind=Float64), intent(in) :: eb

Relative collision energy [keV/amu]

Return Value real(kind=Float64)

Cross Section []


Called by

proc~~p_excit_3_6_janev~~CalledByGraph proc~p_excit_3_6_janev p_excit_3_6_janev proc~p_excit_3_8_janev p_excit_3_8_janev proc~p_excit_3_8_janev->proc~p_excit_3_6_janev proc~p_excit_3_9_janev p_excit_3_9_janev proc~p_excit_3_9_janev->proc~p_excit_3_6_janev proc~p_excit_3_7_janev p_excit_3_7_janev proc~p_excit_3_7_janev->proc~p_excit_3_6_janev proc~p_excit_3_10_janev p_excit_3_10_janev proc~p_excit_3_10_janev->proc~p_excit_3_6_janev proc~p_excit_3_janev p_excit_3_janev proc~p_excit_3_janev->proc~p_excit_3_6_janev proc~p_excit_3_janev->proc~p_excit_3_8_janev proc~p_excit_3_janev->proc~p_excit_3_9_janev proc~p_excit_3_janev->proc~p_excit_3_7_janev proc~p_excit_3_janev->proc~p_excit_3_10_janev proc~p_excit_n p_excit_n proc~p_excit_n->proc~p_excit_3_janev proc~p_excit_n_m p_excit_n_m proc~p_excit_n_m->proc~p_excit_n proc~p_excit p_excit proc~p_excit->proc~p_excit_n proc~write_bb_h_h write_bb_H_H proc~write_bb_h_h->proc~p_excit program~generate_tables generate_tables program~generate_tables->proc~write_bb_h_h

Contents

Source Code


Source Code

function p_excit_3_6_janev(eb) result(sigma)
    !+Calculates cross section for a proton-Hydrogen impact excitation transition from
    !+the \(n=3\) state to the \(m=6\) state at energy `eb`
    !+
    !+###Equation
    !+$$ H^+ + H(3) \rightarrow H^+ + H(6) $$
    !+
    !+###References
    !+* Eq. 35 and Table 7 in Ref. 2 [[atomic_tables(module)]]
    !+
    real(Float64), intent(in) :: eb
        !+ Relative collision energy [keV/amu]
    real(Float64)             :: sigma
        !+ Cross Section [\(cm^2\)]
    real(Float64), dimension(6), parameter :: c = [63.494d0, 0.077953d0, &
                                                   0.53461d0, -1.2881d0, &
                                                   11.507d0,   4.3417d0  ]
        !+ Fitting parameters from Table 7 in Ref. 2

    sigma = 1.d-16 * c(1)*(c(2)*exp(-c(3)*eb)/(eb**c(4)) + &
                     exp(-c(5)/eb)*log(1.+c(6)*eb)/eb)

end function p_excit_3_6_janev