S. H. Suyu1, P. J. Marshall2,3, M. W. Auger3,4, S. Hilbert1,5, R. D. Blandford2, L. V. E. Koopmans6, C. D. Fassnacht4 and T. Treu3,7
1
Argelander Institut für Astronomie, Auf dem Hügel 71, 53121 Bonn, Germany
2
Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, P.O. Box 20450, MS 29, Stanford, CA 94309, USA
3
Department of Physics, University of California, Santa Barbara, CA 93106-9530, USA
4
Department of Physics, University of California at Davis, 1 Shields Avenue, Davis, CA 95616, USA
5
Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany
6
Kapteyn Astronomical Institute, P.O. Box 800, 9700AV Groningen, The Netherlands
7
Sloan Fellow, Packard Fellow.
S. H. Suyu et al 2010 ApJ 711 201
Strong gravitational lens systems with measured time delays between the multiple images provide a method for measuring the "time-delay distance" to the lens, and thus the Hubble constant. We present a Bayesian analysis of the strong gravitational lens system B1608+656, incorporating (1) new, deep Hubble Space Telescope (HST) observations, (2) a new velocity-dispersion measurement of 260 ± 15 km s–1 for the primary lens galaxy, and (3) an updated study of the lens' environment. Our analysis of the HST images takes into account the extended source surface brightness, and the dust extinction and optical emission by the interacting lens galaxies. When modeling the stellar dynamics of the primary lens galaxy, the lensing effect, and the environment of the lens, we explicitly include the total mass distribution profile logarithmic slope γ' and the external convergence κext; we marginalize over these parameters, assigning well-motivated priors for them, and so turn the major systematic errors into statistical ones. The HST images provide one such prior, constraining the lens mass density profile logarithmic slope to be γ' = 2.08 ± 0.03; a combination of numerical simulations and photometric observations of the B1608+656 field provides an estimate of the prior for κext: 0.10+0.08 –0.05. This latter distribution dominates the final uncertainty on H 0. Fixing the cosmological parameters at Ωm = 0.3, ΩΛ = 0.7, and w = –1 in order to compare with previous work on this system, we find H 0 = 70.6+3.1 –3.1 km s–1 Mpc–1. The new data provide an increase in precision of more than a factor of 2, even including the marginalization over κext. Relaxing the prior probability density function for the cosmological parameters to that derived from the Wilkinson Microwave Anisotropy Probe (WMAP) five-year data set, we find that the B1608+656 data set breaks the degeneracy between Ωm and ΩΛ at w = –1 and constrains the curvature parameter to be –0.031 < Ωk < 0.009 (95% CL), a level of precision comparable to that afforded by the current Type Ia SNe sample. Asserting a flat spatial geometry, we find that, in combination with WMAP, H 0 = 69.7+4.9 –5.0 km s–1 Mpc–1 and w = –0.94+0.17 –0.19 (68% CL), suggesting that the observations of B1608+656 constrain w as tightly as the current Baryon Acoustic Oscillation data do.
cosmology: observations; distance scale; galaxies: individual (B1608+656); gravitational lensing: strong; methods: data analysis
98.62.Sb Gravitational lenses and luminous arcs
95.30.Sf Relativity and gravitation
95.30.Tg Thermodynamic processes, conduction, convection, equations of state
Issue 1 (2010 March 1)
Received 2009 October 14 , accepted for publication 2010 January 13
Published 2010 February 9
S. H. Suyu et al 2010 ApJ 711 201
M Tanaka et al 2001 Metrologia 38 301
H Costantini et al 2009 Rep. Prog. Phys. 72 086301
Tuomas Lunttila et al 2009 ApJ 702 L37
M H A Hassan and C J H Watson 1977 Plasma Phys. 19 237
Jes K. Jørgensen and Ewine F. van Dishoeck 2010 ApJ 710 L72
Pierre Gibart 2004 Rep. Prog. Phys. 67 667
A. Kopylov and V. Petukhov JCAP08(2009)006
G S Snider 2007 Nanotechnology 18 365202
A Kopylov and V Petukhov JCAP10(2008)007