A third body in an eclipsing binary system causes regular periodic changes in the observed (O) minus the computed (C) eclipse epochs. Fourth bodies are rarely detected from the O-C data. We apply the new Discrete Chi-square method (DCM) to the O-C data of the eclipsing binary XZ Andromedae. These data contain the periodic signatures of at least ten wide orbit stars (WOSs). Their orbital periods are between 1.6 and 91.7 years. since no changes have been observed in the eclipses of XZ And during the past 127 years, the orbits of all these WOSs are most probably co-planar. We give detailed instructions of how the professional and the amateur astronomers can easily repeat all stages of our DCM analysis with an ordinary PC, as well as apply this method to O-C data of other eclipsing binaries.

Regular changes in the observed (O) and the computed (C) epochs of binary eclipses may reveal the presence of a third or a fourth body. One recent study showed that the probability for detecting a fourth body from the O-C data is only 0.00005. We apply the new Discrete Chi-square Method (DCM) to the O-C data of Algol, and detect three wide orbit stars having orbital periods 1.9 years (Algol C), 18.6 years (Algol D) and 52.5 years (Algol E). Since our estimate for the period of Algol C agrees perfectly with the previous estimates, the signals of Algol D and Algol E are certainly real. The orbits of all these three wide orbit stars are most probably co-planar, because no changes have been observed in eclipses of Algol.

Unambiguous detection of signals superimposed on unknown trends is difficult for unevenly spaced data. Here, we formulate the Discrete Chi-square Method (DCM) that can determine the best model for many signals superimposed on arbitrary polynomial trends. DCM minimizes the Chi-square for the data in the multi-dimensional tested frequency space. The required number of tested frequency combinations remains manageable, because the method test statistic is symmetric in this tested frequency space. With our known tested constant frequency grid values, the non-linear DCM model becomes linear, and all results become unambiguous. We test DCM with simulated data containing different mixtures of signals and trends. DCM gives unambiguous results, if the signal frequencies are not too close to each other, and none of the signals is too weak. It relies on brute computational force, because all possible free parameter combinations for all reasonable linear models are tested. DCM works like winning a lottery by buying all lottery tickets. Anyone can reproduce all our results with the DCM computer code. All files, variables and other program code related items are printed in magenta colour. Our Appendix gives detailed instructions for using dcm.py. We also present one preliminary real use case, where DCM is applied to the observed (O) minus the computed (C) eclipse epochs of a binary star, XZ And. This DCM analysis reveals evidence for the possible presence of a third and a fourth body in this system. One recent study of a very large sample of binary stars indicated that the probability for detecting a fourth body from the O-C data of eclipsing binaries is only about 0.00005.

Previous studies have related surface temperature maps, obtained with the Doppler imaging (DI) technique, of LQ Hya with long-term photometry. We compare surface magnetic field maps, obtained with the Zeeman Doppler imaging (ZDI) technique, with contemporaneous photometry, with the aim of quantifying the star's magnetic cycle characteristics. We inverted Stokes IV spectropolarimetry into magnetic field and surface brightness maps using a tomographic inversion code that models high signal-to-noise ratio mean line profiles produced by the least squares deconvolution (LSD) technique. The magnetic field and surface brightness maps reveal similar patterns to previous DI and ZDI studies: non-axisymmetric polar magnetic field structure, void of fields at mid-latitudes, and a complex structure in the equatorial regions. There is a weak but clear tendency of the polar structures to be linked with strong radial field and the equatorial ones with the azimuthal. We find a polarity reversal in the radial field between 2016 and 2017 coincident with an activity minimum seen in the long-term photometry. The inverted field strengths cannot easily be related with the observed spottedness, but we find that they are partially connected with the retrieved field complexity. Comparing to global magnetoconvection models for rapidly rotating young Suns, this field topology and dominance of the poloidal field component could be explained by a turbulent dynamo, where differential rotation does not play a major role (so called alpha^2 Omega or alpha^2 dynamos), and axi- and non-axisymmetric modes are excited simultaneously. The complex equatorial magnetic field structure could arise from the twisted (helical) wreaths often seen in these simulations, while the polar feature would be connected to the mostly poloidal non-axisymmetric component having a smooth spatial structure.

