In high mass X-ray binaries (HMXBs), an accreting compact object orbits a high mass star which loses mass through a dense and inhomogeneous wind. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass loss rates for high mass stars. We explore the impact of clumpiness on the variability of the column density with a simplified wind model. In particular, we focus on the standard deviation of the column density and the characteristic duration of enhanced absorption episodes, and compare them with analytical predictions based on the porosity length. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a clump in front of the compact object. We provide a recipe to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. The coherence time scale grants direct access to the size of the clumps while their mass can be deduced separately from the amplitude of the variability. If it is due to unaccreted passing-by clumps, the high column density variations in some HMXBs requires high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are hardly compatible with the ones derived from first principles. Alternatively, other components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism still to be identified, or a dense environment in the immediate vicinity of the accretor such as an accretion disk, an outflow or a spherical shell around the magnetosphere of the accreting neutron star.

Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensional. We present a new spectroscopic analysis tool tailored for massive stars that deviate from spherical symmetry. This code (entitled SPAMMS) is a spectroscopic patch model that takes the three-dimensional surface geometry of the system into account to produce spectral profiles at given phases and orientations. In using the Wilson-Devinney-like code PHOEBE in combination with the non local thermodynamic equilibrium (NLTE) radiative transfer code FASTWIND, we created a three-dimensional mesh that represents the surface geometry of our system and we assigned FASTWIND emergent intensity line profiles to each mesh triangle, which take the local parameters such as temperature, surface gravity, and radius into account. These line profiles were then integrated across the visible surface, where their flux contribution and radial velocity are taken into account, thus returning a final line profile for the visible surface of the system at a given phase. We demonstrate that SPAMMS can accurately reproduce the morphology of observed spectral line profiles for overcontact systems. Additionally, we show how line profiles of rapidly-rotating single stars differ when taking rotational distortion into account, and the effects that these can have on the determined parameters. Finally, we demonstrate the code's ability to reproduce the Rossiter-Mclaughlin and Struve-Sahade effects.

Context: Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and corresponding tools need to be developed. Aims: We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering spherically symmetric test models, we discuss potential error sources, and indicate advantages and disadvantages by comparing different solution methods. Further, we analyse the UV resonance line formation in the winds of rapidly rotating O stars. Methods: We consider both a (simplified) continuum model including scattering and thermal sources, and a UV resonance line transition approximated by a two-level-atom. We applied the short-characteristics (SC) method, using linear or monotonic B\'ezier interpolations, to solve the equation of radiative transfer on a non-uniform Cartesian grid. To calculate scattering dominated problems, our solution method is supplemented by an accelerated $\Lambda$-iteration scheme. Results: For the spherical test models, the mean relative error of the source function is on the $5-20\,\%$ level, depending on the applied interpolation technique and the complexity of the considered model. All calculated line profiles are in excellent agreement with corresponding 1D solutions. The predicted line profiles from fast rotating stars show a distinct behaviour as a function of rotational speed and inclination. This behaviour is tightly coupled to the wind structure and the description of gravity darkening and stellar surface distortion. Conclusions: Our SC methods are ready to be used for quantitative analyses of UV resonance line profiles. When calculating optically thick continua, both SC methods give reliable results, in contrast to the alternative finite-volume method.

(Abridged) The behaviour of mass loss across bi-stability jump is a key uncertainty in models of massive stars. While an increase in mass loss is theoretically predicted, this has so far not been observationally confirmed. However, radiation-driven winds of massive stars are known to exhibit clumpy structures triggered by the line-deshadowing instability (LDI). Wind clumping affects empirical mass-loss rates inferred from density square-dependent spectral diagnostics. If clumping properties differ significantly for O and B supergiants across the bi-stability jump, this may help alleviate discrepancies between theory and observations. We investigate with analytical and numerical tools how the onset of clumpy structures behaves in the winds of O supergiants (OSG) and B supergiants (BSG) across the bi-stability jump. We derive a scaling relation for the linear growth rate of the LDI for a single optically thick line and apply it in both regimes. We run 1D time-dependent line-driven instability simulations to study the non-linear evolution of the LDI in clumpy OSG and BSG winds. Linear perturbation analysis for a single line shows that the LDI linear growth rate scales strongly with stellar effective temperature and terminal wind speed. This implies significantly lower growth rates for (cooler, slower) BSG winds than for OSG winds. This is confirmed by the non-linear simulations, which show significant differences in OSG and BSG wind structure formation, with the latter characterized by significantly weaker clumping factors and lower velocity dispersions. This suggests that lower correction factors due to clumping should be employed when deriving empirical mass-loss rates for BSGs on the cool side of the bi-stability jump. Moreover, the non-linear simulations provide a theoretical background toward explaining the general lack of observed intrinsic X-ray emission in (single) B star winds.

