The central kiloparsecs of the Milky Way are known to host an old, spheroidal stellar population, whose spatial and kinematical properties set it apart from the boxy/peanut structure that constitutes most of the central stellar mass. The nature of this spheroidal population, whether a small classical bulge, the innermost stellar halo or a population of disk stars with large initial velocity dispersion, remains unclear. This structure is also a promising candidate to host some of the oldest stars in the Galaxy. Here we address the topic of the inner stellar spheroid age, using spectroscopic and photometric metallicities for a sample of 935 RR Lyrae stars that are constituents of this component. By means of stellar population synthesis, we derive an age-metallicity relation for RR Lyrae populations. We infer, for the RR Lyrae stars in the bulge spheroid, an extremely ancient age of $13.41 \pm 0.54$ Gyr and conclude they were among the first stars to form in what is now the Milky Way galaxy. Our age estimate for the central spheroid shows remarkable agreement with the age profile that has been inferred for the Milky Way stellar halo, suggesting a connection between the two structures. However, we find mild evidence for a transition in the halo properties at $r_{\rm GC} \sim 5$~kpc. We also investigate formation scenarios for metal-rich RR Lyrae stars, such as binarity and helium variations, and whether they can provide alternative explanations for the properties of our sample. We conclude that, within our framework, the only viable alternative is to have younger, slightly helium-rich, RR Lyrae stars, a hypothesis that would open intriguing questions for the formation of the inner stellar spheroid.

Using a sample of 81 galactic, detached eclipsing binary stars we investigated the global zero-point shift of their parallaxes with the Gaia Data Release 2 (DR2) parallaxes. The stars in the sample lay in a distance range of 0.04-2 kpc from the Sun. The photometric parallaxes {\pi}_Phot of the eclipsing binaries were determined by applying a number of empirical surface brightness - color (SBC) relations calibrated on optical-infrared colors. For each SBC relation we calculated the individual differences d{\pi}_i = ({\pi}_Gaia - {\pi}_Phot)_i and then we calculated unweighted and weighted means. As the sample covers the whole sky we interpret the weighted means as the global shifts of the Gaia DR2 parallaxes with respect to our eclipsing binary sample. Depending on the choice of the SBC relation the shifts vary from -0.094 mas to -0.025 mas. The weighted mean of the zero-point shift from all colors and calibrations used is d{\pi} = -0.054 +/- 0.024 mas. However, the SBC relations based on (B-K) and (V-K) colors, which are the least reddening dependent and have the lowest intrinsic dispersions, give a zero-point shift of d{\pi} = -0.031 +/- 0.011 mas in full agreement with results obtained by Lindegren et al. and Arenou et al. Our result confirms the global shift of Gaia DR2 parallaxes of d{\pi} = -0.029 mas reported by the Gaia team, but we do not confirm the larger zero-point shift reported by a number of follow-up papers.

We present the most extended and homogeneous study carried out so far of the main and early shocks in 1485 RR~Lyrae stars in the Galactic bulge observed by the Optical Gravitational Lensing Experiment (OGLE). We selected non-modulated fundamental-mode RR~Lyrae stars with good-quality photometry. Using a self-developed method, we determined the centers and strengths of main and early shock features in the phased light curves. We found that the position of both humps and bumps are highly correlated with the pulsation properties of the studied variables. Pulsators with a pronounced main shock are concentrated in the low-amplitude regime of the period-amplitude diagram, while stars with a strong early shock have average and above-average pulsation amplitudes. A connection between the main and early shocks and the Fourier coefficients is also observed. In the color-magnitude diagram (CMD), we see a separation between stars with strong and weak shocks. Variables with a pronounced main shock cluster close to the fundamental red edge of the instability strip (IS), while stars with a strong early shock tend to clump in the center and near the fundamental blue edge of the IS. The appearance of shocks and their properties seem independent of the direction of evolution estimated from the period change rate of the studied stars. In addition, the differences in the period change rate between the two main Oosterhoff groups found in the Galactic bulge suggest that stars of Oosterhoff type I are located close to the zero-age horizontal branch while Oosterhoff type II variables are on their way toward the fundamental red edge of the instability strip, thus having already left the zero-age horizontal branch.

We present candidates for non-pulsating stars lying in the classical Cepheid instability strip based on OGLE photometric maps combined with Str\"omgren photometry obtained with the 4.1-m SOAR telescope, and Gaia DR2 data in four fields in the Large Magellanic Cloud. We selected 19 candidates in total. After analysis of their light curves from OGLE surveys we found that all these stars appear to be photometrically stable at the level of a few mmag. Our results show that non-pulsating stars might constitute to about 21%-30% of the whole sample of giant stars located in the classical instability strip. Furthermore, we identified potential candidates for classical Cepheids with hot companions based on their Str\"omgren colours.

We analyzed photometry for the up-to-date collection of the first-overtone RR Lyrae stars (RRc; 11415 stars) and double-mode RR Lyrae stars (RRd; 148 stars) towards the Galactic bulge from the Optical Gravitational Lensing Experiment. We analyzed frequency spectra of these stars in search for additional, low-amplitude signals, beyond the radial modes. We focused on stars from two groups: $RR_{0.61}$ and $RR_{0.68}$. In the first group, additional low-amplitude signals have periods shorter than the first-overtone period; period ratios fall in the 0.60-0.64 range. In the second group, additional low-amplitude signals have periods longer than the first-overtone period; period ratios tightly cluster around 0.68. Altogether we have detected 960 and 147 RR Lyrae stars that belong to $RR_{0.61}$ and $RR_{0.68}$ groups, respectively, which yield the incidence rates of 8.3 and 1.3 per cent of the considered sample. We discuss statistical properties of RR Lyrae stars with additional periodicities. For $RR_{0.61}$ group we provide strong arguments that additional periodicities are connected to non-radial pulsation modes of degrees $\ell=8$ and $\ell=9$, as proposed by Dziembowski. We have also detected two double-periodic variables, with two close periodicities, similar to RR Lyrae variable V37 in NGC 6362. Properties of these peculiar variables, which may form a new group of double-mode pulsators, are discussed.

$\varepsilon$ Lupi A is a binary system consisting of two main sequence early B-type stars Aa and Ab in a short period, moderately eccentric orbit. The close binary pair is the only doubly-magnetic massive binary currently known. Using photometric data from the BRITE-Constellation we identify a modest heartbeat variation. Combining the photometry with radial velocities of both components we determine a full orbital solution including empirical masses and radii. These results are compared with stellar evolution models as well as interferometry and the differences discussed. We also find additional photometric variability at several frequencies, finding it unlikely these frequencies can be caused by tidally excited oscillations. We do, however, determine that these signals are consistent with gravity mode pulsations typical for slowly pulsating B stars. Finally we discuss how the evolution of this system will be affected by magnetism, determining that tidal interactions will still be dominant.

We present a kinematical study of RR~Lyrae stars associated with two Oosterhoff groups in the Galactic bulge. We used data published in the first paper of the series, plus proper motions from the {\it Gaia} Data Release 2, and radial velocities from the literature. A 6D kinematical and spatial solution was obtained for 429 RR~Lyrae stars. We use a model of the Galactic gravitational potential to infer stellar orbits. We did not find a difference between the Oosterhoff groups in the individual components of the space velocity. We report that \textit{foreground} and \textit{background} stars with respect to the Galactic bulge stand out in the mean $V$ velocity component, which we interpret as a sign of the Galactic rotation. The movement of the studied stars in the central region of the Galactic bulge is consistent with random motions expected for a classical bulge component. From the orbital integration, we estimate that 8\,\% of the RR~Lyrae stars are halo interlopers currently located in the Galactic bulge. The majority of the stars' orbits are within a 3\,kpc radius from the Galactic bulge. The fraction of Oosterhoff\,II stars increases with increasing Galactic latitude, as well as towards longer orbital periods. We found several RR~Lyrae stars with high space velocities, one of which has an extremely long orbital period of $\sim$1\,Gyr. We conclude that based on their kinematics, the vast majority of the stars in our sample do not seem to contribute to the boxy/peanut component of the Galactic bulge.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

Aim: Our aim is to obtain high-accuracy measurements of the physical and orbital parameters of two evolved eclipsing binary systems, and to use these measurements to study their evolutionary status. We also aim to derive distances to the systems by using a surface brightness - colour relation and compare these distances with the measurements provided by GAIA. Methods: We measured the physical and orbital parameters on both systems based on V-band and I-band photometry from OGLE, near-infrared photometry obtained with the NTT telescope and the instrument SOFI, as well as high-resolution spectra obtained at ESO 3.6m/HARPS and Clay 6.5/MIKE spectrographs. The light curves and radial-velocity curves were analysed with the Wilson-Devinney code. Results: We analysed two double-lined eclipsing binary systems OGLE-BLG-ECL-23903 and OGLE-BLG-ECL-296596 from the Optical Gravitational Lensing Experiment (OGLE) catalogue. Both systems have a configuration of two well-detached giant stars. The masses of the components ofOGLE-BLG-ECL-123903 are M_1= 2.045 $\pm$ 0.027 and M_2=2.074 $\pm$ 0.023 $M_\odot$ and the radii are R_1=9.540 $\pm$ 0.049 and R_2=9.052 $\pm$ 0.060 $R_\odot$. For OGLE-BLG-ECL-296596, the masses are M_1=1.093 $\pm$ 0.015 and M_2=1.125 $\pm$ 0.014 $M_\odot$, while the radii are R_1=18.06 $\pm$ 0.28 and R_2=29.80 $\pm$ 0.33 $R_\odot$. Evolutionary status was discussed based on the isochrones and evolutionary tracks from PARSEC and MESA codes. The ages of the systems were establishes to be around 1.3 Gyr for the OGLE-BLG-ECL-123903 and 7.7 Gyr for the OGLE-BLG-ECL-296596. We also determined the distance to both systems. For OGLE-BLG-ECL-123903 this is equal to d=2.95 $\pm$ 0.06 (stat.) $\pm$ 0.07 (syst.) kpc, while for the OGLE-BLG-ECL-296596 it is d=5.68 $\pm$ 0.07 (stat.) $\pm$ 0.14 (syst.) kpc. This is the first analysis of its kind for these unique evolved eclipsing binary systems.

