Free Access
Issue
A&A
Volume 586, February 2016
Article Number A156
Number of page(s) 11
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/201526403
Published online 11 February 2016

© ESO, 2016

1. Introduction

Spectral lines in stellar atmospheres frequently form in conditions that may deviate significantly from local thermodynamic equilibrium (LTE). This is normally taken into account by solving the statistical equilibrium equations, coupled with non-LTE (NLTE) radiative transfer computations. Fortunately, except for the hottest stars, NLTE effects have relatively minor influence on the thermodynamic structures of the model atmospheres (e.g., Short & Hauschildt 2006). However, in many cases the NLTE treatment in the computation of synthetic spectral line profiles is very important, especially at lower metallicities where NLTE effects may be significant (see, e.g., Asplund 2005, for a review). So far, however, attempts to include NLTE effects have been mostly limited to using 1D hydrostatic model atmospheres.

At the same time, recent studies made with the 3D hydrodynamical model atmospheres in LTE have revealed that various dynamical phenomena may also be very important in the spectral line formation taking place in convective stellar atmospheres (e.g., Asplund et al. 1999; Collet et al. 2007, 2009; Behara et al. 2010; Beeck et al. 2013; Dobrovolskas et al. 2013; Magic et al. 2013). It would therefore be desirable if abundances of chemical elements in stellar atmospheres could be derived using the full 3D NLTE approach, i.e., with NLTE spectral synthesis computations performed using 3D hydrodynamical model atmospheres, which would be the most realistic way to model spectral line formation. This, however, has rarely been possible since the majority of current spectral synthesis codes lack the capability to perform the statistical equilibrium computations in 3D NLTE, even though several successful steps in this direction have already been made, with very likely more to come in the near future (see, for example, Asplund et al. 2003, 2004; Cayrel et al. 2007; Steffen et al. 2010, 2012; Lind et al. 2012, 2013; Holzreuter & Solanki 2013; Prakapavičius et al. 2013; Steffen et al. 2015; Amarsi et al. 2015). Frequently then, various simplifications are still used to substitute for the full 3D NLTE approach. For example, sometimes elemental abundances are derived in 1D NLTE and subsequently corrected for the 3D hydrodynamical effects, with the size of these corrections evaluated under the assumption of LTE (3D+NLTE approach); or where average ⟨3D⟩ model atmospheres are used with the standard 1D NLTE analysis tools. However, since the properties of spectral line formation in full 3D NLTE are still poorly explored, many details of the interaction between hydrodynamical and NLTE effects in 3D are still unknown. Therefore, it is not clear whether such strategies may be generalized to substitute for the full 3D NLTE computations.

Lithium is an important cosmological tracer element, and one way of measuring its primordial (cosmological) abundance is from stellar spectra, mostly of main sequence (MS), main sequence turn-off (TO), and subgiant branch (SGB) stars; however, red giant branch (RGB) stars, may be useful too (see, for example, Mucciarelli et al. 2012). Lithium is known to be prone to both 3D hydrodynamical and NLTE effects, and it is thought that these effects typically tend to work in opposite directions in stellar atmospheres, at least in those of MS to SGB stars, which leads to generally small 3D NLTE abundance corrections (Asplund et al. 2003; Steffen et al. 2010). Nevertheless, it is still not entirely clear whether less sophisticated tools/approaches could be used in the case of lithium to circumvent the full 3D NLTE calculations, which are still computationally expensive and less accessible to wider astronomical community.

In this paper, we therefore investigate whether the 3D+NLTE and/or ⟨3D⟩ NLTE approaches could be used to substitute the full 3D NLTE methodology to obtain 3D NLTE abundances of lithium in the atmospheres of MS, TO, SGB, and RGB stars.

The tests are performed using state-of-the-art 3D hydrodynamical CO5BOLD and 1D hydrostatic LHD model atmospheres, computed at two different metallicities, [ M / H ] = 0.0 and −2.0. The CO5BOLD and LHD model atmospheres share identical atmospheric parameters, chemical composition, equation of state, and opacities, to make a strictly differential analysis of the role of 3D/1D and NLTE/LTE effects in the spectral line formation.

The paper is structured as follows: in Sect. 2 we describe the model atmospheres utilized in our study and outline the methodology of 3D/1D NLTE/LTE spectral synthesis computations. The results obtained using various spectral line synthesis approaches are presented and discussed in Sect. 3, while in Sect. 4 we provide a short summary of the most important results obtained in this work.

2. Methodology

2.1. Model atom of lithium

