The correct interpretation of the large amount of complex data from next-generation (in particular, space-based) observational facilities requires a very strong theoretical underpinning. One can thus predict that in the near future the use of atmospheric models obtained through the use of three-dimensional (3-D) radiation magnetohydrodynamics (RMHD) codes, coupled with advanced radiative transfer treatment including non-local thermodynamic equilibrium (non-LTE) effects and polarisation, will become the norm. In particular, stellar brightness variability in cool (i. e., spectral type F, G, K, and M) stars can be induced by several different effects, besides pulsation. We briefly discuss some literature results and mention some of our recent progress. Finally, we attempt to peek into the future of understanding this important aspect of the life of stars.

Thermal convection is one of the main mechanisms of heat transport and mixing in stars in general and also in the photospheric layers which emit the radiation that we observe with astronomical instruments. The present lecture notes first introduce the role of convection in astrophysics and explain the basic physics of convection. This is followed by an overview on the modelling of convection. Challenges and pitfalls in numerical simulation based modelling are discussed subsequently. Finally, a particular application for the previously introduced concepts is described in more detail: the study of convective overshooting into stably stratified layers around convection zones in stars.

We present a numerical study of double-diffusive convection characterized by a stratification unstable to thermal convection while at the same time a mean molecular weight (or solute concentration) difference between top and bottom counteracts this instability. Convective zones can form in this case either by the stratification being locally unstable to the combined action of both temperature and solute gradients or by another process, the oscillatory double-diffusive convective instability, which is triggered by the faster molecular diffusivity of heat in comparison with that one of the solute. We discuss successive layer formation for this problem in the case of an instantaneously heated bottom (plate) which forms a first layer with an interface that becomes temporarily unstable and triggers the formation of further, secondary layers. We consider both the case of a Prandtl number typical for water (oceanographic scenario) and of a low Prandtl number (giant planet scenario). We discuss the impact of a Couette like shear on the flow and in particular on layer formation for different shear rates. Additional layers form due to the oscillatory double-diffusive convective instability, as is observed for some cases. We also test the physical model underlying our numerical experiments by recovering experimental results of layer formation obtained in laboratory setups.

Space-borne missions CoRoT and Kepler have provided a rich harvest of high-quality photometric data for solar-like pulsators. It is now possible to measure damping rates for hundreds of main-sequence and thousands of red-giant. However, among the seismic parameters, mode damping rates remain poorly understood and thus barely used for inferring the physical properties of stars. Previous approaches to model mode damping rates were based on mixing-length theory or a Reynolds-stress approach to model turbulent convection. While able to grasp the main physics of the problem, those approaches are of little help to provide quantitative estimates as well as a definitive answer on the relative contribution of each physical mechanism. Our aim is thus to assess the ability of 3D hydrodynamical simulations to infer the physical mechanisms responsible for damping of solar-like oscillations. To this end, a solar high-spatial resolution and long-duration hydrodynamical 3D simulation computed with the ANTARES code allows probing the coupling between turbulent convection and the normal modes of the simulated box. Indeed, normal modes of the simulation experience realistic driving and damping in the super-adiabatic layers of the simulation. Therefore, investigating the properties of the normal modes in the simulation provides a unique insight into the mode physics. We demonstrate that such an approach provides constraints on the solar damping rates and is able to disentangle the relative contribution related to the perturbation of the turbulent pressure, the gas pressure, the radiative flux, and the convective flux contributions. Finally, we conclude that using the normal modes of a 3D numerical simulation is possible and is potentially able to unveil the respective role of the different physical mechanisms responsible for mode damping provided the time-duration of the simulation is long enough.

We present results of a large, high resolution 3D hydrodynamical simulation of the surface layers of a DA white dwarf (WD) with $T_{\rm eff}=11800$ K and $\log(g)=8$ using the ANTARES code, the widest and deepest such simulation to date. Our simulations are in good agreement with previous calculations in the Schwarzschild-unstable region and in the overshooting region immediately beneath it. Farther below, in the wave-dominated region, we find that the rms horizontal velocities decay with depth more rapidly than the vertical ones. Since mixing requires both vertical and horizontal displacements, this could have consequences for the size of the region that is well mixed by convection, if this trend is found to hold for deeper layers. We discuss how the size of the mixed region affects the calculated settling times and inferred steady-state accretion rates for WDs with metals observed in their atmospheres.

After an overview of general aspects of modelling the pulsation- convection interaction we present reasons why such simulations (in multidimensions) are needed but, at the same time, pose a considerable challenge. We then discuss, for several topics, what insights multidimensional simulations have either already provided or can be expected to yield in the future. We finally discuss properties of our ANTARES code. Many of these features can be expected to be characteristic of other codes which may possibly be applied to these physical questions in the foreseeable future.