In the era of precision cosmology, it is essential to empirically determine the Hubble constant with an accuracy of one per cent or better. At present, the uncertainty on this constant is dominated by the uncertainty in the calibration of the Cepheid period - luminosity relationship (also known as Leavitt Law). The Large Magellanic Cloud has traditionally served as the best galaxy with which to calibrate Cepheid period-luminosity relations, and as a result has become the best anchor point for the cosmic distance scale. Eclipsing binary systems composed of late-type stars offer the most precise and accurate way to measure the distance to the Large Magellanic Cloud. Currently the limit of the precision attainable with this technique is about two per cent, and is set by the precision of the existing calibrations of the surface brightness - colour relation. Here we report the calibration of the surface brightness-colour relation with a precision of 0.8 per cent. We use this calibration to determine the geometrical distance to the Large Magellanic Cloud that is precise to 1 per cent based on 20 eclipsing binary systems. The final distane is 49.59 +/- 0.09 (statistical) +/- 0.54 (systematic) kiloparsecs.

We report the discovery of the first type II Cepheids (BL Herculis stars) pulsating solely in the first overtone. We found two such objects among tens of millions of stars regularly observed by the OGLE survey in the Large Magellanic Cloud. Our classification and the pulsation mode identification is based on the position of these stars on the period-luminosity and color-magnitude diagrams and on the light curve analysis. We discuss why single-mode first-overtone BL Her pulsators must be very rare. For the two discovered stars we present non-linear models that successfully reproduce their light variation. These models indicate that both first-overtone pulsators should be more massive than it is typically assumed for BL Her stars, i.e. their masses should be above 0.75 Solar masses. However, the higher mass requires higher luminosity to match the observed periods of the stars, which is inconsistent with observations.

We present a study of the Oosterhoff (Oo) dichotomy in the Galactic bulge using 8\,141 fundamental mode RR~Lyrae stars. We used public photometric data from the Optical Gravitational Lensing Experiment (OGLE) and the Vista Variables in the V\'ia L\'actea survey (VVV). We carefully selected fundamental mode stars without modulation and without association with any globular cluster located toward the Galactic bulge. Subsequently, we identified and separated the Oosterhoff groups I and II on the basis of their period-amplitude distribution and using a relation fitted to the Oosterhoff I locus. Both Oosterhoff groups were then compared to observations of two bulge globular clusters and with models of stellar pulsation and evolution. We found that some of the variables classified as Oo\,II belong to a third Oo group. The Oosterhoff II variables are more metal-poor on average, more massive, and cooler than their Oosterhoff I counterparts. The analysis of their spatial distribution shows a systematic difference between \textit{foreground}, central and \textit{background} regions in the occurrence of the Oosterhoff II group. The difference between the Oo\,I and II groups is also seen in their distance distributions with respect to the Galactic bar, but neither group is associated with the bar.

We present the analysis of the Blazhko effect - quasi-periodic modulation of pulsation amplitude and/or phase - in the Galactic bulge first overtone RR Lyrae stars (RRc). We used the data gathered during the fourth phase of the Optical Gravitational Lensing Experiment (OGLE). Out of 10 826 analyzed RRc stars, Blazhko effect was detected in 607 stars which constitute 5.6 percent of the sample. It is the largest and most homogeneous sample of modulated RRc stars analyzed so far. Modulation periods cover a wide range, from slightly above 2 d to nearly 3000 d. Multiperiodic modulation was detected in 47 stars. The appearance of modulation in the frequency domain was studied in detail. Modulation manifests either as close doublets or as equidistant triplets and multiplets centered on radial mode frequency and its harmonics. In a significant fraction (29 percent) of stars, we have detected the modulation frequency itself, which corresponds to the modulation of the mean stellar brightness. Our search for period doubling effect, that was discovered recently in modulated fundamental mode RR Lyrae stars, and triggered development of new model behind the Blazhko modulation, yielded negative result. In 104 stars we detected additional signals that could correspond to both radial and non-radial modes. Statistical properties of modulated stars were analyzed in detail and confronted with properties of non-modulated stars and of modulated fundamental mode RR Lyrae stars. Our analysis provides constraints for the models to explain the Blazhko phenomenon, which still remains a puzzle more than hundred years after its discovery.

We present the analysis of a peculiar W~Virginis (pWVir) type-II Cepheid, OGLE-LMC-T2CEP-211 ($P_{puls}=9.393\,d$), in a double-lined binary system ($P_{orb}=242\,d$), which shed light on virtually unknown evolutionary status and structure of pWVir stars. The dynamical mass of the Cepheid (first ever for a type-II Cepheid) is $0.64\pm{}0.02\,M_\odot$ and the radius $R=25.1\pm{}0.3\,R_\odot$. The companion is a massive ($5.67\,M_\odot$) main-sequence star obscured by a disk. Such configuration suggests a mass transfer in the system history. We found that originally the system ($P_{orb}^{init}=12\,d$) was composed of $3.5$ and $2.8\,M_\odot$ stars, with the current Cepheid being more massive. The system age is now $\sim{}$200 My, and the Cepheid is almost completely stripped of hydrogen, with helium mass of $\sim{}92\%$ of the total mass. It finished transferring the mass 2.5 My ago and is evolving towards lower temperatures passing through the instability strip. Comparison with observations indicate a reasonable $2.7\cdot{}10^{-8}\,M_\odot/y$ mass loss from the Cepheid. The companion is most probably a Be main-sequence star with $T=22000\,K$ and $R=2.5\,R_\odot$. Our results yield a good agreement with a pulsation theory model for a hydrogen-deficient pulsator, confirming the described evolutionary scenario. We detected a two-ring disk ($R_{disk}\sim\,116\,R_{\odot}$) and a shell ($R_{shell}\sim\,9\,R_{\odot}$) around the companion, that is probably a combination of the matter from the past mass transfer, the mass being lost by the Cepheid due to wind and pulsations, and a decretion disk around a rapidly rotating secondary. Our study together with observational properties of pWVir stars suggests that their majority are products of a similar binary evolution interaction.

We present a study of a recently discovered stellar overdensity near the northern edge of the Small Magellanic Cloud (SMCNOD). We exploited variable stars from the fourth release of the Optical Gravitational Lensing Experiment (OGLE). We used mainly pulsating variable stars and investigated their potential association with the SMCNOD using their spatial distribution and distances. We found four rather spatially concentrated anomalous Cepheids and eight evenly dispersed RR~Lyrae stars to be most likely members of this overdensity. The anomalous Cepheids inside the SMCNOD trace possible intermediate-age population with ages ranging between $2-4.5$ Gyr. The age distribution of anomalous Cepheids seems to be in a good agreement with the age distribution of anomalous Cepheids in the SMC. Using empirical relations for RR~Lyrae stars we determined the median metallicity for a possible old population in the SMCNOD to be $\rm [Fe/H]_{\rm SMCNOD}=-1.71\pm0.21$ dex, which is in agreement with median metallicity of the old SMC population. The density profile for anomalous Cepheids shows a small anomaly at the position of the SMCNOD, on the other hand, RR~Lyrae variables show no such deviation. The probability of finding the observed number of variable stars at the location of the SMCNOD by chance is very low for anomalous Cepheids (0.7 %) but high for RR~Lyrae stars (13.0 %). Based on its variable stars content, we thus confirm the presence of a modest overdensity in intermediate-age stars in the SMCNOD and conclude that it probably has its origin in the SMC rather than to be the remnant of an accreted dwarf galaxy.

We have analysed the photometry of 924 type II Cepheids from the Galactic bulge observed by the Optical Gravitational Lensing Experiment. Our goal was to detect and analyse various dynamical phenomena, including, e.g. multi-mode pulsation, period-doubled pulsation, period-k pulsation, or quasi-periodic modulation of pulsation. First examples of double-mode BL Her-type stars, pulsating simultaneously in the fundamental and first overtone modes are reported. Two double-mode pulsators are not isolated cases however, but form an extension of recently identified class of long-period RRd stars. Period doubling effect was detected in three BL Her-type stars: two were identified previously and one is a new discovery. Through identification of numerous period-doubled W Vir stars we show that the appearance of the period doubling effect, and hence the transition towards RV Tau class, is a smooth process that starts at pulsation periods slightly above 15d. Interchange of the deep and shallow pulsation cycles is identified in tens of RV Tau variables. One RV Tau star is a strong candidate for period-4 pulsation. Quasi-periodic modulation of pulsation in all sub-classes of type II Cepheids is reported for the first time. Modulation was detected in 16 BL Her, 9 W Vir and in 7 RV Tau stars. In the last case, modulation is the most likely scenario, but two other scenarios, excitation of non-radial modes, or period-4 pulsation, can also account for the observed pulsation pattern. Irregular changes of pulsation amplitude and period on various time-scales are common in all subgroups of type II Cepheids.

We briefly summarize the BRITE observations of classical Cepheids so far.

Based on new observations and improved modeling techniques, we have reanalyzed seven Cepheids in the Large Magellanic Cloud. Improved physical parameters have been determined for the exotic system OGLE LMC-CEP-1718 composed of two first-overtone Cepheids and a completely new model was obtained for the OGLE LMC-CEP-1812 classical Cepheid. This is now the shortest period Cepheid for which the projection factor is measured. The typical accuracy of our dynamical masses and radii determinations is 1%. The radii of the six classical Cepheids follow period--radius relations in the literature. Our very accurate physical parameter measurements allow us to calculate a purely empirical, tight period--mass--radius relation that agrees well with theoretical relations derived from non-canonical models. This empirical relation is a powerful tool to calculate accurate masses for single Cepheids for which precise radii can be obtained from Baade--Wesselink-type analyses. The mass of the type-II Cepheid $\kappa$ Pav, $0.56 \pm 0.08 M_\odot$, determined using this relation is in a very good agreement with theoretical predictions. We find large differences between the p-factor values derived for the Cepheids in our sample. Evidence is presented that a simple period--p-factor relation shows an intrinsic dispersion, hinting at the relevance of other parameters, such as the masses, radii, and radial velocity variation amplitudes. We also find evidence that the systematic blueshift exhibited by Cepheids, is primarily correlated with their gravity. The companion star of the Cepheid in the OGLE LMC-CEP-4506 system has a very similar temperature and luminosity, and is clearly located inside the Cepheid instability strip, yet it is not pulsating.

