Linear damping rates and modal frequency corrections of radial oscillation modes in selected LEGACY main-sequence stars are estimated by means of a nonadiabatic stability analysis. The selected stellar sample covers stars observed by Kepler with a large range of surface temperatures and surface gravities. A nonlocal, time-dependent convection model is perturbed to assess stability against pulsation modes. The mixing-length parameter is calibrated to the surface-convection-zone depth of a stellar model obtained from fitting adiabatic frequencies to the LEGACY observations, and two of the nonlocal convection parameters are calibrated to the corresponding LEGACY linewidth measurements. The remaining nonlocal convection parameters in the 1D calculations are calibrated so as to reproduce profiles of turbulent pressure and of the anisotropy of the turbulent velocity field of corresponding 3D hydrodynamical simulations. The atmospheric structure in the 1D stability analysis adopts a temperature-optical-depth relation derived from 3D hydrodynamical simulations. Despite the small number of parameters to adjust, we find good agreement with detailed shapes of both turbulent pressure profiles and anisotropy profiles with depth, and with damping rates as a function of frequency. Furthermore, we find the absolute modal frequency corrections, relative to a standard adiabatic pulsation calculation, to increase with surface temperature and surface gravity.

We present the discovery of TOI-197.01, the first transiting planet identified by the Transiting Exoplanet Survey Satellite (TESS) for which asteroseismology of the host star is possible. TOI-197 (HIP116158) is a bright (V=8.2 mag), spectroscopically classified subgiant which oscillates with an average frequency of about 430 muHz and displays a clear signature of mixed modes. The oscillation amplitude confirms that the redder TESS bandpass compared to Kepler has a small effect on the oscillations, supporting the expected yield of thousands of solar-like oscillators with TESS 2-minute cadence observations. Asteroseismic modeling yields a robust determination of the host star radius (2.943+/-0.064 Rsun), mass (1.212 +/- 0.074 Msun) and age (4.9+/-1.1 Gyr), and demonstrates that it has just started ascending the red-giant branch. Combining asteroseismology with transit modeling and radial-velocity observations, we show that the planet is a "hot Saturn" (9.17+/-0.33 Rearth) with an orbital period of ~14.3 days, irradiance of 343+/-24 Fearth, moderate mass (60.5 +/- 5.7 Mearth) and density (0.431+/-0.062 gcc). The properties of TOI-197.01 show that the host-star metallicity - planet mass correlation found in sub-Saturns (4-8 Rearth) does not extend to larger radii, indicating that planets in the transition between sub-Saturns and Jupiters follow a relatively narrow range of densities. With a density measured to ~15%, TOI-197.01 is one of the best characterized Saturn-sized planets to date, augmenting the small number of known transiting planets around evolved stars and demonstrating the power of TESS to characterize exoplanets and their host stars using asteroseismology.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

Availability of reliable atomic and molecular opacity tables is essential in a wide variety of astronomical modeling: the solar and stellar interiors, stellar and planetary atmospheres, stellar evolution, pulsating stars, and protoplanetary disks, to name a few. With the advancement of powerful research techniques such as helio-seismology and asteroseismology, solar neutrino-flux measurements, exoplanet survey satellites, three-dimensional hydrodynamic atmospheric simulations (including non-LTE and granulation effects), high-performance computing of atomic and molecular data, and innovative plasma experiments the accuracy and completeness of opacity tables is being taken to an unprecedented level. The goal of the second Workshop on Astrophysical Opacities was to gather opacity data producers and consumers from both the atomic and molecular sectors to contribute to solving outstanding problems and to develop more effective and integrated interfaces. In this review we attempt to summa- rize the discussion at the workshop and propose future directions for opacity research.

At the Workshop on Astrophysical Opacities, several attendees voiced their interest in a list of absorption data that are missing from or inadequate in current models of astrophysical objects. This wish list by modelers is meant as motivation and inspiration for experimentalists and theoreticians alike.

