In this catalog we compile information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) as well as all other bright ($m_{peak}\leq17$), spectroscopically confirmed supernovae found in 2017, totaling 308 supernovae. We also present UV through near-IR magnitudes gathered from public databases of all host galaxies for the supernovae in the sample. We perform statistical analyses of our full bright supernova sample, which now contains 949 supernovae discovered since 2014 May 1, including supernovae from our previous catalogs. This is the fourth of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team, and this work presents updated data and measurements, including light curves, redshifts, classifications, and host galaxy identifications, that supersede information contained in any previous publications.

We devise a test of nonlinear departures from general relativity (GR) using time delays in strong gravitational lenses. We use a phenomenological model of gravitational screening as a step discontinuity in the measure of curvature per unit mass, at a radius $\Lambda$. The resulting slip between two scalar gravitational potentials leads to a shift in the apparent positions and time delays of lensed sources, relative to the GR predictions, of size $\gamma_{\rm PN}-1$. As a proof of principle, we use measurements of two lenses, RXJ1121-1231 and B1608+656, to constrain deviations from GR to be below $|\gamma_{\rm PN}-1| \leq 0.2 \times (\Lambda/100\, \rm kpc)$. These constraints are complementary to other current probes, and are the tightest in the range $\Lambda=10-200$ kpc, showing that future measurements of strong-lensing time delays have great promise to seek departures from general relativity on kpc-Mpc scales.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models.

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal).

The Dark Ages, probed by the redshifted 21-cm signal, is the ideal epoch for a new rigorous test of the standard LCDM cosmological model. Divergences from that model would indicate new physics, such as dark matter decay (heating) or baryonic cooling beyond that expected from adiabatic expansion of the Universe. In the early Universe, most of the baryonic matter was in the form of neutral hydrogen (HI), detectable via its ground state's spin-flip transition. A measurement of the redshifted 21-cm spectrum maps the history of the HI gas through the Dark Ages and Cosmic Dawn and up to the Epoch of Reionization (EoR). The Experiment to Detect the Global EoR Signature (EDGES) recently reported an absorption trough at 78 MHz (redshift z of 17), similar in frequency to expectations for Cosmic Dawn, but about 3 times deeper than was thought possible from standard cosmology and adiabatic cooling of HI. Interactions between baryons and slightly-charged dark matter particles with electron-like mass provide a potential explanation of this difference but other cooling mechanisms are also being investigated to explain these results. The Cosmic Dawn trough is affected by cosmology and the complex astrophysical history of the first luminous objects. Another trough is expected during the Dark Ages, prior to the formation of the first stars and thus determined entirely by cosmological phenomena (including dark matter). Observations on or in orbit above the Moon's farside can investigate this pristine epoch (15-40 MHz; z=100-35), which is inaccessible from Earth. A single cross-dipole antenna or compact array can measure the amplitude of the 21-cm spectrum to the level required to distinguish (at >5$\sigma$}) the standard cosmological model from that of additional cooling derived from current EDGES results. This observation constitutes a powerful, clean probe of exotic physics in the Dark Ages.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

Recent observations of hydrogen absorption that occurred when the first stars turned on may give insights into the nature of dark matter.

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

We use single-epoch spectroscopy of three gravitationally lensed quasars, HE0435-1223, WFI2033-4723, and HE2149-2745, to study their inner structure (BLR and continuum source). We detect microlensing-induced magnification in the wings of the broad emission lines of two of the systems (HE0435-1223 and WFI2033-4723). In the case of WFI2033-4723, microlensing affects two "bumps" in the spectra which are almost symmetrically arranged on the blue (coincident with an AlIII emission line) and red wings of CIII]. These match the typical double-peaked profile that follows from disk kinematics. The presence of microlensing in the wings of the emission lines indicates the existence of two different regions in the BLR: a relatively small one with kinematics possibly related to an accretion disk, and another one that is substantially more extended and insensitive to microlensing. There is good agreement between the estimated size of the region affected by microlensing in the emission lines, $r_s=10^{+15}_{-7} \sqrt{M/M_{\odot}}$ light-days (red wing of CIV in HE0435-1223) and $r_s=11^{+28}_{-7} \sqrt{M/M_{\odot}}$ light-days (CIII] bumps in WFI2033-4723) with the sizes inferred from the continuum emission, $r_s=13^{+5}_{-4} \sqrt{M/M_{\odot}}$ light-days (HE0435-1223) and $r_s=10^{+3}_{-2} \sqrt{M/M_{\odot}}$ light-days (WFI2033-4723). For HE2149-2745 we measure an accretion disk size $r_s=8^{+11}_{-5} \sqrt{M/M_{\odot}}$ light-days. The estimates of $p$, the exponent of the size vs. wavelength ($r_s\propto\lambda^p$), are $1.2\pm0.6$, $0.8\pm0.2$, and $0.4\pm0.3$ for HE0435-1223, WFI2033-4723, and HE2149-2745, respectively. In conclusion, the continuum microlensing amplitude in the three quasars and chromaticity in WFI2033-4723 and HE2149-2745 are below expectations for the thin disk model. The disks are larger and their temperature gradients are flatter than predicted by this model.

We present new mid-IR observations of the quadruply lensed quasar Q2237+0305 taken with CanariCam on the Gran Telescopio Canarias. Mid-IR emission by hot dust, unlike the optical and near-IR emission from the accretion disk, is unaffected by the interstellar medium (extinction/scattering) or stellar microlensing. We compare these "true" ratios to the (stellar) microlensed flux ratios observed in the optical/near-IR to constrain the structure of the quasar accretion disk. We find a half-light radius of $R_{1/2}=3.4_{-2.1}^{+5.3}\sqrt{\langle M \rangle/0.3\,\rm{M_{\odot}}}$ light-days at $\lambda_{rest}=1736$ {\AA}, and an exponent for the temperature profile $R \propto \lambda^{p}$ of $p=0.79\pm0.55$, where $p=4/3$ for a standard thin-disk model. If we assume that the differences in the mid-IR flux ratios measured over the years are due to microlensing variability, we find a lower limit for the size of the mid-IR-emitting region of $R_{1/2} \gtrsim 200\,\sqrt{\langle M \rangle/0.3\,\rm{M_{\odot}}}$ light-days. We also test for the presence of substructure/satellites by comparing the observed mid-IR flux ratios with those predicted from smooth lens models. We can explain the differences if the surface density fraction in satellites near the lensed images is $\alpha = 0.033_{-0.019}^{+0.046}$ for a singular isothermal ellipsoid plus external shear mass model or $\alpha = 0.013_{-0.008}^{+0.019}$ for a mass model combining ellipsoidal NFW and de Vaucouleurs profiles in an external shear.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

