A small fraction of the radiative flux emitted by hot stars is absorbed by their winds and redistributed towards longer wavelengths. This effect, which leads also to the heating of the stellar photosphere, is termed wind blanketing. For stars with variable winds, the effect of wind blanketing may lead to the photometric variability. We have studied the consequences of line driven wind instability and wind blanketing for the light variability of O stars. We combined the results of wind hydrodynamic simulations and of global wind models to predict the light variability of hot stars due to the wind blanketing and instability. The wind instability causes stochastic light variability with amplitude of the order of tens of millimagnitudes and a typical timescale of the order of hours for spatially coherent wind structure. The amplitude is of the order of millimagnitudes when assuming that the wind consists of large number of independent concentric cones. The variability with such amplitude is observable using present space borne photometers. We show that the simulated light curve is similar to the light curves of O stars obtained using BRITE and CoRoT satellites.

Context. Evolution of massive stars is affected by a significant loss of mass either via (nearly) spherically symmetric stellar winds or by aspherical mass-loss mechanisms, namely the outflowing equatorial disks. However, the scenario that leads to the formation of a disk or rings of gas and dust around massive stars is still under debate. Aims. We study the hydrodynamic and thermal structure of optically thick, dense parts of outflowing circumstellar disks that may be formed around various types of critically rotating massive stars, for example, Be stars, B[e] supergiant (sgB[e]) stars or Pop III stars. Methods. We specify the optical depth of the disk along the line-of-sight from stellar poles. Within the optically thick dense region we calculate the vertical disk thermal structure using the diffusion approximation while for the optically thin outer layers we assume a local thermodynamic equilibrium with the impinging stellar irradiation. We use two of our own types of hydrodynamic codes: two-dimensional operator-split numerical code and unsplit code based on the Roe's method. Results. Our models show the geometric distribution and contribution of viscous heating that begins to dominate in the central part of the disk. In the models of dense viscous disks the viscosity increases the central temperature up to several tens of thousands of Kelvins. The high mass-loss rates and high viscosity lead to instabilities with significant waves or bumps in density and temperature in the very inner disk region. Conclusions. The two-dimensional radial-vertical models of dense outflowing disks including the full Navier-Stokes viscosity terms show very high temperatures that are however limited to only the central disk cores inside the optically thick area, while near the edge of the optically thick region the temperature may be low enough for the existence of neutral hydrogen.

The accretion of the stellar wind material by a compact object represents the main mechanism powering the X-ray emission in classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. In this work we present the first attempt to simulate the accretion process of a fast and dense massive star wind onto a neutron star, taking into account the effects of the centrifugal and magnetic inhibition of accretion ("gating") due to the spin and magnetic field of the compact object. We made use of a radiative hydrodynamical code to model the non-stationary radiatively driven wind of an O-B supergiant star and then place a neutron star characterized by a fixed magnetic field and spin period at a certain distance from the massive companion. Our calculations follow, as a function of time (on a total time scale of several hours), the transition of the system through all different accretion regimes that are triggered by the intrinsic variations in the density and velocity of the non-stationary wind. The X-ray luminosity released by the system is computed at each time step by taking into account the relevant physical processes occurring in the different accretion regimes. Synthetic lightcurves are derived and qualitatively compared with those observed from classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. Although a number of simplifications are assumed in these calculations, we show that taking into account the effects of the centrifugal and magnetic inhibition of accretion significantly reduces the average X-ray luminosity expected for any neutron star wind-fed binary. The present model calculations suggest that long spin periods and stronger magnetic fields are favoured in order to reproduce the peculiar behavior of supergiant fast X-ray transients in the X-ray domain.

The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. Our group investigated the presence of intermediate-mass black holes for a sample of 10 Galactic globular clusters. We measured the inner kinematic profiles with integral-field spectroscopy and determined masses or upper limits of central black holes in each cluster. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M_bh - sigma) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). We ran N-body simulations of globular clusters containing IMBHs in a tidal field and studied their effects on mass-loss rates and remnant fractions and showed that an IMBH in the center prevents core collapse and ejects massive objects more rapidly. These simulations were further used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighbouring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole.

