Free Access
Issue
A&A
Volume 611, March 2018
Article Number A11
Number of page(s) 16
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/201732147
Published online 14 March 2018

© ESO 2018

1 Introduction

The stellar atmosphere is the boundary to the opaque stellar interior, and serves as the link between observations and the models of stellar structure and evolution. The phenomena of stellar evolution manifest themselves in the stellar surface as changes in chemical composition and in fundamental stellar parameters such as radius, surface gravity, effective temperature, and luminosity. The information we use to study distant stars comes from the flux they have emitted. However, the atmospheric layers where this flux forms is the transition region between convective and radiative regime. Thus, the surface structures and dynamics of cool stars are characterised by the presence of convective motions and turbulent flows. Convection manifests in the surface layers as a particular pattern of downflowing cooler plasma and bright areas where hot plasma rises (Nordlund et al. 2009). The size of granules depends on the stellar parameters of the star and, as a consequence, on the extension of their atmosphere (e.g. Magic et al. 2013). Eventually, the convection causes an inhomogeneous stellar surface that changes with time. They affect the atmospheric stratification in the region where the flux forms and also affect the emergent spectral energy distribution (SED), with potential effects on the precise determinations of stellar parameters (e.g. Bigot et al. 2011; Creevey et al. 2012; Chiavassa et al. 2012), radial velocity (e.g. Bigot & Thévenin 2008; Chiavassa et al. 2011; Allende Prieto et al. 2013), chemical abundance (e.g. Asplund et al. 2005, 2009; Caffau et al. 2011), photometric colours (Bonifacio et al. 2017), and on planet detection (Magic et al. 2015; Chiavassa et al. 2017).

Convection is a difficult process to understand because it is non-local, and 3D, and it involves non-linear interactions over many disparate length scales. In this context, the use of numerical 3D radiative hydrodynamical simulations of stellar convection is extremely important. In recent years, with increased computational power, it has been possible to compute grids of 3D simulations that cover a substantial portion of the Hertzsprung–Russell diagram (Magic et al. 2013; Trampedach et al. 2013; Ludwig et al. 2009). With these tools it is possible to predict reliable synthetic spectra for several stellar types.

Photometric systems and filters are designed to be sensitive to temperature, gravity, and metal abundance indicators and thereby to complement spectroscopic determinations of the fundamental properties of stars. In addition, the integrated magnitudes and colours of stars can be used to infer the ages, metallicities, and other properties of the underlying stellar populations (e.g. Casagrande & VandenBerg 2014). For these purposes, several broad-band, or intermediate- and/or narrow-band filters have been designed to probe different regions of stellar spectra sensitive to different atmospheric parameters (Bessell & Murphy 2012; Gunn et al. 2006; Cohen et al. 2003). Additionally, there are the photometric systems used by the Gaia mission.

Gaia (Gaia Collaboration 2016) is an astrometric, photometric, and spectroscopic spaceborne mission of a large part of the Milky Way. Apart from the astrometric instrument, Gaia carries on board two low-resolution spectrophotometers (Blue and Red Prism, BP/RP, Bailer-Jones et al. 2013) and the Radial Velocity Spectrometer (RVS; Katz et al. 2004). The photometric instrument measures the SED over 120 pixels of all detected objects at the same angular resolution and at the same epoch as the astrometric observations. The BP operates in the range 3300–6800 Å, while the RP uses therange 6400–10 500 Å. The main aims of this instrument are to provide proper classifications (e.g. distinguish between stars, galaxies, and quasars) and characterisations (e.g. reddenings and stellar parameters), and to enable chromatic corrections of the astrometric centroid data. Finally, the integral-field spectrograph RVS provides spectra between 8470 and 8710 Å at a spectral resolving power of ≈11 200. The RVS is expected to produce radial velocities through Doppler-shift measurements; interstellar reddening, atmospheric parameters, and projected rotational velocities; and individual element abundances for some elements. In this work, we calculated synthetic stellar spectra and photometry for broad-band and Gaia systems obtained using realistic 3D radiative hydrodynamical simulations of stellar convection from the STAGGER-grid. The predicted photometric 3D colours and the 3D spectra are publicly available for the community through the POLLUX database. The spectra corresponding to the Gaia RVS spectral range will be used for the calibration of the instrument to preserve the measured radial velocities from the convection shift (see forthcoming paper of GaiaDataRelease2/RVS).

2 Methods

2.1 Stellar model atmospheres

Magic et al. (2013) described the STAGGER-grid of realistic 3D radiative hydrodynamical simulations of stellar convection for cool stars using the STAGGER-code (originally developed by Nordlund & Galsgaard 19952, and continuously improved over the years by its user community), a state-of-the-art (magneto) hydrodynamic code that solves the time-dependent hydrodynamic equations for mass-, momentum-, and energy-conservation, coupled with the 3D radiative transfer equation in order to account correctly for the interaction between the radiation field and the plasma. The code uses periodic boundary conditions horizontally and open boundaries vertically. At the bottom of the simulation, the inflows have a constant entropy. The outflows are not tightly constrained and are free to pass through the boundary. The code is based on a sixth-order explicit finite-difference scheme and a fifth-order interpolation. The considered large number over wavelength points is merged into 12 opacity bins (Nordlund 1982; Skartlien 2000; Magic et al. 2013). STAGGER simulations are based on a realistic equation of state that accounts for ionisation, recombination, and dissociation (Mihalas et al. 1988); continuous absorption and scattering coefficients listed in Hayek et al. (2010); and the line opacities listed in Gustafsson et al. (2008). These are in turn based on the VALD-2 database (Stempels et al. 2001) of atomic lines and the SCAN-base (Jørgensen 1997) of molecular lines.

2.2 Three-dimensional radiative transfer

