The magnetic activity of late-type stars is correlated with their rotation rates. Up to a certain limit, stars with smaller Rossby numbers, defined as the rotation period divided by the convective turnover time, have higher activity. A more detailed look at this rotation-activity relation reveals that, rather than being a simple power law relation, the activity scaling has a shallower slope for the low-Rossby stars than for the high-Rossby ones. We find that, for the chromospheric CaII H&K activity, this scaling relation is well modelled by a broken two-piece power law. Furthermore, the knee-point of the relation coincides with the axisymmetry to non-axisymmetry transition seen in both the spot activity and surface magnetic field configuration of active stars. We interpret this knee-point as a dynamo transition between dominating axi- and non-axisymmetric dynamo regimes with a different dependence on rotation and discuss this hypothesis in the light of current numerical dynamo models.

Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other stellar parameters. We also study the differences in the shape of the solar cycle, depending on what activity indicator is used. The observations are also compared to simulated activity cycles. We use the chromospheric Ca II H&K data from the Mount Wilson Observatory HK Project. From this data set we identify 47 individual cycles from 18 stars. We use the statistical skewness of a cycle as a measure of its asymmetry, and compare this to other stellar parameters. A similar analysis has been done to magnetic cycles extracted from direct numerical magnetohydrodynamic simulations of solar-type convection zones. The shape of the solar cycle (fast rise and slower decline) is common in other stars as well, although the Sun has particularly asymmetric cycles. Cycle-to-cycle variations are large, but the average shape of a cycle is still fairly well represented by a sinusoid. We find only slight correlations between the cycle asymmetry and other stellar parameters. There are large differences in the shape of the solar cycle, depending on what activity indicator is used. In the simulated cycles, there is a difference in the symmetry of global simulations covering the full longitudinal range, hence capable of exciting non-axisymmetric dynamo modes, versus wedge simulations covering a partial extent in longitude, where only axisymmetric modes are possible. The former produce preferentially positive skewness, while the latter a negative one.

One interpretation of the activity and magnetism of late-type stars is that these both intensify with decreasing Rossby number up to a saturation level, suggesting that stellar dynamos depend on both rotation and convective turbulence. Some studies have claimed, however, that rotation alone suffices to parametrise this scaling adequately. Here, we tackle the question of the relevance of turbulence to stellar dynamos by including evolved, post main sequence stars in the analysis of the rotation-activity relation. These stars rotate very slowly compared with main sequence stars, but exhibit similar activity levels. We show that the two evolutionary stages fall together in the rotation-activity diagram and form a single sequence in the unsaturated regime in relation only to Rossby numbers derived from stellar models, confirming earlier preliminary results that relied on a more simplistic parametrisation of the convective turnover time. This mirrors recent results of fully convective M dwarfs, which likewise fall on the same rotation-activity sequence as partially convective solar-type stars. Our results demonstrate that turbulence plays a crucial role in driving stellar dynamos and suggest that there is a common turbulence-related dynamo mechanism explaining the magnetic activity of all late-type stars.

Global magnetic fields of active solar-like stars are nowadays routinely detected with spectropolarimetric measurements and are mapped with Zeeman-Doppler imaging (ZDI). However, due to the cancellation of opposite field polarities, polarimetry captures only a tiny fraction of the magnetic flux and cannot assess the overall stellar surface magnetic field if it is dominated by a small-scale component. Analysis of Zeeman broadening in high-resolution intensity spectra can reveal these hidden complex magnetic fields. Historically, there were very few attempts to obtain such measurements for G dwarf stars due to the difficulty of disentangling Zeeman effect from other broadening mechanisms affecting spectral lines. Here we developed a new magnetic field diagnostic method based on relative Zeeman intensification of optical atomic lines with different magnetic sensitivity. Using this technique we obtained 78 field strength measurements for 15 Sun-like stars, including some of the best-studied young solar twins. We find that the average magnetic field strength $Bf$ drops from 1.3-2.0 kG in stars younger than about 120 Myr to 0.2-0.8 kG in older stars. The mean field strength shows a clear correlation with the Rossby number and with the coronal and chromospheric emission indicators. Our results suggest that magnetic regions have roughly the same local field strength $B\approx3.2$ kG in all stars, with the filling factor $f$ of these regions systematically increasing with stellar activity. Comparing our results with the spectropolarimetric analyses of global magnetic fields in the same stars, we find that ZDI recovers about 1% of the total magnetic field energy in the most active stars. This figure drops to just 0.01% for the least active targets.

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.