In this study, we use a model atom of Li i which was originally developed and tested by Cayrel et al. (2007) and Sbordone et al. (2010). For the purposes of the current work, the model atom was updated and now consists of 26 levels and 123 (96 of which are radiative) bound-bound transitions of Li i and the ground level of Li ii, with each level of Li i coupled to the continuum via bound-free transitions. (The ground state of Li ii in the current model atom is always in LTE, since lithium is mostly fully ionized throughout the model atmospheres studied in this work.) This renders the model atom complete up to the principal quantum number n = 6 and spectroscopic term , with additional energy levels up to n = 9 and term 2D (Fig. 1). Data concerning atomic energy levels and transitions (level energies and statistical weights; wavelengths and Einstein coefficients of the bound-bound transitions) were taken from the NIST database. We used electron collisional excitation and ionization rates from the quantum mechanical computations of Osorio et al. (2011) for the energy levels of up to 5s (2S). Elsewhere, collisional excitation by electrons for radiatively permitted transitions was accounted for by using the classical formula of van Regemorter (1962), while the formula of Seaton (1962) was used to compute collisional electron ionization rates. To account for the collisional excitation by hydrogen, we used collisional excitation rates computed by Barklem et al. (2003), while the classical formula of Drawin (in the formulation of Lambert 1993) was used for radiatively permitted transitions when no quantum mechanical data were available. Hydrogen H–Li charge transfer rates were taken from Barklem et al. (2003) for the atomic levels up to 4p inclusive. Bound-free transitions resulting from collisions with hydrogen were expected to be inefficient and thus were ignored. Photoionization cross sections were taken from TOPBASE (Cunto et al. 1993). No scaling of collisional rates was applied in the calculations of bound-free and bound-bound transitions. Information about the energy levels and bound-bound radiative transitions, included in the present version of the Li i model atom, are provided in Tables A.1 and A.2, respectively. Twenty-seven transitions in the model atom are purely collisional. Collisional radiatively-forbidden transitions involving Li i levels beyond 5s were not accounted for since reliable quantum-mechanical data for these transitions are not available. We note that the role of the omitted transitions between the higher levels is minor: when they are taken into account using the formula of Allen (1973), collision strength Ω = 1, the change in the estimated abundance (which directly applies to abundance corrections, too) is always less than 0.05 dex, with typical values being significantly smaller.

thumbnail Fig. 1

Model atom of lithium used in the 3D NLTE spectral line synthesis computations. The thick red line indicates the transition that corresponds to the lithium 670.8 nm resonance doublet. Other radiative bound-bound transitions are shown as grey lines.

2.2. 3D and 1D model atmospheres, spectral line synthesis, and abundance corrections

The 3D hydrodynamical models used in this work were taken from the CIFIST CO5BOLD model atmosphere grid (Ludwig et al. 2009). All simulation runs cover ~13480 convective turnover times, as measured by the Brunt-Vaisälä timescale (see Kučinskas et al. 2013, for the definition). The models were computed using solar-scaled elemental abundances from Asplund et al. (2005), with a constant enhancement of [α/ Fe] = + 0.4 in alpha-element abundances for models at [ M / H ] = −2.0. Monochromatic opacities were taken from the MARCS model atmosphere package (Gustafsson et al. 2008) and were grouped into five and six opacity bins (Ludwig et al. 2009) for models at [ M / H ] = 0.0 and −2.0, respectively. All simulations were performed under the assumption of LTE using a Cartesian model grid; the gravity vector has been assumed constant and antiparallel to z-axis throughout the model box, thus the effects of sphericity were ignored (for more details on the model calculations see Ludwig et al. 2009). For comparison, we also used 1D hydrostatic LHD model atmospheres (Caffau et al. 2008). It is important to note that both CO5BOLD and LHD models shared identical atmospheric parameters, opacities, chemical composition, and equation of state, which allowed us to make a strictly differential comparison of their predictions.

We used 3D hydrodynamical and 1D hydrostatic models of the representative RGB, SGB, TO, and MS stars, in each case computed for the same log g at two different effective temperatures and metallicities (see Table 1). This choice of model parameters allowed us to bracket the range of Teff, log g, and [ M / H ] typical of real stars observed in various Galactic populations. At the same time, this provides an opportunity to obtain a rough estimate of the size of 3D/1D and/or NLTE/LTE effects in real stellar atmospheres that are covered by our parameter range, by interpolating our results obtained at the bracketing values of stellar parameters. Atmospheric parameters of the 3D hydrodynamical CO5BOLD model atmospheres used in our study are provided in Table 1, with their positions in the log gTeff plane shown in Fig. 2. At each point in the log gTeff diagram, calculations were performed with two models of different metallicity, [ M / H ] = 0.0 and −2.0, to assess the differential effect of metallicity on the lithium spectral line formation.

Table 1

3D hydrodynamical CO5BOLD model atmospheres used in this work.

thumbnail Fig. 2

3D hydrodynamical CO5BOLD model atmospheres used in this study, plotted in the log gTeff plane (large filled circles). Several isochrones are shown as lines to indicate the approximate loci of stars in various Galactic populations.

