Issue 
A&A
Volume 530, June 2011



Article Number  A94  
Number of page(s)  9  
Section  Stellar atmospheres  
DOI  https://doi.org/10.1051/00046361/201116745  
Published online  18 May 2011 
On inelastic hydrogen atom collisions in stellar atmospheres
^{1}
Department of Physics and AstronomyUppsala University, Box 515, 75120 Uppsala, Sweden
email: Paul.Barklem@physics.uu.se
^{2}
Department of Theoretical Physics, Herzen University, 191186 St. Petersburg, Russia
^{3}
Université ParisEst, Laboratoire Modélisation et Simulation MultiÉchelle, MSME UMR 8208 CNRS, 5 Bd Descartes, 77454 MarnelaVallée, France
^{4}
LERMA and UMR 8112 of CNRS, Observatoire de ParisMeudon, 92195 Meudon Cedex, France
^{5}
Laboratoire de Chimie et Physique Quantique, UMR 5626 du CNRS, IRSAMC, Université Paul Sabatier, 118 Rte de Narbonne, 31062 Toulouse, France
Received: 18 February 2011
Accepted: 10 April 2011
The influence of inelastic hydrogen atom collisions on nonLTE spectral line formation has been, and remains to be, a significant source of uncertainty for stellar abundance analyses, due to the difficulty in obtaining accurate data for lowenergy atomic collisions either experimentally or theoretically. For lack of a better alternative, the classical “Drawin formula” is often used. Over recent decades, our understanding of these collisions has improved markedly, predominantly through a number of detailed quantum mechanical calculations. In this paper, the Drawin formula is compared with the quantum mechanical calculations both in terms of the underlying physics and the resulting rate coefficients. It is shown that the Drawin formula does not contain the essential physics behind direct excitation by H atom collisions, the important physical mechanism being quantum mechanical in character. Quantitatively, the Drawin formula compares poorly with the results of the available quantum mechanical calculations, usually significantly overestimating the collision rates by amounts that vary markedly between transitions.
Key words: atomic data / line: formation / stars: abundances / stars: atmospheres
© ESO, 2011
1. Introduction
Solution of the nonlocal thermodynamic equilibrium (nonLTE) radiative transfer problem in stellar atmospheres requires detailed and complete knowledge of the radiative and collisional processes that affect the statistical equilibrium of a given atomic species of interest. A difficulty regarding the collisional processes is to determine which, among the almost endless possibilities in a stellar atmosphere, are important. In an early study of the formation of the Na D lines in the solar spectrum, Plaskett (1955) considered the two obvious candidates for the case of the solar atmosphere: inelastic collisions with electrons and with hydrogen atoms in their ground state. The importance of electron collisions in many environments is well known, arising from the fact that electrons are always the most abundant charged particle and have a much higher thermal velocity than atoms, leading directly to a higher collision rate.
Hydrogen atoms are the most abundant perturber by far; however, as pointed out by Plaskett, given that in thermal conditions atoms are much slower than electrons, simple considerations based on the Massey criterion (Massey 1949) indicate H atom collisions are nearly adiabatic and thus should have much smaller inelastic cross sections than electron collisions at the relevant thermal velocities (see Lambert 1993, for a detailed discussion). Plaskett considered experimental estimates for quenching of Na by He and methane, which suggested deexcitation cross sections ranging from nondetectable for He to a small value for methane. Assuming that quenching by H might be similar to that for methane, Plaskett estimated the inelastic H atom collision rate would be of the same order of magnitude as that for electrons, despite n(H)/n(e) ~ 10^{4}. As indicated by Plaskett, if one made the more reasonable assumption that H would have a similar cross section to He, then the rate would probably be negligible.
Based on this rough expectation, and probably also because of the almost total lack of any data for atomatom collisions at low energy due to the increased difficulties in calculating and measuring such cross sections compared to the case of electrons (e.g. Massey 1955), most early nonLTE studies of cool stars during the 1960’s and 1970’s neglected H atom collisions. While some authors mentioned the possible need to include such processes (e.g. van Regemorter 1962; Athay & Canfield 1969), the first to make an attempt at including them quantitatively was the study of solar Na lines by Gehren (1975). That study used the classical formula from Drawin (1969), though details on how these formula were extended from the case of collisions between like atoms to collisions between unlike atoms are not given. Gehren found the H atom collisions to be relatively unimportant, less than half of the electron collision rates.
A later study of Li line formation by Steenbock & Holweger (1984) used a similar approach, though in this case the generalisation of the formula of Drawin to inelastic H atom collisions was now given. The resulting formula is often referred to in stellar astrophysics as the “Drawin formula”. The key development of this study was that it considered the case of a lowmetallicity halo dwarf star, where n(H)/n(e) ~ 10^{6}. They found that if H atom collisions were described to orderofmagnitude accuracy by this formula, then indeed the large ratio n(H)/n(e) would be sufficient to offset the lower efficiency and frequency of H atom collisions relative to electrons, and inelastic collisions with H atoms would be the dominant collisional process in cool metalpoor dwarfs. Thus, the study of Steenbock and Holweger restimulated the idea that H atom collisions might be important in cool stellar atmospheres, and since that study the question has been: does the Drawin formula really provide orderofmagnitude estimates of these processes?
Recent theoretical and experimental studies have put us in a position to answer this question. In this paper we will examine the question from two points of view. First, is there any reasonable expectation that the Drawin formula can provide reasonable results given our present understanding of the physics involved? This is discussed in Sect. 2. Second, how do the results from the Drawin formula compare with those from quantum mechanical calculations? This comparison is made in Sect. 3. Finally, the main conclusions are summarised.
2. Physics of inelastic hydrogen collisions
In this section we discuss the physics behind the classical Drawin formula. We then compare with quantum theory and the physics of inelastic H atom collisions as revealed by experiment and detailed quantum scattering calculations based on quantumchemistry descriptions of the relevant quasimolecules.
2.1. The classical Drawin formula
The Drawin formula is the result of a number of modifications and extensions of the classical formula for ionization of atoms by electron impact due to Thomson (1912). In Thomson’s theory, the bound electron in the target atom is considered as a stationary free classical electron. The Coulomb interaction of an incident classical electron with the atomic electron is then considered and Thomson calculates the impact parameter for which the deflection of the incident electron results in an energy transfer between electrons corresponding precisely to the ionization potential. All collisions inside this impact parameter will have a larger energy transfer and thus lead to ionization and so the cross section can be easily calculated. Elwert (1952) introduced a factor to account for equivalent electrons in the outer shell of the atom, which is usually also applied.
As pointed out by Rudge (1968, Sect. 4), Thomson’s theory is remarkable for being able to make a number of qualitative predictions regarding ionization cross sections and their behaviour with impact energy. However, quantitatively it provides only an orderofmagnitude estimate. For example, Seaton (1962) compared with experimental results for H, He, Ne and Ar, and found the classical theory overestimates by a factor of about 5 at low energy. Seaton explained this from the fact that in quantum mechanical treatments the ionization probability rarely approaches unity even at small impact parameters, while the classical theory assumes the probability is unity for all impact parameters smaller than that leading to a classical energy transfer equal to the ionization potential. It was also observed that the comparison differed markedly for Ne, a difference that is explained by quantum mechanical properties of the atom, and cannot be explained by classical theory.
In Fig. 1 we complement Seaton’s comparisons with new comparisons with modern experimental data from the Queen’s University of Belfast laboratory for the cases of electron impact ionization of ground state hydrogen, oxygen, magnesium and iron (Shah et al. 1987; Thompson et al. 1995; McCallion et al. 1992; Shah et al. 1993, respectively). These experimental data are in good agreement with other experiments, cover a wide range of collision energies including those near threshold and the absolute uncertainty is expected to be less than 20% (see caption to Fig. 1). The general features of the comparison are very similar to those seen in Seaton’s comparison, i.e. the classical formula overestimates the cross sections by around a factor of 5 at low energy. In this comparison, O is seen to be particularly discrepant, differing by a factor of 20 near threshold. Thus, quantitively, the Thomson formula is able to provide only a rough orderofmagnitude estimate of the ionization cross section. Given the difficulties involved in quantum calculations compared to the classical approach, considerable efforts were made to improve the classical models, predominantly through inclusion of the initial velocity of the atomic electron, the effects of the target nucleus, and exchange effects between electrons. These efforts were met with mixed success, and there is an extensive literature on the subject (see Burgess & Percival 1968; Rudge 1968; Bely & van Regemorter 1970; Bates 1978, for some reviews).
Fig. 1 Comparison of the classical Thomson (1912) cross section for ionization of H, O, Mg and Fe atoms by electron impact with experimental results of Shah et al. (1987); Thompson et al. (1995); McCallion et al. (1992); Shah et al. (1993), respectively. The plot shows the ratio of the experimental values with respect to the theoretical cross section, plotted against the collision energy E scaled by the ionization energy I. The errors in the experimental values are in all cases typically 5% or less, with an additional uncertainty of around 7–14% in the absolute scale. 
Drawin (1968, 1969); Drawin & Emard (1973) extended the Thomson theory to the case of ionization and excitation in collisions between like atoms, which was later generalised by Steenbock & Holweger (1984) to the case of H atom collisions with other atoms. The derivation has been redone in considerable detail by Lambert (1993), correcting an error in the cross section expressions due to confusion of centreofmass and laboratory frames; an error that was pointed out earlier by Fleischmann & Dehmel (1972). It is worth noting that if this error was carried through one would expect considerable errors in the rate coefficients due to the incorrect threshold. However, the equations for the rate coefficients generalised by Steenbock & Holweger (1984) seem to be correct, with the final expression for the rate coefficient derived by Lambert (1993) differing from that of Steenbock and Holweger by an essentially negligible factor m_{A}/(m_{A} + m_{H}), of order unity. The resulting expression for the rate coefficient for an inelastic H atom collision with atom A, A_{l} + H → A_{u} + H, where l and u denote lower and upper states respectively, as presented by Lambert is: (1)where f_{lu} is the oscillator strength for the transition l → u, is the ionization potential of H, , , and Ψ is a function arising from the integration over the Maxwell velocity distribution, which for is: (2)All other quantities have their usual meaning.
The main physical assumption in the generalisation of the Thomson theory to ionization in collisions between like atoms, is that for a given collision energy above threshold, the efficiency of the energy transfer in an atomatom collision is the same as in the electronatom collision. This amounts to neglect of the nucleus of the perturbing atom, and thus the model in essence is the classical interaction of the impact atom electron with the stationary free classical atomic electron, allowing for the increased mass of the perturber. Drawin introduces an additional mass factor (m_{A}/m_{H}) without justification, which has been discussed by Fleischmann & Dehmel (1972) and Lambert (1993). While Drawin compared with experiments for ionization in collisions between light atoms and between molecules and found agreement within a factor of two, Fleischmann & Dehmel made comparisons and found only orderofmagnitude agreement becoming worse for heavier atoms. In the case of collisions between atoms of similar mass, Fleischmann & Dehmel found only agreement within a factor of 100. Comparisons with more recent experimental data do not change this conclusion (Kunc & Soon 1991).
Classical descriptions of atomic processes are less well suited to excitation than to ionization, due to the generally smaller energy transfer and the difficulty in choosing a classical final energy state band. In Drawin (1969; see also Drawin & Emard 1973) the ionization formula is extended to the case of excitation by analogy with the ionization formula; that is, the ionization threshold energy is replaced by the excitation threshold energy, and the oscillator strength is inserted into the formula. This amounts to a choice of the final energy band to be all energy transfers greater than the excitation threshold, modified by the oscillator strength. The inclusion of the oscillator strength is almost certainly motivated by analogy with the Bethe approximation for inelastic collisions with electrons (see below). This approximation is valid only for optically allowed transitions and at high collision energies where the momentum transfer is small. The oscillator strength is introduced in some versions of the ionization formula (e.g. Drawin 1969) including the one adopted by Steenbock and Holweger, but not in others (Drawin 1968; Drawin & Emard 1973).
2.2. Quantum collision theory
In certain cases, quantum scattering theory leads to a dependence of the excitation (and ionization) cross section on the (generalized) oscillator strength. This arises within the first Born approximation (e.g. Seaton 1962; Landau & Lifshitz 1965; Bransden & Joachain 2003), which is valid for highenergy collisions and is obtained within perturbation theory. For electronatom collisions, the unperturbed wave functions are taken in the form of atomic initial/final wave functions (r = { r_{α} } denoting a set of coordinates of atomic electrons labeled by α) multiplied by plane wave functions (2π)^{ − 3/2}exp(ik_{l(u)}R) of the impact particle (R being the impact electron coordinate) and the perturbation operator is the interaction potential U between an impact electron and the target atom, which consists of Coulomb interactions with the nucleus U_{en} and the atomic electrons U_{ee}: U = U_{en} + U_{ee}. The direct scattering probability amplitude for excitation of the atom l → u is (3)where q = k_{l} − k_{u} is the momentum transfer, and the integration is taken over both r and R. Integration over the impact electron coordinate R, considering the Coulomb interactions and the orthogonality of the atomic wave functions, removes the interaction with the nucleus and leaves only the interactions between the incident and atomic electrons (4)The cross section is proportional to the square of the scattering amplitude in Eq. (4), and finally the cross section can be shown to be proportional to the quantity (e.g. Bransden & Joachain 2003) (5)which is called the generalized oscillator strength for reasons that will become apparent. For optically allowed transitions, expansion of the exponential in Eq. (5) can be restricted to the dipole term leading to a relation between the excitation cross section and both the transition dipole moment and the corresponding optical oscillator strength f_{lu}; this relation is known as the Bethe approximation (Bethe 1930; Seaton 1962; Landau & Lifshitz 1965; Bransden & Joachain 2003), and leads to (6)Seaton (1962) provides a discussion of the physical interpretation of this relationship between the collisional and radiative process. For optically forbidden transitions, the higherorder expansion terms resulting from Eq. (5) must be considered thus yielding nonzero excitation/ionization cross sections. Thus, the first Born approximation forms a bridge between the quantum approach and the classical approach of Drawin, both being based on the Coulomb interactions of an impact charged particle with atomic electrons and resulting in relations between the excitation (ionization) cross section and the oscillator strength.
The Born approximation is applicable for highenergy electronatom and electronmolecule collisions, but is valid also for charged heavy particle (ion) – atom/molecule collisions. For lowenergy collisions, which are the main interest of this paper, the Born approximation is no longer valid and gives rather poor results (further explaining the poor agreement between the classical calculations and the experimental data at low energies, see Fig. 1).
At low collision energies, quantum atomic collision theory uses a different approach, in particular, the standard adiabatic (BornOppenheimer) approach, which is based on the separation of the nuclear kinetic energy operator and the electronic fixednuclei Hamiltonian in the total Hamiltonian, finally leading to the molecularstate representation for electronic (fixednuclei) wave functions and nonadiabatic nuclear dynamics (e.g. Mott & Massey 1949; Bates 1962; Macías & Riera 1982; Nikitin & Umansky 1984). Total wave functions of the entire collisional system are expressed in terms of products of electronic molecularstate wave functions and nuclear wave functions, e.g., radial nuclear wave functions F_{l(u)}(R)/R multiplied by angular nuclear wave functions; R being nuclear coordinates. This full quantum treatment requires calculations of quantumchemical data (molecular potentials and nonadiabatic couplings) and integration of coupledchannel scattering equations. The nuclear motion induces nonadiabatic transitions between molecular states. In the asymptotic region, where the internuclear distance is large enough, molecular orbitals are determined by corresponding atomic orbitals, but at short and intermediate distances, where nonadiabatic transitions mainly take place, molecular orbitals are usually represented by mixtures of quite different atomic orbitals, for example, covalent and ionic configurations.
Perturbation theory can also be applied to lowenergy atomic collisions under certain conditions. This can be done in the adiabatic or diabatic representations. In both representations the unperturbed functions are molecularstate electronic wave functions, adiabatic or diabatic ones, and perturbation operators appear from the nuclear kinetic energy operator. In the adiabatic representation, the transition probability amplitude reads (Belyaev 2007) (7)where is the nonadiabatic coupling, which is responsible for nonadiabatic transitions between adiabatic molecular states l and u due to the nuclear motion; see Belyaev (2007) for details. It is clearly seen from Eq. (7) that inelastic transitions in lowenergy atomic collisions are determined by the nuclear kinetic energy operator (i.e. ), by the mixing of electronic molecularstate wave functions expressed by nonadiabatic couplings, and by nuclear wave functions, which substantially deviate from plane waves. In the diabatic representation, the corresponding amplitude has the form (Mott & Massey 1949; Nikitin & Umansky 1984) (8)where H_{lu} is an offdiagonal matrix element calculated with diabatic electronic molecularstate wave functions , which are obtained from adiabatic quantumchemical data by the diabatisation procedure under the condition (Smith 1969). Again, the basis for nonadiabatic transitions is the nuclear kinetic energy operator, the molecularstate wave functions and nuclear wave functions, which should be calculated in this case with diabatic potentials.
Comparison of Eqs. (7) and (8) for lowenergy collisions with Eqs. (3) and (4) for highenergy collisions shows that the basic physics behind the slow atomic collision theory and the first Born approximation, which can be related to the classical description, are different. In highenergy collisions the perturbation theory is based on the unperturbed atomicstate wave functions and an interaction between a charged impact particle and atomic electrons, while in lowenergy atomic collisions the perturbation theory is based on the molecularstate wave functions and an operator appearing from the nuclear kinetic energy operator. If one would use unperturbed electronic atomicstate wave functions within perturbation theory for slow atomic collisions, a number of new terms related to other interactions would appear in addition to the operators mentioned above.
Thus, it is not surprising that, as we will see, the Drawin formula and the quantum calculations give substantially different results. In some cases, simplified nonadiabatic transition models, e.g., the LandauZener model (Landau 1932a,b; Zener 1932), the Demkov (1964), and the Nikitin (1962) ones, which are based on the molecularstate representation, can give reliable inelastic cross sections. Nevertheless, the most accurate results are obtained by means of complete quantum treatments.
2.3. Quantum collision calculations
The basic theory of atomic collisions has been reasonably well understood for some time (e.g. Mott & Massey 1949; Bates 1962, and references therein). However, the detailed calculation of accurate quantumchemical data for quasimolecules required as input to the relevant coupledchannel scattering equations, including for excited states, has been a major hurdle. Further, while the solution of the scattering equations in the standard adiabatic approach appears straightforward, it encounters difficulties such as the electron translation problem (Bates & McCarroll 1958; Belyaev et al. 2001; Belyaev 2010).
Over the last 35 years, our understanding of diatomic hydride molecules and lowenergy H atom collisions with atoms, particularly alkalis, has increased markedly. Detailed theoretical studies of alkali hydride molecular structure began with those by Sachs et al. (1975); Meyer & Rosmus (1975); Olson & Liu (1980); Mó et al. (1985), and these works showed the importance of the ionic configuration, as had been expected from experimental studies (e.g. Mulliken 1936). Studies of ionpair production and mutual neutralisation processes in alkalihydrogen collisions followed shortly after (Janev & Radulovic 1978; Olson & Kimura 1985; Errea et al. 1986; Dickinson et al. 1999; Croft et al. 1999).
Though Errea et al. (1986) briefly looked at excitation of sodium by highenergy H atom collisions, the study of excitation processes came somewhat later. Motivated by Steenbock and Holweger’s work, Fleck et al. (1991) performed an experimental study of Na(3s) + H → Na(3p) + H at low energies (15–1500 eV), though due to experimental difficulties, not down to near the threshold (2.1 eV), which is the relevant regime for the temperatures of interest. The results were compared with theory by Mo & Riera (1992) at energies above 200 eV. Revised experimental data, including results down to 10 eV, and detailed quantum scattering calculations down to the threshold, were presented in Belyaev et al. (1999). These studies demonstrated, through comparison with the LandauZener model and detailed quantum scattering calculations, that the experimental data could be reasonably explained by nonadiabatic transitions associated with avoided crossings, also called pseudocrossings, in the NaH molecular potentials. At pseudocrossings, adiabatic potentials approach each other and the transition probabilities may become large even for slow collisions, despite that the Massey criterion based atomic energy level spacings predicts the collisions to be nearly adiabatic (e.g. Bates 1962, pp. 582 and 611).
Fig. 2 Adiabatic potential energies as a function of internuclear distance R for the lowest ten states of the NaH quasimolecule from pseudopotential calculations by Dickinson et al. (1999). The atomic states and energies at dissociation are shown at the righthand side of the figure. The thick, dashed (red) line shows the 1/R interaction corresponding to the interaction for the pure ionic configuration Na^{+}+H^{−} at long range. 
In the particular case of lowlying states of NaH, as discussed above, these pseudocrossings arise from interaction of the covalent Na+H and ionic Na^{+}+H^{−} configurations, and are thus called avoided ionic crossings. The adiabatic potentials showing these pseudocrossings are shown in Fig. 2, where we see clearly the influence of the ionic configuration, which strongly perturbs the potentials from their behaviour at large internuclear distance. Similar effects are seen in other quasimolecules such as LiH (e.g. Croft et al. 1999) and MgH (e.g. Guitou et al. 2011). As shown by diabatic studies, the underlying ionic curve behaves roughly as 1/R and pseudocrossings occur with all neutral states below the ionic limit, which show a much weaker dependence on R (e.g. Boutalib & Gadéa 1992; Berriche & Gadea 1995). Quantum scattering calculations for Li+H (Belyaev & Barklem 2003) and Na+H (Belyaev et al. 2010) including all states below the ionic state asymptotic limit have found this to be an important mechanism for nonadiabatic transitions. Recent calculations for inelastic Mg+H collisions by Guitou et al. (2011), though presently only done for a small number of states, also find this result.
The basic physical interpretation of this mechanism is that, in a given A+H quasimolecule, starting from a covalent state with appropriate symmetry where the electrons are associated to their respective atoms, at avoided ionic crossings the valence electron associated with the atom A has a probability to tunnel to the H atom resulting in a predominantly ionic charge distribution A^{+}+H^{−}. Later during the collision at a different avoided crossing, there is a probability that the electron may tunnel back to a different covalent molecular state leading finally to a different final state of the atom A, and excitation or deexcitation A(nl) + H → A(n′l′) + H. The electron may also stay with the H atom leading to ionpair production A + H → A^{ + } + H^{ − }. Calculations for Li+H and Na+H mentioned above have generally shown that largevalued cross sections for excitation are typically small compared with those for ionpair production from certain states. This can be easily understood: a largevalued cross section for an excitation process results from passing two avoided ionic crossings, one of them (higherlying) passed with a small transition probability, while for an ionpair production process this region is passed with a large nonadiabatic transition probability due to small adiabatic energy splittings at highly excited crossings.
We conclude by emphasizing that the mechanism found to be important in Li+H, Na+H and Mg+H collisions discussed above, that of the avoided ionic crossings, is an example of more general phenomena in atomic collisions, which will be discussed further below. In He+H inelastic collisions Belyaev et al. (2001) couplings arise due to electronic structure properties completely unrelated to the ionic configuration. Systems such as Li+Na, on the other hand, show a similar mechanism (Belyaev 2010).
2.4. Discussion
Lowenergy atomic collisions differ at a fundamental level from electron collisions, especially those at high energies. This difference arises, naturally, from the presence of a complex projectile and thus the effects of nuclear dynamics in the atomic collisions. As put succinctly by Bates (1962): “They are more complicated in that a wider variety of reaction paths may be followed: for example, the projectile as well as the target may have structure and hence may undergo excitation or ionization; and again, in addition to the possibility of electron exchange, the possibility of charge transfer arises”. Put another way, the effect of the nuclear dynamics is that the electronic structure of the atomatom system during the collision can differ significantly from the electronic structure of the separated atoms, leading to interactions and couplings that are very different in character to those found in the separated atoms. The phenomenon of the avoided ionic crossings, where an electron may tunnel from one nucleus to the other to form an ionic electronic structure leading to strong coupling between adiabatic states, is a particular example of the effects of nuclear dynamics in atomic collisions. In contrast, electron collisions with atoms have no such possibility to change the structure of the target atom in the same way.
Thus, we see there is little in common between the physics of the Drawin formula and the physics of lowenergy atomic collisions, especially H atom collisions, as we now understand them. The Drawin formula is based on a classical model and we see that the fundamental mechanism involved in H atom collisions for the cases studied relates to electron tunnelling, which is quantum mechanical in nature. The Drawin formula extends a theory for electrons, neglecting the atomic nuclei other than in scaling the collision energy. Thus, from a physics point of view it seems the grounds for expectations that Drawin formula can provide good estimates are rather weak, despite the fact that Thomson formula on which it is based is capable of providing orderofmagnitude estimates in some situations for ionization by electron impact.
3. Comparison of results
Fig. 3 Comparison of cross sections for the process Na(3s) + H → Na(3p) + H. The dotted line shows the Drawin cross section and the dashed line the LandauZener cross section. The full lines are two quantum scattering calculations using different input quantumchemical data (MRDCI, which is largest near the threshold, and pseudopotential; see Belyaev et al. 2010, for details). The circles show the experimental data of Fleck et al. (1991); Belyaev et al. (1999) with 1σ error bars. 
Detailed quantum scattering calculations of cross sections, and corresponding rate coefficients, have been done for excitation processes between all states below the ionic limit for Li+H (Belyaev & Barklem 2003; Barklem et al. 2003) and Na+H (Belyaev et al. 2010; Barklem et al. 2010). Recently, calculations for the three lowest states of Mg+H have also been performed (Guitou et al. 2011). In this section we compare the results with the predictions of the Drawin formula. This comparison will be done on two levels. Firstly, comparison of the cross sections, which best elucidates the physics. Secondly, comparison of the rate coefficients, which is of most interest for astrophysical applications, and also allows us to focus on more general integrated properties.
3.1. Cross sections
In this section we compare the Drawin cross sections with other theoretical and experimental results. Note, since Steenbock & Holweger (1984) give only a formula for the rate coefficient, and due to above mentioned problems with Drawin’s derivation, the Drawin cross section formula used here is taken from Lambert (1993, Eqs. (A2) and (A8)). As shown there, these expressions give a formula for the rate coefficient differing from that of Steenbock & Holweger (1984) by only a factor of m_{A}/(m_{A} + m_{H}), which is of order unity.
Figure 3 compares cross sections at low energy for Na(3s) + H → Na(3p) + H, the single case where experimental data is available. We see that quantum scattering calculations and the LandauZener model results agree quite well with experiment. Near the threshold, where there is no experimental data, the quantum scattering calculations show substantial differences depending on which quantumchemical data are used. This highlights the sensitivity of the nearthreshold cross sections to the uncertainties in the quantumchemistry data, and gives an estimate of the uncertainties in the calculated cross sections. In this case the uncertainties seem to be around one or two orders of magnitude. However, it should be noted that the uncertainties vary strongly from transition to transition, and the transitions with the largest cross sections have the smallest uncertainties, as low as a factor of 2. Fortunately, these are the most important from an astrophysical perspective (see discussion in Barklem et al. 2010).
In Figs. 4 and 5, theoretical excitation cross sections for inelastic Li+H and Na+H collisions at low energy are compared for a number of transitions. A major problem of the Drawin formula is immediately evident here: it predicts zero cross sections for optically forbidden transitions (or cross sections much smaller than those for optically allowed transitions if one uses the very small fvalues in some cases). The Drawin formula depends on the absorption oscillator strength, while collisional processes do not follow the same selection rules, and no such dependence is seen in the quantum mechanical data. For the optically allowed transitions the Drawin cross section is typically greater than the quantum scattering results by several orders of magnitude. The Drawin formula gives a reasonable description of the cross section behaviour with collision energy in many cases, though in a few, particularly among transitions involving excited initial states, the quantum cross sections have a much flatter behaviour with collision energy.
Fig. 4 Comparison of quantum scattering cross sections for Li+H from Belyaev & Barklem (2003) (full lines) with those of the Drawin formula (dotted lines) for the 10 transitions between the 4 lowest states of Li. 
Fig. 5 Comparison of quantum scattering cross sections for Na+H from Belyaev et al. (2010) (full lines) with those of the Drawin formula (dotted lines) for the 10 transitions between the 4 lowest states of Na. 
3.2. Rate coefficients
In Fig. 6, rate coefficients for excitation at 6000 K by inelastic Li+H, Na+H and Mg+H collisions, found by integration of the cross sections over a Maxwellian velocity distribution, are shown and compared. In addition to the quantum mechanical calculations and the Drawin formula results, we also calculate rate coefficients from the free electron model described by Kaulakys (1991). Cross sections were calculated according to Eq. (18) in Kaulakys, using nonhydrogenic wave functions in momentum space calculated using the methods of Hoang Binh & Van Regemorter (1997). This model is only valid for transitions in Rydberg atoms, and thus is not valid for the lowlying states being examined here; however, we make a comparison in any case to see if there is any possibility it may give reasonable estimates for transitions among lowlying states. The rate coefficients are plotted against the transition energy, ΔE, since this is a main parameter in determining the rate coefficients, as seen in the figure, due to the strong dependence of the integral on the threshold (e.g., see Eq. (1) where the exponential dependence on is easily seen). This allows us to separate this effect from other more intrinsic secondary effects related to the physics of the collisions. The quantum mechanical results correlate strongly with ΔE, though with a significant scatter of several orders of magnitude around the mean relation. One sees, particularly in the case of Na+H, clear secondary trends among transitions with the same initial state, e.g. 3p in Na+H. These trends result from regular behaviour in the properties of the quasimolecule at avoided crossings (potential splitting, radial coupling), leading to trends in the transition probabilities.
Fig. 6 Rate coefficients at 6000 K from various calculations for excitation transitions between the lowest 9 states of Li and Na due to H atom collisions, plotted against the transition energy. The panels in the left column are for Li+H, and in the right column Na+H. In each panel, different symbols are used to denote the initial state of the transition following the key given in the upper panels. The upper panels show the results from detailed quantum scattering calculations (Barklem et al. 2003, 2010), noting there are no results for Li n = 4 initial states as the data for these transitions were deemed unreliable. The dotted line shows a linear fit to the these data, which is repeated in the lower two panels to aid comparison. The upper panel in the right column also shows the limited data now available for Mg+H calculated from cross sections given in Guitou et al. (2011). Note that the data for the two transitions to the 3p ^{1}P° state should be interpreted as upper limits; however, they are probably not more than an order of magnitude or two lower. No data for Mg+H are plotted in the other panels. The middle panels show the results of the Drawin formula, and the insets show the ratio with the quantum scattering calculations. Forbidden transitions, where the Drawin formula is not applicable, are shown in the bottom of the panels. The lower panels show the results of the Kaulakys free electron model (Kaulakys 1991), and the insets show the ratio with the quantum scattering calculations. 
Comparison with the Drawin formula immediately shows the problem already discussed: lack of predictions for forbidden transitions. The Drawin formula data also show a strong correlation with ΔE, as would be expected since for the case of a given element and temperature, the formula has only two parameters, and f_{lu}. Thus, all deviation from the main trend with ΔE is a result of the fvalue and those transitions deviating significantly from the trend are those with very high or low fvalues. These secondary trends are seen not to correlate with those seen in the quantum mechanical data, which is expected as the physical mechanisms are not related. In general, it is seen that the rate coefficients differ from the quantum ones by up to eight orders of magnitude with a general tendency to be significantly larger. There is also a huge scatter spanning around eight orders of magnitude.
Though the Kaulakys free electron model has the advantage of providing data for all transitions, it compares no better with the detailed quantum scattering results in general, except perhaps for transitions with small transition energies ΔE where agreement is within a few orders of magnitude. Comparison shows a clear tendency for increasingly poor agreement with increasing transition energy, with significant scatter. Thus, there is no suggestion that the Kaulakys free electron model could give reasonable estimates for transitions among lowlying states. It should, however, provide reasonable estimates for transitions among Rydberg states above the ionic limit.
4. Conclusions
The Drawin formula does not contain the essential physics behind direct excitation by H atom collisions at low energies, which is quantum mechanical in character. It compares poorly with the results of the available full quantum scattering calculations based on detailed quantumchemical modelling of the quasimolecule. The estimated uncertainties in the rate coefficients from the quantum mechanical data are not sufficient to explain the differences (see, e.g., Barklem et al. 2010). In general, the Drawin formula significantly overestimates the rates by amounts that vary markedly from transition to transition. Thus, we conclude that the Drawin formula does not provide orderofmagnitude estimates for excitation by H atom collisions. We can even more generally conclude that formulae based on the Drawin expressions are also unlikely to be useful for excitation by other atoms.
The Drawin formula is often used in nonLTE stellar atmosphere applications in a semiempirical fashion with a scaling factor, often denoted S_{H}. The problem of forbidden transitions is sometimes avoided with some arbitrary solution to allow data to be estimated, for example an arbitrary constant fvalue or scaled fvalue (e.g. Steenbock 1985). The value of S_{H} is varied either to assess the possible influence of H atom collisions or to calibrate it against some observational criteria. This procedure, in the absence of a better alternative, is not unreasonable in view of the fact that ΔE has a strong influence on the rate coefficient, and this influence resulting from the velocity distribution is accounted for by the Drawin formula. However, as we have seen, there is a large scatter around this basic behaviour of the rate coefficients, which stems from the physics of the atomic collision, and is not reproduced by the Drawin formula. This means such a procedure will not give good estimates of relative rates for collisional transitions, and could lead to modelling errors if more than one transition is of importance. Furthermore, if other atomic collision processes are present but not modelled, effects caused by other processes may be incorrectly attributed to H atom collisions. For example, in modelling of Li and Na nonLTE line formation it has been shown that the ionpair production processes mentioned in Sect. 2.3, often called charge exchange or charge transfer, are far more important than direct excitation (Barklem et al. 2003; Lind et al. 2011). Deficiencies in other aspects of the spectrum modelling might also be masked by such an approach.
It is clear that there is a need for quantum mechanical estimates for direct excitation and ionpair production processes. The most accurate calculations are naturally provided by the full quantum treatment, i.e. detailed quantumchemical data calculation with detailed quantum scattering calculations. However, in order to provide estimates for a wide range of elements of astrophysical interest, simplified models, based on the physical understanding provided by the detailed calculations, would be required. The reasonable success of the LandauZener model, as seen in Fig. 3, suggests this model, when coupled with a method for producing quantumchemical data, could lie at the basis of attempts to produce such simplified models. However, we emphasise that more experience from the full quantum calculations is complementary and still needed.
Acknowledgments
We acknowledge the important contributions of Alan Dickinson, who sadly passed away last year, to parts of this work. We gratefully acknowledge the support of the Royal Swedish Academy of Sciences, the WennerGren foundation, Göran Gustafssons Stiftelse and the Swedish Research Council. P.S.B. is a Royal Swedish Academy of Sciences Research Fellow supported by a grant from the Knut and Alice Wallenberg Foundation. P.S.B. is grateful for hospitality of the Max Planck Institute for Astrophysics, Garching, where parts of this article were written. A.K.B. also acknowledges support from the Russian Foundation for Basic Research (grant No. 100300807a). The authors acknowledge the role of the SAM collaboration (http://www.anst.uu.se/ulhei450/GaiaSAM/) in stimulating this research through regular workshops.
References
 Athay, R. G., & Canfield, R. C. 1969, ApJ, 156, 695 [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., Dickinson, A. S., & Gadéa, F. X. 2010, A&A, 519, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bates, D. R. 1962, in Atomic and Molecular Processes, ed. D. R. Bates, 550 [Google Scholar]
 Bates, D. R. 1978, Phys. Rep., 35, 305 [NASA ADS] [CrossRef] [Google Scholar]
 Bates, D. R., & McCarroll, R. 1958, Roy. Soc. London Proc. Ser. A, 245, 175 [Google Scholar]
 Bely, O., & van Regemorter, H. 1970, ARA&A, 8, 329 [NASA ADS] [CrossRef] [Google Scholar]
 Belyaev, A. K. 2007, Eur. Phys. J. D, 44, 497 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belyaev, A. K. 2010, Phys. Rev. A, 82, 060701 [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., Grosser, J., Hahne, J., & Menzel, T. 1999, Phys. Rev. A, 60, 2151 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Belyaev, A. K., Egorova, D., Grosser, J., & Menzel, T. 2001, Phys. Rev. A, 64, 052701 [NASA ADS] [CrossRef] [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]
 Berriche, H., & Gadea, F. X. 1995, Chem. Phys. Lett., 247, 85 [NASA ADS] [CrossRef] [Google Scholar]
 Bethe, H. 1930, Annal. Phys., 397, 325 [Google Scholar]
 Boutalib, A., & Gadéa, F. X. 1992, J. Chem. Phys., 97, 1144 [NASA ADS] [CrossRef] [Google Scholar]
 Bransden, B. H., & Joachain, C. J. 2003, Physics of Atoms and Molecules (Harlow: Pearson Education) [Google Scholar]
 Burgess, A., & Percival, I. C. 1968, Adv. At. Mol. Phys., 4, 109 [NASA ADS] [CrossRef] [Google Scholar]
 Croft, H., Dickinson, A. S., & Gadéa, F. X. 1999, J. Phys. B At. Mol. Phys., 32, 81 [Google Scholar]
 Demkov, Y. N. 1964, Sov. Phys. JETP, 18, 138 [Google Scholar]
 Dickinson, A. S., Poteau, R., & Gadéa, F. X. 1999, J. Phys. B At. Mol. Phys., 32, 5451 [NASA ADS] [CrossRef] [Google Scholar]
 Drawin, H.W. 1968, Z. Phys., 211, 404 [Google Scholar]
 Drawin, H. W. 1969, Z. Phys., 225, 483 [Google Scholar]
 Drawin, H. W., & Emard, F. 1973, Phys. Lett. A, 43, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Elwert, G. 1952, Z. Naturforschung Teil A, 7, 432 [Google Scholar]
 Errea, L. F., Méndez, L., Mó, O., & Riera, A. 1986, J. Chem. Phys., 84, 147 [NASA ADS] [CrossRef] [Google Scholar]
 Fleck, I., Grosser, J., Schnecke, A., Steen, W., & Voigt, H. 1991, J. Phys. B At. Mol. Phys., 24, 4017 [Google Scholar]
 Fleischmann, H. H., & Dehmel, R. C. 1972, Z. Phys., 252, 435 [Google Scholar]
 Gehren, T. 1975, A&A, 38, 289 [Google Scholar]
 Guitou, M., Belyaev, A. K., Barklem, P. S., Spielfiedel, A., & Feautrier, N. 2011, J. Phys. B: At. Mol. Opt. Phys., 44, 035202 [NASA ADS] [CrossRef] [Google Scholar]
 HoangBinh, D., & Van Regemorter, H. 1997, J. Phys. B At. Mol. Phys., 30, 2403 [Google Scholar]
 Janev, R. K., & Radulovic, Z. M. 1978, Phys. Rev. A, 17, 889 [NASA ADS] [CrossRef] [Google Scholar]
 Kaulakys, B. 1991, J. Phys. B At. Mol. Phys., 24, L127 [Google Scholar]
 Kunc, J. A., & Soon, W. H. 1991, J. Chem. Phys., 95, 5738 [NASA ADS] [CrossRef] [Google Scholar]
 Lambert, D. L. 1993, Phys. Scr. T, 47, 186 [NASA ADS] [CrossRef] [Google Scholar]
 Landau, L. D. 1932a, Phys. Z. Sowiet., 1, 88 [Google Scholar]
 Landau, L. D. 1932b, Phys. Z. Sowiet., 2, 46 [Google Scholar]
 Landau, L. D., & Lifshitz, E. M. 1965, Quantum Mechanics (Oxford: Pergamon Press) [Google Scholar]
 Lind, K., Asplund, M., Barklem, P. S., & Belyaev, A. K. 2011, A&A, 528, A103 [Google Scholar]
 Macías, A., & Riera, A. 1982, Phys. Rep., 90, 299 [NASA ADS] [CrossRef] [Google Scholar]
 Massey, H. S. W. 1949, Rep. Progr. Phys., 12, 248 [Google Scholar]
 Massey, H. S. W. 1955, Vistas in Astronomy, 1, 277 [NASA ADS] [CrossRef] [Google Scholar]
 McCallion, P., Shah, M. B., & Gilbody, H. B. 1992, J. Phys. B At. Mol. Phys., 25, 1051 [NASA ADS] [CrossRef] [Google Scholar]
 Meyer, W., & Rosmus, P. 1975, J. Chem. Phys., 63, 2356 [NASA ADS] [CrossRef] [Google Scholar]
 Mo, O., & Riera, A. 1992, J. Phys. B At. Mol. Phys., 25, L101 [NASA ADS] [CrossRef] [Google Scholar]
 Mó, O., Riera, A., & Yánez, M. 1985, Phys. Rev. A, 31, 3977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Mott, M. F., & Massey, H. S. W. 1949, The theory of atomic collisions (Oxford: Clarendon) [Google Scholar]
 Mulliken, R. S. 1936, Phys. Rev., 50, 1028 [NASA ADS] [CrossRef] [Google Scholar]
 Nikitin, E. E. 1962, Opt. Spectr., 13, 431 [NASA ADS] [Google Scholar]
 Nikitin, E. E., & Umansky, S. Y. 1984, Theory of slow atomic collisions (Berlin: Springer) [Google Scholar]
 Olson, R. E., & Kimura, M. 1985, Phys. Rev. A, 32, 3092 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Olson, R. E., & Liu, B. 1980, J. Chem. Phys., 73, 2817 [NASA ADS] [CrossRef] [Google Scholar]
 Plaskett, H. H. 1955, MNRAS, 115, 256 [NASA ADS] [Google Scholar]
 Rudge, M. R. 1968, Rev. Mod. Phys., 40, 564 [NASA ADS] [CrossRef] [Google Scholar]
 Sachs, E. S., Hinze, J., & Sabelli, N. H. 1975, J. Chem. Phys., 62, 3367 [NASA ADS] [CrossRef] [Google Scholar]
 Seaton, M. J. 1962, in Atomic and Molecular Processes, ed. D. R. Bates, 375 [Google Scholar]
 Shah, M. B., Elliott, D. S., & Gilbody, H. B. 1987, J. Phys. B At. Mol. Phys., 20, 3501 [Google Scholar]
 Shah, M. B., McCallion, P., Okuno, K., & Gilbody, H. B. 1993, J. Phys. B At. Mol. Phys., 26, 2393 [Google Scholar]
 Smith, F. T. 1969, Phys. Rev., 179, 111 [NASA ADS] [CrossRef] [Google Scholar]
 Steenbock, W. 1985, in Cool Stars with Excesses of Heavy Elements, ed. M. Jaschek, & P. C. Keenan (Dordrecht: D. Reidel Publishing Co.,), Astrophys. Space Sci. Lib., 114, 231 [Google Scholar]
 Steenbock, W., & Holweger, H. 1984, A&A, 130, 319 [NASA ADS] [Google Scholar]
 Thomson, J. J. 1912, Phil. Mag., 23, 449 [Google Scholar]
 Thompson, W. R., Shah, M. B., & Gilbody, H. B. 1995, J. Phys. B At. Mol. Phys., 28, 1321 [NASA ADS] [CrossRef] [Google Scholar]
 van Regemorter, H. 1962, ApJ, 136, 906 [NASA ADS] [CrossRef] [Google Scholar]
 Zener, C. 1932, Roy. Soc. London Proc. Ser. A, 137, 696 [Google Scholar]
All Figures
Fig. 1 Comparison of the classical Thomson (1912) cross section for ionization of H, O, Mg and Fe atoms by electron impact with experimental results of Shah et al. (1987); Thompson et al. (1995); McCallion et al. (1992); Shah et al. (1993), respectively. The plot shows the ratio of the experimental values with respect to the theoretical cross section, plotted against the collision energy E scaled by the ionization energy I. The errors in the experimental values are in all cases typically 5% or less, with an additional uncertainty of around 7–14% in the absolute scale. 

In the text 
Fig. 2 Adiabatic potential energies as a function of internuclear distance R for the lowest ten states of the NaH quasimolecule from pseudopotential calculations by Dickinson et al. (1999). The atomic states and energies at dissociation are shown at the righthand side of the figure. The thick, dashed (red) line shows the 1/R interaction corresponding to the interaction for the pure ionic configuration Na^{+}+H^{−} at long range. 

In the text 
Fig. 3 Comparison of cross sections for the process Na(3s) + H → Na(3p) + H. The dotted line shows the Drawin cross section and the dashed line the LandauZener cross section. The full lines are two quantum scattering calculations using different input quantumchemical data (MRDCI, which is largest near the threshold, and pseudopotential; see Belyaev et al. 2010, for details). The circles show the experimental data of Fleck et al. (1991); Belyaev et al. (1999) with 1σ error bars. 

In the text 
Fig. 4 Comparison of quantum scattering cross sections for Li+H from Belyaev & Barklem (2003) (full lines) with those of the Drawin formula (dotted lines) for the 10 transitions between the 4 lowest states of Li. 

In the text 
Fig. 5 Comparison of quantum scattering cross sections for Na+H from Belyaev et al. (2010) (full lines) with those of the Drawin formula (dotted lines) for the 10 transitions between the 4 lowest states of Na. 

In the text 
Fig. 6 Rate coefficients at 6000 K from various calculations for excitation transitions between the lowest 9 states of Li and Na due to H atom collisions, plotted against the transition energy. The panels in the left column are for Li+H, and in the right column Na+H. In each panel, different symbols are used to denote the initial state of the transition following the key given in the upper panels. The upper panels show the results from detailed quantum scattering calculations (Barklem et al. 2003, 2010), noting there are no results for Li n = 4 initial states as the data for these transitions were deemed unreliable. The dotted line shows a linear fit to the these data, which is repeated in the lower two panels to aid comparison. The upper panel in the right column also shows the limited data now available for Mg+H calculated from cross sections given in Guitou et al. (2011). Note that the data for the two transitions to the 3p ^{1}P° state should be interpreted as upper limits; however, they are probably not more than an order of magnitude or two lower. No data for Mg+H are plotted in the other panels. The middle panels show the results of the Drawin formula, and the insets show the ratio with the quantum scattering calculations. Forbidden transitions, where the Drawin formula is not applicable, are shown in the bottom of the panels. The lower panels show the results of the Kaulakys free electron model (Kaulakys 1991), and the insets show the ratio with the quantum scattering calculations. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.