Estimating inelastic heavyparticle – hydrogen collision data
II. Simplified model for ionic collisions and application to bariumhydrogen ionic collisions^{⋆}
^{1} MaxPlanck Institute for Astrophysics, Postfach 1371, 85741 Garching, Germany
email: andrey.k.belyaev@gmail.com
^{2} Department of Theoretical Physics and Astronomy, Herzen University, 191186 St. Petersburg, Russia
Received: 1 September 2017
Accepted: 8 October 2017
Aims. A simplified model is derived for estimating rate coefficients for inelastic processes in lowenergy collisions of heavy particles with hydrogen, in particular, the rate coefficients with high and moderate values. Such processes are important for nonlocal thermodynamic equilibrium modeling of cool stellar atmospheres.
Methods. The derived method is based on the asymptotic approach for electronic structure calculations and the LandauZener model for nonadiabatic transition probability determination.
Results. It is found that the rate coefficients are expressed via statistical probabilities and reduced rate coefficients. It is shown that the reduced rate coefficients for neutralization and ionpair formation processes depend on single electronic bound energies of an atomic particle, while the reduced rate coefficients for excitation and deexcitation processes depend on two electronic bound energies. The reduced rate coefficients are calculated and tabulated as functions of electronic bound energies. The derived model is applied to bariumhydrogen ionic collisions. For the first time, rate coefficients are evaluated for inelastic processes in Ba^{+} + H and Ba^{2+} + H^{−} collisions for all transitions between the states from the ground and up to and including the ionic state.
Key words: atomic data / atomic processes / stars: atmospheres
Tables with calculated data are only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/608/A33
© ESO, 2017
1. Introduction
Nonlocal thermodynamic equilibrium (nonLTE) modelings of different astrophysical phenomena are important for many fundamental problems in modern astrophysics (see, e.g., reviews Asplund 2005; Mashonkina 2014; Barklem 2016a, and references therein). In particular, nonLTE modelings of stellar atmospheres are of importance, for example, for determining absolute and relative abundances of various chemical elements, for the Galactic evolution, and so on. NonLTE modeling requires detailed and complete information about inelastic heavyparticle collision processes, most importantly those in collisions with hydrogen atoms and negative ions. The lack of atomic data on inelastic processes due to collisions with neutral hydrogen atoms and hydrogen negative ions has been a major limitation on nonLTE modeling of stellar spectra. This problem has been well documented (see, e.g., Lambert 1993; Asplund 2005; Mashonkina 2014; Barklem 2016a, and references therein). Recently, significant progress has been achieved through detailed full quantum scattering calculations based on accurate ab initio quantum chemical data for the cases of simple atoms, such as He, Li, Na, Mg, and Ca (Belyaev et al. 1999, 2010, 2012; Croft et al. 1999a,b; Belyaev & Barklem 2003; Guitou et al. 2011, 2015; Belyaev 2015; Mitrushchenkov et al. 2017). This allowed one to perform accurate nonLTE modeling of stellar spectra (Barklem et al. 2003; Lind et al. 2011; Osorio & Barklem 2016; Mashonkina et al. 2017). Nevertheless, full quantum ab initio and scattering calculations are still time consuming and, hence, rather seldom.
Owing to high demand in atomic data for inelastic processes in collisions with hydrogen and a large number of partial processes and atomic states to be taken into account, astrophysicists were forced to use widely the classical Drawin formula (Drawin 1968, 1969; Drawin & Emard 1973; Steenbock & Holweger 1984; Lambert 1993). Later on, it was recognized that the Drawin formula is not reliable (Barklem et al. 2011; see also Asplund 2005), and it was emphasized that it is highly desirable to derive an approximate model approach to inelastic atomic collisions with hydrogen. Low computational cost are also considered important.
In nonLTE modelings many other processes, including radiative processes and other collisional processes, especially those with electrons, are taken into account. A collisional process with hydrogen is only of importance, if the rate coefficient is comparable to those for the other processes (see, e.g., Mashonkina 2014; Barklem 2016a). Thus, a desirable model approach to inelastic collisions with hydrogen should provide rate coefficients with high and moderate values.
Two such approaches for nonadiabatic nuclear dynamics have been proposed: the quantum branching probability current method (Belyaev 2013a) and the quantum multichannel analytic approach (Belyaev 1993; Belyaev & Barklem 2003; Belyaev et al. 2014; Yakovleva et al. 2016). Both approaches are based on electronic structure calculations and the LandauZener model. For electronic structure calculations, in addition to ab initio methods, the approximate methods have been used: the asymptotic method (Belyaev 2013a) and the Linear Combinations of Atomic Orbitals (LCAO) method (Barklem 2016b, 2017). The quantum model approaches have been successfully applied to a number of chemical elements relevant to nonLTE modeling, such as Al, Si, Ca, Be, and Mn, (Belyaev 2013b; Belyaev et al. 2014, 2016; Yakovleva et al. 2016; Barklem 2016b, 2017; Belyaev & Voronov 2017).
Although the approximate model calculations are less timeconsuming than the full quantum calculations, they still require some computation. For this reason, the simplified model for inelastic atomic collisions with hydrogen has been proposed very recently (Belyaev & Yakovleva 2017); hereafter, we refer to this paper to as Paper I. The simplified model practically does not require any computational efforts. In the present paper, the simplified model is extended on inclusion of inelastic ionic collisions with hydrogen and applied to bariumhydrogen ionic collisions, which are of astrophysical interest (Asplund 2005).
2. Simplified model for estimating inelastic collision rate coefficients
2.1. General description
Let us consider a general case of collisions of A^{Z+} + H and A^{(Z+1)+} + H^{−}, where A is a chemical element of interest and Z = 0,1,... is a charge of this element. The case Z = 0 corresponds to atomic collisions with hydrogen, Z ≥ 1 to ionic collisions. The heavyparticle A^{Z+} can be in the ground and in excited states. The main goal of the present model is to estimate highvalued inelastic rate coefficients in collisions with hydrogen, which occur in nonLTE stellar atmosphere modeling.
We are interested in inelastic collision processes that change electronic states of atomic (ionic) particles, nevertheless, theoretical treatments of lowenergy atomic and ionic collisions are performed within the molecular representation, and reaction mechanisms are based on nonadiabatic transitions between molecular states. It has been proven by full quantum calculations (see Belyaev et al. 1999, 2010, 2012; Croft et al. 1999a,b; Belyaev & Barklem 2003) that the rate coefficients with high and moderate values are determined by longrange nonadiabatic transitions due to the interactions between A^{(Z+1)+} + H^{−} ionic and A^{Z +}(j) + H covalent molecular states. The index j denotes an electronic state of A^{Z+}. So, one needs to determine molecular state symmetries that lead to rates with high values. These symmetries are those which correspond to ionic molecular states. Within the LS coupling scheme, molecular states are determined by symmetries such as a total spin quantum number S and an absolute value of a projection quantum number of total orbital electronic angular momentum Λ and other quantum numbers. Since the H^{−} anion has the ^{1}S symmetry, the spin quantum number S of the ionic molecular state has a single value, while Λ can have several values, and both quantum numbers are uniquely determined by the ground state of the A^{(Z+1)+} ion.
As pointed out in Paper I, the rate coefficients K_{if} for an inelastic process i → f can be expressed via the reduced rate coefficients for neutralization and ionpair formation processes (1)or via for (de)excitation processes (2)where a statistical probability for population of the initial molecular state i with a given molecular symmetry Λ from a scattering channel and T is a temperature. The sums in these equations are taken over all possible quantum numbers Λ corresponding to the ground ionic molecular state A^{(Z+1)+} + H^{−}. The energy E_{j} is an electronic bound energy^{1} of an A^{Z +}(j) state to which a molecular state asymptotically corresponds. Although the spin S and other quantum numbers do not appear explicitly, these values are needed to select covalent molecular states that are populated during a collision and, hence, for selecting final atomic (ionic) A^{Z +}(f) states populated.
For estimating rate coefficients K_{if}(T), the reduced rate coefficients and must be known for each molecular symmetry of interest. The statistical probabilities can be readily calculated. It turns out that within some reasonable approximations (see below) the reduced rate coefficient for the neutralization process A^{(Z + 1) +} + H^{−} → A^{Z +}(f) + H depends on a single bound energy E_{f} ≡ E_{AZ +(f)}, while the reduced rate coefficient for the (de) excitation process A^{Z +}(i) + H → A^{Z +}(f) + H depends on two electronic bound energies, E_{i} ≡ E_{AZ +(i)} and E_{f} ≡ E_{AZ +(f)}.
Thus, in order to estimate a rate coefficient for a particular process with participation a particular chemical element, one must know the dependence of a corresponding reduced rate coefficient on bound energies. Such dependences are provided by the present simplified model. The reduced rate coefficients for the case Z = 0 are calculated and tabulated as functions of electronic bound energies in Paper I. The cases Z = 1,2 are treated in the present paper together with the general derivations and an application to bariumhydrogen ionic collisions.
The rate coefficients for exothermic processes (neutralization and deexcitation) depend weakly on the temperature, while the rates for inverse endothermic processes (ionpair formation and excitation) depend strongly on the temperature. Taking this into account, it is better first to calculate the rate coefficients K_{jk}(T) for exothermic processes (j → k transitions) and then to calculate the rate coefficients K_{kj}(T) for the inverse endothermic processes (k → j transitions) by means of the detailed balance relation (3)where ΔE_{jk} = E_{j}−E_{k} is the energy defect of the process. One should note that for exothermic (j → k) processes ΔE_{jk} is positive.
2.2. Reduced rate coefficients for neutralization and ionpair formation processes
Within a molecular symmetry Λ, a reduced rate coefficient for the neutralization process A^{(Z + 1) +} + H^{−} → A^{Z +}(f) + H is calculated by integrating a partial cross section for the process with the Maxwell distribution (4)where E is the collision energy, k_{B} the Boltzmann constant, and μ the reduced nuclear mass. In turn, the neutralization cross section is calculated as a sum over total angular momentum quantum number J(5)where is a probability for a statetostate transition from the initial state i = ionic to the final state f at a given collision energy E and a given quantum number J.
In order to estimate statetostate transition probabilities leading to rates with high values, one needs to know longrange ionic A^{(Z + 1) +} + H^{−} and covalent A^{Z +}(j) + H diabatic molecular potentials H_{ionic ionic} and H_{jj} as functions of the internuclear distance R. These are the Coulomb potential for an ionic molecular state (6)and the flat potentials^{2} for covalent molecular states (7)The electronic bound energy E_{H−} for a hydrogen anion is equal to −0.754 eV. These potential energy curves are plotted in Fig. 1. In a general case, the ionic potential crosses several covalent potentials creating a number of nonadiabatic regions; see Fig. 1a. A center R_{j} of each nonadiabatic region created by the ionic molecular state and the jth covalent molecular state is readily determined by equating potentials (6) and (7), and finally is written as (8)Thus, it is seen that a center of a nonadiabatic region is uniquely determined by a corresponding electronic bound energy.
Fig. 1 Schematic diagrams of a series of ionic (Coulomb) and covalent potentials in the diabatic representation. Panel a: depicts a general multichannel case, panel b: a twochannel case (neutralization), and panel c: a threechannel case (deexcitation). 

