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.

Because follow-up observations of quadruple gravitational lens systems are of extraordinary importance for astrophysics and cosmology, we present single-epoch optical spectra and $r$--band light curves of PS J0147+4630. This recently discovered system mainly consists of four images ABCD of a background quasar around a foreground galaxy G that acts as a gravitational lens. First, we use long--slit spectroscopic data in the Gemini Observatory Archive and a multi-component fitting to accurately resolve the spectra of A, D, and G. The spectral profile of G resembles that of an early-type galaxy at a redshift of 0.678 $\pm$ 0.001, which is about 20\% higher than the previous estimate. Additionally, the stellar velocity dispersion is measured to $\sim$5\% precision. Second, our early $r$-band monitoring with the Liverpool Telescope leads to accurate light curves of the four quasar images. Adopting time delays predicted by the lens model, the new lens redshift, and a standard cosmology, we report the detection of microlensing variations in C and D as large as $\sim$0.1 mag on timescales of a few hundred days. We also estimate an actual delay between A and B of a few days (B is leading), which demonstrates the big potential of optical monitoring campaigns of PS J0147+4630.

Understanding the properties of the continuum radiation and broad emission lines of active galactic nuclei provide significant information not only to model the radiation mechanism and constrain the geometry and kinematics of the broad-line region (BLR) but also to probe the central engine of the sources. Here we investigate the multifractal behaviour of the H$\beta$ emission line and the 5100 {\AA} continuum flux light curves of NGC 5548. The aim is to search for multi-scaling signatures in the light curves and check if there is a possible nonlinear relationship between them. To this end, we use a multifractality analysis technique called Multifractal Detrended Moving Average (MFDMA) analysis. We detect multifractal (nonlinear) signatures in the full monitoring and densely sampled period of the H$\beta$ line and 5100 {\AA} continuum light curves of NGC 5548, possibly indicating the presence of complex and nonlinear interaction in the 5100 {\AA} continuum and H$\beta$ emission line regions. Moreover, the degree of multifractality of H$\beta$ line is found to be about twice that of the 5100 {\AA} continuum. The nonlinearity of both emissions could be generated when the broad-line region reprocesses the radiation from the central compact source. Finally, we found that anti-persistent long-range temporal correlation is the main source of the multifractality detected in both light curves.

We present an $r$--band photometric monitoring of the two images A and B of the gravitationally lensed quasar SDSS J1442+4055 using the Liverpool Telescope (LT). From the LT light curves between 2015 December and 2018 August, we derive at once a time delay of 25.0 $\pm$ 1.5 days (1$\sigma$ confidence interval; A is leading) and microlensing magnification gradients below 10$^{-4}$ mag day$^{-1}$. The delay interval is not expected to be affected by an appreciable microlensing--induced bias, so it can be used to estimate cosmological parameters. This paper also focuses on new Gran Telescopio Canarias (GTC) and LT spectroscopic observations of the lens system. We determine the redshift of two bright galaxies around the doubly imaged quasar using LT spectroscopy, while GTC data lead to low--noise individual spectra of A, B, and the main lensing galaxy G1. The G1 spectral shape is accurately matched to an early--type galaxy template at $z$ = 0.284, and it has potential for further relevant studies. Additionally, the quasar spectra show absorption by metal--rich gas at $z \sim$ 2. This dusty absorber is responsible for an extinction bump at a rest--frame wavelength of 2209 $\pm$ 2 \AA, which has strengths of $\sim$ 0.47 and 0.76 mag $\mu$m$^{-1}$ for A and B, respectively. In such intervening system, the dust--to--gas ratio, gas--phase metallicity indicator [Zn/H], and dust depletion level [Fe/Zn] are relatively high.

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.

In the framework of the Gravitational LENses and DArk MAtter (GLENDAMA) project, we present a database of nine gravitationally lensed quasars (GLQs) that have two or four images brighter than $r$ = 20 mag and are located in the northern hemisphere. This new database consists of a rich variety of follow-up observations included in the GLENDAMA global archive, which is publicly available online and contains 6557 processed astronomical frames of the nine lens systems over the period 1999$-$2016. In addition to the GLQs, our archive also incorporates binary quasars, accretion-dominated radio-loud quasars, and other objects, where about 50% of the non-GLQs were observed as part of a campaign to identify GLQ candidates. Most observations of GLQs correspond to an ongoing long-term macro-programme with 2$-$10 m telescopes at the Roque de los Muchachos Observatory, and these data provide information on the distribution of dark matter at all scales. We outline some previous results from the database, and we additionally obtain new results for several GLQs that update the potential of the tool for astrophysical studies.

