Aiming to extend the census of RR Lyrae stars to highly reddened low-latitude regions of the central Milky Way, we performed a deep near-IR variability search using data from the VISTA Variables in the V\'ia L\'actea (VVV) survey of the bulge, analyzing the photometric time series of over a hundred million point sources. In order to separate fundamental-mode RR Lyrae (RRab) stars from other periodically variable sources, we trained a deep bidirectional long short-term memory recurrent neural network (RNN) classifier using VVV survey data and catalogs of RRab stars discovered and classified by optical surveys. Our classifier attained a ~99% precision and recall for light curves with signal-to-noise ratio above 60, and is comparable to the best-performing classifiers trained on accurate optical data. Using our RNN classifier, we identified over 4300 hitherto unknown bona fide RRab stars toward the inner bulge. We provide their photometric catalog and VVV J,H,Ks photometric time-series.

Prompted by some inconsistencies in the photometry of the VISTA Variables in the V\'ia L\'actea (VVV) survey, we conduct a revision of the standard calibration procedure of VISTA data in the $J$, $H$, and $K_S$ passbands. Two independent sources of bias in the photometric zero-points are identified: First, high sky backgrounds severely affect the $H$-band measurements, but this can mostly be minimized by strict data vetting. Secondly, during the zero-point calibration, stars serving as standards are taken from the 2MASS catalog, which can suffer from high degrees of blending in regions of high stellar density, affecting both the absolute photometric calibration, as well as the scatter of repeated observations. The former affects studies that rely on an absolute magnitude scale, while the latter can also affect the shapes and amount of scatter in the VVV light curves, thus potentially hampering their proper classification. We show that these errors can be effectively eliminated by relatively simple modifications of the standard calibration procedure, and demonstrate the effect of the recalibration on the VVV survey's data quality. We give recommendations for future improvements of the pipeline calibration of VISTA photometry, while also providing preliminary corrections to the VVV $JHK_\mathrm{S}$ observations as a temporary measure.

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 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.

The far side of the Milky Way's disk is one of the most concealed parts of the known Universe due to extremely high interstellar extinction and point source density toward low Galactic latitudes. Large time-domain photometric surveys operating in the near-infrared hold great potential for the exploration of these vast uncharted areas of our Galaxy. We conducted a census of distant classical and type II Cepheids along the southern Galactic mid-plane using near-infrared photometry from the VISTA Variables in the V\'ia L\'actea survey. We performed a machine-learned classification of the Cepheids based on their infrared light curves using a convolutional neural network. We have discovered 640 distant classical Cepheids with up to ~40 magnitudes of visual extinction, and over 500 type II Cepheids, most of them located in the inner bulge. Intrinsic color indices of individual Cepheids were predicted from sparse photometric data using a neural network, allowing their use as accurate reddening tracers. They revealed a steep, spatially varying near-infrared extinction curve toward the inner bulge. Type II Cepheids in the Galactic bulge were also employed to measure robust mean selective-to-absolute extinction ratios. They trace a centrally concentrated spatial distribution of the old bulge population with a slight elongation, consistent with earlier results from RR Lyrae stars. Likewise, the classical Cepheids were utilized to trace the Galactic warp and various substructures of the Galactic disk, and to uncover significant vertical and radial age gradients of the thin disk population at the far side of the Milky Way.

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 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.

This study's objective was to exploit infrared VVV (VISTA Variables in the Via Lactea) photometry for high latitude RRab stars to establish an accurate Galactic Centre distance. RRab candidates were discovered and reaffirmed ($n=4194$) by matching $K_s$ photometry with templates via $\chi^2$ minimization, and contaminants were reduced by ensuring targets adhered to a strict period-amplitude ($\Delta K_s$) trend and passed the Elorietta et al. classifier. The distance to the Galactic Centre was determined from a high latitude Bulge subsample ($|b|>4^{o}$, $R_{GC}=8.30 \pm 0.36$ kpc, random uncertainty is relatively negligible), and importantly, the comparatively low color-excess and uncrowded location mitigated uncertainties tied to the extinction law, the magnitude-limited nature of the analysis, and photometric contamination. Circumventing those problems resulted in a key uncertainty being the $M_{K_s}$ relation, which was derived using LMC RRab stars ($M_{K_s}=-(2.66\pm0.06) \log{P}-(1.03\pm0.06)$, $(J-K_s)_0=(0.31\pm0.04) \log{P} + (0.35\pm0.02)$, assuming $\mu_{0,LMC}=18.43$). The Galactic Centre distance was not corrected for the cone-effect. Lastly, a new distance indicator emerged as brighter overdensities in the period-magnitude-amplitude diagrams analyzed, which arise from blended RRab and red clump stars. Blending may thrust faint extragalactic variables into the range of detectability.