For seven decades, the widely held view has been that the formation, the migration and the decay of short-lived starspots explain the constantly changing light curves of chromospherically active stars. Our hypothesis is that these deceptive observed light curves are interference of two real constant period light curves of long-lived starspots. The slow motion of these long-lived starspots with respect to each other causes the observed light curve changes. This hypothesis contradicts the current views of starspots. Therefore, we subject it to eight reproducible tests. Our new period finding method detects the two real light curves of FK Com. Our hypothesis is a total success: all real light curve parameters are directly connected to the long-lived starspots which are also seen in the Doppler images of FK Com.These parameters are spatially and temporally correlated just like in the Sun, including weak solar-like surface differential rotation. As for other chromospherically active stars, all eight reproducible tests also support our hypothesis. It explains many spurious phenomena: the rapid light curve changes, the short starspot life-times, the rapid rotation period changes, the active longitudes, the starspot migration, the period cycles, the amplitude cycles and the minimum epoch cycles. It also explains why the light curves and the Doppler images give contradicting surface differential rotation estimates even for the same individual star, as well as the abrupt 180 degrees shifts of activity (the flip-flop events) and the long-term mean light curves. We argue that the current views of starspots need to be revised.

Starspots are important manifestations of stellar magnetic activity. By studying their behaviour in young solar analogues, we can unravel the properties of their magnetic cycles. This gives crucial information of the underlying dynamo process. Comparisons with the solar cycle enable us to infer knowledge about how the solar dynamo has evolved during the Sun's lifetime. Here we study the correlation between photometric brightness variations, spottedness and mean temperature in V889 Her, a young solar analogue. Our data covers 18 years of spectroscopic and 25 years of photometric observations. We use Doppler imaging to derive temperature maps from high-resolution spectra. We use the Continuous Period Search method to retrieve mean V-magnitudes from photometric data. Our Doppler imaging maps show a persistent polar spot structure varying in strength. This structure is centered slightly off the rotational pole. The mean temperature derived from the maps shows an overall decreasing trend, as does the photometric mean brightness, until it reaches its minimum around 2017. The filling factor of cool spots, however, shows only a weak tendency to anti-correlate with the decreasing mean brightness. We interpret V889 Her to have entered into a grand maximum in its activity. The clear relation between the mean temperature of the Doppler imaging surface maps and the mean magnitude supports the reliability of the Doppler images. The lack of correlation between the mean magnitude and the spottedness may indicate that bright features in the Doppler images are real.

An ancient Egyptian Calendar of Lucky and Unlucky Days, the Cairo Calendar (CC), assigns luck with the period of 2.850 days. Previous astronomical, astrophysical and statistical analyses of CC support the idea that this was the period of the eclipsing binary Algol three millennia ago. However, next to nothing is known about who recorded Algol's period into CC and especially how. Here, we show that the ancient Egyptian scribes had the possible means and the motives for such astronomical observations. Their principles of describing celestial phenomena as activity of gods reveal why Algol received the title of Horus.

Recently, we presented a general model for the light curves of chromospherically active stars, where the observed light curve is interference of two real constant period light curves of long-lived starspots. In this first paper, we make six specific questions which undermine this argument, because it contradicts the current widely held views about the stellar surface differential rotation and the starspots. Our aim is to answer these six questions. We present evidence that the long-lived starspots of our general model have already been detected in the earlier surface imaging studies. The Lomb-Scargle power spectrum method analysis of the real and the simulated data of FK Com reveals that this method fails to detect the two real constant period light curves of our general model. If our model is valid, this method gives incompatible period, amplitude and minimum epoch estimates telling nothing about the real periods, the real amplitudes and the real minimum epochs of the two real light curves. This would mean that all earlier one-dimensional period analyses of the light curves of chromospherically active stars have given spurious results which have been widely and uncritically accepted since the discovery of the starspots in the year 1947. However, we arrive at a dead end, because we can not solve the real light curves of FK Com. In our second paper, we solve these real light curves with a new two-dimensional period finding method, prove the validity of our general model, and answer all six questions made in this first paper.

