DataSet D108060

\[\mathrm{e}^- + \mathrm{H} \; 1s \rightarrow \mathrm{H} \; n=2 + \mathrm{e}^-\]


Process EEX: Excitation
Data type cross section | uploaded on 2023-02-08
Comment Fit function and data were taken from Janev et al. p6. Error: E< 100 eV -> 15%, E>= 100 eV -> 10%.

Methodexperimental
Columns
  1. E /eV
  2. sigma /cm2
Threshold10.200 eV
Ref
DataDownload (data from fit)

Fitted Data

Fit Function
Details
\[\sigma_{ex}^{n=1 \rightarrow 2} =\begin{cases} \frac{5.984 \times 10^{-16}}{IX} \Bigg[ \sum_{j=1}^5 \frac{A_j}{X^{j-1}} + A_6 \ln X\Bigg]; X=\frac{E}{I}; \; (E > 12.23) \\ A_7; \; (E > 11.56, E < 12.23) \\ 10^{-16} \big[A_8 + A_9(E - I) \big]; \; (E < 11.56). \end{cases}\]
Python
def h_eex_janev_12(E, A1, A2, A3, A4, A5, A6, A7, A8, A9, I):
    """
    This function calculates electron impact excitation cross sections (in cm2) of 
    H n=1 to n=2
    param E: requested electron-impact energy in eV
    type E: np.ndarray
    param Ai: fit coefficient 
    type Ai: float
    param I: threshold energy in eV
    type I: float
    """
    sigma = (5.984e-16/E)*(A1/(E/I)**0 + A2/(E/I)**1 + A3/(E/I)**2 + 
                           A4/(E/I)**3 + A5/(E/I)**4 + A6*np.log(E/I))
    ind_below = np.where(E < 12.23)
    for ind in ind_below[0]:
        if E[ind] < 11.56:
            sigma[ind] = 1e-16*(A8 + A9*(E[ind]-I))
        else:
            sigma[ind] = A7
    return sigma
Fit Precision1.3 %
Fit Coefficients
A1
 1.418e+00
A2
-2.088e+01
A3
 4.973e+01
A4
-4.625e+01
A5
 1.744e+01
A6
 4.498e+00
A7
 5.025e-17
A8
 2.550e-01
A9
 1.865e-01
I
 1.020e+01
x-range 10.2 –