Classical Cepheids are powerful probes of both stellar evolution and near-field cosmology thanks to their high luminosities, pulsations, and that they follow the Leavitt (Period-Luminosity) Law. However, there still exist a number of questions regarding their evolution, such as the role of rotation, convective core overshooting and winds. ln particular, how do these processes impact Cepheid evolution and the predicted fundamental properties such as stellar mass. In this work, we compare a sample of period change that are real-time observations of stellar evolution with new evolution models to test the impact of these first two processes. In our previous study we found that enhanced mass loss is crucial for describing the sample, and here we continue that analysis but for rotational mixing and core overshooting. We show that, while rotation is important for stellar evolution studies, rotation, itself, is insufficient to model the distribution of period change rates from the observed sample. On the other hand, convective core overshooting is needed to explain the magnitude of the rates of period change, but does not explain the number of stars with positive and negative period change rates. In conclusion, we determine that convective core overshooting and stellar rotation alone are not enough to account for the observed distribution of Cepheid rates of period change and another mechanism, such as pulsation-driven mass-loss, may be required.

Cepheid masses derived from pulsations or binary dynamics are generally lower than those derived from stellar evolution models. Recent efforts have been dedicated to investigating the effects of abundances, mass loss, rotation, convection and overshooting prescriptions for modifying the evolution tracks to reduce or remove this Cepheid mass discrepancy. While these approaches are promising, either alone or in combination, more work is required to distinguish between possible solutions. Here we investigate nuclear reaction rate and opacity modifications on Cepheid evolution using the MESA code. We discuss the effects of opacity increases at envelope temperatures of 200,000-400,000 K proposed to explain the pulsation properties of hybrid main-sequence beta Cep/Slowly Pulsating B (SPB) variables which will evolve into Cepheids. We make use of the RSP nonlinear radial pulsation modeling capability in MESA to calculate periods and radial velocity amplitudes of Galactic Cepheids V1334 Cyg, Polaris, and delta Cep.

Polaris is the nearest Cepheid to us and as such holds a special place in our understanding of Cepheids in general and the Leavitt Law. In the past couple of decades, we have learned many new things about the star as a Cepheid and as the primary component of a multiple star system. As such, we are more precisely measuring the mass, radius and evolution of Polaris. However, as we learn more, it is becoming clear that we understand less. There is evidence that Polaris is much less massive than stellar evolution models suggest and that Polaris is a much younger star than its main sequence companion. In this work, we review some of the recent measurements and their connections with past studies. We then present new stellar evolution models and populations synthesis calculations to compare with the new mass measurements by Evans et al. (2018). We find that the mass discrepancy for Polaris is about 50\%. We also find that there is a likely age discrepancy between Polaris and its companion, but that there is also a very small probability that the discrepancy is not real.

V473 Lyr is a classical Cepheid which is unique in having substantial amplitude variations with a period of approximately 3.3 years, thought to be similar to the Blazhko variations in RR Lyrae stars. We obtained an {\it XMM-Newton} observation of this star to followup a previous detection in X-rays. Rather than the X-ray burst and rapid decline near maximum radius seen in $\delta$ Cephei itself, the X-ray flux in V473 Lyr remained constant for a third of the pulsation cycle covered by the observation. Thus the X-rays are most probably not produced by the changes around the pulsation cycle. The X-ray spectrum is soft (kT = 0.6 keV), with X-ray properties which are consistent with a young low mass companion. Previously there was no evidence of a companion in radial velocities or in {\it Gaia} and {\it Hipparcos} proper motions. While this rules out companions which are very close or very distant, a binary companion at a separation between 30 and 300 AU is possible. This is an example of an X-ray observation revealing evidence of a low mass companion, which is important in completing the mass ratio statistics of binary Cepheids. Furthermore, the detection of a young X-ray bright companion is a further indication that the Cepheid (primary) is a Population I star, even though its pulsation behavior differs from other classical Cepheids.

Maunakea, the proposed site of the Thirty Meter Telescope (TMT), is a lightning-rod topic for Native Hawaiians, Hawaii residents, and the international astronomy community. In this paper we, Native Hawaiian natural scientists and allies, identify historical decisions that impact current circumstances on Maunakea and provide approaches to acknowledging their presence. Our aim is to provide an Indigenous viewpoint centered in Native Hawaiian perspectives on the impacts of the TMT project on the Hawaiian community. We summarize the current Maunakea context from the perspective of the authors who are trained in the natural sciences (inclusive of and beyond astronomy and physics), the majority of whom are Native Hawaiian or Indigenous. We highlight three major themes in the conflict surrounding TMT: 1) physical demonstrations and the use of law enforcement against the protectors of Maunakea; 2) an assessment of the benefit of Maunakea astronomy to Native Hawaiians; and 3) the disconnect between astronomers and Native Hawaiians. We close with general short- and long- term recommendations for the astronomy community, which represent steps that can be taken to re-establish trust and engage in meaningful reciprocity and collaboration with Native Hawaiians and other Indigenous communities. Our recommendations are based on established best principles of free, prior, and informed consent and researcher-community interactions that extend beyond transactional exchanges. We emphasize that development of large-scale astronomical instrumentation must be predicated on consensus from the local Indigenous community about whether development is allowed on their homelands. Proactive steps must be taken to center Indigenous voices in the earliest stages of project design.

This white paper explains that professional astronomy has benefited from settler colonial white supremacist patriarchy. We explicate the impact that this has had on communities which are not the beneficiaries of colonialism and white supremacy. We advocate for astronomers to reject these benefits in the future, and we make proposals regarding the steps involved in rejecting colonialist white supremacy's benefits. We center ten recommendations on the timely issue of what to do about the Thirty Meter Telescope (TMT) on Maunakea in Hawaii. This paper is written in solidarity with and support of efforts by Native Hawaiian scientists (e.g. Kahanamoku et al. 2019).

