We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $\mu$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 \deg$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry.

We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk.

We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems.

We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$\sigma$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$).

The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions.

We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30{\deg}-150{\deg} scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed.

The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectroscopic mode with the goal of obtaining precise astrometry and spectroscopy of faint companions to bright stars. We find a clear correlation between residual wavefront error measured by the AO system and the contrast sensitivity by comparing phase errors in observations of the same source, taken on different dates. We find a typical 5-$\sigma$ contrast sensitivity of $2-3~\times~10^{-3}$ at $\sim\lambda/D$. We explore the accuracy of spectral extraction of secondary components of binary systems by recovering the signal from a simulated source injected into several datasets. We outline data reduction procedures unique to GPI's IFS and describe a newly public data pipeline used for the presented analyses. We demonstrate recovery of astrometry and spectroscopy of two known companions to HR 2690 and HD 142527. NRM+polarimetry observations achieve differential visibility precision of $\sigma\sim0.4\%$ in the best case. We discuss its limitations on Gemini-S/GPI for resolving inner regions of protoplanetary disks and prospects for future upgrades. We summarize lessons learned in observing with NRM in spectroscopic and polarimetric modes.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

Ground-based telescopes coupled with adaptive optics (AO) have been playing a leading role in exoplanet direct imaging science and technological development for the past two decades and will continue to have an indispensable role for the next decade and beyond. Over the next decade, extreme AO systems on 8-10m telescopes will 1) mitigate risk for WFIRST-CGI by identifying numerous planets the mission can spectrally characterize, 2) validate performance requirements and motivate improvements to atmosphere models needed to unambiguously characterize solar system-analogues from space, and 3) mature novel technological innovations useful for space. Extremely Large Telescopes can deliver the first thermal infrared (10 $\mu m$) images of rocky planets around Sun-like stars and identify biomarkers. These data provide a future NASA direct imaging flagship mission (i.e. HabEx, LUVOIR) with numerous exo-Earth candidates and critical ancillary information to help clarify whether these planets are habitable.

We report the discovery of a large ($\sim$8500 km diameter) infrared-bright storm at Neptune's equator in June 2017. We tracked the storm over a period of 7 months with high-cadence infrared snapshot imaging, carried out on 14 nights at the 10 meter Keck II telescope and 17 nights at the Shane 120 inch reflector at Lick Observatory. The cloud feature was larger and more persistent than any equatorial clouds seen before on Neptune, remaining intermittently active from at least 10 June to 31 December 2017. Our Keck and Lick observations were augmented by very high-cadence images from the amateur community, which permitted the determination of accurate drift rates for the cloud feature. Its zonal drift speed was variable from 10 June to at least 25 July, but remained a constant $237.4 \pm 0.2$ m s$^{-1}$ from 30 September until at least 15 November. The pressure of the cloud top was determined from radiative transfer calculations to be 0.3-0.6 bar; this value remained constant over the course of the observations. Multiple cloud break-up events, in which a bright cloud band wrapped around Neptune's equator, were observed over the course of our observations. No "dark spot" vortices were seen near the equator in HST imaging on 6 and 7 October. The size and pressure of the storm are consistent with moist convection or a planetary-scale wave as the energy source of convective upwelling, but more modeling is required to determine the driver of this equatorial disturbance as well as the triggers for and dynamics of the observed cloud break-up events.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

Ground-layer adaptive optics (GLAO) systems offer the possibility of improving the "seeing" of large ground-based telescopes and increasing the efficiency and sensitivity of observations over a wide field-of-view. We explore the utility and feasibility of deploying a GLAO system at the W. M. Keck Observatory in order to feed existing and future multi-object spectrographs and wide-field imagers. We also briefly summarize science cases spanning exoplanets to high-redshift galaxy evolution that would benefit from a Keck GLAO system. Initial simulations indicate that a Keck GLAO system would deliver a 1.5x and 2x improvement in FWHM at optical (500 nm) and infrared (1.5 micron), respectively. The infrared instrument, MOSFIRE, is ideally suited for a Keck GLAO feed in that it has excellent image quality and is on the telescope's optical axis. However, it lacks an atmospheric dispersion compensator, which would limit the minimum usable slit size for long-exposure science cases. Similarly, while LRIS and DEIMOS may be able to accept a GLAO feed based on their internal image quality, they lack either an atmospheric dispersion compensator (DEIMOS) or flexure compensation (LRIS) to utilize narrower slits matched to the GLAO image quality. However, some science cases needing shorter exposures may still benefit from Keck GLAO and we will investigate the possibility of installing an ADC.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

The Gemini Planet Imager Exoplanet Survey (GPIES) is a multi-year direct imaging survey of 600 stars to discover and characterize young Jovian exoplanets and their environments. We have developed an automated data architecture to process and index all data related to the survey uniformly. An automated and flexible data processing framework, which we term the Data Cruncher, combines multiple data reduction pipelines together to process all spectroscopic, polarimetric, and calibration data taken with GPIES. With no human intervention, fully reduced and calibrated data products are available less than an hour after the data are taken to expedite follow-up on potential objects of interest. The Data Cruncher can run on a supercomputer to reprocess all GPIES data in a single day as improvements are made to our data reduction pipelines. A backend MySQL database indexes all files, which are synced to the cloud, and a front-end web server allows for easy browsing of all files associated with GPIES. To help observers, quicklook displays show reduced data as they are processed in real-time, and chatbots on Slack post observing information as well as reduced data products. Together, the GPIES automated data processing architecture reduces our workload, provides real-time data reduction, optimizes our observing strategy, and maintains a homogeneously reduced dataset to study planet occurrence and instrument performance.

We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf.