We used the 3D pure-LTE radiative transfer code OPTIM3D (Chiavassa et al. 2009) to compute synthetic spectrum from the snapshots of the radiative-hydrodynamical (RHD) simulations of the STAGGER-grid (see Fig. 1 in Magic et al. 2013). The code takes into account the Doppler shifts due to convective motions. The radiative transfer equation is solved monochromatically using pre-tabulated extinction coefficients as a function of temperature, density, and wavelength. The lookup tables were computed for the same chemical compositions as the RHD simulations using the same extensive atomic and molecular continuum and line opacity data as the latest generation of MARCS models (Gustafsson et al. 2008) with the addition – with respect to Table 2 of Gustafsson’s paper – of the SiS molecule (Cami et al. 2009), which is particularly important for the far-infrared region of the spectrum. While the sources of line opacities used in the RHD simulations of Magic et al. (2013) and in OPTIM3D are the same, the data for the continuum opacities are almost the same: Hayek et al. (2010) reported that the data used in RHD simulations are mostly identical to those used in the MARCS models, but include additional bound-free data from the Opacity Project and the Iron Project (Trampedach et al., priv. comm.) as well as some opacities of the second ionisation stage for many metals.

For the computation of the spectra from RHD simulations, the assumed microturbulence is equal to zero since the velocity fields inherent in RHD models are expected to self-consistently and adequately account for non-thermal Doppler broadening of spectral lines (Asplund 2000). The temperature and density ranges spanned by the tables are optimised for the values encountered in the RHD simulations. The detailed methods used in the code are explained in (2009, 2010). OPTIM3D has already been employed in synergy with the STAGGER simulations in several works (Chiavassa et al. 2010, 2011, 2012, 2014a, 2015, 2017; Magic et al. 2015) either concerning the extraction of synthetic spectra or interferometric observables.

2.3 One-dimensional radiative transfer

For all the following comparisons with 3D simulations, we used plane-parallel, hydrostatic, 1D atmosphere models computed with a similar physical treatment to the MARCS code and the same equation of state and opacities as in the individual 3D simulations (ATMO; Magic et al. 2013). Moreover, we used a 1D version of OPTIM3D, and the chemical compositions, the opacities, and numerics of the radiative transfer calculations for the emergent intensities are the same as used in the 1D and 3D approaches (Fig. 2).

3 Synthetic spectra from 0.1 to 20 μm

The STAGGER-grid includes 3D stellar atmosphere simulations with metallicities [Fe/H] = +0.5, 0.0, −0.5, −1.0, −2.0, −3.0, and −4.0; surface gravity log g between 1.5 and 5.0 in steps of 0.5 dex; and effective temperature Teff from 4000 to 7000 K in steps of 500 K (Fig. 1 of Paper I). In this work we present the synthetic spectra computed for the STAGGER-grid for a total of 181 simulations (Table B.1). The spectra have been calculated with a constant resolving power of λΔλ = 20 000 (nλ = 105 767 wavelength points) from 1000 to 200 000 Å. OPTIM3D computes the emerging intensities for vertical rays cast through the computational box for all required wavelengths. The procedure is repeated after tilting the computational box by an angle θ with respect to the line of sight (vertical axis) and rotating it azimuthally by an angle ϕ. The final result is a spatially resolved intensity spectrum at different angles. In addition, a temporal average is also performed. We performed the calculations for ten snapshots of the 3D RHD simulations of Table B.1, adequately spaced so as to capture several convective turnovers, for ten different inclination angles μ=cos(θ) $\mu=\cos(\theta)$ = [1.00, 0.90, 0.80, 0.70, 0.50, 0.30, 0.20, 0.10, 0.05, 0.01] (see Fig. 4), and four ϕ-angles [0 , 90 , 180 , 270 ]. The strongest decline in the limb darkening is usually found towards the limb; therefore, we decided to resolve with more μ-angles at the limb instead of having an equidistant scale in μ. We tested the discrepancy between the temporal average using a large number of snapshots (e.g. 20) and using only 10 snapshots is lower than 0.3%. The number of ten snapshots was chosen because it represents the best compromise in terms of computational time and accuracy among the whole set of stellar parameters. All things considered, we computed 400 spectra in the range 1000–200 000 Å for every simulation.

Figure 1 displays the set of all synthetic spectra computed. We determined the Teff from the integration of the SED of the spectra from 0.1 to 20 μm. The effective temperature has been computed using Stefan–Boltzmann law as T effspectra = { [ λ 1 λ 2 f ( λ )dλ ]/σ } 0.25 , $T_{\rm{eff}\ \rm{spectra}} = \left\{\left[\int_{\lambda_1}^{\lambda_2}f\left(\lambda\right)\mathrm{d}\lambda\right]/\sigma\right\}^{0.25} ,$(1) where λ1 = 1010Å and λ2 = 199 960 Å, f( λ ) $f\left(\lambda\right)$ is the synthetic flux, and σ is the Stefan–Boltzmann constant. The values of the effective temperatureare listed in Table B.1. Figure 1 shows that increasing Teff returns higher radiated energy per surface area and the peak of the radiation curve movesto shorter wavelengths, as expected by Planck law. Pereira et al. (2013) provided excellent agreement of their 3D solar simulation of the STAGGER-grid with the continuum observation of the Sun. As they did, we used the Kurucz (2005) irradiance3 and normalised flux atlases for the Sun between 3000 and 10 000 Å and found a good agreement (Fig. 3), reinforcing the view that the simulation of thermodynamic structure and post-processing detailed radiative transfer are realistic. This conclusion was also reported by Hayek et al. (2012), who determined that the numerical resolution of the STAGGER 3D RHD models and the spectral resolution for the flux computations are sufficient to predict realistic observables. In particular, some of the RHD simulations presented in this work and for a limited spectral region between 2000 and 10 000 Å have been used in Magic et al. (2015) to provide appropriate coefficients for various bi-parametric and non-linear limb darkening laws.

thumbnail Fig. 1

