Search

Balmer lines serve as important indicators of stellar effective temperatures in late-type stellar spectra. One of their modelling uncertainties is the influence of convective flows on their shape. We aim to characterize the influence of convection on the wings of Balmer lines. We perform a differential comparison of synthetic Balmer line profiles obtained from 3D hydrodynamical model atmospheres and 1D hydrostatic standard ones. The model parameters are appropriate for F,G,K dwarf and subgiant stars of metallicity ranging from solar to 1/1000 solar. The shape of the Balmer lines predicted by 3D models can never be exactly reproduced by a 1D model, irrespective of its effective temperature. We introduce the concept of a 3D temperature correction, as the effective temperature difference between a 3D model and a 1D model which provides the closest match to the 3D profile. The temperature correction is different for the different members of the Balmer series and depends on the adopted mixing-length parameter in the 1D model. Among the investigated models, the 3D correction ranges from -300K to +300K. Horizontal temperature fluctuations tend to reduce the 3D correction. Accurate effective temperatures cannot be derived from the wings of Balmer lines, unless the effects of convection are properly accounted for. The 3D models offer a physically well justified way of doing so. The use of 1D models treating convection with the mixing-length theory do not appear to be suitable for this purpose. In particular, there are indications that it is not possible to determine a single value of the mixing-length parameter which will optimally reproduce the Balmer lines for any choice of atmospheric parameters., Comment: 6 pages, 3 figures, accepted for publication in A&A

Solar-like oscillations have been observed in numerous red giants from ground and from space. An important question arises: could we expect to detect non-radial modes probing the internal structure of these stars? We investigate under what physical circumstances non-radial modes could be observable in red giants; what would be their amplitudes, lifetimes and heights in the power spectrum (PS)? Using a non-radial non-adiabatic pulsation code including a non-local time-dependent treatment of convection, we compute the theoretical lifetimes of radial and non-radial modes in several red giant models. Next, using a stochastic excitation model, we compute the amplitudes of these modes and their heights in the PS. Very distinct results are found depending on the evolutionary status of the star., Comment: Source tex file : 1713.tex Bibliography file : 1713.bbl 24 figures : 1713f1.eps, ..., 1713f24.eps

Local-box hydrodynamical model atmospheres provide statistical information about a star's emergent radiation field which allows one to predict the level of its granulation-related micro-variability. Space-based photometry is now sufficiently accurate to test model predictions. We aim to model the photometric granulation background of HD49933 as well as the Sun, and compare the predictions to the measurements obtained by the COROT and SOHO satellite missions. We construct hydrodynamical model atmospheres representing HD49933 and the Sun, and use a previously developed scaling technique to obtain the observable disk-integrated brightness fluctuations. We further performed exploratory magneto-hydrodynamical simulations to gauge the impact of small scale magnetic fields on the synthetic light-curves. We find that the granulation-related brightness fluctuations depend on metallicity. We obtain a satisfactory correspondence between prediction and observation for the Sun, validating our approach. For HD49933, we arrive at a significant over-estimation by a factor of two to three in total power. Locally generated magnetic fields are unlikely to be responsible, otherwise existing fields would need to be rather strong to sufficiently suppress the granulation signal. Presently suggested updates on the fundamental stellar parameters do not improve the correspondence; however, an ad-hoc increase of the HD49933 surface gravity by about 0.2dex would eliminate most of the discrepancy. We diagnose a puzzling discrepancy between the predicted and observed granulation background in HD49933, with only rather ad-hoc ideas for remedies at hand., Comment: 7 pages, 5 figures, accepted for publication in A&A

The light elements Li, O, Na, Al, and Mg are known to show star to star variations in the globular clusters 47 Tuc and NGC 6752. Such variations are interpreted as due to processing in a previous generation of stars. In this paper we investigate the abundances of the alpha-element sulphur, for which no previous measurements exist. In fact this element has not been investigated in any Galactic globular cluster so far. The only globular cluster for which such measurements are available is Terzan 7, which belongs to the Sgr dSph. We use high resolution spectra of the S I Mult. 1, acquired with the UVES spectrograph at the 8.2m VLT-Kueyen telescope, for turn-off and giant stars in the two globular clusters. The spectra are analyzed making use of ATLAS static plane parallel model atmospheres and SYNTHE spectrum synthesis. We also compute 3D corrections from CO5BOLD hydrodynamic models and apply corrections due to NLTE effects taken from the literature. In the cluster NGC 6752 sulphur has been measured only in four subgiant stars. We find no significant star to star scatter and a mean <[S/Fe]> = +0.49 +- 0.15, consistent with what observed in field stars of the same metallicity. In the cluster 47 Tuc we measured S in 4 turn-off and 5 subgiant stars with a mean <[S/Fe]> = +0.18 +- 0.14. While this result is compatible with no star to star scatter we notice a statistically significant correlation of the sulphur abundance with the sodium abundance and a tentative correlation with the silicon abundance., Comment: 8 pages, 4 figures, accepted for publication in A&A

