David R. Alexander

Low-Temperature Opacities

Previous computations of low-temperature Rosseland and Planck mean opacities from Alexander & Ferguson areupdatedandexpanded.Thenewcomputationsincludeamorecompleteequationofstate(EOS)withmoregrain species and updated optical constants. Grains are now explicitly included in thermal equilibrium in the EOS calculation, which allows for a much wider range of grain compositions to be accurately included than was previously the case. The inclusion of high-temperature condensates such as Al2O3 and CaTiO3 significantly affects the total opacityoveranarrowrangeoftemperaturesbeforetheappearanceofthefirstsilicategrains.Thenewopacitytables are tabulated for temperatures ranging from 30,000 to 500 K with gas densities from 10 � 4 to 10 � 19 gc m � 3 .C omparisons with previous Rosseland mean opacity calculations are discussed. At high temperatures, the agreement with OPAL and Opacity Project is quite good. Comparisons at lower temperatures are more divergent as a result of differences in molecular and grain physics included in different calculations. The computation of Planck mean opacities performed with the opacity sampling method is shown to require a very large number of opacity sampling wavelength points; previously published results obtained with fewer wavelength points are shown to be significantly in error. Methods for requesting or obtaining the new tables are provided. Subject heading gs: atomic data — equation of state — methods: numerical — molecular data

THE LIMITING EFFECTS OF DUST IN BROWN DWARF MODEL ATMOSPHERES

We present opacity sampling model atmospheres, synthetic spectra, and colors for brown dwarfs and very low mass stars in the following two limiting cases of dust grain formation: (1) Inefficient gravitational settling (i.e., the dust is distributed according to the chemical equilibrium predictions) and (2) efficient gravitational settling (i.e., the dust forms and depletes refractory elements from the gas, but their opacity does not affect the thermal structure). The models include the formation of over 600 gas-phase species and 1000 liquids and crystals and the opacities of 30 different types of grains including corundum (Al2O3), the magnesium aluminum spinel MgAl2O4, iron, enstatite (MgSiO3), forsterite (Mg2SiO4), amorphous carbon, SiC, and a number of calcium silicates. The models extend from the beginning of the grain formation regime well into the condensation regime of water ice (Teff = 3000-100 K) and encompass the range of log g = 2.5-6.0 at solar metallicity. We find that silicate dust grains can form abundantly in the outer atmospheric layers of red and brown dwarfs with a spectral type later than M8. The greenhouse effects of dust opacities provide a natural explanation for the peculiarly red spectroscopic distribution of the latest M dwarfs and young brown dwarfs. The grainless (cond) models, on the other hand, correspond closely to methane brown dwarfs such as Gliese 229B. We also discover that the λλ5891, 5897 Na I D and λλ7687, 7701 K I resonance doublets play a critical role in T dwarfs, in which their red wings define the pseudocontinuum from the I to the Z bandpass.

The NEXTGEN Model Atmosphere Grid. II. Spherically Symmetric Model Atmospheres for Giant Stars with Effective Temperatures between 3000 and 6800 K

We present the extension of our NextGen model atmosphere grid to the regime of giant stars. The input physics of the models presented here is nearly identical to that of the NextGen dwarf atmosphere models; however, spherical geometry is used self-consistently in the model calculations (including the radiative transfer). We revisit the discussion of the effects of spherical geometry on the structure of the atmospheres and the emitted spectra and discuss the results of non-LTE calculations for a few selected models.

Models for Old, Metal-poor Stars with Enhanced α-Element Abundances. I. Evolutionary Tracks and ZAHB Loci; Observational Constraints