We analyse new optical observations of the gravitational lens system SDSS J1515+1511. These include a 2.6-year photometric monitoring with the Liverpool Telescope (LT) in the $r$ band, as well as a spectroscopic follow-up with the LT and the Gran Telescopio Canarias (GTC). Our $r$-band LT light curves cover a quiescent microlensing period of the doubly imaged quasar at $z_{\rm s}$ = 2.049, which permits us to robustly estimate the time delay between the two images A and B: 211 $\pm$ 5 days (1$\sigma$ confidence interval; A is leading). Unfortunately, the main lensing galaxy (G1) is so faint and close to the bright quasar that it is not feasible to accurately extract its spectrum through the GTC data. However, assuming the putative redshift $z_{\rm G1}$ = 0.742, the GTC and LT spectra of the distant quasar are used to discuss the macrolens magnification, and the extinction and microlensing effects in G1. The new constraints on the time delay and macrolens magnification ratio essentially do not change previous findings on the mass scale of G1 and external shear, while the redshift of the lensing mass is found to be consistent with the assumed value of $z_{\rm G1}$. This is a clear evidence that G1 is indeed located at $z_{\rm G1}$ = 0.742. From the GTC data we also obtain the redshift of two additional objects (the secondary galaxy G2 and a new absorption system) and discuss their possible role in the lens scenario.

We spectroscopically re-observed the gravitational lens system SDSS J1339+1310 using OSIRIS on the GTC. We also monitored the $r$-band variability of the two quasar images (A and B) with the LT over 143 epochs in the period 2009$-$2016. These new data in both the wavelength and time domains have confirmed that the system is an unusual microlensing factory. The C$\scriptsize{\rm{IV}}$ emission line is remarkably microlensed, since the microlensing magnification of B relative to that for A, $\mu_{\rm{BA}}$, reaches a value of 1.4 ($\sim$ 0.4 mag) for its core. Moreover, the B image shows a red wing enhancement of C$\scriptsize{\rm{IV}}$ flux (relative to A), and $\mu_{\rm{BA}}$ = 2 (0.75 mag) for the C$\scriptsize{\rm{IV}}$ broad-line emission. Regarding the nuclear continuum, we find a chromatic behaviour of $\mu_{\rm{BA}}$, which roughly varies from $\sim$ 5 (1.75 mag) at 7000 \AA\ to $\sim$ 6 (1.95 mag) at 4000 \AA. We also detect significant microlensing variability in the $r$ band, and this includes a number of microlensing events on timescales of 50$-$100 d. Fortunately, the presence of an intrinsic 0.7 mag dip in the light curves of A and B, permitted us to measure the time delay between both quasar images. This delay is $\Delta t_{\rm{AB}}$ = 47$^{+5}_{-6}$ d (1$\sigma$ confidence interval; A is leading), in good agreement with predictions of lens models.

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.

This is the first version (v1) of the Gravitational LENses and DArk MAtter (GLENDAMA) database accessible at http://grupos.unican.es/glendama/database The new database contains more than 6000 ready-to-use (processed) astronomical frames corresponding to 15 objects that fall into three classes: (1) lensed QSO (8 objects), (2) binary QSO (3 objects), and (3) accretion-dominated radio-loud QSO (4 objects). Data are also divided into two categories: freely available and available upon request. The second category includes observations related to our yet unpublished analyses. Although this v1 of the GLENDAMA archive incorporates an X-ray monitoring campaign for a lensed QSO in 2010, the rest of frames (imaging, polarimetry and spectroscopy) were taken with NUV, visible and NIR facilities over the period 1999-2014. The monitorings and follow-up observations of lensed QSOs are key tools for discussing the accretion flow in distant QSOs, the redshift and structure of intervening (lensing) galaxies, and the physical properties of the Universe as a whole.

We are developing a web-based interactive software to simulate crowded-field imaging with the Imaging and Slitless Spectrograph Instrument for Surveys (ISSIS) on board the future World Space Observatory - Ultraviolet (WSO-UV). This new tool is aimed to prepare WSO-UV/ISSIS proposals to observe multicomponent targets and dense fields. For a given combination of UV channel, filters and exposure time, the user creates a set of point-like and extended sources (source model). This source model produces a final image, which takes into account a pixelated field of view, a realistic conversion between physical flux and counts per second, the convolution with the expected point spread function, a sky background and noise fluctuations. The current version of the simulator is available at the Gravitational LENses and DArk MAtter (Glendama) website, and it allows users to specify all relevant parameters of each point-like or extended source, drag-and-drop sources by using a mouse or a fingertip/stylus on a touchscreen, change the frame size or the brightness scale, etc.