This article intends to provide a concise overview, from an observational point-of-view, of the current state of our knowledge of the most relevant properties of the Milky Way's nuclear star cluster (MWNSC). The MWNSC appears to be a typical specimen of nuclear star clusters, which are found at the centers of the majority of all types of galaxies. Nuclear clusters represent the densest and most massive stellar systems in the present-day Universe and frequently coexist with central massive black holes. They are therefore of prime interest for studying stellar dynamics and the MWNSC is the only one that allows us to obtain data on milli-parsec scales. After discussing the main observational constraints, we start with a description of the overall structure and kinematics of the MWNSC, then focus on a comparison to extragalactic systems, summarize the properties of the young, massive stars in the immediate environment of the Milky Way's central black hole, Sagittarius\,A*, and finally focus on the dynamics of stars orbiting the black hole at distances of a few to a few tens of milli parsecs.

Within the central 10pc of our Galaxy lies a dense nuclear star cluster (NSC), and similar NSCs are found in most nearby galaxies. Studying the structure and kinematics of NSCs reveals the history of mass accretion of galaxy nuclei. Because the Milky Way (MW) NSC is at a distance of only 8kpc, we can spatially resolve the MWNSC on sub-pc scales. This makes the MWNSC a reference object for understanding the formation of all NSCs. We have used the NIR long-slit spectrograph ISAAC (VLT) in a drift-scan to construct an integral-field spectroscopic map of the central 9.5 x 8pc of our Galaxy. We use this data set to extract stellar kinematics both of individual stars and from the unresolved integrated light spectrum. We present a velocity and dispersion map from the integrated light and model these kinematics using kinemetry and axisymmetric Jeans models. We also measure CO bandhead strengths of 1,375 spectra from individual stars. We find kinematic complexity in the NSCs radial velocity map including a misalignment of the kinematic position angle by 9 degree counterclockwise relative to the Galactic plane, and indications for a rotating substructure perpendicular to the Galactic plane at a radius of 20" or 0.8pc. We determine the mass of the NSC within r = 4.2pc to 1.4 x 10^7 Msun. We also show that our kinematic data results in a significant underestimation of the supermassive black hole (SMBH) mass. The kinematic substructure and position angle misalignment may hint at distinct accretion events. This indicates that the MWNSC grew at least partly by the mergers of massive star clusters. Compared to other NSCs, the MWNSC is on the compact side of the r_eff - M_NSC relation. The underestimation of the SMBH mass might be caused by the kinematic misalignment and a stellar population gradient. But it is also possible that there is a bias in SMBH mass measurements obtained with integrated light.

During their evolution massive stars can reach the phase of critical rotation when a further increase in rotational speed is no longer possible. Direct centrifugal ejection from a critically or near-critically rotating surface forms a gaseous equatorial decretion disk. Anomalous viscosity provides the efficient mechanism for transporting the angular momentum outwards. The outer part of the disk can extend up to a very large distance from the parent star. We study the evolution of density, radial and azimuthal velocity, and angular momentum loss rate of equatorial decretion disks out to very distant regions. We investigate how the physical characteristics of the disk depend on the distribution of temperature and viscosity. We calculated stationary models using the Newton-Raphson method. For time-dependent hydrodynamic modeling we developed the numerical code based on an explicit finite difference scheme on an Eulerian grid including full Navier-Stokes shear viscosity. The sonic point distance and the maximum angular momentum loss rate strongly depend on the temperature profile and are almost independent of viscosity. The rotational velocity at large radii rapidly drops accordingly to temperature and viscosity distribution. The total amount of disk mass and the disk angular momentum increase with decreasing temperature and viscosity. The time-dependent one-dimensional models basically confirm the results obtained in the stationary models as well as the assumptions of the analytical approximations. Including full Navier-Stokes viscosity we systematically avoid the rotational velocity sign change at large radii. The unphysical drop of the rotational velocity and angular momentum loss at large radii (present in some models) can be avoided in the models with decreasing temperature and viscosity.

