Bias and consistency in time delay estimation methods: case of the double quasar HE 1104-1805

We present a short re-evaluation of a recently published time delay estimate for the gravitational lens system HE 1104-1805 with emphasis on important methodological aspects: bias of the statistics, inconsistency of the methods and use of the purposeful selection of data points(or so-called "cleaning") at the preprocessing stage. We show how the inadequate use of simple analysis methods can lead to too strong conclusions. Our analysis shows that there are indications for the time delay in HE 1104-1805 to be between -0.9 and -0.7 years, but still with a large uncertainty.Comment: 5 pages, 5 figures, accepted as a Letter to the Editor in A&

Three photometric methods tested on ground-based data of Q 2237+0305

The Einstein Cross, Q~2237+0305, has been photometrically observed in four bands on two successive nights at NOT (La Palma, Spain) in October 1995. Three independent algorithms have been used to analyse the data: an automatic image decomposition technique, a CLEAN algorithm and the new MCS deconvolution code. The photometric and astrometric results obtained with the three methods are presented. No photometric variations were found in the four quasar images. Comparison of the photometry from the three techniques shows that both systematic and random errors affect each method. When the seeing is worse than 1.0", the errors from the automatic image decomposition technique and the Clean algorithm tend to be large (0.04-0.1 magnitudes) while the deconvolution code still gives accurate results (1{sigma} error below 0.04) even for frames with seeing as bad as 1.7". Reddening is observed in the quasar images and is found to be compatible with either extinction from the lensing galaxy or colour dependent microlensing. The photometric accuracy depends on the light distribution used to model the lensing galaxy. In particular, using a numerical galaxy model, as done with the MCS algorithm, makes the method less seeing dependent. Another advantage of using a numerical model is that eventual non-homogeneous structures in the galaxy can be modeled. Finally, we propose an observational strategy for a future photometric monitoring of the Einstein Cross.Comment: 9 pages, accepted for publication in A&

A novel approach for extracting time-delays from lightcurves of lensed quasar images

We present a new method to estimate time delays from light curves of lensed quasars. The method is based on chi^2 minimization between the data and a numerical model light curve. A linear variation can be included in order to correct for slow long-term microlensing effects in one of the lensed images. An iterative version of the method can be applied in order to correct for higher order microlensing effects. The method is tested on simulated light curves. When higher order microlensing effects are present the time delay is best constrained with the iterative method. Analysis of a published data set for the lensed double Q0957+561 yields results in agreement with other published estimates.Comment: 6 pages, accepted for publication in A&

Microlensing in the double quasar SBS1520+530

We present the results of a monitoring campaign of the double quasar SBS1520+530 at Maidanak observatory from April 2003 to August 2004. We obtained light curves in V and R filters that show small-amplitude \Delta m~0.1 mag intrinsic variations of the quasar on time scales of about 100 days. The data set is consistent with the previously determined time delay of \Delta t=(130+-3) days by Burud et al. (2002). We find that the time delay corrected magnitude difference between the quasar images is now larger by (0.14+-0.03) mag than during the observations by Burud et al. (2002). This confirms the presence of gravitational microlensing variations in this system.Comment: 6 pages, 7 figures. Accepted for publication in A&

A method for spatial deconvolution of spectra

A method for spatial deconvolution of spectra is presented. It follows the same fundamental principles as the ``MCS image deconvolution algorithm'' (Magain, Courbin, Sohy, 1998) and uses information contained in the spectrum of a reference Point Spread Function (PSF) to spatially deconvolve spectra of very blended sources. An improved resolution rather than an infinite one is aimed at, overcoming the well known problem of ``deconvolution artefacts''. As in the MCS algorithm, the data are decomposed into a sum of analytical point sources and a numerically deconvolved background, so that the spectrum of extended sources in the immediate vicinity of bright point sources may be accurately extracted and sharpened. The algorithm has been tested on simulated data including seeing variation as a function of wavelength and atmospheric refraction. It is shown that the spectra of severely blended point sources can be resolved while fully preserving the spectrophotometric properties of the data. Extended objects ``hidden'' by bright point sources (up to 4-5 magnitudes brighter) can be accurately recovered as well, provided the data have a sufficiently high total signal-to-noise ratio (200-300 per spectral resolution element). Such spectra are relatively easy to obtain, even down to faint magnitudes, within a few hours of integration time with 10m class telescopes.Comment: 18 pages, 6 postscript figures, in press in Ap

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses VII. Time delays and the Hubble constant from WFI J2033-4723

Gravitationally lensed quasars can be used to map the mass distribution in lensing galaxies and to estimate the Hubble constant H0 by measuring the time delays between the quasar images. Here we report the measurement of two independent time delays in the quadruply imaged quasar WFI J2033-4723 (z = 1.66). Our data consist of R-band images obtained with the Swiss 1.2 m EULER telescope located at La Silla and with the 1.3 m SMARTS telescope located at Cerro Tololo. The light curves have 218 independent epochs spanning 3 full years of monitoring between March 2004 and May 2007, with a mean temporal sampling of one observation every 4th day. We measure the time delays using three different techniques, and we obtain Dt(B-A) = 35.5 +- 1.4 days (3.8%) and Dt(B-C) = 62.6 +4.1/-2.3 days (+6.5%/-3.7%), where A is a composite of the close, merging image pair. After correcting for the time delays, we find R-band flux ratios of F_A/F_B = 2.88 +- 0.04, F_A/F_C = 3.38 +- 0.06, and F_A1/F_A2 = 1.37 +- 0.05 with no evidence for microlensing variability over a time scale of three years. However, these flux ratios do not agree with those measured in the quasar emission lines, suggesting that longer term microlensing is present. Our estimate of H0 agrees with the concordance value: non-parametric modeling of the lensing galaxy predicts H0 = 67 +13/-10 km s-1 Mpc-1, while the Single Isothermal Sphere model yields H0 = 63 +7/-3 km s-1 Mpc-1 (68% confidence level). More complex lens models using a composite de Vaucouleurs plus NFW galaxy mass profile show twisting of the mass isocontours in the lensing galaxy, as do the non-parametric models. As all models also require a significant external shear, this suggests that the lens is a member of the group of galaxies seen in field of view of WFI J2033-4723.Comment: 14 pages, 12 figures, published in A&

Chandra X-ray Observations of the Quadruply Lensed Quasar RX J0911.4+0551

