The time evolution of angular momentum and surface rotation of massive stars is strongly influenced by fossil magnetic fields via magnetic braking. We present a new module containing a simple, comprehensive implementation of such a field at the surface of a massive star within the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. We test two limiting scenarios for magnetic braking: distributing the angular momentum loss throughout the star in the first case, and restricting the angular momentum loss to a surface reservoir in the second case. We perform a systematic investigation of the rotational evolution using a grid of OB star models with surface magnetic fields ($M_\star=5-60$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2-1.0$, $B_{\rm p} =1-20$ kG). We then employ a representative grid of B-type star models ($M_\star=5, 10, 15$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2 , 0.5, 0.8$, $B_{\rm p} = 1, 3 ,10, 30$ kG) to compare to the results of a recent self-consistent analysis of the sample of known magnetic B-type stars. We infer that magnetic massive stars arrive at the zero age main sequence with a range of rotation rates, rather than with one common value. In particular, some stars are required to have close-to-critical rotation at the ZAMS. However, magnetic braking yields surface rotation rates converging to a common low value, making it difficult to infer the initial rotation rates of evolved, slowly-rotating 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.

Magnetic confinement of stellar winds leads to the formation of magnetospheres, which can be sculpted into Centrifugal Magnetospheres (CMs) by rotational support of the corotating plasma. The conditions required for the CMs of magnetic early B-type stars to yield detectable emission in H$\alpha$ -- the principal diagnostic of these structures -- are poorly constrained. A key reason is that no detailed study of the magnetic and rotational evolution of this population has yet been performed. Using newly determined rotational periods, modern magnetic measurements, and atmospheric parameters determined via spectroscopic modelling, we have derived fundamental parameters, dipolar oblique rotator models, and magnetospheric parameters for 56 early B-type stars. Comparison to magnetic A- and O-type stars shows that the range of surface magnetic field strength is essentially constant with stellar mass, but that the unsigned surface magnetic flux increases with mass. Both the surface magnetic dipole strength and the total magnetic flux decrease with stellar age, with the rate of flux decay apparently increasing with stellar mass. We find tentative evidence that multipolar magnetic fields may decay more rapidly than dipoles. Rotational periods increase with stellar age, as expected for a magnetic braking scenario. Without exception, all stars with H$\alpha$ emission originating in a CM are 1) rapid rotators, 2) strongly magnetic, and 3) young, with the latter property consistent with the observation that magnetic fields and rotation both decrease over time.

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS.

Motivated by the eruptive mass loss inferred from Luminous Blue Variable (LBV) stars, we present 1D hydrodynamical simulations of the response from sudden energy injection into the interior of a very massive ($100 \, M_\odot$) star. For a fiducial case with total energy addition set to a factor $f=0.5$ of the net stellar binding energy, and applied within the stellar envelope, we detail the dynamical response that leads to ejection of the outermost $7.2 \, M_\odot$. We find that the ejecta's variations in time $t$ and radius $r$ for the velocity $v$, density $\rho$, and temperature $T$ are quite well fit by similarity forms in the variable $r/t \approx v$. Specifically the scaled density follows a simple exponential decline $\rho t^{3} \sim \exp (-r/v_{\rm o} t)$. This `exponential similarity' leads to analytic scaling relations for total ejecta mass $\Delta M$ and kinetic energy $\Delta K$ that agree well with the hydrodynamical simulations, with the specific-energy-averaged speed related to the exponential scale speed $v_{\rm o}$ through ${\bar v} \equiv \sqrt{2 \Delta K/\Delta M} = \sqrt{12} \, v_{\rm o}$, and a value comparable to the star's surface escape speed, $v_{\rm esc}$. Models with energy added in the core develop a surface shock breakout that propels an initial, higher-speed ejecta ($>$5000km s$^{-1}$), but the bulk of the ejected material still follows the same exponential similarity scalings with ${\bar v} \approx v_{\rm esc}$. A broader parameter study examines how the ejected mass and energy depends on the energy-addition factor $f$, for three distinct model series that locate the added energy in either the core, envelope, or near-surface. We conclude by discussing the relevance of these results for understanding LBV outbursts and other eruptive phenomena, such as failed supernovae and pulsational pair instability events.

In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state and high state are the largest observed in any magnetic O-type star, consistent with the expected properties of NGC 1624-2's magnetosphere. We also observe a redshifted absorption component in the low state, a feature not seen in most stars. We present preliminary modelling efforts based on the Analytic Dynamical Magnetosphere (ADM) formalism, demonstrating the necessity of using non-spherically symmetric models to determine wind/magnetospheric properties of magnetic O stars.

In studying the formation of massive stars, it is essential to consider the strong radiative feedback on the stars' natal environments from their high luminosities ($10^4 \sim 10^6 L_\odot$). Given that massive stars contract to main-sequence-like radii before accretion finishes, one form this feedback takes is UV line-acceleration, resulting in outflows much like those expected from main-sequence massive stars. As shown by the prior papers in this series, in addition to driving stellar winds, such line forces also ablate the surface layers off of circumstellar discs within a few stellar radii of the stellar photosphere. This removal of material from an accretion disc in turn results in a decreased accretion rate onto the forming star. Quantifying this, however, requires accounting for the continuum optical depth of the disc along the non-radial rays required for the three-dimensional line-acceleration prescription used in this paper series. We introduce the "thin disc approximation", allowing these continuum optical depths arising from an optically thick but geometrically thin disc to be dynamically treated in the context of radiation-hydrodynamics simulations. Using this approximation in full dynamical simulations, we show that such continuum optical depth effects only reduce the disc ablation by 30 percent or less relative to previous simulations that ignored continuum absorption.

We present results of simulations to assess the feasibility of modeling outflows from massive stars using the Los Alamos 3-D radiation hydrodynamics code Cassio developed for inertial confinement fusion (ICF) applications. We find that a 1-D stellar envelope simulation relaxes into hydrostatic equilibrium using computing resources that would make the simulation tractable in 2-D. We summarize next steps to include more physics fidelity and model the response to a large and abrupt energy deposition at the base of the envelope.

Spectropolarimetric surveys reveal that 8-10\% of OBA stars harbor large-scale magnetic fields, but thus far no such fields have been detected in any classical Be stars. Motivated by this, we present here MHD simulations for how a pre-existing Keplerian disc -- like that inferred to form from decretion of material from rapidly rotating Be stars -- can be disrupted by a rotation-aligned stellar dipole field. For characteristic stellar and disc parameters of a near-critically rotating B2e star, we find that a polar surface field strength of just 10 G can significantly disrupt the disc, while a field of 100 G, near the observational upper limit inferred for most Be stars, completely destroys the disc over just a few days. Our parameter study shows that the efficacy of this magnetic disruption of a disc scales with the characteristic plasma beta {(defined as the ratio between thermal and magnetic pressure)} in the disc, but is surprisingly insensitive to other variations, e.g. in stellar rotation speed, or the mass loss rate of the star's radiatively driven wind. The disc disruption seen here for even a modest field strength suggests that the presumed formation of such Be discs by decretion of material from the star would likely be strongly inhibited by such fields; this provides an attractive explanation for why no large-scale fields are detected from such Be stars.

We analyze recent 2D simulations of the nonlinear evolution of the line-deshadowing instability (LDI) in hot-star winds, to quantify how the associated highly clumped density structure can lead to a "turbulent porosity" reduction in continuum absorption and/or scattering. The basic method is to examine the statistical variations of mass column as a function of path length, and fit these to analytic forms that lead to simple statistical scalings for the associated mean extinction. A key result is that one can characterize porosity effects on continuum transport in terms of a single "turbulent porosity length", found here to scale as $H \approx (f_{\rm cl} - 1) a$, where $f_{\rm cl} \equiv \left < \rho^2 \right >/\left < \rho \right >^2$ is the clumping factor in density $\rho$, and $a$ is the density autocorrelation length. For continuum absorption or scattering in an optically thick layer, we find the associated effective reduction in opacity scales as $\sim 1/\sqrt{1+\tau_{\rm H}}$, where $\tau_{\rm H} \equiv \kappa \rho H$ is the local optical thickness of this porosity length. For these LDI simulations, the inferred porosity lengths are small, only about a couple percent of the stellar radius, $H \approx 0.02 R_\ast$. For continuum processes like bound-free absorption of X-rays that are only marginally optically thick throughout the full stellar wind, this implies $\tau_{\rm H} \ll 1$, and thus that LDI-generated porosity should have little effect on X-ray transport in such winds. The formalism developed here could however be important for understanding the porous regulation of continuum-driven, super-Eddington outflows from luminous blue variables.

In this contribution, we present BRITE observations of the early-B supergiants $\epsilon$ Ori and $\kappa$ Ori. We perform a preliminary analysis of the data acquired over the first two Orion observing runs. We evaluate whether they are compatible with co-rotating bright spots and discuss the challenges of such an approach.

In the subset of luminous, early-type stars with strong, large-scale magnetic fields and moderate to rapid rotation, material from the star's radiatively driven stellar wind outflow becomes trapped by closed magnetic loops, forming a centrifugally supported, co-rotating magnetosphere. We present here a semi-analytic analysis of how this quasi-steady accumulation of wind mass can be balanced by losses associated with a combination of an outward, centrifugally driven drift in the region beyond the Kepler co-rotation radius, and an inward/outward diffusion near this radius. We thereby derive scaling relations for the equilibrium spatial distribution of mass, and the associated emission measure for observational diagnostics like Balmer line emission. We discuss the potential application of these relations for interpreting surveys of the emission line diagnostics for OB stars with centrifugally supported magnetospheres. For a specific model of turbulent field-line-wandering rooted in surface motions associated with the iron opacity bump, we estimate values for the associated diffusion and drift coefficients.

We extend our previous study of radiative ablation of circumstellar disks by line-scattering of the star's radiation, accounting now for the effect of multiple line resonances off the stellar limb. For an analytic, three-dimensional model of the velocity structure of an equatorial Keplerian disk bounded at higher latitudes by a radially accelerating stellar wind outflow, we use root-finding methods to identify multiple resonances from a near-disk circumstellar location along starward rays both on and off the stellar core. Compared to our previous study that accounted only for the effect of on-core resonances in reducing the radiative driving through the scattering of radiation away from a near-disk circumstellar location, including off-limb resonances leads to additional radiative driving from scattering toward this location. Instead of the up-to-50% reduction in line-acceleration previously inferred from multiple resonance effects, we now find a more modest 15-20% net reduction.

Context: Clumping in the radiation-driven winds of hot, massive stars arises naturally due to the strong, intrinsic instability of line-driving (the `LDI'). But LDI wind models have so far mostly been limited to 1D, mainly because of severe computational challenges regarding calculation of the multi-dimensional radiation force. Aims: To simulate and examine the dynamics and multi-dimensional nature of wind structure resulting from the LDI. Methods: We introduce a `pseudo-planar', `box-in-a-wind' method that allows us to efficiently compute the line-force in the radial and lateral directions, and then use this approach to carry out 2D radiation-hydrodynamical simulations of the time-dependent wind. Results: Our 2D simulations show that the LDI first manifests itself by mimicking the typical shell-structure seen in 1D models, but how these shells then quickly break up into complex 2D density and velocity structures, characterized by small-scale density `clumps' embedded in larger regions of fast and rarefied gas. Key results of the simulations are that density-variations in the well-developed wind statistically are quite isotropic and that characteristic length-scales are small; a typical clump size is ~0.01R at 2R, thus resulting also in rather low typical clump-masses ~10^17 g. Overall, our results agree well with the theoretical expectation that the characteristic scale for LDI-generated wind-structure is of order the Sobolev length. We further confirm some earlier results that lateral `filling-in' of radially compressed gas leads to somewhat lower clumping factors in 2D simulations than in comparable 1D models. We conclude by discussing an extension of our method toward rotating LDI wind models that exhibit an intriguing combination of large- and small-scale structure extending down to the wind base.

(abridged) The strong winds of Wolf-Rayet (WR) stars are important for the mechanical and chemical feedback of the most massive stars and determine whether they end their lives as neutron stars or black holes. In this work we investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). We connect stellar structure models for He stars with wind models for optically-thick winds and assess how both types of models can simultaneously fulfill their respective sonic-point conditions. Fixing the outer wind law and terminal wind velocity, we obtain unique solutions for the mass-loss rates of optically-thick, radiatively-driven winds of WR stars in the phase of core He-burning. The resulting mass-loss relations as a function of stellar parameters, agree well with previous empirical relations. Furthermore, we encounter stellar mass limits below which no continuous solutions exist. While these mass limits agree with observations of WR stars in the Galaxy, they are in conflict with observations in the LMC. While our results confirm in particular the slope of oft-used empirical mass-loss relations, they imply that only part of the observed WN population can be understood in the framework of the standard assumptions of a smooth transonic flow and compact stellar core. This means that alternative approaches, such as a clumped and inflated wind structure, or deviations from the diffusion limit at the sonic point may have to be invoked. Qualitatively, the existence of mass limits for the formation of WR-type winds may be relevant for the non-detection of low-mass WR stars in binary systems, which are believed to be progenitors of Type Ib/c supernovae. The sonic-point conditions derived in this work may provide a possibility to include optically-thick winds in stellar evolution models in a more physically motivated form than in current models.

We derive semi-analytic solutions for optically thick, super-Eddington stellar winds, induced by an assumed steady energy addition $\Delta {\dot E}$ concentrated around a near-surface heating radius $R$ in a massive star of central luminosity $L_\ast$. We show that obtaining steady wind solutions requires both that the resulting total luminosity $L_o = L_\ast + \Delta {\dot E}$ exceed the Eddington luminosity, $\Gamma_o \equiv L_o/L_{Edd} > 1$, and that the induced mass loss rate be such that the "photon-tiring" parameter $m \equiv {\dot M} GM/R L_o \le 1-1/\Gamma_o$, ensuring the luminosity is sufficient to overcome the gravitational potential $GM/R$. Our analysis unifies previous super-Eddington wind models that either: (1) assumed a direct radiative flux-driving without accounting for the advection of radiative enthalpy that can become important in such an optically thick flow; or (2) assumed that such super-Eddington outflows are adiabatic, neglecting the effects of the diffusive radiative flux. We show that these distinct models become applicable in the asymptotic limits of small vs. large values of $m \Gamma_o $, respectively. By solving the coupled differential equations for radiative diffusion and wind momentum, we obtain general solutions that effectively bridge the behaviours of these limiting models. Two key scaling results are for the terminal wind speed to escape speed, which is found to vary as $v_\infty^2/v_{esc}^2 = \Gamma_o/(1+m \Gamma_o) -1$, and for the final observed luminosity $L_{ obs}$, which for all allowed steady-solutions with $m < 1 - 1/\Gamma_o$, exceeds the Eddington luminosity, $L_{obs} > L_{Edd}$. Our super-Eddington wind solutions have potential applicability for modeling phases of eruptive mass loss from massive stars, classical novae, and the remnants of stellar mergers.

OB stars exhibit various types of spectral variability historically associated with wind structures, including the apparently ubiquitous discrete absorption components (DACs). These features have been proposed to be caused either by magnetic fields or non-radial pulsations. In this second paper of this series, we revisit the canonical phenomenological hydrodynamical modelling used to explain the formation of DACs by taking into account modern observations and more realistic theoretical predictions. Using constraints on putative bright spots located on the surface of the O giant $\xi$ Persei derived from high precision space-based broadband optical photometry obtained with the Microvariability and Oscillations of STars (MOST) space telescope, we generate two-dimensional hydrodynamical simulations of co-rotating interaction regions in its wind. We then compute synthetic ultraviolet (UV) resonance line profiles using Sobolev Exact Integration and compare them with historical timeseries obtained by the International Ultraviolet Explorer (IUE) to evaluate if the observed behaviour of $\xi$ Persei's DACs is reproduced. Testing three different models of spot size and strength, we find that the classical pattern of variability can be successfully reproduced for two of them: the model with the smallest spots yields absorption features that are incompatible with observations. Furthermore, we test the effect of the radial dependence of ionization levels on line driving, but cannot conclusively assess the importance of this factor. In conclusion, this study self-consistently links optical photometry and UV spectroscopy, paving the way to a better understanding of cyclical wind variability in massive stars in the context of the bright spot paradigm.