Open with DEXTER 
In a multichannel case, a statetostate transition probability can be calculated by the multichannel formula (Belyaev 1993), which requires information about an energy spectrum of A^{Z +}. The spectra are different for different chemical elements. However, the statetostate transition probability for neutralization can be evaluated by a probability calculated in the twochannel approximation depicted in Fig. 1b, that is, one ionic channel and only one covalent f channel. This approximation does not require information about the whole spectrum but only about the energy level f. In this case the probability is expressed via a nonadiabatic transition probability p_{f}, after a singlepass traversal through the nonadiabatic region f formed by the ionic state and the final state f, as follows: (9)The nonadiabatic transition probability p_{j} can be calculated within the LandauZener model as follows: (10)where H_{ionic j} is an offdiagonal matrix element for the (diabatic) ioniccovalent interaction, v is a radial velocity of colliding heavy particles, and primed values are referred to as a derivative with respect to the internuclear distance R. All values are evaluated at the center of the nonadiabatic region R_{j}; see Eq. (8). Substituting the ionic diabatic potential (6) and the covalent diabatic potential (7), and the expression (8) for R_{j} into Eq. (10), one can rewrite the formula for the LandauZener probability (11)The offdiagonal matrix element H_{ionic j} can be evaluated by the semiempirical formula for oneelectron transition charge exchange from Olson et al. (1971), i.e., (12)For the neutralization process the initial state is ionic and the final state f, so the index j in Eqs. (7), (8), (10)–(12) should be replaced by the index f. The radial velocity in the nonadiabatic region formed by the ionic state and the final state f reads (13)It is seen that the statetostate transition probability , Eq. (9), for a neutralization process at a given collision energy E and at given quantum number J is uniquely determined by a single electronic bound energy E_{f}. This means that the reduced rate coefficient for neutralization at a given temperature T also depends on a single electronic bound energy E_{f}.
If the reduced rate coefficient for the exothermic neutralization process (j → k transitions; for neutralization processes j = ionic is known, then the rate coefficient for the inverse endothermic ionpair formation process (k → j transitions) can be calculated by means of the following detailed balance relation: (14)This relation is readily derived from Eqs. (1)–(3).
Fig. 2 General dependence of neutralization reduced rate coefficients on a final electronic bound energy E_{f} obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER 
Thus, Eqs. (4)–(14) allow us to calculate the reduced rate coefficients for neutralization and ionpair formation processes as a function of an electronic bound energy E_{f} at a given temperature T. Within the simplified model the reduced rates do not depend on Λ, but this quantum number selects which final states are populated from the ionic state within this molecular symmetry. The neutralization reduced rate coefficient for the case Z = 0 is calculated in Paper I and for the cases Z = 1,2 in the present paper. The results are published at the CDS for the temperature range T = 1000−10 000 K. For the temperature T = 6000 K the neutralization reduced rate coefficients N_{if}(T = 6000;E_{f}) are shown in Fig. 2.
It is seen in this figure that the neutralization reduced rate coefficient for Z = 0 (dashed curve) reaches the maximum value 6.8 × 10^{8} cm^{3}/ s at the bound energy E_{f} = −2 eV. Also, within the most optimal window^{3} [−2.7, −1.5] eV the reduced rates have high values above 10^{8} cm^{3}/ s, while within the less optimal window (roughly [−4.0, −1.3] eV) the reduced rates have moderate values above 10^{10} cm^{3}/ s, down to nearly three orders of magnitude lower than the maximum value (see Paper I). Outside of these optimal windows the reduced rates have negligible values.
For the case of ionic collisions the behavior of the reduced rate coefficients is different. For Z = 1 the neutralization reduced rate coefficient (solid line) has the broad maximum of 7.6 × 10^{8} cm^{3}/ s at the bound energy E_{f} = −4.4 eV. The most optimal window in these collisions can be estimated by the interval [−2.8, −8.3] eV for the bound energy, while the less optimal window corresponds to the bound energies <−2.1 eV. For the case of Z = 2, the neutralization reduced rate (dashdotted curve) has even broader maximum of 7.7 × 10^{8} cm^{3}/ s at the bound energy E_{f} = −8.2 eV. The most optimal window in these collisions corresponds to the electronic bound energies <−4.5 eV, while the less optimal window <−3 eV. Thus, one may conclude that the optimal windows for neutralization are shifted in the direction toward deeper bound energies with increasing the charge of positive ion in collisions with hydrogen negative ions.
2.3. Reduced rate coefficients for deexcitation and excitation processes
Within a molecular symmetry Λ, a reduced rate coefficient for the deexcitation process A^{Z +}(i) + H → A^{Z +}(f) + H is calculated by integrating a deexcitation cross section with the Maxwell distribution (15)The deexcitation cross section is calculated as a sum over total angular momentum quantum number J,(16)where is a probability for a transition from an initial state i to a final state f.
Within the same approximation as that used in the previous subsection, the deexcitation reaction mechanism is determined by minimum three diabatic molecular states: two covalent states i and f, and one ionic state, as depicted in Fig. 1c. Therefore, there are two nonadiabatic regions that are responsible for the statetostate transition probability; that is, is expressed via two nonadiabatic transition probabilities, p_{i} and p_{f}, as follows: (17)Each probability p_{i} and p_{f} is after a singlepass traversal through a nonadiabatic region formed by the ionic state and the initial or final state i or f, and can be calculated within the LandauZener model by means of Eqs. (11), (12); see also Eq. (8). In these formulas the index j should be replaced by indices i or f for probabilities p_{i} or p_{f}, respectively. The radial velocities for deexcitation in the nonadiabatic regions i and f read, respectively, Finally, one can see that the statetostate transition probability and the reduced rate coefficient for deexcitation depend on two electronic bound energies E_{i} and E_{f}.
Thus, Eqs. (8), (11), (12), and (15)–(19) allow us to calculate the reduced rate coefficients for deexcitation processes at a given temperature T as a function of an electronic bound energies E_{i} and E_{f}. The reduced rates for excitation processes can be calculated from by means of the detailed balance relation similar to Eq. (14) replacing N by D.
The (de)excitation reduced rate coefficient for the case Z = 0 is calculated and discussed in Paper I and for the cases Z = 1,2 in the present paper. The results are published at the CDS for the temperature range T = 1000−10 000 K. For Z = 1 and for the temperature T = 6000 K the (de)excitation reduced rate coefficient D_{if}(T = 6000;E_{i},E_{f}) is shown in Fig. 3.
Fig. 3 Excitation and deexcitation reduced rate coefficient for Z = 1 as a function of the electronic bound energies E_{i} and E_{f} for the temperature T = 6000 K. The bound energies are measured from the ionization limit. These two subfigures present the same surface for the reduced rate D_{if} shown from different view points. 

Open with DEXTER 
The (de)excitation reduced rate coefficient for Z = 1 reaches the maximum value 1.1 × 10^{8} cm^{3}/ s at the bound energies E_{i} ≈ E_{f} ≈ −4 eV. The most optimal window (exceeding the value 10^{8} cm^{3}/ s) is located in a narrow vicinity around this point. The less optimal window (exceeding the value 10^{10} cm^{3}/ s) is roughly [− 8,−2.2] eV for E_{i} and <−2.4 eV for E_{f}.
2.4. Summary of the simplified model
In summary, to estimate the rate coefficients of inelastic processes in collisions with hydrogen by means of the simplified method proposed in the present paper, the following simple steps should be performed:
Step 1. Determining a molecular symmetry for the ground ionic (diabatic) A^{(Z + 1) +} + H^{−} molecular state within the LS scheme, i.e., determining a total electronic spin S, an electronic orbital projection quantum number Λ, and other quantum numbers, if necessary. The values S and Λ for the ground ionic molecular state are uniquely determined by the ground state of a cation A^{(Z+1)+}, since an anion H has both spin and orbital quantum numbers equal to zero; S has a single value, Λ might have several values.
Step 2. Determining A^{Z +}(j) states, which produce covalent A^{Z +}(j) + H molecular states of the same symmetry as the ground ionic state determined in Step 1, and which correspond to oneelectron transitions for the ioniccovalent interaction, i.e., between the A^{(Z + 1) +} + H^{−} and A^{Z +}(j) + H (diabatic) states. The electronic bound energies E_{j} for the selected atomic/ionic states A^{Z +}(j) must be determined, for example, from the NIST database.
Step 3. This step is slightly different for the neutralization and deexcitation processes.
Step 3a for neutralization processes A^{(Z + 1) +} + H^{−} → A^{Z +}(f) + H. For a particular final state A^{Z +}(f), the corresponding electronic bound energy E_{f} is taken and then used to determine the corresponding reduced rate coefficient N_{if} from the tabulated function N_{if}(E_{f}) presented in Fig. 2 for T = 6000 K and provided in the supplementary material for T = 1000−10 000 K; in other words, N_{if} = N_{if}(E_{f}) is determined from the calculated function N_{if}(T;E_{f}) for a given temperature T.
Step 3b for deexcitation processes A^{Z +}(i) + H → A^{Z +}(f) + H. For a particular initial state A^{Z +}(i) and a particular final state A^{Z +}(f), the corresponding electronic bound energies E_{i} and E_{f} are taken and then used to determine the corresponding reduced rate coefficient D_{if} from the tabulated function D_{if}(E_{i},E_{f}) presented in Fig. 3 as the twodimensional (2D) surface for T = 6000 K and provided in the supplementary material for T = 1000−10 000 K; in other words, D_{if} = D_{if}(E_{i},E_{f}) is determined from the calculated function D_{if}(T;E_{i},E_{f}) for a given temperature T. This can be performed in one step by using the 2D function D_{if}(E_{i},E_{f}) or in two steps by first taking a cut of the 2D function D_{if}(E_{i},E_{f}) at a given bound energy E_{i} for the initial state A^{Z +}(i) and then by treating the 1D function D_{if}(E_{f}) at the given bound energy E_{f};
Step 4. In order to get the rate coefficient for the neutralization or deexcitation processes, a statistical probability for the population of the initial molecular state is calculated and the reduced rate is multiplied by this statistical probability according to Eqs. (1) or (2), respectively. If the ground ionic A^{(Z + 1) +} + H^{−} molecular state corresponds to various molecular symmetries (several Λ), then a sum over Λ should be taken. In order to get the rate coefficient for ionpair formation or for excitation processes, the rate coefficient is calculated for an inverse process from the corresponding rate coefficient for the direct process according to the detailed balance relation (3).
3. Application to bariumhydrogen ionic collisions
In this section, the simplified model is applied for estimating rate coefficients of neutralization, ionpair formation, and excitation and deexcitation processes in bariumhydrogen ionic collisions (Z = 1). To the best of our knowledge no bariumhydrogen collision rates have been published.
Step 1: determining a ground ionic Ba^{2 +} + H^{−} (diabatic) molecular state symmetry. For bariumhydrogen collisions, the ground ionic Ba^{2 +}(5p^{6}^{1}S) + H^{−}(1s^{2 1}S) molecular state has the only symmetry, BaH. In this case, each sum in Eqs. (1) and (2) only has the term with Λ = 0 and total spin S = 0.
Step 2: determining the barium ionic states to be treated and the corresponding electronic bound energies. It is necessary to determine bariumion states, which (i) create BaH molecular states and (ii) correspond to singleelectron transitions to the ground ionic molecular state. Only the ionic states with electronic bound energies that have values E_{j} < −1.3 eV are taken into account; above this energy the rate coefficients are negligibly small along with states that correspond to multielectron transitions into the ground ionic molecular state. These barium ionic states with the corresponding scattering channels, electronic bound energies, and statistical probabilities are collected in Table 1.
Step 3: determining reduced rate coefficients. Step 3a. In order to estimate rate coefficients for a neutralization process, a reduced rate N_{if}(T;E_{f}) should be taken at a given temperature as a function of a finalstate electronic bound energy E_{f}. At T = 6000 K this function is plotted in Fig. 2 and shown again in Fig. 4 by solid lines. The reduced rate coefficients for particular neutralization processes in Ba^{2+} + H^{−} collisions are determined by the values of the function N_{if}(T;E_{f}) at the particular final electronic bound energies E_{f} from Table 1; see symbols in Fig. 4. For example, the final state Ba^{+}(6d ^{2}D) has the electronic bound energy E_{f} = −4.29573 eV.
Scattering channels correlated to BaH^{+}() molecular states, their statistical probabilities and the asymptotic energies equal to the electronic bound energies E_{j} ( averaged experimental values taken from NIST, Kramida et al. 2016) with respect to the ionization limit BaS) + H(1s ^{2}S).
Fig. 4 General dependence of neutralization reduced rate coefficient N_{if} (solid line) on a final bound energy E_{f} and the particular reduced rate coefficients for Ba^{2+} + H^{−} → Ba^{+}(f) + H processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER 
The corresponding reduced rate coefficient for the neutralization process Ba^{2 +} + H^{−} → Ba^{+}(6d^{2}D) + H is equal to N_{if}(T = 6000 K,E_{f} = −4.29573 eV) = 7.59 × 10^{8} cm^{3}/ s (see Fig. 4). This rate coefficient has the highest value in these collisions.
Step 3b. In order to estimate rate coefficients for deexcitation and excitation processes, a reduced rate D_{if}(T;E_{i},E_{f}) should be taken at a given temperature as a 2D function of the initial and finalstate electronic bound energies. At T = 6000 K this function is plotted in Fig. 3. The reduced rate coefficients for particular excitation and deexcitation processes Ba^{+}(i) + H → Ba^{+}(f) + H at a given temperature can first be determined as a cut of the 2D surface D_{if}(T;E_{i},E_{f}) at an initial bound energy E_{i}, and then by taking particular values from this cut of D_{if} at the final bound energies E_{f} (see Table 1). Examples of the excitation and deexcitation processes in bariumhydrogen ionic collisions are shown in Figs. 5–7 for the initial states i = 5−7.
Fig. 5 Cut of deexcitation/excitation reduced rate coefficient D_{if}(E_{i},E_{f}) for particular deexcitation and excitation processes Ba^{+}(6d) + H → Ba^{+}(f) + H for the initial bound energy E_{i} = E_{5} as a function of a final bound energy E_{f} (solid line) and the particular reduced rate coefficients for the same processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER 
Fig. 6 Cut of deexcitation/excitation reduced rate coefficient D_{if}(E_{i},E_{f}) for particular deexcitation and excitation processes Ba^{+}(4f) + H → Ba^{+}(f) + H for the initial bound energy E_{i} = E_{6} as a function of a final bound energy E_{f} (solid line) and the particular reduced rate coefficients for the same processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER 
Fig. 7 Cut of deexcitation/excitation reduced rate coefficient D_{if}(E_{i},E_{f}) for particular deexcitation and excitation processes Ba^{+}(7p) + H → Ba^{+}(f) + H for the initial bound energy E_{i} = E_{7} as a function of a final bound energy E_{f} (solid line) and the particular reduced rate coefficients for the same processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER 
Step 4: estimating rate coefficients. Since the statistical probability for the ionic channel is presently equal to 1, the rate coefficients of the neutralization processes K_{if}(T) coincide with the corresponding reduced rate coefficients N_{if}, see Eq. (1). Rate coefficients for the ionpair formation processes are calculated from the neutralization rates by the detailed balance Eq. (3). The rate coefficients K_{if}(T) for the deexcitation processes are calculated from the reduced rate coefficients D_{if} by multiplying the latter by the statistical probability ; see Eq. (2). Although rate coefficients K_{if}(T) for excitation processes can be calculated directly from the reduced rate coefficients D_{if}, it is more practical first to calculate the rate coefficients for the deexcitation processes and then to compute K_{if}(T) for excitation processes from the rate coefficients for the deexcitation processes using the detailed balance Eq. (3).
The rate coefficients K_{if}(T) for all possible processes in bariumhydrogen ionic collisions are estimated by the described simplified model and presented as the supplementary material for the temperature range T = 1000−10 000 K. The graphical representation of the calculated rate coefficients in bariumhydrogen ionic collisions at T = 6000 K is presented in Fig. 8. The highvalued rate coefficients K_{if}(T) are shown in Table 2 for the temperatures T = 2000, 6000, 10 000 K. It is seen that the highest values of the rate coefficients correspond to the neutralization processes into the Ba^{+}(6p)−Ba^{+}(8s) states and the moderate values correspond to the neutralization and the (de)excitation processes from and to these states and the states around them; see Fig. 8 and Table 2, in agreement with the general rules discussed above.
Fig. 8 Graphical representation of the rate coefficients in bariumhydrogen ionic collisions obtained by means of the simplified model at T = 6000 K. 