The starspots on the surface of many chromospherically active binary stars concentrate on long--lived active longitudes separated by 180 degrees. The activity shifts between these two longitudes, the "flip-flop" events, have been observed in single stars like FK Comae and binary stars like $\sigma$ Geminorum. Recently, interferometry has revealed that ellipticity may at least partly explain the "flip-flop" events in $\sigma$ Geminorum. This idea was supported by the double peaked shape of the long--term mean light curve of this star. Here, we show that the long--term mean light curves of fourteen chromospherically active binaries follow a general model which explains the connection betweenm orbital motion, starspot distribution changes, ellipticity and \flip ~events. Surface differential rotation is probably weak in these stars, because the interference of two constant period waves may explain the observed light curve changes. These two constant periods are the active longitude period $(P_{\mathrm{act}})$ and the orbital period $(P_{\mathrm{orb}})$. We also show how to apply the same model to single stars, where only the value of $P_{\mathrm{act}}$ is known. Finally, we present a tentative interference hypothesis about the origin of magnetic fields in all spectral types of stars.

Studying RS CVn binaries is challenging, because in addition to spot activity, other effects such as mass transfer between the components and gravitational distortion of their spherical forms may distort their light curves. Such effects can, however, be removed from the data by subtracting a mean light curve phased with the orbital period. We study a quarter of a century of standard Johnson differential V photometry of the RS CVn binary BM CVn. Our main aims are to determine the activity cycles, the rate of surface differential rotation and the rotation period of the active longitudes of BM CVn. The Continuous Period Search (CPS) is applied to the photometry. The changes of the mean and amplitude of the light curves are used to search for activity cycles. The rotation period changes give an estimate of the rate of surface differential rotation. The Kuiper method is applied to the epochs of the primary and secondary minima to search for active longitudes. The photometry reveals the presence of a stable mean light curve (MLC) connected to the orbital period P_orb = 20d.6252 of this binary. We remove this MLC from the original V magnitudes which gives us the corrected V' magnitudes. These two samples of V and V' data are analysed separately with CPS. The fraction of unreliable CPS models decreases when the MLC is removed. The same significant activity cycle of approximately 12.5 years is detected in both V and V' samples. The estimate for the surface differential rotation coefficient, k >= 0.10, is the same for both samples, but the number of unrealistic period estimates decreases after removing the MLC. The same active longitude period of P_al = 20d.511 +- 0d.005 is detected in the V and V' magnitudes. This long-term regularity in the epochs of primary and secondary minima of the light curves is not caused by the MLC. On the contrary, the MLC hampers the detection of active longitudes.

We study a sample of 21 young and active solar-type stars with spectral types ranging from late F to mid K and characterize the behaviour of their activity. We apply the continuous period search (CPS) time series analysis method on 16 to 27 years of photometry to estimate the surface differential rotation and determine the existence and behaviour of active longitudes and activity cycles on the stars. We supplement the time series results by calculating new $\log{R'_{\rm HK}}$ emission indices for the stars from high resolution spectroscopy. The photometric rotation period variations reveal a positive correlation between the relative differential rotation coefficient and the rotation period, but do not reveal any dependence of the differential rotation on the effective temperature of the stars. Secondary period searches reveal activity cycles in 18 of the stars and temporary or persistent active longitudes in 11 of them. The activity cycles fall into specific activity branches. We find a new split into sub-branches, indicating multiple simultaneously present cycle modes. Active longitudes appear to be present only on the more active stars. There is a sharp break at approximately $\log{R'_{\rm HK}}=-4.46$ separating the less active stars with long-term axisymmetric spot distributions from the more active ones with non-axisymmetric configurations. In seven out of eleven of our stars with clearly detected long-term non-axisymmetric spot activity the estimated active longitude periods are significantly shorter than the mean photometric rotation periods. This systematic trend can be interpreted either as a sign of the active longitudes being sustained from a deeper level in the stellar interior than the individual spots or as azimuthal dynamo waves exhibiting prograde propagation.

