BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

B. Freytag

B. Freytag appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.IMastro-ph.GAastro-ph

Catalog footprint

What is connected

34works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

Probing Red Supergiant dynamics through photo-center displacements measured by Gaia

Red supergiant (RSGs) are cool massive stars in a late phase of their evolution when the stellar envelope becomes fully convective. They are the brightest stars in the universe at infrared light and can be detected in galaxies far beyond the Local Group, allowing for accurate determination of chemical composition of galaxies. The study of their physical properties is extremely important for various phenomena including the final fate of massive stars as type II supernovae and gravitational wave progenitors. We explore the well-studied nearby young stellar cluster chi Per. Using Gaia EDR3 data, we find the distance of the cluster (d = 2.260+-0.020 kpc). We then investigate the variability of the convection-related surface structure as a source for parallax measurement uncertainty. We use state-of-the-art 3D radiative hydrodynamics simulations with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps in the Gaia G photometric system. We calculate the variabiltiy, as a function of time, of the intensity-weighted mean from the synthetic maps. We then select the RSG stars in the cluster and compare their uncertainty on parallaxes to the predictions of photocentre displacements. The synthetic maps of RSG show extremely irregular and temporal variable surfaces due to convection-related dynamics. Consequentially, the position of the photo-center varies during Gaia measurements between 0.033 and 0.130 AU (up to 5% of the corresponding simulation stellar radius). We argue that the variability of the convection-related surface structures accounts for a substantial part of the Gaia EDR3 parallax error of the RSG sample. We suggest that the variation of the uncertainty on Gaia parallax could be exploited quantitatively using appropriate 3D simulations to extract, in a unique way, important information about the stellar dynamics and parameters of RSG stars.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 μm), M (4.5 to 5 μm), and N (8 to 13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics (RHD) simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes Optim3D and Radmc3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of the stellar surface and circumstellar environment across a very large spectral domain for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optical interferometry and Gaia measurement uncertainties reveal the physics of asymptotic giant branch stars

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

VLTI-PIONIER imaging of the red supergiant V602 Carinae

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrmμ$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Enhanced methods for computing spectra from CO5BOLD models using Linfor3D. Molecular bands in metal-poor stars

Molecular features such as the G-band, CN-band and NH-band are important diagnostics for measuring a star's carbon and nitrogen abundances, especially in metal-poor stars where atomic lines are no longer visible in stellar spectra. Unlike atomic transitions, molecular features tend to form in bands, which cover large wavelength regions in a spectrum. While it is a trivial matter to compute carbon and nitrogen molecular bands under the assumption of 1D, it is extremely time consuming in 3D. In this contribution to the 2016 CO5BOLD workshop we review the improvements made to the 3D spectral synthesis code Linfor3D, and discuss the new challenges found when computing molecular features in 3D.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Near-infrared spectro-interferometry of Mira variables and comparisons to 1D dynamic model atmospheres and 3D convection simulations

We obtained a total of 20 near-infrared K-band spectro-interferometric snapshot observations of the Mira variables o Cet, R Leo, R Aqr, X Hya, W Vel, and R Cnc with a spectral resolution of about 1500. We compared observed flux and visibility spectra with predictions by CODEX 1D dynamic model atmospheres and with azimuthally averaged intensities based on CO5BOLD 3D dynamic model atmospheres including convection. Our visibility data confirm the presence of spatially extended molecular atmospheres located above the continuum radii with large-scale inhomogeneities or clumps that contribute a few percent of the total flux. The detailed structure of the inhomogeneities or clumps show a variability on time scales of 3 months and above. Both modeling attempts provided satisfactory fits to our data. In particular, they are both consistent with the observed decrease in the visibility function at molecular bands of water vapor and CO, indicating a spatially extended molecular atmosphere. Observational variability phases are mostly consistent with those of the best-fit CODEX models, except for near-maximum phases, where data are better described by near-minimum models. Rosseland angular diameters derived from the model fits are broadly consistent between those based on the 1D and the 3D models and with earlier observations. We derived fundamental parameters including absolute radii, effective temperatures, and luminosities for our sources. Our results provide a first observational support for theoretical results that shocks induced by convection and pulsation in the 3D CO5BOLD models of AGB stars are roughly spherically expanding and of similar nature to those of self-excited pulsations in 1D CODEX models. Unlike for red supergiants, the pulsation- and shock-induced dynamics can levitate the molecular atmospheres of Mira variables to extensions that are consistent with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

3D Model Atmospheres for Extremely Low-Mass White Dwarfs

We present an extended grid of mean three-dimensional (3D) spectra for low-mass, pure-hydrogen atmosphere DA white dwarfs (WDs). We use CO5BOLD radiation-hydrodynamics 3D simulations covering Teff = 6000-11,500 K and logg = 5-6.5 (cgs units) to derive analytical functions to convert spectroscopically determined 1D temperatures and surface gravities to 3D atmospheric parameters. Along with the previously published 3D models, the 1D to 3D corrections are now available for essentially all known convective DA WDs (i.e., logg = 5-9). For low-mass WDs, the correction in temperature is relatively small (a few per cent at the most), but the surface gravities measured from the 3D models are lower by as much as 0.35 dex. We revisit the spectroscopic analysis of the extremely low-mass (ELM) WDs, and demonstrate that the 3D models largely resolve the discrepancies seen in the radius and mass measurements for relatively cool ELM WDs in eclipsing double WD and WD + milli-second pulsar binary systems. We also use the 3D corrections to revise the boundaries of the ZZ Ceti instability strip, including the recently found ELM pulsators.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On The Evolution of Magnetic White Dwarfs

We present the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarf stars. We demonstrate that convective energy transfer is seriously impeded by magnetic fields when the plasma-beta parameter, the thermal to magnetic pressure ratio, becomes smaller than unity. The critical field strength that inhibits convection in the photosphere of white dwarfs is in the range B = 1-50 kG, which is much smaller than the typical 1-1000 MG field strengths observed in magnetic white dwarfs, implying that these objects have radiative atmospheres. We have then employed evolutionary models to study the cooling process of high-field magnetic white dwarfs, where convection is entirely suppressed during the full evolution (B > 10 MG). We find that the inhibition of convection has no effect on cooling rates until the effective temperature (Teff) reaches a value of around 5500 K. In this regime, the standard convective sequences start to deviate from the ones without convection owing to the convective coupling between the outer layers and the degenerate reservoir of thermal energy. Since no magnetic white dwarfs are currently known at the low temperatures where this coupling significantly changes the evolution, effects of magnetism on cooling rates are not expected to be observed. This result contrasts with a recent suggestion that magnetic white dwarfs with Teff < 10,000 K cool significantly slower than non-magnetic degenerates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pathways for Observing Stellar Surfaces Using 3D Hydrodynamical Simulations of Evolved Stars

Evolved stars are among the largest and brightest stars and they are ideal targets for the new generation of sensitive, high resolution instrumentation that provides spectrophotometric, interferometric, astrometric, and imaging observables. The interpretation of the complex stellar surface images requires numerical simulations of stellar convection that take into account multi-dimensional time-dependent radiation hydrodynamics with realistic input physics. We show how the evolved star simulations are obtained using the radiative hydrodynamics code CO5BOLD and how the accurate observables are computed with the post-processing radiative transfer code Optim3D. The synergy between observations and theoretical work is supported by a proper and quantitative analysis using these simulations, and by strong constraints from the observational side.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

What causes the large extensions of red-supergiant atmospheres? Comparisons of interferometric observations with 1-D hydrostatic, 3-D convection, and 1-D pulsating model atmospheres

We present the atmospheric structure and the fundamental parameters of three red supergiants, increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible mechanisms that may explain the large observed atmospheric extensions of RSGs. We carried out spectro-interferometric observations of 3 RSGs in the near-infrared K-band with the VLTI/AMBER instrument at medium spectral resolution. To comprehend the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D convection, and new 1-D self-excited pulsation models of RSGs. Our near-infrared flux spectra are well reproduced by the PHOENIX model atmospheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, the PHOENIX model visibilities do not predict the large observed extensions of molecular layers, most remarkably in the CO bands. Likewise, the 3-D convection models and the 1-D pulsation models with typical parameters of RSGs lead to compact atmospheric structures as well, which are similar to the structure of the hydrostatic PHOENIX models. They can also not explain the observed decreases in the visibilities and thus the large atmospheric molecular extensions. The full sample of our RSGs indicates increasing observed atmospheric extensions with increasing luminosity and decreasing surface gravity, and no correlation with effective temperature or variability amplitude, which supports a scenario of radiative acceleration on Doppler-shifted molecular lines.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

3D hydrodynamical simulations of evolved stars and observations of stellar surfaces

Evolved stars are among the largest and brightest stars and they are ideal targets for the new generation of sensitive, high resolution instrumentation that pro- vides spectrophotometric, interferometric, astrometric, and imaging observables. The interpretation of the complex stellar surface images requires numerical simulations of stellar convection that take into account multi-dimensional time-dependent radiation hydrodynamics with realistic input physics. We show how the evolved star simulations are obtained using the radiative hydrodynamics code CO5BOLD and how the accurate observables are computed with the post-processing radiative transfer code Optim3D. The synergy between observations and theoretical work is supported by a proper and quantitative analysis using these simulations, and by strong constraints from the obser- vational side.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Calibration of the Mixing-Length Theory for Convective White Dwarf Envelopes

A calibration of the mixing-length parameter in the local mixing-length theory (MLT) is presented for the lower part of the convection zone in pure-hydrogen atmosphere white dwarfs. The parameterization is performed from a comparison of 3D CO5BOLD simulations with a grid of 1D envelopes with a varying mixing-length parameter. In many instances, the 3D simulations are restricted to the upper part of the convection zone. The hydrodynamical calculations suggest, in those cases, that the entropy of the upflows does not change significantly from the bottom of the convection zone to regions immediately below the photosphere. We rely on this asymptotic entropy value, characteristic of the deep and adiabatically stratified layers, to calibrate 1D envelopes. The calibration encompasses the convective hydrogen-line (DA) white dwarfs in the effective temperature range 6000 < Teff (K) < 15,000 and the surface gravity range 7.0 < log g < 9.0. It is established that the local MLT is unable to reproduce simultaneously the thermodynamical, flux, and dynamical properties of the 3D simulations. We therefore propose three different parameterizations for these quantities. The resulting calibration can be applied to structure and envelope calculations, in particular for pulsation, chemical diffusion, and convective mixing studies. On the other hand, convection has no effect on the white dwarf cooling rates until there is a convective coupling with the degenerate core below Teff ~ 5000 K. In this regime, the 1D structures are insensitive to the MLT parameterization and converge to the mean 3D results, hence remain fully appropriate for age determinations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

White Dwarfs In The UKIRT Infrared Deep Sky Survey Data Release 9

We have identified eight to ten new cool white dwarfs from the Large Area Survey (LAS) Data Release 9 of the United Kingdom InfraRed Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS). The data set was paired with the Sloan Digital Sky Survey (SDSS) to obtain proper motions and a broad ugrizYJHK wavelength coverage. Optical spectroscopic observations were secured at Gemini Observatory and confirm the degenerate status for eight of our targets. The final sample includes two additional white dwarf candidates with no spectroscopic observations. We rely on improved 1D model atmospheres and new multi-dimensional simulations with CO5BOLD to review the stellar parameters of the published LAS white dwarf sample along with our additional discoveries. Most of the new objects possess very cool atmospheres with effective temperatures below 5000 K, including two pure-hydrogen remnants with a cooling age between 8.5 and 9.0 Gyr, and tangential velocities in the range 40 km/s < vtan < 60 km/s. They are likely thick disk 10-11 Gyr-old objects. In addition we find a resolved double degenerate system with vtan ~ 155 km/s and a cooling age between 3.0 and 5.0 Gyr. These white dwarfs could be disk remnants with a very high velocity or former halo G stars. We also compare the LAS sample with earlier studies of very cool degenerates and observe a similar deficit of helium-dominated atmospheres in the range 5000 < Teff (K) < 6000. We review the possible explanations for the spectral evolution from helium-dominated towards hydrogen-rich atmospheres at low temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Convective line shifts for the Gaia RVS from the CIFIST 3D model atmosphere grid