Non-local, time-dependent convection models have been used to explain the location of double-mode pulsations in Cepheids in the HR diagram as well as the existence and location of the red edge of the instability strip. These properties are highly sensitive to model parameters. We use 2D radiation hydrodynamical simulations with realistic microphysics and grey radiative-transfer to model a short period Cepheid. The simulations show that the strength of the convection zone varies significantly over the pulsation period and exhibits a phase shift relative to the variations in radius. We evaluate the convective flux and the work integral as predicted by the most common convection models. It turns out that over one pulsation cycle the model parameter $\alpha_{\rm c}$, has to be varied by up to a factor of beyond 2 to match the convective flux obtained from the simulations. To bring convective fluxes integrated over the He II convection zone and the overshoot zone below into agreement, this parameter has to be varied by a factor of up to $\sim 7.5$ (Kuhfu{\ss}). We then present results on the energetics of the convection and overshoot zone by radially symmetric and fluctuating quantities. To successfully model this scenario by a static, one dimensional or even by a simple time-dependent model appears extremely challenging. We conclude that significant improvements are needed to make predictions based on 1D models more robust and to improve the reliability of conclusions on the convection-pulsation coupling drawn from them. Multidimensional simulations can provide guidelines for developing descriptions of convection then applied in traditional 1D modelling.

We investigate the achievable efficiency of both the time and the space discretisation methods used in Antares for mixed parabolic-hyperbolic problems. We show that the fifth order variant of WENO combined with a second order Runge-Kutta scheme is not only more accurate than standard first and second order schemes, but also more efficient taking the computation time into account. Then, we calculate the error decay rates of WENO with several explicit Runge-Kutta schemes for advective and diffusive problems with smooth and non-smooth initial conditions. With this data, we estimate the computational costs of three-dimensional simulations of stellar surface convection and show that SSP RK(3,2) is the most efficient scheme considered in this comparison.

So far, numerical studies of double-diffusive layering in turbulent convective flows have neglected the effects of rotation. We undertake a first step into that direction by investigating how Coriolis forces affect a double-diffusive layer inside a rotating spherical shell. For this purpose we have run simulations in a parameter regime where these layers are expected to form and successively increased the rate of rotation with the result that fast rotation is found to have a similar stabilising effect on the overall convective flux as an increase of the stability ratio $R_{\rho}$ has in a non-rotating setup. We have also studied to what extent the regimes of rotational constraints suggested by King, Stellmach, and Buffett (2013) for rotation in the case of Rayleigh-B\'enard convection are influenced by double-diffusive convection: their classification could also be applicable to the case of double-diffusive convection in a spherical shell if it is extended to be also a function of the stability ratio $R_{\rho}$. Furthermore, we examined the ratio of saline and thermal Nusselt numbers and compared our results with models of Spruit (2013), Rosenblum et al. (2011) and Wood, Garaud, and Stellmach (2013). We find our data to be fitted best by Spruit's model. Our result that fast rotation further decreases the convective transport, which is already lowered by double-diffusive convection, could play a major role for e.g. the modeling of the interior of some rapidly rotating giant planets, as gaseous giant planets have recently been proposed to be influenced by double-diffusive convection.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

We have implemented open boundary conditions into the ANTARES code to increase the realism of our simulations of stellar surface convection. Even though we greatly benefit from the high accuracy of our fifth order numerical scheme (WENO5), the broader stencils needed for the numerical scheme complicate the implementation of boundary conditions. We show that the effective temperature of a numerical simulation cannot be changed by corrections at the lower boundary since the thermal stratification does only change on the Kelvin-Helmholtz time scale. Except for very shallow models, this time scale cannot be covered by multidimensional simulations due to the enormous computational requirements. We demonstrate to what extent numerical simulations of stellar surface convection are sensitive to the initial conditions and the boundary conditions. An ill-conceived choice of parameters for the boundary conditions can have a severe impact. Numerical simulations of stellar surface convection will only be (physically) meaningful and realistic if the initial model, the extent and position of the simulation box, and the parameters from the boundary conditions are chosen adequately.

We investigate the applicability of curvilinear grids in the context of astrophysical simulations and WENO schemes. With the non-smooth mapping functions from Calhoun et al. (2008), we can tackle many astrophysical problems which were out of scope with the standard grids in numerical astrophysics. We describe the difficulties occurring when implementing curvilinear coordinates into our WENO code, and how we overcome them. We illustrate the theoretical results with numerical data. The WENO finite difference scheme works only for high Mach number flows and smooth mapping functions whereas the finite volume scheme gives accurate results even for low Mach number flows and on non-smooth grids.