Is the Sun likely to have a more opaque interior than previously thought? The solar oxygen (or abundance) problem can be solved with higher interior opacities, reconciling abundance analyses based on 3D convective atmospheres with the helioseismic structure of the solar interior. This has been known for more than a decade, but last year we learned that the absorption by just iron may contribute 7\% more to the solar opacity at the bottom of the convection zone than predicted by any opacity calculation so far, and by OP05 in particular. I find that artificial changes to the absorption (calibrated against the iron experiment) by other elements in a solar mixture give an opacity increase of a shape and magnitude that can restore agreement between modern abundance analysis and helioseismology. This suggests that improved opacity calculations will solve the solar oxygen problem.

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

The effect of metallicity on the granulation activity in stars is still poorly understood. Available spectroscopic parameters from the updated APOGEE-\textit{Kepler} catalog, coupled with high-precision photometric observations from NASA's \textit{Kepler} mission spanning more than four years of observation, make oscillating red giant stars in open clusters crucial testbeds. We determine the role of metallicity on the stellar granulation activity by discriminating its effect from that of different stellar properties such as surface gravity, mass, and temperature. We analyze 60 known red giant stars belonging to the open clusters NGC 6791, NGC 6819, and NGC 6811, spanning a metallicity range from [Fe/H] $\simeq -0.09$ to $0.32$. The parameters describing the granulation activity of these stars and their $\nu_\mathrm{max}$, are studied by considering the different masses, metallicities, and stellar evolutionary stages. We derive new scaling relations for the granulation activity, re-calibrate existing ones, and identify the best scaling relations from the available set of observations. We adopted the Bayesian code DIAMONDS for the analysis of the background signal in the Fourier spectra of the stars. We performed a Bayesian parameter estimation and model comparison to test the different model hypotheses proposed in this work and in the literature. Metallicity causes a statistically significant change in the amplitude of the granulation activity, with a dependency stronger than that induced by both stellar mass and surface gravity. We also find that the metallicity has a significant impact on the corresponding time scales of the phenomenon. The effect of metallicity on the time scale is stronger than that of mass. A higher metallicity increases the amplitude of granulation and meso-granulation signals and slows down their characteristic time scales toward longer periods.

We analyse the effect on adiabatic stellar oscillation frequencies of replacing the near-surface layers in 1D stellar structure models with averaged 3D stellar surface convection simulations. The main difference is an expansion of the atmosphere by 3D convection, expected to explain a major part of the asteroseismic surface effect; a systematic overestimation of p-mode frequencies due to inadequate surface physics. We employ pairs of 1D stellar envelope models and 3D simulations from a previous calibration of the mixing-length parameter, alpha. That calibration constitutes the hitherto most consistent matching of 1D models to 3D simulations, ensuring that their differences are not spurious, but entirely due to the 3D nature of convection. The resulting frequency shift is identified as the structural part of the surface effect. The important, typically non-adiabatic, modal components of the surface effect are not included in the present analysis, but relegated to future papers. Evaluating the structural surface effect at the frequency of maximum mode amplitude, $\nu_{\rm max}$, we find shifts from $\delta\nu$=-0.8 microHz for giants at $\log g$=2.2 to -35 microHz for a ($T_{\rm eff}=6901$ K, $\log g$=4.29) dwarf. The fractional effect $\delta\nu(\nu_{\rm max})/\nu_{\rm max}$, ranges from -0.1% for a cool dwarf (4185 K, 4.74) to -6% for a warm giant (4962 K, 2.20).

Stellar evolution codes play a major role in present-day astrophysics, yet they share common issues. In this work we seek to remedy some of those by the use of results from realistic and highly detailed 3D hydrodynamical simulations of stellar atmospheres. We have implemented a new temperature stratification extracted directly from the 3D simulations into the Garching Stellar Evolution Code to replace the simplified atmosphere normally used. Secondly, we have implemented the use of a variable mixing-length parameter, which changes as a function of the stellar surface gravity and temperature -- also derived from the 3D simulations. Furthermore, to make our models consistent, we have calculated new opacity tables to match the atmospheric simulations. Here, we present the modified code and initial results on stellar evolution using it.