We present a determination of precise fundamental physical parameters of twenty detached, double- lined, eclipsing binary stars in the Large Magellanic Cloud (LMC) containing G- or early K-type giant stars. Eleven are new systems, the remaining nine are systems already analyzed by our team for which we present updated parameters. The catalogue results from our long-term survey of eclipsing binaries in the Magellanic Clouds suitable for high-precision determination of distances (the Araucaria project). The V-band brightnesses of the systems range from 15.4 mag to 17.7 mag and their orbital periods range from 49 days to 773 days. Six systems have favorable geometry showing total eclipses. The absolute dimensions of all eclipsing binary components are calculated with a precision of better than 3% and all systems are suitable for a precise distance determination. The measured stellar masses are in the range 1.4 to 4.6 M_sun and comparison with the MESA isochrones gives ages between 0.1 and 2.1 Gyr. The systems show some weak age-metallicity relation. Two systems have components with very different masses: OGLE LMC-ECL-05430 and OGLE LMC-ECL-18365. Neither system can be fitted by single stellar evolution isochrone, explained by a past mass transfer scenario in the case of ECL-18365 and a gravitational capture or a hierarchical binary merger scenario in the case of ECL-05430. The longest period system OGLE LMC SC9 230659 shows a surprising apsidal motion which shifts the apparent position of the eclipses. In one spectrum of OGLE LMC-ECL-12669 we noted a peculiar dimming of one of the components by 65% well outside of the eclipses. We interpret this observation as arising from an extremely rare occultation event as a foreground Galactic object covers only one component of an extragalactic eclipsing binary.

Physical parameters were recently measured for several Cepheids in eclipsing binary systems in the LMC. It is a good opportunity to compare these results with models obtained from pulsation theory. Having well determined physical parameters of the stars, we can calculate the corresponding pulsation periods for a given mode and check for instability. For stars that lack data, unknown or poorly defined parameters can be further constrained. We present the comparison for six Cepheids, showing that, in general, the results are in very good agreement.

We have analysed photometric data for the first overtone Cepheids in the Magellanic Clouds from the OGLE collection. In more than 500 stars, we have detected additional variability with periods shorter than the first overtone period and period ratios in the (0.60, 0.65) range. The sample includes double periodic stars detected previously by the OGLE team as well as new discoveries. In the Petersen diagram, these stars form three well-separated sequences. In some stars, we have found simultaneously two close periodicities corresponding to two sequences in the Petersen diagram. In a significant fraction of stars, we have detected the power excess at half the frequency of the additional variability. Interestingly, most of these stars form the middle sequence in the Petersen diagram.

Time-series photometry from space in the ultraviolet can be presently done with only a few platforms, none of which is able to provide wide-field long-term high-cadence photometry. We present a concept of UVSat, a twin space telescope which will be capable to perform this kind of photometry, filling an observational niche. The satellite will host two telescopes, one for observations in the ultraviolet, the other for observations in the optical band. We also briefly show what science can be done with UVSat.

We present the results of a study of the type II Cepheid ($P_{puls} = 4.974 d$) in the eclipsing binary system OGLE-LMC-T2CEP-098 ($P_{orb} = 397.2 d$). The Cepheid belongs to the peculiar W Vir group, for which the evolutionary status is virtually unknown. It is the first single-lined system with a pulsating component analyzed using the method developed by Pilecki et al. (2013). We show that the presence of a pulsator makes it possible to derive accurate physical parameters of the stars even if radial velocities can be measured for only one of the components. We have used four different methods to limit and estimate the physical parameters, eventually obtaining precise results by combining pulsation theory with the spectroscopic and photometric solutions. The Cepheid radius, mass and temperature are $25.3 \pm 0.2 R_\odot$, $1.51 \pm 0.09 M_\odot$ and $5300 \pm 100 K$, respectively, while its companion has similar size ($26.3 R_\odot$), but is more massive ($6.8 M_\odot$) and hotter ($9500 K$). Our best estimate for the p-factor of the Cepheid is $1.30 \pm 0.03$. The mass, position on the period-luminosity diagram, and pulsation amplitude indicate that the pulsating component is very similar to the Anomalous Cepheids, although it has a much longer period and is redder in color. The very unusual combination of the components suggest that the system has passed through a mass transfer phase in its evolution. More complicated internal structure would then explain its peculiarity.

We report the discovery of periodic modulation of pulsation in 51 fundamental mode classical Cepheids of the Magellanic Clouds observed by the Optical Gravitational Lensing Experiment. Although the overall incidence rate is very low, about 1 per cent in each of the Magellanic Clouds, in the case of the SMC and pulsation periods between 12 and 16d the incidence rate is nearly 40 per cent. On the other hand, in the LMC the highest incidence rate is 5 per cent for pulsation periods between 8 and 14d, and the overall amplitude of the effect is smaller. It indicates that the phenomenon is metallicity dependent. Typical modulation periods are between 70 and 300d. In nearly all stars the mean brightness is modulated, which, in principle, may influence the use of classical Cepheids for distance determination. Fortunately, the modulation of mean brightness does not exceed 0.01 mag in all but one star. Also, the effect averages out in typical observations spanning a long time base. Consequently, the effect of modulation on the determination of the distance moduli is negligible. The relative modulation amplitude of the fundamental mode is also low and, with one exception, it does not exceed 6 per cent. The origin of the modulation is unknown. We draw a hypothesis that the modulation is caused by the 2:1 resonance between the fundamental mode and the second overtone that shapes the famous Hertzsprung bump progression.

Over the recent years, the Petersen diagram for classical pulsators, Cepheids and RR Lyr stars, populated with a few hundreds of new multiperiodic variables. We review our analyses of the OGLE data, which resulted in the significant extension of the known, and in the discovery of a few new and distinct forms of multiperiodic pulsation. The showcase includes not only radial mode pulsators, but also radial-non-radial pulsators and stars with significant modulation observed on top of the beat pulsation. First theoretical models explaining the new forms of stellar variability are briefly discussed.

We present the analysis of the top-quality photometry of RR Lyrae stars in the globular cluster NGC 6362, gathered over 11 observing seasons by the CASE project. 16 stars are fundamental mode pulsators (RRab stars) and 16 are first overtone pulsators (RRc stars). In two stars, previously identified as RRab, V3 and V34, we detect additional periodicity identified as radial first overtone mode. Lower than usual period ratios (0.730 and 0.728), dominant pulsation in the radial fundamental mode and presence of a long-period modulation indicate, that these two variables are not classical RRd stars, but are new members of the recently identified class of anomalous RRd variables. In a significant fraction of RRc stars, 63 per cent, we detect additional shorter-period variability in the (0.60, 0.65)$P_1$ range. This form of double-periodic pulsation must be common in first overtone RR Lyr stars, as space observations indicate. The incidence rate we find in NGC 6362, is the highest in ground-based observations reported so far. We study the properties of these stars in detail; in particular we confirm that in the colour-magnitude diagram, this group is adjacent to the interface between RRab and RRc stars, as first reported in the analysis of M3 observations by Jurcsik et al. The incidence rate of the Blazhko effect is also very high: we observe it in 69 per cent of RRab stars and in 19 per cent of RRc stars. Rare, double-periodic modulation is reported in one RRab and in one RRc star. Finally we discuss V37 - a peculiar variable in which we detect two close high-amplitude periodicities and modulation. Its previous classification as RRc must be treated as tentative.

We report the discovery of a new group of double-periodic stars in the OGLE Galactic bulge photometry. In 38 stars identified as fundamental mode RR~Lyrae and 4 classified as the first-overtone RR~Lyrae, we detected additional shorter periodicity. Periods of the dominant variability in the newly discovered group are $0.28<P_{\rm D}<0.41$\,days. Period ratios $(0.68 - 0.72)$ are smaller than the period ratios of the Galactic bulge RRd stars. The typical amplitude ratio (of the additional to the dominant periodicity) is 20\,\% for the stars identified as fundamental mode RR~Lyrae and 50\,\% for stars classified as the first-overtone RR~Lyrae. Ten stars from our sample exhibit equidistant peaks in the frequency spectrum, that suggest the Blazhko-type modulation of the main pulsation frequency and/or the additional periodicity. The Fourier coefficients $R_{\mathrm{21}}$ and $R_{\mathrm{31}}$ are one of the lowest among fundamental mode RR Lyrae stars, but among the highest for the first-overtone pulsators. For the phase Fourier coefficients $\varphi_{\mathrm{21}}$ and $\varphi_{\mathrm{31}}$, our stars lie in between RRab and RRc stars. Discussed stars were compared with the radial, linear pulsation models. Their position in the Petersen diagram cannot be reproduced assuming that two radial modes are excited and their physical parameters are similar to that characteristic for RR~Lyrae stars. The non-radial mode scenario also faces difficulties. We conclude that the dominant variability is most likely due to pulsation in the radial fundamental mode including stars classified as the first overtone mode pulsators. At this point, we cannot explain the nature of the additional periodicity. Even more, classification of the stars as RR~Lyrae should be treated as tentative.

Single star evolution does not allow extremely low-mass stars to cross the classical instability strip (IS) during the Hubble time. However, within binary evolution framework low-mass stars can appear inside the IS once the mass transfer (MT) is taken into account. Triggered by a discovery of low-mass 0.26 Msun RR Lyrae-like variable in a binary system, OGLE-BLG-RRLYR-02792, we investigate the occurrence of similar binary components in the IS, which set up a new class of low-mass pulsators. They are referred to as Binary Evolution Pulsators (BEPs) to underline the interaction between components, which is crucial for substantial mass loss prior to the IS entrance. We simulate a population of 500 000 metal-rich binaries and report that 28 143 components of binary systems experience severe MT (loosing up to 90% of mass), followed by at least one IS crossing in luminosity range of RR Lyrae (RRL) or Cepheid variables. A half of these systems enter the IS before the age of 4 Gyr. BEPs display a variety of physical and orbital parameters, with the most important being the BEP mass in range 0.2-0.8 Msun, and the orbital period in range 10-2500 d. Based on the light curve only, BEPs can be misclassified as genuine classical pulsators, and as such they would contaminate genuine RRL and classical Cepheid variables at levels of 0.8% and 5%, respectively. We state that the majority of BEPs will remain undetected and we discuss relevant detection limitations.

We briefly summarize the BRITE observations of classical Cepheids. Possible detection of modulation in a fundamental mode Cepheid, T Vul, and of additional non-radial modes in first overtone Cepheids, DT Cyg and V1334 Cyg, are reported.

The BRIght Target Explorer (BRITE) Constellation is the first nanosatellite mission applied to astrophysical research. Five satellites in low-Earth orbits perform precise optical two-colour photometry of the brightest stars in the night sky. BRITE is naturally well suited for variability studies of hot stars. This contribution describes the basic outline of the mission and some initial problems that needed to be overcome. Some information on BRITE data products, how to access them, and how to join their scientific exploration is provided. Finally, a brief summary of the first scientific results obtained by BRITE is given.

We present the first analysis of W Vir stars observed by the Kepler space telescope in the K2 mission. Clear cycle-to-cycle variation were detected in the light curves of KT Sco and the globular cluster member M80-V1. While the variations in the former star seems to be irregular on the short time scale of the K2 data, the latter appears to experience period doubling in its pulsation. Ground-based colour data confirmed that both stars are W Vir-type pulsators, while a comparison with historical photometric time-series data revealed drastic period changes in both stars. For comparison we reexamine ground-based observations of W Vir, the prototype of the class, and conclude that it shows period doubling instead of mode beating. These results support the notion that nonlinear dynamics plays an important role in the pulsation of W Virginis-type stars.

