Open Access
Issue
A&A
Volume 689, September 2024
Article Number A100
Number of page(s) 10
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/202451157
Published online 06 September 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

That the negative hydrogen ion H contributes significantly to the opacity in the atmospheres of solar-type stars was first suggested by Wildt (1939), and it is now well known that it dominates the continuous opacity in the visual and near-infrared region of F, G, and K stars, including in the Sun (Chandrasekhar & Breen 1946, Sect. 6); see also, e.g. Unsöld (1955, Sect. 48) and Aller (1963, Sects. 4-19). Theoretical modelling of the atmospheres and spectra of these stars usually assumes local thermodynamic equilibrium (LTE), including for H(e.g. in 1D with the MARCS (Gustafsson et al. 2008) or ATLAS (Kurucz 2013) codes, or in 3D with the STAGGER (Magic et al. 2013) or COBOLD (Freytag et al. 2012) codes).

Potential departures from LTE, non-LTE effects, can take two forms for H. There is a direct effect on H through the effects of non-local radiation on H in competition with local collisional processes, such that the statistical equilibrium of H differs from LTE. There is also an indirect effect via other species. In particular, the ionisation of electron donor elements in non-LTE may be different than that in LTE, leading to a different electron density and thus a different proportion of H compared to a pure LTE situation. It is also worth noting that these direct and indirect effects are not independent; for example, a changed electron density affects the collisional processes in the statistical equilibrium. Furthermore, there would be feedback effects on the structure of the model atmosphere.

For a few late-type stars, model calculations solving the atmospheric structure self-consistently with a non-LTE treatment of a significant number of important species and their line blanketing have been performed (Short & Hauschildt 2003, 2005, 2009) with the PHOENIX code (Hauschildt & Baron 1999). H itself is treated in LTE, that is, indirect effects are included, and no direct effect. The total impact of non-LTE on the continuum flux for the Sun in the visual region compared to a fully LTE model is about 5% (see Fig. 4 of Short & Hauschildt 2005), although the non-LTE PHOENIX model analysed by Pereira et al. (2013) shows much smaller differences relative to LTE, where small differences in the centre-limb variations are also predicted. While the code PHOENIX might in principle be used to produce non-LTE line-blanketed models for FGK stars, the most recent extensive grid is produced in LTE (Husser et al. 2013). This may reflect the computational and physical difficulties, not least in obtaining complete and accurate collision data for all relevant species, which are often not available or have only become available recently, for example for Fe (Wang et al. 2018; Barklem 2018). Thus, LTE models are used in essentially all analyses of late-type stars.

Semi-empirical modelling of the solar atmosphere such as that of Vernazza et al. (1973, 1981) included a non-LTE treatment of H, but these early studies focussed on the upper atmosphere and did not show significant effects on H in the photosphere. However, more recently, Shapiro et al. (2010) found non-LTE effects of about 10% on the H population in the photosphere, and correspondingly large effects on the continuum flux in the visual region. These models considered both direct and indirect effects. The large non-LTE effects on H in the photosphere found in Shapiro et al. (2010) stem predominantly from the indirect effect, which is due to overionisation of metals with a low ionisation potential such as Fe, Si, and Mg (see their Sect. 4.2.1 and Fig. 8, as well as Fig. 1 of Short & Hauschildt 2005).

The validity of any assumption that the relative abundance of H in cool stellar atmospheres can be predicted by LTE, that is, Saha-Boltzmann statistics, or in other words, that the direct effect is negligible, in practice rests on the study of Lambert & Pagel (1968) (hereafter LP68). This is also the basis of the statistical equilibrium model that was used in semi-empirical studies of the solar atmosphere we discussed in the previous paragraph. LP68 reduced the problem to the three most important processes involving H, H, and H2, enabling them to solve the statistical equilibrium problem analytically. Their work was somewhat hampered by uncertain collision rates for the two important collision processes, namely three-body recombination to form H2 from H, and the reverse process of collisional dissociation, 3HH2+H,$3{\rm{H}} \leftrightarrow {{\rm{H}}_2} + {\rm{H,}}$(1)

and associative detachment to form H2 and destroy H, with the reverse process of dissociative attachment to form H, H+HH2+e.${\rm{H}} + {{\rm{H}}^ - } \leftrightarrow {{\rm{H}}_2} + e.$(2)

Despite the uncertainties in the atomic data used in their study, LP68 made a convincing case that direct non-LTE effects on H in late-type stars are small. However, there are several reasons to revisit this question. First, improved atomic data are available for these and other processes; the improvements are often driven by modelling hydrogen in the early Universe (e.g. Lepp et al. 2002; Galli & Palla 2013) and primordial star formation (e.g. Turk et al. 2010; Kreckel et al. 2010), and in molecular clouds (e.g. Dalgarno & Black 2002; Nesterenok et al. 2019). Secondly, the current modelling of stellar spectra in 3D and non-LTE attempts to obtain accuracies at the percent level across a wide range of stellar parameters (e.g. Amarsi et al. 2019; Giribaldi et al. 2021; Wang et al. 2022) for use in many contexts, including large-scale high-resolution spectro-scopic surveys and studies of exoplanets and their host stars (e.g. Nissen & Gustafsson 2018; Jofré et al. 2019). Thus, even small effects on the percent level are of interest in a modern context. LP68 only used scaling relations for stars beyond the Sun and did not consider bright giants or supergiants (i.e. luminosity class I and II). Finally, LP68 only considered continua and not line formation.

Lites & Mihalas (1984) later formulated a more general statistical equilibrium problem, including other species and internal states. A model like this would require a numerical solution, and even though this is possible with modern computers and numerical methods (e.g. Scharmer & Carlsson 1985; Rybicki & Hummer 1992; Arramy et al. 2024), this has not yet been attempted. It is important that these calculations are made to understand the statistical equilibrium of H in detail, including all possible processes. As a basis for further work, however, including understanding of any direct effects seen in more complex calculations, and since the LP68 model is still used (e.g. Shapiro et al. 2010), we see it as an important first step to revisit the analytic model and to understand any differences compared to LP68. This is the subject of this paper.

2 Theory and model

We employed an analytical model for the statistical equilibrium of H similar to that described in Vernazza et al. (1973) (pg. 619, hereafter VAL73), which is a development of the analytical model of LP68. The model includes H, H2, and H, and is shown schematically in Fig. 1. The difference compared to the LP68 model is that it includes the electron collisional detachment process H + e ↔ H + 2e and a more precise expression for the total radiative recombination, in which stimulated emission is treated using the actual mean intensity Jv, rather than the Planck function Bv. As in LP68 and VAL73, we assumed that H2 can be represented by a single level, and this approximation is discussed further in Sect. 2.2.

The bound-free source function for H is given by Sv=2hv3c21bb1ehv/kT1,${S_v} = {{2h{v^3}} \over {{c^2}}}{1 \over {{{{b_ - }} \over {{b_1}}}{e^{hv/kT}} - 1}},$(3)

where the symbols have their usual meanings (e.g. Rutten 2003) with the subscript “−” referring to H and subscript “1” to the ground state of hydrogen. The departure coefficients b are defined b=nn* and b1=n1n1*${b_ - } = {{{n_ - }} \over {n_ - ^*}}{\rm{ and }}{b_1} = {{{n_1}} \over {n_1^*}}{\rm{, }}$(4)