Adiabatic oscillation frequencies of stellar models, computed with the standard mixing-length formulation for convection, increasingly deviate with radial order from observations in solar-like stars. Standard solar models overestimate adiabatic frequencies by as much as ~ 20 {\mu}Hz. In this letter, we address the physical processes of turbulent convection that are predominantly responsible for the frequency differences between standard models and observations, also called `surface effects'. We compare measured solar frequencies from the MDI instrument on the SOHO spacecraft with frequency calculations that include three-dimensional (3D) hydrodynamical simulation results in the equilibrium model, nonadiabatic effects, and a consistent treatment of the turbulent pressure in both the equilibrium and stability computations. With the consistent inclusion of the above physics in our model computation we are able to reproduce the observed solar frequencies to < 3 {\mu}Hz without the need of any additional ad-hoc functional corrections.

HD 187547 was the first candidate that led to the suggestion that solar-like oscillations are present in delta Scuti stars. Longer observations, however, show that the modes interpreted as solar-like oscillations have either very long mode lifetimes, longer than 960 days, or are coherent. These results are incompatible with the nature of `pure' stochastic excitation as observed in solar-like stars. Nonetheless, one point is certain: the opacity mechanism alone cannot explain the oscillation spectrum of HD 187547. Here we present new theoretical investigations showing that convection dynamics can intrinsically excite coherent pulsations in the chemically peculiar delta Scuti star HD 187547. More precisely, it is the perturbations of the mean Reynold stresses (turbulent pressure) that drives the pulsations and the excitation takes place predominantly in the hydrogen ionization zone.

We perform a calibration of the mixing length of convection in stellar structure models against realistic 3D radiation-coupled hydrodynamics (RHD) simulations of convection in stellar surface layers, determining the adiabat deep in convective stellar envelopes. The mixing-length parameter $\alpha$ is calibrated by matching averages of the 3D simulations to 1D stellar envelope models, ensuring identical atomic physics in the two cases. This is done for a previously published grid of solar-metallicity convection simulations, covering from 4200 K to 6900 K on the main sequence, and 4300-5000 K for giants with logg=2.2. Our calibration results in an $\alpha$ varying from 1.6 for the warmest dwarf, which is just cool enough to admit a convective envelope, and up to 2.05 for the coolest dwarfs in our grid. In between these is a triangular plateau of $\alpha$ ~ 1.76. The Sun is located on this plateau and has seen little change during its evolution so far. When stars ascend the giant branch, they largely do so along tracks of constant $\alpha$, with $\alpha$ decreasing with increasing mass.

Recently the number of main-sequence and subgiant stars exhibiting solar-like oscillations that are resolved into individual mode frequencies has increased dramatically. While only a few such data sets were available for detailed modeling just a decade ago, the Kepler mission has produced suitable observations for hundreds of new targets. This rapid expansion in observational capacity has been accompanied by a shift in analysis and modeling strategies to yield uniform sets of derived stellar properties more quickly and easily. We use previously published asteroseismic and spectroscopic data sets to provide a uniform analysis of 42 solar-type Kepler targets from the Asteroseismic Modeling Portal (AMP). We find that fitting the individual frequencies typically doubles the precision of the asteroseismic radius, mass and age compared to grid-based modeling of the global oscillation properties, and improves the precision of the radius and mass by about a factor of three over empirical scaling relations. We demonstrate the utility of the derived properties with several applications.

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

Relations between temperature, T, and optical depth, tau, are often used for describing the photospheric transition from optically thick to optically thin in stellar structure models. We show that this is well justified, but also that currently used T(tau) relations are often inconsistent with their implementation. As an outer boundary condition on the system of stellar structure equations, T(tau) relations have an undue effect on the overall structure of stars. In this age of precision asteroseismology, we need to re-assess both the method for computing and for implementing T(tau) relations, and the assumptions they rest on. We develop a formulation for proper and consistent evaluation of T(tau) relations from arbitrary 1D or 3D stellar atmospheres, and for their implementation in stellar structure and evolution models. We extract radiative T(tau) relations, as described by our new formulation, from 3D simulations of convection in deep stellar atmospheres of late-type stars from dwarfs to giants. These simulations employ realistic opacities and equation of state, and account for line-blanketing. For comparison, we also extract T(tau) relations from 1D MARCS model atmospheres using the same formulation. T(tau)-relations from our grid of 3D convection simulations display a larger range of behaviours with surface gravity, compared with those of conventional theoretical 1D hydrostatic atmosphere models. Based on this, we recommend no longer to use scaled solar T(tau) relations. Files with T(tau) relations for our grid of simulations are made available to the community, together with routines for interpolating in this irregular grid. We also provide matching tables of atmospheric opacity, for consistent implementation in stellar structure models.

