Next Article in Journal
Measurement of Phthalates in Settled Dust in University Dormitories and Its Implications for Exposure Assessment
Previous Article in Journal
Machine-Learning-Based Downscaling of Hourly ERA5-Land Air Temperature over Mountainous Regions
Previous Article in Special Issue
Seismo Ionospheric Anomalies around and over the Epicenters of Pakistan Earthquakes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling of Energy-Dependent Electron Interactions in the Earth’s Mesosphere

by
Laurence Campbell
1,* and
Michael J. Brunger
1,2,†
1
College of Science and Engineering, Flinders University, Bedford Park, SA 5042, Australia
2
Institute of Actuarial Science and Data Analytics, UCSI University, Kuala Lumpur 56000, Malaysia
*
Author to whom correspondence should be addressed.
Deceased author.
Atmosphere 2023, 14(4), 611; https://doi.org/10.3390/atmos14040611
Submission received: 23 February 2023 / Revised: 9 March 2023 / Accepted: 16 March 2023 / Published: 23 March 2023
(This article belongs to the Special Issue Structure and Dynamics of Mesosphere and Lower Thermosphere)

Abstract

:
Electrons are produced in the Earth’s quiet nighttime mesosphere by ionization by cosmic rays and ionization of NO by Lyman- α  radiation. They are removed by attachment or recombination processes that are usually assumed in modelling to occur at the ambient temperature. However, the electrons have initial energies that are much higher than at thermal equilibrium, and so must have a range of energies as they progress towards equilibrium via interactions with atoms and molecules. As attachment and recombination rates are dependent on the electron energy, it is possible that modelling that considers the actual energy of the electrons will give different results to those based on assuming that the electrons are at the ambient temperature. In this work, starting with electrons at a higher initial energy, the detailed electron interactions (including elastic scattering and vibrational excitation of molecules) are tracked in a time-step simulation. This simulation is implemented by treating electrons in subranges of the electron energy spectrum as chemical species. This allows an investigation of two phenomena in the nighttime mesosphere: the origin of the D-region ledge and the production of radiative emissions from vibrationally excited molecules. It is found that there is negligible difference in the electron densities calculated using the ambient temperature or detailed interaction models, so this study does not provide an explanation for the D-region ledge. However, in the latter model, emissions at various wavelengths are predicted due to reactions involving vibrationally excited molecules. It is also found, using the time-step calculation, that it would take several hours for the predicted electron density to approach equilibrium.

Graphical Abstract

1. Introduction