CONTEXT:The detailed chemical abundances of extremely metal-poor (EMP) stars are key guides to understanding the early chemical evolution of the Galaxy. Most existing data are, however, for giant stars which may have experienced internal mixing later. AIMS: We aim to compare the results for giants with new, accurate abundances for all observable elements in 18 EMP turnoff stars. METHODS:VLT/UVES spectra at R ~45,000 and S/N~ 130 per pixel (330-1000 nm) are analysed with OSMARCS model atmospheres and the TURBOSPECTRUM code to derive abundances for C, Mg, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, Zn, Sr, and Ba. RESULTS: For Ca, Ni, Sr, and Ba, we find excellent consistency with our earlier sample of EMP giants, at all metallicities. However, our abundances of C, Sc, Ti, Cr, Mn and Co are ~0.2 dex larger than in giants of similar metallicity. Mg and Si abundances are ~0.2 dex lower (the giant [Mg/Fe] values are slightly revised), while Zn is again ~0.4 dex higher than in giants of similar [Fe/H] (6 stars only). CONCLUSIONS:For C, the dwarf/giant discrepancy could possibly have an astrophysical cause, but for the other elements it must arise from shortcomings in the analysis. Approximate computations of granulation (3D) effects yield smaller corrections for giants than for dwarfs, but suggest that this is an unlikely explanation, except perhaps for C, Cr, and Mn. NLTE computations for Na and Al provide consistent abundances between dwarfs and giants, unlike the LTE results, and would be highly desirable for the other discrepant elements as well. Meanwhile, we recommend using the giant abundances as reference data for Galactic chemical evolution models., Comment: A&A, Received 16 July 2008/ Accepted 15 March 2009

CONTEXT: In recent years, the solar chemical abundances have been studied in considerable detail because of discrepant values of solar metallicity inferred from different indicators, i.e., on the one hand, the "sub-solar" photospheric abundances resulting from spectroscopic chemical composition analyses with the aid of 3D hydrodynamical models of the solar atmosphere, and, on the other hand, the high metallicity inferred by helioseismology. AIMS: After investigating the solar oxygen abundance using a CO5BOLD 3D hydrodynamical solar model in previous work, we undertake a similar approach studying the solar abundance of nitrogen, since this element accounts for a significant fraction of the overall solar metallicity, Z. METHOD: We used a selection of atomic spectral lines to determine the solar nitrogen abundance, relying mainly on equivalent width measurements in the literature. We investigate the influence on the abundance analysis, of both deviations from local thermodynamic equilibrium ("NLTE effects") and photospheric inhomogeneities ("granulation effects"). RESULTS: We recommend use of a solar nitrogen abundance of A(N)=7.86+-0.12 whose error bar reflects the line-to-line scatter. CONCLUSION: The solar metallicity implied by the CO5BOLD-based nitrogen and oxygen abundances is in the range 0.0145<= Z <= 0.0167. This result is a step towards reconciling photospheric abundances with helioseismic constraints on Z. Our most suitable estimates are Z=0.0156 and Z/X=0.0213., Comment: To be published on A&A

We present 2D local box simulations of near-surface radiative magneto-convection with prescribed magnetic flux, carried out with the MHD version of the CO5BOLD code for the Sun and a solar-like star with a metal-poor chemical composition (metal abundances reduced by a factor 100, [M/H]=-2). The resulting magneto-hydrodynamical models can be used to study the influence of the metallicity on the properties of magnetized stellar atmospheres. A preliminary analysis indicates that the horizontal magnetic field component tends to be significantly stronger in the optically thin layers of metal-poor stellar atmospheres., Comment: Proc. IAU Symposium 259, Cosmic Magnetic Fields: from Planets, to Stars and Galaxies, K.G. Strassmeier, A.G. Kosovichev and J.E. Beckman, eds. (2009) p.233

We present synthetic FeH band spectra in the z-filter range for several M-dwarf models with logg=3.0-5.0 [cgs] and Teff=2800K -3450K. Our aim is to characterize convective velocities in M-dwarfs and to give a rough estimate of the range in which 3D-atmosphere treatment is necessary and where 1D-atmosphere models suffice for the interpretation of molecular spectral features. This is also important in order to distinguish between the velocity-broadening and the rotational- or Zeeman-broadening. The synthetic spectra were calculated using 3D CO5BOLD radiative-hydrodynamic (RHD) models and the line synthesis code LINFOR3D. We used complete 3D-models and high resolution 3D spectral synthesis for the detailed study of some well isolated FeH lines. The FeH line strength shows a dependence on surface gravity and effective temperature and could be employed to measure both quantities in M-type objects. The line width is related to the velocity-field in the model stars, which depends strongly on surface gravity. Furthermore, we investigate the velocity-field in the 3D M-dwarf models together with the related micro- and macro-turbulent velocities in the 1D case. We also search for effects on the lineshapes., Comment: Cool Stars 15 Conference Proceeding, 4 pages

