Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

Nitrogen is an important element in various fields of stellar and Galactic astronomy, and the solar nitrogen abundance is crucial as a yardstick for comparing different objects in the cosmos. In order to obtain a precise and accurate value for this abundance, we carried out N i line formation calculations in a 3D radiative-hydrodynamic STAGGER model solar atmosphere, in full 3D non-local thermodynamic equilibrium (non-LTE), using a model atom that includes physically-motivated descriptions for the inelastic collisions of N i with free electrons and with neutral hydrogen. We selected five N i lines of high excitation energy to study in detail, based on their strengths and on their being relatively free of blends. We found that these lines are slightly strengthened from non-LTE photon losses and from 3D granulation effects, resulting in negative abundance corrections of around $-0.01$ dex and $-0.04$ dex respectively. Our advocated solar nitrogen abundance is $\log\epsilon_{\mathrm{N}} = 7.77$, with the systematic $1\sigma$ uncertainty estimated to be $0.05$ dex. This result is consistent with earlier studies after correcting for differences in line selections and equivalent widths.

The abundances of carbon, oxygen, and iron in late-type stars are important parameters in exoplanetary and stellar physics, as well as key tracers of stellar populations and Galactic chemical evolution. We carried out three-dimensional (3D) non-LTE radiative transfer calculations for CI and OI, and 3D LTE radiative transfer calculations for FeII, across the STAGGER-grid of 3D hydrodynamic model atmospheres. The absolute 3D non-LTE versus 1D LTE abundance corrections can be as severe as $-0.3$ dex for CI lines in low-metallicity F dwarfs, and $-0.6$ dex for OI lines in high-metallicity F dwarfs. The 3D LTE versus 1D LTE abundance corrections for FeII lines are less severe, typically less than $+0.15$ dex. We used the corrections in a re-analysis of carbon, oxygen, and iron in $187$ F and G dwarfs in the Galactic disk and halo. Applying the differential 3D non-LTE corrections to 1D LTE abundances visibly reduces the scatter in the abundance plots. The thick disk and high-$\alpha$ halo population rise in carbon and oxygen with decreasing metallicity, and reach a maximum of [C/Fe]$\approx0.2$ and a plateau of [O/Fe]$\approx0.6$ at [Fe/H]$\approx-1.0$. The low-$\alpha$ halo population is qualitatively similar, albeit offset towards lower metallicities and with larger scatter. Nevertheless, these populations overlap in the [C/O] versus [O/H] plane, decreasing to a plateau of [C/O]$\approx-0.6$ below [O/H]$\approx-1.0$. In the thin-disk, stars having confirmed planet detections tend to have higher values of C/O at given [O/H]; this potential signature of planet formation is only apparent after applying the abundance corrections to the 1D LTE results. Our grids of line-by-line abundance corrections are publicly available and can be readily used to improve the accuracy of spectroscopic analyses of late-type stars.

Older GCE models predict [K/Fe] ratios as much as 1 dex lower than those inferred from stellar observations. Abundances of potassium are mainly based on analyses of the 7698 $\AA$ resonance line, and the discrepancy between models and observations is in part caused by the LTE assumption. We study the statistical equilibrium of KI, focusing on the non-LTE effects on the $7698 \ \AA$ line. We aim to determine how non-LTE abundances of K can improve the analysis of its chemical evolution, and help to constrain the yields of models. We construct a model atom that employs the most up-to-date data. In particular, we calculate and present inelastic e+K collisional excitation cross-sections from the convergent close-coupling and the $B$-Spline $R$-matrix methods, and H+K collisions from the two-electron model. We constructed a fine grid of non-LTE abundance corrections that span $4000<\teff / \rm{K}<8000$, $0.50<\lgg<5.00$, $-5.00<\feh<+0.50$, and applied the corrections to abundances from the literature. In concordance with previous studies, we find severe non-LTE effects in the $7698 \ \AA$ line, which is stronger in non-LTE with abundance corrections that can reach $\sim-0.7\,\dex$. We explore the effects of atmospheric inhomogeneity by computing a full 3D non-LTE stellar spectrum of KI for a test star. We find that 3D is necessary to predict a correct shape of the resonance 7698 $\AA$ line, but the line strength is similar to that found in 1D non-LTE. Our non-LTE abundance corrections reduce the scatter and change the cosmic trends of literature K abundances. In the regime [Fe/H]$\lesssim-1.0$ the non-LTE abundances show a good agreement with the GCE model with yields from rotating massive stars. The reduced scatter of the non-LTE corrected abundances of a sample of solar twins shows that line-by-line differential analysis techniques cannot fully compensate for systematic modelling errors.