Electron densities deduced from reflection of radio waves at night showed a very sharp decline with altitude at heights 83 and 89 km [1], shown in Figure 1 for summer at sunspot maximum (“Deeks a”) and winter at sunspot minimum (“Deeks c”). Direct measurements of the electron densities using rockets [2] showed similar declines (although not as sharp) at about 87 km (“Smith (2)” and “Smith (3)” in Figure 1). The sharp decline was described as a ledge by Thomas and Harrison [3]. Very sharp ledges were observed in some rocket measurements in the 1990s, with examples reproduced in Figure 1 as “Friedrich 01” [4], “Goldberg” [5], “Friedrich 99” and “Fnight 99” [6].
Nearly simultaneous measurements of O-atom and electron densities [7] showed similarities, including ledges at about 85 km. A simultaneous measurement of electron density and O-atoms [8] showed a ledge in each at about 86 km, with the comment “The pronounced ledge in the electron densities can be reproduced in theoretical models by a similar ledge in atomic oxygen, but sufficient concurrent measurements of the two parameters are not available to prove such an assertion”.
In the Earth’s nighttime mesosphere ionization is produced by cosmic rays and by Lyman- α  radiation [9] scattered from the Earth’s geocorona [10]. Sample calculated values of the nighttime electron density from various theories are shown in Figure 2, with the “Fnight 99” measurements included for comparison.
In 1968, Radicella [11] derived a model of the mid-latitude nighttime mesosphere in which atmospheric constituents and reactions were adjusted to reproduce the ledge that was apparent in the measurements by Deeks [1]. Free electrons were assumed to be produced in ionization produced by cosmic rays and precipitating electrons, while Lyman- α  was dismissed as an insignificant source. Electrons produced by galactic cosmic rays (GCRs) were presented as equivalent to “ion-pair” production. It was assumed that the free electrons are rapidly attached to molecules to form negative ions, primarily by
O 2 + O 2 + e O 2 + O 2
and are released by
O + O 2 O 3 + e .
Figure 2. Measured [4] (thick line) and calculated nighttime electron densities as a function of altitude by Friedrich and Torkar [4] (”FIRI”– – –), Baumann et al. [12] (— - —), Thomas and Bowman [13] (- - -) and Radicella [11] (as labelled).
Figure 2. Measured [4] (thick line) and calculated nighttime electron densities as a function of altitude by Friedrich and Torkar [4] (”FIRI”– – –), Baumann et al. [12] (— - —), Thomas and Bowman [13] (- - -) and Radicella [11] (as labelled).
Atmosphere 14 00611 g002
The O 2  ions start a complex chain of negative-ion reactions, with negative charge finally being removed by ion-ion and electron-ion dissociative recombination. Very sharp electron-density ledges were deduced (“Radicella 5h” and “Radicella5b” in Figure 2), but at heights that were too low. (The plotted curves were nonphysical, having two different densities at the same height.) Less sharp ledges closer to the height in the experimental data (“Radicella 11h” and “Radicella 11b” in Figure 2) were produced by including the reaction
N 2 + O 3 NO 2 + NO .
After Strobel et al. [10] determined that Lyman- α  and Lyman- β  radiation were sufficient to produce observed electron densities in the mesosphere, Thomas et al. [13] assumed free electrons were produced by this and by the same GCR model as used by Radicella [11]. Their model included many more reactions (than Radicella’s) and produced a sharp decline in electron density downward from 79 km (curve “Thomas” in Figure 2). A semiempirical model “FIRI” [4] did not predict the ledge.
Using the Sodankylä Ion and Neutral Chemistry (SIC) model [14], which combines the ion chemistry of the D region with the chemistry of neutral species [15], Baumann et al. [12] added meteor smoke particles, giving the “Baumann” curve in Figure 2. This model gives good agreement with the measured electron densities in the altitude range 80–100 km.
It can be seen in Figure 2 that in all models, the very sharp density gradient in the measurements is not predicted. The slope of the ledge in the recent “Baumann” model is the same as in the much earlier “Thomas” model.
The approach in this paper is to investigate whether the sharp ledge can be predicted by considering the detailed energy-dependent processes of electrons that are produced by Lyman- α  radiation. In ionospheric calculations, the rate constants for electron interactions are usually defined as a function of electron temperature  T e , based on the assumption that the electron distribution is a Maxwell–Boltzmann distribution. In this work, the electron distribution is calculated as a function of the production and loss mechanisms for electrons, in conjunction with the photochemical model of the chemical reactions. This calculation is facilitated by the method defined and verified by Campbell et al. [16], in which subranges of the electron energy spectrum are treated as chemical species, allowing the electron interactions to be entered as a series of extra rate equations into existing chemical models. As elastic collisions between electrons normally dominate the emergence of a Maxwell–Boltzmann distribution, this method is only applicable where the electron density is very low compared to the molecular density, so that electron–molecule interactions determine the electron energies.
The method is to run simulations of a minimal model, either with electrons assumed to be at the ambient (neutral) temperature, or with electrons initially having a higher temperature. In the latter case the electrons lose energy in small steps in elastic collisions with molecules and larger steps in inelastic collisions producing vibrational excitation of molecules. Thus, the processes of attachment and recombination, which have rate constants that depend on the electron temperature, may proceed at different rates to those for electrons at the neutral temperature, thus giving a different predicted electron density. In all cases, the approach is to perform a semi-implicit time-step simulation [17] of the populations of the chemical species as a function of time.

2. Materials and Methods