(abridged) In this paper we aim at determining the shape, size, and luminosity/mass of the Milky Way Nuclear Star Cluster (MWNSC). We use Spitzer/IRAC images at 3.6 and 4.5 micrometer, where interstellar extinction is at a minimum but the overall emission is still dominated by stars. We correct the 4.5 micrometer image for PAH emission with the help of the IRAC 8.0 micrometer map and for extinction with the help of a [3.6-4.5] colour map. We present an extinction map for the central ~300x200 pc^2 of the Milky Way, as well as a PAH-emission and extinction corrected image of the stellar emission, with a resolution of about 0.2 pc. We find that the MWNSC appears in projection intrinsically point-symmetric, that it is significantly flattened, with its major axis aligned along the Galactic Plane, and that it is centred on the black hole, Sagittarius A*. Its density follows the well known approximate rho~r^{-2}-law at distances of a few parsecs from Sagittarius A*, but may become as steep as rho~r^{-3} at projected radii around 5 pc. We derive a half light radius of 4.2+-0.4 pc, a total luminosity of L_MWNSC=(4.1+-0.4)x10^{7} L_Sun, and a mass of M_{MWNSC}=(2.5+-0.4)x10^{7} M_Sun. The overall properties of the MWNSC agree well with the ones of its extragalactic counterparts, which underlines its role as a template for these objects. Its flattening agrees well with its previously established rotation parallel to Galactic rotation and suggests that it has formed by accretion of material that fell in preferentially along the Galactic Plane. Our findings support the in situ growth scenario for nuclear clusters and emphasize the need to increase the complexity of theoretical models for their formation and for the interaction between their stars and the central black hole in order to include rotation, axisymmetry, and growth in recurrent episodes.

Massive stars may during their evolution reach the phase of critical rotation when the further increase in rotational speed is no longer possible. The ejection of matter in the equatorial region forms the gaseous outflowing disk, which allows the star to remove the excess of its angular momentum. The outer part of the disk can extend up to large distance from the parent star. We study the evolution of density, radial and azimuthal velocity and angular momentum loss rate of equatorial decretion disk up to quite distant outer regions. We calculate the evolution of density, radial and azimuthal velocity from the initial Keplerian state until the disk approaches the final stationary state.

For galaxies hosting supermassive black holes (SMBHs), it has been observed that the mass of the central black hole (M_BH) tightly correlates with the effective or central velocity dispersion (sigma) of the host galaxy. The origin of this M_BH - sigma scaling relation is assumed to lie in the merging history of the galaxies but many open questions about its origin and the behavior in different mass ranges still need to be addressed. The goal of this work is to study the black-hole scaling relations for low black-hole masses, where the regime of intermediate-mass black holes (IMBHs) in globular clusters (GCs) is entered. We collect all existing reports of dynamical black-hole measurements in globular clusters, providing black-hole masses or upper limits for 14 candidates. We plot the black-hole masses versus different cluster parameters including total mass, velocity dispersion, concentration and half-mass radius. We search for trends and test the correlations in order to quantify their significance using a set of different statistical approaches. For correlations showing a large significance we perform a linear fit, accounting for uncertainties and upper limits. We find a clear correlation between the mass of the IMBH and the velocity dispersion of the globular cluster. As expected, the total mass of the globular cluster then also correlates with the mass of the IMBH. While the slope of the M_BH - sigma correlation differs strongly from the one observed for SMBHs, the other scaling relations M_BH - M_TOT, and M_BH - L are similar to the correlations in galaxies. Significant correlations of black-hole mass with other cluster properties were not found in the present sample.

Intermediate-mass black holes (IMBHs, 10^2-10^5 M_sun) fill the gap between stellar-mass black holes and supermassive black holes (SMBHs). Simulations have shown that IMBHs may form in dense star clusters, and therefore may still be present in these smaller stellar systems. We investigate the Galactic globular cluster NGC 5286 for indications of a central IMBH using spectroscopic data from VLT/FLAMES, velocity measurements from the Rutgers Fabry Perot at CTIO, and photometric data from HST. We run analytic spherical and axisymmetric Jeans models with different central black-hole masses, anisotropy, mass-to-light ratio, and inclination. Further, we compare the data to a grid of N-body simulations without tidal field. Additionally, we use one N-body simulation to check the results of the spherical Jeans models for the total cluster mass. Both the Jeans models and the N-body simulations favor the presence of a central black hole in NGC 5286 and our detection is at the 1- to 1.5-sigma level. From the spherical Jeans models we obtain a best fit with black-hole mass M_BH=(1.5+-1.0)x10^3 M_sun. The error is the 68% confidence limit from Monte Carlo simulations. Axisymmetric models give a consistent result. The best fitting N-body model is found with a black hole of 0.9% of the total cluster mass (4.38+-0.18)x10^5 M_sun, which results in an IMBH mass of M_BH=(3.9+-2.0)x10^3 M_sun. Jeans models give lower values for the total cluster mass. Our test of the Jeans models with N-body simulation data shows that this discrepancy has two reasons: The influence of a radially varying M/L profile, and underestimation of the velocity dispersion as the measurements are limited to bright stars. We conclude that detection of IMBHs in Galactic globular clusters remains a challenging task unless their mass fractions are above those found for SMBHs in nearby galaxies. [abridged]