RR Lyrae stars (RRLs) are tracers of the Milky Way's fossil record, holding valuable information on its formation and early evolution. Owing to the high interstellar extinction endemic to the Galactic plane, distant RRLs lying at low Galactic latitudes have been elusive. We attained a census of 1892 high-confidence RRLs by exploiting the near-infrared photometric database of the VVV survey's disk footprint spanning $\sim$70$^\circ$ of Galactic longitude, using a machine-learned classifier. Novel data-driven methods were employed to accurately characterize their spatial distribution using sparsely sampled multi-band photometry. The RRL metallicity distribution function (MDF) was derived from their $K_s$-band light curve parameters using machine-learning methods. The MDF shows remarkable structural similarities to both the spectroscopic MDF of red clump giants and the MDF of bulge RRLs. We model the MDF with a multi-component density distribution and find that the number density of stars associated with the different model components systematically changes with both the Galactocentric radius and vertical distance from the Galactic plane, equivalent to weak metallicity gradients. Based on the consistency with results from the ARGOS survey, three MDF modes are attributed to the old disk populations, while the most metal-poor RRLs are probably halo interlopers. We propose that the dominant [Fe/H] component with a mean of $-1$ dex might correspond to the outskirts of an ancient Galactic spheroid or classical bulge component residing in the central Milky Way. The physical origins of the RRLs in this study need to be verified by kinematical information.

RR~Lyrae variables are widely used tracers of Galactic halo structure and kinematics, but they can also serve to constrain the distribution of the old stellar population in the Galactic bulge. With the aim of improving their near-infrared photometric characterization, we investigate their near-infrared light curves, as well as the empirical relationships between their light curve and metallicities using machine learning methods. We introduce a new, robust method for the estimation of the light-curve shapes, and hence the average magnitudes of RR~Lyrae variables in the $K_\mathrm{S}$ band, by utilizing the first few principal components (PCs) as basis vectors, obtained from the PC analysis of a training set of light curves. Furthermore, we use the amplitudes of these PCs to predict the light-curve shape of each star in the $J$-band, allowing us to precisely determine their average magnitudes (hence colors), even in cases where only one $J$ measurement is available. Finally, we demonstrate that the $K_\mathrm{S}$-band light-curve parameters of RR~Lyrae variables, together with the period, allow the estimation of the metallicity of individual stars with an accuracy of $\sim 0.2 - 0.25$ dex, providing valuable chemical information about old stellar populations bearing RR~Lyrae variables. The methods presented here can be straightforwardly adopted for other classes of variable stars, bands, or for the estimation of other physical quantities.

We present first direct evidence of modulation in the $K$-band of Blazhko-type RR Lyrae stars that are identified by their secular modulations in the I-band data of OGLE-IV. A method has been developed to decompose the $K$-band light variation into two parts originating from the temperature and the radius changes using synthetic data of atmosphere-model grids. The amplitudes of the temperature and the radius variations derived from the method for non-Blazhko RRab stars are in very good agreement with the results of the Baade-Wesselink analysis of RRab stars in the M3 globular cluster confirming the applicability and correctness of the method. It has been found that the Blazhko modulation is primarily driven by the change in the temperature variation. The radius variation plays a marginal part, moreover it has an opposite sign as if the Blazhko effect was caused by the radii variations. This result reinforces the previous finding based on the Baade-Wesselink analysis of M3 (NGC 5272) RR Lyrae, that significant modulation of the radius variations can only be detected in radial-velocity measurements, which relies on spectral lines that form in the uppermost atmospheric layers. Our result gives the first insight into the energetics and dynamics of the Blazhko phenomenon, hence it puts strong constraints on its possible physical explanations.

$\omega$ Centauri (NGC~5139) contains many variable stars of different types, including the pulsating type II Cepheids, RR Lyrae and SX Phoenicis stars. We carried out a deep, wide-field, near-infrared (IR) variability survey of $\omega$ Cen, using the VISTA telescope. We assembled an unprecedented homogeneous and complete $J$ and $K_{\rm S}$ near-IR catalog of variable stars in the field of $\omega$ Cen. In this paper we compare optical and near-IR light curves of RR Lyrae stars, emphasizing the main differences. Moreover, we discuss the ability of near-IR observations to detect SX Phoenicis stars given the fact that the amplitudes are much smaller in these bands compared to the optical. Finally, we consider the case in which all the pulsating stars in the three different variability types follow a single period-luminosity relation in the near-IR bands.

