Issue 
A&A
Volume 519, September 2010



Article Number  A50  
Number of page(s)  9  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201014111  
Published online  10 September 2010 
Comoving frame models of hot star winds
I. Test of the Sobolev approximation in the case of pure line transitions
J. Krticka^{1}  J. Kubát^{2}
1  Ústav teoretické fyziky a astrofyziky PrF MU, 611 37 Brno, Czech
Republic
2  Astronomický ústav, Akademie ved Ceské republiky, 251 65 Ondrejov,
Czech Republic
Received 21 January 2010 / Accepted 28 April 2010
Abstract
We provide hot star wind models with radiative force calculated using
the solution of comoving frame (CMF) radiative transfer equation. The
wind models are calculated for the first stars, O stars, and
the central stars of planetary nebulae. We show that without line
overlaps and with solely thermal line broadening the pure Sobolev
approximation provides a reliable estimate of the radiative force even
close to the wind sonic point. Consequently, models with the Sobolev
line force provide good approximations to solutions obtained with
nonSobolev transfer. Taking line overlaps into account,
the radiative force becomes slightly lower, leading to a
decrease in the wind massloss rate by roughly 40%. Below the
sonic point, the CMF line force is significantly lower
than the Sobolev one. In the case of pure thermal broadening,
this does not influence the massloss rate, as the wind
massloss rate is set in the supersonic part of the wind. However, when
additional line broadening is present (e.g., the turbulent
one) the region of low CMF line force may extend outwards to the
regions where the massloss rate is set. This results in a
decrease in the wind massloss rate. This effect can at least partly
explain the low wind massloss rates derived from some observational
analyses of luminous O stars.
Key words: stars: winds, outflows  stars: massloss  stars: earlytype  hydrodynamics  radiative transfer
1 Introduction
One of the most important galactic populations consists of massive stars, because these stars dominate the spectra of many galaxies and contribute significantly to the mass and momentum input into the interstellar matter. Moreover, massive stars end their active lives in gigantic explosions such as supernovae or even possibly as the progenitors of gammaray bursts (see Yoon & Langer 2005; Woosley & Heger 2006), producing huge amounts of heavier elements.
An important property of hot stars that significantly influences their final stages is the stellar wind (see, e.g., Puls et al. 2008b; Owocki 2004; Krticka & Kubát 2007a, for reviews dedicated to hot star winds). However, in stellar evolution calculations it is usually unnecessary to know detailed wind properties, but just the amount of mass expelled from the star per unit of time (massloss rate) as a function of stellar parameters (e.g., mass, effective temperature, radius, surface metallicity). However, for many hot stars we simply cannot estimate their true massloss rate with the precision necessary to calculate evolutionary models. The situation may be less problematic for luminous O stars, for which relatively good agreement between theoretical predictions and observational results seems to exist (Vink et al. 2001; Krticka & Kubát 2004; Pauldrach et al. 2001, hereafter nltei).
However, the agreement between theoretically predicted massloss rates and those derived from observations may be an illusion caused by the neglect of some physical effects in the wind, such as clumping (Bouret et al. 2003; Martins et al. 2005). As a result, the true massloss rates of O stars may be a few times lower than the standard wind theory predicts. This seems to be supported by the observations of Fullerton et al. (2006) of weak wind line profiles of P V. Last but not least, the unexpected occurrence of symmetrical Xray line profiles seems to require relatively low wind massloss rates (Kramer et al. 2003).
Possibly unreliable estimates of hot star wind massloss rates are also problematic because altough more realistic evolutionary stellar models can be calculated, by including, e.g., rotation and magnetic fields, the wind massloss rates remain uncertain. Ideally, all observational indicators of massloss rate and theoretical models should find and predict similar massloss rates. From the observational point of view, more detailed models of line formation in inhomogeneous media may be necessary to obtain reliable line profiles, and consequently also estimate massloss rates (Oskinova et al. 2007; Sundqvist et al. 2010).
From the theoretical point of view, disagreement between theory and observations would imply that some of the assumptions used for the hotstar wind modeling are inaccurate. Part of the disagreement may be caused by using incorrect abundances (Krticka & Kubát 2007b), although the reason for a disagreement remains mainly unclear. A thorough inspection of all assumptions involved in the modeling is therefore strongly needed. As a first step in this direction, we studied the influence of Xrays on the wind structure of hot stars. It seems that Xrays alone cannot entirely explain the disagreement between theory and observations (Krticka & Kubát 2009) as their influence on wind massloss rates is small and they do not strongly affect the ionization fraction of many important ions, especially that of P V. On the other hand, the modified ionization equilibrium may affect the Xray line formation (Krticka & Kubát 2009; Oskinova et al. 2006), and too ling cooling time in the postshock region (Krticka & Kubát 2009; Cohen et al. 2008) may cause the socalled ``weak wind problem'' (Martins et al. 2004; Bouret et al. 2003; Marcolino et al. 2009).
One of the most important approximations in the hotstar wind modeling is the Sobolev approximation (Castor 1974; Sobolev 1947), which enables us to solve the line radiation transfer analytically. Some studies confirm its applicability in the supersonic part of smooth linedriven winds (Pauldrach et al. 1986; Hamann 1981; Puls 1987). However, the applicability of the Sobolev approximation is questionable especially in the regions close to the photosphere bacause of the existence of strong source function gradients (Owocki & Puls 1999). On the other hand, some models avoid using the Sobolev approximation and use only the comovingframe (hereafter CMF) method of solving the radiative transfer equation (e.g., Gräfener & Hamann 2005).
We decided to test the applicability of the Sobolev approximation and include the CMF solution of the radiative transfer equation in our wind models. In this first paper of a series, we describe our method, and study the applicability of the Sobolev approximation using models neglecting continuum opacity.
2 Basic model assumptions
The models used in this paper are based on the NLTE wind models of Krticka & Kubát (2004, hereafter <)375#>nltei#. Here we summarise only their basic features and describe the inclusion of CMF line force.
We assume a spherically symmetric stationary stellar wind. The excitation and ionization state of the considered elements is derived from the statistical equilibrium (NLTE) equations. Ionic models are either adopted from the TLUSTY grid of model stellar atmospheres (Lanz & Hubeny 2007,2003) or are prepared by us using the data from the Opacity and Iron Projects (Nahar & Pradhan 1993; Chen & Pradhan 1999; Butler et al. 1993; Seaton et al. 1992; Zhang 1996; Nahar & Pradhan 1996; Luo & Pradhan 1989; Bautista 1996; Seaton 1987; Bautista & Pradhan 1997; Fernley et al. 1987; Sawey & Berrington 1992; Zhang & Pradhan 1997; Hummer et al. 1993). As in nltei, the solution of the radiative transfer equation for NLTE equations is artificially split into two parts, namely the radiative transfer in either the continuum or in lines. The solution to the radiative transfer equation in continuum is based on the Feautrier method in the spherical coordinates (Kubát 1993; Mihalas & Hummer 1974), and the line radiative transfer is solved in the Sobolev approximation (Castor 1974; Rybicki & Hummer 1978) neglecting continuum opacity and line overlaps.
In contrast to our previous models, the radiative transfer in lines used for the calculation of the radiative force is solved in the CMF (see Sect. 3) neglecting the continuum opacity. The line radiative force is calculated directly from the true chemical composition, NLTE ionization and excitation balance, and CMF flux using data from the VALD database (Piskunov et al. 1995; Kupka et al. 1999). We do not use the linestrength distribution function parameterized by force multipliers k, , and .
The flux at the surface (used as the lower boundary condition for the radiative transfer in the wind) is taken from the HHe spherically symmetric NLTE model stellar atmospheres of Kubát (2003, and references therein).
3 CMF calculation of the radiative force
The radiative force is calculated using the solution of the spherically symmetric CMF radiative transfer equation (Mihalas 1978, Eq. (14.99))
where is the intensity seen by the observer moving with the wind at the radial velocity , is the frequency, , is the direction between the given ray and the radial direction, and and are the line emissivity and opacity, respectively, given by
where n_{i}
and n_{j}
are the number densities of individual states with statistical
weights g_{i}
and g_{j}
corresponding to the line transition
with oscillator strength f_{ij}
and the lineprofile
,
and
is the electron mass. Assuming the line profile to be Gaussian produced
by thermal broadening only, is
given by
(3) 
where is the laboratory line frequency, and the line broadening is given by
and is the mass of a given atom. The number densities of individual levels in Eqs. (2) are calculated from statistical equilibrium (NLTE) equations.
Writing Eq. (1), we neglected advection and aberration terms, which is justifiable in nonrelativistic flows (see, e.g., Korcáková & Kubát 2003). We also note that by neglecting the spatial derivatives of intensity in Eq. (1) we obtain the Sobolev approximation (Castor 2004).
Following Mihalas
et al. (1975), we rewrite Eq. (1) for rays with an
impact parameter p
where + and  refers to radiation flowing toward and away from the observer, respectively, r=(p^{2}+z^{2})^{1/2}, and z is the distance along the ray. We transform Eq. (5) using intensitylike and fluxlike variables
(6a)  
(6b) 
to obtain a system of partial differential equations
where
(8)  
(9)  
(10)  
(11) 
The system of equations in Eq. (7) is solved numerically using the long characteristic method of Mihalas et al. (1975), which we modified slightly for the present purpose (see Appendix A). As we are interested in the calculation of the radiative force using the v variable at a particular grid point, in contrast to Mihalas et al. (1975) we specify v at grid points, and u in the middle between them.
In our numerical solution of Eq. (7), we use the same spatial grid as for the solution of hydrodynamical equations. The spacing of the frequency grid is , where is the prespecified expected minimum wind temperature, is the atomic mass of artificial metallic atom, and is the multiplicative factor (see below). The CMF radiative transfer equation is solved only for selected frequencies from the frequency grid that lie close to some line. The selection of frequencies is controlled by two integer numbers , and NCERV (cf. Hillier & Miller 1998). For each line, we select frequencies that lie within line Doppler widths . Redward of the center of each line, we select each NCERV frequency up to the frequency corresponding to the Doppler shift for the wind terminal velocity. The numerical test showed that a sufficiently precise value of the radiative force can be derived for the value of parameters , where is the mass of hydrogen atom, , , NCERV = 30, and typically K.
The radiative force is calculated as an integral
As the calculation of the CMF radiative force is rather timeconsuming, we do not calculate during each iteration of hydrodynamical variables, but adopt a different approach. We calculate the ratio of the CMF and Sobolev line forces
By the Sobolev line force , we mean here the force calculated by assuming the Sobolev approximation for radiative transfer, neglecting line overlaps and using true line opacities and the emergent flux from the underlying stellar atmosphere [Eq. (25) therein]nltei. Unless the base density is known with a precision better than about 30%, we calculate only when the estimate of the base density is changed, and keep fixed during the subsequent iterations of the hydrodynamical structure. When the base density is known with a higher precision, we calculate after each change of the hydrodynamical structure. Moreover, because we solve the hydrodynamical equations using the NewtonRaphson method, we have to calculate the derivatives of with respect to individual hydrodynamical variables. These derivatives are approximated using the derivatives of the Sobolev line force multiplied by . The force term in the critical point condition [see][]nltei is also multiplied by .
We note that direct use of Eq. (13) causes
instability in the model convergence. The reason for these convergence
problems may be numerical, but this behavior may also be connected with
linedriven instability (Feldmeier et al. 1997; Owocki et al.
1988). To avoid this (since we are seeking
stationary solution and not evolution with time) we introduced a weak
smoothing of
,
where is the value of at a given grid point d (as for d1 and d+1) and we use instead of in the models. Our numerical tests showed that the smoothing Eq. (14) does not significantly affect the resulting radiative force.
Table 1: Radius R_{*}, mass M, and the effective temperature of studied model stars.
4 Studied model stars
In our study, we selected three types of stars to study more carefully the Sobolev approximation in different wind environments (see Table 1).
The stellar parameters of the first stars were obtained according to an evolutionary calculation of initially zerometallicity star with initial mass derived by Marigo et al. (2001). For these models, we assumed a stellar wind driven purely by CNO elements (which appear on the stellar surface due to mixing) with a massfraction of CNO Z=10^{3}.
The stellar parameters (effective temperatures and radii) of an O star sample were derived using the model atmospheres with line blanketing (Markova et al. 2004; Martins et al. 2005; Repolust et al. 2004). Stellar masses were obtained using evolutionary tracks either by ourselves (using Schaller et al. 1992 tracks) or by Martins et al. (2005). For these stars, we assumed a solar chemical composition (Asplund et al. 2005).
The stellar parameters of central stars of planetary nebulae were taken from Pauldrach et al. (2004), who derived them from UV spectroscopy. Helium abundance was adopted from Kudritzki et al. (1997), for other elements we assumed a solar chemical composition (after Asplund et al. 2005), which was for some stars slightly modified according to Pauldrach et al. (2004).
5 Comparison of CMF and Sobolev wind models
We calculated wind models with both CMF and Sobolev line forces and compared the final wind structure. The resulting ratio of the CMF to Sobolev line forces for selected stars is shown in Fig. 1. We note that the Sobolev force was calculated using the flux from the stellar atmosphere and by neglecting line overlaps.
Figure 1: The ratio of the CMF to Sobolev line forces given by Eq. (13) as a function of the wind velocity plotted in the terms of the sound speed for two selected stars. The dasheddotted line denotes a model with constant source function and constant level populations equal to their values at the sonic point. Arrows indicate the thermal speed of selected ions. 