Simulations of atmospheric chemical processes usually involve the processing of a set of reactions specified in the form:
A + B k C + D + ,
where A, B, C and D represent atoms or molecules and k is the rate constant. The rate r at which A and B interact is given by:
r   =   k [ A ] [ B ] ,
where  [ x ]  represents the density of species x and the rate constant k has the units l 3 s 1  where “l” is the unit of length in which the density is specified (e.g., cm 3 s 1 ).
A minimal model involving just Lyman- α  ionization, recombination with NO + , attachment to O 2 , charge exchange involving O 3  and associative electronic detachment is implemented by the 7 reactions (Reaction (R2) [18]; Reactions (R3)–(R8) [13]):
NO + γ k Ly NO + + e k Ly = 1.24 × 10 18 cm 3 s 1
NO + + e k Rc N + O k Rc = 6.0 × 10 7 cm 3 s 1
O 2 + O 2 + e k aO 2 O 2 + O 2 k aO 2 = 4.0 × 10 30 e 193 / T n cm 6 s 1
O 2 + N 2 + e k aN 2 O 2 + N 2 k aN 2 = 1.0 × 10 31 cm 6 s 1
O 3 + O 2 k ce O 3 + O 2 k ce = 6.0 × 10 10 cm 3 s 1
O 3 + O k ed 1 e + O 2 + O 2 k ed 1 = 1.0 × 10 10 cm 3 s 1
O + O 2 k ed 2 O 3 + e k ed 2 = 1.5 × 10 10 cm 3 s 1
The rate constants in (R3)–(R5) are for a Maxwellian distribution of electrons with temperature  T e , which for the lower mesosphere is assumed to be equal to the neutral temperature  T n . Hereafter, this will be referred to as the “nT” model.
As an example, consider only Reactions (R2) and (R3) for altitude 85 km where [NO] =  2.02 × 10 6  cm 3  and the Lyman- α  flux  F Ly  is as specified by Thomas and Bowman [13]. (Reactions (R4)–(R8) are only significant at lower altitudes.) An equilibrium calculation, where the gain is set equal to the loss rate, gives the electron density as:
[ e ]   =   F Ly NO k Ly k Rc   =   79 cm 3 .
This is consistent with the values of ∼100 cm 3  of Thomas and Bowman [13], allowing for an extra 25% from cosmic rays in their model.
As the electron density is low, the recombination rate will be low and so it will take time for the electron density to reach an equilibrium value in response to a change in input parameters, particularly the Lyman- α  ionization rate. Thus, a nonequilibrium calculation is required. A way to simulate the evolution of a set of interacting species is to apply Equation (1) to each Reaction (R1) for a time-step  Δ t , so that the change in each species is:
Δ [ A ]   =   Δ [ B ]   =   Δ [ C ]   =   Δ [ D ]   =   k [ A ] [ B ] Δ t .
For each species i, the gains and losses are added up to give the total gain  G i  and total loss  L i  so the new density  n i  of species i after time  Δ t  is:
n i ( t + Δ t )   =   n i ( t ) + G i L i .
The development of the densities of all species can be simulated by iterative application of Equation (4) over the required time, but the size of  Δ t  is limited by the requirement that the density must not go negative and so this “explicit” formula is unsuitable for simulation of systems over long time intervals. This problem is solved in an alternative “semi-implicit” method [17], as detailed by Campbell et al. [16], that transforms Equation (4) to:
n i ( t + Δ t )   =   n i ( t ) + G i 1 + L i n i ( t ) .
To allow the simulation to run to equilibrium with minimum computation time, an adaptive time step  Δ t  is required to implement Equation (5) efficiently. The initial value of  Δ t  is set very small ( 10 8  s) and is then successively increased as:
Δ t j + 1   =   f Δ t min   =   f min n i ( t ) | L i G i | for n i ( t ) > 0 Δ t j ,
where  Δ t min  is the minimum time interval for the density of any of the constituents i to fall to zero in the next time step  Δ t j + 1  at the current rate of change and  f   =   0.0001  is a factor that acts to reduce the cumulative error in the calculations.
To incorporate energy-dependent electron interactions into a time-step simulation, the method [16] is to split the electron energy range  [ E min , E max ]  into N subranges R 1 –R N , with R 0  set as the range  [ 0 , E min ] . All electrons within subrange R n  are assigned the energy E n , where this is the energy at the midpoint of the subrange.
In place of a single reaction  A + e k A * + e * , where A *  represents an excited or ionized state of A and  e *  a change in energy of the electron, a series of reactions:
A + R n k n m A * + R m
is entered into the list of chemical reactions, with the number density of electrons in each energy range R n  treated in the same way as the density of a chemical species. (For example, in the current implementation the number of electrons in the range [0,  E min ] is stored in a variable R0, analogous to the number of N 2  molecules being stored in a variable N2.)
A model (designated “EI”) that includes the detailed electron interactions has the reactions above modified to:
NO + γ k Ly NO + + R n Ly
NO + + R n k Rc n N + O
O 2 + O 2 + R n k O 2 n O 2 + O 2
O 2 + N 2 + R n k N 2 n O 2 + N 2
O 3 + O 2 k ce O 3 + O 2
O 3 + O k ed 1 R 0 + O 2 + O 2
O + O 2 k ed 2 O 3 + R 0
and also includes vibrational excitation of the first level of N 2  and O 2 :
N 2 ( ν = 0 ) + R n k vibN 2 n m N 2 ( ν = 1 ) + R m
O 2 ( ν = 0 ) + R n k vibO 2 n m O 2 ( ν = 1 ) + R m
and elastic collisions of electrons with N 2  and O 2 :
N 2 + R n k elasN 2 n m N 2 + R m ( m = n 1 )
O 2 + R n k elasO 2 n m O 2 + R m ( m = n 1 ) .
In the following subsections, the rate constants for recombination ( k Rc n ), attachment ( k O 2 n  and  k N 2 n ), vibrational excitation ( k vibN 2 n m  and  k vibO 2 n m ) and elastic collsions ( k elasN 2 n m  and  k elasO 2 n m ) are deduced. Due to lack of information about the energy of the electrons released in Reactions (R15) and (R16), it is assumed here that they are at  T n  and so are placed in subrange R 0 . Lyman- α  has an energy of 10.21 eV and NO has an ionization energy of 9.2642 eV, so the electrons produced are assumed to have an energy of 0.944 eV. By setting  E max  to 1.0 eV, R n Ly =R N .

2.1. NO+ Recombination