[Abridged] Context: Radiation-driven mass loss plays a key role in the life-cycles of massive stars. However, basic predictions of such mass loss still suffer from significant quantitative uncertainties. Aims: We develop new radiation-driven, steady-state wind models for massive stars with hot surfaces, suitable for quantitative predictions of global parameters like mass-loss and wind-momentum rates. Methods: The simulations presented here are based on a self-consistent, iterative grid-solution to the spherically symmetric, steady-state equation of motion, using full NLTE radiative transfer solutions in the co-moving frame to derive the radiative acceleration. We do not rely on any distribution functions or parametrization for computation of the line force responsible for the wind driving. Results: In this first paper, we present models representing two prototypical O-stars in the Galaxy, one with a higher stellar mass M/Msun=59 and luminosity log10(L/Lsun)= 5.87 and one with a lower M/Msun= 27 and log10(L/Lsun)= 5.1. For these simulations, basic predictions for global mass-loss rates and velocity laws are given, and the influence from additional parameters like wind clumping and microturbulent speeds discussed. A key result is that our mass-loss rates are significantly lower than those predicted by the mass-loss recipes normally included in models of massive-star evolution. Conclusions: Our results support previous suggestions that Galactic O-star mass-loss rates may be overestimated in present-day stellar evolution models, and that new rates thus might be needed. Indeed, future papers in this series will incorporate our new models into such simulations of stellar evolution, extending the very first simulations presented here toward larger grids covering a range of metallicities, B supergiants across the bistability jump, and possibly also Wolf-Rayet stars.

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.

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.

A small subset of Galactic O-stars possess surface magnetic fields that alter the outflowing stellar wind by magnetically confining it. Key to the magnetic confinement is that it induces rotational modulation of spectral lines over the full EM domain; this allows us to infer basic quantities, e.g., mass-loss rate and magnetic geometry. Here, we present an empirical study of the H$\alpha$ line in Galactic magnetic O-stars to constrain the mass fed from the stellar base into the magnetosphere, using realistic multi-dimensional magnetized wind models, and compare with theoretical predictions. Our results suggest that it may be reasonable to use mass-feeding rates from non-magnetic wind theory if the absolute mass-loss rate is scaled down according to the amount of wind material falling back upon the stellar surface. This provides then some empirical support to the proposal that such magnetic O-stars might evolve into heavy stellar-mass black holes (Petit et al. 2017)

In Supergiant X-ray binaries (SgXB), a compact object captures a fraction of the wind of an O/B supergiant on a close orbit. Proxies exist to evaluate the efficiency of mass and angular momentum accretion but they depend so dramatically on the wind speed that given the current uncertainties, they only set loose constrains. Furthermore, they often bypass the impact of orbital and shock effects on the flow structure. We study the wind dynamics and the angular momentum gained as the flow is accreted. We identify the conditions for the formation of a disc-like structure around the accretor and the observational consequences for SgXB. We use recent results on the wind launching mechanism to compute 3D streamlines, accounting for the gravitational and X-ray ionizing influence of the compact companion on the wind. Once the flow enters the Roche lobe of the accretor, we solve the hydrodynamics equations with cooling. A shocked region forms around the accretor as the flow is beamed. For wind speeds of the order of the orbital speed, the shock is highly asymmetric compared to the axisymmetric bow shock obtained for a purely planar homogeneous flow. With net radiative cooling, the flow always circularizes for wind speeds low enough. Although the donor star does not fill its Roche lobe, the wind can be significantly beamed and bent by the orbital effects. The net angular momentum of the accreted flow is then sufficient to form a persistent disc-like structure. This mechanism could explain the proposed limited outer extension of the accretion disc in Cygnus X-1 and suggests the presence of a disc at the outer rim of the neutron star magnetosphere in Vela X-1, with dramatic consequences on the spinning up of the accretor.