We have extended the ANTARES code to simulate the coupling of pulsation with convection in Cepheid-like variables in an increasingly realistic way, in particular in multidimensions, 2D at this stage. Present days models of radially pulsating stars assume radial symmetry and have the pulsation-convection interaction included via model equations containing ad hoc closures and moreover parameters whose values are barely known. We intend to construct ever more realistic multidimensional models of Cepheids. In the present paper, the first of a series, we describe the basic numerical approach and how it is motivated by physical properties of these objects which are sometimes more, sometimes less obvious. - For the construction of appropriate models a polar grid co-moving with the mean radial velocity has been introduced to optimize radial resolution throughout the different pulsation phases. The grid is radially stretched to account for the change of spatial scales due to vertical stratification and a new grid refinement scheme is introduced to resolve the upper, hydrogen ionisation zone where the gradient of temperature is steepest. We demonstrate that the simulations are not conservative when the original weighted essentially non-oscillatory method implemented in ANTARES is used and derive a new scheme which allows a conservative time evolution. The numerical approximation of diffusion follows the same principles. Moreover, the radiative transfer solver has been modified to improve the efficiency of calculations on parallel computers. We show that with these improvements the ANTARES code can be used for realistic simulations of the convection-pulsation interaction in Cepheids. We discuss the properties of several models which include the upper 42% of a Cepheid along its radial coordinate, assume different opening angles, and are suitable for an in-depth study of convection and pulsation.

The Virtual Atomic and Molecular Data Centre (VAMDC) (M.L. Dubernet et al. 2010, JQSRT 111, 2151) is an EU-FP7 e-infrastructure project devoted to building a common electronic infrastructure for the exchange and distribution of atomic and molecular data. It involves two dozen teams from six EU member states (Austria, France, Germany, Italy, Sweden, United Kingdom) as well as Russia, Serbia, and Venezuela. Within VAMDC scientists from many different disciplines in atomic and molecular physics collaborate with users of their data and also with scientists and engineers from the information and communication technology community. In this presentation an overview of the current status of VAMDC and its capabilities will be provided. In the second part of the presentation I will focus on one of the databases which have become part of the VAMDC platform, the Vienna Atomic Line Data Base (VALD). VALD has developed into a well-known resource of atomic data for spectroscopy particularly in astrophysics. A new release, VALD-3, will provide numerous improvements over its predecessor. This particularly relates to the data contents where new sets of atomic data for both precision spectroscopy (i.e., with data for observed energy levels) as well as opacity calculations (i.e., with data involving predicted energy levels) have been included. Data for selected diatomic molecules have been added and a new system for data distribution and data referencing provides for more convenience in using the upcoming third release of VALD.

We report on modelling in stellar astrophysics with the ANTARES code. First, we describe properties of turbulence in solar granulation as seen in high-resolution calculations. Then, we turn to the first 2D model of pulsation-convection interaction in a cepheid. We discuss properties of the outer and the HEII ionization zone. Thirdly, we report on our work regarding models of semiconvection in the context of stellar physics.

We report on the first asteroseismic analysis of solar-type stars observed by Kepler. Observations of three G-type stars, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation can clearly be distinguished in each star. We discuss the appearance of the oscillation spectra, including the presence of a possible signature of faculae, and the presence of mixed modes in one of the three stars.

We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe.

Convection is the most important physical process that determines the structure of the envelopes of cool stars. It influences the surface radiation flux and the shape of observed spectral line profiles and is responsible for both generating and damping solar-like oscillations, among others. 3D numerical simulations of stellar surface convection have developed into a powerful tool to model and analyse the physical mechanisms operating at the surface of cool stars. This review discusses the main principles of 3D stellar atmospheres used for such applications. The requirements from stellar structure and evolution theory to use them as boundary conditions are analysed as well as the capabilities of using helio- and asteroseismology to reduce modelling uncertainties and probing the consistency and accuracy of 3D stellar atmospheres as part of this process. Simulations for the solar surface made by different teams are compared and some issues concerning the uncertainties of this modelling approach are discussed.

