We present an extensive search for RR Lyrae stars in and around the ultra-faint Milky Way companions Segue 2 and Segue 3. The former (M_V = -2.5, Belokurov et al. 2009) appears to be an extremely faint dwarf galaxy companion of the Milky Way. The latter (M_V = 0.0, Fadely et al. 2011) is among the faintest star clusters known. We use B and V band time-series imaging obtained at the WIYN 0.9 meter telescope at Kitt Peak National Observatory to search for RR Lyrae in these objects. In our Segue 2 observations, we present a previously unknown fundamental mode (RRab) RR Lyrae star with a period of P_ab = 0.748 days. With this measurement, we revisit the inverse correlation between <P_ab> and <[Fe/H]> established in the literature for Milky Way dwarf galaxies and their RR Lyrae. In this context, the long period of Segue 2's RRab star as well as the known significant spread in metallicity in this dwarf galaxy are consistent with the observed trend in <P_ab> and <[Fe/H]>. We derive the first robust distance to Segue 2, using both its RRab star and spectroscopically confirmed blue horizontal branch stars. Using [Fe/H] = -2.16 and -2.44 dex, we find d_RRL = 36.6 +2.5/-2.4 and 37.7 +/- 2.7 kpc; assuming [Fe/H] = -2.257 dex, we find d_BHB = 34.4 +/- 2.6 kpc. Although no RR Lyrae were present in the Segue 3 field, we found a candidate eclipsing binary star system.

Ground-based optical surveys such as PanSTARRS, DES, and LSST, will produce large catalogs to limiting magnitudes of r > 24. Star-galaxy separation poses a major challenge to such surveys because galaxies---even very compact galaxies---outnumber halo stars at these depths. We investigate photometric classification techniques on stars and galaxies with intrinsic FWHM < 0.2 arcsec. We consider unsupervised spectral energy distribution template fitting and supervised, data-driven Support Vector Machines (SVM). For template fitting, we use a Maximum Likelihood (ML) method and a new Hierarchical Bayesian (HB) method, which learns the prior distribution of template probabilities from the data. SVM requires training data to classify unknown sources; ML and HB don't. We consider i.) a best-case scenario (SVM_best) where the training data is (unrealistically) a random sampling of the data in both signal-to-noise and demographics, and ii.) a more realistic scenario where training is done on higher signal-to-noise data (SVM_real) at brighter apparent magnitudes. Testing with COSMOS ugriz data we find that HB outperforms ML, delivering ~80% completeness, with purity of ~60-90% for both stars and galaxies, respectively. We find no algorithm delivers perfect performance, and that studies of metal-poor main-sequence turnoff stars may be challenged by poor star-galaxy separation. Using the Receiver Operating Characteristic curve, we find a best-to-worst ranking of SVM_best, HB, ML, and SVM_real. We conclude, therefore, that a well trained SVM will outperform template-fitting methods. However, a normally trained SVM performs worse. Thus, Hierarchical Bayesian template fitting may prove to be the optimal classification method in future surveys.

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.

We investigate the properties of dark matter substructure in the gravitational lens HE 0435-1223 (z_l=0.455) via its effects on the positions and flux ratios of the quadruply-imaged background quasar (z_s=1.689). We start with a smooth mass model, add individual, truncated isothermal clumps near the lensed images, and use the Bayesian evidence to compare the quality of different models. Compared with smooth models, models with at least one clump near image A are strongly favored. The mass of this clump within its Einstein radius is log(Mein/Msun) = 7.65 +0.87/-0.84. The Bayesian evidence provides weaker support for a second clump near image B, with log(Mein/Msun) = 6.55 +1.01/-1.51. We also examine models with a full population of substructure, and find the mass fraction in substructure at the Einstein radius to be f_sub > 0.00077, assuming the total clump masses follow a mass function dN/dM proportional to M^(-1.9) over the range M = 10^7-10^10 Msun. Few-clump and population models produce similar Bayesian evidence values, so neither type of model is objectively favored.

We investigate the kinematic and photometric properties of the Segue 3 Milky Way companion using Keck/DEIMOS spectroscopy and Magellan/IMACS g and r-band imaging. Using maximum likelihood methods to analyze the photometry, we study the structure and stellar population of Segue 3. We find the half-light radius of Segue 3 is 26" +/- 5" (2.1 +/- 0.4 pc, for a distance of 17 kpc) and the absolute magnitude is a mere M_V = 0.0 +/- 0.8 mag, making Segue 3 the least luminous old stellar system known. We find Segue 3 to be consistent with a single stellar population, with an age of 12.0 +1.5/-0.4 Gyr and an [Fe/H] of -1.7 +0.07/-0.27. Line-of-sight velocities from the spectra are combined with the photometry to determine a sample of 32 stars which are likely associated with Segue 3. The member stars within three half-light radii have a velocity dispersion of 1.2 +/- 2.6 km/s. Photometry of the members indicates the stellar population has a spread in [Fe/H] of <0.3 dex. These facts, together with the small physical size of Segue 3, imply the object is likely an old, faint stellar cluster which contains no significant dark matter. We find tentative evidence for stellar mass loss in Segue 3 through the eleven candidate member stars outside of three half-light radii, as expected from dynamical arguments. Interpretation of the data outside of three half-light radii, is complicated by the object's spatial coincidence with a previously known halo substructure, which may enhance contamination of our member sample.

