We present a refined gravitational lens model of the four-image lens system B1608+656 based on new and improved observational constraints: (i) the three independent time-delays and flux-ratios from VLA observations, (ii) the radio-image positions from VLBA observations, (iii) the shape of the deconvolved Einstein Ring from optical and infrared HST images, (iv) the extinction-corrected lens-galaxy centroids and structural parameters, and (v) a stellar velocity dispersion, sigma_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 which are the key parameters to determine the value of H_0 from lens time delays. This has required the development of a new lens code that is highly optimized for speed. The minimum-chi^2 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_0=75(+7/-6) km/s/Mpc (68 percent CL; Omega_m=0.3, Omega_Lambda=0.7. The non-linear 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_0. The lens galaxy G1 is ~5 times more massive than the secondary lens galaxy (G2), and has a mass density slope of gamma_G1=2.03(+0.14/-0.14) +- 0.03 (68 percent CL) for rho~r^-gamma', very close to isothermal (gamma'=2). (Abridged)

We present a multi-wavelength analysis of high-resolution observations of the quadruple lens B1608+656 from the HST archive, acquired with WFPC2 through filters F606W (V-band) and F814W (I-band), and with NIC1 in filter F160W (H-band). In the three bands, the observations show extended emission from the four images of the source in a ring-like configuration that surrounds the two, resolved, lensing galaxies. B1608+656 was discovered as a double-lobed radio source, and later identified as a post-starburst galaxy in the optical. Based on photometry and optical spectroscopy we estimate that the stellar population of the source has an age of $\sim$ 500 Myr. This provides a model for the spectrum of the source that extends over spectral regions where no observations are available, and is used to generate Tiny Tim PSFs for the filters. Deconvolutions performed with the Lucy-Richardson method are presented, and the limitations of these restorations is discussed. V$-$I and I$-$H color maps show evidence of extinction by dust associated with one of the lensing galaxies, a late type galaxy presumably disrupted after its close encounter with the other lens, an elliptical galaxy. The extinction affects the two lens galaxies and two of the four multiple images. The diagnostic of wavelength-dependent effects in the images shows that corrections for contamination with light from the lenses, extinction, and PSF convolution need to be applied before using the extended structure in the images as a constraint on lens models. We will present the restoration of the images in a subsequent paper.

We apply elliptical isothermal mass models to reproduce the point source properties, i.e. image positions, flux density ratios and time delay ratios, of the quadruple lens B1608+656. A wide set of suitable solutions is found, showing that models that only fit the properties of point sources are under-constrained and can lead to a large uncertainty in the prediction of Ho. We present two examples of models predicting Ho=100 km/(s Mpc) (chi^2=4) and Ho=69 km/(s Mpc) (chi^2=24).

In order to use a gravitational lens to measure the Hubble constant accurately, it is necessary to derive a reliable model of the lens surface potential. If the analysis is restricted to the locations and magnifications of point images, the derived Hubble constant depends upon the class of mass models used to fit the data. However, when there is extended emission from an Einstein ring, it may be possible to derive a potential from the observed surface brightness in a model-independent manner. This procedure is illustrated with reference to B1608+656. The multi-band images are de-reddened, de-convolved and de-contaminated so that the luminous matter and the surface brightness contours in the Einstein ring are both faithfully mapped. This intensity distribution can then be used to reconstruct the potential. Progress in implementing this program is reported. The observed incidence of multiple-imaged galaxies on the Hubble Deep Fields is an order of magnitude smaller than naively predicted on the basis of radio lens surveys, like CLASS, but consistent with the rate computed using surface photometry of candidate lens galaxies assuming standard mass to light ratios. In order to resolve this paradox, it is suggested that most galaxy lenses are located in compact groups.

We estimate the effects of weak lensing by large-scale density inhomogeneities and long-wavelength gravitational waves upon the polarization properties of electromagnetic radiation as it propagates from cosmologically distant sources. Scalar (density) fluctuations do not rotate neither the plane of polarization of the electromagnetic radiation nor the source image. They produce, however, an appreciable shear, which distorts the image shape, leading to an apparent rotation of the image orientation relative to its plane of polarization. In sources with large ellipticity the apparent rotation is rather small, of the order (in radians) of the dimensionless shear. The effect is larger at smaller source eccentricity. A shear of 1% can induce apparent rotations of around 5 degrees in radio sources with the smallest eccentricity among those with a significant degree of integrated linear polarization. We discuss the possibility that weak lensing by shear with rms value around or below 5% may be the cause for the dispersion in the direction of integrated linear polarization of cosmologically distant radio sources away from the perpendicular to their major axis, as expected from models for their magnetic fields. A rms shear larger than 5% would be incompatible with the observed correlation between polarization properties and source orientation in distant radio galaxies and quasars. Gravity waves do rotate both the plane of polarization as well as the source image. Their weak lensing effects, however, are negligible.

We analyse the effects of large-scale inhomogeneities upon the observables of a gravitational lens system, focusing on the issue of whether large-scale structure imperils the program to determine the Hubble parameter through measurements of time delays between multiple images in lens systems. We find that the lens equation in a spatially flat Robertson-Walker cosmology with scalar metric fluctuations is equivalent to that for the same lensing system in the absence of fluctuations, but with a different angular position of the source relative to the lens axis. Since the absolute position of the source is not observable, gravitational lens measurements cannot directly reveal the presence of large scale structure. Large-scale perturbations do not modify the functional relationship between observable lens parameters and the Hubble parameter, and therefore do not seriously affect the determination of $H_0$ from lens time delays.

Using Fermat's principle, we analyze the effects of very long wavelength gravitational waves upon the images of a gravitationally lensed quasar. We show that the lens equation in the presence of gravity waves is equivalent to that of a lens with different alignment between source, deflector, and observer in the absence of gravity waves. Contrary to a recent claim, we conclude that measurements of time delays in gravitational lenses cannot serve as a method to detect or constrain a stochastic background of gravitational waves of cosmological wavelengths, because the wave-induced time delay is observationally indistinguishable from an intrinsic time delay due to the lens geometry.