EDP Sciences
Free Access
Volume 607, November 2017
Article Number L8
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201731902
Published online 17 November 2017

© ESO, 2017

1. Introduction

We present mass-loss predictions for helium (He) stars, including both massive Wolf-Rayet (WR) stars, and low-mass He stars that have presumably lost their hydrogen (H) envelope through binary mass transfer (e.g. Podsiadlowski et al. 1992).

One of the main ingredients for massive star evolution modelling concerns their rates of stellar wind mass loss (Langer 2012; Yoon et al. 2012; Georgy et al. 2013; Limongi & Chieffi 2006). For the H-burning main sequence, we are in the fortunate position that a theoretical recipe such as that by Vink et al. (2000, 2001) can be utilized in a meaningful way, but for more evolved stars, such as WR stars, the community still needs to rely on empirical recipes, such as the recipe by Nugis & lamers (2000). A notable drawback of employing empirical recipes is that they cannot be applied in a parameter space for which they have not been derived. These recipes might fail dramatically outside that parameter space, for example for low-mass He-stripped stars, especially when the recipe has not been constructed from the physically correct stellar parameters (Vink & de Koter 2002; Puls et al. 2008). One of the main aims of this study is to alleviate this shortcoming.

Reliable mass-loss rates are fundamental to our understanding of the late evolution of massive stars towards core collapse, either as supernovae (SNe) Ibc, or direct black hole formation (Eldridge & Vink 2006; Belczynski et al. 2010). Such mass-loss rates are especially relevant for understanding merging black holes in recent gravitational wave events such as GW 150914 (Abbott et al. 2016). These rates are also fundamental for correctly predicting the ionizing radiation of these hot stars with major consequences for interpreting He ii line emission at intermediate and high red shifts as Population iii stars (Cassata et al. 2013; Sobral et al. 2015) or very massive stars (VMS) at low metallicity (Gräfener & Vink 2015; Szécsi et al. 2015).

There has been much recent attention to detecting the elusive binary progenitors of Ibc supernovae. These objects should be stripped stars that in some sense resemble classical WR stars, but because of their weaker winds may either have very little line emission (Götberg et al. 2017) or no line emission at all. Addressing this question is one of the main aims of this paper.

Mass loss from massive stars is driven by radiative forces on spectral lines (Lucy & Solomon 1970; Castor et al. 1975; CAK). The CAK work developed the force multiplier formalism to treat all relevant ionic transitions. This enabled these authors to predict the wind mass-loss rate, , and terminal velocity, ν, simultaneously. Using a “global energy” Monte Carlo approach (Abbott & Lucy 1985) in which the velocity law was adopted – aided by empirical constraints – Vink et al. (2000) predicted mass-loss rates for galactic O-type stars that include multi-line scattering. Müller & Vink (2008) subsequently introduced a new parametrization of the line acceleration, expressing it as a function of radius instead of the velocity gradient (as in the CAK theory). The implementation of this new formalism improved the local consistency of Monte Carlo models that initially assumed a velocity law. Not only do we find reasonably good agreement with observation of wind terminal velocities, but because our method naturally accounts for the physics of multi-line scattering, it is also applicable to denser winds, such as those of WR stars (Vink et al. 2011a; Vink 2015).

This paper is organized as follows. In Sect. 2, we briefly describe the Monte Carlo modelling and parameter space considered for this study. The mass-loss predictions (Sect. 4) are followed by a discussion, before ending with a summary in Sect. 7.

2. Monte Carlo modelling

Helium star mass-loss rates are calculated with a Monte Carlo method, in which we follow the fate of a large number of photon packets (2 × 106) from below the stellar photosphere throughout the entire wind, up to the terminal velocity. The core of the method is linked to the total loss of radiative energy that is coupled to the momentum gain of the outflowing gas (Abbott & Lucy 1985; Vink et al. 1999). As the photon absorptions and scatterings depend on wind density and hence on the mass-loss rate , we are able to obtain a consistent model where the momentum of the gas equals the transferred momentum by photons. Here we shall employ our local dynamical approach (Müller & Vink 2008, 2014; Muijres et al. 2012), in which we predict , ν, and the wind structure parameter, β, simultaneously.

The Monte Carlo code employs the density and temperature stratifications from a model atmosphere calculation performed with isa-wind (de Koter et al. 1993). These model atmospheres are unified, implying continuity between the photosphere and wind region. The models describe the radiative transfer in spectral lines adopting an improved Sobolev approximation that includes the diffuse radiation field. The species that are explicitly calculated in non-LTE are H, He, C, N, O, S, and Si. The iron-group elements, which are crucial for the radiative driving, are treated in a generalized version of the modified nebular approximation (Schmutz 1991). However, we performed test calculations in which we also explicitly computed Fe (iron) in non-LTE. These tests revealed that differences with respect to the assumption of the modified nebular approximation for Fe were negligible, and we therefore decided to keep the treatment of Fe approximate.