$\omega$ Centauri (NGC 5139) hosts hundreds of pulsating variable stars of different types, thus representing a treasure trove for studies of their corresponding period-luminosity (PL) relations. Our goal in this study is to obtain the PL relations for RR Lyrae, and SX Phoenicis stars in the field of the cluster, based on high-quality, well-sampled light curves in the near-infrared (IR). $\omega$ Centauri was observed using VIRCAM mounted on VISTA. A total of 42 epochs in $J$ and 100 epochs in $K_{\rm S}$ were obtained, spanning 352 days. Point-spread function photometry was performed using DoPhot and DAOPHOT in the outer and inner regions of the cluster, respectively. Based on the comprehensive catalogue of near-IR light curves thus secured, PL relations were obtained for the different types of pulsators in the cluster, both in the $J$ and $K_{\rm S}$ bands. This includes the first PL relations in the near-IR for fundamental-mode SX Phoenicis stars. The near-IR magnitudes and periods of Type II Cepheids and RR Lyrae stars were used to derive an updated true distance modulus to the cluster, with a resulting value of $(m-M)_0 = 13.708 \pm 0.035 \pm 0.10$ mag, where the error bars correspond to the adopted statistical and systematic errors, respectively. Adding the errors in quadrature, this is equivalent to a heliocentric distance of $5.52\pm 0.27$ kpc.

We use deep multi-epoch near-IR images of the VISTA Variables in the V\'ia L\'actea (VVV) Survey to search for RR Lyrae stars towards the Southern Galactic plane. Here we report the discovery of a group of RR Lyrae stars close together in VVV tile d025. Inspection of the VVV images and PSF photometry reveals that most of these stars are likely to belong to a globular cluster, that matches the position of the previously known star cluster FSR\,1716. The stellar density map of the field yields a $>100$ sigma detection for this candidate globular cluster, that is centered at equatorial coordinates $RA_{J2000}=$16:10:30.0, $DEC_{J2000}=-$53:44:56; and galactic coordinates $l=$329.77812, $b=-$1.59227. The color-magnitude diagram of this object reveals a well populated red giant branch, with a prominent red clump at $K_s=13.35 \pm 0.05$, and $J-K_s=1.30 \pm 0.05$. We present the cluster RR Lyrae positions, magnitudes, colors, periods and amplitudes. The presence of RR Lyrae indicates an old globular cluster, with age $>10$ Gyr. We classify this object as an Oosterhoff type I globular cluster, based on the mean period of its RR Lyrae type ab, $<P>=0.540$ days, and argue that this is a relatively metal-poor cluster with $[Fe/H] = -1.5 \pm 0.4$ dex. The mean extinction and reddening for this cluster are $A_{K_s}=0.38 \pm 0.02$, and $E(J-K_s)=0.72 \pm 0.02$ mag, respectively, as measured from the RR Lyrae colors and the near-IR color-magnitude diagram. We also measure the cluster distance using the RR Lyrae type ab stars. The cluster mean distance modulus is $(m-M)_0 = 14.38 \pm 0.03$ mag, implying a distance $D = 7.5 \pm 0.2$ kpc, and a Galactocentric distance $R_G=4.3$ kpc.

Deep near-IR images from the VISTA Variables in the V\'ia L\'actea (VVV) Survey were used to search for RR Lyrae stars in the Southern Galactic plane. A sizable sample of 404 RR Lyrae of type ab stars was identified across a thin slice of the 4$^{\rm th}$ Galactic quadrant ($295\deg < l < 350\deg$, $-2.24\deg < b < -1.05\deg$). The sample's distance distribution exhibits a maximum density that occurs at the bulge tangent point, which implies that this primarily Oosterhoff type I population of RRab stars does not trace the bar delineated by their red clump counterparts. The bulge RR Lyrae population does not extend beyond $l \sim340 \deg$, and the sample's spatial distribution presents evidence of density enhancements and substructure that warrants further investigation. Indeed, the sample may be employed to evaluate Galactic evolution models, and is particularly lucrative since half of the discovered RR Lyrae are within reach of Gaia astrometric observations.

Variable stars of RR Lyrae type are a prime tool to obtain distances to old stellar populations in the Milky Way, and one of the main aims of the Vista Variables in the Via Lactea (VVV) near-infrared survey is to use them to map the structure of the Galactic Bulge. Due to the large number of expected sources, this requires an automated mechanism for selecting RR Lyrae,and particularly those of the more easily recognized type ab (i.e., fundamental-mode pulsators), from the 10^6-10^7 variables expected in the VVV survey area. In this work we describe a supervised machine-learned classifier constructed for assigning a score to a K_s-band VVV light curve that indicates its likelihood of being ab-type RR Lyrae. We describe the key steps in the construction of the classifier, which were the choice of features, training set, selection of aperture and family of classifiers. We find that the AdaBoost family of classifiers give consistently the best performance for our problem, and obtain a classifier based on the AdaBoost algorithm that achieves a harmonic mean between false positives and false negatives of ~7% for typical VVV light curve sets. This performance is estimated using cross-validation and through the comparison to two independent datasets that were classified by human experts.