We propose to use the flux variability of lensed quasar images induced by gravitational microlensing to measure the transverse peculiar velocity of lens galaxies over a wide range of redshift. Microlensing variability is caused by the motions of the observer, the lens galaxy (including the motion of the stars within the galaxy), and the source; hence, its frequency is directly related to the galaxy's transverse peculiar velocity. The idea is to count time-event rates (e.g., peak or caustic crossing rates) in the observed microlensing light curves of lensed quasars that can be compared with model predictions for different values of the transverse peculiar velocity. To compensate for the large time-scale of microlensing variability we propose to count and model the number of events in an ensemble of gravitational lenses. We develop the methodology to achieve this goal and apply it to an ensemble of 17 lensed quasar systems. In spite of the shortcomings of the available data, we have obtained tentative estimates of the peculiar velocity dispersion of lens galaxies at $z\sim 0.5$, $\sigma_{\rm pec}(0.53\pm0.18)\simeq(638\pm213)\sqrt{\langle m \rangle/0.3 M_\odot} \, \rm km\, s^{-1}$. Scaling at zero redshift we derive, $\sigma_{\rm pec}(0)\simeq(491\pm164) \sqrt{\langle m \rangle/0.3 M_\odot} \, \rm km\, s^{-1}$, consistent with peculiar motions of nearby galaxies and with recent $N$-body nonlinear reconstructions of the Local Universe based on $\Lambda$CDM. We analyze the different sources of uncertainty of the method and find that for the present ensemble of 17 lensed systems the error is dominated by Poissonian noise, but that for larger ensembles the impact of the uncertainty on the average stellar mass may be significant.

We present estimates for the size and the logarithmic slope of the disk temperature profile of the lensed quasar Q2237+0305, independent of the component velocities. These estimates are based on six epochs of multi-wavelength narrowband images from the Nordic Optical Telescope. For each pair of lensed images and each photometric band, we determine the microlensing amplitude and chromaticity using pre-existing mid-IR photometry to define the baseline for no microlensing magnification. A statistical comparison of the combined microlensing data (6 epochs $\times$ 5 narrow bands $\times$ 6 image pairs) with simulations based on microlensing magnification maps gives Bayesian estimates for the half-light radius of $R_{1/2}=8.5^{+7.5}_{-4.0}\sqrt{ \langle M \rangle/0.3\, M_\odot}$ light-days, and $p=0.95\pm0.33$ for the exponent of the logarithmic temperature profile $T\propto R^{ -1/p}$. This size estimate is in good agreement with most recent studies. Other works based on the study of single microlensing events predict smaller sizes, but could be statistically biased by focusing on high-magnification events.

We investigate a possible method for determining the progenitors of black hole (BH) mergers observed via their gravitational wave (GW) signal. We argue that measurements of the cross-correlation of the GW events with overlapping galaxy catalogs may provide an additional tool in determining if BH mergers trace the stellar mass of the Universe, as would be expected from mergers of the endpoints of stellar evolution. If on the other hand the BHs are of primordial origin, as has been recently suggested, their merging would be preferentially hosted by lower biased objects, and thus have a lower cross-correlation with luminous galaxies. Here we forecast the expected precision of the cross-correlation measurement for current and future GW detectors such as LIGO and the Einstein Telescope. We then predict how well these instruments can distinguish the model that identifies high-mass BH-BH mergers as the merger of primordial black holes that constitute the dark matter in the Universe from more traditional astrophysical sources.

We study three high magnification microlensing events, generally recognized as probable caustic crossings, in the optical light curves of the multiply imaged quasar Q 2237+0305. We model the light curve of each event as the convolution of a standard thin disk luminosity profile with a straight fold caustic. We also allow for a linear gradient that can account for an additional varying background effect of microlensing. This model not only matches noticeably well the global shape of each of the three independent microlensing events but also gives remarkably similar estimates for the disk size parameter. The measured average half-light radius, $R_{1/2}=(3.0\pm 1.5)\sqrt{M/0.3M\odot}$ light-days, agrees with previous estimates. In the three events, the core of the magnification profile exhibits "fine structure" related to the innermost region of the accretion disk (located at a radial distance of $2.7\pm 1.4$ Schwarzschild radii according to our measurement). Relativistic beaming at the internal rim of the accretion disk can explain the shape and size of the fine structure, although alternative explanations are also possible. This is the first direct measurement of the size of a structure, likely the innermost stable circular orbit, at $\sim 3$ Schwarzschild radii in a quasar accretion disk. The monitoring of thousands of lensed quasars with future telescopes will allow the study of the event horizon environment of black holes in hundreds of quasars in a wide range of redshifts $(0.5<z<5)$.

For all exoplanet candidates, the reliability of a claimed detection needs to be assessed through a careful study of systematic errors in the data to minimize the false positives rate. We present a method to investigate such systematics in microlensing datasets using the microlensing event OGLE-2013-BLG-0446 as a case study. The event was observed from multiple sites around the world and its high magnification (A_{max} \sim 3000) allowed us to investigate the effects of terrestrial and annual parallax. Real-time modeling of the event while it was still ongoing suggested the presence of an extremely low-mass companion (\sim 3M_\oplus ) to the lensing star, leading to substantial follow-up coverage of the light curve. We test and compare different models for the light curve and conclude that the data do not favour the planetary interpretation when systematic errors are taken into account.

Constraints to the parameters of inflation models are often derived assuming some plausible range for the number--e.g., $N_k=46$ to $N_k=60$--of $e$-folds of inflation that occurred between the time that our current observable Universe exited the horizon and the end of inflation. However, that number is, for any specific inflaton potential, related to an effective equation-of-state parameter $w_{\mathrm{re}}$ and temperature $T_{\mathrm{re}}$, for reheating. Although the physics of reheating is highly uncertain, there is a finite range of reasonable values for $w_{\mathrm{re}}$. Here we show that by restricting $w_{\mathrm{re}}$ to this range, more stringent constraints to inflation-model parameters can be derived than those obtained from the usual procedure. To do so, we focus in this work in particular on natural inflation and inflation with a Higgs-like potential, and on power law models as limiting cases of those. As one example, we show that the lower limit to the tensor-to-scalar ratio $r$, derived from current measurements of the scalar spectral index, is about 20%-25% higher (depending on the model) with this procedure than with the usual approach.

Although many models have been proposed, the physical mechanisms responsible for the formation of low-mass brown dwarfs are poorly understood. The multiplicity properties and minimum mass of the brown-dwarf mass function provide critical empirical diagnostics of these mechanisms. We present the discovery via gravitational microlensing of two very low-mass, very tight binary systems. These binaries have directly and precisely measured total system masses of 0.025 Msun and 0.034 Msun, and projected separations of 0.31 AU and 0.19 AU, making them the lowest-mass and tightest field brown-dwarf binaries known. The discovery of a population of such binaries indicates that brown dwarf binaries can robustly form at least down to masses of ~0.02 Msun. Future microlensing surveys will measure a mass-selected sample of brown-dwarf binary systems, which can then be directly compared to similar samples of stellar binaries.