The line profiles of the stars with v sin i below a few km/s can reveal direct signatures of local velocity fields (e.g. convection) in stellar atmospheres. This effect is well established in cool main sequence stars, and has been detected and studied in three A stars. This paper reports observations of main sequence B, A and F stars with two goals: (1) to identify additional stars having sufficiently low values of v sin i to search for spectral line profile signatures of local velocity fields, and (2) to explore how the signatures of the local velocity fields in the atmosphere depend on stellar parameters such as effective temperature T_eff and peculiarity type. For stars having T_eff below about 10000 K, we always detect local atmospheric velocity fields indirectly through a non-zero microturbulence parameter, but not for hotter stars. Among the A and F stars in our sample having the sharpest lines, direct tracers of atmospheric velocity fields are found in six new stars. The velocity field signatures identified include asymmetric excess line wing absorption, deeper in the blue line wing than in the red; line profiles of strong lines that are poorly fit by computed profiles; and strong lines that are broader than they should be for the v sin i values deduced from weak lines. These effects are found in both normal and Am stars, but seem stronger in Am stars. These data still have not been satisfactorily explained by models of atmospheric convection, including numerical simulations.

We discuss the general design of the ANTARES code which is intended for simulations in stellar hydrodynamics with radiative transfer and realistic microphysics in 1D, 2D and 3D. We then compare the quality of various numerical methods. We have applied ANTARES in order to obtain high resolution simulations of solar granulation which we describe and analyze. In order to obtain high resolution, we apply grid refinement to a region predominantly occupied by an exploding granule. Strong, rapidly rotating vortex tubes of small diameter (~100 km) generated by the downdrafts and ascending into the photosphere near the granule boundaries evolve, often entering the photosphere from below in an arclike fashion. They essentially contribute to the turbulent velocity field near the granule boundaries.

We present results from 2D radiation-hydrodynamical simulations of fully compressible convection for the surface layers of A-type stars with the ANTARES code. Spectroscopic indicators for photospheric convective velocity fields show a maximum of velocities near Teff ~8000 K. In that range the largest values are measured for the subgroup of Am stars. Thus far, no prognostic model, neither theoretical nor numerical, is able to exactly reproduce the line profiles of sharp line A and Am stars in that temperature range. In general, the helium abundance of A stars is not known from observations. Hence, we have considered two extreme cases for our simulations: a solar helium abundance as an upper limit and zero helium abundance as a lower limit. The simulation for the helium free case is found to differ from the case with solar helium abundance by larger velocities, larger flow structures, and by a sign reversal of the flux of kinetic energy inside the hydrogen ionisation zone. Both simulations show extended shock fronts emerging from the optical surface, as well as mixing far below the region of partial ionisation of hydrogen, and vertical oscillations emerging after initial perturbations have been damped. We discuss problems related to the rapid radiative cooling at the surface of A-type stars such as resolution and efficient relaxation. The present work is considered as a step towards a systematic study of convection in A- to F-type stars, encouraged by the new data becoming available for these objects from both asteroseismological missions and from high resolution spectroscopy.

We infer from different seismic observations the energy supplied per unit of time by turbulent convection to the acoustic modes of Alpha Cen A (HD 128620), a star which is similar but not identical to the Sun. The inferred rates of energy supplied to the modes (i.e. mode excitation rates) are found to be significantly larger than in the Sun. They are compared with those computed with an excitation model that includes two sources of driving, the Reynolds stress contribution and the advection of entropy fluctuations. The model also uses a closure model, the Closure Model with Plumes (CMP hereafter), that takes the asymmetry between the up- and down-flows (i.e. the granules and plumes, respectively) into account. Different prescriptions for the eddy-time correlation function are also confronted to observational data. Calculations based on a Gaussian eddy-time correlation underestimate excitation rates compared with the values derived from observations for Alpha Cen A. On the other hand, calculations based on a Lorentzian eddy-time correlation lie within the observational error bars. This confirms results obtained in the solar case. With respect to the helioseismic data, those obtained for Alpha Cen A constitute an additional support for our model of excitation. We show that mode masses must be computed taking turbulent pressure into account. Finally, we emphasize the need for more accurate seismic measurements in order to discriminate, in the case of Alpha Cen A, between the CMP closure model and the quasi-Normal Approximation as well as to confirm or not the need to include the excitation by the entropy fluctuations.

Using grid refinement, we have simulated solar granulation in 2D. The refined region measures 1.97*2.58 Mm (vertical*horizontal). Grid spacing there is 1.82*2.84 km. The downflows exhibit strong Kelvin-Helmholtz instabilities. Below the photosphere, acoustic pulses are generated. They proceed laterally (in some cases distances of at least the size of our refined domain) and may be enhanced when transversing downflows) as well as upwards where, in the photosphere they contribute significantly to 'turbulence' (velocity gradients, etc.) The acoustic pulses are ubiquitous in that at any time several of them are seen in our high-resolution domain. Their possible contributions to p-mode excitation or heating of the chromosphere needs to be investigated.

