Quasar microlensing offers a unique opportunity to resolve tiny sources in distant active galactic nuclei and study compact object populations in lensing galaxies. We therefore searched for microlensing-induced variability of the gravitationally lensed quasar QSO 2237+0305 (Einstein Cross) using 4374 optical frames taken with the 2.0 m Liverpool Telescope and the 1.5 m Maidanak Telescope. These $gVrRI$ frames over the 2006$-$2019 period were homogeneously processed to generate accurate long-term multi-band light curves of the four quasar images A-D. Through difference light curves, we found strong microlensing signatures. We then focused on the analytical modelling of two putative caustic-crossing events in image C, finding compelling evidence that this image experienced a double caustic crossing. Additionally, our overall results indicate that a standard accretion disc accounts reasonably well for the brightness profile of UV continuum emission sources and for the growth in source radius when the emission wavelength increases: $R_{\lambda} \propto \lambda^{\alpha}$, $\alpha$ = 1.33 $\pm$ 0.09. However, we caution that numerical microlensing simulations are required before firm conclusions can be reached on the UV emission scenario because the $VRI$-band monitoring during the first caustic crossing and one of our two $\alpha$ indicators lead to a few good solutions with $\alpha \approx$ 1.

We report the results of an optical spectroscopic follow-up of four double quasar candidates in the Sloan Digital Sky Survey (SDSS) database. SDSS J1617+3827 is most likely a lensed quasar at z = 2.079, consisting of two images with r ~ 19-21 and separated by ~ 2 arcsec. We identify an extended source northeast of the brightest image as an early-type lensing galaxy at z = 0.602, and detect a candidate for the main deflector in the vicinity of the faintest image. SDSS J2153+2732 consists of two distinct but physically associated quasars at z ~ 2.24 (r ~ 19-20, separation of 3.6 arcsec). Although this system might be a binary quasar, there is evidence of a collision or merger within a galaxy cluster at an early stage. The other two candidates are projected pairs of active galactic nuclei: SDSS J1642+3200 (separation of ~ 3 arcsec) comprises a distant quasar (r ~ 18) at z = 2.263 and the active nucleus (r > 20) of a galaxy at z ~ 0.3, while SDSS J0240-0208 (r ~ 18-19, separation of ~ 1 arcsec) is a pair of quasars at z = 1.687 and z = 1.059. In each of these two systems, the background quasar only suffers a weak gravitational lensing effect by the host galaxy of the foreground active nucleus, so the host galaxy mass is constrained to be less than (2.9-3) x 10^{11} solar masses inside 10 kpc.

Optically bright, wide separation double (gravitationally lensed) quasars can be easily monitored, leading to light curves of great importance in determining the Hubble constant and other cosmological parameters, as well as the structure of active nuclei and halos of galaxies. Searching for new double quasars in the SDSS-III database, we discovered SDSS J1442+4055. This consists of two bright images (18-19 magnitudes in the r band) of the same distant quasar at redshift z = 2.575. The two quasar images are separated by about 2.1 arcsec, show significant parallel flux variations and can be monitored from late 2015. We also found other two double quasar candidates, SDSS J1617+3827 (z = 2.079) and SDSS J1642+3200 (z = 2.264), displaying evidence for the presence of a lensing object and parallel flux variations, but requiring further spectroscopic observations to be confirmed as lensed quasars.

We report the results of our multicolor observations of PG 1115+080 with the 1.5-m telescope of the Maidanak Observatory (Uzbekistan, Central Asia) in 2001-2006. Monitoring data in filter R spanning the 2004, 2005 and 2006 seasons (76 data points) demonstrate distinct brightness variations of the source quasar with the total amplitude of almost 0.4 mag. Our R light curves have shown image C leading B by 16.4d and image (A1+A2) by 12d that is inconsistent with the previous estimates obtained by Schechter et al. in 1997 - 24.7d between B and C and 9.4d between (A1+A2) and C. The new values of time delays in PG 1115+080 must result in larger values for the Hubble constant, thus reducing difference between its estimates taken from the gravitational lenses and with other methods. Also, we analyzed variability of the A2/A1 flux ratio, as well as color changes in the archetypal "fold" lens PG 1115+080. We found the A1/A2 flux ratio to grow during 2001-2006 and to be larger at longer wavelengths. In particular, the A2/A1 flux ratio reached 0.85 in filter I in 2006. We also present evidence that both the A1 and A2 images might have undergone microlensing during 2001-2006, with the descending phase for A1 and initial phase for A2. We find that the A2/A1 flux ratio anomaly in PG 1115 can be well explained both by microlensing and by finite distance of the source quasar from the caustic fold.