We present the analysis of the gravitational microlensing event OGLE-2011-BLG-0251. This anomalous event was observed by several survey and follow-up collaborations conducting microlensing observations towards the Galactic Bulge. Based on detailed modelling of the observed light curve, we find that the lens is composed of two masses with a mass ratio q=1.9 x 10^-3. Thanks to our detection of higher-order effects on the light curve due to the Earth's orbital motion and the finite size of source, we are able to measure the mass and distance to the lens unambiguously. We find that the lens is made up of a planet of mass 0.53 +- 0.21,M_Jup orbiting an M dwarf host star with a mass of 0.26 +- 0.11 M_Sun. The planetary system is located at a distance of 2.57 +- 0.61 kpc towards the Galactic Centre. The projected separation of the planet from its host star is d=1.408 +- 0.019, in units of the Einstein radius, which corresponds to 2.72 +- 0.75 AU in physical units. We also identified a competitive model with similar planet and host star masses, but with a smaller orbital radius of 1.50 +- 0.50 AU. The planet is therefore located beyond the snow line of its host star, which we estimate to be around 1-1.5 AU.

We present an analysis of the anomalous microlensing event, MOA-2010-BLG-073, announced by the Microlensing Observations in Astrophysics survey on 2010-03-18. This event was remarkable because the source was previously known to be photometrically variable. Analyzing the pre-event source lightcurve, we demonstrate that it is an irregular variable over time scales >200d. Its dereddened color, $(V-I)_{S,0}$, is 1.221$\pm$0.051mag and from our lens model we derive a source radius of 14.7$\pm$1.3 $R_{\odot}$, suggesting that it is a red giant star. We initially explored a number of purely microlensing models for the event but found a residual gradient in the data taken prior to and after the event. This is likely to be due to the variability of the source rather than part of the lensing event, so we incorporated a slope parameter in our model in order to derive the true parameters of the lensing system. We find that the lensing system has a mass ratio of q=0.0654$\pm$0.0006. The Einstein crossing time of the event, $T_{\rm{E}}=44.3$\pm$0.1d, was sufficiently long that the lightcurve exhibited parallax effects. In addition, the source trajectory relative to the large caustic structure allowed the orbital motion of the lens system to be detected. Combining the parallax with the Einstein radius, we were able to derive the distance to the lens, $D_L$=2.8$\pm$0.4kpc, and the masses of the lensing objects. The primary of the lens is an M-dwarf with $M_{L,p}$=0.16$\pm0.03M_{\odot}$ while the companion has $M_{L,s}$=11.0$\pm2.0M_{\rm{J}}$ putting it in the boundary zone between planets and brown dwarfs.

The Galactic bulge source MOA-2010-BLG-523S exhibited short-term deviations from a standard microlensing lightcurve near the peak of an Amax ~ 265 high-magnification microlensing event. The deviations originally seemed consistent with expectations for a planetary companion to the principal lens. We combine long-term photometric monitoring with a previously published high-resolution spectrum taken near peak to demonstrate that this is an RS CVn variable, so that planetary microlensing is not required to explain the lightcurve deviations. This is the first spectroscopically confirmed RS CVn star discovered in the Galactic bulge.

We develop a radiation pressure-balanced model for the interstellar medium of high-redshift galaxies that describes many facets of galaxy formation at z>~6, including star formation rates and distributions and gas accretion onto central black holes. We first show that the vertical gravitational force in the disk of such a model is dominated by the disk self-gravity supported by the radiation pressure of ionizing starlight on gas. Constraining our model to reproduce the UV luminosity function of Lyman-break galaxies (LBGs), we limit the available parameter-space to wind mass-loading factors 1--4 times the canonical value for momentum-driven winds. We then focus our study by exploring the effects of different angular momentum transport mechanisms in the galactic disk and find that accretion driven by gravitational torques, such as from linear spiral waves or non-linear orbit crossings, can build up black hole masses by z=6 consistent with the canonical M-sigma relation with a duty cycle of unity, while accretion mediated by a local viscosity such as in an alpha-disk results in negligible BH accretion. Both gravitational torque models produce X-ray emission from active galactic nuclei (AGN) in high-redshift LBGs in excess of the estimated contribution from high-mass X-ray binaries. Using a recent analysis of deep Chandra observations by Cowie et al., we can already begin to rule out the most extreme regions of our parameter-space: the inflow velocity of gas through the disk must either be less than one percent of the disk circular velocity or the X-ray luminosity of the AGN must be substantially obscured. Moderately deeper future observations or larger sample sizes will be able to probe the more reasonable range of angular momentum transport models and obscuring geometries.

We report on a program of spectroscopic observations of gravitationally-lensed QSOs with multiple images. We seek to establish whether microlensing is occurring in each QSO image using only single-epoch observations. We calculate flux ratios for the cores of emission lines in image pairs to set a baseline for no microlensing. The offset of the continuum flux ratios relative to this baseline yields the microlensing magnification free from extinction, as extinction affects the continuum and the lines equally. When we find chromatic microlensing, we attempt to constrain the size of the QSO accretion disk. SDSSJ1004+4112 and HE1104-1805 show chromatic microlensing with amplitudes $0.2< |\Delta m| < 0.6$ and $0.2< |\Delta m| < 0.4$ mag, respectively. Modeling the accretion disk with a Gaussian source ($I\propto \exp(-R^2/2r_s^2)$) of size $r_s\propto \lambda^p$ and using magnification maps to simulate microlensing we find $r_s(\lambda 3363)=7\pm3 light-days (18.1\pm7.8 \times 10^{15} cm$) and $p=1.1\pm 0.4$ for SDSS1004+4112, and $r_s(\lambda 3363)=6\pm2 light-days (15.5\pm5.2 \times 10^{15} cm$) and $p=0.7\pm0.1$ for HE1104-1805. For SDSSJ1029+2623 we find strong chromaticity of $\sim 0.4$ mag in the continuum flux ratio, which probably arises from microlensing although not all the available data fit within this explanation. For Q0957+561 we measure B-A magnitude differences of 0.4 mag, much greater than the $\sim$0.05 mag amplitude usually inferred from lightcurve variability. It may substantially modify the current interpretations of microlensing in this system, likely favoring the hypothesis of smaller sources and/or larger microdeflectors. For HS0818+1227, our data yield posible evidence of microlensing.

Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift.

The gravitational lens system MG J0414+0534 is formed by an elliptical galaxy at redshift ~0.96 and a quasar at z~2.64. The system geometry is typical of lensing by an elliptical galaxy with the QSO close to and inside a fold caustic. It shows 4 images of the background source, and a partial Einstein ring is visible at optical wavelengths. It was observed with a global-VLBI array at 18 cm in June 2008. We present here the imaging results and a preliminary lens model constrained by these observations.