Magnetic massive stars comprise approximately 10% of the total OB star population. Modern spectropolarimetry shows these stars host strong, stable, large-scale, often nearly dipolar surface magnetic fields of 1 kG or more. These global magnetic fields trap and deflect outflowing stellar wind material, forming an anisotropic magnetosphere that can be probed with wind-sensitive UV resonance lines. Recent HST UV spectra of NGC 1624-2, the most magnetic O star observed to date, show atypically unsaturated P-Cygni profiles in the Civ resonant doublet, as well as a distinct variation with rotational phase. We examine the effect of non-radial, magnetically-channeled wind outflow on P-Cygni line formation, using a Sobolev Exact Integration (SEI) approach for direct comparison with HST UV spectra of NGC 1624-2. We demonstrate that the addition of a magnetic field desaturates the absorption trough of the P-Cygni profiles, but further efforts are needed to fully account for the observed line profile variation. Our study thus provides a first step toward a broader understanding of how strong magnetic fields affect mass loss diagnostics from UV lines.

Recent spectropolarimetric surveys of bright, hot stars have found that ~10% of OB-type stars contain strong (mostly dipolar) surface magnetic fields (~kG). The prominent paradigm describing the interaction between the stellar winds and the surface magnetic field is the magnetically confined wind shock (MCWS) model. In this model, the stellar wind plasma is forced to move along the closed field loops of the magnetic field, colliding at the magnetic equator, and creating a shock. As the shocked material cools radiatively it will emit X-rays. Therefore, X-ray spectroscopy is a key tool in detecting and characterizing the hot wind material confined by the magnetic fields of these stars. Some B-type stars are found to have very short rotational periods. The effects of the rapid rotation on the X-ray production within the magnetosphere have yet to be explored in detail. The added centrifugal force due to rapid rotation is predicted to cause faster wind outflows along the field lines, leading to higher shock temperatures and harder X-rays. However, this is not observed in all rapidly rotating magnetic B-type stars. In order to address this from a theoretical point of view, we use the X-ray Analytical Dynamical Magnetosphere (XADM) model, originally developed for slow rotators, with an implementation of new rapid rotational physics. Using X-ray spectroscopy from ESA's XMM-Newton space telescope, we observed 5 rapidly rotating B-type stars to add to the previous list of observations. Comparing the observed X-ray luminosity and hardness ratio to that predicted by the XADM allows us to determine the role the added centrifugal force plays in the magnetospheric X-ray emission of these stars.

The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the wind mass-loss rates of the progenitors were greatly reduced relative to the mass-loss rates of massive stars in the Galaxy, concluding that heavy BHs must form in low-metallicity (Z < 0.25-0.5 Zsun) environments. However, strong surface magnetic fields also provide a powerful mechanism for modifying mass loss and rotation of massive stars, independent of environmental metallicity (ud-Doula & Owocki 2002; ud-Doula et al. 2008). In this paper we explore the hypothesis that some heavy BHs, with masses >25 Msun such as those inferred to compose GW150914, could be the natural end-point of evolution of magnetic massive stars in a solar-metallicity environment. Using the MESA code, we developed a new grid of single, non-rotating, solar metallicity evolutionary models for initial ZAMS masses from 40-80 Msun that include, for the first time, the quenching of the mass loss due to a realistic dipolar surface magnetic field. The new models predict TAMS masses that are significantly greater than those from equivalent non-magnetic models, reducing the total mass lost by a strongly magnetized 80 Msun star during its main sequence evolution by 20 Msun. This corresponds approximately to the mass loss reduction expected from an environment with metallicity Z = 1/30 Zsun.

Slowly rotating magnetic massive stars develop "dynamical magnetospheres" (DM's), characterized by trapping of stellar wind outflow in closed magnetic loops, shock heating from collision of the upflow from opposite loop footpoints, and subsequent gravitational infall of radiatively cooled material. In 2D and 3D magnetohydrodynamic (MHD) simulations the interplay among these three components is spatially complex and temporally variable, making it difficult to derive observational signatures and discern their overall scaling trends.Within a simplified, steady-state analysis based on overall conservation principles, we present here an "analytic dynamical magnetosphere" (ADM) model that provides explicit formulae for density, temperature and flow speed in each of these three components -- wind outflow, hot post-shock gas, and cooled inflow -- as a function of colatitude and radius within the closed (presumed dipole) field lines of the magnetosphere. We compare these scalings with time-averaged results from MHD simulations, and provide initial examples of application of this ADM model for deriving two key observational diagnostics, namely hydrogen H-alpha emission line profiles from the cooled infall, and X-ray emission from the hot post-shock gas. We conclude with a discussion of key issues and advantages in applying this ADM formalism toward derivation of a broader set of observational diagnostics and scaling trends for massive stars with such dynamical magnetospheres.

The detection by Rest et al. (2012) of light echoes from $\eta$ Carinae has provided important new observational constraints on the nature of its 1840's era giant eruption. Spectra of the echoes suggest a relatively cool spectral temperature of about 5500K, lower than the lower limit of about 7000K suggested in the optically thick wind outflow analysis of Davidson (1987). This has lead to a debate about the viability of this steady wind model relative to alternative, explosive scenarios. Here we present an updated analysis of the wind outflow model using newer low-temperature opacity tabulations and accounting for the stronger mass loss implied by the $>$10 Msun mass now inferred for the Homunculus. A major conclusion is that, because of the sharp drop in opacity due to free electron recombination for $T<$6500K, a low temperature of about 5000K is compatible with, and indeed expected from, a wind with the extreme mass loss inferred for the eruption. Within a spherical gray model in radiative equilibrium, we derive spectral energy distributions for various assumptions for the opacity variation of the wind, providing a basis for comparisons with observed light echo spectra. The scaling results here are also potentially relevant for other highly optically thick outflows, including those from classical novae, giant eruptions of LBVs and SN Type IIn precursors. A broader issue therefore remains whether the complex, variable features observed from such eruptions are better understood in terms of a steady or explosive paradigm, or perhaps a balance of these idealizations.

Many types of stars have strong magnetic fields that can dynamically influence the flow of circumstellar matter. In stars with accretion disks, the stellar magnetic field can truncate the inner disk and determine the paths that matter can take to flow onto the star. These paths are different in stars with different magnetospheres and periods of rotation. External field lines of the magnetosphere may inflate and produce favorable conditions for outflows from the disk-magnetosphere boundary. Outflows can be particularly strong in the propeller regime, wherein a star rotates more rapidly than the inner disk. Outflows may also form at the disk-magnetosphere boundary of slowly rotating stars, if the magnetosphere is compressed by the accreting matter. In isolated, strongly magnetized stars, the magnetic field can influence formation and/or propagation of stellar wind outflows. Winds from low-mass, solar-type stars may be either thermally or magnetically driven, while winds from massive, luminous O and B type stars are radiatively driven. In all of these cases, the magnetic field influences matter flow from the stars and determines many observational properties. In this chapter we review recent studies of accretion, outflows, and winds of magnetized stars with a focus on three main topics: (1) accretion onto magnetized stars; (2) outflows from the disk-magnetosphere boundary; and (3) winds from isolated massive magnetized stars. We show results obtained from global magnetohydrodynamic simulations and, in a number of cases compare global simulations with observations.

Magnetically confined winds of early-type stars are expected to be sources of bright and hard X-rays. To clarify the systematics of the observed X-ray properties, we have analyzed a large series of Chandra and XMM observations, corresponding to all available exposures of known massive magnetic stars (over 100 exposures covering ~60% of stars compiled in the catalog of Petit et al. 2013). We show that the X-ray luminosity is strongly correlated with the stellar wind mass-loss-rate, with a power-law form that is slightly steeper than linear for the majority of the less luminous, lower-Mdot B stars and flattens for the more luminous, higher-Mdot O stars. As the winds are radiatively driven, these scalings can be equivalently written as relations with the bolometric luminosity. The observed X-ray luminosities, and their trend with mass-loss rates, are well reproduced by new MHD models, although a few overluminous stars (mostly rapidly rotating objects) exist. No relation is found between other X-ray properties (plasma temperature, absorption) and stellar or magnetic parameters, contrary to expectations (e.g. higher temperature for stronger mass-loss rate). This suggests that the main driver for the plasma properties is different from the main determinant of the X-ray luminosity. Finally, variations of the X-ray hardnesses and luminosities, in phase with the stellar rotation period, are detected for some objects and they suggest some temperature stratification to exist in massive stars' magnetospheres.

The X-ray emission of $\eta$ Carinae shows multiple features at various spatial and temporal scales. The central constant emission (CCE) component is centred on the binary and arises from spatial scales much smaller than the bipolar Homunculus nebula, but likely larger than the central wind--wind collision region between the stars as it does not vary over the $\sim$2-3 month X-ray minimum when it can be observed. Using large-scale 3D smoothed particle hydrodynamics (SPH) simulations, we model both the colliding-wind region between the stars, and the region where the secondary wind collides with primary wind ejected from the previous periastron passage. The simulations extend out to one hundred semimajor axes and make two limiting assumptions (strong coupling and no coupling) about the influence of the primary radiation field on the secondary wind. We perform 3D radiative transfer calculations on the SPH output to synthesize the X-ray emission, with the aim of reproducing the CCE spectrum. For the preferred primary mass-loss rate $\dot{M}_A\approx8.5\times10^{-4}$ M$_\odot$ yr$^{-1}$, the model spectra well reproduce the observation as the strong- and no-coupling spectra bound the CCE observation for longitude of periastron $\omega\approx252^\circ$, and bound/converge on the observation for $\omega\approx90^\circ$. This suggests that $\eta$ Carinae has moderate coupling between the primary radiation and secondary wind, that both the region between the stars and the comoving collision on the backside of the secondary generate the CCE, and that the CCE cannot place constraints on the binary's line of sight. We also discuss comparisons with common X-ray fitting parameters.

The extreme luminosities of hot, massive stars drive strong stellar winds through UV line-scattering. For OB stars with an orbiting circumstellar disk, we explore the effect of such line-scattering in ablating disk material, initially focusing on the marginally optically thin decretion disks of classical Oe and Be stars. For this we apply a multi-dimensional radiation-hydrodynamics code, assuming optically thin ray tracing for the stellar continuum and a multi-ray Sobolev treatment of the line transfer. This accounts for desaturation of line-absorption by Keplerian shear in the disk, and associated driving by non-radial photons. Results show dense, intermediate-speed surface ablation, consistent with the strong, blue-shifted absorption seen in UV wind lines of Be shell stars. The asymptotic ablation rate is typically an order-unity factor times the stellar wind mass loss rate, leading to disk destruction times of order months to years for Be disks, consistent with observations. The much stronger radiative forces of O stars reduce this time to order days, making sustaining a disk difficult, and so providing a natural explanation for the rarity of Galactic Oe stars. Additionally, the weakened line-driving at lower metallicity implies both a reduction in the winds that help spin-down stars from near-critical rotation, and in the ablation of decretion disks, thus providing a natural explanation for the higher fraction of Classical Be stars, and the presence of Oe stars, in the Magellanic Clouds. We conclude with a discussion of future extensions to study line-driven ablation of denser, optically thick, accretion disks around pre-main-sequence massive stars.

The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is Delta kT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion star's pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source.

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.

Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescopes. In wide WR+O binaries, such as WR 140, the X-ray flux is tied to the orbital phase, and is a direct probe of the winds' properties. In the Galactic center, $\sim$30 WRs orbit the super massive black hole (SMBH) within $\sim$10", leading to a smorgasbord of wind-wind collisions. To model the X-ray emission of WR 140 and the Galactic center, we perform 3D hydrodynamic simulations to trace the complex gaseous flows, and then carry out 3D radiative transfer calculations to compute the variable X-ray spectra. The model WR 140 RXTE light curve matches the data well for all phases except the X-ray minimum associated with periastron, while the model spectra agree with the RXTE hardness ratio and the shape of the Suzaku observations throughout the orbit. The Galactic center model of the Chandra flux and spectral shape match well in the region r$<$3", but the model flux falls off too rapidly beyond this radius.

This work explores a Pulsationally Driven Orbital Mass Ejection (PDOME) model for the launching of Classical Be star disks. Under this model, a combination of rapid rotation and non-radial pulsation modes contribute to placing material into the circumstellar environment. Several varieties of non-radial pulsation modes, characterized by their propagation direction and the relative phase of their velocity and density perturbations, are considered. As well, the orbital stability of material launched by such a mechanism is investigated.

For line-driven winds from hot, luminous OB stars, we examine the subtle but important role of diffuse, scattered radiation in determining both the topology of steady-state solutions and intrinsic variability in the transonic wind base. We use a smooth source function formalism to obtain nonlocal, integral expressions for the direct and diffuse components of the line-force that account for deviations from the usual localized, Sobolev forms. As the scattering source function is reduced, we find the solution topology in the transonic region transitions from X-type, with a unique wind solution, to a nodal type, characterized by a degenerate family of solutions. Specifically, in the idealized case of an optically thin source function and a uniformly bright stellar disk, the unique X-type solution proves to be a stable attractor to which time-dependent numerical radiation-hydrodynamical simulations relax. But in models where the scattering strength is only modestly reduced, the topology instead turns nodal, with the associated solution degeneracy now manifest by intrinsic structure and variability that persist down to the photospheric wind base. This highlights the potentially crucial role of diffuse radiation for the dynamics and variability of line-driven winds, and seriously challenges the use of Sobolev theory in the transonic wind region. Since such Sobolev-based models are commonly used in broad applications like stellar evolution and feedback, this prompts development of new wind models, not only for further quantifying the intrinsic variability found here, but also for computing new theoretical predictions of global properties like velocity laws and mass-loss rates.

We observed NGC 1624-2, the O-type star with the largest known magnetic field Bp~20 kG), in X-rays with the ACIS-S camera onboard the Chandra X-ray Observatory. Our two observations were obtained at the minimum and maximum of the periodic Halpha emission cycle, corresponding to the rotational phases where the magnetic field is the closest to equator-on and pole-on, respectively. With these observations, we aim to characterise the star's magnetosphere via the X-ray emission produced by magnetically confined wind shocks. Our main findings are: (i) The observed spectrum of NGC 1624-2 is hard, similar to the magnetic O-type star Theta 1 Ori C, with only a few photons detected below 0.8 keV. The emergent X-ray flux is 30% lower at the Halpha minimum phase. (ii) Our modelling indicated that this seemingly hard spectrum is in fact a consequence of relatively soft intrinsic emission, similar to other magnetic Of?p stars, combined with a large amount of local absorption (~1-3 x 10^22 cm^-2). This combination is necessary to reproduce both the prominent Mg and Si spectral features, and the lack of flux at low energies. NGC 1624-2 is intrinsically luminous in X-rays (log LX emission ~ 33.4) but 70-95% of the X-ray emission produced by magnetically confined wind shocks is absorbed before it escapes the magnetosphere (log LX ISM corrected ~ 32.5). (iii) The high X-ray luminosity, its variation with stellar rotation, and its large attenuation are all consistent with a large dynamical magnetosphere with magnetically confined wind shocks.

Eta Carinae is the nearest example of a supermassive, superluminous, unstable star. Mass loss from the system is critical in shaping its circumstellar medium and in determining its ultimate fate. Eta Car currently loses mass via a dense, slow stellar wind and possesses one of the largest mass loss rates known. It is prone to episodes of extreme mass ejection via eruptions from some as-yet unspecified cause; the best examples of this are the large-scale eruptions which occurred in 19th century. Eta Car is a colliding wind binary in which strong variations in X-ray emission and in other wavebands are driven by the violent collision of the wind of eta Car-A and the fast, less dense wind of an otherwise hidden companion star. X-ray variations are the simplest diagnostic we have to study the wind-wind collision and allow us to measure the state of the stellar mass loss from both stars. We present the X-ray lightcurve over the last 20 years from ROSAT observations and monitoring with the Rossi X-ray Timing Explorer and the X-ray Telescope on the Swift satellite. We compare and contrast the behavior of the X-ray emission from the system over that timespan, including surprising variations during the 2014 X-ray minimum.

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

The initial success of the Rigidly Rotating Magnetosphere (RRM) model application to the B2Vp star sigma OriE by Townsend, Owocki & Groote (2005) triggered a renewed era of observational monitoring of this archetypal object. We utilize high-resolution spectropolarimetry and the magnetic Doppler imaging (MDI) technique to simultaneously determine the magnetic configuration, which is predominately dipolar, with a polar strength Bd = 7.3-7.8 kG and a smaller non-axisymmetric quadrupolar contribution, as well as the surface distribution of abundance of He, Fe, C, and Si. We describe a revised RRM model that now accepts an arbitrary surface magnetic field configuration, with the field topology from the MDI models used as input. The resulting synthetic Ha emission and broadband photometric observations generally agree with observations, however, several features are poorly fit. To explore the possibility of a photospheric contribution to the observed photometric variability, the MDI abundance maps were used to compute a synthetic photospheric light curve to determine the effect of the surface inhomogeneities. Including the computed photospheric brightness modulation fails to improve the agreement between the observed and computed photometry. We conclude that the discrepancies cannot be explained as an effect of inhomogeneous surface abundance. Analysis of the UV light variability shows good agreement between observed variability and computed light curves, supporting the accuracy of the photospheric light variation calculation. We thus conclude that significant additional physics is necessary for the RRM model to acceptably reproduce observations of not only sigma Ori E, but also other similar stars with significant stellar wind-magnetic field interactions.