To derive space velocities of stars along the line of sight from wavelength shifts in stellar spectra requires accounting for a number of second-order effects. For most stars, gravitational redshifts, convective blueshifts, and transverse stellar motion are the dominant contributors. We provide theoretical corrections for the net velocity shifts due to convection expected for the measurements from the Gaia Radial Velocity Spectrometer (RVS). We used a set of three-dimensional time-dependent simulations of stellar surface convection computed with CO5BOLD to calculate spectra of late-type stars in the Gaia RVS range and to infer the net velocity offset that convective motions will induce in radial velocities derived by cross-correlation. The net velocity shifts derived by cross-correlation depend both on the wavelength range and spectral resolution of the observations. Convective shifts for Gaia RVS observations are less than 0.1 km/s for late-K-type stars, and they increase with stellar mass, reaching about 0.3 km/s or more for early F-type dwarfs. This tendency is the result of an increase with effective temperature in both temperature and velocity fluctuations in the line-forming region. Our simulations also indicate that the net RVS convective shifts can be positive (i.e. redshifts) in some cases. Overall, the blueshifts weaken slightly with increasing surface gravity, and are enhanced at low metallicity. Gravitational redshifts amount up to 0.7 km/s and dominate convective blueshifts for dwarfs, but become much weaker for giants.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Granulation properties of giants, dwarfs, and white dwarfs from the CIFIST 3D model atmosphere grid