An updated version of the NLTE3D code1 was used to compute the NLTE departure coefficients, bi(x,y,z), at each geometrical position (x,y,z) in the background model atmosphere for each energy level i, utilizing the model atom of Li i described in the previous section (bi(x,y,z) = ni(x,y,z)NLTE/ni(x,y,z)LTE, where ni(x,y,z)NLTE and ni(x,y,z)NLTE are population densities of a given level i in NLTE and LTE, respectively). NLTE3D computes the departure coefficients by solving the statistical equilibrium equations at each grid cell (x,y,z) using an accelerated Λ-iteration scheme. The line-blanketed mean continuum intensity, Jν(x,y,z), is computed using the BIG opacity distribution functions (ODFs) from the ATLAS9 model atmosphere package, at 600 frequency points, which spans the spectral range from 92.5 to 7300 nm, and accounts for coherent, isotropic scattering in the continuum. Jν was computed as a weighted average of the intensities obtained along the vertical direction, four azimuthal directions, and four inclination angles (i.e., 17 rays in total).

The departure coefficients were computed using three types of model atmospheres: (i) 3D hydrodynamical, calculated with the CO5BOLD code (Freytag et al. 2012); (ii) ⟨3D⟩ average atmospheres, obtained by taking a temporal and spatial average over the sequence of full 3D model structures; and (iii) 1D hydrostatic atmospheres, calculated using the LHD code (Caffau et al. 2008). To calculate ⟨3D⟩ average atmospheres, we averaged the fourth moment of temperature and the first moment of gas pressure over surfaces of constant Rosseland optical depth, according to the prescription provided in Steffen et al. (1995). For detailed reasoning behind the choice of the averaging scheme, see, e.g., Appendix A of Dobrovolskas et al. (2013). While this choice is fully adequate for the purposes of the current study, there are other options for constructing average ⟨3D⟩ model atmospheres, depending on the specific task (see, e.g., Magic et al. 2013). The average ⟨3D⟩ model is in fact a 1D model and, because it is obtained by temporal and horizontal averaging of the full 3D model sequence, it lacks information about the horizontal inhomogeneities that are present in the 3D model structures. Thus, the comparison of the line formation properties in full 3D and ⟨3D⟩ models allows us to assess the importance of horizontal inhomogeneities in the process of line formation.

In most astrophysical situations, the only spectral line available for lithium diagnostics in the UV–IR stellar spectra is the Li i 670.8 nm resonance doublet. We therefore focused on this single tracer by computing its synthetic 3D and 1D profiles with the Linfor3D spectral synthesis code, both in LTE and NLTE. The doublet structure of the Li i 670.8 nm resonance line was accounted for when computing the bound-bound radiative rates, while assuming identical NLTE departure coefficients for the two fine structure sublevels of the 2p level. To synthesize the line profile, the fine structure of the Li i 670.8 nm feature was taken into account. The analysis was done in the cases of (i) weak (equivalent width W = 0.5 pm); and (ii) strong (W = 8 pm) lines. These line strengths bracket the typical EWs of the Li i 670.8 nm line observed in stars with atmospheric parameters similar to those studied in this work. However, as we will see in Sect. 3, there is little sensitivity of the abundance corrections to the line strength for all model atmospheres studied here.

The interplay between 3D hydrodynamical and NLTE effects was studied with the help of abundance corrections (see, e.g., Caffau et al. 2011). The latter were defined as differences in Li abundance that would be determined from the Li i 670.8 nm line of a given strength using different model atmospheres, in NLTE and/or LTE. First, we computed the following abundance corrections: Δ3D NLTE − 1D LTE, Δ1D NLTE − 1D LTE, and Δ3D LTE − 1D LTE. Then, the latter two corrections were used to obtain the Δ3D + NLTE ≡ Δ1D NLTE − 1D LTE + Δ3D LTE − 1D LTE abundance correction, i.e., the total correction expected when the 1D NLTE abundances are corrected for 3D effects by adding the 3D–1D LTE abundance correction. Finally, we also calculated the Δ⟨ 3D ⟩ NLTE − 1D LTE abundance correction obtained by using the average ⟨3D⟩ and 1D model atmospheres. All obtained abundance corrections are listed in Table 2.

Finally, we note that in this work we frequently used 1D LTE abundance estimates as the reference point for computing and assessing various abundance corrections. This choice may be well justified if one is interested in knowing how much the lithium abundances, which were obtained using various more sophisticated approaches, would differ from those determined in 1D LTE. At the same time, this may offer a convenient way to “correct” the 1D LTE abundances for 3D and/or NLTE effects. Nevertheless, it is clearly the full 3D NLTE approach that provides the highest realism in modeling lithium spectral line formation in stellar atmospheres. As such, ideally the 3D NLTE methodology should be applied to obtain the most reliable lithium abundances from the measured equivalent widths of the Li i 670.8 nm lines (e.g., by using the 3D NLTE1D LTE abundance corrections).