The Ancient Egyptians wrote Calendars of Lucky and Unlucky Days that assigned astronomically influenced prognoses for each day of the year. The best preserved of these calendars is the Cairo Calendar (hereafter CC) dated to 1244--1163 B.C. We have presented evidence that the 2.85 days period in the lucky prognoses of CC is equal to that of the eclipsing binary Algol during this historical era. We wanted to find out the vocabulary that represents Algol in the mythological texts of CC. Here we show that Algol was represented as Horus and thus signified both divinity and kingship. The texts describing the actions of Horus are consistent with the course of events witnessed by any naked eye observer of Algol. These descriptions support our claim that CC is the oldest preserved historical document of the discovery of a variable star. The period of the Moon, 29.6 days, has also been discovered in CC. We show that the actions of Seth were connected to this period, which also strongly regulated the times described as lucky for Heaven and for Earth. Now, for the first time, periodicity is discovered in the descriptions of the days in CC. Unlike many previous attempts to uncover the reasoning behind the myths of individual days, we discover the actual rules in the appearance and behaviour of deities during the whole year.

We investigate the spot activity of the young chromospherically active main sequence star LQ Hya. Our aims are to identify possible active longitudes, estimate the differential rotation and study long and short term changes in the activity. Our analysis is based on 24 years of Johnson V-band photometry. We use the previously published Continuous Period Search (CPS) method to model the evolution of the light curve of LQ Hya. The CPS fits a Fourier series model to short overlapping subsets of data. This enables us to monitor the spot configuration of the star with a higher time resolution. We find seasonal variability in the mean level and amplitude of the light curve of LQ Hya. The variability of the light curve amplitude seems not to be cyclic, but the long-term variations in the mean magnitude could be explained by an approximately 13 year cycle. Because of the limited length of the observed time series, it is not yet possible to determine whether this structure really is periodic and represents an activity cycle. We estimate the differential rotation of the star to be small, and the star is potentially very close to a rigid rotator. We search for active longitudes and find that on time scales up to six months there are typically one or two relatively stable active areas on the star with limited phase migration. On time scales longer than one year, no stable active longitudes have been present except for the period between 2003 and 2009 and possibly also some time before 1995. We find any signs of flip-flops with a regular period. The mean time scale of change of the light curve during the observation period is determined to be of the same order of magnitude as the predicted convective turnover time for the star.

We model the photometry of RS CVn star $\sigma$ Geminorum to obtain new information on the changes of the surface starspot distribution, i.e., activity cycles, differential rotation and active longitudes. We use the previously published Continuous Periods Search-method (CPS) to analyse V-band differential photometry obtained between the years 1987 and 2010 with the T3 0.4 m Automated Telescope at the Fairborn Observatory. The CPS-method divides data into short subsets and then models the light curves with Fourier-models of variable orders and provides estimates of the mean magnitude, amplitude, period and light curve minima. These light curve parameters are then analysed for signs of activity cycles, differential rotation and active longitudes. We confirm the presence of two previously found stable active longitudes, synchronised with the orbital period $P_{\rm{orb}}=19.60$d and find eight events where the active longitudes are disrupted. The epochs of the primary light curve minima rotate with a shorter period $P_{\rm{min,1}}=19.47$d than the orbital motion. If the variations in the photometric rotation period were to be caused by differential rotation, this would give a differential rotation coefficient of $\alpha \ge 0.103$. The presence of two slightly different periods of active regions may indicate a superposition of two dynamo modes, one stationary in the orbital frame and the other one propagating in the azimuthal direction. Our estimate of the differential rotation is much higher than previous results. However, simulations show that this can be caused by insufficient sampling in our data.