We study the effects of different descriptions of the solar surface convection on the eigenfrequencies of p-modes. 1-D evolution calculations of the whole Sun and 3-D hydrodynamic and magnetohydrodynamic simulations of the current surface are performed. These calculations rely on realistic physics. Averaged stratifications of the 3-D simulations are introduced in the 1-D solar evolution or in the structure models. The eigenfrequencies obtained are compared to those of 1-D models relying on the usual phenomenologies of convection and to observations of the MDI instrument aboard SoHO. We also investigate how the magnetic activity could change the eigenfrequencies and the solar radius, assuming that, 3 Mm below the surface, the upgoing plasma advects a 1.2 kG horizontal field. All models and observed eigenfrequencies are fairly close below 3 mHz. Above 3 mHz the eigenfrequencies of the phenomenological convection models are above the observed eigenfrequencies. The frequencies of the models based on the 3-D simulations are slightly below the observed frequencies. Their maximum deviation is ~ 3 \mu Hz at 3 mHz but drops below 1 \mu Hz at 4 mHz. Replacing the hydrodynamic by the magnetohydrodynamic simulation increases the eigenfrequencies. The shift is negligible below 2.2 mHz and then increases linearly with frequency to reach ~ 1.7 \mu Hz at 4 mHz. The impact of the simulated activity is a 14 milliarcsecond shrinking of the solar layers near the optical depth unity.

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

Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & M\"uller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & M\"uller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D models are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.

Present grids of stellar atmosphere models are the workhorses in interpreting stellar observations, and determining their fundamental parameters. These models rely on greatly simplified models of convection, however, lending less predictive power to such models of late type stars. We present a grid of improved and more reliable stellar atmosphere models of late type stars, based on deep, 3D, convective, stellar atmosphere simulations. This grid is to be used in general for interpreting observations, and improve stellar and asteroseismic modeling. We solve the Navier Stokes equations in 3D and concurrent with the radiative transfer equation, for a range of atmospheric parameters, covering most of stellar evolution with convection at the surface. We emphasize use of the best available atomic physics for quantitative predictions and comparisons with observations. We present granulation size, convective expansion of the acoustic cavity, asymptotic adiabat, as function of atmospheric parameters. These and other results are also available in electronic form.

It is well known that current spectroscopic determinations of the chemical composition of the Sun are starkly at odds with the metallicity implied by helioseismology. We investigate whether the discrepancy may be due to conversion of photons to a new light boson in the solar photosphere. We examine the impact of particles with axion-like interactions with the photon on the inferred photospheric abundances, showing that resonant axion-photon conversion is not possible in the region of the solar atmosphere in which line-formation occurs. Although non-resonant conversion in the line-forming regions can in principle impact derived abundances, constraints from axion-photon conversion experiments rule out the couplings necessary for these effects to be detectable. We show that this extends to hidden photons and chameleons (which would exhibit similar phenomenological behaviour), ruling out known theories of new light bosons as photospheric solutions to the solar abundance problem.

We obtain lower limits on the amplitude of convective velocities in the deep solar convection zone based only on the observed properties of the differential rotation and meridional circulation together with simple and robust dynamical balances obtained from the fundamental MHD equations. The linchpin of the approach is the concept of gyroscopic pumping whereby the meridional circulation across isosurfaces of specific angular momentum is linked to the angular momentum transport by the convective Reynolds stress. We find that the amplitude of the convective velocity must be at least 30 m s$^{-1}$ in the upper CZ ($r \sim 0.95 R$) and at least 8 m s$^{-1}$ in the lower CZ ($r \sim 0.75 R$) in order to be consistent with the observed mean flows. Using the base of the near-surface shear layer as a probe of the rotational influence, we are further able to show that the characteristic length scale of deep convective motions must be no smaller than 5.5--30 Mm. These results are compatible with convection models but suggest that the efficiency of the turbulent transport assumed in advection-dominated flux-transport dynamo models is generally not consistent with the mean flows they employ.