Ultra-luminous X-ray sources (ULX) have so high X-ray luminosities that they were long thought to be accreting intermediate mass black holes. Yet, some ULX have been shown to display periodic modulations and coherent pulsations, suggestive of a neutron star in orbit around a stellar companion and accreting at super-Eddington rates. In this letter, we propose that the mass transfer in ULX could be qualitatively the same as in Supergiant X-ray binaries (SgXB), with a wind from the donor star highly beamed towards the compact object. Since the star does not fill its Roche lobe, this mass transfer mechanism known as "wind Roche lobe overflow" can remain stable even for large mass ratios. Based on realistic acceleration profiles derived from spectral observations and modeling of the stellar wind, we compute the bulk motion of the wind to evaluate the fraction of the stellar mass outflow captured by the compact object. We identify the orbital and stellar conditions for a SgXB to transfer mass at rates matching the expectations for ULX and show that the transition from SgXB to ULX luminosity levels is progressive. These results indicate that a high stellar Roche lobe filling factor is not necessary to funnel large quantities of material into the Roche lobe of the accretor. Large stellar mass loss rates such as the ones from the Wolf-Rayet star in M101 ULX-1 or the late B9 Supergiant in NGC 7793 P13 are enough to lead to a highly beamed wind and a significantly enhanced mass transfer rate.

Context: State of the art quantitative spectroscopy of OB-stars compares synthetic spectra (calculated by means of 1D, spherically symmetric computer codes) with observations. Certain stellar atmospheres, however, show strong deviations from spherical symmetry, and need to be treated in 3D. Aims: We present a newly developed 3D radiative transfer code, tailored to the solution of the radiation field in rapidly expanding stellar atmospheres. We apply our code to the continuum transfer in wind-ablation models, and to the UV resonance line formation in magnetic winds. Methods: We have used a 3D finite-volume method for the solution of the equation of radiative transfer, to study continuum- and line-scattering problems. Convergence has been accelerated by a non-local approximate Lambda-iteration scheme. Particular emphasis has been put on careful (spherically symmetric) test cases. Results: Typical errors of the source functions, when compared to 1D solutions, are of the order of 10-20 %, and increase for optically thick continua. In circumstellar discs, the radiation temperatures in the (optically thin) transition region from wind to disc are quite similar to corresponding values in the wind. For MHD simulations of dynamical magnetospheres, the line profiles, calculated with our 3D code, agree well with previous solutions using a 3D-SEI method. When compared with profiles resulting from the `analytic dynamical magnetosphere' (ADM) model, significant differences become apparent. Conclusions: Due to similar radiation temperatures in the wind and the transition region to the disc, the same line-strength distribution can be applied within radiation hydrodynamic calculations for circumstellar discs in `accreting high-mass stars'. To properly describe the UV line formation in dynamical magnetospheres, the ADM model needs to be further developed, at least in a large part of the outer wind.

We propose a methodological framework to perform forward asteroseismic modeling of stars with a convective core, based on gravity-mode oscillations. These probe the near-core region in the deep stellar interior. The modeling relies on a set of observed high-precision oscillation frequencies of low-degree coherent gravity modes with long lifetimes and their observational uncertainties. Identification of the mode degree and azimuthal order is assumed to be achieved from rotational splitting and/or from period spacing patterns. This paper has two major outcomes. The first is a comprehensive list and discussion of the major uncertainties of theoretically predicted gravity-mode oscillation frequencies based on linear pulsation theory, caused by fixing choices of the input physics for evolutionary models. Guided by a hierarchy among these uncertainties of theoretical frequencies, we subsequently provide a global methodological scheme to achieve forward asteroseismic modeling. We properly take into account correlations amongst the free parameters included in stellar models. Aside from the stellar mass, metalicity and age, the major parameters to be estimated are the near-core rotation rate, the amount of convective core overshooting, and the level of chemical mixing in the radiative zones. This modeling scheme allows for maximum likelihood estimation of the stellar parameters for fixed input physics of the equilibrium models, followed by stellar model selection considering various choices of the input physics. Our approach uses the Mahalanobis distance instead of the often used $\chi^2$ statistic and includes heteroscedasticity. It provides estimation of the unknown variance of the theoretically predicted oscillation frequencies.

