We present late-time Hubble Space Telescope (HST) images of the site of supernova (SN) 2009ip taken almost 3 yr after its bright 2012 luminosity peak. SN 2009ip is now slightly fainter in broad filters than the progenitor candidate detected by HST in 1999. The current source continues to be dominated by ongoing late-time circumstellar material interaction that produces strong Ha emission and a weak pseudo-continuum, as found previously for 1-2 yr after explosion. The intent of these observations was to search for evidence of recent star formation in the local (similar to 1 kpc; 10 arcsec) environment around SN 2009ip, in the remote outskirts of its host spiral galaxy NGC 7259. We can rule out the presence of any massive star-forming complexes like 30 Dor or the Carina nebula at the SN site or within a few kpc. If the progenitor of SN 2009ip was really a 50-80 M-circle dot star as archival HST images suggested, then it is strange that there is no sign of this type of massive star formation anywhere in the vicinity. A possible explanation is that the progenitor was the product of a merger or binary mass transfer, rejuvenated after a lifetime that was much longer than 4-5 Myr, allowing its natal H II region to have faded. A smaller region like the Orion nebula would be an unresolved but easily detected point source. This is ruled out within similar to 1.5 kpc around SN 2009ip, but a small H II region could be hiding in the glare of SN 2009ip itself. Later images after a few more years have passed are needed to confirm that the progenitor candidate is truly gone and to test for the possibility of a small H II region or cluster at the SN position.

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

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 unequal source functions to reveal that line shapes can be asymmetric about line center. A parameter study is presented and applications discussed.

We present new, contemporaneous Hubble Space Telescope STIS and XMM-Newton observations of the O7 III(n)((f)) star ξ Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of ξ Per. The variable wind of this star is known to have a 2.086-day periodicity. We use a simple, heuristic spot model that 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 yr earlier. We also show that the time variability of the new XMM-Newton fluxes is also consistent with the 2.086-day period. Thus, our new, multiwavelength coordinated observations demonstrate that the mechanism that 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 multiwavelength 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. Furthermore, the high-resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blueshifted. 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°, on the surface of the star. We note that the presence and persistence of two spots separated by 180° suggest that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in ξ Per.

The stellar winds of hot stars have an important impact on both stellar and galactic evolution, yet their structure and internal processes are not fully understood in detail. One of the best nearby laboratories for studying such massive stellar winds is the O4I(n)fp star ζ Pup. After briefly discussing existing X-ray observations from Chandra and XMM, we present a simulation of X-ray emission line profile measurements for the upcoming 840 kilosecond Chandra HETGS observation. This simulation indicates that the increased S/N of this new observation will allow several major steps forward in the understanding of massive stellar winds. By measuring X-ray emission line strengths and profiles, we should be able to differentiate between various stellar wind models and map the entire wind structure in temperature and density. This legacy X-ray spectrum of ζ Pup will be a useful benchmark for future X-ray missions.

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The similar to 3.5mag brightening in the V band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1mag for 6. days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from similar to 10,000K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km s(-1). The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H-2 is detected in emission from vibrational levels as high as v = 11, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, though the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

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 ζ 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/Solar Mass Ejection Imager, 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 multisite spectroscopic monitoring of the star unveiled cyclical modulation of its He IIλ4686 wind emission line with the 1.78-d rotation period, showing signatures of corotating interaction regions 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 λ4686 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 ζ 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.

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

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

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