lir_achem.compute_coefficients_mitra

Here we compute the chesmitry coefficients as presented in Mitra & Rowe, 1974

Functions

B_coeff(M)	Coefficient B from Mitra & Rowe, 1974
O2m_to_Xm(O2, M, O3)
O2m_to_e(N2, O2, O, H, Tn)	O2- -> e
O2p_to_04p(M, O2, Tn)	O2+ -> O4+ From Pavlov, 2014, reaction 36 of Table 6
O2p_to_NOp(NO, N)	O2+ -> NO+ From Pavlov, 2014, reactions 28 and 29 of Table 6
O4p_to_O2p(M, O, Tn)	O4+ -> O2+ From Pavlov, 2014, reactions 37, 44 and 45 of Table 6 We assume that [O2(1Dg)] is 1e10 cm-3 at our altitudes (Crutzen et al., 1970 (Figure 1))
O4p_to_Yp(H2O)
Xm_and_O4p()
Xm_to_e(M)	X- -> e
alpha_i()	Recombination of each positive ion to the negative ones
e_and_04p(Tn)	Recombination coefficient for O4+ and e- From Pavlov, 2014, reaction 110 of Table 6
e_and_NOp(Tn)	Recombination coefficient of NO+ and e- From Pavlov, 2014, reaction 30 of Table 6 Tn is less than 1200 K
e_and_O2p(Tn)	Recombination coefficent of O2+ and e- From Pavlov, 2014, reaction 27 of Table 6 We assume Tn is less than 1200 K
e_and_Yp()	Recombination coefficient
e_to_O2m(Tn, O2, N2)	e- -> O2-
get_mr_coeffs(n_here, e_here)

lir_achem.compute_coefficients_mitra.B_coeff(M)

Coefficient B from Mitra & Rowe, 1974

lir_achem.compute_coefficients_mitra.alpha_i()

Recombination of each positive ion to the negative ones

From Mitra & Rowe, 1972

lir_achem.compute_coefficients_mitra.Xm_to_e(M)

X- -> e

Parameters:

M – [O2] + [N2]

From Breakall, 1982: Private communication in Burns et al., 1991

lir_achem.compute_coefficients_mitra.O2m_to_e(N2, O2, O, H, Tn)

O2- -> e

This comes from reactions 26, 27, 31, 33 and 34 of Table 10 of Pavlov, 2014 We assume that [O2(1Dg)] is 1e10 cm-3 at our altitudes (Crutzen et al., 1970 (Figure 1))

Parameters:

• N2 – Density of N2 (cm-3)

• O2 – Density of N2 (cm-3)

• O – Density of O (cm-3)

• H – Density of H (cm-3)

• Tn – Neutral density (in K)

lir_achem.compute_coefficients_mitra.e_to_O2m(Tn, O2, N2)

e- -> O2-

From reactions 1 and 2 in Table 10 of Pavlov, 2014

Parameters:

• Tn – Neutrals temperature (in K)

• O2 – O2 density (cm-3)

• N2 – N2 density (cm-2)

lir_achem.compute_coefficients_mitra.O2m_to_Xm(O2, M, O3)

Parameters:

• O2 – O2 density (cm-3)

• M – O2 + N2 density (cm-3)

• O3 – O3 density (cm-3)

From Mitra, 1975

lir_achem.compute_coefficients_mitra.e_and_NOp(Tn)

Recombination coefficient of NO+ and e- From Pavlov, 2014, reaction 30 of Table 6 Tn is less than 1200 K

Parameters:

Tn – Neutral temperature (in K)

lir_achem.compute_coefficients_mitra.e_and_O2p(Tn)

Recombination coefficent of O2+ and e- From Pavlov, 2014, reaction 27 of Table 6 We assume Tn is less than 1200 K

Parameters:

Tn – Neutral temperature (in K)

lir_achem.compute_coefficients_mitra.e_and_04p(Tn)

Recombination coefficient for O4+ and e- From Pavlov, 2014, reaction 110 of Table 6

Parameters:

Tn – Neutrals temperature (in K)

lir_achem.compute_coefficients_mitra.e_and_Yp()

Recombination coefficient

lir_achem.compute_coefficients_mitra.O2p_to_04p(M, O2, Tn)

O2+ -> O4+ From Pavlov, 2014, reaction 36 of Table 6

Parameters:

• M – O2 + N2 densities (in cm-3)

• O2 – O2 density (in cm-3)

• Tn – Neutral temperature (in K)

lir_achem.compute_coefficients_mitra.O4p_to_O2p(M, O, Tn)

O4+ -> O2+ From Pavlov, 2014, reactions 37, 44 and 45 of Table 6 We assume that [O2(1Dg)] is 1e10 cm-3 at our altitudes (Crutzen et al., 1970 (Figure 1))

Parameters:

• M – O2 + N2 densities (in cm-3)

• O – O density (in cm-3)

• Tn – Neutrals temperature (in K)

lir_achem.compute_coefficients_mitra.O2p_to_NOp(NO, N)

O2+ -> NO+ From Pavlov, 2014, reactions 28 and 29 of Table 6

Parameters:

• NO – NO density (in cm-3)

• N – N density (in cm-3)