Open with DEXTER 
For very low wind velocities (where ), the CMF force is large, . This is most likely partly connected with the boundary conditions, which are not completely compatible with the wind (cf., Noerdlinger & Rybicki 1974).
For velocities of about one tenth of the sound speed, there is an apparent minimum of . In some cases, the ratio could even be negative, which corresponds to a negative CMF radiative force. The Sobolev approximation is not applicable in this region, but a low value of the radiative force is also connected with positive source function gradients. For subsonic velocities, the Doppler shift is less important, and the line radiative transfer is given basically by the static radiative transfer equation. In the optically thick regions, for frequencies corresponding to line transitions it follows from Eq. (7b) , from Eq. (7a) , and the radiative force is proportional to the negative of the derivative of the source function (see Eq. (12), and Noerdlinger & Rybicki 1974). Because the line source function increases here (see Fig. 2), the line radiative force at low velocities may even be negative. For a constant source function, the minimum of close to the star is significantly weaker (see Fig. 1). The source function minimum below the sonic point is caused by a local temperature minimum, because the line source function of optically thick lines (which are, consequently, in detailed radiative balance) close to the star (asterisk denotes LTE value) depends on temperature. Another source function minimum for nonSobolev source function due to velocity field curvature was also found by Sellmaier et al. (1993), and Owocki & Puls (1999). We note that in the case of the resonance lines plotted in Fig. 2 the line source function at larger radii is roughly proportional to (e.g., Kudritzki & Puls 2000).
The minimum of close to the star is also connected with the velocity gradient changing significantly within the resonance zone. Thus, a given line also picks up the radiation corresponding to a lower velocity gradient leading to a further reduction in the radiative force. For velocities comparable to or higher than the ion thermal speed, the lines are deshadowed because of the Doppler effect, and the radiative force increases. We note that we only include the thermal broadening, hence these effects occur for velocities lower than the sound speed.
Figure 2: The line source function for resonance lines of O IV at 790 Å, S V at 786 Å, and Fe V at 388 Å as a function of relative wind velocity. The source function is plotted relative to its value at the sonic point for the wind model of HD 93204. 

