The structured winds of single massive stars can be classified into two broad groups: stochastic structure and organized structure. While the former is typically identified with clumping, the latter is typically associated with rotational modulations, particularly the paradigm of Co-rotating Interaction Regions (CIRs). While CIRs have been explored extensively in the UV band, and moderately in the X-ray and optical, here we evaluate radio variability from CIR structures assuming free-free opacity in a dense wind. Our goal is to conduct a broad parameter study to assess the observational feasibility, and to this end, we adopt a phenomenological model for a CIR that threads an otherwise spherical wind. We find that under reasonable assumptions, it is possible to obtain radio variability at the 10% level. The detailed structure of the folded light curve depends not only on the curvature of the CIR, the density contrast of the CIR relative to the wind, and viewing inclination, but also on wavelength. Comparing light curves at different wavelengths, we find that the amplitude can change, that there can be phase shifts in the waveform, and the the entire waveform itself can change. These characterstics could be exploited to detect the presence of CIRs in dense, hot winds.

Here, we study the microlensing of radially pulsating stars. Discerning and characterizing the properties of distant, faint pulsating stars is achievable through high-cadence microlensing observations. Combining stellar variability period with microlensing gives the source distance, type, and radius and helps better determine the lens parameters. Considering periodically variations in their radius and surface temperature, their microlensing light curves are resulted from multiplication of the magnification factor with variable finite size effect by the intrinsic brightness curves of pulsing source. The variable finite source size due to pulsation can be significant for transit and single microlensing and while caustic-crossing features. This kind of deviation in the magnification factor is considerable when the ratio of the source radius to the projected lens-source distance is in the range of $\rho_{\star}/u \in[0.4,10]$ and its duration is short and in the same order of the time of crossing the source radius. Other deviations due to variable source intensity and its area make colored and periodic deviations which are asymmetric with respect to the signs of pulsation phase. The positive phases makes deviations with larger amplitude that negative phase. These deviations dominate in filters with short wave lengths (e.g., $B-$band). The position of magnification peaks in microlensing of variable stars varies and this displacement differs in different filters.

We have obtained a very deep exposure (813 ks) of $\zeta\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 \r{A} region, especially well suited for Chandra, which has a significant contribution from continuum emission between well separated emission lines from high-ionization species. These data allow us to study the hottest plasma present through the continuum shape and emission line strengths. Assuming a powerlaw emission measure distribution which has a high-temperature cut-off, we find that the emission is consistent with a thermal spectrum having a maximum temperature of 12 MK. This implies an effective wind shock velocity of $900\,\mathrm{km\,s^{-1}}$, well below the wind terminal speed of $2250\,\mathrm{km\,s^{-1}}$. For X-ray emission which forms close to the star, the speed and X-ray flux are larger than can be easily reconciled with strictly self-excited line-deshadowing-instability models, suggesting a need for a fraction of the wind to be accelerated extremely rapidly right from the base. This is not so much a dynamical instability as a nonlinear response to changing boundary conditions.

We present new ATCA multi-wavelength radio measurements (range 2.1-21.2 GHz) of the early-type magnetic star rho Oph A, performed in March 2019 during 3 different observing sessions. These new ATCA observations evidence a clear rotational modulation of the stellar radio emission and the detection of coherent auroral radio emission from rho Oph A at 2.1 GHz. We collected high-resolution optical spectra of rho Oph A acquired by several instruments over a time span of about ten years. We also report new magnetic field measurements of rho Oph A that, together with the radio light curves and the temporal variation of the equivalent width of the HeI line (lambda=5015 Angstrom), were used to constrain the rotation period and the stellar magnetic field geometry. The above results have been used to model the stellar radio emission, modelling that allowed us to constrain the physical condition of rho Oph A's magnetosphere. Past XMM measurements showed periodic X-ray pulses from rho Oph A. We correlate the X-ray light curve with the magnetic field geometry of rho Oph A. The already published XMM data have been re-analyzed showing that the X-ray spectra of rho Oph A are compatible with the presence of a non-thermal X-ray component. We discuss a scenario where the emission phenomena occurring at the extremes of the electromagnetic spectrum, radio and X-ray, are directly induced by the same plasma process. We interpret the observed X-ray and radio features of rho Oph A as having an auroral origin.

Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have been confirmed due to tight observational constraints. An indirect method of detecting otherwise hidden companions is offered by their effect on the outer parts of Be star disks, which are expected to be disrupted or truncated. In the context of the IR and radio continuum excess radiation originating in the disk, the disk truncation can be revealed by a turndown in the spectral energy distribution due to reduced radio flux levels. In this work we search for signs of spectral turndown in a sample of 57 classical Be stars with radio data, which include new data for 23 stars and the longest wavelength detections so far ($\lambda \approx$ 10\,cm) for 2 stars. We confidently detect the turndown for all 26 stars with sufficient data coverage (20 of which are not known to have close binary companions). For the remaining 31 stars, data are inconclusive as to whether the turndown is present or not. The analysis suggests that many if not all Be stars have close companions influencing their outer disks. If confirmed to be subdwarf companions, the mass transfer spin-up scenario might explain the existence of the vast majority of classical Be stars.

PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.

Among Wolf-Rayet stars, those of subtype WN8 are the intrinsically most variable. We have explored the long-term photometric variability of the brightest known WN8 star, WR 40, through four contiguous months of time-resolved, single-passband optical photometry with the BRIght Target Explorer (BRITE) nanosatellite mission. The Fourier transform of the observed light-curve reveals that the strong light variability exhibited by WR 40 is dominated by many randomly-triggered, transient, low-frequency signals. We establish a model in which the whole wind consists of stochastic clumps following an outflow visibility promptly rising to peak brightness upon clump emergence from the optically thick pseudo-photosphere in the wind, followed by a gradual decay according to the right-half of a Gaussian. Free electrons in each clump scatter continuum light from the star. We explore a scenario where the clump size follows a power-law distribution, and another one with an ensemble of clumps of constant size. Both scenarios yield simulated light curves morphologically resembling the observed light curve remarkably well, indicating that one cannot uniquely constrain the details of clump size distribution with only a photometric light curve. Nevertheless, independent evidence favours a negative-index power law, as seen in many other astrophysical turbulent media.

We use a 3D Monte Carlo radiative transfer code to study the polarimetric and photometric variability from stationary corotating interaction regions (CIR) in the wind of massive stars. Our CIRs are approximated by Archimedean spirals of higher (or lower) density formed in a spherical wind originating from the star and we also made allowance for a bright gaussian spot at the base of the CIR. Comparing results from our code to previous analytical calculations in the optically thin case, we find differences which we attribute mainly to a better estimation of the total unpolarized flux reaching the observer. In the optically thick case, the differences with the analytical calculations are much larger, as multiple scattering introduces extra complexities including occultation effects. The addition of a gaussian spot does not alter the shape of the polarization curve significantly but does create a small excess in polarization. On the other hand, the effect can be larger on the light curve and can become dominant over the resulting CIR, depending on the spot parameters and density of the wind.

High spectral resolution and long exposure times are providing unprecedented levels of data quality of massive stars at X-ray wavelengths. A key diagnostic of the X-ray emitting plasma are the fir lines for He-like triplets. In particular, owing to radiative pumping effects, the forbidden-to-intercombination line luminosity ratio, R=f/i, can be used to determine the proximity of the hot plasma to the UV-bright photospheres of massive stars. Moreover, the era of large observing programs additionally allows for investigation of line variability. This contribution is the second to explore how variability in the line ratio can provide new diagnostic information about distributed X-rays in a massive star wind. While there are many ways to drive variability in the line ratio, we use variable mass loss as an illustrative example. The f/i ratio can be significantly modulated owing to evolving wind properties. We evaluate how variable mass loss might bias measures of f/i.

We present new, contemporaneous HST STIS and XMM observations of the O7 III(n)((f)) star xi Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of xi Pper. The variable wind of this star is known to have a 2.086 day periodicity. We use a simple, heuristic spot model which fits the low velocity (near surface) IUE wind line variability very well, to demonstrate that the low velocity absorption in the new STIS spectra of N IV 1718 and Si IV 1402 vary with the same 2.086 day period. It is remarkable that the period and amplitude of the STIS data agree with those of the IUE spectra obtained 22 years earlier. We also show that the time variability of the new XMM fluxes are also consistent with the 2.086 day period. Thus, our new, multi-wavelength coordinated observations demonstrate that the mechanism which causes the UV wind line variability is also responsible for a significant fraction of the X-rays in single O stars. The sequence of events for the multi-wavelength light curve minima is: Si IV 1402, N IV 1718, and X-ray flux, each separated by a phase of about 0.06 relative to the 2.086 day period. Analysis of the X-ray fluxes shows that they become softer as they weaken. This is contrary to expectations if the variability is caused by periodic excess absorption. Further, the high resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blue shifted. If we interpret the low velocity wind line light curves in terms of our model, it implies that there are two bright regions, i.e., regions with less absorption, separated by 180 deg, on the surface of the star. We note that the presence and persistent of two spots separated by 180 deg suggests that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in xi Per.

Among different types of massive stars in advanced evolutionary stages is the enigmatic WN8h type. There are only a few Wolf-Rayet (WR) stars with this spectral type in our Galaxy. It has long been suggested that WN8h-type stars are the products of binary evolution that may harbor neutron stars (NS). One of the most intriguing WN8h stars is the runaway WR124 surrounded by its magnificent nebula M1-67. We test the presence of an accreting NS companion in WR124 using 100 ks long observations by the Chandra X-ray observatory. The hard X-ray emission from WR124 with a luminosity of $L_\mathrm{X}$= 10$^{31}$ erg s$^{-1}$ is marginally detected. We use the non-LTE stellar atmosphere code PoWR to estimate the WR wind opacity to the X-rays. The wind of a WN8-type star is effectively opaque for X-rays, hence the low X-ray luminosity of WR124 does not rule out the presence of an embedded compact object. We suggest that, in general, high opacity WR winds could prevent X-ray detections of embedded NS, and be an explanation for the apparent lack of WR+NS systems.

