2003ApJ...599...70K

We present a refined gravitational lens model of the four-image lens system B1608+656 based on new and improved observational constraints: (1) the three independent time delays and flux ratios from Very Large Array observations, (2) the radio-image positions from Very Large Baseline Array observations, (3) the shape of the deconvolved Einstein ring from optical and infrared Hubble Space Telescope images, (4) the extinction-corrected lens-galaxy centroids and structural parameters, and (5) a stellar velocity dispersion, σ_ap=247±35 km/s, of the primary lens galaxy (G1), obtained from an echelle spectrum taken with the Keck II Telescope. The lens-mass model consists of two elliptical mass distributions with power-law density profiles and an external shear, totaling 22 free parameters, including the density slopes that are the key parameters for determining the value of H₀from lens time delays. This has required the development of a new lens code that is highly optimized for speed. The minimum-χ² model reproduces all observations very well, including the stellar velocity dispersion and the shape of the Einstein ring. A combined gravitational lens and stellar dynamical analysis leads to a value of the Hubble constant of H₀=75⁺⁷_–6km/s/Mpc (68% CL; Ω_m=0.3, Ω_Λ=0.7). The nonlinear error analysis includes correlations between all free parameters, in particular the density slopes of G1 and G2, yielding an accurate determination of the random error on H₀. The lens galaxy G1 is ∼5 times more massive than the secondary lens galaxy (G2) and has a mass density slope of γ'_G1=2.03^+0.14_–0.14±0.03 (68% CL) for ρ∝r^–γ', very close to isothermal (γ'=2). After extinction correction, G1 exhibits a smooth surface brightness distribution with an R^1/4 profile and no apparent evidence for tidal disruption by interactions with G2. Given the scope of the observational constraints and the gravitational lens models, as well as the careful corrections to the data, we believe this value of H₀to be little affected by known systematic errors (≲5%).