We check the performance of the FRODOSpec integral-field spectrograph for observations of faint sources in crowded fields. Although the standard processing pipeline L2 yields too noisy fibre spectra, we present a new processing software (L2LENS) that gives rise to accurate spectra for the two images of the gravitationally lensed quasar Q0957+561. Among other things, this L2LENS reduction tool accounts for the presence of cosmic-ray events, scattered-light backgrounds, blended sources, and chromatic source displacements due to differential atmospheric refraction. Our non-standard reduction of Q0957+561 data shows the ability of FRODOSpec to provide useful information on a wide variety of targets, and thus, the big potential of integral-field spectrographs on current and future robotic telescopes.

We present three complete seasons and two half-seasons of SDSS r-band photometry of the gravitationally lensed quasar SBS 0909+532 from the U.S. Naval Observatory, as well as two seasons each of SDSS g-band and r-band monitoring from the Liverpool Robotic Telescope. Using Monte Carlo simulations to simultaneously measure the system's time delay and model the r-band microlensing variability, we confirm and significantly refine the precision of the system's time delay to \Delta t_{AB} = 50^{+2}_{-4} days, where the stated uncertainties represent the bounds of the formal 1\sigma\ confidence interval. There may be a conflict between the time delay measurement and a lens consisting of a single galaxy. While models based on the Hubble Space Telescope astrometry and a relatively compact stellar distribution can reproduce the observed delay, the models have somewhat less dark matter than we would typically expect. We also carry out a joint analysis of the microlensing variability in the r- and g-bands to constrain the size of the quasar's continuum source at these wavelengths, obtaining log[(r_{s,r}/cm) [cos{i}/0.5]^{1/2}] = 15.3 \pm 0.3 and log[(r_{s,g}/cm) [cos{i}/0.5]^{1/2}] = 14.8 \pm 0.9, respectively. Our current results do not formally constrain the temperature profile of the accretion disk but are consistent with the expectations of standard thin disk theory.

New light curves of the gravitationally lensed double quasar Q0957+561 in the gr bands during 2008-2010 include densely sampled, sharp intrinsic fluctuations with unprecedentedly high signal-to-noise ratio. These relatively violent flux variations allow us to very accurately measure the g-band and r-band time delays between the two quasar images A and B. Using correlation functions, we obtain that the two time delays are inconsistent with each other at the 2sigma level, with the r-band delay exceeding the 417-day delay in the g band by about 3 days. We also studied the long-term evolution of the delay-corrected flux ratio B/A from our homogeneous two-band monitoring with the Liverpool Robotic Telescope between 2005 and 2010. This ratio B/A slightly increases in periods of violent activity, which seems to be correlated with the flux level in these periods. The presence of the previously reported dense cloud within the cD lensing galaxy, along the line of sight to the A image, could account for the observed time delay and flux ratio anomalies.

Q0957+561 was the first discovered gravitationally lensed quasar. The mirage shows two images of a radio-loud quasar at redshift z = 1.41. The time lag between these two images is well established around one year. We detected a very prominent variation in the optical brightness of Q0957+561A at the beginning of 2009, which allowed us to predict the presence of significant intrinsic variations in multi-wavelength light curves of Q0957+561B over the first semester of 2010. To study the predicted brightness fluctuations of Q0957+561B, we conducted an X-ray, NUV, optical and NIR monitoring campaign using both ground-based and space-based facilities. The continuum NUV-optical light curves revealed evidence of a centrally irradiated, standard accretion disk. In this paper, we focus on the radial structure of the standard accretion disk and the nature of the central irradiating source in the distant radio-loud active galactic nucleus (AGN).

We present the overview and current results of an ongoing optical/NIR monitoring of seven GLQs with the 2-m Liverpool Robotic Telescope. The photometric data over the first seven years of this programme (2005-2011) are leading to high-quality light curves, which in turn are being used as key tools for different standard and novel studies. While brightness records of non-lensed distant quasars may contain unrecognized extrinsic variations, one can disentangle intrinsic from extrinsic signal in certain GLQs. Thus, some GLQs in our sample allow us to assess their extrinsic and intrinsic variations, as well as to discuss the origin of both kinds of fluctuations. We also demonstrate the usefulness of GLQ time-domain data to obtain successful reverberation maps of inner regions of accretion disks around distant supermassive black holes, and to estimate redshifts of distant lensing galaxies.

