A&A 367, 513-520 (2001)
DOI: 10.1051/0004-6361:20000444
A. Frankowski - R. Tylenda
Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Rabianska 8, 87-100 Torun, Poland
Received 3 July 2000 / Accepted 4 December 2000
Abstract
We explore the influence of a companion on the wind from a giant
in a binary system. We have developed a simple method for estimating the
enhancement and intrinsic non-sphericity of the outflow in this case.
Assuming the Roche model (modified to account for the reduced effective
gravity of the giant), local stellar parameters (,
g) are
calculated. This, according to our approach, allows to derive
the local mass loss rate.
The influence of varying the model parameters on the results is studied.
A relation between the ratio of the stellar radius to the Roche lobe radius
and the overall wind enhancement is obtained both numerically and in an
analytical approximation. A similar relation is derived for the
equatorial-to-polar wind intensity contrast.
Key words: binaries: general - stars: mass-loss - stars: late-type - stars: AGB and post-AGB
Mass loss is known to be a key factor for stellar evolution in late stages. Matter escapes easily because of high luminosities and low surface gravities of evolved, cool giants. On the Red Giant Branch (RGB) and especially on the Asymptotic Giant Branch (AGB) the mass loss rate can be larger, sometimes by orders of magnitude, than the nuclear burning rate. Therefore it is mass loss which determines the AGB lifetime and which plays the key role in the stellar initial-final mass relation. Thanks to the stellar winds intermediate mass stars avoid supernova explosions. Instead, due to the copious mass loss terminating the AGB they give origin to planetary nebulae with their amazing shapes.
The stellar winds also affect the evolution of the interstellar medium by influencing its dynamics and spreading elements created in the stellar interiors. It is also believed that interstellar dust has its origin in cool atmospheres and dense winds blowing from the RGB and AGB stars.
The mass loss process becomes even more important, more interesting and more complex when it takes place in a binary system. The most dramatic phenomena occur when an RGB or AGB star happens to fill its Roche lobe. In this case, a dynamically unstable mass transfer to the companion takes place. This leads to the formation of a common envelope and a very rapid evolution of the entire system.
However, even far from the Roche lobe filling the extended atmosphere of an RGB or AGB star can be significantly disturbed by binarity. The gravitational attraction from the companion and the centrifugal force due to the rotation of the system are expected to increase the mass loss rate. Mass loss from the star also means mass and angular momentum loss from the system. Certain mass and angular momentum accretion on the companion is also possible. The binary parameters are therefore expected to evolve accordingly which, in turn, will affect the stellar mass loss. Thus the evolution of the RGB and AGB stars in binary systems deserves thorough studies as it can be significantly different from that of the single stars.
In this paper we try to investigate effects of binarity on the mass loss
rate and the wind structure from cool, evolved stars by analysing
intrinsic conditions under which the mass loss process occurs. Generally
speaking, mass loss is driven by physical conditions in the stellar
atmosphere. Physics involved in mass loss from cool, extended atmospheres
is very complex and difficult for theoretical modelling although a
significant progress has recently been made in its understanding (e.g. Bowen
& Willson 1991; Arndt et al. 1997).
In the classical theory of (plane-parallel) stellar atmospheres
the physical structure of the
atmosphere, for a given chemical compositon, is determined by two parameters:
effective temperature,
,
and gravity, g.
In the case of the RGB and AGB stars the situation is
certainly more complex. Geometrical extension of the atmosphere and phenomena
generated in the stellar interior, e.g. pulsations and convection,
influence
the structure of the atmosphere. The ratio of the photospheric
scale height,
(
), to the stellar radius, R,
can be considered as a measure of
the geometrical extension. If we define
as
then
it seems reasonable to
adopt, at least to a first approximation, that
,
g and
are principal parameters of the RGB and AGB atmospheres and
their winds. This notion finds its support in results of modelling winds from the red
giants. Bowen & Willson (1991) have found that the mass loss rate
from the AGB stars depends first of all on the density scale height in
the atmosphere. Arndt et al. (1997) obtained formulae in which the
AGB wind characteristics are functions of the effective temperature, stellar
mass, and luminosity.
In the famous wind formula for red giants of
Reimers (1975), the mass loss rate depends on the stellar
mass, radius and luminosity.
However, the formulae of Arndt et al. and Reimers can be reduced to
dependencies on local atmospheric parameters,
i.e.
,
g and
(see next section).
This becomes important in the case of
a binary system, where the mass losing star can be noticeably distorted
by tidal forces. Differences in the local conditions over the stellar
surface may then lead to different
local intensities of the outflow. Thus the wind will show an intrinsic
directivity and may be globally enhanced.
Next section presents details of our approach. Results of numerical modelling of the wind structure are given and discussed in Sect. 3. An analytical approximation, which helps to understand and interpret the results, is developed in Sect. 4. General discussion and summary are given in Sect. 5. Preliminary results of this study have been presented in Frankowski & Tylenda (2000).
Principal assumptions made in our model are the following: 1. the binary orbit is
circular, the giant corotates with the orbital motion, and the gravitational
fields of both components can be
treated as originating from point masses -
hence the Roche model for the effective potential applies; 2. the stellar
surface is defined by the Roche equipotential surface; 3. the local mass
loss rate per unit area, ,
is a function of the
local stellar parameters, namely the effective temperature,
,
gravity, g, and the ratio of the photospheric scale height
to the stellar radius,
.
Our formulae for
have been derived from relations
obtained by other authors, either from theoretical modelling or observations.
These relations usually link the global mass loss rate,
,
to the global stellar parameters, i.e. stellar mass,
M, radius, R, and luminosity, L.
They can be transformed to dependencies on the
local atmospheric parameters using standard definitions and relations, i.e.
,
,
g=GM/R2 and
(where
is the
Stefan-Boltzmann constant, G is the gravitational constant, k is the
Boltzmann constant,
is the mean molecular weight, and
is
the hydrogen atom mass). In this way the Reimers' (1975) law can be
written in the form
Now, in general form, we can write:
We have introduced a free parameter, ,
describing the effectiveness
of the Newtonian gravitational force. This has been done to account for
effects of pulsation, radiation pressure, and
convective motions acting to diminish the gravity on the giant's surface.
Thus we assume that the
gravitational acceleration from the giant alone is:
![]() |
(8) |
With the assumptions described above, we have calculated sequences of models,
representing stars filling an
increasing fraction of their Roche lobes for various mass ratios, q.
The filled fraction of the giant's Roche lobe can be represented by the
ratio of the giant's volume radius, R, to the Roche lobe volume radius,
:
.
(The volume radius is defined as
a radius of a sphere of the same volume as the volume of the considered object.)
For each point of the stellar surface we have computed radius, gravity
and effective temperature (relative to the spherical case). This allowed us
to evaluate the local mass loss rate.
Integrating
over the
whole surface leads to the total mass loss rate, which can be compared
with the mass loss rate for a single
(spherical) star of the same intrinsic parameters (L, M, R).
Further, an intensity contrast between the equatorial and polar mass loss rates
can be calculated as a ratio of the local mass loss rate averaged along
the stellar equator,
,
to the polar mass loss rate,
(the latter is usually close to the minimum value of
).
The most important results of our study are presented in Figs. 1-4.
Figures 1a-d show the local mass loss rates
represented by grayscale. For each figure black
denotes maximum and white - minimum in the mass loss rate.
The mass loss rate at the dotted line is equal to that of a single,
undistorted (i.e. spherical) star. It marks the border between the polar
regions where the mass loss rate is lower, and the "equatorial strip'' where
it is higher than for the single star. The calculations have been done using
the local mass loss rate according to
the Arndt et al. prescription (Eq. (4)).
The gravity effectiveness factor
has been adopted.
The effective temperature surface distribution follows from
the assumption of uniform flux distribution over the solid angle (case (B)
described in Sect. 2).
![]() |
Figure 1:
a-d) The local mass loss rate per unit area across
the surface of a giant in a binary system.
Dark and light regions have, respectively, high and low mass loss rate per
unit area. At the dotted line the local mass loss rate is equal to that of
a single star. The mass loss formula of Arndt et al. has been applied.
The mass ratio is q = 0.5. The gravity effectiveness is
![]() ![]() ![]() ![]() |
Open with DEXTER |
Figures 1a-c have been obtained for the mass ratio q=0.5 and
the ratio of the giant's volume radius to the critical Roche surface volume radius
= 0.1, 0.5, 1.0, respectively. The stars are viewed from the
orbital plane, with the hemisphere facing the companion being closer to the
observer. Figure 1d presents for comparison the q = 0.5,
case (the same as in Fig. 1b) viewed from the opposite hemisphere.
![]() |
Figure 2:
a and b) The dependence of the total mass loss rate enhancement,
![]() ![]() ![]() |
Open with DEXTER |
Figure 2 plots the enhancement of the total mass loss rate caused by
the presence of the
binary companion,
,
against the
ratio.
Note the logarithmic scale on both axes. Points represent the numerical results,
lines - the analytical relation derived in Sect. 4, i.e.
![]() |
Figure 3:
a and b)The dependence of the equatorial-to-polar wind intensity contrast,
![]() ![]() |
Open with DEXTER |
Figure 3 presents the wind intensity contrast as a function of the
ratio.
For clarity of the functional dependence,
instead of
is plotted.
Again, note the logarithmic scale on both axes. As in Fig. 2, points
represent the numerical results, lines - the analytical relation derived in
Sect. 4, i.e.
![]() |
Figure 4:
a and b) The dependence of the equatorial-to-polar wind intensity contrast,
![]() ![]() ![]() |
Open with DEXTER |
Figure 4 shows the influence of the temperature distribution
law on the
wind intensity contrast. Results of case (A) with
,
0.08,
0.25 and of case (B) are displayed for q = 1.0 and
.
As in Figs. 2 and 3, points represent
the numerical results, lines - the derived analytical relation.
Figure 4a is for the Arndt et al. mass loss formula,
Fig. 4b - for the Reimers' law. In the case of the Reimers'
formula, the model with
is not shown, because it has a reversed intensity
contrast (the wind is more intense on the poles than at the equator - see
Sect. 5).
We have also made analytical considerations in order to obtain
approximate expressions for
and for
.
As usual,
analytical formulae help to understand and interpret the numerical results.
Our analysis has been done in the limit of small R/A
(which is equivalent to small
but makes the formulae more clear).
Let us define a corotating spherical coordinate system (
), with the origin at the center of mass of the giant, with the
angle measured from the axis parallel to the binary rotation axis,
and with
measured in the orbital plane from the axis pointing to the
companion.
The stellar surface is then defined as a function
.
In order to use Eq. (5) for the mass loss rate, we need expressions for
,
,
and
in the considered limit.
Actually, we need only formulae for r and g, because when R/A is small,
is determined by
(either in (A) or in (B)
case).
For a given point with coordinates (
)
in the corotating
frame, the potential function,
,
is given by:
![]() |
(9) |
![]() |
(14) |
The equation of surface of the giant can be written as:
![]() |
(17) |
![]() |
(20) |
![]() |
(21) |
![]() |
(23) |
Now (again making use of the fact that for small R/A we can regard the
stellar surface as spherical) an integration of Eq. (26) over
the sphere of radius Ryields an expression for the global mass loss rate:
We have constructed a model of a tidally enhanced wind from a
giant in a binary system.
The model has five parameters. Two of them, the
mass ratio q and the Roche lobe filling factor
,
characterize
the binary configuration. The other three - i.e. the gravity effectiveness
factor
,
the temperature distribution parameter
,
and the
mass loss law - describe the properties of the mass losing giant.
Figures 1-4 and the equations derived in
Sect. 4 show
the influence of these parameters quantitatively. A more general description and
interpretation of the results is given below.
At the beginning we note, that
one of the most prominent features of our model (that can be seen from
Sect. 4) is that when the Roche lobe filling factor
is small, the departures of all the local mass loss related quantities from
their values in the spherical case scale with the third power of R/A (and
therefore also with
).
Other model parameters can be thought of as constituting "scaling factors''
for this general dependence.
We shall return to this observation in more detail later on.
As can be seen from Eq. (27), the global mass loss rate
enhancement depends only on one of three exponents in the mass loss
formula (Eq. (5)), namely on the gravity exponent, c.
The reason for this is evident from Eq. (26), as the
term is averaged to 1 when integrating over the
stellar surface. The angular dependence of the stellar radius contributes
only to the wind intensity contrast.
The sign of c decides whether the overall mass loss from the giant in a binary system
will be greater (negative c) or smaller (positive c) than that of the
single giant. As has been noted in Sect. 2 (see comments on
Eq. (6)), as gravity counteracts matter ejection, it seems natural to
assume c<0. Both mass loss formulae considered in this work
(Eqs. (3) and (4)) indeed have c<0, and so the overall mass
loss rate is enhanced compared to the single, non-rotating star case.
This is also true for the AGB wind formulae of Blöcker (1995)
and for the RGB wind formulae reviewed by Catelan (2000 -
with one exception, when a dependence of mass loss on radius only is assumed).
The influence of the mass ratio q is clear when R/A is used as a
variable (Eq. (27)) - the higher the mass of the companion, the
stronger is the wind enhanced. However, things may look a bit confusing when
is used (see Eqs. (29) and Fig. 2).
The results shown in Fig. 2 seem to be - at first glance - counterintuitive:
the lines displaying the wind enhancement for q=0.1 lie always above the lines for
q=1.0. But a greater q means a stronger companion's influence, so why should
the global wind enhancement decrease with increasing q? This is just a
result of using
as the abscissa. The ratio of the Roche lobe
radius to the separation,
,
decreases with
increasing q, and this causes the whole expression
also to decrease. The increasing (1+q) term does not prevail until
.
Thus - in the range of q considered here -
for a given
,
the overall wind
enhancement is in fact greater for a smaller q.
The intrinsic structure of the wind in the presence of a companion, as described by our model, can be outlined as follows. For a star well inside its Roche lobe the local mass loss rate per unit area is governed by Eq. (26). This means that maxima and minima of the local mass loss rate coincide with extrema of the stellar radius. These are easy to find, as the stellar surface can be approximated by an triaxial ellipsoid (Eqs. (18) and (19)). The principal axes of the ellipsoid are defined by: the direction to the companion (semi-axis a1), the direction in the orbital plane that is perpendicular to the latter (semi-axis a2), and the direction perpendicular to the orbital plane (semi-axis a3). Equation (19) implies also that a1>a2>a3 for any q>0.
Magnitudes of the local mass loss rates corresponding to the semi-axes
follow this order when the expression
is
positive, but have a reversed order when it is negative
(see Eqs. (26) and (28)). In the first case
the global maxima of the mass loss rate occur in the directions to/from the
companion, and the global minima occur at the poles. There are also two saddle
points, corresponding to the semi-axis a2. These are the minima in the
coordinate, but maxima in
.
The values of
at the saddle points
are higher than (or equal to) the value for the spherical case,
,
when
,
otherwise they are
lower. When comparing
at the saddle points with the average over the whole
stellar surface,
,
the analogous condition is
.
Thus not always there is a well-defined "equatorial strip'' of enhanced mass loss. However,
for E>0, and provided that
(see above),
the average equatorial mass loss rate,
,
is enhanced
both compared to
,
as well as compared to
,
independently of the mass ratio q. (As a matter of fact, the
requirement on c is even less severe, namely
,
and is needed
only when comparing to
.)
In the case of negative E the mass loss is most intense at
the poles and weakest in the directions to/from the companion. Also the
situation at the saddle points is changed. The condition for
is still
,
but the
one for
is reversed, i.e.
.
Further, if the mass loss from these saddle points is in fact enhanced, then
now a "meridional strip'' of enhanced mass loss is created.
If the opposite is true, there is an equatorial strip of reduced mass
loss.
Finally, if E=0 the local mass loss rate per unit area is constant over the stellar surface.
For a given mass loss law (a given set of a,b,c) it is ,
the
temperature distribution parameter, that determines the sign of E.
In the case of the Reimers' mass loss formula E is positive when
is smaller than 0.1875. For the Arndt et al. prescription
all values of
greater than
result in positive E.
Since
seems unrealistic, all the models with the Arndt et al. mass loss
law have an equatorially enhanced mass loss. However, the case of
Reimers' formula with
falls in the region of negative E.
The wind intensity contrast,
,
is reversed - the mass
loss rate on the poles becomes greater than at the equator.
For this reason Fig. 4b does not show the case of
.
Nevertheless, we may conclude that for the values of
reasonable for
stars with deep convective envelopes, our model always gives the equatorial wind
enhancement, as expected for cool giants.
This holds also for the aforementioned wind formulae of Blöcker
(1995) and Catelan (2000).
When the star approaches to fill its Roche lobe, an asymmetry between the sub-companion and off-companion hemispheres develops. Mass loss towards the companion becomes more intense than in the opposite direction. The asymmetry grows until the gravity at the sub-companion point drops to zero (i.e. when the inner Lagrangian point L1 is reached), and the local mass loss rate per unit area at this point becomes infinite.
The above picture remains qualitatively the same
independently of the value of .
The effect of
is essentially to cause the "effective Roche lobe''
to shrink. As a result, the giant fills its critical equipotential surface
even while it fits comfortably inside its clasically determined Roche lobe.
The mass loss enhancement/reduction phenomena described
above are more pronounced and the deviations from the ellipsoidal symmetry
occur for smaller values of
.
The quantitative
differences are evident from Eq. (29) and
Figs. 2-3. For a given
both the global mass loss enhancement and the
equatorial-to-polar intensity contrast enhancement scale as
.
The use of
in our modelling can be justified in different ways.
As mentioned in Sect. 2 any kind of motions, like pulsations,
waves, turbulence, which are expected
to be important in the atmospheres of giants will effectively act as a
factor reducing the gravity. In the case of AGB stars there are additional
reasons why
significantly smaller than 1 can be expected.
The AGB stars have dusty winds and therefore
it seems reasonable to consider, that the intensity of the wind
should strongly depend on the conditions in the dust condensation region,
which lies well above the photosphere. Replacing in our
equations R by the condensation radius,
,
and assuming, say,
,
pushes the maximal possible stellar radius back to
and
increases 8 times the analytical values for the enhancement/reduction
of mass loss related quantities.
Therefore it is roughly equivalent to setting
.
In reality, however, it seems more likely that it is a certain combination of
conditions in the stellar atmosphere and in the dust condensation region,
that governs the mass loss.
The above considerations are consistent with the results of a study of S stars by Jorissen & Knapp (1998). These authors do not find any significant mass loss enhancement for the so called "extrinsic S stars'' which are believed to be RGB giants in binary systems. However, they find that among the "intrinsic S stars'', which are AGB giants, those which are in binaries tend to have the largest mass loss rates.
The assumption made by Bowen (1988) in his models of dusty
AGB winds, that the ratio of the acceleration due to radiation pressure on
dust to that due to gravity is 0.95, translates into setting
.
The idea of a reduced effective gravity of a mass losing star has also been
remarked by Jorissen et al. (1998) in connection
with short-period Ba stars.
An interesting result has been obtained by Mürset & Schmid (1999)
from their analysis of symbiotic systems. These authors
conclude that the radius of the cool giant in the symbiotic binaries
is limited to
,
where
is the distance from the center of
the giant to the inner Lagrangian point.
This translates into
-
,
depending on q. Mürset & Schmid suggest that when the giant exceeds this limit, a dynamical
mass transfer ending the symbiotic phenomenon occurs.
All the systems involved are s-type symbiotics so they are expected
to contain an RGB star. Thus Reimers' formula for the mass loss rate can be applied.
The intensity of the outflow from the sub-companion point in our models
with the Reimers' wind exhibits a steep rise
at
when
-0.35is adopted, depending on q.
In theoretical modelling of the evolution of binary systems one has to cope
with the fact that
the stellar mass loss rate can be affected by the presence
of the companion. Tout & Eggleton (1988) have suggested
on a heuristic basis that the wind from a star in a binary system is
enhanced by means of tidal forces by a factor of
,
where B is a free parameter. These authors have proposed B=104.
Since then, this prescription has often been invoked in various studies
of the binary evolution (e.g. Miko
ajewska et al. 1989;
Tout & Hall 1991; Soker 1994, 1998;
Han et al. 1995; Skopal et al. 1997).
However, attempts of establishing a reliable estimate of Bhave not brought a conclusive result. Han et al. (1995) have
preferred B=500, whereas in a later work Han (1998) has adopted B=1000.
Soker et al. (1998) have written the enhancement
factor in a more general form, i.e.
with
and
being free parameters. From their modelling of
off-center planetary nebula nuclei they argue for
and
.
The numerical and analytical results presented in the present work suggest
a different form of the enhancement factor.
As has already been mentioned, when the Roche lobe filling factor
is
small, the enhancement/reduction of all the local mass loss related quantities
scales as
.
Also the formula for the global wind
enhancement is of the form
(see Eq. (29)). The value of the scaling factor B' is again uncertain,
but the
pattern is common for all the
mass loss related quantities. In fact, it stems from the Roche model geometry
and thus should not depend on details of any particular account to the
mass loss description. In this paper
we have concentrated on the winds from red giants of RGB
and AGB, but as we allow for various
exponents in the
vs. g law and for various wind formulae, our approach is not limited to
the late type stars only.
Our general results should apply to all stars losing mass
in the form of wind, irrespectively of their spectral type.
It is worth noting, that the applicability of the above pattern also extends
onto single, rotating stars.
Although single red giants are generally expected to be slow rotators, and
therefore the effects of rotation on the wind should be negligible, recently
Heger & Langer (1998) and García-Segura et al. (1999)
conjectured mechanisms that allow single giants at the tip of the AGB to
rotate at considerable rates.
Our results can be applied also in this case.
Indeed, the expression for the global
mass loss enhancement (Eq. (27)) with q = 0 becomes identical
to that for a single star rotating with the angular velocity
.
For a single star, however, A has a different
physical meaning - it is a break-up radius.
If one defines the break-up velocity as
then Eq. (27) can be rewritten in a more
convenient form for single, rotating stars, i.e.
![]() |
(30) |
![]() |
(31) |
In this paper we have dealt with the intrinsic structure of the stellar wind, caused by the local differences in physical conditions at the origin of the outflow. Further fate of the ejected matter is beyond the scope of our work. Certainly, the giant's orbital motion and the still present gravitational influence of the companion can change the equatorial-to-polar contrast and complicate the structure of the outflow. Scenarios for such processes have been proposed by Soker (1994, 1997) and - along with a detailed numerical modelling - by Mastrodemos & Morris (1998, 1999). The equatorial compression of the already ejected matter is possible also in winds from single, rotating stars, as shown by Bjorkman & Cassinelli (1993) and Asida & Tuchman (1995). The picture can be further enriched by introducing magnetic effects (García-Segura et al. 1999). Our results could be used as inner boundary conditions for considerations of these types.
Acknowledgements
This work has been supported from the grant No. 2.P03D.020.17 of the Polish State Committee for Scientific Research.