Context. Studying the spots of late-type stars is crucial for distinguishing between the various proposed dynamo mechanisms believed to be the main cause of starspot activity. For this research it is important to collect observation time series that are long enough to unravel both long- and short-term spot evolution. Doppler imaging is a very efficient method for studying spots of stars that cannot be angularly resolved. Aims. High-resolution spectral observations during 1994-2017 are analysed in order to reveal long- and short-term changes in the spot activity of the FK Comae-type subgiant HD 199178. Methods. Most of the observations were collected with the Nordic Optical Telescope. The Doppler imaging temperature maps were calculated using an inversion technique based on Tikhonov regularisation and utilising multiple spectral lines. Results. We present a unique series of 41 temperature maps spanning more than 23 years. All reliable images show a large cool spot region centred near the visible rotation pole. Some lower latitude cool features are also recovered, although the reliability of these is questionable. There is an expected anti-correlation between the mean surface temperature and the spot coverage. Using the Doppler images, we construct the equivalent of a solar butterfly diagram for HD 199178. Conclusions. HD 199178 clearly has a long-term large and cool spot structure at the rotational pole. This spot structure dominated the spot activity during the years 1994-2017. The size and position of the structure has evolved with time, with a gradual increase during the last years. The lack of lower latitude features prevents the determination of a possible differential rotation.

Context. LQ Hya is one of the most frequently studied young solar analogue stars. Recently, it has been observed to show intriguing behaviour by analysing long-term photometry: during 2003--2009, a coherent spot structure migrating in the rotational frame has been reported by various authors, but since that time the star has entered a chaotic state where coherent structures seem to have disappeared and rapid phase jumps of the photometric minima occur irregularly over time. Aims. LQ Hya is one of the stars included in the SOFIN/FIES long-term monitoring campaign extending over 25 years. Here we publish new temperature maps for the star during 2006--2017, covering the chaotic state of the star. Methods. We use a Doppler imaging technique to derive surface temperature maps from high-resolution spectra. Results. From the mean temperatures of the Doppler maps we see a weak but systematic increase in the surface temperature of the star. This is consistent with the simultaneously increasing photometric magnitude. During nearly all observing seasons we see a high-latitude spot structure which is clearly nonaxisymmetric. The phase behaviour of this structure is very chaotic but agrees reasonably well with the photometry. Equatorial spots are also frequently seen, but many of them we interpret to be artefacts due to the poor to moderate phase coverage. Conclusions. Even during the chaotic phase of the star, the spot topology has remained very similar to the higher activity epochs with more coherent and long-lived spot structures: we see high-latitude and equatorial spot activity, the mid latitude range still being most often void of spots. We interpret the erratic jumps and drifts in phase of the photometric minima to be caused by changes in the high-latitude spot structure rather than the equatorial spots.

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.

Debate over the existence of branches in the stellar activity-rotation diagrams continues. Application of modern time series analysis tools to study the mean cycle periods in chromospheric activity index is lacking. We develop such models, based on Gaussian processes, for one-dimensional time series and apply it to the extended Mount Wilson Ca H&K sample. Our main aim is to study how the previously commonly used assumption of strict harmonicity of the stellar cycles as well as handling of the linear trends affects the results. We introduce three methods of different complexity, starting with the simple Bayesian harmonic model and followed by Gaussian Process models with periodic and quasi-periodic covariance functions. We confirm the existence of two populations in the activity-period diagram. We find only one significant trend in the inactive population, namely that the cycle periods get shorter with increasing rotation. This is in contrast with earlier studies, that postulate the existence of trends in both of the populations. In terms of rotation to cycle period ratio, our data is consistent with only two activity branches such that the active branch merges together with the transitional one. The retrieved stellar cycles are uniformly distributed over the R'HK activity index, indicating that the operation of stellar large-scale dynamos carries smoothly over the Vaughan-Preston gap. At around the solar activity index, however, indications of a disruption in the cyclic dynamo action are seen. Our study shows that stellar cycle estimates depend significantly on the model applied. Such model-dependent aspects include the improper treatment of linear trends, while the assumption of strict harmonicity can result in the appearance of double cyclicities that seem more likely to be explained by the quasi-periodicity of the cycles.

