CollisionDB: Fit Function: h_hex_janev

Fit Function: `h_hex_janev_1n`

Fit Function	$σ_{e x}^{1 \to n > 6} (E) = (\frac{6}{A_{9}})^{3} 10^{- 16} A_{1} [\frac{e^{- \frac{A_{2}}{E}} \ln (1 + A_{3} E)}{E} + \frac{A_{4} e^{- A_{5} E}}{E^{A_{6}} + A_{7} E^{A_{8}}}]$
Comments	Python code requires NumPy imported as `np`.

Python

Arguments

name	description	units	type(s)
`E`	impact energy	keV u^-1	`float`, `np.ndarray`
`A1`	fit coefficient	keV cm² u^-1	`float`
`A2`	fit coefficient	keV u^-1	`float`
`A3`	fit coefficient	keV^-1 u	`float`
`A4`	fit coefficient	undef	`float`
`A5`	fit coefficient	keV^-1 u	`float`
`A6`	fit coefficient		`float`
`A7`	fit coefficient	undef	`float`
`A8`	fit coefficient		`float`
`A9`	main quantum number		`float`

Return values

name	description	units	type(s)
`sigma`	cross section	cm²	`float`, `np.ndarray`

Code

def h_hex_janev_1n(E, A1, A2, A3, A4, A5, A6, A7, A8, A9):
    """
    This function calculates proton impact excitation cross sections (in cm2) of 
    H n=1 to n>6 and He+2 of H n=1 to n>6
    param E: requested impact energy in keV/u
    type E: float, np.ndarray
    param Ai: fit coefficient 
    type Ai: float
    """
    sigma = (6/A9)**3 * 1e-16 * A1 * (np.exp(-A2/E) * np.log(1 + A3*E) / E + 
                    A4*np.exp(-A5*E)/ (E**A6 + A7*E**A8))
    return sigma

Fit Function: h_hex_janev_1n

Python

Fit Function: `h_hex_janev_1n`