Abstract available through American Astronomical Society.

The Chandra and XMM-Newton X-ray telescopes have led to numerous advances in the study and understanding of astrophysical X-ray sources. Particularly important has been the much increased spectral resolution of modern X-ray instrumentation. Wind-broadened emission lines have been spectroscopically resolved for many massive stars. This contribution reviews approaches to the modeling of X-ray emission line profile shapes from single stars, including smooth winds, winds with clumping, optically thin versus thick lines, and the effect of a radius-dependent photoabsorption coefficient.

Theoretical models predict that the compressed interstellar medium around runaway O stars can produce high-energy non-thermal diffuse emission, in particular, non-thermal X-ray and γ-ray emission. So far, detection of non-thermal X-ray emission was claimed for only one runaway star, AE Aur. We present a search for non-thermal diffuse X-ray emission from bow shocks using archived XMM-Newton observations for a clean sample of six well-determined runaway O stars. We find that none of these objects present diffuse X-ray emission associated with their bow shocks, similarly to previous X-ray studies toward ζ Oph and BD+43°3654. We carefully investigated multi-wavelength observations of AE Aur and could not confirm previous findings of non-thermal X-rays. We conclude that so far there is no clear evidence of non-thermal extended emission in bow shocks around runaway O stars.

The β Cephei stars represent an important class of massive star pulsators that probe 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, including 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 the theoretical rates of period change, that 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 low to infer any information from stellar evolution models, whereas for three other stars the rates are too high. We infer stellar parameters, such as mass and age, for two β Cephei stars: ξ1 CMa and δ Cet, which agree well with independent measurements. We explore ideas for why models may not predict the higher rates of period change. In particular, period drifts in β Cep stars can artificially lead to overestimated rates of secular period change.

The intriguing WN4b star WR 6 has been known to display epoch-dependent spectroscopic, photometric and polarimetric variability for several decades. In this paper, we set out to verify if a simplified analytical model in which corotating interaction regions (CIRs) threading an otherwise spherical wind is able to reproduce the many broad-band continuum light curves from the literature with a reasonable set of parameters. We modified the optically thin model developed by Ignace, St-Louis & Proulx-Giraldeau to approximately account for multiple scattering and used it to fit 13 separate data sets of this star. By including two CIRs in the wind, we obtained reasonable fits for all data sets with coherent values for the inclination of the rotation axis (i0 = 166°) and for its orientation in the plane of the sky, although in the latter case we obtained two equally acceptable values (ψ = 63° and 152°) from the polarimetry. Additional line profile variation simulations using the Sobolev approximation for the line transfer allowed us to eliminate the ψ = 152° solution. With the adopted configuration (i0 = 166° and ψ = 63°), we were able to reproduce all data sets relatively well with two CIRs located near the stellar equator and always separated by ∼90° in longitude. The epoch dependence comes from the fact that these CIRs migrate along the surface of the star. Density contrasts smaller than a factor of 2 and large opening angles for the CIR (β⪆35∘) were found to best reproduce the type of spectroscopic variability reported in the literature.

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-shock nebulae 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 ionization 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 for 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 implications 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

With a deep Chandra/HETGS exposure of WR 6, we have resolved emission lines whose profiles show that the X-rays originate from a uniformly expanding spherical wind of high X-ray-continuum optical depth. The presence of strong helium-like forbidden lines places the source of X-ray emission at tens to hundreds of stellar radii from the photosphere. Variability was present in X-rays and simultaneous optical photometry, but neither were correlated with the known period of the system or with each other. An enhanced abundance of sodium revealed nuclear-processed material, a quantity related to the evolutionary state of the star. The characterization of the extent and nature of the hot plasma in WR 6 will help to pave the way to a more fundamental theoretical understanding of the winds and evolution of massive stars.

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, James Webb Space Telescope (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.

Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations toward the runaway stars ζ Oph by Chandra and Suzaku and of BD+43°3654 by XMM-Newton to search for the presence of non-thermal X-ray emission. We found no evidence of non-thermal emission spatially coincident with the bow shocks; nonetheless, diffuse emission was detected in the vicinity of ζ Oph. After a careful analysis of its spectral characteristics, we conclude that this emission has a thermal nature with a plasma temperature of T ≈ 2 × 106 K. The cometary shape of this emission seems to be in line with recent predictions of radiation-hydrodynamic models of runaway stars. The case of BD+43°3654 is puzzling, as non-thermal emission has been reported in a previous work for this source.

To this date ψ Per is the only classical Be star that was angularly resolved in radio (by the VLA at λ = 2 cm). Gaussian fit to the azimuthally averaged visibility data indicates a disk size (FWHM) of ~500 stellar radii (Dougherty & Taylor 1992). Recently, we obtained new multi-band cm flux density measurements of ψ Per from the enhanced VLA. We modeled the observed spectral energy distribution (SED) covering the interval from ultraviolet to radio using the Monte Carlo radiative transfer code HDUST (Carciofi & Bjorkman 2006). An SED turndown, that occurs between far-IR and radio wavelengths, is explained by a truncated viscous decretion disk (VDD), although the shallow slope of the radio SED suggests that the disk is not simply cut off, as is assumed in our model. The best-fit size of a truncated disk derived from the modeling of the radio SED is 100+5 −15 stellar radii, which is in striking contrast with the result of Dougherty & Taylor (1992). The reasons for this discrepancy are under investigation.