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.

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.

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

Broad absorption line quasars (BALQSOs) are key objects for studying the structure and emission/absorption properties of AGN. However, despite their fundamental importance, the properties of BALQSOs are still not well understood. In order to investigate the X-ray nature of these sources, as well as the correlations between X-ray and rest-frame UV properties, we compile a large sample of 88 BALQSOs observed by XMM-Newton. We performed a full X-ray spectral analysis on a sample of 39 sources with higher X-ray spectral quality, and an approximate HR analysis on the remaining sources. Using available optical spectra, we calculate the BALnicity index and investigate the dependence between this optical parameter and different X-ray properties. Using the neutral absorption model, we found that 36% of our BALQSOs have NH < 5x10^21 cm^-2, lower than the expected X-ray absorption for such objects. However, when we used a physically-motivated model for the X-ray absorption in BALQSOs, i.e. ionized absorption, \sim 90% of the objects are absorbed. The absorption properties also suggest that LoBALs may be physically different objects from HiBALs. In addition, we report on a correlation between the ionized absorption column density and BAL parameters. There is evidence (at 98% level) that the amount of X-ray absorption is correlated with the strength of high-ionization UV absorption. This correlation, not previously reported, can be naturally understood in virtually all BALQSO models, as driven by the total amount of gas mass flowing towards the observer.

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.

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.

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.

We present the results of an AAT wide-field camera survey of the stars in the Monoceros Ring and purported Canis Major overdensity in the Galactic longitudes of {\it l} = (193 - 276)$^\circ$. Current numerical simulations suggest that both of these structures are the result of a single on-going accretion event, although an alternative solution is that the warped and flared disc of the Galaxy can explain the origin of both of these structures. Our results show that, with regards the Monoceros Ring, the warped and flared disc is unable to reproduce the locations and strengths of the detections observed around the Galaxy. This supports a non-Galactic origin for this structure. We report 8 new detections and 2 tentative detections of the Monoceros Ring in this survey. The exact nature of the Canis Major overdensity is still unresolved although this survey provides evidence that invoking the Galactic warp is not a sufficient solution when compared with observation. Several fields in this survey are highly inconsistent with the current Galactic disc models that include a warp and flare, to such an extent that explaining their origins with these structures is problematic. (Abridged)

The influence of microlensing in the profiles of the emission lines generated in a biconical geometry is discussed. Microlensing amplification in this anisotropic model is not directly related to the bicone's intrinsic size but depends on the orientation of the bicone axis and on the cone aperture. The orientation of the projected bicone with respect to the shear of the magnification pattern can induce very interesting effects, like the quasi-periodic enhancements of the red/blue part of the emission line profile or the lack of correlation between the broad line region (BLR) and continuum light curves of QSOs. The emission line profiles of a BLR moving in a high caustic concentration exhibit sharp features that are well defined in wavelength. These features (spikes) correspond to the scanning of the kinematics of the BLR by the caustic clusters. The biconical model can qualitatively reproduce with a transversal (with respect to the shear) movement of the BLR, the recurrent blue-wing enhancement detected in the emission line profile of the A image of the quasar lensed system SDSS J1004+4112. The probability of observing this repetitive event is almost a 2% for a fraction of matter in stars of a 5%. This result would make plausible the detection of the spectral variability in SDSS J1004+4112 under the hypothesis of microlensing of a bicone.

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 magnification induced by gravitational microlensing is sensitive to the size of a source relative to the Einstein radius, the natural microlensing scale length. This paper investigates the effect of source size in the case where the microlensing masses are distributed with a bimodal mass function, with solar mass stars representing the normal stellar masses, and smaller masses (down to $8.5\times 10^{-5}$M$_\odot$) representing a dark matter component. It is found that there exists a critical regime where the dark matter is initially seen as individual compact masses, but with an increasing source size the compact dark matter acts as a smooth mass component. This study reveals that interpretation of microlensing light curves, especially claims of small mass dark matter lenses embedded in an overall stellar population, must consider the important influence of the size of the source.

We qualitatively study the effects of gravitational microlensing on our view of unresolved extragalactic star formation regions. Using a general gravitational microlensing configuration, we perform a number of simulations that reveal that specific imprints of the star forming region are imprinted, both photometrically and spectroscopically, upon observations. Such observations have the potential to reveal the nature and size of these star forming regions, through the degree of variability observed in a monitoring campaign, and hence resolve the star formation regions in distant galaxies which are too small to be probed via more standard techniques.

The multiply imaged, gravitationally lensed quasar, QSO 2237+0305, has been the subject of recent optical monitoring campaigns, with its light curves displaying uncorrelated variability attributed to gravitational microlensing by masses in the foreground galaxy. Based on these light curves, it has been claimed that the dominant microlensing population must be a population of free-floating Jupiter-like objects; such a conclusion is not new, with several similar claims in the literature. Historically, however, it has been shown that such conclusions are flawed, with an incorrect interpretation of the complex caustic network that arises at significant optical depth. This paper examines this more recent claim, illustrating that it too is flawed.

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 present CHANDRA X-ray and NTT optical observations of the distant z=0.52 galaxy cluster RBS380 -- the most distant cluster of the ROSAT Bright Source (RBS) catalogue. We find diffuse, non-spherically symmetric X-ray emission with a X-ray luminosity of L_X(0.3-10 keV)=1.6 10^(44) erg/s, which is lower than expected from the RBS. The reason is a bright AGN in the centre of the cluster contributing considerably to the X-ray flux. This AGN could not be resolved with ROSAT. In optical wavelength we identify several galaxies belonging to the cluster. The galaxy density is at least 2 times higher than expected for such a X-ray faint cluster, which is another confirmation of the weak correlation between X-ray luminosity and optical richness. The example of the source confusion in this cluster shows how important high-resolution X-ray imaging is for cosmological research.

We present X-ray and optical observations of the z=0.52 galaxy cluster RBS380. This is the most distant cluster in the ROSAT Bright Source catalog. The cluster was observed with the CHANDRA satellite in September 2000. The optical observations were carried out with the NTT-SUSI2 camara in filters V and R in August and September 2001. The preliminary conclusions are that we see a very rich optical galaxy cluster but with a relative low X-ray luminosity. We also compare our results to other clusters with similar properties.

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.

We present a new determination of the time delay of the gravitational lens system HE1104-1805 ('Double Hamburger') based on a previously unpublished dataset. We argue that the previously published value of dt_(A-B)=0.73 years was affected by a bias of the employed method. We determine a new value of dt_(A-B)=0.85+/-0.05 years (2 sigma confidence level), using six different techniques based on non interpolation methods in the time domain. The result demonstrates that even in the case of poorly sampled lightcurves, useful information can be obtained with regard to the time delay. The error estimates were calculated through Monte Carlo simulations. With two already existing models for the lens and using its recently determined redshift, we infer a range of values of the Hubble parameter: Ho=48+/-4 km/s Mpc^-1 (2 sigma) for a singular isothermal ellipsoid (SIE) and Ho=62+/-4 km/s Mpc^-1 (2 sigma) for a constant mass-to-light ratio plus shear model (M/L+gamma). The possibly much larger errors due to systematic uncertainties in modeling the lens potential are not included in this error estimate.

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