We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the expected number of quasars lensed by early-type galaxies and their image separation distribution assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant model (w=-1) we obtain \Omega_\Lambda=0.74^{+0.11}_{-0.15}(stat.)^{+0.13}_{-0.06}(syst.). Allowing w to be a free parameter we find \Omega_M=0.26^{+0.07}_{-0.06}(stat.)^{+0.03}_{-0.05}(syst.) and w=-1.1\pm0.6(stat.)^{+0.3}_{-0.5}(syst.) when combined with the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, and provide additional confirmation of the presence of dark energy consistent with a cosmological constant, derived independently of type Ia supernovae.

We report the first results of our systematic search for strongly lensed quasars using the spectroscopically confirmed quasars in the Sloan Digital Sky Survey (SDSS). Among 46,420 quasars from the SDSS Data Release 3 (~4188 deg^2), we select a subsample of 22,683 quasars that are located at redshifts between 0.6 and 2.2 and are brighter than the Galactic extinction corrected i-band magnitude of 19.1. We identify 220 lens candidates from the quasar subsample, for which we conduct extensive and systematic follow-up observations in optical and near-infrared wavebands, in order to construct a complete lensed quasar sample at image separations between 1'' and 20'' and flux ratios of faint to bright lensed images larger than 10^{-0.5}. We construct a statistical sample of 11 lensed quasars. Ten of these are galaxy-scale lenses with small image separations (~1''-2'') and one is a large separation (15'') system which is produced by a massive cluster of galaxies, representing the first statistical sample of lensed quasars including both galaxy- and cluster-scale lenses. The Data Release 3 spectroscopic quasars contain an additional 11 lensed quasars outside the statistical sample.

We report the discovery of the two-image gravitationally lensed quasar SDSS J133222.62+034739.9 (SDSS J1332+0347) with an image separation of Delta_theta=1.14". This system consists of a source quasar at z_s=1.445 and a lens galaxy at z_l=0.191. The agreement of the luminosity, ellipticity and position angle of the lens galaxy with those expected from lens model confirms the lensing hypothesis.

We investigate the dependence of the time delays for the large-separation gravitationally lensed quasar SDSS J1004+4112 on the inner mass profile of the lensing cluster. Adopting the mass model whose innermost density profile is parameterized as \rho \propto r^{-\alpha}, we derive a series of mass models which can fit observational data and then compute the probability distribution functions of time delays. We find that larger \alpha has longer time delays, longer tails at the higher end of the probability distribution, and larger model uncertainties. The ratios of time delays slightly depend on the slope \alpha. Among others, time delays between images C and A (or B) have little dependence on the inner slope, particularly when the time delays are short. The dependence of time delays on \alpha is well fitted by a linear form, which reflects well-known degeneracy between the mass profile and time delays. We perform a Monte-Carlo simulation to illustrate how well the inner slope can be constrained from measurements of time delays. We find that measurements of more than one time delays result in reasonably tight constraints on the inner slope (\sigma_{\alpha} \lesssim 0.25), while only one time delay cannot determine the inner slope very well. Our result indicates that time delays indeed serve as a powerful tool to determine the mass profile, despite the complexity of the lensing cluster.

We report the discoveries of two, two-image gravitationally lensed quasars selected from the Sloan Digital Sky Survey: SDSS J0806+2006 at z_s=1.540 and SDSS J1353+1138 at z_s=1.629 with image separations of 1.40" and 1.41" respectively. Spectroscopic and optical/near-infrared imaging follow-up observations show that the quasar images have identical redshifts and possess extended objects between the images that are likely to be lens galaxies at z_l~0.6 in SDSS J0806+2006 and z_l~0.3 in SDSS J1353+1138. The field of SDSS J0806+2006 contains several nearby galaxies that may significantly perturb the system, and SDSS J1353+1138 has an extra component near its Einstein ring that is probably a foreground star. Simple mass models with reasonable parameters reproduce the quasar positions and fluxes of both systems.

In modeling strong gravitational lens systems, one often adopts simple models, such as singular isothermal elliptical density plus lowest-order external perturbation. However, such simple models may mislead us if the real mass distribution is more complicated than that in the assumed models. In particular, assumptions on mass models are crucial in studying flux ratio anomalies that have been suggested as evidence for a cold dark-matter substructure. We reinvestigated four quadruple lens systems using power-law Fourier models, which have advantages of clear physical meanings and the applicability of a linear method, as well as a simple singular isothermal elliptical density model. We also investigated the effect of external perturbations, including a singular isothermal sphere, lowest-order expansion, and next-order expansion. We have found that the $\cos3\theta$ terms of the primary galaxy and/or of external perturbation significantly reduce the \chi^2 in PG 1115+080 and B1422+231. In particular, we could reproduce the flux ratios of B1422+231 with next-order external perturbation assuming 5% flux uncertainties, suggesting that external perturbation cannot be described by a simple singular isothermal sphere approximation. On the other hand, we could not fit B0712+472 and B2045+265 very well even with our models, although the \chi^2 were reduced compared with the case of using simple models. Our results clearly demonstrate that both the primary lens galaxy and the external perturbation are often more complicated than we usually assume.

An attempt to measure the Hubble constant with gravitational lens time delays is often limited by the strong degeneracy between radial mass profiles of lens galaxies and the Hubble constant. We show that strong gravitational lensing of type Ia supernovae breaks this degeneracy; the standard candle nature of type Ia supernova luminosity function allows one to measure the magnification factor directly, and this information is essential to constrain radial mass profiles and the Hubble constant separately. Our numerical simulation demonstrates that the Hubble constant can be determined with \sim 5% accuracy from only several lens events if magnification factors are used as constraints. Therefore, distant supernova survey is a promising way to measure the global Hubble constant independently with the local estimates.