We present an abundance analysis of the low-metallicity benchmark red giant star HD 122563 based on realistic, state-of-the-art, high-resolution, three-dimensional (3D) model stellar atmospheres including non-grey radiative transfer through opacity binning with four, twelve, and 48 bins. The 48-bin 3D simulation reaches temperatures lower by ~ 300 - 500 K than the corresponding 1D model in the upper atmosphere. Small variations in the opacity binning, adopted line opacities, or chemical mixture can cool the photospheric layers by a further ~ 100 - 300 K and alter the effective temperature by ~ 100 K. A 3D local thermodynamic equilibrium (LTE) spectroscopic analysis of Fe I and Fe II lines gives discrepant results in terms of derived Fe abundance, which we ascribe to non-LTE effects and systematic errors on the stellar parameters. We also determine C, N, and O abundances by simultaneously fitting CH, OH, NH, and CN molecular bands and lines in the ultraviolet, visible, and infrared. We find a small positive 3D-1D abundance correction for carbon (+0.03 dex) and negative ones for nitrogen (-0.07 dex) and oxygen (-0.34 dex). From the analysis of the [O I] line at 6300.3 {\AA}, we derive a significantly higher oxygen abundance than from molecular lines (+0.46 dex in 3D and +0.15 dex in 1D). We rule out important OH photodissociation effects as possible explanation for the discrepancy and note that lowering the surface gravity would reduce the oxygen abundance difference between molecular and atomic indicators.