The line list used for the calculations consists of over 105 of the strongest transitions of the elements H - Zn extracted from the line list constructed by Kurucz & Bell (1995). The wind was divided into 90 concentric shells, with many narrow shells in the subsonic region and wider shells in supersonic layers.

Other assumptions involve wind stationarity and sphericity. The latter seems to be a reasonable approximation, given that the bulk of WR stars are unpolarized (Vink & Harries 2017). However, there are notable exceptions, involving those that are surrounded by ejecta nebulae (Vink et al. 2011b). Furthermore, polarization variability due to small-scale structure or wind clumping has been encountered in luminous blue variables (Davies et al. 2005, 2007). It has been well established that small-scale clumping of the outflowing gas has a pronounced effect on the ionization structure of both O-star and WR atmospheres (e.g. Hillier 1991). This has lead to a downward revision of empirical mass-loss rates, by a factor of ~23 (Moffat & Robert 1994; Puls et al. 2008; Hamann et al. 2008; Vink & Gräfener 2012; Ramírez-Agudelo et al. 2017).

Moreover, small-scale clumping might have a quantitative effect on the radiative driving itself as well, thereby affecting predicted mass-loss rates. The subtle issues of clumping and porosity on the mass-loss rate predictions were studied by Muijres et al. (2011), who found that whilst the impact on can be large for certain specific clumping prescriptions, it was concluded overall that moderate clumping and porosity does not appear to change the mass-loss rate predictions substantially. In the present set of computations we therefore do not account for the effects of clumping. Nevertheless, it should be kept in mind that clumping might affect the predicted values.

3. Parameter space

We wish to predict mass-loss rates for He stars over a wide mass range, focusing on the 220 M range, but also going up to 60 M and down to 0.6 M into the regime of the subdwarf O (sdO) stars. The main focus of this study involves the stripped He stars of order 4 M, as this is the mass range of stripped stars that fall just below that of the classical WR stars (in the range 520 M; Nugis & Lamers 2000, which should represent absolute values fairly well in the WR range). The luminosities in our study follow from the mass-luminosity relation of Gräfener et al. (2011). This leads to similar values for the low-mass stripped He stars as the models by Götberg et al. (2017).

We do not express as a function of the Eddington parameter, as this would depend on the amount of ionized He versus H, whilst Vink & de Koter (2002) and Vink et al. (2011a) showed that should hardly depend of the amount of He. However, the effective Eddington factor remains of key importance for higher mass models in the classical WR range.

The effective temperature sets the ionization stratification in the atmosphere and determines which lines are most active in driving the wind. As a result, Teff affects the predicted mass-loss rate. Here, we fix Teff to 50 000 K, which is characteristic for the best observed example of a stripped He star HD 45166 (Groh et al. 2008). We do not express the mass-loss rates of classical WR stars as a function of Teff, as radii are highly uncertain because of envelope inflation and the possibly associated clumped nature of stellar envelopes (Gräfener et al. 2012).

Our grid was constructed to predict the mass-loss behaviour as a function of L (or M via the ML relation) and Z, which are scaled to the solar values (Anders & Grevesse 1989).

4. Results

Table 1

Mass-loss predictions for the He main sequence. Teff is kept constant at 50 000.

The main aim of the paper is to get a handle on the mass-loss rates for stripped He stars with masses just below those of classical WR stars, i.e. at about 4 M. We predict a rate of log  = −7.43 and a terminal wind velocity of 2871 km s-1 for our 4 M model. This mass-loss rate is an order of magnitude smaller than that assumed by Götberg et al. (2017) through extrapolation of empirical WR stars by Nugis & Lamers (2000).

Table 1 lists the rest of our mass-loss predictions. The initial stellar parameter columns are self-explanatory. The resulting predicted wind terminal velocities, mass-loss rates, and wind acceleration parameter β are given in Cols. (4)(6), respectively. The predicted mass-loss rates (Col. 5) are shown in Fig. 1. Various symbols are used to identify various Z ranges. For visual clarity, not all Z models feature here, but only a subset – with alternating Z values – from Table 1 were plotted.

thumbnail Fig. 1

Predicted mass-loss rates vs. log (L/L) for solar metallicity models (open squares), models of 10% Z (open triangles), and 1% Z (asterisks). The high-mass models of 60 M follow the bulk of the models very well; the very low-mass models with M = 0.6M (log (L/L) = 1.95) fall somewhat below the general behaviour.

Open with DEXTER