where n are particle densities, and n* denotes particle densities in LTE. We used a definition for the departure coefficients, b = n/n*, that is standard today (e.g. Rutten 2003, Sect. 2.6), rather than the definition used in LP68, where the relation to the modern notation, for example, is bLP68 =b/b1$b_ - ^{{\rm{LP68 }}} = {b_ - }/{b_1}$ (see Rutten’s lecture notes for further details).

Writing the statistical equilibrium equations for H and H2 including the processes shown in Fig. 1, which are described further below, employing the detailed balance relations for inverse processes and solving the equations, we find the departure coefficient for H as b=b1R+rR+r,${b_ - } = {b_1}{{{R^\dag } + r} \over {R + r}},$(5)

where R is the photodetachment rate per H ion R=4πv0σ(v)Jvdvhv,$R = 4\pi \int_{{v_0}}^\infty \sigma (v){J_v}{{{\rm{d}}v} \over {hv}},$(6)

with σ(v) the photodetachment cross section with the frequency threshold v0, and R is R=4πv0σ(v)(2hv3c2+Jv)ehv/kTdvhv,${R^\dag } = 4\pi \int_{{v_0}}^\infty \sigma (v)\left( {{{2h{v^3}} \over {{c^2}}} + {J_v}} \right){e^{ - hv/kT}}{{{\rm{d}}v} \over {hv}},$(7)

such that the total radiative recombination rate (including spontaneous and induced recombination) is n1(n*/n1*)R${n_1}\left( {n_ - ^*/n_1^*} \right){R^\dag }$ (e.g. Mihalas 1978; Rutten 2003; Hubeny & Mihalas 2014). The quantity r is defined as r=neke+n1kH1+Ω,$r = {n_e}{k_e} + {{{n_1}{k_{\rm{H}}}} \over {1 + \Omega }},$(8)

with Ω=n*kHn1*n1k3H,$\Omega = {{n_ - ^*{k_{\rm{H}}}} \over {n_1^*{n_1}{k_{3{\rm{H}}}}}},$(9)

where ke is the rate coefficient (L3T−1) for electron detachment, kH is the rate coefficient (L3T−1) for associative detachment, and k3H is the rate coefficient (L6T−1) for three-body recombination; these rate coefficients are discussed below. The quantity r may be interpreted as an effective total collision rate per H ion, coupling it to H. Thus, b_ as given by Eq. (5), shows the expected behaviour that when b1 = 1, and collisions dominate such that rR and rR, then b_ → 1. When b1 = 1 and the radiative process dominates such that rR and rR, then b_ → R/R. In LTE (Jv = Bv), again denoting LTE values with an asterisk superscript, we have (R)*=4πv0σ(v)Bvdvhv=R*,${\left( {{R^\dag }} \right)^*} = 4\pi \int_{{v_0}}^\infty \sigma (v){B_v}{{{\rm{d}}v} \over {hv}} = {R^*},$(10)

and b_ = b1.

We note that our results are very similar to but not precisely the same as those found in VAL73; their equations seem only to be strictly valid if b1 = 1. In our derivation we distinguished between n1 , the density for ground-state hydrogen atoms, and nH, the total hydrogen atom density, although in practice at the temperatures of interest, they will often be very similar.

Finally, the departure coefficient for H2, denoted bH2${b_{{{\rm{H}}_2}}}$, can be found to be bH2=nH2nH2*=1+(b/b1)Ω1+Ω${b_{{{\rm{H}}_2}}} = {{{n_{{{\rm{H}}_2}}}} \over {n_{{{\rm{H}}_2}}^*}} = {{1 + \left( {{b_ - }/{b_1}} \right)\Omega } \over {1 + \Omega }}{\rm{. }}$(11)

When Ω = 0, and therefore, no coupling of H2 to H, then bH2=1${b_{{{\rm{H}}_2}}} = 1$. In practice, Ω > 0, and it can be seen that bH2${b_{{{\rm{H}}_2}}}$ is always closer to 1 than b_/b1 is to 1, and that if b_/b1 > 1, then bH2>1${b_{{{\rm{H}}_2}}} > 1$, and if b_/b1 < 1, then bH2<1${b_{{{\rm{H}}_2}}} < 1$.

thumbnail Fig. 1

Bubble diagram of the model and the included processes. Collision processes are drawn in blue, and the radiative processes (involving a photon) are drawn in red, where the other particles involved in the process are shown along the relevant line. The corresponding rate coefficients for processes that appear in the model are shown in brackets: R for photodetachment (see Eq. (6) and Sect. 2.1.1), kH for associative detachment (see Sect. 2.1.2), k3H for three-body recombination of hydrogen (see Sect. 2.1.3), and ke for electron collisional detachment (see Sect. 2.1.4).

Table 1

Comparison of rate coefficients for collision processes at temperatures of interest.

2.1 Processes and rates

We now discuss the processes we employed in the model and the data we used for the rates in each case. In Table 1, we compare the rate coefficients for collision processes we adopted with those used by LP68. Even though the data used in LP68 were uncertain, the data they used for the important cases of kH and k3H agree reasonably well with modern data at the temperatures of interest. Electron collisional detachment ke was discussed by LP68, but was not included in the calculations, and Table 1 shows that modern data suggest that this process is even less efficient than the estimate discussed in LP68.

2.1.1 Photodetachment and recombination: H + hv ↔ H + e [σ(v)]

The rate for photodetachment given by Eq. (6) requires the cross section σ(v). The most commonly used data are from Wishart (1979) and cover 16 300 Å (threshold, 0.751 eV) to 1250 Å (almost 10 eV). McLaughlin et al. (2017) have presented new calculations including resonances at high energy (> 10 eV) in the far-UV (first resonance at roughly 1120 Å). The differences compared to the Wishart (1979) data are very small in the visual and UV (photon energies up to 7 eV, roughly redward of 1700 Å; McLaughlin et al. 2017, see Fig. 1), and there are numerous theoretical and experimental results that agree well in the region up to 3 eV (McLaughlin et al. 2017, see Fig. 2). At higher energies, the results depart, not least due to the resonances. For the stars of interest, the radiative rate is dominated by the contribution from the visual, and we used the cross sections from Wishart (1979).

2.1.2 Associative detachment: H + H ↔ H2 + e [kH]

The rate for associative detachment is n-n1kH. We employed the total rate coefficient from Kreckel et al. (2010), the fit given in their Table S1, which is based on the theoretical method developed in Cízek et al. (1998), and shown by Kreckel et al. (2010) to agree very well with their results from a merged-beams experiment. State-resolved data for the reverse process of dissociative attachment were also presented in Launay et al. (1991). These give a lower total associative detachment rate than the data used here (see Fig. S3 of Kreckel et al. (2010)).

LP68 employed a temperature-independent estimate of this rate coefficient based on calculations by Dalgarno & Browne (1967). As listed in Table 1, this value is within 40% of the data adopted in this work, although modern data suggest a significant temperature dependence.

2.1.3 Three-body recombination collisions: 3H ↔ H2 + H [k3H]

The rate for three-body recombination of hydrogen forming H2 is (n1)3k3H. Forrey (2013a) have performed calculations based on a Sturmian theory of three-body recombination (Forrey 2013b), and reported a fit to the rate coefficients valid over the range 300–10000 K, k3H=6×1032T1/4+2×1031T1/2[ cm6 s1 ],${k_{3H}} = 6 \times {10^{ - 32}}{T^{ - 1/4}} + 2 \times {10^{ - 31}}{T^{ - 1/2}}\left[ {{\rm{c}}{{\rm{m}}^6}{{\rm{s}}^{ - 1}}} \right],$(12)