3D model atmospheres for giants, dwarfs, and white dwarfs, computed with the CO5BOLD code and part of the CIFIST grid, have been used for spectroscopic and asteroseismic studies. Unlike existing plane-parallel 1D structures, these simulations predict the spatially and temporally resolved emergent intensity so that granulation can be analysed, which provides insights on how convective energy transfer operates in stars. The wide range of atmospheric parameters of the CIFIST 3D simulations (3600 < Teff (K) < 13,000 and 1 < log g < 9) allows the comparison of convective processes in significantly different environments. We show that the relative intensity contrast is correlated with both the Mach and Peclet numbers in the photosphere. The horizontal size of granules varies between 3 and 10 times the local pressure scale height, with a tight correlation between the factor and the Mach number of the flow. Given that convective giants, dwarfs, and white dwarfs cover the same range of Mach and Peclet numbers, we conclude that photospheric convection operates in a very similar way in those objects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Oxygen spectral line synthesis: 3D non-LTE with CO5BOLD hydrodynamical model atmospheres

In this work we present first results of our current project aimed at combining the 3D hydrodynamical stellar atmosphere approach with non-LTE (NLTE) spectral line synthesis for a number of key chemical species. We carried out a full 3D-NLTE spectrum synthesis of the oxygen IR 777 nm triplet, using a modified and improved version of our NLTE3D package to calculate departure coefficients for the atomic levels of oxygen in a CO5BOLD 3D hydrodynamical solar model atmosphere. Spectral line synthesis was subsequently performed with the Linfor 3D code. In agreement with previous studies, we find that the lines of the oxygen triplet produce deeper cores under NLTE conditions, due to the diminished line source function in the line forming region. This means that the solar oxygen IR 777 nm lines should be stronger in NLTE, leading to negative 3D NLTE-LTE abundance corrections. Qualitatively this result would support previous claims for a relatively low solar oxygen abundance. Finally, we outline several further steps that need to be taken in order to improve the physical realism and numerical accuracy of our current 3D-NLTE calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Progress in Modeling Very Low Mass Stars, Brown Dwarfs, and Planetary Mass Objects

