Rovibrationally Resolved Time-Dependent Collisional-Radiative Model of Molecular Hydrogen and Its Application to a Fusion Detached Plasma

Abstract: A novel rovibrationally resolved collisional-radiative model of molecular hydrogen that includes 4,133 rovibrational levels for electronic states whose united atom principal quantum number is below six is developed. The rovibrational XΣg population distribution in a SlimCS fusion demo detached divertor plasma is investigated by solving the model time dependently with an initial 300 K Boltzmann distribution. The effective reaction rate coefficients of molecular assisted recombination and of other processes in which atomic hydrogen is produced are calculated using the obtained time-dependent population distribution.


Introduction
The collisional-radiative model (CR model) [1] provides population distributions of atoms and molecules.Using this model, the effective reaction rate coefficients of various atomic and molecular processes including the contributions of the excited states can be calculated based on the population distribution.Emission intensities of atoms and molecules are also calculated using the population distribution.
We have been developing a CR model of molecular hydrogen, H 2 .Previously, a CR model for molecular hydrogen in which only the electronic states are considered was constructed by Sawada and Fujimoto [2].This model was used to calculate effective reaction rate coefficients of dissociation and ionization of the molecule in plasmas.This CR model was extended by Greenland and Reiter [3] and Fantz and Wündrlich [4], and the rate coefficients calculated by these models were included in the neutral transport code EIRENE [5].
Molecular-assisted recombination (MAR) was then proposed to understand the detached recombining plasmas in fusion edge plasmas [6].MAR comprises two series of processes.The first starts with the production of H + 2 from the collision of H 2 and H + , and the second starts with the dissociative attachment of electrons.To evaluate the contribution of the molecular vibrational distribution to the rate coefficient of MAR, we constructed a model in which the vibrational states were included in addition to the electronic states [7].
It has recently been recognized that the production of state-resolved cross-sections is essential for the detailed analysis of fusion divertor plasmas [8].In this study, we develop a CR model in which the electronic, vibrational and rotational states are considered in order to model molecular processes whose cross-sections strongly depend on the initial vibrational and rotational states, e.g., the dissociative attachment of electrons in MAR.
As for the emission spectroscopy of molecular hydrogen, a corona model [1] for the Fulcher transition (d 3 Π u → a 3 Σ + g ) that includes vibrational states was constructed to determine the vibrational temperature in the electronic ground state X 1 Σ + g [9].Although the model was developed to also include rotational states in order to determine both rotational and vibrational temperatures in X 1 Σ + g [10], corona models may not be suited to high electron density fusion detached plasmas because they neglect electron impact transition among the excited states.To describe high electron density plasmas, a rovibrationally-resolved CR model is necessary.In addition to the Fulcher band, this model provides emission line intensities for many bands, which may be used to determine the electron temperature and density as in helium atom spectroscopy [11] and to determine the rovibrational population distribution in the X 1 Σ + g without assuming the vibrational and rotational Boltzmann distributions.The purpose of this paper is to introduce a newly-developed CR model and to provide examples of its calculation.In this study, we investigated the rovibrational population distribution in a fusion detached divertor plasma.Because the quasi-steady-state solution [1] for rovibrational states in X 1 Σ + g cannot be applied to the rate equations in the CR model owing to the large population relaxation time, we solved the rate equations time dependently for these states.
In this paper, we first introduce our newly-developed state-resolved CR model of molecular hydrogen and then list considered energy levels and relevant data as well as spontaneous transition probabilities, cross-sections and rate coefficients used in the model.Finally, we report on the modeling of a SlimCS fusion demo detached divertor plasma [12,13] in which effective reaction rate coefficients including MAR are calculated using the rovibrational population distribution.An example of emission spectra is also shown.

