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 report an increase of more than 0.2 mag in the optical brightness of the leading image (A) of the gravitational lens Q0957+561, detected during the 09/2000 -- 06/2001 monitoring campaign (2001 observing season). The brightening is similar to or even greater than the largest change ever detected during the 20 years of monitoring of this system. We discuss two different provisional explanations to this event: intrinsic source variability or microlensing (either short timescale microlensing or cessation of the historical microlensing). An exhaustive photometric monitoring of Q0957+561 is needed until summer of 2002 and during 2003 to discriminate between these possibilities.

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.

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.

Photometric optical data of QSO 0957+561 covering the period 1984-99 are analyzed to discern between the two values of the time delay (417 and 424 days) mostly accepted in the recent literature. The observations, performed by groups from three different institutions-Princeton University, Harvard-Smithsonian Center for Astrophysics, and Instituto de Astrofisica de Canarias-and including new unpublished 1998-9 data from the IAC80 Telescope, were obtained in five filters (V, R, I, g, and r). The different light curves have been divided into observational seasons and two restriction have been applied to better calculate the time delay: (i) points with a strange photometric behavior have been removed; and (ii) data sets without large gaps have been selected. Simulated data were generated to test several numerical methods intended to compute the time delay. The methods giving the best results-the discrete correlation function, delta-square, z-transformed discrete correlation function, and linear interpolation-were then applied to real data. A first analysis of the 23 different time delays derived from each technique shows that the time delay must be into the interval 420-424 days. From our statistical study, a most probable value of 422.6 +/- 0.6 days is inferred.

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