Model Atmospheres For Stars Research Articles

The BT-Settl Model Atmospheres for Stars, Brown Dwarfs and Planets

AbstractWe present a grid of stellar and substellar atmosphere models covering the range from solar-mass stars to the latest-type T and Y dwarfs with a single setup. For the first time our synthetic spectra and photometry reproduce the formation of clouds and in particular their clearing at the L/T transition. The BT-Settl models also naturally explain the dustier infrared properties of planets as an effect of low surface gravity.

Theoretical backgrounds for interpretation of spectroscopic observations

The advantage of analysis of observed data using their fit with theoretical model of the directly observed quantities is shown as well as the need for simultaneous solution of all available data. Some particular problems of disentangling of stellar spectra and model atmospheres of component stars of multiple systems are discussed.

Three-dimensional hydrodynamical CO5BOLD model atmospheres of red giant stars

We investigate the role of convection in the formation of atomic and molecular lines in the atmosphere of a red giant star. For this purpose we study the formation properties of spectral lines that belong to a number of astrophysically important tracer elements, including neutral and singly ionized atoms, and molecules. We focus our investigation on a prototypical red giant located close to the red giant branch (RGB) tip (Teff=3660K, logg=1.0, [M/H]=0.0). We used two types of model atmospheres, 3D hydrodynamical and classical 1D, calculated with the CO5BOLD and LHD stellar atmosphere codes, respectively. Both codes share the same atmospheric parameters, chemical composition, equation of state, and opacities, which allowed us to make a strictly differential comparison between the line formation properties predicted in 3D and 1D. The influence of convection on the spectral line formation was assessed with the aid of 3D-1D abundance corrections, which measure the difference between the abundances of chemical species derived with the 3D hydrodynamical and 1D classical model atmospheres. We find that convection plays a significant role in the spectral line formation in this particular red giant. The derived 3D-1D abundance corrections rarely exceed \pm0.1 dex when lines of neutral atoms and molecules are considered, which is in line with the previous findings for solar-metallicity red giants located on the lower RGB. The situation is different with lines that belong to ionized atoms, or to neutral atoms with high ionization potential. In both cases, the corrections for high-excitation lines (\chi>8 eV) may amount to \Delta_3D-1D ~ -0.4 dex. The 3D--1D abundance corrections generally show a significant wavelength dependence; in most cases they are smaller in the near-infrared, at 1600-2500 nm.

X-ray bursting neutron star atmosphere models using an exact relativistic kinetic equation for Compton scattering

