Direct detection of galaxy stellar halos : NGC 3957 as a test case

We present a direct detection of the stellar halo of the edge-on S0 galaxy NGC 3957, using ultra-deep VLT/VIMOS V and R images. This is achieved with a sky subtraction strategy based on infrared techniques. These observations allow us to reach unprecedented high signal-to-noise ratios up to 15 kpc away from the galaxy center, rendering photon-noise negligible. The 1 sigma detection limits are R = 30.6 mag/arcsec^2 and V = 31.4 mag/arcsec^2. We conduct a thorough analysis of the possible sources of systematic errors that could affect the data: flat-fielding, differences in CCD responses, scaling of the sky background, the extended halo itself, and PSF wings. We conclude that the V-R colour of the NGC 3957 halo, calculated between 5 and 8 kpc above the disc plane where the systematic errors are modest, is consistent with an old and preferentially metal-poor normal stellar population, like that revealed in nearby galaxy halos from studies of their resolved stellar content. We do not find support for the extremely red colours found in earlier studies of diffuse halo emission, which we suggest might have been due to residual systematic errors.Comment: Accepted for publication in A&A - "language edited

Firedec: a two-channel finite-resolution image deconvolution algorithm

We present a two-channel deconvolution method that decomposes images into a parametric point-source channel and a pixelized extended-source channel. Based on the central idea of the deconvolution algorithm proposed by Magain, Courbin & Sohy (1998), the method aims at improving the resolution of the data rather than at completely removing the point spread function (PSF). Improvements over the original method include a better regularization of the pixel channel of the image, based on wavelet filtering and multiscale analysis, and a better controlled separation of the point source vs. the extended source. In addition, the method is able to simultaneously deconvolve many individual frames of the same object taken with different instruments under different PSF conditions. For this purpose, we introduce a general geometric transformation between individual images. This transformation allows the combination of the images without having to interpolate them. We illustrate the capability of our algorithm using real and simulated images with complex diffraction-limited PSF.Comment: Accepted in A&A. An application of the technique to real data is available in Cantale et al. http://arxiv.org/abs/1601.05192v

The Spatial Structure of An Accretion Disk

Based on the microlensing variability of the two-image gravitational lens HE1104-1805 observed between 0.4 and 8 microns, we have measured the size and wavelength-dependent structure of the quasar accretion disk. Modeled as a power law in temperature, T proportional to R^-beta, we measure a B-band (0.13 microns in the rest frame) half-light radius of R_{1/2,B} = 6.7 (+6.2 -3.2) x 10^15 cm (68% CL) and a logarithmic slope of beta=0.61 (+0.21 -0.17) for our standard model with a logarithmic prior on the disk size. Both the scale and the slope are consistent with simple thin disk models where beta=3/4 and R_{1/2,B} = 5.9 x 10^15 cm for a Shakura-Sunyaev disk radiating at the Eddington limit with 10% efficiency. The observed fluxes favor a slightly shallower slope, beta=0.55 (+0.03 -0.02), and a significantly smaller size for beta=3/4.Comment: 5 pages, 4 figures, submitted to Ap

A fast empirical method for galaxy shape measurements in weak lensing surveys

We describe a simple and fast method to correct ellipticity measurements of galaxies from the distortion by the instrumental and atmospheric point spread function (PSF), in view of weak lensing shear measurements. The method performs a classification of galaxies and associated PSFs according to measured shape parameters, and corrects the measured galaxy ellipticites by querying a large lookup table (LUT), built by supervised learning. We have applied this new method to the GREAT10 image analysis challenge, and present in this paper a refined solution that obtains the competitive quality factor of Q = 104, without any shear power spectrum denoising or training. Of particular interest is the efficiency of the method, with a processing time below 3 ms per galaxy on an ordinary CPU.Comment: 8 pages, 6 figures. Metric values updated according to the final GREAT10 analysis software (Kitching et al. 2012, MNRAS 423, 3163-3208), no qualitative changes. Associated code available at http://lastro.epfl.ch/megalu

Discovery of a Galaxy Cluster in the Foreground of the Wide-Separation Quasar Pair UM425