Open with DEXTER 
As the wind accelerates, the ratio of the CMF to Sobolev line force increases and reaches a value close to one for velocities higher than the thermal speed of the wind driving ions (Fig. 1). This is unsurprising, because the Sobolev approximation is applicable to regions with a large velocity gradient, which exist already close to the sonic point . Owing to line overlaps, is less than one in the outer wind regions, where it reaches only 0.70.8.
Figure 3: The radial variation in the ratio of the CMF to Sobolev line forces in the wind model of Tc 1 central star with and without line overlaps. 

Open with DEXTER 
To test the influence of line overlaps, we calculated the radiative force with only 50 carefully selected optically thick lines that do not overlap (see Fig. 3). The pronounced minimum for velocities lower than the sound speed is still present here, but in the outer regions the value of is approximately one, supporting the validity of the Sobolev approximation for supersonic velocities.
Figure 4: The ratio of the CMF to Sobolev line forces at radius r=5 R_{*}in the Tc1 wind model in the dependence on the line shift. The value of hydrogen thermal speed is denoted in the figure. 

Open with DEXTER 
To understand more clearly the influence of line overlaps on the radiative force, we constructed another artificial line list using our set of nonoverlapping lines. Each line in this new line list is counted twice with all parameters being completely the same, however with a line center shifted by , where is a free parameter. For , all twin lines completely overlap leading to a significant decrease in the radiative force with respect to the Sobolev one that does not account for the line overlaps (see Fig. 4). For , the lines at a given point do not overlap, but one of the twin lines ``sees'' the flux absorbed by the second line, leading to a reduction in the radiative force even in this case. For , one of the lines is affected by the emission from the second one, leading to an increase in the CMF radiative force relative to the Sobolev one.
A similar reduction in the line force by multiline effects was found by Puls (1987). We note that the multiline effects were also studied with respect to the multiple radiative momentum deposition in WolfRayet star winds (Gayley et al. 1995). However, these effects are probably of minor importance here due to the low density of the studied winds.
Table 2: Comparison of calculated wind parameters derived using CMF and Sobolev line forces.
The CMF radiative force, which is lower than the Sobolev one because of line overlaps, causes a decrease in the massloss rate of CMF models with respect to Sobolev ones (see Table 2). The ratio of CMF to Sobolev massloss rates is about 0.58. The only exception is the model M5001, for which the CMF massloss rate is nearly the same as the Sobolev one. The reason is that the star is so hot, that the wind is accelerated mainly by a dozen O V and O VI lines. For a critical point velocity, these lines do not overlap, hence , and the CMF and Sobolev massloss rates are nearly the same.
Figure 5: Comparison of our derived massloss rates ( upper panel) and terminal velocities ( lower panel) of the central stars of planetary nebulae with those derived by Pauldrach et al. (2004). 