Recent BRITE-Constellation space photometry of the slowly rotating, magnetic $\beta$ Cep pulsator $\xi^1$ CMa permits a new analysis of its pulsation properties. Analysis of the two-colour BRITE data reveals the well-known single pulsation period of $0.209$ d, along with its first and second harmonics. A similar analysis of SMEI and TESS observations yields compatible results, with the higher precision TESS observations also revealing several low-amplitude modes with frequencies below 5 d$^{-1}$; some of these are likely $g$ modes. The phase lag between photometric and radial velocity maxima - equal to 0.334 cycles - is significantly larger than the typical value of $1/4$ observed in other large-amplitude $\beta$ Cep stars. The phase lag, as well as the strong dependence of phase of maximum light on wavelength, can be reconciled with seismic models only if the dominant mode is the fundamental radial mode. We employ all published photometric and radial velocity measurements, spanning over a century, to evaluate the stability of the pulsation period. The $O-C$ diagram exhibits a clear parabolic shape consistent with a mean rate of period change $\dot P=0.34\pm 0.02$ s/cen. The residuals from the best-fit parabola exhibit scatter that is substantially larger than the uncertainties. In particular, dense sampling obtained during the past $\sim$20 years suggests more complex and rapid period variations. Those data cannot be coherently phased with the mean rate of period change, and instead require $\dot P\sim0.9$ s/cen. We examine the potential contributions of binarity, stellar evolution, and stellar rotation and magnetism to understand the apparent period evolution.

(Abridged) Canadian astronomy has, for decades, benefited from access to observatories and participating in international consortia on one of the best astronomical sites in the world: Maunakea. However, Maunakea is part of the unceded territory of the Native Hawaiian peoples and has always been of special significance to Hawaiian culture. The use of the summit and its science reserve has created tensions in the past decade, particularly with the development of the Thirty Meter Telescope. A meaningful and respectful response from the International astronomy community is still lacking. It is expected that the LRP 2020 will continue to support Canadian astronomy on Maunakea so a better official statement on the position and involvement of CASCA should be prepared. In this paper we present recommendations, based on the United Nation Declaration for the Rights of Indigenous Peoples, for the Canadian astronomical community to better support Indigenous rights on Maunakea and Hawaii while providing clear guidelines for the astronomical community to participate in activities conducted on Indigenous land. This framework is designed to motivate conversations with Indigenous communities regarding our place on Indigenous lands and our roles, and responsibilities toward the communities we are working with. Furthermore, we propose this framework as a basis for engaging with communities around the world regarding consent for astronomical facilities.

(Abridged) The Truth and Reconciliation Commission of Canada published its calls to action in 2015 with 94 recommendations. Many of these 94 recommendations are directly related to education, language, and culture, some of which the Canadian Astronomy community can address and contribute to as part of reconciliation. The Canadian Astronomy community has an additional obligation since it benefits from facilities on Indigenous territories across Canada and the world. Furthermore, Indigenous people are still underrepresented at all levels in Canadian astronomy. The purpose of this Community Paper is to develop recommendations for the Canadian astronomy community to support Indigenous inclusion in the science community, support Indigenous learning by developing Indigenous-based learning materials and facilitate access to professionals and science activities, and to recognize and acknowledge the great contributions of Indigenous communities to our science activities. As part of this work we propose the ten following recommendations for CASCA as an organization and throughout this Community Paper we will include additional recommendations for individuals: astronomers, students and academics.

We introduce an observational tool based on visibility nulls in optical spectro-interferometry fringe data to probe the structure of stellar atmospheres. In a preliminary demonstration, we use both Navy Precision Optical Interferometer (NPOI) data and stellar atmosphere models to show that this tool can be used, for example, to investigate limb darkening. Using bootstrapping with either multiple linked baselines or multiple wavelengths in optical and infrared spectro-interferometric observations of stars makes it possible to measure the spatial frequency $u_0$ at which the real part of the fringe visibility ${\rm Re}(V)$ vanishes. That spatial frequency is determined by $u_0 = B_\perp/\lambda_0$, where $B_\perp$ is the projected baseline length, and $\lambda_0$ is the wavelength at which the null is observed. Since $B_\perp$ changes with the Earth's rotation, $\lambda_0$ also changes. If $u_0$ is constant with wavelength, $\lambda_0$ varies in direct proportion to $B_\perp$. Any departure from that proportionality indicates that the brightness distribution across the stellar disk varies with wavelength via variations in limb darkening, in the angular size of the disk, or both. In this paper, we introduce the use of variations of $u_0$ with $\lambda$ as a means of probing the structure of stellar atmospheres. Using the equivalent uniform disk diameter $\theta_{\rm UD, 0}(\lambda_0)$, given by $\theta_{\rm UD, 0} = 1.22/u_0(\lambda_0)$, as a convenient and intuitive parameterization of $u_0(\lambda_0)$, we demonstrate this concept by using model atmospheres to calculate the brightness distribution for $\nu$ Ophiuchi and predict $\theta_{\rm UD, 0}(\lambda_0)$, and then comparing the predictions to coherently averaged data from observations taken with the NPOI.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

