p_ioniz Function

public function p_ioniz(eb, n_max) result(sigma)

Calculates an array of cross sections for proton-Hydrogen impact ionization interactions from the n = 1..n_max state at energy eb

Equation

References

Arguments

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

Relative collision energy [keV/amu]

integer, intent(in) :: n_max

Number of initial atomic energy level/state

Return Value real(kind=Float64), dimension(n_max)

Array of cross sections where the index refers to the n'th state []


Calls

proc~~p_ioniz~~CallsGraph proc~p_ioniz p_ioniz proc~p_ioniz_n p_ioniz_n proc~p_ioniz->proc~p_ioniz_n proc~p_ioniz_5_omullane p_ioniz_5_omullane proc~p_ioniz_n->proc~p_ioniz_5_omullane proc~p_ioniz_4_omullane p_ioniz_4_omullane proc~p_ioniz_n->proc~p_ioniz_4_omullane proc~p_ioniz_3_omullane p_ioniz_3_omullane proc~p_ioniz_n->proc~p_ioniz_3_omullane proc~p_ioniz_2_omullane p_ioniz_2_omullane proc~p_ioniz_n->proc~p_ioniz_2_omullane proc~p_ioniz_1_janev p_ioniz_1_janev proc~p_ioniz_n->proc~p_ioniz_1_janev

Called by

proc~~p_ioniz~~CalledByGraph proc~p_ioniz p_ioniz 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


Source Code

function p_ioniz(eb,n_max) result(sigma)
    !+Calculates an array of cross sections for proton-Hydrogen impact ionization interactions
    !+from the n = 1..`n_max` state at energy `eb`
    !+
    !+###Equation
    !+ $$H^+ + H(n=1..n_{max}) \rightarrow H^+ + H^+ + e$$
    !+###References
    !+* Eq. 40 and Table 8 in Ref. 2 for \(n=1\) [[atomic_tables(module)]]
    !+* Eq. 5 and Table 1 in Ref. 3 for \(n \geq 2\) [[atomic_tables(module)]]
    real(Float64), intent(in)       :: eb
        !+ Relative collision energy [keV/amu]
    integer, intent(in)             :: n_max
        !+ Number of initial atomic energy level/state
    real(Float64), dimension(n_max) :: sigma
        !+ Array of cross sections where the index refers to the `n`'th state [\(cm^2\)]

    integer :: i

    do i=1,n_max
        sigma(i) = p_ioniz_n(eb,i)
    enddo

end function p_ioniz