We compare the microlensing-based continuum emission region size measurements in a sample of 15 gravitationally lensed quasars with estimates of luminosity-based thin disk sizes to constrain the temperature profile of the quasar continuum accretion region. If we adopt the standard thin disk model, we find a significant discrepancy between sizes estimated using the luminosity and those measured by microlensing of $\log(r_{L}/r_{\mu})=-0.57\pm0.08\,\text{dex}$. If quasar continuum sources are simple, optically thick accretion disks with a generalized temperature profile $T(r) \propto r^{-\beta}$, the discrepancy between the microlensing measurements and the luminosity-based size estimates can be resolved by a temperature profile slope $0.37 < \beta < 0.56$ at $1\,\sigma$ confidence. This is shallower than the standard thin disk model ($\beta=0.75$) at $3\,\sigma$ significance. We consider alternate accretion disk models that could produce such a temperature profile and reproduce the empirical continuum size scaling with black hole mass, including disk winds or disks with non-blackbody atmospheres.

We use thirteen seasons of R-band photometry from the 1.2m Leonard Euler Swiss Telescope at La Silla to examine microlensing variability in the quadruply-imaged lensed quasar WFI 2026-4536. The lightcurves exhibit ${\sim}\,0.2\,\text{mag}$ of uncorrelated variability across all epochs and a prominent single feature of ${\sim}\,0.1\,\text{mag}$ within a single season. We analyze this variability to constrain the size of the quasar's accretion disk. Adopting a nominal inclination of 60$^\text{o}$, we find an accretion disk scale radius of $\log(r_s/\text{cm}) = 15.74^{+0.34}_{-0.29}$ at a rest-frame wavelength of $2043\,\unicode{xC5}$, and we estimate a black hole mass of $\log(M_{\text{BH}}/M_{\odot}) = 9.18^{+0.39}_{-0.34}$, based on the CIV line in VLT spectra. This size measurement is fully consistent with the Quasar Accretion Disk Size - Black Hole Mass relation, providing another system in which the accretion disk is larger than predicted by thin disk theory.

In this work, we calculated the sizes of unresolved X-ray emission regions in three gravitationally lensed radio-loud quasars, B\,1422+231, MG\,J0414+0534 and Q\,0957+561, using a combination of imaging and spectral analysis on the X-ray data taken from the \textit{Chandra X-Ray Observatory}. We tentatively detected FeK$\alpha$ emission lines in MG\,J0414+0534 and Q\,0957+561 with over 95\% significance, whereas, we did not significantly detect FeK$\alpha$ emission in B\,1422+231. We constructed differential microlensing light curves from absorption corrected count rates. We subsequently performed a microlensing analysis on the X-ray microlensing light curves to measure the X-ray source sizes in soft (0.83--3.6 keV), hard (3.6--21.8 keV), and full (0.83--21.8 keV) bands, based on either Bayesian or maximum likelihood probabilities. For B\,1422+231, sizes from the two methods are consistent with each other, e.g. $R_X^{hard}/R_G = 6.17\pm5.48 \text{ (Bayesian), } 11.81\pm3.75 \text{ (maximum likelihood)}$, where $R_G=GM_{BH}/c^2)$. However, for MG\,J0414+0534 and Q\,0957+561, the two methods yield completely different results suggesting that more frequently sampled data with better signal-to-noise ratio are needed to measure the source size for these two objects. Comparing the acquired size values with the radio-quiet sample in the literature we found that our results are consistent with X-ray source size scaling approximately as $R_X \propto M_{BH}$ with the mass of the central supermassive black hole. Our results also indicate that radio-loud quasars tend to have larger unresolved X-ray emission sizes compared to the radio-quiet ones.

We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.

