Non-Periodic Phenomena in Variable Stars
                                                IAU Colloquium, Budapest, 1968

             
         
                  SUDDEN CHANGES IN THE PERIOD OF ALGOL

                              T. HERCZEG
                         Hamburg - Bergedorf 
                        (Read by W. Seitter)


  Orbital motion in binary systems was considered, for a long time, a true 
paradigm of periodic phenomena. Even its occasional changes were 
supposed to be of strictly regular, i.e. periodic nature like rotation 
of the apsidal line or light time effects in multiple systems. 
Exceptional objects like beta Lyrae, with a secular increase of the period, 
could have been recognized as cases still showing some regular, 
predictable character. Later on, however, after longer series of 
observations became available, a sort of eclipsing systems (like XZ And) 
were found showing obviously non-periodic changes of the period: an 
erratic but continuous up and down in the O-C diagrams for the times of 
minima. Other well observed photometric binaries exhibited a different 
type of non-periodic changes: the eclipsing period was, from time to 
time, subject to sharp sudden changes ("jumps") while between two 
consecutive jumps it remained practically constant. AR Lacertae or AH 
Virginis are important examples of this type.
  Schneller (1962) in a lecture at the Budapest Observatory advocated his 
thesis that, while "regular" changes of the eclipsing period are 
relatively infrequent, these sudden, erratic changes seem to be the 
rule, at least in the case of the contact and semi-detached systems. 
Among other typical objects he mentioned Algol (beta Persei) itself 
- certainly one of the best observed variables - as indicating period 
changes of this nature; a similar proposal concerning Algol's period, as 
a matter of fact, was already put forward by Sterne.
  In the case of Algol, Schneller tried to represent the O-C curve of the 
times of minima by a polygon having two angles, two "sharp bends". They 
correspond to sudden changes, discontinuities of the period, having 
occurred in 1840 and 1924. This proposal has been essentially 
substantiated by a recent analysis of all photoelectrically determined 
epochs of minima, carried out at the Hamburg observatory (cf. a forthcoming 
paper by Frieboes-Conde, Herczeg and Hog). In particular, the occurrence of 
"period jumps" seems to be definitely established now though numerical details 
(like the data of these changes) had to be revised.
  This interpretation of the O-C curve offers us an explanation of the 
much discussed "great inequality" of Algol, too. This is a hypothetical 
long period light time effect requiring the existence of an additional 
unseen companion (Algol D) in the system. Up to now there are three observed 
and, at least, two hypothetical members listed as belonging to the complicated 
Algol system. Those observed, definitely recognized components are:
  Algol AB the eclipsing pair, a semi-detached system;
  Algol C, a distant companion with a period of 1.862a. This orbital motion
in the triple system gives rise to a light time effect of the times of 
minima; having a small amplitude (~~ 6 min.), it is rather difficult to 
handle. On the other hand, there exists an accurate and reliable set of 
spectroscopic elements for this 1.9 year orbit, derived by Ebbighausen (1958).
  The hypothetical members are:
  Algol E, introduced for explaining by light time effect a clearly 
distinguishable 32-year periodicity of the O-C values by light time effect. 
This periodicity can be, however, much more convincingly explained by the aid 
of apsidal motion in the eclipsing system;
  Algol D, a component proposed already in the last century in order to 
find an explication for the above mentioned great inequality, dominant 
feature of the O-C diagram, with a suggested period of about 170-180 
years. But an interpretation of the great inequality as a light time 
effect faces also serious difficulties. Until now, the periodicity 
itself defied exact representation: predictions based upon a recent, 
very elaborate discussion of the period changes by Kopal, Plavec and 
Reilly (1960) yielded residuals amounting to half an hour. Besides, the 
mass of this hypothetical unseen companion turns out to be unexpectedly 
high, about 4 M_Sun. A further important question is the following: a 
possible light time effect would require an orbit large enough to cause 
detectable changes (> 1 sec of arc) in the position of Algol AB, a much 
observed 2nd magnitude star.
Our discussion of all meridian observation reaching back to about A. D. 
1750 indicated, however, that Algol's proper motion remained sensibly 
constant - no traces of the orbital motion in the hypothetical quadruple 
system Algol ABC-D could be found. This seems to be a direct and decisive 
argument against the existence of the component Algol D.
  We are now obliged to find an alternative explanation for the great 
inequality. Let us make the basic assumption that - besides the eclipsing 
period, the 1.9-year light time effect and the 32-year period - no further 
periodicities exist in the O-C diagram. (This assumption, which might be 
opposed by some observers, has turned out a very useful and also successful 
working hypothesis.) Then a rather sensitive test can be carried out, based on 
the best determined, photoelectrically observed epochs of minima, about 60 
in number since 1920. We treated this rather accurate observational material 
in the following way:
  We subtracted the effect of 32-years periodicity. This correction was 
based on data by Hellerich (1919);
  We accepted as a definite representation of the 1.862a orbit the spectroscopic 
elements given by Ebbighausen (1958, 1962). Thereafter, the corresponding light 
time effect can easily be taken into consideration.
  Then the residuals have shown with considerable accuracy a polygon indicating 
