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.

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.

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.

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.