thumbnail Fig. 3

Δ3D NLTE − 1D LTE (shown in all panels), Δ3D + NLTE − 1D LTE (top), Δ⟨ 3D ⟩ NLTE − 1D LTE (middle) and Δ1D NLTE − 1D LTE (bottom) abundance corrections for the lithium 670.8 nm line in the atmospheres of RGB, SGB, TO, and MS stars, at [ M / H ] = 0.0 and −2.0, in the order they appear in Table 2.

Table 2

Abundance corrections for the lithium 670.8 nm line in the atmospheres of RGB, SGB, TO, and MS stars, at [ M / H ] = 0.0 and −2.0.

thumbnail Fig. 4

Thermodynamic structures and Li 670.8 nm line formation properties in two model atmospheres of red giant stars (Teff = 4480 K, log g = 2.5), at [ M / H ] = 0.0 (left) and [ M / H ] = −2.0 (right). From top to bottom: (1) temperature profiles in the 3D hydrodynamical (probability density, grayscale), average ⟨3D⟩ (dashed red line), and 1D LHD (solid red line) model atmospheres; (2) horizontal temperature fluctuations in the 3D model (ΔTrms, solid black line) and differences between the temperature profiles of the ⟨3D⟩ and 1D LHD model atmospheres (blue dashed lines); (3–5) contribution functions of the Li 670.8 nm line (i.e., rates of the line equivalent width growth, dEW/ dlog τRoss, as a function of optical depth) in 3D NLTE/LTE (black), ⟨3D⟩ NLTE/LTE (red), and 1D NLTE/LTE (blue). For each metallicity, the lithium abundance is fixed such that W(1D LTE) = 0.5 pm.

thumbnail Fig. 5

Same as Fig. 4, but for TO stars.

3. Results and discussion