We estimate the fraction of mass that is composed of compact objects in gravitational lens galaxies. This study is based on microlensing measurements (obtained from the literature) of a sample of 29 quasar image pairs seen through 20 lens galaxies. We determine the baseline for no microlensing magnification between two images from the ratios of emission line fluxes. Relative to this baseline, the ratio between the continua of the two images gives the difference in microlensing magnification. The histogram of observed microlensing events peaks close to no magnification and is concentrated below 0.6 magnitudes, although two events of high magnification, $\Delta m \sim 1.5$, are also present. We study the likelihood of the microlensing measurements using frequency distributions obtained from simulated microlensing magnification maps for different values of the fraction of mass in compact objects, $\alpha$. The concentration of microlensing measurements close to $\Delta m \sim 0$ can be explained only by simulations corresponding to very low values of $\alpha$ (10% or less). A maximum likelihood test yields $\alpha=0.05_{-0.03}^{+0.09}$ (90% confidence interval) for a quasar continuum source of intrinsic size $r_{s_0}\sim 2.6 \cdot 10^{15} \rm cm$. Regarding the current controversy about Milky Way/LMC and M31 microlensing studies, our work supports the hypothesis of a very low content in MACHOS (Massive Compact Halo Objects).

The influence of microlensing in the profiles of the emission lines generated in a biconical geometry is discussed. Microlensing amplification in this anisotropic model is not directly related to the bicone's intrinsic size but depends on the orientation of the bicone axis and on the cone aperture. The orientation of the projected bicone with respect to the shear of the magnification pattern can induce very interesting effects, like the quasi-periodic enhancements of the red/blue part of the emission line profile or the lack of correlation between the broad line region (BLR) and continuum light curves of QSOs. The emission line profiles of a BLR moving in a high caustic concentration exhibit sharp features that are well defined in wavelength. These features (spikes) correspond to the scanning of the kinematics of the BLR by the caustic clusters. The biconical model can qualitatively reproduce with a transversal (with respect to the shear) movement of the BLR, the recurrent blue-wing enhancement detected in the emission line profile of the A image of the quasar lensed system SDSS J1004+4112. The probability of observing this repetitive event is almost a 2% for a fraction of matter in stars of a 5%. This result would make plausible the detection of the spectral variability in SDSS J1004+4112 under the hypothesis of microlensing of a bicone.

We present integral field spectroscopic observations of the quadruple-lensed QSO SDSS 1004+4112 taken with the fiber system INTEGRAL at the William Herschel Telescope on 2004 January 19. In May 2003 a blueward enhancement in the high ionization lines of SDSS 1004+4112A was detected and then faded. Our observations are the first to note a second event of similar characteristics less than one year after. Although initially attributed to microlensing, the resemblance among the spectra of both events and the absence of microlensing-induced changes in the continuum of component A are puzzling. The lack of a convincing explanation under the microlensing or intrinsic variability hypotheses makes the observed enhancements particularly relevant, calling for close monitoring of this object.

We consider whether the variability properties of Gamma-Ray Bursts (GRBs) that produce bright optical and longer wavelength transient afterglows (A-GRBs) are the same as a larger, inclusive sample of bright, long-duration GRBs, selected only by their gamma-ray emission. This sample may include a significant population of physically distinct ``dark'' or ``faint-afterglow'' GRBs with different variability properties, or may be composed of a single population, some of which lack afterglows only because of observational selection effects. We argue that the structure function is the most appropriate method for measuring the variability of bursts because of their transient and aperiodic nature. We define a simple statistic: the ratio of the integrated structure function from 0.1 to 1 s compared to that from 0.1 to 10 s, as measured in the observer frame. To avoid instrumental effects we restrict our analysis to GRBs with BATSE data. Comparing 10 A-GRBs to a ``main'' sample of about 500 bursts, we find there is a probability of only 0.03 of the samples being drawn from the same population, with the A-GRBs tending to have relatively less power on sub-second timescales. We conclude that this result is tentative evidence for variations in the properties of GRB progenitors that affect both the gamma-ray and afterglow properties of long-duration GRBs. In addition, our method of analyzing variability identifies a characteristic timescale of ~1 s, below which variability is suppressed, and finds a trend of increased short timescale variability at higher gamma-ray energies. The long-duration GRBs that we identify as having the most sub-second timescale variability, may be relatively bright examples of short-duration GRBs.

We present VR magnification ratios of QSO 0957+561, which are inferred from the GLITP light curves of Q0957+561A and new frames taken with the 2.56m Nordic Optical Telescope about 14 months after the GLITP monitoring. From two photometric approaches and a reasonable range for the time delay in the system (415-430 days), we do not obtain achromatic optical continuum ratios, but ratios depending on the wavelength. These new measurements are consistent with differential extinction in the lens galaxy, the Lyman limit system, the damped Ly-alpha system, or the host galaxy of the QSO. The possible values for the differential extinction and the ratio of total to selective extinction in the V band are reasonable. Moreover, crude probability arguments suggest that the ray paths of the two components cross a similar dusty environment, including a network of compact dust clouds and compact dust voids. As an alternative (in fact, the usual interpretation of the old ratios), we also try to explain the new ratios as caused by gravitational microlensing in the deflector. From magnification maps for each of the gravitationally lensed images, using different fractions of the surface mass density represented by the microlenses, as well as different sizes and profiles of the V-band and R-band sources, several synthetic distributions of V-band and R-band ratios are derived. In some gravitational scenarios, there is an apparent disagreement between the observed pair of ratios and the simulated distributions. However, several microlensing pictures work well. To decide between either extinction, or microlensing, or a mixed scenario (extinction + microlensing), new observational and interpretation efforts are required.

We estimate the dust extinction laws in two intermediate redshift galaxies. The dust in the lens galaxy of LBQS1009-0252, which has an estimated lens redshift of zl~0.88, appears to be similar to that of the SMC with no significant feature at 2175 A. Only if the lens galaxy is at a redshift of zl~0.3, completely inconsistent with the galaxy colors, luminosity or location on the fundamental plane, can the data be fit with a normal Galactic extinction curve. The dust in the zl=0.68 lens galaxy for B0218+357, whose reddened image lies behind a molecular cloud, requires a very flat ultraviolet extinction curve with (formally) R(V)=12 +- 2. Both lens systems seem to have unusual extinction curves by Galactic standards.

The GLITP collaboration observed the first gravitational lens system (QSO 0957+561) from 2000 February 3 to 2000 March 31. The daily VR observations were made with the 2.56-m Nordic Optical Telescope at Roque de los Muchachos Observatory, La Palma (Spain). We have derived detailed and robust VR light curves of the two components Q0957+561A and Q0957+561B. In spite of the excellent sampling rate, we have not found evidence in favor of true daily variability. With respect to variability on time-scales of several weeks, we measure VR gradients of about -0.8 mmag/day in Q0957+561A and + 0.3 mmag/day in Q0957+561B. The gradients are very probably originated in the far source, thus adopting this reasonable hypothesis (intrinsic variability), we compare them to the expected gradients during the evolution of a compact supernova remnant at the redshift of the source quasar. The starburst scenario is roughly consistent with some former events, but the new gradients do not seem to be caused by supernova remnant activity.