Low-energy inelastic collisions with neutral hydrogen atoms are important processes in stellar atmospheres, and a persistent source of uncertainty in non-LTE modelling of stellar spectra. We have calculated and studied excitation and charge transfer of C i and of N i due to such collisions. We used a previously presented method that is based on an asymptotic two-electron linear combination of atomic orbitals (LCAO) model of ionic-covalent interactions for the adiabatic potential energies, combined with the multichannel Landau-Zener model for the collision dynamics. We find that charge transfer processes typically lead to much larger rate coefficients than excitation processes do, consistent with studies of other atomic species. Two-electron processes were considered and lead to non-zero rate coefficients that can potentially impact statistical equilibrium calculations. However, they were included in the model in an approximate way, via an estimate for the two-electron coupling that was presented earlier in the literature: the validity of these data should be checked in a future work.

Oxygen abundance measurements are important for understanding stellar structure and evolution. Measured in Cepheids, they further provide clues on the metallicity gradient and chemo-dynamical evolution in the Galaxy. However, most of the abundance analyses of Cepheids to date have been based on one-dimensional (1D) hydrostatic model atmospheres. Here, we test the validity of this approach for the key oxygen abundance diagnostic, the \ion{O}{I} $777\,\mathrm{nm}$~triplet lines. We carry out 2D non-LTE radiative transfer clculations across two different 2D radiation hydrodynamics simulations of Cepheid atmospheres, having stellar parameters of $T_\mathrm{eff}= 5600$ K, solar chemical compositions, and $\log\,g= 1.5$ and $2.0$, corresponding to pulsation periods of 9 and 3 days, respectively. We find that the 2D non-LTE versus 1D LTE abundance differences range from $-1.0$~dex to $-0.25$~dex depending on pulsational phase. The 2D non-LTE versus 1D non-LTE abundance differences range from $-0.2$~dex to $0.8$~dex. The abundance differences are smallest when the Cepheid atmospheres are closest to hydrostatic equilibrium, corresponding to phases of around $0.3$ to $0.8$, and we recommend these phases for observers deriving the oxygen abundance from \ion{O}{I} $777\,\mathrm{nm}$ triplet with 1D hydrostatic models.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Carbon abundances in late-type stars are important in a variety of astrophysical contexts. However C i lines, one of the main abundance diagnostics, are sensitive to departures from local thermodynamic equilibrium (LTE). We present a model atom for non-LTE analyses of C i lines, that uses a new, physically-motivated recipe for the rates of neutral hydrogen impact excitation. We analyse C i lines in the solar spectrum, employing a three-dimensional (3D) hydrodynamic model solar atmosphere and 3D non-LTE radiative transfer. We find negative non-LTE abundance corrections for C i lines in the solar photosphere, in accordance with previous studies, reaching up to around 0.1 dex in the disk-integrated flux. We also present the first fully consistent 3D non-LTE solar carbon abundance determination: we infer log $\epsilon_{\text{C}}$ = $8.44\pm0.02$, in good agreement with the current standard value. Our models reproduce the observed solar centre-to-limb variations of various C i lines, without any adjustments to the rates of neutral hydrogen impact excitation, suggesting that the proposed recipe may be a solution to the long-standing problem of how to reliably model inelastic collisions with neutral hydrogen in late-type stellar atmospheres.

