The ability to observe extrasolar planets transiting their stars has profoundly changed our understanding of these planetary systems. However, these measurements depend on how well we understand the properties of the host star, such as radius, luminosity and limb darkening. Traditionally, limb darkening is treated as a parameterization in the analysis, but these simple parameterizations are not accurate representations of actual center-to-limb intensity variations (CLIV) to the precision needed for interpreting these transit observations. This effect leads to systematic errors for the measured planetary radii and corresponding measured spectral features. We compute synthetic planetary transits using model stellar atmosphere CLIV and corresponding best-fit limb-darkening laws for a grid spherically symmetric model stellar atmospheres. From these light curves we measure the differences in flux as a function of the star's effective temperature, gravity, mass, and the inclination of the planet's orbit.

The transit method, employed by MOST, \emph{Kepler}, and various ground-based surveys has enabled the characterization of extrasolar planets to unprecedented precision. These results are precise enough to begin to measure planet atmosphere composition, planetary oblateness, star spots, and other phenomena at the level of a few hundred parts-per-million. However, these results depend on our understanding of stellar limb darkening, that is, the intensity distribution across the stellar disk that is sequentially blocked as the planet transits. Typically, stellar limb darkening is assumed to be a simple parameterization with two coefficients that are derived from stellar atmosphere models or fit directly. In this work, we revisit this assumption and compute synthetic planetary transit light curves directly from model stellar atmosphere center-to-limb intensity variations (CLIV) using the plane-parallel \textsc{Atlas} and spherically symmetric \textsc{SAtlas} codes. We compare these light curves to those constructed using best-fit limb-darkening parameterizations. We find that adopting parametric stellar limb-darkening laws lead to systematic differences from the more geometrically realistic model stellar atmosphere CLIV of about 50 -- 100 ppm at the transit center and up to 300 ppm at ingress/egress. While these errors are small they are systematic, and appear to limit the precision necessary to measure secondary effects. Our results may also have a significant impact on transit spectra.

One of the great challenges in understanding stars is measuring their masses. The best methods for measuring stellar masses include binary interaction, asteroseismology and stellar evolution models, but these methods are not ideal for red giant and supergiant stars. In this work, we propose a novel method for inferring stellar masses of evolved red giant and supergiant stars using interferometric and spectrophotometric observations combined with spherical model stellar atmospheres to measure what we call the stellar mass index, defined as the ratio between the stellar radius and mass. The method is based on the correlation between different measurements of angular diameter, used as a proxy for atmospheric extension, and fundamental stellar parameters. For a given star, spectrophotometry measures the Rosseland angular diameter while interferometric observations generally probe a larger limb-darkened angular diameter. The ratio of these two angular diameters is proportional to the relative extension of the stellar atmosphere, which is strongly correlated to the star's effective temperature, radius and mass. We show that these correlations are strong and can lead to precise measurements of stellar masses.

Limb darkening is a fundamental ingredient for interpreting observations of planetary transits, eclipsing binaries, optical/infrared interferometry and microlensing events. However, this modeling traditionally represents limb darkening by a simple law having one or two coefficients that have been derived from plane-parallel model stellar atmospheres, which has been done by many researchers. More recently, researchers have gone beyond plane-parallel models and considered other geometries. We previously studied the limb-darkening coefficients from spherically symmetric and plane-parallel model stellar atmospheres for cool giant and supergiant stars, and in this investigation we apply the same techniques to FGK dwarf stars. We present limb-darkening coefficients, gravity-darkening coefficients and interferometric angular diameter corrections from Atlas and SAtlas model stellar atmospheres. We find that sphericity is important even for dwarf model atmospheres, leading to significant differences in the predicted coefficients.

Model stellar atmospheres are fundamental tools for understanding stellar observations from interferometry, microlensing, eclipsing binaries and planetary transits. However, the calculations also include assumptions, such as the geometry of the model. We use intensity profiles computed for both plane-parallel and spherically symmetric model atmospheres to determine fitting coefficients in the BVRIHK, CoRot and Kepler wavebands for limb darkening using several different fitting laws, for gravity-darkening and for interferometric angular diameter corrections. Comparing predicted variables for each geometry, we find that the spherically symmetric model geometry leads to different predictions for surface gravities log g < 3. In particular, the most commonly used limb-darkening laws produce poor fits to the intensity profiles of spherically symmetric model atmospheres, which indicates the need for more sophisticated laws. Angular diameter corrections for spherically symmetric models range from 0.67 to 1, compared to the much smaller range from 0.95 to 1 for plane-parallel models.