Cross-sections as a function of energy for Reaction (R11) were measured by Walls and Dunn [19]. From these, they calculated rate constants for Maxwellian electron distributions as a function of electron temperature  T e , shown by the solid line in the lower panel of Figure 3.
The rate constant  k Rc  is related to the cross-section  σ  by:
k Rc ( E ) = v σ ( E ) = 2 E m e σ ( E )
where v is the speed of the electron,  m e  is the electron mass and  σ ( E )  is the cross section for electrons of energy E. Applying this to the measured cross sections (digitised from a graph of Walls and Dunn) gives the curve in the upper panel of Figure 3.
To calculate the rate constant for a Maxwellian distribution, an integration is made of the product of  k Rc ( E ) E  and the Electron Energy Distribution Function  F ( E )  at each energy E. These calculations are repeated here to verify the current implementation, giving the dashed curve in Figure 3. These values are slightly less than those of Walls and Dunn, but this discrepancy could be due to inaccuracy in the original publication and/or digitisation errors.
To demonstrate the application of the set of Reaction (R9) and investigate the accuracy of the approximation, the electron energy range  0.0001 10  eV is divided into 11 subranges R 0 –R 10  and the rate constants are averaged in each subrange (as shown by the horizontal segments in the upper panel of Figure 3). The rate constants for Maxwellian distributions at each temperature  T e  are calculated by applying the 11 average rate constants to the Maxwellian electron density functions at the centre of each subrange:
k Rc ( T e ) = n = 0 10 k Rc n F ( E n ) ( E n , n + 1 E n 1 , n )
where  E n , n + 1  is the energy at the boundary of subranges R n  and R n + 1 . Plotting the rate constants calculated using Equation (8) in the lower panel of Figure 3 shows that there is little difference between the precise values and the approximation, showing that the representation of the electron spectrum by just 11 average values is a viable approach.
For the application here  [ E min , E max ]   is   set   to   [ 0.02 , 1.0 ] , so that the 0.944-eV electrons produced by Lyman- α  ionization are in the highest-energy subrange and the electron energy of 0.019 eV, corresponding to the neutral temperature of 220 K, is in the lowest subrange R 0 . As shown in Figure 4 R 0  is defined as the range [0.002–0.02], for which the averaged value of the rate constant is  6.0 × 10 7 . This is the value used by Thomas and Bowman [13] in the study that is to be partially emulated here. In this case, 21 energy subranges are used to be consistent with the resolution needed for vibrational excitation (to be described below).

2.2. Attachment

Measured rate constants for attachment of electrons for O 2  as a function of electron energy [20] are reproduced in Figure 5 for Reaction (R12) at neutral temperatures of 77 K and 300 K and for Reaction (R13) at 300 K. For comparison, the rate constants for these reactions for a Maxwellian distribution of electrons as a function of the electron temperature are shown. A vertical line corresponding to the neutral temperature of 220 K shows that the rate constants for electrons at ∼0.1 eV are substantially higher than that used in the nT model (labelled “O 2  Thomas”) for Reaction (R12), i.e., about  2 × 10 30  cm 6 s 1  at T n   =   220  K.
To obtain the rate constants for Reaction (R12) at any neutral temperature, an interpolation is performed between the values at 77 K and 300 K. The shape of the curve at 300 K is copied to extend the measured values for Reaction (R13) down through lower energies, then the entire curve is mapped to another neutral temperature by comparison with the interpolation for Reaction (R12). The results of these operations are plotted in Figure 5 for a neutral temperature of 220 K. (i.e., “O 2  220 K fit” and “N 2  220 K fit”).
To obtain the rate coefficients  k O 2 n  and  k N 2 n , the functions plotted are then averaged over subranges in the same way as for the recombination rate constants.

2.3. Vibrational Excitation and Elastic Collisions

The 0.944-eV electrons produced in Reaction (R10) can excite the first vibrational of N 2 , losing 0.289 eV in the process. The inelastic cross sections  σ 01 ( E )  for this excitation (the “Flinders” set of Campbell et al. [21]) are shown as a function of electron-impact energy in the lower panel of Figure 6. Representative rate constants  k vibN 2 ( E ) , calculated using Equation (7), are shown by lines drawn from the initial energy of the electron to the final energy, at the level of the rate constant on the vertical axis in the upper panel.
The averaged rate constants for all transitions from one electron-energy subrange to another are shown by arrows (with solid lines) drawn from the centre of the initial subrange to the centre of the final subrange. These arrows represent the set of Reaction (R17). (The bullets show minor adjustments, where the electron remains in the same subrange, to implement energy conservation for logarithmically-spaced energy intervals [16].) For vibrational excitation of O 2 , the cross-sections published by Jones et al. [22] were used to derive the averaged rate constants for Reaction (R18).
The average energy transferred in an elastic collision between an electron with energy E and a molecule of mass M is given by [23]
Δ E   =   2 m e M σ MT ( E ) E v
where  σ MT  is the elastic momentum transfer cross-section. These cross-sections were obtained from Table 6.3.5.3 for e–N 2  collisions and Table 6.3.5.6 for e–O 2  collisions in a compilation by Elford et al. [24].
The rate constant for elastic transitions of electrons from subrange m to  m 1  is then that required to transfer energy  E m E m 1  i.e.,:
k elas m , m 1   =   Δ E E m E m 1 .
These averaged rate constants are shown by the arrows with dotted lines in Figure 6.

2.4. Emissions Arising from Vibrational Excitation

