• A majority part of magnetic chemically peculiar (mCP) stars of the upper main sequence exhibits strictly periodic light, magnetic, radio, and spectral variations that can be fully explained by the model of a rigidly rotating main-sequence star with persistent surface structures and stable global magnetic field frozen into the body of the star. Nevertheless, there is an inhomogeneous group consisting of a few mCP stars whose rotation periods vary on timescales of decades, while the shapes of their phase curves remain nearly unchanged. Alternations in the rotational period variations, proven in the case of some of them, offer new insight on this theoretically unpredicted phenomenon. We present a novel and generally applicable method of period analysis based on the simultaneous exploitation of all available observational data containing phase information. This phenomenological method can monitor gradual changes in the observed instantaneous period very efficiently and reliably. We present up to date results of the period monitoring of V901~Ori, CU~Vir, $\sigma$\,Ori\,E, and BS~Cir, known to be mCP stars changing their observed periods and discuss the physics of this unusual behaviour. To compare the period behavior of those stars, we treated their data with an orthogonal polynomial model, which was specifically developed for this purpose. We confirmed period variations in all stars and showed that they reflect real changes in the angular velocity of outer layers of the stars, fastened by their global magnetic fields. However, the nature of the observed rotational instabilities has remained elusive up to now. The discussed group of mCP stars is inhomogeneous to such extent that each of the stars may experience a different cause for its period variations.
  • The magnetic chemically peculiar (mCP) stars of the upper main sequence exhibit periodic light, magnetic, radio, and spectroscopic variations that can be adequately explained by a model of a rigidly rotating magnetized star with persistent surface structures. The majority of mCP stars rotate at strictly constant periods. However, there are a few mCP stars whose rotation periods vary on timescales of decades while the shape of their phase curves remains unchanged. In the case of CU Vir and V901 Ori, we have detected cyclic period variations. We demonstrate that the period oscillations of CU Vir may be a consequence of the interaction of the internal magnetic field and differential rotation.
  • This paper presents an extensive overview of known and proposed RR Lyrae stars in binaries. The aim is to revise and extend the list with new Galactic field systems. We utilized maxima timings for eleven RRab type stars with suspicious behaviour from the GEOS database, and determined maxima timings from data of sky surveys and our own observations. This significantly extended the number of suitable maxima timings. We modelled the proposed Light Time Effect (LiTE) in O-C diagrams to determine orbital parameters for these systems. In contrast to recent studies, our analysis focused on decades-long periods instead of periods in the order of years. Secondary components were found to be predominantly low-mass objects. However, for RZ Cet and AT Ser the mass of the suspected companion of more than one solar mass suggests that it is a massive white dwarf, a neutron star or even a black hole. We found that the semi-major axes of the proposed orbits are between 1 and 20 au. Because the studied stars belong to the closest RR Lyraes, maximal angular distances between components during orbit should at least be between 1 and 13 mas and this improves the chance to detect both stars using current telescopes. However, our interpretation of the O-C diagrams as a consequence of the LiTE should be considered as preliminary without reliable spectroscopic measurements. On the other hand our models give a prediction of the period and radial velocity evolution which should be sufficient for plausible proof of binarity.
  • We present a simple phenomenological method for detrending of raw Kepler and MOST photometry, which is illustrated by means of photometric data processing of two periodically variable chemically peculiar stars, KIC 6950556 and HD 37633. In principle, this method may be applied to any type of periodically variable objects and satellite or ground based photometries. As a by product, we have identified KIC 6950556 as a magnetic chemically peculiar star with an ACV type variability.
  • We present an overview and current status of research on RR Lyrae stars in binary systems. In present days the number of binary candidates has steeply increased and suggested that multiple stellar systems with an RR Lyrae component is much higher than previously thought. We discuss the probability of their detection using various observing methods, compare recent results regarding selection effects, period distribution, proposed orbital parameters and the Blazhko effect.
  • We investigate in detail the LiTE of the most probable binary candidate TU UMa. We used a new code applied on previously available and newly determined maxima timings to confirm binarity and refine parameters of the orbit of the RRab component in the binary system. The binary hypothesis is also tested using RV measurements. We used new approach to determine brightness maxima timings based on template fitting. This can also be used on sparse or scattered data. This approach was successfully applied on measurements from different sources. To determine the orbital parameters of the double star TU UMa, we developed a new code to analyse LiTE that also includes secular variation in the pulsation period. Its usability was successfully tested on CL Aur, an eclipsing binary with mass-transfer in a triple system that shows similar changes in the O-C diagram. Since orbital motion would cause systematic shifts in mean RVs (dominated by pulsations), we computed and compared our model with centre-of-mass velocities. They were determined using high-quality templates of RV curves of RRab stars. Maxima timings adopted from the GEOS database (168) together with those newly determined from sky surveys and new measurements (85) were used to construct an O-C diagram spanning almost five proposed orbital cycles. This data set is three times larger than data sets used by previous authors. Modelling of the O-C dependence resulted in 23.3-year orbital period, which translates into a minimum mass of the second component of about 0.33 M_s. Secular changes in the pulsation period of TU UMa over the whole O-C diagram were satisfactorily approximated by a parabolic trend with a rate of -2.2 ms/yr. To confirm binarity, we used RV measurements from nine independent sources. Although our results are convincing, additional long-term monitoring is necessary to unambiguously confirm the binarity of TU UMa.
  • Context: The presence of heavier chemical elements in stellar atmospheres influences the spectral energy distribution (SED) of stars. An uneven surface distribution of these elements, together with flux redistribution and stellar rotation, are commonly believed to be the primary causes of the variability of chemically peculiar (CP) stars. Aims: We aim to model the photometric variability of the CP star PHI Dra based on the assumption of inhomogeneous surface distribution of heavier elements and compare it to the observed variability of the star. We also intend to identify the processes that contribute most significantly to its photometric variability. Methods: We use a grid of TLUSTY model atmospheres and the SYNSPEC code to model the radiative flux emerging from the individual surface elements of PHI Dra with different chemical compositions. We integrate the emerging flux over the visible surface of the star at different phases throughout the entire rotational period to synthesise theoretical light curves of the star in several spectral bands. Results: The synthetic light curves in the visible and in the near-UV regions are in very good agreement with the observed variability of the star. The lack of usable far-UV measurements of the star precludes making any conclusions about the correctness of our model in this spectral region. We also obtained 194 new BVRI observations of PHI Dra and improved its rotational period to P=1.716500(2). Conclusions: We show that the inhomogeneous distribution of elements, flux redistribution, and rotation of the star are fully capable of explaining the stellar variability in the visible and the near-UV regions. The flux redistribution is mainly caused by bound--free transitions of silicon and bound--bound transitions of iron.
  • The observed light curves of most eclipsing binaries and stars with transiting planets can be well described and interpreted by current advanced physical models which also allow for the determination of many physical parameters of eclipsing systems. However, for several common practical tasks there is no need to know the detailed physics of a variable star, but only the shapes of their light curves or other phase curves. We present a set of phenomenological models for the light curves of eclipsing systems. We express the observed light curves of eclipsing binaries and stars, transited by their exoplanets orbiting in circular trajectories, by a sum of special, analytical, few-parameter functions that enable fitting their light curves with an accuracy of better than 1%. The proposed set of phenomenological models of eclipsing variable light curves were then tested on several real systems. For XY Bootis, we also compare in details the results obtained using our phenomenological modelling with those found using available physical models. We demonstrate that the proposed phenomenological models of transiting exoplanet and eclipsing binary light curves applied to ground-based photometric observations yields results compatible with those obtained by the application of more complex physical models. The suggested phenomenological modelling appears useful to solve a number of common tasks in the field of eclipsing variable research.
  • BS Cir is a representative of moderately cool magnetic chemically peculiar stars which displays very strong light variations in Stroemgren index c1 indicating large changes in the height of the Balmer jump. We present two-spot model of light variations fitting successfully all of nine light curves obtained in the spectral region 335-750 nm. We also discuss the nature of the observed variations of intensities of Fe, Cr, Ti, Si, Mg and RE spectral lines and possible mechanisms matching the observed light variations. It was confirmed that the observed period of BS Cir 2.204 d is rising with the rate of dP/dt=5.4(4)x10^-9. The found minor secular changes in the shape of light curve should be compatible with the period changes caused by precessional motion due to magnetic distortion of the star.
  • Our paper presents new methods for finding and testing of weak periodic variability of stellar objects developed for the purpose of detecting expected regular light variations of magnetic chemically peculiar (mCP) candidates in the Large Magellanic Cloud. We introduce two new periodograms of the mCP star, BS Cir (HD 125630), appropriate for rotating spotted variables and compare the results with those obtained by the well-known Lomb-Scargle periodogram. The usage of periodograms and the testing of the significance of the found period candidates are demonstrated with two examples: the observed and simulated observations of the magnetic field of the mCP star CQ UMa (HD 119213) and the mCP candidate OGLE LMC136.7 16501. Three newly developed tests of the periodic variability - the shuffling, bootstrap and subsidiary ones, are presented. We demonstrate that the found periodic variations known with Signal-to-Noise ratio larger than 6 can be approved as real.
  • We present a new few-parameter phenomenological model of light curves of eclipsing binaries and stars with transiting planets that is able to fit the observed light curves with the accuracy better than 1\% of their amplitudes. The model can be used namely for appropriate descriptions of light curve shapes, classification, mid-eclipse time determination, and fine period analyses.
  • The trustworthiness of orbital period analyses of eclipsing binaries strictly depends on the correctness of the observed mid-eclipse time determination, as well the reliability of its uncertainty estimation. The majority of them has been determined by means of the Kwee-van Woerden method (KWM). There are also other possibilities: to use physical models of eclipsing binaries or light curve templates and to determine mid-eclipse times using the least square method (LSM). We compared results yielded by both methods by means of a computer simulation on the synthetic model of AR\,Aur primary minimum. Minima times determined by the KWM and the exact LSM approach are nearly the same, while the scatter of LSM times is always smaller than the scatter of KWM times. KWM uncertainties are systematically underestimated. We think that the time is ripe for the Kwee-van Woerden method to retire.
  • In-depth analysis of eclipsing binary (EB) observational data collected for several decades can inform us about a lot of astrophysically interesting processes taking place in the systems. We have developed a wide-ranging method for the phenomenological modelling of eclipsing binary phase curves that enables us to combine even very disparate sources of phase information. This approach is appropriate for the processing of both standard photometric series of eclipses and data from photometric surveys of all kind. We conclude that mid-eclipse times, determined using the latest version of our 'hi-fi' phenomenological light curve models, as well as their accuracy, are nearly the same as the values obtained using much more complex standard physical EB models.
  • The majority of magnetic chemically peculiar (mCP) stars exhibit periodic light, magnetic, radio, and spectroscopic variations that can be adequately modelled as a rigidly-rotating main-sequence star with persistent surface structures. Nevertheless, there is a small sample of diverse mCP stars whose rotation periods vary on timescales of decades while the shapes of their phase curves remain unchanged. Alternating period increases and decreases have been suspected in the hot CP stars CU Vir and V901 Ori, while rotation in the moderately cool star BS Cir has been decelerating. These examples bring new insight into this theoretically unpredicted phenomenon. We discuss possible causes of such behaviour and propose that dynamic interactions between a thin, outer, magnetically-confined envelope braked by the stellar wind, and an inner faster-rotating stellar body are able to explain the observed rotational variability
  • Plenty of mid-eclipse timings of short-periodic eclipsing binaries derived from series of visual observations appear to be an acceptable source of information for orbital period analyses, namely if they were done in time-intervals not covered by other types of observations. However, our thorough period analysis of the nearly contact eclipsing binary BS Vulpeculae proves that visually determined timings done in 1979--2003 were undoubtedly biased to accommodate the existing linear ephemeris. The heavily subjective character of visual observations disqualifies them as a source of true phase information apt for fine eclipsing binary period analyses. Consequently we warn against the use of visual timings without a preceding careful verification.
  • Models of period variations are basic tools for period analyzes of variable stars. We introduce phase function and instant period and formulate basic relations and equations among them. Some simple period models are also presented.
  • We present the first catalogue of known variable stars in open cluster regions and with up to two times the given cluster radius. This gives basic information about the distribution of variable stars in cluster fields for the complete sky. Knowledge of the variable star contents in open clusters is a significant advantage in their study. Analysing variability of cluster members and fields stars as well, allows us to study the characteristics of stars and clusters together. This catalogue of variable stars in open cluster fields is the first step in supporting such studies. We took all variable and suspected variable stars into account from the most complete collection, "The AAVSO Variable Star Index", and did a cross-match of these stars with the most complete catalogue of galactic open clusters named DAML02. Our on-line catalogue presently contains 18 065 variable stars. We present the basic statistical distribution according to types of variability.
  • Two methods for the determination of extrema timings and their uncertainties appropriate for the analysis of time series of variable stars using matrix calculus are presented. The method I is suitable for determination of times of extrema of non-periodical variables or objects, whose light curves vary. The method II is apt for O-C analyses of objects whose light curves are more or less repeating.
  • CU Vir and V901 Ori belong among these few magnetic chemically peculiar stars whose rotation periods vary on timescales of decades. We aim to study the stability of the periods in CU Vir and V901 Ori using all accessible observational data containing phase information. We collected all available relevant archived observations supplemented with our new measurements of these stars and analysed the period variations of the stars using a novel method that allows for the combination of data of diverse sorts. We found that the shapes of their phase curves were constant, while the periods were changing. Both stars exhibit alternating intervals of rotational braking and acceleration. The rotation period of CU Vir was gradually shortening until the year 1968, when it reached its local minimum of 0.52067198 d. The period then started increasing, reaching its local maximum of 0.5207163 d in the year 2005. Since that time the rotation has begun to accelerate again. We also found much smaller period changes in CU Vir on a timescale of several years. The rotation period of V901 Ori was increasing for the past quarter-century, reaching a maximum of 1.538771 d in the year 2003, when the rotation period began to decrease. A theoretically unexpected alternating variability of rotation periods in these stars would remove the spin-down time paradox and brings a new insight into structure and evolution of magnetic upper-main-sequence stars.
  • Photometrically and spectroscopically variable chemically peculiar (CP) stars are the optimum laboratories for testing the rotational evolution of main sequence stars. A vast majority of well-studied CP stars have quite steady rotational periods (e.g. SrCrEu star CQ UMa). However, there are several CP stars that exhibit apparent period variations. The origin of the period variations is unclear in many cases. We describe the observed period variations of several individual CP stars, especially V901 Ori, sigma Ori E, HR 7355, CU Vir, SX Ari, and EE Dra. The CP stars with unsteady periods now represent a very diverse group with dissimilar O-C diagrams and time scales. We also discuss the causes of the period changes found and a possible cyclicity or chaoticity of them.
  • We present two diagnostic methods based on ideas of Principal Component Analysis and demonstrate their efficiency for sophisticated processing of multicolour photometric observations of variable objects.
  • Magnetic chemically peculiar (mCP) stars of the upper part of the main sequence usually exhibit strictly periodic variations of their light, spectrum and magnetic field. These changes can be well explained using the oblique rotator model of the rigidly rotating star with persistent photometric spots, uneven horizontal abundance of chemical elements and dipole-like magnetic field. Long-term observations of mCP stars variability enable to determine their rotational periods with an extraordinary accuracy and reliability. We compare rotational periods of mCP stars with those of normal main sequence stars of the same spectral type and discuss the found discrepancies. We also discuss the cases of mCP stars V901 Ori, SX Ari, CU Vir, (and possibly also EE Dra), which have been reported to display an increase of their rotational periods.
  • We calculate the uvby light curves of the helium-strong chemically peculiar (CP) star HD 37776. Assuming that the chemical peculiarity influences the monochromatic radiative flux due to mainly the bound-free processes and using the model of the surface distribution of helium and silicon on HD 37776, as derived from spectroscopy, we calculate the emergent fluxes for a series of its rotational phases. We show that it is possible to consistently simulate light curves of a chemically peculiar star and that the basic properties of variability of helium CP stars can be understood in terms of the model of spots with a peculiar chemical composition.
  • We present an analysis of light variations in UBVRI of the anomalous object in the center of planetary nebula Sh 2-71. We refined the linear ephemeris of the light maxima to JD{max} = 2449862.0 + 68.101 (E-96), but also identified long-term, obviously non-periodic variations. The latter manifest themselves in large O-C shifts, a variable profile of light curves (hereafter LC) and changes in the mean brightness of the object. Our spectroscopic observations suggested the presence of a superdense nebula in the center of Sh 2-71.
  • We present a principle component analysis of multicolour light variations of the He-strong magnetic chemically peculiar star HD 37776 and its (O-C) diagram. The period increase rate found, (dP/dt) / P = 5.5 x 10^{-6}/year, is consistent with the conception of angular momentum loss through magnetically confined stellar wind.