We give a progress report about the activities within the CIFIST Team related to the search for extremely metal-poor stars in the Sloan Digital Sky Survey's spectroscopic catalog. So far the search has provided 25 candidates with metallicities around or smaller -3. For 15 candidates high resolution spectroscopy with UVES at the VLT has confirmed their extremely metal-poor status. Work is under way to extend the search to the SDSS's photometric catalog by augmenting the SDSS photometry, and by gauging the capabilities of X-shooter when going to significantly fainter targets., Comment: 6 pages, 6 figures, Proceedings paper of the conference "A stellar journey: A symposium in celebration of Bengt Gustafsson's 65th birthday"

Context: The stable element hafnium (Hf) and the radioactive element thorium (Th) were recently suggested as a suitable pair for radioactive dating of stars. The applicability of this elemental pair needs to be established for stellar spectroscopy. Aims: We aim at a spectroscopic determination of the abundance of Hf and Th in the solar photosphere based on a \cobold 3D hydrodynamical model atmosphere. We put this into a wider context by investigating 3D abundance corrections for a set of G- and F-type dwarfs. Method: High-resolution, high signal-to-noise solar spectra were compared to line synthesis calculations performed on a solar CO5BOLD model. For the other atmospheres, we compared synthetic spectra of CO5BOLD 3D and associated 1D models. Results: For Hf we find a photospheric abundance A(Hf)=0.87+-0.04, in good agreement with a previous analysis, based on 1D model atmospheres. The weak Th ii 401.9 nm line constitutes the only Th abundance indicator available in the solar spectrum. It lies in the red wing of an Ni-Fe blend exhibiting a non-negligible convective asymmetry. Accounting for the asymmetry-related additional absorption, we obtain A(Th)=0.09+-0.03, consistent with the meteoritic abundance, and about 0.1 dex lower than obtained in previous photospheric abundance determinations. Conclusions: Only for the second time, to our knowledge, has am non-negligible effect of convective line asymmetries on an abundance derivation been highlighted. Three-dimensional hydrodynamical simulations should be employed to measure Th abundances in dwarfs if similar blending is present, as in the solar case. In contrast, 3D effects on Hf abundances are small in G- to mid F-type dwarfs and sub-giants, and 1D model atmospheres can be conveniently used., Comment: A&A, in press

Context. Europium is an almost pure r-process element, which may be useful as a reference in nucleocosmochronology. Aims. To determine the photospheric solar abundance using CO5BOLD 3-D hydrodynamical model atmospheres. Methods. Disc-centre and integrated-flux observed solar spectra are used. The europium abundance is derived from the equivalent width measurements. As a reference 1D model atmospheres have been used, in addition. Results. The europium photospheric solar abundance is 0.52 +- 0.02 in agreement with previous determinations. We also determine the photospheric isotopic fraction of Eu(151) to be 49 % +- 2.3 % from the intensity spectra and 50% +-2.3 from the flux spectra. This compares well to the the meteoritic isotopic fraction 47.8%. We explore the 3D corrections also for dwarfs and sub-giants in the temperature range ~5000 K to ~6500 K and solar and 1/10--solar metallicities and find them to be negligible for all the models investigated. Conclusions. Our photospheric Eu abundance is in good agreement with previous determinations based on 1D models. This is in line with our conclusion that 3D effects for this element are negligible in the case of the Sun., Comment: Accepted for publication on A&A

In ESO period 65 (April-September 2000) the large programme 165.N-0276, led by Roger Cayrel, began making use of UVES at the Kueyen VLT telescope. Known within the Team and outside as ``First Stars'', it was aimed at obtaining high resolution, high signal-to-noise ratio spectra in the range 320 nm -- 1000 nm for a large sample of extremely metal-poor (EMP) stars identified from the HK objective prism survey.In this contribution we highlight the main results of the large programme., Comment: to be published in the proceedings of the Workshop "Science with VLT in the ELT era" 8-12 October 2007, Garching, ed. A. Moorwood

We have used high-resolution, high-S/N ratio spectra from the UVES spectrograph at the ESO VLT telescope. Long-term radial-velocity measurements and broad-band photometry allow us to determine improved orbital elements and stellar parameters for both components. We use OSMARCS 1D models and the {{\scshape turbospectrum}} spectral synthesis code to determine the abundances of Li, O, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Fe, Co and Ni. We also use the CO$^5$BOLD model atmosphere code to compute the 3D abundance corrections, especially for Li and O. We find a metallicity of [Fe/H]$\sim -3.6$ for both stars, using 1D models with 3D corrections of $\sim -0.1$ dex from averaged 3D models. We determine the oxygen abundance from the near-UV OH bands; the 3D corrections are large, -1 and -1.5 dex for the secondary and primary respectively, and yield [O/Fe] $\sim 0.8$, close to the high-quality results obtained from the [OI] 630 nm line in metal-poor giants. Other [$\alpha$/Fe] ratios are consistent with those measured in other dwarfs and giants with similar [Fe/H], although Ca and Si are somewhat low ([X/Fe]$\la 0$). Other element ratios follow those of other halo stars. The Li abundance of the primary star is consistent with the Spite plateau, but the secondary shows a lower abundance; 3D corrections are small. The Li abundance in the primary star supports the extension of the {{\em Spite Plateau}} value at the lowest metallicities, without any decrease. The low abundance in the secondary star could be explained by endogenic Li depletion, due to its cooler temperature. If this is not the case, another, yet unknown mechanism may be causing increased scatter in A(Li) at the lowest metallicities., Comment: Accepted for publication in Astronomy & Astrophysics