We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition coincides with a change in the rotation profile from antisolar- to solar-like differential rotation with a faster equator and slow poles. In the solar-like rotation regime, the field configuration consists of an axisymmetric oscillatory field accompanied by an m=1 azimuthal mode (two active longitudes), which also shows temporal variability. At slow (rapid) rotation, the axisymmetric (nonaxisymmetric) mode dominates. The axisymmetric mode produces latitudinal dynamo waves with polarity reversals, while the nonaxisymmetric mode often exhibits a drift in the rotating reference frame and the strength of the active longitudes changes cyclically over time between the different hemispheres. Most of the obtained dynamo solutions exhibit cyclic variability either caused by latitudinal or azimuthal dynamo waves. In an activity-period diagram, the cycle lengths normalized by the rotation period form two different populations as a function of rotation rate or magnetic activity level. The slowly rotating axisymmetric population lies close to what is called the inactive branch in observations, while the rapidly rotating models are close to the superactive branch with a declining cycle to rotation frequency ratio with increasing rotation rate. We can successfully reproduce the transition from axi- to nonaxisymmetric dynamo solutions for high rotation rates, but high-resolution simulations are required to limit the effect of rotational quenching of convection at rotation rates above 20 times the solar value.

Period estimation is one of the central topics in astronomical time series analysis, where data is often unevenly sampled. Especially challenging are studies of stellar magnetic cycles, as there the periods looked for are of the order of the same length than the datasets themselves. The datasets often contain trends, the origin of which is either a real long-term cycle or an instrumental effect, but these effects cannot be reliably separated, while they can lead to erroneous period determinations if not properly handled. In this study we aim at developing a method that can handle the trends properly, and by performing extensive set of testing, we show that this is the optimal procedure when contrasted with methods that do not include the trend directly to the model. The effect of the form of the noise (whether constant or heteroscedastic) on the results is also investigated. We introduce a Bayesian Generalised Lomb-Scargle Periodogram with Trend (BGLST), which is a probabilistic linear regression model using Gaussian priors for the coefficients and uniform prior for the frequency parameter. We show, using synthetic data, that when there is no prior information on whether and to what extent the true model of the data contains a linear trend, the introduced BGLST method is preferable to the methods which either detrend the data or leave the data untrended before fitting the periodic model. Whether to use noise with different than constant variance in the model depends on the density of the data sampling as well as on the true noise type of the process.

For flares to be generated, stars have to have a sufficiently deep outer convection zone (F5 and later), strong large--scale magnetic fields (Ap/Bp-type stars) or strong, radiatively driven winds (B5 and earlier). Normal A-type stars possess none of these and therefore should not flare. Nevertheless, flares have previously been detected in the Kepler lightcurves of 33 A-type stars and interpreted to be intrinsic to the stars. Here we present new and detailed analyses of these 33 stars, imposing very strict criteria for the flare detection. We confirm the presence of flare-like features in 27 of the 33 A-type stars. A study of the pixel data and the surrounding field-of-view (FOV) reveals that 14 of these 27 flaring objects have overlapping neighbouring stars and 5 stars show clear contamination in the pixel data. We have obtained high-resolution spectra for 2/3 of the entire sample and confirm that our targets are indeed A-type stars. Detailed analyses revealed that 11 out of 19 stars with multiple epochs of observations are spectroscopic binaries. Furthermore, and contrary to previous studies, we find that the flares can originate from a cooler, unresolved companion. We note the presence of H$\alpha$ emission in eight stars. Whether this emission is circumstellar or magnetic in origin is unknown. In summary, we find possible alternative explanations for the observed flares for at least 19 of the 33 A-type stars, but find no truly convincing target to support the hypothesis of flaring A-type stars.

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.

The goal of this work is to study the magnetic fields of six young solar-analogue stars both individually, and collectively, to search for possible magnetic field trends with age. If such trends are found, they can be used to understand magnetism in the context of stellar evolution of solar-like stars and, the past of the Sun and the solar system. This is also important for the atmospheric evolution of the inner planets, Earth in particular. We used Stokes IV data from two different spectropolarimeters, NARVAL and HARPSpol. The least-squares deconvolution multi-line technique was used to increase the signal-to-noise ratio of the data. We then applied a modern Zeeman-Doppler imaging code in order to reconstruct the magnetic topology of all stars and the brightness distribution of one of our studied stars. Our results show a significant decrease in the magnetic field strength and energy as the stellar age increases from 100Myr to 250Myr while there is no significant age dependence of the mean magnetic field strength for stars with ages 250-650Myr. The spread in the mean field strength between different stars is comparable to the scatter between different observations of individual stars. The meridional field component has the weakest strength compared to the radial and azimuthal field components in 15 out of the 16 magnetic maps. It turns out that 89-97% of the magnetic field energy is contained in l=1-3. There is also no clear trend with age and distribution of field energy into poloidal/toroidal and axisymmetric/non-axisymmetric components within the sample. The two oldest stars in this study do show a twice as strong octupole component compared to the quadrupole component. This is only seen in one out of 13 maps of the younger stars. One star, chi1 Ori displays two field polarity switches during almost 5 years of observations suggesting a magnetic cycle length of either 2, 6 or 8 years.

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.