two sudden changes of the period (in 1944 and 1952) and three different values 
of a practically constant eclipsing period outside these "jumps" (Fig. 1). 
The sudden changes amount to Delta P = +3.5s and Delta P = -2.0s, respectively. 
The constancy of the period within the time intervals, say 1920-1942 and 
1955-1965 can be judged from the good representation of the 1.9-year light 
time effect, the remaining scatter of the O-C values being but slightly higher 
than the usual error of photoelectrically determined epochs of minima. It seems 
very improbable that any further periodicities could be accomodated within 
this small margin.




                              Fig. 1
  

  The reliability of Hellerich's numerical representation of the 22-year 
period is, of course, a basic question. Though small systematic deviations are 
not to be entirely excluded, no serious error can, - in my opinion - arise 
hereby. Especially the fact that the two jumps are rather close together, 
separated by a quarter of the period only, means that the existence of 
two sudden changes remains incontestable.
  This somewhat perhaps surprising representation of the photoelectric 
times of minima encourages me to an experiment: to apply a similar 
procedure to the earlier measurements, too. Minima were observed since 
1784 in a very great number. The overwhelming majority of these times of 
minima is based on visual estimates, and their rather modest accuracy 
(of the order of +-0.02d), doesn't permit a detection of the 1.9-year 
light time effect. Besides, the in fluence of minor systematic errors in 
the representation of the 32-year period will be considerably enhanced 
by going back to the 18th century observations. This makes our proposed 
representation to a tentative one; on the other hand the deviations from 
a reasonable picture turn out relatively small, thus suggesting the 
proposal I am going to discuss is not completely arbitrary.
  Because of the large scatter of the individual epochs I used the normal 
points meticulously derived by Ferrari (1934). Again, I subtracted Hellerich's 
formula for the 32-year term. The residuals not only allow, they clearly 
suggest a representation given by a polygon. The representation derived earlier 
from the photoelectric minima is a continuation of this new polygon - the long 
interval of constant period 1920-1944 can well be extended back to 1915. 
There is no indication of a sudden period change in 1924 as it was proposed 
by Schneller; the short sharp change around 1914/15 seems to be real (Fig. 2). 
  The following table summarizes this possible sequence of "events": 


        P_1 = 2.8673442d       (1784)- 1835 (approx.) Delta P ~~-3.7s 
	P_2 = 2.8673012d	1835 - 1854                   ~~-2.3s
	P_3 = 2.8672775d	1854 - 1901		      ~~+1.7s		
	P_4 = 2.8672967d	1901 - 1913		      ~~-4.0s		
	P_5 = 2.8672506d	1913 - 1915		      ~~+5.0s		
	P_6 = 2.86730807d	1915 - 1944		      ~~+3.5s		
	P_7 = 2.86734862d	1944 - 1952		      ~~-2.0s		
	P_8 = 2.86732442d	1952 -(1965)				




                              Fig. 2


  This interpretation decomposes the great inequality into a series of 
sudden period changes. Apparently they are forming a random set, the 
changes occurring on an average once in 25-30 years. This is certainly a 
non-periodic phenomenon, probably caused by discontinuities of the mass 
exchange between the components or in the mass loss from the whole 
system. Perhaps, the "overflowing" of the secondary component of its 
Roche-limit is not a smooth phenomenon but it causes from time to time 
directional outbursts of mass from the secondary, hereby changing the 
eclipsing period in a sudden way.
  The obviously spurious periodicity of about 160-180 years length was an 
understandable suggestion of the early observers: it was favoured by the 
existence of two relatively long intervals of constant period. P_1 and P_3.
This led to the use of a mean period in the linear ephemeris formula 
which, in its turn, gave rise to the characteristic wedge-shaped figure 
in the O-C diagram, simulating the "first half" of a sine curve (Fig. 3).




                              Fig. 3


                              DISCUSSION

Detre: Algol's O-C diagram looks rather like to one resulting from 
     cumulative effects of random period variations: on long cycles there are 
     superposed shorter cycles. How far is the 32 year period established?
Seitter: Irrespective of smaller errors having an effect on the period 
     length and the zero phase (and Hellerich's representation of the 32 year 
     period seems to be a very reliable one!), most of the jumps, especially 
     those rather thoroughly observed in 1944 and 1952, are quite well 
     established. This is mainly based on the fact that these two jumps 
     occurred within a time interval short enough not to be seriously 
     affected by any erroneous assumption concerning the 32 year periodicity. 
     Numerical data could certainly be shifted to some extent but the 
     existence of sudden period changes remains beyond doubt.


                              REFERENCES

Ebbighausen, E. G., 1958, Astrophys. J. 128, 598.
Ebbighausen, E. G., and Gange, J. J., 1962. Publ. Dominion Astrophys. 
                    Obs. Victoria 12, 151.
Ferrari, K., 1934, Astr. Nachr. 253, 225.
Hellerich, J., 1919, Astr. Nachr. 209, 227.
Kopal, Z., Plavec, M. and Reilly, Edith, 1960, Jodrell Bank Ann. 1, 374. 
Schneller, H., 1962, Mitt. Sternwarte Budapest No. 53.(CoKon N°.53)