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.

We collected over 6000 high-resolution spectra of four dozen field RR Lyrae (RRL) variables pulsating either in the fundamental (39 RRab) or in the first overtone (9 RRc) mode. We measured radial velocities (RVs) of four strong metallic and four Balmer lines along the entire pulsational cycle and derived RV amplitudes with accuracies better than 1$-$2~\kmsec. The new amplitudes were combined with literature data for 23~RRab and 3~RRc stars (total sample 74 RRLs) which allowed us to investigate the variation of the Bailey diagram (photometric amplitude versus period) when moving from optical to mid-infrared bands and to re-cast the Bailey diagram in terms of RV amplitudes. We found that RV amplitudes for RRab are minimally affected by nonlinear phenomena (shocks) and multi-periodicity (Blazhko effect). The RV slope ($\log P$--A(V$_r$)) when compared with the visual slope ($\log P$--A($V$)) is shallower and the dispersion, at fixed period, decreases by a factor of two. We constructed homogeneous sets of Horizontal Branch evolutionary models and nonlinear, convective pulsation models of RRLs to constrain the impact of evolutionary effects on their pulsation properties. Evolution causes, on the Bailey diagram based on RV amplitudes, a modest variation in pulsation period and a large dispersion in amplitude. The broad dispersion in period of the Bailey diagram is mainly caused by variation in RRL intrinsic parameters (stellar mass, chemical composition). Empirical evidence indicates that RV amplitudes are an optimal diagnostic for tracing the mean effective temperature across the RRab instability strip.

We present a kinematical study of 314 RR~Lyrae stars in the solar neighbourhood using the publicly available photometric, spectroscopic, and {\it Gaia} DR2 astrometric data to explore their distribution in the Milky Way. We report an overdensity of 22 RR~Lyrae stars in the solar neighbourhood at a pericenter distance of between 5--9\,kpc from the Galactic center. Their orbital parameters and their chemistry indicate that these 22 variables share the kinematics and the [Fe/H] values of the Galactic disc, with an average metallicity and tangential velocity of [Fe/H]=$-0.60$\,dex and $v_{\theta} = 241$\,km\,s$^{-1}$, respectively. From the distribution of the Galactocentric spherical velocity components, we find that these 22 disc-like RR~Lyrae variables are not consistent with the {\it Gaia} Sausage ({\it Gaia}-Enceladus), unlike almost half of the local RR~Lyrae stars. Chemical information from the literature shows that the majority of the selected pericenter peak RR~Lyrae variables are $\alpha$-poor, a property shared by typically much younger stars in the thin disc. Using the available photometry we rule out a possible misclassification with the known classical and anomalous Cepheids. The similar kinematic, chemical, and pulsation properties of these disc RR~Lyrae stars suggest they share a common origin. In contrast, we find the RR~Lyrae stars associated with the {\it Gaia}-Enceladus based on their kinematics and chemical composition show a considerable metallicity spread in the old population ($\sim$~1\,dex).

We present chemodynamical links between the present-day Milky Way halo field star population and Galactic globular clusters (GCs) using a dataset that combines information from the $\rm{\it Gaia}$ space mission and the Sloan Digital Sky Survey (SDSS-IV, DR14). Moreover, we incorporated a sample of halo giant stars with a distinct chemical signature (strong CN bandheads) that resembles the light-elements anomaly otherwise only seen in the second generation of GC stellar populations. Using three different tagging techniques, we could successfully establish unique associations between 151 extratidal stars in the neighborhood of eight GCs. In addition, we traced the possible origin of about $62\%$ of the sample of CN-strong giants to their potential host clusters. Several of the involved GCs have been brought into connection with the Gaia-Enceladus and Sequoia merger events. By establishing kinematic and chemical connections between 17 CN-strong stars and their surrounding fields, we could identify co-moving groups of stars at the same [Fe/H] with a possible cluster origin. We found strong evidence that four CN-strong stars and their associates are connected to the Sagittarius stream whilst their tightly confined [Fe/H] may hint to a birth site in M 54. Finally, we provide tentative estimates for the fraction of first-generation cluster stars among all stars lost to the halo. In the immediate cluster vicinity, this value amounts to $50.0\pm16.7\%$ while the associations in the halo field rather imply $80.2_{-5.2}^{+4.9}\%$. We speculate that -- if proven real by spectroscopic follow-up -- the disparity between these numbers could indicate a major contribution of low-mass clusters to the overall number of stars escaped to the halo or, alternatively, point toward a strong mass loss from the first generation during early cluster dissolution. [abridged]