for T in K, and we used this rate coefficient in our modelling. This is estimated to be accurate to roughly a factor of two. This result gives a much flatter temperature dependence than previous results of Jacobs et al. (1967), Palla et al. (1983), Abel et al. (2002), and Flower & Harris (2007). In addition, it is much lower than that from Flower & Harris (by roughly a factor 2–3 at T ~ 5000–7000 K) and significantly higher than that from Abel et al. (slightly more than a factor of 10 at T ~ 5000– 7000 K); see Fig. 1 of Forrey (2013a). State-resolved data for the reverse process of three-body collisional dissociation have also been presented in Bossion et al. (2018).

LP68 employed a semi-empirical estimate of this rate coefficient based on experiments with a large scatter. As listed in Table 1, their estimate generally agrees reasonably well with the data adopted in this work, although modern theoretical data suggest a stronger temperature dependence.

2.1.4 Electron collisional detachment and recombination: H + e ↔ H + 2e [ke]

The rate for collisional detachment of H by electron collisions is n_neke. Various early measurements of the cross section for this process were made at energies well above the threshold (e.g. Tisone & Branscomb 1966; Dance et al. 1967; Tisone & Branscomb 1968; Peart et al. 1970). The first measurement near the threshold was by Walton et al. (1971). Later storage ring experiments probed the low-energy regime in detail, first for D(Andersen et al. 1995; Vejby-Christensen et al. 1996), and later for H (Fritioff et al. 2004).

The experiments at low energy show a very small, possibly zero, cross section at energies just above the electron affinity of H, about 0.75 eV, until an observed threshold at about 2 eV, with the cross section then increasing rapidly to roughly 100 a02$a_0^2$ at about 20 eV. This general behaviour is reasonably explained with a classical over-the-barrier model. The negligible cross section between 0.75 eV and roughly 2 eV due to the Coulomb repulsion means that the low-energy electron is unable to approach the Hion sufficiently closely for the detachment process to occur.

For solar temperature stellar atmospheres, average collision energies are kT ~ 0.2 to 0.6 eV, and they are therefore below threshold for this process. We must therefore expect this process to be inefficient. A reasonable estimate for the cross section σ as a function of collision energy E can be made from the classical model (Andersen et al. 1995; Vejby-Christensen et al. 1996), σ(E)=pπR2max[ 0,(1Uc/E) ],$\sigma (E) = p\pi {R^2}\max \left[ {0,\left( {1 - {U_c}/E} \right)} \right],$(13)

where Uc = 27.2/(R/a0) eV, and the best-fit parameters are p = 0.2 and R = 14.5a0. Using this cross section, we find that the rate coefficient integrated over a Maxwellian velocity distribution can be well represented over the range T = 40–105 K by ke=σv=a1(T/300)a2exp(a3/T)${k_e} = \langle \sigma v\rangle = {a_1}{(T/300)^{{a_2}}}\exp \left( { - {a_3}/T} \right)$(14)

with parameters a1 = 9.28372 × 10−10 cm3s−1, a2 = 0.5, and a3 = 21 787.8 K.

LP68 considered an estimate of this rate coefficient from Thomas & Athay (1961) based on cross-section calculations by Geltman (1960), and concluded that this process was likely to be unimportant. As listed in Table 1, the rate coefficients adopted in this work are lower by roughly two to three orders of magnitude at the temperatures of interest.

2.2 Treatment of H2

LP68, VAL73, and our model all assume that the internal structure of H2, specifically, the vibrational and rotational levels of the ground electronic state, does not need to be considered. Thus, any departure from equilibrium for H2 is captured in a single overall population nH2${n_{{{\rm{H}}_2}}}$ and departure coefficient bH2${b_{{{\rm{H}}_2}}}$.

The basis of this assumption is that H2 is likely to be in, or very close to, internal equilibrium, since the dense system of vibrational levels with quantum numbers v and rotational levels with quantum numbers J are expected to be strongly coupled through collisions. LP68 discussed this assumption at length, and presented a simple two-level model of H2 based on the understanding at the time of which vibrational states are involved in the three-body recombination and associative detachment processes. The discussion was hindered by the fact that no theoretical or laboratory estimates of de-excitation of vibra-tional levels by hydrogen atoms were available. Now, theoretical state-resolved data are available for hydrogen collisions on H2 (Bossion et al. 2018), including excitation and de-excitation, and three-body collisional dissociation (reverse process described in Sect. 2.1.3). We also have state-resolved data for associative detachment (Launay et al. 1991). Based on these data, we extract a somewhat different and more complex picture of what happens in H2 in connection to H formation than found by LP68, which we now discuss.

Based on energy considerations, LP68 expected H2 to be preferentially formed in or destroyed from high vibrational states υ ≥ 8, the highest fully bound vibrational states having υ = 14 (Komasa et al. 2011). This has since been given support from calculations (Orel 1987; Esposito & Capitelli 2009; Bossion et al. 2018; Yuen & Kokoouline 2020). Figure 2 shows data from Bossion et al. (2018) at 5000 K, and it shows that capture occurs most efficiently into states with high υ and low J, which are all clustered in energy close to the dissociation limit.

Furthermore, LP68 expected associative detachment to preferentially form υ ≈ 8, based on examination of the potential curves for H2 and H2${\rm{H}}_2^ - $ in Chen & Peacher (1968). The crossing between the potential curves lies approximately 3.7 eV above the bottom of the H2 potential well, and the (υ, J) = (0,0) state of H2 lies approximately 0.27 eV above the bottom of the well, and thus the crossing could be expected to populate states of H2 with excitation about 3.4 eV above (0,0), that is, roughly E/k ≈ 40 000 K. This led LP68 to expect the υ ≈ 8 state to be most strongly populated. The modern calculations change this picture. Extensive ro-vibrationally resolved calculations were performed by Launay et al. (1991), and the data are shown at 5000 K in Fig. 2. The state with the highest rate for population is (3,19), which lies at roughly this level of excitation. The production in a wide spread of v and J states, and thus energy, is also significant, however. For example, there is significant capture into (2,27) and (3,27), with (3,27) lying very close to the dissociation limit, and thus similar in energy to states with high three-body dissociation rates. In general, Fig. 2 shows that at the temperatures of interest, associative detachment is now expected to form a wide range of states with υ ≈ 1–7 in rotational states J ≈ 7–27, and energies spanning from roughly 2.6 eV (~30 000 K) to the dissociation limit, ~4.48 eV (~52 000 K). Very similar results were reported by Bieniek & Dalgarno (1979), but the temperature grid was much more sparse.

As regards coupling between rotational and vibrational states due to collisions with hydrogen, calculations by Lique (2015) and by Bossion et al. (2018) showed that rotational de-excitation (or in some cases, excitation) rate coefficients for transitions to the nearest-lying state with ∆J = ±2 are about 10−10 cm3 s−1 at temperatures of interest around 5000 K; this transition is favoured as it conserves nuclear spin, that is, ortho→ortho or para→para. Bossion et al. (2018) also showed that vibrational transitions with ∆υ = −1 are also efficient, for example (5,5) to (4,7) has a rate coefficient of about 8 × 10−11 cm3s−1 at around 5000 K. Figure 2 also shows for each state the transition with the highest rate coefficient in the data of Bossion et al. (2018), and this behaviour can be seen.

In summary, modern data suggest that the main collisional pathway for the formation of H is formation of H2 in states with high υ and low J in three-body recombination collisions, 3HH2( high v, low J)+H,$3{\rm{H}} \to {{\rm{H}}_2}({\rm{ high }}\v {\rm{, low }}J) + {\rm{H}},$(15)

