Apsidal advance in SS 433?
M. G. Bowler
University of Oxford, Department of Physics, Keble Road, Oxford, OX1 3RH, UK
Received 14 October 2009, Accepted 15 November 2009
Context. The Galactic microquasar SS 433 launches oppositely-directed jets at speeds approximately a quarter of the speed of light. Both the speed and direction of the jets exhibit small fluctuations. A component of the speed variation has 13 day periodicity and the orbital phase at which its maximum speed occurs has advanced approximately 90 in 25 years.
Aims. To examine the possibility that these variations are associated with a mildly eccentric orbit and conditions necessary to achieve this apsidal advance.
Methods. The advance of the orbital phase for maximum speed is taken to be advance of the apses of the putative elliptical orbit. It is compared with calculations of the effects of tides induced in the companion and also with gravitational perturbations from the circumbinary disc. These calculations are made in the light of recent results on the SS 433 system.
Results. The 13-day periodicity in the speed of the jets of SS 433 might be attributed to a mildly elliptical orbit, through periodic approaches of the donor and the compact object. Advance of the apses of such an elliptical orbit due to tidal effects induced in a normal companion looks to be too small; if caused by the circumbinary disc the mass of the inner regions of that disc is .
Key words: stars: individual: SS 433 - binaries: eclipsing
The relativistic jets of SS 433 have a mean speed of about 0.26 c but that speed is not constant; it fluctuates with rms deviation . These fluctuations are rather symmetric between the jets and a component of the fluctuations has a period of about 13 days, the period of the binary orbit (Blundell & Bowler 2004, 2005; Blundell et al. 2007). The orbital phase of maximum speed of this component has advanced by in 25 years (Blundell et al. 2007). The mechanisms whereby the jet speed is modulated have not as yet been elucidated, but the periodic component suggests that variation of separation between the compact object and the companion is implicated; an elliptical orbit. Such an elliptical orbit has also been suggested as a reason for flares being associated with particular orbital phases claimed by Fabrika & Irsmambetova (2002).
If the orbit of the binary is elliptical then the volume of the Roche lobe of the companion will reach a minimum once each orbit for eccentricity e as small as 0.01 (Bland & Grindlay 1984). Periodic approaches might increase the rate of transfer of material to the accretion disc or otherwise disturb its outer regions. It may be relevant that a period of higher than usual jet speeds set in on or just before JD 2453294 (Blundell et al. 2007), coinciding with an optical outburst and a radio flare (Schmidtobreick & Blundell 2006). The ratio of the putative period of apsidal advance to the period of the system, 2796, is not atypical of massive close binaries where apsidal advance is induced by tidal effects (J. I. Katz private communication); for Y Cygni the ratio is 5791. The mass ratio of the SS 433 system and the masses of the compact object and of the companion have been much better established as a result of detailed studies of the H structure which revealed the circumbinary disc (Blundell et al. 2008); the system is now known to be massive (roughly 40 ). In the light of these new data, it seems worth while to examine conditions under which an elliptical orbit would precess at a rate of 90 in 25 years. Tidal distortion of the companion star by the compact object is one possible mechanism, another is the effect of gravitational perturbations due to the circumbinary disc itself. These are addressed in turn.
The classic work on apsidal motion in binary systems as a result of tidal distortions is Sterne (1939).
In that treatment distortion is along the line joining the centres
(equilibrium tides) and account is taken of the variation of tidal
deformations with time. In this note his Eq. (14a) is applied,
making the following assumptions. First, that the compact object (plus
accretion disc) can be treated as a point mass (but in the context of
SS 433 see Collins & Newsom 1988).
Secondly, that the companion co-rotates with the binary system and
finally that terms in his Eq. (14) which involve the eccentricity e (which enters as the square) can be ignored. Equation (14a) of Sterne (1939) then reads
Here is the frequency of the induced apsidal advance and the orbital frequency. The structure constant is for the companion star (it would be equal to 0.0144 for a polytrope of n=3). The orbital speed of the circumbinary disc has revealed that q is greater than 0.7 (Blundell et al. 2008); it probably lies in the range 0.7 - 1.0. The quantity a is the radius of the companion and A is the semi-major axis of the binary orbit, equal to the separation of the two components for negligible ellipticity. For an orbit of period 13.08 days and periapsis advance of 90 in 25 years, the left hand side of Eq. (1) is .
The ratio a/A is given by the length of time of
total eclipse by the companion of an ideal luminous point in the orbit
of the compact object. The length of time
for which such an object is eclipsed when viewed in the orbital plane is
where is measured in days. The orbital plane is in fact inclined at 12 degrees to the line of sight and the actual eclipse does not last so long, rather
The primary eclipses in SS 433 have been studied in both the optical and X-ray regions of the spectrum (see for example Goranskii et al. 1998; Cherepashchuk et al. 2005) but in neither case is a pointlike object being eclipsed. The ideal eclipse time P must however be close to a single day; for P = 1 day, a/A = 0.32. Another indication of the value of a/A is obtained from the value of the mass ratio q and the assumption that the companion fills its Roche lobe. The Roche lobe radius of the companion, in units of A, is a function only of q (Eggleton 1983) and hence . Thus Eq. (1) is easily satisfied for log , provided that this value is compatible with the companion mass of and the radius km. The latter corresponds to a surface gravity given by log (g in units of cm s-2). Here there seems to be a problem - a normal star is too small if log or too concentrated after expanding to the required radius. The close OB stars in Y Cygni (Hill & Holmgren 1995) provide an example where the masses and radii are known and log are -1.94, inferred from the apsidal precession. The stars each have mass , radii and hence log g - they are too small to match the properties of the companion in SS 433. The Y Cygni system is fairly young, but log becomes increasingly negative as the star slowly expands. Stellar models have been evolved by Claret & Gimenez (1992) and later by Claret (2004) and their tables list both log k2 and log g. The calculated structure constants k2 which determine the effects of equilibrium tides on the precession of the apses of a binary orbit survive very well comparison with close binary systems (Claret & Gimenez 1993). These evolved models have, for suitable masses, log as their radii approach the companion radius (which is close to that of its Roche lobe). Such a structure is too concentrated to drive the suggested apsidal precession in SS 433. The tables of Claret (2004) also list structures with log computed for old bloated stars, but for isolated stars the surface gravity is then log g < 0, corresponding to a surface radius much larger than the separation of the two components of the SS 433 binary; diffuse outer regions would have been stripped off. Thus the companion must have a distinctly peculiar structure if it is responsible for the precession of the apses at a rate of one complete rotation in 100 years. It has certainly endured a peculiar history.
It is very simple to estimate the rate at which the apses of an
elliptical orbit advance as a result of a circumbinary disc. It is not
necessary to consider tidal distortions because the presence of a ring
of matter, external to the binary itself, augments the attractive
inverse square radial force on each component with a linear term
directed outward. Very little is known about the circumbinary disc; the
calculations below are for the simplest case of a ring of circumbinary
material coplanar with the binary orbit. If the ring of matter orbits
the centre of the binary at a radius R and has (linear) mass density ,
then the radial force per unit mass F on a component at radius r is given by
for q=1 and component masses M, ignoring terms in r3/R3. The ratio r/R is about 0.2; for present purposes Eq. (4) is adequate.
In the presence of the linear term the frequency of small radial oscillations about a radius
is smaller than the rotational frequency at that radius and hence the apses advance. The vibrational frequency
is given by
and the circular frequency by
where the small quantity x is given by
and one complete rotation of the ellipse (100 years) takes 1/6x 13.08 day periods. Hence if M is 20 the mass of the inner part of the circumbinary disc must be 0.15 if it is responsible for an advance of the periapsis of in 25 years.
This work is speculative in the sense that it is by no means clear that the observed oscillation of jet speeds with the period of the orbit is associated with the approach of the two components in an elliptical orbit. The advance of the phase corresponding to maximum speed might however be explained by apsidal advance of the orbit.
Both tidal distortion of the companion and the presence of the circumbinary disc will advance the periapsis; if the effects of tidal distortion of the companion are responsible then log must be accompanied by log . Such a structure would be remarkable.
The rate of advance of periapsis induced by the circumbinary disc would match that of the jet speed maximum for a rather modest mass in the disc - approximately 0.15 , which corresponds to about 1000 years of mass transfer.Acknowledgements
I thank J. I. Katz for comments.
- Band, D. L., & Grindlay, J. E. 1984, ApJ, 285, 702 [NASA ADS] [CrossRef] (In the text)
- Blundell, K. M., & Bowler, M. G. 2004, ApJ, 616, L159 [NASA ADS] [CrossRef] (In the text)
- Blundell, K. M., & Bowler, M. G. 2005, ApJ, 622, L129 [NASA ADS] [CrossRef] (In the text)
- Blundell, K. M., Bowler, M. G., & Schmidtobreick, L. 2007, A&A, 474, 903 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Blundell, K. M., Bowler, M. G., & Schmidtobreick, L. 2008, ApJ, 678, L47 [NASA ADS] [CrossRef] (In the text)
- Cherepashchuk, A. M., Sunyaev, R. A., Fabrika, S. N., et al. 2005, A&A, 437, 561 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Claret, A. 2004, A&A 424, 919 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Claret, A., & Gimenez, A. 1992 A&AS 96, 255 (In the text)
- Claret, A., & Gimenez, A. 1993 A&A 277, 487 (In the text)
- Collins, G. W., & Newsom, G. H. 1988 ApJ. 331, 486 (In the text)
- Eggleton, P. P. 1983 ApJ, 268, 368 (In the text)
- Fabrika, S., & Irsmambetova, T. 2002, in New Views on Microquasars, 4th Microquasars Workshop, May 27-June 1, Cargèse, Corsica, France, ed. Ph. Durouchoux, Y. Fuchs & J. Rodriguez (Kolkata, India: Center for Space Physics), 276, [arXiv:astro-ph/0207254] (In the text)
- Goranskii, V. P., Esipov, V. F., & Cherepashchuk, A. M. 1998, Astron. Rep., 42, 209 [NASA ADS] (In the text)
- Hill, G., & Holmgren, D. E. 1995, A&A, 297, 127 [NASA ADS] (In the text)
- Schmidtobreick, L., & Blundell, K. M. 2006, in Proc. VI Microquasar Workshop, 2006 September 18-22, Como, Italy (In the text)
- Sterne, T. E. 1939, MNRAS, 99, 451 [NASA ADS] (In the text)
Copyright ESO 2010