[Abridged] Clumping in the radiation-driven winds of hot, massive stars affects the derivation of synthetic observables across the electromagnetic spectrum. We implement a formalism for treating wind clumping - in particular the light-leakage effects associated with a medium that is porous in physical and velocity space - into the global (photosphere+wind) NLTE model atmosphere code FASTWIND. We assume a stochastic, two-component wind consisting of a mixture of optically thick and thin clumps embedded in a rarefied inter-clump medium. We account fully for the reductions in opacity associated with porosity in physical and velocity-space, and for the well-known effect that opacities depending on rho^2 are higher in clumpy winds than in smooth ones of equal mass-loss rate. By formulating our method in terms of suitable mean and effective opacities for the clumpy wind, we are able to compute models with the same speed (~15 min. on a modern laptop) as in previous code-generations. Some first, generic results of the new models include: i) Confirming earlier results that velocity-space porosity is critical for analysis of UV wind lines in O-stars; ii) for the optical Halpha line, optically thick clumping effects are small for O-stars, but potentially very important for late B and A-supergiants; iii) spatial porosity is a marginal effect for absorption of high-energy X-rays in O-stars, as long as the mean-free path between clumps are kept at realistic values; iv) porosity is negligible at typical O-star radio-photosphere radii; v) regarding the wind ionization balance, a general trend is that increased rates of recombination in simulations with optically thin clumps lead to overall lower degrees of ionization than in corresponding smooth models, but that this effect now is counteracted by the increased levels of light-leakage associated with porosity in physical and velocity space.

(abridged) In the inner regions of AGB outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium (TE) chemical models. The presence of the majority of these species can be explained by shock-induced non-TE chemical models, where shocks caused by the pulsating star take the chemistry out of TE in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. The effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps, are examined. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming TE chemistry, such as HCN in O-rich outflows, H$_2$O in C-rich outflows, and NH$_3$ in both. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O- and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H$_2$O in C-rich outflows and of NH$_3$ in O- and C-rich outflows are however underpredicted.

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.

Supergiant X-ray Binaries (SgXB) host a compact object, often a neutron star (NS), orbiting an evolved O/B star. Mass transfer proceeds through the intense line-driven wind of the stellar donor, a fraction of which is captured by the gravitational field of the NS. The subsequent accretion process onto the NS is responsible for the abundant X-ray emission from SgXB. They also display peak-to-peak variability of the X-ray flux by a factor of a few 10 to 100, along with changes in the hardness ratios possibly due to varying absorption along the line-of-sight. We use recent radiation-hydrodynamic simulations of inhomogeneities (aka clumps) in the non-stationary wind of massive hot stars to evaluate their impact on the time-variable accretion process. For this, we run 3D hydrodynamic simulations of the wind in the vicinity of the accretor to investigate the formation of the bow shock and follow the inhomogeneous flow over several spatial orders of magnitude, down to the NS magnetosphere. In particular, we show that the impact of the wind clumps on the time-variability of the intrinsic mass accretion rate is severely tempered by the crossing of the shock, compared to the purely ballistic Bondi-Hoyle-Lyttleton estimation. We also account for the variable absorption due to clumps passing by the line-of-sight and estimate the final effective variability of the column density and mass accretion rate for different orbital separations. Finally, we compare our results to the most recent analysis of the X-ray flux and the hardness ratio in Vela X-1.

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.

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.

Massive stars, at least $\sim$ 10 times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after its Dark Ages. Another important property of massive stars are the strong stellar winds and outflows they produce. This mass loss, and finally the explosion of a massive star as a supernova or a gamma-ray burst, provide a significant input of mechanical and radiative energy into the interstellar space. These two properties together make massive stars one of the most important cosmic engines: they trigger the star formation and enrich the interstellar medium with heavy elements, that ultimately leads to formation of Earth-like rocky planets and the development of complex life. The study of massive star winds is thus a truly multidisciplinary field and has a wide impact on different areas of astronomy. [...] This detailed review summarises the current knowledge on the theory and observations of winds from massive stars, as well as on observations and accretion processes in wind-fed high mass X-ray binaries. The aim is to combine in the near future all available theoretical diagnostics and observational measurements to achieve a unified picture of massive star winds in isolated objects and in binary systems.

Strong winds from massive stars are a topic of interest to a wide range of astrophysical fields. In High-Mass X-ray Binaries the presence of an accreting compact object on the one side allows to infer wind parameters from studies of the varying properties of the emitted X-rays; but on the other side the accretor's gravity and ionizing radiation can strongly influence the wind flow. Based on a collaborative effort of astronomers both from the stellar wind and the X-ray community, this presentation attempts to review our current state of knowledge and indicate avenues for future progress.