Oscillations of stellar p modes, excited by turbulent convection, are investigated. We take into account the asymmetry of the up and downflows created by turbulent plumes through an adapted closure model. In a companion paper, we apply it to the formalism of excitation of solar p modes developed by Samadi & Goupil 2001. Using results from 3D numerical simulations of the upper most part of the solar convection zone, we show that the two-scale-mass-flux model (TFM) is valid only for quasi-laminar or highly skewed flows (Gryanik & Hartmann 2002). We build a generalized-Two-scale-Mass-Flux Model (GTFM) model which takes into account both the skew introduced by the presence of two flows and the effects of turbulence in each flow. In order to apply the GTFM to the solar case, we introduce the plume dynamics as modelled by Rieutord & Zahn (1995) and construct a Closure Model with Plumes (CMP). When comparing with 3D simulation results, the CMP improves the agreement for the fourth order moments, by approximatively a factor of two compared with the use of the quasi-normal approximation or a skewness computed with the classical TFM. The asymmetry of turbulent convection in the solar case has an important impact on the vertical-velocity fourth-order moment which has to be accounted for by models. The CMP is a significant improvement and is expected to improve the modelling of solar p-mode excitation.

Excitation of solar-like oscillations is attributed to turbulent convection and takes place at the upper-most part of the outer convective zones. Amplitudes of these oscillations depend on the efficiency of the excitation processes as well as on the properties of turbulent convection. We present past and recent improvements on the modeling of those processes. We show how the mode amplitudes and mode line-widths can bring information about the turbulence in the specific cases of the Sun and Alpha Cen A.

Our goal is to improve the theoretical modelling of stochastic excitation of p modes by turbulent convection. With the help of the closure model with plume (CMP) developed in a companion paper, we refine the theoretical description of the excitation by the turbulent Reynolds stress term. The CMP is generalized for two-point correlation products so as to apply it to the formalism developed by Samadi & Goupil (2001). The excitation source terms are then computed with this improvement, and a comparison with solar data from the GOLF instrument is performed. The present model provides a significant improvement when comparing absolute values of theoretical ampplitudes with observational data. It gives rise to a frequency dependence of the power supplied to solar p modes, which agrees with GOLF observations. It is shown that the asymmetry of the turbulent convection zone (up- and downflows) plays a major role in the excitation processes. Despite an increase in the Reynolds stress term contribution due to our improved description, an additional source of excitation, identified as the entropy source term, is still necessary for reproducing the observational data. Theoretical excitation rates in the frequency range [2.5 mHz, 4 mHz] now are in agreement with the observational data from the GOLF instrument. However, at lower frequencies, it exhibits small discrepancies at the maximum level of a few per cent. Improvements are likely to come from a better physical description of the excitation by entropy fluctuations in the superadiabatic zone.

For applications in population synthesis, libraries of theoretical stellar spectra are often considered an alternative to template libraries of observed spectra, because they allow a complete sampling of stellar parameters. Most attention in published theoretical spectral libraries has been devoted to the visual wavelength range. We present a detailed comparison of theoretical spectra in the range 1.57-1.67$\mu$m, for spectral types from A to early M and for giants and dwarf stars, with observed stellar spectra at resolutions around 3000, which would be sufficient to disentangle the different groups of late type stars. We have selected the NeMo grids of stellar atmospheres to perform such a comparison. We first demonstrate that after combining atomic and molecular line lists, it is possible to match observed spectral flux distributions with theoretical ones very well for almost the entire parameter range covered by the NeMo grids at moderate resolution in the visual range. In the infrared range, although the overall shape of the observed flux distributions is still matched reasonably well, the individual spectral features are reproduced by the theoretical spectra only for stars earlier than mid F type. For later spectral types the differences increase and theoretical spectra of K type stars have systematically weaker line features than those found in observations. These discrepancies are traced back to stem primarily from incomplete data on neutral atomic lines, although some of them are also related to molecules. Improving atomic data in the near infrared is a key element in making the construction of reliable libraries of stellar spectra in the infrared feasible.

We compute the rates P at which acoustic energy is injected into the solar radial p modes for several solar models. The solar models are computed with two different local treatments of convection: the classical mixing-length theory (MLT hereafter) and Canuto et al (1996)'s formulation (CGM hereafter). Among the models investigated here, our best models reproduce both the solar radius and the solar luminosity at solar age and the observed Balmer line profiles. For the MLT treatment, the rates P do depend significantly on the properties of the atmosphere whereas for the CGM's treatment the dependence of P on the properties of the atmosphere is found smaller than the error bars attached to the seismic measurements. The excitation rates P for modes associated with the MLT models are significantly underestimated compared with the solar seismic constraints. The CGM models yield values for P closer to the seismic data than the MLT models. We conclude that the solar p-mode excitation rates provide valuable constraints and according to the present investigation clearly favor the CGM treatment with respect to the MLT, although neither of them yields values of P as close to the observations as recently found for 3D numerical simulations.