Emissions can arise from vibrationally excited O 2  and N 2 , by Reactions (R21)–(R23) [25] and (R24) [26]:
N 2 ( 1 ) + O k r 1 N 2 ( 0 ) + O + 2331 cm 1 k r 1 = 2.45 × 10 22 T 2.87
O 2 ( 1 ) + N 2 k r 2 O 2 ( 0 ) + N 2 + 1556 cm 1 k r 2 = 1.364 × 10 12 T e 137.93 / T 1 / 3
O 2 ( 1 ) + O 2 k r 3 O 2 + O 2 + 1556 cm 1 k r 3 = 1.364 × 10 12 T e 137.93 / T 1 / 3
N 2 ( 1 ) + CO 2 ( 0 ) k r 4 N 2 ( 0 ) + CO 2 ( ν 3 ) k r 4 = 8.91 × 10 12 T 1
In the last case, the excited CO 2  can produce 4.3 μm radiation [26] and also emissions at 2.7 μm, 2.0 μm, etc. [25].

3. Results

In Figure 7, the “nT” model (Reactions (R2)–(R8)) is applied with rate constants based on a Maxwellian distribution of electrons at a neutral temperature  T n   = 220 K. The electron density as a function of altitude (shown by the solid red line, labelled “e  nT”) is calculated by time-step simulation applied for midnight conditions for a time of 10 h, starting with zero electron density.
The simulation is then repeated for the electron interaction model without elastic scattering i.e., Reactions (R10)–(R18), with the densities after 10 h shown by symbols. At 77 km and below the electron density is less than in the “nT” model, and substantially less at 72 km and below. Most of the electrons are in subrange R 0 , but some are in subranges R 1 –R 13 , with the proportion increasing with altitude. Densities of vibrationally excited molecules N 2 (1) and O 2 (1) are higher than the electron densities.
The “EI” model is repeated in Figure 8 with elastic scattering (Reactions (R19) and (R20)) and emission-producing Reactions (R21)–(R24) included. There is now no apparent difference in the electron densities for the “nT” and “EI” cases. Here (for the same density range as before), there are electrons in subranges R 1  and R 2  only at the highest altitudes, with all in R 0  at lower altitudes. The total production of CO 2 ( ν 3 )  (labelled “CO 2 ( 3 ) ”) is plotted as a density, as are the number of photons/cc emitted in 10 h for 2331 cm 1  and 1556 cm 1 .
In Figure 9, the densities of both background species and those produced by the EI model are shown as a function of time, at an altitude of 80 km. The electron density reaches the equilibrium value of ∼80 cm 3  after about 6 h.

4. Discussion

In Figure 7, there is little difference between the electron densities calculated by the nT and EI models, except for a reduction at altitudes 71–76 km. The reduction is as expected due to the higher rate of the Reaction (R12) at higher electron temperatures. However, as the reduction is small, it does not provide an explanation for the D-region ledge. This is more evident in Figure 8 when the electrons also lose energy in elastic collisions with N 2  and O 2 , with the result that there is no discernible difference in the predicted electron densities. In the nT model, the electrons do not have sufficient energy to produce vibrational excitation of molecules, but, as seen in Figure 8, there is some production of vibrationally excited N 2 (1) and O 2 (1), leading to radiative emissions at 1556 cm 1  and 2331 cm 1  and from CO 2 ( ν 3 ) . However, these emissions may be very small compared to what is expected due to thermal production. For example, if the gas is at thermal equilibrium, at 80 km [N 2 ] =  3.8 × 10 14  cm 3 , so [N 2 (1)] is expected to be about  9 × 10 7  cm 3 , which is much larger than the predicted density of ∼0.3 cm 3 . However, this depends on the rate of translational–vibrational (TV) transitions, which would need to be included in the model to make an accurate prediction of the thermal emissions.
In Figure 9, it is shown that it takes about 6 h for the electron density to reach the equilibrium value. While starting from zero electron density is not physically realistic, it implies that there will be a substantial time lag in the reduction in electron density from daytime values, compared to instantaneous equilibrium. Thus time-dependent modelling is required for calculations of the nighttime mesosphere. This suggests that it could be worthwhile to implement the above calculation (which applied midnight parameters over 10 h) to a fully time-dependent simulation, with the Lyman- α  flux varied according to time of day or night. As this flux and also the rates of chemical reactions vary with altitude, it is likely that the rate of progress towards equilibrium will also vary with height, producing a different altitude profile for the electron density relative to that produced by assuming equilibrium is attained.
The loss of energy by higher-energy electrons due to elastic collisions with thermal electrons is not included in this study. While it is not expected to be significant, given the very high ratio of molecules to electrons, it does need to be quantified so that the method can be applied with confidence at higher electron densities.
The EI approach could be useful in predicting emissions produced by higher-energy electrons (>5 eV) that could not be produced by thermal excitation. This depends on the secondary electron spectrum produced by cosmic rays—in particular whether cosmic rays produce ions and electrons with high energy, or ion-pairs from which electrons are then released with low energy. Thus, the EI model could be used to predict emissions that could be used as diagnostics of the secondary electron spectrum produced by cosmic rays.

5. Conclusions