For a sample of nine Galactic globular clusters we measured the inner kinematic profiles with integral-field spectroscopy that we combined with existing outer kinematic measurements and HST luminosity profiles. With this information we are able to detect the crucial rise in the velocity-dispersion profile which indicates the presence of a central black hole. In addition, N-body simulations compared to our data will give us a deeper insight in the properties of clusters with black holes and stronger selection criteria for further studies. For the first time, we obtain a homogeneous sample of globular cluster integral- field spectroscopy which allows a direct comparison between clusters with and without an intermediate-mass black hole.

The formation of supermassive black holes at high redshift still remains a puzzle to astronomers. Their growth becomes reasonable only when starting from a massive seed black hole with mass of the order of 10^2 - 10^5 M_SUN. Intermediate-mass black holes (IMBHs) are therefore an important field of research. Especially the possibility of finding them in the centers of globular clusters has recently drawn attention. The search for IMBHs in the centers of globular clusters could therefore shed light on the process of black-hole formation and cluster evolution. We are investigating six galactic globular clusters for the presence of an IMBH at their centers. Based on their kinematic and photometric properties, we selected the globular clusters NGC 1851, NGC 1904 (M79), NGC 5694, NGC 5824, NGC 6093 (M80) and NGC 6266 (M62). We use integral field spectroscopy in order to obtain the central velocity-dispersion profile of each cluster. We compute the cluster photometric center and the surface brightness profile using HST data. After combining these datasets we compare them to analytic Jeans models. We use varying M/L_V profiles for clusters with enough data points in order to reproduce their kinematic profiles in an optimal way. Finally, we vary the mass of the central black hole and test whether the cluster is better fitted with or without an IMBH. We present the statistical significance, including upper limits, of the black-hole mass for each cluster. NGC 1904 and NGC 6266 provide the highest significance for a black hole. Jeans models in combination with a M/L_V profile obtained from N-body simulations (in the case of NGC 6266) predict a central black hole of M_BH = (3 +- 1) x 10^3 M_SUN for NGC 1904 and M_BH = (2 +- 1) x 10^3 M_SUN for NGC 6266. Furthermore, we discuss the possible influence of dark remnants and mass segregation at the center of the cluster on the detection of an IMBH.

The clumping of massive star winds is an established paradigm, which is confirmed by multiple lines of evidence and is supported by stellar wind theory. We use the results from time-dependent hydrodynamical models of the instability in the line-driven wind of a massive supergiant star to derive the time-dependent accretion rate on to a compact object in the Bondi-Hoyle-Lyttleton approximation. The strong density and velocity fluctuations in the wind result in strong variability of the synthetic X-ray light curves. Photoionization of inhomogeneous winds is different from the photoinization of smooth winds. The degree of ionization is affected by the wind clumping. The wind clumping must also be taken into account when comparing the observed and model spectra of the photoionized stellar wind.

The true mass-loss rates from massive stars are important for many branches of astrophysics. For the correct modeling of the resonance lines, which are among the key diagnostics of stellar mass-loss, the stellar wind clumping turned out to be very important. In order to incorporate clumping into radiative transfer calculation, 3-D models are required. Various properties of the clumps may have strong impact on the resonance line formation and, therefore, on the determination of empirical mass-loss rates. We incorporate the 3-D nature of the stellar wind clumping into radiative transfer calculations and investigate how different model parameters influence the resonance line formation. We develop a full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. The number density of clumps follows the mass conservation. For the first time, realistic 3-D models that describe the dense as well as the tenuous wind components are used to model the formation of resonance lines in a clumped stellar wind. At the same time, non-monotonic velocity fields are accounted for. The 3-D density and velocity wind inhomogeneities show very strong impact on the resonance line formation. The different parameters describing the clumping and the velocity field results in different line strengths and profiles. We present a set of representative models for various sets of model parameters and investigate how the resonance lines are affected. Our 3-D models show that the line opacity is reduced for larger clump separation and for more shallow velocity gradients within the clumps. Our new model demonstrates that to obtain empirically correct mass-loss rates from the UV resonance lines, the wind clumping and its 3-D nature must be taken into account.