We review recent advancements in modeling the stellar to substellar transition. The revised molecular opacities, solar oxygen abundances and cloud models allow to reproduce the photometric and spectroscopic properties of this transition to a degree never achieved before, but problems remain in the important M-L transition characteristic of the effective temperature range of characterizable exoplanets. We discuss of the validity of these classical models. We also present new preliminary global Radiation HydroDynamical M dwarfs simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pure-hydrogen 3D model atmospheres of cool white dwarfs

A sequence of pure-hydrogen CO5BOLD 3D model atmospheres of DA white dwarfs is presented for a surface gravity of log g = 8 and effective temperatures from 6000 to 13,000 K. We show that convective properties, such as flow velocities, characteristic granulation size and intensity contrast of the granulation patterns, change significantly over this range. We demonstrate that these 3D simulations are not sensitive to numerical parameters unlike the 1D structures that considerably depend on the mixing-length parameters. We conclude that 3D spectra can be used directly in the spectroscopic analyses of DA white dwarfs. We confirm the result of an earlier preliminary study that 3D model spectra provide a much better characterization of the mass distribution of white dwarfs and that shortcomings of the 1D mixing-length theory are responsible for the spurious high-log g determinations of cool white dwarfs. In particular, the 1D theory is unable to account for the cooling effect of the convective overshoot in the upper atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectroscopic analysis of DA white dwarfs with 3D model atmospheres