The dynamical structure of the atmosphere of Cepheids has been well studied in the optical. Several authors have found very interesting spectral features in the J band, but little data have been secured beyond 1.6 um. However, such observations can probe different radial velocities and line asymmetry regimes, and are able to provide crucial insights into stellar physics. Our goal was to investigate the infrared line-forming region in the K-band domain, and its impact on the projection factor and the k-term of Cepheids. We secured CRIRES observations for the long-period Cepheid l Car, with a focus on the unblended spectral line NaI2208.969 nm. We measured the corresponding radial velocities (by using the first moment method) and the line asymmetries (by using the bi-Gaussian method). These quantities are compared to the HARPS visible spectra we previously obtained on l Car. The optical and infrared radial velocity curves show the same amplitude (only about 3% of difference), with a slight radial velocity shift of about 0.5 +/- 0.3 km s^-1 between the two curves. Around the minimum radius (phase ~ 0.9) the visible radial velocity curve is found in advance compared to the infrared one (phase lag), which is consistent with an infrared line forming higher in the atmosphere (compared to the visible line) and with a compression wave moving from the bottom to the top of the atmosphere during maximum outward velocity. The asymmetry of the K-band line is also found to be significantly different from that of the optical line.

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.

Bow shocks and related density enhancements produced by the winds of massive stars moving through the interstellar medium provide important information regarding the motions of the stars, the properties of their stellar winds, and the characteristics of the local medium. Since bow shocks are aspherical structures, light scattering within them produces a net polarization signal even if the region is spatially unresolved. Scattering opacity arising from free electrons and dust leads to a distribution of polarized intensity across the bow shock structure. That polarization encodes information about the shape, composition, opacity, density, and ionisation state of the material within the structure. In this paper we use the Monte Carlo radiative transfer code SLIP to investigate the polarization created when photons scatter in a bow shock-shaped region of enhanced density surrounding a stellar source. We present results assuming electron scattering, and investigate the polarization behaviour as a function of optical depth, temperature, and source of photons for two different cases: pure scattering and scattering with absorption. In both regimes we consider resolved and unresolved cases. We discuss the implication of these results as well as their possible use along with observational data to constrain the properties of observed bow shock systems. In different situations and under certain assumptions, our simulations can constrain viewing angle, optical depth and temperature of the scattering region, and the relative luminosities of the star and shock.

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.

We present the discovery of a rapidly evolving transient by the Korean Microlensing Telescope Network Supernova Program (KSP). KSP is a novel high-cadence supernova survey that offers deep ($\sim21.5$ mag in $BVI$ bands) nearly continuous wide-field monitoring for the discovery of early and/or fast optical transients. KSP-OT-201509a, reported here, was discovered on 2015 September 27 during the KSP commissioning run in the direction of the nearby galaxy NGC~300, and stayed above detection limit for $\sim$ 22 days. We use our $BVI$ light-curves to constrain the ascent rate, $-3.7(7)$ mag day$^{-1}$ in $V$, decay time scale, $t^{V}_{2}=1.7(6)$ days, and peak absolute magnitude, $-9.65\leq M_{V}\leq -9.25$ mag. We also find evidence for a short-lived pre-maximum halt in all bands. The peak luminosity and lightcurve evolution make KSP-OT-201509a consistent with a bright, rapidly decaying nova outburst. We discuss constraints on the nature of the progenitor and its environment using archival HST/ACS images and conclude with a broad discussion on the nature of the system.

The projection factor p is the key quantity used in the Baade-Wesselink (BW) method for distance determination; it converts radial velocities into pulsation velocities. Several methods are used to determine p, such as geometrical and hydrodynamical models or the inverse BW approach when the distance is known. We analyze new HARPS-N spectra of delta Cep to measure its cycle-averaged atmospheric velocity gradient in order to better constrain the projection factor. We first apply the inverse BW method to derive p directly from observations. The projection factor can be divided into three subconcepts: (1) a geometrical effect (p0); (2) the velocity gradient within the atmosphere (fgrad); and (3) the relative motion of the optical pulsating photosphere with respect to the corresponding mass elements (fo-g). We then measure the fgrad value of delta Cep for the first time. When the HARPS-N mean cross-correlated line-profiles are fitted with a Gaussian profile, the projection factor is pcc-g = 1.239 +/- 0.034(stat) +/- 0.023(syst). When we consider the different amplitudes of the radial velocity curves that are associated with 17 selected spectral lines, we measure projection factors ranging from 1.273 to 1.329. We find a relation between fgrad and the line depth measured when the Cepheid is at minimum radius. This relation is consistent with that obtained from our best hydrodynamical model of delta Cep and with our projection factor decomposition. Using the observational values of p and fgrad found for the 17 spectral lines, we derive a semi-theoretical value of fo-g. We alternatively obtain fo-g = 0.975+/-0.002 or 1.006+/-0.002 assuming models using radiative transfer in plane-parallel or spherically symmetric geometries, respectively. The new HARPS-N observations of delta Cep are consistent with our decomposition of the projection factor.

We review ways in which observations of the changing periods of pulsating variable stars can be used to detect and directly measure their evolution. We briefly describe the two main techniques of analysis -- (O-C) analysis and wavelet analysis -- and results for pulsating variable star types which are reasonably periodic: type I and II Cepheids, RR Lyrae stars, beta Cephei stars, and Mira stars. We comment briefly on delta Scuti stars and pulsating white dwarfs. For some of these variable star types, observations agree approximately with the predictions of evolutionary models, but there still exist significant areas of disagreement that challenge future models of stellar evolution. There may be a need, for instance, to include processes such as rotation, mass loss, and magnetic fields. There may also be non-evolutionary processes which are contributing to the period changes.