Mass is the most important stellar parameter, but it is not directly observable for a single star. Spherical model stellar atmospheres are explicitly characterized by their luminosity ($L_\star$), mass ($M_\star$) and radius ($R_\star$), and observations can now determine directly $L_\star$ and $R_\star$. We computed spherical model atmospheres for red giants and for red supergiants holding $L_\star$ and $R_\star$ constant at characteristic values for each type of star but varying $M_\star$, and we searched the predicted flux spectra and surface-brightness distributions for features that changed with mass. For both stellar classes we found similar signatures of the star's mass in both the surface-brightness distribution and the flux spectrum. The spectral features have been use previously to determine $\log_{10} (g)$, and now that the luminosity and radius of a non-binary red giant or red supergiant can be observed, spherical model stellar atmospheres can be used to determine the star's mass from currently achievable spectroscopy. The surface-brightness variations with mass are slightly smaller than can be resolved by current stellar imaging, but they offer the advantage of being less sensitive to the detailed chemical composition of the atmosphere.

Context. Limb darkening is an important tool for understanding stellar atmospheres, but most observations measuring limb darkening assume various parameterizations that yield no significant information about the structure of stellar atmospheres. Aims. We use a specific limb-darkening relation to study how the best-fit coefficients relate to fundamental stellar parameters from spherically symmetric model stellar atmospheres. Methods. Using a grid of spherically symmetric Atlas model atmospheres, we compute limb-darkening coefficients, and develop a novel method to predict fundamental stellar parameters. Results. We find our proposed method predicts the mass of stellar atmosphere models given only the radius and limb-darkening coefficients, suggesting that microlensing, interferometric, transit and eclipse observations can constrain stellar masses. Conclusions. This novel method demonstrates that limb-darkening parameterizations contain important information about the structure of stellar atmospheres, with the potential to be a valuable tool for measuring stellar masses.

One of the largest uncertainties for using the Baade-Wesselink method to measure Cepheid distances is the value of the projection factor (p-factor). However, p-factors measured using the IRSB technique and from hydrodynamic models disagree. In this work, we compute spherically-symmetric static model stellar atmospheres and predict a period p-factor relation.

Stellar limb-darkening is an important tool for constraining the properties of a stellar atmosphere. We present a novel method for relating the fundamental stellar parameters mass and radius to limb-darkening laws using grids of spherical model stellar atmospheres. This method is applied to interferometric observations of the red supergiant Betelgeuse, where an unique measure of the stellar mass is determined.

(Abridged) Context. Stellar limb darkening, I({\mu} = cos{\theta}), is an important constraint for microlensing, eclipsing binary, planetary transit, and interferometric observations, but is generally treated as a parameterized curve, such as a linear-plus-square-root law. Many analyses assume limb-darkening coefficients computed from model stellar atmospheres. However, previous studies, using I({\mu}) from plane- parallel models, have found that fits to the flux-normalized curves pass through a fixed point, a common {\mu} location on the stellar disk, for all values of T eff, log g and wavelength. Aims. We study this fixed {\mu}-point to determine if it is a property of the model stellar atmospheres or a property of the limb-darkening laws. Furthermore, we use this limb-darkening law as a tool to probe properties of stellar atmospheres for comparison to limb- darkening observations. Methods. Intensities computed with plane-parallel and spherically-symmetric Atlas models (characterized by the three fundamental parameters L\star, M\star and R\star) are used to reexamine the existence of the fixed {\mu}-point for the parametrized curves. Results. We find that the intensities from our spherical models do not have a fixed point, although the curves do have a minimum spread at a {\mu}-value similar to the parametrized curves. We also find that the parametrized curves have two fixed points, {\mu}1 and {\mu}2, although {\mu}2 is so close to the edge of the disk that it is missed using plane-parallel atmospheres. We also find that the spherically- symmetric models appear to agree better with published microlensing observations relative to plane-parallel models.

We present a novel method for measuring the masses of evolved stars from their limb-darkening observations parameterized as a linear-plus-square-root function with two coefficients. The coefficients of the law are related to integrated moments of the intensity and carry information about the extension of the stellar atmosphere, which is correlated to the ratio of the stellar radius to mass, R/M. Here, we show why the limb-darkening law is related to the R/M, and apply this result to limb-darkening observations of the microlensing event EROS-BLG-2000-5.

In this article, we test the hypothesis that Cepheids have infrared excesses due to mass loss. We fit a model using the mass-loss rate and the stellar radius as free parameters to optical observations from the OGLE-III survey and infrared observations from the 2MASS and SAGE data sets. The sample of Cepheids have predicted minimum mass-loss rates ranging from zero to $10^{-8}M_\odot$ $yr^{-1}$, where the rates depend on the chosen dust properties. We use the predicted radii to compute the Period-Radius relation for LMC Cepheids, and to estimate the uncertainty caused by the presence of infrared excess for determining angular diameters with the infrared surface brightness technique. Finally, we calculate the linear and non-linear Period-Luminosity (P-L) relations for the LMC Cepheids at VIJHK + IRAC wavelengths and we find that the P-L relations are consistent with being non-linear at infrared wavelengths, contrary to previous results.