The continuum intensity at wavelengths around 1 mm provides an excellent way to probe the solar chromosphere. Future high-resolution millimetre arrays, such as the Atacama Large Millimeter Array (ALMA), will thus produce valuable input for the ongoing controversy on the thermal structure and the dynamics of this layer. Synthetic brightness temperature maps are calculated on basis of three-dimensional radiation (magneto-)hydrodynamic (MHD) simulations. While the millimetre continuum at 0.3mm originates mainly from the upper photosphere, the longer wavelengths considered here map the low and middle chromosphere. The effective formation height increases generally with wavelength and also from disk-centre towards the solar limb. The average intensity contribution functions are usually rather broad and in some cases they are even double-peaked as there are contributions from hot shock waves and cool post-shock regions in the model chromosphere. Taking into account the deviations from ionisation equilibrium for hydrogen gives a less strong variation of the electron density and with it of the optical depth. The result is a narrower formation height range. The average brightness temperature increases with wavelength and towards the limb. The relative contrast depends on wavelength in the same way as the average intensity but decreases towards the limb. The dependence of the brightness temperature distribution on wavelength and disk-position can be explained with the differences in formation height and the variation of temperature fluctuations with height in the model atmospheres., Comment: 15 pages, 10 figures, accepted for publication in A&A (15.05.07)

We constructed hydrodynamical model atmospheres for mid M-type main-, as well as pre-main-sequence (PMS) objects. Despite the complex chemistry encountered in these cool atmospheres a reasonably accurate representation of the radiative transfer is possible, even in the context of time-dependent and three-dimensional models. The models provide detailed information about the morphology of M-type granulation and statistical properties of the convective surface flows. In particular, we determined the efficiency of the convective energy transport, and the efficiency of mixing by convective overshoot. The convective transport efficiency was expressed in terms of an equivalent mixing-length parameter alpha in the formulation of mixing-length theory (MLT) given by Mihalas (1978). Alpha amounts to values around 2 for matching the entropy of the deep, adiabatically stratified regions of the convective envelope, and lies between 2.5 and 3.0 for matching the thermal structure of the deep photosphere. For current spectral analysis of PMS objects this implies that MLT models based on alpha=2.0 overestimate the effective temperature by 100 K and surface gravities by 0.25 dex. The average thermal structure of the formally convectively stable layers is little affected by convective overshoot and wave heating, i.e., stays close to radiative equilibrium conditions. Our models suggest that the rate of mixing by convective overshoot declines exponentially with geometrical distance to the Schwarzschild stability boundary. It increases at given effective temperature with decreasing gravitational acceleration., Comment: 17 pages, 23 figures, accepted for publication in A&A

We show that contrary to what is expected from 1D stationary model atmospheres, 3D hydrodynamical modeling predicts a considerable influence of convection on the spectral properties of late-type giants. This is due to the fact that convection overshoots into the formally stable outer atmospheric layers producing a notable granulation pattern in the 3D hydrodynamical models, which has a direct influence on the observable spectra and colors. Within the framework of standard 1D model atmospheres the average thermal stratification of the 3D hydro model can not be reproduced with any reasonable choice of the mixing length parameter and formulation of the turbulent pressure. The differences in individual photometric colors -- in terms of 3D versus 1D -- reach up to ~0.2 mag, or \Delta Teff~70K. We discuss the impact of full 3D hydrodynamical models on the interpretation of observable properties of late-type giants, briefly mentioning problems and challenges which need to be solved for bringing these models to a routine use within the astronomical community in 5-10 years from now., Comment: 4 pages, 3 figures. Proceedings of the IAU Symposium 232 "The Scientific Requirements for Extremely Large Telescopes", eds. P. Whitelock, B. Leibundgut, and M. Dennefeld