We report the discovery of a cluster of galaxies in the field of UM425, a pair of quasars separated by 6.5arcsec. Based on this finding, we revisit the long-standing question of whether this quasar pair is a binary quasar or a wide-separation lens. Previous work has shown that both quasars are at z=1.465 and show broad absorption lines. No evidence for a lensing galaxy has been found between the quasars, but there were two hints of a foreground cluster: diffuse X-ray emission observed with Chandra, and an excess of faint galaxies observed with the Hubble Space Telescope. Here we show, via VLT spectroscopy, that there is a spike in the redshift histogram of galaxies at z=0.77. We estimate the chance of finding a random velocity structure of such significance to be about 5%, and thereby interpret the diffuse X-ray emission as originating from z=0.77, rather than the quasar redshift. The mass of the cluster, as estimated from either the velocity dispersion of the z=0.77 galaxies or the X-ray luminosity of the diffuse emission, would be consistent with the theoretical mass required for gravitational lensing. The positional offset between the X-ray centroid and the expected location of the mass centroid is about 40kpc, which is not too different from offsets observed in lower redshift clusters. However, UM425 would be an unusual gravitational lens, by virtue of the absence of a bright primary lensing galaxy. Unless the mass-to-light ratio of the galaxy is at least 80 times larger than usual, the lensing hypothesis requires that the galaxy group or cluster plays a uniquely important role in producing the observed deflections. Based on observations performed with the Very Large Telescope at the European Southern Observatory, Paranal, Chile.Comment: 12 pages, accepted by ApJ 2005, May 1

X-ray Observations of Gravitationally Lensed Quasars; Evidence for a Hidden Quasar Population

The large magnification factors of gravitationally lensed (GL) systems allow us to investigate the properties of quasars with X-ray luminosities that are substantially lower than those of unlensed ones and also provide an independent means of estimating the contribution of faint quasars to the hard X-ray component of the cosmic X-ray background. Our spectral analysis indicate a flattening of the spectral index in the hard band (4-20keV restframe) for 2 radio-loud quasars in the GL quasar sample for which the data have moderate signal-to-noise ratio. We have identified a large fraction of Broad Absorption Line (BAL) quasars amongst the GL quasar population. We find that approximately 35% of radio-quiet GL quasars contain BAL features which is significantly larger than the 10% fraction of BAL quasars presently found in optically selected flux limited quasar samples. We present a simple model that estimates the effects of attenuation and lens magnification on the luminosity function of quasars and that explains the observed fraction of GL BAL quasars. These observations suggest that a large fraction of BAL quasars are missed from flux limited optical surveys. Modeling of several X-ray observations of the GL BAL quasar PG1115+080 suggests that the observed large X-ray variability may be caused in part by a variable intrinsic absorber consistent with previously observed variability of the BAL troughs in the UV band. The observed large X-ray flux variations in PG1115+080 offer the prospect of considerably reducing errors in determining the time delay with future X-ray monitoring of this system and hence constraining the Hubble constant H$_{0}$.Comment: 19 pages, 10 figures, 9 Tables, accepted for publication in Ap

Deep Imaging of AXJ2019+112: The Luminosity of a ``Dark Cluster''

We detect a distant cluster of galaxies centered on the QSO lens and luminous X-ray source AXJ2019+112, a.k.a. ``The Dark Cluster'' (Hattori et al 1997). Using deep V,I Keck images and wide-field K_s imaging from the NTT, a tight red sequence of galaxies is identified within a radius of 0.2 h^{-1} Mpc of the known z=1.01 elliptical lensing galaxy. The sequence, which includes the central elliptical galaxy, has a slope in good agreement with the model predictions of Kodama et al (1998) for z~1. We estimate the integrated rest-frame luminosity of the cluster to be L_V > 3.2 x 10^{11}h^{-2}L_{\sun} (after accounting for significant extinction at the low latitude of this field), more than an order of magnitude higher than previous estimates. The central region of the cluster is deconvolved using the technique of Magain, Courbin & Sohy (1998), revealing a thick central arc coincident with an extended radio source. All the observed lensing features are readily explained by differential magnification of a radio loud AGN by a shallow elliptical potential. The QSO must lie just outside the diamond caustic, producing two images, and the arc is a highly magnified image formed from a region close to the center of the host galaxy, projecting inside the caustic. The mass--to--light ratio within an aperture of 0.4 h ^{-1} Mpc is M_x/L_V= 224^{+112}_{-78}h(M/L_V)_{\sun}, using the X-ray temperature. The strong lens model yields a compatible value, M/L_V= 372^{+94}_{-94}h(M/L_V)_{\sun}, whereas an independent weak lensing analysis sets an upper limit of M/L_V <520 h(M/L_V)_{\sun}, typical of massive clusters.Comment: AAS Latex format, 24 pages, 9 figures. Fig 1a,b available at http://astro.berkeley.edu/~benitezn/cluster.html . Submitted to ApJ on August 15t