Open with DEXTER 
The resulting wind parameters of the central stars of planetary nebulae can be compared with those derived from observations by Pauldrach et al. (2004, see Fig. 5). There is reasonable agreement between the wind parameters predicted by ourselves and those derived by Pauldrach et al. (2004, see Fig. 5). The massloss rates of Pauldrach et al. (2004) are on average a factor of about 1.6 higher than those derived by ourselves. This is most likely partly because of the simplifications included in our code, e.g., the neglect of continuum opacity sources, and partly by the different abundances adopted.
6 Models with base turbulence
Figure 6: The massloss rate of HD 209975 in the models with additional turbulent line broadening relative to the models with zero turbulent velocity. 

Open with DEXTER 
The existence of a region close to the stellar surface where the CMF line force is low compared to the Sobolev one (see Fig. 1) is partly caused by the source function gradients at the wind base and partly by the Sobolev approximation not being applicable to the subsonic regions. The CMF line force increases at the moment when the line starts to absorb the radiation that has not been absorbed yet, i.e., the radiation from the line wing. Because up to now we have assumed pure thermal line broadening, the velocity width of low CMF line force is of the order of the metallic thermal speed (which is roughly 0.13 a in the case of iron). The wind massloss rate in our models is determined in the region of supersonic wind, close to the critical point where the wind velocity approaches the speed of radiativeacoustic waves (Feldmeier et al. 2008; Abbott 1980). Thus, the region of low CMF line force close to the star does not significantly affect the wind massloss rate.
However, if the line broadening were larger (due to
surface turbulence), then the region of low
CMF line force could spread out to large velocities comparable
to the turbulent one. When the turbulent velocity is comparable to the
critical point velocity, below which the wind massloss rate
is set, this could cause a significant decrease in the wind massloss
rate. To test this, we calculated wind models with additional
line broadening, which we attributed to the turbulent one.
In this case, the line profile width is given not by
Eq. (4),
but by
where is the adopted turbulent velocity.
The results of numerical models indicate that with increasing turbulent broadening the velocity width of low CMF line force increases leading to a lower massloss rate (see Fig. 6, cf. Lucy 2007). Hence, in the presence of turbulence the wind parameters may not depend only on the basic stellar parameters (effective temperature, radius, mass) but also on the line turbulent broadening. Moreover, this effect can possibly be one of the reasons why the massloss rates derived from observational analyses that take the clumping into account (Bouret et al. 2003; Martins et al. 2005) are systematically lower than the predicted ones.
For velocities higher than a few times the turbulent one, the Sobolev approximation should be applicable. At these high velocities, one expects that the line force becomes close to the Sobolev one. Because now the same force (as in the model with zero turbulent broadening) accelerates the wind of lower density, one expects the terminal velocity to increase (e.g., Gayley 2001), becoming much higher than the observed one. However, our models do not predict a significant increase in the terminal velocity , which is in the range for the wind models of HD 209975 with different turbulent broadening. This is caused by the stronger blocking of stellar radiation by increased line overlaps mainly in the region with . We note that lines broadened by turbulent motions are able to block the flux more efficiently than lines broadened purely thermally.
Observational studies consider the turbulence already present in the photospheres of O stars (e.g., Bouret et al. 2005; Martins et al. 2004; Bouret et al. 2003; Martins et al. 2005) with turbulent velocities of about . Macroturbulent velocities in B supergiants may be even higher, about (Markova & Puls 2008; Howarth et al. 1997). Convective layers and surface pulsational motions are also expected theoretically (Cantiello et al. 2009; Aerts et al. 2009). Turbulence can spread in the wind (Feldmeier et al. 1997), leading to a decrease in the wind massloss rate, as shown here. We also note that many O stars exhibit turbulent velocities in the range , where we expect a high sensitivity of the predicted massloss rate to the turbulent velocity (see Fig. 6).
The basic results presented here will, in the future, be tested in more detail using models that also account for the continuum opacity and CMF line source function in a separate study.
7 The solution topology
Figure 7: The dependence of the radial velocity on radius for solutions with different boundary densities (massloss rates) close to the stellar surface for the model star 5001. Each consecutive model (from down to up) differs by a factor of 1.5 in the massloss rate (thin lines). The thick line denotes the unique solution that smoothly intercepts the critical point. The dashed line denotes solution with the Sobolev line force. 