Theoretical spectra of X-ray bursting neutron star (NS) model atmospheres are widely used to determine the basic NS parameters such as their masses and radii. We construct accurate NS atmosphere models using for the first time an exact treatment of Compton scattering via the integral relativistic kinetic equation. We also compare the results with the previous calculations based on the Kompaneets operator. We solve the radiation transfer equation together with the hydrostatic equilibrium equation accounting exactly for the radiation pressure by electron scattering. We thus construct a new set of plane-parallel atmosphere models in LTE for hot NSs. The models were computed for six chemical compositions (pure H, pure He, solar H/He mix with various heavy elements abundances Z = 1, 0.3, 0.1, and 0.01 Z_sun, and three log g = 14.0, 14.3, and 14.6. For each chemical composition and log g, we compute more than 26 model atmospheres with various luminosities relative to the Eddington luminosity L_Edd computed for the Thomson cross-section. The maximum relative luminosities L/L_Edd reach values of up to 1.1 for high gravity models. The emergent spectra of all models are redshifted and fitted by diluted blackbody spectra in the 3--20 keV energy range appropriate for the RXTE/PCA. We also compute the color correction factors f_c. The radiative acceleration g_rad in our luminous, hot-atmosphere models is significantly smaller than in corresponding models based on the Kompaneets operator, because of the Klein-Nishina reduction of the electron scattering cross-section, and therefore formally "super-Eddington" model atmospheres do exist. The differences between the new and old model atmospheres are small for L / L_Edd < 0.8. For the same g_rad / g, the new f_c are slightly larger (by approximately 1%) than the old values.

3D LTE spectral line formation with scattering in red giant stars

We investigate the effects of coherent isotropic continuum scattering on the formation of spectral lines in local thermodynamic equilibrium (LTE) using 3D hydrodynamical and 1D hydrostatic model atmospheres of red giant stars. Continuum flux levels, spectral line profiles and curves of growth for different species are compared with calculations that treat scattering as absorption. Photons may escape from deeper, hotter layers through scattering, resulting in significantly higher continuum flux levels beneath a wavelength of 5000 A. The magnitude of the effect is determined by the importance of scattering opacity with respect to absorption opacity; we observe the largest changes in continuum flux at the shortest wavelengths and lowest metallicities; intergranular lanes of 3D models are more strongly affected than granules. Continuum scattering acts to increase the profile depth of LTE lines: continua gain more brightness than line cores due to their larger thermalization depth in hotter layers. We thus observe the strongest changes in line depth for high-excitation species and ionized species, which contribute significantly to photon thermalization through their absorption opacity near the continuum optical surface. Scattering desaturates the line profiles, leading to larger abundance corrections for stronger lines, which reach -0.5 dex at 3000 A for Fe II lines in 3D with excitation potential 2 eV at [Fe/H]=-3.0. The corrections are less severe for low-excitation lines, longer wavelengths, and higher metallicity. Velocity fields increase the effects of scattering by separating emission from granules and intergranular lanes in wavelength. 1D calculations exhibit similar scattering abundance corrections for weak lines, but those for strong lines are generally smaller compared to 3D models and depend on the choice of microturbulence.

Three-dimensional surface convection simulations of metal-poor stars

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

Galactic evolution of oxygen

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.

VLTI/AMBER spectro-interferometric imaging of VX Sagittarii'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.

Spectral Evolution Models for the Next Decade

AbstractSpectral evolution models are a widely used tool for determining the stellar content of galaxies. I provide a review of the latest developments in stellar atmosphere and evolution models, with an emphasis on massive stars. In contrast to the situation for low- and intermediate-mass stars, the current main challenge for spectral synthesis models are the uncertainties and rapid revision of current stellar evolution models. Spectral libraries, in particular those drawn from theoretical model atmospheres for hot stars, are relatively mature and can complement empirical templates for larger parameter space coverage. I introduce a new ultraviolet spectral library based on theoretical radiation-hydrodynamic atmospheres for hot massive stars. Application of this library to star-forming galaxies at high redshift, i.e., Lyman-break galaxies, will provide new insights into the abundances, initial mass function and ages of stars in the very early universe.

Interferometric properties of pulsating C-rich AGB stars

We present the first theoretical study on center-to-limb variation (CLV) properties and relative radius interpretation for narrow and broad-band filters, on the basis of a set of dynamic model atmospheres of C-rich AGB stars. We computed visibility profiles and the equivalent uniform disc radii (UD-radii) in order to investigate the dependence of these quantities upon the wavelength and pulsation phase. After an accurate morphological analysis of the visibility and intensity profiles determined in narrow and broad-band filter, we fitted our visibility profiles with a UD function simulating the observational approach. UD-radii have been computed using three different fitting-methods to investigate the influence of the sampling of the visibility profile: single point, two points and least square method. The intensity and visibility profiles of models characterized by mass loss show a behaviour very different from a UD. We found that UD-radii are wavelength dependent and this dependence is stronger if mass loss is present. Strong opacity contributions from C2H2 affect all radius measurements at 3 micron and in the N-band, resulting in higher values for the UD-radii. The predicted behaviour of UD-radii versus phase is complicated in the case of models with mass loss, while the radial changes are almost sinusoidal for the models without mass loss. Compared to the M-type stars, for the C-stars no windows for measuring the pure continuum are available.

Multiwavelength study of M33's giant H ii regions NGC 588 and NGC 592

We present a detailed multiwavelength photometric study of giant H II regions NGC 592 and NGC 588 in the nearby small spiral galaxy M33. We use data taken with the Wide Field and Planetary Camera 2 (WFPC2) on board the Hubble Space Telescope (HST). We detect several massive stars in both ionizing clusters. Six Wolf–Rayet (WR) stars are known to exist within those regions and we are able to constrain their physical properties by comparing their photometry to the latest grid of model atmospheres for WR stars of the nitrogen sequence (WN subclass). We estimate the age and mass of both regions by fitting our photometry to models of integrated stellar populations.

Fundamental parameters of BE UMa revised

We have determined a complete set of parameters for the young pre-cataclysmic variable BE UMa from a comprehensive photometric and spectroscopic analysis using model atmospheres. Our precise photometry and spectroscopy were acquired with the 6-m telescope and Zeiss-1000 telescope of the Special Astrophysical Observatory and the 1.5-m Russian-Turkish telescope at a wide range of orbital phases, including times of primary eclipses.We performed a detailed identification of emission lines of ten elements. At phases of minimum brightness, the spectra reveal absorption lines and molecular bands formed in the secondary’s atmosphere, whose effective temperature was determined to be Teff(2) = 4750 ± 150 K. We have studied the radial-velocity curves of the cool star using lines of various elements. All the curves exhibit the previously predicted distortions due to reflection effects in the close binary. The derived component-mass ratio is q = 0.43 ± 0.09, and the component masses are M1 = 0.59 ± 0.07 M⊙ and M2 = 0.25 ± 0.08 M⊙. We analyzed the light curves using model atmospheres for irradiated stars; all the parameters of BE UMa were refined. We demonstrate the validity of our modeling of the binary’s spectra at phases of brightness maximum, which provides a good description of the observed intensities of most lines of heavy elements. The abundances of helium and several light elements (C, N, O, Ne, Mg) in the atmosphere of the cool star are probably higher than the solar values. We conclude that the physical characteristics of the primary are in good agreement with evolutionary tracks for planetary-nebula nuclei, and that the secondary is overluminous by a factor of 30 compared to main-sequence stars of the same mass.

A grid of MARCS model atmospheres for late-type stars

A grid of MARCS model atmospheres for late-type stars

Spurious Eccentricities of Distorted Binary Components

I discuss the effect of physical distortion on the velocities of close binary components and how we may use the resulting distortion of velocity curves to constrain some properties of binary systems, such as inclination and mass ratio. Precise new velocities for 5 Cet convincingly detect these distortions with their theoretically predicted phase dependence. We can even use such distortions of velocity curves to test Lucy's theory of convective gravity darkening. The observed distortions for TT Hya and 5 Cet require the contact components of those systems to be gravity darkened, probably somewhat more than predicted by Lucy's theory, but clearly not as much as expected for a radiative star. These results imply there is no credible evidence for eccentric orbits in binaries with contact components. I also present some speculative analyses of the observed properties of a binary encased in a nonrotating common envelope, if such an object could actually exist, and discuss how the limb darkening of some recently calculated model atmospheres for giant stars may bias my results for velocity-curve distortions, as well as other results from a wide range of analyses of binary stars.

Model Atmospheres for Irradiated Giant Stars: Implications for the Galactic Center

Irradiation of a stellar atmosphere by an external source (e.g. an AGN) changes its structure and therefore its spectrum. Using a state-of-the-art stellar atmosphere code, we calculate the infrared spectra of such irradiated and transformed stars. We show that the original spectrum of the star, which is dominated by molecular bands, changes dramatically when irradiated even by a low-luminosity AGN ($L_{\rm X} = 10^{33}$ erg s$^{-1}$), becoming dominated by atomic lines in absorption. We study the changes in the spectrum of low-mass carbon- and oxygen-rich giant stars as they are irradiated by a modest AGN, similar to the one at the Galactic center (GC). The resulting spectra are similar to those of the faintest S-cluster stars observed in the GC. The spectrum of a star irradiated by a much brighter AGN, like that powered by a tidally disrupted star, is very different from that of any star currently observed near the GC. For the first time we have discovered that the structure of the atmosphere of an irradiated giant changes dramatically and induces a double inversion layer. We show that irradiation at the current level can explain the observed trend of CO band intensities decreasing as a function of increasing proximity to Sg $A^{*}$. This may indicate that (contrary to previous claims) there is no paucity of old giants in the GC, which coexist simultaneously with young massive stars.

2MASS J05162881+2607387: A New Low‐mass Double‐lined Eclipsing Binary

We show that the star known as 2MASS J05162881+2607387 (hereafter J0516) is a double-lined eclipsing binary with nearly identical low-mass components. The spectroscopic elements derived from 18 spectra obtained with the High Resolution Spectrograph on the Hobby-Eberly Telescope during the fall of 2005 are K1 = 88.45 ± 0.48 and K2 = 90.43 ± 0.60 km s-1, resulting in a mass ratio of q = K1/K2 = 0.978 ± 0.018 and minimum masses of M1 i = 0.775 ± 0.016 M☉ and M2 i = 0.759 ± 0.012 M☉, respectively. We have extensive differential photometry of J0516 obtained over several nights between 2004 January and March (epoch 1) and between 2004 October and 2005 January plus 2006 January (epoch 2) using the 1 m telescope at the Mount Laguna Observatory. The source was roughly 0.1 mag brighter in all three bandpasses during epoch 1 when compared to epoch 2. Also, phased light curves from epoch 1 show considerable out-of-eclipse variability, presumably due to bright spots on one or both stars. In contrast, the phased light curves from epoch 2 show little out-of-eclipse variability. The light curves from epoch 2 and the radial velocity curves were analyzed using our ELC code with updated model atmospheres for low-mass stars. We find the following: M1 = 0.787 ± 0.012 M☉, R1 = 0.788 ± 0.015 R☉, M2 = 0.770 ± 0.009 M☉, and R2 = 0.817 ± 0.010 R☉. The stars in J0516 have radii that are significantly larger than model predictions for their masses, similar to what is seen in a handful of other well-studied low-mass double-lined eclipsing binaries. We compiled all recent mass and radius determinations from low-mass binaries and determine an empirical mass-radius relation of the form R(R☉) = 0.0324 + 0.9343M(M☉) + 0.0374M2(M☉).

Geometry of giant star model atmospheres: a consistency test

We investigate the effect of a geometric inconsistency in the calculation of synthetic spectra of giant stars. Spectra computed with model atmospheres calculated in spherical geometry while using the plane-parallel approximation for line formation calculations (s_p), as well as the fully plane-parallel case (p_p), are compared to the consistently spherical case (s_s). We present abundance differences for solar metallicity models with Teff ranging from 4000 to 6500 K and logg from 0.5 to 3.0 [cgs]. The effects are smaller for s_p calculations (-0.1 dex in the worst case) than for the p_p case (up to +0.35 dex for minority species and at most -0.04 dex for majority species), both with respect to the s_s case. In the s_p case the differences increase slightly with temperature, while in the p_p case they show a more complex behaviour. In both cases the effects decrease with increasing logg and increase with equivalent width. Thus, within the parameter range of F, G and K giants, consistency seems to be less important than using a spherical model atmosphere. The abundance differences due to sphericity effects presented here can be used for error estimation in abundance studies relying on plane-parallel modelling.

Atmospheric NLTE-models for the spectroscopic analysis of blue stars with winds

We present new or improved methods for calculating NLTE, line-blanketed model atmospheres for hot stars with winds (spectral types A to O), with particular emphasis on a fast performance. These methods have been implemented into a previous, more simple version of the model atmosphere code FASTWIND (Santolaya-Rey et al.1997) and allow to spectroscopically analyze rather large samples of massive stars in a reasonable time-scale, using state-of-the-art physics. We describe our (partly approximate) approach to solve the equations of statistical equilibrium for those elements which are primarily responsible for line-blocking and blanketing, as well as an approximate treatment of the line-blocking itself, which is based on a simple statistical approach using suitable means for line opacities and emissivities. Furthermore, we comment on our implementation of a consistent temperature structure. In the second part, we concentrate on a detailed comparison with results from those two codes which have been used in alternative spectroscopical investigations, namely CMFGEN (Hillier & Miller 1998) and WM-Basic (Pauldrach et al. 2001). All three codes predict almost identical temperature structures and fluxes for lambda > 400 A, whereas at lower wavelengths a number of discrepancies are found. Optical H/He lines as synthesized by FASTWIND are compared with results from CMFGEN, obtaining a remarkable coincidence, except for the HeI singlets in the temperature range between 36,000 to 41,000 K for dwarfs and between 31,000 to 35,000 K for supergiants, where CMFGEN predicts much weaker lines. Consequences due to these discrepancies are discussed.

Hydrodynamic model atmospheres for WR stars

We present the first non-LTE atmosphere models for WR stars that incorporate a self-consistent solution of the hydrodynamic equations. The models take iron-group line-blanketing and clumping into account, and compute the hydrodynamic structure of a radiatively driven wind consistently with the non-LTE radiation transport in the co-moving frame. We construct a self-consistent wind model that reproduces all observed properties of an early-type WC star (WC5). We find that the WR-type mass-loss is initiated at high optical depth by the so-called “Hot Iron Bump” opacities (–xvi). The acceleration of the outer wind regions is due to iron-group ions of lower excitation in combination with C and O. Consequently, the wind structure shows two acceleration regions, one close to the hydrostatic wind base in the optically thick part of the atmosphere, and another farther out in the wind. In addition to the radiative acceleration, the “Iron Bump” opacities are responsible for an intense heating of deep atmospheric layers. We find that the observed narrow emission lines in the optical spectra of WC stars originate from this region. From their dependence on the clumping factor we gain important information about the location where the density inhomogeneities in WR-winds start to develop.

Grids of model spectra for WN stars, ready for use

Grids of model atmospheres for Wolf-Rayet stars of the nitrogen sequence (WN subclass) are presented. The calculations account for the expansion of the atmosphere, non-LTE, clumping, and line blanketing from iron-group elements. Observed spectra of single Galactic WN stars can in general be reproduced consistently by this generation of models. The parameters of the presented model grids cover the whole relevant range of stellar temperatures and mass-loss rates. We point out that there is a degeneracy of parameters for very thick winds; their spectra tend to depend only on the ratio $L/{dot M}^{4/3}$. Abundances of the calculated grids are for Galactic WN stars without hydrogen and with 20% hydrogen (by mass), respectively. Model spectra and fluxes are available via internet (http://www.astro.physik.uni- potsdam.de/PoWR.html).