followed by collisional de-excitation (or excitation) by hydrogen to states with lower υ and larger J via ∆J = ±2 and ∆υ = −1, that is, H2( high v, low J)+HH2( low v, high J)+H. ${{\rm{H}}_2}({\rm{ high }}\v {\rm{, low }}J) + {\rm{H}} \to {{\rm{H}}_2}({\rm{ low }}\v {\rm{, high }}J) + {\rm{H}}{\rm{. }}$(16)

This is followed by dissociative attachment from a large set of vibrational and rotational states with a wide span in energy, namely H2( low v,highJ)+eH+H,${{\rm{H}}_2}({\rm{ low }}\v ,{\mathop{\rm high}\nolimits} J) + e \to {{\rm{H}}^ - } + {\rm{H,}}$(17)

in total turning a hydrogen atom and an electron into a H ion.

Given this more complex picture, we do not immediately see a way to test the assumption of neglecting the internal structure of H2 other than detailed calculations. Simplification of the system would involve approximations that are difficult to justify in that they must involve assuming at least partial internal equilibrium, which is what we would be attempting to test. We may simply note that LP68 concluded that if the rate coefficient for de-excitation of a collapsed level populated by three-body recombination collisions were k > 10−10 cm3s−1, then this approximation would be valid, within their model, and the individual rates are usually of this magnitude. However, the picture seems more complicated as there are myriad possible internal paths through H2 between H and H.

thumbnail Fig. 2

Graphical representation of the efficiency of collision processes on H2 states at 5000 K, shown in terms of state energy on the y-axis and vibrational quantum number v on the x-axis. The left (blue) vertical bar for each state has a height proportional to the associative detachment rate coefficient (H + H → H2(υ, J) + e) from Launay et al. (1991). The right (orange-red) vertical bar for each state has a height proportional to the collisional dissociation rate coefficient (H2(υ, J) + H → 3H) from Bossion et al. (2018). The length scales are arbitrarily chosen in each case so that the relative proportions can be seen. For each state (υ, J), the de-excitation (most cases) or excitation transition (H2(υ, J) + H → H2(υ′, J′) + H) with the highest rate coefficient k is shown with a grey line; data from Bossion et al. (2018).

3 Calculations

Using this model for the statistical equilibrium, we have calculated departures from LTE for H and attempted to quantify the impact of the departures on stellar spectra. This was done via post-processing of 1D LTE model atmospheres. That is, we only studied the first-order impact on the populations, synthetic spectral lines, and synthetic continua; the feedback on the model atmosphere itself was not taken into account (the so-called restricted non-LTE problem; e.g. Hummer & Rybicki 1971). It should also be emphasised that since LTE model atmospheres were used, the population of neutral hydrogen in the ground state was assumed to be in LTE, b1 = 1; this assumption in the photosphere is strongly supported by calculations (e.g. Fig. 4 of Barklem 2007).

We selected a subset of model atmospheres, covering Teff = [4000, 5000, 6000, 7000] K, log 𝑔 = [1, 2, 3, 4, 5] cm s−2, and [Fe/H] = [−3, −2, −1,0], drawn from a grid of standard 1D LTE MARCS model atmospheres (Gustafsson et al. 2008). The models of FGK-type stars used here are the same as those in, for example, Amarsi et al. (2020), and full details can be found therein. We note that plane-parallel models, computed with microturbu-lence of 1 km s−1, were used for stars with log g ≥ 4.0, while spherical models of solar mass, computed with microturbulence of 2 km s−1 , were used otherwise. Models with log 𝑔 = 1 for Teff = 6000 and 7000 K are not covered, and there is no model with Teff/ log 𝑔/[Fe/H] = 7000/3.0/−1 in this grid, leaving a total of 71 models. We also note that they have a standard composition: namely, solar-scaled models based on the solar abundance mixture of Grevesse et al. 2007, with enhancements to α-elements for metal-poor stars of [α/Fe] = −0.4 × [Fe/H] for −1.0 ≤ [Fe/H] ≤ 0.0 and [α/Fe] = +0.4 for [Fe/H] ≤ −1.0.

The radiative transfer post-processing calculations were carried out using Balder (Amarsi et al. 2018, 2022), a modified version of Multi3D (Botnen & Carlsson 1999; Leenaarts & Carlsson 2009). Background opacities were determined using Blue (Zhou et al. 2023). Two sets of calculations were carried out: the first calculation was to determine the angle-averaged intensity Jν as a function of depth, which was then used to compute departure coefficients for H. In a second step, these departure coefficients were read back into Balder and used in the calculation of synthetic spectra.

For the first set of calculations (to determine Jv), background line opacities were pre-computed with Blue on a grid from 91 nm to 1644 nm with a uniform spacing in log λ such that R = λ/∆λ = 50 000 (corresponding to 144 703 wavelength points). The angle-averaged intensity Jv was then computed with Balder, with Rayleigh scattering in the red wing of Lyman-α included and other processes treated in pure absorption. The calculated Jv were confirmed to agree well with those given by the MARCS code. Given Jv, departure coefficients were then computed following Sect. 2.

For the second set of calculations (to determine synthetic spectra), the departure coefficients b_ determined in the previous step were read back into Balder. Specifically, the contribution to the true H bound-free opacity was modified from nH*σH$n_{{{\rm{H}}^ - }}^*{\sigma _{{{\rm{H}}^ - }}}$ to b×nH*σH${b_ - } \times n_{{{\rm{H}}^ - }}^*{\sigma _{{{\rm{H}}^ - }}}$. The stimulated emission and the emissivity were not corrected (i.e. atomic hydrogen was assumed to be in LTE). With this modification, we carried out a spectrum synthesis of the visual region, roughly 400–800 nm, including 142 Fe II lines following Amarsi et al. (2019), with the atomic data originating from Meléndez & Barbuy (2009). The choice of these lines was arbitrary: They are convenient because they are spread across the optical, and they apparently show no significant deviations from LTE in 1D or 3D for much of the parameter space (Amarsi et al. 2022). The background line opacity (used in the calculation of Jv) was excluded from the final spectrum synthesis for clarity. The theoretical continuum was computed in the same fashion, but without any Fe II line opacity.

4 Results

We now present the results for the departure coefficients and spectra. The spectrum synthesis calculations described above allowed us to study the effects on continua and spectral lines in the visual across our grid of models, and thus, to extract some general tendencies and behaviours.

4.1 Departure coefficients

The results for b_ across our grid are shown in Fig. 3. We focused our attention on the region around optical depth unity (log τ500 = 0) where the visual continuum forms, roughly −1 ≳ log τ500 ≳ 0.5 (e.g. Gray 2022, Chapter 9). While the departure coefficients tend to unity in the very deepest layers, models with high Teff and/or low log 𝑔, show departures from LTE populations with b_ > 1 in this region. Generally, the departures consist of a sharp peak just below optical depth unity and a broader peak higher in the atmosphere. The behaviour at these depths is not very sensitive to metallicity, but there is a tendency for departures to be slightly smaller at lower metallicity.

The reason for these departures can be seen from examination of other quantities shown for two examples in Fig. 4. In particular, the behaviour (but not the magnitude) of the departure coefficient follows R/R, which is greater than unity in this region. This indicates that radiative recombination of electrons with hydrogen to form H at the local temperature exceeds the photodetachment resulting from the non-local radiation field, the mechanism known as overrecombination (Bruls et al. 1992). This is a result of the fact that in this part of the atmosphere, Jv < Bv at infrared wavelengths, combined with the photodetachment cross section peaking at around 850 nm in the infrared, leading to R > R.