We study the effects of optically thick clumps, non-void inter-clump medium, variation of the onset of clumping, and velocity dispersion inside clumps on the formation of resonance lines. For this purpose we developed a full 3-D Monte Carlo Radiative Transfer (MCRT) code that is able to handle 3-D shapes of clumps and arbitrary 3-D velocity fields. The method we developed allows us to take into account contributions from density and velocity wind inhomogeneities to the total opacity very precisely. The first comparison with observation shows that 3-D density and velocity wind inhomogeneities have a very strong influence on the resonance line formation, and that they have to be accounted for in order to obtain reliable mass-loss rate determinations.

The clumping of massive star winds is an established paradigm confirmed by multiple lines of evidence and supported by stellar wind theory. The purpose of this paper is to bridge the gap between detailed models of inhomogeneous stellar winds in single stars and the phenomenological description of donor winds in supergiant high-mass X-ray binaries (HMXBs). We use results from time-dependent hydrodynamical models of the instability in the line-driven wind of a massive supergiant star to derive the time-dependent accretion rate onto a compact object in the Bondi-Hoyle-Lyttleton approximation. The strong density and velocity fluctuations in the wind result in strong variability of the synthetic X-ray light curves. The model predicts a large scale X-ray variability, up to eight orders of magnitude, on relatively short timescales. The apparent lack of evidence for such strong variability in the observed HMXBs indicates that the details of accretion process act to reduce the variability due to the stellar wind velocity and density jumps. We, also, study the absorption of X-rays in the clumped stellar wind by means of a 2-D stochastic wind model and find that absorption of X-rays changes strongly at different orbital phases. Furthermore, we address the photoionization in the clumped wind, and show that the degree of ionization is affected by the wind clumping. A correction factor for the photoionization parameter is derived. It is shown that the photoionization parameter is reduced by a factor Xi compared to the smooth wind models with the same mass-loss rate, where Xi is the wind inhomogeneity parameter. We conclude that wind clumping must also be taken into account when comparing the observed and model spectra of the photoionized stellar 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)

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

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.

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

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.

Hot stars are sources of X-ray emission originating in their winds. Although hydrodynamical simulations that are able to predict this X-ray emission are available, the inclusion of X-rays in stationary wind models is usually based on simplifying approximations. To improve this, we use results from time-dependent hydrodynamical simulations of the line-driven wind instability (seeded by the base perturbation) to derive the analytical approximation of X-ray emission in the stellar wind. We use this approximation in our non-LTE wind models and find that an improved inclusion of X-rays leads to a better agreement between model ionization fractions and those derived from servations. Furthermore, the slope of the L_x-L relation is in better agreement with observations, however the X-ray luminosity is underestimated by a factor of three. We propose a possible solution for this discrepancy.

Analyzing the spectra from Wolf-Rayet stars requires adequate non-LTE modeling of their expanding atmosphere. The numerical schemes for solving the radiative transfer in the co-moving frame of reference have been developed by Mihalas and co-workers 30 years ago. The most elaborate codes can cope today with many hundred explicit non-LTE levels or super-levels and account for metal-line blanketing. The limited agreement with observed spectra indicates that the model simplifications are still severe. One approximation that has to be blamed is homogeneity. Stellar-wind clumping on small scales was easily implemented, while "macro-clumping" is still a big challenge. First studies showed that macro-clumping can reduce the strength of predicted P-Cygni line profiles in O-star spectra, and largely affects the X-ray line spectra from stellar winds. The classical model for radiation-driven winds by Castor, Abbot and Klein fails to explain the very dense winds from Wolf-Rayet stars. Only when we solved the detailed non-LTE radiative transfer consistently with the hydrodynamic equations, mass-loss rates above the single-scattering limit have been obtained.

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