Recombination line profile shapes are derived for ionized spherical stellar winds at radio wavelengths. It is assumed that the wind is optically thick owing to free-free opacity. Emission lines of arbitrary optical depth are obtained assuming that the free-free photosphere forms in the outer, constant expansion portion of the wind. Previous works have derived analytic results for isothermal winds when the line and continuum source functions are equal. Here, semi-analytic results are derived for when the source functions are not equal to reveal that line shapes can be asymmetric about line center. A parameter study is presented and applications discussed.

We report the detection of the auroral radio emission from the early-type magnetic star HD142301. New VLA observations of HD142301 detected highly polarized amplified emission occurring at fixed stellar orientations. The coherent emission mechanism responsible for the stellar auroral radio emission amplifies the radiation within a narrow beam, making the star where this phenomenon occurs similar to a radio lighthouse. The elementary emission process responsible for the auroral radiation mainly amplifies one of the two magneto-ionic modes of the electromagnetic wave. This explains why the auroral pulses are highly circularly polarized. The auroral radio emission of HD142301 is characterized by a reversal of the sense of polarization as the star rotates. The effective magnetic field curve of HD142301 is also available making it possible to correlate the transition from the left to the right-hand circular polarization sense (and vice-versa) of the auroral pulses with the known orientation of the stellar magnetic field. The results presented in this letter have implications for the estimation of the dominant magneto-ionic mode amplified within the HD142301 magnetosphere.

We examine over 160 archival H$\alpha$ spectra from the Ritter Observatory for the interacting binary $\beta$~Lyr obtained between 1996 and 2000. The emission is characteristically double-peaked, but asymmetric, and with an absorption feature that is persistently blueshifted. Using a set of simplifying assumptions, phase varying emission line profiles are calculated for H$\alpha$ formed entirely in a Keplerian disk, and separately for the line formed entirely from an off-center bipolar flow. However, a dynamic spectrum of the data indicate the blueshifted feature is not always present, and the data are even suggestive of a drift of the feature in velocity shift. We explore whether a circumbinary envelope, hot spot on the accretion disk, or accretion stream could explain the observations. While none are satisfactory, an accretion stream explanation is somewhat promising.

Oxygen sequence Wolf-Rayet stars (WO) are thought to be the final evolution phase of some high mass stars, as such they may be the progenitors of type Ic SNe as well as potential progenitors of broad-lined Ic and long gamma-ray bursts. We present the first spectropolarimetric observations of the Galactic WO stars WR93b and WR102 obtained with FORS1 on the VLT. We find no sign of a line effect, which could be expected if these stars were rapid rotators. We also place constraints on the amplitude of a potentially undetected line effect. This allows us to derive upper limits on the possible intrinsic continuum polarisation, and find P$_{\rm cont}$ < 0.077 percent and P$_{\rm cont}$ < 0.057 percent for WR93b and WR102, respectively. Furthermore, we derive upper limits on the rotation of our WO stars by considering our results in the context of the wind compression effect. We estimate that for an edge-on case the rotational velocity of WR93b is v$_{\rm rot}$ < 324 km/s while for WR102 v$_{\rm rot}$ < 234 km/s. These correspond to values of v$_{\rm rot}$/v$_{\rm crit}$ <19 percent and <10 percent, respectively, and values of log(j)<18.0 cm$^2$/s for WR93b and <17.6 cm^2 /s for WR102. The upper limits found on v$_{\rm rot}$/v$_{\rm crit}$ and log(j) for our WO stars are therefore similar to the estimates calculated for Galactic WR stars that do show a line effect. Therefore, although the presence of a line effect in single WR stars is indicative of fast rotation, the absence of a line effect does not rule out significant rotation, even when considering the edge-on scenario.

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.

We present new radio/millimeter measurements of the hot magnetic star HR5907 obtained with the VLA and ALMA interferometers. We find that HR5907 is the most radio luminous early type star in the cm-mm band among those presently known. Its multi-wavelength radio light curves are strongly variable with an amplitude that increases with radio frequency. The radio emission can be explained by the populations of the non-thermal electrons accelerated in the current sheets on the outer border of the magnetosphere of this fast rotating magnetic star. We classify HR5907 as another member of the growing class of strongly magnetic fast rotating hot stars where the gyro-synchrotron emission mechanism efficiently operates in their magnetospheres. The new radio observations of HR5907 are combined with archival X-ray data to study the physical condition of its magnetosphere. The X-ray spectra of HR5907 show tentative evidence for the presence of non-thermal spectral component. We suggest that non-thermal X-rays originate a stellar X-ray aurora due to streams of non-thermal electrons impacting on the stellar surface. Taking advantage of the relation between the spectral indices of the X-ray power-law spectrum and the non-thermal electron energy distributions, we perform 3-D modeling of the radio emission for HR5907. The wavelength-dependent radio light-curves probe magnetospheric layers at different heights above the stellar surface. A detailed comparison between simulated and observed radio light-curves leads us to conclude that the stellar magnetic field of HR5907 is likely non-dipolar, providing further indirect evidence of the complex magnetic field topology of HR5907.

The intriguing WN4b star WR6 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 broadband continuum light curves from the literature with a reasonable set of parameters. We modified the optically thin model we developed in Ignace, St-Louis & Proulx-Giraldeau (2015) to approximately account for multiple scattering and used it to fit 13 separate datasets of this star. By including two CIRs in the wind, we obtained reasonable fits for all datasets with coherent values for the inclination of the rotation axis ($i_0=166^{\circ}$) and for its orientation in the plane of the sky, although in the latter case we obtained two equally acceptable values ($\psi=63^{\circ}$ and $\psi=152^{\circ}$) from the polarimetry. Additional line profile variation simulations using the Sobolev approximation for the line transfer allowed us to eliminate the $\psi=152^{\circ}$ solution. With the adopted configuration ($i_0=166^{\circ}$ and $\psi=63^{\circ}$), we were able to reproduce all datasets relatively well with two CIRs located near the stellar equator and always separated by $\sim 90^{\circ}$ in longitude. The epoch-dependency comes from the fact that these CIRs migrate along the surface of the star. Density contrasts smaller than a factor of two and large opening angles for the CIR ($\beta \geq 35^{\circ}$) were found to best reproduce the type of spectroscopic variability reported in the literature.

From $5.5$ months of dual-band optical photometric monitoring at the $1$ mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star $\zeta$ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the $1.78$ d signal previously detected by Coriolis/SMEI, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light curve inversion algorithm we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II $\lambda4686$ line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in $\zeta$ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.

AH Cep and CW Cep are both early B-type binaries with short orbital periods of 1.8~d and 2.7~d, respectively. All four components are B0.5V types. The binaries are also double-lined spectroscopic and eclipsing. Consequently, solutions for orbital and stellar parameters make the pair of binaries ideal targets for a study of the colliding winds between two B~stars. {\em Chandra} ACIS-I observations were obtained to determine X-ray luminosities. AH~Cep was detected with an unabsorbed X-ray luminosity at a 90\% confidence interval of $(9-33)\times 10^{30}$ erg s$^{-1}$, or $(0.5-1.7)\times 10^{-7} L_{\rm Bol}$, relative to the combined Bolometric luminosities of the two components. While formally consistent with expectations for embedded wind shocks, or binary wind collision, the near-twin system of CW~Cep was a surprising non-detection. For CW~Cep, an upper limit was determined with $L_X/L_{\rm Bol} < 10^{-8}$, again for the combined components. One difference between these two systems is that AH~Cep is part of a multiple system. The X-rays from AH~Cep may not arise from standard wind shocks nor wind collision, but perhaps instead from magnetism in any one of the four components of the system. The possibility could be tested by searching for cyclic X-ray variability in AH~Cep on the short orbital period of the inner B~stars.

The blue hypergiant Cyg OB2-12 (B3Ia+) is a representative member of the class of very massive stars in a poorly understood evolutionary stage. We obtained its high-resolution X-ray spectrum using Chandra observatory. PoWR model atmospheres were calculated to provide realistic wind opacities and to establish the wind density structure. We find that collisional de-excitation is the dominant mechanism de-populating the metastable upper levels of the forbidden lines of the He-like ions SiXIV and MgXII. Comparison between the model and observations reveals that X-ray emission is produced in a dense plasma, which could reside only at the photosphere or in a colliding wind zone between binary components. The observed X-ray spectra are well fitted by thermal plasma models, with average temperatures in excess of 10 MK. The wind speed in Cyg OB2-12 is not high enough to power such high temperatures, but the collision of two winds in a binary system can be sufficient. We used archival data to investigate the X-ray properties of other blue hypergiants. In general, stars of this class are not detected as X-rays sources. We suggest that our new Chandra observations of Cyg OB2-12 can be best explained if Cyg OB2-12 is a colliding wind binary possessing a late O-type companion. This makes Cyg OB2-12 only the second binary system among the 16 known Galactic hypergiants. This low binary fraction indicates that the blue hypergiants are likely products of massive binary evolution during which they either accreted a significant amount of mass or already merged with their companion.

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.

Context. The structure of the inner parts of Be star disks (20 stellar radii) is well explained by the viscous decretion disk (VDD) model, which is able to reproduce the observable properties of most of the objects studied so far. The outer parts, on the ther hand, are not observationally well-explored, as they are observable only at radio wavelengths. A steepening of the spectral slope somewhere between infrared and radio wavelengths was reported for several Be stars that were previously detected in the radio, but a convincing physical explanation for this trend has not yet been provided. Aims. We test the VDD model predictions for the extended parts of a sample of six Be disks that have been observed in the radio to address the question of whether the observed turndown in the spectral energy distribution (SED) can be explained in the framework of the VDD model, including recent theoretical development for truncated Be disks in binary systems. Methods. We combine new multi-wavelength radio observations from the Karl. G. Jansky Very Large Array (JVLA) and Atacama Pathfinder Experiment (APEX) with previously published radio data and archival SED measurements at ultraviolet, visual, and infrared wavelengths. The density structure of the disks, including their outer parts, is constrained by radiative transfer modeling of the observed spectrum using VDD model predictions. In the VDD model we include the presumed effects of possible tidal influence from faint binary companions. Results. For 5 out of 6 studied stars, the observed SED shows strong signs of SED turndown between far-IR and radio wavelengths. A VDD model that extends to large distances closely reproduces the observed SEDs up to far IR wavelengths, but fails to reproduce the radio SED. ... (abstract continues but did not fit here)

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 $\gamma$-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 6 well-determined runaway O stars. We find that none of these objects present diffuse X-ray emission associated to their bow shocks, similarly to previous X-ray studies toward $\zeta$ Oph and BD$+$43$^{\circ}$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.