To assess the current status in the theoretical modeling of the spectral properties of late-type giants, we provide a comparison of synthetic photometric colors of late-type giants (calculated with PHOENIX, MARCS and ATLAS model atmospheres) with observations, at [M/H]=0.0 and -2.0. Overall, there is a good agreement between observed and synthetic colors, and synthetic colors and published Teff-color relations, both at [M/H]=0.0 and -2.0. Deviations from the observed trends in Teff-color planes are generally within \pm 150K (or less) in the effective temperature range Teff=3500-4800K. Synthetic colors calculated with different stellar atmosphere models typically agree to ~100K, within a large range of effective temperatures and gravities. Some discrepancies are seen in the Teff-(B-V) plane below Teff~3800K at [M/H]=0.0, due to difficulties in reproducing the 'turn-off' to the bluer colors which is seen in the observed data at Teff~3600K. Note that at [M/H]=-2.0 effective temperatures given by the scale of Alonso et al. (1999) are generally lower than those resulting from other Teff-color relations based both on observed and synthetic colors., Comment: 2 pages, 1 figure. Proceedings of the IAU Symposium 232 "The Scientific Requirements for Extremely Large Telescopes", eds. P. Whitelock, B. Leibundgut, and M. Dennefeld

We present new synthetic broad-band photometric colors for late-type giants based on synthetic spectra calculated with the PHOENIX model atmosphere code. The grid covers effective temperatures T_eff=3000...5000K, gravities log g=-0.5...+3.5, and metallicities [M/H]=+0.5...-4.0. We show that individual broad-band photometric colors are strongly affected by model parameters such as molecular opacities, gravity, microturbulent velocity, and stellar mass. Our exploratory 3D modeling of a prototypical late-type giant shows that convection has a noticeable effect on the photometric colors too, as it alters significantly both the vertical and horizontal thermal structures in the outer atmosphere. The differences between colors calculated with full 3D hydrodynamical and 1D model atmospheres are significant (e.g., \Delta(V-K)~0.2 mag), translating into offsets in effective temperature of up to ~70K. For a sample of 74 late-type giants in the Solar neighborhood, with interferometric effective temperatures and broad-band photometry available in the literature, we compare observed colors with a new PHOENIX grid of synthetic photometric colors, as well as with photometric colors calculated with the MARCS and ATLAS model atmosphere codes. (abridged), Comment: 30 pages, 21 figures, A&A in press. Table 2 can be obtained from the CDS or directly from the authors

Recent three-dimensional radiation hydrodynamic simulations by Wedemeyer et al. (2004) suggest that the solar chromosphere is highly structured in space and time on scales of only 1000 km and 20-25 sec, resp.. The resulting pattern consists of a network of hot gas and enclosed cool regions which are due to the propagation and interaction of shock fronts. In contrast to many other diagnostics, the radio continuum at millimeter wavelengths is formed in LTE, and provides a rather direct measure of the thermal structure. It thus facilitates the comparison between numerical model and observation. While the involved time and length scales are not accessible with todays equipment for that wavelength range, the next generation of instruments, such as the Atacama Large Millimeter Array (ALMA), will provide a big step towards the required resolution. Here we present results of radiative transfer calculations at mm and sub-mm wavelengths with emphasis on spatial and temporal resolution which are crucial for the ongoing discussion about the chromospheric temperature structure., Comment: 4 pages, 3 figures (1 colour)

We present first realistic numerical simulations of 3D radiative convection in the surface layers of main sequence A-type stars with Teff = 8000 K and 8500 K, log g = 4.4 and 4.0, recently performed with the CO5BOLD radiation hydrodynamics code. The resulting models are used to investigate the structure of the H+HeI and the HeII convection zones in comparison with the predictions of local and non-local convection theories, and to determine the amount of "overshoot" into the stable layers below the HeII convection zone. The simulations also predict how the topology of the photospheric granulation pattern changes from solar to A-type star convection. The influence of the photospheric temperature fluctuations and velocity fields on the shape of spectral lines is demonstrated by computing synthetic line profiles and line bisectors for some representative examples, allowing us to confront the 3D model results with observations., Comment: 5 pages, 6 figures (17 figure files), 1 Table

We have fit synthetic visibilities from 3-D (CO5BOLD + PHOENIX) and 1-D (PHOENIX, ATLAS 12) model stellar atmospheres of Procyon (F5 IV) to high-precision interferometric data from the VLTI Interferometer (K-band) and from the Mark III interferometer (500 nm and 800 nm). These data sets provide a test of theoretical wavelength dependent limb-darkening predictions. The work of Allende Prieto et al. has shown that the temperature structure from a spatially and temporally averaged 3-D hydrodynamical model produces significantly less limb darkening at 500 nm relative to the temperature structure of a 1-D MARCS model atmosphere with a standard mixing-length approximation for convection. Our direct fits to the interferometric data confirm this prediction. A 1-D ATLAS 12 model with ``approximate overshooting'' provides the required temperature gradient. We show, however, that 1-D models cannot reproduce the ultraviolet spectrophotometry below 160 nm with effective temperatures in the range constrained by the measured bolometric flux and angular diameter. We find that a good match to the full spectral energy distribution can be obtained with a composite model consisting of a weighted average of twelve 1-D model atmospheres based on the surface intensity distribution of a 3-D granulation simulation. We emphasize that 1-D models with overshooting may realistically represent the mean temperature structure of F-type stars like Procyon, but the same models will predict redder colors than observed because they lack the multicomponent temperature distribution expected for the surfaces of these stars., Comment: 24 pages, 8 figures, accepted for publication in the Astrophysical Journal