The B-W method is used to determine the distance of Cepheids and consists in combining the angular size variations of the star, as derived from infrared surface-brightness relations or interferometry, with its linear size variation, as deduced from visible spectroscopy using the projection factor. While many Cepheids have been intensively observed by infrared beam combiners, only a few have been observed in the visible. This paper is part of a project to observe Cepheids in the visible with interferometry as a counterpart to infrared observations already in hand. Observations of delta Cep itself were secured with the VEGA/CHARA instrument over the full pulsation cycle of the star. These visible interferometric data are consistent in first approximation with a quasi-hydrostatic model of pulsation surrounded by a static circumstellar environment (CSE) with a size of theta_cse=8.9 +/- 3.0 mas and a relative flux contribution of f_cse=0.07+/-0.01. A model of visible nebula (a background source filling the field of view of the interferometer) with the same relative flux contribution is also consistent with our data at small spatial frequencies. However, in both cases, we find discrepancies in the squared visibilities at high spatial frequencies (maximum 2sigma) with two different regimes over the pulsation cycle of the star, phi=0.0-0.8 and phi=0.8-1.0. We provide several hypotheses to explain these discrepancies, but more observations and theoretical investigations are necessary before a firm conclusion can be drawn. For the first time we have been able to detect in the visible domain a resolved structure around delta~Cep. We have also shown that a simple model cannot explain the observations, and more work will be necessary in the future, both on observations and modelling.

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.

The classical Cepheid $l$ Carinae is an essential calibrator of the Cepheid Leavitt Law as a rare long-period Galactic Cepheid. Understanding the properties of this star will also constrain the physics and evolution of massive ($M \ge 8$ $M_\odot$) Cepheids. The challenge, however, is precisely measuring the star's pulsation period and its rate of period change. The former is important for calibrating the Leavitt Law and the latter for stellar evolution modeling. In this work, we combine previous time-series observations spanning more than a century with new observations to remeasure the pulsation period and compute the rate of period change. We compare our new rate of period change with stellar evolution models to measure the properties of $l$ Car, but find models and observations are, at best, marginally consistent. The results imply that $l$ Car does not have significantly enhanced mass-loss rates like that measured for $\delta$ Cephei. We find that the mass of $l$ Car is about 8 - 10 $M_\odot$. We present Hubble Space Telescope COS observations that also differ from measurements for $\delta$ Cep, and $\beta$ Dor. These measurements further add to the challenge of understanding the physics of Cepheids, but do hint at the possible relation between enhanced mass loss and ultraviolet emission, perhaps both due to the strength of shocks propagating in the atmospheres of Cepheids.

The $\beta$ Cephei stars represent an important class of massive star pulsators probing the evolution of B-type stars and the transition from main sequence to hydrogen-shell burning evolution. By understanding $\beta$ Cep stars, we gain insights into the detailed physics of massive star evolution such as rotational mixing, convective core overshooting, magnetic fields and stellar winds, all of which play important roles. Similarly, modeling their pulsation provides additional information into their interior structures. Furthermore, measurements of the rate of change of pulsation period offer a direct measure of $\beta$ Cephei stellar evolution. In this work, we compute state-of-the-art stellar evolution models assuming different amounts of initial rotation and convective core overshoot and measure theoretical rates of period change for which we compare to rates previously measured for a sample of $\beta$ Cephei stars. The results of this comparison are mixed. For three stars, the rates are too small to infer any information from stellar evolution models, whereas for three other stars the rates are too large. We infer stellar parameters, such as mass and age, for two $\beta$ Cephei stars: $\xi^1$ CMa and $\delta$ Cet, that agree well with independent measurements. We explore ideas for why models may not predict the larger rates of period change. In particular, period drifts in $\beta$ Cep stars can artificially lead to overestimated rates of secular period change.

Classical Cepheids are key probes of both stellar astrophysics and cosmology as standard candles and pulsating variable stars. It is important to understand Cepheids in unprecedented detail in preparation for upcoming GAIA, JWST and extremely-large telescope observations. Cepheid eclipsing binary stars are ideal tools for achieving this goal, however there are currently only three known systems. One of those systems, OGLE-LMC-CEP1812, raises new questions about the evolution of classical Cepheids because of an apparent age discrepancy between the Cepheid and its red giant companion. We show that the Cepheid component is actually the product of a stellar merger of two main sequence stars that has since evolved across the Hertzsprung gap of the HR diagram. This post-merger product appears younger than the companion, hence the apparent age discrepancy is resolved. We discuss this idea and consequences for understanding Cepheid evolution.

Classical Cepheids, like binary stars, are laboratories for stellar evolution and Cepheids in binary systems are especially powerful ones. About one-third of Galactic Cepheids are known to have companions and Cepheids in eclipsing binary systems have recently been discovered in the Large Magellanic Cloud. However, there are no known Galactic binary Cepheids with orbital periods less than one year. We compute population synthesis models of binary Cepheids to compare to the observed period and eccentricity distributions of Galactic Cepheids as well as to the number of observed eclipsing binary Cepheids in the LMC. We find that our population synthesis models are consistent with observed binary properties of Cepheids. Furthermore, we show that binary interaction on the red giant branch prevents some red giant stars from becoming classical Cepheids. Such interactions suggest that the binary fraction of Cepheids should be significantly less than that of their main-sequence progenitors, and that almost all binary Cepheids have orbital periods longer than one year. If the Galactic Cepheid spectroscopic binary fraction is about 35%, then the spectroscopic binary fraction of their intermediate mass main sequence progenitors is about 40-45%.

We present results from a study of the eclipsing, colliding-wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMM-Newton observations. In addition, we have characterized the intrinsic, polarimetric phase-dependent behavior of the strongest optical emission lines using data obtained with the University of Wisconsin's Half-Wave Spectropolarimeter. We have detected evidence of the Coriolis distortion of the wind-wind collision in the X-ray regime, which manifests itself through asymmetric behavior around the eclipses in the system's X-ray light curves. The large opening angle of the X-ray emitting region, as well as its location (i.e. the WN wind does not collide with the O star, but rather its wind) are evidence of radiative braking/inhibition occurring within the system. Additionally, the polarimetric results show evidence of the cavity the wind-wind collision region carves out of the Wolf-Rayet star's wind.