We present and analyse new R-band frames of the gravitationally lensed double quasar FBQ 0951+2635. These images were obtained with the 1.5m AZT-22 Telescope at Maidanak (Uzbekistan) in the 2001-2006 period. Previous results in the R band (1999-2001 period) and the new data allow us to discuss the dominant kind of microlensing variability in FBQ 0951+2635. The time evolution of the flux ratio A/B does not favour the continuous production of short-timescale (months) flares in the faintest quasar component B (crossing the central region of the lensing galaxy). Instead of a rapid variability scenario, the observations are consistent with the existence of a long-timescale fluctuation. The flux ratio shows a bump in the 2003-2004 period and a quasi-flat trend in more recent epochs. Apart from the global behaviour of A/B, we study the intra-year variability over the first semester of 2004, which is reasonably well sampled. Short-timescale microlensing is not detected in that period. Additional data in the i band (from new i-band images taken in 2007 with the 2m Liverpool Robotic Telescope at La Palma, Canary Islands) also indicate the absence of short-timescale events in 2007.

To go into the details about the variability of the double quasar SBS 0909+532, we designed a monitoring programme with the 2 m Liverpool Robotic Telescope in the r Sloan filter, spanning 1.5 years from 2005 January to 2006 June. The r-band light curves of the A and B components, several cross-correlation techniques and a large number of simulations (synthetic light curves) lead to a robust delay of 49 +/- 6 days (1-sigma interval) that agrees with our previous results (the B component is leading). Once the time delay and the magnitude offset are known, the magnitude- and time-shifted light curve of image A is subtracted from the light curve of image B. This difference light curve of SBS 0909+532 is consistent with zero, so any possible extrinsic signal must be very weak, i.e., the observed variability in A and B is basically due to observational noise and intrinsic signal. We then make the combined light curve and analyse its statistical properties (structure functions). The structure function of the intrinsic luminosity is fitted to predictions of simple models of two physical scenarios: accretion disc instabilities and nuclear starbursts. Although no simple model is able to accurately reproduce the observed trend, symmetric triangular flares in an accretion disc seems to be the best option to account for it.

The time delays between the components of a lensed quasar are basic tools to analyze the expansion of the Universe and the structure of the main lens galaxy halo. In this paper, we focus on the variability and time delay of the double system SBS 0909+532A,B as well as the time behaviour of the field stars. We use VR optical observations of SBS 0909+532A,B and the field stars in 2003. The frames were taken at Calar Alto, Maidanak and Wise observatories, and the VR light curves of the field stars and quasar components are derived from aperture and point-spread function fitting methods. We measure the R-band time delay of the system from the chi-square and dispersion techniques and 1000 synthetic light curves based on the observed records. One nearby field star (SBS 0909+532c) is found to be variable, and the other two nearby field stars are non-variable sources. With respect to the quasar components, the R-band records seem more reliable and are more densely populated than the V-band ones. The observed R-band fluctuations permit a pre-conditioned measurement of the time delay. From the chi-square minimization, if we assume that the quasar emission is observed first in B and afterwards in A (in agreement with basic observations of the system and the corresponding predictions), we obtain a delay of - 45 (+ 1)/(- 11) days (95% confidence interval). The dispersion technique leads to a similar delay range. A by-product of the analysis is the determination of a totally corrected flux ratio in the R band (corrected by the time delay and the contamination due to the galaxy light). Our 95% measurement of this ratio (0.575 +/- 0.014 mag) is in excellent agreement with previous results from contaminated fluxes at the same time of observation.