Observations with powerful X-ray telescopes, such as XMM-Newton and Chandra, significantly advance our understanding of massive stars. Nearly all early-type stars are X-ray sources. Studies of their X-ray emission provide important diagnostics of stellar winds. High-resolution X-ray spectra of O-type stars are well explained when stellar wind clumping is taking into account, providing further support to a modern picture of stellar winds as non-stationary, inhomogeneous outflows. X-ray variability is detected from such winds, on time scales likely associated with stellar rotation. High-resolution X-ray spectroscopy indicates that the winds of late O-type stars are predominantly in a hot phase. Consequently, X-rays provide the best observational window to study these winds. X-ray spectroscopy of evolved, Wolf-Rayet type, stars allows to probe their powerful metal enhanced winds, while the mechanisms responsible for the X-ray emission of these stars are not yet understood.

In this paper we investigate the multiwavelengths properties of the magnetic early B-type star HR7355. We present its radio light curves at several frequencies, taken with the Jansky Very Large Array, and X-ray spectra, taken with the XMM X-ray telescope. Modeling of the radio light curves for the Stokes I and V provides a quantitative analysis of the HR7355 magnetosphere. A comparison between HR7355 and a similar analysis for the Ap star CUVir, allows us to study how the different physical parameters of the two stars affect the structure of the respective magnetospheres where the non-thermal electrons originate. Our model includes a cold thermal plasma component that accumulates at high magnetic latitudes that influences the radio regime, but does not give rise to X-ray emission. Instead, the thermal X-ray emission arises from shocks generated by wind stream collisions close to the magnetic equatorial plane. The analysis of the X-ray spectrum of HR7355 also suggests the presence of a non-thermal radiation. Comparison between the spectral index of the power-law X-ray energy distribution with the non-thermal electron energy distribution indicates that the non-thermal X-ray component could be the auroral signature of the non-thermal electrons that impact the stellar surface, the same non-thermal electrons that are responsible for the observed radio emission. On the basis of our analysis, we suggest a novel model that simultaneously explains the X-ray and the radio features of HR7355 and is likely relevant for magnetospheres of other magnetic early type stars.