The number of stars observed by the Optical Gravitational Lensing Experiment (OGLE) project in the Galactic bulge offers an invaluable chance to study RR Lyrae stars in a statistical manner. We used data of 3141 fundamental-mode RR Lyrae stars showing the Blazhko effect observed in OGLE-IV to investigate a possible connection between modulation amplitudes and periods, light curve and pulsation characteristics. We found that there is no simple monotonic correlation between any combination of two parameters concerning the Blazhko and pulsation amplitudes, periods and the shape of the light curves. There are only systematic limits. There is a bottom limit of the modulation period with respect to the pulsation period. We also found that the possible range of modulation amplitudes decreases with increasing pulsation period which could point towards that the Blazhko effect is suppressed in cooler, larger, more luminous and less metal abundant bulge RR Lyrae stars. Our investigation revealed that the distribution of the modulation periods can be described with two populations of stars with the mean modulation periods of 48 and 186 days. There is a certain region with a low density of the modulated stars, which we call the Blazhko valley, in the pulsation period-modulation period plane. Based on the similarity of the modulation envelopes, basically every star can be assigned to one of six morphological classes. Double modulation was found in 25 per cent of the studied stars. Only 6.3 per cent of modulated stars belongs to the Oosterhoff group II.

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 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 present 20 newly discovered candidates for binary systems with an RR~Lyrae companion. Using the photometric data from the Optical Gravitational Lensing Experiment (OGLE) and Korea Microlensing Telescope Network (KMTNet) for the Galactic bulge, we searched for binary candidates among non-modulated fundamental-mode RR~Lyrae variables. We constructed and investigated over 9000 $O-C$ diagrams for individual pulsators, in order to find signs of the light-travel time effect. We found 20 variables exhibiting cyclic variation in the $O-C$ (time of the observed minus calculated brightness maximum) diagram, which can be associated with a second component in the system, but for confirmation of binarity, a long-term radial velocity study is necessary. The modeling of the $O-C$ diagrams yields orbital parameters, estimates of the semi-amplitude of the radial velocity curve, and the mass function. The orbital periods of our candidates range from 3 to 15 years. Two of the binary candidates display a minimum mass higher than the assumed mass of the RR~Lyrae component, which points towards an evolved companion that can under some circumstances contribute significantly to the total light of the system.

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 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 present our results of searching for differences in light curves of modulated and non-modulated RRab stars in the Galactic bulge. We examined a sample of more than 8000 stars. The most important results are that Blazhko stars have shorter pulsation periods, less skewed mean light curves, lower mean amplitudes, larger rise-time, no difference in spatial distribution and metallicity.

We present the first release of the Czech Variable star catalogue that currently contains 1228 stars whose variability was discovered by 60 Czech observers. The catalogue contains confirmed variable stars of various types, but also candidates. We give precise coordinates, cross identification with other catalogues, information about constellation, variability type, brightness, light elements, name of the discoverer and year of discovery. In eighty-eight percent of stars the variability type is estimated, for more than 60 % of the stars the light ephemerides are given.

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 first paper of a series focused on the Blazhko effect in RR Lyrae type stars pulsating in the fundamental mode, that are located in the Galactic bulge. A~comprehensive overview about the incidence rate and light-curve characteristics of the Blazhko stars is given. We analysed 8\,282 stars having the best quality data in the OGLE-IV survey, and found that at least $40.3$\,\% of stars show modulation of their light curves. The number of Blazhko stars we identified is 3\,341, which is the largest sample ever studied implying the most relevant statistical results currently available. Using combined data sets with OGLE-III observations, we found that 50\,\% of stars that show unresolved close peaks to the main component in OGLE-IV are actually Blazhko stars with extremely long periods. Blazhko stars with modulation occur preferentially among RR Lyrae stars with shorter pulsation periods in the Galactic bulge. Fourier amplitude and phase coefficients based on the mean light curves appear to be substantially lower for Blazhko stars than for stars with unmodulated light curve in average. We derived new relations for the compatibility parameter $D_{m}$ in $I$ passband and relations that allow for differentiating modulated and non-modulated stars easily on the basis of $R_{31}$, $\phi_{21}$ and $\phi_{31}$. Photometric metallicities, intrinsic colours and absolute magnitudes computed using empirical relations are the same for Blazhko and non-modulated stars in the Galactic bulge suggesting no correlation between the occurrence of the Blazhko effect and these parameters.

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.

We investigated 1234 fundamental mode RR Lyrae stars observed by the All Sky Automated Survey (ASAS) to identify the Blazhko (BL) effect. A sample of 1547 BL stars from the literature was collected to compare the modulation-period distribution with stars newly identified in our sample. A classical frequency spectra analysis was performed using Period04 software. Data points from each star from the ASAS database were analysed individually to avoid confusion with artificial peaks and aliases. Statistical methods were used in the investigation of the modulation-period distribution. Altogether we identified 87 BL stars (48 new detections), 7 candidate stars, and 22 stars showing long-term period variations. The distribution of modulation periods of newly identified BL stars corresponds well to the distribution of modulation periods of stars located in the Galactic field, Galactic bulge, Large Magellanic Cloud, and globular cluster M5 collected from the literature. As a very important by-product of this comparison, we found that pulsation periods of BL stars follow Gaussian distribution with the mean period of $0.54\pm0.07$ d, while the modulation periods show log-normal distribution with centre at $\log(P_{\rm m}~{\rm [d]})=1.78\pm0.30$ dex. This means that 99.7 % of all known modulated stars have BL periods between 7.6 and 478 days. We discuss the identification of long modulation periods and show, that a significant percentage of stars showing long-term period variations could be classified as BL stars.

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.

A paper about variable stars with possible multiple occurrence in the VSX catalogue is presented. Our main criteria for identification of such duplicities were the angular distance among stars (below 1 arcmin) and close periods of objects. In our approach, we also considered double or half values of periods to reveal possible misclassification among stars with similar light curve shapes. The probability of false identification is expressed by the parameter R giving the relative difference between periods. We found 1487 pairs of stars in angular distance lower than 1 arcmin with period difference R lower than 0.1 %, which are high-probable candidates on duplicates. From this sample, 354 pairs have exactly the same periods (R = 0.0 %) and should be considered as definite duplicates. The main contribution of certain duplicates comes from the Catalina Sky Survey (73 pairs have two names with CSS acronym) and from the BEST projects (71 pairs). Distribution of identified duplicates on the sky is not homogeneous but contains surprising depression in Galactic plane.