Photometry of the Q2237+0305gravitational lens in VRI spectral bands with the 1.5-m telescope of the high-altitude Maidanak observatory in 1995-2000 is presented. Monitoring of Q2237+0305 in July-October 2000, made at nearly daily basis, did not reveal rapid (night-to-night and intranight) variations of brightness of the components during this time period. Rather slow changes of magnitudes of the components were observed, such as 0.08 mag fading of B and C components and 0.05 mag brightening of D in R band during July 23 - October 7, 2000. By good luck three nights in 1999 were almost at the time of the strong brightness peak of image C, and approximately in the middle of the ascending slope of the image A brightness peak. The C component was the most blue one in the system in 1998 and 1999, having changed its (V-I) color from 0.56 mag to 0.12 mag since August 1997, while its brightness increased almost 1.2 mag during this time period. The A component behaved similarly between August 1998 and August 2000, having become 0.47 mag brighter in R, and at the same time, 0.15 mag bluer. A correlation between the color variations and variations of magnitudes of the components is demonstrated to be significant and reaches 0.75, with a regression line slope of 0.33. A color (V-I) vrs color (V-R) plot shows the components settled in a cluster, stretched along a line with a slope of 1.31. Both slopes are noticeably smaller than those expected if a standard galactic interstellar reddening law were responsible for the differences between the colors of images and their variations over time. We attribute the brightness and color changes to microlensing of the quasar's structure, which we conclude is more compact at shorter wavelengths, as predicted by most quasar models featuring an energizing central source.

We obtained a series of more than two hundred R-band CCD images for the crowded central (115"x77") region of the metal-poor globular cluster M15 with an angular resolution of 0.5"-0.9" in most images. Optimal image subtraction was used to identify variable stars. Brightness variations were found in 83 stars, 55 of which were identified with known cluster variables and the remaining 28 are candidates for new variables. Two of them are most likely SX Phe variables. The variability type of two more stars is incertain. The remaining stars were tentatively classified as RR Lyrae variables. A preliminary analysis of published data and our results shows that the characteristics of RR Lyrae variables in the densest part (r<35") of the cluster probably change. More specifically, the maximum of the period distribution of first- and second-overtone (RR1, RR2) pulsating stars shifts toward shorter periods; i.e., there is an increase in the fraction of stars pulsating with periods < 0.3d and a deficiency of stars with 0.35d-0.40d. The ratio of the number of these short-period RR Lyrae variables to the number of fundamental-tone (RR0) pulsating variables changes appreciably. We found and corrected the error of transforming the coordinates of variables V128-155 in M15 to the coordinate system used in the catalog of variable stars in globular clusters.

We report on an observing campaign in March 2001 to monitor the brightness of the later arriving Q0957+561 B image in order to compare with the previously published brightness observations of the (first arriving) A image. The 12 participating observatories provided 3543 image frames which we have analyzed for brightness fluctuations. From our classical methods for time delay determination, we find a 417.09 +/- 0.07 day time delay which should be free of effects due to incomplete sampling. During the campaign period, the quasar brightness was relatively constant and only small fluctuations were found; we compare the structure function for the new data with structure function estimates for the 1995--6 epoch, and show that the structure function is statistically non-stationary. We also examine the data for any evidence of correlated fluctuations at zero lag. We discuss the limits to our ability to measure the cosmological time delay if the quasar's emitting surface is time resolved, as seems likely.

An observing campaign with 10 participating observatories has undertaken to monitor the optical brightness of the Q0957 gravitationally lensed quasar for 10 consecutive nights in January 2000. The resulting A image brightness curve has significant brightness fluctuations and makes a photometric prediction for the B image light curve for a second campaign planned for 12-21 March 2001. The ultimate purpose is to determine the gravitational lens time delay to a fraction of an hour, and to seek evidence for rapid microlensing.