A time-step simulation that includes energy-dependent electron interactions has been applied to simulate electron interactions in the Earth’s quiet nighttime mesosphere. A minimal model with just a few processes (Lyman- α  production of electrons and NO + , attachment of electrons to O 2 , charge transfer to O 3 , associative electron detachment and dissociative recombination of electrons with NO + ) was run for both the assumption that the electrons are at the neutral temperature and that they have higher energy for some finite time as they progressively lose it in elastic and inelastic collisions with molecules. It was found that there is little difference in the electron densities predicted by the two approaches, so consideration of energy-dependent reaction rates does not provide an explanation for the D-region ledge. It was found that the presence of higher-energy electrons leads to the production of vibrationally excited N 2  and O 2  with subsequent radiative emissions, but that the densities are small compared to those expected from thermal excitation. However, the EI model could be used to predict emissions produced by cosmic rays, if these produce higher-energy electrons. It was also found that it takes of the order of 6 h for the calculated electron densities to reach equilibrium at 80 km, showing that only time-dependent models are adequate for calculations of the lower mesosphere.

Author Contributions

Conceptualization, L.C. and M.J.B.; methodology, L.C. and M.J.B.; software, L.C.; validation, L.C.; formal analysis, L.C.; investigation, L.C.; resources, M.J.B.; data curation, L.C.; writing—original draft preparation, L.C. and M.J.B.; writing—review and editing, L.C.; visualization, L.C; supervision, M.J.B.; project administration, M.J.B.; funding acquisition, M.J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Australian Research Council (Projects No. DP190100696 and No. DP180101655).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used in this study was obtained from the literature as quoted or is available from the LXCAT website.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Deeks, D.G. D-region electron distributions in middle latitudes deduced from the reflexion of long radio waves. Proc. R. Soc. Lond. A 1966, 291, 413–437. [Google Scholar]
  2. Smith, L.G. A sequence of rocket observations of night-time sporadic-E. J. Atmos. Terr. Phys. 1970, 32, 1247–1257. [Google Scholar] [CrossRef]
  3. Thomas, L.; Harrison, M.D. The electron density distributions in the D-region during the night and pre-sunrise period. J. Atmos. Terr. Phys. 1970, 32, 1–14. [Google Scholar] [CrossRef]
  4. Friedrich, M.; Torkar, K.M. FIRI: A semiempirical model of the lower ionosphere. J. Geophys. Res. 2001, 106, 21409–21418. [Google Scholar] [CrossRef]
  5. Goldberg, R.A.; Lehmacher, G.A.; Schmidlin, F.J.; Fritts, D.C.; Mitchell, J.D.; Croskey, C.L.; Friedrich, M.; Swartz, W.E. Equatorial dynamics observed by rocket, radar, and satellite during the CADRE/MALTED campaign 1. Programmatics and small-scale fluctuations. J. Geophys. Res. 1997, 102, 26179–26190. [Google Scholar] [CrossRef] [Green Version]
  6. Friedrich, M.; Gumbel, J.; Pilgram, R. Atomic oxygen in the mesosphere and its relevance for the ionosphere: A summary of empirical evidence. In Proceedings of the 14th ESA Symposium on European Rocket and Balloon Programmes and Related Research, Potsdam, Germany, 31 May–3 June 1999; (ESA SP-437, September 1999). pp. 287–290. [Google Scholar]
  7. Dickinson, P.H.G.; Bain, W.C.; Thomas, L.; Williams, E.R.; Jenkins, D.B.; Twiddy, N.D. The Determination of the Atomic Oxygen Concentration and Associated Parameters in the Lower Ionosphere. Proc. R. Soc. Lond. A 1980, 369, 379–480. [Google Scholar]
  8. Friedrich, M. The lower ionosphere: Abandoned by communication, to be re-discovered for aeronomy. Radio Sci. Bull. 2004, 38–46. [Google Scholar] [CrossRef]
  9. Mironova, I.A.; Aplin, K.L.; Arnold, F.; Bazilevskaya, G.A.; Harrison, R.G.; Krivolutsky, A.A.; Nicoll, K.A.; Rozanov, E.V.; Turunen, E.; Usoskin, I.G. Energetic Particle Influence on the Earth’s Atmosphere. Space Sci. Rev. 2015, 194, 1. [Google Scholar] [CrossRef] [Green Version]
  10. Strobel, D.F.; Young, T.R.; Meier, R.R.; Coffey, T.P.; Ali, A.W. The Nighttime Ionosphere: E Region and Lower F Region. J. Geophys. Res. 1974, 79, 3171–3178. [Google Scholar] [CrossRef]
  11. Radicella, S.M. Theoretical models of electron and ion density in the night-time D-region. J. Atmos. Terr. Phys. 1968, 30, 1745–1760. [Google Scholar] [CrossRef]
  12. Baumann, C.; Rapp, M.; Anttila, M.; Kero, A.; Verronen, P.T. Effects of meteoric smoke particles on the D region ion chemistry. J. Geophys. Res. Space Phys. 2015, 120, 10823–10839. [Google Scholar] [CrossRef] [Green Version]
  13. Thomas, L.; Bowman, M.R. Model studies of the D-region negative-ion composition during day-time and night-time. J. Atmos. Terr. Phys. 1985, 47, 547–556. [Google Scholar] [CrossRef]
  14. Verronen, P.T.; Seppälä, A.; Clilverd, M.A.; Rodger, C.J.; Kyrölä, E.; Enell, C.F.; Ulich, T.; Turunen, E. Diurnal variatin of ozone depletion during the October-November 2003 solar proton events. J. Geophys. Res. 2005, 110, e2004JA010932. [Google Scholar] [CrossRef] [Green Version]
  15. Verronen, P.T.; Santee, M.L.; Manney, G.L.; Lehmann, R.; Salmi, S.M.; Seppälä, A. Nitric acid enhancements in the mesosphere during the January 2005 and December 2006 solar proton events. J. Geophys. Res. 2011, 116, e2011JD016075. [Google Scholar] [CrossRef] [Green Version]
  16. Campbell, L.; Muccignat, D.L.; Brunger, M.J. Inclusion of Electron Interactions by Rate Equations in Chemical Models. Atoms 2022, 10, 62. [Google Scholar] [CrossRef]
  17. Cariolle, D.; Moinat, P.; Teyssèdre, H.; Giraud, L.; Josse, B.; Lefèvre, F. ASIS v1.0: An adaptive solver for the simulation of atmospheric chemistry. Geosci. Model Dev. 2017, 10, 1467–1485. [Google Scholar] [CrossRef] [Green Version]
  18. Watanabe, K.; Marmo, F.F.; Inn, E.C.Y. Photoionization cross section of nitric oxide. Phys. Rev. 1953, 91, 1155–1158. [Google Scholar] [CrossRef]
  19. Walls, F.L.; Dunn, G.H. Measurement of total cross sections for electron recombination with NO+ and O 2 + using ion storage techniques. J. Geophys. Res. 1974, 79, 1911–1915. [Google Scholar] [CrossRef]
  20. Chanin, L.M.; Phelps, A.V.; Biondi, M.A. Measurements of the attachment of low-energy electrons to oxygen molecules. Phys. Rev. 1962, 128, 219–230. [Google Scholar] [CrossRef]
  21. Campbell, L.; Brunger, M.J.; Cartwright, D.C.; Teubner, P.J.O. Production of vibrationally excited N2 by electron impact. Planet. Space Sci. 2004, 52, 815–822. [Google Scholar] [CrossRef]
  22. Jones, D.B.; Campbell, L.; Bottema, M.J.; Teubner, P.J.O.; Cartwright, D.C.; Newell, W.R.; Brunger, M.J. Electron-driven excitation of O2 under night-time auroral conditions: Excited state densities and band emissions. Planet. Space Sci. 2006, 54, 45–59. [Google Scholar] [CrossRef]
  23. Winstead, C.; McKoy, V. Electron-molecule collisions in low-temperature plasmas: The role of theory. In Advances in Atomic, Molecular and Optical Physics; Bederson, B., Walther, H., Eds.; Academic Press: San Diego, CA, USA, 2000; Volume 43. [Google Scholar]
  24. Elford, M.T.; Buckman, S.J.; Brunger, M.J. Elastic momentum transfer cross sections. In Photon and Electron Interactions with Atoms, Molecules and Ions; Itikawa, Y., Ed.; Landolt-Börnstein; Springer: Berlin/Heidelberg, Germany, 2003; Volume 17/C, Chapter 6.3. [Google Scholar]
  25. Shved, G.M.; Kutepov, A.A.; Ogivalov, V.P. Non-local thermodynamic equilibrium in CO2 in the middle atmosphere. I. Input data and populations of the ν3 mode manifold states. J. Atmos. Solar-Terr. Phys. 1998, 60, 289–314. [Google Scholar] [CrossRef]
  26. Panka, P.A.; Kutepov, A.A.; Kalogerakis, K.S.; Janches, D.; Russell, J.M.; Rezac, L.; Feofilov, A.G.; Mlynczak, M.G.; Yiğit, E. Resolving the mesospheric nighttime 4.3 μm emission puzzle: Comparison of the CO2(ν3) and OH(ν) emission models. Atmos. Chem. Phys. 2017, 17, 9751–9760. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Nighttime electron densities as a function of altitude, measured by radio reflection (Deeks a,c) [1] and by rockets: (Friedrich 01) [4], (Goldberg) [5], (Smith (2,3)) [2] and (Fnight 99, Friedrich 99)) [6].