Open with DEXTER 
Rate coefficients, in units of cm^{3}/s, for selected temperatures in the range T = 2000−10 000 K for the excitation, deexcitation, ionpair formation, and mutual neutralization processes in bariumhydrogen ionic collisions.
4. Conclusions
The simplified model is derived and described in detail in the present paper to estimate rate coefficients of inelastic processes in lowenergy collisions of hydrogen atoms and negative ions with atoms and positive ions of other chemical elements; i.e., A^{(Z + 1) +} + H^{−} and A^{Z +}(j) + H collisions, Z = 0,1,... The main goal of the model is to estimate rate coefficients with high and moderate values; that is, the rates that are important for nonLTE stellar atmosphere modeling. The model is based on the main reaction mechanism of the processes under consideration, which is determined by nonadiabatic transitions in longrange regions created by ioniccovalent molecularstate interactions.
Within the simplified model, the rate coefficients are expressed via reduced rate coefficients that are different for neutralization and ionpair formation processes A^{(Z + 1) +} + H^{−} ⇌ A^{Z +}(f) + H and for excitation/deexcitation processes A^{Z +}(i) + H → A^{Z +}(f) + H. It is shown that the reduced rate coefficients for neutralization and ionpair formation processes depend on a single electronic bound energy of an atomic/ionic state A^{Z +}(f), which correlates to a corresponding covalent molecular state, while the reduced rate coefficients for excitation/deexcitation processes depend on two bound energies: electronic bound energies for an initial and a final atomic/ionic states. The reduced rate coefficients are calculated and tabulated for the case Z = 0 in Paper I, while these coefficients are calculated for cases Z = 1,2 in the present paper. Optimal windows for neutralization and ionpair formation and for excitation/deexcitation processes are determined as well.
The simplified method does not have any adjustable parameters. Moreover, one does not need to perform any complicated calculations; all calculations of reduced rate coefficients have been accomplished and tabulated in Paper I and in the present paper, and finally the rate coefficients can be estimated by the simple Eqs. (1)–(3). Without any additional calculation, the method indicates which processes have high and moderate rate coefficients, and hence are expected to be important for nonLTE stellar atmosphere modeling.
Furthermore, the simplified method helps in understanding a general distribution of inelastic rate coefficients in collisions of various chemical elements with hydrogen. In particular, the method selects atomic/ionic states that are efficiently involved in the inelastic Hcollision processes. This is carried out by determining the optimal windows of atomic/ionic states; that is, by determining the electronic bound energy intervals for atomic/ionic states whose participation leads to high and moderate rate coefficients. It is shown that these intervals are the same for all chemical elements. For Z = 0, they are located in the vicinity of the electronic bound energy −2 eV, for Z = 1 around −4.4 eV, and for Z = 2 around −8.2 eV.
The simplified method allows rate coefficients to be estimated for neutralization, ionpair formation, excitation, and deexcitation processes in collisions of hydrogen with many chemical elements. In the present paper, the simplified model is applied
to bariumhydrogen ionic collisions, Z = 1. The rate coefficients for the most important neutralization, ionpair formation, excitation, and deexcitation processes are estimated and presented at the CDS.
An electronic bound energy E_{j} is measured from an ionization limit of A^{Z+}, and, hence, it is negative. If an electronic energy is measured from a groundstate energy of A^{Z+} as a zeroenergy level, then the electronic bound energy should be recalculated as , where I_{ionization} is an ionization potential of A^{Z+}.
By “an optimal window” we mean an electronic energy interval for atomic or ionic states whose participation in inelastic collision processes leads to relatively high rate coefficients. We also define “the most optimal window” as rate coefficient values exceeding 10^{8} cm^{3}/ s and “the least optimal window” as rate coefficient values exceeding 10^{10} cm^{3}/ s.
Acknowledgments
The authors gratefully acknowledge partial support from the Ministry for Education and Science (Russian Federation), projects # 3.1738.2017/4.6, 3.5042.2017/6.7.
References
 Asplund, M. 2005, ARA&A, 43, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Barklem, P. S. 2016a, A&ARv, 24, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Barklem, P. S. 2016b, Phys. Rev. A, 93, 042705 [NASA ADS] [CrossRef] [Google Scholar]
 Barklem, P. S. 2017, Phys. Rev. A, 95, 069906(E) [NASA ADS] [CrossRef] [Google Scholar]
 Barklem, P. S., Belyaev, A. K., & Asplund, M. 2003, A&A, 409, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barklem, P. S., Belyaev, A. K., Guitou, M., et al. 2011, A&A, 530, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belyaev, A. K. 1993, Phys. Rev. A, 48, 4299 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Belyaev, A. K. 2013a, Phys. Rev. A, 88, 052704 [NASA ADS] [CrossRef] [Google Scholar]
 Belyaev, A. K. 2013b, A&A, 560, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belyaev, A. K. 2015, Phys. Rev. A, 91, 062709 [NASA ADS] [CrossRef] [Google Scholar]
 Belyaev, A. K., & Barklem, P. S. 2003, Phys. Rev. A, 68, 062703 [NASA ADS] [CrossRef] [Google Scholar]
 Belyaev, A. K., & Voronov, Y. V. 2017, A&A, 606, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belyaev, A. K., & Yakovleva, S. A. 2017, A&A, 606, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belyaev, A. K., Grosser, J., Hahne, J., & Menzel, T. 1999, Phys. Rev. A, 60, 2151 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Belyaev, A. K., Barklem, P. S., Dickinson, A. S., & Gadéa, F. X. 2010, Phys. Rev. A, 81, 032706 [NASA ADS] [CrossRef] [Google Scholar]
 Belyaev, A. K., Barklem, P., Spielfiedel, A., et al. 2012, Phys. Rev. A, 85, 032704 [NASA ADS] [CrossRef] [Google Scholar]
 Belyaev, A. K., Yakovleva, S. A., & Barklem, P. S. 2014, A&A, 572, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belyaev, A. K., Yakovleva, S., Guitou, M., et al. 2016, A&A, 587, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Croft, H., Dickinson, A. S., & Gadéa, F. X. 1999a, J. Phys. B: At. Mol. Opt. Phys., 32, 81 [NASA ADS] [CrossRef] [Google Scholar]
 Croft, H., Dickinson, A. S., & Gadéa, F. X. 1999b, MNRAS, 304, 327 [NASA ADS] [CrossRef] [Google Scholar]
 Drawin, H.W. 1968, Z. Phys., 211, 404 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Drawin, H.W. 1969, Z. Phys., 225, 483 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Drawin, H.W., & Emard, F. 1973, Phys. Lett. A, 43, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Guitou, M., Belyaev, A. K., Barklem, P. S., Spielfiedel, A., & Feautrier, N. 2011, J. Phys. B, 44, 035202 [NASA ADS] [CrossRef] [Google Scholar]
 Guitou, M., Spielfiedel, A., Rodionov, D. S., et al. 2015, Chem. Phys., 462, 94 [NASA ADS] [CrossRef] [Google Scholar]
 Kramida, A., Ralchenko, Y., Reader, J., & NIST ASD Team 2016, NIST Atomic Spectra Database (Version 5.4), http://physics.nist.gov/asd [Google Scholar]
 Lambert, D. L. 1993, Phys. Scr. T, 47, 186 [NASA ADS] [CrossRef] [Google Scholar]
 Lind, K., Asplund, M., Barklem, P., & Belyaev, A. 2011, A&A, 528, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mashonkina, L. 2014, IAU Symp., 298, 355 [NASA ADS] [Google Scholar]
 Mashonkina, L., Sitnova, T., & Belyaev, A. K. 2017, A&A, 605, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mitrushchenkov, A., Guitou, M., Belyaev, A. K., et al. 2017, J. Chem. Phys., 146, 014304 [NASA ADS] [CrossRef] [Google Scholar]
 Olson, R. E., Smith, F. T., & Bauer, E. 1971, Appl. Opt., 10, 1848 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Osorio, Y., & Barklem, P. S. 2016, A&A, 586, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Steenbock, W., & Holweger, H. 1984, A&A, 130, 319 [NASA ADS] [Google Scholar]
 Yakovleva, S., Voronov, Y. V., & Belyaev, A. K. 2016, A&A, 593, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Scattering channels correlated to BaH^{+}() molecular states, their statistical probabilities and the asymptotic energies equal to the electronic bound energies E_{j} ( averaged experimental values taken from NIST, Kramida et al. 2016) with respect to the ionization limit BaS) + H(1s ^{2}S).
