Search

We fit the high-resolution Chandra X-ray spectra of the O supergiant ζ Puppis using the variable boundary condition (VBC) line model to test the stability of its mass-loss rate between two epochs of observation: 2000 March and 2018 July – 2019 August. At issue is whether the observed variations are induced by global changes in the cool (unshocked) wind itself or are isolated to the local pockets of hot gas (i.e. changes in the frequency and location of the shocks). Evidence in the literature favoured the possibility of a 40 per cent increase in the mass flux of the entire stellar wind, based on X-ray reabsorption from a line-deshadowing-instability-inspired parametrization, whereas our fit parameters are consistent with a constant mass flux with a change in the velocity variations that determine the locations where shocks form. Our results suggest the shocks in the more recent data are formed at somewhat larger radii, mimicking the enhanced blueshifts and increased line fluxes interpreted in the previous analysis as being due to increases in both the X-ray generation and reabsorption from an overall stronger wind. [ABSTRACT FROM AUTHOR]

We use the line ratios of resolved X-ray resonance and satellite line pairs from Si  xiii and Mg  xi ions to test the widely held assumption that the shocked plasma in the wind of hot stars is in collisional ionization equilibrium. We specifically apply this test to the prototypical O supergiant ζ Puppis, as its line ratios are well resolved in its total 813 ks Chandra observation. This satellite-to-resonance line ratio is known to be a reliable diagnostic of plasma temperature but has yet to be applied to hot stars due to needing a long exposure time to see both lines above the continuum. In testing the assumed plasma state, we also provide evidence against significant ionization and excitation from proton collisions. Since dielectronic recombination (DR) allows somewhat lower energy free electrons to excite resonance lines than by the usual collisional pathways, we also use the DR flux ratios to extend diagnostics of the ionization fraction of the parent ion to lower temperatures than any other method can achieve. Also, the DR ratio supports a concept of an average free electron temperature, for each target ion used, which can constrain heating models. [ABSTRACT FROM AUTHOR]

We fit a new line shape model to Chandra X-ray spectra of the O supergiant ζ Puppis to test the robustness of mass-loss rates derived from X-ray wind line profiles against different assumed heating models. Our goal is to track the hot gas by replacing the common assumption that it is proportional to the cool gas emission measure. Instead of assuming a turn-on radius for the hot gas (as appropriate for the line-deshadowing instability internal to the wind), we parametrize the hot gas in terms of a mean-free path for accelerated low-density gas to encounter slower high-density material. This alternative model is equally successful as previous approaches at fitting X-ray spectral lines in the 5–17 Å wavelength range. We find that the characteristic radius where the hottest gas appears is inversely proportional to line-formation temperature, suggesting that stronger shocks appear generally closer to the surface. This picture is more consistent with pockets of low density, rapid acceleration at the lower boundary than with an internally generated wind instability. We also infer an overall wind mass-loss rate from the profile shapes with a technique used previously in the literature. In doing so, we find evidence that the mass-loss rate derived from X-ray wind line profiles is not robust with respect to changes in the specific heating picture used. [ABSTRACT FROM AUTHOR]

We attempt to determine the driver for clumping in hot-star winds by extending the measure of the spectral variability level of Galactic Wolf–Rayet stars to by far the hottest known among them, the WN2 star WR 2 and the WO2 stars WR 102 and WR 142. These three stars have T⋆ = 140 and 200 kK, the last two being well above the bulk of WR stars with T⋆ ∼ 40–120 kK. This full temperature range for WR stars is much broader than that of their O-star progenitors (∼30–50 kK), so is better suited to look for any temperature dependence of wind clumping. We have obtained multiple observations with high signal-to-noise, moderate-resolution spectroscopy in search of small-scale variability in the strong emission lines from the dense winds of these three extreme stars, and find a very low-level of variability in both stars. Temperature and terminal velocity are correlated, so faster winds show a lower variability, though this trend goes against any predictions made involving Line Deshadowing Instability (LDI) only, implying that instabilities intrinsic to LDI are not the main source of wind clumping. Instead, it could be taken as support for the suggestion that clumps are caused by a subsurface convection zone (SSCZ) at T ∼ 170 kK, since such an SSCZ would have little opportunity to operate under the hydrostatic surface of these hottest WR stars. It is still possible, however, that an SSCZ-related driver could interact with nonlinear line instability effects to enhance or possibly even produce clumps. [ABSTRACT FROM AUTHOR]

We present a new method for using measured X-ray emission line fluxes from O stars to determine the shock-heating rate due to instabilities in their radiation-driven winds. The high densities of these winds means that their embedded shocks quickly cool by local radiative emission, while cooling by expansion should be negligible. Ignoring for simplicity any non-radiative mixing or conductive cooling, the method presented here exploits the idea that the cooling post-shock plasma systematically passes through the temperature characteristic of distinct emission lines in the X-ray spectrum. In this way, the observed flux distribution among these X-ray lines can be used to construct the cumulative probability distribution of shock strengths that a typical wind parcel encounters as it advects through the wind. We apply this new method to Chandra grating spectra from five O stars with X-ray emission indicative of embedded wind shocks in effectively single massive stars. The results for all the stars are quite similar: the average wind mass element passes through roughly one shock that heats it to at least 106 K as it advects through the wind, and the cumulative distribution of shock strengths is a strongly decreasing function of temperature, consistent with a negative power law of index n ≈ 3, implying a marginal distribution of shock strengths that scales as T−4, and with hints of an even steeper decline or cut-off above 107 K. [ABSTRACT FROM AUTHOR]

In collisional ionization equilibrium, the X-ray spectrum from a plasma depends on the differential emission measure (DEM), distributed over temperature. Due to the well-known ill conditioning problem, no precisely resolved DEM can be inverted directly from the spectrum, so often only a gross parameterization of the DEM is used to approximate the data, in hopes that the parameterization can provide useful model-independent constraints on the heating process. However, ill conditioning also introduces ambiguity into the various different parameterizations that could approximate the data, which may spoil the perceived advantages of model independence. Thus, this paper instead suggests a single parameterization for both the heating mechanism and the X-ray sources, based on a model of impulsive heating followed by radiative cooling. This approach is similar to a “cooling flow” approach but allows injection at multiple initial temperatures and applies even when the steady state is a distribution of different shock strengths, as for a standing shock with a range of obliquities, or for embedded stochastic shocks that are only steady in a statistical sense. This produces an alternative parameterization for X-ray spectra that is especially streamlined for higher density plasmas with efficient radiative cooling and minimal thermal conduction and mixing. The method also provides some internal consistency checks on the validity of its assumptions. A heuristic general version is then applied over a wide range of astrophysical applications to schematically explore potential alternative models for these phenomena. [ABSTRACT FROM AUTHOR]