We present an analysis of Chandra spectra of five gravitationally lensed active galactic nuclei. We confirm the previous detections of FeK$\alpha$ emission lines in most images of these objects with high significance. The line energies range from 5.8 to 6.8 keV with widths from unresolved to 0.6 keV, consistent with emission close to spinning black holes viewed at different inclination angles. We also confirm the positive offset from the Iwasawa-Taniguchi effect, the inverse correlation between the FeK$\alpha$ equivalent width and the X-ray luminosity in AGN, where our measured equivalent widths are larger in lensed quasars. We attribute this effect to microlensing, and perform a microlensing likelihood analysis to constrain the emission size of the relativistic reflection region and the spin of supermassive black holes, assuming that the X-ray corona and the reflection region, responsible for the iron emission line, both follow power-law emissivity profiles. The microlensing analysis yields strong constraints on the spin and emissivity index of the reflection component for Q2237+0305, with $a > 0.92$ and $n > 5.4$. For the remaining four targets, we jointly constrain the two parameters, yielding $a=0.8\pm0.16$ and an emissivity index of $n=4.0\pm 0.8$, suggesting that the relativistic X-ray reflection region is ultra-compact and very close to the innermost stable circular orbits of black holes, which are spinning at close to the maximal value. We successfully constrain the half light radius of the emission region to $< 2.4$ $r_g$ ($r_g = GM/c^2$) for Q2237+0305 and in the range 5.9--7.4 $r_g$ for the joint sample.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$

We present 13 seasons of $R$-band photometry of the quadruply-lensed quasar WFI 2033-4723 from the 1.3m SMARTS telescope at CTIO and the 1.2m Euler Swiss Telescope at La Silla, in which we detect microlensing variability of $\sim0.2$ mags on a timescale of $\sim$6 years. Using a Bayesian Monte Carlo technique, we analyze the microlensing signal to obtain a measurement of the size of this system's accretion disk of $\log (r_s/{\rm cm}) = 15.86^{+0.25}_{-0.27}$ at $\lambda_{rest} = 2481{\rm \AA}$, assuming a $60^\circ$ inclination angle. We confirm previous measurements of the BC and AB time delays, and we obtain a tentative measurement of the delay between the closely spaced A1 and A2 images of $\Delta t_{A1A2} = t_{A1} - t_{A2} = -3.9^{+3.4}_{-2.2}$ days. We conclude with an update to the Quasar Accretion Disk Size - Black Hole Mass Relation, in which we confirm that the accretion disk size predictions from simple thin disk theory are too small.

The requirements-driven OSIRIS-REx Camera Suite (OCAMS) acquires images essential to collecting a sample from the surface of Bennu. During proximity operations, these images document the presence of satellites and plumes, record spin state, enable an accurate digital terrain model of the shape of the asteroid and identify any surface hazards. They confirm the presence of sampleable regolith on the surface, observe the sampling event itself, and image the sample head in order to verify its readiness to be stowed. They document the history of Bennu as an example of early solar system material, as a microgravity body with a planetesimal size-scale, and as a carbonaceous object. OCAMS is fitted with three cameras. The MapCam records point-source color images on approach to the asteroid in order to connect ground-based point-source observations of Bennu to later higher-resolution surface spectral imaging. The SamCam documents the sample site before, during, and after it is disturbed by the sample mechanism. The PolyCam, using its focus mechanism, observes the sample site at sub-centimeter resolutions, revealing surface texture and morphology. While their imaging requirements divide naturally between the three cameras, they preserve a strong degree of functional overlap. OCAMS and the other spacecraft instruments allow the OSIRIS-REx mission to collect a sample from a microgravity body on the same visit during which it was first optically acquired from long range, a useful capability as humanity explores near-Earth, Main-Belt and Jupiter Trojan asteroids.

We present a promising new technique, the g-distribution method, for measuring the inclination angle (i), the innermost stable circular orbit (ISCO), and the spin of a supermassive black hole. The g-distribution method uses measurements of the energy shifts in the relativistic iron line emitted by the accretion disk of a supermassive black hole due to microlensing by stars in a foreground galaxy relative to the g-distribution shifts predicted from microlensing caustic calculations. We apply the method to the gravitationally lensed quasars RX J1131-1231 (z_s=0.658, z_l=0.295), QJ 0158-4325 (z_s=1.294, z_l=0.317), and SDSS 1004+4112 (z_s=1.734, z_l=0.68). For RX J1131-1231 our initial results indicate that r_ISCO<8.5 gravitational radii (r_g) and i > 76 degrees. We detect two shifted Fe lines, in several observations, as predicted in our numerical simulations of caustic crossings. The current DeltaE-distribution of RX J1131-1231 is sparsely sampled but further X-ray monitoring of RX J1131-1231 and other lensed quasars will provide improved constraints on the inclination angles, ISCO radii and spins of the black holes of distant quasars.

