A cyclic bipolar wind in the interacting binary V 393 Scorpii

V393 Scorpii is a double periodic variable characterized by a relatively stable non-orbital photometric cycle of 253 d. Mennickent et al. argue for the presence of a massive optically thick disc around the more massive B-type component and describe the evolutionary stage of the system. In this paper, we analyse the behaviour of the main spectroscopic optical lines during the long non-orbital photometric cycle. We study the radial velocity of the donor determining its orbital elements and find a small but significant orbital eccentricity (e = 0.04). The donor spectral features are modelled and removed from the spectrum at every observing epoch using the light-curve model given by Mennickent et al. We find that the line emission is larger during eclipses and mostly comes from a bipolar wind. We also find that the long cycle is explained in terms of a modulation of the wind strength; the wind has a larger line and continuum emissivity at the high state. We report the discovery of highly variable chromospheric emission in the donor, as revealed by the Doppler maps of the emission lines Mg II 4481 and C I 6588. We discuss notable and some novel spectroscopic features like discrete absorption components, especially visible at blue depressed O I 7773 absorption wings during the second half-cycle, Balmer double emission with V/R curves showing 'Z-type' and 'S-type' excursions around secondary and main eclipses, respectively, and H beta emission wings extending up to +/- 2000 km s(-1). We also discuss possible causes for these phenomena and for their modulations with the long cycle.

The relationship between gamma Cassiopeiae's X-ray emission and its circumstellar environment

gamma Cas is the prototypical classical Be star and is recently best known for its variable hard X-ray emission. To elucidate the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around four XMM-Newton observations. The observational techniques included long baseline optical interferometry (LBOI) from two instruments at CHARA, photometry carried out by an automated photometric telescope and H alpha observations. Because gamma Cas is also known to be in a binary, we measured radial velocities from the H alpha line and redetermined its period as 203.55 +/- 0.20 days and its eccentricity as near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an system inclination angle near 45 degrees, and a larger radius than previously reported. In addition, the Be star began an "outburst" at the beginning of our campaign, made visible by a brightening and reddening of the disk during our campaign and beyond. Our analyses of the new high resolution spectra disclosed many attributes also found from spectra obtained in 2001 (Chandra) and 2004 (XMM-Newton). As well as a dominant hot (approximate to 14 keV) thermal component, the familiar attributes included: (i) a fluorescent feature of Fe K even stronger than observed at previous times; (ii) strong lines of N VII and Ne XI lines indicative of overabundances; and (iii) a subsolar Fe abundance from K-shell lines but a solar abundance from L-shell ions. We also found that two absorption columns are required to fit the continuum. While the first one maintained its historical average of 1 x 10(21) cm(-2), the second was very large and doubled to 7.4 x 10(23) cm(-2) during our X-ray observations. Although we found no clear relation between this column density and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind matter ejected by gamma Cas into our line of sight.

On the nature of candidate luminous blue variables in M 33

Context. Luminous blue variables (LBVs) are a class of highly unstable stars that have been proposed to play a critical role in massive stellar evolution as well as being the progenitors of some of the most luminous supernovae known. However the physical processes underlying their characteristic instabilities are currently unknown. Aims. In order to provide observational constraints on this behaviour we have initiated a pilot study of the population of (candidate) LBVs in the Local Group galaxy M 33. Methods. To accomplish this we have obtained new spectra of 18 examples within M 33. These provide a baseline of ≥ 4 yr with respect to previous observations, which is well suited to identifying LBV outbursts. We also employed existing multi-epoch optical and mid-IR surveys of M 33 to further constrain the variability of the sample and search for the presence of dusty ejecta. Results. Combining the datasets reveals that spectroscopic and photometric variability appears common, although in the majority of cases further observations will be needed to distinguish between an origin for this behavour in short lived stochastic wind structure and low level photospheric pulsations or coherent long term LBV excursions. Of the known LBVs we report a hitherto unidentified excursion of M 33 Var C between 2001-5, while the transition of the WNLh star B517 to a cooler B supergiant phase between 1993−2010 implies an LBV classification. Proof-of-concept quantitative model atmosphere analysis is provided for Romano’s star; the resultant stellar parameters being consistent with the finding that the LBV excursions of this star are accompanied by changes in bolometric luminosity. The combination of temperature and luminosity of two stars, the BHG [HS80] 110A and the cool hypergiant B324, appear to be in violation of the empirical Humphreys-Davidson limit. Mid-IR observations demonstrate that a number of candidates appear associated with hot circumstellar dust, although no objects as extreme as η Car are identified. The combined dataset suggests that the criteria employed to identify candidate LBVs results in a heterogeneous sample, also containing stars demonstrating the B[e] phenomenon. Of these, a subset of optically faint, low luminosity stars associated with hot dust are of particular interest since they appear similar to the likely progenitor of SN 2008S and the 2008 NGC 300 transient (albeit suffering less intrinsic extinction). Conclusions. The results of such a multiwavelength observational approach, employing multiplexing spectrographs and supplemented with quantitative model atmosphere analysis, appears to show considerable promise in both identifying and characterising the physical properties of LBVs as well as other short lived phases of massive stellar evolution.