The Quintuplet cluster is one of the most massive star clusters in the Milky Way, situated very close to the Galactic center. We present a new search for variable stars in the vicinity of the cluster, using the five-year database of the Vista Variables in the Via Lactea (VVV) ESO Public Survey in the near-infrared. A total of 7586 objects were identified in the zone around $2'$ from the cluster center, using 55 $K_S$-band epochs. Thirty-three stars show $K_S$-band variability, 24 of them being previously undiscovered. Most of the variable stars found are slow/semiregular variables, long-period variables of the Mira type, and OH/IR stars. In addition, a good number of our candidates show variations in a rather short timescale. We also propose four Young Stellar Object (YSO) candidates, which could be cluster members.

This project is a massive near-infrared (NIR) search for variable stars in highly reddened and obscured open cluster (OC) fields projected on regions of the Galactic bulge and disk. The search is performed using photometric NIR data in the $J$-, $H$- and $K_s$- bands obtained from the Vista Variables in the V\'ia L\'actea (VVV) Survey. We performed in each cluster field a variability search using Stetson's variability statistics to select the variable candidates. Later, those candidates were subjected to a frequency analysis using the Generalized Lomb-Scargle and the Phase Dispersion Minimization algorithms. The number of independent observations range between 63 and 73. The newly discovered variables in this study, 157 in total in three different known OCs, are classified based on their light curve shapes, periods, amplitudes and their location in the corresponding color-magnitude $(J-K_s,K_s)$ and color-color $(H-K_s,J-H)$ diagrams. We found 5 possible Cepheid stars which, based on the period-luminosity relation, are very likely type II Cepheids located behind the bulge. Among the newly discovered variables, there are eclipsing binaries, $\delta$ Scuti, as well as background RR Lyrae stars. Using the new version of the Wilson & Devinney code as well as the "Physics Of Eclipsing Binaries" (PHOEBE) code, we analyzed some of the best eclipsing binaries we discovered. Our results show that these studied systems turn out to be ranging from detached to double-contact binaries, with low eccentricities and high inclinations of approximately $80^{\circ}$. Their surface temperatures range between $3500$K and $8000$K.

Solid insight into the physics of the inner Milky Way is key to understanding our Galaxy's evolution, but extreme dust obscuration has historically hindered efforts to map the area along the Galactic mid-plane. New comprehensive near-infrared time-series photometry from the VVV Survey has revealed 35 classical Cepheids, tracing a previously unobserved component of the inner Galaxy, namely a ubiquitous inner thin disk of young stars along the Galactic mid-plane, traversing across the bulge. The discovered period (age) spread of these classical Cepheids implies a continuous supply of newly formed stars in the central region of the Galaxy over the last 100 million years.

Despite their importance, very few RR Lyrae (RRL) stars have been known to reside in binary systems. We report on a search for binary RRL in the OGLE-III Galactic bulge data. Our approach consists in the search for evidence of the light-travel time effect in so-called observed minus calculated ($O-C$) diagrams. Analysis of 1952 well-observed fundamental-mode RRL in the OGLE-III data revealed an initial sample of 29 candidates. We used the recently released OGLE-IV data to extend the baselines up to 17 years, leading to a final sample of 12 firm binary candidates. We provide $O-C$ diagrams and binary parameters for this final sample, and also discuss the properties of 8 additional candidate binaries whose parameters cannot be firmly determined at present. We also estimate that $\gtrsim 4$ per cent of the RRL reside in binary systems.

[ABRIDGED] $\omega$ Centauri (NGC 5139) contains large numbers of variable stars of different types and, in particular, more than a hundred RR Lyrae stars. We have conducted a variability survey of $\omega$ Cen in the NIR, using ESO's 4.1m Visible and Infrared Survey Telescope for Astronomy (VISTA). This is the first paper of a series describing our results. $\omega$ Cen was observed using VIRCAM mounted on VISTA. A total of 42 and 100 epochs in $J$ and $K_{\rm S}$, respectively, were obtained, distributed over a total timespan of 352 days. PSF photometry was performed, and periods of the known variable stars were improved when necessary using an ANOVA analysis. An unprecedented homogeneous and complete NIR catalogue of RR Lyrae stars in the field of $\omega$ Cen was collected, allowing us to study, for the first time, all the RR Lyrae stars associated to the cluster, except 4 located far away from the cluster center. Membership status, subclassifications between RRab and RRc subtypes, periods, amplitudes, and mean magnitudes were derived for all the stars in our sample. Additionally, 4 new RR Lyrae stars were discovered, 2 of them with high probability of being cluster members. The distribution of $\omega$ Cen stars in the Bailey (period-amplitude) diagram is also discussed. Reference lines in this plane, for both Oosterhoff type I (OoI) and II (OoII) components, are provided. In the present paper, we clarify the status of many (candidate) RR Lyrae stars that had been unclear in previous studies. This includes stars with anomalous positions in the color-magnitude diagram, uncertain periods or/and variability types, and possible field interlopers. We conclude that $\omega$ Cen hosts a total of 88 RRab and 101 RRc stars, for a grand total of 189 likely members. We confirm that most RRab stars in the cluster belong to an OoII component, as previously found using visual data.