In the course of a project to study eclipsing binary stars in vinicity of the Sun, we found that the cooler component of LL Aqr is a solar twin candidate. This is the first known star with properties of a solar twin existing in a non-interacting eclipsing binary, offering an excellent opportunity to fully characterise its physical properties with very high precision. We used extensive multi-band, archival photometry and the Super-WASP project and high-resolution spectroscopy obtained from the HARPS and CORALIE spectrographs. The spectra of both components were decomposed and a detailed LTE abundance analysis was performed. The light and radial velocity curves were simultanously analysed with the Wilson-Devinney code. The resulting highly precise stellar parameters were used for a detailed comparison with PARSEC, MESA, and GARSTEC stellar evolution models. LL Aqr consists of two main-sequence stars (F9 V + G3 V) with masses of M1 = 1.1949$\pm$0.0007 and M2=1.0337$\pm$0.0007 $M_\odot$, radii R1 = 1.321$\pm$0.006 and R2 = 1.002$\pm$0.005 $R_\odot$, temperatures T1=6080$\pm$45 K and T2=5703$\pm$50 K and solar chemical composition [M/H]=0.02$\pm$0.05 dex. The absolute dimensions, radiative and photometric properties, and atmospheric abundances of the secondary are all fully consistent with being a solar twin. Both stars are cooler by about 3.5 $\sigma$ or less metal abundant by 5$\sigma$ than predicted by standard sets of stellar evolution models. When advanced modelling was performed, we found that full agreement with observations can only be obtained for values of the mixing length and envelope overshooting parameters that are hard to accept. The most reasonable and physically justified model fits found with MESA and GARSTEC codes still have discrepancies with observations but only at the level of 1$\sigma$.

We report the discovery of a new subclass of double-mode RR Lyrae stars in the Large and Small Magellanic Clouds. The sample of 22 pulsating stars have been extracted from the latest edition of the OGLE collection of RR Lyrae variables in the Magellanic System. The stars pulsating simultaneously in the fundamental (F) and first-overtone (1O) modes have distinctly different properties than regular double-mode RR Lyrae variables (RRd stars). The $P_{1O}/P_F$ period ratios of our anomalous RRd stars are within a range 0.725-0.738, while "classical" double-mode RR Lyrae variables have period ratios in the range 0.742-0.748. In contrast to the typical RRd stars, in the majority of the anomalous pulsators the F-mode amplitudes are higher than the 1O-mode amplitudes. The light curves associated with the F-mode in the anomalous RRd stars show different morphology than the light curves of, both, regular RRd stars and single-mode RRab stars. Most of the anomalous double-mode stars show long-term modulations of the amplitudes (Blazhko-like effect). Translating the period ratios into the abundance parameter, Z, we find for our stars Z in the range (0.002,0.005) - an order of magnitude higher values than typical for RR Lyrae stars. The mass range of the RRd stars inferred from the $W_I$ vs. $P_F$ diagram is (0.55-0.75) of the solar mass. These parameters cannot be accounted for with single star evolution assuming a Reimers-like mass loss. Much greater mass loss caused by interaction with other stars is postulated. We blame the peculiar pulsation properties of our stars to the parametric resonance instability of the 1O-mode to excitation of the F- and 2O-modes as with the inferred parameters of the stars $2\omega_{\rm 1O}\approx\omega_{\rm F}+\omega_{\rm 2O}$.

We present the discovery of a new, peculiar form of double-periodic pulsation in RR Lyrae stars. In four, long-period ($P>0.6$ d) stars observed by the Optical Gravitational Lensing Experiment, and classified as fundamental mode pulsators (RRab), we detect additional, low-amplitude variability, with period shorter than fundamental mode period. The period ratios fall in a range similar to double-mode fundamental and first overtone RR Lyrae stars (RRd), with the exception of one star, in which the period ratio is significantly lower and nearly exactly equals 0.7. Although period ratios are fairly different for the four stars, the light curve shapes corresponding to the dominant, fundamental mode are very similar. The peak-to-peak amplitudes and amplitude ratios (Fourier parameters $R_{21}$ and $R_{31}$) are among the highest observed in RRab stars of similar period, while Fourier phases ($\varphi_{21}$ and $\varphi_{31}$) are among the lowest observed in RRab stars. If the additional variability is interpreted as due to radial first overtone, then, the four stars are the most extreme RRd variables of the longest pulsation periods known. Indeed, the observed period ratios can be well modelled with high metallicity pulsation models. However, at such long pulsation periods, first overtone is typically damped. Five other candidates, with weak signature of additional variability, sharing the same characteristics, were also detected and are briefly discussed.

RR Lyrae stars are classical pulsating stars. They pulsate mostly in the radial fundamental mode (RRab stars), in the radial first overtone mode (RRc stars), or in both modes simultaneously (RRd stars). Collection of variable stars from the Optical Gravitational Lensing Experiment (OGLE) contains more than 38 000 RR Lyrae stars from the Galactic bulge. We analysed these data for RRc and RRd stars. We have found new members of radial-non-radial double-mode RR Lyrae stars, with characteristic period ratio of the two modes around 0.61. We increased the number of known RR Lyrae stars of this type by a factor of 8. We have also discovered another group of double-mode RR Lyrae stars. They pulsate in the first overtone and in another, unidentified mode, which has period longer than period of the undetected fundamental mode. The period ratios tightly cluster around 0.686. These proceedings are focused on this puzzling group. In particular, we report eight new members of the group.

We analyse photometry for 138 first overtone Cepheids from the Small Magellanic Cloud, in which Optical Gravitational Lensing Experiment (OGLE) team discovered additional variability with period shorter than first overtone period, and period ratios in the (0.60, 0.65) range. In the Petersen diagram these stars form three well separated sequences. The additional variability cannot correspond to other radial mode. This form of pulsation is still puzzling. We find that amplitude of the additional variability is small, typically 2-4 per cent of the first overtone amplitude, which corresponds to 2-5 mmag. In some stars we find simultaneously two close periodicities corresponding to two sequences in the Petersen diagram. The most important finding is the detection of power excess at half the frequency of the additional variability (at subharmonic) in 35 per cent of the analysed stars. Interestingly, power excess at subharmonic frequency is detected mostly for stars of the middle sequence in the Petersen diagram (74 per cent), incidence rate is much lower for stars of the top sequence (31 per cent), and phenomenon is not detected for stars of the bottom sequence. The amplitude and/or phase of the additional periodicities strongly vary in time. Similar form of pulsation is observed in first overtone RR Lyrae stars. Our results indicate that the nature and cause of this form of pulsation is the same in the two groups of classical pulsators; consequently, a common model explaining this form of pulsation should be searched for. Our results favour the theory of the excitation of non-radial modes of angular degrees 7, 8 and 9, proposed recently by Dziembowski.

We review the recent progress in our understanding of the Blazhko effect in RR Lyrae stars.

The Optical Gravitational Lensing Experiment (OGLE) is a great source of top-quality photometry of classical pulsators. Collection of variable stars from the fourth part of the project contains more than 38 000 RR Lyrae stars. These stars pulsate mostly in the radial fundamental mode (RRab), in radial first overtone (RRc) or in both modes simultaneously (RRd). Analysis of the OGLE data allowed to detect additional non-radial modes in RRc and in RRd stars. We have found more than 260 double-mode stars with characteristic period ratio of the additional (shorter) period to first overtone period around 0.61, increasing the number of known stars of this type by factor of 10. Stars from the OGLE sample form three nearly parallel sequences in the Petersen diagram. Some stars show more than one non-radial mode simultaneously. These modes belong to different sequences.

We briefly discuss three cases of RR~Lyrae variable stars with Blazhko modulation from the OGLE Galactic bulge collection. In one of the discussed stars modulation with two periods is detected, in other star, we detect three Blazhko modulations.

With the help of linear pulsation models we briefly discuss the peculiar period ratios characteristic for modulated double-mode RR Lyrae stars.

We present a grid on non-linear convective type-II Cepheid models. The dense model grids are computed for 0.6M_Sun and a range of metallicities ([Fe/H]=-2.0,-1.5,-1.0), and for 0.8M_Sun ([Fe/H]=-1.5). Two sets of convective parameters are considered. The models cover the full temperature extent of the classical instability strip, but are limited in luminosity; for the most luminous models violent pulsation leads to the decoupling of the outermost model shell. Hence, our survey reaches only the shortest period RV Tau domain. In the Hertzsprung-Russel diagram we detect two domains in which period doubled pulsation is possible. The first extends through the BL Her domain and low luminosity W Vir domain (pulsation periods ~2-6.5 d). The second domain extends at higher luminosities (W Vir domain; periods >9.5d). Some models within these domains display period-4 pulsation. We also detect very narrow domains (~10 K wide) in which modulation of pulsation is possible. Another interesting phenomenon we detect is double-mode pulsation in the fundamental mode and in the fourth radial overtone. Fourth overtone is a surface mode, trapped in the outer model layers. Single-mode pulsation in the fourth overtone is also possible on the hot side of the classical instability strip. The origin of the above phenomena is discussed. In particular, the role of resonances in driving different pulsation dynamics as well as in shaping the morphology of the radius variation curves is analysed.

We present a detailed study of the classical Cepheid in the double-lined, highly eccentric eclipsing binary system OGLE-LMC562.05.9009. The Cepheid is a fundamental mode pulsator with a period of 2.988 days. The orbital period of the system is 1550 days. Using spectroscopic data from three 4-8-m telescopes and photometry spanning 22 years, we were able to derive the dynamical masses and radii of both stars with exquisite accuracy. Both stars in the system are very similar in mass, radius and color, but the companion is a stable, non-pulsating star. The Cepheid is slightly more massive and bigger (M_1 = 3.70 +/- 0.03M_sun, R_1 = 28.6 +/- 0.2R_sun) than its companion (M_2 = 3.60 +/- 0.03M_sun, R_2 = 26.6 +/- 0.2R_sun). Within the observational uncertainties both stars have the same effective temperature of 6030 +/- 150K. Evolutionary tracks place both stars inside the classical Cepheid instability strip, but it is likely that future improved temperature estimates will move the stable giant companion just beyond the red edge of the instability strip. Within current observational and theoretical uncertainties, both stars fit on a 205 Myr isochrone arguing for their common age. From our model, we determine a value of the projection factor of p = 1.37 +/- 0.07 for the Cepheid in the OGLE-LMC562.05.9009 system. This is the second Cepheid for which we could measure its p-factor with high precision directly from the analysis of an eclipsing binary system, which represents an important contribution towards a better calibration of Baade-Wesselink methods of distance determination for Cepheids.