Rate coefficients, in units of cm^{3}/s, for selected temperatures in the range T = 2000−10 000 K for the excitation, deexcitation, ionpair formation, and mutual neutralization processes in bariumhydrogen ionic collisions.
All Figures
Fig. 1 Schematic diagrams of a series of ionic (Coulomb) and covalent potentials in the diabatic representation. Panel a: depicts a general multichannel case, panel b: a twochannel case (neutralization), and panel c: a threechannel case (deexcitation). 

Open with DEXTER  
In the text 
Fig. 2 General dependence of neutralization reduced rate coefficients on a final electronic bound energy E_{f} obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER  
In the text 
Fig. 3 Excitation and deexcitation reduced rate coefficient for Z = 1 as a function of the electronic bound energies E_{i} and E_{f} for the temperature T = 6000 K. The bound energies are measured from the ionization limit. These two subfigures present the same surface for the reduced rate D_{if} shown from different view points. 

Open with DEXTER  
In the text 
Fig. 4 General dependence of neutralization reduced rate coefficient N_{if} (solid line) on a final bound energy E_{f} and the particular reduced rate coefficients for Ba^{2+} + H^{−} → Ba^{+}(f) + H processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER  
In the text 
Fig. 5 Cut of deexcitation/excitation reduced rate coefficient D_{if}(E_{i},E_{f}) for particular deexcitation and excitation processes Ba^{+}(6d) + H → Ba^{+}(f) + H for the initial bound energy E_{i} = E_{5} as a function of a final bound energy E_{f} (solid line) and the particular reduced rate coefficients for the same processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER  
In the text 
Fig. 6 Cut of deexcitation/excitation reduced rate coefficient D_{if}(E_{i},E_{f}) for particular deexcitation and excitation processes Ba^{+}(4f) + H → Ba^{+}(f) + H for the initial bound energy E_{i} = E_{6} as a function of a final bound energy E_{f} (solid line) and the particular reduced rate coefficients for the same processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER  
In the text 
Fig. 7 Cut of deexcitation/excitation reduced rate coefficient D_{if}(E_{i},E_{f}) for particular deexcitation and excitation processes Ba^{+}(7p) + H → Ba^{+}(f) + H for the initial bound energy E_{i} = E_{7} as a function of a final bound energy E_{f} (solid line) and the particular reduced rate coefficients for the same processes (symbols) obtained by means of the present simplified model at T = 6000 K. 

Open with DEXTER  
In the text 
Fig. 8 Graphical representation of the rate coefficients in bariumhydrogen ionic collisions obtained by means of the simplified model at T = 6000 K. 

Open with DEXTER  
In the text 