The VISTA Variables in the V\'ia L\'actea (VVV) Survey is one of the six ESO public surveys currently ongoing at the VISTA telescope on Cerro Paranal, Chile. VVV uses near-IR ($ZYJHK_{\rm s}$) filters that at present provide photometry to a depth of $K_{\rm s} \sim 17.0$ mag in up to 36 epochs spanning over four years, and aim at discovering more than 10$^6$ variable sources as well as trace the structure of the Galactic bulge and part of the southern disk. A variability search was performed to find RR Lyrae variable stars. The low stellar density of the VVV tile $\textit{b201}$, which is centered at ($\ell, b$) $\sim$ ($-9^\circ, -9^\circ$), makes it suitable to search for variable stars. Previous studies have identified some RR Lyrae stars using optical bands that served to test our search procedure. The main goal is to measure the reddening, interstellar extinction, and distances of the RR Lyrae stars and to study their distribution on the Milky Way bulge. A total of 1.5 sq deg were analyzed, and we found 39 RR Lyrae stars, 27 of which belong to the ab-type and 12 to the c-type. Our analysis recovers all the previously identified RR Lyrae variables in the field and discovers 29 new RR Lyrae stars. The reddening and extinction toward all the RRab stars in this tile were derived, and distance estimations were obtained through the period--luminosity relation. Despite the limited amount of RR Lyrae stars studied, our results are consistent with a spheroidal or central distribution around $\sim 8.1$ and $\sim 8.5$ kpc. for either the Cardelli or Nishiyama extinction law.

We report the discovery of a pair of extremely reddened classical Cepheid variable stars located in the Galactic plane behind the bulge, using near-infrared time-series photometry from the VVV Survey. This is the first time that such objects have ever been found in the opposite side of the Galactic plane. The Cepheids have almost identical periods, apparent brightnesses and colors. From the near-infrared Leavitt law, we determine their distances with ~1.5% precision and ~8% accuracy. We find that they have a same total extinction of A(V)~32 mag, and are located at the same heliocentric distance of <d>=11.4+/-0.9 kpc, and less than 1 pc from the true Galactic plane. Their similar periods indicate that the Cepheids are also coeval, with an age of ~48+/-3 Myr, according to theoretical models. They are separated by an angular distance of only 18.3", corresponding to a projected separation of ~1 pc. Their position coincides with the expected location of the Far 3 kpc Arm behind the bulge. Such a tight pair of similar classical Cepheids indicates the presence of an underlying young open cluster, that is both hidden behind heavy extinction and disguised by the dense stellar field of the bulge. All our attempts to directly detect this "invisible cluster" have failed, and deeper observations are needed.

The VISTA Variables in the Via Lactea (VVV) ESO Public Survey is opening a new window to study the inner Galactic globular clusters using their variable stars. These globular clusters have been neglected in the past due to the difficulties caused by the presence of an elevated extinction and high field stellar densities in their lines of sight. However, the discovery and study of any present variables in these clusters, especially RRLyrae stars, can help to greatly improve the accuracy of their physical parameters. It can also help to shed some light on the interrogations brought by the intriguing Oosterhoff dichotomy in the Galactic globular cluster system. In a series of papers we plan to explore the variable stars in the globular clusters falling inside the field of the VVV survey. In this first paper we search and study the variables present in two highly-reddened, moderately metal-poor, faint, inner Galactic globular clusters: 2MASS-GC02 and Terzan10. We report the discovery of sizable populations of RR Lyrae stars in both globular clusters. We use near-infrared period-luminosity relations to determine the color excess of each RR Lyrae star, from which we obtain both accurate distances to the globular clusters and the ratios of the selective to total extinction in their directions. We find the extinction towards both clusters to be elevated, non-standard, and highly differential. We also find both clusters to be closer to the Galactic center than previously thought, with Terzan10 being on the far side of the Galactic bulge. Finally, we discuss their Oosterhoff properties, and conclude that both clusters stand out from the dichotomy followed by most Galactic globular clusters.