We analyse the OGLE-IV photometry of the first overtone and double-mode RR Lyrae stars (RRc/RRd) in the two fields towards the Galactic bulge observed with high cadence. In 27 per cent of RRc stars we find additional non-radial mode, with characteristic period ratio, P x /P 1O \in (0.6, 0.64). It strongly corroborates the conclusion arising from the analysis of space photometry of RRc stars, that this form of pulsation must be common. In the Petersen diagram the stars form three sequences. In 20 stars we find two or three close secondary modes simultaneously. The additional modes are clearly non-stationary. Their amplitude and/or phase vary in time. As a result, the patterns observed in the frequency spectra of these stars may be very complex. In some stars the additional modes split into doublets, triplets or appear as a more complex bands of increased power. Subharmonics of additional modes are detected in 20 per cent of stars. They also display a complex structure. Including our previous study of the OGLE-III Galactic bulge data, we have discovered 260 RRc and 2 RRd stars with the additional non-radial mode, which is the largest sample of these stars so far. The additional mode is also detected in two Blazhko RRc stars, which shows that the modulation and additional non-radial mode are not exclusive.

A detailed analysis is presented of 33 RR Lyrae stars in Pisces observed with the Kepler space telescope over the 8.9-day long K2 Two-Wheel Concept Engineering Test. The sample includes not only fundamental-mode and first overtone (RRab and RRc) stars but the first two double-mode (RRd) stars that Kepler detected and the only modulated first-overtone star ever observed from space so far. The precision of the extracted K2 light curves made it possible to detect low-amplitude additional modes in all subtypes. All RRd and non-modulated RRc stars show the additional mode at P_X/P_1~0.61 that was detected in previous space-based photometric measurements. A periodicity longer than the fundamental mode was tentatively identified in one RRab star that might belong to a gravity mode. We determined the photometric [Fe/H] values for all fundamental-mode stars and provide the preliminary results of our efforts to fit the double-mode stars with non-linear hydrodynamic pulsation models. The results from this short test run indicate that the K2 mission will be, and has started to be, an ideal tool to expand our knowledge about RR Lyrae stars. As a by-product of the target search and analysis, we identified 165 bona-fide double-mode RR Lyrae stars from the Catalina Sky Survey observations throughout the sky, 130 of which are new discoveries.

We have analyzed the double-lined eclipsing binary system ASAS J180057-2333.8 from the All Sky Automated Survey (ASAS) catalogue . We measure absolute physical and orbital parameters for this system based on archival $V$-band and $I$-band ASAS photometry, as well as on high-resolution spectroscopic data obtained with ESO 3.6m/HARPS and CORALIE spectrographs. The physical and orbital parameters of the system were derived with an accuracy of about 0.5 - 3%. The system is a very rare configuration of two bright well-detached giants of spectral types K1 and K4 and luminosity class II. The radii of the stars are $R_1$ = 52.12 $\pm$ 1.38 and $R_2$ = 67.63 $\pm$ 1.40 R$_\odot$ and their masses are $M_1$ = 4.914 $\pm$ 0.021 and $M_2$ = 4.875$\pm$ 0.021 M$_\odot$ . The exquisite accuracy of 0.5% obtained for the masses of the components is one of the best mass determinations for giants. We derived a precise distance to the system of 2.14 $\pm$ 0.06 kpc (stat.) $\pm$ 0.05 (syst.) which places the star in the Sagittarius-Carina arm. The Galactic rotational velocity of the star is $\Theta_s=258 \pm 26$ km s$^{-1}$ assuming $\Theta_0=238$ km s$^{-1}$. A comparison with PARSEC isochrones places the system at the early phase of core helium burning with an age of slightly larger than 100 million years. The effect of overshooting on stellar evolutionary tracks was explored using the MESA star code.

We report the discovery of a new group of double-periodic RR Lyrae stars from the analysis of the OGLE-IV Galactic bulge photometry. In 11 stars identified in the OGLE catalog as first overtone pulsators (RRc stars) we detect additional longer period variability of low amplitude, in the mmag regime. One additional star of the same type is identified in a published analysis of the Kepler space photometry. The period ratio between the shorter first overtone period and a new, longer period lies in a narrow range around 0.686. Thus, the additional period is longer than the expected period of the undetected radial fundamental mode. The obvious conclusion that addition periodicity corresponds to a gravity or a mixed mode faces difficulties, however.

We present here the first spectroscopic and photometric analysis of the double-lined eclipsing binary containing the classical, first-overtone Cepheid OGLE-LMC-CEP-2532 (MACHO 81.8997.87). The system has an orbital period of 800 days and the Cepheid is pulsating with a period of 2.035 days. Using spectroscopic data from three high-class telescopes and photometry from three surveys spanning 7500 days we are able to derive the dynamical masses for both stars with an accuracy better than 3%. This makes the Cepheid in this system one of a few classical Cepheids with an accurate dynamical mass determination (M_1=3.90 +/- 0.10 M_sun). The companion is probably slightly less massive (3.82 +/- 0.10 M_sun), but may have the same mass within errors (M_2/M_1= 0.981 +/- 0.015). The system has an age of about 185 million years and the Cepheid is in a more advanced evolutionary stage. For the first time precise parameters are derived for both stars in this system. Due to the lack of the secondary eclipse for many years not much was known about the Cepheid's companion. In our analysis we used extra information from the pulsations and the orbital solution from the radial velocity curve. The best model predicts a grazing secondary eclipse shallower than 1 mmag, hence undetectable in the data, about 370 days after the primary eclipse. The dynamical mass obtained here is the most accurate known for a first-overtone Cepheid and will contribute to the solution of the Cepheid mass discrepancy problem.

We report the discovery of an intriguing triple-mode RR Lyrae star found in the OGLE Galactic bulge collection, OGLE-BLG-RRLYR-24137. In the OGLE catalog the star was identified as RRd star - double-mode pulsator, pulsating simultaneously in the fundamental and in the first overtone modes. We find that third mode is excited and firmly detect its period doubling. Period ratios are not far from that expected for triple-mode - fundamental, first and third overtone - pulsation. Unfortunately, we cannot reproduce period ratios of the three modes with a consistent set of pulsation models. Therefore the other interpretation, that additional mode is non-radial, is also likely.

We present the analysis of Blazhko-type modulation in double-mode RR Lyrae (RRd) stars found in the OGLE photometry of the Galactic bulge. Modulation is detected in 15 stars. Most of them have non-typical period ratio of the radial modes. In the Petersen diagram, at a given period of the fundamental mode, they are located significantly below or above the sequence formed by the majority of RRd stars. Multi-periodic modulation is very frequent; two or three modulation periods are detected in 8 stars. Modulation periods vary from ~20 to more than 300 days. Radial mode amplitudes can be modulated by a few to nearly hundred per cent. Both radial modes may be modulated with the same period. More commonly however, dominant modulation for the fundamental mode has different period than dominant modulation for the first overtone. Quite often modulation of only one mode is detected in the data. We find a clear feedback between pulsation amplitude of the dominant mode and mean stellar brightness: lower the pulsation amplitude, brighter the star. At phases of low pulsation amplitude, the mode periods are prone to fast changes. All the stars share the common feature: their pulsation properties are non-stationary. Amplitudes and phases of the radial modes vary irregularly on a long time-scale of a few hundred or thousand days. The short-term modulations are also irregular. One of the stars has switched the pulsation mode recently: from single-mode fundamental mode pulsation to RRd state. In other star the non-radial mode with characteristic ~0.61 period ratio to the first overtone is detected. This non-radial mode is likely modulated with the same period as the radial modes.

We present the analysis of four first overtone RR Lyrae stars observed with the Kepler space telescope, based on data obtained over nearly 2.5yr. All four stars are found to be multiperiodic. The strongest secondary mode with frequency f_2 has an amplitude of a few mmag, 20 - 45 times lower than the main radial mode with frequency f_1. The two oscillations have a period ratio of P_2/P_1 = 0.612 - 0.632 that cannot be reproduced by any two radial modes. Thus, the secondary mode is nonradial. Modes yielding similar period ratios have also recently been discovered in other variables of the RRc and RRd types. These objects form a homogenous group and constitute a new class of multimode RR Lyrae pulsators, analogous to a similar class of multimode classical Cepheids in the Magellanic Clouds. Because a secondary mode with P_2/P_1 ~ 0.61 is found in almost every RRc and RRd star observed from space, this form of multiperiodicity must be common. In all four Kepler RRc stars studied, we find subharmonics of f_2 at ~1/2 f_2 and at ~3/2 f_2. This is a signature of period doubling of the secondary oscillation, and is the first detection of period doubling in RRc stars. The amplitudes and phases of f_2 and its subharmonics are variable on a timescale of 10 - 200d. The dominant radial mode also shows variations on the same timescale, but with much smaller amplitude. In three Kepler RRc stars we detect additional periodicities, with amplitudes below 1mmag, that must correspond to nonradial g-modes. Such modes never before have been observed in RR Lyrae variables.

Non-radial modes are excited in classical pulsators, both in Cepheids and in RR Lyrae stars. Firm evidence come from the first overtone pulsators, in which additional shorter period mode is detected with characteristic period ratio falling in between 0.60 and 0.65. In the case of first overtone Cepheids three separate sequences populated by nearly 200 stars are formed in the Petersen diagram, i.e. the diagram of period ratio versus longer period. In the case of first overtone RR Lyrae stars (RRc stars) situation is less clear. A dozen or so such stars are known which form a clump in the Petersen diagram without any obvious structure. Interestingly, all first overtone RR Lyrae stars for which precise space-borne photometry is available show the additional mode, which suggests that its excitation is common. Motivated by these results we searched for non-radial modes in the OGLE-III photometry of RRc stars from the Galactic bulge. We report the discovery of 147 stars, members of a new group of double-mode, radial-non-radial mode pulsators. They form a clear and tight sequence in the Petersen diagram, with period ratios clustering around 0.613 with a signature of possible second sequence with higher period ratio (0.631). The scatter in period ratios of the already known stars is explained as due to population effects. Judging from the results of space observations this still mysterious form of pulsation must be common among RRc stars and with our analysis of the OGLE data we just touch the tip of the iceberg.

