Research Note
Improved velocity law parameterization for hot star winds
^{1}
Ústav teoretické fyziky a astrofyziky PřF MU, 611 37 Brno, Czech Republic
email: krticka@physics.muni.cz
^{2}
Astronomický ústav, Akademie věd České republiky, 251 65 Ondřejov, Czech Republic
Received: 9 February 2011
Accepted: 29 August 2011
The velocity law of hot star winds is usually parameterized via the socalled beta velocity law. Although this parameterization stems from theoretical considerations, it is not the most accurate description of the wind velocity law that follows from hydrodynamical calculations. We show that the velocity profile of our hydrodynamical wind models is described much better by polynomial approximation. This approximation provides a better fit than the beta velocity law already for the same number of free parameters.
Key words: stars: winds, outflows / stars: massloss / stars: earlytype / hydrodynamics / stars: atmospheres
© ESO, 2011
Hot stars lose an important part of their mass via linedriven winds. Consequently, massloss rate prescriptions are an important part of evolutionary codes. Such prescriptions can be derived from theoretical modelling of hot star winds. However, any theoretical prescription should be tested against the massloss rates derived from observations. Although there are several possibilities for deriving wind massloss rates (e.g., from Hα line profiles, from ultraviolet P Cygni line profiles, and from radio continua), knowing of the wind velocity law is needed to derive the massloss rates from observations. Although there exist sophisticated hydrodynamic calculations giving the wind velocity law in a tabular form (e.g., Gräfener & Hamann 2005; Krtička & Kubát 2010), analytic formulae are useful in many analyses. Consequently, the question of proper velocity law parameterization is important for analysing observed hot star wind spectra. Here we develop a parameterization that is as close as possible to the results of hydrodynamical calculations.
The most frequently used parameterization of the wind velocity law is the socalled beta velocity law, (1)where v_{∞,β} is the wind terminal velocity, β a free parameter describing the steepness of the velocity law, and R_{∗} the stellar radius. Thus, for parameterization of the wind velocity law using Eq. (1) we need two parameters, namely β and the wind terminal velocity v_{∞,β}. The popularity of the beta velocity law Eq. (1) stems from its simplicity and from the belief that the velocity profile described by Eq. (1) is close to the real one. Also simple analytical wind solutions have the form of Eq. (1). For example, in the case of accelerating force proportional to ~ 1/r^{2}, the solution of the momentum equation (Milne 1926; Chandrasekhar 1934; Rublev 1965) leads to , implying β = 0.5 in Eq. (1). Also the velocity law of a classical Castor et al. (1975) linedriven wind model can be approximated with β = 0.5 (cf., Gayley 2000). Motivated by this, observational studies widely use Eq. (1) and usually consider β as a free parameter, yielding a typical value of β = 0.7 − 1 for O star winds (e.g., Puls et al. 1996). However, for some stars (e.g., early B supergiants), higher β up to 3 may be determined (e.g. Crowther et al. 2006).
Bestfit parameters to the radial velocity from hydrodynamic models (Krtička & Kubát 2010) for selected stars and individual velocity parameterizations.
Several variants of the simple formula Eq. (1) can be found in the literature and generally aim to improve the fits of wind line profiles. For example, Hillier (1988) introduced the velocity law that also encompasses the wind part close to the hydrostatic photosphere. Subsequently, Hillier & Miller (1999) introduced a second component in the beta law, which provides better fits to the blue edges of P Cygni line profiles. This was also supported by Gräfener & Hamann (2005), who found that different values of β describe the wind velocity law in different regions. The failure of the simple βvelocity law to describe the wind velocity law was also mentioned by Gayley et al. (1995).
Although the original beta velocity law was derived theoretically, there is no reason to believe that it provides the best approximation to the real situation. It is well known from hydrodynamical simulations (Feldmeier et al. 1997; Owocki & Puls 1999) that the wind velocity field is variable with time, and consequently the beta law may at most describe mean flow. But even in stationary models, the wind velocity law may deviate from Eq. (1), at least because the wind is driven by different ions close to the star and at large radii (Vink et al. 2001; Krtička 2006).
It would be highly desirable to develop formulae that are able to describe the velocity law in a more accurate form than Eq. (1). It is not clear how to develop a more general formula based on Eq. (1). Luckily, the problem of functional fitting is fully settled in mathematics, where the concept of orthogonal functions was introduced. Consequently, we propose to use the set of polynomials that may form an orthogonal basis in a given functional space, e.g., (2a)from which the terminal velocity follows as . Equation (2a) can be equivalently rewritten in terms of Legendre polynomials. This guarantees that the derived coefficients decrease with increasing order i and that the coefficients do not significantly change with increasing n (e.g., Mikulášek 2007). In this case instead of Eq. (2a), the parameterizations is (2b)and the terminal velocity . To ensure the orthogonality in the interval [0,1 ] we use shifted Legendre polynomials here. The first four are (e.g., Koornwinder et al. 2010) (3)The relations between coefficients and v_{i} follow from Eqs. (2) and (3) (4)Only three of these () coefficients are independent.
Fig. 1 Comparison of the fit to the calculated hydrodynamical model radial velocity of HD 24912 (crosses) by the beta velocity law (solid line), and by the polynomial fit Eq. (2b) for n = 3. 

