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.

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.

A systematic search for multi-mode Cepheids using the database of the ASAS-SN survey has led to the detection of thirteen new double-mode and two triple-mode Cepheids in the Galactic disk. These discoveries have increased the number of Galactic disk multi-mode Cepheids by 33%. One of the new triple-mode variables pulsates simultaneously in the fundamental and in the first and the second radial overtone modes and the other in the first three radial overtone modes. Overtone triple-mode Cepheids were identified only in the Galactic bulge and in the Large and Small Magellanic Clouds previously.

The ratio of the first overtone (1O) / fundamental (F) periods of mixed-mode Cepheids that pulsate simultaneously in these two modes (F/1O) is metallicity-dependent. It can therefore be used to characterize the systems that host such variable stars. We want to take advantage of the F/1O double-mode Cepheids listed in the Gaia DR2 catalogue to derive the metallicity gradient in the Milky Way disk. The metallicity is derived from the ratio of the first overtone and fundamental periods provided by Gaia DR2 while the Gaia DR2 parallaxes are used to determine the Galactocentric distances of the stars. From a visual inspection of the light curves, it turns out that a large fraction (77%) of the Galactic F/1O double-mode Cepheids in Gaia DR2 are spurious detections. Gaia DR2 provides 3 new bona fide F/1O Cepheids. Combining them with the currently known F/1O Cepheids and using the Gaia DR2 parallaxes for the entire sample, we can derive the metallicity gradient in the Milky Way disk. We find a slope of -0.045$\pm$0.007 dex/kpc using a bootstrap method, and of -0.040$\pm$0.002 dex/kpc using a total least squares method. These results are in good agreement with previous determinations of the [Fe/H] gradient in the disk based on canonical Cepheids. The period ratio of F/1O Cepheids allows for a reliable determination of the metallicity gradient in the Milky Way, and in turn, in other systems that would be difficult to reach via classical spectroscopic methods.

We use deep multi-epoch near-IR images of the VISTA Variables in the Via Lactea (VVV) Survey to measure proper motions (PMs) of stars in the Milky Way globular cluster FSR1716 = VVV-GC05. The color-magnitude diagram of this object, made using PM selected members, shows an extended horizontal branch, nine confirmed RR Lyrae members in the instability strip, and possibly several hotter stars extending to the blue. Based on the fundamental-mode (ab-type) RR Lyrae stars that move coherently with the cluster, we confirmed that FSR1716 is an Oosterhoff I globular cluster with a mean period Pab = 0.574 days. Intriguingly, we detect tidal extensions to both sides of this cluster in the spatial distribution of PM selected member stars. Also, one of the confirmed RRabs is located 11 arcmin in projection from the cluster center, suggesting that FSR1716 may be losing stars due to the gravitational interaction with the Galaxy. We also measure radial velocities (RVs) for five cluster red giants selected using the PMs. The combination of RVs and PMs allow us to compute for the first time the orbit of this globular cluster, using an updated Galactic potential. The orbit results to be confined within |Zmax | < 2.0 kpc, and has eccentricity 0.4 < e < 0.6, with perigalactic distance 1.5 < Rperi (kpc) < 2.3, and apogalactic distance 5.3 < Rapo (kpc) < 6.4. We conclude that, in agreement with its relatively low metallicity ([Fe/H]= -1.4 dex), this is an inner halo globular cluster plunging into the disk of the Galaxy. As such, this is a unique object to test the dynamical processes that contribute to the disruption of Galactic globular clusters.

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

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.

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.