We present non-linear, convective, BL Her-type hydrodynamic models that show complex variability characteristic for deterministic chaos. The bifurcation diagram reveals a rich structure, with many phenomena detected for the first time in hydrodynamic models of pulsating stars. The phenomena include not only period doubling cascades en route to chaos (detected in earlier studies) but also periodic windows within chaotic band, type-I and type-III intermittent behaviour, interior crisis bifurcation and others. Such phenomena are known in many textbook chaotic systems, from the simplest discrete logistic map, to more complex systems like Lorenz equations. We discuss the physical relevance of our models. Although except of period doubling such phenomena were not detected in any BL Her star, chaotic variability was claimed in several higher luminosity siblings of BL Her stars - RV Tau variables, and also in longer-period, luminous irregular pulsators. Our models may help to understand these poorly studied stars. Particularly interesting are periodic windows which are intrinsic property of chaotic systems and are not necessarily caused by resonances between pulsation modes, as sometimes claimed in the literature.

We have obtained extensive high-quality spectroscopic observations of the OGLE-LMC-CEP-1718 eclipsing binary system in the Large Magellanic Cloud which Soszynski et al. (2008) had identified as a candidate system for containing two classical Cepheids in orbit. Our spectroscopic data clearly demonstrate binary motion of the Cepheids in a 413-day eccentric orbit, rendering this eclipsing binary system the first ever known to consist of two classical Cepheid variables. After disentangling the four different radial velocity variations in the system we present the orbital solution and the individual pulsational radial velocity curves of the Cepheids. We show that both Cepheids are extremely likely to be first overtone pulsators and determine their respective dynamical masses, which turn out to be equal to within 1.5 %. Since the secondary eclipse is not observed in the orbital light curve we cannot derive the individual radii of the Cepheids, but the sum of their radii derived from the photometry is consistent with overtone pulsation for both variables. The existence of two equal-mass Cepheids in a binary system having different pulsation periods (1.96 and 2.48 days, respectively) may pose an interesting challenge to stellar evolution and pulsation theories, and a more detailed study of this system using additional datasets should yield deeper insight about the physics of stellar evolution of Cepheid variables. Future analysis of the system using additional near-infrared photometry might also lead to a better understanding of the systematic uncertainties in current Baade-Wesselink techniques of distance determinations to Cepheid variables.

Many dynamical systems of different complexity, e.g. 1D logistic map, the Lorentz equations, or real phenomena, like turbulent convection, show chaotic behaviour. Despite huge differences, the dynamical scenarios for these systems are strikingly similar: chaotic bands are born through the series of period doubling bifurcations and merge through interior crises. Within chaotic bands periodic windows are born through the tangent bifurcations, preceded by the intermittent behaviour. This is a universal behaviour of dynamical systems (Feigenbaum 1983). We demonstrate such behaviour in models of pulsating stars.

In this review we focus on non-linear phenomena in pulsating stars the mode selection and amplitude limitation. Of many linearly excited modes only a fraction is detected in pulsating stars. Which of them and why (the problem of mode selection) and to what amplitude (the problem of amplitude limitation) are intrinsically non-linear and still unsolved problems. Tools for studying these problems are briefly discussed and our understanding of mode selection and amplitude limitation in selected groups of self-excited pulsators is presented. Focus is put on classical pulsators (Cepheids and RR Lyrae stars) and main sequence variables (delta Scuti and beta Cephei stars). Directions of future studies are briefly discussed.

As a response to the white paper call, we propose to turn Kepler to the South Ecliptic Pole (SEP) and observe thousands of large amplitude variables for years with high cadence in the frame of the Kepler-SEP Mission. The degraded pointing stability will still allow observing these stars with reasonable (probably better than mmag) accuracy. Long-term continuous monitoring already proved to be extremely helpful to investigate several areas of stellar astrophysics. Space-based missions opened a new window to the dynamics of pulsation in several class of pulsating variable stars and facilitated detailed studies of eclipsing binaries. The main aim of this mission is to better understand the fascinating dynamics behind various stellar pulsational phenomena (resonances, mode coupling, chaos, mode selection) and interior physics (turbulent convection, opacities). This will also improve the applicability of these astrophysical tools for distance measurements, population and stellar evolution studies. We investigated the pragmatic details of such a mission and found a number of advantages: minimal reprogramming of the flight software, a favorable field of view, access to both galactic and LMC objects. However, the main advantage of the SEP field comes from the large sample of well classified targets, mainly through OGLE. Synergies and significant overlap (spatial, temporal and in brightness) with both ground- (OGLE, LSST) and space-based missions (GAIA, TESS) will greatly enhance the scientific value of the Kepler-SEP mission. GAIA will allow full characterization of the distance indicators. TESS will continuously monitor this field for at least one year, and together with the proposed mission provide long time series that cannot be obtained by other means. If Kepler-SEP program is successful, there is a possibility to place one of the so-called LSST "deep-drilling" fields in this region.

A novel method of analysis of double-lined eclipsing binaries containing a radially pulsating star is presented. The combined pulsating-eclipsing light curve is built up from a purely eclipsing light curve grid created using an existing modeling tool. For every pulsation phase the instantaneous radius and surface brightness are taken into account, being calculated from the disentangled radial velocity curve of the pulsating star and from its out-of-eclipse pulsational light curve and the light ratio of the components, respectively. The best model is found using the Markov Chain Monte Carlo method. The method is applied to the eclipsing binary Cepheid OGLE-LMC-CEP-0227 (P_puls = 3.80 d, P_orb = 309 d). We analyze a set of new spectroscopic and photometric observations for this binary, simultaneously fitting OGLE V-band, I-band and Spitzer 3.6 {\mu}m photometry. We derive a set of fundamental parameters of the system significantly improving the precision comparing to the previous results obtained by our group. The Cepheid mass and radius are M_1 = 4.165 +/- 0.032 M_solar and R_1 = 34.92 +/- 0.34 R_solar, respectively. For the first time a direct, geometrical and distance-independent determination of the Cepheid projection factor is presented. The value p = 1.21 +/- 0.03(stat.) +/- 0.04(syst.) is consistent with theoretical expectations for a short period Cepheid and interferometric measurements for {\delta} Cep. We also find a very high value of the optical limb darkening coefficients for the Cepheid component, in strong disagreement with theoretical predictions for static atmospheres at a given surface temperature and gravity.

We performed a new and accurate fit of light and radial velocity curves of the Large Magellanic Cloud (LMC) Cepheid --OGLE-LMC-CEP-0227-- belonging to a detached double-lined eclipsing binary system. We computed several sets of nonlinear, convective models covering a broad range in stellar mass, effective temperature and in chemical composition. The comparison between theory and observations indicates that current theoretical framework accounts for luminosity --V and I band-- and radial velocity variations over the entire pulsation cycle. Predicted pulsation mass --M=4.14+-0.06 Mo-- and mean effective temperature --Te=6100+-50 K-- do agree with observed estimates with an accuracy better than 1 sigma. The same outcome applies, on average, to the luminosity amplitudes and to the mean radius. We find that the best fit solution requires a chemical composition that is more metal--poor than typical LMC Cepheids (Z=0.004 vs 0.008) and slightly helium enhanced (Y=0.27 vs 0.25), but the sensitivity to He abundance is quite limited. Finally, the best fit model reddening --E(V-I)=0.171+-0.015 mag-- and the true distance modulus corrected for the barycenter of the LMC --mu_{0,LMC}=18.50+-0.02+-0.10 (syst) mag--, agree quite well with similar estimates in the recent literature.

In the era of precision cosmology it is essential to determine the Hubble Constant with an accuracy of 3% or better. Currently, its uncertainty is dominated by the uncertainty in the distance to the Large Magellanic Cloud (LMC) which as the second nearest galaxy serves as the best anchor point of the cosmic distance scale. Observations of eclipsing binaries offer a unique opportunity to precisely and accurately measure stellar parameters and distances. The eclipsing binary method was previously applied to the LMC but the accuracy of the distance results was hampered by the need to model the bright, early-type systems used in these studies. Here, we present distance determinations to eight long-period, late- type eclipsing systems in the LMC composed of cool giant stars. For such systems we can accurately measure both the linear and angular sizes of their components and avoid the most important problems related to the hot early-type systems. Our LMC distance derived from these systems is demonstrably accurate to 2.2 % (49.97 +/- 0.19 (statistical) +/- 1.11 (systematic) kpc) providing a firm base for a 3 % determination of the Hubble Constant, with prospects for improvement to 2 % in the future.

We present non-linear hydrodynamic pulsation models for OGLE-BLG-RRLYR-02792 - a 0.26M_sun pulsator, component of the eclipsing binary system, analysed recently by Pietrzynski et al. The star's light and radial velocity curves mimic that of classical RR Lyrae stars, except for the bump in the middle of the ascending branch of the radial velocity curve. We show that the bump is caused by the 2:1 resonance between the fundamental mode and the second overtone - the same mechanism that causes the Hertzsprung bump progression in classical Cepheids. The models allow to constrain the parameters of the star, in particular to estimate its absolute luminosity (approx 33L_sun) and effective temperature (approx 6970K, close to the blue edge of the instability strip). We conduct a model survey for the new class of low mass pulsators similar to OGLE-BLG-RRLYR-02792 - products of evolution in the binary systems. We compute a grid of models with masses corresponding to half (and less) of the typical mass of RR Lyrae variable, 0.20M_sun<=M<=0.30M_sun, and discuss the properties of the resulting light and radial velocity curves. Resonant bump progression is clear and may be used to distinguish such stars from classical RR Lyrae stars. We present the Fourier decomposition parameters for the modelled light and radial velocity curves. The expected values of the phi_31 Fourier phase for the light curves differ significantly from that observed in RR Lyrae stars, which is another discriminant of the new class.

We present a detailed period analysis of the bright Cepheid-type variable star V1154 Cygni (V =9.1 mag, P~4.9 d) based on almost 600 days of continuous observations by the Kepler space telescope. The data reveal significant cycle-to-cycle fluctuations in the pulsation period, indicating that classical Cepheids may not be as accurate astrophysical clocks as commonly believed: regardless of the specific points used to determine the O-C values, the cycle lengths show a scatter of 0.015-0.02 days over the 120 cycles covered by the observations. A very slight correlation between the individual Fourier parameters and the O-C values was found, suggesting that the O - C variations might be due to the instability of the light curve shape. Random fluctuation tests revealed a linear trend up to a cycle difference 15, but for long term, the period remains around the mean value. We compare the measurements with simulated light curves that were constructed to mimic V1154 Cyg as a perfect pulsator modulated only by the light travel time effect caused by low-mass companions. We show that the observed period jitter in V1154 Cyg represents a serious limitation in the search for binary companions. While the Kepler data are accurate enough to allow the detection of planetary bodies in close orbits around a Cepheid, the astrophysical noise can easily hide the signal of the light-time effect.