Open with DEXTER 
Fig. 2 The mean relative difference (see text) between the velocity law (upper panel) and its derivative v′ = dv/dr (bottom panel) taken from hydrodynamical models and its approximations averaged over stars from Table 1 as a function of radius. 

Open with DEXTER 
To demonstrate the convenience of the formulas Eq. (2) we fit the radial velocity structure obtained from hydrodynamical NLTE wind models with CMF line force of Krtička & Kubát (2010, see Fig. 1. The bestfit parameters of the beta velocity law and polynomial fit Eq. (2b) with n = 3 for individual stars are given in Table 1. To test the global accuracy of individual orders of the velocity approximation Eq. (2), we plot the relative difference between the velocity from the numerical models v^{mod} and its approximations v^{app} (Eq. (2b) for n = 1,2,3 and Eq. (1)) averaged over individual stars in Table 1, ⟨  (v^{mod} − v^{app})/v^{app}  ⟩ , as a function of radius in Fig. 2.
Already the first term in Eq. (2a), i.e. the approximation v(r) = v_{1}^{(}1 − R_{∗}/r^{)} (in fact Eq. (1) for β = 1), provides reasonable approximation to the wind velocity law suitable for simplified analysis (see Fig. 2). Including the second term in Eq. (2a) we obtain better approximation than the usual beta velocity law both for the velocity and its derivative using the same number of free parameters (β and v_{∞,β} for the beta velocity law, v_{1} and v_{2} for Eq. (2a)). The agreement in the outer parts of the wind (r/R_{∗} ≳ 1.1) can be improved further by adding the third term in Eq. (2a) (see Fig. 2).
It has not escaped our attention that even relatively high order polynomials do not provide reasonable fits in the region close to the star, for r/R_{∗} ≲ 1.1. We do not aim to fit the velocity in this region, because here the wind density approaches the hydrostatic density stratification. Consequently, it would be better to derive the velocity from the continuity equation using density from the hydrostatic equilibrium equation instead of a polynomial.
Many studies fit the observed line profiles of particular stars using sophisticated radiative transfer calculations (e.g., Puls et al. 1996; Hillier & Miller 1999) by assuming the velocity law of type Eq. (1) and trying to find the most suitable value of the parameter β in combination with the value of the terminal velocity v_{∞}. Although our formulae Eq. (2) were derived to fit the results of numerical hydrodynamical calculations, it might be interesting to use them to also fit the beta velocity profiles. The test showed that the polynomial approximation Eq. (2) with n = 2 provides very good fits to beta laws (Eq. (1)) both with small β < 1 and large β > 1 and to the combined law of Hillier & Miller (1999). Only for large β ≳ 2.3 is it better to increase the polynomial degree (n = 3) to assure the positivity of the fit close to the star. The modified beta law with additional parameter b < 1 (Puls & Hummer 1988) is frequently used to improve the fit close to the star. It can be satisfactorily reproduced by either introducing a fixed low value of v_{0} in Eq. (2a) or by a
substitution x = 1 − bR_{ ∗ }/r in Eq. (2b). Moreover, our experiments with fitting the trial P Cygni line profiles have shown that the fits using Eq. (2) and the beta velocity law require a comparable number of line profile evaluations. A reasonable fit strategy could be to derive v_{∞} from observations and then keep the sum fixed.
We conclude that the beta velocity law provides reasonable approximation to the hydrodynamical wind velocity law and its derivative. The approximation can be improved with polynomial fits Eq. (2a) (or, equivalently Eq. (2b)) using the same number of free parameters as in the commonly used beta velocity law.
Acknowledgments
We thank Prof. Achim Feldmeier and Dr. Joachim Puls for discussions of this subject. This work was supported by grant GA ČR 205/08/0003. The Astronomical Institute Ondřejov is supported by the project AV0 Z10030501.
References
 Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157 (CAK) [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1934, MNRAS, 94, 522 [NASA ADS] [Google Scholar]
 Crowther, P. A., Lennon, D. J., & Walborn, N. R. 2006, A&A, 446, 279 [Google Scholar]
 Feldmeier, A., Puls, J., & Pauldrach, A. W. A. 1997, A&A, 322, 878 [NASA ADS] [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]
 Hillier, D. J. 1988, ApJ, 327, 822 [NASA ADS] [CrossRef] [Google Scholar]
 Hillier, D. J., & Miller, D. L. 1999, ApJ, 519, 354 [NASA ADS] [CrossRef] [Google Scholar]
 Koornwinder, T. H., Wong, R. S. C., Koekoek, R., & Swarttouw, R. F. 2010, in NIST Handbook of Mathematical Functions, ed. F. W. J. Olver, D. M. Lozier, R. F. Boisvertet et al. (Cambridge: Cambridge University Press) [Google Scholar]
 Krtička, J. 2006, MNRAS, 367, 1282 [NASA ADS] [CrossRef] [Google Scholar]
 Krtička, J., & Kubát, J. 2010, A&A, 519, A50 [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., et al. 2005, A&A, 441, 735 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mikulášek, Z. 2007, Odessa Astronomical Publications, 20, 138 [NASA ADS] [Google Scholar]
 Milne, E. A. 1926, MNRAS, 86, 459 [NASA ADS] [Google Scholar]
 Owocki, S. P., & Puls, J. 1999, ApJ, 510, 355 [NASA ADS] [CrossRef] [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]
 Puls, J., & Hummer, D. G. 1988, A&A, 191, 87 [NASA ADS] [Google Scholar]
 Puls, J., Kudritzki, R.P., Herrero, A., et al. 1996, A&A, 305, 171 [NASA ADS] [Google Scholar]
 Repolust, T., Puls, J., & Herrero, A. 2004, A&A, 415, 349 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rublev, S. V. 1965, AZh, 42, 718 (SvA, 9, 555) [NASA ADS] [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]
All Tables
Bestfit parameters to the radial velocity from hydrodynamic models (Krtička & Kubát 2010) for selected stars and individual velocity parameterizations.
All Figures
Fig. 1 Comparison of the fit to the calculated hydrodynamical model radial velocity of HD 24912 (crosses) by the beta velocity law (solid line), and by the polynomial fit Eq. (2b) for n = 3. 

Open with DEXTER  
In the text 
Fig. 2 The mean relative difference (see text) between the velocity law (upper panel) and its derivative v′ = dv/dr (bottom panel) taken from hydrodynamical models and its approximations averaged over stars from Table 1 as a function of radius. 

Open with DEXTER  
In the text 