Over the past decade, the Secret Lives of Cepheids (SLiC) program has been carried out at Villanova University to study aspects and behaviors of classical Cepheids that are still not well-understood. In this, the first of several planned papers on program Cepheids, we report the current results for delta Cep, the Cepheid prototype. Ongoing photometry has been obtained to search for changes in the pulsation period, light curve morphology and amplitude. Combining our photometry with the times of maximum light compilation by Berdnikov 2000 returns a small period change of dP/dt ~ -0.1006 +/- 0.0002 sec yr^-1. There is also evidence for a gradual light amplitude increase of ~0.011-mag (V-band) and ~0.012-mag (B-band) per decade over the last ~50 years. In addition, HST-COS UV spectrophotometry and XMM-Newton X-ray data were carried out to investigate the high-temperature plasmas present above the Cepheid photospheres. In total, from the five visits (eight exposures) with XMM-Newton, delta Cep is found to be a soft X-ray source (Lx [0.3--2 keV] ~ 4.5-13x10^28 erg sec^-1) with peak flux at kT = 0.6-0.9 keV. The X-ray activity is found to vary, possibly in phase with the stellar pulsations. From 2010-2013, nine observations of delta Cep were carried out with HST-COS. The UV emissions are also variable, and well-phased with the stellar pulsations. Maximum UV line emissions occur near, or slightly before, maximum optical light, varying by as much as twenty times. This variability shows that pulsation-induced shock-heating plays a significant role in Cepheid atmospheres, possibly in addition to a quiescent, magnetic heating. The results of this study show Cepheid atmospheres to be rather complex and dynamic.

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

We discuss two important effects for the astrospheres of runaway stars: the propagation of ionizing photons far beyond the astropause, and the rapid evolution of massive stars (and their winds) near the end of their lives. Hot stars emit ionizing photons with associated photoheating that has a significant dynamical effect on their surroundings. 3D simulations show that HII regions around runaway O stars drive expanding conical shells and leave underdense wakes in the medium they pass through. For late O stars this feedback to the interstellar medium is more important than that from stellar winds. Late in life, O stars evolve to cool red supergiants more rapidly than their environment can react, producing transient circumstellar structures such as double bow shocks. This provides an explanation for the bow shock and linear bar-shaped structure observed around Betelgeuse.

Measurements of secular period change probe real-time stellar evolution of classical Cepheids making these measurements powerful constraints for stellar evolution models, especially when coupled with interferometric measurements. In this work, we present stellar evolution models and measured rates of period change for two Galactic Cepheids: Polaris and l Carinae, both important Cepheids for anchoring the Cepheid Leavitt law (period-luminosity relation). The combination of previously-measured parallaxes, interferometric angular diameters and rates of period change allows for predictions of Cepheid mass loss and stellar mass. Using the stellar evolution models, We find that l Car has a mass of about 9 $M_\odot$ consistent with stellar pulsation models, but is not undergoing enhanced stellar mass loss. Conversely, the rate of period change for Polaris requires including enhanced mass-loss rates. We discuss what these different results imply for Cepheid evolution and the mass-loss mechanism on the Cepheid instability strip.

Polarization is a powerful tool for understanding stellar atmospheres and circumstellar environments. Mira and semi-regular variable stars have been observed for decades and some are known to be polarimetrically variable, however, the semi-regular variable V Canes Venatici displays an unusually large, unexplained amount of polarization. We present ten years of optical polarization observations obtained with the HPOL instrument, supplemented by published observations spanning a total interval of about forty years for V CVn. We find that V CVn shows large polarization variations ranging from 1 - 6%. We also find that for the past forty years the position angle measured for V CVn has been virtually constant suggesting a long-term, stable, asymmetric structure about the star. We suggest that this asymmetry is caused by the presence of a stellar wind bow shock and tail, consistent with the star's large space velocity.

Two red supergiants of the Per OB1 association, RS Per and T Per, have been observed in H band using the MIRC instrument at the CHARA array. The data show clear evidence of departure from circular symmetry. We present here new techniques specially developed to analyze such cases, based on state-of-the-art statistical frameworks. The stellar surfaces are first modeled as limb-darkened discs based on SATLAS models that fit both MIRC interferometric data and publicly available spectrophotometric data. Bayesian model selection is then used to determine the most probable number of spots. The effective surface temperatures are also determined and give further support to the recently derived hotter temperature scales of red su- pergiants. The stellar surfaces are reconstructed by our model-independent imaging code SQUEEZE, making use of its novel regularizer based on Compressed Sensing theory. We find excellent agreement between the model-selection results and the reconstructions. Our results provide evidence for the presence of near-infrared spots representing about 3-5% of the stellar flux.

The circumstellar medium around massive stars is strongly impacted by stellar winds, radiation, and explosions. We use numerical simulations of these interactions to constrain the current properties and evolutionary history of various stars by comparison with observed circumstellar structures. Two- and three-dimensional simulations of bow shocks around red supergiant stars have shown that Betelgeuse has probably only recently evolved from a blue supergiant to a red supergiant, and hence its bow shock is very young and has not yet reached a steady state. We have also for the first time investigated the magnetohydrodynamics of the photoionised H II region around the nearby runaway O star Zeta Oph. Finally, we have calculated a grid of models of bow shocks around main sequence and evolved massive stars that has general application to many observed bow shocks, and which forms the basis of future work to model the explosions of these stars into their pre-shaped circumstellar medium.