In this paper we report 23 magnetic field measurements of the B3IV star HD 23478: 12 obtained from high resolution Stokes $V$ spectra using the ESPaDOnS (CFHT) and Narval (TBL) spectropolarimeters, and 11 from medium resolution Stokes $V$ spectra obtained with the DimaPol spectropolarimeter (DAO). HD 23478 was one of two rapidly rotating stars identified as potential "centrifugal magnetosphere" hosts based on IR observations from the Apache Point Observatory Galactic Evolution Experiment survey. We derive basic physical properties of this star including its mass ($M=6.1^{+0.8}_{-0.7}\,M_\odot$), effective temperature ($T_{\rm eff}=20\pm2\,$kK), radius ($R=2.7^{+1.6}_{-0.9}\,R_\odot$), and age ($\tau_{\rm age}=3^{+37}_{-1}\,$Myr). We repeatedly detect weakly-variable Zeeman signatures in metal, He and H lines in all our observations corresponding to a longitudinal magnetic field of $\langle B_z\rangle\approx-2.0\,$kG. The rotational period is inferred from Hipparcos photometry ($P_{\rm rot}=1.0498(4)\,$d). Under the assumption of the Oblique Rotator Model, our obsevations yield a surface dipole magnetic field of strength $B_d\geq9.5\,$kG that is approximately aligned with the stellar rotation axis. We confirm the presence of strong and broad H$\alpha$ emission and gauge the volume of this star's centrifugal magnetosphere to be consistent with those of other H$\alpha$ emitting centrifugal magnetosphere stars based on the large inferred Alfv\'en to Kepler radius ratio.

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.

We present results from a study of the eclipsing, colliding-wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMM-Newton observations. In addition, we have characterized the intrinsic, polarimetric phase-dependent behavior of the strongest optical emission lines using data obtained with the University of Wisconsin's Half-Wave Spectropolarimeter. We have detected evidence of the Coriolis distortion of the wind-wind collision in the X-ray regime, which manifests itself through asymmetric behavior around the eclipses in the system's X-ray light curves. The large opening angle of the X-ray emitting region, as well as its location (i.e. the WN wind does not collide with the O star, but rather its wind) are evidence of radiative braking/inhibition occurring within the system. Additionally, the polarimetric results show evidence of the cavity the wind-wind collision region carves out of the Wolf-Rayet star's wind.

The high luminosities of massive stars drive strong stellar winds, through line scattering of the star's continuum radiation. This paper reviews the dynamics of such line driving, building first upon the standard CAK model for steady winds, and deriving the associated analytic scalings for the mass loss rate and wind velocity law. It next summarizes the origin and nature of the strong Line Deshadowing Instability (LDI) intrinsic to such line-driving, including also the role of a diffuse-line-drag effect that stabilizes the wind base, and then describes how both instability and drag are incorporated in the Smooth Source Function (SSF) method for time-dependent simulations of the nonlinear evolution of the resulting wind structure. The review concludes with a discussion of the effect of the resulting extensive structure in temperature, density and velocity for interpreting observational diagnostics. In addition to the usual clumping effect on density-squared diagnostics, the spatial porosity of optically thick clumps can reduce single-density continuum absorption, and a kind of velocity porosity, or vorocity, can reduce the absorption strength of spectral lines. An overall goal is to illuminate the rich physics of radiative driving and the challenges that lie ahead in developing dynamical models for the often complex structure and variability of hot-star winds.

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 (Gayley 2014) to Chandra grating spectra from five O stars with X-ray emission indicative of embedded wind shocks in effectively single massive stars. Correcting for wind absorption of the X-ray line emission is a crucial component of our analysis, and we use wind optical depth values derived from X-ray line-profile fitting (Cohen et al. 2014) in order to make that correction. The shock-heating rate results we derive for all the stars are quite similar: the average wind mass element passes through roughly one shock that heats it to at least $10^6$ 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 \approx 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 $10^7$ K.

For 30 years, cyclical wind variability in OB stars has puzzled the astronomical community. Phenomenological models involving co-rotating bright spots provide a potential explanation for the observed variations, but the underlying physics remains unknown. We present recent results from hydrodynamical simulations constraining bright spot properties and compare them to what can be inferred from space-based photometry. We also explore the possibility that these spots are caused by magnetic fields and discuss the detectability of such fields.

[Abridged] We develop and benchmark a fast and easy-to-use effective-opacity formalism for line and continuum radiative transfer in an accelerating two-component clumpy medium. The formalism bridges the limits of optically thin and thick clumps, and is here used to i) design a simple vorosity-modified Sobolev with exact integration (vmSEI) method for analyzing UV wind resonance lines in hot, massive stars, and ii) derive simple correction factors to the line force driving the outflows of such stars. We show that (for a given ionization factor) UV resonance doublets may be used to analytically predict the upward corrections in empirically inferred mass-loss rates associated with porosity in velocity space (a.k.a. velocity-porosity, or vorosity), but that severe solution degeneracies exist. For an inter-clump density set to 1 % of the mean density, we for O and B supergiants derive upward empirical mass-loss corrections of typically factors of either ~5 or ~50, depending on which of the two applicable solutions is chosen. Overall, our results indicate this solution dichotomy severely limits the use of UV resonance lines as direct mass-loss indicators of clumped hot stellar winds. We next apply the effective-opacity formalism to the standard CAK theory of line-driven winds. By analytic and numerical hydrodynamics calculations, we show that in cases where vorosity is important at the critical point setting the mass-loss rate, the reduced line-force leads to a lower theoretical mass loss, by a factor scaling with the normalized velocity filling factor fvel. On the other hand, if vorosity is important only above this critical point, the predicted mass loss is not affected, but the wind terminal speed is reduced. This shows that porosity in velocity space can have a significant impact not only on the diagnostics, but also on the dynamics and theory of radiatively driven winds.

Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to $\sim$3000 km s$^{-1}$, and their collision leads to hot plasma (up to $\sim10^8$K) that emit thermal X-rays (up to $\sim$10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes' response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems $\eta$ Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only $\eta$ Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.

We use 2D MHD simulations to examine the effects of radiative cooling and inverse Compton (IC) cooling on X-ray emission from magnetically confined wind shocks (MCWS) in magnetic massive stars with radiatively driven stellar winds. For the standard dependence of mass loss rate on luminosity $\Mdot \sim L^{1.7} $, the scaling of IC cooling with $L$ and radiative cooling with $\Mdot$ means that IC cooling become formally more important for lower luminosity stars. However, because the sense of the trends is similar, we find the overall effect of including IC cooling is quite modest. More significantly, for stars with high enough mass loss to keep the shocks radiative, the MHD simulations indicate a linear scaling of X-ray luminosity with mass loss rate; but for lower luminosity stars with weak winds, X-ray emission is reduced and softened by a {\em shock retreat} resulting from the larger post-shock cooling length, which within the fixed length of a closed magnetic loop forces the shock back to lower pre-shock wind speeds. A semi-analytic scaling analysis that accounts both for the wind magnetic confinement and this shock retreat yields X-ray luminosities that have a similar scaling trend, but a factor few higher values, compared to time-averages computed from the MHD simulations. The simulation and scaling results here thus provide a good basis for interpreting available X-ray observations from the growing list of massive stars with confirmed large-scale magnetic fields.

The high luminosity of Very Massive Stars (VMS) means that radiative forces play an important, dynamical role both in the structure and stability of their stellar envelope, and in driving strong stellar-wind mass loss. Focusing on the interplay of radiative flux and opacity, with emphasis on key distinctions between continuum vs. line opacity, this chapter reviews instabilities in the envelopes and winds of VMS. Specifically, we discuss how: 1) the iron opacity bump can induce an extensive inflation of the stellar envelope; 2) the density dependence of mean opacity leads to strange mode instabilities in the outer envelope; 3) desaturation of line-opacity by acceleration of near-surface layers initiates and sustains a line-driven stellar wind outflow; 4) an associated line-deshadowing instability leads to extensive small-scale structure in the outer regions of such line-driven winds; 5) a star with super-Eddington luminosity can develop extensive atmospheric structure from photon bubble instabilities, or from stagnation of flow that exceeds the "photon tiring" limit; 6) the associated porosity leads to a reduction in opacity that can regulate the extreme mass loss of such continuum-driven winds. Two overall themes are the potential links of such instabilities to Luminous Blue Variable (LBV) stars, and the potential role of radiation forces in establishing the upper mass limit of VMS.

We quantitatively investigate the extent of wind absorption signatures in the X-ray grating spectra of all non-magnetic, effectively single O stars in the Chandra archive via line profile fitting. Under the usual assumption of a spherically symmetric wind with embedded shocks, we confirm previous claims that some objects show little or no wind absorption. However, many other objects do show asymmetric and blue shifted line profiles, indicative of wind absorption. For these stars, we are able to derive wind mass-loss rates from the ensemble of line profiles, and find values lower by an average factor of 3 than those predicted by current theoretical models, and consistent with H-alpha if clumping factors of f_cl ~ 20 are assumed. The same profile fitting indicates an onset radius of X-rays typically at r ~ 1.5 R_star, and terminal velocities for the X-ray emitting wind component that are consistent with that of the bulk wind. We explore the likelihood that the stars in the sample that do not show significant wind absorption signatures in their line profiles have at least some X-ray emission that arises from colliding wind shocks with a close binary companion. The one clear exception is zeta Oph, a weak-wind star that appears to simply have a very low mass-loss rate. We also reanalyse the results from the canonical O supergiant zeta Pup, using a solar-metallicity wind opacity model and find Mdot = 1.8 \times 10^{-6} M_sun/yr, consistent with recent multi-wavelength determinations.

We examine X-rays from radiatively cooled shocks, focusing on how their thin-shell instability reduces X-ray emission. For 2D simulations of collision between equal expanding winds, we carry out a parameter study of such instability as a function of the ratio of radiative vs. adiabatic-expansion cooling lengths. In the adiabatic regime, the extended cooling layer suppresses instability, leading to planar shock compression with X-ray luminosity that follows closely the expected (Lx ~ M^2) quadratic scaling with mass-loss rate M . In the strongly radiative limit, the X-ray emission now follows an expected linear scaling with mass loss (Lx ~ M), but the instability deforms the shock compression into extended shear layers with oblique shocks along fingers of cooled, dense material. The spatial dispersion of shock thermalization limits strong X-ray emission to the tips and troughs of the fingers, and so reduces the X-ray emission (here by about a factor 1/50) below what is expected from analytic radiative-shock models without unstable structure. Between these two limits, X-ray emission can switch between a high-state associated with extended shock compression, and a low-state characterized by extensive shear. Further study is needed to clarify the origin of this "shear mixing reduction factor" in X-ray emission, and its dependence on parameters like the shock Mach number.

The Magnetism in Massive Stars (MiMeS) project represents the largest systematic survey of stellar magnetism ever undertaken. Based on a sample of over 550 Galactic B and O-type stars, the MiMeS project has derived the basic characteristics of magnetism in hot, massive stars. Herein we report preliminary results.

Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling, and radiative forces. We find that Eta A's mass-loss rate greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star (Eta B) switches from the adiabatic to the radiative-cooling regime during periastron passage. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line-of-sight to the system lies close to the dense walls of the WWC zone. A modest decrease in Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

Massive, hot OB-stars show clear evidence of strong macroscopic broadening (in addition to rotation) in their photospheric spectral lines. This paper examines the occurrence of such "macro-turbulence" in slowly rotating O-stars with strong, organised surface magnetic fields. Focusing on the CIV 5811A line, we find evidence for significant macro-turbulent broadening in all stars except NGC1624-2, which also has (by far) the strongest magnetic field. Instead, the very sharp CIV lines in NGC1624-2 are dominated by magnetic Zeeman broadening, from which we estimate a dipolar field of approximately 20 kG. By contrast, magnetic broadening is negligible in the other stars (due to their weaker field strengths, on order 1 kG), and their CIV profiles are typically very broad and similar to corresponding lines observed in non-magnetic O-stars. Quantifying this by an isotropic, Gaussian macro-turbulence, we derive vmac = 2.2 (+- 0.9/2.2) km/s for NGC1624-2, and vmac = 20-65 km/s for the rest of the magnetic sample. We use these observational results to test the hypothesis that the field can stabilise the atmosphere and suppress the generation of macro-turbulence down to stellar layers where the magnetic pressure PB and the gas pressure Pg are comparable. Using a simple grey atmosphere to estimate the temperature T0 at which PB = Pg, we find that T0 > Teff for all investigated magnetic stars, but that T0 reaches the ~160000 K layers associated with the iron opacity-bump in hot stars only for NGC1624-2. This is consistent with the view that the responsible physical mechanism for photospheric O-star macro-turbulence may be stellar gravity-mode oscillations excited by sub-surface convection zones, and suggests that a sufficiently strong magnetic field can suppress such iron-bump generated convection and associated pulsational excitation.

We fit X-ray emission line profiles in high resolution XMM-Newton and Chandra grating spectra of the early O supergiant Zeta Pup with models that include the effects of porosity in the stellar wind. We explore the effects of porosity due to both spherical and flattened clumps. We find that porosity models with flattened clumps oriented parallel to the photosphere provide poor fits to observed line shapes. However, porosity models with isotropic clumps can provide acceptable fits to observed line shapes, but only if the porosity effect is moderate. We quantify the degeneracy between porosity effects from isotropic clumps and the mass-loss rate inferred from the X-ray line shapes, and we show that only modest increases in the mass-loss rate (<~ 40%) are allowed if moderate porosity effects (h_infinity <~ R_*) are assumed to be important. Large porosity lengths, and thus strong porosity effects, are ruled out regardless of assumptions about clump shape. Thus, X-ray mass-loss rate estimates are relatively insensitive to both optically thin and optically thick clumping. This supports the use of X-ray spectroscopy as a mass-loss rate calibration for bright, nearby O stars.

Even though it is broadly accepted that single Be stars are rapidly rotating stars surrounded by a flat rotating circumstellar disk, there is still a debate about how fast these stars rotate and also about the mechanisms involved in the angular-momentum and mass input in the disk. We study the properties of stars that rotate near their critical-rotation rate and investigate the properties of the disks formed by equatorial mass ejections. We used the most recent Geneva stellar evolutionary tracks for rapidly rotating stars that reach the critical limit and used a simple model for the disk structure. We obtain that for a 9 Msun star at solar metallicity, the minimum average velocity during the Main Sequence phase to reach the critical velocity is around 330 km/s, whereas it would be 390 km/s at the metallicity of the Small Magellanic Cloud (SMC). Red giants or supergiants originating from very rapid rotators rotate six times faster and show N/C ratios three times higher than those originating from slowly rotating stars. This difference becomes stronger at lower metallicity. It might therefore be very interesting to study the red giants in clusters that show a large number of Be stars on the MS band. On the basis of our single-star models, we show that the observed Be-star fraction with cluster age is compatible with the existence of a temperature-dependent lower limit in the velocity rate required for a star to become a Be star. The mass, extension, and diffusion time of the disks produced when the star is losing mass at the critical velocity, obtained from simple parametrized expressions, are not too far from those estimated for disks around Be-type stars. At a given metallicity, the mass and the extension of the disk increase with the initial mass and with age on the MS phase. Denser disks are expected in low-metallicity regions.

X-ray satellites since Einstein have empirically established that the X-ray luminosity from single O-stars scales linearly with bolometric luminosity, Lx ~ 10^{-7} Lbol. But straightforward forms of the most favored model, in which X-rays arise from instability-generated shocks embedded in the stellar wind, predict a steeper scaling, either with mass loss rate Lx ~ Mdot ~ Lbol^{1.7} if the shocks are radiative, or with Lx ~ Mdot^{2} ~ Lbol^{3.4} if they are adiabatic. This paper presents a generalized formalism that bridges these radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to the local radius. Noting that the thin-shell instability of radiative shocks should lead to extensive mixing of hot and cool material, we propose that the associated softening and weakening of the X-ray emission can be parametrized as scaling with the cooling length ratio raised to a power m$, the "mixing exponent". For physically reasonable values m ~= 0.4, this leads to an X-ray luminosity Lx ~ Mdot^{0.6} ~ Lbol that matches the empirical scaling. To fit observed X-ray line profiles, we find such radiative-shock-mixing models require the number of shocks to drop sharply above the initial shock onset radius. This in turn implies that the X-ray luminosity should saturate and even decrease for optically thick winds with very high mass-loss rates. In the opposite limit of adiabatic shocks in low-density winds (e.g., from B-stars), the X-ray luminosity should drop steeply with Mdot^2. Future numerical simulation studies will be needed to test the general thin-shell mixing ansatz for X-ray emission.