By studying young magnetically active late-type stars, i.e. analogues to the young Sun, one can draw conclusions on the evolution of the solar dynamo. We determine the topology of the surface magnetic field and study the relation between the magnetic field and cool photospheric spots in three young late-type stars. High-resolution spectropolarimetry of the targets were obtained with the HARPSpol instrument mounted at the ESO 3.6 m telescope. The signal-to-noise ratio of the Stokes IV measurements were boosted by combining the signal from a large number of spectroscopic absorption lines through the least squares deconvolution technique. Surface brightness and magnetic field maps were calculated using the Zeeman-Doppler imaging technique. All the three targets show clear signs of both magnetic fields and cool spots. Only one of the targets, namely V1358 Ori, shows evidence of the dominance of non-axisymmetric modes. In two of the targets, the poloidal field is significantly stronger than the toroidal one, indicative of an $\alpha^2$-type of a dynamo, in which convective turbulence effects dominate over the weak differential rotation. In two of the cases there is a slight anti-correlation between the cool spots and the strength of the radial magnetic field. However, even in these cases the correlation is much weaker than in the case of sunspots. The weak correlation between the measured radial magnetic field and cool spots may indicate a more complex magnetic field structure in the spots or spot groups involving mixed magnetic polarities. Comparison with a previously published magnetic field map shows that on one of the stars, HD 29615, the underlying magnetic field has changed its polarity between 2009 and 2013.

Over the duration of the Kepler mission, KIC8462852 was observed to undergo irregularly shaped, aperiodic dips in flux of up to $\sim 20$\%. The dipping activity can last for between 5 and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy distribution fitting, radial velocity measurements, high-resolution imaging, and Fourier analyses of the Kepler light curve. We determine that KIC8462852 is a typical main-sequence F3 V star that exhibits no significant IR excess, and has no very close interacting companions. In this paper, we describe various scenarios to explain the dipping events observed in the Kepler light curve. We confirm that the dipping signals in the data are not caused by any instrumental or data processing artifact, and thus are astrophysical in origin. We construct scenario-independent constraints on the size and location of a body in the system that is needed to reproduce the observations. We deliberate over several assorted stellar and circumstellar astrophysical scenarios, most of which have problems explaining the data in hand. By considering the observational constraints on dust clumps in orbit around a normal main-sequence star, we conclude that the scenario most consistent with the data in hand is the passage of a family of exocomet or planetesimal fragments, all of which are associated with a single previous break-up event, possibly caused by tidal disruption or thermal processing. The minimum total mass associated with these fragments likely exceeds $10^{-6}$~\mearth, corresponding to an original rocky body of $>100$~km in diameter. We discuss the necessity of future observations to help interpret the system.

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 study LQ Hya photometry for 1982-2014 with the carrier fit (CF) -method and compare our results to earlier photometric analysis and recent Doppler imaging maps. We first utilize different types of statistical methods to estimate various candidates for the carrier period for the CF method. Secondly, a global fit to the whole data set and local fits to shorter segments are computed with the period that is found to be the optimal one. The harmonic least-squares analysis of all the available data reveals a short period close to 1.6 days as a limiting value for a set of significant frequencies. We interpret this as the rotation period of the spots near the equatorial region. In addition, the distribution of the significant periods is found to be bimodal, hinting of a longer-term modulating period, which we set out to study with a two-harmonic CF model. Weak modulation signal is, indeed retrieved, with a period of roughly 6.9 years. The phase dispersion analysis gives a clear symmetric minimum for coherence times lower than and around 100 days. We interpret this as the mean rotation period of the spots (1.60514 days), and this value is chosen to be used as the carrier period for the CF analysis. With the CF method we seek for any systematic trends in the spot distribution in the global time frame, and locally look for abrupt phase changes earlier reported in rapidly rotating objects. During 2005-2008 the global CF reveals a coherent structure rotating with a period of 1.6037 days, while during most other times the spot distribution appears rather random in phase. The evolution of the spot distribution of the object is found to be very chaotic, with no clear signs of an azimuthal dynamo wave that would persist over longer time scales, although the short-lived coherent structures observed occasionally do not rotate with the same speed as the mean spot distribution.

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.

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.