Stellar models generally use simple parametrizations to treat convection. The most widely used parametrization is the so-called "Mixing Length Theory" where the convective eddy sizes are described using a single number, \alpha, the mixing-length parameter. This is a free parameter, and the general practice is to calibrate \alpha using the known properties of the Sun and apply that to all stars. Using data from NASA's Kepler mission we show that using the solar-calibrated \alpha is not always appropriate, and that in many cases it would lead to estimates of initial helium abundances that are lower than the primordial helium abundance. Kepler data allow us to calibrate \alpha for many other stars and we show that for the sample of stars we have studied, the mixing-length parameter is generally lower than the solar value. We studied the correlation between \alpha and stellar properties, and we find that \alpha increases with metallicity. We therefore conclude that results obtained by fitting stellar models or by using population-synthesis models constructed with solar values of \alpha are likely to have large systematic errors. Our results also confirm theoretical expectations that the mixing-length parameter should vary with stellar properties.

Asteroseismology with the Kepler space telescope is providing not only an improved characterization of exoplanets and their host stars, but also a new window on stellar structure and evolution for the large sample of solar-type stars in the field. We perform a uniform analysis of 22 of the brightest asteroseismic targets with the highest signal-to-noise ratio observed for 1 month each during the first year of the mission, and we quantify the precision and relative accuracy of asteroseismic determinations of the stellar radius, mass, and age that are possible using various methods. We present the properties of each star in the sample derived from an automated analysis of the individual oscillation frequencies and other observational constraints using the Asteroseismic Modeling Portal (AMP), and we compare them to the results of model-grid-based methods that fit the global oscillation properties. We find that fitting the individual frequencies typically yields asteroseismic radii and masses to \sim1% precision, and ages to \sim2.5% precision (respectively 2, 5, and 8 times better than fitting the global oscillation properties). The absolute level of agreement between the results from different approaches is also encouraging, with model-grid-based methods yielding slightly smaller estimates of the radius and mass and slightly older values for the stellar age relative to AMP, which computes a large number of dedicated models for each star. The sample of targets for which this type of analysis is possible will grow as longer data sets are obtained during the remainder of the mission.

More than 1000 red giants have been observed by NASA/Kepler mission during a nearly continuous period of ~ 13 months. The resulting high-frequency resolution (< 0.03 muHz) allows us to study the granulation parameters of these stars. The granulation pattern results from the convection motions leading to upward flows of hot plasma and downward flows of cooler plasma. We fitted Harvey-like functions to the power spectra, to retrieve the timescale and amplitude of granulation. We show that there is an anti-correlation between both of these parameters and the position of maximum power of acoustic modes, while we also find a correlation with the radius, which agrees with the theory. We finally compare our results with 3D models of the convection.

The granulation pattern that we observe on the surface of the Sun is due to hot plasma from the interior rising to the photosphere where it cools down, and descends back into the interior at the edges of granules. This is the visible manifestation of convection taking place in the outer part of the solar convection zone. Because red giants have deeper convection zones and more extended atmospheres than the Sun, we cannot a priori assume that granulation in red giants is a scaled version of solar granulation. Until now, neither observations nor 1D analytical convection models could put constraints on granulation in red giants. However, thanks to asteroseismology, this study can now be performed. The resulting parameters yield physical information about the granulation. We analyze \sim1000 red giants that have been observed by Kepler during 13 months. We fit the power spectra with Harvey-like profiles to retrieve the characteristics of the granulation (time scale tau_gran and power P_gran). We also introduce a new time scale, tau_eff, which takes into account that different slopes are used in the Harvey functions. We search for a correlation between these parameters and the global acoustic-mode parameter (the position of maximum power, nu_max) as well as with stellar parameters (mass, radius, surface gravity (log g) and effective temperature (T_eff)). We show that tau_eff nu_max^{-0.89} and P_gran nu_max^{-1.90}, which is consistent with the theoretical predictions. We find that the granulation time scales of stars that belong to the red clump have similar values while the time scales of stars in the red-giant branch are spread in a wider range. Finally, we show that realistic 3D simulations of the surface convection in stars, spanning the (T_eff, log g)-range of our sample of red giants, match the Kepler observations well in terms of trends.