Stellar evolutionary tracks have been computed for 17 [Fe/H] values from -2.31 to -0.30 assuming, in each case, [?/Fe] = 0.0, 0.3, and 0.6. The helium abundance was assumed to vary from Y = 0.2352 at [Fe/H] = -2.31 to Y = 0.2550 at [Fe/H] = -0.30 and held constant for the different choices of [?/Fe] at a fixed iron content. Masses in the range 0.5 ? ? ? 1.0, in 0.1 ? steps, were generally considered, though sequences for higher mass values were computed, as necessary, to ensure that isochrones as young as 8 Gyr could be generated for each grid. All of the stellar models are based on an equation of state that treats nonideal effects, the latest nuclear reaction and neutrino cooling rates, and opacities that were computed specifically for the adopted chemical mixtures. The tracks were extended to the tip of the giant branch or to an age of 30 Gyr, whichever came first, and zero-age horizontal-branch (ZAHB) loci were constructed using the helium core masses and chemical profiles from appropriate red giant precursors. Selected models have been compared with those computed by A. V. Sweigart, for the same masses and chemical compositions, to demonstrate that the results obtained from two entirely independent stellar evolution codes agree well with one another when very similar input physics is assumed. In the case of extremely metal-deficient stars, an enhancement in the abundance of the ?-elements causes a single, fairly significant bump in the opacity at a temperature just above 106 K, which is caused by absorption processes involving the K shell of oxygen. This peak becomes steadily more pronounced as the overall metallicity increases and a second bump, arising from the L edges of Ne, Mg, and Si, eventually appears near log T = 5.6. As far as the tracks and isochrones are concerned, we find that, as already reported by others, it is possible to mimic the computations for [?/Fe] > 0 remarkably well by those for scaled-solar mixes simply by requiring the total mass-fraction abundance of the heavy elements, Z, to be the same. However, this result holds only for metallicities significantly less than solar. Above [Fe/H] -0.8, tracks and isochrones for enhanced ?-element mixtures begin to have systematically hotter/bluer turnoffs and red giant branches than those for scaled-solar mixtures of the heavy elements. Also addressed is the extent to which our models satisfy the constraints posed by the local subdwarfs, the distances of which are based on Hipparcos parallax measurements. Our analysis suggests that the predicted metallicity dependence of the location of the lower main sequence on the C-M diagram is in good agreement with the observed dependence. In fact, we do not find any compelling evidence from the local Population II calibrators that the colors of our models require significant adjustments. In further support of our calculations, we find that, both in zero point and slope, the computed giant branches on the (Mbol, log Teff)-plane agree well with those inferred for globular clusters from observations in the infrared. Moreover, our ZAHB models have luminosities that are just outside the 1 ? error bars of the mean MVs inferred for RR Lyrae stars from Baade-Wesselink, statistical parallax, and trigonometric parallax studies. Lower helium contents or higher ?-element abundances or an increase in the conductive opacities are among the possible ways of reducing the differences that remain. To facilitate comparisons with observations, the tracks/ZAHBs are provided with predicted BV(RI)C photometry.

The YY Isochrones for alpha-element Enhanced Mixtures

We present a new set of isochrones in which the effect of the alpha-element enhancement is fully incorporated. These isochrones are an extension of the already published set of YY Isochrones (Yi et al. 2001: Paper 1), constructed for the scaled-solar mixture. As in Paper 1, helium diffusion and convective core overshoot have been taken into account.The range of chemical compositions covered is 0.00001 < Z < 0.08. The models were evolved from the pre-main-sequence stellar birthline to the onset of helium burning in the core. The age range of the full isochrone set is 0.1 -- 20 Gyr, while younger isochrones of age 1 -- 80 Myr are also presented up to the main-sequence turn-off. Combining this set with that of Paper 1 for scaled-solar mixture isochrones, we provide a consistent set of isochrones which can be used to investigate populations of any value of alpha-enhancement. We confirm the earlier results of Paper 1 that inclusion of alpha-enhancement effects further reduces the age estimates of globular clusters by approximately 8 percent if [alpha/Fe]=+0.3. It is important to note the metallicity dependence of the change in age estimates (larger age reductions in lower metallicities). This reduces the age gap between the oldest metal-rich and metal-poor Galactic stellar populations and between the halo and the disk populations. The isochrone tables, together with interpolation routines have been made available via internet; this http URL this http URL this http URL

The Y2 Isochrones for α-Element Enhanced Mixtures