X-ray spectroscopy is a sensitive probe of stellar winds. X-rays originate from optically thin shock-heated plasma deep inside the wind and propagate outwards throughout absorbing cool material. Recent analyses of the line ratios from He-like ions in the X-ray spectra of O-stars highlighted problems with this general paradigm: the measured line ratios of highest ions are consistent with the location of the hottest X-ray emitting plasma very close to the base of the wind, perhaps indicating the presence of a corona, while measurements from lower ions conform with the wind-embedded shock model. Generally, to correctly model the emerging X-ray spectra, a detailed knowledge of the cool wind opacities based on stellar atmosphere models is prerequisite. A nearly grey stellar wind opacity for the X-rays is deduced from the analyses of high-resolution X-ray spectra. This indicates that the stellar winds are strongly clumped. Furthermore, the nearly symmetric shape of X-ray emission line profiles can be explained if the wind clumps are radially compressed. In massive binaries the orbital variations of X-ray emission allow to probe the opacity of the stellar wind; results support the picture of strong wind clumping. In high-mass X-ray binaries, the stochastic X-ray variability and the extend of the stellar-wind part photoionized by X-rays provide further strong evidence that stellar winds consist of dense clumps.

Context: The mass-loss rate is a key parameter of massive stars. Adequate stellar atmosphere models are required for spectral analyses and mass-loss determinations. Present models can only account for the inhomogeneity of stellar winds in the approximation of small-scale structures that are optically thin. This treatment of ``microclumping'' has led to reducing empirical mass-loss rates by factors of two and more. Aims: Stellar wind clumps can be optically thick in spectral lines. We investigate how this ``macroclumping'' impacts on empirical mass-loss rates. Methods: The Potsdam Wolf-Rayet (PoWR) model atmosphere code is generalized in the ``formal integral'' to account for clumps that are not necessarily optically thin. Results: Optically thick clumps reduce the effective opacity. This has a pronounced effect on the emergent spectrum. Our modeling for the O-type supergiant zeta Puppis reveals that the optically thin H-alpha line is not affected by wind porosity, but that the PV resonance doublet becomes significantly weaker when macroclumping is taken into account. The reported discrepancies between resonance-line and recombination-line diagnostics can be resolved entirely with the macroclumping modeling without downward revision of the mass-loss rate. Conclusions: Mass-loss rates inferred from optically thin emission, such as the H-alpha line in O stars, are not influenced by macroclumping. The strength of optically thick lines, however, is reduced because of the porosity effects. Therefore, neglecting the porosity in stellar wind modeling can lead to underestimating empirical mass-loss rates.

We have developed a time-dependent two-component hydrodynamics code to simulate radiatively-driven stellar winds from hot stars. We use a time-explicit van Leer scheme to solve the hydrodynamic equations of a two-component stellar wind. Dynamical friction due to Coulomb collisions between the passive bulk plasma and the line-scattering ions is treated by a time-implicit, semi-analytic method using a polynomial fit to the Chandrasekhar function. This gives stable results despite the stiffness of the problem. This method was applied to model stars with winds that are both poorly and well-coupled. While for the former case we reproduce the mCAK solution, for the latter case our solution leads to wind decoupling.

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.

Archival X-ray spectra of the four prominent single, non-magnetic O stars Zeta Pup, Zeta Ori, Ksi Per and Zeta Oph, obtained in high resolution with Chandra HETGS/MEG have been studied. The resolved X-ray emission line profiles provide information about the shocked, hot gas which emits the X-radiation, and about the bulk of comparably cool stellar wind material which partly absorbs this radiation. In this paper, we synthesize X-ray line profiles with a model of a clumpy stellar wind. We find that the geometrical shape of the wind inhomogeneities is important: better agreement with the observations can be achieved with radially compressed clumps than with spherical clumps. The parameters of the model, i.e. chemical abundances, stellar radius, mass-loss rate and terminal wind velocity, are taken from existing analyses of UV and optical spectra of the programme stars. On this basis, we also calculate the continuum-absorption coefficient of the cool-wind material, using the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The radial location of X-ray emitting gas is restricted from analysing the fir line ratios of helium-like ions. The only remaining free parameter of our model is the typical distance between the clumps; here, we assume that at any point in the wind there is one clump passing by per one dynamical time-scale of the wind. The total emission in a model line is scaled to the observation. There is a good agreement between synthetic and observed line profiles. We conclude that the X-ray emission line profiles in O stars can be explained by hot plasma embedded in a cool wind which is highly clumped in the form of radially compressed shell fragments.

