Modeling the behavior of hot oxygen ions

Photochemical processes in the upper atmosphere are known to create significant amounts of energetic oxygen atoms or “hot O”. In this research we simulate the effects of ionized hot oxygen, hot O , on the ionosphere. We find that hot O is not able to maintain a temperature substantially above the ambient ion temperature at most altitudes, the exception being around the F-region ion density peak. However, the thermalization of hot O , due to Coulomb collisions, represents an important heating process for the ambient ions. A time-dependent, fluid-kinetic model of the ionosphere (TRANSCAR) is used to self-consistently simulate hot O by considering it to be a separate species from O . A Maxwellian neutral hot O population having characteristics consistent with current knowledge is added to TRANSCAR. The production of the hot O + is then computed by considering ion charge exchange with the neutral hot O population that we have assumed. Loss of hot O + results from these charge exchange reactions and from reactions with molecular atoms.


Introduction
It has been observed that the ion temperature in the ionosphere significantly exceeds what theoretical calculations say the temperature should be, (Oliver and Glotfelty, 1996;Oliver and Schoendorf, 1999).A neglected heat source could account for this difference.This fact led researchers to propose the existence of a non-thermal population of oxygen atoms, known as hot O.There is considerable in situ evidence for the existence of hot O (Hedin, 1989;Cotton et al., 1993).
Correspondence to: M. Zettergren (mattz@bu.edu)While there are alternative explanations to Hedin's conclusion of a population of hot O (Keating et al., 1997), there is also radar and optical evidence to support the existence of hot O. Oliver (1997) analyzed radar data and concluded that an unknown heat source for the ions, such as hot O, is needed to explain the observed ion temperature.Optical observations made by Hernandez (1971) show a broadened oxygen spectrum, indicating the presence a high-temperature population of oxygen.Most work estimates the effective temperature of hot O to be in the range of 4000 K−6000 K, and the concentration to be about 0.1%−1% of the cold O density at 400 km altitude, (Rohrbaugh and Nisbet, 1973;Shematovich et al., 1994;Oliver, 1997;Litvin and Oliver, 2000).The profile shape of the hot O density has not been determined conclusively; hot O could form a layer shape (Cotton et al., 1993;Schoendorf et al., 2000), it could approximate diffusive equilibrium (Oliver, 1997), or it could form some other type of exponential shape (Yee et al., 1980;Shematovich et al., 1994).Research has shown that hot O can be produced in numerous reactions (Richards et al., 1994;Hickey et al., 1995) and can have a potentially large effect on the energy dynamics of the ionosphere (Oliver, 1997;Alcaydé et al., 2001).Hot O ionizes in the same ways as thermal oxygen: via photoionization and charge exchange chemistry.Airglow observations have confirmed the fact that hot O + can exist at least long enough to undergo radiative decay and emission (Yee et al., 1980).However, the production, loss, and heating effects of hot O + have not yet been investigated in any detail, despite this evidence that supports its existence.
The purpose of our work is to make a first attempt to characterize hot O + behavior and its effects on the upper atmosphere.We want to establish two main things.First, can hot O + exist in appreciable densities?Second, can hot O + maintain an elevated temperature significantly above the ambient ion temperature, or does it cool to the same temperature?We approach this goal in two ways.First we present a simplified analysis of the density and temperature of hot O + .Then Published by Copernicus GmbH on behalf of the European Geosciences Union.
Table 1.TRANSCAR chemistry model for hot O + and the associated reactions rates.k 1 and k 1 are taken from Stancil et al. (1999); k 3 and k 4 are taken from Diloy et al. (1996).We have assumed that k 2 =k 1 since there are no available reaction rates for O + -O charge exchange.