The scale length over which convection mixes mass in a star can be calculated as the inverse of the vertical derivative of the unidirectional (up or down) mass flux. This is related to the mixing length in the mixing length theory of stellar convection. We give the ratio of mass mixing length to pressure scale height for a grid of 3D surface convection simulations, covering from 4300\,K to 6900\,K on the main-sequence, and up to giants at $\log g = 2.2$, all for solar composition. These simulations also confirm what is already known from solar simulations, that convection doesn't proceed by discrete convective elements, but rather as a continuous, slow, smooth, warm upflow and turbulent, entropy deficient, fast down drafts. This convective topology also results in mixing on a scale as that of the classic mixing length formulation, and is simply a consequence of mass conservation on flows in a stratified atmosphere.

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

As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary conditions and minimization of sub-grid scale diffusive processes become increasingly crucial. We here assess changes in the dynamics and the energy flux balance of the flows established in rotating spherical shell segments that capture much of the NSSL with the Curved Spherical Segment (CSS) code using two different diffusion schemes. The CSS code is a new massively parallel modeling tool capable of simulating 3-D compressible MHD convection with a realistic solar stratification in rotating spherical shell segments.

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

We 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.

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

We present the results of realistic, 3D, hydrodynamical, simulations of surface convection in red giant stars with varying effective temperatures and metallicities. We use the convection simulations as time-dependent, hydrodynamical, model atmospheres to compute spectral line profiles for a number of ions and molecules under the assumption of local thermodynamic equilibrium (LTE). We compare the results with the predictions of line formation calculations based on 1D, hydrostatic, model stellar atmospheres in order to estimate the impact of 3D models on the derivation of elemental abundances. We find large negative 3D-1D LTE abundance corrections (typically -0.5 to -1 dex) for weak low-excitation lines from molecules and neutral species in the very low metallicity cases. Finally, we discuss the extent of departures from LTE in the case of neutral iron spectral line formation.

We investigate the impact of 3D hydrodynamical model atmospheres of red giant stars at different metallicities on the formation of spectral lines of a number of ions and molecules. We carry out realistic 3D simulations of surface convection in red giant stars with varying stellar parameters. We use the simulations as time-dependent hydrodynamical model stellar atmospheres to compute atomic (Li, O, Na, Mg, Ca, Fe) and molecular (CH, NH, OH) spectral lines under the assumption of local thermodynamic equilibrium (LTE). We compare the line strengths computed in 3D with the results of analogous line formation calculations for 1D, hydrostatic, plane-parallel MARCS model atmospheres in order to estimate the impact of 3D models on the derivation of elemental abundances. The temperature and density inhomogeneities and correlated velocities in 3D models, as well as the differences between the 1D and mean 3D structures significantly affect the predicted line strengths. Under the assumption of LTE, the low atmospheric temperatures of very metal-poor 3D model atmospheres cause the lines from neutral species and molecules to appear stronger than in 1D. Therefore, elemental abundances derived from these lines using 3D models are significantly lower than according to 1D analyses. Differences between 3D and 1D abundances of C, N, and O derived from CH, NH, and OH weak low-excitation lines are found to be in the range -0.5 dex to -1.0 dex for the the red giant stars at [Fe/H]=-3 considered here. At this metallicity, large negative corrections (about -0.8 dex) are also found for weak low-excitation Fe I lines. We caution, however, that departures from LTE might be significant for these and other elements and comparable to the effects due to stellar granulation.

We investigate the impact of realistic 3D hydrodynamical model stellar atmospheres on the determination of elemental abundances in the carbon-rich, hyper iron-poor stars HE0107-5240 and HE1327-2326. We derive the chemical compositions of the two stars by means of a detailed 3D analysis of spectral lines under the assumption of local thermodynamic equilibrium (LTE). The lower temperatures of the line-forming regions of the hydrodynamical models cause changes in the predicted spectral line strengths. In particular we find the 3D abundances of C, N, and O to be lower by ~ -0.8 dex (or more) than estimated from a 1D analysis. The 3D abundances of iron-peak elements are also decreased but by smaller factors (~ -0.2 dex). We caution however that the neglected non-LTE effects might actually be substantial for these metals. We finally discuss possible implications for studies of early Galactic chemical evolution.