V582 Aur is an FU Ori-type young eruptive star in outburst since $\sim$1985. The eruption is currently in a relatively constant plateau phase, with photometric and spectroscopic variability superimposed. Here we will characterize the progenitor of the outbursting object, explore its environment, and analyse the temporal evolution of the eruption. We are particularly interested in the physical origin of the two deep photometric dips, one occurred in 2012, and one is ongoing since 2016. We collected archival photographic plates, and carried out new optical, infrared, and millimeter wave photometric and spectroscopic observations between 2010 and 2017, with high sampling rate during the current minimum. Beside analysing the color changes during fading, we compiled multiepoch spectral energy distributions, and fitted them with a simple accretion disk model. Based on pre-outburst data and a millimeter continuum measurement, we suggest that the progenitor of the V582 Aur outburst is a low-mass T Tauri star with average properties. The mass of an unresolved circumstellar structure, probably a disk, is 0.04 M$_{\odot}$. The optical and near-infrared spectra demonstrate the presence of hydrogen and metallic lines, show the CO bandhead in absorption, and exhibit a variable H$\alpha$ profile. The color variations strongly indicate that both the $\sim$1 year long brightness dip in 2012, and the current minimum since 2016 are caused by increased extinction along the line of sight. According to our accretion disk models, the reddening changed from $A_V$=4.5 mag to 12.5 mag, while the accretion rate remained practically constant. Similarly to the models of the UXor phenomenon of intermediate and low-mass young stars, orbiting disk structures could be responsible for the eclipses.

The simultaneous photometric and spectroscopic observations of the RR Lyrae variables in the globular cluster, M3, published in Jurcsik et al. (2017, Paper I) made it possible to perform Baade-Wesselink (BW) analysis of a large sample of Blazhko stars for the first time. The BW distances of Blazhko stars turned out to be unreliable, as significantly different distances were obtained for the stars of the Blazhko sample and also for the same star in different modulation phases. Even the results of small modulation-amplitude Blazhko stars may be doubtful. This result warns that the application of the BW method to Blazhko stars is not trustworthy. Keeping the distance fixed for each Blazhko star in each modulation phase, a significant difference between the spectroscopic and the photometric radius ($R_{\textrm{sp}}$, $R_{\textrm{ph}}$) variations is detected. The phase and amplitude variations of $R_{\textrm{sp}}$ follow the changes of the light curve during the Blazhko cycle but the $R_{\textrm{ph}}$ curve seems to be not (or only marginally) affected by the modulation. The asynchronous behaviour of $R_{\textrm{sp}}$ and $R_{\textrm{ph}}$ supports the interpretation of the Blazhko effect as a depth-dependent phenomenon, as the spectroscopic radius variation reflects the radial displacement of the line-forming region high in the atmosphere, while the photospheric radius variation is derived from the information of the observed visual-band light emitted mostly by the lower photosphere. The stability of $R_{\textrm{ph}}$ may be interpreted as a strong argument against the non-radial-mode explanation of the Blazhko phenomenon.

We present the first simultaneous photometric and spectroscopic investigation of a large set of RR Lyrae variables in a globular cluster. The radial-velocity data presented comprise the largest sample of RVs of RR Lyrae stars ever obtained. The target is M3; $BVI_{\mathrm{C}}$ time-series of 111 and $b$ flux data of further 64 RRab stars, and RV data of 79 RR Lyrae stars are published. Blazhko modulation of the light curves of 47 percent of the RRab stars are detected. The mean value of the center-of-mass velocities of RR Lyrae stars is $-146.8$ km s$^{-1}$ with 4.52 km s$^{-1}$ standard deviation, which is in good agreement with the results obtained for the red giants of the cluster. The ${\Phi_{21}}^{\mathrm RV}$ phase difference of the RV curves of RRab stars is found to be uniformly constant both for the M3 and for Galactic field RRab stars; no period or metallicity dependence of the ${\Phi_{21}}^{\mathrm RV}$ is detected. The Baade-Wesselink distances of 26 non-Blazhko variables with the best phase-coverage radial-velocity curves are determined; the corresponding distance of the cluster, $10480\pm210$ pc, agrees with the previous literature information. A quadratic formula for the $A_{\mathrm{puls}}-A_V$ relation of RRab stars is given, which is valid for both OoI and OoII variables. We also show that the $(V-I)_0$ of RRab stars measured at light minimum is period dependent, there is at least 0.1 mag difference between the colours at minimum light of the shortest- and longest-period variables.

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.

