A&A 395, 899-905 (2002)
Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany
Received 3 April 2002 / Accepted 13 September 2002
The structure equations of contact binaries are discussed. An equation for the transfer of mass between the components is derived. Serious uncertainties concern only the transfer of energy. They are expressed as unknown functions. An evolution code for contact binaries is presented. The restrictions imposed by a spherically averaged treatment of the components are discussed. Apart from these restrictions the code can be adapted to any choice of the functions describing the energy transfer. A simple tentative choice for these functions is proposed.
Key words: stars: binaries: close
Lucy (1968) had the key idea of a common convective envelope in which the entropy is constant and energy is transferred from the primary to the secondary. Thermal equilibrium turned out to be usually not possible. Lucy (1968) found age zero equilibrium models (violating the period-colour relation) only in a limited mass range and for extreme Population I composition. Moss & Whelan (1970) and Whelan (1972) encountered similar difficulties. Indeed, in unevolved systems with typical masses and typical composition thermal equilibrium was found to be in conflict with the equal entropy condition (Kähler 1995). If observed systems are assumed to be evolved, thermal equilibrium, when possible, turned out to be unstable (Kähler 1997a).
All these results are based on the assumption that the energy transfer between the components occurs in the adiabatic part of the common convective envelope (adiabatic transfer). The modification due to superadiabatic transfer is small (Whelan 1972; Biermann & Thomas 1973). Extreme superadiabatic transfer (Moss & Whelan 1973) gives more freedom, but Hazlehurst (1974) showed that this case must be excluded since the heat capacity of the secondary's subphotospheric layers is too small.
Abandoning the restriction to thermal equilibrium, Lucy (1976), Flannery (1976), and Robertson & Eggleton (1977) obtained contact binary solutions evolving in thermal cycles about a state of marginal contact. The solutions are in the required mass range and in agreement with the period-colour relation, but they have bad light curves since for a considerable part of the time the temperature difference between the components is too large. In a series of papers Rahunen and Vilhu encountered the same difficulty, the so-called light curve paradox, also for systems evolving in cycles without loss of contact. They concluded that the TRO (thermal relaxation oscillation) theory is not able to give a correct description of contact binaries (Rahunen 1983). The cycles investigated by Kähler et al. (1986) have indeed again bad light curves. In a recent paper Hazlehurst (2001) confirms the conclusion that the present theory of cyclic contact binaries is not in a position to resolve the light curve paradox.
As long as a proper hydrodynamic treatment of contact binaries is not possible, assumptions are necessary in any discussion. The results on thermal cycles obtained so far are based on few and apparently weak assumptions. The deficiencies of the resulting solutions have led to a number of investigations based on stronger assumptions in order to find contact binary models in thermal equilibrium which are compatible with the observations and sound from a theoretical viewpoint. All these attempts failed. In particular, the contact discontinuity hypothesis (Shu et al. 1976) violates the second law of thermodynamics (Hazlehurst & Refsdal 1978; Papaloizou & Pringle 1979), the assumption of an extended turbulent envelope in the primary (Kähler 1989; Kähler & Fehlberg 1991) is also in conflict with the second law (Hazlehurst 1993), and the assumption of sufficiently rapid internal mass motions in the primary's envelope (Kähler 1995) is in conflict with the observations (Kähler 1997a). For additional arguments against the possibility of thermal equlibrium see Hazlehurst (1993).
Collecting these results, thermal equilibrium appears to be as impossible as thermal disequilibrium. Nevertheless contact binaries are observed. We are thus led to the conclusion that the assumptions made so far in the theoretical discussion are not weak but too restrictive. This casts a doubt upon results which have been regarded as reliable for decades. In particular, we expect that if the restrictions are removed the light curve paradox for thermal cycles can be resolved and/or thermal equilibrium can be shown to be possible. The purpose of the present paper is to collect the structure equations of contact binaries on the basis of truly weak assumptions, to discuss the remaining uncertainties, and to present an evolution code for contact binaries which is flexible enough to be an appropriate tool when investigating the theoretical problems and trying to resolve the conflict between theory and observations.
We adopt a self-consistent spherically averaged treatment of the
components (Kähler 1986, 1997b) and assume uniform rotation
and conservation of total mass M=M1+M2 and angular momentum J.
Each component is assumed to be in hydrostatic equation, but the Roche equipotential condition is not assumed to be satisfied.
The system may be either in contact or semi-detached or detached.
|Figure 1: The function describing the area of the neck in systems in shallow contact (see text).|
|Open with DEXTER|
Since mass M and angular momentum J are conserved, the global balance of energy (Kähler 1997b) requires that
The area of the neck can be used can be used to derive constraints
for the exchange of mass and energy. We shall need the function Bidefined by
If the system is semi-detached, matter in the neck is streaming from the loser to the gainer, approximately with the velocity of sound, and through an area which is smaller than
(cf. Fig. 1.6 of Pringle 1985). Accordingly we have
We need an equation for the rate of mass exchange which depends continuously on the structure of the system and which is physically plausible in all cases. Equation (33) has these properties. If matter is streaming almost freely from the loser to the gainer, the rate can be written as the product of density, sound velocity and effective area in the neck and thus in the form of Eq. (33). This is the case not only in semi-detached systems but also in systems in contact as long as . In detached systems without mass exchange the equation is trivially valid.
The general case of systems in contact is more complex since matter is streaming from the loser to the gainer (at the rate , say) as well as from the gainer to the loser (at the rate ). Each rate separately can be approximated by a product of characteristic values for density, velocity (usually smaller than the sound velocity) and effective area. For the net rate however an expression of this form is not obvious, except for the case (already mentioned) that ans thus . Nevertheless Eq. (33) is plausible as an approximation since providing fulfilment of the Roche equipotential condition in a close approximation, with the possible exception of short phases during a cycle just before contact is broken or just after contact has been reestablished. As long as in these phases, Eq. (33) is certainly adequate. For these reasons Eq. (33), with taken in the range , is expected to be reasonable in any system.
With the rate of mass exchange also the accretion luminosity is known from Eq. (8). The luminosity transferred by circulation currents however is still unknown.
It is convenient to use the variable
xi=mi/Mi and to write
the instantaneous distribution of the energy sources and sinks in the form
So far it has almost always been assumed that the energy sources and sinks occur only in the outer layers. For this reason it is useful to introduce a simple and smooth normalised function which is zero in the interior and positive in the outer layers. The function
Sources and sinks caused by circulation currents occur
only in systems in contact. Since the
fitting mass is very close to unity, the contact condition
The uncertainties concerning the energy transfer by circulation currents have
been expressed as uncertainties in the quantities
We shall assume that these quantities
are (in way so far unknown) determined by the instantaneous spherically averaged structure of both components and thus by the set of functions
The equations derived in the preceding subsections contain the parameters
and the unknown functions listed in Eq. (45).
The presence of these parameters and functions reflects the difficulties
of the hydrodynamical problem. Concerning the choice of the parameters there
is little freedom. We decided to use as standard values
When studying the thermal evolution of unevolved systems, however, we are interested only in the long-time behaviour, which is expected to be either a stable limit-cycle (i.e. a thermal cycle) or a fixed point (i.e. a configuration in stable thermal equilibrium). More complex attractors cannot be excluded, but previous results (Kähler et al. 1986) suggest that they are not realised. Under these circumstances the problem is less difficult and a spherically averaged treatment of the components is reasonable.
To see this, consider first a given system evolving in relaxed thermal cycles. All properties (e.g. ) are periodic functions of the time t, and t can certainly be eliminated by inverting functions describing the spherically averaged structure of the components as a function of t. For example, if the primary's degree of contact F1(t) is a monotonic function in a certain time interval, then t=t(F1) and thus in this interval. Accordingly, for a given system evolving in cycles the functions (45) can be defined in a variety of ways. We are interested in a physical definition which is (approximately) valid not only for the particular system under consideration but for all systems evolving in cycles as well for all systems in (stable or unstable) thermal equilibrium. Such a definition can hopefully be obtained in the future from 3-D hydrodynamic calculations. In the meantime we have to make plausible assumptions and to introduce free parameters. Equation (46) is an example.
Suppose next that thermal equilibrium can be achieved. In an approximate discussion the uncertainties in the parameters listed in Eq. (47) can be neglected. The structure of a configuration is then determined by the values of , i.e. by few numbers, not functions. The investigation of possible equilibrium configurations is therefore much less difficult than the investigation of thermal cycles. In fact, it is possible to investigate all candidates for reasonable equilibrium configurations. We are interested in configurations which are dynamically as well as thermally stable. The dynamical stability problem is well determined. More difficult is the thermal stability problem since it involves the change of the quantities in the course of thermal perturbations, and thus, if these quantities can be treated as functions of , the derivatives of these functions.
Some choice of the functions in Eq. (45) being necessary,
here we define a simple model for late-type contact binaries involving
only few parameters in addition to those in Eq. (47).
We assume that the energy flux in the neck is not larger than
In view of the persistent difficulties in the theoretical treatment of contact binaries we cannot expect that this model gives results which are fully compatible with the observations as well as with theoretical constraints. We do, however, expect that the deficiencies of the results can be used to reduce the inconsistencies and to improve the model. In other words, the model is intended as a starting point in the attempt to find a consistent description of contact binaries.
The non-Lagrangian mass variables xi=mi/Mi are used. Since the components exchange mass, this choice minimises the changes in pressure and temperature at the meshpoints near the surface during a time step.
As expected, numerical stability turned out to be sometimes a serious problem. While in quiet evolutionary phases time steps of more than 103 years can be used, in some difficult phases steps of 10-4 years are required. We found it necessary to treat the difference equations as proposed by Sugimoto (1970).
A problem with the temperature gradient concerns only the secondary's envelope. In some phases during a cycle the extent of the convective layers may not converge in the iteration process. Some layers oscillate between radiative transport and convection. This difficulty is removed if the extent of radiative and of turbulent layers is artificially kept fixed after a certain number of iterations.
Concerning the choice of an initial configuration, one of the equations in the Henyey scheme deals with the exchange of mass between the components. If this equation is replaced by an equation prescribing the mass ratio, and if artificial energy sources/sinks are introduced in the secondary's/primary's interior, systems in thermal equilibrium with any given mass ratio can easily be constructed. The energy sources and sinks can be adjusted to obtain a configuration in contact satisfying the Roche equipotential condition. Switching then off the artificial energy sources and sinks we obtain a contact configuration in thermal disequilibrium which, if dynamically stable, can be used as initial configuration.
We have collected the equations for zero age contact systems, assuming conservation of angular momentum. (In our opinion it is premature to study the effects of nuclear evolution and loss of angular momentum as long as the structure of unevolved systems is not well understood.) Are these equations based on truly weak assumptions? The continuity equation is beyond doubt. The equation of hydrostatic equilibrium (2) is a very close appoximation in the interior, and a close approximation also in the outer layers if (i) the velocity of the internal mass motions (in a rotating frame) is not larger than the the sound velocity and thus much smaller than the orbital velocity and (ii) the turbulent pressure is unimportant. If the velocity of the internal mass motions cannot be neglected, modifications of the equation of hydrostatic equilibrium and of the Roche equipotential condition will be necessary as described by Kähler (1995), and this requires a modification of Eq. (33). There are also uncertainties in the transport equation which can be estimated comparing results for different assumptions (i.e. different values of ). Serious uncertainties concern only the energy balance. Some quantities decribing the energy transfer between the components are unknown. When investigating thermal cycles, these quantities can be treated as (so far unknown) functions depending on properties of the spherically averaged structure of the components. An evolution code has been presented which can be adapted to any choice of these functions. With this freedom and with the possibility to introduce other modifications (concerning the temperature gradient and departures from hydrostatic equilibrium in the outer layers) the treatment of contact binaries is free from restrictive assumptions.
The transfer of mass and the exchange of energy between the components have been treated separately. In reality these processes are coupled since determined by the internal mass motions. A combined treatment will be possible only when we get more insight into these motions. A simple model has been proposed as a starting point for more detailed investigations, and in an accompanying paper first results are presented.