We compute, for several solar models, the rates P at which the solar radial p modes are expected to be excited. The solar models are computed with two different local treatments of convection : the classical mixing-length theory (MLT hereafter) and Canuto, Goldmann and Mazzitelli(1996, CGM hereafter)'s formulation. For one set of solar models (EMLT and ECGM models), the atmosphere is gray and assumes Eddington's approximation. For a second set of models (KMLT and KCGM models), the atmosphere is built using a T(tau) law which has been obtained from a Kurucz's model atmosphere computed with the same local treatment of convection. The mixing-length parameter in the model atmosphere is chosen so as to provide a good agreement between synthetic and observed Balmer line profiles, while the mixing-length parameter in the interior model is calibrated so that the model reproduces the solar radius at solar age. For the MLT treatment, the rates P do depend significantly on the properties of the atmosphere. On the other hand, for the CGM treatment, differences in P between the ECGM and the KCGM models are very small compared to the error bars attached to the seismic measurements. The excitation rates P for modes from the EMLT model are significantly under-estimated compared with the solar seismic constraints. The KMLT model results in intermediate values for P and shows also an important discontinuity in the temperature gradient and the convective velocity. On the other hand, the KCGM model and the ECGM model yield values for P closer to the seismic data than the EMLT and KMLT models. We conclude that the solar p-mode excitation rates provide valuable constraints and according to the present investigation cleary favor the CGM treatment with respect to the MLT.

NEMO.2003 (vienna NEw MOdel grid of stellar atmospheres) is a DVD with 91,520 stellar model atmospheres representing a 5D grid of modified ATLAS9 atmospheres; the purpose of the modifications was to include different treatments of convection and higher vertical resolution. In addition, for every model fluxes are provided and color indices for 14 different photometric systems. Because the model grid contained gaps due to non-converging models, we developed and applied a set of 4D interpolation routines to complete the grid. All the data, which will continuously be updated, can be found on the NEMO homepage (http://ams.astro.univie.ac.at/nemo/) and are available via DVD.

We present results of a fully non-local model of convection for white dwarf envelopes. We show that this model is able to reproduce the results of numerical simulations for convective efficiencies ranging from very inefficient to moderately efficient; this agreement is made more impressive given that no closure parameters have been adjusted in going from the previously reported case of A-stars to the present case of white dwarfs; for comparison, in order to match the peak convective flux found in numerical simulations for both the white dwarf envelopes discussed in this paper and the A-star envelopes discussed in our previous work requires changing the mixing length parameter of commonly used local models by a factor of 4. We also examine in detail the overshooting at the base of the convection zone, both in terms of the convective flux and in terms of the velocity field: we find that the flux overshoots by approximately 1.25 H_P and the velocity by approximately 2.5 H_P. Due to the large amount of overshooting found at the base of the convection zone the new model predicts the mixed region of white dwarf envelopes to contain at least 10 times more mass than local mixing length theory (MLT) models having similar photospheric temperature structures. This result is consistent with the upper limit given by numerical simulations which predict an even larger amount of mass to be mixed by convective overshooting. Finally, we attempt to parametrise some of our results in terms of local MLT-based models, insofar as is possible given the limitations of MLT

The Teff location of Pre-Main Sequence (PMS) evolutionary tracks depends on the treatment of over-adiabaticity. Since the convection penetrates into the stellar atmosphere, also the treatment of convection in the modeling of stellar atmospheres will affect the location of the Hayashi tracks. We present new non-grey PMS tracks for Teff>4000 K. We compute several grids of evolutionary tracks varying: 1. the treatment of convection: either the Mixing Length Theory (MLT) or Canuto et al. (1996, CGM) treatment; 2.the atmospheric boundary conditions: we use the new Vienna grids of ATLAS9 atmospheres (Heiter et al. 2002a), which were computed using either MLT (with alpha=Lambda/Hp=0.5) or CGM treatments. We compute as well grids of models with the NextGen (Allard & Hauschildt 1997) atmosphere models, and a 1 Msun grey MLT evolutionary track using the alpha calibration based on 2D-hydrodynamical models (Ludwig et al. 1999)}. These different grids of models allow us to analyze the effects of convection modeling on the non--grey PMS evolutionary tracks. We conclude that: 1. In spite of the solar calibration, if MLT convection is adopted a large uncertainty results in the shape and location of PMS tracks. 2. As long as the model of convection is not the same in the interior and in the atmosphere the optical depth at which we take the boundary conditions is an additional parameter of the models. 3. The comparison between NextGen and ATLAS9 based models shows that, in the Teff domain they have in common (4000--10000K), the PMS location is mainly determined by the treatment of the over-adiabatic convection. 4. Comparison with observations indicates that, for any convection model and for any of the atmosphere grids, convection in the PMS phase appears to be less efficient than what is necessary in order to fit the Sun.