Approximation of hot O + profile shape
We propose that, in the steady-state, hot O + forms a layer dominated by chemistry at lower altitudes and by diffusion at higher altitudes.A layer of this type has maximum density where the chemical loss time constant equals the diffusion time constant (Rishbeth and Garriot, 1969).The chemical reactions that we consider for hot O + are shown in Table 1.
The loss processes for hot O + are the "primed" reactions, R1 -R4 .The chemical loss time constant for these reactions, β −1 , is given by: k i is the reaction constant for reaction R i , and n N i (z) is the number density of the neutral reactant of reaction R i as a function of altitude z.The diffusion time constant is given by Rishbeth and Garriot (1969): . (2) Here D(z) is the diffusion coefficient, H s is the scale height for s=O + hot , k b is Boltzmann's constant, T s is the temperature for species s=O + hot , ν sj is the collision frequency between species s and species j , m s is the mass of species s, and g is the acceleration of gravity.The collision frequency that we are concerned with in this approximation is that of the dominant collision process involving hot O + , collisions with O + .With these considerations the hot O + layer density peak occurs where the chemical and diffusion time constants are equal.
For computation we adopt profiles for neutral densities and collision frequencies from our TRANSCAR control simulation.In this control simulation we use TRANSCAR to simulate the ionosphere without hot O or hot O + included.Using T s =4000 K the peak of the layer occurs at z pk ≈840 km.
It is assumed that there is zero flux divergence at the hot O + density peak, and we consider only chemical sources on hot O + .Under these assumptions the density is given by the ratio of the production rate of hot O + to the loss frequency of hot O + , evaluated at the peak altitude z pk .
k 1 is the reaction constant for reaction R 1 of Table 1, and k 2 is the reaction constant for reaction R 2 .Our simulations of hot O + show that the assumption of zero flux divergence at the peak is accurate and that photoionization can also be ignored.Again using values from the control simulation of TRANSCAR for the O + density and assuming hot O to be in diffusive equilibrium with a reference density of 1% of O density at 400 km, we find the hot O + peak density to be n s (z pk )≈5.62×10 10 m −3 .Below this peak the hot O + density drops off with altitude as a result of increasingly rapid loss chemistry.Far above the hot O + peak, collisions, and thus chemical reactions, become scarce and diffusion dominates.

Approximation of hot O + temperature behavior
To analyze the behavior of the hot O + temperature, we first consider the case where a sample population of hot O + is injected into the ionosphere and allowed to "relax" and thermalize down to the temperature of the ambient ions.
The 1-D equation describing the evolution of the hot O + temperature with time is given by Blelly et al. (1996b): For the case of the hot O + temperature s=O + hot , z is the independent spatial variable (distance along magnetic field line), t is the independent time variable, and A is the cross section area of the magnetic flux tube as a function of z.L s and Q s are energy loss and source rates, respectively.(δT s /δt) * is a source term, explained later, that includes energy input due to production processes.For our simplified calculations we assume that Q s =L s =0 and (δT s /δt) * =0.
We make a few approximations in order to solve Eq. ( 5).All drift velocities and the hot O + heat flow are assumed to be zero so that transport of energy becomes unimportant.This enables us to consider the temperature at a given altitude to be independent of the temperature at other altitudes, and to approximate how the hot O + temperature behaves as a function of time.With these approximations and noting that s=O + hot , the temperature equation becomes: There is an equation of this form for the temperature at each altitude.Since Coulomb collisions dominate the dynamics of the F-region, the summation can be reduced to a sum over the major ions, j 1 =O + at lower altitudes and j 2 =H + at high altitudes: If the collision frequencies and temperatures vary slowly with time as compared to the hot O + temperature, then the solution of this equation is: From the solution for T s it can be seen that T s (∞)=C 1 .For most altitudes the O + and H + temperatures will be almost the same so the approximation T j 1 =T j 2 =T i is appropriate.Inserting these values for C 1 yields: This leads us to predict that hot O + thermalizes to the ambient ion temperature, T i with an approximate time constant of 1/ν , as seen from the solution for T s .Using our control simulation to get representative values for the collision frequencies involved in these calculations, we obtain a value for the decay time constant τ =1/ν as a function of height.A plot of this time constant is shown in Fig. 1 along with a plot of the time constant found from our simulation of hot O + , which is discussed in detail later.As expected the time constant increases with altitude as collisions become increasingly infrequent.
The energy source, (δT s /δt) * , needs to be considered for determining the steady-state equilibrium temperature of hot O + .(δT s /δt) * represents the energy released into the hot O + medium during its production, given by: Here P s represents the total production rate per unit volume of hot O + .Intuitively, (δT s /δt) * is the rate of temperature change due to production of s=O + hot .The main reactions that produce hot O + are charge exchange reactions involving neutral hot O, in which neutral hot O loses an electron to become hot O + .At the instant of creation, a hot O + ion has kinetic energy similar to that of the neutral hot O atom from which it is produced.The fluid description of ions that is used by TRANSCAR does not track individual particle energy, but, rather the average energy and temperature of each ion species population.Any hot O + produced by charge exchange is assumed to have the temperature of the rest of the population at the time of creation.To model the fact that the newly created hot O + has energy comparable to that of its  this value for (δE s /δt) * into equation 14 yields equation 13.The inclusion of a heat source due to chemical creation, (δT s /δt) * , can be applied to ions other than hot O + , as well.However, as seen from equation 13, this source term is significant only when the ion and its parent have a large temperature difference.
If the (δT s /δt) * heat source is included, then the temperature equation can be solved for the steady state equilibrium temperature by setting ∂T s /∂t = 0: T s,eq is now a weighted sum of the ion and hot O temperatures.In regimes where the production of hot O + is large relative to its density, the equilibrium temperature is higher than T i .However, at altitudes where the production is relatively low, this heat source can be ignored and the temperature of hot O + is equal to that of the other ions.Our simulations will answer the question of at what altitudes P s /n s is large and the temperature of hot O + remains elevated, and at what altitudes this heating effect is insignificant and the hot O + thermalizes.It should also be noted that if equation 5 is solved with the (δT s /δt) * term from equation 13 included, then the solution retains the form of equation 9 with