Stellar variability in the near-infrared (NIR) remains largely unexplored. The exploitation of public science archives with data-mining methods offers a perspective for the time-domain exploration of the NIR sky. We perform a comprehensive search for stellar variability using the optical-NIR multi-band photometric data in the public Calibration Database of the WFCAM Science Archive (WSA), with the aim of contributing to the general census of variable stars, and to extend the current scarce inventory of accurate NIR light curves for a number of variable star classes. We introduce new variability indices designed for multi-band data with correlated sampling, and apply them for pre-selecting variable star candidates, i.e., light curves that are dominated by correlated variations, from noise-dominated ones. Pre-selection criteria are established by robust numerical tests for evaluating the response of variability indices to colored noise characteristic to the data. We find 275 periodic variable stars and an additional 44 objects with suspected variability with uncertain periods or apparently aperiodic variation. Only 44 of these objects had been previously known, including 11 RR~Lyrae stars in the outskirts of the globular cluster M3 (NGC~5272). We provide a preliminary classification of the new variable stars that have well-measured light curves, but the variability types of a large number of objects remain ambiguous. We classify most of the new variables as contact binary stars, but we also find several pulsating stars, among which 34 are probably new field RR~Lyrae and 3 are likely Cepheids. We also identify 32 highly reddened variable objects close to previously known dark nebulae, suggesting that these are embedded young stellar objects. We publish our results and all light-curve data as the WFCAM Variable Star Catalog.

The VISTA Variables in the V\'ia L\'actea (VVV) survey is one of six public ESO surveys, and is now in its 4th year of observing. Although far from being complete, the VVV survey has already delivered many results, some directly connected to the intended science goals (detection of variables stars, microlensing events, new star clusters), others concerning more exotic objects, e.g. novae. Now, at the end of the fourth observing period, and comprising roughly 50% of the proposed observations, the actual status of the survey, as well some of the results based on the VVV data, are presented.

Context. The Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey is a variability survey of the Milky Way bulge and an adjacent section of the disk carried out from 2010 on ESO Visible and Infrared Survey Telescope for Astronomy (VISTA). VVV will eventually deliver a deep near-IR atlas with photometry and positions in five passbands (ZYJHK_S) and a catalogue of 1-10 million variable point sources - mostly unknown - which require classifications. Aims. The main goal of the VVV Templates Project, that we introduce in this work, is to develop and test the machine-learning algorithms for the automated classification of the VVV light-curves. As VVV is the first massive, multi-epoch survey of stellar variability in the near-infrared, the template light-curves that are required for training the classification algorithms are not available. In the first paper of the series we describe the construction of this comprehensive database of infrared stellar variability. Methods. First we performed a systematic search in the literature and public data archives, second, we coordinated a worldwide observational campaign, and third we exploited the VVV variability database itself on (optically) well-known stars to gather high-quality infrared light-curves of several hundreds of variable stars. Results. We have now collected a significant (and still increasing) number of infrared template light-curves. This database will be used as a training-set for the machine-learning algorithms that will automatically classify the light-curves produced by VVV. The results of such an automated classification will be covered in forthcoming papers of the series.

The Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey is an ongoing time-series, near-infrared (IR) survey of the Galactic bulge and an adjacent portion of the inner disk, covering 562 square degrees of the sky, using ESO's VISTA telescope. The survey has provided superb multi-color photometry in 5 broadband filters ($Z$, $Y$, $J$, $H$, and $K_s$), leading to the best map of the inner Milky Way ever obtained, particularly in the near-IR. The main variability part of the survey, which is focused on $K_s$-band observations, is currently underway, with bulge fields having been observed between 31 and 70 times, and disk fields between 17 and 36 times. When the survey is complete, bulge (disk) fields will have been observed up to a total of 100 (60) times, providing unprecedented depth and time coverage. Here we provide a first overview of stellar variability in the VVV data, including examples of the light curves that have been collected thus far, scientific applications, and our efforts towards the automated classification of VVV light curves.

We have combined optical and near-infrared data of known RR Lyrae (RRL) stars in the bulge in order to study the spatial distribution of its metal-poor component by measuring precise reddening values and distances of 7663 fundamental-mode RRL stars with high-quality photometry. We obtain a distance to the Galactic center of R0 = 8.33 +/- 0.05 +/- 0.14 kpc. We find that the spatial distribution of the RRL stars differs from the structures traced by the predominantly metal-rich red clump (RC) stars. Unlike the RC stars, the RRL stars do not trace a strong bar, but have a more spheroidal, centrally concentrated distribution, showing only a slight elongation in its very center. We find a hint of bimodality in the density distribution at high southern latitudes (b < -5), which needs to be confirmed by extending the areal coverage of the current census. The different spatial distributions of the metal-rich and metal-poor stellar populations suggest that the Milky Way has a composite bulge.

