Issue 
A&A
Volume 612, April 2018



Article Number  A90  
Number of page(s)  6  
Section  Atomic, molecular, and nuclear data  
DOI  https://doi.org/10.1051/00046361/201732365  
Published online  03 May 2018 
Excitation and charge transfer in lowenergy hydrogen atom collisions with neutral iron^{★}
Theoretical Astrophysics, Department of Physics and Astronomy, Uppsala University,
Box 516,
751 20
Uppsala, Sweden
email: Paul.Barklem@physics.uu.se
Received:
27
November
2017
Accepted:
26
December
2017
Data for inelastic processes due to hydrogen atom collisions with iron are needed for accurate modelling of the iron spectrum in latetype stars. Excitation and charge transfer in lowenergy Fe+H collisions is studied theoretically using a previously presented method based on an asymptotic twoelectron linear combination of atomic orbitals model of ioniccovalent interactions in the neutral atomhydrogenatom system, together with the multichannel Landau–Zener model. An extensive calculation including 166 covalent states and 25 ionic states is presented and rate coefficients are calculated for temperatures in the range 1000–20 000 K. The largest rates are found for charge transfer processes to and from two clusters of states around 6.3 and 6.6 eV excitation, corresponding in both cases to active 4d and 5p electrons undergoing transfer. Excitation and deexcitation processes among these two sets of states are also significant.
Key words: atomic data / atomic processes / line: formation / Sun: abundances / stars: abundances
Full Tables and rate coefficient 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/612/A90
© ESO 2018
1 Introduction
Spectral lines of neutral and singly ionised iron are so common in latetype stellar spectra that the iron abundance is used as a standard measure of overall metallicity and lines often used as a diagnostic of stellar effective temperature and surface gravity via excitation and ionisation equilibria. Thus, the iron spectrum can be used to derive the three most important parameters for latetype stellar atmospheres if its formation is well understood. However, it has long been realised from both theoretical and observational evidence that the spectrum of neutral iron, Fe I, departs from local thermodynamical equilibrium (LTE; e.g. Athay & Lites 1972; Rutten & van der Zalm 1984; Ruland et al. 1980; Steenbock 1985; Drake & Smith 1991). The precise magnitude of the corrections due to departures from LTE has been a matter of significant debate, the efficiency of inelastic collision processes due to hydrogen leading to significant uncertainty in line formation calculations (e.g. Holweger 1996; Thévenin & Idiart 1999; Gehren et al. 2001; Collet et al. 2005; Mashonkina et al. 2011). Many studies have employed the Drawin formula, which is an extension of the classical Thomson model for ionisation by electrons to heavy particle collisions between like atoms (Drawin 1968), including excitation (Drawin 1969; Drawin & Emard 1973), and is further extended to unlike atoms by Steenbock & Holweger (1984). This approach has been shown to be unsatisfactory for simple atoms, having poor agreement with experiment and full quantum calculations (e.g. Lambert 1993; Barklem et al. 2011) and not capturing the dominant process, namely charge transfer. However, for lack of better alternatives, this approach is still extensively used in the case of complex atoms such as iron, often with a scaling factor calibrated on spectra of standard stars. This procedure is undesirable, as it is unclear that all discrepancies between theory and observation can be attributed to hydrogen collisions in the form described by the Drawin formula, and therefore doubtful that extrapolation of the calibration to other stars is accurate.
Detailed full quantum structure and scattering calculations of the type that have been employed for simple atoms (e.g. Belyaev et al. 1999; Belyaev & Barklem 2003; Barklem et al. 2003; Belyaev et al. 2010; Barklem et al. 2010; Guitou et al. 2011; Belyaev et al. 2012; Barklem et al. 2012), at the scale required for astrophysical applications, are currently beyond existing methods for complex atoms such as Fe. In this paper, calculations for Fe + H inelastic collision processes are performed using a recently proposed asymptotic model approach (Barklem 2016, 2017). This approach employs a theoretical twoelectron linear combination of atomic orbitals (LCAO) asymptotic model for the molecular structure, together with the multichannel Landau–Zener model of nonadiabatic collision dynamics. The model captures the ionic curve crossing mechanism seen to be important in experiment (Fleck et al. 1991) and in the detailed calculations for simple atoms mentioned. The method has been demonstrated to give reasonable estimates for the processes with large cross sections and rates when compared with detailed calculations,especially charge transfer processes. Thus, these calculations aim to provide a more physically sound basis for understanding the role of collisions with hydrogen atoms in the formation of the iron spectrum in latetype stellar atmospheres.
2 Calculations
The calculations employ the method presented in Barklem (2016, 2017, hereafter B16) with updates as detailed in Barklem (2018). These papers should be consulted for details of the model, including the notation for atomic parameters. The aim was to include all states below the first ionic limit, i.e. below the asymptotic limit of the molecular state dissociating to Fe ^{+} (a^{6}D) + H^{−}, which lies 57 252 cm^{−1} (~7.1 eV, ~2.6 au) above the ground term of Fe + H. A total of 166 states were included, up to and including ^{5} F^{o} at 56 719 cm^{−1}, which has a crossing with the ionic state at over 400 a_{0} and thus the crossing is practically diabatic. The calculations are carried out in LS coupling and states abovethis state might require an alternate coupling scheme, but in any case the crossings involving these states are practically diabatic.
The energy level data employed were taken from National Institute of Standards and Technology (NIST) Atomic Spectra Database (Kramida et al. 2016) and are presented in Table 1 (full table available at the CDS). The required coefficients of fractional parentage are taken from standard tabulations (e.g. Sobelman 1979), and in the case of mixed configurations are combined using the method of Kelly & Armstrong (1959). In order to cover all ionic states appearing as core states among the included configurations, 25 ionic states were included as detailed in Table 1. Many of these cores lead to crossings with covalent states only at very short internuclear distance, and could be excluded from the calculations without significant changes to the results; however, we choose to include these cores for completeness and this provides a check on the coefficients of fractional parentage (correct normalisation).
The resulting molecular term symmetries are presented in Table 2, including those 18 symmetries for which there are both ionic and covalent states, and thus are calculated. The large number of ionic states (25 cores) and symmetries (18) lead to a large number of individual potential and dynamical calculations, although it should be noted that each ionic state has a limited number of symmetries. As for previous calculations, cross sections are computed from thresholds to 100 eV collision energies, and rate coefficients are then calculated and finally summed over all symmetries and cores. Final results are obtained for temperatures in the range 1000–20 000 K with steps of 1000 K.
Input data for the calculations.
Possible symmetries for Fe + H molecular states arising from various asymptotic atomic states and the total statistical weights.
3 Results and discussion
In simpler atoms, the core corresponding to the ground state of the ion is usually completely dominant because it leads to curve crossings at a range of internuclear distances, including those optimal for large transition probabilities at intermediate internuclear distances (~30 to 50 au). This is not the case for complex atoms such as iron, since there are lowlying states of Fe II, and thus other lowlying cores also make contributions. However, the importance certainly decreases as the excitation of the core increases, thus leading more and more only to crossings with lowlying covalent states at shorter and shorter range. For this reason we take examples from the two lowest cores, Fe ^{+} (a^{6}D) and Fe ^{+} (a^{4}F), which lead to the two lowest lying ionic states. The Fe^{+}(a^{6}D) core leads to three symmetries, ^{6}Σ^{+}, ^{6} Π, and ^{6} Δ, and the Fe ^{+} (a^{4}F) core leads to four symmetries, ^{4}Σ^{−}, ^{4} Π, ^{4} Δ, and ^{4} Φ. Both cores have significant statistical weight, thus strengthening their importance.
Figure 1 shows the ^{6}Σ^{+} potentials resulting from the Fe^{+}(a^{6}D) core and the ^{4}Σ^{−} potentials resulting from the Fe^{+}(a^{4}F) core, where the series of ionic curve crossings are clearly seen. Figure 2 shows the corresponding derived ioniccovalent couplings H_{1j} at the avoided crossings, derived in the adiabatic and diabatic representations (see §B of B16), which are in good agreement, except for a few cases at very short range. These cases are unimportant as they cannot lead to large cross sections, and the adopted value from the adiabatic representation is used as it is expected to be more reliable; see B16. For the ^{6} Σ ^{+} Fe ^{+} (a^{6}D) core case it is seen that there are essentially two clusters of crossings: the lowest lying covalent states with crossings at short distance (~10 to 20 au) correspond to Fe states with excitations of 0.1–0.2 au (~24 000 cm^{−1} to 40 000 cm^{−1} ~3 to 5 eV) arising from configuration , and anotherlarge, tightly bunched group of somewhat excited states with crossings at intermediate internuclear distances (~ 30 to 40 au) correspond to Fe states with excitations of 0.22–0.24 (~50 000 cm^{−1} ~6.3 eV) arising from configurations and . Similarly, for the ^{4}Σ^{−} Fe ^{+} (a^{4}F) core case, two clusters are also seen, but do not correspond to the same groups of states. The lowest lying covalent states with crossings at short distance (~ 10 to 20 au) correspond to Fe states with excitations of 0.14–0.22 au (~31 000 cm^{−1} to 48 000 cm^{−1} ~3.8 to 6 eV) arising predominantly from , and an additional tightly bunched group of many excited states with crossings at intermediate internuclear distances (~30 to 40 au) correspond to Fe states with excitations of 0.24–0.25 (~53 000 cm^{−1} ~6.6 eV) arising from and .
Example final results for the rate coefficients are shown in two forms in Figs. 3 and 4. Figure 3 is particularly useful for identifying the dominant cores, as these have large rates for mutual neutralisation. As expected, and discussed above, the core corresponding to the ground state of Fe ^{+} (a^{6}D) gives rise to the most efficient charge transfer (mutual neutralisation and ion pair production) processes of the form Fe ^{+} (a^{6}D) + H^{−}⇄Fe + H, as these show the largest rate coefficients. The processes are to and from the Fe states corresponding to the cluster of states at around ~0.23 au, ~50 000 cm^{−1} ~6.3 eV. The first excited core, Fe^{+}(a^{4}F), also provides some significant rate coefficients, but these coefficients are smaller than those for the ground core and are shifted towards the more excited states (around 0.24 au, 53 000 cm^{−1}, 6.6 eV). The remaining cores show some contributions to charge transfer processes, but are typically orders of magnitude lower.
Figure 4 shows the endothermic rates as a function of change in energy state for the system, ΔE. Earlier calculations for simple atoms, have all shown a rather regular behaviour of the rate coefficients for charge transfer processes with ΔE, that is the asymptotic energy difference between initial and final molecular states. The rate coefficients for ionpair production usually show an arch, peaking at transitions around ΔE = 1 eV (e.g. Barklem et al. 2012), corresponding to transitions resulting from crossings at optimal internuclear distance, i.e. ~30 to 50 au, corresponding to crossings about 1 eV (~0.04 au) below the asymptotic limit for the ionic state. Ezzeddine et al. (2016) have exploited this regular behaviour by fitting the data for simple atoms to permit extrapolation to complex atoms such as Fe. The same structure for the ionpair production processes is seen in Fig. 4, although additional structure is seen at ΔE = 3–5 eV, where rather than a monotonic decrease with ΔE, two branches are seen. This is due to the complicated interplay of the influence of more than one core and the clustering of states correlating to different core configurations. This structure only appears for processes with rather small rates.
As in previous work (B16, Barklem 2018) calculations were also performed with alternate models for the couplings, for example the SEMIEMP and LHJ models of B16, which provide fluctuations and an indication of the uncertainties, at least for the largest rates for which the ionic crossing mechanism is dominant. These fluctuations are not shown in Fig. 4, as the plot becomes too dense. However, the fluctuations of the rate coefficients show very similar behaviour to that for other atoms studied in earlier work, namely, roughly one order of magnitude fluctuations for the largest rates and larger fluctuations for smaller rates.
Fig. 1
Example potentials energies for Fe + H from the LCAO model. Upper panel: ^{6} Σ ^{+} with Fe ^{+}(a^{6}D) core. Lower panel: ^{4}Σ^{−} with Fe ^{+}(a^{4}F) core. 
Fig. 2
Example couplings H_{1j} from the LCAO model plotted against the crossing distance R_{c} for Fe + H states. Upper panel: ^{6}Σ^{+} with Fe ^{+}(a^{6}D) core. Lower panel: ^{4}Σ^{−} with Fe ^{+}(a^{4}F) core. Results are shown for the adiabatic and diabatic models and the final adopted values. 
Fig. 3
Graphical representation of the rate coefficient matrix ⟨σ⟩ (in cm^{3} s^{−1}) for inelastic Fe + H and Fe^{+} + H^{−} collisions at temperature T = 6000 K. Results are from the LCAO asymptotic model. The logarithms in the legend are to base 10. Only every fifth state is labelled. The charge transfer processes, involving initial or final ionic states, are outlined with (coloured) boxes, for ionpair production at the upper right (red) and for mutual neutralisation at the lower left (blue). 
Fig. 4
Rate coefficients ⟨σv⟩, for Fe + Hcollision processes at 6000 K, plotted against the asymptotic energy difference between initial and final molecular states, ΔE. The data are shown for endothermic processes. i.e. excitation and ionpair production. The legend labels the initial state of the transition and processes leading to a final ionic state (ionpair production) are circled. 
4 Concluding remarks
Largescale model calculations have been performed for inelastic Fe + H collision processes using a model approach accounting for the ionic crossing mechanism associated with electron transfer between ionic and covalent configurations. The calculations are expected to give reasonable estimates for the processes with the largest rates, particularly charge transfer and excitation processes between nearlying states. The results show general behaviour very similar to that seen for earlier calculations with simple atoms, although some additional complexity is seen owing to the influence of more than a single core configuration for the iron atom. In general it is found that the largest rates are for charge transfer processes toand from the two clusters of states around 6.3 and 6.6 eV, corresponding in both cases to active 4d and 5p electrons undergoing transfer. The excitation and deexcitation transition rates are largest for transitions between nearby states; these clusters are extensively coupled.
Uncertainties for the largest rates are expected to be about one order of magnitude. The uncertainties for processes with smaller rates are larger, both because of the uncertainties in the present model and uncertainties from mechanisms not captured by the present model. It should be noted that recent work on comparing the centretolimb variation of the solar 777 nm oxygen lines, perhaps indicates the need for larger efficiency of hydrogen collision processes for excitation among lowlying states, and it seems possible that this could be provided by mechanisms other than the ionic crossing mechanism (Amarsi et al., in prep.). This may well be the case for other complex atoms such as iron, and will be the subject of future research, both in terms of comparison with centretolimb variation in the sun and theoretical calculations. Previous work on the solar centretolimb variation for Fe by Lind et al. (2017) made use of preliminary calculations using this method. The present calculation corrects an earlier error (Barklem 2017) and includes a number of other improvements (Barklem 2018) and many more states and cores.
Acknowledgements
This work received financial support from the Swedish Research Council and the project grant “The New Milky Way” from the Knut and Alice Wallenberg Foundation.
References
 Athay, R. G., & Lites, B. W. 1972, ApJ, 176, 809 [NASA ADS] [CrossRef] [Google Scholar]
 Barklem, P. S. 2016, Phys. Rev. A, 93, 042705 [NASA ADS] [CrossRef] [Google Scholar]
 Barklem, P. S. 2017, Phys. Rev. A, 95, 069906 [NASA ADS] [CrossRef] [Google Scholar]
 Barklem, P. S. 2018, A&A, 610, A57 [NASA ADS] [CrossRef] [EDP Sciences] [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]
 Barklem, P. S., Belyaev, A. K., Guitou, M., et al. 2011, A&A, 530, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barklem, P. S., Belyaev, A. K., Spielfiedel, A., Guitou, M., & Feautrier, N. 2012, A&A, 541, A80 [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., 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. S., Spielfiedel, A., et al. 2012, Phys. Rev. A, 85, 32704 [Google Scholar]
 Collet, R., Asplund, M., & Thévenin, F. 2005, A&A, 442, 643 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Drake, J. J., & Smith, G. 1991, MNRAS, 250, 89 [NASA ADS] [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]
 Ezzeddine, R., Plez, B., Merle, T., Gebran, M., & Thévenin, F. 2016, A&A, submitted [arXiv:1612.09302] [Google Scholar]
 Fleck, I., Grosser, J., Schnecke, A., Steen, W., & Voigt, H. 1991, J. Phys. B: At. Mol. Opt. Phys., 24, 4017 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Gehren, T., Korn, A. J., & Shi, J. 2001, A&A, 380, 645 [NASA ADS] [CrossRef] [EDP Sciences] [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]
 Holweger, H. 1996, Phys. Scr. T, 65, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Kelly, P., & Armstrong, B. H. 1959, ApJ, 129, 786 [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., Amarsi, A. M., Asplund, M., et al. 2017, MNRAS, 468, 4311 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mashonkina, L., Gehren, T., Shi, J.R., Korn, A. J., & Grupp, F. 2011, A&A, 528, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruland, F., Biehl, D., Holweger, H., Griffin, R., & Griffin, R. 1980, A&A, 92, 70 [NASA ADS] [Google Scholar]
 Rutten, R. J., & van der Zalm, E. B. J. 1984, A&AS, 55, 143 [NASA ADS] [Google Scholar]
 Sobelman, I. I. 1979, Atomic Spectra and Radiative Transitions (Springer) [CrossRef] [Google Scholar]
 Steenbock, W. 1985, in Cool Stars with Excesses of Heavy Elements, 114, 231 [NASA ADS] [CrossRef] [Google Scholar]
 Steenbock, W., & Holweger, H. 1984, A&A, 130, 319 [NASA ADS] [Google Scholar]
 Thévenin, F., & Idiart, T. P. 1999, ApJ, 521, 753 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Possible symmetries for Fe + H molecular states arising from various asymptotic atomic states and the total statistical weights.
All Figures
Fig. 1
Example potentials energies for Fe + H from the LCAO model. Upper panel: ^{6} Σ ^{+} with Fe ^{+}(a^{6}D) core. Lower panel: ^{4}Σ^{−} with Fe ^{+}(a^{4}F) core. 

In the text 
Fig. 2
Example couplings H_{1j} from the LCAO model plotted against the crossing distance R_{c} for Fe + H states. Upper panel: ^{6}Σ^{+} with Fe ^{+}(a^{6}D) core. Lower panel: ^{4}Σ^{−} with Fe ^{+}(a^{4}F) core. Results are shown for the adiabatic and diabatic models and the final adopted values. 

In the text 
Fig. 3
Graphical representation of the rate coefficient matrix ⟨σ⟩ (in cm^{3} s^{−1}) for inelastic Fe + H and Fe^{+} + H^{−} collisions at temperature T = 6000 K. Results are from the LCAO asymptotic model. The logarithms in the legend are to base 10. Only every fifth state is labelled. The charge transfer processes, involving initial or final ionic states, are outlined with (coloured) boxes, for ionpair production at the upper right (red) and for mutual neutralisation at the lower left (blue). 

In the text 
Fig. 4
Rate coefficients ⟨σv⟩, for Fe + Hcollision processes at 6000 K, plotted against the asymptotic energy difference between initial and final molecular states, ΔE. The data are shown for endothermic processes. i.e. excitation and ionpair production. The legend labels the initial state of the transition and processes leading to a final ionic state (ionpair production) are circled. 

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.