Polaris the Cepheid has been observed for centuries, presenting surprises and changing our view of Cepheids and stellar astrophysics, in general. Specifically, understanding Polaris helps anchor the Cepheid Leavitt law, but the distance must be measured precisely. The recent debate regarding the distance to Polaris has raised questions about its role in calibrating the Leavitt law and even its evolutionary status. In this work, I present new stellar evolution models of Cepheids to compare with previously measured CNO abundances, period change and angular diameter. Based on the comparison, I show that Polaris cannot be evolving along the first crossing of the Cepheid instability strip and cannot have evolved from a rapidly-rotating main sequence star. As such, Polaris must also be at least 118 pc away and pulsates in the first overtone, disagreeing with the recent results of Turner et al. (2013).

Analyses of recent observations of the sole classical Cepheid in the Kepler field, V1154 Cygni, found random changes of about 30 minutes in the pulsation period. These period changes challenge standard theories of pulsation and evolution because the period change is non-secular, and explaining this period jitter is necessary for understanding stellar evolution and the role of Cepheids as precise standard candles. We suggest that convection and convective hot spots can explain the observed period jitter. Convective hot spots alter the timing of flux maximum and minimum in the Cepheid light curve, hence change the measured pulsation period. We present a model of random hot spots that generate a localized flux excess that perturbs the Cepheid light curve and consequently the pulsation period which is consistent with the observed jitter. This result demonstrates how important understanding convection is for modeling Cepheid stellar structure and evolution, how convection determines the red edge of the instability strip, and just how sensitive Cepheid light curves are to atmospheric physics.

Both pulsation and mass loss are commonly observed in stars and are important ingredients for understanding stellar evolution and structure, especially for massive stars. There is a growing body of evidence that pulsation can also drive and enhance mass loss in massive stars and that pulsation-driven mass loss is important for stellar evolution. In this review, I will discuss recent advances in understanding pulsation driven mass loss in massive main sequence stars, classical Cepheids and red supergiants and present some challenges remaining.

Mira and semi-regular variable stars have been studied for centuries but continue to be enigmatic. One unsolved mystery is the presence of polarization from these stars. In particular, we present 40 years of polarization measurements for the prototype o Ceti and V CVn and find very different phenomena for each star. The polarization fraction and position angle for Mira is found to be small and highly variable. On the other hand, the polarization fraction for V CVn is large and variable, from 2 - 7 %, and its position angle is approximately constant, suggesting a long-term asymmetric structure. We suggest a number of potential scenarios to explain these observations.

Classical Cepheids form one of the foundations of modern cosmology and the extragalactic distance scale, however, cosmic microwave background observations measure cosmological parameters and indirectly the Hubble Constant, H0, to unparalleled precision. The coming decade will provide opportunities to measure H0 to 2% uncertainty thanks to the GAIA satellite, JWST, ELTs and other telescopes using Cepheids and other standard candles. In this work, we discuss the upcoming role for variable stars and asteroseismology in calibrating the distance scale and measuring H0 and what problems exist in understanding these stars that will feedback on these measurements.

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.

The nearby red supergiant (RSG) Betelgeuse has a complex circumstellar medium out to at least 0.5 parsecs from its surface, shaped by its mass-loss history within the past 0.1 Myr, its environment, and its motion through the interstellar medium (ISM). In principle its mass-loss history can be constrained by comparing hydrodynamic models with observations. Observations and numerical simulations indicate that Betelgeuse has a very young bow shock, hence the star may have only recently become a RSG. To test this possibility we calculated a stellar evolution model for a single star with properties consistent with Betelgeuse. We incorporated the resulting evolving stellar wind into 2D hydrodynamic simulations to model a runaway blue supergiant (BSG) undergoing the transition to a RSG near the end of its life. The collapsing BSG wind bubble induces a bow shock-shaped inner shell which at least superficially resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We investigate whether this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

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.

Measurements of rates of period change of Classical Cepheids probe stellar physics and evolution. Additionally, better understanding of Cepheid structure and evolution provides greater insight into their use as standard candles and tools for measuring the Hubble constant. Our recent study of the period change of the nearest Cepheid, Polaris, suggested that it is undergoing enhanced mass loss when compared to canonical stellar evolution model predictions. In this work, we expand the analysis to rates of period change measured for about 200 Galactic Cepheids and compare them to population synthesis models of Cepheids including convective core overshooting and enhanced mass loss. Rates of period change predicted from stellar evolution models without mass loss do not agree with observed rates whereas including enhanced mass loss yields predicted rates in better agreement with observations. This is the first evidence that enhanced mass loss as suggested previously for Polaris and delta Cephei must be a ubiquitous property of Classical Cepheids.

Galactic Cepheids are necessary tools for calibrating the period-luminosity relation, but distances to individual Galactic Cepheids are difficult to precisely measure and are limited to a small number of techniques such as direct parallax, main-sequence fitting to open clusters that host Cepheids and Baade-Wesselink (BW) methods. In this work, we re-examine the application of Wesenheit function in determining distances to more than 300 Galactic Cepheids, by taking advantage of the fact that the Wesenheit function is extinction free by definition. The Wesenheit distances were used to calibrate the projection factor (p-factor) for Galactic Cepheids that also possess BW distances. Based on ~70 Cepheids, we found that a period-p factor relation may exhibit a non-linear trend with a considerable scatter. During our investigation, discrepant p factors for delta Cephei were found in the literature. This may be due to inconsistent measurements of the angular diameters using different empirical techniques. We discuss the reason for the inconsistency of angular diameter measurements and offer a possible remedy for this problem.

Our ability to extract information from the spectra of stars depends on reliable models of stellar atmospheres and appropriate techniques for spectral synthesis. Various model codes and strategies for the analysis of stellar spectra are available today. We aim to compare the results of deriving stellar parameters using different atmosphere models and different analysis strategies. The focus is set on high-resolution spectroscopy of cool giant stars. Spectra representing four cool giant stars were made available to various groups and individuals working in the area of spectral synthesis, asking them to derive stellar parameters from the data provided. The results were discussed at a workshop in Vienna in 2010. Most of the major codes currently used in the astronomical community for analyses of stellar spectra were included in this experiment. We present the results from the different groups, as well as an additional experiment comparing the synthetic spectra produced by various codes for a given set of stellar parameters. Similarities and differences of the results are discussed. Several valid approaches to analyze a given spectrum of a star result in quite a wide range of solutions. The main causes for the differences in parameters derived by different groups seem to lie in the physical input data and in the details of the analysis method. This clearly shows how far from a definitive abundance analysis we still are.

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.