Late type giants (RGB/AGB stars) will be important tracers of the Galactic morphology and evolution in the framework of Gaia, as they are intrinsically bright and thus can probe distant stellar populations or those obscured by interstellar extinction. A realistic representation of their atmospheres and spectra with stellar atmosphere models is thus of crucial importance, both for the design and optimization of Gaia instruments, as well as the interpretation of provided astrophysical data. Our analysis of synthetic photometric colors of late-type giants based on PHOENIX, MARCS and ATLAS model atmospheres indicates a general agreement between the current theoretical predictions and observations in the framework of stationary 1-D model atmospheres. Presently available models allow temperature determinations of RGB/AGB stars to an accuracy of ~\pm100 K. In an exploratory study we try to quantify possible residual systematic effects due to the approximations made in 1-D models using full 3-D hydrodynamical models. We find that differences in broad-band photometric colors calculated with 1-D and 3-D models are significant, translating to the offsets in effective temperature of up to ~70 K. Clearly, full 3-D hydrodynamical models will help to alleviate such ambiguities in current theoretical modeling. Additionally, they will allow to study new phenomena, to open qualitatively new windows for stellar astrophysics in the Gaia-era., Comment: Proceedings of "The Three Dimensional Universe With Gaia", Paris, October 4-7, 2004, 4 pages, 3 figures

Since about 25 years ago work has been dedicated to the development of hydrodynamical model atmospheres for cool stars (of A to T spectral type). Despite their obviously sounder physical foundation in comparison with standard hydrostatic models, their general application has been rather limited. In order to understand why this is, and how to progress, we review the present status of hydrodynamical modelling of cool star atmospheres. The development efforts were and are motivated by the theoretical interest of understanding the dynamical processes operating in stellar atmospheres. To show the observational impact, we discuss examples in the fields of spectroscopy and stellar structure where hydrodynamical modelling provided results on a level qualitatively beyond standard models. We stress present modelling challenges, and highlight presently possible and future observations that would be particularly valuable in the interplay between model validation and interpretation of observables, to eventually widen the usage of hydrodynamical model atmospheres within the astronomical community., Comment: Proc. 13th Cool Stars Workshop, Hamburg, 5-9 July 2004, F.Favata et al., eds.,7 pages, 5 figures

We used a series of COBOLD hydrodynamical model atmospheres covering stellar objects from white dwarfs to red giants to derive theoretical estimates of the photometric and photocentric stellar variability in wavelength-integrated light across the Hertzsprung-Russell diagram. We validated our models against solar measurements from the SOHO/VIRGO instrument. Within our set of models we find a systematic increase of the photometric as well as photocentric variability - which turn out to be closely connected - with decreasing surface gravity. The estimated absolute levels of the photocentric variability do not affect astrometric observations on a precision level expected to be achieved by the GAIA mission - with the exception of close-by giants. The case of supergiants remains to be investigated. In view of the ongoing debate about the photometric non-detection of p-modes in Procyon by the Canadian MOST satellite we remark that we obtain a factor of ca. 3 in amplitude between the granular background noise in the Sun and Procyon. This statement refers to a particular representation of temporal power spectra as discussed in Sect. 5., Comment: Proc. 13th Cool Stars Workshop, Hamburg, 5-9 July 2004, F.Favata et al., eds., 5 pages, 8 figures

Magnetic activity in M-dwarfs present enigmatic questions: On the one hand they have higher field strengths and larger filling factors than the magnetic field on the Sun, on the other hand, they are fully convective and their atmospheres are more neutral, hence they do not have an undershoot layer for magnetic flux storage and as we show here, cannot have small-scale dynamo action in their photospheres either. We present a discussion of these facts and propose a new numerical model to investigate M-dwarf magnetism., Comment: 3 pages, 4 figures, submitted to the proceedings of Cool Stars 13, to appear in an ESA-SP volume

Spectral-line displacements away from the wavelengths naively expected from the Doppler shift caused by stellar radial motion may originate as convective shifts (correlated velocity and brightness patterns in the photosphere), as gravitational redshifts, or perhaps be induced by wave motions. Absolute lineshifts, in the past studied only for the Sun, are now accessible also for other stars thanks to astrometric determination of stellar radial motion, and spectrometers with accurate wavelength calibration. Comparisons between spectroscopic apparent radial velocities and astrometrically determined radial motions reveal greater spectral blueshifts in F-type stars than in the Sun (as theoretically expected from their more vigorous convection), further increasing in A-type stars (possibly due to atmospheric shockwaves). An important near-future development to enable a further analysis of stellar surface structure will be the study of wavelength variations across spatially resolved stellar disks, e.g., the center-to-limb wavelength changes along a stellar diameter, and their spatially resolved time variability., Comment: Proc. 13th Cool Stars Workshop, Hamburg, 5-9 July 2004, F.Favata et al., eds. 7 pages, 8 figures