An XMM-Newton view of FeKα in high-mass X-ray binaries

We present a comprehensive analysis of the whole sample of available XMM-Newton observations of high-mass X-ray binaries (HMXBs) until August 2013, focusing on the FeKα emission line. This line is key to better understanding the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We collected observations from 46 HMXBs and detected FeKα in 21 of them. We used the standard classification of HMXBs to divide the sample into different groups. We find that (1) different classes of HMXBs display different qualitative behaviours in the FeKα spectral region. This is visible especially in SGXBs (showing ubiquitous Fe fluorescence but not recombination Fe lines) and in γ Cass analogues (showing both fluorescent and recombination Fe lines). (2) FeKα is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states that are lower than Fe xviii. (3) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (4) We observe an inverse correlation between the X-ray luminosity and the equivalent width of FeKα (EW). This phenomenon is known as the X-ray Baldwin effect. (5) FeKα is narrow (σline< 0.15 keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering, and Doppler shifts (with velocities of the reprocessing material V ~ 1000 km s-1). (6) The equivalent hydrogen column (NH) directly correlates to the EW of FeKα, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. (7) The observed NH in supergiant X-ray binaries (SGXBs) is in general higher than in supergiant fast X-ray transients (SFXTs). We suggest two possible explanations: different orbital configurations or a different interaction compact object – wind. (8) Finally, we analysed the sources IGR J16320-4751 and 4U 1700-37 in more detail, covering several orbital phases. The observed variation in NH between phases is compatible with the absorption produced by the wind of their optical companions. The results clearly point to a very important contribution of the donor’s wind in the FeKα emission and the absorption when the donor is a supergiant massive star.

The donor star of the X-ray pulsar X1908+075

High-mass X-ray binaries consist of a massive donor star and a compact object. While several of those systems have been well studied in X-rays, little is known for most of the donor stars as they are often heavily obscured in the optical and ultraviolet regime. There is an opportunity to observe them at infrared wavelengths, however. The goal of this study is to obtain the stellar and wind parameters of the donor star in the X1908+075 high-mass X-ray binary system with a stellar atmosphere model to check whether previous studies from X-ray observations and spectral morphology lead to a sufficient description of the donor star. We obtained H- and K-band spectra of X1908+075 and analysed them with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. For the first time, we calculated a stellar atmosphere model for the donor star, whose main parameters are: Mspec = 15 ± 6 M⊙, T∗ = 23-3+6 kK, log geff = 3.0 ± 0.2 and log L/L⊙ = 4.81 ± 0.25. The obtained parameters point towards an early B-type (B0–B3) star, probably in a supergiant phase. Moreover we determined a more accurate distance to the system of 4.85 ± 0.50 kpc than the previously reported value.

On the Radial Onset of Clumping in the Wind of the B0I Massive Star QV Nor

We present an analysis of a 78 ks Chandra high-energy transmission gratings observation of the B0I star QV Nor, the massive donor of the wind-accreting pulsar 4U1538−52. The neutron star (NS) orbits its companion in a very close orbit (r < 1.4R*, in units of the stellar radii), thereby allowing probing of the innermost wind regions. The flux of the Fe Kα line during eclipse reduces to only ∼30% of the flux measured out of eclipse. This indicates that the majority of Fe fluorescence must be produced in regions close to the NS, at distances smaller than 1R* from its surface. The fact that the flux of the continuum decreases to only ∼3% during eclipse allows for a high contrast of the Fe Kα line fluorescence during eclipse. The line is not resolved and centered at 1.9368 0.0018 l = 0.0032 - + Å. From the inferred plasma speed limit of v < c l < 800 l D km s−1 and range of ionization parameters of log 1, 2 x = [- ], together with the stellar density profile, we constrain the location of the cold, dense material in the stellar wind of QV Nor using simple geometrical considerations. We then use the Fe Kα line fluorescence as a tracer of wind clumps and determine that these clumps in the stellar wind of QV Nor (B0I) must already be present at radii r < 1.25R*, close to the photosphere of the star.

Measuring the stellar wind parameters in IGR J17544-2619 and Vela X-1 constrains the accretion physics in supergiant fast X-ray transient and classical supergiant X-ray binaries

Context. Classical supergiant X-ray binaries (SGXBs) and supergiant fast X-ray transients (SFXTs) are two types of high-mass X-ray binaries (HMXBs) that present similar donors but, at the same time, show very different behavior in the X-rays. The reason for this dichotomy of wind-fed HMXBs is still a matter of debate. Among the several explanations that have been proposed, some of them invoke specific stellar wind properties of the donor stars. Only dedicated empiric analysis of the donors’ stellar wind can provide the required information to accomplish an adequate test of these theories. However, such analyses are scarce. Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J17544-2619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 ± 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e< 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (ν∞ = 1500 km s-1 in IGR J17544-2619 and ν∞ = 700 km s-1 in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.

Chemical abundances and winds of massive stars in M31: a B-type supergiant and a WC star in OB 10

We present high quality spectroscopic data for two massive stars in the OB 10 association of M31, OB 10-64 (B0 la) and OB 10-WRI (WC6). Medium resolution spectra of both stars were obtained using the ISIS spectrograph on the William Herschel Telescope. This is supplemented with Hubble Space Telescope STIS UV spectroscopy and Keck I HIRES data for OB 10-64. A non- local thermodynamic equilibrium (LTE) model atmosphere and abundance analysis for OB 10-64 is presented, indicating that this star has similar photospheric CNO, Mg and Si abundances to solar neighbourhood massive stars. A wind analysis of this early B-type supergiant reveals a mass-loss rate of (M)over dot = 1.6 x 10(-6) M-circle dot yr(-1), and v(infinity) = 1650 km s(-1). The corresponding wind momentum is in good agreement with the wind momentum-luminosity relationship found for Galactic early-B supergiants. Observations of OB 10-WRI are analysed using a non-LTE, line-blanketed code, to reveal approximate stellar parameters of log L/L-circle dot similar to 5.7, T-* - 75 kK, v(infinity) similar to 3000 km s(-1), (M)over dot/(M-circle dot yr(-1)) similar to 10(-4.3) adopting a clumped wind with a filling factor of 10 per cent. Quantitative comparisons are made with the Galactic WC6 star HD 92809 (WR23) revealing that OB 10-WR1 is 0.4 dex more luminous, though it has a much lower C/He ratio (similar to0.1 versus 0.3 for HD 92809). Our study represents the first detailed, chemical model atmosphere analysis for either a B-type supergiant or a Wolf- Rayet (WR) star in Andromeda, and shows the potential of how such studies can provide new information on the chemical evolution of galaxies and the evolution of massive stars in the local Universe.

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

The VLT-FLAMES survey of massive stars: stellar parameters and rotational velocities in NGC 3293, NGC 4755 and NGC 6611

An analysis is presented of VLT-FLAMES spectroscopy for three Galactic clusters, NGC3293, NGC4755 and NGC6611. Non-LTE model atmosphere calculations have been used to estimate effective temperatures (from either the helium spectrum or the silicon ionization equilibrium) and gravities (from the hydrogen spectrum). Projected rotational velocities have been deduced from the helium spectrum (for fast and moderate rotators) or the metal line spectrum (for slow rotators). The origin of the low gravity estimates for apparently near main sequence objects is discussed and is related to the stellar rotational velocity. The atmospheric parameters have been used to estimate cluster distances (which are generally in good agreement with previous determinations) and these have been used to estimate stellar luminosities and evolutionary masses. The observed Hertzsprung-Russell diagrams are compared with theoretical predictions and some discrepancies including differences in the main sequence luminosities are discussed. Cluster ages have been deduced and evidence for non-coeval star formation is found for all three of the clusters. Projected rotational velocities for targets in the older clusters, NGC3293 and NGC4755, have been found to be systematically larger than those for the field, confirming recent results in other similar age clusters. The distribution of projected rotational velocities are consistent with a Gaussian distribution of intrinsic rotational velocities. For the relatively unevolved targets in the older clusters, NGC3293 and NGC4755, the peak of the velocity distribution would be 250 km s(-1) with a full-width-half-maximum of approximately 180 km s(-1). For NGC6611, the sample size is relatively small but implies a lower mean rotational velocity. This may be evidence for the spin-down effect due to angular momentum loss through stellar winds, although our results are consistent with those found for very young high mass stars. For all three clusters we deduce present day mass functions with Gamma-values in the range of -1.5 to -1.8, which are similar to other young stellar clusters in the Milky Way.

The VLT-FLAMES survey of massive stars: mass loss and rotation of early-type stars in the SMC

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method (Mokiem et al. 2005, A&A, 441, 711), which combines the stellar atmosphere code FASTWIND (Puls et al. 2005, A&A, 435, 669) with the genetic algorithm based optimisation routine PIKAIA (Charbonneau 1995, ApJS, 101, 309). Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The automated method allows for a detailed determination of the projected rotational velocities. The present day v(r) sin i distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v(r)) distribution that can be characterised by a mean velocity of about 160-190 km s(-1) and an effective half width of 100-150 km s(-1). The vr distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct (Maeder & Meynet 2001, A&A, 373, 555), the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring vr. The mass loss rates of the SMC objects having luminosities of log L-star/L-circle dot greater than or similar to 5.4 are in excellent agreement with predictions by Vink et al. (2001, A&A, 369, 574). However, for lower luminosity stars the winds are too weak to determine. M accurately from the optical spectrum. Three targets were classified as Vz stars, two of which are located close to the theoretical zero-age main sequence. Three lower luminosity targets that were not classified as Vz stars are also found to lie near the ZAMS. We argue that this is related to a temperature effect inhibiting cooler from displaying the spectral features required for the Vz luminosity class.

The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational velocity distributions for B-type stars in the Magellanic Clouds

Aims.We aim to provide the atmospheric parameters and rotational velocities for a large sample of O- and early B-type stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. Methods: Atmospheric parameters, stellar masses, and rotational velocities have been estimated for approximately 250 early B-type stars in the Large (LMC) and Small (SMC) Magellanic Clouds from high-resolution VLT-FLAMES data using the non-LTE TLUSTY model atmosphere code. This data set has been supplemented with our previous analyses of some 50 O-type stars (Mokiem et al. 2006, 2007) and 100 narrow-lined early B-type stars (Hunter et al. 2006; Trundle et al. 2007) from the same survey, providing a sample of ~400 early-type objects. Results: Comparison of the rotational velocities with evolutionary tracks suggests that the end of core hydrogen burning occurs later than currently predicted and we argue for an extension of the evolutionary tracks. We also show that the large number of the luminous blue supergiants observed in the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational velocities nor their position on the H-R diagram are predicted. We suggest that blue loops or mass-transfer binary systems may populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to a similar Galactic sample, we find that (at 3s confidence level) massive stars (above 8 M?) in the SMC rotate faster than those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity distributions in the Galaxy and the LMC. We find that the v sin i distributions in the SMC and LMC can modelled with an intrinsic rotational velocity distribution that is a Gaussian peaking at 175 km s-1 (SMC) and 100 km s-1 (LMC) with a 1/e half width of 150 km s-1. We find that in NGC 346 in the SMC, the 10-25 M? main-sequence stars appear to rotate faster than their higher mass counterparts. It is not expected that O-type stars spin down significantly through angular momentum loss via stellar winds at SMC metallicity, hence this could be a reflection of mass dependent birth spin rates. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star progenitors. Our measured rotational velocity distribution for the 10-25 M? stars is peaked at slightly higher velocities than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14 M? at low metallicities.

The empirical metallicity dependence of the mass-loss rate of O- and early B-type stars

We present a comprehensive study of the observational dependence of the mass-loss rate in stationary stellar winds of hot massive stars on the metal content of their atmospheres. The metal content of stars in the Magellanic Clouds is discussed, and a critical assessment is given of state-of-the-art mass-loss determinations of OB stars in these two satellite systems and the Milky-Way. Assuming a power-law dependence of mass loss on metal content,. M. Z(m), and adopting a theoretical relation between the terminal flow velocity and metal content, v(infinity). Z(0.13) (Leitherer et al. 1992, ApJ, 401, 596), we find m = 0.83 +/- 0.16 for non-clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than 105.2 L circle dot. Within the errors, this result is in agreement with the prediction m = 0.69 +/- 0.10 by Vink et al. (2001, A& A, 369, 574). Absolute empirical values for the mass loss, based on Ha and ultraviolet (UV) wind lines, are found to be a factor of two higher than predictions in this high luminosity regime. If this difference is attributed to inhomogeneities in the wind, and this clumping does not impact the predictions, this would imply that luminous O and early-B stars have clumping factors in their Ha and UV line forming regions of about a factor of four. For lower luminosity stars, the winds are so weak that their strengths can generally no longer be derived from optical spectral lines (essentially Ha) and one must currently rely on the analysis of UV lines. We confirm that in this low-luminosity domain the observed Galactic WLR is found to be much steeper than expected from theory (although the specific sample is rather small), leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L similar to 10(4.75) L circle dot, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars would turn out to be overestimated as a result of wind clumping, but the degree of clumping would be rather independent of metallicity, the scalings derived in this study are expected to remain correct.

Blue luminous stars in nearby galaxies: Quantitative spectral analysis of M33 B-type supergiant stars

We present the detailed spectral analysis of a sample of M33 B-type supergiant stars, aimed at the determination of their fundamental parameters and chemical composition. The analysis is based on a grid of non-LTE metal line-blanketed model atmospheres including the effects of stellar winds and spherical extension computed with the code FASTWIND. Surface abundance ratios of C, N, and O are used to discuss the chemical evolutionary status of each individual star. The comparison of observed stellar properties with theoretical predictions of massive star evolutionary models shows good agreement within the uncertainties of the analysis. The spatial distribution of the sample allows us to investigate the existence of radial abundance gradients in the disk of M33. The comparison of stellar and H II region O abundances ( based on direct determinations of the electron temperature of the nebulae) shows good agreement. Using a simple linear radial representation, the stellar oxygen abundances result in a gradient of -0.0145 +/- 0.005 dex arcmin(-1) (or -0.06 +/- 0.02 dex kpc(-1)) up to a distance equal to similar to 1.1 times the isophotal radius of the galaxy. A more complex representation cannot be completely discarded by our stellar sample. The stellar Mg and Si abundances follow the trend displayed by O abundances, although with shallower gradients. These differences in gradient slope cannot be explained at this point. The derived abundances of the three alpha-elements yield solar metallicity in the central regions of the disk of M33. A comparison with recent planetary nebula data from Magrini and coworkers indicates that the disk of M33 has not suffered from a significant O enrichment in the last 3 Gyr.

Wolf-Rayet optically thick winds with Alfven waves

The Wolf-Rayet (WR) stars are hot luminous objects which are suffering an extreme mass loss via a continuous stellar wind. The high values of mass loss rates and high terminal velocities of the WR stellar winds constitute a challenge to the theories of radiation driven winds. Several authors incorporated magnetic forces to the line driven mechanism in order to explain these characteristics of the wind. Observations indicate that the WR stellar winds may reach, at the photosphere, velocities of the order of the terminal values, which means that an important part of the wind acceleration occurs at the optically thick region. The aim of this study is to analyze a model in which the wind in a WR star begins to be accelerated in the optically thick part of the wind. We used as initial conditions stellar parameters taken from the literature and solved the energy, mass and momentum equations. We demonstrate that the acceleration only by radiative forces is prevented by the general behavior of the opacities. Combining radiative forces plus a flux of Alfven waves, we found in the simulations a fast drop in the wind density profile which strongly reduces the extension of the optically thick region and the wind becomes optically thin too close its base. The understanding how the WR wind initiate is still an open issue. (C) 2010 COSPAR. Published by Elsevier Ltd. All rights reserved.