We examine the wavelength dependence of flux ratios for six gravitationally lensed quasars using K and L' images obtained at the Gemini North 8m telescope. We select lenses with source redshifts z_s < 2.8 so that K-band images probe rest-frame optical emission from accretion disks, while L'-band images probe rest-frame near-infrared flux emitted (in part) from the more extended surrounding torus. Since the observations correspond to different source sizes, the K and L' flux ratios are sensitive to structure on different scales and may be useful for studying small-structure in the lens galaxies. Four of the six lenses show differences between K and L' flux ratios. In HE 0435$-1223, SDSS 0246-0825, and HE 2149-2745 the differences may be attributable to known microlensing and/or intrinsic variability. In SDSS 0806+2006 the wavelength dependence is not easily attributed to known variations, and may indicate the presence of substructure. By contrast, in Q0142-100 and SBS 0909+523 the K and L' flux ratios are consistent within the uncertainties. We discuss the utility of the current data for studying chromatic effects related to microlensing, dust extinction, and dark matter substructure.

We present low-resolution, rest-frame ~ 5 - 12 micron Spitzer/IRS spectra of two lensed z ~ 2 UV-bright star-forming galaxies, SDSS J120602.09+514229.5 and SDSS J090122.37+181432.3. Using the magnification boost from lensing, we are able to study the physical properties of these objects in greater detail than is possible for unlensed systems. In both targets, we detect strong PAH emission at 6.2, 7.7, and 11.3 microns, indicating the presence of vigorous star formation. For J1206, we find a steeply rising continuum and significant [S IV] emission, suggesting that a moderately hard radiation field is powering continuum emission from small dust grains. The strength of the [S IV] emission also implies a sub-solar metallicity of ~ 0.5 Z_{Sun}, confirming published rest-frame optical measurements. In J0901, the PAH lines have large rest-frame equivalent widths (> 1 micron) and the continuum rises slowly with wavelength, suggesting that any AGN contribution to L_{IR} is insignificant, in contrast to the implications of optical emission-line diagnostics. Using [O III] line flux as a proxy for AGN strength, we estimate that the AGN in J0901 provides only a small fraction of its mid-infrared continuum flux. By combining the detection of [Ar II] with an upper limit on [Ar III] emission, we infer a metallicity of > 1.3 Z_{Sun}. This work highlights the importance of combining rest-frame optical and mid-IR spectroscopy in order to understand the detailed properties of star-forming galaxies at high redshift.

We present a detailed strong lensing analysis of an HST/ACS legacy dataset for the first gravitational lens, Q0957+561. With deep imaging we identify 24 new strongly lensed features, which we use to constrain mass models. We model the stellar component of the lens galaxy using the observed luminosity distribution, and the dark matter halo using several different density profiles. We draw on the weak lensing analysis by Nakajima et al. (2009) to constrain the mass sheet and environmental terms in the lens potential. Adopting the well-measured time delay, we find H_0 = 85 (+14/-13) km/s/Mpc (68% CL) using lensing constraints alone. The principal uncertainties in H_0 are tied to the stellar mass-to-light ratio (a variant of the radial profile degeneracy in lens models). Adding constraints from stellar population synthesis models, we obtain H_0 = 79.3 (+6.7/-8.5) km/s/Mpc (68% CL). We infer that the lens galaxy has a rising rotation curve and a dark matter distribution with an inner core. Intriguingly, we find the quasar flux ratios predicted by our models to be inconsistent with existing radio measurements, suggesting the presence of substructure in the lens.

Attempts to constrain the Hubble constant using the strong gravitational lens system Q0957+561 are limited by systematic uncertainties in the mass model, since the time delay is known very precisely. One important systematic effect is the mass sheet degeneracy, which arises because strong lens modeling cannot constrain the presence or absence of a uniform mass sheet $\kappa$, which rescales $H_0$ by the factor $(1-\kappa)$. In this paper we present new constraints on the mass sheet derived from a weak-lensing analysis of the Hubble Space Telescope imaging of a 6 arcmin square region surrounding the lensed quasar. The average mass sheet within a circular aperture (the strong lens model region) is constrained by integrating the tangential weak gravitational shear over the surrounding area. We find the average convergence within a $30"$ radius around the lens galaxy to be $\kappa(<30") = 0.166\pm0.056$ ($1\sigma$ confidence level), normalized to the quasar redshift. This includes contributions from both the lens galaxy and the surrounding cluster. We also constrain a few other low-order terms in the lens potential by applying a multipole aperture mass formalism to the gravitational shear in an annulus around the strong lensing region. Implications for strong lens models and the Hubble constant are discussed in an accompanying paper.