With a large, unique spectroscopic survey in the fields of 28 galaxy-scale strong gravitational lenses, we identify groups of galaxies in the 26 adequately-sampled fields. Using a group finding algorithm, we find 210 groups with at least five member galaxies; the median number of members is eight. Our sample spans redshifts of 0.04 $\le z_{grp} \le$ 0.76 with a median of 0.31, including 174 groups with $0.1 < z_{grp} < 0.6$. Groups have radial velocity dispersions of 60 $\le \sigma_{grp} \le$ 1200 km s$^{-1}$ with a median of 350 km s$^{-1}$. We also discover a supergroup in field B0712+472 at $z =$ 0.29 consisting of three main groups. We recover groups similar to $\sim$ 85% of those previously reported in these fields within our redshift range of sensitivity and find 187 new groups with at least five members. The properties of our group catalog, specifically 1) the distribution of $\sigma_{grp}$, 2) the fraction of all sample galaxies that are group members, and 3) the fraction of groups with significant substructure, are consistent with those for other catalogs. The distribution of group virial masses agrees well with theoretical expectations. Of the lens galaxies, 12 of 26 (46%) (B1422+231, B1600+434, B2114+022, FBQS J0951+2635, HE0435-1223, HST J14113+5211, MG0751+2716, MGJ1654+1346, PG 1115+080, Q ER 0047-2808, RXJ1131-1231, and WFI J2033-4723) are members of groups with at least five galaxies, and one more (B0712+472) belongs to an additional, visually identified group candidate. There are groups not associated with the lens that still are likely to affect the lens model; in six of 25 (24%) fields (excluding the supergroup), there is at least one massive ($\sigma_{grp} \ge$ 500 km s$^{-1}$) group or group candidate projected within 2$^{\prime}$ of the lens.

Since 2013, observations of Neptune with small telescopes have resulted in several detections of long-lived bright atmospheric features that have also been observed by large telescopes such as Keck II or Hubble. The combination of both types of images allows the study of the long term evolution of major cloud systems in the planet. In 2013 and 2014 two bright features were present on the planet at southern mid latitudes. These may have merged in late 2014, possibly leading to the formation of a single bright feature observed during 2015 at the same latitude. This cloud system was first observed in January 2015 and nearly continuously from July to December 2015 in observations with telescopes in the 2 to 10 meter class and in images from amateur astronomers. These images show the bright spot as a compact feature at 40.1 deg South planetographic latitude well resolved from a nearby bright zonal band that extended from 42 deg South to 20 deg South. Tracking its motion from July to November 2015 suggests a longitudinal oscillation of 16 deg in amplitude with a 90 day period, typical of dark spots on Neptune and similar to the Great Red Spot oscillation in Jupiter. The limited time covered by high-resolution observations only covers one full oscillation and other interpretations of the changing motions could be possible. HST images in September 2015 show the presence of a dark spot at short wavelengths in the southern flank of the bright cloud observed throughout 2015.

At a distance of 2~pc, our nearest brown dwarf neighbor, Luhman 16 AB, has been extensively studied since its discovery 3 years ago, yet its most fundamental parameter -- the masses of the individual dwarfs -- has not been constrained with precision. In this work we present the full astrometric orbit and barycentric motion of Luhman 16 AB and the first precision measurements of the individual component masses. We draw upon archival observations spanning 31 years from the European Southern Observatory (ESO) Schmidt Telescope, the Deep Near-Infrared Survey of the Southern Sky (DENIS), public FORS2 data on the Very Large Telescope (VLT), and new astrometry from the Gemini South Multiconjugate Adaptive Optics System (GeMS). Finally, we include three radial velocity measurements of the two components from VLT/CRIRES, spanning one year. With this new data sampling a full period of the orbit, we use a Markov Chain Monte Carlo algorithm to fit a 16-parameter model incorporating mutual orbit and barycentric motion parameters and constrain the individual masses to be~$27.9^{+1.1}_{-1.0}$~$M_{J}$ for the T dwarf and~$34.2^{+1.3}_{-1.1}$~$M_{J}$ for the L dwarf. Our measurements of Luhman 16 AB's mass ratio and barycentric motion parameters are consistent with previous estimates in the literature utilizing recent astrometry only. The GeMS-derived measurements of the Luhman 16 AB separation in 2014-2015 agree closely with Hubble Space Telescope (HST) measurements made during the same epoch Bedin et al. 2017, and the derived mutual orbit agrees with those measurements to within the HST uncertainties of $0.3 - 0.4$ milliarcseconds.

We present a combined strong and weak lensing analysis of the J085007.6+360428 (J0850) field, which was selected by its high projected concentration of luminous red galaxies and contains the massive cluster Zwicky 1953. Using Subaru/Suprime-Cam $BVR_{c}I_{c}i^{\prime}z^{\prime}$ imaging and MMT/Hectospec spectroscopy, we first perform a weak lensing shear analysis to constrain the mass distribution in this field, including the cluster at $z = 0.3774$ and a smaller foreground halo at $z = 0.2713$. We then add a strong lensing constraint from a multiply-imaged galaxy in the imaging data with a photometric redshift of $z \approx 5.03$. Unlike previous cluster-scale lens analyses, our technique accounts for the full three-dimensional mass structure in the beam, including galaxies along the line of sight. In contrast with past cluster analyses that use only lensed image positions as constraints, we use the full surface brightness distribution of the images. This method predicts that the source galaxy crosses a lensing caustic such that one image is a highly-magnified "fold arc", which could be used to probe the source galaxy's structure at ultra-high spatial resolution ($< 30$ pc). We calculate the mass of the primary cluster to be $\mathrm{M_{vir}} = 2.93_{-0.65}^{+0.71} \times 10^{15}~\mathrm{M_{\odot}}$ with a concentration of $\mathrm{c_{vir}} = 3.46_{-0.59}^{+0.70}$, consistent with the mass-concentration relation of massive clusters at a similar redshift. The large mass of this cluster makes J0850 an excellent field for leveraging lensing magnification to search for high-redshift galaxies, competitive with and complementary to that of well-studied clusters such as the HST Frontier Fields.

We present spectro-photometry spanning 1-5 $\mu$m of 51 Eridani b, a 2-10 M$_\text{Jup}$ planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new $K1$ (1.90-2.19 $\mu$m) and $K2$ (2.10-2.40 $\mu$m) spectra taken with the Gemini Planet Imager as well as an updated $L_P$ (3.76 $\mu$m) and new $M_S$ (4.67 $\mu$m) photometry from the NIRC2 Narrow camera. The new data were combined with $J$ (1.13-1.35 $\mu$m) and $H$ (1.50-1.80 $\mu$m) spectra from the discovery epoch with the goal of better characterizing the planet properties. 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4-T8) and we propose that 51 Eri b might be in the process of undergoing the transition from L-type to T-type. We used two complementary atmosphere model grids including either deep iron/silicate clouds or sulfide/salt clouds in the photosphere, spanning a range of cloud properties, including fully cloudy, cloud free and patchy/intermediate opacity clouds. Model fits suggest that 51 Eri b has an effective temperature ranging between 605-737 K, a solar metallicity, a surface gravity of $\log$(g) = 3.5-4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the SED. From the model atmospheres, we infer a luminosity for the planet of -5.83 to -5.93 ($\log L/L_{\odot}$), leaving 51 Eri b in the unique position as being one of the only directly imaged planet consistent with having formed via cold-start scenario. Comparisons of the planet SED against warm-start models indicates that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M$_{\oplus}$.