Figure 1. Nighttime electron densities as a function of altitude, measured by radio reflection (Deeks a,c) [1] and by rockets: (Friedrich 01) [4], (Goldberg) [5], (Smith (2,3)) [2] and (Fnight 99, Friedrich 99)) [6].
Atmosphere 14 00611 g001
Figure 3. (upper panel) Rate constants (cm 3 s 1 ) for the recombination NO +  +  e  as a function of electron energy (eV), as calculated from the cross sections of Walls and Dunn [19] (curve) and averaged over energy subranges (horizontal segments). (lower panel) Rate constants for a Maxwellian distribution of electrons as a function of the electron temperature, as given by Walls and Dunn [19] (—–), a current emulation of their calculations (- - - -) and an approximate calculation using the 11 averaged rate constants (– – –).
Figure 3. (upper panel) Rate constants (cm 3 s 1 ) for the recombination NO +  +  e  as a function of electron energy (eV), as calculated from the cross sections of Walls and Dunn [19] (curve) and averaged over energy subranges (horizontal segments). (lower panel) Rate constants for a Maxwellian distribution of electrons as a function of the electron temperature, as given by Walls and Dunn [19] (—–), a current emulation of their calculations (- - - -) and an approximate calculation using the 11 averaged rate constants (– – –).
Atmosphere 14 00611 g003
Figure 4. Rate constants (cm 3 s 1 ) for the recombination NO +  +  e  as a function of electron energy (eV), as calculated from the cross sections of Walls and Dunn [19] (curve), with this function averaged over energy subranges (horizontal segments) for the parameters of the EI calculation.
Figure 4. Rate constants (cm 3 s 1 ) for the recombination NO +  +  e  as a function of electron energy (eV), as calculated from the cross sections of Walls and Dunn [19] (curve), with this function averaged over energy subranges (horizontal segments) for the parameters of the EI calculation.
Atmosphere 14 00611 g004
Figure 5. Rate constants (cm 6 s 1 ) as a function of electron energy (eV) for attachment Reaction (R12) (labelled with “O 2 ”) measured [20] at 77 K (――) and 300 K (– – – –), calculated for 220 K (— —) and as given [13] for a Maxwelllian distribution (— - — - —). Similarly, for (R13), (labelled with “N 2 ”) measured [20] at 300 K (— — —), calculated for 220 K (- - -) and for a Maxwellian distribution [13] (―― – ――).
Figure 5. Rate constants (cm 6 s 1 ) as a function of electron energy (eV) for attachment Reaction (R12) (labelled with “O 2 ”) measured [20] at 77 K (――) and 300 K (– – – –), calculated for 220 K (— —) and as given [13] for a Maxwelllian distribution (— - — - —). Similarly, for (R13), (labelled with “N 2 ”) measured [20] at 300 K (— — —), calculated for 220 K (- - -) and for a Maxwellian distribution [13] (―― – ――).
Atmosphere 14 00611 g005
Figure 6. (Lower panel) Cross-section (cm 2 ) as a function of electron energy (eV) for electron-impact excitation of the first vibrational level in N 2 . (upper panel) Rate constants (cm 3 s 1 ) as a function of electron energy (eV) derived from the cross-sections, represented by horizontal lines drawn from the initial energy of the electron to its final energy. Averaged rate constants for all transitions from one energy subrange to another, for vibrational excitation (arrows with solid lines) and elastic collisions (arrows with dotted lines).
Figure 6. (Lower panel) Cross-section (cm 2 ) as a function of electron energy (eV) for electron-impact excitation of the first vibrational level in N 2 . (upper panel) Rate constants (cm 3 s 1 ) as a function of electron energy (eV) derived from the cross-sections, represented by horizontal lines drawn from the initial energy of the electron to its final energy. Averaged rate constants for all transitions from one energy subrange to another, for vibrational excitation (arrows with solid lines) and elastic collisions (arrows with dotted lines).
Atmosphere 14 00611 g006
Figure 7. Calculated densities as a function of altitude: electron (—) and NO +  (- - -) densities for the nT method; densities for electrons (⋄), NO +  (□), molecules and electron subranges as labelled for the EI method without elastic collisions included.
Figure 7. Calculated densities as a function of altitude: electron (—) and NO +  (- - -) densities for the nT method; densities for electrons (⋄), NO +  (□), molecules and electron subranges as labelled for the EI method without elastic collisions included.
Atmosphere 14 00611 g007
Figure 8. Calculated densities as a function of altitude: electron (—) and NO +  (- - -) densities for the nT method; densities for electrons (⋄), NO +  (□), molecules, electron subranges and emitted photons as labelled for the EI method with elastic scattering included.
Figure 8. Calculated densities as a function of altitude: electron (—) and NO +  (- - -) densities for the nT method; densities for electrons (⋄), NO +  (□), molecules, electron subranges and emitted photons as labelled for the EI method with elastic scattering included.
Atmosphere 14 00611 g008
Figure 9. Densities of atmospheric constituents and cumulative emitted photons (as labelled) as a function of time (h) for the EI simulation.
Figure 9. Densities of atmospheric constituents and cumulative emitted photons (as labelled) as a function of time (h) for the EI simulation.
Atmosphere 14 00611 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Campbell, L.; Brunger, M.J. Modelling of Energy-Dependent Electron Interactions in the Earth’s Mesosphere. Atmosphere 2023, 14, 611. https://doi.org/10.3390/atmos14040611

AMA Style

Campbell L, Brunger MJ. Modelling of Energy-Dependent Electron Interactions in the Earth’s Mesosphere. Atmosphere. 2023; 14(4):611. https://doi.org/10.3390/atmos14040611

Chicago/Turabian Style

Campbell, Laurence, and Michael J. Brunger. 2023. "Modelling of Energy-Dependent Electron Interactions in the Earth’s Mesosphere" Atmosphere 14, no. 4: 611. https://doi.org/10.3390/atmos14040611

APA Style

Campbell, L., & Brunger, M. J. (2023). Modelling of Energy-Dependent Electron Interactions in the Earth’s Mesosphere. Atmosphere, 14(4), 611. https://doi.org/10.3390/atmos14040611

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop