110 research outputs found
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
- …