We present H band spectroscopic and Halpha photometric observations of HD 100546 obtained with GPI and MagAO. We detect H band emission at the location of the protoplanet HD 100546b, but show that choice of data processing parameters strongly affects the morphology of this source. It appears point-like in some aggressive reductions, but rejoins an extended disk structure in the majority of the others. Furthermore, we demonstrate that this emission appears stationary on a timescale of 4.6 yrs, inconsistent at the 2sigma level with a Keplerian clockwise orbit at 59 au in the disk plane. The H band spectrum of the emission is inconsistent with any type of low effective temperature object or accreting protoplanetary disk. It strongly suggests a scattered light origin, as it is consistent with the spectrum of the star and the spectra extracted at other locations in the disk. A non detection at the 5sigma level of HD 100546b in differential Halpha imaging places an upper limit, assuming the protoplanet lies in a gap free of extinction, on the accretion luminosity and accretion rate of 1.7E-4 Lsun and MMdot<6.4E-7Mjup^2/yr for 1Rjup. These limits are comparable to the accretion luminosity and rate of TTauri-stars or LkCa 15b. Taken together, these lines of evidence suggest that the H band source at the location of HD 100546b is not emitted by a planetary photosphere or an accreting circumplanetary disk but is a disk feature enhanced by the PSF subtraction process. This non-detection is consistent with the non-detection in the K band reported in an earlier study but does not exclude the possibility that HD 100546b is deeply embedded.

We present optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at H{\alpha}. The new GPI H -band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple datasets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 c varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.

We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12--1.3 micron) and H (1.50--1.80 micron) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 AU (70 year) orbit with a 68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119 degrees with a 68% confidence interval between 114 degrees and 125 degrees. To address considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare to known spectral types using field brown dwarf and low-mass star spectra and find a best fit match of a field gravity M6.5+/-1.5 spectral type with a corresponding temperature of 2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun) = -2.88+/-0.07 dex, using DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34+/-1 to 95+/-4 M_Jup. These results are consistent with previous measurements of the object.

The Gemini Planet Imager is a high-contrast near-infrared instrument specifically designed to image exoplanets and circumstellar disks over a narrow field of view. We use science data and AO telemetry taken during the first 1.5 yr of the GPI Exoplanet Survey to quantify the performance of the AO system. In a typical 60 sec H-band exposure, GPI achieves a 5$\sigma$ raw contrast of 10$^{-4}$ at 0.4"; typical final 5$\sigma$ contrasts for full 1 hr sequences are more than 10 times better than raw contrasts. We find that contrast is limited by bandwidth wavefront error over much of the PSF. Preliminary exploratory factor analysis can explain 60-70% of the variance in raw contrasts with combinations of seeing and wavefront error metrics. We also examine the effect of higher loop gains on contrast by comparing wavefront error maps reconstructed from AO telemetry to concurrent IFS images. These results point to several ways that GPI performance could be improved in software or hardware.

The HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2-2.3 microns that further constrains its outer morphology (projected separations of 27-135 AU). We also present complementary Gemini Planet Imager 1.6 micron total intensity and polarized light detections that probe down to projected separations less than 10 AU. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40-52 AU and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 AU to a Jupiter mass at 5 AU.