We analyzed the Long Cadence photometry of 4 first overtone RR Lyr-type stars (RRc stars) observed by the KEPLER telescope. All studied variables are multiperiodic. The strongest secondary peak appears for f_2/f_1 = 1.58-1.63, or P_2/P_1 = 0.61-0.63. In each star we detect at least one subharmonic of f_2, either at ~1/2 f_2 or at ~3/2 f_2. The presence of subharmonics is a characteristic signature of a period doubling.

We present the hydrodynamic BL Herculis-type models which display a long-term modulation of pulsation amplitudes and phases. The modulation is either strictly periodic or it is quasi-periodic, with the modulation period and modulation pattern varying from one cycle to the other. Such behaviour has not been observed in any BL Her variable so far, however, it is a common property of their lower luminosity siblings - RR Lyrae variables showing the Blazhko effect. These models provide a support for the recent mechanism proposed by Buchler & Kollath to explain this still mysterious phenomenon. In their model, a half-integer resonance that causes the period doubling effect, discovered recently in the Blazhko RR Lyrae stars, is responsible for the modulation of the pulsation as well. Although our models are more luminous than is appropriate for RR Lyrae stars, they clearly demonstrate, through direct hydrodynamic computation, that the mechanism can indeed be operational. Of great importance are models which show quasi-periodic modulation - a phenomenon observed in Blazhko RR Lyrae stars. Our models coupled with the analysis of the amplitude equations show that such behaviour may be caused by the dynamical evolution occurring in the close proximity of the unstable single periodic saddle point.

Rapid and strong changes in the Blazhko modulation of RR Lyrae stars, as they have recently been detected in high precision satellite data, have become a crucial topic in finding an explanation of the long-standing mystery of the Blazhko effect. We present here an analysis of the most extreme case detected so far, the RRab star V445 Lyr (KIC 6186029) which was observed with the Kepler space mission. V445 Lyr shows very strong cycle-to-cycle changes in its Blazhko modulation, which are caused both by a secondary long-term modulation period as well as irregular variations. In addition to the complex Blazhko modulation, V445 Lyr also shows a rich spectrum of additional peaks in the frequency range between the fundamental pulsation and the first harmonic. Among those peaks, the second radial overtone could be identified, which, combined with a metallicity estimate of [Fe/H]=-2.0 dex from spectroscopy, allowed to constrain the mass (0.55-0.65 M_sun) and luminosity (40-50 L_sun) of V445 Lyr through theoretical Petersen diagrams. A non-radial mode as well as possibly the first overtone are also excited. Furthermore, V445 Lyr shows signs of the period doubling phenomenon and a long term period change. A detailed Fourier analysis along with a study of the O-C variation of V445 Lyr is presented, and the origin of the additional peaks and possible causes of the changes in the Blazhko modulation are discussed. The results are then put into context with those of the only other star with a variable Blazhko effect for which a long enough set of high precision continuous satellite data has been published so far, the CoRoT star 105288363.

We report on the discovery of the first BL Herculis star displaying period-doubling behaviour as predicted by the hydrodynamic models of Buchler & Moskalik. The star, with P=2.4d, is located in the Galactic bulge and was discovered with OGLE-III photometry. We present new nonlinear convective models, which, together with recent evolutionary tracks, put constraints on the stellar parameters. In particular, we estimate the mass and metallicity of the object.

We revisit a phenomenological description of turbulent thermal convection along the lines proposed by Gough (1977) in which eddies grow solely by extracting energy from the unstably stratified mean state and are subsequently destroyed by internal shear instability. This work is part of an ongoing investigation for finding a procedure to calculate the turbulent fluxes of heat and momentum in the presence of a shearing background flow in stars.

We report the discovery of a period doubling behaviour in a 2.4d BL Herculis-type variable of the Galactic bulge. Another bulge BL Her-type star (P=2.25d) is a strong period doubling candidate. Both objects have been identified with the OGLE-III photometry. Possibility of period doubling in this type of pulsators has been predicted almost twenty years ago by Buchler & Moskalik. Our finding is the first observational confirmation of their theoretical results. Discovery of the first BL Herculis star showing the period doubling effect motivates a new theoretical investigation with the state-of-the-art convective pulsation codes. We present the results of initial model survey, specifically aimed at studying the observed period-doubled BL Her variable. All of our non-linear models have P=2.4d. The computations confirm that the period doubling effect is caused by the 3:2 resonance between the fundamental mode and the first overtone, as indicated by earlier radiative models of Buchler & Moskalik. Comparison of the computed and the observed light curves allows to constrain the parameters of the star, in particular its metallicity, which appears to be high, Z is approximately 0.01. The recent evolutionary tracks put also constraint on the mass of the star, which is close to M=0.50 solar masses.

This paper summarizes the main results of our recent study of the non-Blazhko RR Lyrae stars observed with the Kepler space telescope. These stars offer the opportunity for studying the stability of the pulsations of RR Lyrae stars and for providing a reference against which the Blazhko RR Lyrae stars can be compared. Of particular interest is the stability of the low-dispersion (sigma < 1mmag) light curves constructed from ~18,000 long-cadence (30-min) and (for FN Lyr and AW Dra) the ~150,000 short-cadence (1-min) photometric data points. Fourier-based [Fe/H] values and other physical characteristics are also derived. When the observed periods are compared with periods computed with the Warsaw non-linear convective pulsation code better agreement is achieved assuming pulsational L and M values rather than the (higher) evolutionary L and M values.