We are using OSIRIS at the 10.4-m Gran Telescopio CANARIAS (GTC) to obtain spectra of faint galaxies close to multiply imaged quasars. We initially focused on the fields of HE 1413+117 (Cloverleaf quasar) and SDSS 1116+4118. In this contribution, we present long-slit spectroscopy of two galaxies in the southwest of the Cloverleaf, and show that one of them makes a negligible contribution to the external shear of the gravitational lens system. Spectra of the main lensing galaxy candidate in SDSS 1116+4118 are also analysed and discused. If gravitational lensing is causing the quasar image splitting, our spectra reveal that the main lens can not consist of only one dominant galaxy.

We present the results of our X-ray, UV and optical monitoring campaign of the first gravitationally lensed AGN from late 2009 to mid 2010. The trailing (B) image of the AGN 0957+561 shows the intrinsic continuum variations that were predicted in advance based on observations of the leading (A) image in the gr optical bands. This multiwavelength variability of the B image allows us to carry out a reverberation mapping analysis in the radio-loud AGN 0957+561 at redshift z = 1.41. We find that the U-band and r-band light curves are highly correlated with the g-band record, leading and trailing it by 3 +/- 1 days (U band) and 4 +/- 1 days (r band). These 1-sigma measurements are consistent with a scenario in which flares originated in the immediate vicinity of the supermassive black hole are thermally reprocessed in a standard accretion disk at about 10-20 Schwarzschild radii from the central dark object. We also report that the light curve for the X-ray emission with power-law spectrum is delayed with respect to those in the Ugr bands by about 32 days. Hence, the central driving source can not be a standard corona emitting the observed power-law X-rays. This result is also supported by X-ray reprocessing simulations and the absence of X-ray reflection features in the spectrum of 0957+561. We plausibly interpret the lack of reflection and the 32-day delay as evidence for a power-law X-ray source in the base of the jet at a typical height of about 200 Schwarzschild radii. A central EUV source would drive the variability of 0957+561.

Gravitationally lensed QSOs (GLQs) at redshift z = 1-2 play a key role in understanding the cosmic evolution of the innermost parts of active galaxies (black holes, accretion disks, coronas and internal jets), as well as the structure of galaxies at intermediate redshifts. With respect to studies of normal QSOs, GLQ programmes have several advantages. For example, a monitoring of GLQs may lead to unambiguous detections of intrinsic and extrinsic variations. Both kinds of variations can be used to discuss central engines in distant QSOs, and mass distributions and compositions of lensing galaxies. In this context, UV data are of particular interest, since they correspond to emissions from the immediate surroundings of the supermassive black hole. We describe some observation strategies to analyse optically bright GLQs at z of about 1.5, using ISSIS (CfS) on board World Space Observatory-Ultraviolet.

The quadruple quasar H1413+117 (z_s = 2.56) has been monitored with the 2.0 m Liverpool Telescope in the r Sloan band from 2008 February to July. This optical follow-up leads to accurate light curves of the four quasar images (A-D), which are defined by 33 epochs of observation and an average photometric error of \sim 15 mmag. We then use the observed (intrinsic) variations of \sim 50-100 mmag to measure the three time delays for the lens system for the first time (1\sigma confidence intervals): \Delta \tau_{AB} = -17 +/- 3, \Delta \tau_{AC} = -20 +/- 4, and \Delta \tau_{AD} = 23 +/- 4 days (\Delta \tau_{ij} = \tau_j - \tau_i; B and C are leading, while D is trailing). Although time delays for lens systems are often used to obtain the Hubble constant (H_0), the unavailability of the spectroscopic lens redshift (z_l) in the system H1413+117 prevents a determination of H_0 from the measured delays. In this paper, the new time delay constraints and a concordance expansion rate (H_0 = 70 km s^{-1} Mpc^{-1}) allow us to improve the lens model and to estimate the previously unknown z_l. Our 1\sigma estimate z_l = 1.88^{+0.09}_{-0.11} is an example of how to infer the redshift of very distant galaxies via gravitational lensing.

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.