A detailed comparison is carried out between two popular equations of state (EOS), the Mihalas-Hummer-Dappen (MHD) and the OPAL equations of state, which have found widespread use in solar and stellar modeling during the past two decades. They are parts of two independent efforts to recalculate stellar opacities; the international Opacity Project (OP) and the Livermore-based OPAL project. We examine the difference between the two equations of state in a broad sense, over the whole applicable rho-T range, and for three different chemical mixtures. Such a global comparison highlights both their differences and their similarities. We find that omitting a questionable hard-sphere correction, tau, to the Coulomb interaction in the MHD formulation, greatly improves the agreement between the MHD and OPAL EOS. We also find signs of differences that could stem from quantum effects not yet included in the MHD EOS, and differences in the ionization zones that are probably caused by differences in the mechanisms for pressure ionization. Our analysis do not only give a clearer perception of the limitations of each equation of state for astrophysical applications, but also serve as guidance for future work on the physical issues behind the differences. The outcome should be an improvement of both equations of state.

I present an overview of an equation of state, being developed in the chemical picture, and based on the very successful MHD equation of state. The flexibility of the chemical picture combined with the free-energy minimization procedure, makes it rather straight-forward, albeit laborious, to include new effects in the model free-energy, simply by adding new terms. The most notable additions to the original MHD equation of state, are relativistic effects, quantum effects, improved higher order Coulomb terms and a long list of molecules other than the H2 and H2+ treated so far.

The convection code of Nordlund & Stein has been used to evaluate the 3D, radiation-coupled convection in a stellar atmosphere with Teff=7300K, logg=4.3 and [Fe/H]=0.0, corresponding to a main-sequence A9-star. I will present preliminary comparisons between the 3D-simulation and a conventional 1D stellar structure calculation, and elaborate on the consequences of the differences.

We compute the rates P at which energy is injected into the p modes for a set of 3D simulations of outer layers of stars. We found that Pmax - the maximum in P - scales as (L/M)^s where s is the slope of the power law, L and M are the luminosity and the mass of the 1D stellar models associated with the simulations. The slope is found to depend significantly on the adopted representation for the turbulent eddy-time correlation function, chi_k. According to the expected performances of COROT, it will likely be possible to measure Pmax as a function of L/M and to constrain the properties of stellar turbulence as the turbulent eddy time-correlation.

Hydrodynamical, 3D simulations of the outer layers of the Sun and Alpha Cen A are used to obtain constraints on the properties of turbulent convection in such stars. These constraints enable us to compute - on the base of a theoretical model of stochastic excitation - the rate P at which p modes are excited by turbulent convection in those two stars. Results are then compared with solar seismic observations and recent observations of Alpha Cen A. For the Sun, a good agreement between observations and computed P is obtained. For Alpha Cen A a large discrepancy is obtained which origin cannot be yet identified: it can either be caused by the present data quality which is not sufficient for our purpose or by the way the intrinsic amplitudes and the life-times of the modes are determined or finally attributed to our present modelling. Nevertheless, data with higher quality or/and more adapted data reductions will likely provide constraints on the p-mode excitation mechanism in Alpha Cen A.

The solar photospheric Fe abundance has been determined using realistic ab initio 3D, time-dependent, hydrodynamical model atmospheres. The study is based on the excellent agreement between the predicted and observed line profiles directly rather than equivalent width, since the intrinsic Doppler broadening from the convective motions and oscillations provide the necessary non-thermal broadening. Thus, three of the four hotly debated parameters (equivalent widths, microturbulence and damping enhancement factors) in the center of the recent solar Fe abundance dispute regarding FeI lines no longer enter the analysis, leaving the transition probabilities as the main uncertainty. Both FeI (using the samples of lines of both the Oxford and Kiel studies) and FeII lines have been investigated, which give consistent results: log FeI = 7.44 +- 0.05 and log FeII = 7.45 +- 0.10. Also the wings of strong FeI lines return consistent abundances, log FeII = 7.42 +- 0.03, but due to the uncertainties inherent in analyses of strong lines we give this determination lower weight than the results from weak and intermediate strong lines. In view of the recent slight downward revision of the meteoritic Fe abundance log Fe = 7.46 +- 0.01, the agreement between the meteoritic and photospheric values is very good, thus appearingly settling the debate over the photospheric Fe abundance from FeI lines.

