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.
Bisectors of the HARPS cross-correlation function (CCF) can discern between planetary radial-velocity (RV) signals and spurious RV signals from stellar magnetic activity variations. However, little is known about the effects of the stellar atmosphere
on CCF bisectors or how these effects vary with spectral type and luminosity class. Here we investigate the variations in the shapes of HARPS CCF bisectors across the HR diagram in order to relate these to the basic stellar parameters, surface gravity and temperature. We use archive spectra of 67 well studied stars observed with HARPS and extract mean CCF bisectors. We derive previously defined bisector measures (BIS, v_bot, c_b) and we define and derive a new measure called the CCF Bisector Span (CBS) from the minimum radius of curvature on direct fits to the CCF bisector. We show that the bisector measures correlate differently, and non-linearly with log g and T_eff. The resulting correlations allow for the estimation of log g and T_eff from the bisector measures. We compare our results with 3D stellar atmosphere models and show that we can reproduce the shape of the CCF bisector for the Sun.
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 de
rived 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.
One of the most important factors in determining the stellar lithium abundance is the effective temperature. In a previous study by the authors, new effective temperatures Teff for sixteen metal-poor halo dwarfs were derived using a local thermodynam
ic equilibrium (LTE) description of the formation of Fe lines. This new Teff scale reinforced the discrepancy. For six of the stars from our previous study we calculate revised temperatures using a non-local thermodynamic equilibrium (NLTE) approach. These are then used to derive a new mean primordial lithium abundance in an attempt to solve the lithium discrepancy. Using the code MULTI we calculate NLTE corrections to the LTE abundances for the Fe I lines measured in the six stars, and determine new Teffs. We keep other physical parameters, i.e. log g, [Fe/H] and xi, constant at the values calculated in Paper I. With the revised Teff scale we derive new Li abundances. We compare the NLTE values of Teff with the photometric temperatures of Ryan et al. (1999, ApJ, 523, 654), the infrared flux method (IRFM) temperatures of Melendez & Ramirez (2004, ApJ, 615, 33), and the Balmer line wing temperatures of Asplund et al. (2006, ApJ, 644, 229). We find that our temperatures are hotter than both the Ryan et al. and Asplund et al. temperatures by typically ~ 110 K - 160 K, but are still cooler than the temperatures of Melendez & Ramirez by typically ~ 190 K. The temperatures imply a primordial Li abundance of 2.19 dex or 2.21 dex, depending on the magnitude of collisions with hydrogen in the calculations, still well below the value of 2.72 dex inferred from WMAP + BBN. We discuss the effects of collisions on trends of 7Li abundances with [Fe/H] and Teff, as well as the NLTE effects on the determination of log g through ionization equilibrium, which imply a collisional scaling factor SH > 1 for collisions between Fe and H atoms.