The composition of the Sun is an essential piece of reference data for astronomy, cosmology, astroparticle, space and geo-physics. This article, dealing with the intermediate-mass elements Na to Ca, is the first in a series describing the comprehensive re-determination of the solar composition. In this series we severely scrutinise all ingredients of the analysis across all elements, to obtain the most accurate, homogeneous and reliable results possible. We employ a highly realistic 3D hydrodynamic solar photospheric model, which has successfully passed an arsenal of observational diagnostics. To quantify systematic errors, we repeat the analysis with three 1D hydrostatic model atmospheres (MARCS, MISS and Holweger & M\"{u}ller 1974) and a horizontally and temporally-averaged version of the 3D model ($\langle$3D$\rangle$). We account for departures from LTE wherever possible. We have scoured the literature for the best transition probabilities, partition functions, hyperfine and other data, and stringently checked all observed profiles for blends. Our final 3D+NLTE abundances are: $\log\epsilon_{\mathrm{Na}}=6.21\pm0.04$, $\log\epsilon_{\mathrm{Mg}}=7.59\pm0.04$, $\log\epsilon_{\mathrm{Al}}=6.43\pm0.04$, $\log\epsilon_{\mathrm{Si}}=7.51\pm0.03$, $\log\epsilon_{\mathrm{P}}=5.41\pm0.03$, $\log \epsilon_{\mathrm{S}}=7.13\pm0.03$, $\log\epsilon_{\mathrm{K}}=5.04\pm0.05$ and $\log\epsilon_{\mathrm{Ca}}=6.32\pm0.03$. The uncertainties include both statistical and systematic errors. Our results are systematically smaller than most previous ones with the 1D semi-empirical Holweger & M\"uller model. The $\langle$3D$\rangle$ model returns abundances very similar to the full 3D calculations. This analysis provides a complete description and a slight update of the Na to Ca results presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), with full details of all lines and input data.

The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution.

We study the implications of averaging methods with different reference depth scales for 3D hydrodynamical model atmospheres computed with the Stagger-code. The temporally and spatially averaged (hereafter denoted as <3D>) models are explored in the light of local thermodynamic equilibrium (LTE) spectral line formation by comparing spectrum calculations using full 3D atmosphere structures with those from <3D> averages. We explore methods for computing mean <3D> stratifications from the Stagger-grid time-dependent 3D radiative hydro- dynamical atmosphere models by considering four different reference depth scales (geometrical depth, column-mass density, and two optical depth scales). Furthermore, we investigate the influence of alternative averages (logarithmic or enforced hydrostatic equilibrium, flux-weighted temperatures). For the line formation we compute curves of growth for Fe i and Fe ii lines in LTE . The resulting <3D> stratifications for the four reference depth scales can be considerably different. We find typically that in the upper atmosphere and in the superadiabatic region just below the optical surface, where the temperature and density fluctuations are highest, the differences become considerable and increase for higher Teff, lower logg, and lower [Fe/H]. The differential comparison of spectral line formation shows distinctive differences depending on which <3D> model is applied. The averages over layers of constant column-mass density yield the best mean <3D> representation for LTE line formation, while the averages on layers at constant geometrical height are the least appropriate. Unexpectedly, the usually preferred averages over layers of constant optical depth are prone to the increasing interference of the reversed granulation towards higher effective temperature, in particular at low metallicity.

We present the Stagger-grid, a comprehensive grid of time-dependent, 3D hydrodynamic model atmospheres for late-type stars with realistic treatment of radiative transfer, covering a wide range in stellar parameters. This grid of 3D models is intended for various applications like stellar spectroscopy, asteroseismology and the study of stellar convection. In this introductory paper, we describe the methods used for the computation of the grid and discuss the general properties of the 3D models as well as their temporal and spatial averages (<3D>). All our models were generated with the Stagger-code, using realistic input physics for the equation of state (EOS) and for continuous and line opacities. Our ~220 grid models range in Teff from 4000 to 7000K in steps of 500K, in log g from 1.5 to 5.0 in steps of 0.5 dex, and [Fe/H] from -4.0 to +0.5 in steps of 0.5 and 1.0 dex. We find a tight scaling relation between the vertical velocity and the surface entropy jump, which itself correlates with the constant entropy value of the adiabatic convection zone. The range in intensity contrast is enhanced at lower metallicity. The granule size correlates closely with the pressure scale height sampled at the depth of maximum velocity. We compare the <3D> models with widely applied 1D models, as well as with theoretical 1D hydrostatic models generated with the same EOS and opacity tables as the 3D models, in order to isolate the effects of using self-consistent and hydrodynamic modeling of convection, rather than the classical mixing length theory approach. For the first time, we are able to quantify systematically over a broad range of stellar parameters the uncertainties of 1D models arising from the simplified treatment of physics, in particular convective energy transport. In agreement with previous findings, we find that the differences can be significant, especially for metal-poor stars.

We compare limb darkening laws derived from 3D hydrodynamical model atmospheres and 1D hydrostatic MARCS models for the host stars of the two transiting exoplanet systems HD 209458 and HD 189733. The surface brightness distribution of the stellar disks is calculated using 3D LTE spectrum formation and opacity sampling. We test our predictions using least-squares fits of model light curves to primary eclipses that were observed with the Hubble Space Telescope (HST). The limb darkening law derived from the 3D model of HD 209458 between 2900 A and 5700 A produces significantly better fits to the HST data, removing systematic residuals that were previously observed for model light curves based on 1D predictions. This difference arises mainly from the shallower mean temperature structure of the 3D model, which is a consequence of the explicit simulation of surface granulation. In the case of HD 189733, the model atmospheres produce practically equivalent limb darkening curves between 2900 A and 5700 A, partly due to obstruction by spectral lines, and the data are not sufficient to distinguish between the light curves. We also analyze HST observations between 5350 A and 10500 A for this star; the 3D model leads to a better fit compared to 1D predictions. The significant improvement of fit quality for the HD 209458 system demonstrates the higher degree of realism of 3D models and the importance of surface granulation for the formation of the atmospheric radiation field of late-type stars. This result agrees well with recent investigations of limb darkening in the solar continuum and other observational tests. The case of HD 189733 is no contradiction as the model light curves are less sensitive to the temperature stratification of the atmosphere and the observed data in the 2900 A - 5700 A region are not sufficient to distinguish more clearly between 3D and 1D limb darkening.

Radiative hydrodynamic simulations of solar and stellar surface convection have become an important tool for exploring the structure and gas dynamics in the envelopes and atmospheres of late-type stars and for improving our understanding of the formation of stellar spectra. We quantitatively compare results from three-dimensional, radiative hydrodynamic simulations of convection near the solar surface generated with three numerical codes CO5BOLD, MURaM, and STAGGER and different simulation setups in order to investigate the level of similarity and to cross-validate the simulations. For all three simulations, we considered the average stratifications of various quantities (temperature, pressure, flow velocity, etc.) on surfaces of constant geometrical or optical depth, as well as their temporal and spatial fluctuations. We also compared observables, such as the spatially resolved patterns of the emerging intensity and of the vertical velocity at the solar optical surface as well as the center-to-limb variation of the continuum intensity at various wavelengths. The depth profiles of the thermodynamical quantities and of the convective velocities as well as their spatial fluctuations agree quite well. Slight deviations can be understood in terms of differences in box size, spatial resolution and in the treatment of non-gray radiative transfer between the simulations. The results give confidence in the reliability of the results from comprehensive radiative hydrodynamic simulations.

Motivated by the work of Guillot (2010), we present a semi-analytical formalism for calculating the temperature-pressure profiles in hot Jovian atmospheres which includes the effects of clouds/hazes and collision-induced absorption. Using the dual-band approximation, we assume that stellar irradiation and thermal emission from the hot Jupiter occur at distinct wavelengths ("shortwave" versus "longwave"). For a purely absorbing cloud/haze, we demonstrate its dual effect of cooling and warming the upper and lower atmosphere, respectively, which modifies, in a non-trivial manner, the condition for whether a temperature inversion is present in the upper atmosphere. The warming effect becomes more pronounced as the cloud/haze deck resides at greater depths. If it sits below the shortwave photosphere, the warming effect becomes either more subdued or ceases altogether. If shortwave scattering is present, its dual effect is to warm and cool the upper and lower atmosphere, respectively, thus counteracting the effects of enhanced longwave absorption by the cloud/haze. We make a tentative comparison of a 4-parameter model to the temperature-pressure data points inferred from the observations of HD 189733b and estimate that its Bond albedo is approximately 10%. Besides their utility in developing physical intuition, our semi-analytical models are a guide for the parameter space exploration of hot Jovian atmospheres via three-dimensional simulations of atmospheric circulation.

We investigate the effects of coherent isotropic continuum scattering on the formation of spectral lines in local thermodynamic equilibrium (LTE) using 3D hydrodynamical and 1D hydrostatic model atmospheres of red giant stars. Continuum flux levels, spectral line profiles and curves of growth for different species are compared with calculations that treat scattering as absorption. Photons may escape from deeper, hotter layers through scattering, resulting in significantly higher continuum flux levels beneath a wavelength of 5000 A. The magnitude of the effect is determined by the importance of scattering opacity with respect to absorption opacity; we observe the largest changes in continuum flux at the shortest wavelengths and lowest metallicities; intergranular lanes of 3D models are more strongly affected than granules. Continuum scattering acts to increase the profile depth of LTE lines: continua gain more brightness than line cores due to their larger thermalization depth in hotter layers. We thus observe the strongest changes in line depth for high-excitation species and ionized species, which contribute significantly to photon thermalization through their absorption opacity near the continuum optical surface. Scattering desaturates the line profiles, leading to larger abundance corrections for stronger lines, which reach -0.5 dex at 3000 A for Fe II lines in 3D with excitation potential 2 eV at [Fe/H]=-3.0. The corrections are less severe for low-excitation lines, longer wavelengths, and higher metallicity. Velocity fields increase the effects of scattering by separating emission from granules and intergranular lanes in wavelength. 1D calculations exhibit similar scattering abundance corrections for weak lines, but those for strong lines are generally smaller compared to 3D models and depend on the choice of microturbulence.

We present Hubble Space Telescope optical and near-ultraviolet transmission spectra of the transiting hot-Jupiter HD189733b, taken with the repaired Space Telescope Imaging Spectrograph (STIS) instrument. The resulting spectra cover the range 2900-5700 Ang and reach per-exposure signal-to-noise levels greater than 11,000 within a 500 Ang bandwidth. We used time series spectra obtained during two transit events to determine the wavelength dependance of the planetary radius and measure the exoplanet's atmospheric transmission spectrum for the first time over this wavelength range. Our measurements, in conjunction with existing HST spectra, now provide a broadband transmission spectrum covering the full optical regime. The STIS data also shows unambiguous evidence of a large occulted stellar spot during one of our transit events, which we use to place constraints on the characteristics of the K dwarf's stellar spots, estimating spot temperatures around Teff~4250 K. With contemporaneous ground-based photometric monitoring of the stellar variability, we also measure the correlation between the stellar activity level and transit-measured planet-to-star radius contrast, which is in good agreement with predictions. We find a planetary transmission spectrum in good agreement with that of Rayleigh scattering from a high-altitude atmospheric haze as previously found from HST ACS camera. The high-altitude haze is now found to cover the entire optical regime and is well characterised by Rayleigh scattering. These findings suggest that haze may be a globally dominant atmospheric feature of the planet which would result in a high optical albedo at shorter optical wavelengths.

Context: Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very different physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the different layers interact strongly. These physical regimes are very diverse and it takes a highly efficient massively parallel numerical code to solve the associated equations. Aims: The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods: The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results: Bifrost completes the tests with good results and shows near linear efficiency scaling to thousands of computing cores.

Three-dimensional (3D) radiative hydrodynamic model atmospheres of metal-poor late-type stars are characterized by cooler upper photospheric layers than their 1D counterparts. This property of 3D models can dramatically affect elemental abundances derived from temperature-sensitive spectral lines. We investigate whether the cool surface temperatures predicted by metal-poor 3D models can be ascribed to the approximated treatment of scattering in the radiative transfer. We use the Bifrost code to test three different ways to handle scattering in 3D model atmospheres of metal-poor stars. First, we solve self-consistently the radiative transfer equation for a source function with a coherent scattering term. Second, we solve the radiative transfer equation for a Planckian source function, neglecting the contribution of continuum scattering to extinction in the optically thin layers; this has been the default mode in previous models of ours. Third, we treat scattering as pure absorption everywhere, which is the standard case in CO5BOLD models. We find that the second approach produces temperature structures with cool upper photospheric layers very similar to the correct coherent scattering solution. In contrast, treating scattering as pure absorption leads to significantly hotter and shallower temperature stratifications. The main differences in temperature structure between our published models and those generated with the CO5BOLD code can be traced to the different treatments of scattering. Neglecting the contribution of continuum scattering to extinction in optically thin layers provides a good approximation to the full radiative transfer solution for metal-poor stars. Our results demonstrate that the cool temperature stratifications predicted for metal-poor late-type stellar atmospheres by previous models of ours are not an artifact of the approximated treatment of scattering.

We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350\,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.

We report on a detailed abundance analysis of two strongly r-process enhanced, very metal-poor stars newly discovered in the HERES project, CS 29491-069 ([Fe/H]=-2.51, [r/Fe]=+1.1) and HE 1219-0312 ([Fe/H]=-2.96, [r/Fe]=+1.5). The analysis is based on high-quality VLT/UVES spectra and MARCS model atmospheres. We detect lines of 15 heavy elements in the spectrum of CS 29491-069, and 18 in HE 1219-0312; in both cases including the Th II 4019 {\AA} line. The heavy-element abundance patterns of these two stars are mostly well-matched to scaled solar residual abundances not formed by the s-process. We also compare the observed pattern with recent high-entropy wind (HEW) calculations, which assume core-collapse supernovae of massive stars as the astrophysical environment for the r-process, and find good agreement for most lanthanides. The abundance ratios of the lighter elements strontium, yttrium, and zirconium, which are presumably not formed by the main r-process, are reproduced well by the model. Radioactive dating for CS 29491-069 with the observed thorium and rare-earth element abundance pairs results in an average age of 9.5 Gyr, when based on solar r-process residuals, and 17.6 Gyr, when using HEW model predictions. Chronometry seems to fail in the case of HE 1219-0312, resulting in a negative age due to its high thorium abundance. HE 1219-0312 could therefore exhibit an overabundance of the heaviest elements, which is sometimes called an "actinide boost".

We present a preliminary local thermodynamic equilibrium (LTE) abundance analysis of the template halo red giant HD122563 based on a realistic, three-dimensional (3D), time-dependent, hydrodynamical model atmosphere of the very metal-poor star. We compare the results of the 3D analysis with the abundances derived by means of a standard LTE analysis based on a classical, 1D, hydrostatic model atmosphere of the star. Due to the different upper photospheric temperature stratifications predicted by 1D and 3D models, we find large, negative, 3D-1D LTE abundance differences for low-excitation OH and Fe I lines. We also find trends with lower excitation potential in the derived Fe LTE abundances from Fe I lines, in both the 1D and 3D analyses. Such trends may be attributed to the neglected departures from LTE in the spectral line formation calculations.