We conducted a survey of seven magnetic O and eleven B-type stars with masses above $8M_{\odot}$ using the Very Large Array in the 1cm, 3cm and 13cm bands. The survey resulted in a detection of two O and two B-type stars. While the detected O-type stars - HD 37742 and HD 47129 - are in binary systems, the detected B-type stars, HD 156424 and ALS 9522, are not known to be in binaries. All four stars were detected at 3cm, whereas three were detected at 1cm and only one star was detected at 13cm. The detected B-type stars are significantly more radio luminous than the non-detected ones, which is not the case for O-type stars. The non-detections at 13cm are interpreted as due to thermal free-free absorption. Mass-loss rates were estimated using 3cm flux densities and were compared with theoretical mass-loss rates, which assume free-free emission. For HD 37742, the two values of the mass-loss rates were in good agreement, possibly suggesting that the radio emission for this star is mainly thermal. For the other three stars, the estimated mass-loss rates from radio observations were much higher than those expected from theory, suggesting either a possible contribution from non- thermal emission from the magnetic star or thermal or non-thermal emission due to interacting winds of the binary system, especially for HD 47129. All the detected stars are predicted to host centrifugal magnetospheres except HD 37742, which is likely to host a dynamical magnetosphere. This suggests that non-thermal radio emission is favoured in stars with centrifugal magnetospheres.

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.

X-rays/EUV radiation emitted from wind-embedded shocks in hot, massive stars can affect the ionization balance in their outer atmospheres, and can be the mechanism responsible for the production of highly ionized species. To allow for these processes in the context of spectral analysis, we have implemented such emission into our unified, NLTE model atmosphere/spectrum synthesis code FASTWIND. The shock structure and corresponding emission is calculated as a function of user-supplied parameters. We account for a temperature and density stratification inside the post-shock cooling zones, calculated for radiative and adiabatic cooling in the inner and outer wind, respectively. The high-energy absorption of the cool wind is considered by adding important K-shell opacities, and corresponding Auger ionization rates have been included into the NLTE network. We tested and verified our implementation carefully against corresponding results from various alternative model atmosphere codes, and studied the effects from shock emission for important ions from He, C, N, O, Si, and P. Surprisingly, dielectronic recombination turned out to play an essential role for the ionization balance of OIV/OV around Teff = 45,000 K. Finally, we investigated the behavior of the mass absorption coefficient, kappa_nu(r), important in the context of X-ray line formation in massive star winds. In almost all considered cases, direct ionization is of major influence, and Auger ionization significantly affects only NVI and OVI. The approximation of a radially constant kappa_nu is justified for r > 1.2 Rstar and lambda < 18 A, and also for many models at longer wavelengths. To estimate the actual value of this quantity, however, the HeII opacities need to be calculated from detailed NLTE modeling, at least for wavelengths longer than 18 to 20 A, and information on the individual CNO abundances has to be present.

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

The Of?p star CPD -28 2561 was monitored at high energies with XMM-Newton and HST. In X-rays, this magnetic oblique rotator displays bright and hard emission that varies by ~55% with rotational phase. These changes occur in phase with optical variations, as expected for magnetically confined winds; there are two maxima and two minima in X-rays during the 73d rotational period of CPD -28 2561. However, contrary to previously studied cases, no significant hardness variation is detected between minima and maxima, with the exception of the second minimum which is slightly distinct from the first one. In the UV domain, broad-band fluxes remain stable while line profiles display large variations. Stronger absorptions at low velocities are observed when the magnetic equator is seen edge-on, which can be reproduced by a detailed 3D model. However, a difference in absorption at high velocities in the CIV and NV lines is also detected for the two phases where the confined wind is seen nearly pole-on. This suggests the presence of strong asymmetries about the magnetic equator, mostly in the free-flowing wind (rather than in the confined dynamical magnetosphere).

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.

We have carried out a study of radio emission from a small sample of magnetic O- and B-type stars using the Giant Metrewave Radio Telescope, with the goal of investigating their magnetospheres at low frequencies. These are the lowest frequency radio measurements ever obtained of hot magnetic stars. The observations were taken at random rotational phases in the 1390 and the 610 MHz bands. Out of the 8 stars, we detect five B-type stars in both the 1390 and the 610 MHz bands. The O-type stars were observed only in the 1390 MHz band, and no detections were obtained. We explain this result as a consequence of free-free absorption by the free-flowing stellar wind exterior to the closed magnetosphere. We also study the variability of individual stars. One star - HD 133880 - exhibits remarkably strong and rapid variability of its low frequency flux density. We discuss the possibility of this emission being coherent emission as reported for CU Vir by Trigilio et al. (2000).