We present results from X-ray observations of the quadruply lensed quasar RX J0911.4+0551 using data obtained with the Advanced CCD Imaging Spectrometer (ACIS) on board the Chandra X-ray Observatory. The 29 ks observation detects a total of ~404 X-ray photons (0.3 to 7.0 keV) from the four images of the lensed quasar. Deconvolution of the aspect corrected data resolves all four lensed images, with relative positions in good agreement with optical measurements. When compared to contemporaneous optical data, one of the lensed images (component A3) is dimmer by a factor of ~6 in X-rays with respect to the 2 brighter images (components A1 and A2). Spectral fitting for the combined images shows significant intrinsic absorption in the soft (0.2 to 2.4 keV) energy band, consistent with the mini-BAL nature of this quasar, while a comparison with ROSAT PSPC observations from 1990 shows a drop of ~6.5 in the total soft bandpass flux. The observations also detect ~157 X-ray photons arising from extended emission of the nearby cluster (peaked ~42" SW of RXJ0911.4+0551) responsible for the large external shear present in the system. The Chandra observation reveals the cluster emission to be complex and non-spherical, and yields a cluster temperature of kT = 2.3^{+1.8}_{-0.8} keV and a 2.0 to 10 keV cluster luminosity within a 1 Mpc radius of L_X = 7.6_{-0.2}^{+0.6} x 10^{43} ergs/s (error bars denote 90% confidence limits). Our mass estimate of the cluster within its virial radius is 2.3^{+1.8}_{-0.7} x 10^{14} solar, and is a factor of 2 smaller than, although consistent with, previous mass estimates based on the observed cluster velocity dispersion.Comment: 16 pages, 3 figures (figure 1 is color ps). Accepted by Ap

Three photometric methods tested on ground-based data of Q 2237+0305

The Einstein Cross, Q 2237+0305, has been photometrically observed in four bands on two successive nights at NOT (La Palma, Spain) in October 1995. Three independent algorithms have been used to analyse the data: an automatic image decomposition technique, a CLEAN algorithm and the new MCS deconvolution code. The photometric and astrometric results obtained with the three methods are presented. No photometric variations were found in the four quasar images. Comparison of the photometry from the three techniques shows that both systematic and random errors affect each method. When the seeing is worse than 1farcs0 , the errors from the automatic image decomposition technique and the Clean algorithm tend to be large (0.04-0.1 magnitudes) while the deconvolution code still gives accurate results (1sigma error below 0.04) even for frames with seeing as bad as 1farcs7 . Reddening is observed in the quasar images and is found to be compatible with either extinction from the lensing galaxy or colour dependent microlensing. The photometric accuracy depends on the light distribution used to model the lensing galaxy. In particular, using a numerical galaxy model, as done with the MCS algorithm, makes the method less seeing dependent. Another advantage of using a numerical model is that eventual non-homogeneous structures in the galaxy can be modeled. Finally, we propose an observational strategy for a future photometric monitoring of the Einstein Cross. Based on observations obtained at NOT, La Palma

Bayesian approach to gravitational lens model selection: constraining H_0 with a selected sample of strong lenses

Bayesian model selection methods provide a self-consistent probabilistic framework to test the validity of competing scenarios given a set of data. We present a case study application to strong gravitational lens parametric models. Our goal is to select a homogeneous lens subsample suitable for cosmological parameter inference. To this end we apply a Bayes factor analysis to a synthetic catalog of 500 lenses with power-law potential and external shear. For simplicity we focus on double-image lenses (the largest fraction of lens in the simulated sample) and select a subsample for which astrometry and time-delays provide strong evidence for a simple power-law model description. Through a likelihood analysis we recover the input value of the Hubble constant to within 3\sigma statistical uncertainty. We apply this methodology to a sample of double image lensed quasars. In the case of B1600+434, SBS 1520+530 and SDSS J1650+4251 the Bayes' factor analysis favors a simple power-law model description with high statistical significance. Assuming a flat \LambdaCDM cosmology, the combined likelihood data analysis of such systems gives the Hubble constant H_0=76+15-5 km/s/Mpc having marginalized over the lens model parameters, the cosmic matter density and consistently propagated the observational errors on the angular position of the images. The next generation of cosmic structure surveys will provide larger lens datasets and the method described here can be particularly useful to select homogeneous lens subsamples adapted to perform unbiased cosmological parameter inferenceComment: 13 pages; 13 figures; includes Bayesian analysis of a synthetic lens catalog generated with GRAVLENS, several additional results; matches MNRAS accepted versio

Extinction Curves of Lensing Galaxies out to z=1

We present a survey of the extinction properties of ten lensing galaxies, in the redshift range z = 0.04 - 1.01, using multiply lensed quasars imaged with the ESO VLT in the optical and near infrared. The multiple images act as 'standard light sources' shining through different parts of the lensing galaxy, allowing for extinction studies by comparison of pairs of images. We explore the effects of systematics in the extinction curve analysis, including extinction along both lines of sight and microlensing, using theoretical analysis and simulations. In the sample, we see variation in both the amount and type of extinction. Of the ten systems, seven are consistent with extinction along at least one line of sight. The mean differential extinction for the most extinguished image pair for each lens is A(V) = 0.56 +- 0.04, using Galactic extinction law parametrization. The corresponding mean R_V = 2.8 +- 0.4 is consistent with that of the Milky Way at R_V = 3.1, where R_V = A(V)/E(B-V). We do not see any strong evidence for evolution of extinction properties with redshift. Of the ten systems, B1152+199 shows the strongest extinction signal of A(V) = 2.43 +- 0.09 and is consistent with a Galactic extinction law with R_V = 2.1 +- 0.1. Given the similar redshift distribution of SN Ia hosts and lensing galaxies, a large space based study of multiply imaged quasars would be a useful complement to future dark energy SN Ia surveys, providing independent constraints on the statistical extinction properties of galaxies up to z~1