We constructed hydrodynamical model atmospheres for mid M-type main-, as well as pre-main-sequence objects. Despite the complex chemistry encountered in such cool atmospheres a reasonably accurate representation of the radiative transfer is possible. The detailed treatment of the interplay between radiation and convection in the hydrodynamical models allows to study processes usually not accessible within the framework conventional model atmospheres. In particular, we determined the efficiency of the convective energy transport, and the efficiency of mixing by convective overshoot. The convective transport efficiency expressed in terms of an equivalent mixing-length parameter amounts to values around ~2 in the optically thick, and ~2.8 in the optically thin regime. The thermal structure of the formally convectively stable layers is little affected by convective overshoot and wave heating, i.e. stays close to radiative equilibrium. Mixing by convective overshoot shows an exponential decline with geometrical distance from the Schwarzschild stability boundary. The scale height of the decline varies with gravitational acceleration roughly as g^(-1/2), with 0.5 pressure scale heights at log(g)=5.0., Comment: Proceedings of IAU Symp. 210. 12 pages, 8 figures (2 color)

Based on detailed 2D and 3D numerical radiation-hydrodynamics (RHD) simulations of time-dependent compressible convection, we have studied the dynamics and thermal structure of the convective surface layers of a prototypical late-type M-dwarf (Teff~2800K log(g)=5.0, solar chemical composition). The RHD models predict stellar granulation qualitatively similar to the familiar solar pattern. Quantitatively, the granular cells show a convective turn-over time scale of ~100s, and a horizontal scale of 80km; the relative intensity contrast of the granular pattern amounts to 1.1%, and root-mean-square vertical velocities reach 240m/s at maximum. Deviations from radiative equilibrium in the higher, formally convectively stable atmospheric layers are found to be insignificant allowing a reliable modeling of the atmosphere with 1D standard model atmospheres. A mixing-length parameter of alpha=2.1 provides the best representation of the average thermal structure of the RHD model atmosphere while alternative values are found when fitting the asymptotic entropy encountered in deeper layers of the stellar envelope alpha=1.5, or when matching the vertical velocity field alpha=3.5. The close correspondence between RHD and standard model atmospheres implies that presently existing discrepancies between observed and predicted stellar colors in the M-dwarf regime cannot be traced back to an inadequate treatment of convection in the 1D standard models. The RHD models predict a modest extension of the convectively mixed region beyond the formal Schwarzschild stability boundary which provides hints for the distribution of dust grains in cooler (brown dwarf) atmospheres., Comment: 19 pages, 16 figures, accepted for publication in A&A

We performed kinematic studies of the evolution of small-scale magnetic fields in the surface laters of M-dwarfs. We solved the induction equation for a prescribed velocity field, magnetic Reynolds number, and boundary conditions in a Cartesian box, representing a volume comprising the optically thin stellar atmosphere and the uppermost part of the optically thick convective envelope. The velocity field is spatially and temporally variable, and stems from detailed radiation-hydrodynamics simulations of convectione flows in a proto-typical late-type M-dwarf. We find dynamo action for large magnetic Reynolds numbers. Growth time scales of the magnetic field is comparable to the convective turn-over time scale (approximately 150 seconds). The convective velocity field concentrates the magnetic field in sheets and tubular structures in the inter-granular down-flows. Scaling from solar conditions suggests that field strengths as high as 20 kG might be reached locally. Perhps surprisingly, the magnetic Reynolds number is of the order unity in the surface layers of cooler M-dwarfs, rendering the dynamo inoperative. In all studied cases we find a rather low spatial filling factor of the magnetic field., Comment: 5 pages, 7 figures, submitted to Astron. Nachr. (from 1st Potsdam Thinkshop on Sunspots and Starspots)

We present here numerical simulations of surface solar convection which cover a box of 30$\times30\times$3.2 Mm$^3$ with a resolution of 315$\times315\times$82, which is used to investigate the dynamics of scales larger than granulation. No structure resembling supergranulation is present; possibly higher Reynolds numbers (i.e. higher numerical resolution), or magnetic fields, or greater depth are necessary. The results also show interesting aspects of granular dynamics which are briefly presented, like extensive p-mode ridges in the k-$\omega$ diagram and a ringlike distribution of horizontal vorticity around granules. At large scales, the horizontal velocity is much larger than the vertical velocity and the vertical motion is dominated by p-mode oscillations., Comment: Contribution to the proceedings of the workshop entitled "THEMIS and the new frontiers of solar atmosphere dynamics" (March 2001), 6 pages, to appear in Nuovo Cimento C

Using a numerical simulation of compressible convection with radiative transfer mimicking the solar photosphere, we compare the velocity field derived from granule motions to the actual velocity field of the plasma. We thus test the idea that granules may be used to trace large-scale velocity fields at the sun's surface. Our results show that this is indeed the case provided the scale separation is sufficient. We thus estimate that neither velocity fields at scales less than 2500 km nor time evolution at scales shorter than 0.5 hr can be faithfully described by granules. At larger scales the granular motions correlate linearly with the underlying fluid motions with a slope of ~< 2 reaching correlation coefficients up to ~0.9., Comment: 4 pages - accepted in Astronomy and Astrophysics