Carbon and oxygen are key tracers of the Galactic chemical evolution; in particular, a reported upturn in [C/O] towards decreasing [O/H] in metal-poor halo stars could be a signature of nucleosynthesis by massive Population III stars. We reanalyse carbon, oxygen, and iron abundances in thirty-nine metal-poor turn-off stars. For the first time, we take into account three-dimensional (3D) hydrodynamic effects together with departures from local thermodynamic equilibrium (LTE) when determining both the stellar parameters and the elemental abundances, by deriving effective temperatures from 3D non-LTE H$\beta$ profiles, surface gravities from Gaia parallaxes, iron abundances from 3D LTE Feii equivalent widths, and carbon and oxygen abundances from 3D non-LTE Ci and Oi equivalent widths. We find that [C/Fe] stays flat with [Fe/H], whereas [O/Fe] increases linearly up to $0.75$ dex with decreasing [Fe/H] down to $-3.0$ dex. As such [C/O] monotonically decreases towards decreasing [O/H], in contrast to previous findings, mainly by virtue of less severe non-LTE effects for Oi at low [Fe/H] with our improved calculations.

The OI 777 nm triplet is a key diagnostic of oxygen abundances in the atmospheres of FGK-type stars; however it is sensitive to departures from local thermodynamic equilibrium (LTE). The accuracy of non-LTE line formation calculations has hitherto been limited by errors in the inelastic O+H collisional rate coefficients: several recent studies have used the so-called Drawin recipe, albeit with a correction factor $\mathrm{S_{H}}$ that is calibrated to the solar centre-to-limb variation of the triplet. We present a new model oxygen atom that incorporates inelastic O+H collisional rate coefficients using an asymptotic two-electron model based on linear combinations of atomic orbitals, combined with a free electron model, based on the impulse approximation. Using a 3D hydrodynamic stagger model solar atmosphere and 3D non-LTE line formation calculations, we demonstrate that this physically-motivated approach is able to reproduce the solar centre-to-limb variation of the triplet to 0.02 dex, without any calibration of the inelastic collisional rate coefficients or other free parameters. We infer $\log\epsilon_{\mathrm{O}}=8.69\pm0.03$ from the triplet alone, strengthening the case for a low solar oxygen abundance.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

Several recent spectroscopic investigations have presented conflicting results on the existence of Na-rich asymptotic giant branch (AGB) stars in the Galactic globular cluster M4 (NGC6121). The studies disagree on whether or not Na-rich red giant branch (RGB) stars evolve to the AGB. For a sample of previously published HER- MES/AAT AGB and RGB stellar spectra we present a re-analysis of O, Na, and Fe abundances, and a new analysis of Mg and Al abundances; we also present CN band strengths for this sample, derived from low-resolution AAOmega spectra. Following a detailed literature comparison, we find that the AGB samples of all studies consistently show lower abundances of Na and Al, and are weaker in CN, than RGB stars in the cluster. This is similar to recent observations of AGB stars in NGC 6752 and M 62. In an attempt to explain this result, we present new theoretical stellar evolutionary models for M 4; however, these predict that all stars, including Na-rich RGB stars, evolve onto the AGB. We test the robustness of our abundance results using a variety of atmospheric models and spectroscopic methods; however, we do not find evidence that systematic modelling uncertainties can explain the apparent lack of Na- rich AGB stars in M4. We conclude that an unexplained, but robust, discordance between observations and theory remains for the AGB stars in M 4.

Hydrogen Balmer lines are commonly used as spectroscopic effective temperature diagnostics of late-type stars. However, the absolute accuracy of classical methods that are based on one-dimensional (1D) hydrostatic model atmospheres and local thermodynamic equilibrium (LTE) is still unclear. To investigate this, we carry out 3D non-LTE calculations for the Balmer lines, performed, for the first time, over an extensive grid of 3D hydrodynamic STAGGER model atmospheres. For H$\alpha$, H$\beta$, and H$\gamma$, we find significant 1D non-LTE versus 3D non-LTE differences (3D effects): the outer wings tend to be stronger in 3D models, particularly for H$\gamma$, while the inner wings can be weaker in 3D models, particularly for H$\alpha$. For H$\alpha$, we also find significant 3D LTE versus 3D non-LTE differences (non-LTE effects): in warmer stars ($T_{\text{eff}}\approx6500$K) the inner wings tend to be weaker in non-LTE models, while at lower effective temperatures ($T_{\text{eff}}\approx4500$K) the inner wings can be stronger in non-LTE models; the non-LTE effects are more severe at lower metallicities. We test our 3D non-LTE models against observations of well-studied benchmark stars. For the Sun, we infer concordant effective temperatures from H$\alpha$, H$\beta$, and H$\gamma$; however the value is too low by around 50K which could signal residual modelling shortcomings. For other benchmark stars, our 3D non-LTE models generally reproduce the effective temperatures to within $1\sigma$ uncertainties. For H$\alpha$, the absolute 3D effects and non-LTE effects can separately reach around 100K, in terms of inferred effective temperatures. For metal-poor turn-off stars, 1D LTE models of H$\alpha$ can underestimate effective temperatures by around 150K. Our 3D non-LTE model spectra are publicly available, and can be used for more reliable spectroscopic effective temperature determinations.

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope.