We report magnetic and spectroscopic observations and modeling of the Of?p star CPD -28 2561. Using more than 75 new spectra, we have measured the equivalent width variations and examined the dynamic spectra of photospheric and wind-sensitive spectral lines. A period search results in an unambiguous 73.41 d variability period. High resolution spectropolarimetric data analyzed using Least-Squares Deconvolution yield a Zeeman signature detected in the mean Stokes V profile corresponding to phase 0.5 of the spectral ephemeris. Interpreting the 73.41 d period as the stellar rotational period, we have phased the equivalent widths and inferred longitudinal field measurements. The phased magnetic data exhibit a weak sinusoidal variation, with maximum of about 565 G at phase 0.5, and a minimum of about -335 G at phase 0.0, with extrema approximately in phase with the (double-wave) Halpha equivalent width variation. Modeling of the Halpha equivalent width variation assuming a quasi-3D magnetospheric model produces a unique solution for the ambiguous couplet of inclination and magnetic obliquity angles: (i, beta) or (beta, i)=(35 deg,90 deg). Adopting either geometry, the longitudinal field variation yields a dipole polar intensity Bd=2.6\pm 0.9~kG, consistent with that obtained from direct modelling of the Stokes V profiles. We derive a wind magnetic confinement parameter eta*\simeq 100, leading to an Alfv\'en radius RA\simeq 3-5~R*, and a Kepler radius RK\simeq 20~R*. This supports a physical scenario in which the Halpha emission and other line variability have their origin in an oblique, co-rotating 'dynamical magnetosphere' structure resulting from a magnetically channeled wind. Nevertheless, the details of the formation of spectral lines and their variability within this framework remain generally poorly understood.

Binary systems with an accreting compact object are a unique chance to investigate the strong, clumpy, line-driven winds of early type supergiants by using the compact object's X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, the O9.7Iab giant companion of the black hole Cyg X-1 using 4.77 Ms of RXTE observations of the system taken over the course of 16 years. Absorption changes strongly over the 5.6 d binary orbit, but also shows a large scatter at a given orbital phase, especially at superior conjunction. The orbital variability is most prominent when the black hole is in the hard X-ray state. Our data are poorer for the intermediate and soft state, but show signs for orbital variability of the absorption column in the intermediate state. We quantitatively compare the data in the hard state to a toy model of a focussed Castor-Abbott-Klein-wind: as it does not incorporate clumping, the model does not describe the observations well. A qualitative comparison to a simplified simulation of clumpy winds with spherical clumps shows good agreement in the distribution of the equivalent hydrogen column density for models with a porosity length on the order of the stellar radius at inferior conjunction; we conjecture that the deviations between data and model at superior conjunction could be either due to lack of a focussed wind component in the model or a more complicated clump structure.

We review potential mass-loss mechanisms in the various evolutionary stages of massive stars, from the well-known line-driven winds of O-stars and BA-supergiants to the less-understood winds from Red Supergiants. We discuss optically thick winds from Wolf-Rayet stars and Very Massive Stars, and the hypothesis of porosity-moderated, continuum-driven mass loss from stars formally exceeding the Eddington limit, which might explain the giant outbursts from Luminous Blue Variables. We finish this review with a glance on the impact of rapid rotation, magnetic fields and small-scale inhomogeneities in line-driven 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.

A quite remarkable aspect of non-interacting O-stars with detected surface magnetic fields is that they all are very slow rotators. This paper uses this unique property to first demonstrate that the projected rotational speeds of massive, hot stars, as derived using current standard spectroscopic techniques, can be severely overestimated when significant "macroturbulent" line-broadening is present. This may, for example, have consequences for deriving the statistical distribution of rotation rates in massive-star populations, and for the use of these rates in stellar evolution models. It is next shown how such macroturbulence (seemingly a universal feature of hot, massive stars) is present in all but one of the magnetic O-stars, namely NGC 1624-2. Assuming then a simple model in which NGC 1624-2's exceptionally strong, large-scale magnetic field suppresses atmospheric motions down to layers where the magnetic and gas pressures are comparable, first empirical constraints on the formation depth of this enigmatic hot-star macroturbulence are derived. The results suggest an origin in the thin sub-surface convection zone of massive stars, consistent with a physical origin due to, e.g., stellar pulsations excited by the convective motions.