Recent observations of Cepheids using infrared interferometry and Spitzer photometry have detected the presence of circumstellar envelopes (CSE) of dust and it has been hypothesized that the CSE's are due to dust forming in a Cepheid wind. Here we use a modified Castor, Abbott & Klein formalism to produce a Cepheid wind, and this is used to estimate the contribution of mass loss to the Cepheid mass discrepancy Furthermore, we test the OGLE-III Classical Cepheids using the IR fluxes from the SAGE survey to determine if Large Magellanic Cloud Cepheids have CSE's. It is found that IR excess is a common phenomenon for LMC Cepheids and that the resulting mass-loss rates can explain at least a fraction of the Cepheid mass discrepancy, depending on the assumed dust-to-gas ratio in the wind.

It has been claimed that Period-Luminosity relations derived from infrared observations of Large Magellanic Cloud (LMC) Cepheids are less dependent on the metallicity of the Cepheids. In this work, infrared observations of LMC Cepheids from the SAGE survey are combined with OGLE II optical observation to model and predict mass-loss rates. The mass-loss rates are fit to the data and are predicted to range from about 10^{-12} to 10^{-7}M_\odot/yr; however, the rates depend on the assumed value of the dust-to-gas ratio. By comparing the relations derived from observations to the relations derived from predicted infrared stellar luminosities from the mass-loss model, it is shown mass loss affects the structure and scatter of the infrared Period-Luminosity relation. Mass loss produces shallower slopes of the infrared relations and a lower zero point. There is also evidence for non-linearity in the predicted Period-Luminosity relations, and it is argued that mass loss produces larger infrared excess at lower periods, which affects the slope and zero point, making the PL relations more linear in the wavelength range of 3.6 to 5.8 micron. Because the dust-to-gas ratio is metallicity dependent and mass loss may have a metallicity dependence, infrared Period-Luminosity relations have additional uncertainty due to metallicity.

It has been observed that Cepheids in the Magellanic Clouds have lower masses for the same luminosity than those in the Milky Way. The model, from Neilson & Lester (2008), of pulsation-driven mass loss for Cepheids is applied to theoretical models of Cepheids with metallicity consistent with the Milky Way and Large and Small Magellanic Clouds. The mass-loss model is analyzed using the metallicity correction of the Period-Luminosity relation to compare the ratio of mass loss of Cepheids with lower metallicity to that of Cepheids with solar metallicity. It is determined that mass loss may be larger for the lower metallicity Cepheids, counterintuitive to radiative driving estimates. Also the mass-loss rates of theoretical Cepheid models are found to be up to $5\times 10^{-9}$ for Galactic Cepheids, $5\times 10^{-8}$ for Large Magellanic Cloud Cepheids, and $2\times 10^{-7}M_\odot /yr$ for Small Magellanic Cloud Cepheids. It is argued that mass loss increases as metallicity decreases for Cepheids with periods less than 20 days and that mass loss decreases for longer periods. Assuming dust forms in the wind of a Cepheid at some distance, the infrared excess of the models is computed, finding the infrared brightness is approximately a magnitude larger due to mass loss. The infrared magnitudes are compared to recently published Period-Luminosity relations as a test of our predictions.

Context: Optical interferometry is a powerful tool for observing the intensity structure and angular diameter of stars. When combined with spectroscopy and/or spectrophotometry, interferometry provides a powerful constraint for model stellar atmospheres. Aims: The purpose of this work is to test the robustness of the spherically symmetric version of the Atlas stellar atmosphere program, SAtlas, using interferometric and spectrophotometric observations. Methods: Cubes (three dimensional grids) of model stellar atmospheres, with dimensions of luminosity, mass, and radius, are computed to fit observations for three evolved giant stars, \psi Phoenicis, \gamma Sagittae, and \alpha Ceti. The best-fit parameters are compared with previous results. Results: The best-fit angular diameters and values of \chi^2 are consistent with predictions using Phoenix and plane-parallel Atlas models. The predicted effective temperatures, using SAtlas, are about 100 to 200 K lower, and the predicted luminosities are also lower due to the differences in effective temperatures. Conclusions: It is shown that the SAtlas program is a robust tool for computing models of extended stellar atmospheres that are consistent with observations. The best-fit parameters are consistent with predictions using Phoenix models, and the fit to the interferometric data for \psi Phe differs slightly, although both agree within the uncertainty of the interferometric observations.