Building on results from the Magnetism in Massive Stars (MiMeS) project, this paper shows how a two-parameter classification of massive-star magnetospheres in terms of the magnetic wind confinement (which sets the Alfv\'en radius RA) and stellar rotation (which sets the Kepler co-rotation radius RK) provides a useful organisation of both observational signatures and theoretical predictions. We compile the first comprehensive study of inferred and observed values for relevant stellar and magnetic parameters of 64 confirmed magnetic OB stars with Teff > 16 kK. Using these parameters, we locate the stars in the magnetic confinement-rotation diagram, a log-log plot of RK vs. RA. This diagram can be subdivided into regimes of centrifugal magnetospheres (CM), with RA > RK, vs. dynamical magnetospheres (DM), with RK > RA. We show how key observational diagnostics, like the presence and characteristics of Halpha emission, depend on a star's position within the diagram, as well as other parameters, especially the expected wind mass-loss rates. In particular, we identify two distinct populations of magnetic stars with Halpha emission: namely, slowly rotating O-type stars with narrow emission consistent with a DM, and more rapidly rotating B-type stars with broader emission associated with a CM. For O-type stars, the high mass-loss rates are sufficient to accumulate enough material for line emission even within the relatively short free-fall timescale associated with a DM: this high mass-loss rate also leads to a rapid magnetic spindown of the stellar rotation. For the B-type stars, the longer confinement of a CM is required to accumulate sufficient emitting material from their relatively weak winds, which also lead to much longer spindown timescales. [abbreviated]

We present the first fully 3D MHD simulation for magnetic channeling and confinement of a radiatively driven, massive-star wind. The specific parameters are chosen to represent the prototypical slowly rotating magnetic O star \theta^1 Ori C, for which centrifugal and other dynamical effects of rotation are negligible. The computed global structure in latitude and radius resembles that found in previous 2D simulations, with unimpeded outflow along open field lines near the magnetic poles, and a complex equatorial belt of inner wind trapping by closed loops near the stellar surface, giving way to outflow above the Alfv\'{e}n radius. In contrast to this previous 2D work, the 3D simulation described here now also shows how this complex structure fragments in azimuth, forming distinct clumps of closed loop infall within the Alfv\'{e}n radius, transitioning in the outer wind to radial spokes of enhanced density with characteristic azimuthal separation of $15-20 \degr$. Applying these results in a 3D code for line radiative transfer, we show that emission from the associated 3D `dynamical magnetosphere' matches well the observed H\alpha emission seen from \theta^1 Ori C, fitting both its dynamic spectrum over rotational phase, as well as the observed level of cycle to cycle stochastic variation. Comparison with previously developed 2D models for Balmer emission from a dynamical magnetosphere generally confirms that time-averaging over 2D snapshots can be a good proxy for the spatial averaging over 3D azimuthal wind structure. Nevertheless, fully 3D simulations will still be needed to model the emission from magnetospheres with non-dipole field components, such as suggested by asymmetric features seen in the H\alpha equivalent-width curve of \theta^1 Ori C.

We investigate the effects of stellar limb-darkening and photospheric perturbations for the onset of wind structure arising from the strong, intrinsic line-deshadowing instability (LDI) of a line-driven stellar wind. A linear perturbation analysis shows that including limb-darkening reduces the stabilizing effect of the diffuse radiation, leading to a net instability growth rate even at the wind base. Numerical radiation-hydrodynamics simulations of the non-linear evolution of this instability then show that, in comparison with previous models assuming a uniformly bright star without base perturbations, wind structure now develops much closer ($r \la 1.1 R_\star$) to the photosphere. This is in much better agreement with observations of O-type stars, which typically indicate the presence of strong clumping quite near the wind base.

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.

We have compiled a list of 36 O+O and 89 Wolf-Rayet binary candidates in the Milky Way and Magellanic clouds detected with the Chandra, XMM-Newton and ROSAT satellites to probe the connection between their X-ray properties and their system characteristics. Of the WR binaries with published parameters, all but two have kT > 0.9 keV. The most X-ray luminous WR binaries are typically very long period systems. The WR binaries show a nearly four-order of magnitude spread in X-ray luminosity, even among among systems with very similar WR primaries. Among the O+O binaries, short-period systems generally have soft X-ray spectra and longer period systems show harder X-ray spectra, again with a large spread in LX/Lbol.

It is observationally as well as theoretically well established that the winds of hot, massive OB-stars are highly structured on a broad range of spatial scales. This paper first discusses consequences of the small-scale structures associated with the strong instability inherent to the line-driving of these winds. We demonstrate the importance of a proper treatment of such wind clumping to obtain reliable estimates of mass-loss rates, and also show that instability simulations that are perturbed at the lower boundary indeed display significant clumping quite close to the wind base, in general agreement with observations. But a growing subset of massive stars has also been found to possess strong surface magnetic fields, which may channel the star's outflow and induce also large-scale wind structures and cyclic behavior of spectral diagnostics. The paper concludes by showing that multi-dimensional, magneto-hydrodynamical wind simulations, together with detailed radiative-transfer modeling, can reproduce remarkably well the periodic Balmer line emission observed in slowly rotating magnetic O stars like HD191612.

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

The magnetic O-star HD191612 exhibits strongly variable, cyclic Balmer line emission on a 538-day period. We show here that its variable Halpha emission can be well reproduced by the rotational phase variation of synthetic spectra computed directly from full radiation magneto-hydrodynamical simulations of a magnetically confined wind. In slow rotators such as HD191612, wind material on closed magnetic field loops falls back to the star, but the transient suspension of material within the loops leads to a statistically overdense, low velocity region around the magnetic equator, causing the spectral variations. We contrast such "dynamical magnetospheres" (DMs) with the more steady-state "centrifugal magnetospheres" of stars with rapid rotation, and discuss the prospects of using this DM paradigm to explain periodic line emission from also other non-rapidly rotating magnetic massive stars.

(shortened) It has been proposed that the envelopes of luminous stars may be subject to substantial radius inflation. The inflation effect has been discussed in relation to the radius problem of WR stars, but has yet failed to explain the large observed radii of Galactic WR stars. We wish to obtain a physical perspective of the inflation effect, and study the consequences for the radii of WR stars, and LBVs. For WR stars the observed radii are up to an order of magnitude larger than predicted by theory, whilst S Doradus-type LBVs are subject to humongous radius variations, which remain as yet ill-explained. We use a dual approach to investigate the envelope inflation, based on numerical models for stars near the Eddington limit, and a new analytic formalism to describe the effect. An additional new aspect is that we take the effect of density inhomogeneities (clumping) within the outer stellar envelopes into account. Due to the effect of clumping we are able to bring the observed WR radii in agreement with theory. Based on our new formalism, we find that the radial inflation is a function of a dimensionless parameter W, which largely depends on the topology of the Fe-opacity peak, i.e., on material properties. For W>1, we discover an instability limit, for which the stellar envelope becomes gravitationally unbound, i.e. there no longer exists a static solution. Within this framework we are also able to explain the S Doradus-type instabilities for LBVs like AG Car, with a possible triggering due to changes in stellar rotation. The stellar effective temperatures in the upper HR diagram are potentially strongly affected by the inflation effect. This may have particularly strong effects on the evolved massive LBV and WR stars just prior to their final collapse, as the progenitors of SNe Ibc, SNe II, and long GRBs.

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA {\theta} that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38 degrees, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i \approx 130 to 145 degrees, {\theta} \approx -15 to +30 degrees, and an orbital axis projected on the sky at a PA \approx 302 to 327 degrees east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta B, thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.

We present a generalised formalism for treating the porosity-associated reduction in continuum opacity that occurs when individual clumps in a stochastic medium become optically thick. We consider geometries resulting in either isotropic or anisotropic effective opacity, and, in addition to an idealised model in which all clumps have the same local overdensity and scale, we also treat an ensemble of clumps with optical depths set by Markovian statistics. This formalism is then applied to the specific case of bound-free absorption of X-rays in hot star winds, a process not directly affected by clumping in the optically thin limit. We find that the Markov model gives surprisingly similar results to those found previously for the single clump model, suggesting that porous opacity is not very sensitive to details of the assumed clump distribution function. Further, an anisotropic effective opacity favours escape of X-rays emitted in the tangential direction (the `venetian blind' effect), resulting in a 'bump' of higher flux close to line centre as compared to profiles computed from isotropic porosity models. We demonstrate how this characteristic line shape may be used to diagnose the clump geometry, and we confirm previous results that for optically thick clumping to significantly influence X-ray line profiles, very large porosity lengths, defined as the mean free path between clumps, are required. Moreover, we present the first X-ray line profiles computed directly from line-driven instability simulations using a 3-D patch method, and find that porosity effects from such models also are very small. This further supports the view that porosity has, at most, a marginal effect on X-ray line diagnostics in O stars, and therefore that these diagnostics do indeed provide a good `clumping insensitive' method for deriving O star mass-loss rates.

Magnetic massive stars -- which are being discovered with increasing frequency -- represent a new category of wind-shaping mechanism for O and B stars. Magnetic channeling of these stars' radiation-driven winds, the Magnetically Confined Wind Shock paradigm, leads to the formation of a shock-heated magnetosphere, which can radiate X-rays, modify UV resonance lines, and create disks of Halpha emitting material. The dynamical properties of these magnetospheres are well understood from a theoretical point of view as an interplay between the magnetic wind confinement and rotation. However, the manifestations of magnetospheres across the spectrum may be more complex and various than first anticipated. On the other hand, recent advances in modeling these magnetospheres provide a key to better understand massive star winds in general. We will summarize the coordinated observational, theoretical, and modeling efforts from the Magnetism in Massive Star Project, addressing key outstanding questions regarding magnetosphere manifestations across the spectral domain.

The colliding wind binary (CWB) systems \eta\ Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these X-ray light curves and spectra, we apply 3D hydrodynamic simulations of the wind-wind collision using smoothed particle hydrodynamics (SPH), with the recent improvements of radiative cooling and the acceleration of the stellar winds according to a \beta\ law. For both systems, the 2-10 keV RXTE light curves are well-reproduced in absolute units for most phases, but the light curve dips associated with the periastron passages are not well matched. In WR140, the dip is too weak, and in \eta\ Carinae, the large difference in wind speeds of the two stars leads to a hot, post-periastron bubble that produces excess emission toward the end of the X-ray minimum.

The colliding wind binary (CWB) systems \eta\ Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these RXTE X-ray light curves, we model the wind-wind collision using 3D smoothed particle hydrodynamics (SPH) simulations. Adiabatic simulations that account for the absorption of X-rays from an assumed point source at the apex of the wind-collision shock cone by the distorted winds can closely match the observed 2-10keV RXTE light curves of both \eta\ Car and WR140. This point-source model can also explain the early recovery of \eta\ Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of \eta\ Car. Our more recent models relax the point-source approximation and account for the spatially extended emission along the wind-wind interaction shock front. For WR140, the computed X-ray light curve again matches the RXTE observations quite well. But for \eta\ Car, a hot, post-periastron bubble leads to an emission level that does not match the extended X-ray minimum observed by RXTE. Initial results from incorporating radiative cooling and radiatively-driven wind acceleration via a new anti-gravity approach into the SPH code are also discussed.

X-ray satellites since Einstein have empirically established that the X-ray luminosity from single O-stars scales linearly with bolometric luminosity, Lx ~ 10^-7 Lbol. But straightforward forms of the most favored model, in which X-rays arise from instability-generated shocks embedded in the stellar wind, predict a steeper scaling, either with mass loss rate Lx ~ Mdot ~ Lbol^1.7 if the shocks are radiative, or with Lx ~ Lx ~ Mdot^2 ~ Lbol^3.4 if they are adiabatic. We present here a generalized formalism that bridges these radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to the local radius. Noting that the thin-shell instability of radiative shocks should lead to extensive mixing of hot and cool material, we then propose that the associated softening and weakening of the X-ray emission can be parametrized by the cooling length ratio raised to a power m, the "mixing exponent." For physically reasonable values m ~= 0.4, this leads to an X-ray luminosity Lx ~ Mdot^0.6 ~ Lbol that matches the empirical scaling. We conclude by noting that such thin-shell mixing may also be important for X-rays from colliding wind binaries, and that future numerical simulation studies will be needed to test this thin-shell mixing ansatz for X-ray emission.

This review describes the evidence for small-scale structure, `clumping', in the radiation line-driven winds of hot, massive stars. In particular, we focus on examining to what extent simulations of the strong instability inherent to line-driving can explain the multitude of observational evidence for wind clumping, as well as on how to properly account for extensive structures in density and velocity when interpreting the various wind diagnostics used to derive mass-loss rates.

We have obtained 18 new high-resolution spectropolarimetric observations of the B2Vp star sigma Ori E with both the Narval and ESPaDOnS spectropolarimeters. The aim of these observations is to test, with modern data, the assumptions of the Rigidly Rotating Magnetosphere (RRM) model of Townsend & Owocki (2005), applied to the specific case of sigma Ori E by Townsend et al. (2005). This model includes a substantially offset dipole magnetic field configuration, and approximately reproduces previous observational variations in longitudinal field strength, photometric brightness, and Halpha emission. We analyze new spectroscopy, including H I, He I, C II, Si III and Fe III lines, confirming the diversity of variability in photospheric lines, as well as the double S-wave variation of circumstellar hydrogen. Using the multiline analysis method of Least-Squares Deconvolution (LSD), new, more precise longitudinal magnetic field measurements reveal a substantial variance between the shapes of the observed and RRM model time-varying field. The phase resolved Stokes V profiles of He I 5876 A and 6678 A lines are fit poorly by synthetic profiles computed from the magnetic topology assumed by Townsend et al. (2005). These results challenge the offset dipole field configuration assumed in the application of the RRM model to sigma Ori E, and indicate that future models of its magnetic field should also include complex, higher-order components.

We study the interaction between the Be star and the pulsar in the TeV binary PSR B1259-63/LS 2883, using 3-D SPH simulations of the tidal and wind interactions in this Be-pulsar system. We first run a simulation without pulsar wind nor Be wind, taking into account only the gravitational effect of the pulsar on the Be disk. In this simulation, the gas particles are ejected at a constant rate from the equatorial surface of the Be star, which is tilted in a direction consistent with multi-waveband observations. We run the simulation until the Be disk is fully developed and starts to repeat a regular tidal interaction with the pulsar. Then, we turn on the pulsar wind and the Be wind. We run two simulations with different wind mass-loss rates for the Be star, one for a B2V type and the other for a significantly earlier spectral type. Although the global shape of the interaction surface between the pulsar wind and the Be wind agrees with the analytical solution, the effect of the pulsar wind on the Be disk is profound. The pulsar wind strips off an outer part of the Be disk, truncating the disk at a radius significantly smaller than the pulsar orbit. Our results, therefore, rule out the idea that the pulsar passes through the Be disk around periastron, which has been assumed in the previous studies. It also turns out that the location of the contact discontinuity can be significantly different between phases when the pulsar wind directly hits the Be disk and those when the pulsar wind collides with the Be wind. It is thus important to adequately take into account the circumstellar environment of the Be star, in order to construct a satisfactory model for this prototypical TeV binary.

We present analysis of both the resolved X-ray emission line profiles and the broadband X-ray spectrum of the O2 If* star HD 93129A, measured with the Chandra HETGS. This star is among the earliest and most massive stars in the Galaxy, and provides a test of the embedded wind shock scenario in a very dense and powerful wind. A major new result is that continuum absorption by the dense wind is the primary cause of the hardness of the observed X-ray spectrum, while intrinsically hard emission from colliding wind shocks contributes less than 10% of the X-ray flux. We find results consistent with the predictions of numerical simulations of the line-driving instability, including line broadening indicating an onset radius of X-ray emission of several tenths Rstar. Helium-like forbidden-to-intercombination line ratios are consistent with this onset radius, and inconsistent with being formed in a wind-collision interface with the star's closest visual companion at a distance of ~100 AU. The broadband X-ray spectrum is fit with a dominant emission temperature of just kT = 0.6 keV along with significant wind absorption. The broadband wind absorption and the line profiles provide two independent measurements of the wind mass-loss rate: Mdot = 5.2_{-1.5}^{+1.8} \times 10^{-6} Msun/yr and Mdot = 6.8_{-2.2}^{+2.8} \times 10^{-6} Msun/yr, respectively. This is the first consistent modeling of the X-ray line profile shapes and broadband X-ray spectral energy distribution in a massive star, and represents a reduction of a factor of 3 to 4 compared to the standard H-alpha mass-loss rate that assumes a smooth wind.

The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent and ongoing star formation and to probe Carina's regions of bright diffuse X-ray emission. This study has yielded a catalog of properties of >14,000 X-ray point sources; >9800 of them have multiwavelength counterparts. Using Chandra's unsurpassed X-ray spatial resolution, we have separated these point sources from the extensive, spatially-complex diffuse emission that pervades the region; X-ray properties of this diffuse emission suggest that it traces feedback from Carina's massive stars. In this introductory paper, we motivate the survey design, describe the Chandra observations, and present some simple results, providing a foundation for the 15 papers that follow in this Special Issue and that present detailed catalogs, methods, and science results.

Mass-loss rates currently in use for hot, massive stars have recently been seriously questioned, mainly because of the effects of wind clumping. We investigate the impact of clumping on diagnostic ultraviolet resonance and optical recombination lines. Optically thick clumps, a non-void interclump medium, and a non-monotonic velocity field are all accounted for in a single model. We used 2D and 3D stochastic and radiation-hydrodynamic (RH) wind models, constructed by assembling 1D snapshots in radially independent slices. To compute synthetic spectra, we developed and used detailed radiative transfer codes for both recombination lines (solving the "formal integral") and resonance lines (using a Monte-Carlo approach). In addition, we propose an analytic method to model these lines in clumpy winds, which does not rely on optically thin clumping. Results: Synthetic spectra calculated directly from current RH wind models of the line-driven instability are unable to in parallel reproduce strategic optical and ultraviolet lines for the Galactic O-supergiant LCep. Using our stochastic wind models, we obtain consistent fits essentially by increasing the clumping in the inner wind. A mass-loss rate is derived that is approximately two times lower than predicted by the line-driven wind theory, but much higher than the corresponding rate derived from spectra when assuming optically thin clumps. Our analytic formulation for line formation is used to demonstrate the potential impact of optically thick clumping in weak-winded stars and to confirm recent results that resonance doublets may be used as tracers of wind structure and optically thick clumping. (Abridged)

We examine the nature and role of mass loss via an equatorial decretion disk in massive stars with near-critical rotation induced by evolution of the stellar interior. In contrast to the usual stellar wind mass loss set by exterior driving from the stellar luminosity, such decretion-disk mass loss stems from the angular momentum loss needed to keep the star near and below critical rotation, given the interior evolution and decline in the star's moment of inertia. Because the specific angular momentum in a Keplerian disk increases with the square root of the radius, the decretion mass loss associated with a required level of angular momentum loss depends crucially on the outer radius for viscous coupling of the disk, and can be significantly less than the spherical mass loss the spherical, wind-like mass loss commonly assumed in evolutionary calculations. We discuss the physical processes that affect the outer disk radius, including thermal disk outflow, and ablation of the disk material via a line-driven wind induced by the star's radiation. We present parameterized scaling laws for taking account of decretion-disk mass loss in stellar evolution codes, including how these are affected by metallicity, or by presence within a close binary and/or a dense cluster. Effects similar to those discussed here should also be present in accretion disks during star formation, and may play an important role in shaping the distribution of rotation speeds on the ZAMS.

There are currently two optical interferometry recombiners that can provide spectral resolutions better than 10000, AMBER/VLTI operating in the H-K bands, and VEGA/CHARA, recently commissioned, operating in the visible. These instruments are well suited to study the wind activity of the brightest AB supergiants in our vicinity, in lines such as H$\alpha$ or BrGamma. We present here the first observations of this kind, performed on Rigel (B8Ia) and Deneb (A2Ia). Rigel was monitored by AMBER in two campaigns, in 2006-2007 and 2009-2010, and observed in 2009 by VEGA; whereas Deneb was monitored in 2008-2009 by VEGA. The extension of the Halpha and BrGamma line forming regions were accurately measured and compared with CMFGEN models of both stars. Moreover, clear signs of activity were observed in the differential visibility and phases. These pioneer observations are still limited, but show the path for a better understanding of the spatial structure and temporal evolution of localized ejections using optical interferometry.

Small-scale inhomogeneities, or `clumping', in the winds of hot, massive stars are conventionally included in spectral analyses by assuming optically thin clumps. To reconcile investigations of different diagnostics using this microclumping technique, very low mass-loss rates must be invoked for O stars. Recently it has been suggested that by using the microclumping approximation one may actually drastically underestimate the mass-loss rates. Here we demonstrate this, present a new, improved description of clumpy winds, and show how corresponding models, in a combined UV and optical analysis, can alleviate discrepancies between previously derived rates and those predicted by the line-driven wind theory. Furthermore, we show that the structures obtained in time-dependent, radiation-hydrodynamic simulations of the intrinsic line-driven instability of such winds, which are the basis to our current understanding of clumping, in their present-day form seem unable to provide a fully self-consistent, simultaneous fit to both UV and optical lines. The reasons for this are discussed.

Perhaps the most extreme examples of "Active OB stars" are the subset of high-mass X-ray binaries -- consisting of an OB star plus compact companion -- that have recently been observed by Fermi and ground-based Cerenkov telescopes like HESS to be sources of very high energy (VHE; up to 30 TeV) gamma-rays. This paper focuses on the prominent gamma-ray source, LS5039, which consists of a massive O6.5V star in a 3.9-day-period, mildly elliptical (e = 0.24) orbit with its companion, assumed here to be a black-hole or unmagnetized neutron star. Using 3-D SPH simulations of the Bondi-Hoyle accretion of the O-star wind onto the companion, we find that the orbital phase variation of the accretion follows very closely the simple Bondi-Hoyle-Lyttleton (BHL) rate for the local radius and wind speed. Moreover, a simple model, wherein intrinsic emission of gamma-rays is assumed to track this accretion rate, reproduces quite well Fermi observations of the phase variation of gamma-rays in the energy range 0.1-10 GeV. However for the VHE (0.1-30 TeV) radiation observed by the HESS Cerenkov telescope, it is important to account also for photon-photon interactions between the gamma-rays and the stellar optical/UV radiation, which effectively attenuates much of the strong emission near periastron. When this is included, we find that this simple BHL accretion model also quite naturally fits the HESS light curve, thus making it a strong alternative to the pulsar-wind-shock models commonly invoked to explain such VHE gamma-ray emission in massive-star binaries.

This paper presents results obtained from Stokes I and V spectra of the B2Vp star sigma Ori E, observed by both the Narval and ESPaDOnS spectropolarimeters. Using Least- Squares Deconvolution, we investigate the longitudinal magnetic field at the current epoch, including period analysis exploiting current and historical data. sigma Ori E is the prototypical helium-strong star that has been shown to harbor a strong magnetic field, as well as a magnetosphere, consisting of two clouds of plasma forced by magnetic and centrifugal forces to co-rotate with the star on its 1.19 day period. The Rigidly Rotating Magnetosphere (RRM) model of Townsend & Owocki (2005) approximately reproduces the observed variations in longitudinal field strength, photometric brightness, Halpha emission, and various other observables. There are, however, small discrepancies between the observations and model in the photometric light curve, which we propose arise from inhomogeneous chemical abundances on the star's surface. Using Magnetic Doppler Imaging (MDI), future work will attempt to identify the contributions to the photometric variation due to abundance spots and due to circumstellar material.

The Magnetism in Massive Stars (MiMeS) Project is a consensus collaboration among many of the foremost international researchers of the physics of hot, massive stars, with the basic aim of understanding the origin, evolution and impact of magnetic fields in these objects. At the time of writing, MiMeS Large Programs have acquired over 950 high-resolution polarised spectra of about 150 individual stars with spectral types from B5-O4, discovering new magnetic fields in a dozen hot, massive stars. The quality of this spectral and magnetic mat\'eriel is very high, and the Collaboration is keen to connect with colleagues capable of exploiting the data in new or unforeseen ways. In this paper we review the structure of the MiMeS observing programs and report the status of observations, data modeling and development of related theory.

We present a method for computing the net transmission of X-rays emitted by shock-heated plasma distributed throughout a partially optically thick stellar wind from a massive star. We find the transmission by an exact integration of the formal solution, assuming that the emitting plasma and absorbing plasma are mixed at a constant mass ratio above some minimum radius, below which there is assumed to be no emission. This model is more realistic than either the slab absorption associated with a corona at the base of the wind or the exospheric approximation that assumes that all observed X-rays are emitted without attenuation from above the radius of optical depth unity. Our model is implemented in XSPEC as a pre-calculated table that can be coupled to a user-defined table of the wavelength dependent wind opacity. We provide a default wind opacity model that is more representative of real wind opacities than the commonly used neutral interstellar medium (ISM) tabulation. Preliminary modeling of \textit{Chandra} grating data indicates that the X-ray hardness trend of OB stars with spectral subtype can largely be understood as a wind absorption effect.

We analyze interferometric measurements of the Luminous Blue Variable Eta Carinae with the goal of constraining the rotational velocity of the primary star and probing the influence of the companion. Using 2-D radiative transfer models of latitude-dependent stellar winds, we find that prolate wind models with a ratio of the rotational velocity (vrot) to the critical velocity (vcrit) of W=0.77-0.92, inclination angle of i=60-90 degrees, and position angle PA=108-142 degrees reproduce simultaneously K-band continuum visibilities from VLTI/VINCI and closure phase measurements from VLTI/AMBER. Interestingly, oblate models with W=0.73-0.90 and i=80-90 degrees produce similar fits to the interferometric data, but require PA=210-230 degrees. Therefore, both prolate and oblate models suggest that the rotation axis of the primary star is not aligned with the Homunculus polar axis. We also compute radiative transfer models of the primary star allowing for the presence of a cavity and dense wind-wind interaction region created by the companion star. We find that the wind-wind interaction has a significant effect on the K-band image mainly via free-free emission from the compressed walls and, for reasonable model parameters, can reproduce the VLTI/VINCI visibilities taken at phase 0.92-0.93. We conclude that the density structure of the primary wind can be sufficiently disturbed by the companion, thus mimicking the effects of fast rotation in the interferometric observables. Therefore, fast rotation may not be the only explanation for the interferometric observations. Intense temporal monitoring and 3-D modeling are needed to resolve these issues.

Recent observations of luminous Type IIn supernovae (SNe) provide compelling evidence that massive circumstellar shells surround their progenitors. In this paper we investigate how the properties of such shells influence the SN lightcurve by conducting numerical simulations of the interaction between an expanding SN and a circumstellar shell ejected a few years prior to core collapse. Our parameter study explores how the emergent luminosity depends on a range of circumstellar shell masses, velocities, geometries, and wind mass-loss rates, as well as variations in the SN mass and energy. We find that the shell mass is the most important parameter, in the sense that higher shell masses (or higher ratios of M_shell/M_SN) lead to higher peak luminosities and higher efficiencies in converting shock energy into visual light. Lower mass shells can also cause high peak luminosities if the shell is slow or if the SN ejecta are very fast, but only for a short time. Sustaining a high luminosity for durations of more than 100 days requires massive circumstellar shells of order 10 M_sun or more. This reaffirms previous comparisons between pre-SN shells and shells produced by giant eruptions of luminous blue variables (LBVs), although the physical mechanism responsible for these outbursts remains uncertain. The lightcurve shape and observed shell velocity can help diagnose the approximate size and density of the circumstellar shell, and it may be possible to distinguish between spherical and bipolar shells with multi-wavelength lightcurves. These models are merely illustrative. One can, of course, achieve even higher luminosities and longer duration light curves from interaction by increasing the explosion energy and shell mass beyond values adopted here.

We present new U-band photometry of the magnetic Helium-strong star Sigma Ori E, obtained over 2004-2009 using the SMARTS 0.9-m telescope at Cerro Tololo Inter-American Observatory. When combined with historical measurements, these data constrain the evolution of the star's 1.19 d rotation period over the past three decades. We are able to rule out a constant period at the p_null = 0.05% level, and instead find that the data are well described (p_null = 99.3%) by a period increasing linearly at a rate of 77 ms per year. This corresponds to a characteristic spin-down time of 1.34 Myr, in good agreement with theoretical predictions based on magnetohydrodynamical simulations of angular momentum loss from magnetic massive stars. We therefore conclude that the observations are consistent with Sigma Ori E undergoing rotational braking due to its magnetized line-driven wind.

We report near-IR spectroscopic observations of the Eta Carinae massive binary system during 2008-2009 using VLT/CRIRES. We detect a strong, broad absorption wing in He I 10833 extending up to -1900 km/s across the 2009.0 spectroscopic event. Archival HST/STIS ultraviolet and optical data shows a similar high-velocity absorption (up to -2100 km/s) in the UV resonance lines of Si IV 1394, 1403 across the 2003.5 event. UV lines from low-ionization species, such as Si II 1527, 1533 and C II 1334, 1335, show absorption up to -1200 km/s, indicating that the absorption with v from -1200 to -2100 km/s originates in a region markedly faster and more ionized than the nominal wind of the primary star. Observations obtained at the OPD/LNA during the last 4 spectroscopic cycles (1989-2009) also display high-velocity absorption in He I 10833 during periastron. Based on the OPD/LNA dataset, we determine that material with v < -900 km/s is present in the phase range 0.976 < phi < 1.023 of the spectroscopic cycle, but absent in spectra taken at phi < 0.947 and phi > 1.049. Therefore, we constrain the duration of the high-velocity absorption to be 95 to 206 days (or 0.047 to 0.102 in phase). We suggest that the high-velocity absorption originates from shocked gas in the wind-wind collision zone, at distances of 15 to 45 AU in the line-of-sight to the primary star. Using 3-D hydrodynamical simulations of the wind-wind collision zone, we find that the dense high-velocity gas is in the line-of-sight to the primary star only if the binary system is oriented in the sky so that the companion is behind the primary star during periastron, corresponding to a longitude of periastron of omega ~ 240 to 270 degrees. We study a possible tilt of the orbital plane relative to the Homunculus equatorial plane and conclude that our data are broadly consistent with orbital inclinations in the range i=40 to 60 degrees.

We fit every emission line in the high-resolution Chandra grating spectrum of zeta Pup with an empirical line profile model that accounts for the effects of Doppler broadening and attenuation by the bulk wind. For each of sixteen lines or line complexes that can be reliably measured, we determine a best-fitting fiducial optical depth, tau_* = kappa*Mdot/4{pi}R_{\ast}v_{\infty}, and place confidence limits on this parameter. These sixteen lines include seven that have not previously been reported on in the literature. The extended wavelength range of these lines allows us to infer, for the first time, a clear increase in tau_* with line wavelength, as expected from the wavelength increase of bound-free absorption opacity. The small overall values of tau_*, reflected in the rather modest asymmetry in the line profiles, can moreover all be fit simultaneously by simply assuming a moderate mass-loss rate of 3.5 \pm 0.3 \times 10^{-6} Msun/yr, without any need to invoke porosity effects in the wind. The quoted uncertainty is statistical, but the largest source of uncertainty in the derived mass-loss rate is due to the uncertainty in the elemental abundances of zeta Pup, which affects the continuum opacity of the wind, and which we estimate to be a factor of two. Even so, the mass-loss rate we find is significantly below the most recent smooth-wind H-alpha mass-loss rate determinations for zeta Pup, but is in line with newer determinations that account for small-scale wind clumping. If zeta Pup is representative of other massive stars, these results will have important implications for stellar and galactic evolution.

Dimtri Mihalas' textbooks in the 70's and 80's on "Stellar Atmospheres" and "Foundations of Radiation Hydrodynamics" helped lay the early groundwork for understanding the moving atmospheres and winds of massive, luminous stars. Indeed, the central role of the momentum of stellar radiation in driving the mass outflow makes such massive-star winds key prototypes for radiation hydrodynamical processes. This paper reviews the dynamics of such radiative driving, building first upon the standard CAK model, and then discussing subtleties associated with the development and saturation of instabilities, and wind initiation near the sonic point base. An overall goal is to illuminate the rich physics of radiative driving and the challenges that lie ahead in developing dynamical models that can explain the broad scaling of mass loss rate and flow speed with stellar properties, as well as the often complex structure and variability observed in massive-star outflows.

In the subclass of high-mass X-ray binaries known as "microquasars", relativistic hadrons in the jets launched by the compact object can interact with cold protons from the star's radiatively driven wind, producing pions that then quickly decay into gamma rays. Since the resulting gamma-ray emissivity depends on the target density, the detection of rapid variability in microquasars with GLAST and the new generation of Cherenkov imaging arrays could be used to probe the clumped structure of the stellar wind. We show here that the fluctuation in gamma rays can be modeled using a "porosity length" formalism, usually applied to characterize clumping effects. In particular, for a porosity length defined by h=l/f, i.e. as the ratio of the characteristic size l of clumps to their volume filling factor f, we find that the relative fluctuation in gamma-ray emission in a binary with orbital separation a scales as sqrt(h/pi a) in the "thin-jet" limit, and is reduced by a factor 1/sqrt(1 + phi a/(2 l)) for a jet with a finite opening angle phi. For a thin jet and quite moderate porosity length h ~ 0.03 a, this implies a ca. 10 % variation in the gamma-ray emission. Moreover, the illumination of individual large clumps might result in isolated flares, as has been recently observed in some massive gamma-ray binaries.

We examine the angular momentum loss and associated rotational spin-down for magnetic hot stars with a line-driven stellar wind and a rotation-aligned dipole magnetic field. Our analysis here is based on our previous 2-D numerical MHD simulation study that examines the interplay among wind, field, and rotation as a function of two dimensionless parameters, W(=Vrot/Vorb) and 'wind magnetic confinement', $\eta_\ast$ defined below. We compare and contrast the 2-D, time variable angular momentum loss of this dipole model of a hot-star wind with the classical 1-D steady-state analysis by Weber and Davis (WD), who used an idealized monopole field to model the angular momentum loss in the solar wind. Despite the differences, we find that the total angular momentum loss averaged over both solid angle and time follows closely the general WD scaling $\dot {J} \sim \dot {M} \Omega R_A^2$. The key distinction is that for a dipole field Alfv\`en radius $R_A$ is significantly smaller than for the monopole field WD used in their analyses. This leads to a slower stellar spin-down for the dipole field with typical spin-down times of order 1 Myr for several known magnetic massive stars.

Stars can produce steady-state winds through radiative driving as long as the mechanical luminosity of the wind does not exceed the radiative luminosity at its base. This upper bound on the mass loss rate is known as the photon-tiring limit. Once above this limit, the radiation field is unable to lift all the material out of the gravitational potential of the star, such that only part of it can escape and reach infinity. The rest stalls and falls back toward the stellar surface, making a steady-state wind impossible. Photon-tiring is not an issue for line-driven winds since they cannot achieve sufficiently high mass loss rates. It can however become important if the star exceeds the Eddington limit and continuum interaction becomes the dominant driving mechanism. This paper investigates the time-dependent behaviour of stellar winds that exceed the photon-tiring limit, using 1-D numerical simulations of a porosity moderated, continuum-driven stellar wind. We find that the regions close to the star show a hierarchical pattern of high density shells moving back and forth, unable to escape the gravitational potential of the star. At larger distances, the flow eventually becomes uniformly outward, though still quite variable. Typically, these winds have a very high density but a terminal flow speed well below the escape speed at the stellar surface. Since most of the radiative luminosity of the star is used to drive the stellar wind, such stars would appear much dimmer than expected from the super-Eddington energy generation at their core. The visible luminosity typically constitutes less then half of the total energy flow and can become as low as ten percent or less for those stars that exceed the photon-tiring limit by a large margin.

LS 5039 is a TeV gamma-ray binary with extended radio emission. It consists of a compact object in the mildly eccentric (e=0.35), 3.9-day orbit around a massive O star. The nature of the compact object is not yet established. In this paper, assuming that the compact object is a black hole, we study the accretion of O-star wind by the black hole, by performing three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations. In order to roughly emulate the effect of the stellar radiation effectively canceling the stellar gravity, we assume that the O star's gravity does not exert on the wind. The wind particles are ejected with half the observed terminal velocity in a narrow range of azimuthal and vertical angles toward the black hole, in order to emulate the wind significantly slower than the terminal speed, and optimize the resolution and computational efficiency of simulations. We find that the mass-accretion rate closely follows the classical Bondi-Hoyle-Littleton accretion rate, which is of the order of 10^{16}g/s around periastron. The accretion rate at this level would provide jets enough power to produce the gamma-rays detected by HESS. Since the accretion peak occurs near the periastron passage, we need a strong gamma-ray absorption around periastron in order for the microquasar scenario to be consistent with the observed orbital modulation of the TeV gamma-ray flux.

We examine the angular momentum loss and associated rotational spindown for magnetic hot stars with a line-driven stellar wind and a rotation-aligned dipole magnetic field. Our analysis here is based on our previous 2-D numerical MHD simulation study that examines the interplay among wind, field, and rotation as a function of two dimensionless parameters, one characterizing the wind magnetic confinement ($\eta_{\ast} \equiv B_{eq}^{2} R_{\ast}^{2}/{\dot M} v_{\infty}$), and the other the ratio ($W \equiv V_{rot}/V_{orb}$) of stellar rotation to critical (orbital) speed. We compare and contrast the 2-D, time variable angular momentum loss of this dipole model of a hot-star wind with the classical 1-D steady-state analysis by Weber and Davis (WD), who used an idealized monopole field to model the angular momentum loss in the solar wind. Despite the differences, we find that the total angular momentum loss ${\dot J}$ averaged over both solid angle and time follows closely the general WD scaling ${\dot J} = (2/3) {\dot M} \Omega R_{A}^{2}$, where ${\dot M}$ is the mass loss rate, $\Omega$ is the stellar angular velocity, and $R_{A}$ is a characteristic Alfv\'{e}n radius. However, a key distinction here is that for a dipole field, this Alfv\'{e}n radius has a strong-field scaling $R_{A}/R_{\ast} \approx \eta_{\ast}^{1/4}$, instead of the scaling $R_{A}/R_{\ast} \sim \sqrt{\eta_{\ast}}$ for a monopole field. This leads to a slower stellar spindown time that in the dipole case scales as $\tau_{spin} = \tau_{mass} 1.5k/\sqrt{\eta_{\ast}}$, where $\tau_{mass} \equiv M/{\dot M}$ is the characteristic mass loss time, and $k$ is the dimensionless factor for stellar moment of inertia. The full numerical scaling relation we cite gives typical spindown times of order 1 Myr for several known magnetic massive stars.

We present the results of numerical simulations of continuum-driven winds of stars that exceed the Eddington limit and compare these against predictions from earlier analytical solutions. Our models are based on the assumption that the stellar atmosphere consists of clumped matter, where the individual clumps have a much larger optical thickness than the matter between the clumps. This `porosity' of the stellar atmosphere reduces the coupling between radiation and matter, since photons tend to escape through the more tenuous gas between the clumps. This allows a star that formally exceeds the Eddington limit to remain stable, yet produce a steady outflow from the region where the clumps become optically thin. We have made a parameter study of wind models for a variety of input conditions in order to explore the properties of continuum-driven winds. The results show that the numerical simulations reproduce quite closely the analytical scalings. The mass loss rates produced in our models are much larger than can be achieved by line driving. This makes continuum driving a good mechanism to explain the large mass loss and flow speeds of giant outbursts, as observed in eta Carinae and other luminous blue variable (LBV) stars. Continuum driving may also be important in population III stars, since line driving becomes ineffective at low metalicities. We also explore the effect of photon tiring and the limits it places on the wind parameters.

The very massive star system $\eta$ Carinae exhibits regular 5.54-year (2024-day) period disruptive events in wavebands ranging from the radio to X-ray. There is a growing consensus that these events likely stem from periastron passage of an (as yet) unseen companion in a highly eccentric ($\epsilon \sim 0.9$) orbit. This paper presents three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the orbital variation of the binary wind-wind collision, and applies these to modeling the X-ray light curve observed by the Rossi X-ray Timing Explorer (RXTE). By providing a global 3-D model of the phase variation of the density of the interacting winds, the simulations allow computation of the associated variation in X-ray absorption, presumed here to originate from near the apex of the wind-wind interaction cone. We find that the observed RXTE light curve can be readily fit if the observer's line of sight is within this cone along the general direction of apastron. Specifically, the data are well fit by an assumed inclination $i = 45^{\circ}$ for the orbit's polar axis, which is thus consistent with orbital angular momentum being along the inferred polar axis of the Homunculus nebula. The fits also constrain the position angle $\phi$ that an orbital-plane projection makes with the apastron side of the semi-major axis, strongly excluding positions $\phi < 9^{\circ}$ along or to the retrograde side of the axis, with the best fit position given by $\phi = 27^{\circ}$. Overall the results demonstrate the utility of a fully 3-D dynamical model for constraining the geometric and physical properties of this complex colliding-wind binary system.

We present results of a 90 ksec Chandra X-ray observation of the young sigma Orionis cluster (age ~3 Myr) obtained with the High Energy Transmission Grating Spectrometer. We use the high resolution grating spectrum and moderate resolution CCD spectrum of the massive central star sigma Ori AB (O9.5V + B0.5V) to test wind shock theories of X-ray emission and also analyze the high spatial resolution zero-order ACIS-S image of the central cluster region. Chandra detected 42 X-ray sources on the primary CCD (ACIS-S3). All but five have near-IR or optical counterparts and about one-fourth are variable. Notable high-mass stellar detections are sigma Ori AB, the magnetic B star sigma Ori E, and the B5V binary HD 37525. Most of the other detections have properties consistent with lower mass K or M-type stars. We present the first X-ray spectrum of the unusual infrared source IRS1 located 3.3 arc-sec north of sigma Ori AB, which is likely an embedded T Tauri star whose disk/envelope is being photoevaporated by sigma Ori AB. We focus on the radiative wind shock interpretation of the soft luminous X-ray emission from sigma Ori AB, but also consider possible alternatives including magnetically-confined wind shocks and colliding wind shocks. Its emission lines show no significant asymmetries or centroid shifts and are moderately broadened to HWHM ~ 264 km/s, or one-fourth the terminal wind speed. Forbidden lines in He-like ions are formally undetected, implying strong UV suppression. The Mg XI triplet forms in the wind acceleration zone within one stellar radius above the surface. These X-ray properties are consistent in several respects with the predictions of radiative wind shock theory for an optically thin wind, but explaining the narrow line widths presents a challenge to the theory.

We study colliding winds in the superluminous binary eta Carinae by performing three-dimensional, Smoothed Particle Hydrodynamics (SPH) simulations. For simplicity, we assume both winds to be isothermal. We also assume that wind particles coast without any net external forces. We find that the lower density, faster wind from the secondary carves out a spiral cavity in the higher density, slower wind from the primary. Because of the phase-dependent orbital motion, the cavity is very thin on the periastron side, whereas it occupies a large volume on the apastron side. The model X-ray light curve using the simulated density structure fits very well with the observed light curve for a viewing angle of i=54 degrees and phi=36 degrees, where i is the inclination angle and phi is the azimuth from apastron.

During the course of their evolution, massive stars lose a substantial fraction of their initial mass, both through steady winds and through relatively brief eruptions during their Luminous Blue Variable (LBV) phase. This talk reviews the dynamical driving of this mass loss, contrasting the line-driving of steady winds to the potential role of continuum driving for eruptions during LBV episodes when the star exceeds the Eddington limit. A key theme is to emphasize the inherent limits that self-shadowing places on line-driven mass loss rates, whereas continuum driving can in principle drive mass up to the "photon-tiring" limit, for which the energy to lift the wind becomes equal to the stellar luminosity. We review how the "porosity" of a highly clumped atmosphere can regulate continuum-driven mass loss, but also discuss recent time-dependent simulations of how base mass flux that exceeds the tiring limit can lead to flow stagnation and a complex, time-dependent combination of inflow and outflow regions. A general result is thus that porosity-mediated continuum driving in super-Eddington phases can explain the large, near tiring-limit mass loss inferred for LBV giant eruptions.

Building upon our previous MHD simulation study of magnetic channeling in radiatively driven stellar winds, we examine here the additional dynamical effects of stellar {\em rotation} in the (still) 2-D axisymmetric case of an aligned dipole surface field. In addition to the magnetic confinement parameter $\eta_{\ast}$ introduced in Paper I, we characterize the stellar rotation in terms of a parameter $W \equiv V_{\rm{rot}}/V_{\rm{orb}}$ (the ratio of the equatorial surface rotation speed to orbital speed), examining specifically models with moderately strong rotation $W =$ 0.25 and 0.5, and comparing these to analogous non-rotating cases. Defining the associated Alfv\'{e}n radius $R_{\rm{A}} \approx \eta_{\ast}^{1/4} \Rstar$ and Kepler corotation radius $R_{\rm{K}} \approx W^{-2/3} \Rstar$, we find rotation effects are weak for models with $R_{\rm{A}} < R_{\rm{K}}$, but can be substantial and even dominant for models with $R_{\rm{A}} \gtwig R_{\rm{K}}$. In particular, by extending our simulations to magnetic confinement parameters (up to $\eta_{\ast} = 1000$) that are well above those ($\eta_{\ast} = 10$) considered in Paper I, we are able to study cases with $R_{\rm{A}} \gg R_{\rm{K}}$; we find that these do indeed show clear formation of the {\em rigid-body} disk predicted in previous analytic models, with however a rather complex, dynamic behavior characterized by both episodes of downward infall and outward breakout that limit the buildup of disk mass. Overall, the results provide an intriguing glimpse into the complex interplay between rotation and magnetic confinement, and form the basis for a full MHD description of the rigid-body disks expected in strongly magnetic Bp stars like $\sigma$ Ori E.

Results from the most extensive study of the time-evolving dust structure around the prototype "Pinwheel" nebula WR 104 are presented. Encompassing 11 epochs in three near-infrared filter bandpasses, a homogeneous imaging data set spanning more than 6 years (or 10 orbits) is presented. Data were obtained from the highly successful Keck Aperture Masking Experiment, which can recover high fidelity images at extremely high angular resolutions, revealing the geometry of the plume with unprecedented precision. Inferred properties for the (unresolved) underlying binary and wind system are orbital period 241.5 +/- 0.5 days and angular outflow velocity of 0.28 +/- 0.02 mas/day. An optically thin cavity of angular size 13.3 +/- 1.4 mas was found to lie between the central binary and the onset of the spiral dust plume. Rotational motion of the wind system induced by the binary orbit is found to have important ramifications: entanglement of the winds results in strong shock activity far downstream from the nose of the bowshock. The far greater fraction of the winds participating in the collision may play a key role in gas compression and the nucleation of dust at large radii from the central binary and shock stagnation point. Investigation of the effects of radiative braking pointed towards significant modifications of the simple hydrostatic colliding wind geometry, extending the relevance of this phenomena to wider binary systems than previously considered. Limits placed on the maximum allowed orbital eccentricity of e < 0.06 argue strongly for a prehistory of tidal circularization in this system. Finally we discuss the implications of Earth's polar (i < 16 deg) vantage point onto a system likely to host supernova explosions at future epochs.

LS I +61 303 is a puzzling object detected from radio up to high-energy gamma-rays. Variability has recently been observed in its high-energy emission. The object is a binary system, with a compact object and a Be star as primary. The nature of the secondary and the origin of the gamma-ray emission are not clearly established at present. Recent VLBA radio data have been used to claim that the system is a Be/neutron star colliding wind binary, instead of a microquasar. We review the main views on the nature of LS I +61 303 and present results of 3D SPH simulations that can shed some light on the nature of the system. Our results support an accretion powered source, compatible with a microquasar interpretation.

We introduce a new Rigid-Field Hydrodynamics approach to modeling the magnetospheres of massive stars in the limit of very-strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the 1-dimensional flow along each, subject to pressure, radiative, gravitational, and centrifugal forces. We solve these equations numerically for a large ensemble of field lines, to build up a 3-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star sigma Ori E. Since the flow along each field line can be solved for independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment. The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disk defined by the locus of minima of the effective (gravitational plus centrifugal) potential. But a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to ~3,000 km/s) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disk can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.

Continuum driving is an effective method to drive a strong stellar wind. It is governed by two limits: the Eddington limit and the photon-tiring limit. A star must exceed the effective Eddington limit for continuum driving to overcome the stellar gravity. The photon-tiring limit places an upper limit on the mass loss rate that can be driven to infinity, given the energy available in the radiation field of the star. Since continuum driving does not require the presence of metals in the stellar atmosphere it is particularly suited to removing mass from low- and zero-metallicity stars and can play a crucial part in their evolution. Using a porosity length formalism we compute numerical simulations of super-Eddington, continuum driven winds to explore their behaviour for stars both below and above the photon-tiring limit. We find that below the photon tiring limit, continuum driving can produce a large, steady mass loss rate at velocities on the order of the escape velocity. If the star exceeds the photon-tiring limit, a steady solution is no longer possible. While the effective mass loss rate is still very large, the wind velocity is much smaller

LS I +61 303 is a puzzling Be/X-ray binary with variable gamma-ray emission at up TeV energies. The nature of the compact object and the origin of the high-energy emission are unclear. One family of models invokes particle acceleration in shocks from the collision between the B-star wind and a relativistic pulsar wind, while another centers on a relativistic jet powered by accretion. Recent high-resolution radio observations showing a putative "cometary tail" pointing away from the Be star near periastron have been cited as support for the pulsar-wind model. We wish here to carry out a quantitative assessment of these competing models for this extraordinary source. We apply a 3D SPH code for dynamical simulations of both the pulsar-wind-interaction and accretion-jet models. The former yields a description of the shape of the wind-wind interaction surface. The latter provides an estimation of the accretion rate. The results allow critical evaluation of how the two distinct models confront the data in various wavebands under a range of conditions. When one accounts for the 3D dynamical wind interaction under realistic constraints for the relative strength of the B-star and pulsar winds, the resulting form of the interaction front does not match the putative "cometary tail" claimed from radio observations. On the other hand, dynamical simulations of the accretion-jet model indicate that the orbital phase variation of accretion power includes a secondary broad peak well away from periastron, thus providing a plausible way to explain the observed TeV gamma ray emission toward apastron. We conclude that the colliding-wind model is not clearly established for LS I +61 303, while the accretion-jet model can reproduce many key characteristics of the observed TeV gamma-ray emission.

Gamma-rays can be produced by the interaction of a relativistic jet and the matter of the stellar wind in the subclass of massive X-ray binaries known as "microquasars". The relativistic jet is ejected from the surroundings of the compact object and interacts with cold protons from the stellar wind, producing pions that then quickly decay into gamma-rays. Since the resulting gamma-ray emissivity depends on the target density, the detection of rapid variability in microquasars with GLAST and the new generation of Cherenkov imaging arrays could be used to probe the clumped structure of the stellar wind. In particular, we show here that the relative fluctuation in gamma rays may scale with the square root of the ratio of porosity length to binary separation, \sqrt{h/a}, implying for example a ca. 10 % variation in gamma ray emission for a quite moderate porosity, h/a ~ 0.01.

We present XMM-Newton Reflection Grating Spectrometer observations of pairs of X-ray emission line profiles from the O star Zeta Pup that originate from the same He-like ion. The two profiles in each pair have different shapes and cannot both be consistently fit by models assuming the same wind parameters. We show that the differences in profile shape can be accounted for in a model including the effects of resonance scattering, which affects the resonance line in the pair but not the intercombination line. This implies that resonance scattering is also important in single resonance lines, where its effect is difficult to distinguish from a low effective continuum optical depth in the wind. Thus, resonance scattering may help reconcile X-ray line profile shapes with literature mass-loss rates.

We fit the Doppler profiles of the He-like triplet complexes of \ion{O}{7} and \ion{N}{6} in the X-ray spectrum of the O star $\zeta$ Puppis, using XMM-Newton RGS data collected over $\sim 400$ ks of exposure. We find that they cannot be well fit if the resonance and intercombination lines are constrained to have the same profile shape. However, a significantly better fit is achieved with a model incorporating the effects of resonance scattering, which causes the resonance line to become more symmetric than the intercombination line for a given characteristic continuum optical depth $\tau_*$. We discuss the plausibility of this hypothesis, as well as its significance for our understanding of Doppler profiles of X-ray emission lines in O stars.

We analyze the steady 1D flow equations for a rotating stellar wind based on a ``nozzle'' analogy for terms that constrain the local mass flux. For low rotation, we find the nozzle minimum occurs near the stellar surface, allowing a transition to a standard, CAK-type steep-acceleration solution; but for rotations > 75% of the critical rate, this inner nozzle minimum exceeds the global minimum, implying near-surface supercritical solutions would have an overloaded mass loss rate. In steady, analytic models in which the acceleration is assumed to be monotonically positive, this leads the solution to switch to a slow acceleration mode. However, time-dependent simulations using a numerical hydrodynamics code show that, for rotation rates 75 - 85% of critical, the flow can develop abrupt "kink" transitions from a steep acceleration to a decelerating solution. For rotations above 85% of critical, the hydrodynamic simulations confirm the slow acceleration, with the lower flow speed implying densities 5 - 30 times higher than the polar (or a nonrotating) wind. Still, when gravity darkening and 2D flow effects are accounted for, it seems unlikely that rotationally modified equatorial wind outflows could account for the very large densities inferred for the equatorial regions around B[e] supergiants.

We suggest that the mass lost during the evolution of very massive stars may be dominated by optically thick, continuum-driven outbursts or explosions, instead of by steady line-driven winds. In order for a massive star to become a WR star, it must shed its H envelope, but new estimates of the effects of clumping in winds indicate that line driving is vastly insufficient. We discuss massive stars above roughly 40-50 Msun, for which the best alternative is mass loss during brief eruptions of luminous blue variables (LBVs). Our clearest example of this phenomenon is the 19th century outburst of eta Car, when the star shed 12-20 Msun or more in less than a decade. Other examples are circumstellar nebulae of LBVs, extragalactic eta Car analogs (``supernova impostors''), and massive shells around SNe and GRBs. We do not yet fully understand what triggers LBV outbursts, but they occur nonetheless, and present a fundamental mystery in stellar astrophysics. Since line opacity from metals becomes too saturated, the extreme mass loss probably arises from a continuum-driven wind or a hydrodynamic explosion, both of which are insensitive to metallicity. As such, eruptive mass loss could have played a pivotal role in the evolution and fate of massive metal-poor stars in the early universe. If they occur in these Population III stars, such eruptions would profoundly affect the chemical yield and types of remnants from early SNe and hypernovae.

X-ray line profile analysis has proved to be the most direct diagnostic of the kinematics and spatial distribution of the very hot plasma around O stars. In this paper we apply several analysis techniques to the emission lines in the Chandra HETGS spectrum of the late-O supergiant zeta Ori (O9.7 Ib), including the fitting of a simple line-profile model. We show that there is distinct evidence for blue shifts and profile asymmetry, as well as broadening in the X-ray emission lines of zeta Ori. These are the observational hallmarks of a wind-shock X-ray source, and the results for zeta Ori are very similar to those for the earlier O star, zeta Pup, which we have previously shown to be well-fit by the same wind-shock line-profile model. The more subtle effects on the line-profile morphologies in zeta Ori, as compared to zeta Pup, are consistent with the somewhat lower density wind in this later O supergiant. In both stars, the wind optical depths required to explain the mildly asymmetric X-ray line profiles imply reductions in the effective opacity of nearly an order of magnitude, which may be explained by some combination of mass-loss rate reduction and large-scale clumping, with its associated porosity-based effects on radiation transfer. In the context of the recent reanalysis of the helium-like line intensity ratios in both zeta Ori and zeta Pup, and also in light of recent work questioning the published mass-loss rates in OB stars, these new results indicate that the X-ray emission from zeta Ori can be understood within the framework of the standard wind-shock scenario for hot stars.

We present results of numerical simulations of wind variability in the magnetic B1 IVe star $\beta$ Cephei. 2D-MHD simulations are used to determine the structure of the wind. From these wind models we calculate line profiles for different aspect angles to simulate rotation. The results are compared with the observed UV wind line profiles.

We present 2D MHD simulations of the radiatively driven outflow from a rotating hot star with a dipole magnetic field aligned with the star's rotation axis. We focus primarily on a model with moderately rapid rotation (half the critical value), and also a large magnetic confinement parameter, $\eta_{\ast} \equiv B_{\ast}^2 R_{\ast}^{2} / \dot{M} V_{\infty} = 600$. The magnetic field channels and torques the wind outflow into an equatorial, rigidly rotating disk extending from near the Kepler corotation radius outwards. Even with fine-tuning at lower magnetic confinement, none of the MHD models produce a stable Keplerian disk. Instead, material below the Kepler radius falls back on to the stellar surface, while the strong centrifugal force on material beyond the corotation escape radius stretches the magnetic loops outwards, leading to episodic breakout of mass when the field reconnects. The associated dissipation of magnetic energy heats material to temperatures of nearly $10^{8}$K, high enough to emit hard (several keV) X-rays. Such \emph{centrifugal mass ejection} represents a novel mechanism for driving magnetic reconnection, and seems a very promising basis for modeling X-ray flares recently observed in rotating magnetic Bp stars like $\sigma$ Ori E.

We investigate the degree to which the nearly symmetric form of X-ray emission lines seen in Chandra spectra of early-type supergiant stars could be explained by a possibly porous nature of their spatially structured stellar winds. Such porosity could effectively reduce the bound-free absorption of X-rays emitted by embedded wind shocks, and thus allow a more similar transmission of red- vs. blue-shifted emission from the back vs. front hemispheres. For a medium consisting of clumps of size l and volume filling factor f, in which the `porosity length' h=l/f increases with local radius as h = h' r, we find that a substantial reduction in wind absorption requires a quite large porosity scale factor h' > 1, implying large porosity lengths h > r. The associated wind structure must thus have either a relatively large scale l~ r, or a small volume filling factor f ~ l/r << 1, or some combination of these. The relatively small-scale, moderate compressions generated by intrinsic instabilities in line-driving seem unlikely to give such large porosity lengths, leaving again the prospect of instead having to invoke a substantial (ca. factor 5) downward revision in assumed mass-loss rates.

We attempt to characterize the observed variability of the magnetic helium-strong star sigma Ori E in terms of a recently developed rigidly rotating magnetosphere model. This model predicts the accumulation of circumstellar plasma in two co-rotating clouds, situated in magnetohydrostatic equilibrium at the intersection between magnetic and rotational equators. We find that the model can reproduce well the periodic modulations observed in the star's light curve, H alpha emission-line profile, and longitudinal field strength, confirming that it furnishes an essentially correct, quantitative description of the star's magnetically controlled circumstellar environment.

We report on four Chandra grating observations of the oblique magnetic rotator theta^1 Ori C (O5.5 V) covering a wide range of viewing angles with respect to the star's 1060 G dipole magnetic field. We employ line-width and centroid analyses to study the dynamics of the X-ray emitting plasma in the circumstellar environment, as well as line-ratio diagnostics to constrain the spatial location, and global spectral modeling to constrain the temperature distribution and abundances of the very hot plasma. We investigate these diagnostics as a function of viewing angle and analyze them in conjunction with new MHD simulations of the magnetically channeled wind shock mechanism on theta^1 Ori C. This model fits all the data surprisingly well, predicting the temperature, luminosity, and occultation of the X-ray emitting plasma with rotation phase.

We present initial attempts to include the multi-dimensional nature of radiation transport in hydrodynamical simulations of the small-scale structure that arises from the line-driven instability in hot-star winds. Compared to previous 1D or 2D models that assume a purely radial radiation force, we seek additionally to treat the lateral momentum and transport of diffuse line-radiation, initially here within a 2D context. A key incentive is to study the damping effect of the associated diffuse line-drag on the dynamical properties of the flow, focusing particularly on whether this might prevent lateral break-up of shell structures at scales near the lateral Sobolev angle of ca. $1^{\rm o}$. We first explore nonlinear simulations that cast the lateral diffuse force in the simple, local form of a parallel viscosity. Second, to account for the lateral mixing of radiation associated with the radial driving, we next explore models in which the radial force is azimuthally smoothed over a chosen scale. Third, to account for both the lateral line-drag and the lateral mixing in a more self-consistent way, we explore further a method first proposed by Owocki (1999), which uses a restricted 3-ray approach that combines a radial ray with two oblique rays set to have an impact parameter $p < R_{\ast}$ within the stellar core. From numerical simulations, we find that, compared to equivalent 1-ray simulations, the high-resolution 3-ray models show systematically a much higher lateral coherence.... (Full abstract in paper)

The migration of profile sub-peaks identified in time-monitored optical emission lines of Wolf-Rayet star spectra provides a direct diagnostic of the dynamics of their stellar winds via a measured line-of-sight velocity change per unit time. Inferring the associated wind acceleration scale from such an apparent acceleration then relies on the adopted intrinsic velocity of the wind material at the origin of this variable pattern. We develop radiative transfer tools to characterise the Line Emission Region (LER) of the program lines, including turbulence and line-optical depth effects. We find that monitored-lines can be fitted well with a pure optically thin formation mechanism, line-broadening resulting from the finite velocity extent of the LER rather than turbulence.Our new estimates of LER velocity centroids are systematically shifted outwards closer to terminal velocity compared to previous determinations, now suggesting WR-wind acceleration length scales of ca. 10-20Rsun, a factor of a few smaller than previously inferred. Based on radiation-hydrodynamics simulations of the line-driven-instability mechanism, we compute synthetic line profile variability for CIII5696A, for WR111 The results match well the measured observed migration of 20-30 m/s2, However, our model stellar radius of 19 Rsun, typical of an O-type supergiant, is a factor 2--10 larger than generally expected for WR core radii.Such small radii leave inferred acceleration scales to be more extended than expected from dynamical models of line driving, but the severity of the discrepancy is substantially reduced compared to previous analyses. We conclude with a discussion of how using lines formed deeper in the wind would provide a stronger constraint on the key wind dynamics as well as information on the onset of wind clumping.

We present a semi-analytical approach for modeling circumstellar emission from rotating hot stars with a strong dipole magnetic field tilted at an arbitrary angle to the rotation axis. By assuming the rigid-field limit in which material driven (e.g., in a wind outflow) from the star is forced to remain in strict rigid-body co-rotation, we are able to solve for the effective centrifugal-plus-gravitational potential along each field line, and thereby identify the location of potential minima where material is prone to accumulate. Applying basic scalings for the surface mass flux of a radiatively driven stellar wind, we calculate the circumstellar density distribution that obtains once ejected plasma settles into hydrostatic stratification along field lines. The resulting accumulation surface resembles a rigidly rotating, warped disk, tilted such that its average surface normal lies between the rotation and magnetic axes. Using a simple model of the plasma emissivity, we calculate time-resolved synthetic line spectra for the disk. Initial comparisons show an encouraging level of correspondence with the observed rotational phase variations of Balmer-line emission profiles from magnetic Bp stars like sigma Ori E.

New constraints on the large-scale wind structures responsible for discrete absorption components are obtained from far-ultraviolet time-series observations of O supergiants in the Large Magellanic Cloud.

We examine radiatively driven mass loss from stars near and above the Eddington limit (Ledd). We begin by reviewing the instabilities that are expected to form extensive structure near Ledd. We investigate how this "porosity" can reduce the effective coupling between the matter and radiation. Introducing a new "porosity-length'' formalism, we derive a simple scaling for the reduced effective opacity, and use this to derive an associated scaling for the porosity-moderated, continuum-driven mass loss rate from stars that formally exceed Ledd. For a simple super-Eddington model with a single porosity length that is assumed to be on the order of the gravitational scale height, the overall mass loss is similar to that derived in previous porosity work. This is much higher than is typical of line-driven winds, but is still only a few percent of the photon tiring limit--for which the luminosity becomes insufficient to carry the flow out of the gravitational potential. To obtain still stronger mass loss that approaches observationally inferred values near this limit, we introduce a power-law-porosity model in which the associated structure has a broad range of scales. We show that the mass loss rate can be enhanced by a factor that increases with the Eddington parameter Gamma, such that for moderately large Gamma (> 3-4), mass loss rates could approach the photon tiring limit. Together with the ability to drive quite fast outflow speeds, the derived mass loss could explain the large inferred mass loss and flow speeds of giant outbursts in eta Carinae and other LBV stars.

We present an efficient technique to study the 1D evolution of instability-generated structure in winds of hot stars out to very large distances (more than 1000 stellar radii). This technique makes use of our previous finding that external forces play little role in the outer evolution of structure. Rather than evolving the entire wind, as is traditionally done, the technique focuses on a representative portion of the structure and follows it as it moves out with the flow. This requires the problem to be formulated in a moving reference frame. The lack of Galilean invariance of the spherical equations of hydrodynamics is circumvented by recasting them in a pseudo-planar form. By applying the technique to a number of problems we show that it is fast and accurate, and has considerable conceptual advantages. It is particularly useful to test the dependence of solutions on the Galilean frame in which they were obtained. As an illustration, we show that, in a one-dimensional approximation, the wind can remain structured out to distances of more than 1300 stellar radii from the central star.

We report the analysis of an XMM-Newton observation of the close binary HD159176 (O7V + O7V). The observed L_X/L_bol ratio reveals an X-ray luminosity exceeding by a factor ~ 7 the expected value for X-ray emission from single O-stars, therefore suggesting a wind-wind interaction scenario. EPIC and RGS spectra are fitted consistently with a two temperature mekal optically thin thermal plasma model, with temperatures ranging from ~ 2 to 6 10^6 K. At first sight, these rather low temperatures are consistent with the expectations for a close binary system where the winds collide well before reaching their terminal velocities. We also investigate the variability of the X-ray light curve of HD159176 on various short time scales. No significant variability is found and we conclude that if hydrodynamical instabilities exist in the wind interaction region of HD159176, they are not sufficient to produce an observable signature in the X-ray emission. Hydrodynamic simulations using wind parameters from the literature reveal some puzzling discrepancies. The most striking one concerns the predicted X-ray luminosity which is one or more orders of magnitude larger than the observed one. A significant reduction of the mass loss rate of the components compared to the values quoted in the literature alleviates the discrepancy but is not sufficient to fully account for the observed luminosity. Because hydrodynamical models are best for the adiabatic case whereas the colliding winds in HD159176 are most likely highly radiative, a totally new approach has been envisaged... (see paper for complete abstract)

The Discrete Absorption Components (DACs) commonly observed in the ultraviolet lines of hot stars have previously been modelled by dynamical simulations of Corotating Interaction Regions (CIRs) in their line-driven stellar winds. Here we apply the kinematic DAC inversion method of Brown et al. to the hydrodynamical CIR models and test the reliability of the results obtained. We conclude that the inversion method is able to recover valuable information on the velocity structure of the mean wind and to trace movement of velocity plateaux in the hydrodynamical data, though the recovered density profile of the stream is correct only very near to the stellar surface.

We argue that, in general, observational studies of Be-star rotation have paid insufficient attention to the effects of equatorial gravity darkening. We present new line-profile calculations that emphasize the insensitivity of line width to rotation for fast rotators. Coupled with a critical review of observational procedures, these calculations suggest that the observational parameter v sin(i) may systematically underestimate the true projected equatorial rotation velocity, ve sin(i), by some tens of per cent for rapid rotators. The crucial implication of this work is that Be stars may be rotating much closer to their critical velocities than is generally supposed, bringing a range of new processes into contention for the elusive physical mechanism responsible for the circumstellar disk thought to be central to the Be phenomenon.

We perform an analytic investigation of the stability of line-driven disk winds, independent of hydrodynamic simulations. Our motive is to determine whether or not line-driven disk winds can account for the wide/broad UV resonance absorption lines seen in cataclysmic variables (CVs) and quasi-stellar objects (QSOs). In both CVs and QSOs observations generally indicate that the absorption arising in the outflowing winds has a steady velocity structure on time scales exceeding years (for CVs) and decades (for QSOs). However, published results from hydrodynamic simulations of line-driven disk winds are mixed, with some researchers claiming that the models are inherently unsteady, while other models produce steady winds. The analytic investigation presented here shows that if the accretion disk is steady, then the line-driven disk wind emanating from it can also be steady. In particular, we show that a gravitational force initially increasing along the wind streamline, which is characteristic of disk winds, does not imply an unsteady wind. The steady nature of line-driven disk winds is consistent with the 1D streamline disk-wind models of Murray and collaborators and the 2.5D time-dependent models of Pereyra and collaborators. This paper emphasizes the underlying physics behind the steady nature of line-driven disk winds using mathematically simple models that mimic the disk environment

This talk summarizes results from recent MHD simulations of the role of a dipole magnetic field in inducing large-scale structure in the line-driven stellar winds of hot, luminous stars. Unlike previous fixed-field analyses, the MHD simulations here take full account of the dynamical competition between the field and the flow. A key result is that the overall degree to which the wind is influenced by the field depends largely on a single, dimensionless `wind magnetic confinement parameter', $\eta_\ast (= B_{eq}^2 R_{\ast}^2/\dot{M} v_\infty$), which characterizes the ratio between magnetic field energy density and kinetic energy density of the wind. For weak confinement, $\eta_\ast \le 1$, the field is fully opened by wind outflow, but nonetheless, for confinement as small as $\eta_\ast=1/10$ it can have significant back-influence in enhancing the density and reducing the flow speed near the magnetic equator. For stronger confinement, $\eta_\ast > 1$, the magnetic field remains closed over limited range of latitude and height above the equatorial surface, but eventually is opened into nearly radial configuration at large radii. Within the closed loops, the flow is channeled toward loop tops into shock collisions that are strong enough to produce hard X-rays. Within the open field region, the equatorial channeling leads to oblique shocks that are again strong enough to produce X-rays and also lead to a thin, dense, slowly outflowing ``disk'' at the magnetic equator.

We carry out an extended analytic study of how the tilt and faster-than-radial expansion from a magnetic field affect the mass flux and flow speed of a line-driven stellar wind. A key motivation is to reconcile results of numerical MHD simulations with previous analyses that had predicted non-spherical expansion would lead to a strong speed enhancement. By including finite-disk correction effects, a dynamically more consistent form for the non-spherical expansion, and a moderate value of the line-driving power index $\alpha$, we infer more modest speed enhancements that are in good quantitative agreement with MHD simulations, and also are more consistent with observational results. Our analysis also explains simulation results that show the latitudinal variation of the surface mass flux scales with the square of the cosine of the local tilt angle between the magnetic field and the radial direction. Finally, we present a perturbation analysis of the effects of a finite gas pressure on the wind mass loss rate and flow speed in both spherical and magnetic wind models, showing that these scale with the ratio of the sound speed to surface escape speed, $a/v_{esc}$, and are typically 10-20% compared to an idealized, zero-gas-pressure model.

We summarize recent 2D MHD simulations of line-driven stellar winds from rotating hot-stars with a dipole magnetic field aligned to the star's rotation axis. For moderate to strong fields, much wind outflow is initially along closed magnetic loops that nearly corotate as a solid body with the underlying star, thus providing a torque that results in an effective angular momentum spin-up of the outflowing material. But instead of forming the ``magnetically torqued disk'' (MTD) postulated in previous phenemenological analyses, the dynamical simulations here show that material trapped near the tops of such closed loops tends either to fall back or break out, depending on whether it is below or above the Keplerian corotation radius. Overall the results raise serious questions about whether magnetic torquing of a wind outflow could naturally result in a Keplerian circumstellar disk. However, for very strong fields, it does still seem possible to form a %``magnetically confined, centrifugally supported, rigid-body disk'', centrifugally supported, ``magnetically rigid disk'' (MRD), in which the field not only forces material to maintain a rigid-body rotation, but for some extended period also holds it down against the outward centrifugal force at the loop tops. We argue that such rigid-body disks seem ill-suited to explain the disk emission from Be stars, but could provide a quite attractive paradigm for circumstellar emission from the magnetically strong Bp and Ap stars.

X-ray emission line profiles provide the most direct insight into the dynamics and spatial distribution of the hot, X-ray-emitting plasma above the surfaces of OB stars. The O supergiant zeta Puppis shows broad, blueshifted, and asymmetric line profiles, generally consistent with the wind-shock picture. We model the profiles of eight lines in the Chandra HETGS spectrum of this prototypical hot star. The fitted lines indicate that the plasma is distributed throughout the wind starting close to the photosphere, that there is significantly less attenuation of the X-rays by the overlying wind than is generally supposed, and that there is not a strong trend in wind absorption with wavelength.

We calculate multicomponent line-driven wind models of stars at extremely low metallicity suitable for massive first generation stars. For most of the models we find that the multicomponent wind nature is not important for either wind dynamics or for wind temperature stratification. However, for stars with the lowest metallicities we find that multicomponent effects influence the wind structure. These effects range from pure heating to possible fallback of the nonabsorbing wind component. We present a simple formula for the calculation of metallicity for which the multicomponent effects become important. We show that the importance of the multicomponent nature of winds of low metallicity stars is characterised not only by the low density of driving ions, but also by lower mass-loss rate.

Long known to be an unusual early-type star by virtue of its hard and strong X-ray emission, tau Scorpii poses a severe challenge to the standard picture of O-star wind-shock X-ray emission. The Chandra HETGS spectrum now provides significant direct evidence that this B0.2 star does not fit this standard wind-shock framework. The many emission lines detected with the Chandra gratings are significantly narrower than what would be expected from a star with the known wind properties of tau Sco, although they are broader than the corresponding lines seen in late-type coronal sources. While line ratios are consistent with the hot plasma on this star being within a few stellar radii of the photosphere, from at least one He-like complex there is evidence that the X-ray-emitting plasma is located more than a stellar radius above the photosphere. The Chandra spectrum of tau Sco is harder and more variable than those of other hot stars, with the exception of the young magnetized O star theta Ori C. We discuss these new results in the context of wind, coronal, and hybrid wind-magnetic models of hot-star X-ray emission.

The launch of high-spectral-resolution x-ray telescopes (Chandra, XMM) has provided a host of new spectral line diagnostics for the astrophysics community. In this paper we discuss Doppler-broadened emission line profiles from highly supersonic outflows of massive stars. These outflows, or winds, are driven by radiation pressure and carry a tremendous amount of kinetic energy, which can be converted to x rays by shock-heating even a small fraction of the wind plasma. The unshocked, cold wind is a source of continuum opacity to the x rays generated in the shock-heated portion of the wind. Thus the emergent line profiles are affected by transport through a two-component, moving, optically thick medium. While complicated, the interactions among these physical effects can provide quantitative information about the spatial distribution and velocity of the x-ray-emitting and absorbing plasma in stellar winds. We present quantitative models of both a spherically-symmetric wind and a wind with hot plasma confined in an equatorial disk by a dipole magnetic field.

We present numerical magnetohydrodynamic (MHD) simulations of the effect of stellar dipole magnetic fields on line-driven wind outflows from hot, luminous stars. Unlike previous fixed-field analyses, the simulations here take full account of the dynamical competition between field and flow, and thus apply to a full range of magnetic field strength, and within both closed and open magnetic topologies. A key result is that the overall degree to which the wind is influenced by the field depends largely on a single, dimensionless, `wind magnetic confinement parameter', $\eta_{\ast}$ ($ = B_{eq}^2 R_\ast^2/{\dot M} v_\infty$), which characterizes the ratio between magnetic field energy density and kinetic energy density of the wind. For weak confinement $\eta_{\ast} \le 1$, the field is fully opened by the wind outflow, but nonetheless for confinements as small as $\eta_{\ast}=1/10$ can have a significant back-influence in enhancing the density and reducing the flow speed near the magnetic equator. For stronger confinement $\eta_{\ast} > 1$, the magnetic field remains closed over a limited range of latitude and height about the equatorial surface, but eventually is opened into a nearly radial configuration at large radii.

High-resolution X-ray spectra of high-mass stars and low-mass T-Tauri stars obtained during the first year of the Chandra mission are providing important clues about the mechanisms which produce X-rays on very young stars. For zeta Puppis (O4 If) and zeta Ori (O9.5 I), the broad, blue-shifted line profiles, line ratios, and derived temperature distribution suggest that the X-rays are produced throughout the wind via instability-driven wind shocks. For some less luminous OB stars, like theta^1 Ori C (O7 V) and tau Sco (B0 V), the line profiles are symmetric and narrower. The presence of time-variable emission and very high-temperature lines in theta^1 Ori C and tau Sco suggest that magnetically confined wind shocks may be at work. The grating spectrum of the classical T-Tauri star TW Hya is remarkable because the forbidden-line emission of He-like Ne IX and O VII is very weak, implying that the X-ray emitting region is very dense, n = 6E+12 cgs, or that the X-rays are produced very close to the ultraviolet hotspot at the base of an accretion funnel. ACIS light curves and spectra of flares and low-mass and high-mass young stellar objects in Orion and rho Ophiuchus further suggest that extreme magnetic activity is a general property of many very young stars.

Motivated by recent detections by the XMM and Chandra satellites of X-ray line emission from hot, luminous stars, we present synthetic line profiles for X-rays emitted within parameterized models of a hot-star wind. The X-ray line emission is taken to occur at a sharply defined co-moving-frame resonance wavelength, which is Doppler-shifted by a stellar wind outflow parameterized by a `beta' velocity law, $v(r)=v_{\infty} (1-\R_{\ast}/r)^\beta$. Above some initial onset radius $R_o$ for X-ray emission, the radial variation of the emission filling factor is assumed to decline as a power-law in radius, $f(r) \sim r^{-q}$. The computed emission profiles also account for continuum absorption within the wind, with the overall strength characterized by a cumulative optical depth $\tau_\ast$. In terms of a wavelength shift from line-center scaled in units of the wind terminal speed $v_{\infty}$, we present normalized X-ray line profiles for various combinations of the parameters $\beta$, $\tau_\ast$, $q$ and $R_o$, and including also the effect of instrumental broadening as characterized by a Gaussian with a parameterized width $\sigma$. We discuss the implications for interpreting observed hot-star X-ray spectra, with emphasis on signatures for discriminating between ``coronal'' and ``wind-shock'' scenarios. In particular, we note that in profiles observed so far the substantial amount of emission longward of line center will be difficult to reconcile with the expected attenuation by the wind and stellar core in either a wind-shock or coronal model.

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 discuss the EUV and X-ray properties of B stars, focusing on $\epsilon$ CMa (B2 II) which is the only star with both emission lines and a photospheric continuum detected with EUVE. We explore the modest effects of the photospheric EUV continua on the wind, as well as the much stronger effects of the short-wavelength EUV and soft X-ray emission lines. Attenuation of the EUV and soft X-ray emission by the wind plays an important role, and leads to the reprocessing of X-rays via He$^+$ ionization and the Bowen mechanism in the wind. Finally, we explore some of the new diagnostics that will shortly become available with the next generation of high spectral resolution X-ray telescopes. All of this analysis is presented in the context of a two component stellar wind--a dense component (clumps) that contains most of the mass but fills a negligible fraction of the volume, and a rarefied component that fills most of the volume but accounts for only a small fraction of the mass.

We present a simple analysis of X-ray emission and absorption for hot-star winds, designed to explore the natural scalings of the observed X-ray luminosity with wind and sstellar properties. We show that an exospheric approximation, in which all of the emission above the optical depth unity radius escapes the wind, reproduces very well the detailed expression for radiation transport through a spherically symmetric wind. Using this approximation we find that the X-ray luminosity $L_x$ scales naturally with the wind density parameter $\Mdot/\vinf$, obtaining $L_x \sim (\Mdot/\vinf)^2$ for optically thin winds, and $L_x \sim (\Mdot/\vinf)^{1+s}$ for optically thick winds with an X-ray filling factor that varies in radius as $f \sim r^s$. These scalings with wind density contrast with the commonly inferred empirical scalings of X-ray luminosity $L_x $ with bolometric luminosity $L_B$. The empirically derived linear scaling of $L_x \sim L_B$ for thick winds can however be reproduced, through a delicate cancellation of emission and absorption, if one assumes modest radial fall-off in the X-ray filling factor ($s \approx -0.25$ or $s \approx -0.4$, depending on details of the secondary scaling of wind density with luminosity). We also explore the nature of the X-ray spectral energy distribution in the context of this model, and find that the spectrum is divided into a soft, optically thick part and a hard, optically thin part. Finally, we conclude that the energy-dependent emissivity must have a high-energy cut-off, corresponding to the maximum shock energy, in order to reproduce the general trends seen in X-ray spectral energy distributions of hot stars.

We present two-dimensional hydrodynamical simulations of corotating stream structure in the wind from a rotating O star, together with resulting synthetic line profiles showing discrete absorption components (DACs). An azimuthal variation is induced by a local increase or decrease in the radiative driving force, as would arise from a bright or dark ``star spot'' in the equatorial plane. Since much of the emergent wind structure seems independent of the exact method of perturbation, we expect similar morphology in winds perturbed by localized magnetic fields or nonradial pulsations, as well as by either rotationally-modulated structure or transient mass ejections. We find that bright spots with enhanced driving generate high-density, low-speed streams, while dark spots generate low-density, high-speed streams. Corotating interaction regions (CIRs) form where fast material collides with slow material -- e.g. at the leading (trailing) edge of a stream from a dark (bright) spot, often steepening into shocks. The unperturbed supersonic wind obliquely impacts the high-density CIR and sends back a nonlinear signal which takes the form of a sharp propagating discontinuity (``kink'' or ``plateau'') in the radial velocity gradient. These features travel inward in the co-moving frame at the radiative-acoustic characteristic speed, and thus slowly outward in the star's frame. We find that these slow kinks, rather than the CIRs themselves, are more likely to result in high-opacity DACs in the absorption troughs of unsaturated P Cygni line profiles.