Highlight
Free Access

 Issue A&A Volume 509, January 2010 A14 13 Interstellar and circumstellar matter https://doi.org/10.1051/0004-6361/200912084 12 January 2010
A&A 509, A14 (2010)

I. A grid of solar-metallicity wind models

L. Mattsson - R. Wahlin - S. Höfner

Dept. of Physics and Astronomy, Div. of Astronomy and Space Physics, Uppsala University, Box 515, 751 20 Uppsala, Sweden

Received 16 March 2009 / Accepted 27 August 2009

Abstract
Context. Knowing how the mass loss of carbon-rich AGB stars depends on stellar parameters is crucial for stellar evolution modelling, as well as for the understanding of when and how circumstellar structures emerge around these stars, e.g., dust shells and so-called detached shells of expelled gas.
Aims. The purpose of this paper is to explore the stellar parameter space using a numerical radiation hydrodynamic (RHD) model of carbon-star atmospheres, including a detailed description of dust formation and frequency-dependent radiative transfer, in order to determine how the mass loss of carbon stars changes with stellar parameters.
Methods. We have computed a grid of 900 numeric dynamic model atmospheres (DMAs) using a well-tested computer code. This grid of models covers most of the expected combinations of stellar parameters, which are made up of the stellar temperature, the stellar luminosity, the stellar mass, the abundance of condensible carbon, and the velocity amplitude of the pulsation.
Results. The resultant mass-loss rates and wind speeds are clearly affected by the choice of stellar temperature, mass, luminosity and the abundance of available carbon. In certain parts of the parameter space there is also an inevitable mass-loss threshold, below which a dust-driven wind is not possible. Contrary to some previous studies, we find a strong dependence on the abundance of free carbon, which turns out to be a critical parameter. Furthermore, we have found that the dust grains that form in the atmosphere may grow too large for the commonly used small particle approximation of the dust opacity to be strictly valid. This may have some bearing on the wind properties, although further study of this problem is needed before quantitative conclusions can be drawn.
Conclusions. The wind properties show relatively simple dependences on stellar parameters above the mass-loss threshold, while the threshold itself is of a more complicated nature. Hence, we chose not to derive any simplistic mass-loss formula, but rather provide a mass-loss prescription in the form of an easy-to-use FORTRAN routine (available at http://coolstars.astro.uu.se). Since this mass-loss routine is based on data coming from an essentially self-consistent model of mass loss, it may therefore serve as a better mass-loss prescription for stellar evolution calculations than empirical formulae. Furthermore, we conclude that there are still some issues that need to be investigated, such as the rôle of grain-sizes.

Key words: stars: AGB and post-AGB - stars: atmospheres - stars: carbon - circumstellar matter - stars: evolution - stars: mass-loss

1 Introduction

How the mass loss of carbon-rich AGB stars (carbon stars) depends on stellar parameters is not very well known. It is crucial, however, to have that kind of information in many contexts, such as for stellar evolution modelling. A variety of empirical (e.g., Groenewegen 1995,1993; Vassiliadis & Wood 1993) as well as theoretical (e.g., Arndt et al. 1997; Wachter et al. 2008,2002; Blöcker 1995) formulae with different parameterisations are available in the literature and the agreement between them is rather poor. These simple mass-loss prescriptions are nevertheless used in many models of stellar evolution, sometimes without enough consideration of their applicability in specific cases.

Determining the mass-loss rates observationally is complicated and riddled with many possible sources of uncertainty.Also, the fundamental stellar parameters, such as effective temperature and mass, can be hard to determine observationally with high precision, which makes it very difficult to empirically determine the mass-loss rate as a function of stellar parameters with sufficient accuracy. Empirical studies as such can only provide a limited amount of information, since one cannot easily disentangle the effects of individual stellar parameters. The reason for this is of course that stellar evolution and observed stellar parameters are connected, and we can only observe each star at one evolutionary stage. Observations can tell us about the mass loss for a given set of stellar parameters, but the exact dependences on these parameters can unfortunately not be determined easily due to selection effects. Existing empirical relations suffer from selection effects, with a very steep dependence of mass loss rates on stellar parameters.

The fundamental idea behind a so-called dust driven wind is that stellar photons, incident on dust grains, will lead to a transfer of momentum from the radiation to the atmospheric gas, which is dragged along by the dust grains. Models of this process have existed for quite some time and have become increasingly more complex. Wood (1979) and Bowen (1988) made dynamic atmosphere models introducing a parameterised opacity to describe the effects of dust formation in the circumstellar envelope. Bowen & Willson (1991) further elaborated on the implications of their models for the mass loss during AGB evolution, i.e., the development of a so-called superwind'' phase. Blöcker (1995) then derived a modified version of Reimers' law based on the models of Bowen (1988) to prescribe mass loss in his calculations of stellar evolution on the AGB. This mass-loss prescription is still widely used, although it is entirely based on the early work by Bowen (1988).

Fleischer et al. (1992) presented self-consistent dynamic models of carbon rich AGB stars that include a time-dependent description of dust formation. Their work revealed phenomena like a discrete spatial structure of the circumstellar dust shell and multiperiodicity. A picture has emerged with dust-driven winds as a complex phenomenon, where the onset of mass loss requires particular circumstances. For about a decade, time-dependent dynamical models for stellar winds (e.g., Höfner & Dorfi 1997; Winters et al. 2000) were based on rather crude descriptions of radiative transfer and simplified micro-physics, which led to unrealistic atmospheric density-temperature structures.

Höfner et al. (2003) presented dynamic model atmospheres which couple time-dependent dynamics and frequency-dependent radiative transfer, including a non-grey description of the dust component. Comparing the observable properties calculated from these state-of-the-art dynamical models (spectra and their variations over time) with mboxobservations has shown drastic improvements compared to the synthetic spectra based on previous (grey) generations of models (see, e.g., Gautschy-Loidl et al. 2004; Nowotny et al. 2005b,a). While further improvement may become possible in the future, a reasonable level of realism has now been reached in the modelling to allow application to stellar evolution. The model used in the present study is identical to that used in Mattsson et al. (2007a) and contains only minor modifications of the Höfner et al. (2003) model. It is probably the most advanced theoretical tool available at present for studying how dust driven winds of carbon stars depend on stellar parameters.

Detailed numerical modelling provides valuable insights into the mass loss problem, since it can be used to constrain the actual physics involved in wind formation. In this paper, we explore the stellar parameter space using a numerical radiation hydrodynamic (RHD) model of carbon-star atmospheres, including time-dependent dust formation and frequency-dependent radiative transfer, in order to determine how the mass-loss of carbon stars depends on stellar parameters. The results we give here reflect primarily the intense mass loss during the so-called thermal pulse phase, i.e., the very late stages of carbon-rich AGB evolution.

2 Theory and methods

2.1 RHD atmosphere models

The model includes frequency-dependent radiative transfer and non-equilibrium dust formation, i.e., we solve the coupled system of frequency-dependent radiation hydrodynamics (RHD) and time-dependent dust formation employing an implicit numerical method and an adaptive grid (cf. Höfner et al. 2003).

In our model, the stellar atmosphere and circumstellar envelope are described in terms of conservation laws for the gas, the dust and the radiation field, expressed by the following set of coupled, nonlinear partial differential equations (PDEs):

• The three equations describing conservation of mass, momentum and energy for the gas.
• The 0th and 1st moment equation of radiative transfer.
• Four moment equations of dust formation (cf. Sect. 2.5).
• The Poisson equation (self gravity).
To our system of nonlinear partial differential equations describing the physics of a C-star atmosphere we add a so called grid equation'' which determines the locations of the grid points according to accuracy considerations (Dorfi & Drury 1987) and an equation keeping track of the amount of condensible carbon, leaving us with a total of 12 partial differential equations (PDEs). This system of PDEs is then solved implicitly using a Newton-Raphson scheme. All equations are discretised in a volume-integrated conservation form on a staggered mesh. The spatial discretisation of the advection term is a monotonic second-order advection scheme (van Leer 1977). The same order of numerical precision is used for all PDEs. Details of the numerical method are discussed by Dorfi & Feuchtinger (1995) and in several previous papers about dust-driven wind models (cf. Höfner et al. 1995, and references therein).

2.2 General procedure

The RHD computations are started from hydrostatic dust-free initial models providing the atmospheric structure at t=0. When the dust equations are switched on, dust condensation starts and the resulting radiative acceleration creates an outward motion of the dust and the gas. In the first computational phase, the expansion of the atmospheric layers is followed by the grid to about (usually around to 1015 cm). At this distance, the location of the outer boundary is fixed, allowing outflow. The outflow model then evolves typically for more than 100 periods. To avoid a significant depletion of mass inside the computational domain, the model calculation is stopped after about 105 time steps. This is important since we cannot (due to the computational method) allow for material to flow over the inner boundary. Due to this fact, a large enough mass loss will eventually lead to a depletion of mass in the considered part of the atmosphere and circumstellar envelope.

2.3 Inner boundary: stellar pulsations

The dynamic atmosphere model does not include a physical model for the pulsation mechanism. The pulsation is modelled by a piston boundary condition'' (Bowen 1988) located at in order to be in the optically thick part of the atmosphere but outside the zone where the pulsations are excited. It should be noted that all results obtained in this paper are, to some extent, subject to this inner boundary condition, and as shown by Mattsson et al. (2008), the effect on the mass-loss rate is not negligible in the non-saturated wind regime where the degree of dust condensation is affected by the input of kinetic energy by pulsations.

We use a harmonic piston boundary given by the expressions

 (1)

and

 (2)

where is the period and is the velocity amplitude of the pulsations. To restrict the number of free parameters, we also employ an empirical period-luminosity relation (Feast et al. 1989) and keep the period tied to the luminosity for all models. The piston amplitude , or the kinetic-energy injection by pulsation (cf. the q-parameter in Mattsson et al. 2008)

 (3)

is essentially a free parameter in the model, but previous experience (e.g., Höfner & Dorfi 1997; Mattsson et al. 2007a) shows that only a limited range can be considered as reasonable. Nonetheless, we must consider the effects of picking different amplitudes, and we therefore consider the cases and 6.0 km s-1. The motion of the piston is accompanied by luminosity variations, since the radiative flux through the (moving) inner boundary is kept fixed, as described in previous papers, (e.g., Höfner & Dorfi 1997). In order to have luminosity variations that are consistent with observations, the velocity amplitudes given above turn out to be a good choice.

The inner boundary condition, i.e., the location of , was chosen so that all models have a similar kinetic-energy input due to pulsation for a given piston amplitude. This is critical in order to obtain a consistent picture of how the wind properties may vary with stellar parameters (Mattsson et al. 2008).

2.4 Outer boundary: mass loss

During the initial phase of the computation the expansion of the model is followed by the grid. When the outer-most point of the adaptive grid reaches some given distance from the inner boundary, the program automatically switches to a fixed outer boundary located at this point, allowing outflow over the outer boundary. There are no specific conditions posed on the location of the outer boundary, except that it must be at a large enough distance to ensure that the wind velocity at the outer boundary is close to the terminal velocity. At this point, the conditions

 (4)

are applied, where reflects the angular intensity distribution of the radiation field.

From the models described here we obtain the density and the wind velocity, both as functions of radius and time. The mass loss rate is then a function of time given by

 (5)

where and are the density and wind velocity at , i.e. the radius where the outer boundary is fixed (usually around 1015 cm, or typically 20-30 stellar radii, see above). Since practically all momentum transfer from radiation to matter is due to the interaction with dust, we expect the wind velocity to be correlated with the degree of dust condensation.

2.5 Dust formation

The models presented here include a time-dependent description of dust grain growth and evaporation using the moment method by Gail & Sedlmayr (1988) and Gauger et al. (1990). The dust component is described in terms of moments of the grain size distribution function weighted with a power j of the grain radius. The zeroth moment is the total number density of grains (simply the integral of the size distribution function over all grain sizes), while the third moment is proportional to the average volume of the grains. The moment equations are solved simultaneously with the RHD equations.

In order to calculate how much momentum is transferred from photons to dust grains, we need to know the frequency-dependent opacity of these grains. This can be expressed in terms of the extinction efficiency, which is the ratio of the extinction cross section to the geometrical cross section of the grains. In the small-particle limit (which is used here), the dust grain opacity becomes a simple function of the grain radius. The dependence of the opacity on wavelength and grain size can thus be separated into two independent factors, which greatly simplifies the calculations. However, as we shall see, this approximation may need to be relaxed, since the carbon grains tend to grow quite rapidly. The models in this paper are calculated using the refractive index data of Rouleau & Martin (1991) to obtain the dust extinction (see Höfner et al. 2003; Andersen et al. 2003, for further discussion). The intrinsic density of the grain material used in the model is set to  g cm-3, which matches the material in Rouleau & Martin (1991).

We assume that dust grains can be considered to be spherical, consisting of amorphous carbon only, since we are modelling carbon stars (see Andersen et al. 2003, for details about these assumptions). The nucleation, growth and evaporation of grains is assumed to proceed by reactions involving , , and . In this model of grain growth, a so called sticking coefficient (sometimes referred to as the reaction efficiency factor) is used that enters into the net growth rate of the dust grains. This parameter, , is not definitely known unless we know the exact sequence of chemical reactions responsible for dust formation. However, Gail & Sedlmayr (1988) argued that the sticking coefficient must be of the order of unity, mainly because it is expected that neutral radical reactions play a major role in the formation of carbon grains. We set , but note that a set of test models with give similar results.

Dust grains in a stellar atmosphere influence both its energy and momentum balance. For simplicity we assume complete momentum and position coupling of gas and dust, i.e. the momentum gained by the dust from the radiation field is directly transferred to the gas and there is no drift between dust and gas. However, this strong coupling between the dust and the gas phase is not obvious. In a previous attempt to relax this phase coupling approximation, Sandin & Höfner (2003, 2004) found that the effects of decoupling the phases might be quite significant. The most striking feature is that the dust formation may increase significantly, but this does not necessarily increase the predicted mass loss rates for a given set of stellar parameters. A coupled solution of the detailed frequency-dependent equations of RHD, including dust formation and drift would increase the computing time for each model to an extent that would make it very laborious to compute a grid of the kind presented here. Thus, we chose to give frequency-dependent radiative transfer priority over relaxing the phase coupling approximation, even if this approximation cannot hold for very low gas densities, i.e. at very large distances from the star or for wind models with very low mass loss rates. Furthermore, since transfer of internal energy between gas and dust is negligible compared to the interaction of each component with the radiative field (Gauger et al. 1990), we assume radiative equilibrium for the dust. This allows us to estimate the grain temperatures from the radiation temperature, as we know the dust opacities.

3 Definition of the grid

Since one purpose of the grid is (among other things) to provide information about how the mass loss rate is affected by changes in various stellar parameters as the star evolves on the AGB, it is wise to consider which values of the fundamental parameters ( ) never appear in evolutionary tracks. Systematic studies by, e.g., Blöcker (1995) show for example that for carbon stars at solar metallicity we should not expect average luminosities above or masses above on the AGB. The temperature of a carbon star lies approximately in the range 2000-3500 K according to both observations (e.g. Lambert et al. 1986; Bergeat & Chevallier 2005) and evolutionary models. Thus, the region in the plane (HR-diagram) that we are interested in is quite small for a star of a given mass and initial metallicity. It is important, however, to stress that during its evolution an AGB star may temporarily have parameter configurations that do not agree'' with typical observations. This is due to the fact that we are extremely unlikely to observe an AGB star as it undergoes, e.g., a thermal pulse, where the luminosity is very high and the temperature rather low (see Fig. 1 in Mattsson et al. 2007a). Therefore, defining the grid according to parameter combinations found in evolutionary tracks will require a somewhat larger parameter space than observations may suggest is necessary (see Table 1).

 Figure 1: Mass-loss rates as a function of wind speeds for models with  km s-1 and various values of . Open with DEXTER

Table 1:   Definition of the grid. denotes the grid step.

As shown by Mattsson et al. (2008), the kinetic-energy injection by pulsation can affect the wind properties significantly. Placing the inner boundary at a fixed optical or geometric depth led to different energy injections due to different density structures, resulting from different chemical compositions. This fact needs to be emphasised - theoretical work on AGB mass loss does not always include such energy control'' of the piston boundary, which may affect the results. Here we have tried to keep the equilibrium density/pressure at the inner boundary roughly the same for all models and used three different velocity amplitudes for the piston ( ).

 Figure 2: Mass-loss rate vs. wind speed for models with different , but the same carbon excess, . Open with DEXTER

Initially, parts of the grid were supposed to be computed using the solar composition by Grevesse & Sauval (1998) as well as that of Asplund et al. (2005), which were first thought to give slightly different results. However, after some scrutiny, it became clear that the differences were due to a poorly constrained inner boundary condition (see Mattsson et al. 2008). The grid of models presented here is computed using the Asplund et al. (2005) solar composition only.

4 Results and discussion

Table 2:   Input parameters (, , (C-O), P) and the resulting avergage mass loss rate, average wind speed and the mean degree of dust condensation at the outer boundary for a subset of models with and .

4.1 General trends

Wind velocities and mass loss rates are obviously correlated with some stellar parameters, such as luminosity, temperature and mass. The presence of trends with these parameters can be seen in observations, but in general the information is degenerate since the different stellar parameters are very much coupled throughout the evolution of a star. This is why, for example, one finds a simple mass-loss-period law (see, e.g., Groenewegen 1995,1993; Vassiliadis & Wood 1993) which of course does not contain all the information needed for a proper treatment of mass loss in, e.g., stellar evolution modelling. Empirical results do not tell us directly how the mass loss rate is affected by changes in individual stellar parameters.

 Figure 3: Wind speed ( left) and mass-loss rate ( right) as functions of the acceleration parameter for models with  km s-1. The dashed and dotted lines represent analytic models as described in Lamers & Cassinelli (1999) and Groenewegen et al. (1998), assuming different . Open with DEXTER

The mass-loss rate and the wind speed is, from a theoretical point of view, often related to the ratio of radiative to gravitational acceleration,

 (6)

where is the total flux-mean opacity and c is the speed of light. Here we define a similar quantity,

 (7)

where

 (8)

in which , i.e., the carbon excess (condensible carbon) and is the degree of condensation. is in principle proportional to and can be used in its place for most purposes. The reason we use instead of  is that only contains the physical quantities studied in this paper. In Fig. 3 we show the wind speed and the mass-loss rate as functions of for all models with . The wind speed is relatively well-correlated with , which is expected from simple analytical arguments (see Lamers & Cassinelli 1999, Chaps. 4 & 7). It is straight forward to show that, for a stationary polytropic wind,

 (9)

where is the escape speed, is the sound speed, is the condensation radius and  is the polytropic index in a polytropic equation of state. For the mass-loss rate, the correlation with  (or ) is weaker, which is due to variations of the gas-density profiles of the wind regions between different models. However, the mass-loss rate can (in the stationary case) be expressed analytically as (see Chap. 7 in Lamers & Cassinelli 1999, and the Appendix in Groenewegen et al. 1998)

 (10)

where is the flux-mean optical depth far out in the wind, and all other symbols are as previously defined. By analysing the time development of several wind models, it becomes obvious that the density profiles are not directly a consequence of the momentum transfer, as the wind speed profiles are, but of a combination of momentum transfer and pulsation dynamics. Thus, one should not expect a strong correlation between mass-loss rate and  for the non-stationary case.

 Figure 4: Wind speed ( left) and mass-loss rate ( right) as functions of the dust-to-gas ratio for models with  km s-1. Grey dots shows modelling results by Winters et al. (2000) and the dashed line marks the division between A- and B-type models according to Winters et al. (2000). Open with DEXTER

Previous efforts to model the winds of carbon stars are in many cases not directly comparable to the present results. For instance, stationary wind models with detailed dust formation should not be compared with DMAs, nor should models including drift but having grey radiative transfer be compared with models not including drift but featuring frequency-dependent radiative transfer instead (as in the present work). In comparison with existing theoretical work similar to the present (Wachter et al. 2002; Winters et al. 2000), our mass-loss rates show over-all trends that are comparable, although individual combinations of stellar parameters in some cases lead to quite different results compared to their models. More recent work by the Berlin group (Wachter et al. 2008) has been focusing on sub-solar metallicities and is not adequate for comparison. The wind speeds that we obtain for the most carbon rich models are higher than the highest observed (cf. Schöier et al. 2001), but consistent with the results by Winters (2000) and Wachter et al. (2002) for similar dust-to-gas ratios (see Fig. 4).

4.2 Wind velocity

The wind velocity depends strongly on the carbon excess , which appears as a steep slope in Fig. 6. This is expected, since the acceleration of the gas depends on the dust opacity, which in turn is proportional to the dust-to-gas ratio. The strong correlation between dust abundance and wind velocity (see Fig. 4) demonstrates how the wind speed is a direct consequence of momentum transfer from radiation to gas via dust grains. In addition, that the abrupt decline in the wind speeds at high temperatures (see Fig. 6) is mostly due to inefficient dust formation. At high the condensation radius is moved out to a distance too far away from the star for any significant amount of material to reach this radius. Consequently, almost no dust is formed and therefore no wind either.

4.3 Mass-loss rate

The classical Reimers (1975) formula for the mass loss rate of red giants,

 (11)

where we have used to eliminate , can be obtained by relating the kinetic energy of the outflow to the radiative energy flux and the gravitational potential energy. The Reimers law was derived strictly for red giants and cannot be applied to AGB stars without significant modifications (Kudritzki & Reimers 1978). Fitting a simple function to a set of wind models reveals a stronger dependence on luminosity and effective temperature for carbon stars (see, e.g., Arndt et al. 1997; Wachter et al. 2002). In Fig. 8 one may note a weak correlation between the mass-loss rates obtained from our detailed modelling and the Reimers formula, but since the Reimers formula (and commonly used modified versions of it, e.g., the formula suggested by Blöcker 1995) does not include any dependence on the free carbon abundance, it cannot easily be modified so that it would reproduce our results. Using a modified Reimers law to account for the mass-loss during the carbon star phase in stellar evolution models is therefore not to be recommended.

In general, the mass loss rate follows a slightly different pattern than the wind velocity. First of all there is an obvious anti-correlation between effective temperature and the mass loss rate (see Fig. 6), which partly can be interpreted as an effect of the correlation. But this anti-correlation is not directly due to the effect of surface gravity, , on the wind dynamics. With increasing , the star becomes more compact and the outer atmospheric layers become thinner, which means that there will be less material available in the dust-formation zone around the condensation radius. With lower density where the dust forms, the mass loss rate may drop, since there will be less gas available in the wind forming region.

The parameter configurations we have used here give mass-loss rates mostly within the interval  yr-1, which is in agreement with mass loss rates derived from observations (see, e.g., Schöier et al. 2001, and Fig. 1 in this paper), apart from the lack of models with mass loss rates below and the correlation between mass-loss rate and wind speed obtained from observations. This correlation could be an effect of the coupling between stellar parameters as the stars evolve on the AGB. But a detailed discussion about how to compare the present model grid with observations so that biases from sparse grid sampling and observational selection effects are reduced to a minimum, is beyond the scope of this paper. However, we may note that stars will stay much longer in evolutionary states that correspond to low mass-loss rates than those that correspond the superwind phase (typically the tip of the AGB) during which the mass-loss rate is very high and the evolution is rapid. Hence, one is more likely to find a star with a low mass loss rate. The low mass-loss rates may correspond to the transition regions where the mass-loss thresholds appear, and these transitions from no (or very inefficient) mass loss to strong, dust driven winds are perhaps not adequately resolved, the way the grid is defined. We have explored some of these transition regions in more detail by taking smaller steps in effective temperature and luminosity, which showed that the wind appears to be switched on'' at more or less a specific parameter value (see Fig. 9). Whether this steep threshold is a real property of dust driven winds or a model artefact is not clear at present, although mass-loss thresholds should exist for physical reasons (see Sect. 4.4). It is possible that a smooth transition is more realistic and that the observed stars with low mass-loss rates are sitting in these transition regions, which may correspond to very restricted parameter intervals. We will return to this important issue in a forthcoming paper.

The energy injection by pulsations, or more precisely, the kinetic energy density at the inner boundary of the model atmosphere, also plays an important rôle in the wind-formation process. A correlation between this quantity and the mass loss rate obtained for specific sets of stellar parameters was pointed out by Mattsson et al. (2008) and stems from the fact that the pulsations are needed to levitate the atmosphere such that dust formation can take place. The efficiency of dust formation is therefore more or less affected by the pulsation energy. Hence, for AGB stars, it is necessary to consider the kinetic energy from the pulsations as a part of the wind-driving mechanism, although it is the momentum transfer from radiation to dust that maintains the outflow. We have tried different piston-velocity amplitudes and, as expected, stronger pulsations seem to favour wind formation (see Mattsson et al. 2008, for further details), although such effects of are only important in the critical wind regime (near the mass-loss threshold). The most obvious effect of changing is that reducing the strength of the pulsations means that a smaller region in stellar parameter space will correspond to stable outflows and vice versa, i.e., changing can shift the locations of the mass-loss thresholds (see Fig. 2 in Mattsson et al. 2007b), but the typical mass-loss rate as well as the range of mass-loss rates in the wind-forming cases are not affected significantly (see Fig. 2).

Comparing the mass-loss rates we obtain with the single-scattering limit'' mass loss, i.e.,

 (12)

where c is the speed of light, we find that a fair number of models where the mass-loss rate is higher than , which can be seen in Fig. 5. The quantity is, however, not a proper upper limit for single scattering'' due to the assumptions made in deriving it (see, e.g., Lamers & Cassinelli 1999). The optical depth of the wind may very well be larger than unity around the flux peak in some cases, even if the total optical depth is not. Furthermore, it has been shown by Gail & Sedlmayr (1986) that heavily obscured objects must have mean optical depths that exceed unity, so the result in Fig. 5 is hardly surprising. Furthermore, the time-dependent luminosity L may not be represented properly by the luminosity  of the static start model used in Eq. (12) and in Fig. 5, since L is phase dependent in the dynamical models and wind acceleration may occur within a limited range of phases in a given model, thus changing the effective luminosity''.

 Figure 5: Modelled mass-loss rates vs. the single-scattering limit'' (see Eq. (12)). The dashed line marks the case of a one-to-one correspondance. Open with DEXTER

 Figure 6: Mean mass-loss rate ( upper panels), mean dust-loss rate ( middle panels) and wind speed ( lower panels) as functions of effective temperature and the abundance of condensible carbon for models with and  km s-1. Open with DEXTER

 Figure 7: Mean dust-to-gas ratio ( upper panels) and degree of dust condensation ( lower panels) as functions of effective temperature and the abundance of condensable carbon for models with and  km s-1. Open with DEXTER

4.3.1 Dust-to-gas ratios

The dust opacity scales with the relative dust abundance . Thus, we should expect the wind velocity to behave in a way similar to , which is actually the case if one compares Fig. 6 with Fig. 7. According to observations (e.g., Bergeat & Chevallier 2005; Schöier et al. 2001), the -ratio is nearly independent of . This fairly consistent with the weak -dependence that we obtain. However, the clearly increasing trend with the abundance of condensible carbon is not seen in observations (see, e.g., Bergeat & Chevallier 2005) although the existence of such a trend should be rather obvious since .

Because the stellar luminosity and mass is kept constant in Fig. 7, the plotted dust-to-gas ratio will also tell us how the acceleration ratio  varies with and the abundance of free carbon, . Below a certain , the dust-to-gas ratio depends mainly on , and we may thus conclude that is to a large extent set by for given a mass and luminosity.

4.3.2 Average grain size

We return now to the grain-size issue and whether the small-particle limit actually applies to the dusty atmospheres of carbon stars. The characteristic grain radius can be computed as

 (13)

where r0 is the monomer radius and K0, K1 are the zeroth and first moment of the grain size distribution function, respectively. Using this definition of grain radius, we obtain a typical grain radius of cm in the set of models with ,  km s-1, that produce winds. The small particle approximation is assumed to be applicable as long as the true radiation pressure efficiency factor does not deviate from that of the approximation by more than 10% at 1 m. Thus, the small-particle limit requires that  cm, i.e., not larger than a few times 10-6 cm (see Fig.
11), which indicates that the grains forming in our models are on average a little too large for to be well-approximated by the small-particle limit. The wind speed is anti-correlated with dust-grain size (see Fig. 10), which we interpret as a simple consequence of the fact that for a slow wind the grains stay longer in the dust-formation zone and thus they have time to grow bigger.

Since (which describes the net effect of momentum transfer by absorption and scattering for dust) will level-out to a constant above a certain grain radius, i.e., is only proportional to the grain size for small grains, our model (which is based on the small particle approximation of , irrespective of grain size) will tend overestimate the dust opacity for really large particles. However, in the transition region, where goes from being proportional to the grain size to being almost independent of the grain size, the small-particle approximation will instead underestimate the extinction by as much as a factor of five (again, see Fig. 11). Since there is a competition between nucleation and grain growth, we expect this underestimation to be important only in critical cases with slow winds. The majority of the grid models, however, show grain sizes that correspond to the transition region where is likely underestimated. Exactly how this may affect the mass-loss rates and other results in this paper is difficult to estimate and needs further investigation.

4.4 Mass-loss thresholds

Mattsson et al. (2007b) demonstrated the necessity of considering thresholds in mass-loss prescriptions for dust-driven winds. As we mentioned above, a mass-loss threshold appears, as one would expect, and we find that for and/or  K a dust driven does not seem to form regardless of the combination of other stellar parameters (see Figs. 6 and 7, as well as the additional online material). Since the surface gravity scales with and dust formation is a process highly sensitive to temperature, it is expected that the mass-loss rate is highly dependent on too. But, as shown above, there is a rather strong dependence on condensible carbon for both the wind velocity and the mass-loss rate, which is quite interesting in comparison with previous studies of this kind. Arndt et al. (1997) as well as Wachter et al. (2002) find a rather weak dependence on (which is their choice of parameter). Their results stand in quite sharp contrast to the results presented here. However, our findings here are, qualitatively speaking, hardly a new discovery. Höfner & Dorfi (1997) and Winters (2000) have already pointed out the C/O-dependence, especially in the critical wind regime, although this has not been widely recognised.

The mass-loss thresholds originate from the simple fact that there exists a critical acceleration ratio and that there consequently also exist critical values for , and for which . The principles behind this may easily be understood from basic physics, but none of the existing parametric formulae for dust-driven mass loss (e.g. Arndt et al. 1997; Wachter et al. 2002) contain any thresholds. These transition regions are perhaps not adequately covered due to the grid spacing we have chosen and the sudden onset of an outflow discussed earlier may partly be due to model assumptions. But, the existence of mass-loss thresholds as such is a fact, and should definitely be included in a mass-loss prescription for stellar evolution modelling, even if it is unlikely that one can find a sufficiently simple functional form that will cover all aspects of how mass loss depends on stellar parameters.

The threshold due to the amount of carbon excess implies that Galactic carbon stars with low observed C/O ratios should typically show no or very little mass loss. Some support for this idea comes from the work by Lambert et al. (1986), where about 50% of their sample of carbon stars showed C/O <1.1 (all stars in their sample have roughly solar oxygen abundances), and a fair number of those stars do not have any detected outflows. A puzzeling fact, however, is that many stars with C/O  do have rather strong winds (Ramstedt et al. 2006), observations which cannot be properly explained at present.

4.5 No simple formula!

 Figure 8: Comparison between our mass-loss models (with  km s-1) and the Reimers (1975) formula with . The dashed line marks a one-to-one relation and the dotted lines indicate a 0.5 dex spread. Open with DEXTER

 Figure 9: Wind speed ( left) and mass-loss rate ( right) as functions of the effective temperature near a threshold point in the grid. Note that linear extrapolation of the mass-loss rate may be adequate, since the steep threshold could be an effect of the approximations and limitations of the model. Open with DEXTER

Consider two stars with the same mass, luminosity and effective temperature. Increasing the amount of condensible carbon by a factor of 5 may increase the mass loss rate by almost an order of magnitude (see Fig. 6). This may not appear as a huge problem in a log-log-plot, as in Fig. 8 where the grid models are compared to a scaled Reimers (1975) law, even if is not a parameter in Reimers' formula. But note that there is a 0.5 dex scatter in Fig. 8 (dotted lines) which actually makes a considerable difference for, e.g., stellar evolution modelling, where changes during the carbon-star phase, causing a different mass-loss evolution and the occurrence of a superwind. A direct fit to the whole set of data, including both and as fitting parameters does not really fix the problem, since statistically preferred correlations (formulae obtained through, e.g., least-square fitting or -minimisation) only show a rather weak dependence on these two parameters (see Arndt et al. 1997; Wachter et al. 2002). On the other hand, if one keeps all other parameters fixed and then varies or , the mass-loss rate is definitely strongly affected and it therefore seems dangerous to conclude that these two parameters can be omitted in the prescription as Arndt et al. (1997) and Wachter et al. (2002) suggested.

For the reasons given above, we do not recommend that formulae derived from the data given in this paper are used to replace the actual grid. Instead, using the tables (on-line material) to create a look-up matrix for the mass loss (combined with multi-dimensional interpolation) in stellar evolution modelling will guarantee that detailed features (such as the increase of the mass loss rate as the amount of condensible carbon increases) will be properly included, while simplified formulae will not. As an example, Mattsson et al. (2007a) showed that, as the amount of condensible carbon changes rapidly (relatively speaking, often more rapidly than luminosity, effective temperature etc.) during the TP-AGB phase, simple prescriptions of mass loss may therefore not be correct, since that parameter (condensible carbon) is usually not included. Note, however, that the mass-loss prescription we provide here best describes the superwind of carbon stars (at the tip of the AGB) and winds corresponding to the early/mid AGB (where mass-loss rates are much lower) may not always be accurately covered (see Sect. 4.4 and below).

4.6 Uncertainties

Generally, numerical models are subject to various uncertainties due to simplifying assumptions, numerical limitations and the fact that some relevant aspects of the physical reality may not be included in the model. Several of the models presented here show intermittent, chaotic behaviour that may, or may not, reflect non-linearities in the physics of real dynamic stellar atmospheres. The existence of such non-linear features can affect the time averages of mass-loss rates, wind speeds and other wind properties that we have derived here. By necessity, the time series from which these averages are computed must be limited and the time series that we obtain for a specific set of stellar parameters over the 100-400 pulsation periods the models cover may not always be representative. For example, the deviating'' mass-loss rate seen in Fig. 6 for , T = 3000 K and is probably an effect of dynamical non-linearities, which turn out to be favourable for dust formation. But it is hard to know for certain without a very long time series, and the models cannot be evolved for more than a few hundred periods, since it would lead to significant mass depletion in the modelled region. For the more well-behaved'' cases, the relative error in the time average due to time variations of the mass-loss rate, wind speed etc., is typically no bigger than 10% (usually much less).

 Figure 10: Wind speed ( left) and mass-loss rate ( right) versus the characteristic grain radii for models with  km s-1 and . Note the absence of any correlatiopn with mass-loss rate. Open with DEXTER

Another source of uncertainty (or bias) is our choices of stellar parameters, i.e., the definition of the grid. As previously discussed, the threshold regions may not be adequately covered, which can affect the interpolations in the mass-loss prescription described above. As shown in Fig. 9, linear interpolation is probably sufficiently accurate in most cases. But the very steep thresholds we find when the grid spacing is decreased may not be reflecting the true nature of these critical regions. It cannot be excluded that one possible explanation for the steepness might be the resolution of the radiative transfer, which is relatively low (64 frequency points) in the present study. In critical cases, a higher resolution may slightly affect the momentum transfer efficiency (from radiation to dust). We have recently begun to study the effects of frequency resolution (i.e., computing models with high-resolution radiative transfer), but it is yet too early to draw general conclusions. It is also possible that the effects of dust grain sizes can affect the radiative transfer significantly in models near the thresholds. These effects are currently being investigated as well, and will be discussed in a forthcoming paper.

Finally, there are two, possibly important, ingredients that are not included the present DMAs: dynamical gas/dust decoupling (i.e., drift) and the effects of the size distribution of dust grains on the momentum transfer efficiency. Grain size appears to be a critical issue for M-type stars (Höfner 2008) and it may be of some importance also for the winds of carbon stars. It is known that drift does have an effect on the wind properties (see, e.g., Sandin & Höfner 2003, 2004), but it has not been studied for the frequency-dependent case together with a detailed description of dust condensation in time-dependent wind models.

5 Summary and conclusions

The present grid of wind models has shown that mass-loss rates and wind speeds are clearly affected by the choice of stellar temperature, mass, luminosity and the abundance of available carbon. Furthermore, there are also inevitable mass-loss thresholds below which a dust-driven wind cannot form, in certain parts of the parameter space. Contrary to some previous studies, we find a strong dependence on the abundance of free carbon, which turns out to be a critical parameter. In a stellar evolution context, we expect that this dependence on the carbon excess will naturally lead to the development of a superwind after a few dredge-up events. Hence, the AGB would be terminated as soon as the atmosphere has become sufficiently carbon rich, and the number of thermal pulses will be limited by a self-regulating mechanism: each thermal pulse increases the carbon excess, which in turn increases the mass loss.

 Figure 11: The absorption, scattering and extinction efficiency, and the resultant efficiency factor for radiative pressure as functions of grain radius for amorphous carbon dust evaluated at m. The vertical, thin, dashed line shows where the small particle approximation breaks down. Open with DEXTER

The set of models with  km s-1 show mass-loss rates that are typically a few times  yr-1 and only one model has a really low mass-loss rate . With  km s-1 and  km s-1 the region in the stellar parameter space where winds occur is somewhat smaller or larger, respectively. The typical mass-loss rates are otherwise similar to that of the  km s-1 subset. Models with a steady mass loss rate appear to be difficult to produce regardless of the combination of stellar parameters and piston amplitude. The difficulty of producing really low mass-loss rates may represent a shortcoming of the present wind model. It is yet unclear whether this is due to, e.g., the low-resolution radiative transfer of the models, or the phase coupling between gas and dust (i.e., absence of drift). A smaller number of test models with high resolution (440 frequency points instead of 64) has been computed, but no dramatic differences were found. Further study of this matter is needed, however.

Winds of carbon stars are known to be relatively fast, compared to other long-period variables. The wind speeds that we obtain lie in the range , where the fastest winds are found in the most carbon rich cases, i.e., carbon stars on the tip of the AGB. The lack of observed carbon stars with very high wind speeds ( ) may be a selection effect, since the chance of finding a carbon star at the very final stage of its evolution must be quite small. The very low wind speeds seen in a few models ( ), but not in observed carbon stars, is an interesting prediction of the models. Such slow outflows may exist, since a wind can, in theory, be maintained as long as the flow surpasses the escape speed at some point in the velocity profile.

The mean degree of dust condensation in the wind-producing models varies between and , while the typical value is around 0.3. In the light of the recent models of M-type stars by Höfner (2008), we derived the characteristic grain sizes for a subset of models. The characteristic grain radius in the the models that produce winds is  cm, which is a somewhat too large for the small-particle limit to be a good approximation. The dust-grain size shows an anti-correlation with the wind speed, while the mass-loss rate seems totally uncorrelated. Large grains tend to form in slow winds, since they stay longer in the dust-formation zone and therefore have extra time to grow. However, the anti-correlation seen in Fig. 10 may not be entirely correct due to the fact that we compute the dust extinction in the small-particle limit, which is not fully consistent with the grain sizes we obtain from the model of dust formation. The grain size is also anti-correlated with the abundance of condensible carbon, which is likely an effect of the competition between nucleation and grain growth.

Despite some shortcomings in reproducing low mass-loss rates, the present grid is still a significant step forward. This is the first systematic study of mass loss and other wind properties as functions of stellar parameters that is based on a dynamic atmosphere model including both frequency-dependent radiative transfer and detailed, time-dependent dust formation, and is probably the most advanced theoretical constraint on the super-wind of carbon stars at present. Our results have shown that strong, dust driven winds cannot develop under arbitrary physical circumstances, i.e., in certain regions of (stellar) parameter space, carbon stars do not experience any significant dust-driven mass loss.

The thresholds that appear when, e.g., the effective temperature becomes too high or the abundance of condensible carbon too low, must be taken into account when constructing a mass-loss prescription to be used with other types of stellar models, such as stellar evolution and nucleosynthesis. Simple mass-loss formulae cannot reproduce all aspects of these thresholds and we have therefore not derived any such formula. Instead, we provide an easy-to-use FORTRAN-code that finds the corresponding mass-loss rate from interpolated data cubes. The code will soon be made available on-line.

Acknowledgements
The authors wish to thank the referee, Martin Groenewegen, for his careful reading of the manuscript and constructive criticism that helped to improve the clarity of this paper. B. Gustafsson, K. Eriksson and A. Wachter are thanked for their reading of the manuscript draft and for giving valuable suggestions. This work was partly supported by the Swedish Research Council, Vetenskapsrådet.

Footnotes

... models
Tables 3-14 and FORTAN code are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/509/A14
... stationary
The time dependent case will, at least on average, resemble the stationary case in the sense that the wind speed is mainly due to the amount of dust opacity and how strong the gravitational field is. The main differences, both in snap-shots and average structures, occur in the dust-free pulsating atmosphere, influencing the conditions in the dust-formation zone, therefore also .
... loss
A FORTRAN routine for this purpose is available on-line at http://coolstars.astro.uu.se

All Tables

Table 1:   Definition of the grid. denotes the grid step.

Table 2:   Input parameters (, , (C-O), P) and the resulting avergage mass loss rate, average wind speed and the mean degree of dust condensation at the outer boundary for a subset of models with and .

All Figures

 Figure 1: Mass-loss rates as a function of wind speeds for models with  km s-1 and various values of . Open with DEXTER In the text

 Figure 2: Mass-loss rate vs. wind speed for models with different , but the same carbon excess, . Open with DEXTER In the text

 Figure 3: Wind speed ( left) and mass-loss rate ( right) as functions of the acceleration parameter for models with  km s-1. The dashed and dotted lines represent analytic models as described in Lamers & Cassinelli (1999) and Groenewegen et al. (1998), assuming different . Open with DEXTER In the text

 Figure 4: Wind speed ( left) and mass-loss rate ( right) as functions of the dust-to-gas ratio for models with  km s-1. Grey dots shows modelling results by Winters et al. (2000) and the dashed line marks the division between A- and B-type models according to Winters et al. (2000). Open with DEXTER In the text

 Figure 5: Modelled mass-loss rates vs. the single-scattering limit'' (see Eq. (12)). The dashed line marks the case of a one-to-one correspondance. Open with DEXTER In the text

 Figure 6: Mean mass-loss rate ( upper panels), mean dust-loss rate ( middle panels) and wind speed ( lower panels) as functions of effective temperature and the abundance of condensible carbon for models with and  km s-1. Open with DEXTER In the text

 Figure 7: Mean dust-to-gas ratio ( upper panels) and degree of dust condensation ( lower panels) as functions of effective temperature and the abundance of condensable carbon for models with and  km s-1. Open with DEXTER In the text

 Figure 8: Comparison between our mass-loss models (with  km s-1) and the Reimers (1975) formula with . The dashed line marks a one-to-one relation and the dotted lines indicate a 0.5 dex spread. Open with DEXTER In the text

 Figure 9: Wind speed ( left) and mass-loss rate ( right) as functions of the effective temperature near a threshold point in the grid. Note that linear extrapolation of the mass-loss rate may be adequate, since the steep threshold could be an effect of the approximations and limitations of the model. Open with DEXTER In the text

 Figure 10: Wind speed ( left) and mass-loss rate ( right) versus the characteristic grain radii for models with  km s-1 and . Note the absence of any correlatiopn with mass-loss rate. Open with DEXTER In the text

 Figure 11: The absorption, scattering and extinction efficiency, and the resultant efficiency factor for radiative pressure as functions of grain radius for amorphous carbon dust evaluated at m. The vertical, thin, dashed line shows where the small particle approximation breaks down. Open with DEXTER In the text