Source author record

Hilding R. Neilson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

42 published item(s)

preprint2020arXiv

Phase-modulated X-ray Emission from Cepheids due to Pulsation-Driven Shocks

Cepheids are pulsating variable stars with a periodic chromospheric response at UV wavelengths close to their minimum radius phase. Recently, an X-ray variable signature was captured in observations during the \emph{maximum} radius phase. This X-ray emission came as a surprise and is not understood. In this work, we use the modern astrophysical code, PLUTO, to investigate the effects of pulsations on Cepheid X-ray emission. We run a number of hydrodynamic numerical simulations with a variety of initial and boundary conditions in order to explore the capability of shocks to produce the observed phase-dependent X-ray behavior. Finally we use the Simulated Observations of X-ray Sources (SOXS) package to create synthetic spectra for each simulation case and link our simulations to observables. We show that, for certain conditions, we can reproduce observed X-ray fluxes at phases 0.4--0.8 when the Cepheid is at maximum radius. Our results span a wide range of mass-loss rates, $2\times10^{-13}$--$3\times10^{-8}$ $M_\odot$ yr$^{-1}$, and peak X-ray luminosities, $5\times10^{-17}$--$1.4\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$. We conclude that Cepheids exhibit \textit{two component emission} with (a) shock waves being responsible for the phase dependent variable emission (phases 0.2 - 0.6), and (b) a separate quiescent mechanism being the dominant emission mechanism for the remaining phases.

preprint2020arXiv

Reframing astronomical research through an anticolonial lens -- for TMT and beyond

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

preprint2020arXiv

Rotation, convective core overshooting, and period changes in classical Cepheid stellar evolution models

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.

preprint2020arXiv

The Curious Case of the North Star: the continuing tension between evolution models and measurements of Polaris

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.

preprint2019arXiv

Canadian Astronomy on Maunakea: On Respecting Indigenous Rights

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

preprint2018arXiv

Limb Darkening and Planetary Transits II: Intensity profile correction factors for a grid of model stellar atmospheres

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.

preprint2016arXiv

Period Changes and Evolution in Pulsating Variable Stars

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.

preprint2016arXiv

Stellar atmospheres, atmospheric extension and fundamental parameters: weighing stars using the stellar mass index

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.

preprint2016arXiv

The Secret Lives of Cepheids: Evolution, Mass Loss, and Ultraviolet Emission of the Long-Period Classical Cepheid $l$ Carinae

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 $δ$ 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 $δ$ Cep, and $β$ 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.

preprint2015arXiv

Period change and stellar evolution of $β$ Cephei stars

The $β$ 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 $β$ 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 $β$ 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 $β$ 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 $β$ Cephei stars: $ξ^1$ CMa and $δ$ 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 $β$ Cep stars can artificially lead to overestimated rates of secular period change.

preprint2015arXiv

The strange evolution of the Large Magellanic Cloud Cepheid OGLE-LMC-CEP1812

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.

preprint2014arXiv

Convection, granulation and period jitter in classical Cepheids

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.

preprint2014arXiv

Effects of stellar evolution and ionizing radiation on the environments of massive stars

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.

preprint2014arXiv

Evidence of a Mira-like tail and bow shock about the semi-regular variable V CVn from four decades of polarization measurements

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.

preprint2014arXiv

Interacting supernovae from photoionization-confined shells around red supergiant stars

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.

preprint2014arXiv

Pulsation Period Change & Classical Cepheids: Probing the Details of Stellar Evolution

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.

preprint2014arXiv

Revisiting the fundamental properties of Cepheid Polaris using detailed stellar evolution models

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

preprint2014arXiv

The Circumstellar Medium of Massive Stars in Motion

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.

preprint2014arXiv

The occurrence of classical Cepheids in binary systems

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

preprint2014arXiv

The Secret Lives of Cepheids: Evolutionary Changes and Pulsation-Induced Shock Heating in the Prototype Classical Cepheid δ Cep

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.

preprint2014arXiv

V444 Cyg X-ray and polarimetric variability: Radiative and Coriolis forces shape the wind collision region

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.

preprint2013arXiv

Indicators of Mass in Spherical Stellar Atmospheres

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.

preprint2013arXiv

Long-term polarization observations of Mira variable stars suggest asymmetric structures

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.

preprint2013arXiv

Numerical models for the circumstellar medium around Betelgeuse

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.

preprint2013arXiv

Pulsation and Mass Loss Across the HR Diagram: From OB stars to Cepheids to Red Supergiants

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.

preprint2013arXiv

Spherically symmetric model stellar atmospheres and limb darkening II: limb-darkening laws, gravity-darkening coefficients and angular diameter corrections for FGK dwarf stars

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.

preprint2013arXiv

Spherically-symmetric model stellar atmospheres and limb darkening I: limb-darkening laws, gravity-darkening coefficients and angular diameter corrections for red giant stars

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.

preprint2012arXiv

Calibrating the projection factor for Galactic Cepheids

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.

preprint2012arXiv

Cepheid limb darkening, angular diameter corrections, and projection factor from static spherical model stellar atmospheres

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.

preprint2012arXiv

Classical Cepheids Require Enhanced Mass Loss

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.

preprint2012arXiv

Constructing a Cepheid Period p-Factor Relation from Static Model Stellar Atmospheres

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.

preprint2012arXiv

Double bow shocks around young, runaway red supergiants: application to Betelgeuse

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.

preprint2012arXiv

The period change of the Cepheid Polaris suggests enhanced mass loss

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.

preprint2012arXiv

Using limb darkening to measure fundamental parameters of stars

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.

preprint2011arXiv

Is there a mass discrepancy in the Cepheid binary OGLE-LMC-CEP0227?

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.

preprint2011arXiv

Limb darkening in spherical stellar atmospheres

(Abridged) Context. Stellar limb darkening, I(μ = cosθ), 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(μ) from plane- parallel models, have found that fits to the flux-normalized curves pass through a fixed point, a common μ location on the stellar disk, for all values of T eff, log g and wavelength. Aims. We study this fixed μ-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 μ-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 μ-value similar to the parametrized curves. We also find that the parametrized curves have two fixed points, μ1 and μ2, although μ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.

preprint2011arXiv

The Cepheid mass discrepancy and pulsation-driven mass loss

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.

preprint2010arXiv

Convective Core Overshoot & Mass Loss in Classical Cepheids: A Solution to the Mass Discrepancy?

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.

preprint2010arXiv

Testing Mass Loss in Large Magellanic Cloud Cepheids using Infrared and Optical Observations II. Predictions and Tests of the OGLE-III Fundamental-Mode Cepheids

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.

preprint2010arXiv

Using Limb-Darkening to Measure the Masses of Red Giants

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.

preprint2008arXiv

On the Enhancement of Mass Loss in Cepheids Due to Radial Pulsation. II. The Effect of 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.

preprint2008arXiv

Testing Mass Loss in Large Magellanic Cloud Cepheids using Infrared and Optical Observations

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.