Framework of the Collisional-Radiative Model
To provide a framework of the CR model itself, we describe a model of atomic hydrogen [1,14,15] for simplicity.Here we assume that only the atomic hydrogen, electrons and protons exist in a plasma.The CR model gives the population of an excited state p, n(p).The temporal variation of the population density of an excited state p can be described by a rate equation as: where n e and n H + are the electron and proton densities, respectively, C(p, q) is the excitation rate coefficient for electron collisions from state p to q and F(q, p) is the inverse de-excitation rate coefficient.
A(p, q) is the spontaneous transition probability from p to q. S(p) is the ionization rate coefficient for state p, and α(p) and β(p) are three-body and radiative recombination rate coefficients, respectively.Each rate coefficient is a function of the electron temperature T e .In many cases, the quasi-steady-state solution (QSS) [1] can be used to approximate the left-hand side of Equation (1) as zero for all states except the ground state: ( By solving Equation (1) with Equation (2), this is a set of coupled linear equations, the population density of an excited state p in the form of: is obtained, where R 0 (p) and R 1 (p) are the population coefficients, each of which is a function of n e and T e .The first and the second terms are referred to as the recombining and the ionizing components, respectively [1].
The population of the ground state n(1) should be solved time dependently using the population coefficients for the excited state in Equation (3) [15]: where the first and the second terms in the right-hand side of Equation ( 4) denote inflow to and outflow from the ground state, respectively.Equation ( 4) is simply rewritten as by defining the effective recombination and ionization rate coefficients, α CR and S CR , respectively.Both α CR and S CR are functions of n e and T e .We define τ(p) as the reciprocal of the outflow rate in Equation (1): The value of τ(p) gives the relaxation time after the inflow to the level p becomes steady.The relaxation time of all excited levels is given approximately by the maximum value of τ(p) in the excited levels [15], which we define as τ max .In applying the QSS for the excited states, the environment of the excited states, n(1), n e , n H + , and T e , should be steady on the time scale of τ max .In many cases, this condition is satisfied; the ground state τ( 1) is much larger than τ max because there is no spontaneous transition.For more details, see [15].

Collisional-Radiative Model for Molecular Hydrogen
In the case of molecular hydrogen, the values of τ(p) for the rovibrational states in the electronic ground state X 1 Σ + g are nearly the same, and the populations of all of these states should therefore be solved time dependently.In this study, the populations of c 3 Π − u (v = 0), which is known as the metastable state, and the a 3 Σ + g , which couples strongly to the c 3 Π − u (v = 0) state, are solved time dependently in addition to the population of X 1 Σ + g .In the application of this study, as shown later, we deal with the population of molecular hydrogen originating from the electron impact excitation of the ground state X 1 Σ + g , which corresponds to the ionizing plasma component in Equation (3).

Elementary Processes
Here, we introduce elementary processes included in the molecular hydrogen CR model.The projectiles considered are listed in Table 1.

Target Projectile
H 2 e, H + , H, H 2 In the following, references are given for the data used; for discussion of their respective reliabilities, the original papers should be consulted.We will provide comments regarding data with large uncertainty.

Energy Levels of H 2
In this study, the energy levels of molecular hydrogen are labeled e, v, N, and J based on Hund's coupling Case (b) [16,17].Here, e is the set of quantum numbers characterizing the electronic state; v is the vibrational quantum number; N is the rotational quantum number of the total angular momentum excluding electron and nuclear spins; and J is the quantum number of the total angular momentum excluding nuclear spin.For the statistical weights and allowed transition in the rotational levels, see [16].For electronic states with Λ > 0, the two levels associated with the two possible values ±|Λ| are considered.Here, Λ is a projection of total electronic orbital angular momentum about the internuclear axis [16].
The energy level data are taken from [18][19][20].Table 2 lists references adopted in the present model.For triplet data in [18], 149.6 cm −1 is subtracted according to [21,22]. Figure 1 shows the electronic states whose energy levels are given in [18][19][20].4,133 rovibrational states with the united atom principal quantum number n ≤ 6 is included in the present model.

State
United Atom Energy Level Potential [25] [29] ) ) ) The energy level of the rovibrational ground state for each electronic state in [18][19][20] is shown.
For the repulsive state b 3 Σ + u , the energy level at the nuclear distance of the potential minimum of X 1 Σ + g is shown.In the present model, states up to n = 6 are considered.
We calculate the rovibrational energy levels for which data are not available using the electronic potential as follows.