The redshift distribution of flat-spectrum radio sources with 5 GHz flux densities S>5 mJy is a key component in using current radio lens surveys to probe the cosmological model. We have constructed the first flat-spectrum radio sample in the flux density range 3-20 mJy. Our new sample has 33 sources; we have determined the redshifts of 14 of these (42% complete). The low mean redshift, <z>~0.75, of our faintest sample needs to be confirmed by further observations to improve the sample completeness. We also increased the redshift completeness of several surveys of brighter flat-spectrum sources. While the mean redshift, <z>~1.1 of flat-spectrum samples fainter than 1 Jy is nearly constant, the fraction of the sources identifiable as quasars steadily drops from ~80% to ~10% as the flux density of the sources decreases.

We investigate the evolution of mass-selected early-type field galaxies using a sample of 28 gravitational lenses spanning the redshift range 0 < z < 1. Based on the redshift-dependent intercept of the fundamental plane in the rest frame B band, we measure an evolution rate of d log (M/L)_B / dz = -0.56 +/- 0.04 (all errors are 1 sigma unless noted) if we directly compare to the local intercept measured from the Coma cluster. Re-fitting the local intercept helps minimize potential systematic errors, and yields an evolution rate of d log (M/L)_B / dz = -0.54 +/- 0.09. An evolution analysis of properly-corrected aperture mass-to-light ratios (defined by the lensed image separations) is closely related to the Faber-Jackson relation. In rest frame B band we find an evolution rate of d log (M/L)_B / dz = -0.41 +/- 0.21, a present-day characteristic magnitude of M_{*0} = -19.70 + 5 log h +/- 0.29 (assuming a characteristic velocity dispersion of sigma_{DM*} = 225 km/s), and a Faber-Jackson slope of gamma_{FJ} = 3.29 +/- 0.58. The measured evolution rates favor old stellar populations (mean formation redshift z_f > 1.8 at 2 sigma confidence for a Salpeter initial mass function and a flat Omega_m =0.3 cosmology) among early-type field galaxies, and argue against significant episodes of star formation at z < 1.