We present an X-ray photometric analysis of six gravitationally lensed quasars, with observation campaigns spanning from 5 to 14 years, measuring the total (0.83 - 21.8 keV restframe), soft (0.83 - 3.6 keV), and hard (3.6 - 21.8 keV) band image flux ratios or each epoch. Using the ratios of the model-predicted macro-magnifications as baselines, we build differential microlensing light curves and obtain joint likelihood functions for the average X-ray emission region sizes. Our analysis yields a Probability Distribution Function for the average half-light radius of the X-Ray emission region in the sample that peaks slightly above 1 gravitational radius and with nearly indistinguishable 68% confidence (one-sided) upper limits of 17.8 and 18.9 gravitational radii for the soft and hard X-ray emitting regions, assuming a mean stellar mass of 0.3 solar masses. We see hints of energy dependent microlensing between the soft and hard bands in two of the objects. In a separate analysis on the root-mean-square of the microlensing variability, we find significant differences between the soft and hard bands but the sign of the difference is not consistent across the sample. This suggests the existence of some kind of spatial structure to the X-ray emission in an otherwise extremely compact source. We also discover a correlation between the RMS microlensing variability and the average microlensing amplitude.

We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing analysis are significantly larger than those predicted by the theoretical-thin-disk estimate. In a few cases we also constrain the slope of the size-wavelength relation. Our size constraints of the soft and hard X-ray emission regions of quasars indicate that in some objects of our sample the hard X-ray emission region is more compact than the soft and in others the soft emission region is smaller. This difference may be the result of the relative strengths of the disk-reflected (harder and extended) versus corona-direct (softer and compact) components in the quasars of our sample. Finally, we present the analysis of several strong microlensing events where we detect an evolution of the relativistic Fe line profile as the magnification caustic traverses the accretion disk. These caustic crossings are used to provide constraints on the innermost stable circular orbit (ISCO) radius and the accretion disk inclination angle of the black hole in quasar RX J1131-1231.

We analyze the optical, UV, and X-ray microlensing variability of the lensed quasar SDSS J0924+0219 using six epochs of Chandra data in two energy bands (spanning 0.4-8.0 keV, or 1-20 keV in the quasar rest frame), 10 epochs of F275W (rest-frame 1089A) Hubble Space Telescope data, and high-cadence R-band (rest-frame 2770A) monitoring spanning eleven years. Our joint analysis provides robust constraints on the extent of the X-ray continuum emission region and the projected area of the accretion disk. The best-fit half-light radius of the soft X-ray continuum emission region is between 5x10^13 and 10^15 cm, and we find an upper limit of 10^15 cm for the hard X-rays. The best-fit soft-band size is about 13 times smaller than the optical size, and roughly 7 GM_BH/c^2 for a 2.8x10^8 M_sol black hole, similar to the results for other systems. We find that the UV emitting region falls in between the optical and X-ray emitting regions at 10^14 cm < r_1/2,UV < 3x10^15 cm. Finally, the optical size is significantly larger, by 1.5*sigma, than the theoretical thin-disk estimate based on the observed, magnification-corrected I-band flux, suggesting a shallower temperature profile than expected for a standard disk.