HBC 722 (V2493 Cyg) is a young eruptive star in outburst since 2010. It is an FU Orionis-type object with an atypically low outburst luminosity. Because it was well characterized in the pre-outburst phase, HBC 722 is one of the few FUors where we can learn about the physical changes and processes associated with the eruption. We monitored the source in the BVRIJHKs bands from the ground, and at 3.6 and 4.5 $\mu$m from space with the Spitzer Space Telescope. We analyzed the light curves and the spectral energy distribution by fitting a series of steady accretion disk models at many epochs. We also analyzed the spectral properties of the source based on new optical and infrared spectra. We also mapped HBC 722 and its surroundings at millimeter wavelengths. From the light curve analysis we concluded that the first peak of the outburst in 2010 September was due to an abrupt increase of the accretion rate in the innermost part of the system. This was followed by a long term process, when the brightening was mainly due to a gradual increase of the accretion rate and the emitting area. Our new observations show that the source is currently in a constant plateau phase. We found that around the peak the continuum was bluer and the H$\alpha$ profile changed significantly between 2012 and 2013. The source was not detected in the millimeter continuum, but we discovered a flattened molecular gas structure with a diameter of 1700 au and mass of 0.3 M$_{\odot}$ centered on HBC 722. While the first brightness peak could be interpreted as a rapid fall of piled-up material from the inner disk onto the star, the later monotonic flux rise suggests the outward expansion of a hot component according to the theory of Bell & Lin (1994). Our study of HBC 722 demonstrated that accretion-related outbursts can occur in young stellar objects even with very low mass disks, in the late Class II phase.

We present an update of a spectroscopic follow-up survey at low-resolution of a large number of RR Lyrae halo overdensity candidates found in the southern sky. The substructure candidates were identified in the RR Lyrae catalog of Torrealba et al. (2015) using Catalina Real-time Transient Survey (CRTS) data. Radial velocities and mean metallicities have been estimated for target stars in almost half of the original overdensities to assess their potential membership to coherent halo features.

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.

The overtone and multi-mode RR Lyrae stars in the globular cluster M3 are studied using a 200-d long, $B,V$ and $I_{\mathrm C}$ time-series photometry obtained in 2012. 70\% of the 52 overtone variables observed show some kind of multi-periodicity (additional frequency at ${f_{0.61}}={f_{\mathrm {1O}}}/0.61$ frequency ratio, Blazhko effect, double/multi-mode pulsation, period doubling). A signal at 0.587 frequency ratio to the fundamental-mode frequency is detected in the double-mode star, V13, which may be identified as the second radial overtone mode. If this mode-identification is correct, than V13 is the first RR Lyrae star showing triple-mode pulsation of the first three radial modes. Either the Blazhko effect or the ${f_{0.61}}$ frequency (or both of these phenomena) appear in 7 double-mode stars. The $P_{\mathrm{1O}}/P_{\mathrm{F}}$ period ratio of RRd stars showing the Blazhko effect are anomalous. A displacement of the main frequency component at the fundamental-mode with the value of modulation frequency (or its half) is detected in three Blazhko RRd stars parallel with the appearance of the overtone-mode pulsation. The ${f_{0.61}}$ frequency appears in RRc stars that lie at the blue side of the double-mode region and in RRd stars, raising the suspicion that its occurrence may be connected to double-mode pulsation. The changes of the Blazhko and double-mode properties of the stars are also reviewed using the recent and archive photometric data.

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.

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.

A new, extended time-series photometry of M3 RR Lyrae stars has revealed that four of the ten double-mode stars show large-amplitude Blazhko modulation of both radial modes. The first, detailed analysis of the peculiar behavior of the unique, Blazhko RRd stars is given. While the P1/P0 period ratio is normal, and the overtone mode is dominant in the other RRd stars of the cluster, the period ratio is anomalous and the fundamental mode has a larger (or similar) mean amplitude than the overtone has in Blazhko RRd stars. The modulations of the fundamental and overtone modes are synchronized only in one of the Blazhko RRd stars. No evidence of any connection between the modulations of the modes in the other three stars is found. The Blazhko modulation accounts, at least partly, for the previously reported amplitude and period changes of these stars. Contrary to the $\sim50$$%$ Blazhko statistics of RRab and RRd stars, Blazhko modulation occurs only in 10$%$ of the overtone variables in M3. Four of the five Blazhko RRc stars are bright, evolved objects, and one has a similar period and brightness as Blazhko RRd stars have. The regions of the instability strip with high and low occurrence rate of the Blazhko modulation overlap with the regions populated by first- and second-generation stars according to theoretical and observational studies, raising up the possibility that the Blazhko modulation occurs preferentially in first-generation RR Lyrae stars.

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.