The projection factor (p), which converts the radial velocity to pulsational velocity, is an important parameter in the Baade-Wesselink (BW) type analysis and distance scale work. The p-factor is either adopted as a constant or linearly depending on the logarithmic of pulsating periods. The aim of this work is to calibrate the p-factor if a Cepheid has both the BW distance and an independent distance measurement, and examine the p-factor for delta Cephei -- the prototype of classical Cepheids. We calibrated the p-factor for several Galactic Cepheids that have both the latest BW distances and independent distances either from Hipparcos parallaxes or main-sequence fitting distances to Cepheid-hosted stellar clusters. Based on 25 Cepheids, the calibrated p-factor relation is consistent with latest p-factor relation in literature. The calibrated p-factor relation also indicates that this relation may not be linear and may exhibit an intrinsic scatter. We also examined the discrepancy of empirical p-factors for delta Cephei, and found that the reasons for this discrepancy include the disagreement of angular diameters, the treatment of radial velocity data, and the phase interval adopted during the fitting procedure. Finally, we investigated the impact of the input p-factor in two BW methodologies for delta Cephei, and found that different p-factors can be adopted in these BW methodologies and yet result in the same angular diameters.

A significant fraction of massive stars are moving supersonically through the interstellar medium (ISM), either due to disruption of a binary system or ejection from their parent star cluster. The interaction of their wind with the ISM produces a bow shock. In late evolutionary stages these stars may undergo rapid transitions from red to blue and vice versa on the Hertzsprung-Russell diagram, with accompanying rapid changes to their stellar winds and bow shocks. Recent 3D simulations of the bow shock produced by the nearby runaway red supergiant (RSG) Betelgeuse, under the assumption of a constant wind, indicate that the bow shock is very young (<30000 years old), hence Betelgeuse may have only recently become a RSG. To test this possibility, we have calculated stellar evolution models for single stars which match the observed properties of Betelgeuse in the RSG phase. The resulting evolving stellar wind is incorporated into 2D hydrodynamic simulations in which we model a runaway blue supergiant (BSG) as it undergoes the transition to a RSG near the end of its life. We find that the collapsing BSG wind bubble induces a bow shock-shaped inner shell around the RSG wind that resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We suggest that this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

Context. One challenge for measuring the Hubble constant using Classical Cepheids is the calibration of the Leavitt Law or period-luminosity relationship. The Baade-Wesselink method for distance determination to Cepheids relies on the ratio of the measured radial velocity and pulsation velocity, the so-called projection factor and the ability to measure the stellar angular diameters. Aims. We use spherically-symmetric model stellar atmospheres to explore the dependence of the p-factor and angular diameter corrections as a function of pulsation period. Methods. Intensity profiles are computed from a grid of plane-parallel and spherically-symmetric model stellar atmospheres using the SAtlas code. Projection factors and angular diameter corrections are determined from these intensity profiles and compared to previous results. Results. Our predicted geometric period-projection factor relation including previously published state-of-the-art hydrodynamical predictions is not with recent observational constraints. We suggest a number of potential resolutions to this discrepancy. The model atmosphere geometry also affects predictions for angular diameter corrections used to interpret interferometric observations, suggesting corrections used in the past underestimated Cepheid angular diameters by 3 - 5%. Conclusions. While spherically-symmetric hydrostatic model atmospheres cannot resolve differences between projection factors from theory and observations, they do help constrain underlying physics that must be included, including chromospheres and mass loss. The models also predict more physically-based limb-darkening corrections for interferometric observations.

Polaris is one of the most observed stars in the night sky, with recorded observations spanning more than 200 years. From these observations, one can study the real-time evolution of Polaris via the secular rate of change of the pulsation period. However, the measurements of the rate of period change do not agree with predictions from state-of-the-art stellar evolution models. We show that this may imply that Polaris is currently losing mass at a rate of $\dot{M} \approx 10^{-6} M_\odot$ yr$^{-1}$ based on the difference between modeled and observed rates of period change, consistent with pulsation-enhanced Cepheid mass loss. A relation between the rate of period change and mass loss has important implications for understanding stellar evolution and pulsation, and provides insight into the current Cepheid mass discrepancy.

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.

Context. The Cepheid mass discrepancy, the difference between masses predicted from stellar evolution and stellar pulsation calculations, is a challenge for the understanding of stellar astrophysics. Recent models of the eclipsing binary Cepheid OGLE-LMC-CEP-0227 have suggested that the discrepancy may be resolved. Aims. We explore for what physical parameters do stellar evolution models agree with the measured properties of OGLE-LMC-CEP0227 and compare to canonical stellar evolution models assuming no convective core overshooting. Methods. We construct state-of-the-art stellar evolution models for varying mass, metallicity, and convective core overshooting and compare the stellar evolution predictions with the observed properties. Results. The observed mass, effective temperature, and radius of the two stars in the binary system are well fit by numerous combinations of physical parameters, suggesting a Cepheid mass discrepancy of 10-20% relative to canonical stellar evolution models. Conclusions. The properties of the observed binary Cepheid suggest that the Cepheid mass discrepancy is still a challenge and requires more specific observations, such as the rate of period change, to better constrain and understand the necessary physics for stellar evolution models to resolve the discrepancy.