We determine Mg abundances in 6 Gaia benchmark stars using theoretical one-dimensional (1D) hydrostatic model atmospheres, as well as temporally- and spatially-averaged 3D model atmospheres (<3D>). The stars cover a range of Teff from 4700 to 6500 K, log g from 1.6 to 4.4 dex, and [Fe/H] from -3.0 dex to solar. Spectrum synthesis calculations are performed in local thermodynamic equilibrium (LTE) and in non-LTE (NLTE) using the oscillator strengths recently published by Pehlivan Rhodin et al. We find that: a) Mg abundances determined from the infrared spectra are as accurate as the optical diagnostics, b) the NLTE effects on Mg I line strengths and abundances in this sample of stars are minor (although for a few Mg I lines the NLTE effects on abundance exceed 0.6 dex in <3D> and 0.1 dex in 1D, c) the solar Mg abundance is 7.56 +/- 0.05 dex (total error), in the excellent agreement with the Mg abundance measured in CI chondritic meteorites, d) the 1D NLTE and <3D> NLTE approach can be used with confidence to analyse optical Mg I lines in spectra of dwarfs and sub-giants, but for red giants the Mg I 5711 A line should be preferred, e) low-excitation Mg I lines are sensitive to the atmospheric structure; for these lines, LTE calculations with <3D> models lead to significant systematic abundance errors. The methods developed in this work will be used to study Mg abundances of a large sample of stars in the next paper in the series.

Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral hydrogen; the latter effectively remove a free parameter that has hampered all previous line formation studies of Fe in non-local thermodynamic equilibrium (NLTE). For the first time, we use realistic 3D NLTE calculations of Fe for a quantitative comparison to solar observations. We confront our theoretical line profiles with observations taken at different viewing angles across the solar disk with the Swedish 1-m Solar Telescope. We find that 3D modelling well reproduces the observed centre-to-limb behaviour of spectral lines overall, but highlight aspects that may require further work, especially cross-sections for inelastic collisions with electrons. Our inferred solar iron abundance is log(eps(Fe))=7.48+-0.04.

We present three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer calculations for silicon in the solar photosphere, using an extensive model atom that includes recent, realistic neutral hydrogen collisional cross-sections. We find that photon losses in the SiI lines give rise to slightly negative non-LTE abundance corrections of the order -0.01 dex. We infer a 3D non-LTE based solar silicon abundance of 7.51 dex. With silicon commonly chosen to be the anchor between the photospheric and meteoritic abundances, we find that the meteoritic abundance scale remains unchanged compared with the Asplund et al. (2009) and Lodders et al. (2009) results.