Deconvolution with correct sampling

A new method for improving the resolution of astronomical images is presented. It is based on the principle that sampled data cannot be fully deconvolved without violating the sampling theorem. Thus, the sampled image should not be deconvolved by the total Point Spread Function, but by a narrower function chosen so that the resolution of the deconvolved image is compatible with the adopted sampling. Our deconvolution method gives results which are, in at least some cases, superior to those of other commonly used techniques: in particular, it does not produce ringing around point sources superimposed on a smooth background. Moreover, it allows to perform accurate astrometry and photometry of crowded fields. These improvements are a consequence of both the correct treatment of sampling and the recognition that the most probable astronomical image is not a flat one. The method is also well adapted to the optimal combination of different images of the same object, as can be obtained, e.g., from infrared observations or via adaptive optics techniques.Comment: 22 pages, LaTex file + 10 color jpg and postscript figures. To be published in ApJ, Vol 484 (1997 Feb.

Keck Spectroscopy of Three Gravitational Lens Systems Discovered in the JVAS and CLASS Surveys

We present spectra of three gravitational lens systems taken with the Low Resolution Imaging Spectrograph on the W. M. Keck Telescopes. All of the systems were discovered in the JVAS and CLASS radio surveys, which were designed to find lenses suitable for measuring $H_0$. Previous spectra of these systems had low signal-to-noise ratios, and only one of the source redshifts was secure. Our observations give unambiguous lens and source redshifts for all of the systems, with ($z_l$, $z_s$) = (0.4060,1.339), (0.5990,1.535) and (0.4144,1.589) for B0712+472, B1030+074 and B1600+434, respectively. The observed image splittings in the systems imply that the masses of the lensing galaxies within their Einstein rings are 5.4$\times 10^{10}$, 1.2$\times 10^{11}$, and 6.3\times 10^{10} h^{-1} M_{\sun}. The resulting V-band mass-to-light ratios for B0712+472 and B1030+074, measured inside their Einstein ring radii, are \sim 10h (M/L)_{\sun, V}, slightly higher than values observed in nearby ellipticals. For B1600+434, the mass-to-light ratio is 48h (M/L)_{\sun, V}. This high value can be explained, at least in part, by the prominent dust lane running through the galaxy. Two of the three lens systems show evidence of variability, so monitoring may yield a time delay and thus a measurement of $H_0$.Comment: 8 pages, 5 Postscript Figures, uses aastex. To appear in A.

Simultaneous Estimation of Time Delays and Quasar Structure

We expand our Bayesian Monte Carlo method for analyzing the light curves of gravitationally lensed quasars to simultaneously estimate time delays and quasar structure including their mutual uncertainties. We apply the method to HE1104-1805 and QJ0158-4325, two doubly-imaged quasars with microlensing and intrinsic variability on comparable time scales. For HE1104-1805 the resulting time delay of (Delta t_AB) = t_A - t_B = 162.2 -5.9/+6.3 days and accretion disk size estimate of log(r_s/cm) = 15.7 -0.5/+0.4 at 0.2 micron in the rest frame are consistent with earlier estimates but suggest that existing methods for estimating time delays in the presence of microlensing underestimate the uncertainties. We are unable to measure a time delay for QJ0158-4325, but the accretion disk size is log(r_s/cm) = 14.9 +/- 0.3 at 0.3 micron in the rest frame.Comment: 21 pages, 6 figures, submitted to Ap