Radiative transfer in a clumped winds is used to describe X-ray emission line profiles observed in the XMM-Newton RGS spectrum of the OI star Zeta Puppis. It is shown that this X-ray spectrum can be explained as originating from a multi-temperature collisional plasma located in the wind acceleration zone. The X-rays are attenuated in the clumped stellar wind, which gives characteristic profiles to the emergent lines. We specifically study the N VII emission line in the spectrum of Zeta Pup. Long RGS exposures reveal that the N VII line profile is structured. On the basis of our Zeta Pup atmosphere model, we rule out the presence of N VII in the cool wind component. We suggest that the detailed N VII line structure is due to self-absorption in the hot plasma. Wind clumping also affects the transfer of ionizing radiation in high-mass X-ray binaries (HMXBs). We derive analytical formulae for the ionisation parameter in dependence on the parameters of wind clumping.

It is commonly adopted that X-rays from O stars are produced deep inside the stellar wind, and transported outwards through the bulk of the expanding matter which attenuates the radiation and affects the shape of emission line profiles. The ability of Chandra and XMM-Newton to resolve these lines spectroscopically provided a stringent test for the theory of X-ray production. It turned out that none of the existing models was able to reproduce the observations consistently. The major caveat of these models was the underlying assumption of a smooth stellar wind. Motivated by the various observational evidence that the stellar winds are in fact structured, we present a 2-D model of a stochastic, inhomogeneous wind. The X-ray radiative transfer is derived for such media. It is shown that profiles from a clumped wind differ drastically from those predicted by conventional homogeneous models. We review the up-to-date observations of X-ray line profiles from stellar winds and present line fits obtained from the inhomogeneous wind model. The necessity to account for inhomogeneities in calculating the X-ray transport in massive star winds, including for HMXB is highlighted.

There is strong observational evidence of shocks and clumping in radiation-driven stellar winds from hot, luminous stars. The resulting non monotonic velocity law allows for radiative coupling between distant locations, which is so far not accounted for in hydrodynamic wind simulations. In the present paper, we determine the Sobolev source function and radiative line force in the presence of radiative coupling in spherically symmetric flows, extending the geometry-free formalism of Rybicki and Hummer (1978) to the case of three-point coupling, which can result from, e.g., corotating interaction regions, wind shocks, or mass overloading. For a simple model of an overloaded wind, we find that, surprisingly, the flow decelerates at all radii above a certain height when nonlocal radiative coupling is accounted for. We discuss whether radiation-driven winds might in general not be able to re-accelerate after a non monotonicity has occurred in the velocity law.

It is commonly adopted that X-rays from O stars are produced deep inside the stellar wind, and transported outwards through the bulk of the expanding matter which attenuates the radiation and affects the shape of emission line profiles. None of the existing models was able to fit the observationed lines consistently. The possible caveat of these models was the assumption of a smooth stellar wind. Motivated by the various evidences that the stellar winds are in fact structured, we present a 2-D numerical model of a stochastic, inhomogeneous wind. Small parcels of hot, X-ray emitting gas are permeated with cool, absorbing wind material which is compressed into thin shell fragments. Wind fragmentation alters the radiative transfer drastically, compared to homogeneous models of the same mass-loss rate. X-rays produced deep inside the wind, which would be totally absorbed in a homogeneous flow, can effectively escape from a fragmented wind. The wind absorption becomes wavelength independent if the individual fragments are optically thick. The X-ray line profiles are flat-topped in the blue part and decline steeply in the red part for the winds with short acceleration zone. For the winds where the acceleration extends over significant distances, the lines can appear nearly symmetric and only slightly blueshifted, in contrast to the skewed, triangular line profiles typically obtained from homogeneous wind models of high optical depth. We show that profiles from a fragmented wind model can reproduce the observed line profiles from zeta Orionis.

We discuss X-ray line formation in dense O star winds. A random distribution of wind shocks is assumed to emit X-rays that are partially absorbed by cooler wind gas. The cool gas resides in highly compressed fragments oriented perpendicular to the radial flow direction. For fully opaque fragments, we find that the blueshifted part of X-ray line profiles remains flat-topped even after severe wind attenuation, whereas the red part shows a steep decline. These box-type, blueshifted profiles resemble recent Chandra observations of the O3 star zeta Pup. For partially transparent fragments, the emission lines become similar to those from a homogeneous wind.