To investigate if magnetic fields are present in Wolf-Rayet stars, we selected a few stars in the Galaxy and one in the Large Magellanic Cloud (LMC). We acquired low-resolution spectropolarimetric observations with the ESO FORS2 instrument during two different observing runs. During the first run in visitor mode, we observed the LMC Wolf-Rayet star BAT99 7 and the stars WR6, WR7, WR18, and WR23 in our Galaxy. The second run in service mode was focused on monitoring the star WR6. Linear polarization was recorded immediately after the observations of circular polarization. During our visitor observing run, the magnetic field for the cyclically variable star WR6 was measured at a significance level of 3.3sigma (<B_z> = 258+-78G). Among the other targets, the highest value for the longitudinal magnetic field, <B_z> = 327+-141G, was measured in the LMC star BAT99 7. Spectropolarimetric monitoring of the star WR6 revealed a sinusoidal nature of the <B_z> variations with the known rotation period of 3.77d, significantly adding to the confidence in the detection. The presence of the rotation-modulated magnetic variability is also indicated in our frequency periodogram. The reported field magnitude suffers from significant systematic uncertainties at the factor 2 level, in addition to the quoted statistical uncertainties, owing to the theoretical approach used to characterize it. Linear polarization measurements showed no line effect in the stars, apart from WR6. BAT99 7, WR7, and WR,23 do not show variability of the linear polarization over two nights.

Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations towards the runaway stars $\zeta$ Oph (Chandra and Suzaku) and BD+43$^{\circ}$3654 (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 is detected in the vicinity of $\zeta$ Oph. After a careful analysis of its spectral characteristics we conclude that this emission has a thermal nature with a plasma temperature of $T \approx 2 \times10^{6}$ 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$^{\circ}$3654 is puzzling as non-thermal emission has been reported in a previous work for this source.

The influence of macroclumps for free-free spectral energy distributions (SEDs) of ionized winds is considered. The goal is to emphasize distinctions between microclumping and macroclumping effects. Microclumping can alter SED slopes and flux levels if the volume filling factor of the clumps varies with radius; however, the modifications are independent of the clump geometry. To what extent does macroclumping alter SED slopes and flux levels? In addressing the question, two specific types of macroclump geometries are explored: shell fragments ("pancake"-shaped) and spherical clumps. Analytic and semi-analytic results are derived in the limiting case that clumps never obscure one another. Numerical calculations based on a porosity formalism is used when clumps do overlap. Under the assumptions of a constant expansion, isothermal, and fixed ionization wind, the fragment model leads to results that are essentially identical to the microclumping result. Mass-loss rate determinations are not affected by porosity effects for shell fragments. By contrast, spherical clumps can lead to a reduction in long-wavelength fluxes, but the reductions are only significant for extreme volume filling factors.

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.

The MiMeS project is a large-scale, high resolution, sensitive spectropolarimetric investigation of the magnetic properties of O and early B type stars. Initiated in 2008 and completed in 2013, the project was supported by 3 Large Program allocations, as well as various programs initiated by independent PIs and archival resources. Ultimately, over 4800 circularly polarized spectra of 560 O and B stars were collected with the instruments ESPaDOnS at the Canada-France-Hawaii Telescope, Narval at the T\'elescope Bernard Lyot, and HARPSpol at the European Southern Observatory La Silla 3.6m telescope, making MiMeS by far the largest systematic investigation of massive star magnetism ever undertaken. In this paper, the first in a series reporting the general results of the survey, we introduce the scientific motivation and goals, describe the sample of targets, review the instrumentation and observational techniques used, explain the exposure time calculation designed to provide sensitivity to surface dipole fields above approximately 100 G, discuss the polarimetric performance, stability and uncertainty of the instrumentation, and summarize the previous and forthcoming publications.

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.

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.

We present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, Delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary, Delta Ori A provides a unique system with which to test the spatial distribution of the X-ray emitting gas around Delta Ori Aa1 via occultation by the photosphere of and wind cavity around the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ksec and covering nearly the entire binary orbit. Companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities simultaneous with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectrum. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about $0.3-0.5\times$ the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of \ion{Fe}{17} and \ion{Ne}{10} are inconsistent with model predictions, which may be an effect of resonance scattering

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral dataset of the $\delta$ Orionis Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of ~479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 $\AA$ is confirmed, with maximum amplitude of about +/-15% within a single ~125 ks observation. Periods of 4.76d and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phase=0.0 when the secondary $\delta$ Orionis Aa2 is at inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.

Beta Cephei-type variables are early B-type stars that are characterized by oscillations observable in their optical light curves. At least one Beta Cep-variable also shows periodic variability in X-rays. Here we study the X-ray light curves in a sample of beta Cep-variables to investigate how common X-ray pulsations are for this type of stars. We searched the Chandra and XMM-Newton X-ray archives and selected stars that were observed by these telescopes for at least three optical pulsational periods. We retrieved and analyzed the X-ray data for kappa Sco, beta Cru, and alpha Vir. The X-ray light curves of these objects were studied to test for their variability and periodicity. While there is a weak indication for X-ray variability in beta Cru, we find no statistically significant evidence of X-ray pulsations in any of our sample stars. This might be due either to the insufficient data quality or to the physical lack of modulations. New, more sensitive observations should settle this question.

Massive star winds are complex radiation-hydrodynamic (sometimes magnetohydrodynamic) outflows that are propelled by their enormously strong luminosities. The winds are often found to be structured and variable, but can also display periodic or quasi-periodic behavior in a variety of wind diagnostics. The regular variations observed in putatively single stars, especially in UV wind lines, have often been attributed to corotating interaction regions (CIRs) like those seen in the solar wind. We present light curves for variable polarization from winds with CIR structures. We develop a model for a time-independent CIR based on a kinematical description. Assuming optically thin electron scattering, we explore the range of polarimetric light curves that result as the curvature, latitude, and number of CIRs are varied. We find that a diverse array of variable polarizations result from an exploration of cases. The net polarization from an unresolved source is weighted more toward the inner radii of the wind. Given that most massive stars have relatively fast winds compared to their rotation speeds, CIRs tend to be conical at inner radii, transitioning to a spiral shape at a few to several stellar radii in the wind. Winds with a single CIR structure lead to easily identifiable polarization signatures. By contrast allowing for multiple CIRs, all emerging from a range of azimuth and latitude positions at the star, can yield complex polarimetric behavior. Although our model is based on some simplifying assumptions, it produces qualitative behavior that we expect to be robust, and this has allowed us to explore a wide range of CIR configurations that will prove useful for interpreting polarimetric data.

A rather strong mean longitudinal magnetic field of the order of a few hundred Gauss was detected a few years ago in the Of?p star CPD -28 2561 using FORS2 low-resolution spectropolarimetric observations. In this work we present additional low-resolution spectropolarimetric observations obtained during several weeks in 2013 December using FORS2 (FOcal Reducer low dispersion Spectrograph) mounted at the 8-m Antu telescope of the VLT. These observations cover a little less than half of the stellar rotation period of 73.41d mentioned in the literature. The behaviour of the mean longitudinal magnetic field is consistent with the assumption of a single-wave variation during the stellar rotation cycle, indicating a dominant dipolar contribution to the magnetic field topology. The estimated polar strength of the surface dipole B_d is larger than 1.15kG. Further, we compared the behaviour of the line profiles of various elements at different rotation phases associated with different magnetic field strengths. The strongest contribution of the emission component is observed at the phases when the magnetic field shows a negative or positive extremum. The comparison of the spectral behaviour of CPD -28 2561 with that of another Of?p star, HD148937 of similar spectral type, reveals remarkable differences in the degree of variability between both stars. Finally, we present new X-ray observations obtained with the Suzaku X-ray Observatory. We report that the star is X-ray bright with log L_X/L_bol ~= -5.7. The low resolution X-ray spectra reveal the presence of a plasma heated up to 24MK. We associate the 24MK plasma in CPD -28 2561 with the presence of a kG strong magnetic field capable to confine stellar wind.

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.

X-ray emission from stars much more massive than the Sun was discovered only 35 years ago. Such stars drive fast stellar winds where shocks can develop, and it is commonly assumed that the X-rays emerge from the shock-heated plasma. Many massive stars additionally pulsate. However, hitherto it was neither theoretically predicted nor observed that these pulsations would affect their X-ray emission. Here we report the discovery of pulsating X-rays from the massive B-type star Xi1 Canis Majoris. This star is a variable of beta Cephei type and has a strong magnetic field. Our observations with the XMM-Newton telescope reveal X-ray pulsations with the same period as the fundamental stellar pulsation. This discovery challenges our understanding of stellar winds from massive stars, their X-ray emission, and their magnetism.

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.

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.

We obtained four pointings of over 100 ks each of the well-studied Wolf-Rayet star WR 6 with the XMM-Newton satellite. With a first paper emphasizing the results of spectral analysis, this follow-up highlights the X-ray variability clearly detected in all four pointings. However, phased light curves fail to confirm obvious cyclic behavior on the well-established 3.766 d period widely found at longer wavelengths. The data are of such quality that we were able to conduct a search for "event clustering" in the arrival times of X-ray photons. However, we fail to detect any such clustering. One possibility is that X-rays are generated in a stationary shock structure. In this context we favor a co-rotating interaction region (CIR) and present a phenomenological model for X-rays from a CIR structure. We show that a CIR has the potential to account simultaneously for the X-ray variability and constraints provided by the spectral analysis. Ultimately, the viability of the CIR model will require both intermittent long-term X-ray monitoring of WR 6 and better physical models of CIR X-ray production at large radii in stellar winds.

Faraday rotation is a valuable tool for detecting magnetic fields. Here the technique is considered in relation to wind-blow bubbles. In the context of spherical winds with azimuthal or split monopole stellar magnetic field geometries, we derive maps of the distribution of position angle (PA) rotation of linearly polarized radiation across projected bubbles. We show that the morphology of maps for split monopole fields are distinct from those produced by the toroidal field topology; however, the toroidal case is the one most likely to be detectable because of its slower decline in field strength with distance from the star. We also consider the important case of a bubble with a spherical sub-volume that is field-free to approximate crudely a "swept-up" wind interaction between a fast wind (or possibly a supernova ejecta shell) overtaking a slower magnetized wind from a prior state of stellar evolution. With an azimuthal field, the resultant PA map displays two arc-like features of opposite rotation measure, similar to observations of the supernova remnant G296.5+10.0. We illustrate how PA maps can be used to disentangle Faraday rotation contributions made by the interstellar medium versus the bubble. Although our models involve simplifying assumptions, their consideration leads to a number of general robust conclusions for use in the analysis of radio mapping datasets.

An increasing number of OB stars have been shown to possess magnetic fields. Although the sample remains small, it is surprising that the magnetic and X-ray properties of these stars appear to be far less correlated than expected. This contradicts model predictions, which generally indicate that the X-rays from magnetic stars to be harder and more luminous than their non-magnetic counterparts. Instead, the X-ray properties of magnetic OB stars are quite diverse. $\tau$ Sco is one example where the expectations are better met. This bright main sequence, early B star has been studied extensively in a variety of wavebands. It has a surface magnetic field of around 500 G, and Zeeman Doppler tomography has revealed an unusual field configuration. Furthermore, tau Sco displays an unusually hard X-ray spectrum, much harder than similar, non-magnetic OB stars. In addition, the profiles of its UV P Cygni wind lines have long been known to possess a peculiar morphology. Recently, two stars, HD 66665 and HD 63425, whose spectral types and UV wind line profiles are similar to those of $\tau$ Sco, have also been determined to be magnetic. In the hope of establishing a magnetic field - X-ray connection for at least a sub-set of the magnetic stars, we obtained XMM-Newton EPIC spectra of these two objects. Our results for HD 66665 are somewhat inconclusive. No especially strong hard component is detected; however, the number of source counts is insufficient to rule out hard emission. longer exposure is needed to assess the nature of the X-rays from this star. On the other hand, we do find that HD 63425 has a substantial hard X-ray component, thereby bolstering its close similarity to tau Sco.

We report the detection of a strong, organized magnetic field in the secondary component of the massive O8III/I+O7.5V/III double-lined spectroscopic binary system HD 47129 (Plaskett's star), in the context of the Magnetism in Massive Stars (MiMeS) survey. Eight independent Stokes $V$ observations were acquired using the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope and the Narval spectropolarimeter at the T\'elescope Bernard Lyot. Using Least-Squares Deconvolution we obtain definite detections of signal in Stokes $V$ in 3 observations. No significant signal is detected in the diagnostic null ($N$) spectra. The Zeeman signatures are broad and track the radial velocity of the secondary component; we therefore conclude that the rapidly-rotating secondary component is the magnetized star. Correcting the polarized spectra for the line and continuum of the (sharp-lined) primary, we measured the longitudinal magnetic field from each observation. The longitudinal field of the secondary is variable and exhibits extreme values of $-810\pm 150$ G and $+680\pm 190$ G, implying a minimum surface dipole polar strength of $2850\pm 500$ G. In contrast, we derive an upper limit ($3\sigma$) to the primary's surface magnetic field of 230 G. The combination of a strong magnetic field and rapid rotation leads us to conclude that the secondary hosts a centrifugal magnetosphere fed through a magnetically confined wind. We revisit the properties of the optical line profiles and X-ray emission - previously interpreted as a consequence of colliding stellar winds - in this context. We conclude that HD 47129 represents a heretofore unique stellar system - a close, massive binary with a rapidly rotating, magnetized component - that will be a rich target for further study.

\mu\ Columbae is a prototypical weak-wind O-star for which we have obtained a high-resolution X-ray spectrum with the Chandra LETG/ACIS-S instrument and a low resolution spectrum with Suzaku. This allows us, for the first time, to investigate the role of X-rays on the wind structure in a bona fide weak-wind system and to determine whether there actually is a massive, hot wind. The X-ray emission measure indicates that the outflow is an order of magnitude greater than that derived from UV lines and is commensurate with the nominal wind-luminosity relationship for O-stars. Therefore, the ``weak-wind problem''---identified from cool wind UV/optical spectra---is largely resolved by accounting for the hot wind seen in X-rays. From X-ray line profiles, Doppler shifts, and relative strengths, we find that this weak-wind star is typical of other late O dwarfs. The X-ray spectra do not suggest a magnetically confined plasma---the spectrum is soft and lines are broadened; Suzaku spectra confirm the lack of emission above 2 keV. Nor do the relative line shifts and widths suggest any wind decoupling by ions. The He-like triplets indicate that the bulk of the X-ray emission is formed rather close to the star, within 5 stellar radii. Our results challenge the idea that some OB stars are ``weak-wind'' stars that deviate from the standard wind-luminosity relationship. The wind is not weak, but it is hot and its bulk is only detectable in X-rays.

The O9IV star HD 57682, discovered to be magnetic within the context of the MiMeS survey in 2009, is one of only eight convincingly detected magnetic O-type stars. Among this select group, it stands out due to its sharp-lined photospheric spectrum. Since its discovery, the MiMeS Collaboration has continued to obtain spectroscopic and magnetic observations in order to refine our knowledge of its magnetic field strength and geometry, rotational period, and spectral properties and variability. In this paper we report new ESPaDOnS spectropolarimetric observations of HD 57682, which are combined with previously published ESPaDOnS data and archival H{\alpha} spectroscopy. This dataset is used to determine the rotational period (63.5708 \pm 0.0057 d), refine the longitudinal magnetic field variation and magnetic geometry (dipole surface field strength of 880\pm50 G and magnetic obliquity of 79\pm4\circ as measured from the magnetic longitudinal field variations, assuming an inclination of 60\circ), and examine the phase variation of various lines. In particular, we demonstrate that the H{\alpha} equivalent width undergoes a double-wave variation during a single rotation of the star, consistent with the derived magnetic geometry. We group the variable lines into two classes: those that, like H{\alpha}, exhibit non-sinusoidal variability, often with multiple maxima during the rotation cycle, and those that vary essentially sinusoidally. Based on our modelling of the H{\alpha} emission, we show that the variability is consistent with emission being generated from an optically thick, flattened distribution of magnetically-confined plasma that is roughly distributed about the magnetic equator. Finally, we discuss our findings in the magnetospheric framework proposed in our earlier study.

We present the first high-resolution X-ray spectrum of a putatively single Wolf-Rayet star. 400 ks observations of WR 6 by the XMM-Newton-telescope resulted in a superb quality high-resolution X-ray spectrum. Spectral analysis reveals that the X-rays originate far out in the stellar wind, more than 30 stellar radii from the photosphere, and thus outside the wind acceleration zone where the line-driving instability could create shocks. The X-ray emitting plasma reaches temperatures up to 50\,MK, and is embedded within the un-shocked, "cool" stellar wind as revealed by characteristic spectral signatures. We detect a fluorescent Fe line at approx 6.4 keV. The presence of fluorescence is consistent with a two-component medium, where the cool wind is permeated with the hot X-ray emitting plasma. The wind must have a very porous structure to allow the observed amount of X-rays to escape. We find that neither the line-driving instability nor any alternative binary scenario can explain the data. We suggest a scenario where X-rays are produced when the fast wind rams into slow "sticky clumps" that resist acceleration. Our new data show that the X-rays in single WR-star are generated by some special mechanism different from the one operating in the O-star winds.

The consequences of structured flows continue to be a pressing topic in relating spectral data to physical processes occurring in massive star winds. In a preceding paper, our group reported on hydrodynamic simulations of hypersonic flow past a rigid spherical clump to explore the structure of bow shocks that can form around wind clumps. Here we report on profiles of emission lines that arise from such bow shock morphologies. To compute emission line profiles, we adopt a two component flow structure of wind and clumps using two "beta" velocity laws. While individual bow shocks tend to generate double horned emission line profiles, a group of bow shocks can lead to line profiles with a range of shapes with blueshifted peak emission that depends on the degree of X-ray photoabsorption by the interclump wind medium, the number of clump structures in the flow, and the radial distribution of the clumps. Using the two beta law prescription, the theoretical emission measure and temperature distribution throughout the wind can be derived. The emission measure tends to be power law, and the temperature distribution broad in terms of wind velocity. Although restricted to the case of adiabatic cooling, our models highlight the influence of bow shock effects for hot plasma temperature and emission measure distributions in stellar winds and their impact on X-ray line profile shapes. Previous models have focused on geometrical considerations of the clumps and their distribution in the wind. Our results represent the first time that the temperature distribution of wind clump structures are explicitly and self-consistently accounted in modeling X-ray line profile shapes for massive stars.

Line ratios in "fir" triplets of helium-like ions have proven to be a powerful diagnostic of conditions in X-ray emitting gas surrounding massive stars. Recent observations indicate that these ratios can be variable with time. The possible causes of variation in line ratios are limited: changes in the radiation field or changes in density, and changes in mass-loss or geometry. In this paper, we investigate the ability of changes in the radiation field to induce variability in the ratio R=f/i. To isolate the radiative effect, we use a heuristic model of temperature and radius changes in variable stars in the B and O range with low-density, steady-state winds. We model the changes in emissivity of X-ray emitting gas close to the star due to differences in level-pumping from available UV photons at the location of the gas. We find that under these conditions, variability in R is dominated by the stellar temperature. Although the relative amplitude of variability is roughly comparable for most lines at most temperatures, detectable variations are limited to a few lines for each spectral type. We predict that variable values in R due to stellar variability must follow predictable trends found in our simulations. Our model uses radial pulsations as a mode of stellar variability that maximizes the amplitude of variation in R. This model is robust enough to show which ions will provide the best opportunity for observing variability in the f/i ratio at different stellar temperatures, and the correlation of that variability with other observable parameters. In real systems, the effects would be more complex than in our model, with differences in phase and suppressed amplitude in the presence of non-radial pulsations. This suggests that changes in R across many lines concurrently are not likely to be produced by a variable radiation field.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and complex magnetic field. Recently, Petit et al. 2011 presented the discovery of the first two tau Sco analogues -- HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS and Narval spectropolarimetric observations were obtained by the Magnetism in Massive Stars CFHT and TBL Large Programs, in order to characterize the stellar and magnetic properties of these stars. A magnetic field of similar surface strength was found on both stars, reinforcing the connection between the presence of a magnetic field and wind peculiarities. We present additional phase-resolved observations secured by the MiMeS collaboration for HD 66665 in order to measure its magnetic geometry, and correlate that geometry with diagnostics of mass-loss.

We present a comprehensive study of X-ray emission and wind properties of massive magnetic early B-type stars. Dedicated XMM-Newton observations were obtained for three stars xi1 CMa, V2052 Oph, and zeta Cas. We report the first detection of X-ray emission from V2052 Oph and zeta Cas. The observations show that the X-ray spectra of our program stars are quite soft. We compile the complete sample of early B-type stars with detected magnetic fields to date and existing X-ray measurements, in order to study whether the X-ray emission can be used as a general proxy for stellar magnetism. We find that hard and strong X-ray emission does not necessarily correlate with the presence of a magnetic field. We analyze the UV spectra of five non-supergiant B stars with magnetic fields by means of non-LTE iron-blanketed model atmospheres. The latter are calculated with the PoWR code, which treats the photosphere as well as the the wind, and also accounts for X-rays. Our models accurately fit the stellar photospheric spectra in the optical and the UV. The parameters of X-ray emission, temperature and flux are included in the model in accordance with observations. We confirm the earlier findings that the filling factors of X-ray emitting material are very high. Our analysis reveals that the magnetic early type B stars studied here have weak winds. The mass-loss rates are significantly lower than predicted by hydrodynamically consistent models. We find that, although the X-rays strongly affect the ionization structure of the wind, this effect is not sufficient in reducing the total radiative acceleration. When the X-rays are accounted for at the intensity and temperatures observed, there is still sufficient radiative acceleration to drive stronger mass-loss than we empirically infer from the UV spectral lines. (abridged)

Aims: This paper extends previous studies in modeling time varying linear polarization due to axisymmetric magnetic fields in rotating stars. We use the Hanle effect to predict variations in net line polarization, and use geometric arguments to generalize these results to linear polarization due to other mechanisms. Methods: Building on the work of Lopez Ariste et al., we use simple analytic models of rotating stars that are symmetric except for an axisymmetric magnetic field to predict the polarization lightcurve due to the Hanle effect. We highlight the effects for the variable line polarization as a function of viewing inclination and field axis obliquity. Finally, we use geometric arguments to generalize our results to linear polarization from the weak transverse Zeeman effect. Results: We derive analytic expressions to demonstrate that the variable polarization lightcurve for an oblique magnetic rotator is symmetric. This holds for any axisymmetric field distribution and arbitrary viewing inclination to the rotation axis. Conclusions: For the situation under consideration, the amplitude of the polarization variation is set by the Hanle effect, but the shape of the variation in polarization with phase depends largely on geometrical projection effects. Our work generalizes the applicability of results described in Lopez Ariste et al., inasmuch as the assumptions of a spherical star and an axisymmetric field are true, and provides a strategy for separating the effects of perspective from the Hanle effect itself for interpreting polarimetric lightcurves.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and highly complex magnetic field. This paper presents the discovery of the first two tau Sco analogues - HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS spectropolarimetric observations were secured by the Magnetism in Massive Stars CFHT Large Program, in order to characterize the stellar and magnetic properties of these stars. CMFGEN modelling of optical ESPaDOnS spectra and archived IUE UV spectra showed that these stars have stellar parameters similar to those of tau Sco. A magnetic field of similar surface strength is found on both stars, reinforcing the connection between the presence of a magnetic field and wind peculiarities. However, additional phase-resolved observations will be required in order to assess the potential complexity of the magnetic fields, and verify if the wind anomalies are linked to this property.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its high X-ray activity and remarkable wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star that still defies interpretation in terms of a physically coherent model. Recently, two early B-type stars were discovered as tau Sco analogues, identified by the striking similarity of their UV spectra to that of tau Sco, which was - until now - unique among OB stars. We present the recent detection of their magnetic fields by the MiMeS collaboration, reinforcing the connection between the presence of a magnetic field and wind anomalies (Petit et al. 2010). We will also present ongoing observational efforts undertaken to establish the precise magnetic topology, in order to provide additional constrains for existing models attempting to reproduce the unique wind structure of tau Sco-like stars.

X-ray emission is ubiquitous among massive stars. In the last decade, X-ray observations revolutionized our perception of stellar winds but opened a Pandora's box of urgent problems. X-rays penetrating stellar winds suffer mainly continuum absorption, which greatly simplifies the radiative transfer treatment. The small and large scale structures in stellar winds must be accounted for to understand the X-ray emission from massive stars. The analysis of X-ray spectral lines can help to infer the parameters of wind clumping, which is prerequisite for obtaining empirically correct stellar mass-loss rates. The imprint of large scale structures, such as CIRs and equatorial disks, on the X-ray emission is predicted, and new observations are testing theoretical expectations. The X-ray emission from magnetic stars proves to be more diverse than anticipated from the direct application of the magnetically-confined wind model. Many outstanding questions about X-rays from massive stars will be answered when the models and the observations advance.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its remarkable, superionized wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star that still defies interpretation in terms of a physically coherent model. Recently, two early B-type stars were discovered as tau Sco analogues, identified by the striking similarity of their UV spectra to that of tau Sco, which was - until now - unique among OB stars. We present the recent detection of their magnetic fields by the MiMeS collaboration, reinforcing the connection between the presence of a magnetic field and a superionized wind. We will also present ongoing observational efforts undertaken to establish the precise magnetic topology, in order to provide additional constrains for existing models attempting to reproduce the unique wind structure of tau Sco-like stars.

This paper focuses on the polarized profiles of resonance scattering lines that form in magnetized disks. Optically thin lines from Keplerian planar disks are considered. Model line profiles are calculated for simple field topologies of axial fields (i.e., vertical to the disk plane) and toroidal fields (i.e., purely azimuthal). A scheme for discerning field strengths and geometries in disks is developed based on Stokes Q-U diagrams for the run of polarization across line profiles that are Doppler broadened by the disk rotation. A discussion of the Hanle effect for magnetized disks in which the magnetorotational instability (MRI) is operating is also presented. Given that the MRI has a tendency to mix the vector field orientation, it may be difficult to detect the disk fields with the longitudinal Zeeman effect, since the amplitude of the circularly polarized signal scales with the net magnetic flux in the direction of the observer. The Hanle effect does not suffer from this impediment, and so a multi-line analysis could be used to constrain field strengths in disks dominated by the MRI.

We obtained relatively high signal-to-noise X-ray spectral data of the early massive star tau Sco (B0.2V) with the Suzaku XIS instrument. This source displays several unusual features that motivated our study: (a) redshifted absorption in UV P Cygni lines to approximately +250 km/s suggestive of infalling gas, (b) unusually hard X-ray emission requiring hot plasma at temperatures in excess of 10 MK whereas most massive stars show relatively soft X-rays at a few MK, and (c) a complex photospheric magnetic field of open and closed field lines. In an attempt to understand the hard component better, X-ray data were obtained at six roughly equally spaced phases within the same epoch of tau Sco's 41 day rotation period. The XIS instrument has three operable detectors: XIS1 is back illuminated with sensitivity down to 0.2 keV; XIS0 and XIS2 are front illuminated with sensivitity only down to 0.4 keV and have overall less effective area than XIS1. The XIS0 and XIS3 detectors show relatively little variability. In contrast, there is a 4sigma detection of a 4% drop in the count rate of the XIS1 detector at one rotational phase. In addition, all three detectors show a 3% increase in count rate at the same phase. The most optimistic prediction of X-ray variability allows for a 40% change in the count rate, particularly near phases where we have pointings. Observed modulations in the X-ray light curve on the rotation cycle is an order of magnitude smaller than this, which places new stringent constraints on future modeling of this interesting magnetic massive star.

We obtained X-ray XMM-Newton observations of the open cluster Berkeley 87 and the massive star-forming region (SFR) ON 2. In addition, archival infrared Spitzer Space Telescope observations were used. It is likely that the SFR ON 2 and Berkeley 87 are at the same distance, 1.23 kpc, and hence are associated. The XMM-Newton observations detected X-rays from massive stars in Berkeley 87 as well as diffuse emission from the SFR ON 2. The two patches of diffuse X-ray emission are encompassed in the shell-like H II region GAL 75.84+0.40 in the northern part of ON 2 and in the ON 2S region in the southern part of ON 2. The diffuse emission from GAL 75.84+0.40 suffers an absorption column equivalent to A_V approx. 28 mag. Its spectrum can be fitted either with a thermal plasma model at T < 30 MK or by an absorbed power-law model with gamma; approx. -2.6. The X-ray luminosity of GAL 75.84+0.40 is L_X approx. 6 10^31 erg/s. The diffuse emission from ON 2S is adjacent to the ultra-compact H II (UCHII) region Cygnus 2N, but does not coincide with it or with any other known UCHII region. It has a luminosity of L_X approx. 4 10^31 erg/s. The spectrum can be fitted with an absorbed power-law model with gamma; approx.-1.4. We adopt the view of Turner and Forbes (1982) that the SFR ON 2 is physically associated with the massive star cluster Berkeley 87 hosting the WO type star WR 142. We discuss different explanations for the apparently diffuse X-ray emission in these SFRs. These include synchrotron radiation, invoked by the co-existence of strongly shocked stellar winds and turbulent magnetic fields in the star-forming complex, cluster wind emission, or an unresolved population of discrete sources.

A number of substellar companions to evolved cool stars have now been reported. Cool giants are distinct from their progenitor Main Sequence (MS) low-mass stars in a number of ways. First, the mass loss rates of cool giant stars are orders of magnitude greater than for the late-type MS stars. Second, on the cool side of the Linsky-Haisch "dividing line", K and M giant stars are not X-ray sources, although they do show evidence for chromospheres. As a result, cool star winds are largely neutral for those spectral types, suggesting that planetary or brown dwarf magnetospheres will not be effective in standing off the stellar wind. In this case one expects the formation of a bow shock morphology at the companion, deep inside its magnetosphere. We explore radio emissions from substellar companions to giant stars with ionised winds or neutral winds.

In dense hot star winds, the infrared and radio continua are dominated by free-free opacity and recombination emission line spectra. In the case of a spherically symmetric outflow that is isothermal and expanding at constant radial speed, the radiative transfer for the continuum emission from a dense wind is analytic. Even the emission profile shape for a recombination line can be derived. Key to these derivations is that the opacity scales with only the square of the density. These results are well-known. Here an extension of the derivation is developed that also allows for line blends and the inclusion of an additional power-law dependence beyond just the density dependence. The additional power-law is promoted as a representation of a radius dependent clumping factor. It is shown that differences in the line widths and equivalent widths of the emission lines depend on the steepness of the clumping power-law. Assuming relative level populations in LTE in the upper levels of He {\sc ii}, an illustrative application of the model to {\em Spitzer}/IRS spectral data of the carbon-rich star WR 90 is given.

The polarimetric and photometric variability of Wolf-Rayet (WR) stars as caused by clumps in the winds, is revisited. In the model which is improved from Li et al. 2000, the radial expansion of the thickness is accounted for, but we retain the dependence on the beta velocity law, stellar occultation effects. We again search for parameters that can yield results consistent with observations in regards to the mean polarization, the ratio of polarimetric to photometric variability, and the volume filling factor. Clump generation and spatial distribution are randomized by the Monte Carlo method so as to produce clumps which are, in the mean, distributed uniformly in space and have time intervals with a Gaussian distribution. The generated clumps move radially outward with a velocity law determined by a beta index, and the angular size of the clumps is assumed to keep fixed. By fitting the observational results and the volume filling factor, the clump velocity law index beta and clump ejection rate are inferred, and are found to be well constrained. In addition, the subpeak features on broad emission lines seem to support the clump ejection rate. Meanwhile, the fraction of the total mass loss rate that is contained in the clumps is obtained by fitting the observed polarization. We conclude that this picture for the clump properties produces a valuable diagnostic of WR wind structure.

This paper presents calculations for forbidden emission line profile shapes arising from colliding wind binaries. The main application is for systems involving a Wolf-Rayet (WR) star and an OB star companion. The WR wind is assumed to dominate the forbidden line emission. The colliding wind interaction is treated as an archimedean spiral with an inner boundary. Under the assumptions of the model, the major findings are as follows. (a) The redistribution of the WR wind as a result of the wind collision is not flux conservative but typically produces an excess of line emission; however, this excess is modest at around the 10% level. (b) Deviations from a flat-top profile shape for a spherical wind are greatest for viewing inclinations that are more nearly face-on to the orbital plane. At intermediate viewing inclinations, profiles display only mild deviations from a flat-top shape. (c) The profile shape can be used to constrain the colliding wind bow shock opening angle. (d) Structure in the line profile tends to be suppressed in binaries of shorter periods. (e) Obtaining data for multiple forbidden lines is important since different lines probe different characteristic radial scales. Our models are discussed in relation to ISO data for WR 147 and gamma Vel (WR11). The lines for WR 147 are probably not accurate enough to draw firm conclusions. For gamma Vel, individual line morphologies are broadly reproducible but not simultaneously so for the claimed wind and orbital parameters. Overall, the effort demonstrates how lines that are sensitive to the large-scale wind can help to deduce binary system properties and provide new tests of numerical simulations.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of stars and stellar evolution. Community-based White Paper to Astro2010 in response to the call for such papers.

We report the discovery of weak yet hard X-ray emission from the Wolf-Rayet (WR) star WR142 with the XMM-Newton X-ray telescope. Being of spectral subtype WO2, WR142 is a massive star in a very advanced evolutionary stage, short before its explosion as a supernova or gamma-ray burst. This is the first detection of X-ray emission from a WO-type star. We rule out any serendipitous X-ray sources within approx 1" of WR142. WR142 has an X-ray luminosity of L_X=7\times10^{30} erg/s, which constitutes only $\lsim 10^{-8}$ of its bolometric luminosity. The hard X-ray spectrum suggests a plasma temperature of about 100 MK. Commonly, X-ray emission from stellar winds is attributed to embedded shocks due to the intrinsic instability of the radiation driving. From qualitative considerations we conclude that this mechanism cannot account for the hardness of the observed radiation. There are no hints for a binary companion. Therefore the only remaining, albeit speculative explanation must refer to magnetic activity. Possibly related, WR142 seems to rotate extremely fast, as indicated by the unusually round profiles of its optical emission lines. Our detection implies that the wind of WR142 must be relatively transparent to X-rays, which can be due to strong wind ionization, wind clumping, or non-spherical geometry from rapid rotation.

The presence of magnetic fields in O-type stars has been suspected for a long time. The discovery of such fields would explain a wide range of well documented enigmatic phenomena in massive stars, in particular cyclical wind variability, Halpha emission variations, chemical peculiarity, narrow X-ray emission lines and non-thermal radio/X-ray emission. Here we present the results of our studies of magnetic fields in O-type stars, carried out over the last years.

We consider the polarization arising from scattering in an envelope illuminated by a central anisotropic source. Spherical harmonics are used to describe both the light source anisotropy and the envelope density distribution functions of the scattering particles. This framework demonstrates how the net resultant polarization arises from a superposition of three basic "shape" functions: the distribution of source illumination, the distribution of envelope scatterers, and the phase function for dipole scattering. Specific expressions for the Stokes parameters and scattered flux are derived for the case of an ellipsoidal light source inside an ellipsoidal envelope, with principal axes that are generally not aligned. Two illustrative examples are considered: (a) axisymmetric mass loss from a rapidly rotating star, such as may apply to some Luminous Blue Variables, and (b) a Roche-lobe filling star in a binary system with a circumstellar envelope. As a general conclusion, the combination of source anisotropy with distorted scattering envelopes leads to more complex polarimetric behavior such that the source characteristics should be carefully considered when interpreting polarimetric data.

We report on low resolution (R~3000) spectropolarimetry of the A0 supergiant star HD 92207. This star is well-known for significant spectral variability. The source was observed on seven different nights spanning approximately 3 months in time. With a rotation period of approximately 1 year, our data covers approximately a quarter of the star's rotational phase. Variability in the continuum polarization level is observed over this period of time. The polarization across the Halpha line on any given night is typically different from the degree and position angle of the polarization in the continuum. Interestingly, Hbeta is not in emission and does not show polarimetric variability. We associate the changes at Halpha as arising in the wind, which is in accord with the observed changes in the profile shape and equivalent width of Halpha along with the polarimetric variability. For the continuum polarization, we explore a spiral shaped wind density enhancement in the equatorial plane of the star, in keeping with the suggestion of Kaufer etal (1997). Variable polarization signatures across Halpha are too complex to be explained by this simple model and will require a more intensive polarimetric follow-up study to interpret properly.

We present preliminary calculations of electron scattering polarizations from models of structured cool star envelopes. We note that net polarizations from unresolved sources can result from non-spherical scattering envelopes and/or anisotropic illumination from a photosphere that has brightness variations. The resultant polarizations are quite small (hundreths of a percent); however, Rayleigh scattering from molecular opacity and/or dust scattering from the more extended envelope under similar considerations may produce higher polarizations.

The presence of magnetic fields in O-type stars has been suspected for a long time. The discovery of such fields would explain a wide range of well documented enigmatic phenomena in massive stars, in particular cyclical wind variability, H alpha emission variations, chemical peculiarity, narrow X-ray emission lines and non-thermal radio/X-ray emission. To investigate the incidence of magnetic fields in O stars, we acquired 38 new spectropolarimetric observations with FORS1 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Kueyen telescope of the VLT. Spectropolarimetric observations have been obtained at different phases for a sample of 13 O stars. 10 stars were observed in the spectral range 348-589nm, HD36879 and HD148937 were observed in the spectral region 325-621nm, and HD155806 was observed in both settings. To prove the feasibility of the FORS1 spectropolarimetric mode for the measurements of magnetic fields in hot stars, we present in addition 12 FORS1 observations of the mean longitudinal magnetic field in theta^1 Ori C and compare them with measurements obtained with the MuSiCoS, ESPaDOnS and Narval spectropolarimeters. Most stars in our sample which have been observed on different nights show a change of the magnetic field polarity, but a field at a significance level of 3sigma has been detected only in four stars, HD36879, HD148937, HD152408, and HD164794. The largest longitudinal magnetic field, <B_z> = -276+-88G, was detected in the Of?p star HD148937. We conclude that large-scale organised magnetic fields with polar field strengths larger than 1kG are not widespread among O-type stars.

It is now well established that stellar winds of hot stars are fragmentary and that the X-ray emission from stellar winds has a strong contribution from shocks in winds. Chandra high spectral resolution observations of line profiles of O and B stars have shown numerous properties that had not been expected. Here we suggest explanations by considering the X-rays as arising from bow shocks that occur where the stellar wind impacts on spherical clumps in the winds. We use an accurate and stable numerical hydrodynamical code to obtain steady-state physical conditions for the temperature and density structure in a bow shock. We use these solutions plus analytic approximations to interpret some major X-ray features: the simple power-law distribution of the observed emission measure derived from many hot star X-ray spectra and the wide range of ionization stages that appear to be present in X-ray sources throughout the winds. Also associated with the adiabatic cooling of the gas around a clump is a significant transverse velocity for the hot plasma flow around the clumps, and this can help to understand anomalies associated with observed line widths, and the differences in widths seen in stars with high and low mass-loss rates. The differences between bow shocks and the planar shocks that are often used for hot stars are discussed. We introduce an ``on the shock'' (OTSh) approximation that is useful for interpreting the X-rays and the consequences of clumps in hot star winds and elsewhere in astronomy.

Microlensing surveys have proven to be tremendously fruitful in providing valuable data products for many fields of astrophysics, from eclipse lightcurves for substellar candidates to limb darkening in stellar atmospheres. We report on a program of modeling observables from microlensing of circumstellar envelopes, particularly those of red giant stars that are the most likely to show finite source effects. We will summarize work for how polarization light curves can be used to infer envelope properties and will describe recent modeling of the time dependent spectral energy distributions (SEDs) for microlensing of dusty winds. One of the most exciting developments is the possibility of measuring variable polarization from microlensing in a suitable source using the RINGO polarimeter at La Palma. Also quite interesting is the possibility of probing a dusty wind using IRAC data for a suitable source in the event that Spitzer has a ``warm'' cycle.

We review the recent discoveries of magnetic fields in different types of massive stars and briefly discuss strategies for spectropolarimetric observations to be carried out in the future.

We report on observations of the eclipsing and interacting binary beta Lyrae from the Suzaku X-ray telescope. This system involves an early B star embedded in an optically and geometrically thick disk that is siphoning atmospheric gases from a less massive late B II companion. Motivated by an unpublished X-ray spectrum from the Einstein X-ray telescope suggesting unusually hard emission, we obtained time with Suzaku for pointings at three different phases within a single orbit. From the XIS detectors, the softer X-ray emission appears typical of an early-type star. What is surprising is the remarkably unchanging character of this emission, both in luminosity and in spectral shape, despite the highly asymmetric geometry of the system. We see no eclipse effect below 10 keV. The constancy of the soft emission is plausibly related to the wind of the embedded B star and Thomson scattering of X-rays in the system, although it might be due to extended shock structures arising near the accretion disk as a result of the unusually high mass-transfer rate. There is some evidence from the PIN instrument for hard emission in the 10-60 keV range. Follow-up observations with the RXTE satellite will confirm this preliminary detection.

Given that dynamically significant magnetic fields in at least some massive stars have now been measured, our contribution addresses the question, to what extent can fields be directly detected in circumstellar gas? The question speaks directly to the very interesting topic of line-driving physics coupled with magnetized plasmas, and how this coupling produces structure in the wind flow. We focus our attention on weak-field diagnostics. These come in two main types: the Hanle effect, which pertains to coherence effects for linear polarization from line scattering, and the weak longitudinal Zeeman effect, which pertains to circular polarization in lines.

We report on neon abundances derived from {\it Spitzer} high resolution spectral data of eight Wolf-Rayet (WR) stars using the forbidden line of [\ion{Ne}{3}] 15.56 microns. Our targets include four WN stars of subtypes 4--7, and four WC stars of subtypes 4--7. We derive ion fraction abundances $\gamma$ of Ne^{2+} for the winds of each star. The ion fraction abundance is a product of the ionization fraction $Q_{\rm i}$ in stage i and the abundance by number ${\cal A}_E$ of element E relative to all nuclei. Values generally consistent with solar are obtained for the WN stars, and values in excess of solar are obtained for the WC stars.

A number of strong infrared forbidden lines have been observed in several evolved Wolf-Rayet star winds, and these are important for deriving metal abundances and testing stellar evolution models. In addition, because these optically thin lines form at large radius in the wind, their resolved profiles carry an imprint of the asymptotic structure of the wind flow. This work presents model forbidden line profile shapes formed in axisymmetric winds. It is well-known that an optically thin emission line formed in a spherical wind expanding at constant velocity yields a flat-topped emission profile shape. Simulated forbidden lines are produced for a model stellar wind with an axisymmetric density distribution that treats the latitudinal ionization self-consistently and examines the influence of the ion stage on the profile shape. The resulting line profiles are symmetric about line centre. Within a given atomic species, profile shapes can vary between centrally peaked, doubly peaked, and approximately flat-topped in appearance depending on the ion stage (relative to the dominant ion) and viewing inclination. Although application to Wolf-Rayet star winds is emphasized, the concepts are also relevant to other classes of hot stars such as luminous blue variables and Be/B[e] stars.

This paper examines line profile evolution due to the linear expansion of circumstellar material obsverved during a microlensing event. This work extends our previous papers on emission line profile evolution from radial and azimuthal flow during point mass lens events and fold caustic crossings. Both ''flavours'' of microlensing were shown to provide effective diagnostics of bulk motion in circumstellar envelopes. In this work a different genre of flow is studied, namely linear homologous expansion, for both point mass lenses and fold caustic crossings. Linear expansion is of particular relevance to the effects of microlensing on supernovae at cosmological distances. We derive line profiles and equivalent widths for the illustrative cases of pure resonance and pure recombination lines, modelled under the Sobolev approximation. The efficacy of microlensing as a diagnostic probe of the stellar environs is demonstrated and discussed.

The high gradient of magnification across the source during both small impact parameter events and caustic crossings offers a unique opportunity for determining the surface brightness profile of the source. Furthermore, models indicate that these events can also provide an appreciable polarisation signal -- arising from differential magnification across the otherwise symmetric source. In a previous paper, polarimetric variations from point lensing of a circumstellar envelope were considered, as would be suitable for an extended envelope around a red giant. In this work we examine the polarisation in the context of caustic crossing events, the scenario which represents the most easily accessible situation for actually observing a polarisation signal in Galactic microlensing.

The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0), where l_0 is insensitive to limiting magnitude, band choice, and side of Galactic center: l_0= 17-30 degrees with a best fit value in the the 4.5 micron band of l_0=24 +/- 4 degrees. Modeling the source distribution as an exponential disk yields a radial scale length of H= 3.9 +/- 0.6 kpc. There is a pronounced north-south asymmetry in source counts for |l| < 30 degrees, with ~ 25% more stars in the north. For l=10-30 degrees, there is a strong enhancement of stars of m= 11.5-13.5 mag. A linear bar passing through the Galactic center with half-length R_bar=4.4 +/- 0.5 kpc, tilted by phi=44 +/- 10 degrees to the Sun-Galactic Center line, provides the simplest interpretation of this data. We examine the possibility that enhanced source counts at l=26-28 degrees, 31.5-34 degrees, and 306-309 degrees are related to Galactic spiral structure. Total source counts are depressed in regions where the counts of red objects (m_K-m_[8.0] >3) peak. In these areas, the counts are reduced by extinction due to molecular gas and/or high diffuse backgrounds associated with star formation.

The hot star wind momentum problem \eta = mdot*vinfty/(L/c) >> 1 is revisited, and it is shown that the conventional belief, that it can be solved by a combination of clumping of the wind and multiple scattering of photons, is not self-consistent for optically thick clumps. Clumping does reduce the mass loss rate mdot, and hence the momentum supply, required to generate a specified radio emission measure \varepsilon, while multiple scattering increases the delivery of momentum from a specified stellar luminosity L. However, in the case of thick clumps, when combined the two effects act in opposition rather than in unison since clumping reduces multiple scattering. Thus the ratio of momentum deliverable to momentum required is maximal for a smooth wind and the momentum problem remains for the thick clump case. In the case of thin clumps, all of the benefit of clumping in reducing \eta lies in reducing \mdot for a given \varepsilon so that extremely small filling factors f = 10e-4 are needed. It is also shown that clumping affects the inference of \mdot from radio \varepsilon not only by changing the emission measure per unit mass but also by changing the radio optical depth unity radius R_{rad}, and hence the observed wind volume, at radio wavelengths.

The Hanle effect has been proposed as a new diagnostic of circumstellar magnetic fields for early-type stars, for which it is sensitive to field strengths in the 1-300 G range. In this paper we compute the polarized P-Cygni line profiles that result from the Hanle effect. For modeling the polarization, we employ a variant of the ``last scattering approximation''. For cases in which the Sobolev optical depths are greater than unity, the emergent line intensity is assumed to be unpolarized; while for smaller optical depths, the Stokes source functions for the Hanle effect with optically thin line scattering are used. For a typical P Cygni line, the polarized emission forms in the outer wind, because the Sobolev optical depth is large at the inner wind. For low surface field strengths, weak P Cygni lines are needed to measure the circumstellar field. For high values of the surface fields, both the Zeeman and Hanle diagnostics can be used, with the Zeeman effect probing the photospheric magnetic fields, and the Hanle effect measuring the magnetic field in the wind flow. Polarized line profiles are calculated for a self-consistent structure of the flow and the magnetic geometry based on the WCFields model, which is applicable to slowly rotating stellar winds with magnetic fields drawn out by the gas flow. For surface fields of a few hundred Gauss, we find that the Hanle effect can produce line polarizations in the range of a few tenths of a percent up to about 2 percent.

Observations with the Infrared Space Observatory (ISO) SWS spectrometer are used to constrain the velocity law and wind clumping of the well-studied Wolf-Rayet (WR) star WR 136 (HD 192163) (WN6). Because the free-free continuum opacity in WR winds increases steadily with wavelength in the IR, each point in the continuous spectrum may be regarded as forming in a pseudo-photosphere of larger radius for longer wavelength. Using this idea in combination with an analysis of the Doppler-broadened widths of several \ion{He}{2} recombination lines, we can derive information about the velocity law and clumpiness of the stellar wind of WR 136. Adopting a $\beta$-law distribution for the outflow velocity law, we compute the continuous energy distribution and line profiles. Interestingly,the continuum data are found to constrain the wind to have fairly low clumping values for the IR-forming region.

The physics of the Hanle effect is briefly reviewed, and its application as a diagnostic of hot star magnetic fields is described. Emphasis is given to the practicalities of using spectropolarimetry of resolved wind emission lines to infer information about the circumstellar magnetic field strength and its geometry. A model for a weakly magnetized wind from "WCFields" theory is used as the backdrop for investigating polarized line profile effects for P Cygni resonance lines using a kind of "last scattering approximation". Model results are presented for a typical P Cygni line that forms in a spherical wind. Significant line polarizations of a few tenths of a percent can result for circumstellar fields of order 100G. Information about the field topology and surface field strength is gleaned from the Stokes-Q and U-profiles. Simplistically, the Q-profile polarization is governed by the field strength, and the U-profile symmetry (whether symmetric or anti-symmetric or even zero) is governed by the field geometry.

The carbon-rich WC5 star WR 114 was not detected during a 15.9 ksec XMM-Newton observation, implying an upper limit to the X-ray luminosity of Lx < 2.5*10^{30} erg/s and to the X-ray to bolometric luminosity ratio of Lx/Lbol < 4*10^{-9}. This confirms indications from earlier less sensitive measurements that there has been no convincing X-ray detection of any single WC star. This lack of detections is reinforced by XMM-Newton and Chandra observations of WC stars. Thus the conclusion has to be drawn that the stars with radiatively-driven stellar winds of this particular class are insignificant X-ray sources. We attribute this to photoelectronic absorption by the stellar wind. The high opacity of the metal-rich and dense winds from WC stars puts the radius of optical depth unity at hundreds or thousands of stellar radii for much of the X-ray band. We believe that the essential absence of hot plasma so far out in the wind exacerbated by the large distances and correspondingly high ISM column densities makes the WC stars too faint to be detectable with current technology. The result also applies to many WC stars in binary systems, of which only about 20% are identified X-ray sources, presumably due to colliding winds.

We consider the consequences of appreciable line optical depth for the profile shape of X-ray emission lines formed in stellar winds. The hot gas is thought to arise in distributed wind shocks, and the line formation is predominantly via collisional excitation followed by radiative decay. Such lines are often modelled as optically thin, but the theory has difficulty matching resolved X-ray line profiles. We suggest that for strong lines of abundant metals, newly created photons may undergo resonance scattering, modifying the emergent profile. Using Sobolev theory in a spherically symmetric wind, we show that thick-line resonance scattering leads to emission profiles that still have blueshifted centroids like the thin lines, but which are considerably less asymmetric in appearance. We focus on winds in the constant-expansion domain, and derive an analytic form for the profile shape in the limit of large line and photoabsorptive optical depths. Our theory is applied to published {\it Chandra} observations of the O star $\zeta$ Pup.

The problem of detecting Jovian-sized planets orbiting White Dwarf stars is considered. Significant IR excesses result from warm Jupiters orbiting a White Dwarf of $T_{\rm eff}=10000$ K at a distance of $\sim 10^3$ White Dwarf radii (corresponding to $\sim 10^2$ Jupiter radii or a few tenths of an AU) with an orbital period of $\sim 100$ days. Such a planet will have a 10 micron flux density at its Wien peak that is comparable to the emission of the White Dwarf at that wavelength. Although the White Dwarf is much hotter than the planet, the planet will have peak brightness at the IR, well into the Rayleigh-Jeans tail of the White Dwarf, plus Jovians are about 10 times larger than White Dwarfs, so there is a substantial gain in the planet to star brightness contrast as compared to planets around Main Sequence stars. In the solar neighborhood, there are 51 White Dwarf stars within 13 pc of the Sun. At 10 pc, the IR flux density of ``warm'' Jupiters (a few hundred Kelvin) will fall in the range 10--100 micro-Jansky which should be observable with {\it SIRTF}.

Resonance lines that are sensitive to the Hanle effect are prominent in the UV spectra of early-type stars. To understand the differences from the solar application of the Hanle effect, we focus on the formation of P-Cygni lines both as a scattering process, and as one that allows a spectral isolation of sectors in the wind. Some complications occuring in the solar case are found to be absent for the Hanle effect for hot stars. Rocket observations from the Far Ultraviolet SpectroPolarimeter (FUSP) experiment should allow for a determination of fields in the dynamically interesting range from 1 to 300 Gauss.

We describe X-ray production in the atmospheres of hot, early-type stars in the framework of a ``stochastic shock model''. The extended envelope of a star is assumed to possess numerous X-ray emitting ``hot'' zones that are produced by shocks and embedded in the ambient ``cold'' medium in dynamical equilibrium. It is shown that the apparent lack of X-ray variability on short (order of hours) timescales do not contradict a shock model for X-ray production. The character of the X-ray variability is found to depend on the frequency with which hot zones are generated, the cool wind opacity to X-rays, and the wind flow parameters, such as mass loss rate and terminal speed.

We examine the effects of optically thick line forces on orbiting circumstellar disks, such as occur around Be stars. For radially streaming radiation, line forces are only effective if there is a strong radial velocity gradient, as occurs, for example, in a line-driven stellar wind. However, within an orbiting disk, the radial shear of the azimuthal velocity leads to strong line-of-sight velocity velocity gradients along nonradial directions. As such, in the proximity of a stellar surface extending over a substantial cone angle, the nonradial stellar radiation can impart a significant line force, even in the case of purely circular orbits. Given the highly supersonic nature of orbital velocity variations, we use the Sobolev approximation, thereby extending to the disk case the standard CAK formalism developed for line-driven winds. We delineate the parameter regimes for which radiative forces might alter disk properties; but even when radiative forces are small, we analytically quantify higher-order effects in the linear limit, including the precession of weakly elliptical orbits. We find that optically thick line forces can have observable implications for the dynamics of disks around Be stars, including the generation of either prograde or retrograde precession in slightly eccentric orbits. However, our analysis suggests a net retrograde effect, in apparent contradiction with observed long-term variations of violet/red line profile asymmetries from Be stars, which are generally thought to result from prograde propagation of a so-called ``one arm mode''. We also conclude that radiative forces may alter the dynamical properties at the surface of the disk where disk winds originate, and may even make low-density disks vulnerable to being blown away.

We derive a method for inverting emission line profiles formed in supernova ejecta. The derivation assumes spherical symmetry and homologous expansion (i.e., $v(r) \propto r$), is analytic, and even takes account of occultation by a pseudo-photosphere. Previous inversion methods have been developed which are restricted to optically thin lines, but the particular case of homologous expansion permits an analytic result for lines of {\it arbitrary} optical depth. In fact, we show that the quantity that is generically retrieved is the run of line intensity $I_\lambda$ with radius in the ejecta. This result is quite general, and so could be applied to resonance lines, recombination lines, etc. As a specific example, we show how to derive the run of (Sobolev) optical depth $\tau_\lambda$ with radius in the case of a pure resonance scattering emission line.

This paper presents first results on the line profile shapes from a circumstellar envelope in bulk motion as modified by a microlensing event. Only geometrically and optically thin spherical shells in uniform expansion or rotation are considered here so as to emphasise the information content available in the profile shapes. In particular it is demonstrated that for the case of expansion, the line emission can increase by significant factors and the time variation of the profile shape is symmetric about line centre. For uniform rotation the line emission also increases significantly, but the time evolution of the profile shape is distinctly asymmetric. Thus, microlensing is seen to yield information about the velocity field in the extended envelope. We elaborate on (a) the observational advantages of tailoring microlensing programs toward detecting extended circumstellar envelopes, (b) the use of multiline observations to infer other properties of the envelopes, such as the ionization stratification, (c) the use of the continuum excess emission at infrared wavelengths as a means of probing the envelope structure, and (d) the use of polarisation for constraining the properties of ``clumpy'' winds.