We present the first grid of mean three-dimensional (3D) spectra for pure-hydrogen (DA) white dwarfs based on 3D model atmospheres. We use CO5BOLD radiation-hydrodynamics 3D simulations instead of the mixing-length theory for the treatment of convection. The simulations cover the effective temperature range of 6000 < Teff (K) < 15,000 and the surface gravity range of 7 < log g < 9 where the large majority of DAs with a convective atmosphere are located. We rely on horizontally averaged 3D structures (over constant Rosseland optical depth) to compute <3D> spectra. It is demonstrated that our <3D> spectra can be smoothly connected to their 1D counterparts at higher and lower Teff where the 3D effects are small. Analytical functions are provided in order to convert spectroscopically determined 1D effective temperatures and surface gravities to 3D atmospheric parameters. We apply our improved models to well studied spectroscopic data sets from the Sloan Digital Sky Survey and the White Dwarf Catalog. We confirm that the so-called high-log g problem is not present when employing <3D> spectra and that the issue was caused by inaccuracies in the 1D mixing-length approach. The white dwarfs with a radiative and a convective atmosphere have derived mean masses that are the same within ~0.01 Msun, in much better agreement with our understanding of stellar evolution. Furthermore, the 3D atmospheric parameters are in better agreement with independent Teff and log g values from photometric and parallax measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Imaging the heart of astrophysical objects with optical long-baseline interferometry

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was "blind" in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of the solar near-surface layers with the CO5BOLD, MURaM, and Stagger codes

Radiative hydrodynamic simulations of solar and stellar surface convection have become an important tool for exploring the structure and gas dynamics in the envelopes and atmospheres of late-type stars and for improving our understanding of the formation of stellar spectra. We quantitatively compare results from three-dimensional, radiative hydrodynamic simulations of convection near the solar surface generated with three numerical codes CO5BOLD, MURaM, and STAGGER and different simulation setups in order to investigate the level of similarity and to cross-validate the simulations. For all three simulations, we considered the average stratifications of various quantities (temperature, pressure, flow velocity, etc.) on surfaces of constant geometrical or optical depth, as well as their temporal and spatial fluctuations. We also compared observables, such as the spatially resolved patterns of the emerging intensity and of the vertical velocity at the solar optical surface as well as the center-to-limb variation of the continuum intensity at various wavelengths. The depth profiles of the thermodynamical quantities and of the convective velocities as well as their spatial fluctuations agree quite well. Slight deviations can be understood in terms of differences in box size, spatial resolution and in the treatment of non-gray radiative transfer between the simulations. The results give confidence in the reliability of the results from comprehensive radiative hydrodynamic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiative hydrodynamics simulations of red supergiant stars. IV gray versus non-gray opacities

Red supergiants are massive evolved stars that contribute extensively to the chemical enrichment of our Galaxy. It has been shown that convection in those stars gives rise to large granules that cause surface inhomogeneities and shock waves in the photosphere. The understanding of their dynamics is crucial to unveil the unknown mass-loss mechanism, their chemical composition and stellar parameters. We present a new generation of red supergiants simulations with a more sophisticated opacity treatment done with 3D radiative- hydrodynamics CO5BOLD. In the code, the coupled equations of compressible hydrodynamics and non-local radiation transport are solved in the presence of a spherical potential. The stellar core is replaced by a special spherical inner boundary condition, where the gravitational potential is smoothed and the energy production by fusion is mimicked by a simply producing heat corresponding to the stellar luminosity. The post-processing radiative transfer code OPTIM3D is used to extract spectroscopic and interferometric observables. We show that the relaxation of the assumption of frequency-independent opacities shows a steeper mean thermal gradient in the optical thin region that affect strongly the atomic strengths and the spectral energy distribution. Moreover, the weaker temperature fluctuations reduce the incertitude on the radius determination with interferometry. We show that 1D models of red supergiants must include a turbulent velocity calibrated on 3D simulations to obtain the effective surface gravity that mimic the effect of turbulent pressure on the stellar atmosphere. We provide an empirical calibration of the ad-hoc micro- and macroturbulence parameters for 1D models using the 3D simulations: we find that there is not a clear distinction between the different macroturbulent profiles needed in 1D models to fit 3D synthetic lines.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Solution to the problem of the surface gravity distribution of cool DA white dwarfs from improved 3D model atmospheres