Numerical Modeling
The computer simulation that we use to model the effects of hot O on the ionosphere and the behavior of hot O + is TRANSCAR (described in Blelly et al. (1996b)).TRAN-SCAR is a time-dependent, 8-moment fluid model of the ionosphere that solves for the density n s , drift velocity u s , temperature T s , and heat flow q s for 6 different ion species (s = O + , H + , N + , N + 2 , NO + , and O + 2 ).TRANSCAR computes a numerical solution to the conservation equations for each of these moments in 1 dimension, along the geomagnetic field lines.
TRANSCAR assumes the distribution function of the gases that it models to be adequately described with 8 moments.An 8-moment distribution function can describe gases that are somewhat non-Maxwellian through the heat flow vector, but a gas that is highly non-Maxwellian, such as a mixture of hot and cold oxygen, cannot be adequately accounted for in a 8-moment model.In the case of oxygen, the distribution function has a Maxwellian center that represents thermal, or cold O, with a pronounced high-energy tail that represents hot O.We have chosen to break apart the oxygen distribution function into ambient oxygen and hot O and model them as separate species with different distribution functions.
In this paper, we perform 3 simulations with TRANSCAR.One simulation is a control and does not include hot O or hot O + .The second simulation includes just neutral hot O. parent, hot O atom and not the same energy as the "old" hot O + , the energy source of Eq. ( 13) is needed.
The expression for (δT s /δt) * is derived by considering the energy exchange between neutral hot O and hot O + during production processes involving these species.(δT s /δt) * can be expressed in the following form: (δn s /δt) * , the density change due to production, is P s .By definition: δε/δt is the rate of energy released into the hot O + medium due to production and c s is the random velocity of the gas.
If we assume that all energy is released at the same level 3/2 k b T O hot , then this expression becomes: This expression assumes that a newly created hot O + ion inherits the kinetic energy of its parent atom.Substituting this value for (δE s /δt) * into Eq.( 14) yields Eq. ( 13).The inclusion of a heat source due to chemical creation, (δT s /δt) * , can be applied to ions other than hot O + , as well.However, as seen from Eq. ( 13), this source term is significant only when the ion and its parent have a large temperature difference.
If the (δT s /δt) * heat source is included, then the temperature equation can be solved for the steady state equilibrium temperature by setting ∂T s /∂t=0: T s,eq is now a weighted sum of the ion and hot O temperatures.In regimes where the production of hot O + is large relative to its density, the equilibrium temperature is higher than T i .However, at altitudes where the production is relatively low, this heat source can be ignored and the temperature of hot O + is equal to that of the other ions.Our simulations will answer the question of at what altitudes P s /n s is large and the temperature of hot O + remains elevated, and at what altitudes this heating effect is insignificant and the hot O + thermalizes.It should also be noted that if Eq. ( 5) is solved with the (δT s /δt) * term from Eq. ( 13) included, then the solution retains the form of Eq. ( 9  to the weighted sum of the collision frequencies O + and H + .To see how accurate these calculations are, we have run a simplified hot O + simulation that omits the effects of productions, energy transport, and drag on the hot O + .Using this simulation, we obtain from the time-dependent temperature output a time constant which may be compared with that estimated analytically from equation 10. Figure 1 shows that these two time constants follow the same general trend with altitude.τ ranges from ≈ 0.5 min at lower altitudes to ≈ 2 min at the highest altitudes that we simulate.
Hot O + does not thermalize completely near the F-region peak.This is a result of the heat source of equation 13.At the F-region peak the hot O + density is low due to loss processes involving H, O, O 2 , and N 2 .However, the O + density is at its peak so there is still a large hot O + production rate from reaction R 2 of table 1.This results in a significant heat source for the hot O + population and alters its equilibrium temperature to be greater than that of the ambient ions.
The introduction of hot O + into TRANSCAR has a considerable effect on the rest of the simulated ionosphere.Specifically, hot O + is a heat source for the ambient ions.Figure 3 shows that hot O + has the effect of increasing the ambient ion temperature, in addition to the increase imposed by the neutral hot O.At higher altitudes the presence of hot O + has the effect of raising the ambient ion temperature another ≈ 150 K above the temperature of the simulation with neutral hot O.This is a direct result of the thermalization of hot O + .Hot O + collides with and imparts its energy to the ambient ion population.Consequently, hot O + thermalizes to the ambient ion temperature, but, in the process, causes the ambient temperature to increase.A plot of the heat input to the O + population due to collisions with hot O and hot O + is shown in figure 5.This figure shows that hot O and hot O + are comparably important heat sources for the ambient ions.The simulation with hot O + also shows small but noticeable effects on the densities of O + and H + (due to reactions R 1 , R 1 , R 2 , and R 2 in table 1, and changes in scale height) when compared to the control simulation (the simulation without hot O).
Hot O + represents a sink for neutral hot O through chemical reactions R 1 and R 2 in table 1.We can use our simulation to evaluate the loss rate of hot O due to the chemical production of hot O + .The results of this evaluation are shown in figure 6.The percentage loss of hot O due to production of

Numerical modeling
The computer simulation that we use to model the effects of hot O on the ionosphere and the behavior of hot O + is TRANSCAR (described in Blelly et al., 1996b).TRAN-SCAR is a time-dependent, 8-moment fluid model of the ionosphere that solves for the density n s , drift velocity u s , temperature T s , and heat flow q s for 6 different ion species (s=O + , H + , N + , N + 2 , NO + , and O + 2 ).TRANSCAR computes a numerical solution to the conservation equations for each of these moments in 1-D, along the geomagnetic field lines.
TRANSCAR assumes the distribution function of the gases that it models to be adequately described with 8 moments.An 8-moment distribution function can describe gases that are somewhat non-Maxwellian through the heat flow vector, but a gas that is highly non-Maxwellian, such as a mixture of hot and cold oxygen, cannot be adequately accounted for in a 8-moment model.In the case of oxygen, the distribution function has a Maxwellian center that rep-resents thermal, or cold O, with a pronounced high-energy tail that represents hot O.We have chosen to break apart the oxygen distribution function into ambient oxygen and hot O and model them as separate species with different distribution functions.
In this paper, we perform 3 simulations with TRANSCAR.One simulation is a control and does not include hot O or hot O + .The second simulation includes just neutral hot O.The third includes both hot O and hot O + .All simulations are done during the month of October, for a latitude of 45 • , f 10.7 =187, and K p =0.

Simulation with neutral hot O
We adopt a diffusive equilibrium density profile for the neutral hot O.The distribution function for the collective O + species is given by: for s 1 = O + and s 2 = O + hot .The moments of the O + distribution function can be calculated in terms of the moments of the separate species by using this distribution function: The heat flow is more complicated but can be found in the same manner.
Analysis of this expression simplifies if we introduce the following variables: Substituting these variables for u, the heat flow equation becomes: where The reference n(z 0 ) is taken to be 1% of the cold O density at z 0 =400 km: The temperature for hot O is chosen to be 4000 K at all altitudes.
Other neutral species are obtained by TRANSCAR from the MSIS model (Hedin, 1991).The high temperature of hot O gives it a large scale height and causes it to become progressively more important than ambient O at higher altitudes.An example of the neutral atmosphere calculated in TRAN-SCAR is shown in Fig. 2. Hot O affects the ions via collisions.These collision terms are added into the momentum, energy, and heat flow equations for all ions in TRANSCAR.
When we run TRANSCAR with the hot O in the simulation, the most pronounced effect is on the ion temperature.Figure 3 shows that even the trace amount (1% at 400 km) of hot O included in the simulation has significant effect on the ion temperature, especially at high altitudes.At lower altitudes, hot O is present only in small concentrations relative to the cooler neutral species densities (Fig. 2), so any energy imparted to an ion by a collision with a hot O atom is quickly lost due to the relatively high collision frequency of the ion with the cooler neutral species.At higher altitudes hot O becomes the most common neutral species, and, thus, has a more significant effect on the ion temperature (Alcaydé et al., 2001).At these altitudes ion collisions with hot O atoms are more frequent than with any other neutral species.Thus, energy imparted to the ions is not lost as quickly via collisions with the cooler neutrals.The result is that, at higher altitudes, the ions are hotter by about ≈300 K in the simulation with hot O than they are in the control simulation.
These simulations exhibit the non-monotonic electron temperature profile that is characteristic of solar max.This decrease in electron temperature near the F-region peak is a result of the very high electron density there during solar max.Since the electron temperature is inversely related to density (Rishbeth and Garriot, 1969), this causes the temperature profile to "dip" at altitudes where the electron density is extremely high.

Simulation of hot O +
In order to self-consistently model hot O + , we now alter TRANSCAR to solve its continuity, momentum, energy and heat flow equations.Hot O + is produced from hot O as a result of photoionization and charge exchange chemistry.The Ann. Geophys., 24, 1625Geophys., 24, -1637Geophys., 24, , 2006 www  chemical model that we use to describe hot O + is shown in Table 1.Reactions R 1 and R 2 represent sources of hot O + and reactions R 1 -R 4 represent losses.Our simulations show that charge exchange production reactions represent a much larger portion of the overall production rate for hot O + than does photoionization.This is not surprising since the ionized reactants for the charge exchange reactions, O + and H + , are the dominant ions at the F-region peak and topside altitudes, respectively.Also, the ionization potentials for O and H are similar and allow for very fast charge exchange.
In addition to the production and loss effects of hot O + on the continuity equations, hot O + affects (and is affected by) other ions through collisions.The terms describing these collisions are now added to the appropriate equations, in much the same manner as were the terms for hot O.The collision cross sections for ion-ion interactions are larger than those of ion-neutral interactions, since the potential for Coulomb interaction (ion-ion) is larger than that for Maxwell molecule interactions (ion-neutral).From these considerations, we predict that hot O + likely has a significant effect on the ambient ion temperature due to frequent collisions.However, the collisional heat exchange process is two-way.The ambient ions also have an equally significant heating effect on hot O + .
A sample plot of the output from TRANSCAR for the simulation including hot O + is shown in Fig. 4. As expected the ambient ion population does have a large effect on the energy dynamics of hot O + .In fact, the ambient ions cause the hot O + to thermalize to the ambient ion temperature at most altitudes, the exception being around the F-region density peak (see Fig. 4).Here ion-ion collision frequencies are the greatest, but production of hot O + is also large.The high production rate results in a net heat input, via the source of Eq. ( 13), into the population and an elevated temperature at the F-region peak.Hot O + forms a layer with a peak density of n s (z pk )≈7.46×10 10 m −3 at an altitude of z pk ≈940 km.These simulated values are reasonably close to the values predicted in section 2, n s (z pk )≈5.62×10 10 m −3 and z pk ≈840 km.According to the simulation, hot O + becomes the major ion at altitudes of 1000 km−2000 km.However, at these altitudes the temperature of hot O + is essentially equal to that of O + .Thus, the two species are indistinguishable in this altitude regime.
Hot O + does not maintain a temperature above that of the rest of the ion population at most altitudes, instead thermalizing to the ambient ion temperature quickly.This result was predicted earlier using an approximate analytic solution to the temperature equation for hot O + .From our simulation Ann.Geophys., 24, 1625Geophys., 24, -1637Geophys., 24, , 2006 www.ann-geophys.net/24/1625/2006/results, we conclude that the ratio P s /n s from Eq. ( 17) is small for hot O + at most altitudes.Earlier we calculated the thermalization time constant to be inversely proportional to the weighted sum of the collision frequencies O + and H + .To see how these calculations are, we have run a simplified hot O + simulation that omits the effects of productions, energy transport, and drag on the hot O + .Using this simulation, we obtain from the time-dependent temperature output a time constant which may be compared with that estimated analytically from Eq. ( 10). Figure 1 shows that these two time constants follow the same general trend with altitude.τ ranges from ≈0.5 min at lower altitudes to ≈2 min at the highest altitudes that we simulate.Hot O + does not thermalize completely near the F-region peak.This is a result of the heat source of Eq. ( 13).At the Fregion peak the hot O + density is low due to loss processes involving H, O, O 2 , and N 2 .However, the O + density is at its peak so there is still a large hot O + production rate from reaction R 2 of Table 1.This results in a significant heat source for the hot O + population and alters its equilibrium temperature to be greater than that of the ambient ions.
The introduction of hot O + into TRANSCAR has a considerable effect on the rest of the simulated ionosphere.Specifically, hot O + is a heat source for the ambient ions.Figure 3 shows that hot O + has the effect of increasing the ambient ion temperature, in addition to the increase imposed by the neutral hot O.At higher altitudes the presence of hot O + has the effect of raising the ambient ion temperature another ≈150 K above the temperature of the simulation with neutral hot O.This is a direct result of the thermalization of hot O + .Hot O + collides with and imparts its energy to the ambient ion population.Consequently, hot O + thermalizes to the ambient ion temperature, but, in the process, causes the ambient temperature to increase.A plot of the heat input to the O + population due to collisions with hot O and hot O + is shown in Fig. 5.This figure shows that hot O and hot O + are comparably important heat sources for the ambient ions.The simulation with hot O + also shows small but noticeable effects on the densities of O + and H + (due to reactions R 1 , R 1 , R 2 , and R 2 in   ion will contain a significant contribution from the differential drift.The simulations show that while the hot O + particles may not constitute a separate species, they may manifest themselves through an enhanced heat flow in the overall O + population, via the differential drift velocity.However, hot O + and ambient O + do not maintain enough drift velocity difference (even during sunrise and sunset) to cause the overall heat flow to deviate substantially from the sum of the individual heat flows.

Thermalization as a loss process
To this point the hot O + has been treated as a completely separate species from O + .However, when hot O + thermalizes to the ambient O + temperature there is no feature that differentiates the hot O + from the ambient O + .The thermalization can be considered a loss process for the hot O + .If hot O + is modeled this way, then the density will be lower than in our simulations and this will affect the equilibrium temperature of the hot O + , since the ratio P s /n s will be larger (Eq.17).
We present here a simulation of this alternative model of hot O + .We model the effects of thermalization loss by adding in two cooling reactions, R th,1 and R th,2 of Table 2.
The total loss frequency for these reactions, β th , is taken to be the inverse time constant for thermalization that was derived in section 3.
To fully model the thermalization reactions of Table 2, we must also include the energy input into the O + and H + populations due to these reactions.The energy input has the form of Eq. (13): For this case, s=O + , H + .Figure 7 shows the TRANSCAR results for the ionosphere with hot O + thermalization losses included.Figure 7 shows that the new thermalization loss reactions cause the hot O + density to be much lower.As a consequence the heat source of Eq. ( 13), which is inversely proportional to density, is large enough to keep the hot O + at a significantly higher temperature than the ambient ions (Fig. 7).However, the hot O + would still not be detectable, since it has such a low density.The heating effects of hot O + on the ambient O + now include two components: collisional heating, δT s /δt, and a production heating (δT s /δt) * .Our simulations show that if these two components are added together, they equal the hot O + heating rate computed from our original model of hot O + , Fig. 5.

Conclusions
Our findings indicate that hot O + does not maintain a temperature above that of the ambient ions, but is an important heat source for the ambient ions.
We have made approximate calculations describing the behavior of hot O + .Also we have made additions to an existing simulation of the ionosphere (TRANSCAR) in order to investigate the possible effects of a hypothesized neutral hot O population and a self-consistently modeled hot O + population.Since it has a considerably higher temperature than the ambient ions, neutral hot O serves as a significant heat source to the ionosphere at all altitudes (Fig. 3).However, hot O + cannot maintain an elevated temperature at most altitudes due to thermalization by the other ions (Fig. 4).Through this process of thermalization, hot O + represents a heat source (in addition to the hot O heat source) for the ambient ions (Fig. 5). Figure 3 shows that the inclusion of hot O + causes the ambient ion temperature to be higher than in the control simulation (no hot O, no hot O + ) and the simulation with only hot O. From Figs. 3 and 5 it can be seen that the combined effect of the hot O and hot O + on the energy dynamics of the ionosphere is considerable.
Figure 4, which is representative of the noontime, steadystate ionosphere, shows that hot O + is the dominant ion at 1000 km−2000 km altitude.The temperature of hot O + is approximately equal to the ambient ion temperature (Fig. 4) at these altitudes, so that there is really no distinction between hot O + and O + .The ionization of hot O and the subsequent thermalization of hot O + appears to be an efficient heating process in the ionosphere.
We have also considered an alternative way to model hot O + that includes the thermalization process as a loss process in the continuity equation.This alternate model gives a significantly different density and temperature profile.However, the heating effects of hot O + on the ambient ions remains the same as before.Since detecting hot O + does not appear to be feasible, the heating effects, which are the same under both models, are the most important aspect of hot O + .

Fig. 1 .
Fig. 1.Comparison of calculated and simulated thermalization time constants for hot O + .The chemical time constant (equation 1) is also shown for comparison

Fig. 1 .
Fig. 1.Comparison of calculated and simulated thermalization time constants for hot O + .The chemical time constant (Eq. 1) is also shown for comparison.

Fig. 2 .
Fig. 2. The Neutral Atmosphere From TRANSCAR.n s is the number density for species s.T O hot is the hot O temperature and T N eutral is the other neutral particle temperatures from the MSIS model.

Fig. 2 .
Fig. 2. The neutral atmosphere from TRANSCAR.n s is the number density for species s.T O hot is the hot O temperature and T Neutral is the other neutral particle temperatures from the MSIS model.

Fig. 3 .
Fig. 3. Comparison of Temperatures Calculated by TRANSCAR for 3 different simulations: control sim., sim.with hot O, and sim.with hot O and hot O + .T N eutral is the cold neutral temperature, T O + is the O + temperature, and T e is the electron temperature.

Fig. 3 .
Fig. 3. Comparison of temperatures calculated by TRANSCAR for 3 different control sim., sim.with hot O, and sim.with hot O and hot O + .T Neutral is the cold neutral temperature, T O + is the O + temperature, and T e is the electron temperature.

Fig. 4 .Fig. 5 .Fig. 6 .
Fig. 4. Ion densities n s and temperatures T s as simulated in TRANSCAR.This simulation includes our self-consistently modeled population of hot O + .Fig. 4. Ion densities n s and temperatures T s as simulated in TRANSCAR.This simulation includes our self-consistently modeled population of hot O + .

Fig. 5 .
Fig. 5. O + noontime heating rates due to collisions with hot O, hot O + , and e − .The total steady-state cooling rate (from collisions with neutral is equal to the total heating rate shown above.

Fig. 7 .
Fig.7.Ion densities n s and temperatures T s as simulated in TRANSCAR.This simulation includes our self-consistently modeled population of hot O + , with thermalization losses included.Fig.7.Ion densities n s and temperatures T s as simulated in TRANSCAR.This simulation includes our self-consistently modeled population of hot O + , with thermalization losses included.

Fig. 7 .
Fig.7.Ion densities n s and temperatures T s as simulated in TRANSCAR.This simulation includes our self-consistently modeled population of hot O + , with thermalization losses included.Fig.7.Ion densities n s and temperatures T s as simulated in TRANSCAR.This simulation includes our self-consistently modeled population of hot O + , with thermalization losses included.
Temperature rate of change, δT O +/δt, for different heating processes (K s −1 ) Fig.5.O + noontime heating rates due to collisions with hot O, hot O + , and e − .The total steady-state cooling rate (from collisions with neutral atoms) is equal to the total heating rate shown above.

Table 1 ,
and changes in scale height) when compared to the control simulation (the simulation without hot O).Hot O + represents a sink for neutral hot O through chemical reactions R 1 and R 2 in Table1.We can use our simulation to evaluate the loss rate of hot O due to the chemical production of hot O + .The results of this evaluation are www.ann-geophys.net/24/1625/2006/Ann.Geophys., 24, 1625-1637, 2006 Matthew Zettergren, et al: Modeling hot oxygen ions 11