Surface rendering for all the synthetic spectra computed for the 3D RHD simulations in Table B.1. The vertical bar on the right displays the colour scale for the emerging flux in erg s−1 cm−2−1. For clarity, the wavelength range has been reduced to 1000–25 000 Å.

thumbnail Fig. 2

Synthetic spectra of the solar simulation in the spectral range 2000–25 000 Å and for the different μ=cos(θ) $\mu=\cos(\theta)$ inclination angles used in the computation, where θ is the angle with respect to the line of sight (vertical axis).

thumbnail Fig. 3

Comparison of the solar simulation (black) with the observed flux of the Sun (red, Kurucz 2005). The solarirradiance is converted to flux at the solar surface using the multiplicative factor of [(1AU )R]2=46202$\left[\left(1 \rm{AU}\right)/R_{\odot}\right]^2 = 46\,202$. For clarity, the spectra have been resampled to a lower spectral resolution with fewer frequency points (nλ = 2115).

thumbnail Fig. 4

Synthetic spectrum of the solar simulation at full spectral resolution in the spectral range 2000–25 000 Å (grey, top panel) and 2000–12 000 Å (central and bottom panels). Several system response functions (Table 1), from which synthetic colours have been computed, are overplotted. Johnson-Cousins system response functions (U, B, V, Rc, Ic) are plotted in green and 2MASS in pink-violet (top left panel); SDSS (u, g, r, i, z) in yellow-red and Gaia (BP, RP, G) (top right panel); Strömgren (uvby) in red-blue and SkyMapper (us , vs , gs, rs, is, zs) (bottom left panel); and the 15 filters of the HST-WFC3 (bottom right panel). For clarity, SDSS and SkyMapper functions are normalised to 0.5.

4 Photometric synthetic observables

Photometric systems and filters are designed to probe fundamental physical parameters, such as the effective temperature, surface gravity, and metallicity of stars. Colour and magnitude relations are used for a variety of purposes: interpreting the observed distribution of stars in colour–colour and colour–magnitude diagrams, deriving distances to stars and star clusters, and testing stellar evolutionary theory by comparing with observations to name just a few. Thus, it is important to have realistic model fluxes to generate colours which match the observed values. In addition to synthetic model fluxes, details on the photometric standardisation are also a part of this quest.

In essence, photometry condenses the information encoded in a spectrum f( λ ) $f\left(\lambda\right)$ over a system response function T( λ ) $T\left(\lambda\right)$, i.e. f ( λ )T( λ )dλ $\int f\left(\lambda\right) T\left(\lambda\right) \mathrm{d}\lambda$. Each existing photometric system then varies in the details. Most notably, T( λ ) $T\left(\lambda\right)$ will depend on the filter under consideration and the response function of the detector. This means that a distinction must be made between photo-counting and energy-integration detectors, meaning that a measurement of energy f ( λ )T( λ )dλ $\int f\left(\lambda\right) T\left(\lambda\right) \mathrm{d}\lambda$ will correspond to (hc) 1 f ( λ )λT( λ )dλ $(hc)^{-1} \int f\left(\lambda\right) \lambda T\left(\lambda\right) \mathrm{d}\lambda$ photons (see e.g. Bessell 2000). Another aspect that often varies among different photometric systems is how their standardisation (zero-point and absolute calibration) is achieved. Here, for all systems but Gaia we adopted the exact same procedure as used by Casagrande & VandenBerg (2014), where details on the adopted filter transmission curves, the photo-counting and energy-integration formalism, and zero-points and absolute calibration can be found4. We computed synthetic colours in the Johnson-Cousins, SDSS, 2MASS, Gaia, SkyMapper, Strömgren, HST-WFC3, and Gaia systems (Table 1 and Fig. 4 for a comparison of the solar spectrum with the filter transmission curves studied here). For the HST-WFC3 systems our tables are provided in the VEGA, ST, and AB systems.

A full characterisation of the Gaia photometric system, including zero-points and standardisation is expected to be released in 2018. In this work, we used the transmission curves available from the ESA website5, and computed Gaia colours following Jordi et al. (2010). We fixed Vega’s magnitudes to be G = BP = RP = 0.03, and for the absolute calibration used a Kurucz synthetic Vega spectrum rescaled to the measured flux value at 5556 Å from Megessier (1995).

Similarly to Casagrande & VandenBerg (2014), instead of colour indices we provide bolometric corrections in different bands (Table B.1 and 4–7, available at the CDS) because they are more versatile and can be rearranged in any colour combination, as follows from Equations. (3) and (4). The bolometric magnitude is defined as M Bol =2.5log L L + M Bol, , $ M_{\rm{Bol}} = -2.5 \log\frac{L}{L_{\odot}} + M_{\rm{Bol},\odot}, $(2) where we adopt MBol,⊙ = 4.74. It follows that the bolometric correction in a given band BCζ is B C ζ = m Bol m ζ = M Bol M ζ , $BC_{\zeta} = m_{\mathrm{Bol}}-m_{\zeta}=M_{\mathrm{Bol}}-M_{\zeta}, $(3)

where the lower and upper cases refer to apparent and absolute magnitudes, respectively. From this it follows that colour indices can be obtained from the difference in bolometric corrections, where ζ and η are two given bands: ζη= m ζ m η =B C η B C ζ . $\zeta - \eta = m_{\zeta} - m_{\eta} = BC_{\eta} - BC_{\zeta}. $(4)

Thus, in the rest of the paper when we talk about synthetic colours, these have been obtained as differences in bolometric corrections from our tables.

Table 1

Photometric systems used in this work and overplotted on the synthetic spectra in Fig. 4.

4.1 Microturbulence