We study the influence of gravitational microlensing on the AGN Fe K-alpha line confirming that unexpected enhancements recently detected in the iron line of some AGNs can be produced by this effect. We use a ray tracing method to study the influence of microlensing in the emission coming from a compact accretion disc considering both geometries, Schwarzschild and Kerr. Thanks to the small dimensions of the region producing the AGN Fe K-alpha line, the Einstein Ring Radii associated to even very small compact objects have size comparable to the accretion disc hence producing noticeable changes in the line profiles. Asymmetrical enhancements contributing differently to the peaks or to the core of the line are produced by a microlens, off-centered with respect to the accretion disc. In the standard configuration of microlensing by a compact object in an intervening galaxy, we found that the effects on the iron line are two orders of magnitude larger than those expected in the optical or UV emission lines. In particular, microlensing can satisfactorily explain the excess in the iron line emission found very recently in two gravitational lens systems, H 1413+117 and MG J0414+0534. Exploring other physical {scenario} for microlensing, we found that compact objects (of the order of one Solar mass) which belong to {the bulge or the halo} of the host galaxy can also produce significant changes in the Fe K$_\alpha$ line profile of an AGN. However, the optical depth estimated for this type of microlensing is {very small, $\tau\sim 0.001$, even in a favorable case.

We study a recently observed gravitational microlensing peak in the V-band light curve of Q2237+0305A using a relatively simple, but highly consistent with the data (the best-fit reduced \chi^2 is very close to 1), physical model. The source quasar is assumed to be a Newtonian geometrically-thin and optically-thick accretion disk. The disk has an arbitrary orientation, and both blackbody and greybody emission spectra are considered. When the electron-photon scattering plays a role, the greybody spectrum will be a simplified version of the exact one. In our model the microlensing variability result from the source crossing a caustic straight line. The main goal is to estimate the black hole mass and the mass accretion rate in QSO 2237+0305 as well as to discuss the power and the weakness of the technique, some possible improvements, and the future prospects from multifrequency monitoring of new microlensing peaks. We also put into perspective the new methodology and the results on the central engine in QSO 2237+0305. From the fitted microlensing parameters and reasonable dynamical/cosmological constraints, it is concluded that QSO 2237+0305 harbours a central massive black hole: 10^7 M_Sun < M < 6 10^8 M_Sun. While the information about the central dark mass is very interesting, the mass accretion rate is not so well constrained. The typical values of the disk luminosity/Eddington luminosity ratio are in the (1-20)*\epsilon range, where \epsilon \leq 1 is the emissivity relative to a blackbody and the highest L/L_{Edd} ratio corresponds to the largest deflector motion. Therefore, in order to verify L/L_{Edd} \leq 1, a relatively small projected peculiar motion of the lens galaxy and a greybody emission seem to be favored.

From the peak of a gravitational microlensing high-magnification event in the A component of QSO 2237+0305, which was accurately monitored by the GLITP collaboration, we derived new information on the nature and size of the optical V-band and R-band sources in the far quasar. If the microlensing peak is caused by a microcaustic crossing, we firstly obtained that the standard accretion disk is a scenario more reliable/feasible than other usual axially symmetric models. Moreover, the standard scenario fits both the V-band and R-band observations with reduced chi-square values very close to one. Taking into account all these results, a standard accretion disk around a supermassive black hole is a good candidate to be the optical continuum main source in QSO 2237+0305. Secondly, using the standard source model and a robust upper limit on the transverse galactic velocity, we inferred that 90 per cent of the V-band and R-band luminosities are emitted from a region with radial size less than 1.2 10^{-2} pc (= 3.7 10^{16} cm, at 2*sigma confidence level).

We discuss the effects of microlensing on the broad emission lines (BELs) of QSOs in the light of recent determinations of the size of the broad line region (BLR) and its scaling with QSO luminosity. Microlensing by star-sized objects can produce significant amplifications in the BEL of some multiple-imaged QSOs, and could be very relevant for high-ionization lines. We have identified a group of ten gravitational lens systems (30% of the selected sample) in which microlensing could be observed. Using standard kinematic models for AGNs, we have studied the changes induced in the line profile by a microlens located at different positions with respect to the center of the BLR. We found that microlensing could produce important effects such as the relative enhancement of different parts of the line profile or the displacement of the peak of the line. The study of BEL profiles of different ionization in a microlensed QSO image could be an alternative method for probing the BLR structure and size.

We report an increase of more than 0.2 mag in the optical brightness of the leading image (A) of the gravitational lens Q0957+561, detected during the 09/2000 -- 06/2001 monitoring campaign (2001 observing season). The brightening is similar to or even greater than the largest change ever detected during the 20 years of monitoring of this system. We discuss two different provisional explanations to this event: intrinsic source variability or microlensing (either short timescale microlensing or cessation of the historical microlensing). An exhaustive photometric monitoring of Q0957+561 is needed until summer of 2002 and during 2003 to discriminate between these possibilities.

We present VR observations of QSO 2237+0305 conducted by the GLITP collaboration from 1999 October 1 to 2000 February 3. The observations were made with the 2.56 m Nordic Optical Telescope at Roque de los Muchachos Observatory, La Palma (Spain). The PSF fitting method and an adapted version of the ISIS subtraction method have been used to derive the VR light curves of the four components (A-D) of the quasar. The mean errors range in the intervals 0.01-0.04 mag (PSF fitting) and 0.01-0.02 mag (ISIS subtraction), with the faintest component (D) having the largest uncertainties. We address the relatively good agreement between the A-D light curves derived using different filters, photometric techniques, and telescopes. The new VR light curves of component A extend the time coverage of a high magnification microlensing peak, which was discovered by the OGLE team.

We determine the extinction curve in the z_l=0.83 lens galaxy of the gravitational lens SBS0909+532 from the wavelength dependence of the flux ratio between the lensed quasar images (z_s=1.38) from 3400 to 9200\AA. It is the first measurement of an extinction curve at a cosmological distance of comparable quality to those obtained within the Galaxy. The extinction curve has a strong 2175\AA feature, a noteworthy fact because it has been weak or non-existent in most estimates of extinction curves outside the Galaxy. The extinction curve is fitted well by a standard $R_V=2.1\pm0.9$ Galactic extinction curve. If we assume standard Galactic extinction laws, the estimated dust redshift of $z=0.88\pm0.02$ is in good agreement with the spectroscopic redshift of the lens galaxy. The widespread assumption that SMC extinction curves are more appropriate models for cosmological dust may be incorrect.

We report on a program to obtain HST observations of galaxy-mass gravitational lens systems at optical and infrared wavelengths. Here we discuss the properties of 10 two-image gravitational lens systems (Q0142-100=UM673, B0218+357, SBS0909+532, BRI0952-0115, LBQS1009-0252, Q1017-207=J03.13, B1030+074, HE1104-1805, Q1208+1011, and PKS1830-211). We grouped these 10 systems because they have limited lens model constraints and often show poor contrast between the images and the lens galaxy. Of the 10 lens galaxies, 7 are probably early-type galaxies, 2 are probably late-type galaxies (B0218+357 and PKS1830-211), and one was not detected (Q1208+1011). We detect the host galaxies of the z_s=4.50 lensed quasar in BRI0952-0115, the z_s=2.32 lensed quasar in HE1104-1805, and the unlensed z=1.63 quasar near LBQS1009-0252. We fit a set of four standard lens models to each lens that had sufficient constraints to compare isothermal dark matter and constant mass-to-light lens models, and to explore the effects of local tidal shears.

An observing campaign with 10 participating observatories has undertaken to monitor the optical brightness of the Q0957 gravitationally lensed quasar for 10 consecutive nights in January 2000. The resulting A image brightness curve has significant brightness fluctuations and makes a photometric prediction for the B image light curve for a second campaign planned for 12-21 March 2001. The ultimate purpose is to determine the gravitational lens time delay to a fraction of an hour, and to seek evidence for rapid microlensing.

Recent observations of galaxy luminosity profiles and dark matter simulations find luminosity and mass distributions characterized by central cusps rather than finite core radii. We introduce and implement a set of cusped ellipsoidal lens models which include limits similar to the Jaffe, Hernquist, \eta and NFW models and apply them to the gravitational lenses APM 08279+5255 and B 1933+503. A successful model of APM 08279+5255 with its central, odd image requires a very shallow cusp, $\gamma \ltorder 0.4$ where $\rho \propto r^{-\gamma}$ as $r\to 0$, which is similar to a core rather than the favored $1 \ltorder \gamma \ltorder 2$ cusps. B~1933+503, by contrast, is well modeled with a steep density cusp, $1.6 \ltorder \gamma \ltorder 2.0$.

Photometric optical data of QSO 0957+561 covering the period 1984-99 are analyzed to discern between the two values of the time delay (417 and 424 days) mostly accepted in the recent literature. The observations, performed by groups from three different institutions-Princeton University, Harvard-Smithsonian Center for Astrophysics, and Instituto de Astrofisica de Canarias-and including new unpublished 1998-9 data from the IAC80 Telescope, were obtained in five filters (V, R, I, g, and r). The different light curves have been divided into observational seasons and two restriction have been applied to better calculate the time delay: (i) points with a strange photometric behavior have been removed; and (ii) data sets without large gaps have been selected. Simulated data were generated to test several numerical methods intended to compute the time delay. The methods giving the best results-the discrete correlation function, delta-square, z-transformed discrete correlation function, and linear interpolation-were then applied to real data. A first analysis of the 23 different time delays derived from each technique shows that the time delay must be into the interval 420-424 days. From our statistical study, a most probable value of 422.6 +/- 0.6 days is inferred.

HST V and I-band observations show that the gravitational lens B1359+154 consists of six images of a single z_s=3.235 radio source and its star-forming host galaxy, produced by a compact group of galaxies at z_l = 1. VLBA observations at 1.7 GHz strongly support this conclusion, showing six compact cores with similar low-frequency radio spectra. B1359+154 is the first example of galaxy-scale gravitational lensing in which more than four images are observed of the same background source. The configuration is due to the unique lensing mass distribution: three primary lens galaxies lying on the vertices of a triangle separated by 0.7 arcsec (4/h kpc), inside the 1.7 arcsec diameter Einstein ring defined by the radio images. The gravitational potential has additional extrema within this triangle, creating a pair of central images that supplement the ``standard'' four-image geometry of the outer components. Simple mass models consisting of three lens galaxies constrained by HST and VLBA astrometry naturally reproduce the observed image positions but must be finely-tuned to fit the flux densities.

We present wide-field images of the quadruple gravitational lenses B1422+231 and MGJ0414+0534 obtained from global Very Long Baseline Interferometry (VLBI) observations at 8.4 GHz on 23 November 1997. We present also a lens model for MGJ0414+0534, which reproduces the core positions and flux densities of the VLBI images, combining a singular isothermal ellipsoid with external shear, and a singular isothermal sphere to represent, respectively, the main lens galaxy and its neighbor, a faint galaxy near one of the images.

We observed the quadruple gravitationally lensed image of MGJ0414+0534 on 23 November 1997 with a global VLBI array at 8.4 GHz. We report wide-field imaging results of its four components at submilliarcsecond resolution, displaying complex core-like and jet-like extended structures. A simple model combining a singular isothermal ellipsoid to represent the main lens galaxy, external shear, and a singular isothermal sphere to represent an additional, nearby object accounts well for the core positions and flux densities of the VLBI images. This model predicts delays between the different lensed images of several weeks.

We report the discovery of a new double-image quasar that was found during a search for gravitational lenses in the southern sky. Radio source PMN J1838-3427 is composed of two flat-spectrum components with separation 1", flux density ratio 14:1 and matching spectral indices, in VLA and VLBA images. Ground-based BRI images show the optical counterpart (total I=18.6) is also double with the same separation and position angle as the radio components. An HST/WFPC2 image reveals the lens galaxy. The optical flux ratio (27:1) is higher than the radio value probably due to differential extinction of the components by the lens galaxy. An optical spectrum of the bright component contains quasar emission lines at z=2.78 and several absorption features, including prominent Ly-alpha absorption. The lens galaxy redshift could not be measured but is estimated to be z=0.36 +/- 0.08. The image configuration is consistent with the simplest plausible models for the lens potential. The flat radio spectrum and observed variability of PMN J1838-3427 suggest the time delay between flux variations of the components is measurable, and could thus provide an independent measurement of H_0.

We report the discovery of a new two-image gravitational lens candidate. The system RX J0921+4529 contains two z_s=1.66 quasars separated by 6."93 with an H band magnitude difference of $\Delta m=1.39$. The HST NIC2 H band images reveal an H=18.2 spiral galaxy between the quasar images, which is probably a member of a z_l=0.32 X-ray cluster centered on the field. We detect an extended source near the fainter quasar image but not in the brighter image. If this extended source is the host galaxy of the fainter quasar, then the system is a binary quasar rather than a gravitational lens. VLA observations at 3.6 cm reveal emission from the lens galaxy at the flux level of 1 mJy and a marginal detection of the brighter quasar.

Most gravitational lenses are early-type galaxies in relatively low density environments -- a ``field'' rather than a ``cluster'' population. We show that field early-type galaxies with 0 < z < 1, as represented by the lens galaxies, lie on the same fundamental plane as those in rich clusters at similar redshifts. We then use the fundamental plane to measure the combined evolutionary and K-corrections for early-type galaxies in the V, I and H bands. Only for passively evolving stellar populations formed at z > 2 (H_0=65 km/s Mpc, Omega_0=0.3, Lambda_0=0.7) can the lens galaxies be matched to the local fundamental plane. The high formation epoch and the lack of significant differences between the field and cluster populations contradict many current models of the formation history of early-type galaxies. Lens galaxy colors and the fundamental plane provide good photometric redshift estimates with an empirical accuracy of -0.03 +/- 0.11 for the 17 lenses with known redshifts. A mass model dominated by dark matter is more consistent with the data than either an isotropic or radially anisotropic constant M/L mass model, and a radially anisotropic model is better than an isotropic model.

Infrared images of the Q 0957+561 gravitational lens obtained with the Hubble Space Telescope show two large (\sim 5 arcsec) lensed images of the z_s=1.41 quasar host galaxy. Parts of the host galaxy are doubly-imaged like the quasar, while other parts are quadruply-imaged. The distortions of the host galaxy offer the best probe yet of the global structure of the lensing potential, which is essential for determining the Hubble constant from the measured time delay. The distortions are inconsistent with the predictions of previously published lens models, which invalidates those models and their implications for H_0. New models show that the distortions finally break the long-standing degeneracy between the shape of the lens galaxy and the tidal shear contributed by the cluster containing the lens galaxy. The shape of the lens galaxy's mass distribution must be remarkably similar to the shape of its luminosity distribution, and most models that produce reasonable values for the Hubble constant roughly match the observed ellipticity gradient and isophote twist of the lens galaxy. Also, the cluster must be non-spherical and produce a relatively small tidal shear. Although there are still degeneracies in the lens models that lead to a 25% uncertainty in the derived value of the Hubble constant, there are also strong prospects for new observations to further improve the constraints and reduce the uncertainties.

We present H-band observations of gravitationally lensed QSO host galaxies obtained with NICMOS on HST as part of the CfA-Arizona-Gravitational-Lens-Survey (CASTLES). The detections are greatly facilitated by the lensing magnification in these systems; we find that most hosts of radio-quiet QSOs (RQQ) at z~2 are of modest luminosity (L<L_*). They are 2-5 times fainter than the hosts of radio-loud QSOs at the same epoch. Compared to low redshifts, RQQ hosts at z~2 also support higher nuclear luminosities at given stellar host mass. This suggests that the supermassive black holes at their centers grew faster at early epochs than the stellar body of their surrounding host galaxies.

Most gravitational lens galaxies are early-type galaxies in relatively low density environments. We show that they lie on the same fundamental plane as early-type galaxies in both local and distant rich clusters. Their surface brightness evolution requires a typical star formation epoch of z=2-3, almost indistinguishable from that of rich cluster galaxies at comparable redshifts. The restricted galaxy type range of the lenses means that photometric redshifts work well even with only 1-3 filter photometry. We make preliminary measurements of the mass and luminosity functions of the lens galaxies, and find they are consistent with the standard model used for deriving cosmological limits using lens statistics. As expected for a mass-weighted sample, they are more massive and more luminous than the overall early-type galaxy population.

We use observations from the CASTLES survey of gravitational lenses to study extinction in 23 lens galaxies with $0 < z_l < 1$. The median differential extinction between lensed images is $\Delta E(\bv) = 0.05$ mag, and the directly measured extinctions agree with the amount needed to explain the differences between the statistics of radio and (optical) quasar lens surveys. We also measure the first extinction laws outside the local universe, including an $R_V=7.2$ curve for a molecular cloud at $z_l=0.68$ and an $R_V=1.5$ curve for the dust in a redshift $z_l=0.96$ elliptical galaxy.

We describe a series of new applications of gravitational lenses as astrophysical and cosmological tools. Such applications are becoming possible thanks to advances in the quality and quantity of observations. CASTLES (CfA-Arizona-Space-Telescope-LEns-Survey) {cfa-www.harvard.edu/castles} is an ongoing project that exploits the sensitivity and resolution of the Hubble Space Telescope (HST) at optical and infrared wavelengths to study the sample of over 50 known gravitational lenses. The observational goal of CASTLES is a uniform sample of multi-band images of all known galaxy-mass lens systems, to derive precise photometry and astrometry for the lens galaxies, all the known images, and any source or lens components that might have escaped detection. With these measurements we are investigating: (1) the properties of dust and of dark matter in lens galaxies out to z\sim1; (2) the dark matter in lens galaxies and in their environments; (3) the evolution of lens galaxies; and (4) the cosmological model, for instance by refining constraints on the Hubble constant H_0.

The lens candidate RXJ 0921+4529 consists of two z_s=1.66 quasar separated by 6."93 with an H band magnitude difference of \Delta m=1.39. The lens appears to be a z_l=0.31 X-ray cluster, including a m_H=18.5 late-type galaxy lying between the quasar images. We detect an extended source overlapping the faint quasar but not the bright quasar. If this extended source is the host galaxy of the fainter quasar, then the system is a quasar binary rather than a gravitational lens.

There are $\sim 1$, 0.1 and 0.01 gravitationally lensed X-ray sources per square degree with soft X-ray fluxes exceeding $10^{-15}, 10^{-14}$ and $10^{-13} ergs/s cm^{-2}$ respectively. These sources will be detected serendipitously with the Chandra X-ray Observatory at a rate of 1--3 lenses per year of high resolution imaging. The low detection rate is due to the small area over which the HRC and ACIS cameras have the $<1\farcs5$ FWHM resolution necessary to find gravitational lenses produced by galaxies. Deep images of rich clusters at intermediate redshifts should yield one wide separation ($\Delta\theta \gtorder 5\farcs0$) multiply-imaged background X-ray source for every $\sim 10$, 30 and 300 clusters imaged to the same flux limits.

Deep HST/NICMOS H (F160W) band observations of the z=1.96 quasar pair Q 1634+267A,B reveal no signs of a lens galaxy to a 1 sigma threshold of approximately 22.5 mag. The minimum luminosity for a normal lens galaxy would be a 6L_* galaxy at z > 0.5, which is 650 times greater than our detection threshold. Our observation constrains the infrared mass-to-light ratio of any putative, early-type, lens galaxy to (M/L)_H > 690h_65 (1200h_65) for Omega_0=0.1 (1.0) and H_0=65h_65 km/s/Mpc. We would expect to detect a galaxy somewhere in the field because of the very strong Mg II absorption lines at z=1.1262 in the Q 1634+267 A spectrum, but the HST H-band, I-band (F785LP) and V-band (F555W) images require that any associated galaxy be very under-luminous less than 0.1 L^*_H (1.0 L^*_I) if it lies within less than 40 h^{-1} (100 h^{-1}) kpc from Q 1634+267 A,B. While the large image separation (3.86 arcsec) and the lack of a lens galaxy strongly favor interpreting Q 1634+267A,B as a binary quasar system, the spectral similarity remains a puzzle. We estimate that at most 0.06% of randomly selected quasar pairs would have spectra as similar to each other as the spectra of Q 1634+267 A and B. Moreover, spectral similarities observed for the 14 quasar pairs are significantly greater than would be expected for an equivalent sample of randomly selected field quasars. Depending on how strictly we define similarity, we estimate that only 0.01--3% of randomly drawn samples of 14 quasar pairs would have as many similar pairs as the observational sample.

We describe the goals of the CASTLES (CfA-Arizona-Space-Telescope-LEns-Survey) project including a sample of NICMOS images of gravitational lenses and a brief list of the preliminary findings.

We determine 37 differential extinctions in 23 gravitational lens galaxies over the range 0 < z_l < 1. Only 7 of the 23 systems have spectral differences consistent with no differential extinction. The median differential extinction for the optically-selected (radio-selected) subsample is E(B-V)=0.04 (0.06) mag. The extinction is patchy and shows no correlation with impact parameter. The median total extinction of the bluest images is E(B-V)=0.08 mag, although the total extinction distribution is dominated by the uncertainties in the intrinsic colors of quasars. The directly measured extinction distributions are consistent with the mean extinction estimated by comparing the statistics of quasar and radio lens surveys, thereby confirming the need for extinction corrections when using the statistics of lensed quasars to estimate the cosmological model. A disjoint subsample of two face-on, radio-selected spiral lenses shows both high differential and total extinctions, but standard dust-to-gas ratios combined with the observed molecular gas column densities overpredict the amount of extinction by factors of 2-5. For several systems we can estimate the extinction law, ranging from R_V=1.5+/-0.2 for a z_l=0.96 elliptical, to R_V=7.2+/-0.1 for a z_l=0.68 spiral. For the four radio lenses where we can construct non-parametric extinction curves we find no evidence for gray dust over the IR-UV wavelength range. The dust can be used to estimate lens redshifts with reasonable accuracy, although we sometimes find two degenerate redshift solutions.

We use simple comparisons of the optical and radio properties of the wide separation (3'' to 10'') quasar pairs to demonstrate that they are binary quasars rather than gravitational lenses. The most likely model is that all the pairs are binary quasars, with a one-sided 2--sigma (1--sigma) upper limit of 22% (8%) on the lens fraction. Simple models for the expected enhancement of quasar activity during galaxy mergers that are consistent with the enhancement observed at low redshift can explain the incidence, separations, redshifts, velocity differences, and radio properties of the binary quasar population. Only a modest fraction (< 5%) of all quasar activity need be associated with galaxy mergers to explain the binary quasars.

We obtained MMT spectra of the gravitational lens candidate UM425 to compare the redshifts and line profile structures of lens components A and B, which are separated by ~6.5". The CIV lambda 1550 A emission in both A and B exhibits Broad Absorption Line (BAL) structure, consistent with the earlier detection of BAL structure in OVI lambda 1033 A and NV lambda 1240 A that was found with the International Ultraviolet Explorer (IUE) in component A. Cross-correlation of the spectra of A and B using emission lines of CIV lambda 1550 A, HeII lambda 1640 A, NIII lambda 1750 A, CIII] lambda 1909 A and MgII lambda 2800 A reveals a difference in the redshifts of A & B. However, the detailed BAL line profile structure found in spectra of A & B are strikingly similar to one another, which suggests the system is lensed. The spectra of A & B also indicate significant dust extinction, which we base on the presence of the lambda 2175 A absorption feature in the rest frame of the QSO (z_{QSO}=1.47). This feature is commonly seen in galactic sources, but is not generally observed in QSO spectra. Our spectra show the presence of the lambda 2175 A absorption feature in spectra of both images associated with the gravitational lens UM425. Based upon the strong similarity of BAL line profile structure exhibited by UM425A&B, particularly the presence of the lambda 2175 dust absorption feature in spectra of both images, we conclude that UM425 is a gravitational lens.

We are conducting a redshift survey of 177 flat-spectrum radio sources in 3 samples covering the 5GHz flux ranges 50-100, 100-200 and 200-250 mJy. So far, we have measured 124 redshifts with completenesses of 80%, 68% and 58% for the bright, intermediate, and faint flux ranges. Using the newly determined redshift distribution we can derive cosmological limits from the statistics of the 6 gravitational lenses in the JVAS sample of 2500 flat-spectrum radio sources brighter than 200 mJy at 5GHz. For flat cosmological models with a cosmological constant, the limit using only radio data is Omega_0 > 0.27 at 2-sigma (0.47 < Omega_0 < 1.38 at 1-sigma). The limits are statistically consistent with those for lensed quasars, and the combined radio + optical sample requires Omega_0 > 0.38 at 2-sigma (0.64 < Omega_0 < 1.66 at 1-sigma) for our most conservative redshift completeness model and assuming that there are no quasar lenses produced by spiral galaxies. Our best fit model improves by approximately 1-sigma if extinction in the early-type galaxies makes the lensed quasars fainter by Delta m = 0.58 +- 0.45 mag, but we still find a limit of Omega_0 > 0.26 at 2-sigma in flat cosmologies. The increasing fraction of radio galaxies as compared to quasars at fainter radio fluxes (rising from ~10% at 1 Jy to ~50% at 0.1 Jy) explains why lensed optical emission is common for radio lenses and partly explains the red color of radio-selected lenses.