The surface gravities of cool (Teff < 13,000 K) hydrogen-atmosphere DA white dwarfs, determined from spectroscopic analyses, are found to be significantly higher than the canonical value of log g ~ 8 expected for these stars. It was recently concluded that a problem with the treatment of convective energy transport within the framework of the mixing-length theory was the most plausible explanation for this high-log g problem. We pursue the investigation of this discrepancy by computing model spectra of cool convective white dwarfs from a small sequence (11,300 K < Teff < 12,800 K) of 3D hydrodynamical model atmospheres, which feature a sophisticated treatment of convection and radiative transfer. Our approach is to proceed with a differential analysis between 3D and standard 1D models. We find that the 3D spectra predict significantly lower surface gravities, with corrections of the right amplitude as a function of effective temperature to obtain values of log g ~ 8 on average. We conclude that the surface gravity distribution of cool convective DA white dwarfs is much closer to that of hotter radiative objects when using, for the treatment of the convection, 3D models instead of the mixing-length framework.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galactic Evolution of Oxygen: OH lines in 3D hydrodynamical model atmospheres

The measurement of oxygen lines in metal-poor unevolved stars, in particular near-UV OH lines, can provide invaluable information on the properties of the Early Galaxy. Near-UV OH lines constitute an important tool to derive oxygen abundances in metal-poor dwarf stars. Therefore, it is important to correctly model the line formation of OH lines, especially in metal-poor stars, where 3D hydrodynamical models commonly predict cooler temperatures than plane-parallel hydrostatic models in the upper photosphere. We have made use of a grid of 52 3D hydrodynamical model atmospheres for dwarf stars computed using the code CO5BOLD, extracted from the more extended CIFIST grid. The 52 models cover the effective temperature range 5000-6500K, the surface gravity range 3.5-4.5 and the metallicity range -3<[Fe/H]<0. We determine 3D-LTE abundance corrections in all the 52 3D models for several OH lines and FeI lines of different excitation potentials. These 3D-LTE corrections are generally negative reaching values of roughly -1 dex (for the OH 3167 with excitation potential of approximately 1 eV) for the higher temperatures and surface gravities. We apply these 3D-LTE corrections to the individual O abundances, derived from OH lines, of a sample the metal-poor dwarf stars reported in Israelian et al.(1998, 2001) and Boesgaard et al.(1999), by interpolating the stellar parameters of the dwarfs in the grid of 3D-LTE corrections. The new 3D-LTE [O/Fe] ratio still keeps a similar trend as the 1D-LTE, i.e, increasing towards lower [Fe/H] values. We applied 1D-NLTE corrections to 3D FeI abundances and we still see an increasing [O/Fe] ratio towards lower metallicites. However, the Galactic [O/Fe] ratio must be revisited once 3D-NLTE corrections become available for OH and Fe lines for a grid of 3D hydrodynamical model atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Model Atmospheres From Very Low Mass Stars to Brown Dwarfs

Since the discovery of brown dwarfs in 1994, and the discovery of dust cloud formation in the latest Very Low Mass Stars (VLMs) and Brown Dwarfs (BDs) in 1996, the most important challenge in modeling their atmospheres as become the understanding of cloud formation and advective mixing. For this purpose, we have developed radiation hydrodynamic 2D model atmosphere simulations to study the formation of forsterite dust in presence of advection, condensation, and sedimentation across the M-L-T VLMs to BDs sequence (Teff = 2800 K to 900 K, Freytag et al. 2010). We discovered the formation of gravity waves as a driving mechanism for the formation of clouds in these atmospheres, and derived a rule for the velocity field versus atmospheric depth and Teff , which is relatively insensitive to gravity. This rule has been used in the construction of the new model atmosphere grid, BT-Settl, to determine the microturbulence velocity, the diffusion coefficient, and the advective mixing of molecules as a function of depth. This new model grid of atmospheres and synthetic spectra has beencomputedfor100,000K>Teff >400K,5.5>logg>-0.5,and[M/H]=+0.5to -1.5, and the reference solar abundances of Asplund et al. (2009). We found that the new solar abundances allow an improved (close to perfect) reproduction of the photo- metric and spectroscopic VLMs properties, and, for the first time, a smooth transition between stellar and substellar regimes -- unlike the transition between the NextGen models from Hauschildt et al. 1999a,b, and the AMES-Dusty models from Allard et al. 2001). In the BDs regime, the BT-Settl models propose an improved explanation for the M-L-T spectral transition. In this paper, we therefore present the new BT-Settl model atmosphere grid, which explains the entire transition from the stellar to planetary mass regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiative hydrodynamics simulations of red supergiant stars: II. simulations of convection on Betelgeuse match interferometric observations