We present three complete seasons and two half-seasons of SDSS r-band photometry of the gravitationally lensed quasar SBS 0909+532 from the U.S. Naval Observatory, as well as two seasons each of SDSS g-band and r-band monitoring from the Liverpool Robotic Telescope. Using Monte Carlo simulations to simultaneously measure the system's time delay and model the r-band microlensing variability, we confirm and significantly refine the precision of the system's time delay to \Delta t_{AB} = 50^{+2}_{-4} days, where the stated uncertainties represent the bounds of the formal 1\sigma\ confidence interval. There may be a conflict between the time delay measurement and a lens consisting of a single galaxy. While models based on the Hubble Space Telescope astrometry and a relatively compact stellar distribution can reproduce the observed delay, the models have somewhat less dark matter than we would typically expect. We also carry out a joint analysis of the microlensing variability in the r- and g-bands to constrain the size of the quasar's continuum source at these wavelengths, obtaining log[(r_{s,r}/cm) [cos{i}/0.5]^{1/2}] = 15.3 \pm 0.3 and log[(r_{s,g}/cm) [cos{i}/0.5]^{1/2}] = 14.8 \pm 0.9, respectively. Our current results do not formally constrain the temperature profile of the accretion disk but are consistent with the expectations of standard thin disk theory.

We present four new seasons of optical monitoring data and six epochs of X-ray photometry for the doubly-imaged lensed quasar Q J0158-4325. The high-amplitude, short-period microlensing variability for which this system is known has historically precluded a time delay measurement by conventional methods. We attempt to circumvent this limitation by application of a Monte Carlo microlensing analysis technique, but we are only able to prove that the delay must have the expected sign (image A leads image B). Despite our failure to robustly measure the time delay, we successfully model the microlensing at optical and X-ray wavelengths to find a half light radius for soft X-ray emission log(r_{1/2,X,soft}/cm) = 14.3^{+0.4}_{-0.5}, an upper limit on the half-light radius for hard X-ray emission log(r_{1/2,X,hard}/cm) <= 14.6 and a refined estimate of the inclination-corrected scale radius of the optical R-band (rest frame 3100 Angstrom) continuum emission region of log(r_s/cm) = 15.6+-0.3.

We present \emph{Chandra} monitoring data for six gravitationally lensed quasars: QJ 0158$-$4325, HE 0435$-$1223, HE 1104$-$1805, SDSS 0924+0219, SDSS 1004+4112, and Q 2237+0305. We detect X-ray microlensing variability in all six lenses with high confidence. We detect energy dependent microlensing in HE 0435$-$1223, SDSS 1004+4112, SDSS 0924+0219 and Q 2237+0305. We present a detailed spectral analysis for each lens, and find that simple power-law models plus Gaussian emission lines give good fits to the spectra. We detect intrinsic spectral variability in two epochs of Q 2237+0305. We detect differential absorption between images in four lenses. We also detect the \feka\ emission line in all six lenses, and the Ni XXVII K$\alpha$ line in two images of Q 2237+0305. The rest frame equivalent widths of the \feka\ lines are measured to be 0.4--1.2 keV, significantly higher than those measured in typical active galactic nuclei of similar X-ray luminosities. This suggests that the \feka\ emission region is more compact or centrally concentrated than the continuum emission region.

We present evidence for ultraviolet/optical microlensing in the gravitationally lensed quasar Q0957+561. We combine new measurements from our optical monitoring campaign at the United States Naval Observatory, Flagstaff (USNO) with measurements from the literature and find that the time-delay-corrected r-band flux ratio m_A - m_B has increased by ~0.1 magnitudes over a period of five years beginning in the fall of 2005. We apply our Monte Carlo microlensing analysis procedure to the composite light curves, obtaining a measurement of the optical accretion disk size, log {(r_s/cm)[cos(i)/0.5]^{1/2}} = 16.2^{+0.5}_{-0.6}, that is consistent with the quasar accretion disk size - black hole mass relation.

In this paper, we use cubic polynomials to approximate RR Lyrae light curves and apply the method to HST data of RR Lyraes in the halo of M31. We compare our method to the standard method of Fourier decomposition and find that the method of cubic polynomials eliminates virtually all ringing effects and does so with significantly fewer parameters than the Fourier technique. Further, for RRc stars the parameters in the fit are all physical. Our study also reveals a number of additional periodicites in this data not found previously: we find 23 RRc stars, 29 RRab stars and 3 multiperiodic stars.