Open with DEXTER 
In Fig. 7, we plot solutions with different base densities (massloss rates). In generall, with increasing base density the wind velocity increases until the density reaches a maximum value. There is no solution that is smooth out to large radii for the densities higher than the maximum one. The solution with maximum density is very similar to the critical solution of Sobolev models (see Fig. 7). Moreover, there are many solutions that smoothly pass through the sonic point v=a for different massloss rates.
This indicates that the critical point of nonSobolev models is close to the CAK critical point (Castor et al. 1975) and that the sonic point is not a point where the wind massloss rate is determined. The reason is that even in the nonSobolev models the radiative force is not given locally by wind density and velocity, but depends on the wind properties in a close neighborhood of a studied point. This dependence on the nonlocal properties at its limit approaches the Sobolev approximation for very thin resonance layers (for very large velocity gradients).
8 Conclusions
We have presented hot star wind models in which the radiative force is calculated using the solution of the comoving frame (CMF) radiative transfer equation. The wind models were calculated for three different groups of stellar parameters (corresponding to evolved first stars, O stars, and the central stars of planetary nebulae) to compare the CMF and Sobolev radiative forces for a broader range of stellar parameters.
The comparison of the CMF radiative force with an approximate one calculated by assuming the Sobolev approximation showed that the Sobolev line force is slightly higher due to the neglect of line overlaps. Thus, the massloss rate of wind models that include the Sobolev line force is on average a factor of about 1.7 higher than a more realistic one calculated using CMF wind models. However, we note that the simple Sobolev approximation applied here for reference can be improved to account for line overlaps and continuum absorption (Olson 1982; Hummer & Rybicki 1985; Pavlakis & Kylafis 1996; Puls & Hummer 1988). We emphasize that modern hot star wind models include line overlaps (e.g., Vink et al. 2001; Gräfener & Hamann 2005; Pauldrach et al. 2001).
Without line overlaps in the case of purely thermal line broadening, the Sobolev approximation provides a reliable estimate of the radiative force even close to the wind sonic point. The CAK model therefore provides a good approximation for a solution obtained with nonSobolev transfer. Below the sonic point, the CMF line force is significantly lower than the Sobolev one partly because of the strong gradients in the source function. This does not influence the massloss rate, as the wind massloss rate is set in the supersonic part of the wind below the critical point. However, when additional line broadening is present (e.g., the turbulent one) then the region of low CMF line force may extend outwards to the regions where the massloss rate is set. This results in a significant decrease in the wind massloss rate. We note that this is not a shortcoming of the Sobolev approximation because the Sobolev approximation is applicable to velocities higher than the turbulent velocity in this case.
The influence of turbulent line broadening may cause a dependence of the wind massloss rate on the atmospheric turbulent motions. Because theoretical models are not yet able to predict the atmospheric turbulent motions in hot stars in detail, we are unable to provide reliable wind massloss rate predictions until the theory of the atmospheric turbulence develops considerably (however see Aerts et al. 2009; Cantiello et al. 2009). Nowadays, hot star evolution seems to be a deterministic one depending only on the initial stellar parameters, i.e., mass, metallicity, and rotational rate. However, as the properties of atmospheric turbulent motions seem to be a nontrivial function of stellar parameters, the evolution of hot stars may become less deterministic, becoming instead dependent on free parameters describing the role of surface turbulence.
AcknowledgementsThis work was supported by grant GA CR 205/07/ 0031. The Astronomical Institute Ondrejov is supported by the project AV0 Z10030501.
References
 Abbott, D. C. 1980, ApJ, 242, 1183 [NASA ADS] [CrossRef] [Google Scholar]
 Aerts, C., Puls, J., Godart, M., & Dupret, M.A. 2009, A&A, 508, 409 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Anderson, E., Bai, Z., Bischof, C., et al. 1999, LAPACK Users' Guide, 3rd edn. (Philadelphia: SIAM) [Google Scholar]
 Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Cosmic Abundances as Records of Stellar Evolution and in Nucleosynthesis, ed. T. G. Barnes III, & F. N. Bash (San Francisco: ASP), ASP Conf. Ser., 336, 25 [Google Scholar]
 Bautista, M. A. 1996, A&AS, 119, 105 [Google Scholar]
 Bautista, M. A., & Pradhan, A. K. 1997, A&AS, 126, 365 [Google Scholar]
 Bouret, J.C., Lanz, T., Hillier, D. J., et al. 2003, ApJ, 595, 1182 [NASA ADS] [CrossRef] [Google Scholar]
 Bouret, J.C., Lanz, T., & Hillier, D. J. 2005, A&A, 438, 301 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Butler, K., Mendoza, C., & Zeippen, C. J. 1993, J. Phys. B., 26, 4409 [NASA ADS] [CrossRef] [Google Scholar]
 Cantiello, M., Langer, N., Brott, I., et al. 2009, A&A, 499, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Castor, J. I. 1974, MNRAS, 169, 279 [NASA ADS] [CrossRef] [Google Scholar]
 Castor, J. I. 2004, Radiation Hydrodynamics (Cambridge: Cambridge University Press) [Google Scholar]
 Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157 (CAK) [NASA ADS] [CrossRef] [Google Scholar]
 Chen, G. X., & Pradhan, A. K. 1999, A&AS, 136, 395 [Google Scholar]
 Cohen, D. H., Kuhn, M. A., Gagné, M., Jensen, E. L. N., & Miller, N. A. 2008, MNRAS, 386, 1855 [NASA ADS] [CrossRef] [Google Scholar]
 Feldmeier, A., Puls, J., & Pauldrach, A. W. A. 1997, A&A, 322, 878 [NASA ADS] [Google Scholar]
 Feldmeier, A., Rätzel, D., & Owocki, S. P. 2008, ApJ, 679, 704 [NASA ADS] [CrossRef] [Google Scholar]
 Fernley, J. A., Taylor, K. T., & Seaton, M. J. 1987, J. Phys. B, 20, 6457 [NASA ADS] [CrossRef] [Google Scholar]
 Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 [NASA ADS] [CrossRef] [Google Scholar]
 Gayley, K. G. 2000, ApJ, 529, 1019 [NASA ADS] [CrossRef] [Google Scholar]
 Gayley, K. G., Owocki, S. P., & Cranmer, S. R. 1995, ApJ, 442, 296 [NASA ADS] [CrossRef] [Google Scholar]
 Gräfener, G., & Hamann, W.R. 2005, A&A, 432, 633 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hamann, W.R. 1981, A&A, 93, 353 [NASA ADS] [Google Scholar]
 Hillier, D. J., & Miller, D. L. 1998, ApJ, 496, 407 [NASA ADS] [CrossRef] [Google Scholar]
 Howarth, I. D., Siebert, K. W., Hussain, G. A. J., & Prinja, R. K. 1997, MNRAS, 284, 265 [NASA ADS] [CrossRef] [Google Scholar]
 Hummer, D. G., & Rybicki, G. B. 1985, ApJ, 293, 258 [NASA ADS] [CrossRef] [Google Scholar]
 Hummer, D. G., Berrington, K. A., Eissner, W., et al. 1993, A&A, 279, 298 [NASA ADS] [Google Scholar]
 Korcáková, D., & Kubát, J. 2003, A&A, 401, 419 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kramer, R. H., Cohen, D. H., & Owocki, S. P. 2003, ApJ, 592, 532 [NASA ADS] [CrossRef] [Google Scholar]
 Krticka, J., & Kubát, J. 2001, A&A, 369, 222 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Krticka, J., & Kubát, J. 2004, A&A, 417, 1003 nltei [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Krticka, J., & Kubát, J. 2007a, in Active OBStars: Laboratories for Stellar & Circumstellar Physics, ed. S. Stefl, S. P. Owocki, & A. T. Okazaki (San Francisco: ASP), 153 [Google Scholar]
 Krticka, J., & Kubát, J. 2007b, A&A, 464, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Krticka, J., & Kubát, J. 2009, MNRAS, 394, 2065 [NASA ADS] [CrossRef] [Google Scholar]
 Kubát, J. 1993, Ph.D. Thesis, Astronomický ústav AV CR, Ondrejov [Google Scholar]
 Kubát, J. 2003, in Modelling of Stellar Atmospheres, ed. N. E. Piskunov, W. W. Weiss, & D. F. Gray (San Francisco: ASP), IAU Symp., 210, A8 [Google Scholar]
 Kudritzki, R. P., & Puls, J. 2000, ARA&A, 38, 613 [NASA ADS] [CrossRef] [Google Scholar]
 Kudritzki, R.P., Méndez, R. H., Puls, J., & McCarthy, J. K. 1997, in Planetary Nebulae, ed. H. J. Habing, & H. J. G. L. M. Lamers, Proc. IAU Symp., 180, 64 [Google Scholar]
 Kupka, F., Piskunov, N. E., Ryabchikova, T. A., Stempels, H. C., & Weiss, W. W. 1999, A&AS, 138, 119 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Lanz, T., & Hubeny, I. 2003, ApJS, 146, 417 [NASA ADS] [CrossRef] [Google Scholar]
 Lanz, T., & Hubeny, I. 2007, ApJS, 169, 83 [Google Scholar]
 Lucy, L. B. 2007, A&A, 468, 649 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Luo, D., & Pradhan, A. K. 1989, J. Phys. B, 22, 3377 [Google Scholar]
 Marcolino, W. L. F., Bouret, J.C., Martins, F., et al. 2009, A&A, 498, 837 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Marigo, P., Girardi, L., Chiosi, C., & Wood, P. R. 2001, A&A, 371, 152 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Markova, N., & Puls, J. 2008, A&A, 478, 823 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Markova, N., Puls, J., Repolust, T., & Markov, H. 2004, A&A, 413, 693 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Martins, F., Schaerer, D., Hillier, D. J., & HeydariMalayeri, M. 2004, A&A, 420, 1087 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Martins, F., Schaerer, D., Hillier, D. J., et al. 2005, A&A, 441, 735 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mihalas, D. 1978, Stellar Atmospheres (San Francisco: Freeman & Co.) [Google Scholar]
 Mihalas, D., & Hummer, D. G. 1974, ApJS, 28, 343 [NASA ADS] [CrossRef] [Google Scholar]
 Mihalas, D., Kunasz, P. B., & Hummer, D. G. 1975, ApJ, 202, 465 [NASA ADS] [CrossRef] [Google Scholar]
 Nahar, S. N., & Pradhan, A. K. 1993, J. Phys. B, 26, 1109 [NASA ADS] [CrossRef] [Google Scholar]
 Nahar, S. N., & Pradhan, A. K. 1996, A&AS, 119, 509 [Google Scholar]
 Noerdlinger, P. D., & Rybicki, G. B. 1974, ApJ, 193, 651 [NASA ADS] [CrossRef] [Google Scholar]
 Olson, G. L. 1982, ApJ, 255, 267 [NASA ADS] [CrossRef] [Google Scholar]
 Oskinova, L. M., Hamann, W.R., & Feldmeier, A. 2006, in High Resolution Xray Spectroscopy, ed. G. BranduardiRaymont, 27 [Google Scholar]
 Oskinova, L. M., Hamann, W.R., & Feldmeier, A. 2007, A&A, 476, 1331 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Owocki, S. P. 2004, in Evolution of Massive Stars, Mass Loss and Winds, ed. M. HeydariMalayeri, Ph. Stee, & J.P. Zahn, EAS PS, 13, 163 [Google Scholar]
 Owocki, S. P., & Puls, J. 1999, ApJ, 510, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Owocki, S. P., Castor, J. I., & Rybicki, G. B. 1988, ApJ, 335, 914 [NASA ADS] [CrossRef] [Google Scholar]
 Pauldrach, A., Puls, J., & Kudritzki, R. P. 1986, A&A, 164, 86 [NASA ADS] [Google Scholar]
 Pauldrach, A. W. A., Hoffmann, T. L., & Lennon, M. 2001, A&A, 375, 161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pauldrach, A. W. A., Hoffmann, T. L., & Méndez, R. H. 2004, A&A, 419, 1111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pavlakis, K. G., & Kylafis, N. D. 1996, ApJ, 467, 292 [NASA ADS] [CrossRef] [Google Scholar]
 Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525 [Google Scholar]
 Puls, J. 1987, A&A, 184, 227 [NASA ADS] [Google Scholar]
 Puls, J., & Hummer, D. G. 1988, A&A, 191, 87 [NASA ADS] [Google Scholar]
 Puls, J., Markova, N., & Scuderi, S. 2008a, in Mass Loss from Stars and the Evolution of Stellar Clusters, ed. A. de Koter, L. Smith, & R. Waters (San Francisco: ASP), 101 [Google Scholar]
 Puls, J., Vink, J. S., & Najarro, F. 2008b, A&ARv, 16, 209 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Repolust, T., Puls, J., & Herrero, A. 2004, A&A, 415, 349 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rybicki, G. B., & Hummer, D. G. 1978, ApJ, 219, 654 [NASA ADS] [CrossRef] [Google Scholar]
 Sawey, P. M. J., & Berrington, K. A. 1992, J. Phys. B, 25, 1451 [NASA ADS] [CrossRef] [Google Scholar]
 Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 [Google Scholar]
 Seaton, M. J. 1987, J. Phys. B, 20, 6363 [Google Scholar]
 Seaton, M. J., Zeippen, C. J., Tully, J. A., et al. 1992, Rev. Mex. Astron. Astrofis., 23, 19 [NASA ADS] [EDP Sciences] [Google Scholar]
 Sellmaier, F., Puls, J., Kudritzki, R. P., et al. 1993, A&A, 273, 533 [NASA ADS] [Google Scholar]
 Sobolev, V. V. 1947, Dvizhushchiesia obolochki zvedz, Leningr. Gos. Univ., Leningrad [Google Scholar]
 Sundqvist, J. O., Puls, J., & Feldmeier, A. 2010, A&A, 510, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Woosley, S., & Heger, A. 2006, ApJ, 637, 914 [NASA ADS] [CrossRef] [Google Scholar]
 Yoon, S.C., & Langer, N. 2005, A&A, 443, 643 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zhang, H. L. 1996, A&AS, 119, 523 [Google Scholar]
 Zhang, H. L., & Pradhan, A. K. 1997, A&AS, 126, 373 [Google Scholar]
Appendix A: The solution of CMF radiative transfer equation
To calculate the radiative force, the Mihalas et al. (1975) method for the solution of the CMF radiative transfer equation is modified in such a way that the v variable is specified on the spatial grid, and u on the intermediate one.
Following the notation of Mihalas et al. (1975), the depth index d increases inward, , where is the radius of the outer model boundary. The impact parameters p are labeled in order of increasing size by index j, , where NC is the number of rays intersecting the core. Along each ray with impact parameter p_{j}, we define grid in z and optical depth , , where for , and for , . The frequencies are labeled by index k in order of decreasing values, .
We assume that v is specified on the depth grid, and u is specified at intermediate grid points labeled by . Suppressing the ray index j in the following, we define on each ray p_{j}
(A.1)  
(A.2)  
(A.3) 
The difference form of the system of the equations in Eq. (7) is
where , and
(A.5)  
(A.6)  
(A.7) 
Solving Eq. (A.4b) for , we obtain
where
(A.9) 
Substituting Eqs. (A.8) into (A.4a), we derive a linear system of equations for
This system should be supplemented by equations corresponding to the boundary and initial conditions. At the outer spatial boundary , we assume no infalling radiation, consequently u=v and we derive from Eq. (7b)
(A.11) 
or, in a difference form
(A.12) 
The infalling radiation at the inner boundary is taken from the model atmospheres. The initial solution for is derived using the solution of the radiative transfer equation neglecting the velocity fields (Kubát 1993; Mihalas & Hummer 1974).
The velocity derivatives at the grid points are approximated as in the hydrodynamical code (see Krticka & Kubát 2001, Eq. (A.4a) therein). The derivatives in the middle points between grid points are calculated as the average of the derivatives at the grid points.
The system of algebraic equations Eq. (A.10) with boundary conditions is solved using the LAPACK package (http://www.cs.colorado.edu/ lapack, Anderson et al. 1999).
All Tables
Table 1: Radius R_{*}, mass M, and the effective temperature of studied model stars.
Table 2: Comparison of calculated wind parameters derived using CMF and Sobolev line forces.
All Figures
Figure 1: The ratio of the CMF to Sobolev line forces given by Eq. (13) as a function of the wind velocity plotted in the terms of the sound speed for two selected stars. The dasheddotted line denotes a model with constant source function and constant level populations equal to their values at the sonic point. Arrows indicate the thermal speed of selected ions. 

Open with DEXTER  
In the text 
Figure 2: The line source function for resonance lines of O IV at 790 Å, S V at 786 Å, and Fe V at 388 Å as a function of relative wind velocity. The source function is plotted relative to its value at the sonic point for the wind model of HD 93204. 

Open with DEXTER  
In the text 
Figure 3: The radial variation in the ratio of the CMF to Sobolev line forces in the wind model of Tc 1 central star with and without line overlaps. 

Open with DEXTER  
In the text 
Figure 4: The ratio of the CMF to Sobolev line forces at radius r=5 R_{*}in the Tc1 wind model in the dependence on the line shift. The value of hydrogen thermal speed is denoted in the figure. 

Open with DEXTER  
In the text 
Figure 5: Comparison of our derived massloss rates ( upper panel) and terminal velocities ( lower panel) of the central stars of planetary nebulae with those derived by Pauldrach et al. (2004). 

Open with DEXTER  
In the text 
Figure 6: The massloss rate of HD 209975 in the models with additional turbulent line broadening relative to the models with zero turbulent velocity. 

Open with DEXTER  
In the text 
Figure 7: The dependence of the radial velocity on radius for solutions with different boundary densities (massloss rates) close to the stellar surface for the model star 5001. Each consecutive model (from down to up) differs by a factor of 1.5 in the massloss rate (thin lines). The thick line denotes the unique solution that smoothly intercepts the critical point. The dashed line denotes solution with the Sobolev line force. 

Open with DEXTER  
In the text 
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.