The Delta a photometric system provides an efficient observational method to identify and distinguish magnetic and several other types of chemically peculiar (CP) stars of spectral types B to F from other classes of stars in the same range of effective temperatures. We have developed a synthetic photometric system which can be used to explore the capability of model atmospheres with individual element abundances to predict photometric Delta a magnitudes which measure the extent of the flux depression around 5200A found in different types of CP stars. In this first paper, we confirm the observed dependency of the a-index as a function of various colour indices sensitive to the effective temperature of stars as well as its average scatter expected from surface gravity variations within the main sequence band. The behaviour of the so-called normality line of Delta a systems used in photometric observations of CP stars is well reproduced. The metallicity dependence of the normality line of the Delta a system was computed for several grids of model atmospheres where the abundances of elements heavier than He had been scaled +-0.5 dex with respect to the solar value. We estimate a lowering of Delta a magnitudes for CP stars within the Magellanic Clouds by approximately -3 mmag relative to those in the solar neighbourhood assuming an average metallicity of [Fe/H]= -0.5 dex. Using these results on the metallicity bias of the Delta a system we find the observational systems in use suitable to identify CP stars in other galaxies or distant regions of our own galaxy and capable to provide data samples on a statistically meaningful basis. In turn, the synthetic system is suitable to test the performance of model atmospheres for CP stars. This work will be presented in follow-up papers of this series.

In previous work, Gardiner et al. (1999) found evidence for a discrepancy between the Teff obtained from Balmer lines with that from photometry and fundamental values for A-type stars. An investigation into this anomaly is presented using Balmer line profiles of stars in binary system with fundamental values of both Teff and log g. A revision of the fundamental parameters for binary systems given by Smalley & Dworetsky (1995) is also presented. The Teff obtained by fitting Halpha and Hbeta line profiles is compared to the fundamental values and those obtained from uvby photometry. We find that the discrepancy found by Gardiner et al. (1999) for stars in the range 7000 K < Teff < 9000 K is no longer evident.

Balmer lines are an important diagnostic of stellar atmospheric structure, since they are formed at a wide range of depths within the atmosphere. The different Balmer lines are formed at slightly different depths making them useful atmospheric diagnostics. The low sensitivity to surface gravity for stars cooler than ~8000 K makes them excellent diagnostics in the treatment of atmospheric convection. For hotter stars Balmer profiles are sensitive to both effective temperature and surface gravity. Provided we know the surface gravity of these stars from some other method (e.g. from eclipsing binary systems), we can use them to determine effective temperature. In previous work, we have found no significant systematic problems with using uvby photometry to determine atmospheric parameters of fundamental (and standard) stars. In fact, uvby was found to be very good for obtaining both Teff and log g. Using Halpha and Hbeta profiles, we have found that both the Canuto & Mazzitteli and standard Kurucz mixing-length theory without approximate overshooting are both in agreement to within the uncertainties of the fundamental stars. Overshooting models were always clearly discrepant. Some evidence was found for significant disagreement between all treatments of convection and fundamental values around 8000~9000K, but these results were for fundamental stars without fundamental surface gravities. We have used stars with fundamental values of both Teff and log g to explore this region in more detail.

We present several new sets of grids of model stellar atmospheres computed with modified versions of the ATLAS9 code. Each individual set consists of several grids of models with different metallicities ranging from [M/H] = -2.0 to +1.0 dex. The grids range from 4000 to 10000 K in T_eff and from 2.0 to 5.0 dex in logg. The individual sets differ from each other and from previous ones essentially in the physics used for the treatment of the convective energy transport, in the higher vertical resolution of the atmospheres and in a finer grid in the (T_eff, logg) plane. These improvements enable the computation of derivatives of color indices accurate enough for pulsation mode identification. In addition, we show that the chosen vertical resolution is necessary and sufficient for the purpose of stellar interior modelling. To explain the physical differences between the model grids we provide a description of the currently available modifications of ATLAS9 according to their treatment of convection. Our critical analysis of the dependence of the atmospheric structure and observable quantities on convection treatment, vertical resolution and metallicity reveals that spectroscopic and photometric observations are best represented when using an inefficient convection treatment. This conclusion holds whatever convection formulation investigated here is used, i.e. MLT(alpha=0.5), CM and CGM are equivalent. We also find that changing the convection treatment can lead to a change in the effective temperature estimated from Stroemgren color indices from 200 to 400 K.