Context: The current stellar atmosphere programs still cannot match some fundamental observations of the brightest stars, and with new techniques, such as optical interferometry, providing new data for these stars, additional development of stellar atmosphere codes is required. Aims: To modify the open-source model atmosphere program Atlas to treat spherical geometry, creating a test-bed stellar atmosphere code for stars with extended atmospheres. Methods: The plane-parallel Atlas has been changed by introducing the necessary spherical modifications in the pressure structure, in the radiative transfer and in the temperature correction. Results: Several test models show that the spherical program matches the plane-parallel models in the high surface gravity regime, and matches spherical models computed by Phoenix and by MARCS in the low gravity case.

One of the current opportunities for stellar atmospheric modeling is the interpretation of optical interferometric data of stars. Starting from the robust, open source ATLAS atmospheric code (Kurucz, 1979), we have developed a spherically symmetric code, SATLAS, as a new option for modeling stellar atmospheres of low gravity stars. The SATLAS code is tested against both interferometric observations of M giants by Wittkowski and collaborators, and spherically symmetric M giant NextGen models from the PHOENIX code. The SATLAS models predict interferometric visibilities that agree with the observed visibilities and with predicted visibilities, and the SATLAS atmospheric structures also agree with those from spherical PHOENIX models, with just small differences in temperature and pressure at large depths in the atmospheres.

The purpose of this preliminary work is to determine if Large Magellanic Cloud (LMC) Cepheids have stellar winds. If a Cepheid undergoes mass loss then at some distance from the star, a fraction of the gas becomes dust, which causes an infrared excess. Mass loss is tested using OGLE II optical observations and SAGE infrared observations for a sample of 488 Cepheids. The resultant mass-loss rates range from $10^{-12}$ to $10^{-7}$ $M_\odot/yr$. Using the mass--loss model we compute infrared stellar luminosities for the sample of Cepheids and compare predicted infrared PL relations with observed relations. The predicted relations not only vary from the observed relations, implying mass loss plays a significant role, but also show evidence for non-linearity. It is determined that mass loss is important for LMC Cepheids.

An analytical derivation is presented for computing mass-loss rates of Cepheids by using the method of Castor, Abbott, & Klein (1975) modified to include a term for momentum input from pulsation and shocks generated in the atmosphere. Using this derivation, mass-loss rates of Cepheids are determined as a function of stellar parameters. When applied to a set of known Cepheids, the calculated mass-loss rates range from 10^{-10} to 10^{-7}M_{Sun}/yr, larger than if the winds were driven by radiation alone. Infrared excesses based on the predicted mass-loss rates are compared to observations from optical interferometry and IRAS, and predictions are made for Spitzer observations. The mass-loss rates are consistent with the observations, within the uncertainties of each. The rate of period change of Cepheids is discussed and shown to relate to mass loss, albeit the dependence is very weak. There is also a correlation between the large mass-loss rates and the Cepheids with slowest absolute rate of period change due to evolution through the instability strip. The enhanced mass loss helps illuminate the issue of infrared excess and the mass discrepancy found in Cepheids.

Long-baseline interferometric observations obtained with the Navy Prototype Optical Interferometer of the H-alpha-emitting envelopes of the Be stars eta Tauri and beta Canis Minoris are presented. For compatibility with the previously published interferometric results in the literature of other Be stars, circularly symmetric and elliptical Gaussian models were fitted to the calibrated H-alpha observations. The models are sufficient in characterizing the angular distribution of the H-alpha-emitting circumstellar material associated with these Be stars. To study the correlations between the various model parameters and the stellar properties, the model parameters for eta Tau and beta CMi were combined with data for other Be stars from the literature. After accounting for the different distances to the sources and stellar continuum flux levels, it was possible to study the relationship between the net H-alpha emission and the physical extent of the H-alpha-emitting circumstellar region. A clear dependence of the net H-alpha emission on the linear size of the emitting region is demonstrated and these results are consistent with an optically thick line emission that is directly proportional to the effective area of the emitting disk. Within the small sample of stars considered in this analysis, no clear dependence on the spectral type or stellar rotation is found, although the results do suggest that hotter stars might have more extended H-alpha-emitting regions.

Synthetic spectra of Arcturus computed with the latest theoretical stellar atmospheres disagree strongly with both the spectral energy distribution and with the spectroscopic data for $\lambda < 4000$ \AA , even though the stellar parameters are well constrained. We find that the discrepancy can be removed by adding a continuous absorption opacity that is approximately equal to that which is already included in the {\sc atlas9} code. Spectroscopy of other K2 III stars shows that the spectrum of Arcturus in not anomalous. Therefore, missing opacity may be a common feature of the models for stars in this temperature and luminosity range.