Our aim was to study the light curve periodicity of the solar analogue V352 CMa and in particular show that the presence or absence of periodicity in low amplitude light curves can be modelled with the Continuous Period Search (CPS) method. We applied the CPS method to 14 years of V-band photometry of V352 CMa and obtained estimates for the mean, amplitude, period and minima of the light curves in the selected datasets. We also applied the Power Spectrum Method (PSM) to these datasets and compared the performance of this frequently applied method to that of CPS. We detected signs of a $11.7 \pm 0.5$ year cycle in in the mean brightness. The long-term average photometric rotation period was $7.24 \pm 0.22$ days. The lower limit for the differential rotation coefficient would be $|k| > 0.12$, if the law of solar surface differential rotation were valid for V352 CMa and the period changes traced this phenomenon. Signs of stable active longitudes rotating with a period of $7.157 \pm 0.002$ days were detected from the epochs of the light minima with the Kuiper method. CPS performed better than the traditional PSM, because the latter always used a sinusoidal model for the data even when this was clearly not the correct model.

The eclipses in binary stars give precise information of orbital period changes. Goodricke discovered the 2.867 days period in the eclipses of Algol in the year 1783. The irregular orbital period changes of this longest known eclipsing binary continue to puzzle astronomers. The mass transfer between the two members of this binary should cause a long-term increase of the orbital period, but observations over two centuries have not confirmed this effect. Here, we present evidence indicating that the period of Algol was 2.850 days three millenia ago. For religious reasons, the ancient Egyptians have recorded this period into the Cairo Calendar, which describes the repetitive changes of the Raging one. Cairo Calendar may be the oldest preserved historical document of the discovery of a variable star.

According to earlier Doppler images of the magnetically active primary giant component of the RS CVn binary II Peg, the surface of the star was dominated by one single active longitude that was clearly drifting in the rotational frame of the binary system during 1994-2002; later imaging for 2004-2010, however, showed decreased and chaotic spot activity, with no signs of the drift pattern. Here we set out to investigate from a more extensive photometric dataset whether such a drift is a persistent phenomenon, in which case it could be due to either an azimuthal dynamo wave or an indication of the binary system orbital synchronization still being incomplete. We analyse the datasets using the Carrier Fit method (hereafter CF), especially suitable for analyzing time series in which a fast clocking frequency (such as the rotation of the star) is modulated with a slower process (such as the stellar activity cycle). We combine all collected photometric data into one single data set, and analyze it with the CF method. As a result, we confirm the earlier results of the spot activity having been dominated by one primary spotted region almost through the entire data set, and the existence of a persistent, nearly linear drift. Disruptions of the linear trend and complicated phase behavior are also seen, but the period analysis reveals a rather stable periodicity with P(spot)=6.71054d plus/minus 0.00005d. After the linear trend is removed from the data, we identify several abrupt phase jumps, three of which are analyzed in more detail with the CF method. These phase jumps closely resemble what is called flip-flop event, but the new spot configurations do not, in most cases, persist for longer than a few months.

FK Comae is a rapidly rotating magnetically active star, the light curve of which is modulated by cool spots on its surface. It was the first star where the "flip-flop" phenomenon was discovered. Since then, flip-flops in the spot activity have been reported in many other stars. Therefore, it is of interest to perform a more thorough study of the evolution of the spot activity in FK Com. In this study, we analyse 15 years of photometric observations with two different time series analysis methods, with a special emphasis on detecting flip-flop type events from the data. We apply the continuous period search and carrier fit methods on long-term standard Johnson-Cousins V-observations from the years 1995--2010. The observations were carried out with two automated photometric telescopes, Phoenix-10 and Amadeus T7 located in Arizona. We identify complex phase behaviour in 6 of the 15 analysed data segments. We identify five flip-flop events and two cases of phase jumps, where the phase shift is \Delta \phi < 0.4. In addition we see two mergers of spot regions and two cases where the apparent phase shifts are caused by spot regions drifting with respect to each other. Furthermore we detect variations in the rotation period corresponding to a differential rotation coefficient of |k|>0.031. The flip-flop cannot be interpreted as a single phenomenon, where the main activity jumps from one active longitude to another. In some of our cases the phase shifts can be explained by differential rotation: Two spot regions move with different angular velocity and even pass each other. Comparison between the methods show that the carrier fit utility is better in retrieving slow evolution especially from a low amplitude light curve, while the continuous period search is more sensitive in case of rapid changes.