The stellar surface convection produces a velocity field where the emerging spectral lines form. The Doppler broadening of these lines is a direct consequence of the velocity field in these crucial layers (Asplund et al. 2000b; Nordlund et al. 2009). In traditional 1D models, this effect can be accounted for by the use of arbitrary micro- and macroturbulence parameters. Full 3D line formation calculations using 3D RHD simulations have demonstrated that in late-type stars the required non-thermal Doppler line broadening is fully included in the convection-related motions of the stellar atmosphere (e.g. Collet et al. 2007). One-dimensional microturbulence represents the small-scale end of turbulent motions and is applied to the spectral line absorption coefficient. It affects the strong lines to a greater extent, reducing their saturation, and to a lesser extent the widths of weak lines. For 1D-based SEDs, microturbulence partly redistributes the flux in spectral regions probed by the photometricsystems, in particular in regions crowded with lines towards the blue and the ultraviolet, and in filters with smaller wavelength coverages (Casagrande & VandenBerg 2014).

The values of the microturbulence parameters are usually determined by comparing synthetic and observed spectral line profiles and line strengths and often using a depth-independent value. For reference, a typical value for dwarfs and subgiants is around 1−1.5 km s−1, which increases to 2−2.5 km s−1 for stars on the red giant branch (e.g. Gray et al. 2001). A constant value of 2 km s−1 is usually assumed in large grids of synthetic stellar spectra (Castelli & Kurucz 2004; Brott & Hauschildt 2005). To compute our 1D hydrostatic comparison models, we explored different values of microturbulence: 0, 1, and 2 km s−1 . We found that there is no clear and no unique relation between microturbulence and the stellar parameters, as reported by Casagrande & VandenBerg (2014). For clarity we adopted, as a guiding example, a value of 1 km s−1 when performing 1D calculations.

4.2 Three-dimensional versus one-dimensional bolometric differences in correction

The figures in the Appendix display the bolometric corrections between 3D simulations and the corresponding 1D hydrostatic models with microturbulence = 1 km s−1. The values of BC are reported in Table B.1 to Table 8 for all the filters and, to retrieve the absolute colours, Equations. (3) and (4) should be used. Considering the SDSS, SkyMapper, 2MASS, and HST-WFC3 systems, the overall deviations are limited to small fraction (less than 5%) from BCr toBCKs, but they increase to 10% for BCu and in BCg (SDSS) and for BCgsky, BCvsky, and BCusky (SkyMapper) where the optical and line crowded region of the spectrum is probed with rather narrow filters (Fig. 4, bottom); these differences decline with increasing effective temperature. This behaviour is more or less visible for all the photometric systems; there is no clear correspondence of ΔBC with the other stellar parameters. On a broad scale, the wide infrared photometric systems like 2MASS (BCH and BCKs) and optical Gaia (BCG) display a noticeable offset with respect to ΔBC = 0.0. This is due to a redistribution of the spectral energy flux among the different filters and is a direct effect of the impact of 3D dynamics and thermodynamic structure on spectral line formation.

Gaia photometric systems return 3D and 1D deviations of less than 3% with higher values for the bluer system (BP). This effect may be not negligible and should become important for future releases of Gaia data. For this purpose, part of the spectra presented in this work have already been provided to Gaia consortium-CU86.

5 Convective velocity shifts for RVS

Measurements of stellar radial velocities are fundamental in order to determine stellar space velocities. This is needed, for example to investigate the kinematic structure of stellar populations in the Galaxy or to monitor for radial velocity variations, either of which would point to the presence of unseen companion(s). Convection plays a crucial role in the formation of spectral lines and deeply influences the shape, shift, and asymmetries of lines in late-type stars (e.g. Asplund et al. 2000a). These stars represent most of the objects that will be observed during the Gaia mission. Absorption lines may be blueshifted as a result of convective movements in the stellar atmosphere: bright and rising convective elements contribute more photons than the cool dark shrinking gas, and as a consequence, the absorption lines appear blueshifted (Dravins 1982). However, the convective line shift is not the same for all the spectral lines. Each line has a unique fingerprint in the spectrum that depends on line strength, depth, shift, width, and asymmetry across the granulation pattern depending on their height of formation and sensitivity to the atmospheric conditions. In this context, the line strengths play a major role (Asplund et al. 2000c).

The aim of the present section is to derive the overall convective shift for 3D simulations. First, we computed the 1D and 3D spectra with a constant resolving power of λ∕Δ λ = 300 000 from 8470 to 8710 Å for a limited number of 3D simulations (see Table B.1) covering stellar parameters observed by RVS (i.e. [Fe/H] ≥ −2.0). Then, from our spectra we selected only a series of non-blended Fe I lines and masked the others (Fig. 5). The oscillator strengths of these Fe I lines (Table 2) have been accurately determined by Bigot & Thévenin (2006) using 3D RHD simulation where the Fe I and Ca II lines are indicated. It should be noted that the synthetic spectra, when compared to the observations, have to be gravitationally redshifted (e.g. Pasquini et al. 2011)by a certain amount corresponding to the type of star considered (e.g. for the Sun it is 636.486 ± 0.024 m s−1 Lindegren & Dravins 2003). Gravitational shifts for late-type dwarfs (log g ≈ 4.5) range between 0.7 and 0.8 km s−1 and they dramatically decrease with surface gravity down to 0.02–0.03 kms−1 for K giant stars with log g ≈ 1.5 (Allende Prieto et al. 2013).

The velocity gradient through the photosphere sets the basic shape of the absorption lines in terms of asymmetry and position of the emerging intensity. One way to detect the asymmetries in the line is the bisector7 . A symmetric profile has a straight vertical bisector (i.e. in the case of hydrostatic 1D spectra). The spectral lines with C-shape bisectors are formed mostly in the upflows (granules) and therefore blueshifted. The reverse C-shapes are generally formed in downflows (Dravins et al. 1981). Reversed C-shape bisectors can be explained by a combination of a steep decline in velocities with height with a flux deficit spanning only a fraction of the red wing of the line profiles (Gray 2010). Different articles show the presence of bisectors revealing asymmetriesand wavelength shifts that indicate the presence of granulation for several kinds of stars (e.g. Ramírez et al. 2008; Gray 2009).