The computationally demanding nature of radiative-hydrodynamical simulations of stellar surface convection warrants an investigation of the sensitivity of the convective structure and spectral synthesis to the numerical resolution and dimension of the simulations, which is presented here. With too coarse a resolution the predicted spectral lines tend to be too narrow, reflecting insufficient Doppler broadening from the convective motions, while at the currently highest affordable resolution the line shapes have converged essentially perfectly to the observed profiles. Similar conclusions are drawn from the line asymmetries and shifts. In terms of abundances, weak FeI and FeII lines show a very small dependence (~0.02 dex) while for intermediate strong lines with significant non-thermal broadening the sensitivity increases (~0.10 dex). Problems arise when using 2D convection simulations to describe an inherent 3D phenomenon, which translates to inaccurate atmospheric velocity fields and temperature and pressure structures. In 2D the theoretical line profiles tend to be too shallow and broad compared with the 3D calculations and observations, in particular for intermediate strong lines. In terms of abundances, the 2D results are systematically about 0.1 dex lower than for the 3D case for FeI lines. Furthermore, the predicted line asymmetries and shifts are much inferior in 2D. Given these shortcomings and computing time considerations it is better to use 3D simulations of even modest resolution than high-resolution 2D simulations., Comment: Accepted for A&A

Based on detailed 2D numerical radiation hydrodynamics (RHD) calculations of time-dependent compressible convection, we have studied the dynamics and thermal structure of the convective surface layers of solar-type stars. The RHD models provide information about the convective efficiency in the superadiabatic region at the top of convective envelopes and predict the asymptotic value of the entropy of the deep, adiabatically stratified layers. This information is translated into an effective mixing-length parameter alpha suitable to construct standard stellar structure models. We validate the approach by a detailed comparison to helioseismic data. The grid of RHD models for solar metallicity comprises 58 simulation runs with a helium abundance of Y=0.28 in the range of effective temperatures 4300K < Teff < 7100K and gravities 2.54 < logg < 4.74. We find a moderate, nevertheless significant variation of alpha between about 1.3 for F-dwarfs and 1.75 for K-subgiants with a dominant dependence on Teff. In the close neighbourhood of the Sun we find a plateau where alpha remains almost constant. The internal accuracy of the calibration of alpha is estimated to be +/- 0.05 with a possible systematic bias towards lower values. An analogous calibration of the convection theory of Canuto & Mazzitelli gives a different temperature dependence but a similar variation of the free parameter. For the first time, values for the gravity-darkening exponent beta are derived independently of mixing-length theory: beta=0.07...0.10. We show that our findings are consistent with constraints from stellar stability considerations and provide compact fitting formulae for the calibrations., Comment: 14 pages, 6 figures, 4 tables

A CO(1-0) on-source spectrum of the well-known carbon star IRC+10216 taken with the IRAM 30m telescope in June 1996 shows excess emission between -18.3 and -14.3 km/s at the red wing of the underlying profile. The excess emission is confirmed in January 1997 but is absent in April 1997 and June 1998 IRAM spectra. Such a transient feature has not been seen before in this star or any other AGB star. In April 1997 we mapped the circumstellar shell out to 110". Both the J = 1-0 and J = 2-1 spectra show ``spikes'' or components which vary in strength with position in the envelope. One of these components corresponds to the velocity interval mentioned above. An immediate conclusion is that the circumstellar shell is not spherically symmetric, contrary to what was believed based on lower spectral resolution data. We are probably seeing emission from a complex geometrical structure. Neither a disk structure nor a double-wind structure seem to be able to explain the observations. The on-source transient behaviour of the red excess emission can reasonably well be explained by a single large (about 5 x 10^{13} cm) blob, that expands due to internal motion., Comment: Astronomy & Astrophysics, accepted. Incorrect figure 2 replaced

Based on radiation hydrodynamics modeling of stellar convection zones, a diffusion scheme has been devised describing the downward penetration of convective motions beyond the Schwarzschild boundary (overshoot) into the radiative interior. This scheme of exponential diffusive overshoot has already been successfully applied to AGB stars. Here we present an application to the Sun in order to determine the time scale and depth extent of this additional mixing, i.e. diffusive overshoot at the base of the convective envelope. We calculated the associated destruction of lithium during the evolution towards and on the main-sequence. We found that the slow-mixing processes induced by the diffusive overshoot may lead to a substantial depletion of lithium during the Sun's main-sequence evolution., Comment: 8 pages including 4 PostScript figures, using spackap.sty, Proc. of the workshop "Solar Composition and its Evolution - From Core to Corona, to appear in Space Science Reviews, also available from http://www.speckle.mpifr-bonn.mpg.de/publications.html