We present $H$-band near-infrared polarimetric imaging observations of the F5V star HD~157587 obtained with the Gemini Planet Imager (GPI) that reveal the debris disk as a bright ring structure at a separation of $\sim$80$-$100~AU. The new GPI data complement recent HST/STIS observations that show the disk extending out to over 500~AU. The GPI image displays a strong asymmetry along the projected minor axis as well as a fainter asymmetry along the projected major axis. We associate the minor and major axis asymmetries with polarized forward scattering and a possible stellocentric offset, respectively. To constrain the disk geometry we fit two separate disk models to the polarized image, each using a different scattering phase function. Both models favor a disk inclination of $\sim 70\degr$ and a $1.5\pm0.6$ AU stellar offset in the plane of the sky along the projected major axis of the disk. We find that the stellar offset in the disk plane, perpendicular to the projected major axis is degenerate with the form of the scattering phase function and remains poorly constrained. The disk is not recovered in total intensity due in part to strong adaptive optics residuals, but we recover three point sources. Considering the system's proximity to the galactic plane and the point sources' positions relative to the disk, we consider it likely that they are background objects and unrelated to the disk's offset from the star.

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Lo\`eve Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $\beta$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$\sigma$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $\beta$ Pic b firmly rule out the possibility of a transit of the planet at 10-$\sigma$ significance. However, we confirm that the Hill sphere of $\beta$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.

Using the Gemini Planet Imager (GPI), we have resolved the circumstellar debris disk around HD 111520 at a projected range of ~30-100 AU in both total and polarized $H$-band intensity. The disk is seen edge-on at a position angle of ~165$^{\circ}$ along the spine of emission. A slight inclination or asymmetric warping are covariant and alters the interpretation of the observed disk emission. We employ 3 point spread function (PSF) subtraction methods to reduce the stellar glare and instrumental artifacts to confirm that there is a roughly 2:1 brightness asymmetry between the NW and SE extension. This specific feature makes HD 111520 the most extreme examples of asymmetric debris disks observed in scattered light among similar highly inclined systems, such as HD 15115 and HD 106906. We further identify a tentative localized brightness enhancement and scale height enhancement associated with the disk at ~40 AU away from the star on the SE extension. We also find that the fractional polarization rises from 10 to 40% from 0.5" to 0.8" from the star. The combination of large brightness asymmetry and symmetric polarization fraction leads us to believe that an azimuthal dust density variation is causing the observed asymmetry.

We present astrometric monitoring of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager between 2013 and 2016. A small but significant position angle change is detected at constant separation; the orbital motion is confirmed with literature measurements. Efficient Monte Carlo techniques place preliminary constraints on the orbital parameters of HD 95086 b. With 68% confidence, a semimajor axis of 61.7^{+20.7}_{-8.4} au and an inclination of 153.0^{+9.7}_{-13.5} deg are favored, with eccentricity less than 0.21. Under the assumption of a co-planar planet-disk system, the periastron of HD 95086 b is beyond 51 au with 68% confidence. Therefore HD 95086 b cannot carve the entire gap inferred from the measured infrared excess in the SED of HD 95086. We use our sensitivity to additional planets to discuss specific scenarios presented in the literature to explain the geometry of the debris belts. We suggest that either two planets on moderately eccentric orbits or three to four planets with inhomogeneous masses and orbital properties are possible. The sensitivity to additional planetary companions within the observations presented in this study can be used to help further constrain future dynamical simulations of the planet-disk system.

We present new $H$ (1.5-1.8 $\mu$m) photometric and $K_1$ (1.9-2.2 $\mu$m) spectroscopic observations of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager. The $H$-band magnitude has been significantly improved relative to previous measurements, whereas the low resolution $K_1$ ($\lambda/\delta\lambda \approx 66$) spectrum is featureless within the measurement uncertainties, and presents a monotonically increasing pseudo-continuum consistent with a cloudy atmosphere. By combining these new measurements with literature $L^{\prime}$ photometry, we compare the spectral energy distribution of the planet to other young planetary-mass companions, field brown dwarfs, and to the predictions of grids of model atmospheres. HD 95086 b is over a magnitude redder in $K_1-L^{\prime}$ color than 2MASS J12073346-3932539 b and HR 8799 c and d, despite having a similar $L^{\prime}$ magnitude. Considering only the near-infrared measurements, HD 95086 b is most analogous to the brown dwarfs 2MASS J2244316+204343 and 2MASS J21481633+4003594, both of which are thought to have dusty atmospheres. Morphologically, the spectral energy distribution of HD 95086 b is best fit by low temperature ($T_{\rm eff} =$ 800-1300 K), low surface gravity spectra from models which simulate high photospheric dust content. This range of effective temperatures is consistent with field L/T transition objects, but the spectral type of HD 95086 b is poorly constrained between early L and late T due to its unusual position the color-magnitude diagram, demonstrating the difficulty in spectral typing young, low surface gravity substellar objects. As one of the reddest such objects, HD 95086 b represents an important empirical benchmark against which our current understanding of the atmospheric properties of young extrasolar planets can be tested.

Detecting polarized light from self-luminous exoplanets has the potential to provide key information about rotation, surface gravity, cloud grain size, and cloud coverage. While field brown dwarfs with detected polarized emission are common, no exoplanet or substellar companion has yet been detected in polarized light. With the advent of high contrast imaging spectro-polarimeters such as GPI and SPHERE, such a detection may now be possible with careful treatment of instrumental polarization. In this paper, we present 28 minutes of $H$-band GPI polarimetric observations of the benchmark T5.5 companion HD 19467 B. We detect no polarization signal from the target, and place an upper limit on the degree of linear polarization of $p_{\text{CL}99.73\%} \leq 2.4\%$. We discuss our results in the context of T dwarf cloud models and photometric variability.

We present the first scattered-light image of the debris disk around HD 131835 in $H$ band using the Gemini Planet Imager. HD 131835 is a $\sim$15 Myr old A2IV star at a distance of $\sim$120 pc in the Sco-Cen OB association. We detect the disk only in polarized light and place an upper limit on the peak total intensity. No point sources resembling exoplanets were identified. Compared to its mid-infrared thermal emission, the disk in scattered light shows similar orientation but different morphology. The scattered-light disk extends from $\sim$75 to $\sim$210 AU in the disk plane with roughly flat surface density. Our Monte Carlo radiative transfer model can well describe the observations with a model disk composed of a mixture of silicates and amorphous carbon. In addition to the obvious brightness asymmetry due to stronger forward scattering, we discover a weak brightness asymmetry along the major axis with the northeast side being 1.3 times brighter than the southwest side at a 3-{\sigma} level.

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to $2\times10^{-7}$, an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of $14^{+7}_{-3}$ AU, corresponding to a period of $41^{+35}_{-12}$ years (assuming a mass of $1.75$ $M_{\odot}$ for the central star), and an inclination of $138^{+15}_{-13}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit which does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.

We present the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius $\sim$50 AU, and an outer extent $>$500 AU. The HST data show the outer regions are highly asymmetric, resembling the ''needle'' morphology seen for the HD 15115 debris disk. The planet candidate is oriented $\sim$21$\deg$ away from the position angle of the primary's debris disk, strongly suggesting non-coplanarity with the system. We hypothesize that HD 106906b could be dynamically involved in the perturbation of the primary's disk, and investigate whether or not there is evidence for a circumplanetary dust disk or cloud that is either primordial or captured from the primary. We show that both the existing optical properties and near-infrared colors of HD 106906b are weakly consistent with this possibility, motivating future work to test for the observational signatures of dust surrounding the planet.

SN 2006gy was the most luminous SN ever observed at the time of its discovery and the first of the newly defined class of superluminous supernovae (SLSNe). The extraordinary energetics of SN 2006gy and all SLSNe (>10^51 erg) require either atypically large explosion energies (e.g., pair-instability explosion) or the efficient conversion of kinetic into radiative energy (e.g., shock interaction). The mass-loss characteristics can therefore offer important clues regarding the progenitor system. For the case of SN 2006gy, both a scattered and thermal light echo from circumstellar material (CSM) have been reported at later epochs (day ~800), ruling out the likelihood of a pair-instability event and leading to constraints on the characteristics of the CSM. Owing to the proximity of the SN to the bright host-galaxy nucleus, continued monitoring of the light echo has not been trivial, requiring the high resolution offered by the Hubble Space Telescope (HST) or ground-based adaptive optics (AO). Here we report detections of SN 2006gy using HST and Keck AO at ~3000 days post-explosion and consider the emission mechanism for the very late-time light curve. While the optical light curve and optical spectral energy distribution are consistent with a continued scattered-light echo, a thermal echo is insufficient to power the K'-band emission by day 3000. Instead, we present evidence for late-time infrared emission from dust that is radiatively heated by CSM interaction within an extremely dense dust shell, and we consider the implications on the CSM characteristics and progenitor system.

We present $H$-band observations of $\beta$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $\beta$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($\phi=85.27{\deg}^{+0.26}_{-0.19}$) with a position angle $\theta_{PA}=30.35{\deg}^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $\theta_{PA}\approx29{\deg}$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $\beta$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $\Omega=31.75{\deg}\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $\beta$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $\beta$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $\beta$ Pic b; 2) the mutual inclination of the inner disk and $\beta$ Pic $b$ is $4{\deg}$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $\beta$ Pic $b$ or self-stirring by the disk's parent bodies.

We present Gemini Planet Imager (GPI) observations of AU Microscopii, a young M dwarf with an edge-on, dusty debris disk. Integral field spectroscopy and broadband imaging polarimetry were obtained during the commissioning of GPI. In our broadband imaging polarimetry observations, we detect the disk only in total intensity and find asymmetries in the morphology of the disk between the southeast and northwest sides. The southeast side of the disk exhibits a bump at 1$''$ (10 AU projected separation) that is three times more vertically extended and three times fainter in peak surface brightness than the northwest side at similar separations. This part of the disk is also vertically offset by 69$\pm$30 mas to the northeast at 1$''$ when compared to the established disk mid-plane and consistent with prior ALMA and Hubble Space Telescope/STIS observations. We see hints that the southeast bump might be a result of detecting a horizontal sliver feature above the main disk that could be the disk backside. Alternatively when including the morphology of the northwest side, where the disk mid-plane is offset in the opposite direction $\sim$50 mas between 0$.''$4 and 1$.''$2, the asymmetries suggest a warp-like feature. Using our integral field spectroscopy data to search for planets, we are 50% complete for $\sim$4 $M_\mathrm{Jup}$ planets at 4 AU. We detect a source, resolved only along the disk plane, that could either be a candidate planetary mass companion or a compact clump in the disk.

We search the complete Hubble Frontier Field dataset of Abell 2744 and its parallel field for z~10 sources to further refine the evolution of the cosmic star-formation rate density (SFRD) at z>8. We independently confirm two images of the recently discovered triply-imaged z~9.8 source by Zitrin et al. (2014) and set an upper limit for similar z~10 galaxies with red colors of J_125-H_160>1.2 in the parallel field of Abell 2744. We utilize extensive simulations to derive the effective selection volume of Lyman-break galaxies at z~10, both in the lensed cluster field and in the adjacent parallel field. Particular care is taken to include position-dependent lensing shear to accurately account for the expected sizes and morphologies of highly-magnified sources. We show that both source blending and shear reduce the completeness at a given observed magnitude in the cluster, particularly near the critical curves. These effects have a significant, but largely overlooked, impact on the detectability of high-redshift sources behind clusters, and substantially reduce the expected number of highly-magnified sources. The detections and limits from both pointings result in a SFRD which is higher by 0.4+-0.4 dex than previous estimates at z~10 from blank fields. Nevertheless, the combination of these new results with all other estimates remain consistent with a rapidly declining SFRD in the 170 Myr from z~8 to z~10 as predicted by cosmological simulations and dark-matter halo evolution in LambdaCDM. Once biases introduced by magnification-dependent completeness are accounted for, the full six cluster and parallel Frontier Field program will be an extremely powerful new dataset to probe the evolution of the galaxy population at z>8 before the advent of the JWST.

We present the laboratory verification of a method for re- moving the effects of frozen-flow atmospheric turbulence using a Linear Quadratic Gaussian (LQG) controller, also known as a Kalman Filter. This method, which we term "Predictive Fourier Control," can identify correlated atmospheric motions due to layers of frozen flow turbulence, and can predictively remove the effects of these correlated motions in real-time. Our laboratory verification suggests a factor of 3 improvement in the RMS residual wavefront error and a 10% improvement in measured Strehl of the system. We found that the RMS residual wavefront error was suppressed from 35.0 nm to 11.2 nm due to the use of Predictive Fourier Control, and that the far field Strehl improved from 0.479 to 0.520.

The Gemini Multi-conjugate adaptive optics System (GeMS) is a facility instrument for the Gemini-South telescope. It delivers uniform, near-diffraction-limited image quality at near-infrared wavelengths over a 2 arcminute field of view. Together with the Gemini South Adaptive Optics Imager (GSAOI), a near-infrared wide field camera, GeMS/GSAOI's combination of high spatial resolution and a large field of view will make it a premier facility for precision astrometry. Potential astrometric science cases cover a broad range of topics including exo-planets, star formation, stellar evolution, star clusters, nearby galaxies, black holes and neutron stars, and the Galactic center. In this paper, we assess the astrometric performance and limitations of GeMS/GSAOI. In particular, we analyze deep, mono-epoch images, multi-epoch data and distortion calibration. We find that for single-epoch, un-dithered data, an astrometric error below 0.2 mas can be achieved for exposure times exceeding one minute, provided enough stars are available to remove high-order distortions. We show however that such performance is not reproducible for multi-epoch observations, and an additional systematic error of ~0.4 mas is evidenced. This systematic multi-epoch error is the dominant error term in the GeMS/GSAOI astrometric error budget, and it is thought to be due to time-variable distortion induced by gravity flexure.

Certain configurations of massive structures projected along the line of sight maximize the number of detections of gravitationally lensed $z\sim10$ galaxies. We characterize such lines of sight with the \'etendue $\sigma_\mu$, the area in the source plane magnified over some threshold $\mu$. We use the Millennium I and Millennium XXL cosmological simulations to determine the frequency of high $\sigma_\mu$ beams on the sky, their properties, and efficient selection criteria. We define the best beams as having $\sigma_{\mu>3} >2000$ arcsec$^2$, for a $z\sim10$ source plane, and predict $477 \pm 21$ such beams on the sky. The total mass in the beam and $\sigma_{\mu>3}$ are strongly correlated. After controlling for total mass, we find a significant residual correlation between $\sigma_{\mu>3}$ and the number of cluster-scale halos ($>10^{14} M_\odot h^{-1}$) in the beam. Beams with $\sigma_{\mu>3} >2000$ arcsec$^2$, which should be best at lensing $z\sim10$ galaxies, are ten times more likely to contain multiple cluster-scale halos than a single cluster-scale halo. Beams containing an Abell 1689-like massive cluster halo often have additional structures along the line of sight, including at least one additional cluster-scale ($M_{200}>10^{14}M_\odot h^{-1}$) halo 28% of the time. Selecting beams with multiple, massive structures will lead to enhanced detection of the most distant and intrinsically faint galaxies.

Lines of sight with multiple, projected, cluster-scale halos have high total masses and complex lens plane interactions that can boost the area of magnification, or etendue, making detection of faint background sources more likely than elsewhere. To identify these new "compound" cosmic telescopes, we have found lines-of-sight with the highest integrated mass densities, as traced by the projected concentrations of Luminous Red Galaxies (LRGs). We use 1151 MMT Hectospec spectra to derive preliminary magnification maps for two such lines of sight with total mass exceeding ~ 3 x 10$^{15}$ Msun -- J0850+3604 (0850) and J1306+4632 (1306). We identify 2-3 group- and cluster-scale halos in each beam over 0.1 < z < 0.7, all of which are well-traced by LRGs. In Subaru Suprime-Cam imaging of beam 0850, we discover serendipitously a candidate multiply-imaged V-dropout source at z = 5.03, whose location is consistent with the critical curves for a source plane of $z_s$ = 5.03 predicted by our mass model. Incorporating the position of the candidate multiply-imaged galaxy as a constraint on the critical curve location in 0850 narrows the 68% confidence band on lens plane area with mu > 10 for a source plane of $z_s$ = 10 to [1.8, 4.2] square arcminutes, comparable to that of MACS 0717+3745 and El Gordo, two of the most powerful known single cluster lenses. The 68% confidence intervals on the lens plane area with mu > 10 for 1306 are [2.3, 6.7] square arcminutes. The significant lensing power of our beams makes them powerful probes of reionization and galaxy formation in the early Universe.

The best gravitational lenses for detecting distant galaxies are those with the largest mass concentrations and the most advantageous configurations of that mass along the line of sight. Our new method for finding such gravitational telescopes uses optical data to identify projected concentrations of luminous red galaxies (LRGs). LRGs are biased tracers of the underlying mass distribution, so lines of sight with the highest total luminosity in LRGs are likely to contain the largest total mass. We apply this selection technique to the Sloan Digital Sky Survey and identify the 200 fields with the highest total LRG luminosities projected within a 3.5' radius over the redshift range 0.1 < z < 0.7. The redshift and angular distributions of LRGs in these fields trace the concentrations of non-LRG galaxies. These fields are diverse; 22.5% contain one known galaxy cluster and 56.0% contain multiple known clusters previously identified in the literature. Thus, our results confirm that these LRGs trace massive structures and that our selection technique identifies fields with large total masses. These fields contain 2-3 times higher total LRG luminosities than most known strong-lensing clusters and will be among the best gravitational lensing fields for the purpose of detecting the highest redshift galaxies.

High-precision astrometry can identify exoplanets and measure their orbits and masses, while coronagraphic imaging enables detailed characterization of their physical properties and atmospheric compositions through spectroscopy. In a previous paper, we showed that a diffractive pupil telescope (DPT) in space can enable sub-microarcsecond accuracy astrometric measurements from wide-field images by creating faint but sharp diffraction spikes around the bright target star. The DPT allows simultaneous astrometric measurement and coronagraphic imaging, and we discuss and quantify in this paper the scientific benefits of this combination for exoplanet science investigations: identification of exoplanets with increased sensitivity and robustness, and ability to measure planetary masses to high accuracy. We show how using both measurements to identify planets and measure their masses offers greater sensitivity and provides more reliable measurements than possible with separate missions, and therefore results in a large gain in mission efficiency. The combined measurements reliably identify potentially habitable planets in multiple systems with a few observations, while astrometry or imaging alone would require many measurements over a long time baseline. In addition, the combined measurement allows direct determination of stellar masses to percent-level accuracy, using planets as test particles. We also show that the DPT maintains the full sensitivity of the telescope for deep wide-field imaging, and is therefore compatible with simultaneous scientific observations unrelated to exoplanets. We conclude that astrometry, coronagraphy, and deep wide-field imaging can be performed simultaneously on a single telescope without significant negative impact on the performance of any of the three techniques.

Astrometric detection and mass determination of Earth-mass exoplanets requires sub-microarcsec accuracy, which is theoretically possible with an imaging space telescope using field stars as an astrometric reference. The measurement must however overcome astrometric distortions which are much larger than the photon noise limit. To address this issue, we propose to generate faint stellar diffraction spikes using a two-dimensional grid of regularly spaced small dark spots added to the surface of the primary mirror (PM). Accurate astrometric motion of the host star is obtained by comparing the position of the spikes to the background field stars. The spikes do not contribute to scattered light in the central part of the field and therefore allow unperturbed coronagraphic observation of the star's immediate surrounding. Because the diffraction spikes are created on the PM and imaged on the same focal plane detector as the background stars, astrometric distortions affect equally the diffraction spikes and the background stars, and are therefore calibrated. We describe the technique, detail how the data collected by the wide-field camera are used to derive astrometric motion, and identify the main sources of astrometric error using numerical simulations and analytical derivations. We find that the 1.4 m diameter telescope, 0.3 sq.deg field we adopt as a baseline design achieves 0.2 microarcsec single measurement astrometric accuracy. The diffractive pupil concept thus enables sub-microarcsec astrometry without relying on the accurate pointing, external metrology or high stability hardware required with previously proposed high precision astrometry concepts.

We investigate the gravitational lensing properties of lines of sight containing multiple cluster-scale halos, motivated by their ability to lens very high-redshift (z ~ 10) sources into detectability. We control for the total mass along the line of sight, isolating the effects of distributing the mass among multiple halos and of varying the physical properties of the halos. Our results show that multiple-halo lines of sight can increase the magnified source-plane region compared to the single cluster lenses typically targeted for lensing studies, and thus are generally better fields for detecting very high-redshift sources. The configurations that result in optimal lensing cross sections benefit from interactions between the lens potentials of the halos when they overlap somewhat on the sky, creating regions of high magnification in the source plane not present when the halos are considered individually. The effect of these interactions on the lensing cross section can even be comparable to changing the total mass of the lens from 10^15 M_sun to 3x10^15 M_sun. The gain in lensing cross section increases as the mass is split into more halos, provided that the lens potentials are projected close enough to interact with each other. A nonzero projected halo angular separation, equal halo mass ratio, and high projected halo concentration are the best mass configurations, whereas projected halo ellipticity, halo triaxiality, and the relative orientations of the halos are less important. Such high mass, multiple-halo lines of sight exist in the SDSS.

NGC 6240 is a pair of colliding disk galaxies, each with a black hole in its core. We have used laser guide star adaptive optics on the Keck II telescope to obtain high-resolution ($\sim 0.06$") near-infrared integral-field spectra of the region surrounding the supermassive black hole in the south nucleus of this galaxy merger. We use the K-band CO absorption bandheads to trace stellar kinematics. We obtain a spatial resolution of about 20 pc and thus directly resolve the sphere of gravitational influence of the massive black hole. We explore two different methods to measure the black hole mass. Using a Jeans Axisymmetric Multi-Gaussian mass model, we investigate the limit that a relaxed mass distribution produces all of the measured velocity dispersion, and find an upper limit on the black hole mass at $2.0 \pm 0.2 \times 10^9 M_{\sun}$. When assuming the young stars whose spectra we observe remain in a thin disk, we compare Keplerian velocity fields to the measured two-dimensional velocity field measured and fit for a mass profile containing a black hole point mass plus a radially-varying spherical component, which suggests a lower limit for the black hole mass of $8.7 \pm 0.3 \times 10^8 M_{\sun}$. Our measurements of the stellar velocity dispersion place this AGN within the scatter of the $M_{BH}$-$\sigma_{*}$ relation. As NGC 6240 is a merging system, this may indicate that the relation is preserved during a merger at least until the final coalescence of the two nuclei.

We present uncontaminated rest-frame u - R colors of 78 X-ray-selected AGN hosts at 0.5 < z < 1.5 in the Chandra Deep Fields measured with HST/ACS/NICMOS and VLT/ISAAC imaging. We also present spatially-resolved NUV - R color gradients for a subsample of AGN hosts imaged by HST/WFC3. Integrated, uncorrected photometry is not reliable for comparing the mean properties of soft and hard AGN host galaxies at z ~ 1 due to color contamination from point-source AGN emission. We use a cloning simulation to develop a calibration between concentration and this color contamination and use this to correct host galaxy colors. The mean u - R color of the unobscured/soft hosts beyond ~6 kpc is statistically equivalent to that of the obscured/hard hosts (the soft sources are 0.09 +/- 0.16 magnitudes bluer). Furthermore, the rest-frame V - J colors of the obscured and unobscured hosts beyond ~6 kpc are statistically equivalent, suggesting that the two populations have similar distributions of dust extinction. For the WFC3/IR sample, the mean NUV - R color gradients of unobscured and obscured sources differ by less than ~0.5 magnitudes for r > 1.1 kpc. These three observations imply that AGN obscuration is uncorrelated with the star formation rate beyond ~1 kpc. These observations favor a unification scenario for intermediate-luminosity AGNs in which obscuration is determined geometrically. Scenarios in which the majority of intermediate-luminosity AGN at z ~ 1 are undergoing rapid, galaxy-wide quenching due to AGN-driven feedback processes are disfavored.

Active galaxies hosting two accreting and merging super-massive black holes (SMBHs) -- dual Active Galactic Nuclei (AGN) -- are predicted by many current and popular models of black hole-galaxy co-evolution. We present here the results of a program that has identified a set of probable dual AGN candidates based on near Infra-red (NIR) Laser Guide-Star Adaptive Optics (LGS AO) imaging with the Keck II telescope. These candidates are selected from a complete sample of radio-quiet Quasi-stellar Objects (QSOs) drawn from the Sloan Digital Sky Survey (SDSS), which show double-peaked narrow AGN emission lines. Of the twelve AGNs imaged, we find six with double galaxy structure, of which four are in galaxy mergers. We measure the ionization of the two velocity components in the narrow AGN lines to test the hypothesis that both velocity components come from an active nucleus. The combination of a well-defined parent sample and high-quality imaging allows us to place constraints on the fraction of SDSS QSOs that host dual accreting black holes separated on kiloparsec (kpc) scales: ~0.3%-0.65%. We derive from this fraction the time spent in a QSO phase during a typical merger and find a value that is much lower than estimates that arise from QSO space densities and galaxy merger statistics. We discuss possible reasons for this difference. Finally, we compare the SMBH mass distributions of single and dual AGN and find little difference between the two within the limited statistics of our program, hinting that most SMBH growth happens in the later stages of a merger process.

We investigate the utility of the asymptotic giant branch (AGB) and the red giant branch (RGB) as probes of the star formation history (SFH) of the nearby (D=2.5 Mpc) dwarf irregular galaxy, KKH 98. Near-infrared (IR) Keck Laser Guide Star Adaptive Optics (AO) images resolve 592 IR bright stars reaching over 1 magnitude below the Tip of the Red Giant Branch. Significantly deeper optical (F475W and F814W) Hubble Space Telescope images of the same field contain over 2500 stars, reaching to the Red Clump and the Main Sequence turn-off for 0.5 Gyr old populations. Compared to the optical color magnitude diagram (CMD), the near-IR CMD shows significantly tighter AGB sequences, providing a good probe of the intermediate age (0.5 - 5 Gyr) populations. We match observed CMDs with stellar evolution models to recover the SFH of KKH 98. On average, the galaxy has experienced relatively constant low-level star formation (5 x 10^-4 Mo yr^-1) for much of cosmic time. Except for the youngest main sequence populations (age < 0.1 Gyr), which are typically fainter than the AO data flux limit, the SFH estimated from the the 592 IR bright stars is a reasonable match to that derived from the much larger optical data set. Differences between the optical and IR derived SFHs for 0.1 - 1 Gyr populations suggest that current stellar evolution models may be over-producing the AGB by as much as a factor of three in this galaxy. At the depth of the AO data, the IR luminous stars are not crowded. Therefore these techniques can potentially be used to determine the stellar populations of galaxies at significantly further distances.

One important frontier for astronomical adaptive optics (AO) involves methods such as Multi-Object AO and Multi-Conjugate AO that have the potential to give a significantly larger field of view than conventional AO techniques. A second key emphasis over the next decade will be to push astronomical AO to visible wavelengths. We have conducted the first laboratory simulations of wide-field, laser guide star adaptive optics at visible wavelengths on a 10-meter-class telescope. These experiments, utilizing the UCO/Lick Observatory's Multi-Object / Laser Tomography Adaptive Optics (MOAO/LTAO) testbed, demonstrate new techniques in wavefront sensing and control that are crucial to future on-sky MOAO systems. We (1) test and confirm the feasibility of highly accurate atmospheric tomography with laser guide stars, (2) demonstrate key innovations allowing open-loop operation of Shack-Hartmann wavefront sensors (with errors of ~30 nm) as will be needed for MOAO, and (3) build a complete error budget model describing system performance. The AO system maintains a performance of 32.4% Strehl on-axis, with 24.5% and 22.6% at 10" and 15", respectively, at a science wavelength of 710 nm (R-band) over the equivalent of 0.8 seconds of simulation. The MOAO-corrected field of view is ~25 times larger in area than that limited by anisoplanatism at R-band. Our error budget is composed of terms verified through independent, empirical experiments. Error terms arising from calibration inaccuracies and optical drift are comparable in magnitude to traditional terms like fitting error and tomographic error. This makes a strong case for implementing additional calibration facilities in future AO systems, including accelerometers on powered optics, 3D turbulators, telescope and LGS simulators, and external calibration ports for deformable mirrors.

We present a pilot study of the stellar populations of 8 AGN hosts at z~1 and compare to (1) lower redshift samples and (2) a sample of nonactive galaxies of similar redshift. We utilize K' images in the GOODS South field obtained with the laser guide star adaptive optics (LGSAO) system at Keck Observatory. We combine this K' data with B, V, i, and z imaging from the ACS on HST to give multi-color photometry at a matched spatial resolution better than 100 mas in all bands. The hosts harbor AGN as inferred from their high X-ray luminosities (L_X > 10^42 ergs/s) or mid-IR colors. We find a correlation between the presence of younger stellar populations and the strength of the AGN, as measured with [OIII] line luminosity or X-ray (2-10 keV) luminosity. This finding is consistent with similar studies at lower redshift. Of the three Type II galaxies, two are disk galaxies and one is of irregular type, while in the Type I sample there only one disk-like source and four sources with smooth, elliptical/spheroidal morphologies. In addition, the mid-IR SEDs of the strong Type II AGN indicate that they are excited to LIRG (Luminous InfraRed Galaxy) status via galactic starbursting, while the strong Type I AGN are excited to LIRG status via hot dust surrounding the central AGN. This supports the notion that the obscured nature of Type II AGN at z~1 is connected with global starbursting and that they may be extincted by kpc-scale dusty features that are byproducts of this starbursting.

We combine high-resolution images in four optical/infra-red bands, obtained with the laser guide star adaptive optics system on the Keck Telescope and with the Hubble Space Telescope, to study the gravitational lens system SDSSJ0737+3216 (lens redshift 0.3223, source redshift 0.5812). We show that (under favorable observing conditions) ground-based images are comparable to those obtained with HST in terms of precision in the determination of the parameters of both the lens mass distribution and the background source. We also quantify the systematic errors associated with both the incomplete knowledge of the PSF, and the uncertain process of lens galaxy light removal, and find that similar accuracy can be achieved with Keck LGSAO as with HST. We then exploit this well-calibrated combination of optical and gravitational telescopes to perform a multi-wavelength study of the source galaxy at 0.01" effective resolution. We find the Sersic index to be indicative of a disk-like object, but the measured half-light radius (0.59+-0.007+-0.1 kpc) and stellar mass (2.0+-1.0+-0.8e9Msun) place it more than three sigma away from the local disk size-mass relation. The SDSSJ0737+3216 source has the characteristics of the most compact faint blue galaxies studied, and has comparable size and mass to dwarf early-type galaxies in the local universe. With the aid of gravitational telescopes to measure individual objects' brightness profiles to 10% accuracy, the study of the high-redshift size-mass relation may be extended by an order of magnitude or more beyond existing surveys at the low-mass end, thus providing a new observational test of galaxy formation models.

We report high-precision atmospheric parameters for 1907 stars in the N2K low-resolution spectroscopic survey, designed to identify metal-rich FGK dwarfs likely to harbor detectable planets. 284 of these stars are in the ideal temperature range for planet searches, T_eff <= 6000K, and have a 10% or greater probability of hosting planets based on their metallicities. The stars in the low-resolution spectroscopic survey should eventually yield >60 new planets, including 8-9 hot Jupiters. Short-period planets have already been discovered orbiting the survey targets HIP 14801 and HD 149143.

We have created a framework to facilitate the construction of specialized target lists for radial velocity surveys that are biased toward stars that (1) possess planets and (2) are easiest to observe with current detection techniques. We use a procedure that uniformly estimates fundamental stellar properties of Tycho 2 stars, with errors, using spline functions of broadband photometry and proper motion found in Hipparcos/Tycho 2 and 2MASS. We provide estimates of temperature and distance for 2.4 million Tycho 2 stars that lack trigonometric distances. For stars that appear to be FGK dwarfs according to estimated temperature and absolute magnitude, we also derive [Fe/H] and identify unresolved binary systems with mass ratios between 1.25 and 3. Our spline function models are trained on the unique Valenti & Fischer (2005) set, composed of 1000 dwarfs with precise stellar parameters estimated from HIRES spectroscopy. For FGK dwarfs with V photometric error less than 0.05 magnitudes, or V < 9, our temperature model gives a one-sigma error of +58.7/-65.9 K and our metallicity model gives a one-sigma error of +0.13/-0.14 dex. Our estimates of distance and spectral type enable us to isolate 354,822 Tycho 2 dwarfs, 321,996 of which are absent from Hipparcos, with giant and subgiant contamination at 2.6% and 7.2%, respectively. 2,500 of these FGK dwarfs are bright (V < 9.0) and metal-rich ([Fe/H] > 0.2). Our metallicity estimates have been used to identify targets for N2K (Fischer et al. 2005), a large-scale radial velocity search for Hot Jupiters, which has published the detection of 4 Hot Jupiters with one transit. The broadband filtering outlined here is the first screening tier for N2K; the second tier is a low-resolution spectroscopy program headed by S.E. Robinson (astro-ph/0510150).

We present empirical calibrations that estimate stellar metallicity, effective temperature and surface gravity as a function of Lick/IDS indices. These calibrations have been derived from a training set of 261 stars for which (1) high-precision measurements of [Fe/H], T_eff and log g have been made using spectral-synthesis analysis of HIRES spectra, and (2) Lick indices have also been measured. Our [Fe/H] calibration, which has precision 0.07 dex, has identified a number of bright (V < 9) metal-rich stars which are now being screened for hot Jupiter-type planets. Using the Yonsei-Yale stellar models, we show that the calibrations provide distance estimates accurate to 20% for nearby stars. This paper outlines the second tier of the screening of planet-search targets by the N2K Consortium, a project designed to identify the stars most likely to harbor extrasolar planets. Discoveries by the N2K Consortium include the transiting hot Saturn HD 149026 b (Sato et al. 2005, astro-ph/0507009) and HD 88133 b (Fischer et al. 2005). See Ammons et al. (2005, In Press) for a description of the first tier of N2K metallicity screening, calibrations using broadband photometry.