Figure 6 shows the line bisectors for the Fe I and Ca II triplet lines for stars with different Teff and log g, but with the same metallicity. The gas is strongly horizontally divergent due to mass conservation and its velocities diminish with height. Weak lines (with typically high excitation potential), which form in deeper layers, are more blueshifted than strong lines whose core and part of the wings are formed in higher layers. This effect is particularly visible in Figure 6 when comparing the solar bisector with the hottest Teff = 7000 K simulations. In addition, the velocity field in 3D simulations of STAGGER-grid largely affects the overall shape of the iron lines in the range of RVS and for all the stars with the strongest effects for Ca II. This has already been shown for other spectral regions (e.g. Asplund et al. 2000b; Allende Prieto et al. 2002; Ramírez et al. 2009; Pereira et al. 2013; Magic et al. 2014).

We determined the convective shift considering only Fe I and only Ca II triplet lines, we cross-correlated each 3D spectrum with the corresponding 1D by using a lag vector corresponding to radial velocities (RV) in the range −10 < v < +10 km s−1 for Fe I, and −150 < v < +150 km s−1 for Ca II triplet lines, in steps of 0.3 km s−1. These velocity ranges were chosen to largely cover the wavelength frequency points of all the single lines. For each RV value, we Doppler-shifted the 1D spectrum and computed its cross-correlation function (CCF) with the 3D spectrum. The final step is to compute the weighted average to obtain the location of the CCF maximum, which corresponds to the actual 3D convective shifts (CS) with respect to 1D models: CS= 10 +10 R V( v )CCF( v )dv 10 +10 C CF( v )dv $ CS = \frac{\int^{&#x002B;10}_{-10}RV\left(v\right)\cdot CCF\left(v\right)\textrm{d}v}{\int^{&#x002B;10}_{-10}CCF\left(v\right)\textrm{d}v} $(5)

Figure 7 displays the convective shifts for all the simulations either in the HR diagram (top left panel) or as a function of the Teff for the Fe I of Table 2 (top right panel), and for the Ca II triplet lines (bottom left panel). We found that surface gravity and metallicity have a small effect on the convective shifts, as already noticed by Allende Prieto et al. (2013). The values for the Fe I are in the range between −0.235 and +0.361 km s−1. The convective shifts of Ca II lines are strongly redshifted (as shown by red bisectors in Fig. 6) and are between −0.023 and +0.698 km s−1. In Fig. 7 (top right panel), there is a net correlation of the convective shifts with the effective temperature: Teff ≾ 4500 K denotes redshifts, while Teff ≿ 5000 K denotes blueshifts (except for the hottest Teff ≈ 7000 K). This result is in agreement with Allende Prieto et al. (2013), who had performed the calculations for a different set of iron lines and found a milder correlation, where Teff with warmer stars tend to exhibit larger blueshifts.

To quantify the differences in the convective shifts, we selected 26 simulations from Table B.1 with the same surface gravity, metallicity, and Δ Teff  < 50 K with respect to a subset of CIFIST-grid simulations from Allende Prieto et al. (2013). Convective shifts, as a function of metallicity (Table 3), from RHD simulations in this work are CSStagger,[Fe∕H]=0 = [−0.135, 0.142], CSStagger,[Fe∕H]=−1 = [−0.175, 0.277], and CSStagger,[Fe∕H]=−2 = [−0.114, 0.119] km s−1, and from CIFIST-grid simulations CSCIFIST,[Fe∕H]=0 = [−0.050, 0.054], CSCIFIST,[Fe∕H]=−1 = [−0.170, 0.198], and CSCIFIST,[Fe∕H]=−2 = [−0.316, 0.119] km s−1. The spanned shift values from both grids are similar (Fig. 7, bottom right panel), and show smaller deviations at solar metallicity (ΔCS ≤ 0.195 km s−1) and slightly larger deviations at [Fe/H] = −1 (Δ CS ≤ 0.370 km s−1) and [Fe/H] = −2 (Δ CS ≤ 0.221 km s−1). Apart from the possible numerical differences in the simulations and in the radiative transfer, the shift deviations may also be due to the set of spectral lines considered.

The extraction of accurate radial velocities from RVS needs an appropriate wavelength calibration from convective shifts. This is directly processed in RVS pipeline using the synthetic spectra presented in this work and provided to Gaia consortium-CU68.

thumbnail Fig. 5

Flux normalised to the continuum for the solar simulation (Table B.1) in the RVS range (8400–8900 Å). The Fe I lines (black arrows) are from Bigot & Thévenin (2006), while the Ca II triplet are indicated with red arrows.

Table 2

Central wavelength position, oscillator strength (log gf), and excitation potential (χ) for the 20 Fe I in the spectral domain of RVS (Bigot & Thévenin 2006).

thumbnail Fig. 6

Line bisectors of the 20 Fe I lines (black) from Bigot & Thévenin (2006) and Ca II triplet (red) for five 3D simulations in the grid.

thumbnail Fig. 7

Top left panel: convective shifts predicted by the 3D hydrodynamical simulations for the Gaia RVS spectral range and the Fe I lines (see text for details). Top right panel: convective shifts from Fe I lines as a function of the effective temperature of the 3D simulations in Table B.1. Bottom left panel: convective shifts from Ca II lines. Bottom right panel: comparison of convective shifts for a selected number of RHD simulations from this work (stars) with simulations with equivalent stellar parameters from Allende Prieto et al. (2013) (circles). See Table 3 for details.

Table 3

Selected convective shifts for Fe I only (CSStagger in km s−1 ) from RHD simulations in this work andfrom CIFIST-grid simulations (CSCIFIST in km s−1 ) with equivalent stellar parameters from Allende Prieto et al. (2013).

6 Conclusions

We provided synthetic spectra from the STAGGER-grid:

  • low-resolution spectra from 1000 to 200 000 Å with a constant resolving power of λ∕Δλ = 20 000;

  • high-resolution spectra from 8400 to 8900 Å (Gaia RVS spectral range), with a constant resolving power of λ∕Δλ = 300 000.

We used the low-resolution spectra to compute synthetic colours in the Johnson-Cousins UBV (RI)C, SDSS, 2MASS, Gaia systems, SkyMapper, Strömgren, and HST-WFC3. We extracted the bolometric corrections for the 3D simulations and the corresponding 1D hydrostatic models. We probed that 1D versus 3D deviations are limited to small fraction (less than 5%) except for BCu and to a lesser extent BCg, where the differences are larger than 10%. Systems u and g span the optical and line crowded region of the spectrum. Moreover, we showed that there is a clear correlation between effective temperature and 3D and 1D deviations (Δ BC): it decreases with increasing effective temperature.

The Gaia photometric system return 3D and 1D deviations of less than 3% with higher values for the bluer filter (BP). This effect should become important for future releases of Gaia data. For this purpose, part of the spectra presented in this work have already been provided to Gaia consortium-CU8.

We used the high-resolution spectra to denote the impact of the velocity gradient through the photosphere on the basic shape of the absorption lines. For this purpose, we reported the line bisectors of non-blended Fe I and Ca II triplet lines for different stars. We showed that weak lines (high excitation potential), which form in deeper layers, are more blueshifted than strong lines (low excitation potential), whose core and part of the wings are formed in higher layers.

As a final step to derive the overall convective shift for 3D simulations with respect to the reference 1D hydrostatic models, we cross-correlated each 3D spectrum with the corresponding 1D spectrum. The spanned values are between −0.235 and +0.361 km s−1. We showed a net correlation of the convective shifts with the effective temperature: lower Teff denotes redshifts and higher Teff blueshifts; this result is in agreement with Allende Prieto et al. (2013). In addition, we quantified the differences in the convective shifts between a subset of the RHD simulations in this work and the corresponding CIFIST-grid simulations. The spanned shift values from the two grids are similar, and show smaller deviations at solar metallicity. The extraction of accurate radial velocities from RVS spectra need an appropriate wavelength calibration from convection shifts. The spectra presented in this work have been provided to Gaia consortium-CU6 to directly process the observed spectra in RVS pipeline.

We have made all the spectra publicly available for the community through the POLLUX database (Palacios et al. 2010). POLLUX9 is a stellar spectra database proposing access to theoretical data including high-resolution synthetic spectra and spectral energy distributions from several model atmospheres. Continuous development either of the STAGGER-grid simulations or of the spectral synthesis calculations will be uploaded there in the future.

Acknowledgements

L.C. gratefully acknowledges support from the Australian Research Council (grants DP150100250, FT160100402). This work was granted access to the HPC resources of Observatoire de la Côte d’Azur – Mésocentre SIGAMM.

Appendix A One-dimensional versus three-dimensional bolometric corrections

thumbnail Fig. A.1

Bolometric correction (BC) differences computed for the photometric filters Johnson-Cousins and 2MASS with 3D simulations (Table B.1) and the corresponding 1D hydrostatic models with microturbulence = 1 km s−1: Δ BC = BC1DBC3D. The red dashed line indicates ΔBC = 0.0.

thumbnail Fig. A.2

Bolometric correction differences for the SkyMapper and Gaia photometric filters (Table 1). The notation is the same as in Fig. A.1.

thumbnail Fig. A.3

Bolometric correction differences for the HST-WFC3 provided in the VEGA system (for the ST and AB systems the differences are similar); first set of filters (Table 1). The notation is the same as in Fig. A.1.

thumbnail Fig. A.4

Bolometric correction differences for the HST-WFC3 provided in the VEGA system (for ST and AB systems the differences are similar); second set of filters (Table 1). The notation is the same as in Fig. A.1.

Appendix

Table B.1

RHD simulations’ steAppendix Bllar parameters (first three cols.), bolometric Magnitude (Mbol ), and bolometric correction (BC) for Johnson-Cousins, 2MASS, SDSS, and Gaia systems (Table 1).

References

  1. Allende Prieto, C., Asplund, M., García López, R. J., & Lambert, D. L. 2002, ApJ, 567, 544 [NASA ADS] [CrossRef] [Google Scholar]
  2. Allende Prieto, C., Koesterke, L., Ludwig, H.-G., Freytag, B., & Caffau, E. 2013, A&A, 550, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Asplund, M. 2000, A&A, 359, 755 [NASA ADS] [Google Scholar]
  4. Asplund, M., Ludwig, H., Nordlund, Å., & Stein, R. F. 2000a, A&A, 359, 669 [Google Scholar]
  5. Asplund, M., Nordlund, Å., Trampedach, R., Allende Prieto, C., & Stein, R. F. 2000b, A&A, 359, 729 [Google Scholar]
  6. Asplund, M., Nordlund, Å., Trampedach, R., & Stein, R. F. 2000c, A&A, 359, 743 [Google Scholar]
  7. Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Cosmic Abundances as Records of Stellar Evolution and Nucleosynthesis, eds. T. G. Barnes, III, & F. N. Bash, ASP Conf. Ser., 336, 25 [Google Scholar]
  8. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bailer-Jones, C. A. L., Andrae, R., Arcay, B., et al. 2013, A&A, 559, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bessell, M. S. 2000, PASP, 112, 961 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bessell M. S. 2011, PASP, 123, 1442 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bessell, M., & Murphy, S. 2012, PASP, 124, 140 [Google Scholar]
  13. Bessell, M., Bloxham, G., Schmidt, B., et al. 2011, PASP, 123, 789 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bigot, L., & Thévenin, F. 2006, MNRAS, 372, 609 [NASA ADS] [CrossRef] [Google Scholar]
  15. Bigot, L., & Thévenin, F. 2008, in SF2A-2008, 3, eds. C. Charbonnel, F. Combes, & R. Samadi [Google Scholar]
  16. Bigot, L., Mourard, D., Berio, P., et al. 2011, A&A, 534, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Bonifacio, P., Caffau, E., Ludwig, H.-G., et al. 2017, Mem. Soc. Astron. It., 88, 90 [NASA ADS] [Google Scholar]
  18. Brott, I., & Hauschildt, P. H. 2005, in The Three-Dimensional Universe with Gaia, ed. C. Turon, K. S. O’Flaherty, & M. A. C. Perryman, ESA SP, 576, 565 [Google Scholar]
  19. Caffau, E., Ludwig, H.-G., Steffen, M., Freytag, B., & Bonifacio, P. 2011, Sol. Phys., 268, 255 [NASA ADS] [CrossRef] [Google Scholar]
  20. Cami, J., Sloan, G. C., Markwick-Kemper, A. J., et al. 2009, ApJ, 690, L122 [NASA ADS] [CrossRef] [Google Scholar]
  21. Casagrande, L., & VandenBerg, D. A. 2014, MNRAS, 444, 392 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  22. Castelli, F., & Kurucz, R. L. 2004, ArXiv e-prints [arXiv:0405087] [Google Scholar]
  23. Chiavassa, A., Plez, B., Josselin, E., & Freytag, B. 2009, A&A, 506, 1351 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
  24. Chiavassa, A., Collet, R., Casagrande, L., & Asplund, M. 2010, A&A, 524, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Chiavassa, A., Bigot, L., Thévenin, F., et al. 2011, J. Phys. Conf. Ser., 328, 012012 [NASA ADS] [CrossRef] [Google Scholar]
  26. Chiavassa, A., Bigot, L., Kervella, P., et al. 2012, A&A, 540, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Chiavassa, A., Ligi, R., Magic, Z., et al. 2014a, A&A, 567, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Chiavassa, A., Thévenin, F., Magic, Z., Collet, R., & Asplund, M. 2014b, 3D spectral library for RVS radial velocities, GAIA-C8-TN-OCA-AC-001-1, Tech. Rep. [Google Scholar]
  29. Chiavassa, A., Thévenin, F., Magic, Z., Collet, R., & Asplund, M. 2014c, The 3D spectral library for BP/RP, GAIA-C8-TN-OCA-AC-002-1, Tech. Rep. [Google Scholar]
  30. Chiavassa, A., Pere, C., Faurobert, M., et al. 2015, A&A, 576, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Chiavassa, A., Caldas, A., Selsis, F., et al. 2017, A&A, 597, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Cohen, M., Wheaton, W. A., & Megeath, S. T. 2003, AJ, 126, 1090 [NASA ADS] [CrossRef] [Google Scholar]
  33. Collet, R., Asplund, M., & Trampedach, R. 2007, A&A, 469, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Creevey, O. L., Thévenin, F., Boyajian, T. S., et al. 2012, A&A, 545, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Deustua, S., Baggett, S., Brammer, G., et al. 2016, WFC3 Data Handbook. Version 3.0 (Baltimore: STScI) Tech. Rep. [Google Scholar]
  36. Doi, M., Tanaka, M., Fukugita, M., et al. 2010, AJ, 139, 1628 [NASA ADS] [CrossRef] [Google Scholar]
  37. Dravins, D. 1982, ARA&A, 20, 61 [NASA ADS] [CrossRef] [Google Scholar]
  38. Dravins, D., Lindegren, L., & Nordlund, A. 1981, A&A, 96, 345 [NASA ADS] [Google Scholar]
  39. Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Gray, D. F. 2009, ApJ, 697, 1032 [NASA ADS] [CrossRef] [Google Scholar]
  41. Gray, D. F. 2010, ApJ, 721, 670 [NASA ADS] [CrossRef] [Google Scholar]
  42. Gray, R. O., Graham, P. W., & Hoyt, S. R. 2001, AJ, 121, 2159 [NASA ADS] [CrossRef] [Google Scholar]
  43. Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332 [NASA ADS] [CrossRef] [Google Scholar]
  44. Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Hayek, W., Asplund, M., Carlsson, M., et al. 2010, A&A, 517, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Hayek, W., Sing, D., Pont, F., & Asplund, M. 2012, A&A, 539, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Jordi, C., Gebran, M., Carrasco, J. M., et al. 2010, A&A, 523, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Jørgensen, U. G. 1997, IAU Symp., 178, 441 [Google Scholar]
  49. Katz, D., Munari, U., Cropper, M., et al. 2004, MNRAS, 354, 1223 [NASA ADS] [CrossRef] [Google Scholar]
  50. Kurucz, R. L. 2005, Mem. Soc. Astron. It. Suppl., 8, 189 [Google Scholar]
  51. Lindegren, L., & Dravins, D. 2003, A&A, 401, 1185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Ludwig, H., Caffau, E., Steffen, M., et al. 2009, Mem. Soc. Astron. It., 80, 711 [NASA ADS] [Google Scholar]
  53. Magic, Z., Collet, R., Asplund, M., et al. 2013, A&A, 557, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Magic, Z., Collet, R., & Asplund, M. 2014, A&A, submitted [arXiv:1403.6245] [Google Scholar]
  55. Magic, Z., Chiavassa, A., Collet, R., & Asplund, M. 2015, A&A, 573, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Megessier, C. 1995, A&A, 296, 771 [NASA ADS] [Google Scholar]
  57. Mihalas, D., Dappen, W., & Hummer, D. G. 1988, ApJ, 331, 815 [NASA ADS] [CrossRef] [Google Scholar]
  58. Nordlund, A. 1982, A&A, 107, 1 [Google Scholar]
  59. Nordlund, Å., Stein, R. F., & Asplund, M. 2009, Liv. Rev. Sol. Phys., 6, 2 [Google Scholar]
  60. Palacios, A., Gebran, M., Josselin, E., et al. 2010, A&A, 516, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Pasquini, L., Melo, C., Chavero, C., et al. 2011, A&A, 526, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Pereira, T. M. D., Asplund, M., Collet, R., et al. 2013, A&A, 554, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Ramírez, I., Allende Prieto, C., & Lambert, D. L. 2008, A&A, 492, 841 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Ramírez, I., Allende Prieto, C., Koesterke, L., Lambert, D. L., & Asplund, M. 2009, A&A, 501, 1087 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Skartlien R. 2000, ApJ, 536, 465 [NASA ADS] [CrossRef] [Google Scholar]
  66. Stempels, H. C., Piskunov, N., & Barklem, P. S. 2001, in 11th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun, eds. R. J. Garcia Lopez, R. Rebolo, & M. R. Zapaterio Osorio, ASP Conf. Ser., 223, 878 [Google Scholar]
  67. Trampedach, R., Asplund, M., Collet, R., Nordlund, Å., & Stein, R. F. 2013, ApJ, 769, 18 [NASA ADS] [CrossRef] [Google Scholar]

4

The only difference with respect to Casagrande & VandenBerg (2014) is that here we have adopted MBol = 4.74.

6

For more details, see the technical note “The 3D spectral library for BP/RP” (Chiavassa et al. 2014c).

7

Defined as the locus of the midpoints of the spectral line.

8

For more details, see the technical note “3D spectral library for RVS radial velocities” (Chiavassa et al. 2014b) and the paper on CU6 design and performance (Sartoretti et al., in prep.).

9

Available at http://pollux.oreme.org

All Tables

Table 1

Photometric systems used in this work and overplotted on the synthetic spectra in Fig. 4.

Table 2

Central wavelength position, oscillator strength (log gf), and excitation potential (χ) for the 20 Fe I in the spectral domain of RVS (Bigot & Thévenin 2006).

Table 3

Selected convective shifts for Fe I only (CSStagger in km s−1 ) from RHD simulations in this work andfrom CIFIST-grid simulations (CSCIFIST in km s−1 ) with equivalent stellar parameters from Allende Prieto et al. (2013).

Table B.1

RHD simulations’ steAppendix Bllar parameters (first three cols.), bolometric Magnitude (Mbol ), and bolometric correction (BC) for Johnson-Cousins, 2MASS, SDSS, and Gaia systems (Table 1).

All Figures

thumbnail Fig. 1

Surface rendering for all the synthetic spectra computed for the 3D RHD simulations in Table B.1. The vertical bar on the right displays the colour scale for the emerging flux in erg s−1 cm−2−1. For clarity, the wavelength range has been reduced to 1000–25 000 Å.

In the text
thumbnail Fig. 2

Synthetic spectra of the solar simulation in the spectral range 2000–25 000 Å and for the different μ=cos(θ) $\mu=\cos(\theta)$ inclination angles used in the computation, where θ is the angle with respect to the line of sight (vertical axis).

In the text
thumbnail Fig. 3

Comparison of the solar simulation (black) with the observed flux of the Sun (red, Kurucz 2005). The solarirradiance is converted to flux at the solar surface using the multiplicative factor of [(1AU )R]2=46202$\left[\left(1 \rm{AU}\right)/R_{\odot}\right]^2 = 46\,202$. For clarity, the spectra have been resampled to a lower spectral resolution with fewer frequency points (nλ = 2115).

In the text
thumbnail Fig. 4

Synthetic spectrum of the solar simulation at full spectral resolution in the spectral range 2000–25 000 Å (grey, top panel) and 2000–12 000 Å (central and bottom panels). Several system response functions (Table 1), from which synthetic colours have been computed, are overplotted. Johnson-Cousins system response functions (U, B, V, Rc, Ic) are plotted in green and 2MASS in pink-violet (top left panel); SDSS (u, g, r, i, z) in yellow-red and Gaia (BP, RP, G) (top right panel); Strömgren (uvby) in red-blue and SkyMapper (us , vs , gs, rs, is, zs) (bottom left panel); and the 15 filters of the HST-WFC3 (bottom right panel). For clarity, SDSS and SkyMapper functions are normalised to 0.5.

In the text
thumbnail Fig. 5

Flux normalised to the continuum for the solar simulation (Table B.1) in the RVS range (8400–8900 Å). The Fe I lines (black arrows) are from Bigot & Thévenin (2006), while the Ca II triplet are indicated with red arrows.

In the text
thumbnail Fig. 6

Line bisectors of the 20 Fe I lines (black) from Bigot & Thévenin (2006) and Ca II triplet (red) for five 3D simulations in the grid.

In the text
thumbnail Fig. 7

Top left panel: convective shifts predicted by the 3D hydrodynamical simulations for the Gaia RVS spectral range and the Fe I lines (see text for details). Top right panel: convective shifts from Fe I lines as a function of the effective temperature of the 3D simulations in Table B.1. Bottom left panel: convective shifts from Ca II lines. Bottom right panel: comparison of convective shifts for a selected number of RHD simulations from this work (stars) with simulations with equivalent stellar parameters from Allende Prieto et al. (2013) (circles). See Table 3 for details.

In the text
thumbnail Fig. A.1

Bolometric correction (BC) differences computed for the photometric filters Johnson-Cousins and 2MASS with 3D simulations (Table B.1) and the corresponding 1D hydrostatic models with microturbulence = 1 km s−1: Δ BC = BC1DBC3D. The red dashed line indicates ΔBC = 0.0.

In the text
thumbnail Fig. A.2

Bolometric correction differences for the SkyMapper and Gaia photometric filters (Table 1). The notation is the same as in Fig. A.1.

In the text
thumbnail Fig. A.3

Bolometric correction differences for the HST-WFC3 provided in the VEGA system (for the ST and AB systems the differences are similar); first set of filters (Table 1). The notation is the same as in Fig. A.1.

In the text
thumbnail Fig. A.4

Bolometric correction differences for the HST-WFC3 provided in the VEGA system (for ST and AB systems the differences are similar); second set of filters (Table 1). The notation is the same as in Fig. A.1.

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.