We present results of a fully non-local, compressible model of convection for A-star envelopes. This model quite naturally reproduces a variety of results from observations and numerical simulations which local models based on a mixing length do not. Our principal results, which are for models with Teff between 7200 K and 8500 K, are the following: First, the photospheric velocities and filling factors are in qualitative agreement with those derived from observations of line profiles of A-type stars. Second, the He II and H I convection zones are separated in terms of convective flux and thermal interaction, but joined in terms of the convective velocity field, in agreement with numerical simulations. In addition, we attempt to quantify the amount of overshooting in our models at the base of the He II convection zone.

The non-local hydrodynamic moment equations for compressible convection are compared to numerical simulations. Convective and radiative flux typically deviate less than 20% from the 3D simulations, while mean thermodynamic quantities are accurate to at least 2% for the cases we have investigated. The moment equations are solved in minutes rather than days on standard workstations. We conclude that this convection model has the potential to considerably improve the modelling of convection zones in stellar envelopes and cores, in particular of A and F stars.

We report the first identification of the Eu III lambda 6666.317 line in optical spectra of CP stars. This line is clearly present in the spectra of HR 4816, 73 Dra, HR 7575, and beta CrB, while it is marginally present or absent in spectra of the roAp stars alpha Cir, gamma Equ, BI Mic, 33 Lib, and HD 24712.

Some of the rapidly oscillating (CP2) stars, have frequencies which are larger than the theoretical acoustic cut-off frequency. As the cut-off frequency depends on the T(tau) relation in the atmosphere, we have computed models and adiabatic frequencies for pulsating Ap stars with T(tau) laws based on Kurucz model atmospheres and on Hopf's purely radiative relation. The frequency-dependent treatment of radiative transfer as well as an improved calculation of the radiative pressure in Kurucz model atmospheres increase the theoretical acoustic cut-off frequency by about 200 microHz, which is closer to the observations. For alpha Cir we find models with Kurucz atmospheres which have indeed a cut-off frequency beyond the largest observed frequency and which are well within the Teff - L error box. For HD 24712 only models which are hotter by about 100 K and less luminous by nearly 10% than what is actually the most probable value would have an acoustic cut-off frequency large enough. One may thus speculate that the old controversy about a mismatch between observed largest frequencies and theoretical cut-off frequencies of roAp star models is resolved. However, the observational errors for the astrophysical fundamental parameters have to be reduced further and the model atmospheres refined. Further details can be found in Audard et al. (1997)

We present a discussion on the effects of convection on the uvby colours of A and F stars. The mixing-length theory used in ATLAS9 is compared to the turbulent convection theory of Canuto & Mazzitelli (1991, 1992). Comparison with fundamental stars reveals that colours calculated using the Canuto & Mazzitelli convection theory are generally in better agreement than those obtained using mixing-length theory.

We describe a new method for the computation of opacity distribution functions (ODFs) useful to calculate one-dimensional model atmospheres in local thermal equilibrium (LTE). The new method is fast enough to be applied on current workstations and allows the computation of model atmospheres which deviate significantly from (scaled) solar chemical composition. It has reproduced existing ODFs and model atmospheres for solar abundances. Depending on the type of chemical peculiarity the ``individual'' model atmosphere may have a structure and surface fluxes similar to atmospheres based on (scaled) solar abundances or deviate in a way that cannot be reproduced by any of the conventional models. Examples are given to illustrate this behavior. The availability of models with ``individualized'' abundances is crucial for abundance analyses and Doppler imaging of extreme CP stars.

Some of the rapidly oscillating (roAp) stars, have frequencies which are larger than the theoretical acoustic cut-off frequency. We have computed models and adiabatic frequencies for pulsating Ap stars with T(tau) laws based on Kurucz model atmospheres and on the Hopf's purely radiative relation. We compare the values of the cut-off frequency derived from expressions for the potential from Vorontsov & Zarkhov (1989), from Gough (1986), and from the approximation of an isothermal atmosphere. Kurucz model atmospheres increase the theoretical acoustic cut-off frequency for roAp stars by about 200 mu$Hz, which is closer to the observations. We compare our computations with observations of the 3 roAp stars HD 24712, alphaCir, and HD 134214. For alpha Cir we find models with Kurucz atmospheres which have indeed a cut-off frequency beyond the largest observed frequency. For HD 24712 only models which are hotter by about 100 K and less luminous by nearly 10% than what is actually the most probable value would have an acoustic cut-off frequency large enough. One may thus speculate that the old controversy about a mismatch between observed largest frequencies and theoretical cut-off frequencies of roAp star models is resolved. However, the model atmospheres have to be refined by investigating NLTE effects, among others, before a definite conclusion can be drawn.