Nineteen of the ~40 RR Lyr stars in the Kepler field have been identified as candidate non-Blazhko (or unmodulated) stars. In this paper we present the results of Fourier decomposition of the time-series photometry of these stars acquired during the first 417 days of operation (Q0-Q5) of the Kepler telescope. Fourier parameters based on ~18400 long-cadence observations per star (and ~150000 short-cadence observations for FN Lyr and for AW Dra) are derived. None of the stars shows the recently discovered `period-doubling' effect seen in Blazhko variables; however, KIC 7021124 has been found to pulsate simultaneously in the fundamental and second overtone modes with a period ratio P2/P0 ~ 0.59305 and is similar to the double-mode star V350 Lyr. Period change rates are derived from O-C diagrams spanning, in some cases, over 100 years; these are compared with high-precision periods derived from the Kepler data alone. Extant Fourier correlations by Kovacs, Jurcsik et al. (with minor transformations from the V to the Kp passband) have been used to derive underlying physical characteristics for all the stars. This procedure seems to be validated through comparisons of the Kepler variables with galactic and LMC RR Lyr stars. The most metal-poor star in the sample is NR Lyr, with [Fe/H]=-2.3 dex; and the four most metal-rich stars have [Fe/H] ranging from -0.6 to +0.1 dex. Pulsational luminosities and masses are found to be systematically smaller than L and mass values derived from stellar evolution models, and are favoured over the evolutionary values when periods are computed with the Warsaw linear hydrodynamics code. Finally, the Fourier parameters are compared with theoretical values derived using the Warsaw non-linear convective pulsation code.

The amplitude and phase modulation observed in a significant fraction of the RR Lyrae variables - the Blazhko effect - represents a long-standing enigma in stellar pulsation theory. No satisfactory explanation for the Blazhko effect has been proposed so far. In this paper we focus on the Stothers (2006) idea, in which modulation is caused by changes in the structure of the outer convective zone, caused by a quasi-periodically changing magnetic field. However, up to this date no quantitative estimates were made to investigate whether such a mechanism can be operational and whether it is capable of reproducing the light variation we observe in Blazhko variables. We address the latter problem. We use a simplified model, in which the variation of turbulent convection is introduced into the non-linear hydrodynamic models in an ad hoc way, neglecting interaction with the magnetic field. We study the light curve variation through the modulation cycle and properties of the resulting frequency spectra. Our results are compared with Kepler observations of RR Lyr. We find that reproducing the light curve variation, as is observed in RR Lyr, requires a huge modulation of the mixing length, of the order of +/-50 per cent, on a relatively short time-scale of less than 40 days. Even then, we are not able to reproduce neither all the observed relations between modulation components present in the frequency spectrum, nor the relations between Fourier parameters describing the shape of the instantaneous light curves.

We report results of initial work done on selected candidate Cepheids to be observed with the Kepler space telescope. Prior to the launch 40 candidates were selected from previous surveys and databases. The analysis of the first 322 days of Kepler photometry, and recent ground-based follow-up multicolour photometry and spectroscopy allowed us to confirm that one of these stars, V1154 Cyg (KIC 7548061), is indeed a 4.9-d Cepheid. Using the phase lag method we show that this star pulsates in the fundamental mode. New radial velocity data are consistent with previous measurements, suggesting that a long-period binary component is unlikely. No evidence is seen in the ultra-precise, nearly uninterrupted Kepler photometry for nonradial or stochastically excited modes at the micromagnitude level. The other candidates are not Cepheids but an interesting mix of possible spotted stars, eclipsing systems and flare stars.

We present our analysis of the long cadence Kepler data for the well-studied Blazhko star RR Lyr, gathered during the first two quarters of the satellite's observations and covering a total of 127d. Besides being of great importance for our understanding of RR Lyrae stars in general, these RR Lyr data can be regarded as a case study for observations of bright stars with Kepler. Kepler can perform high-precision photometry on targets like RR Lyr, as the saturated flux is conserved to a very high degree. The Kepler data on RR Lyr are revolutionary in several respects. Even with long-cadence sampling (one measurement per 29.4 min), the unprecedented precision (< mmag) of the Kepler photometry allows the study of the star's extreme light curve variations in detail. The multiplet structures at the main frequency and its harmonics, typical for Blazhko stars, are clearly detected up to the quintuplets. For the first time, photometric data of RR Lyr reveal the presence of half-integer frequencies, linked to a period doubling effect. This phenomenon may be connected to the still unexplained Blazhko modulation. Moreover, with three observed Blazhko cycles at our disposal, we observe that there is no exact repetition in the light curve changes from one modulation cycle to the next for RR Lyr. This may be due to additional periodicities in the star, or to transient or quasi-periodic changes.

We revisit a phenomenological description of turbulent thermal convection along the lines proposed originally by Gough (1965) in which eddies grow solely by extracting energy from the unstably stratified mean state and are subsequently destroyed by internal shear instability. This work is part of an ongoing investigation for finding a procedure to calculate the turbulent fluxes of heat and momentum in the presence of a shearing background flow in stars.

The phenomenon of double-periodic Cepheid pulsation is still poorly understood. Recently we rediscussed the problem of modelling the double-periodic pulsation with non-linear hydrocodes. We showed that the published non-resonant double-mode models are incorrect, because they exclude the negative buoyancy effects. Aims. We continue our efforts to verify whether the Kuhfuss one-equation convection model with negative buoyancy included can reproduce the double-periodic Cepheid pulsation. Methods. Using the direct time integration hydrocode, which implements the Kuhfuss convection model, we search for stable double-periodic Cepheid models. We search for models pulsating in both fundamental and first overtone modes (F+1O), as well as in the two lowest order overtones (1O+2O). In the latter case, we focus on reproducing double-overtone Cepheids of the Large Magellanic Cloud (LMC). Results. We have found full amplitude non-linear beat Cepheid models of both types, F+1O and 1O+2O. In the case of F+1O models, the beat pulsation is most likely caused by the three-mode resonance, 2omega_1=omega_0+omega_2, while in the double-overtone models the underlying mechanism (resonant or non-resonant) cannot be identified beyond doubt. Double-periodic models found in our survey exist, however, only in narrow period ranges and cannot explain the majority of the observed double-periodic objects. Conclusions. With only little doubt left, we conclude that current one dimensional one-equation convection models are incapable of reproducing the majority of the observed beat Cepheids. Among the shortcomings of current pulsation hydrocodes, the simple treatment of convection seems to be the most severe one. Growing evidence for the presence of non-radial modes in Cepheids suggests that the interaction between radial and non-radial modes should also be investigated.

We present our analysis of Kepler observations of 29 RR Lyrae stars, based on 138-d of observation. We report precise pulsation periods for all stars. Nine of these stars had incorrect or unknown periods in the literature. Fourteen of the stars exhibit both amplitude and phase Blazhko modulations, with Blazhko periods ranging from 27.7 to more than 200 days. For V445 Lyr, a longer secondary variation is also observed in addition to its 53.2-d Blazhko period. The unprecedented precision of the Kepler photometry has led to the discovery of the the smallest modulations detected so far. Moreover, additional frequencies beyond the well-known harmonics and Blazhko multiplets have been found. These frequencies are located around the half-integer multiples of the main pulsation frequency for at least three stars. In four stars, these frequencies are close to the first and/or second overtone modes. The amplitudes of these periodicities seem to vary over the Blazhko cycle. V350 Lyr, a non-Blazhko star in our sample, is the first example of a double mode RR Lyrae star that pulsates in its fundamental and second overtone modes.

Double-mode Cepheid behaviour may arise either from non-resonant or resonant mode interaction. Inclusion of turbulent convection into pulsation codes by Kollath et al. (1998) led to stable double-mode F+1O Cepheid pulsation. However, our recent computations (Smolec & Moskalik 2008) indicate, that these models resulted from incorrect neglect of negative buoyancy effects. Once these effects are taken into account, only some resonant double-mode models could be found. These models do not solve the puzzle of double-mode phenomenon, as they are restricted to narrow parameter ranges. For majority of the observed double-mode Cepheids non-resonant mechanism has to be operational.

Data on multimode Cepheids from OGLE-III catalog of the LMC Cepheids are confronted with results from model calculations. Models whose radial mode periods are consistent with observation are not always in agreement with published evolutionary models. Nonradial mode interpretation is considered for the cases of unusual period ratios. The greatest challenge for stellar pulsation theory is explanation of double-mode pulsators with period ratios near 0.6.

V440 Per is a Population I Cepheid with the period of 7.57 day and low amplitude, almost sinusoidal light and radial velocity curves. With no reliable data on the 1st harmonic, its pulsation mode identification remained controversial. We obtained a radial velocity curve of V440 Per with our new high precision and high throughput Poznan Spectroscopic Telescope. Our data reach the accuracy of 130 m/s per individual measurement and yield a secure detection of the 1st harmonic with the amplitude of A_2= 140+/- 15 m/s. The velocity Fourier phase \phi_21 of V440 Per is inconsistent at the 7.25 \sigma level with those of the fundamental mode Cepheids, implying that the star must be an overtone Cepheid, as originally proposed by Kienzle et al.(1999). Thus, V440 Per becomes the longest period Cepheid with the securely established overtone pulsations. We show, that the convective nonlinear pulsation hydrocode can reproduce the Fourier parameters of V440 Per very well. Requirement to match the observed properties of V440 Per constrains free parameters of the dynamical convection model used in the pulsation calculations, in particular the radiative losses parameter.

One of the most interesting results from the OGLE-III study of the LMC Cepheids is the large number of objects that pulsate simultaneously in the first and second overtone (denoted 1O/2O). Double-mode Cepheids yield important constraint on stellar evolution models. We show that great majority of the LMC 1O/2O Cepheids have masses M=3.0+/-0.5 Msun. According to current stellar evolution calculations, these masses are lower than needed for the blue loop in the helium burning phase to reach the instability strip. On the other hand, we found most of these stars significantly overluminous if they are crossing the instability before helium ignition. A possible solution of this discrepancy is to allow for a large overshooting from the convective core in the main sequence phase. We also discuss origin of double-mode pulsation. At the short period range we find two types of resonances that are conducive to this form of pulsation. However, at longer periods, it has a different (non-resonant) origin.

The longstanding problem of modeling double-mode behaviour of classical pulsators was solved with the incorporation of turbulent convection into pulsation hydrocodes. However, the reasons for the computed double-mode behaviour were never clearly identified. In our recent papers (Smolec & Moskalik 2008a,b) we showed that the double-mode behaviour results from the neglect of negative buoyancy effects in some of the hydrocodes. If these effects are taken into account, no stable non-resonant double-mode behaviour can be found. In these proceedings we focus our attention on the role of negative buoyancy effects in classical Cepheid models.

For many years modeling of double-mode pulsation of classical pulsators was a challenging problem. Inclusion of turbulent convection into pulsation hydrocodes finally led to stable double-mode models. However, it was never analysed, which factor of turbulent convection is crucial. We show that the double-mode behaviour displayed in the computed models results from incorrect assumptions adopted in some of the pulsation hydrocodes, namely from the neglect of buoyant forces in convectively stable layers. This leads to significant turbulent energies and consequently to strong eddy-viscous damping in deep, convectively stable layers of the model. Resulting differential reduction of fundamental and first overtone amplitudes favours the occurrence of double-mode pulsation. Once buoyant forces in convectively stable regions are taken into account (as they should), no stable double-mode behaviour is found. The problem of modeling double-mode behaviour of classical pulsators remains open.

In this paper we describe our convective hydrocodes for radial stellar pulsation. We adopt the Kuhfuss (1986) model of convection, reformulated for the use in stellar pulsation hydrocodes. Physical as well as numerical assumptions of the code are described in detail. Described tests show, that our models are numerically robust and reproduce basic observational constraints. We discuss the effects of different treatment of some quantities in other pulsation hydrocodes. Our most important finding concerns the treatment of the turbulent source function in convectively stable regions. In our code we allow for negative values of source function in convectively stable zones, which reflects negative buoyancy. However, some authors restrict the source term to non-negative values. We show that this assumption leads to very high turbulent energies in convectively stable regions. The effect looks like overshooting, but it is not, because turbulence is generated by pulsations. Also, turbulent elements do not carry kinetic nor thermal energy, into convectively stable layers. The range of this artificial overshooting (as we shall call it) is as large as 6 local pressure scale heights, leading to unphysical internal damping through the eddy-viscous forces, in deep, convectively stable parts of the star.

As a by-product of beta Cephei and delta Cephei radiative model surveys we have found interesting cases of multimode pulsation. In case of beta Cephei stars, we found two multimode domains with two or three modes being involved. The origin of the multimode pulsation can be traced to one of the two different mechanisms: either to the non-resonant coupling of the two excited modes (double-mode pulsation) or to the 2omega_1 = omega_0 + omega_2 parametric resonance (triple-mode pulsation). In case of delta Cephei models, a triple-mode pulsation domain, connected with the 2:1 resonance, 2omega_0 = omega_2, is found. These models are of theoretical interest only, since they do not obey a proper mass-luminosity relation.

Although the driving mechanism acting in beta Cephei pulsators is well known (e.g. Dziembowski & Pamyatnykh 1993), problems concerning identification of amplitude limitation mechanism and non-uniform filling of the theoretical instability strip, remain to be solved. In the present analysis, these problems are addressed by non-linear modelling of radial pulsations of these stars. In this approach radial modes are treated as representative for all acoustic oscillations. Several models of different masses and metallicities were converged to limit cycles through Stellingwerf (1974) relaxation technique. Resulting peak-to-peak amplitudes are of order of DeltaV=0.3 mag. Such amplitudes are significantly larger than those observed in beta Cephei pulsators. Assuming that all acoustic modes are similar, we show that collective saturation of the driving mechanism by several acoustic modes can easily lower predicted saturation amplitudes to the observed level. Our calculations predict significant decrease of saturation amplitudes as we go to high mass/high luminosity models. However, this effect is not strong enough to explain scarcity of high mass beta Cephei variables. We also discuss robust double-mode behaviour, encountered in our radiative models. On a single evolutionary track we identify two double-mode domains with two different mechanisms resposible for double-mode behaviour. The non-resonant double-mode domain separates first overtone and fundamental mode pulsation domains. The resonant double-mode domain appears in the middle of the first overtone pulsation domain. Its origin can be traced to the 2omega_1=omega_0+omega_2 parametric resonance, which destabilizes the first overtone limit cycle.

Using nonlinear hydrocodes we have calculated single mode saturation amplitudes for beta Cephei models. Predicted saturation amplitudes are systematically higher than amplitudes observed in beta Cephei variables, even in monoperiodic ones. We argue that collective saturation of the pulsation instability by a dozen or so acoustic modes brings the theoretical amplitudes to the observed level.

The microlensing surveys, such as OGLE or MACHO, have led to the discovery of thousands of RR Lyrae stars in the Galactic Bulge and in the Magellanic Clouds, allowing for detailed investigation of these stars, specially the still mysterious Blazhko phenomenon. Higher incidence rate of Blazhko (BL) variables in the more metal-rich Galactic Bulge than in the LMC, suggests that occurance of Blazhko effect correlates with metallicity (Moskalik and Poretti 2003). To investigate this problem, we calibrate the photometric method of determining the metallicity of the RRab star (Kovacs and Zsoldos 1995) to the I-band and apply it to the OGLE Galactic Bulge and LMC data. In both systems, metallicities of non Blazhko and Blazhko variables are close to each other. The LMC Blazhko pulsators prefer slightly lower metallicities. The different metallicities of the Galactic Bulge and the LMC, can't explain the observed incidence rates. As a by-product of our metallicity estimates, we investigate the luminosity-metallicity relation, finding a steep dependence of the luminosity on [Fe/H]