OB stars exhibit various types of spectral variability associated with wind structures, including the apparently ubiquitous discrete absorption components (DACs). These are proposed to be caused by either magnetic fields or non-radial pulsations (NRPs). In this paper, we evaluate the possible relation between large-scale, dipolar magnetic fields and the DAC phenomenon by investigating the magnetic properties of a sample of 13 OB stars exhibiting well-documented DAC behaviour. Using high-precision spectropolarimetric data acquired in part in the context of the Magnetism in Massive Stars (MiMeS) project, we find no evidence for surface dipolar magnetic fields in any of these stars. Using Bayesian inference, we compute upper limits on the strengths of the fields and use these limits to assess two potential mechanisms by which the field may influence wind outflow: magnetic wind confinement and local photospheric brightness enhancements. Within the limits we derive, both mechanisms fail to provide a systematic process capable of producing DACs in all of the stars of our sample. Therefore, this implies that dipolar fields are highly unlikely to be responsible for these structures in all massive stars, meaning that some other mechanism must come into play.

Important insights into the formation, structure, evolution and environment of all types of stars can be obtained through the measurement of their winds and possible magnetospheres. However, this has hardly been done up to now mainly because of the lack of UV instrumentation available for long periods of time. To reach this aim, we have designed UVMag, an M-size space mission equipped with a high-resolution spectropolarimeter working in the UV and visible spectral range. The UV domain is crucial in stellar physics as it is very rich in atomic and molecular lines and contains most of the flux of hot stars. Moreover, covering the UV and visible spectral domains at the same time will allow us to study the star and its environment simultaneously. Adding polarimetric power to the spectrograph will multiply tenfold the capabilities of extracting information on stellar magnetospheres, winds, disks, and magnetic fields. Examples of science objectives that can be reached with UVMag are presented for pre-main sequence, main sequence and evolved stars. They will cast new light onto stellar physics by addressing many exciting and important questions. UVMag is currently undergoing a Research and Technology study and will be proposed at the forthcoming ESA call for M-size missions. This spectropolarimeter could also be installed on a large UV and visible observatory (e.g. NASA's LUVOIR project) within a suite of instruments.

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

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.

Context: Rotation plays a key role in the life cycles of stars with masses above 8 Msun. Hence, accurate knowledge of the rotation rates of such massive stars is critical for understanding their properties and for constraining models of their evolution. Aims: This paper investigates the reliability of current methods used to derive projected rotation speeds v sin i from line-broadening signatures in the photospheric spectra of massive stars, focusing on stars that are not rapidly rotating. Methods: We use slowly rotating magnetic O-stars with well-determined rotation periods to test the Fourier transform (FT) and goodness-of-fit (GOF) methods typically used to infer projected rotation rates of massive stars. Results: For our two magnetic test stars with measured rotation periods longer than one year, i.e., with v sin i < 1 km/s, we derive v sin i ~ 40-50 km/s from both the FT and GOF methods. These severe overestimates are most likely caused by an insufficient treatment of the competing broadening mechanisms referred to as microturbulence and macroturbulence. Conclusions: These findings warn us not to rely uncritically on results from current standard techniques to derive projected rotation speeds of massive stars in the presence of significant additional line broadening, at least when v sin i <~ 50 km/s. This may, for example, be crucial for i) determining the statistical distribution of observed rotation rates of massive stars, ii) interpreting the evolutionary status and spin-down histories of rotationally braked B-supergiants, and iii) explaining the deficiency of observed O-stars with spectroscopically inferred v sin i ~ 0 km/s. Further investigations of potential shortcomings of the above techniques are presently under way.

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.

We present for the first time phase-resolved UV spectroscopy of an Of?p star, namely, HD 191612. The observations were acquired with the Space Telescope Imaging Spectrograph (STIS) on-board the Hubble Space Telescope (HST). We report the variability observed in the main photospheric and wind features and compare the results with previous findings for the Of?p star HD 108. We show that UV line strengths, H(alpha), and longitudinal magnetic field, vary coherently according to the rotational period (P = 537.6d), providing additional support for the magnetic oblique rotator scenario. The stellar and wind parameters of HD 191612 are obtained based on NLTE expanding atmosphere models. The peculiar wind line profile variations revealed by the new STIS data - not reproduced by 1D atmosphere models - are addressed through non-spherical MHD simulations coupled with radiative transfer. The basic aspects of the UV variability observed are explained and the structure of the dynamical magnetosphere of HD 191612 is 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 ~ 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.

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

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.

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.

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.

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

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)

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.