We use the microlensing variability observed for eleven gravitationally lensed quasars to show that the accretion disk size at a rest-frame wavelength of 2500 Angstroms is related to the black hole mass by log(R_{2500}/cm)=(15.78\pm0.12) + (0.80\pm0.17)\log(M_BH/10^9M_sun). This scaling is consistent with the expectation from thin disk theory (R ~ M_BH^{2/3}), but when interpreted in terms of the standard thin disk model (T ~ R^{-3/4}), it implies that black holes radiate with very low efficiency, log(eta) = -1.77\pm0.29 + log(L/L_E) where eta=L/(Mdot*c^2). Only by making the maximum reasonable shifts in the average inclination, Eddington factors and black hole masses can we raise the efficiency estimate to be marginally consistent with typical efficiency estimates (eta ~ 10%). With one exception, these sizes are larger by a factor of ~4 than the size needed to produce the observed 0.8 micron quasar flux by thermal radiation from a thin disk with the same T ~ R^{-3/4} temperature profile. While scattering a significant fraction of the disk emission on large scales or including a large fraction of contaminating line emission can reduce the size discrepancy, resolving it also appears to require that accretion disks have flatter temperature/surface brightness profiles.

We use gravitational microlensing of the four images of the z=0.658 quasar RXJ1131-1231 to measure the sizes of the optical and X-ray emission regions of the quasar. The (face-on) scale length of the optical disk at rest frame 400 nm is 1.3 10^15cm, while the half-light radius of the rest frame 0.3-17 keV X-ray emission is 2.3 10^14cm. The formal uncertainties are factors of 1.6 and 2.0, respectively. With the exception of the lower limit on the X-ray size, the results are very stable against any changes in the priors used in the analysis. Based on the Hbeta line-width, we estimate that the black hole mass is ~10^8 Msun, which corresponds to a gravitational radius of r_g~2 10^13 cm. Thus, the X-ray emission is emerging on scales of ~10r_g and the 400 nm emission on scales of ~70 r_g. A standard thin disk of this size should be significantly brighter than observed. Possible solutions are to have a flatter temperature profile or to scatter a large fraction of the optical flux on larger scales after it is emitted. While our calculations were not optimized to constrain the dark matter fraction in the lens galaxy, dark matter dominated models are favored. With well-sampled optical and X-ray light curves over a broad range of frequencies there will be no difficulty in extending our analysis to completely map the structure of the accretion disk as a function of wavelength.

Gravitationally lensed quasars can be used to map the mass distribution in lensing galaxies and to estimate the Hubble constant H0 by measuring the time delays between the quasar images. Here we report the measurement of two independent time delays in the quadruply imaged quasar WFI J2033-4723 (z = 1.66). Our data consist of R-band images obtained with the Swiss 1.2 m EULER telescope located at La Silla and with the 1.3 m SMARTS telescope located at Cerro Tololo. The light curves have 218 independent epochs spanning 3 full years of monitoring between March 2004 and May 2007, with a mean temporal sampling of one observation every 4th day. We measure the time delays using three different techniques, and we obtain Dt(B-A) = 35.5 +- 1.4 days (3.8%) and Dt(B-C) = 62.6 +4.1/-2.3 days (+6.5%/-3.7%), where A is a composite of the close, merging image pair. After correcting for the time delays, we find R-band flux ratios of F_A/F_B = 2.88 +- 0.04, F_A/F_C = 3.38 +- 0.06, and F_A1/F_A2 = 1.37 +- 0.05 with no evidence for microlensing variability over a time scale of three years. However, these flux ratios do not agree with those measured in the quasar emission lines, suggesting that longer term microlensing is present. Our estimate of H0 agrees with the concordance value: non-parametric modeling of the lensing galaxy predicts H0 = 67 +13/-10 km s-1 Mpc-1, while the Single Isothermal Sphere model yields H0 = 63 +7/-3 km s-1 Mpc-1 (68% confidence level). More complex lens models using a composite de Vaucouleurs plus NFW galaxy mass profile show twisting of the mass isocontours in the lensing galaxy, as do the non-parametric models. As all models also require a significant external shear, this suggests that the lens is a member of the group of galaxies seen in field of view of WFI J2033-4723.