We extend the gr-band time coverage of the gravitationally lensed double quasar Q0957+561. New gr light curves permit us to detect significant intrinsic fluctuations, to determine new time delays, and thus to gain perspective on the mechanism of intrinsic variability in Q0957+561. We use new optical frames of Q0957+561 in the g and r passbands from January 2005 to July 2007. These frames are part of an ongoing long-term monitoring with the Liverpool robotic telescope. We also introduce two photometric pipelines that are applied to the new gr frames of Q0957+561. The transformation pipeline incorporates zero-point, colour, and inhomogeneity corrections to the instrumental magnitudes, so final photometry to the 1-2% level is achieved for both quasar components. The two-colour final records are then used to measure time delays. The gr light curves of Q0957+561 show several prominent events and gradients, and some of them (in the g band) lead to a time delay between components of 417 +/- 2 d (1 sigma). We do not find evidence of extrinsic variability in the light curves of Q0957+561. We also explore the possibility of a delay between a large event in the g band and the corresponding event in the r band. The gr cross-correlation reveals a time lag of 4.0 +/- 2.0 d (1 sigma; the g-band event is leading) that confirms a previous claim of the existence of a delay between the g and r band in this lensed quasar. The time delays (between quasar components and between optical bands) from the new records and previous ones in similar bands indicate that most observed variations in Q0957+561 (amplitudes of about 100 mmag and timescales of about 100 d) are very probably due to reverberation within the gas disc around the supermassive black hole.