Line-driven winds from stars and accretion disks are accelerated by scattering in numerous line transitions. The wind is believed to adopt a unique critical solution, out of the infinite variety of shallow and steep solutions. We study the inherent dynamics of the transition towards the critical wind. A new runaway wind mechanism is analyzed in terms of radiative-acoustic (Abbott) waves which are responsible for shaping the wind velocity law and fixing the mass loss. Three different flow types result, depending on the location of perturbations. First, if the shallow solution is perturbed sufficiently far downstream, a single critical point forms in the flow, which is a barrier for Abbott waves, and the solution tends to the critical one. Second, if the shallow solution is perturbed upstream from this critical point, mass overloading results, and the critical point is shifted inwards. This wind exhibits a broad, stationary region of decelerating flow and its velocity law has kinks. Third, for perturbations even further upstream, the overloaded wind becomes time-dependent, and develops shocks and dense shells.

Observations and theory suggest that line driven winds from hot stars and luminous accretion disks adopt a unique, critical solution which corresponds to maximum mass loss rate. We analyze the numerical stability of the infinite family of shallow wind solutions, which resemble solar wind breezes, and their transition to the critical wind. Shallow solutions are sub-critical with respect to radiative (or Abbott) waves. These waves can propagate upstream through shallow winds at high speeds. If the waves are not accounted for in the Courant time step, numerical runaway results. The outer boundary condition is equally important for wind stability. Assuming pure outflow conditions, as is done in the literature, triggers runaway of shallow winds to the critical solution or to accretion flow.

Line-driven winds from hot stars and accretion disks are thought to adopt a unique, critical solution which corresponds to maximum mass loss rate and a particular velocity law. We show that in the presence of negative velocity gradients, radiative-acoustic (Abbott) waves can drive shallow wind solutions towards larger velocities and mass loss rates. Perturbations introduced downstream from the wind critical point lead to convergence towards the critical solution. By contrast, low-lying perturbations cause evolution towards a mass-overloaded solution, developing a broad deceleration region in the wind. Such a wind differs fundamentally from the critical solution. For sufficiently deep-seated perturbations, overloaded solutions become time-dependent and develop shocks and shells.

We analyze the dynamics of 2D stationary line-driven winds from accretion disks in cataclysmic variables (CVs), by generalizing the Castor, Abbott and Klein theory. In paper 1, we have solved the wind Euler equation, derived its two eigenvalues, and addressed the solution topology and wind geometry. Here, we focus on mass loss and velocity laws. We find that disk winds, even in luminous novalike variables, have low optical depth, even in the strongest driving lines. This suggests that thick-to-thin transitions in these lines occur. For disks with a realistic radial temperature, the mass loss is dominated by gas emanating from the inner decade in r. The total mass loss rate associated with a luminosity 10 Lsun is 10^{-12} Msun/yr, or 10^{-4} of the mass accretion rate. This is one order of magnitude below the lower limit obtained from P Cygni lines, when the ionizing flux shortwards of the Lyman edge is supressed. The difficulties with such small mass loss rates in CVs are principal, and confirm our previous work. We conjecture that this issue may be resolved by detailed nonLTE calculations of the line force within the context of CV disk winds, and/or better accounting for the disk energy distribution and wind ionization structure. We find that the wind velocity profile is well approximated by the empirical law used in kinematical modeling. The acceleration length scale is given by the footpoint radius of the wind streamline in the disk. This suggests an upper limit of 10 Rwd to the acceleration scale, which is smaller by factors of a few as compared to values derived from line fitting.

We analyze the dynamics of 2-D stationary, line-driven winds from accretion disks in cataclysmic variable stars. The driving force is that of line radiation pressure, in the formalism developed by Castor, Abbott & Klein for O stars. Our main assumption is that wind helical streamlines lie on straight cones. We find that the Euler equation for the disk wind has two eigenvalues, the mass loss rate and the flow tilt angle with the disk. Both are calculated self-consistently. The wind is characterized by two distinct regions, an outer wind launched beyond four white dwarf radii from the rotation axis, and an inner wind launched within this radius. The inner wind is very steep, up to 80 degrees with the disk plane, while the outer wind has a typical tilt of 60 degrees. In both cases the ray dispersion is small. We, therefore, confirm the bi-conical geometry of disk winds as suggested by observations and kinematical modeling. The wind collimation angle appears to be robust and depends only on the disk temperature stratification. The flow critical points lie high above the disk for the inner wind, but close to the disk photosphere for the outer wind. Comparison with existing kinematical and dynamical models is provided. Mass loss rates from the disk as well as wind velocity laws are discussed in a subsequent paper.