We show the results for models with a standard treatment of convection, mixing length efficiency parameter α = 1.5, and also without convection, α = 0. The sharp peak below optical depth unity can be attributed to the significant flattening in temperature gradient deep in the atmosphere due to convection. We note that even when convection is turned off completely, a small residual bump persists near this region that stems from a decline in the effective collision rate r, which is caused by H becoming ionised deeper in the atmosphere.

Finally, we note that these results are not very sensitive to the chosen atomic data within reasonable limits. For example, tests using the data from LP68 found comparable results.

4.2 Spectrum: Continuum and spectral lines

An example of the flux spectrum calculated as described in Sect. 3 for a model with Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0 is shown in Fig. 5. The general behaviour is that the increased abundance of H around optical depth unity increases the opacity, and thus, the continuum is depressed, by about 1% here.

Figure 6 shows the relative changes in the spectral line equivalent widths for this spectrum. We first note the strong trend with wavelength: The relative change increases with increasing wavelength. Secondly, at any given wavelength, strong lines are least affected in relative terms, while weak lines are more affected.

The basic character of these effects, direction and rough magnitude, can be explained as opacity effects via the Milne-Eddington model (see, e.g. Mihalas 1978, Sect. 10.2, Rutten 2003, Sect. 9.1.2, Hubeny & Mihalas 2014, Sect. 17.5). This model assumes that line shape parameters and the ratio of the line and continuum extinction, ηλ=αλl/αλc${\eta _\lambda } = \alpha _\lambda ^l/\alpha _\lambda ^c$, are constant throughout the atmosphere. The model assumes LTE and that the Planck function is a linear function with the continuum optical depth Bλ=B0+B1τλc${B_\lambda } = {B_0} + {B_1}\tau _\lambda ^c$, corresponding to a temperature structure T(τλc)$T\left( {\tau _\lambda ^c} \right)$, where τλc${\tau _\lambda ^c}$ is the continuum optical depth in LTE. The model can now be adjusted in a way analogous to the calculations presented in Sect. 3; that is, the temperature structure is held fixed, while the continuum extinction is adjusted for the non-LTE effect on H with departure coefficient b_, also constant with depth. Then the continuum optical depth in the non-LTE case is bτλc${b_ - }\tau _\lambda ^c$, and thus Bλ=B0+B1bτλc${B_\lambda } = {B_0} + {{{B_1}} \over {{b_ - }}}\tau _\lambda ^c$.

In this model, the flux is then given by Fλ=λ/π=B0+23B11b+ηλ,${F_\lambda } = {{\cal F}_\lambda }/\pi = {B_0} + {2 \over 3}{B_1}{1 \over {{b_ - } + {\eta _\lambda }}},$(18)