The results provided in Table 2 and Fig. 3 show that, in general, the full Δ3D NLTE − 1D LTE correction (Col. 5) is different from the Δ3D + NLTE (Δ1D NLTE − 1D LTE + Δ3D LTE − 1D LTE) correction (Col. 4) at all metallicities and in all model atmospheres studied here. We note, however, that while these differences are generally small at solar metallicity (MS #2 is a clear exception), they do indeed become significant at [ M / H ] = −2.0. In the latter case, the full 3D NLTE correction is small and, depending on the line strength, amounts to Δ3D NLTE − 1D LTE = −0.06··· + 0.12 dex. At the same time, the Δ3D + NLTE correction is always large and negative, and may reach, for example, to ~−0.18···−0.23 dex and ~−0.31···−0.50 dex for RGB and TO stars, respectively. On the other hand, the Δ⟨ 3D ⟩ NLTE − 1D LTE and the Δ1D NLTE − 1D LTE abundance corrections are in many cases very similar to the full Δ3D NLTE − 1D LTE correction: the differences are | Δ3D NLTE − 1D LTE − Δ⟨ 3D ⟩ NLTE − 1D LTE | ≤ 0.04 and | Δ3D NLTE − 1D LTE − Δ1D NLTE − 1D LTE | ≤ 0.08. Although these discrepancies are small, one should note that full Δ3D NLTE − 1D LTE corrections are small too, especially at [ M / H ] = −2.0. As a consequence, at lower metallicities the differences Δ3D NLTE − 1D LTE − Δ⟨ 3D ⟩ NLTE − 1D LTE and Δ3D NLTE − 1D LTE − Δ1D NLTE − 1D LTE may become comparable to (or even exceed) the size of the full Δ3D NLTE − 1D LTE corrections.

The behavior of different abundance corrections shown in Table 2 is relatively easy to understand. As seen from Figs. 4 and 5, the metal content plays a significant role in shaping the structure of the photospheric layers where the Li i line forms in RGB and TO stars. At solar metallicity, the mean temperature of the 3D models of RGB and TO stars is close to the 1D radiative equilibrium solution, and the amplitude of the horizontal temperature fluctuations in both cases is moderate, ΔTrms ≲ 350 K (, where the angled brackets indicate temporal and horizontal averaging on surfaces of equal optical depth, and T0 = ⟨ Tx,y,t is the depth-dependent average temperature). In the metal-poor TO model atmosphere, the mean temperature of the 3D model is significantly lower than that predicted in 1D (up to 700 K), and at the same time the horizontal temperature fluctuations are substantial, increasing with height to ΔTrms ≳ 800 K (Fig. 5, second panels from top). Assuming LTE, both effects lead to significant line strengthening in 3D with respect to 1D in the metal-poor case. This happens because line opacity is a very sensitive non-linear function of temperature. Therefore, even if temperature had a symmetric distribution around the mean value at a given optical depth, the contribution from low-T regions (e.g., inter-granular lanes) towards the total opacity would be more important than that from high-T regions (e.g., granules). As a consequence, the net increase of line opacity in the former would outweigh the net decrease of line opacity in the latter. This would make lines appear stronger in 3D than in ⟨3D⟩ or 1D and would result in negative Δ3D LTE − ⟨ 3D ⟩ LTE and Δ3D LTE − 1D LTE abundance corrections (for details see, e.g., Steffen & Holweger 2002; Kučinskas et al. 2013). The situation is similar in the atmosphere of a RGB star, although the line strengthening in 3D is slightly less in this case, because the difference between the temperature profiles of the average ⟨3D⟩ and 1D model atmospheres is small, and only the temperature fluctuations contribute to the line strengthening. Qualitatively, the effect is similar in all other model atmospheres studied here. Note, however, that the situation is different in NLTE (see below).

Comparison of the line contribution functions obtained in LTE (see Figs. 4 and 5) reveals their close similarity at [ M / H ] = 0.0. Consequently, one may infer that the effect of the horizontal temperature fluctuations is small at solar metallicity. However, the line contribution functions in 3D LTE and ⟨3D⟩ LTE are significantly different at [ M / H ] = −2.0, both in RGB and TO stars, which shows that the role of horizontal fluctuations now becomes substantial and dominates over the effect of the reduced mean temperature. In the metal-poor 3D hydrodynamical models, the fluctuations produce regions in the atmosphere where temperature drops significantly below that predicted by the 1D model. This leads to larger concentration of Li i in the 3D models, and thus, to stronger lines in 3D LTE (i.e., with respect to those in 1D LTE) and negative Δ3D LTE − 1D LTE abundance corrections. In NLTE (both 1D and 3D), however, atomic level population numbers are more sensitive to the average radiation field than to local temperature (e.g., Cayrel & Steffen 2000; Asplund et al. 2003). As a consequence, Li i gets significantly over-ionized in the outer atmosphere. Therefore, the mean concentration of Li i with respect to what would be expected in LTE is reduced. This leads to weaker spectral lines in 3D NLTE than in 3D LTE, and thus, significantly reduced Δ3D NLTE − 1D LTE abundance corrections which are similar in magnitude to Δ1D NLTE − 1D LTE corrections.

It is important to note that, in general, the choice of the microtubulence velocity used in the spectral line synthesis with the average ⟨3D⟩ and 1D model atmospheres may influence the corresponding abundance corrections. Fortunately, the resonance line of lithium observed in stars covered by our atmospheric parameter range is relatively weak, and so the influence of microturbulence on the line strength is, in fact, minor. Our tests show that the difference in the abundance correction obtained at microturbulent velocities of 1.0 km s-1 and 5.0 km s-1 is small, 0.05 dex. The difference in abundance corrections computed at more typical values of 1.0 and 2.0 km s-1 was less than 0.02 dex, irrespective of whether the line formation was treated in NLTE or LTE, even for lines as strong as W ≈ 8 pm. This illustrates that saturation of the stronger line of the Li i 670.8 nm doublet is not significant in the range of Ws investigated. However, in the Li-rich stars, this line may become saturated and thus experience an enhanced sensitivity to microturbulence and larger NLTE effects due to photon losses (e.g., Lind et al. 2009).

To check the importance of continuum scattering for the Li i statistical equilibrium solution, we computed an additional set of synthetic spectra where the departure coefficients were computed by treating scattering as true absorption. The effect of this change on the estimates of Li abundance in 3D, ⟨3D⟩, and 1D is relatively small2, as shown in Fig. 6: at solar metallicity, the differences in abundance increase from about 0.001 dex in the dwarfs to up to 0.012 dex in the giants, while a more pronounced effect of scattering is seen in metal-poor stars, where the abundance differences range from 0.005 dex in the dwarfs to 0.028 dex in the case of the metal-poor red giants. This is because, especially in the metal-poor giants, the opacity coefficients of continuum scattering and true absorption become comparable in magnitude. However, the effect on the abundance corrections (e.g., Δ3D NLTE − 1D NLTE) is much smaller and does not exceed 0.002 dex, except for the two metal-poor giants. Qualitatively, these findings are similar to those obtained by Hayek et al. (2011) who found that the impact of scattering on the formation of fictitious Fe i resonance lines in the atmosphere of red giant star is negligible at 500 nm ([ M / H ] = −2.0 and 0.0), both in 3D and 1D, with the resulting differences in the abundance corrections of <0.015 dex. In general, we find that 3D spectral line formation is slightly more sensitive to scattering effects compared to 1D, such that the abundance corrections are slightly more positive when continuum scattering is treated consistently.

thumbnail Fig. 6

Changes in A(Li) when different treatments of continuum scattering are applied in calculations of the departure coefficients (coherent scattering versus true absorption). Differently shaded areas represent different model atmospheres as they are listed in Table 2. Symbols on the left side of each shaded region represent spectral lines with W = 0.5 pm, while symbols on the right represent those with W = 8.0 pm.

General uncertainties involved in 3D NLTE spectral synthesis (excluding uncertainties in atomic data) can be illustrated by comparing Li abundance corrections, Δ3D NLTE − 1D LTE, computed in this work with those obtained using the interpolation formula from Sbordone et al. (2010, Eq. (B.1)). This formula is based on the 3D NLTE and 1D LTE spectral synthesis computations obtained using a grid of 3D hydrodynamical CO5BOLD and 1D hydrostatic LHD model atmospheres of dwarfs and subgiant stars (Teff ≈ 5470...6560 K, log g = 3.5...4.5, [ M / H ] = −2.0). The spectral synthesis computations of Sbordone et al. (2010) utilize departure coefficients computed with a simpler Li model atom and an older version of the NLTE3D code, treating continuum scattering as true absorption. Their formula is supposed to have an internal precision of ± 0.005 dex as long as extrapolation towards the values of equivalent widths and atmospheric parameters that are not covered in their simulations is avoided (see Figs. 6 and 7 in Sbordone et al. 2010). For this comparison we used models #12, 13, 14, and 15 from our grid, because their atmospheric parameters fall in the range covered by the simulations of Sbordone et al. (2010). We find that corrections obtained in this work differ from those computed using the formula of Sbordone et al. (2010) by −0.027 ± 0.012 dex on average for lines with W = 5 pm and by −0.055 ± 0.017 dex for lines with W = 8 pm, without a notable systematic dependence on atmospheric parameters. The larger disagreement in the case of stronger lines is not unexpected since, in the present calculations of the departure coefficients, we have taken the doublet fine structure of the 670.8 nm line into account, while the single component approximation used by Sbordone et al. (2010) is only valid for weak lines, and consequently Sbordone et al. (2010) does not recommend using their interpolation formula when saturation effects become important.

4. Summary and conclusions

We studied whether various simplified approaches (such as deriving abundances in 1D NLTE and applying to them 3D–1D LTE abundance corrections; or using ⟨3D⟩ models with the 1D NLTE analysis tools) could be used to substitute for the full 3D NLTE methodology to derive lithium abundances in the atmospheres of RGB, SGB, TO, and MS stars. Our results clearly show that the application of the combined Δ3D + NLTE correction may lead to erroneous results in the case of lithium at sub-solar metallicities with all model atmospheres studied here. On the other hand, both ⟨3D⟩ NLTE and 1D NLTE approaches provide a good approximation to the full 3D NLTE solution.

Our results suggest that, in most cases, ⟨3D⟩ NLTE or 1D NLTE modeling should be sufficient for deriving lithium abundances accurate to ~0.1 dex, i.e., with respect to those obtained in 3D NLTE. Nevertheless, since the full Δ3D NLTE − 1D LTE correction, in certain cases, is only larger than 0.1 dex, application of Δ1D NLTE − 1D LTE or Δ⟨ 3D ⟩ NLTE − 1D LTE corrections may lead to systematical bias that is comparable to the size of the Δ3D NLTE − 1D LTE correction itself. This may be obviously unacceptable when high accuracy/reliability in lithium abundances is needed. Ideally, the full 3D NLTE corrections should be preferred over any of the approximate approaches.

The use of ⟨3D⟩ NLTE or 1D NLTE approaches, moreover, cannot be extended in a straightforward fashion to other chemical elements: the importance of different physical processes leading to the departures from LTE may vary from one chemical element to another. As we have seen from the results obtained in our study (which are in line with the earlier findings of Asplund et al. 2003), the 3D–1D and NLTE–LTE abundance corrections generally tend to compensate each other, so that the total Δ3D NLTE − 1D LTE abundance correction is small. This, however, seems not to be the case with other elements where this kind of compensation does not occur. For example, in their study of NLTE Fe i spectral line formation based on a comparison with the observed spectra, Mashonkina et al. (2013) find the ⟨3D⟩ NLTE approach is only applicable to a limited subset of Fe i lines. Therefore, verification with the full 3D NLTE approach is needed on a case-by-case basis before any simplified approach may be applied with confidence to the abundance analysis of other elements in the atmospheres of real stars.


1

For details on the NLTE3D code see Sbordone et al. (2010), Prakapavičius et al. (2013).

2

Convergence limits for the computations of Jν and W were set at 0.01% and 0.001%, respectively, to ensure the precision of 0.001 dex or better in the computed abundance corrections.

Acknowledgments

We thank K. Lind for valuable comments and suggestions that helped to improve the paper. We are grateful to P. Bonifacio, V. Dobrovolskas, and D. Prakapavičius for their help during various stages of the paper’s preparation. This work was supported by grants from the Research Council of Lithuania (MIP-065/2013) and the bilateral French-Lithuanian program Gilibert (TAP LZ 06/2013, Research Council of Lithuania; 28471NE, Campus France). H.G.L. and M.S. acknowledge funding from the Research Council of Lithuania for research visits to Vilnius. H.G.L. acknowledges financial support by the Sonderforschungsbereich SFB 881 “The Milky Way System” (subprojects A4, A5) of the German Research Foundation (DFG). E.C. is grateful to the FONDATION MERAC for funding her fellowship. Some model computations were performed using resources at the High Performance Computing Center, HPC Sauletekis, of the Faculty of Physics, Vilnius University.

References

  1. Allen, C. W. 1973, Astrophysical quantities (London: University of London, Athlone Press) [Google Scholar]
  2. Amarsi, A. M., Asplund, M., Collet, R., & Leenaarts, J. 2015, MNRAS, 454, L11 [NASA ADS] [CrossRef] [Google Scholar]
  3. Asplund, M. 2005, ARA&A, 43, 481 [NASA ADS] [CrossRef] [Google Scholar]
  4. Asplund, M., Nordlund, Å., Trampedach, R., & Stein, R. F. 1999, A&A, 346, L11 [Google Scholar]
  5. Asplund, M., Carlsson, M., & Botnen, A. V. 2003, A&A, 399, L31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Asplund, M., Grevesse, N., Sauval, A. J., Allen de Prieto, C., & Kiselman, D. 2004, A&A, 417, 751 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Asplund, M., Grevesse, N., Sauval, A. J., Allen de Prieto, C., & Kiselman, D. 2005, A&A, 435, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Barklem, P. S., Belyaev, A. K., & Asplund, M. 2003, A&A, 409, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Beeck, B., Cameron, R. H., Reiners, A., & Schüssler, M. 2013, A&A, 558, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Behara, N. T., Bonifacio, P., Ludwig, H.-G., et al. 2010, A&A, 513, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Caffau, E., Ludwig, H.-G., Steffen, M., et al. 2008, A&A, 488, 1031 [CrossRef] [EDP Sciences] [Google Scholar]
  12. Caffau, E., Ludwig, H.-G., Steffen, M., Freytag, B., & Bonifacio, P. 2011, Sol. Phys., 268, 255 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cayrel, R., & Steffen, M. 2000, Proc. IAU Symp., 198, 437 [Google Scholar]
  14. Cayrel, R., Steffen, M., Chand, H., et al. 2007, A&A, 473, L37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Collet, R., Asplund, M., & Trampedach, R. 2007, A&A, 469, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Collet, R., Asplund, M., & Nissen, P. E. 2009, PASA, 26, 330 [NASA ADS] [CrossRef] [Google Scholar]
  17. Cunto, W., Mendoza, C., Ochsenbein, F., & Zeipen, C. J. 1993, A&A, 275, L5 [NASA ADS] [Google Scholar]
  18. Dobrovolskas, V., Kučinskas, A., Steffen, M., et al. 2013, A&A, 559, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Freytag, B., Steffen, M., Ludwig, H.-G., et al. 2012, J. Comp. Phys., 231, 919 [Google Scholar]
  20. Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Hayek, W., Asplund, M., Collet, R., & Nordlund, A. 2011, A&A, 529, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Holzreuter, R., & Solanki, S. K. 2013, A&A, 558, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Kučinskas, A., Steffen, M., Ludwig, H.-G., et al. 2013, A&A, 549, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Lambert, D. L. 1993, Phys. Scr., 47, 186 [NASA ADS] [CrossRef] [Google Scholar]
  25. Lind, K., Asplund, M., & Barklem, P. S. 2009, A&A, 503, 541 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Lind, K., Asplund, M., Collet, R., & Melendez, J. 2012, Mem. Soc. Astron. It. Suppl., 22, 124 [Google Scholar]
  27. Lind, K., Melendez, J., Asplund, M., Collet, R., & Magic, Z. 2013, A&A, 554, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Ludwig, H.-G., Caffau, E., Steffen, M., et al. 2009, Mem. Soc. Astron. Italiana, 80, 711 [Google Scholar]
  29. Magic, Z., Collet, R., Hayek, W., & Asplund, M. 2013, A&A, 560, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Mashonkina, L., Ludwig, H.-G., Korn, A., Sitnova, T., & Caffau, E. 2013, Mem. Soc. Astron. It. Suppl., 24, 120 [Google Scholar]
  31. Mucciarelli, A., Salaris, M., & Bonifacio, P. 2012, MNRAS, 419, 2195 [NASA ADS] [CrossRef] [Google Scholar]
  32. Osorio, Y., Barklem, P. S., Lind, K., & Asplund, M. 2011, A&A, 529, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2006, ApJ, 642, 797 [NASA ADS] [CrossRef] [Google Scholar]
  34. Prakapavičius, D., Steffen, M., Kučinskas, A., et al. 2013, Mem. Soc. Astron. It. Suppl., 24, 111 [Google Scholar]
  35. Sbordone, L., Bonifacio, P., Caffau, E., et al. 2010, A&A, 522, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Seaton, M. J. 1962, Proc. Phys. Soc., 79, 1105 [NASA ADS] [CrossRef] [Google Scholar]
  37. Short, C. I., & Hauschildt, P. H. 2006, ApJ, 641, 494 [NASA ADS] [CrossRef] [Google Scholar]
  38. Steffen, M., & Holweger, H. 2002, A&A, 387, 258 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Steffen, M., Ludwig, H.-G., & Freytag, B. 1995, A&A, 300, 473 [Google Scholar]
  40. Steffen, M., Cayrel, R., Bonifacio, P., Ludwig, H.-G., & Caffau, E. 2010, Proc. IAU Symp., 268, 215 [NASA ADS] [Google Scholar]
  41. Steffen, M., Cayrel, R., Caffau, E., et al. 2012, Mem. Soc. Astron. It., 22, 152 [Google Scholar]
  42. Steffen, M., Prakapavičius, D., Caffau, E., et al. 2015, A&A, 583, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. van Regemorter, H. 1962, ApJ, 136, 906 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Model atom of lithium

The data we used to update the model atom of Li i are provided in Tables A.1 and A.2. We list atomic parameters of the energy levels and/or transitions, as well as the sources from which the data were taken.

Table A.1

Energy levels of the Li i model atom.

Table A.2

Radiative bound-bound transitions of the Li i model atom.

All Tables

Table 1

3D hydrodynamical CO5BOLD model atmospheres used in this work.

Table 2

Abundance corrections for the lithium 670.8 nm line in the atmospheres of RGB, SGB, TO, and MS stars, at [ M / H ] = 0.0 and −2.0.

Table A.1

Energy levels of the Li i model atom.

Table A.2

Radiative bound-bound transitions of the Li i model atom.

All Figures

thumbnail Fig. 1

Model atom of lithium used in the 3D NLTE spectral line synthesis computations. The thick red line indicates the transition that corresponds to the lithium 670.8 nm resonance doublet. Other radiative bound-bound transitions are shown as grey lines.

In the text
thumbnail Fig. 2

3D hydrodynamical CO5BOLD model atmospheres used in this study, plotted in the log gTeff plane (large filled circles). Several isochrones are shown as lines to indicate the approximate loci of stars in various Galactic populations.

In the text
thumbnail Fig. 3

Δ3D NLTE − 1D LTE (shown in all panels), Δ3D + NLTE − 1D LTE (top), Δ⟨ 3D ⟩ NLTE − 1D LTE (middle) and Δ1D NLTE − 1D LTE (bottom) abundance corrections for the lithium 670.8 nm line in the atmospheres of RGB, SGB, TO, and MS stars, at [ M / H ] = 0.0 and −2.0, in the order they appear in Table 2.

In the text
thumbnail Fig. 4

Thermodynamic structures and Li 670.8 nm line formation properties in two model atmospheres of red giant stars (Teff = 4480 K, log g = 2.5), at [ M / H ] = 0.0 (left) and [ M / H ] = −2.0 (right). From top to bottom: (1) temperature profiles in the 3D hydrodynamical (probability density, grayscale), average ⟨3D⟩ (dashed red line), and 1D LHD (solid red line) model atmospheres; (2) horizontal temperature fluctuations in the 3D model (ΔTrms, solid black line) and differences between the temperature profiles of the ⟨3D⟩ and 1D LHD model atmospheres (blue dashed lines); (3–5) contribution functions of the Li 670.8 nm line (i.e., rates of the line equivalent width growth, dEW/ dlog τRoss, as a function of optical depth) in 3D NLTE/LTE (black), ⟨3D⟩ NLTE/LTE (red), and 1D NLTE/LTE (blue). For each metallicity, the lithium abundance is fixed such that W(1D LTE) = 0.5 pm.

In the text
thumbnail Fig. 5

Same as Fig. 4, but for TO stars.

In the text
thumbnail Fig. 6

Changes in A(Li) when different treatments of continuum scattering are applied in calculations of the departure coefficients (coherent scattering versus true absorption). Differently shaded areas represent different model atmospheres as they are listed in Table 2. Symbols on the left side of each shaded region represent spectral lines with W = 0.5 pm, while symbols on the right represent those with W = 8.0 pm.

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.