We analyzed the microlensing of the X-ray and optical emission of the lensed quasar PG 1115+080. We find that the effective radius of the X-ray emission is 1.3(+1.1 -0.5) dex smaller than that of the optical emission. Viewed as a thin disk observed at inclination angle i, the optical accretion disk has a scale length, defined by the point where the disk temperature matches the rest frame energy of the monitoring band (kT=hc/lambda_rest with lambda_rest=0.3 micron), of log[(r_{s,opt}/cm)(cos(i) / 0.5)^{1/2}] = 16.6 \pm 0.4. The X-ray emission region (1.4-21.8 keV in the rest frame) has an effective half-light radius of log[r_{1/2,X}/cm] = 15.6 (+0.6-0.9}. Given an estimated black hole mass of 1.2 * 10^9 M_sun, corresponding to a gravitational radius of log[r_g/cm] = 14.3, the X-ray emission is generated near the inner edge of the disk while the optical emission comes from scales slightly larger than those expected for an Eddington-limited thin disk. We find a weak trend supporting models with low stellar mass fractions near the lensed images, in mild contradiction to inferences from the stellar velocity dispersion and the time delays.

We report the serendipitous detection of a very bright, very nearby microlensing event. In late October 2006, an otherwise unremarkable A0 star at a distance ~1 kpc (GSC 3656-1328) brightened achromatically by a factor of nearly 40 over the span of several days and then decayed in an apparently symmetrical way. We present a light curve of the event based on optical photometry from the Center for Backyard Astrophysics and the All Sky Automated Survey, as well as near-infrared photometry from the Peters Automated Infrared Imaging Telescope. This light curve is well-fit by a generic microlensing model. We also report optical spectra, and Swift X-ray and UV observations that are consistent with the microlensing interpretation. We discuss and reject alternative explanations for this variability. The lens star is probably a low-mass star or brown dwarf, with a relatively high proper motion of >20 mas/yr, and may be visible using precise optical/infrared imaging taken several years from now. A modest, all-sky survey telescope could detect ~10 such events per year, which would enable searches for very low-mass planetary companions to relatively nearby stars.

This paper presents spectroscopy of supernovae discovered in the first season of the Sloan Digital Sky Survey-II Supernova Survey. This program searches for and measures multi-band light curves of supernovae in the redshift range z = 0.05 - 0.4, complementing existing surveys at lower and higher redshifts. Our goal is to better characterize the supernova population, with a particular focus on SNe Ia, improving their utility as cosmological distance indicators and as probes of dark energy. Our supernova spectroscopy program features rapid-response observations using telescopes of a range of apertures, and provides confirmation of the supernova and host-galaxy types as well as precise redshifts. We describe here the target identification and prioritization, data reduction, redshift measurement, and classification of 129 SNe Ia, 16 spectroscopically probable SNe Ia, 7 SNe Ib/c, and 11 SNe II from the first season. We also describe our efforts to measure and remove the substantial host galaxy contamination existing in the majority of our SN spectra.

We expand our Bayesian Monte Carlo method for analyzing the light curves of gravitationally lensed quasars to simultaneously estimate time delays and quasar structure including their mutual uncertainties. We apply the method to HE1104-1805 and QJ0158-4325, two doubly-imaged quasars with microlensing and intrinsic variability on comparable time scales. For HE1104-1805 the resulting time delay of (Delta t_AB) = t_A - t_B = 162.2 -5.9/+6.3 days and accretion disk size estimate of log(r_s/cm) = 15.7 -0.5/+0.4 at 0.2 micron in the rest frame are consistent with earlier estimates but suggest that existing methods for estimating time delays in the presence of microlensing underestimate the uncertainties. We are unable to measure a time delay for QJ0158-4325, but the accretion disk size is log(r_s/cm) = 14.9 +/- 0.3 at 0.3 micron in the rest frame.