where Fλ is the so-called astrophysical flux, while Fλ is the emergent flux at the stellar surface (see Rutten 2003, Sect. 2.1.1). The ratio compared to the case of LTE, b_ = 1, has the following limiting cases: Fλ/FλLTE={ 1, for B0B11+ηλb+ηλ for B0B1. ${F_\lambda }/F_\lambda ^{{\rm{LTE}}} = \left\{ {\matrix{ {1,} \hfill & {{\rm{ for }}{B_0} \gg {B_1}} \hfill \cr {{{1 + {\eta _\lambda }} \over {{b_ - } + {\eta _\lambda }}}} \hfill & {{\rm{ for }}{B_0} \ll {B_1}.} \hfill \cr } } \right.$(19)

The first limiting case, B0B1 is trivial, an isothermal atmosphere. For the second limiting case, closer to reality, B0B1, the behaviour in the deep lines (ηλ → ∞) and continuum (ηλ → 0) gives limiting cases Fλ/FλLTE={ 1, for ηλ1/bfor ηλ0. ${F_\lambda }/F_\lambda ^{{\rm{LTE}}} = \left\{ {\matrix{ {1,} \hfill & {{\rm{ for }}{\eta _\lambda } \to \infty } \hfill \cr {1/{b_ - }} \hfill & {{\rm{for }}{\eta _\lambda } \to 0.} \hfill \cr } } \right.$(20)

Thus, the flux ratio relative to LTE varies between 1 and 1/b_. In the cores of deep lines, the flux is unchanged because the change in continuum extinction has no effect. In the continuum, the flux is changed by the factor 1/b_, and thus, for example if b_ > 1 in the continuum-forming region near optical depth unity, then the flux is reduced and the continuum lowered. This behaviour is shown in Fig. 5, and the approximate magnitude of the effect can be explained. In Figs. 3 and 4, for Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0 around optical depth unity, b_ varies between 1 and about 1.02. Based on the simple model presented above, the continuum is expected to be lowered by a factor in the range 1 to 1/1.02 ≈ 0.98. Fig. 5 shows that the continuum is lower by a factor ≈0.99, roughly 1%.

Following the analysis given in the above references, assuming that the line parameters are constant with depth, and retaining the modification of a departure coefficient for the continuum extinction, we can derive expressions for the equivalent widths of spectral lines and the ratio compared to LTE. Again, it is most instructive to consider limiting cases. For weak lines with η0 ≪ 1, where the subscript 0 denotes the line centre, we obtain WλWλLIE={ 1/(b2) for B0B11/bfor B1B0, ${{{W_\lambda }} \over {W_\lambda ^{{\rm{LIE}}}}} = \left\{ {\matrix{ {1/\left( {b_ - ^2} \right)} \hfill & {{\rm{ for }}{B_0} \gg {B_1}} \hfill \cr {1/{b_ - }} \hfill & {{\rm{for }}{B_1} \gg {B_0},} \hfill \cr } } \right.$(21)

and for strong lines η0 ≫ 1, WλWλLTE={ 1/(b3/2) for B0B11/bfor B1B0, ${{{W_\lambda }} \over {W_\lambda ^{{\rm{LTE}}}}} = \left\{ {\matrix{ {1/\left( {b_ - ^{3/2}} \right)} \hfill & {{\rm{ for }}{B_0} \gg {B_1}} \hfill \cr {1/\sqrt {{b_ - }} } \hfill & {{\rm{for }}{B_1} \gg {B_0},} \hfill \cr } } \right.$(22)

noting that in the first limiting case, the equivalent widths tend towards zero. These results, especially those for the more realistic second limiting case of B1B0, again explain the main effects and their magnitudes as shown in Fig. 6. In particular, the relative effects on the equivalent widths of weak lines are greater than for strong lines. In strong lines, only the continuum formation depth is changed, while for weak lines, both continuum and line formation depths are affected. The dependence on wavelength, with stronger effects on equivalent widths in the red, is explained by the fact that the H extinction increases towards its maximum near 850 nm, leading to a higher (cooler) formation depth and to a reduction of the temperature gradient in the region where the spectrum forms, and thus a reduction in B1.

Finally, in Fig. 7, we attempt to summarise the effects on visual continua across our model grid. In general, the continuum is lower by about 1% for hotter and/or low-gravity models above the upward-sloping diagonal of the figure. Furthermore, from solar metallicity to metal-poor atmospheres, the effects are either similar for Teff = 7000 K models or are reduced for cooler models. The reasons for these behaviours were explained by LP68; in hotter stars, electrons are mainly supplied by hydrogen, and the metallicity therefore matters little. In cooler stars, electrons are supplied by metals, and in metal-poor stars, at a given optical depth in order to retain the same total pressure, the partial pressure due to hydrogen must be higher to compensate for the lower electron pressure. This increases the efficiency of the collision processes involving H (via H2), and thus brings the system closer to LTE. Finally, we note that for some log 𝑔 = 1 models, we have very small positive effects, that is, increase of the continuum, because the departure coefficient dips below unity in the region above optical depth unity (see e.g. the Teff/ log 𝑔/[Fe/H] = 5000/1.0/ − 3.0 model in Fig. 3).

thumbnail Fig. 3

Results for departure coefficients b_ with the logarithm of the optical depth at 500 nm, log τ500, across a grid of MARCS model atmospheres. In the upper right corner of each panel Teff/ log 𝑔 is given, and each panel shows plots for different metallicity [Fe/H] = 0, −1, −2, and −3. The thickest (purple) line shows 0, and the thinnest (orange) line shows -3. This grid does not include a model with Teff/ log 𝑔/[Fe/H] = 7000/3.0/ − 1.

thumbnail Fig. 4

Non-LTE mechanism for a model at roughly solar parameters (Teff = 5750 K, log 𝑔 = 4.5 cm s−2, [Fe/H]= 0) in the left panels and for a hotter star (Teff = 7000 K, log 𝑔 = 4.0 cm s−2, [Fe/H]= 0) showing larger departures from LTE in the right panels. In each case, the upper panel shows the run of the departure coefficient b_ with the logarithm of the optical depth at 500 nm, log τ500. The second panel shows the ratio R/R. The third panel shows the photodetachment rate per H ion R and effective collisional rate per H ion r. The lowest panel shows the temperature T. In each case, models with standard convection for MARCS models (mixing length parameter α = 1.5, thinner blue lines) and without convection (α = 0, thicker orange lines) are shown.

thumbnail Fig. 5

Effect on the spectrum (continuum and Fe II lines) in the visual, roughly 400–800 nm, for a model with Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0. The upper plot shows the emergent flux at the stellar surface 𝓕λ. The blue line shows H treated in non-LTE (𝓕λ), and the dotted red line shows it in LTE (λLTE)$\left( {{\cal F}_\lambda ^{{\rm{LTE}}}} \right)$. The lower panel shows the ratio λ/λLTE${{\cal F}_\lambda }/{\cal F}_\lambda ^{{\rm{LTE}}}$.

thumbnail Fig. 6

Effect on the equivalent widths of Fe II lines in the visual, roughly 400 to 800 nm, for a model with Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0. Each circle shows the ratio of the equivalent width for H treated in non-LTE (Wλ) to that in LTE (WλLTE)$\left( {W_\lambda ^{{\rm{LTE}}}} \right)$, i.e. Wλ/WλLTE${W_\lambda }/W_\lambda ^{{\rm{LTE}}}$. The circle diameter is proportional to the equivalent width.

thumbnail Fig. 7

Summary of the effects of departures from LTE for H on the visual continuum, roughly 400-800 nm, across the grid of models. For each Teff, log 𝑔 combination, the results are shown for four [Fe/H] each slightly offset in the vertical direction such that the highest metallicity is highest in the plot. For each case, a circle marks the position of no change, and then an asterisk displaced in the horizontal direction gives the maximum change in the continuum relative to LTE. The scale is shown in the upper left corner, and the position of a −1% change is marked in each case with a dashed red line. For example, for the case given in Fig. 5, Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0, the uppermost line in this group shows a relative change of just below −1 %, as shown in the lower panel of Fig. 5.

5 Conclusions

We have revisited the statistical equilibrium of H in late-type stellar atmospheres, a question that has been largely untouched since the pioneering work of LP68. Using a slightly modified version of their analytical model including H, H2 and H, together with modern atomic data and a grid of 1D LTE theoretical stellar atmosphere models, we find direct non-LTE effects on populations, continua, and spectral lines of about 1–2% in stars with higher Teff and/or lower log g within the grid. The effects in models for solar parameters are smaller by a factor of 10, about 0.1–0.2%, and for models with lower Teff and/or higher log g, they are practically absent. The non-LTE effects found in the spectrum-forming regions in our calculations originate from overrecombination, that is, radiative recombination of electrons with hydrogen to form H exceeding photodetachment.

Modern atomic data are not a source of significant differences in the results compared to LP68; the atomic data are now significantly more secure than those used by LP68, however, which places the calculations on a firm footing. In addition, the detailed data for processes on H2 resolved with vibrational and rotational states paint a more complete and complex picture of the role of H2 in the formation and equilibrium of H than was found in LP68.

We emphasise that this work only considered the direct effect, is only for a limited reaction network where H is assumed to be in LTE, and was furthermore conducted for 1D models. Lites & Mihalas (1984) formulated a more complete reaction network and statistical equilibrium problem including H2+${\rm{H}}_2^ + $, H+, and internal states of H out of LTE, but no calculations were attempted. This should be done, ideally with the inclusion of internal states of H2, and possibly also H3+${\rm{H}}_3^ + $. In addition, mutual neutralisation processes of H with metal ions, H + X+ ↔ H + X* have in recent years been demonstrated to be important in non-LTE studies of neutral metal atoms, (e.g. Barklem et al. 2003, 2011; Belyaev et al. 2014), and they might affect the statistical equilibrium of H.

Furthermore, the 1D modelling presented here suggests a connection to the temperature structure as modified by convection, and investigation in 3D models may be important. Finally, modelling including both direct and indirect effects, including feedback on the atmospheric structure, should be attempted. Based on the results of Short & Hauschildt (2005) and Shapiro et al. (2010), this requires simultaneous accurate non-LTE modelling of important electron donors such as Fe, Si, and Mg. This is a very significant challenge, even in 1D. Consideration of all these additional factors may significantly change or remove the direct effects seen here, but it is important to understand the physics in detail, both in its own right for understanding stellar spectra and for applications in other astrophysical problems.

Acknowledgements

We thank Bengt Edvardsson for calculating MARCS models with different convection parameters, and doing the first calculations of Jν used in preliminary studies, and Bengt Gustafsson for pointing out the role of convection in determining the atmospheric temperature structures deep in the atmosphere. We thank Yohann Scribano for providing the H2 collision data. Maríla Carlos and Sema Caliskan gave comments that improved the figures. P.S.B. and A.M.A. acknowledge support from the Swedish Research Council (VR 2020-03404 and 2020-03940). This research was supported by computational resources provided by the Australian Government through the National Computational Infrastructure (NCI) under the National Computational Merit Allocation Scheme and the ANU Merit Allocation Scheme (project y89).

References

  1. Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 [CrossRef] [Google Scholar]
  2. Aller, L. H. 1963, Astrophysics. The Atmospheres of the Sun and Stars (New York: The Ronald Press Company) [Google Scholar]
  3. Amarsi, A. M., Nordlander, T., Barklem, P. S., et al. 2018, A&A, 615, A139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Amarsi, A. M., Nissen, P. E., & Skúladóttir, Á. 2019, A&A, 630, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Amarsi, A. M., Lind, K., Osorio, Y., et al. 2020, A&A, 642, A62 [EDP Sciences] [Google Scholar]
  6. Amarsi, A. M., Liljegren, S., & Nissen, P. E. 2022, A&A, 668, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Andersen, L. H., Mathur, D., Schmidt, H. T., & Vejby-Christensen, L. 1995, Phys. Rev. Lett., 74, 892 [NASA ADS] [CrossRef] [Google Scholar]
  8. Arramy, D., Rodríguez, J. d. 1. C., & Leenaarts, J. 2024, A&A, in press https://doi.org/10.1051/0004-6361/202449963 [Google Scholar]
  9. Barklem, P. S. 2007, A&A, 466, 327 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Barklem, P. S. 2018, A&A, 612, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Barklem, P. S., Belyaev, A. K., & Asplund, M. 2003, A&A, 409, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Barklem, P. S., Belyaev, A. K., Guitou, M., et al. 2011, A&A, 530, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Belyaev, A. K., Yakovleva, S. A., & Barklem, P. S. 2014, A&A, 572, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Bieniek, R. J., & Dalgarno, A. 1979, ApJ, 228, 635 [NASA ADS] [CrossRef] [Google Scholar]
  15. Bossion, D., Scribano, Y., Lique, F., & Parlant, G. 2018, MNRAS, 480, 3718 [CrossRef] [Google Scholar]
  16. Botnen, A., & Carlsson, M. 1999, Astrophys. Space Sci. Lib., 240, 379 [NASA ADS] [CrossRef] [Google Scholar]
  17. Bruls, J. H. M. J., Rutten, R. J., & Shchukina, N. G. 1992, A&A, 265, 237 [NASA ADS] [Google Scholar]
  18. Chandrasekhar, S., & Breen, F. H. 1946, ApJ, 104, 430 [CrossRef] [Google Scholar]
  19. Chen, J. C., & Peacher, J. L. 1968, Phys. Rev., 168, 56 [NASA ADS] [CrossRef] [Google Scholar]
  20. Cízek, M., Horácek, J., & Domcke, W. 1998, J. Phys. B: At. Mol. Opt. Phys., 31, 2571 [CrossRef] [Google Scholar]
  21. Dalgarno, A., & Black, J. H. 2002, Rep. Prog. Phys., 39, 573 [Google Scholar]
  22. Dalgarno, A., & Browne, J. C. 1967, ApJ, 149, 231 [NASA ADS] [CrossRef] [Google Scholar]
  23. Dance, D. F., Harrison, M. F. A., & Rundel, R. D. 1967, Proc. R. Soc. Lond. A, 299, 525 [NASA ADS] [CrossRef] [Google Scholar]
  24. Esposito, F., & Capitelli, M. 2009, J. Phys. Chem. A, 113, 15307 [CrossRef] [Google Scholar]
  25. Flower, D. R., & Harris, G. J. 2007, MNRAS, 377, 705 [NASA ADS] [CrossRef] [Google Scholar]
  26. Forrey, R. C. 2013a, ApJ, 773, L25 [NASA ADS] [CrossRef] [Google Scholar]
  27. Forrey, R. C. 2013b, Phys. Rev. A, 88, 052709 [CrossRef] [Google Scholar]
  28. Freytag, B., Steffen, M., Ludwig, H.-G., et al. 2012, J. Comp. Phys., 231, 919 [NASA ADS] [CrossRef] [Google Scholar]
  29. Fritioff, K., Sandström, J., Andersson, P., et al. 2004, Phys. Rev. A, 69, 042707 [CrossRef] [Google Scholar]
  30. Galli, D., & Palla, F. 2013, ARA&A, 51, 163 [NASA ADS] [CrossRef] [Google Scholar]
  31. Geltman, S. 1960, Proc. Phys. Soc., 75, 67 [NASA ADS] [CrossRef] [Google Scholar]
  32. Giribaldi, R. E., da Silva, A. R., Smiljanic, R., & Cornejo Espinoza, D. 2021, A&A, 650, A194 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Gray, D. F. 2022, The Observation and Analysis of Stellar Photospheres (Cambridge: Cambridge University Press) [Google Scholar]
  34. Grevesse, N., Asplund, M., & Sauval, A. J. 2007, Space Sci. Rev., 130, 105 [Google Scholar]
  35. Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Hauschildt, P. H., & Baron, E. 1999, J. Comp. App. Math., 109, 41 [NASA ADS] [CrossRef] [Google Scholar]
  37. Hubeny, I. & Mihalas, D. 2014, Theory of Stellar Atmospheres (Princeton, NJ: Princeton University Press) [Google Scholar]
  38. Hummer, D. G., & Rybicki, G. 1971, ARA&A, 9, 237 [NASA ADS] [CrossRef] [Google Scholar]
  39. Husser, T. O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Jacobs, T. A., Giedt, R. R., & Cohen, N. 1967, J. Chem. Phys., 47, 54 [NASA ADS] [CrossRef] [Google Scholar]
  41. Jofré, P., Heiter, U., & Soubiran, C. 2019, ARA&A, 57, 571 [Google Scholar]
  42. Komasa, J., Piszczatowski, K., Lach, G., et al. 2011, J. Chem. Theory Comput., 7, 3105 [NASA ADS] [CrossRef] [Google Scholar]
  43. Kreckel, H., Bruhns, H., Cížek, M., et al. 2010, Science, 329, 69 [NASA ADS] [CrossRef] [Google Scholar]
  44. Kurucz, R. L. 2013, Astrophysics Source Code Library [record ascl:1303.024] [Google Scholar]
  45. Lambert, D. L., & Pagel, B. E. J. 1968, MNRAS, 141, 299 [NASA ADS] [Google Scholar]
  46. Launay, J. M., Le Dourneuf, M., & Zeippen, C. J. 1991, A&A, 252, 842 [NASA ADS] [Google Scholar]
  47. Leenaarts, J., & Carlsson, M. 2009, ASP Conf. Ser., 415, 87 [Google Scholar]
  48. Lepp, S., Stancil, P. C., & Dalgarno, A. 2002, J. Phys. B At. Mol. Opt. Phys., 35, 57 [Google Scholar]
  49. Lique, F. 2015, MNRAS, 453, 810 [NASA ADS] [CrossRef] [Google Scholar]
  50. Lites, B. W., & Mihalas, D. 1984, Sol. Phys., 93, 23 [NASA ADS] [Google Scholar]
  51. Magic, Z., Collet, R., Asplund, M., et al. 2013, A&A, 557, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. McLaughlin, B. M., Stancil, P. C., Sadeghpour, H. R., & Forrey, R. C. 2017, J. Phys. B At. Mol. Opt. Phys., 50, 114001 [CrossRef] [Google Scholar]
  53. Meléndez, J., & Barbuy, B. 2009, A&A, 497, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Mihalas, D. 1978, Stellar Atmospheres, 2nd edn. (San Francisco: Freeman) [Google Scholar]
  55. Nesterenok, A. V., Bossion, D., Scribano, Y., & Lique, F. 2019, MNRAS, 489, 4520 [NASA ADS] [CrossRef] [Google Scholar]
  56. Nissen, P. E., & Gustafsson, B. 2018, A&ARv, 26, 6 [Google Scholar]
  57. Orel, A. E. 1987, J. Chem. Phys., 87, 314 [NASA ADS] [CrossRef] [Google Scholar]
  58. Palla, F., Salpeter, E. E., & Stahler, S. W. 1983, ApJ, 271, 632 [NASA ADS] [CrossRef] [Google Scholar]
  59. Peart, B., Walton, D. S., & Dolder, K. T. 1970, J. Phys. B At. Mol. Phys., 3, 1346 [NASA ADS] [CrossRef] [Google Scholar]
  60. Pereira, T. M. D., Asplund, M., Collet, R., et al. 2013, A&A, 554, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Rutten, R. J. 2003, Radiative : Transfer in Stellar Atmospheres (UtrechtLecture Notes Utrecht University) [Google Scholar]
  62. Rybicki, G. B., & Hummer, D. G. 1992, A&A, 262, 209 [NASA ADS] [Google Scholar]
  63. Scharmer, G. B., & Carlsson, M. 1985, J. Comp. Phys., 59, 56 [NASA ADS] [CrossRef] [Google Scholar]
  64. Shapiro, A. I., Schmutz, W., Schoell, M., Haberreiter, M., & Rozanov, E. 2010, A&A, 517, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Short, C. I., & Hauschildt, P. H. 2003, ApJ, 596, 501 [Google Scholar]
  66. Short, C. I., & Hauschildt, P. H. 2005, ApJ, 618, 926 [NASA ADS] [CrossRef] [Google Scholar]
  67. Short, C. I., & Hauschildt, P. H. 2009, ApJ, 691, 1634 [CrossRef] [Google Scholar]
  68. Thomas, R. N., & Athay, R. G. 1961, Physics of the Solar Chromosphere (New York: Interscience) [Google Scholar]
  69. Tisone, G., & Branscomb, L. M. 1966, Phys. Rev. Lett., 17, 236 [NASA ADS] [CrossRef] [Google Scholar]
  70. Tisone, G. C., & Branscomb, L. M. 1968, Phys. Rev., 170, 169 [NASA ADS] [CrossRef] [Google Scholar]
  71. Turk, M. J., Clark, P., Glover, S. C. O., et al. 2010, ApJ, 726, 55 [Google Scholar]
  72. Unsöld, A. 1955, Physik Der Sternatmospharen, Mit Besonderer Berucksichti-gung Der Sonne (Berlin: Springer) [Google Scholar]
  73. Vejby-Christensen, L., Kella, D., Mathur, D., et al. 1996, Phys. Rev. A, 53, 2371 [CrossRef] [Google Scholar]
  74. Vernazza, J. E., Avrett, E. H., & Loeser, R. 1973, ApJ, 184, 605 [NASA ADS] [CrossRef] [Google Scholar]
  75. Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 [Google Scholar]
  76. Walton, D. S., Peart, B., & Dolder, K. T. 1971, J. Phys. B At. Mol. Phys., 4, 1343 [NASA ADS] [CrossRef] [Google Scholar]
  77. Wang, K., Bartschat, K., & Zatsarinny, O. 2018, ApJ, 867, 63 [NASA ADS] [CrossRef] [Google Scholar]
  78. Wang, H. S., Quanz, S. P., Yong, D., et al. 2022, MNRAS, 513, 5829 [NASA ADS] [Google Scholar]
  79. Wildt, R. 1939, ApJ, 90, 611 [NASA ADS] [CrossRef] [Google Scholar]
  80. Wishart, A. W. 1979, MNRAS, 187, 59P [NASA ADS] [CrossRef] [Google Scholar]
  81. Yuen, C. H., & Kokoouline, V. 2020, Phys. Rev. A, 101, 042709 [NASA ADS] [CrossRef] [Google Scholar]
  82. Zhou, Y., Amarsi, A. M., Aguirre Børsen-Koch, V., et al. 2023, A&A, 677, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Comparison of rate coefficients for collision processes at temperatures of interest.

All Figures

thumbnail Fig. 1

Bubble diagram of the model and the included processes. Collision processes are drawn in blue, and the radiative processes (involving a photon) are drawn in red, where the other particles involved in the process are shown along the relevant line. The corresponding rate coefficients for processes that appear in the model are shown in brackets: R for photodetachment (see Eq. (6) and Sect. 2.1.1), kH for associative detachment (see Sect. 2.1.2), k3H for three-body recombination of hydrogen (see Sect. 2.1.3), and ke for electron collisional detachment (see Sect. 2.1.4).

In the text
thumbnail Fig. 2

Graphical representation of the efficiency of collision processes on H2 states at 5000 K, shown in terms of state energy on the y-axis and vibrational quantum number v on the x-axis. The left (blue) vertical bar for each state has a height proportional to the associative detachment rate coefficient (H + H → H2(υ, J) + e) from Launay et al. (1991). The right (orange-red) vertical bar for each state has a height proportional to the collisional dissociation rate coefficient (H2(υ, J) + H → 3H) from Bossion et al. (2018). The length scales are arbitrarily chosen in each case so that the relative proportions can be seen. For each state (υ, J), the de-excitation (most cases) or excitation transition (H2(υ, J) + H → H2(υ′, J′) + H) with the highest rate coefficient k is shown with a grey line; data from Bossion et al. (2018).

In the text
thumbnail Fig. 3

Results for departure coefficients b_ with the logarithm of the optical depth at 500 nm, log τ500, across a grid of MARCS model atmospheres. In the upper right corner of each panel Teff/ log 𝑔 is given, and each panel shows plots for different metallicity [Fe/H] = 0, −1, −2, and −3. The thickest (purple) line shows 0, and the thinnest (orange) line shows -3. This grid does not include a model with Teff/ log 𝑔/[Fe/H] = 7000/3.0/ − 1.

In the text
thumbnail Fig. 4

Non-LTE mechanism for a model at roughly solar parameters (Teff = 5750 K, log 𝑔 = 4.5 cm s−2, [Fe/H]= 0) in the left panels and for a hotter star (Teff = 7000 K, log 𝑔 = 4.0 cm s−2, [Fe/H]= 0) showing larger departures from LTE in the right panels. In each case, the upper panel shows the run of the departure coefficient b_ with the logarithm of the optical depth at 500 nm, log τ500. The second panel shows the ratio R/R. The third panel shows the photodetachment rate per H ion R and effective collisional rate per H ion r. The lowest panel shows the temperature T. In each case, models with standard convection for MARCS models (mixing length parameter α = 1.5, thinner blue lines) and without convection (α = 0, thicker orange lines) are shown.

In the text
thumbnail Fig. 5

Effect on the spectrum (continuum and Fe II lines) in the visual, roughly 400–800 nm, for a model with Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0. The upper plot shows the emergent flux at the stellar surface 𝓕λ. The blue line shows H treated in non-LTE (𝓕λ), and the dotted red line shows it in LTE (λLTE)$\left( {{\cal F}_\lambda ^{{\rm{LTE}}}} \right)$. The lower panel shows the ratio λ/λLTE${{\cal F}_\lambda }/{\cal F}_\lambda ^{{\rm{LTE}}}$.

In the text
thumbnail Fig. 6

Effect on the equivalent widths of Fe II lines in the visual, roughly 400 to 800 nm, for a model with Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0. Each circle shows the ratio of the equivalent width for H treated in non-LTE (Wλ) to that in LTE (WλLTE)$\left( {W_\lambda ^{{\rm{LTE}}}} \right)$, i.e. Wλ/WλLTE${W_\lambda }/W_\lambda ^{{\rm{LTE}}}$. The circle diameter is proportional to the equivalent width.

In the text
thumbnail Fig. 7

Summary of the effects of departures from LTE for H on the visual continuum, roughly 400-800 nm, across the grid of models. For each Teff, log 𝑔 combination, the results are shown for four [Fe/H] each slightly offset in the vertical direction such that the highest metallicity is highest in the plot. For each case, a circle marks the position of no change, and then an asterisk displaced in the horizontal direction gives the maximum change in the continuum relative to LTE. The scale is shown in the upper left corner, and the position of a −1% change is marked in each case with a dashed red line. For example, for the case given in Fig. 5, Teff = 7000 K, log 𝑔 = 4.0, and [Fe/H] = 0, the uppermost line in this group shows a relative change of just below −1 %, as shown in the lower panel of Fig. 5.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.