Massive stars are inherently extreme objects, in terms of radiation, mass loss, rotation, and sometimes also magnetic fields. Concentrating on a (personally biased) subset of processes related to pulsations, rapid rotation and its interplay with mass-loss, and the bi-stability mechanism, we will discuss how active (and normal) OB stars can serve as appropriate laboratories to provide further clues.

Small-scale clumping in the winds of hot, massive stars is conventionally included in spectral analyses by assuming optically thin clumps, a void inter-clump medium, and a smooth velocity field. To reconcile investigations of different diagnostics within such models, a highly clumped wind with very low mass-loss rates needs to be invoked. Particularly, unsaturated UV resonance lines seem to indicate rates an order of magnitude (or even more) lower than previously accepted values. We investigate resonance line formation in inhomogeneous hot star winds with non-monotonic velocity fields by means of 2D stochastic and pseudo-2D radiation-hydrodynamic wind models. A Monte-Carlo radiative transfer code is presented and used to produce synthetic line spectra. Results: The optically thin clumping limit is only valid for very weak lines. For intermediate strong lines, the velocity spans of the clumps are of central importance. Current hydrodynamical models predict spans that are too large to reproduce observed profiles unless a very low mass-loss rate is invoked. By simulating lower spans in 2D stochastic models, the profile strengths become drastically reduced, and are consistent with higher mass-loss rates. To simultaneously meet the constraints from strong lines, the inter-clump medium must be non-void. A first comparison to the observed PV doublet in the O6 supergiant lam Cep confirms that a stochastic 2D model reproduces observations with a mass-loss rate roughly ten times higher than that derived from the same lines but assuming optically thin clumping. Tentatively this may resolve discrepancies between theoretical predictions, evolutionary constraints, and recent derived mass-loss rates, and suggests a re-investigation of the structure predicted by current hydrodynamical models.

We review recent developments regarding radiation driven mass loss from OB-stars. We first summarize the fundamental theoretical predictions, and then compare these to observational results (including the VLT-FLAMES survey of massive stars). Especially we focus on the mass loss-metallicity dependence and on the so-called bi-stability jump. Subsequently we concentrate on two urgent problems, weak winds and wind clumping, that have been identified from various diagnostics and that challenge our present understanding of radiation driven winds. We discuss the problems of ``measuring'' mass-loss rates from weak winds and the potential of the near infrared, Br-alpha line as a tool to enable a more precise quantification, and comment on physical explanations for mass-loss rates that are much lower than predicted by the standard model. Wind clumping, conventionally interpreted as the consequence of a strong instability inherent to radiative line-driving, has severe implications for the interpretation of observational diagnostics, since derived mass-loss rates are usually overestimated when clumping is present but ignored in the analyses. Simplified techniques to account for clumping indicate overestimates by factors of 2 to 10, or even more. If actually true, these results would have a dramatic impact on the evolution of, and the feedback from, massive stars. We discuss ongoing attempts (including own work) to interprete the corresponding observations in terms of more sophisticated models. By allowing for porosity in density and velocity space, and for a non-void inter-clump medium, such models might require only moderate reductions of mass-loss rates.

The solar Mg I emission lines at 12 micrometer have already been observed and analyzed well. Previous modeling attempts for other stars have, however, been made only for Procyon and two cool evolved stars, with unsatisfactory results for the latter. We present high-resolution observational spectra for the K giants Pollux, Arcturus, and Aldebaran, which show strong Mg I emission lines at 12 micrometer as compared to the Sun. We also present the first observed stellar emission lines from Mg I at 18 micrometer and from Al I, Si I, and presumably Ca I at 12 micrometer. To produce synthetic line spectra, we employ standard non-LTE modeling for trace elements in cool stellar photospheres. We compute model atmospheres with the MARCS code, apply a comprehensive magnesium model atom, and use the radiative transfer code MULTI to solve for the magnesium occupation numbers in statistical equilibrium. We successfully reproduce the observed Mg I emission lines simultaneously in the giants and in the Sun, but show how the computed line profiles depend critically on atomic input data and how the inclusion of energy levels with n > 9 and collisions with neutral hydrogen are necessary to obtain reasonable fits.