In an extended photometric campaign of RR Lyrae variables of the globular cluster M3, an aberrant light-curve, non-Blazhko RRab star, V123, was detected. Based on its brightness, colors and radial velocity curve, V123 is a bona fide member of M3. The light curve of V123 exhibits neither a bump preceding light minimum, nor a hump on the rising branch, and has a longer than normal rise time, with a convex shape. Similar shape characterizes the mean light curves of some large-modulation-amplitude Blazhko stars, but none of the regular RRab variables with similar pulsation periods. This peculiar object thus mimics Blazhko variables without showing any evidence of periodic amplitude and/or phase modulation. We cannot find any fully convincing answer to the peculiar behavior of V123, however, the phenomenon raises again the possibility that rotation and aspect angle might play a role in the explanation of the Blazhko phenomenon, and some source of inhomogeneity acts (magnetic field, chemical inhomogeneity) that deforms the radial pulsation of Blazhko stars during the modulation.

We present the analysis of the multicolour CCD observations in the 2010 season of the Blazhko RRc star TV Boo. TV Boo shows a complex Blazhko modulation dominated by two independent modulations with P_{fm1}=9.74 d and P_{fm2}=21.43 d long periods. Both modulation components appear in the frequency spectra as multiplet structures around the harmonics of the pulsation. The positive value of the asymmetry parameter (Q=+0.51) of the primary modulation suggests that it is similar in its character to the Blazhko effect of most of the modulated RRab stars. Interestingly, the secondary, lower-amplitude modulation exhibits a negative asymmetry parameter (Q=-0.22), which is an unusually low value when compared to the Blazhko-modulated RRab stars. Apart from the two modulation frequencies, the spectra also show an additional frequency f' and its linear combinations with the pulsation and the primary modulation f_{1}+f' and f_{1}-f'-f_{m1}. We conclude that the additional frequency most probably belongs to a non-radial mode.

A thorough analysis of multicolour CCD observations of two modulated RRab-type variables, XY And and UZ Vir is presented. These Blazhko stars show relatively simple light-curve modulation with the usual multiplet structures in their Fourier spectra. One additional, independent frequency with linear-combination terms of the pulsation frequency is also detected in the residual spectrum of each of the two stars. The amplitude and phase relations of the triplet components are studied in detail. Most of the epoch-independent phase differences show a slight, systematic colour dependence, however, these trends have the opposite sign in the two stars. The mean values of the global physical parameters and their changes with Blazhko phase are determined utilizing the Inverse Photometric Method (IPM). The modulation properties and the IPM results are compared for the two variables. The pulsation period of XY And is the shortest when its pulsation amplitude is the highest, while UZ Vir has the longest pulsation period at this phase of the modulation. Despite this opposite behaviour, the phase relations of their mean-physical-parameter variations are similar. These results are not in accord with the predictions of the Blazhko model of Stothers (2006, ApJ, 652, 643).

The correct amplitude and phase modulation formalism of the Blazhko modulation is given. The harmonic order dependent amplitude and phase modulation form is equivalent with the Fourier decomposition of multiplets. The amplitude and phase modulation formalism used in electronic transmission technique as introduced by Benk\H{o}, Szab\'o and Papar\'o (2011, MNRAS 417, 974) for Blazhko stars oversimplifies the amplitude and phase modulation functions thus it does not describe the light variation in full detail. The results of the different formalisms are compared and documented by fitting the light curve of a real Blazhko star, CM UMa.

The analysis of recent, extended multicolour CCD and archive photoelectric, photographic and visual observations has revealed several important properties of RZ Lyr, an RRab-type variable exhibiting large-amplitude Blazhko modulation. On the time-base of \sim110 yr, a strict anticorrelation between the pulsation and modulation period changes is established. The light curve of RZ Lyr shows a remarkable bump on the descending branch in the small-amplitude phase of the modulation, similarly to the light curves of bump Cepheids. We speculate that the stellar structure temporally suits a 4:1 resonance between the periods of the fundamental and one of the higher-order radial modes in this modulation phase. The light-curve variation of RZ Lyr can be correctly fitted with a two-modulation-component solution; the 121 d period of the main modulation is nearly but not exactly four times longer than the period of the secondary modulation component. Using the inverse photometric method, the variations in the pulsation-averaged values of the physical parameters in different phases of both modulation components are determined.