We present 107 new epochs of optical monitoring data for the four brightest images of the gravitational lens SDSS J1004+4112 observed between October 2006 and June 2007. Combining this data with the previously obtained light curves, we determine the time delays between images A, B and C. We confirm our previous measurement finding that A leads B by dt_BA=40.6+-1.8 days, and find that image C leads image A by dt_CA=821.6+-2.1 days. The lower limit on the remaining delay is that image D lags image A by dt_AD>1250 days. Based on the microlensing of images A and B we estimate that the accretion disk size at a rest wavelength of 2300 angstrom is 10^{14.8+-0.3} cm for a disk inclination of cos{i}=1/2, which is consistent with the microlensing disk size-black hole mass correlation function given our estimate of the black hole mass from the MgII line width of logM_BH/M_sun=8.44+-0.14. The long delays allow us to fill in the seasonal gaps and assemble a continuous, densely sampled light curve spanning 5.7 years whose variability implies a structure function with a logarithmic slope of gamma = 0.35+-0.02. As C is the leading image, sharp features in the C light curve can be intensively studied 2.3 years later in the A/B pair, potentially allowing detailed reverberation mapping studies of a quasar at minimal cost.

The Sloan Digital Sky Survey-II (SDSS-II) has embarked on a multi-year project to identify and measure light curves for intermediate-redshift (0.05 < z < 0.35) Type Ia supernovae (SNe Ia) using repeated five-band (ugriz) imaging over an area of 300 sq. deg. The survey region is a stripe 2.5 degrees wide centered on the celestial equator in the Southern Galactic Cap that has been imaged numerous times in earlier years, enabling construction of a deep reference image for discovery of new objects. Supernova imaging observations are being acquired between 1 September and 30 November of 2005-7. During the first two seasons, each region was imaged on average every five nights. Spectroscopic follow-up observations to determine supernova type and redshift are carried out on a large number of telescopes. In its first two three-month seasons, the survey has discovered and measured light curves for 327 spectroscopically confirmed SNe Ia, 30 probable SNe Ia, 14 confirmed SNe Ib/c, 32 confirmed SNe II, plus a large number of photometrically identified SNe Ia, 94 of which have host-galaxy spectra taken so far. This paper provides an overview of the project and briefly describes the observations completed during the first two seasons of operation.

We analyze V, I and H band HST images and two seasons of R-band monitoring data for the gravitationally lensed quasar SDSS0924+0219. We clearly see that image D is a point-source image of the quasar at the center of its host galaxy. We can easily track the host galaxy of the quasar close to image D because microlensing has provided a natural coronograph that suppresses the flux of the quasar image by roughly an order of magnitude. We observe low amplitude, uncorrelated variability between the four quasar images due to microlensing, but no correlated variations that could be used to measure a time delay. Monte Carlo models of the microlensing variability provide estimates of the mean stellar mass in the lens galaxy (0.02 Msun < M < 1.0 Msun), the accretion disk size (the disk temperature is 5 x 10^4 K at 3.0 x 10^14 cm < rs < 1.4 x 10^15 cm), and the black hole mass (2.0 x 10^7 Msun < MBH \eta_{0.1}^{-1/2} (L/LE)^{1/2} < 3.3 x 10^8 Msun), all at 68% confidence. The black hole mass estimate based on microlensing is consistent with an estimate of MBH = 7.3 +- 2.4 x 10^7 Msun from the MgII emission line width. If we extrapolate the best-fitting light curve models into the future, we expect the the flux of images A and B to remain relatively stable and images C and D to brighten. In particular, we estimate that image D has a roughly 12% probability of brightening by a factor of two during the next year and a 45% probability of brightening by an order of magnitude over the next decade.

We use the microlensing variability observed for nine gravitationally lensed quasars to show that the accretion disk size at 2500 Angstroms is related to the black hole mass by log(R_2500/cm) = (15.6+-0.2) + (0.54+-0.28)log(M_BH/10^9M_sun). This scaling is consistent with the expectation from thin disk theory (R ~ M_BH^(2/3)), but it implies that black holes radiate with relatively low efficiency, log(eta) = -1.29+-0.44 + log(L/L_E) where eta=L/(Mdot c^2). These sizes are also larger, by a factor of ~3, than the size needed to produce the observed 0.8 micron quasar flux by thermal radiation from a thin disk with the same T ~ R^(-3/4) temperature profile. More sophisticated disk models are clearly required, particularly as our continuing observations improve the precision of the measurements and yield estimates of the scaling with wavelength and accretion rate.