Models of star formation in the early universe depend on the details of accretion, fragmentation and radiative feedback. Different simulations predict different initial mass functions of the first stars, ranging from predominantly low mass (0.1-10 Msol), to massive (10-100 Msol), or even supermassive (100-1000 Msol). The mass distribution of the first stars should lead to unique chemical imprints on the low-mass second and later generation metal-poor stars still in existence. The chemical composition of SMSS0313-6708, which has the lowest abundances of Ca and Fe of any star known, indicates it was enriched by a single massive supernova. However, even weak spectral lines may be affected by strong 3D and NLTE effects in metal-poor stars. If these effects are ignored or treated incorrectly, errors in the inferred abundances may significantly bias the inferred properties of the polluting supernovae. We redetermine the chemical composition of SMSS0313-6708 using 3D NLTE radiative transfer to obtain accurate abundances for Li, Na, Mg, Al, Ca and Fe. The model atoms employ realistic collisional rates, with no calibrated free parameters. We find significantly higher abundances in 3D NLTE than 1D LTE by 0.8 dex for Fe, and 0.5 dex for Mg, Al and Ca, while Li and Na are unaffected to within 0.03 dex. In particular, our upper limit for [Fe/H] is now a factor ten larger, at [Fe/H] < -6.53 (3 sigma), than previous estimates based on <3D> NLTE (i.e. using averaged 3D models). This higher estimate is due to a conservative upper limit estimation, updated NLTE data, and 3D-<3D> NLTE differences, all of which lead to a higher abundance determination. We find that the revised chemical composition of SMSS0313-6708 matches supernova yields for massive progenitors of 20-60 Msol exploding with low energies (1-2 x 10^51 erg), as well as progenitors of 10 Msol with very low explosion energies (< 10^51 erg).

As one of the most important elements in astronomy, iron abundance determinations need to be as accurate as possible. We investigate the accuracy of spectroscopic iron abundance analyses using archetypal metal-poor stars. We perform detailed 3D non-LTE radiative transfer calculations based on 3D hydrodynamic Stagger model atmospheres, and employ a new model atom that includes new quantum-mechanical neutral hydrogen collisional rate coefficients. With the exception of the red giant HD122563, we find that the 3D non-LTE models achieve Fe i/Fe ii excitation and ionization balance as well as not having any trends with equivalent width to within modelling uncertainties of 0.05 dex, all without having to invoke any microturbulent broadening; for HD122563 we predict that the current best parallax-based surface gravity is over-estimated by 0.5 dex. Using a 3D non-LTE analysis, we infer iron abundances from the 3D model atmospheres that are roughly 0.1 dex higher than corresponding abundances from 1D MARCS model atmospheres; these differences go in the same direction as the non-LTE effects themselves. We make available grids of departure coefficients, equivalent widths and abundance corrections, calculated on 1D MARCS model atmospheres and horizontally- and temporally-averaged 3D Stagger model atmospheres.

The OI 777 nm lines are among the most commonly used diagnostics for the oxygen abundances in the atmospheres of FGK-type stars. However, they form in conditions that are far from local thermodynamic equilibrium (LTE). We explore the departures from LTE of atomic oxygen, and their impact on OI lines, across the Stagger-grid of three-dimensional hydrodynamic model atmospheres. For the OI 777 nm triplet we find significant departures from LTE. These departures are larger in stars with larger effective temperatures, smaller surface gravities, and larger oxygen abundances. We present grids of predicted 3D non-LTE based equivalent widths for the OI616nm, [OI] 630 nm, [OI] 636 nm, and OI 777 nm lines, as well as abundance corrections to 1D LTE based results.

We revisit the Galactic chemical evolution of oxygen, addressing the systematic errors inherent in classical determinations of the oxygen abundance that arise from the use of one dimensional hydrostatic (1D) model atmospheres and from the assumption of local thermodynamic equilibrium (LTE). We perform detailed 3D non-LTE radiative transfer calculations for atomic oxygen lines across a grid of 3D hydrodynamic stag- ger model atmospheres for dwarfs and subgiants. We apply our grid of predicted line strengths of the [OI] 630 nm and OI 777 nm lines using accurate stellar parameters from the literature. We infer a steep decay in [O/Fe] for [Fe/H] $\gtrsim$ -1.0, a plateau [O/Fe] $\approx$ 0.5 down to [Fe/H] $\approx$ -2.5 and an increasing trend for [Fe/H] $\lesssim$ -2.5. Our 3D non-LTE calculations yield overall concordant results from the two oxygen abundance diagnostics.

Line contribution functions are useful diagnostics for studying spectral line formation in stellar atmospheres. I derive an expression for the contribution function to the abso- lute flux depression that emerges from three-dimensional box-in-a-star model stellar atmospheres. I illustrate the result by comparing the local thermodynamic equilibrium (LTE) spectral line formation of the high-excitation permitted OI777nm lines with the non-LTE case.