Limb-darkening is a fundamental constraint for modeling eclipsing binary and planetary transit light curves. As observations, for example from Kepler, CoRot, and Most, become more precise then a greater understanding of limb-darkening is necessary. However, limb-darkening is typically modeled as simple parameterizations fit to plane-parallel model stellar atmospheres that ignores stellar atmospheric extension. In this work, I compute linear, quadratic and four-parameter limb-darkening laws from grids of plane-parallel and spherically-symmetric model stellar atmospheres in a temperature and gravity range representing stars evolving on the Red Giant branch. The limb-darkening relations for each geometry are compared and are found to fit plane-parallel models much better than the spherically-symmetric models. Assuming that limb-darkening from spherically-symmetry model atmospheres are more physically representative of actual stellar limb-darkening than plane-parallel models then these limb-darkening laws will not fit the limb of a stellar disk leading to errors in a light curve fit. This error will increase with a star's atmospheric extension.

The Cepheid mass discrepancy, the difference between mass estimates from stellar pulsation and stellar evolution models, is a long standing challenge for the understanding of stellar astrophysics. We discuss the current state of the mass discrepancy and test the role of pulsation-driven mass loss in state-of-the-art stellar evolution calculations of Large Magellanic Cepheids in resolving it. We find that Cepheid mass loss is a significant contributor to the mass discrepancy, but it is not clear if the metallicity dependence of Cepheid mass loss is consistent with the measured metallicity dependence of the mass discrepancy.

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.

Context. A longstanding challenge for understanding classical Cepheids is the Cepheid mass discrepancy, where theoretical mass estimates using stellar evolution and stellar pulsation calculations have been found to differ by approximately 10 - 20%. Aims. We study the role of pulsation-driven mass loss during the Cepheid stage of evolution as a possible solution to this mass discrepancy. Methods. We computed stellar evolution models with a Cepheid mass-loss prescription and various amounts of convective core overshooting. The contribution of mass loss towards the mass discrepancy is determined using these models, Results. Pulsation-driven mass loss is found to trap Cepheid evolution on the instability strip, allowing them to lose about 5 - 10% of their total mass when moderate convective core overshooting, an amount consistent with observations of other stars, is included in the stellar models. Conclusions. We find that the combination of moderate convective core overshooting and pulsation-driven mass loss can solve the Cepheid mass discrepancy.

We used WIRCam on CFHT to observe four transits of the super-Earth GJ 1214b in the near-infrared. For each transit we observed in two bands nearly-simultaneously by rapidly switching the WIRCam filter wheel back and forth for the duration of the observations. By combining all our J-band (~1.25 microns) observations we find a transit depth in this band of 1.338\pm0.013% - a value consistent with the optical transit depth reported by Charbonneau and collaborators. However, our best-fit combined Ks-band (~2.15 microns) transit depth is deeper: 1.438\pm0.019%. Formally our Ks-band transits are deeper than the J-band transits observed simultaneously by a factor of 1.072\pm0.018 - a 4-sigma discrepancy. The most straightforward explanation for our deeper Ks-band depth is a spectral absorption feature from the limb of the atmosphere of the planet; for the spectral absorption feature to be this prominent the atmosphere of GJ 1214b must have a large scale height and a low mean molecular weight. That is, it would have to be hydrogen/helium dominated and this planet would be better described as a mini-Neptune. However, recently published observations from 0.78 - 1.0 microns, by Bean and collaborators, show a lack of spectral features and transit depths consistent with those obtained by Charbonneau and collaborators. The most likely atmospheric composition for GJ 1214b that arises from combining all these observations is less clear; if the atmosphere of GJ 1214b is hydrogen/helium dominated then it must have either a haze layer that is obscuring transit depth differences at shorter wavelengths, or significantly different spectral features than current models predict. Our observations disfavour a water-world composition, but such a composition will remain a possibility until observations reconfirm our deeper Ks-band transit depth or detect features at other wavelengths. [Abridged]

We explore the role of mass loss and convective core overshoot in the evolution of Classical Cepheids. Stellar evolution models are computed with a recipe for pulsation-driven mass loss and it is found that mass loss alone is unable to account for the long-standing Cepheid mass discrepancy. However, the combination of mass loss and moderate convective core overshooting does provide a solution, bringing the amount of convective core overshooting in Cepheids closer to that found for other stars.

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 paper, we derive the period-luminosity (P-L) relation for Large Magellanic Cloud (LMC) Cepheids based on mid-infrared AKARI observations. AKARI's IRC sources were matched to the OGLE-III LMC Cepheid catalog. Together with the available I band light curves from the OGLE-III catalog, potential false matches were removed from the sample. This procedure excluded most of the sources in the S7 and S11 bands: hence only the P-L relation in the N3 band was derived in this paper. Random-phase corrections were included in deriving the P-L relation for the single epoch AKARI data, even though the derived P-L relation is consistent with the P-L relation without random-phase correction, though there is a \sim 7 per-cent improvement in the dispersion of the P-L relation. The final adopted N3 band P-L relation is N3 = -3.246 log(P) + 15.844, with a dispersion of 0.149.

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.

In this Paper, we have derived Cepheid period-luminosity (P-L) relations for the Large Magellanic Cloud (LMC) fundamental mode Cepheids, based on the data released from OGLE-III. We have applied an extinction map to correct for the extinction of these Cepheids. In addition to the VIW band P-L relations, we also include JHK and four Spitzer IRAC band P-L relations, derived by matching the OGLE-III Cepheids to the 2MASS and SAGE datasets, respectively. We also test the non-linearity of the Cepheid P-L relations based on extinction-corrected data. Our results (again) show that the LMC P-L relations are non-linear in VIJH bands and linear in KW and the four IRAC bands, respectively.

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: 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.

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.

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.