Electronic Potential Data
A review of electronic state potentials is available in [32], and Table 2 lists the potential references adopted in our model.Figure 2a-c shows the potential curves from the references in Table 2.
To calculate the energy level and the Franck-Condon factor [16], we solve the Schrödinger equation for Hund's (b) case [16]: where ψ evN (R) is the wave function.This is solved using the Numerov method [33].

Spontaneous Transition and Predissociation
The bound-bound transition probability data in Table 3 are included in the model.We calculated other transition probabilities using the following equation [34]: where the labels e, v, J and e , v , J denote the initial and final electronic, vibrational and rotational states, respectively, and S J J is the Hönl-London factor [16], which is given in [35].In the present model, the transition moments [25,27,29,[36][37][38] in Table 4 are used, while the united helium atom transition moments are used for other transitions.
Table 3. References for the transition probability.
Transition Present Table 4. References for the transition moments.

Transition
Present We included the transitions in Table 5, which are accompanied by the continuum emission.The transition probability from a discrete vibrational state v to a continuum vibrational state with wavenumber k is calculated as [45][46][47]: where R is the nuclear distance and ν vk is the photon frequency in cm −1 .The wave functions ψ v (R) and ψ k (R) pertain to the upper and lower states, respectively.D(R) is the transition moment between the electronic states; ν vk is the frequency of the emitted light; and g is the statistical weight.Table 5 lists references for D(R) used in the present model, and ψ k (R) is calculated by applying the Numerov method [33] to Equation (7), and ψ k (R) is normalized at large R according to [46,47].
Table 5. References for the transition moment of the bound-continuum.
Transition Present Predissociation lifetimes from the metastable state c 3 Π − u (v = 0) to the b 3 Σ + u , as given in [48], are included in the model with values on the order of 10 −3 s.

Electron Impact Electronic Excitation
3.4.1.Excitation from the Electronic Ground State X 1 Σ + g Tables 6 and 7 show the references for the electron impact cross-sections from the X 1 Σ + g to the electronic excited states.Cross-sections to the B 1 Σ + u and C 1 Π u are shown in Figure 3 as an example.Aside from some data from [49], we use semi-empirical data from [50,51], although the reliability of these data has not been closely examined.Excited State e→e' ev→e' ev→e'v' Table 7. References for electron impact excitation cross-section from X 1 Σ + g to various triplet states.
Excited State e→e' ev→e'

Electron Impact Excitation between Electronic Excited States
Table 8 shows available data on electron impact excitation between excited electronic states.As there are no available data for other transitions, we used cross-sections for the united atom helium (see Table 2) calculated using the convergent close-coupling (CCC) method [56].The helium atom cross-sections corresponding to the transitions in Table 8 differ from those for the molecular transitions within a factor of 0.1-10, this should be improved in the future.

Rovibrationally-Resolved Excitation
The state-resolved electron impact rate coefficient from (e, v, N) to (e , v , N ), R e v N evN subject to the Born-Oppenheimer approximation, is written as [17]: where q e v N evN is the Franck-Condon factor.And σv e e is the rate coefficient for the transition between the electronic states.δ as is one or zero according to the selection rule from the symmetry of nuclear spin.The factor a N N dominantly determines the rotational dependence.For the transition from the X 1 Σ + g to the d 3 Π u , a N N is given by [17] as: where Q r (r = 1, • • • , 4) have experimental values 0.76, 0.122, 0.1 and 0.014, respectively.The rate coefficient value is distributed to allow J states according to their statistical weights.For other transitions, we distribute the rate coefficient into the final J states within |N − N | = 0, 1 instead of a N N in Equation (10).For transitions from the X 1 Σ + g to the B 1 Σ + u and C 1 Π u , we use the rate coefficients, which are vibrationally resolved for the upper and lower electronic states (see Table 6).

Processes Involving H + 2 Table 9 lists processes involving H +
2 included in the present model.For ionization from the excited electronic states, we use the united atom helium data [56].Figure 4 shows the cross-sections for the second process in Table 9. Transition ev→e' ev→e'v' Data are taken from [62,63].
In [54], the cross-sections for g state are also given.Their rate coefficients are included in the present model.Rate coefficients for other (v, transitions are obtained by assuming that their ratios as the same as their cross-sections at 4 eV given in [54].We assume that the rate coefficients for (v, J = 0) → (v + 1, J = 0) (v > 6) are the same as that for (v = 6, J = 0) → (v = 7, J = 0).We include the dissociation process from In the transitions in [54], the rotational state is fixed.We include (v = 0, J) → (v = 1 ∼ 14, J = J, J ± 2) and (v, J) → (v = v + 1, J = J, J ± 2) in the model assuming the cross-sections for ∆J = 0, ±2 transition are the same.
Vibrational excitation cross-sections for all combinations of (v) → (v ) are calculated in [66] by considering H − 2 (X 2 Σ + g ) resonance.∆J = 0, ±2 transition is included neglecting the rotational dependence of the cross-section.
3.8.Vibrational and Rotational Excitation in the H 2 (X 1 Σ + g , v, J) by Proton Impact For a pure rotational excitation, we include cross-sections for the (v = 0, J = 0) → (v = 0, J = 2) and (v = 0, J = 1) → (v = 0, J = 3) transitions from [55] (see Figure 3) in the H 2 CR model.Transitions for other ∆J = 2 in v = 0 are also considered under the assumption that their cross-sections are the same as that for (v, = 0, J = 1) → (v = 0, J = 3).For other vibrational states, we assume that rate coefficients for the pure rotational excitations are the same as those for v = 0.
3.9.Rotational Excitation in the H 2 (X 1 Σ + g , v = 0) by Molecular Hydrogen Impact As shown in Figure 5, a rate coefficient dataset for the rotational excitation in H 2 (X 1 Σ + g , v = 0) arising from the collision of two molecules is given in [19]: the rate coefficients for collisions between a molecule with J ≤ 8 and a molecule with J = 0 or 1 are given for 2 K ≤ T H 2 ≤ 10,000 K. Data in [19] are included in the present model.We use the data from the H 2 (X 1 Σ + g , v = 0) state for other vibrational states because no data are available.

by Atomic Hydrogen Impact
A cross-sectional dataset for the rovibrational excitation in X 1 Σ + g owing to the collision of atomic hydrogen is given in [67] (see Figure 6).Rate coefficients in the present model are calculated using this dataset for transitions among v = 0-6, J = 0-23 in the barycentric collision energy range 200-60,000 K. 81-)9:,, Data are taken from [67].

by Atomic Hydrogen Impact
Vibrationally-resolved rate coefficients for the dissociation of H 2 (X 1 Σ + g , v) by atomic hydrogen impact given in [68] (Figure 7) are included in the present model.Here, rotational dependence is neglected because no relevant data are available.

Dissociative Attachment by Electron Impact
As shown in Figure 8, the dissociative attachment rate coefficient is calculated in [69] by considering H − 2 (X 2 Σ + u ) resonance.The complete dataset from all rovibrational states in the H 2 (X 1 Σ + g , v, J) is included in the CR model.The dissociative attachment cross-section from all vibrational states in H 2 (X 1 Σ + g , v) calculated in [70] by considering H − 2 (X 2 Σ + g ) resonance is also included in the CR model.In this case, the rotational dependence is neglected because of a lack of data.

by Electron Impact
The dissociative processes in Table 10 are included in the model in order to determine the population of the H 2 (X 1 Σ + g , v, J); once again, rotational dependence is neglected.
Table 10.References for the electron impact dissociative processes starting from H 2 (X 1 Σ + g ).

Results and Discussion
We applied the H 2 CR model to a fusion detached plasma as an example of the application of the present model using the following parameters taken from a SlimCS DEMO detached plasma simulation [12,13]: T e = T H + = T H = 2 eV, n e = n H + = 10 16 cm −3 and n H = 10 15 cm −3 .In the present calculation, background molecular hydrogen, which contributes to the rotational excitation, is limited to H 2 (X 1 Σ + g , v = 0, J = 0, 1), and the kinetic temperature and total density are fixed at T H 2 = 1 eV and n H 2 = 10 15 cm −3 , respectively.
The present code solves the rovibrational populations of the X 1 Σ + g , a 3 Σ + g and c 3 Π − u time dependently.The values of τ(p) from Equation ( 6) for the rovibrational states of X 1 Σ + g are 1 × 10 −9 ∼ 1 × 10 −8 s.The values of τ max in the electronic excited states is 2 × 10 −9 s for a level in the a 3 Σ + g .Except for the a 3 Σ + g , τ is smaller than 1 × 10 −10 s.In this case, the c 3 Π − u (v = 0), known as the metastable state, does not have a significantly larger τ(p) than the other levels owing to the large value of n e .The Runge-Kutta method was used with the time step of ∆t = 10 −11 s.The QSS was applied to the rovibrational states in the other electronic excited states.
Figure 9 shows the population of the rotational states in H 2 (X 1 Σ + g , v = 0) as a function of time.The initial population of 300 K Boltzmann distribution is given at t = 0 s.Initial total H 2 density was set at 1 cm −3 .In the present case, the distribution may be understood to be that of a gas puffed molecule.In modeling the recycling of a molecule at a wall, a proper rovibrational distribution should be given.
Population divided by the statistical weight of the rotational states in Figure 10a-c shows the rovibrational population distributions at t = 10 −10 s, t = 10 −8 s and t = 10 −7 s, respectively.Figure 10b,c suggests a different populating mechanism from that in Figure 10a.
Figure 11a-c shows the inflow to each rovibrational level at t = 10 −10 s, t = 10 −8 s and t = 10 −7 s, respectively.At t = 10 −10 s, the dominant inflow to each level of 0 ≤ J ≤ 8 is the electron impact transition in H 2 (X 1 Σ + g , v, J).The origin of the shoulder-like profile for J > 8 in Figure 11a is the atomic hydrogen impact transition from lower lying J = 0, 1 states with the large population.The J dependence of the contributions of molecular and atomic collisions can be attributed to the differences in the J dependence of the rate coefficients or cross-sections in Figures 5 and 6.The transition with photoemission to the H 2 (X 1 Σ + g , v, J) is limited to J = 14 because the upper limit of J in the electronic excited states is J = 13 in the B 1 Σ + u (v = 2) state.As shown in Figure 11b,c, increasing the excited rovibrational state population leads to an increase in electron impact transition among these states.
Figure 12a-c shows the outflow from each rovibrational level at t = 10 −10 s, t = 10 −8 s and t = 10 −7 s, respectively.As the highly excited rovibrational state population increases, the dissociative attachment increases.Finally, as shown in Figures 11c and 12c, the population distribution is dominated by the large flow in the X 1 Σ + g from electron impact and quenching processes including the dissociative attachment.After t = 10 −7 s, the total population decreases with the relative rovibrational population in X 1 Σ + g kept nearly constant.g , (open triangle) proton collision transition in X 1 Σ + g , (closed circle) atomic collision transition in X 1 Σ + g , (closed square) molecular collision transition in  Figure 14 shows the effective rate coefficients for MAR (P1-P3) and (P4-P7) in Table 11 and other processes in which atomic hydrogen is produced (P8-P12) and (P14-P15).Here, the effective rate coefficients are defined by: dn H dt where P H H 2 (P1-P12) and P H H + (P14-P15) are the effective rate coefficients for the production of atomic hydrogen from H 2 and H + , respectively, including the contributions of their excited states, n H is the population density of the ground state atom, n(1) in Equation (4), n H 2 is total population density of H 2 including the ground and excited states, and P H H + is the same as α CR in Equation ( 5).The effective rate coefficients are calculated using cross-section data and models in the references in Table 11, in addition to the present H 2 CR model.In (P3, P7, P12, P15), the atomic hydrogen CR model [14,15] is used to calculate population flows into the ground state from excited states.The vibrationally-resolved CR model for H + 2 [7] that includes processes (P4-P6) and (P16) in Table 11 is used to calculate the population distribution of H 2 + (X 2 Σ + g , v) assuming the QSS.Rotational levels are not considered in the H 2 + CR model.For processes (P11) and (P12), we use dissociative excitation data derived from experimental emission cross-sections given in [2,71] instead of those in the references in Table 10 in order to include the contribution from all of the excited states of the H 2 molecule.The rovibrational dependence is not included in the present calculation owing to a lack of relevant information.Figure 14 shows that the effective rate coefficients for (P1-P3) and (P4-P7) increase with the excited rovibrational population in X 1 Σ + g .For t > 10 −7 s, the sum of the effective rate coefficients of the two types of MAR is larger by a factor of approximately 25 than that at t = 10 −11 s.Furthermore, MAR cannot be neglected relative to the well-known dissociative process (P8) as a source of atomic hydrogen after time t > 10 −8 s.
In this study, we used the rate coefficient of dissociative attachment based on the nonlocal resonance model [69].According to this model, it is pointed out that non-resonant scattering enhances the rate coefficient [72] and, although a complete dataset including it is not available, this effect may enhance the effective rate coefficient of (P1-P3) by an approximate factor of two beyond the value calculated here.
The photoemission transition from triplet electronic excited states (P10) (see Table 5) is also important for atom production.The populations of the upper states of the photoemission transitions in (P9) and (P10) increase up to t ∼ 10 −8 s.In a later time, these populations decrease because the relative population of lower rotational states in X 1 Σ + g , which are the sources of the populations of the upper states of the photoemission, decreases; the electron impact excitation for ∆J = 0, ±2 is considered and the number of the rotational states of the upper state is smaller than that of the X 1 Σ + g .In this calculation, as the quenching from higher lying rotational states, only the dissociative attachment is effectively considered, because the rotational dependence of the cross-section is not included for other quenching processes.Furthermore, the electron impact dissociation through the X 1 Σ + g continuum, which may play an important role, is not included here because no relevant data are available.If the dissociative attachment process is not dominant, the population density becomes smaller, which in turn reduces the effective rate coefficient for the MAR starting from the attachment process.The quenching cross-section from the higher rotational states should be studied in detail in future research.For precise calculation of the effective rates, precise cross-section data for the various transitions in X 1 Σ + g will also be indispensable.

Label
Process Reference  Figure 15a-c shows an example of the emission spectra of molecular hydrogen calculated from the obtained rovibrational population distribution of the excited electronic states by multiplying the spontaneous transition coefficient and photon energy.At present, the uncertainty in the electron impact excitation cross-section between the electronic states is large.For qualitative spectroscopic diagnostic, a more reliable dataset is necessary.
Here, we would like to discuss the application of the model presented in this paper to more general cases.In this study, we addressed processes starting from the H 2 (X 1 Σ + g , v, J).In future work, a model for the recombining component starting from H 2 + , (P19-P20) in Table 11, should be constructed for lower electron temperature plasmas.Furthermore, recombination starting from atomic hydrogen, (P21) in Table 11, may contribute to the rovibrational population of the H 2 (X 1 Σ + g , v, J) in high atomic density plasmas [77].In considering the recycling of molecules at a wall (P22), a proper initial population distribution at the wall should be given to the model.As for the present calculation, interpreting it as that for the history of a gas puffed molecule, the result is valid regardless of the negligence of the above processes.Figure 16 shows the rate coefficients of (P6) + (P16) and (P17) in Table 11.When the electron density is much smaller than the molecular density, H + 3 may be produced by (P17), which is a competing process of (P6) + (P16).In this case, the dissociative process from H + 3 (P18) may contribute to the population of H 2 (X 1 Σ + g , v, J).Production of the rovibrationally-resolved cross-section for this process is indispensable for such plasmas.

Conclusions
In this study, we constructed a rovibrationally-resolved collisional-radiative model of molecular hydrogen and applied it to a demo fusion detached divertor plasma.Using a time dependently solved population density distribution, the effective rate coefficients for MAR and the production of atomic hydrogen processes were calculated.The present model can provide rate coefficients for the purpose of tracing molecules in the rovibrational states in the X 1 Σ + g as different species in neutral transport codes.Molecular spectra intensity is a function of n H2 , n e , T e and the rovibrational population distribution of the X 1 Σ + g .It was shown that the present CR model has the potential to determine these parameters from observed spectra.For this purpose, the population distribution of the excited electronic states should be calculated precisely by obtaining a reliable dataset for the electron impact excitation cross-sections between the electronic states.In this paper, we listed the used elementary process data.This will give some information on the necessary data.The present rovibrationally-resolved CR model will be utilized as a storage box for the rovibrationally-resolved elementary process data to be provided successively in the future.In addition to the above uses, the rovibrationally-resolved H 2 CR model can also be applied to the precise determination of the energy loss of projectile particles from collision with H 2 in plasmas.We are constructing related models for isotope species D 2 , T 2 , HD and DT that will be reported elsewhere. '-(78+'*-,9):&*(8',,,0+1,9:/+.10;-<*.,:+0004=6

Figure 3 .
Figure 3.Some examples of the cross-sections for electron and proton impacts included in the present model.Data are taken from [49,53-55].

Figure 11 .
Figure 11.Inflow to each rovibrational level divided by statistical weight.The total inflow from all initial states is shown.(a) t = 10 −10 s; (b) t = 10 −8 s; and (c) t = 10 −7 s. (Open circle) Emission from the electronic excited states, (open square) electron collision transition in X 1 Σ +g , (open triangle) proton collision transition in X 1 Σ + g , (closed circle) atomic collision transition in X 1 Σ + g , (closed square) molecular collision transition in X 1 Σ + g .

Figure 12 .
Figure 12.Outflow from each rovibrational level divided by statistical weight.The total outflow to all final states is shown.(a) t = 10 −10 s; (b) t = 10 −8 s; and (c) t = 10 −7 s. (Open square) Electron collision transition in X 1 Σ + g , (open triangle) proton collision transition in X 1 Σ + g , (closed square) molecular collision transition in X 1 Σ + g , (closed circle) atomic collision transition in X 1 Σ + g , (closed triangle) electron collision to the excited electronic states, (plus) dissociative attachment, (open circle) dissociation through the excited electronic states including b 3 Σ + u , (cross) charge transfer to H + , (open diamond) dissociation by atomic hydrogen impact.

Figure 13 Figure 13 .
Figure13shows the vibrational population distributions obtained by summing all rotational populations.The distribution is not expressed as a simple Boltzmann distribution.

Table 1 .
Target and projectiles in this study.

Table 2 .
References for energy levels and potentials.

Table 6 .
References for electron impact excitation cross-section from X 1 Σ + g to various singlet states.

Table 8 .
References for electron impact excitation cross-sections between electronic excited states.

Table 11 .
Molecular-assisted recombination (P1-P3 and P4-P7) and some other processes.p is the principal quantum number of the H atom. H * and H * 2 denote the excited atom and molecule, respectively.

Figure 14 .
(12)ils of P H H 2 in Equation(12).Effective rate coefficients for molecular-assisted recombination (MAR) and other processes for the production of atomic hydrogen (P1-P12) and (P14-P15) in Table 11 are shown.The P H H 2 for (P13), defined as dn H dt = P H H 2 n H 2 n H , is also shown.The values for (P8), (P9), (P10) and (P13) in this figure should be multiplied by two because two atoms are produced.