Realistic ab-initio 3D, radiative-hydrodynamical convection simulations of the solar granulation have been applied to FeI and FeII line formation. In contrast to classical analyses based on hydrostatic 1D model atmospheres the procedure contains no adjustable free parameters but the treatment of the numerical viscosity in the construction of the 3D, time-dependent, inhomogeneous model atmosphere and the elemental abundance in the 3D spectral synthesis. However, the numerical viscosity is introduced purely for numerical stability purposes and is determined from standard hydrodynamical test cases with no adjustments allowed to improve the agreement with the observational constraints from the solar granulation. The non-thermal line broadening is mainly provided by the Doppler shifts arising from the convective flows in the solar photosphere and the solar oscillations. The almost perfect agreement between the predicted temporally and spatially averaged line profiles for weak Fe lines with the observed profiles and the absence of trends in derived abundances with line strengths, seem to imply that the micro- and macroturbulence concepts are obsolete in these 3D analyses. Furthermore, the theoretical line asymmetries and shifts show a very satisfactory agreement with observations with an accuracy of typically 50-100 m/s on an absolute velocity scale. The remaining minor discrepancies point to how the convection simulations can be refined further.

Realistic 3-dimensional (3D), radiative hydrodynamical surface convection simulations of the metal-poor halo stars HD 140283 and HD 84937 have been performed. Due to the dominance of adiabatic cooling over radiative heating very low atmospheric temperatures are encountered. The lack of spectral lines in these metal-poor stars thus causes much steeper temperature gradients than in classical 1D hydrostatic model atmospheres where the temperature of the optically thin layers is determined by radiative equilibrium. The modified atmospheric structures cause changes in the emergent stellar spectra. In particular, the primordial Li abundances may have been overestimated by 0.2-0.35 dex with 1D model atmospheres. However, we caution that our result assumes local thermodynamic equilibrium (LTE), while the steep temperature gradients may be prone to e.g. over-ionization.

Differences between observed and theoretical eigenfrequencies of the Sun have characteristics which identify them as arising predominantly from properties of the oscillations in the vicinity of the solar surface: in the super-adiabatic, convective boundary layer and above. These frequency differences may therefore provide useful information about the structure of these regions, precisely where the theory of solar structure is most uncertain. In the present work we use numerical simulations of the outer part of the Sun to quantify the influence of turbulent convection on solar oscillation frequencies. Separating the influence into effects on the mean model and effects on the physics of the modes, we find that the main model effects are due to the turbulent pressure that provides additional support against gravity, and thermal differences between average 3-D models and 1-D models. Surfaces of constant pressure in the visible photosphere are elevated by about 150 km, relative to a standard envelope model. As a result, the turning points of high-frequency modes are raised, while those of the low-frequency modes remain essentially unaffected. The corresponding gradual lowering of the mode frequencies accounts for most of the frequency difference between observations and standard solar models. Additional effects are expected to come primarily from changes in the physics of the modes, in particular from the modulation of the turbulent pressure by the oscillations.

Advances in hydrodynamical simulations have provided new insights into the effects of convection on the frequencies of solar oscillations. As more accurate observations become available, this may lead to an improved understanding of the dynamics of convection and the interaction between convection and pulsation (\cite{Rosenthal+99}). Recent high-resolution observations from the SOI/MDI instrument on the SOHO spacecraft have provided the so-far most-detailed observations of high-degree modes of solar oscillations, which are particularly sensitive to the near-surface properties of the Sun. Here we present preliminary results of a comparison between these observations and frequencies computed for models based on realistic simulations of near-surface convection. Such comparisons may be expected to help in identifying the causes of the remaining differences between the observed frequencies and those of solar models.

Following the report of solar-like oscillations in the G0 V star eta Boo (Kjeldsen et al. 1995, AJ 109, 1313), a first attempt to model the observed frequencies was made by Christensen-Dalsgaard et al. (1995, ApJ Letters, in press). This attempt succeeded in reproducing the observed frequency separations, although there remained a difference of about 10 microHz between observed and computed frequencies. In those models, the near-surface region of the star was treated rather crudely. Here we consider more sophisticated models that include non-local mixing-length theory, turbulent pressure and nonadiabatic oscillations.