We present a detailed structure function analysis of the UV variability of Q0957+561. From new optical observations, we constructed normalized structure functions of the quasar luminosity at restframe wavelengths of 2100 and 2600 \AA. Old optical records also allow the structure function to be obtained at 2100 \AA, but 10 years ago in the observer's frame. These three structure functions are then compared to predictions of both simple and relatively sophisticated (incorporating two independent variable components) Poissonian models. We do not find clear evidence of a chromatic mechanism of variability. From the recent data, 100-d time-symmetric and 170-d time-asymmetric flares are produced at both restframe wavelengths. Taking into account measurements of time delays and the existence of an EUV/radio jet, reverberation is probably the main mechanism of variability. Thus, two types of EUV/X-ray fluctuations would be generated within or close to the jet and later reprocessed by the disc gas in the two emission rings. The 100-d time-symmetric shots are also responsible for most of the 2100 \AA variability detected in the old experiment. However, there is no evidence of asymmetric shots in the old UV variability. If reverberation is the involved mechanism of variability, this could mean an intermittent production of high-energy asymmetric fluctuations. The old records are also consistent with the presence of very short-lifetime (10 d) symmetric flares, which may represent additional evidence of time evolution. We also discuss the quasar structure that emerges from the variability scenario.

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 aim of this contribution is to present the two first phases of the optical monitoring programme of the Gravitational Lenses group at the Universidad de Cantabria (GLUC, http://grupos.unican.es/glendama/). In an initial stage (2003 March-June), the Estacion de Observacion de Calar Alto (EOCA) was used to obtain VR frames of SBS 0909+532 and QSO 0957+561. These observations in 2003 led to accurate fluxes of the two components of both double QSOs, which are being compared and complemented with data from other 1-1.5 m telescopes located in the North Hemisphere: Fred Lawrence Whipple Observatory (USA), Maidanak Observatory (Uzbekistan) and Wise Observatory (Israel). On the other hand, the GLUC started the second phase of its monitoring programme in 2005 January. In this second phase, they are using the 2 m fully robotic Liverpool Telescope (LT). The key idea is the two-band photometric follow-up of four lensed QSOs with different main lensing galaxies: SBS 0909+532 (elliptical), QSO 0957+561 (giant cD), B1600+434 (edge-on spiral) and QSO 2237+0305 (face-on spiral). Thus, the light rays associated with the components of the four gravitational mirages cross different galaxy environments, and the corresponding light curves could unveil the content of these environments. While SBS 0909+532 and QSO 0957+561 are the targets for the two first years with the LT (2005-2006), the rest of targets (B1600+434 and QSO 2237+0305) will be monitored starting from 2007.

We re-investigate the gravitationally lensed system Q2237+0305 data record to quantify the probability of having a caustic crossing in the A component. Several works assume that this is the case, but no quantitative analysis is available in the literature. We combine the datasets from the OGLE and GLITP collaborations to accurately trace the prominent event in the lightcurve for the A component of the system. Then the observed event is compared with synthetic light curves derived from trajectories in magnification maps. These maps are generated using a ray-tracing technique. We take more than 10^9 trajectories and test a wide range of different physical properties of the lensing galaxy and the source quasar (lens transverse velocity, microlens mass, source intensity profile and source size). We found that around 75% of our good trajectories (i.e. that are consistent with the observations) are caustic crossings. In addition, a high transverse velocity exceeding 300 km/s, a microlens mass of about 0.1 M_sun and a small standard accretion disk is the best parameter combination. The results justify the interpretation of the OGLE-GLITP event in Q2237+0305A as a caustic crossing. Moreover, the physical properties of the lens and source are in very good agreement with previous works. We also remark that a standard accretion disk is prefered to those simpler approaches, and that the former should be used in subsequent simulations.

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.

We present an analysis of infrared/optical/ultraviolet spectra of the two images of the first gravitationally lensed quasar Q0957+561A, B. The Hubble Space Telescope observations of Q0957+561A and Q0957+561B are separated in time by the known time delay in this system, so we can directly deduce the flux ratios. These flux ratios of images lead to important information on the dark halo of the main lens galaxy (a giant elliptical at redshift z = 0.36). Our measurements for the continuum are in good agreement with extinction in the elliptical galaxy and a small fraction of mass in collapsed objects (no need for gravitational microlensing). From the continuum and emission line ratios, we also show evidence in favour of the existence of a network of compact dusty clouds.

We present a method for the determination of upper limits on the transverse velocity of the lensing galaxy in the quadruple quasar system Q2237+0305, based on the lack of strong microlensing signatures in the quasar lightcurves. The limits we derive here are based on four months of high quality monitoring data, by comparing the low amplitudes of the lightcurves of the four components with extensive numerical simulations. We make use of the absence of strong variability of the components (especially components B and D) to infer that a "flat" time interval of such a length is only compatible with an effective transverse velocity of the lensing galaxy of v_bulk <=630 km/s for typical microlenses masses of M_microlens=0.1 M_solar (or v_bulk <=2160 km/s for M_microlens=1.0 M_solar) at the 90% confidence level. This method may be applicable in the future to other systems.

We present VR magnification ratios of QSO 0957+561, which are inferred from the GLITP light curves of Q0957+561A and new frames taken with the 2.56m Nordic Optical Telescope about 14 months after the GLITP monitoring. From two photometric approaches and a reasonable range for the time delay in the system (415-430 days), we do not obtain achromatic optical continuum ratios, but ratios depending on the wavelength. These new measurements are consistent with differential extinction in the lens galaxy, the Lyman limit system, the damped Ly-alpha system, or the host galaxy of the QSO. The possible values for the differential extinction and the ratio of total to selective extinction in the V band are reasonable. Moreover, crude probability arguments suggest that the ray paths of the two components cross a similar dusty environment, including a network of compact dust clouds and compact dust voids. As an alternative (in fact, the usual interpretation of the old ratios), we also try to explain the new ratios as caused by gravitational microlensing in the deflector. From magnification maps for each of the gravitationally lensed images, using different fractions of the surface mass density represented by the microlenses, as well as different sizes and profiles of the V-band and R-band sources, several synthetic distributions of V-band and R-band ratios are derived. In some gravitational scenarios, there is an apparent disagreement between the observed pair of ratios and the simulated distributions. However, several microlensing pictures work well. To decide between either extinction, or microlensing, or a mixed scenario (extinction + microlensing), new observational and interpretation efforts are required.

We introduce a new scheme to study the nature of the central engine in a lensed QSO. The compact emission regions could have different sizes in different optical wavelengths, and our framework permits to obtain the source size ratios when a microlensing special high-magnification event (e.g., a caustic crossing event, a two-dimensional maximum crossing event and so on) is produced in one of the QSO components. To infer the source size ratios, only cross-correlations between the brightness records in different optical bands are required. While the deconvolution method leads to a richer information (1D intrinsic luminosity profiles), the new approach is free of the technical problems with complex inversion procedures. Using simulations related to recent VR data of Q2237+0305A, we discuss the ability of the scheme in the determination of the visible-to-red ratio q = $R_V/R_R$. We conclude that extremely accurate fluxes (with a few microJy uncertainties, or equivalently, a few milli-magnitudes errors) can lead to ~10% measurements of q. Taking into account the errors in the fluxes of Q2237+0305A from a normal ground-based telescope, ~10 microJy (~10 mmag), it must be possible the achievement of smaller errors from the current superb-telescopes, and thus, an accurate determination of q. Obviously, to measure the visible-to-red ratio, the light curves cannot be contaminated by an intrinsic event or an important high-frequency intrinsic signal, i.e., exceeding the microJy (mmag) level. For an arbitrary lensed QSO, we finally remark that the framework seems to work better with very fast microlensing events.

The GLITP collaboration observed the first gravitational lens system (QSO 0957+561) from 2000 February 3 to 2000 March 31. The daily VR observations were made with the 2.56-m Nordic Optical Telescope at Roque de los Muchachos Observatory, La Palma (Spain). We have derived detailed and robust VR light curves of the two components Q0957+561A and Q0957+561B. In spite of the excellent sampling rate, we have not found evidence in favor of true daily variability. With respect to variability on time-scales of several weeks, we measure VR gradients of about -0.8 mmag/day in Q0957+561A and + 0.3 mmag/day in Q0957+561B. The gradients are very probably originated in the far source, thus adopting this reasonable hypothesis (intrinsic variability), we compare them to the expected gradients during the evolution of a compact supernova remnant at the redshift of the source quasar. The starburst scenario is roughly consistent with some former events, but the new gradients do not seem to be caused by supernova remnant activity.

We review some techniques to study the nature and size of the optical continuum sources in multiple QSOs. We focus on the source originating the events with several months timescale (the rapid variability source) as well as the source that is responsible for the non-variable background component (the background source). The techniques are used to study both the rapid variability source in Q0957+561 and the main (compact) background source in Q2237+0305.

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.

We study a recently observed gravitational microlensing peak in the V-band light curve of Q2237+0305A using a relatively simple, but highly consistent with the data (the best-fit reduced \chi^2 is very close to 1), physical model. The source quasar is assumed to be a Newtonian geometrically-thin and optically-thick accretion disk. The disk has an arbitrary orientation, and both blackbody and greybody emission spectra are considered. When the electron-photon scattering plays a role, the greybody spectrum will be a simplified version of the exact one. In our model the microlensing variability result from the source crossing a caustic straight line. The main goal is to estimate the black hole mass and the mass accretion rate in QSO 2237+0305 as well as to discuss the power and the weakness of the technique, some possible improvements, and the future prospects from multifrequency monitoring of new microlensing peaks. We also put into perspective the new methodology and the results on the central engine in QSO 2237+0305. From the fitted microlensing parameters and reasonable dynamical/cosmological constraints, it is concluded that QSO 2237+0305 harbours a central massive black hole: 10^7 M_Sun < M < 6 10^8 M_Sun. While the information about the central dark mass is very interesting, the mass accretion rate is not so well constrained. The typical values of the disk luminosity/Eddington luminosity ratio are in the (1-20)*\epsilon range, where \epsilon \leq 1 is the emissivity relative to a blackbody and the highest L/L_{Edd} ratio corresponds to the largest deflector motion. Therefore, in order to verify L/L_{Edd} \leq 1, a relatively small projected peculiar motion of the lens galaxy and a greybody emission seem to be favored.

From the peak of a gravitational microlensing high-magnification event in the A component of QSO 2237+0305, which was accurately monitored by the GLITP collaboration, we derived new information on the nature and size of the optical V-band and R-band sources in the far quasar. If the microlensing peak is caused by a microcaustic crossing, we firstly obtained that the standard accretion disk is a scenario more reliable/feasible than other usual axially symmetric models. Moreover, the standard scenario fits both the V-band and R-band observations with reduced chi-square values very close to one. Taking into account all these results, a standard accretion disk around a supermassive black hole is a good candidate to be the optical continuum main source in QSO 2237+0305. Secondly, using the standard source model and a robust upper limit on the transverse galactic velocity, we inferred that 90 per cent of the V-band and R-band luminosities are emitted from a region with radial size less than 1.2 10^{-2} pc (= 3.7 10^{16} cm, at 2*sigma confidence level).

We present VR observations of QSO 2237+0305 conducted by the GLITP collaboration from 1999 October 1 to 2000 February 3. The observations were made with the 2.56 m Nordic Optical Telescope at Roque de los Muchachos Observatory, La Palma (Spain). The PSF fitting method and an adapted version of the ISIS subtraction method have been used to derive the VR light curves of the four components (A-D) of the quasar. The mean errors range in the intervals 0.01-0.04 mag (PSF fitting) and 0.01-0.02 mag (ISIS subtraction), with the faintest component (D) having the largest uncertainties. We address the relatively good agreement between the A-D light curves derived using different filters, photometric techniques, and telescopes. The new VR light curves of component A extend the time coverage of a high magnification microlensing peak, which was discovered by the OGLE team.

Q0957+561A,B is a double-imaged quasar that has been intensively observed during the last 10 years in different optical bands and with several telescopes, and we concentrated on recent public data obtained at the Apache Point Observatory (APO) and the Teide Observatory (TO). When an intrinsic event appears in the light curve of Q0957+561A, its twin event (a similar feature) is seen in the brightness record of Q0957+561B, and thus, one can measure the corresponding time delay. The TO dataset includes two prominent twin events, which were detected with a time separation of 425 +/- 4 days. On the other hand, from the APO dataset, we found a clear evidence for two different time delays associated with two pairs of twin events: 417.0 +/- 0.6 (APO main twin events) and 432.0 +/- 1.9 days (APO secondary twin events), where the APO(main)-APO(secondary) difference delay is of - 15 +/- 2 days. In agreement with the Yonehara's idea, if the three pairs of twin events are originated inside a standard hybrid source (accretion disk and circumnuclear stellar region), the three measured time delays indicate that they do not come from a common zone in the source. Therefore, we can consider that the prominent features are caused by flares in a standard hybrid source and discuss on the size and nature of the region of flares. In this paper it is showed that the more plausible interpretation is that two of the three flares are generated at distances (from the central black hole) larger than 90 pc. Some stellar scenarios can explain the two flares far away from the black hole, while phenomena in the accretion disk cannot cause them.

From optical R band data of the double quasar QSO 0957+561A,B, we made two new difference light curves (about 330 days of overlap between the time-shifted light curve for the A image and the magnitude-shifted light curve for the B image). We observed noisy behaviours around the zero line and no short-timescale events (with a duration of months), where the term event refers to a prominent feature that may be due to microlensing or another source of variability. Only one event lasting two weeks and rising - 33 mmag was found. Measured constraints on the possible microlensing variability can be used to obtain information on the granularity of the dark matter in the main lensing galaxy and the size of the source. In addition, one can also test the ability of the observational noise to cause the rms averages and the local features of the difference signals. We focused on this last issue. The combined photometries were related to a process consisting of an intrinsic signal plus a Gaussian observational noise. The intrinsic signal has been assumed to be either a smooth function (polynomial) or a smooth function plus a stationary noise process or a correlated stationary process. Using these three pictures without microlensing, we derived some models totally consistent with the observations. We finally discussed the sensitivity of our telescope (at Teide Observatory) to several classes of microlensing variability.

CCD observations of the gravitational lens system Q0957+561A,B in the BVRI bands are presented in this paper. The observations, taken with the 82 cm IAC-80 telescope, at Teide Observatory, Spain, were made from the beginning of 1996 February to 1998 July, as part of an on-going lens monitoring program. Accurate photometry was obtained by simultaneously fitting a stellar two-dimensional profile on each component by means of DAOPHOT software. This alternative method equals and even improves the results obtained with previous techniques. The final dataset is characterized by its high degree of homogeneity as it was obtained using the same telescope and instrumentation during a period of almost 3 years. The resulting delay, obtained with a new method, the delta2-test, is of 425 +/- 4 days, slightly higher than the value previously accepted (417 days), but concordant with the results obtained by Pelt et al. (1996); Oscoz et al. (1997); Pijpers (1997); and Goicoechea et al. (1999).

We present new long-slit William Herschel Telescope spectroscopic observations of the lens galaxy G1 associated with the double-imaged QSO 0957+561A,B. The obtained central stellar velocity dispersion, sigma_l = 310 +/- 20 km/s, is in reasonable agreement with other measurements of this dynamical parameter. Using all updated measurements of the stellar velocity dispersion in the internal region of the galaxy (at angular separations < 1".5) and a simple isotropic model, we discuss the mass of a possible central massive dark object. It is found that the data of Falco et al. (1997) suggest the existence of an extremely massive object of (0.5-2.1) x 10E10/h M_\odot (80% confidence level), whereas the inclusion of very recent data (Tonry & Franx 1998, and this paper) substantially changes the results: the compact central mass must be $\le$ 6 x10E9/h M_\odot at the 90% confidence level. We note that, taking into account all the available dynamical data, a compact nucleus with a mass of 10E9/h M_\odot (best fit) cannot be ruled out.

The best studied gravitational lens is the Twin QSO Q0957+561, and as a result of recent advances in observations and modelling, a ten percent measurement of the Hubble constant is attainable. We discuss in detail this system and review the perspectives in a near future.