The VVV is currently surveying the central regions of the Milky Way at near-infrared wavelengths, including 36 known Galactic globular clusters. We already have data in several filters to generate their color-magnitude diagrams. We also have enough epochs to begin producing the light-curves of the cluster stars and look for any possible member variable stars. We are especially interested in RR Lyrae stars, since they are abundant in Galactic globular clusters, and the period-luminosity relation they show at near-infrared wavelengths, can help enormously in providing accurate distances and reddenings for the most extincted and poorly studied inner Galactic globular clusters. We center our attention here in the preliminary analysis of three of these clusters: NGC 6441, Terzan 10 and 2MASS-GC02.

The putative association between the 10.8 d classical Cepheid TW Nor and the open cluster Lynga 6 has generated considerable debate in the literature. New JHKs photometry in tandem with existing radial velocities for Lynga 6 stars imply cluster membership for TW Nor, and establish the variable as a high-weight calibrator for classical Cepheid relations. Fundamental mean parameters determined for Lynga 6 are: d=1.91+-0.10 kpc, E(J-H)=0.38+-0.02, and log(t)=7.9+-0.1. The Benedict et al.(2007)/Turner (2010) Galactic VIc Wesenheit function was revised using TW Nor's new parameters: W_VI,0=(-3.37+-0.08)log(P_0)-2.48+-0.08. TW Nor/Lynga 6 lie beyond the Sagittarius-Carina spiral arm and occupy the Centaurus arm, along with innumerable young Cepheids and clusters (e.g., VW Cen & VVV CL070).

Thorough analysis of the multicolour CCD observations of the RRab-type variable, CZ Lacertae is presented. The observations were carried out in two consecutive observing seasons in 2004 and 2005 within the framework of the Konkoly Blazhko Survey of bright, northern, short-period RRab variables. The O-C variation of CZ Lac indicated that a significant period decrease took place just around the time of the CCD observations. Our data gave a unique opportunity to study the related changes in the pulsation and modulation properties of a Blazhko star in detail. Two different period components ($\approx$14.6 d and $\approx$18.6 d) of the Blazhko modulation were identified. Both modulation components had similar strength. The periods and amplitudes of the modulations changed significantly from the first season to the next, meanwhile, the mean pulsation amplitude slightly decreased. The modulation frequencies were in a 5:4 resonance ratio in the first observing season then the frequencies shifted in opposite directions, and their ratio was close to the 4:3 resonance in the next season. The interaction of the two modulations caused beating with a period of 74~d in the first season, which resembled the 4-yr-long cycle of the $\approx$40-d modulation of RR Lyr. The mean values of the global physical parameters and their changes with Blazhko phase of both modulation components were determined by the Inverse Photometric Method.

We present more than 1000-day long photometry of EY Draconis in BV(RI)C passbands. The changes in the light curve are caused by the spottedness of the rotating surface. Modelling of the spotted surface shows that there are two large active regions present on the star on the opposite hemispheres. The evolution of the surface patterns suggests a flip-flop phenomenon. Using Fourier analysis, we detect a rotation period of P_rot=0.45875d, and an activity cycle with P~350d, similar to the 11-year long cycle of the Sun. This cycle with its year-long period is the shortest one ever detected on active stars. Two bright flares are also detected and analysed.

We describe the public ESO near-IR variability survey (VVV) scanning the Milky Way bulge and an adjacent section of the mid-plane where star formation activity is high. The survey will take 1929 hours of observations with the 4-metre VISTA telescope during five years (2010-2014), covering ~10^9 point sources across an area of 520 deg^2, including 33 known globular clusters and ~350 open clusters. The final product will be a deep near-IR atlas in five passbands (0.9-2.5 microns) and a catalogue of more than 10^6 variable point sources. Unlike single-epoch surveys that, in most cases, only produce 2-D maps, the VVV variable star survey will enable the construction of a 3-D map of the surveyed region using well-understood distance indicators such as RR Lyrae stars, and Cepheids. It will yield important information on the ages of the populations. The observations will be combined with data from MACHO, OGLE, EROS, VST, Spitzer, HST, Chandra, INTEGRAL, WISE, Fermi LAT, XMM-Newton, GAIA and ALMA for a complete understanding of the variable sources in the inner Milky Way. This public survey will provide data available to the whole community and therefore will enable further studies of the history of the Milky Way, its globular cluster evolution, and the population census of the Galactic Bulge and center, as well as the investigations of the star forming regions in the disk. The combined variable star catalogues will have important implications for theoretical investigations of pulsation properties of stars.