We formulate an improved time series analysis method for the analysis of photometry of active stars. This new Continuous Period Search (CPS) method is applied to 12 years of V band photometry of the young solar analogue HD 116956 (NQ UMa). The new method is developed from the previous Three Stage Period Analysis (TSPA) method. Our improvements are the use of a sliding window in choosing the modelled datasets, a criterion applied to select the best model for each dataset and the computation of the time scale of change of the light curve. We test the performance of CPS with simulated and real data. The CPS has a much improved time resolution which allows us to better investigate fast evolution of stellar light curves. We can also separate between the cases when the data is best described by periodic and aperiodic (e.g. constant brightness) models. We find, however, that the performance of the CPS has certain limitations. It does not determine the correct model complexity in all cases, especially when the underlying light curve is constant and the number of observations too small. Also the sensitivity in detecting two close light curve minima is limited and it has a certain amount of intrinsic instability in its period estimation. Using the CPS, we find persistent active longitudes in the star HD 116956 and a "flip-flop" event that occurred during the year 1999. Assuming that the surface differential rotation of the star causes observable period variations in the stellar light curve, we determine the differential rotation coefficient to be |k|>0.11. The mean timescale of change of the light curve during the whole 12 year observing period was T_C=44.1 d, which is of the same order as the predicted convective turnover time of the star. We also investigate the presence of activity cycles on the star, but do not find any conclusive evidence supporting them.

Aims: In this study, we aim at investigating the properties of the active longitude system of the young solar analogue HD 116956 in detail, especially concentrating on determining the rotation period of the spot-generating mechanism with respect to the photometric rotation period of the star itself. Because the nonparametric approach, like the Kuiper method, can only give the period of active longitudes, we formulate a new method that can determine the parameters the active longitude distribution uniquely. Methods: For this purpose, we have developed an analysis method, based on Bayesian statistics using Markov chain Monte Carlo, presented in this manuscript. One of the advantages of this method is that an estimate of the active longitude system rotation period, as well as the parameters of the shape and location of the active longitudes together with their respective error estimates. This allows us to compare the active longitude and mean photospheric rotation periods of the star. Results: Our analysis confirms previous results of the object having two stable active longitudes with a phase difference of Delta phi = 0.5, the other longitude having dominated over the other one during almost the entire span of the time series. Our method gives the rotation period of the active longitude system P_al = 7.8412 +- 0.00002 d, which is significantly different from the mean photospheric rotation period of the star P_rot = 7.817 +- 0.003 d. Conclusions: Our analysis indicates that the spot-generating mechanism, manifesting itself as a system of two active longitudes, is lagging behind the overall rotation of the star. This behaviour may be interpreted as a nonaxisymmetric dynamo wave propagating in the rotational reference frame of the stellar surface.

We present surface temperature maps for the FK Comae-type star HD 199178 (V1794 Cygni) calculated from high resolution spectra obtained in 1994 and 1995. The spot pattern evolves, but all maps reveal a large cool spot remaining nearly at the same high latitude. The main spot is 1200 -- 1600 K cooler than the mean surface temperature. The observed slightly flat bottomed absorption lines would usually be interpreted as evidence for a large cool polar spot. We argue that antisolar surface differential rotation offers a better explanation for the box like shape of the line profiles. However, we do not find conclusive evidence for antisolar differential rotation and note that there are still other possible explanations for the slightly flat bottomed line profiles.

The most significant periodicities in the terrestrial impact crater record are due to the human-signal: the bias of assigning integer values for the crater ages. This bias seems to have eluded the proponents and opponents of real periodicity in the occurrence of these events, as well as the theorists searching for an extraterrestrial explanation for such periodicity. The human-signal should be seriously considered by scientists in astronomy, geology and paleontology when searching for a connection between terrestrial major comet or asteroid impacts and mass extinctions of species.