Recent observations support the idea that nuclear black holes grew by gas accretion while shining as luminous quasars at high redshift, and they establish a relation of the black hole mass with the host galaxy's spheroidal stellar system. We develop an analytic model to calculate the expected impact of mergers on the masses of black holes in massive clusters of galaxies. We use the extended Press-Schechter formalism to generate Monte Carlo merger histories of halos with a mass 10^{15} h^{-1} Msun. We assume that the black hole mass function at z=2 is similar to that inferred from observations at z=0 (since quasar activity declines markedly at z<2), and we assign black holes to the progenitor halos assuming a monotonic relation between halo mass and black hole mass. We follow the dynamical evolution of subhalos within larger halos, allowing for tidal stripping, the loss of orbital energy by dynamical friction, and random orbital perturbations in gravitational encounters with subhalos, and we assume that mergers of subhalos are followed by mergers of their central black holes. Our analytic model reproduces numerical estimates of the subhalo mass function. We find that the most massive black holes in massive clusters typically grow by a factor ~ 2 by mergers after gas accretion has stopped. In our ten realizations of 10^{15} h^{-1} Msun clusters, the highest initial (z=2) black hole masses are 5-7 x 10^9 Msun, but four of the clusters contain black holes in the range 1-1.5 x 10^{10} Msun at z=0. Satellite galaxies may host black holes whose mass is comparable to, or even greater than, that of the central galaxy. Thus, black hole mergers can significantly extend the very high end of the black hole mass function.

Microlensing of gravitationally lensed quasars by the stars in the foreground lens galaxy can be used to probe the nature of dark matter, to determine the mean stellar mass in the lens galaxy, and to measure the internal structure of quasar accretion disks. Until recently, little progress has been made toward using microlensing for these purposes because of the difficulty in obtaining the necessary data and the lack of good analysis methods. In the last few years, both problems have been solved. In particular, Bayesian analysis methods provide a general approach to measuring quantities of physical interest and their uncertainties from microlensing light curves. We discuss the data and the analysis methods and show preliminary results for all three astrophysical applications.

We present 426 epochs of optical monitoring data spanning 1000 days from December 2003 to June 2006 for the gravitationally lensed quasar SDSS J1004+4112. The time delay between the A and B images is 38.4+/-2.0 days in the expected sense that B leads A and the overall time ordering is C-B-A-D-E. The measured delay invalidates all published models. The models failed because they neglected the perturbations from cluster member galaxies. Models including the galaxies can fit the data well, but strong conclusions about the cluster mass distribution should await the measurement of the longer, and less substructure sensitive, delays of the C and D images. For these images, a CB delay of 681+/-15 days is plausible but requires confirmation, while CB and AD delays of >560 days and > 800 days are required. We clearly detect microlensing of the A/B images, with the delay-corrected flux ratios changing from B-A=0.44+/-0.01 mag in the first season to 0.29+/-0.01 mag in the second season and 0.32+/-0.01 mag in the third season.

We present RETROCAM, an auxiliary CCD camera that can be rapidly inserted into the optical beam of the MDM 2.4m telescope. The speed and ease of reconfiguring the telescope to use the imager and a straightforward user interface permit the camera to be used during the course of other observing programs. This in turn encourages RETROCAM's use for a variety of monitoring projects.

In searches for planetary transits in the field, well over half of the survey stars are typically giants or other stars that are too large to permit straightforward detection of planets. For all-sky searches of bright V<~11 stars, the fraction is ~90%. We show that the great majority of these contaminants can be removed from the sample by analyzing their reduced proper motions (RPMs): giants have much lower RPMs than dwarfs of the same color. We use Hipparcos data to design a RPM selection function that eliminates most evolved stars, while rejecting only 9% of viable transit targets. Our method can be applied using existing or soon-to-be-released all-sky data to stars V<12.5 in the northern hemisphere and V<12 in the south. The method degrades at fainter magnitudes, but does so gracefully. For example, at V=14 it can still be used to eliminate giants redward of V-I~0.95, that is, the blue edge of the red giant clump.