Context. The red supergiant (RSG) Betelgeuse is an irregular variable star. Convection may play an important role in understanding this variability. Interferometric observations can be interpreted using sophisticated simulations of stellar convection. Aims. We compare the visibility curves and closure phases obtained from our 3D simulation of RSG convection with CO5BOLD to various interferometric observations of Betelgeuse from the optical to the H band in order to characterize and measure the convection pattern on this star. Methods. We use 3D radiative-hydrodynamics (RHD) simulation to compute intensity maps in different filters and we thus derive interferometric observables using the post-processing radiative transfer code OPTIM3D. The synthetic visibility curves and closure phases are compared to observations. Results. We provide a robust detection of the granulation pattern on the surface of Betelgeuse in the optical and in the H band based on excellent fits to the observed visibility points and closure phases. Moreover, we determine that the Betelgeuse surface in the H band is covered by small to medium scale (5-15 mas) convection-related surface structures and a large (30 mas) convective cell. In this spectral region, H2O molecules are the main absorbers and contribute to the small structures and to the position of the first null of the visibility curve (i.e. the apparent stellar radius).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiative hydrodynamics simulations of red supergiant stars. III. Spectro-photocentric variability, photometric variability, and consequences on Gaia measurements

Context. It has been shown that convection in red supergiant stars gives rise to large granules causing surface inhomogeneities together with shock waves in the photosphere. The resulting motion of the photocenter (on time scales ranging from months to years) could possibly have adverse effects on the parallax determination with Gaia. Aims. We explore the impact of the granulation on the photocentric and photometric variability. We quantify these effects in order to better characterize the error possibly altering the parallax. Methods. We use 3D radiative-hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps and spectra in the Gaia G band [325-1030 nm]. Results. We provide astrometric and photometric predictions from 3D simulations of RSGs that are used to evaluate the degradation of the astrometric parameters of evolved stars derived by Gaia. We show from RHD simulations that a supergiant like Betelgeuse exhibits a photocentric noise characterised by a standard deviation of the order of 0.1 AU. The number of bright giant and supergiant stars whose Gaia parallaxes will be altered by the photocentric noise ranges from a few tens to several thousandths. The degradation of the astrometric fit due to the presence of this photocentric noise will be noticeable up to about 5 kpc for the brightest supergiants. Moreover, parallaxes of supergiants are affected by a error of the order of a few percents. We show that the photocentric noise, as predicted by the 3D simulation, does account for a substantial part of the supplementary 'cosmic noise' that affects Hipparcos measurements of Betelgeuse and Antares.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The metal-poor end of the Spite plateau. 1: Stellar parameters, metallicities and lithium abundances

We present VLT-UVES Li abundances for 28 halo dwarf stars between [Fe/H]=-2.5 and -3.5, 10 of which have [Fe/H]<-3. Four different T_eff scales have been used. Direct Infrared Flux Method (IRFM) has been used on the basis of 2MASS infrared photometry. H_alpha wings have been fitted against synthetic grids computed by means of 1D LTE atmosphere models, assuming different self-broadening theories. Finally, a grid of H_alpha profiles has been computed by means of 3D hydrodynamical atmosphere models. The Li I doublet at 670.8 nm has been used to measure A(Li) by means of 3D hydrodynamical NLTE spectral syntheses. An analytical fit of A(Li)(3D, NLTE) as a function of equivalent width, T_eff, log g, and [Fe/H] has been derived and is made available. A(Li) does not exhibit a plateau below [Fe/H]=-3. A strong positive correlation with [Fe/H] appears, not influenced by the choice of the T_eff estimator. From a linear fit, we obtain a strong slope of about 0.30 dex in A(Li) per dex in [Fe/H], significant to 2-3 sigma, and consistent among all the four T_eff estimators. A significant slope is also detected in the A(Li)--T_eff plane, driven mainly by the coolest stars in the sample which appear Li-poor. Removing such stars does not alter the behavior in the A(Li)-[Fe/H] plane. The scatter in A(Li) increases by a factor of 2 towards lower metallicities, while the plateau appears very thin above [Fe/H]=-2.8. The meltdown of the Spite plateau below [Fe/H]\sim-3 is established, but its cause is unclear. If the primordial A(Li) is the one derived from WMAP, it appears difficult to envision a single depletion phenomenon producing a thin, metallicity independent plateau above [Fe/H]=-2.8, and a highly scattered, metallicity dependent distribution below. The fact that no star below [Fe/H]=-3 lies above the plateau suggests that they formed at plateau level and underwent subsequent depletion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The solar photospheric abundance of carbon.Analysis of atomic carbon lines with the CO5BOLD solar model