We present a new set of isochrones in which the effect of the alpha-element enhancement is fully incorporated. These isochrones are an extension of the already published set of YY Isochrones (Yi et al. 2001: Paper 1), constructed for the scaled-solar mixture. As in Paper 1, helium diffusion and convective core overshoot have been taken into account.The range of chemical compositions covered is 0.00001 < Z < 0.08. The models were evolved from the pre-main-sequence stellar birthline to the onset of helium burning in the core. The age range of the full isochrone set is 0.1 -- 20 Gyr, while younger isochrones of age 1 -- 80 Myr are also presented up to the main-sequence turn-off. Combining this set with that of Paper 1 for scaled-solar mixture isochrones, we provide a consistent set of isochrones which can be used to investigate populations of any value of alpha-enhancement. We confirm the earlier results of Paper 1 that inclusion of alpha-enhancement effects further reduces the age estimates of globular clusters by approximately 8 percent if [alpha/Fe]=+0.3. It is important to note the metallicity dependence of the change in age estimates (larger age reductions in lower metallicities). This reduces the age gap between the oldest metal-rich and metal-poor Galactic stellar populations and between the halo and the disk populations. The isochrone tables, together with interpolation routines have been made available via internet; this http URL this http URL this http URL

An opacity-sampled treatment of water vapor

Although the bands of H{sub 2}O are strong in the spectra of cool stars and calculations have repeatedly demonstrated their significance as opacity sources, only approximate opacities are currently available, due both to the difficulty of accounting for the millions of lines involved and to the inadequacy of laboratory and theoretical data. To overcome these obstacles, a new treatment is presented, based upon a statistical representation of the water vapor spectrum derived from available laboratory data. This statistical spectrum of water vapor employs an exponential distribution of line strengths and random positions of lines whose overall properties are forced to reproduce the mean opacities observed in the laboratory. The resultant data set is then treated by the opacity-sampling method exactly as are all other lines, both molecular and atomic. Significant differences are found between the results of this improved treatment and the results obtained with previous treatments of water-vapor opacity. 36 refs.

Non-Local Thermodynamic Equilibrium Effects of Ti I in M Dwarfs and Giants

We present detailed NLTE Ti I calculations in model atmospheres of cool dwarf and giant stars. A fully self-consistent NLTE treatment for a Ti I model atom with 395 levels and 5279 primary bound-bound transitions is included, and we discuss the implication of departures from LTE in this atom for the strengths of Ti I lines and TiO molecular bands in cool star spectra. We show that in the atmospheric parameter range investigated, LTE is a poor approximation to Ti I line formation, as expected from the low collisional rates in cool stars. The secondary effects of Ti I overionization on the TiO number density and the TiO molecular opacities, however, are found to be negligible in the molecular line-forming region for the relatively small parameter range studied in this paper.

Model atmospheres for K and M giants

To aid in the interpretations of observations of K and M giant stars, model atmospheres have been calculated, covering the range in effective temperature between 3000 and 4000 K and in surface gravity-log g between 0.0 and 2.0. These models are computed under the usual assumptions of hydrostatic equilibrium, constancy of total (radiative plus convective) flux, and LTE in plane-parallel, horizontally homogeneous geometry. A significant development in these models is the incorporation of an opacity-sampled treatment of absorption by water vapor, which becomes dominant at lower temperatures. The thermal structure and scaling properties of the models are discussed, as are the effects of certain composition changes expected in the course of stellar evolution. 51 refs.

Outer layers of a carbon star: The view from the Hubble Space Telescope

To advance our understanding of the relationship between stellar chromospheres and mass loss, which is a common property of carbon stars and other asymptotic giant branch stars, we have obtained ultraviolet spectra of the nearby N-type carbon star UU Aur using the Hubble Space Telescope (HST). In this paper we describe the HST observations, identify spectral features in both absorption and emission, and attempt to infer the velocity field in the chromosphere, upper troposphere, and circumstellar envelope from spectral line shifts. A mechanism for producing fluoresced emission to explain a previously unobserved emission line is proposed. Some related ground-based observations are also described.