Figure 1 shows that increases with L and Z, as expected. Whilst the very high-mass model (with M = 60M; log (L/L) = 6.4) fits the relationship formed by the bulk of models in the 220 M range (log (L/L) = 3.49−5.69 range) very well, the very low-mass model (M = 0.6M; log (L/L) = 1.95) with log  = −11.12 falls below the general behaviour. It will therefore not be included in the mass loss versus L power-law fitting.

When we compare the M = 0.6 M model to the predictions of Vink & Cassisi (2002) for horizontal branch stars, we find log  = −10.40. However, the new dynamically consistent rate has ν a factor 2 higher, and for a constant wind momentum one would expect an value a factor ~2 lower, i.e. log  ≃ −10.70. This is within a factor of two from the current value. Our value for the 0.6 M model is a factor of ~5 larger than the SdO mass-loss predictions by Krtička et al. (2016).

We next turn our attention to the wind velocity structure parameter, β, which describes how rapidly the wind accelerates. The predicted values of β are shown in Col. (6) of the table. β shows a slight dependence on stellar mass, but mostly β values are of order unity, in accordance with the O-star models of Pauldrach et al. (1986), Müller & Vink (2008), and Muijres et al. (2012).

5. Mass-loss recipe for He stars

In order to determine the dependence of the mass-loss rate on L and Z simultaneously, we perform multiple linear regression, and we find (1)with a fitting error of σ = 0.11. The formula was derived for the Z range (Z/Z) = 1−10-2 and the mass range 260 M, but for reasons given earlier application of this formula is recommended only for lower-mass stripped stars, and not for classical WR stars.

The mass-loss versus mass relationship (as explored by e.g. Langer 1989) can be transformed using relevant mass-luminosity relationships (e.g. Gräfener et al. 2011).

6. Discussion

Comparing our mass-loss predictions against observed mass-loss rates is a non-trivial undertaking, as low-mass stripped He stars are as yet elusive.

Arguably the best example of a stripped star is HD 45166, which consists of a 4 M quasi-WR (qWR) star orbiting a B7V companion (Groh et al. 2008), which may have evolved from a 12 M primary (Götberg et al. 2017). HD 45166 has a surprisingly low (equatorial)) terminal wind velocity of 350 km s-1. For such a low ν, and employing our global approach (Vink et al. 2000), we would arrive at a mass-loss rate that is an order of magnitude larger than when we use a higher value, such as for the presumed polar wind component. The agreement between the observed mass-loss rate from HD 45166 from Groh et al. (2008) and the mass-loss rates given by extrapolation of the empirical Nugis & Lamers (2000) WR recipe is probably the main reason why Götberg et al. (2017) assumed that the Nugis & Lamers (2000) WR mass-loss rates provide rates in the expected range for low-mass He stars, whilst we find that they are actually overestimated by an order of magnitude.

Regarding the metallicity dependence, the exponent of 0.61 is slightly less than the predicted value of 0.85 in Vink & de Koter (2005). This is to be expected, as the earlier predictions were based on the global approach, and the difference may simply be attributed to the dependence of ν on Z (Vink et al. 2001).

Our -Z exponent is also smaller than that derived in the recent empirical study by Hainich et al. (2015), who suggest an -Z exponent larger than unity. Given that, firstly, the Hainich et al. empirical rates were derived for a data set showing a significant scatter when expressed as a function of L, which suggests the Hainich et al. (2015) study misses a relevant parameter and, secondly, that our exponent was derived over a far wider range in Z than the Hainich et al. data, we consider our lower exponent of 0.61 to be more likely close-to-correct.

Note that these -Z exponents are predominately determined by the Fe abundance, and distinct from the Z-dependence in the empirical Nugis & Lamers (2000) recipe. We do not consider the Nugis & Lamers (2000) formulation physically astute, but it should still provide reasonably accurate estimates within the actual classical WR range.

7. Summary

We presented mass-loss predictions from Monte Carlo radiative transfer models for stripped He stars by stellar winds or mass transfer through a companion.

  • The mass-loss rate is expressed as a function of L and Z through multiple linear regression.

  • We find log =13.3 + 1.36 log (L/L)+ 0.61 log (Z/Z), derived for (Z/Z)=1−10-2.

  • Our mass-loss rates for low-mass stripped stars are an order of magnitude lower than would be found by extrapolating empirical WR recipes.

This has major implications for the search for low-mass stripped He stars as the progenitors of SNe Ibc and ionizing sources in the Universe.


All Tables

Table 1

Mass-loss predictions for the He main sequence. Teff is kept constant at 50 000.

All Figures

thumbnail Fig. 1

Predicted mass-loss rates vs. log (L/L) for solar metallicity models (open squares), models of 10% Z (open triangles), and 1% Z (asterisks). The high-mass models of 60 M follow the bulk of the models very well; the very low-mass models with M = 0.6M (log (L/L) = 1.95) fall somewhat below the general behaviour.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.