The analysis of the multicolour photometric observations of MW Lyr, a large modulation amplitude Blazhko variable, shows for the first time how the mean global physical parameters vary during the Blazhko cycle. About 1-2 percent changes in the mean radius, luminosity and surface effective temperature are detected. The mean radius and temperature changes are in good accordance with pulsation model results, which show that these parameters do indeed vary within this order of magnitude if the amplitude of the pulsation changes significantly. We interpret the phase modulation of the pulsation to be a consequence of period changes. Its magnitude corresponds exactly what one expects from the detected changes of the mean radius assuming that the pulsation constant remains the same during the modulation. Our results indicate that during the modulation the pulsation remains purely radial, and the underlying mechanism is most probably a periodic perturbation of the stellar luminosity with the modulation period.

We have obtained the most extensive and most accurate photometric data of a Blazhko variable MW Lyr during the 2006-2007 observing seasons. The data within each 0.05 phase bin of the modulation period (P_m=1/f_m) cover the entire light cycle of the primary pulsation period (P_0=1/f_0), making possible a very rigorous and complete analysis. The modulation period is found to be 16.5462 d, which is about half of that was reported earlier from visual observations. Previously unknown features of the modulation have been detected. Besides the main modulation frequency f_m, sidelobe modulation frequencies around the pulsation frequency and its harmonics appear at +/- 2 f_m, +/- 4 f_m, and +/- 12.5 f_m separations as well. Residual signals in the prewhitened light curve larger than the observational noise appear at the minimum-rising branch-maximum phase of the pulsation, which most probably arise from some stochastic/chaotic behaviour of the pulsation/modulation. The Fourier parameters of the mean light curve differ significantly from the averages of the Fourier parameters of the observed light curves in the different phases of the Blazhko cycle. Consequently, the mean light curve of MW Lyrae never matches its actual light variation. The Phi_21, Phi_31 phase differences in different phases of the modulation show unexpected stability during the Blazhko cycle. A new phenomenological description of the light curve variation is defined that separates the amplitude and phase (period) modulations utilising the phase coherency of the lower order Fourier phases.

We present the frequency analysis of the multicolour time series photometry of the field RRd variable BS Comae. The large number of data points in each of the BV(RI)_c bands and the ~0.01 magnitude accuracy of the individual measurements allow us a high precision analysis of the properties of the combination frequencies due to nonlinear coupling. Through the combination of the frequency spectra in different colors we show that except for the components corresponding to the linear combinations of the two pulsation modes, there are no other components present above the millimagnitude amplitude level.

Extended BV(RI)c CCD observations of SS Cnc, a short period RRab star are presented. Nearly 1400 data points in each band have been obtained spanning over 79 days during the spring of 2005. The star exhibits light curve modulation, the so called Blazhko effect with small amplitude (B maximum brightness varies 0.1 mag) and with the shortest modulation period (5.309 d) ever observed. In the Fourier spectrum of the V light curve the pulsation frequency components are detected up to the 24th harmonic order, and modulation side lobe frequencies with significantly asymmetric amplitudes are seen up to the 15th and 9th orders for the lower and higher frequency components, respectively. Detailed comparison of the modulation behavior of SS Cnc and RR Gem, the two recently discovered small amplitude, short modulation period Blazhko stars is presented. The modulation frequency (f_m) appears in the Fourier spectrum of both stars with similar amplitude. We also demonstrate that the modulation frequencies have basically different properties as the pulsation and modulation side lobe frequencies have, indicating that the physics behind these frequency components are not the same. The discovery of small amplitude modulations of RRab stars cautions that the large photometric surveys (MACHO, OGLE) may seriously underestimate the number of modulated RR Lyrae stars.

Extended CCD monitoring of RR Gem revealed that it is a Blazhko type RRab star with the shortest Blazhko period (7.23d) and smallest modulation amplitude (Delta Mmax<0.1 mag) currently known. The short period of the modulation cycle enabled us to obtain complete phase coverage of the pulsation at each phase of the modulation. This is the first multicolour observation of a Blazhko star which is extended enough to define accurate mean magnitudes and colours of the variable at different Blazhko phases. Small, but real, changes in the intensity mean colours at different Blazhko phases have been detected. The Fourier analysis of the light curves shows that, in spite of the mmag and smaller order of the amplitudes, the triplet structure is noticeable up to about the 14th harmonic. The modulation is concentrated to a very narrow, 0.2 phase range of the pulsation, centred on the supposed onset of the H emission during rising light. These observational results raise further complications for theoretical explanation of the long known but poorly understood Blazhko phenomenon.