The use of hydrodynamical simulations, the selection of atomic data, and the computation of deviations from local thermodynamical equilibrium for the analysis of the solar spectra have implied a downward revision of the solar metallicity. We are in the process of using the latest simulations computed with the CO5BOLD code to reassess the solar chemical composition. We determine the solar photospheric carbon abundance by using a radiation-hydrodynamical CO5BOLD model, and compute the departures from local thermodynamical equilibrium by using the Kiel code. We measure equivalent widths of atomic CI lines on high resolution, high signal-to-noise ratio solar atlases. Deviations from local thermodynamic equilibrium are computed in 1D with the Kiel code. Our recommended value for the solar carbon abundance, relies on 98 independent measurements of observed lines and is A(C)=8.50+-0.06, the quoted error is the sum of statistical and systematic error. Combined with our recent results for the solar oxygen and nitrogen abundances this implies a solar metallicity of Z=0.0154 and Z/X=0.0211. Our analysis implies a solar carbon abundance which is about 0.1 dex higher than what was found in previous analysis based on different 3D hydrodynamical computations. The difference is partly driven by our equivalent width measurements (we measure, on average, larger equivalent widths with respect to the other work based on a 3D model), in part it is likely due to the different properties of the hydrodynamical simulations and the spectrum synthesis code. The solar metallicity we obtain from the CO5BOLD analyses is in slightly better agreement with the constraints of helioseismology than the previous 3D abundance results. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Radiative hydrodynamics simulations of red supergiant stars: I. interpretation of interferometric observations

Context. It has been suggested that convection in Red Supergiant (RSG) stars gives rise to large-scale granules causing observable surface inhomogeneities. This convection is also extremely vigorous, and suspected to be one of the causes of mass-loss in RSGs. It must thus be understood in details. Evidence has been accumulated that there are asymmetries in the photospheres of RSGs, but detailedstudies of granulation are still lacking. Interferometric observations offer an exciting possibility to tackle this question, but they are still often interpreted using smooth symmetrical limb-darkened intensity distributions, or very simple spotted ad hoc models. Aims. We explore the impact of the granulation on visibility curves and closure phases using the radiative transfer code OPTIM3D. We simultaneously assess how 3D simulations of convection in RSG with CO5BOLD can be tested against these observations. Methods. We use 3D radiative-hydrodynamics (RHD) simulations of convection to compute intensity maps at various wavelengths and time, from which we derive interferometric visibility amplitudes and phases. We study their behaviour with time, position angle, and wavelength, and compare them to observations of the RSG alpha Ori Results. We provide average limb-darkening coefficients for RSGs. We detail the prospects for the detection and characterization of granulation (contrast, size) on RSGs. We demonstrate that our RHD simulations provide an excellent fit to existing interferometric observation of alpha Ori, contrary to limb darkened disks. This confirms the existence of large convective cells on the surface of Betelgeuse.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Solar abundances and 3D model atmospheres

We present solar photospheric abundances for 12 elements from optical and near-infrared spectroscopy. The abundance analysis was conducted employing 3D hydrodynamical (CO5BOLD) as well as standard 1D hydrostatic model atmospheres. We compare our results to others with emphasis on discrepancies and still lingering problems, in particular exemplified by the pivotal abundance of oxygen. We argue that the thermal structure of the lower solar photosphere is very well represented by our 3D model. We obtain an excellent match of the observed center-to-limb variation of the line-blanketed continuum intensity, also at wavelengths shortward of the Balmer jump.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

VLTI/AMBER spectro-interferometric imaging of VX Sgr's inhomogenous outer atmosphere

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

3-D hydrodynamic simulations of convection in A stars

Broadening and asymmetry of spectral lines in slowly rotating late A-type stars provide evidence for high-amplitude convective motions. The properties of turbulence observed in the A-star atmospheres are not understood theoretically and contradict results of previous numerical simulations of convection. Here we describe an ongoing effort to understand the puzzling convection signatures of A stars with the help of 3-D hydrodynamic simulations. Our approach combines realistic spectrum synthesis and non-grey hydrodynamic models computed with the CO5BOLD code. We discuss these theoretical predictions and confront them with high-resolution spectra of A stars. Our models have, for the first time, succeeded in reproducing the observed profiles of weak spectral lines without introducing fudge broadening parameters.

0 citations0 reviewsNo code linkedOpen access