Next: Conclusions Up: An analysis of the Previous: The geometry of the
The derivative of the orbital period and the mass loss rate
In Sect. 4 it was found that the orbital ephemeris has a significant negative quadratic term, indicating a decreasing period. The derivative of the period is given by dP/dt = 2b/P, where b is the coefficient of E2 of the quadratic ephemeris. Thus, dP/dt=-(1.30+-0.04) 10-10 if P=1 and dP/dt=-(2.60+-0.09) 10-10 if P=P2 [note that b(P2=4 b(P1)]. The relative period decrease is then (dP/dt)/P = -(1.15+-0.04) 10-9 d11.
The period derivative is three orders of magnitudes larger than predicted by GB83 who assumed gravitational radiation to be the sole agent for angular momentum (AM) loss. Qualitatively, this difference is readily explained if the mass loss which gave origin to the formation of the planetary nebula is still ongoing.
In order to quantify this idea we assume that one or both binary components lose mass via a stellar wind. The ejected mass takes away rotational AM of the stars. In a close binary system such as MT Ser bound rotation may be assumed. Thus, rotational AM loss will result in a decrease of the binary orbit. In order to derive constraints on the mass loss from the system, me may therefore, in a reasonable approximation, equate the orbital angular momentum loss as indicated by the negative period derivative to the rotational AM loss of the component stars due to a stellar wind. This is an approximate approach because regarding the AM balance both, orbital and rotational contributions have to be considered, whereas equating the above mentioned quantities only the orbital angular momentum is regarded. While the orbital motion is certainly responsible for the major part of the total AM of the system, in view of the relative sizes of the components as derived in Sect. 5, and depending on their internal structure, the rotational angular momentum may not be negligible altogether. However, for an approximate determination of the mass loss from the system it will suffice to regard only the orbital AM loss here.
In a strongly idealized scenario we assume an isotropic wind to emanate
from the stellar surface and to take away an amount of angular momentum
corresponding to the AM of the wind mass due to the stellar rotation at
the location where it left the star. The quantification of this concept
for the individual stellar surface elements and the integration over the
entire surface yields a rate of angular momentum loss of
where dot(M)=dM/dt is the entire mass loss rate from the star,
Rerr is
the effective stellar radius (assuming a spherical star) which - for the time
being, but see below - is taken to be the geometrical radius, and P is
the rotational period which is considered here as identical to the orbital
period.
In an inertial reference frame the orbital AM, Jorb, (assuming a circular orbit) can be expressed as a function of the period, P, the masses, M1 and M2, of the components (assumed to be point-like) and the orbital radius, a1, of the primary star. The temporal derivative of the orbital AM, dJorb/dt, is then the sum of the partial derivatives with respect to P, M1, M2 and a1. Of course, in a Keplerian orbit these are not independent. Since the derivative of the orbital period, dP/dt, is a measured quantity, and because we are interested in the mass loss rates (i.e. dM1/dt and dM2/dt) we use Kepler's third law to express a1 in terms of P, M1 and M2. This yields da1?dt in terms of dP/dt, dM1/dt and dM2/dt. Using the corresponding expression to eleminate da1/dt from the equation for dJorb/dt, we finally get a (somewhat lengthy) expression for the derivative of the orbital angular momentum which depends on P, M1, M2, a1, dP/dt, dM1/dt and dM2/dt. Except for the latter two quantities all are known either from direct observations or were derived within a specific model (see Table 6). Thus, choosing values for dM1/dt and dM2/dt such that dJrot,1/dt + dJrot,2/dt = dJorb/dt gives us an estimate (within the numerous approximations involved) for the mass loss of the system.
For Model 1 we assume that only the primary (evolved) star loses mass. Then, dMtot/dt = dM1/dt +dM2/dt = dM1/dt. For Model 2 we investigate the cases dMtot/dt=dM1/dt; dMtot/dt=dM2/dt and dM1/dt = dM2/dt = 0.5 dMtot/dt. Within Model 1 agreement between dJrot/dt and dJorb/dt is achieved for mass loss rates of 1.8 10-7 Moy-1, 7.3 10-8 Moy-1 and 2.4 10-8 Moy-1 in the H, I and L cases, respectively.
In contrast, no solution exists for Model 2. For all mass loss rates the rotational angular momentum loss is less than the predicted orbital angular momentum loss. This is only different, if an ``efficiency factor'' f for the rotational AM loss is introduced, assuming that Reff = fR, where R is the true geometrical radius of the star. This may be thought of as roughly simulating a magnetically coupled wind which dynamically decouples from the stellar rotation only at the radius Reff. It is then found that solutions for dM/dt such that dJrot/dt = dJorb/dt exist for efficiency factors f>=3.5 for Model 2.2.1, and - in the case of Model 2.2.2 - for f>=2.5 (dMtot/dt = dM1/dt), f>=3.5 (dM1/dt = dM2/dt) and f>=5 (dMtot/dt = dM2/dt). The functional dependence between the efficiency factor and the total mass loss rate is shown in Fig. 4 for Model 1.1 (upper frame), Model 2.2.1 (case dM1/dt = dM2/dt; middle frame) and Model 2.2.2 (case (dM1/dt = dM2/dt); lower frame).