The two parts of the Konkoly Blazhko Survey (KBS I and II) are introduced. The most important preliminary findings of the second part are presented in comparison to the results of the first part. Two interesting cases of very strong modulation from the KBS II are also shown.

The period-change behaviour of 134 RR Lyrae stars in the globular cluster Messier 3 (M3) is investigated on the ~120-year time base of the photometric observations. The mean period-change rates (\beta \approx 0.01 d Myr^-1) of the subsamples of variables exhibiting the most regular behaviour are in good agreement with theoretical expectations based on Horizontal-Branch stellar evolution models. However, a large fraction of variables show period changes that contradict the evolutionary expectations. Among the 134 stars studied, the period-change behaviour of only 54 variables is regular (constant or linearly changing), slight irregularities are superimposed on the regular variations in 23 cases and the remaining 57 stars display irregular period variations. The light curve of ~50 per cent of the RRab stars is not stable, i.e., these variables exhibit Blazhko modulation. The large fraction of variables with peculiar behaviour (showing light-curve modulation and/or irregular O-C variation) indicate that, probably, variables with regular period changes incompatible with their evolutionary stages also could display some kind of instability of the pulsation light curve and/or period, but the available observations have not disclosed it yet. The temporal appearence of the Blazhko effect in some stars, and the 70-90 years long regular changes preceded or followed by irregular, rapid changes of the pulsation period in some cases support this hypothesis. [...] Abstract truncated due to the limitations of astroph. See full abstract in the paper.

A systematic survey to establish the true incidence rate of the Blazhko modulation among short-period, fundamental-mode, Galactic field RR Lyrae stars has been accomplished. The Konkoly Blazhko Survey (KBS) was initiated in 2004. Since then more than 750 nights of observation have been devoted to this project. A sample of 30 RRab stars was extensively observed, and light-curve modulation was detected in 14 cases. The 47% occurrence rate of the modulation is much larger than any previous estimate. The significant increase of the detected incidence rate is mostly due to the discovery of small-amplitude modulation. Half of the Blazhko variables in our sample show modulation with so small amplitude that definitely have been missed in the previous surveys. We have found that the modulation can be very unstable in some cases, e.g. RY Com showed regular modulation only during one part of the observations while during two seasons it had stable light curve with abrupt, small changes in the pulsation amplitude. This type of light-curve variability is also hard to detect in other Survey's data. The larger frequency of the light-curve modulation of RRab stars makes it even more important to find the still lacking explanation of the Blazhko phenomenon. The validity of the [Fe/H](P,phi_{31}) relation using the mean light curves of Blazhko variables is checked in our sample. We have found that the formula gives accurate result for small-modulation-amplitude Blazhko stars, and this is also the case for large-modulation-amplitude stars if the light curve has complete phase coverage. However, if the data of large-modulation-amplitude Blazhko stars are not extended enough (e.g. < 500 data points from < 15 nights), the formula may give false result due to the distorted shape of the mean light curve used.

DM Cyg, a fundamental mode RRab star was observed in the 2007 and 2008 seasons in the frame of the Konkoly Blazhko Survey. Very small amplitude light curve modulation was detected with 10.57 d modulation period. The maximum brightness and phase variations do not exceed 0.07 mag and 7 min, respectively. In spite of the very small amplitude of the modulation, beside the frequency triplets characterizing the Fourier spectrum of the light curve two quintuplet components were also identified. The accuracy and the good phase coverage of our observations made it possible to analyse the light curves at different phases of the modulation separately. Utilizing the IP method (S\'odor, Jurcsik and Szeidl, 2009) we could detect very small systematic changes in the global mean physical parameters of DM Cyg during its Blazhko cycle. The detected changes are similar to what we have already found for a large modulation amplitude Blazhko variable MW Lyrae. The amplitudes of the detected changes in the physical parameters of DM Cyg are only about 10% of that what have been found in MW Lyr. This is in accordance with its small modulation amplitude being about one tenth of the modulation amplitude of MW Lyr.