Comparing chemical abundances of a planet and the host star reveals the origin and formation path. Stellar abundance is measured with high-resolution spectroscopy. Planet abundance, on the other hand, is usually inferred from low-resolution data. For directly imaged exoplanets, the data are available from a slew of high-contrast imaging/spectroscopy instruments. Here, we study the chemical abundance of HR 8799 and its planet c. We measure stellar abundance using LBT/PEPSI (R=120,000) and archival HARPS data: stellar [C/H], [O/H], and C/O are 0.11$\pm$0.12, 0.12$\pm$0.14, and 0.54$^{+0.12}_{-0.09}$, all consistent with solar values. We conduct atmospheric retrieval using newly obtained Subaru/CHARIS data together with archival Gemini/GPI and Keck/OSIRIS data. We model the planet spectrum with petitRADTRANS and conduct retrieval using PyMultiNest. Retrieved planetary abundance can vary by $\sim$0.5 dex, from sub-stellar to stellar C and O abundances. The variation depends on whether strong priors are chosen to ensure a reasonable planet mass. Moreover, comparison with previous works also reveals inconsistency in abundance measurements. We discuss potential issues that can cause the inconsistency, e.g., systematics in individual data sets and different assumptions in the physics and chemistry in retrieval. We conclude that no robust retrieval can be obtained unless the issues are fully resolved.

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.

Age is a fundamental parameter of stars, yet in many cases ages of individual stars are presented without robust estimates of the uncertainty. We have developed a Bayesian framework, BAFFLES, to produce the age posterior for a star from its calcium emission strength (log($R'_{HK}$)) or lithium abundance (Li EW) and $B-V$ color. We empirically determine the likelihood functions for calcium and lithium as functions of age from literature measurements of stars in benchmark clusters with well-determined ages. We use a uniform prior on age which reflects a uniform star formation rate. The age posteriors we derive for several test cases are consistent with literature ages found from other methods. BAFFLES represents a robust method to determine the age posterior probability distribution for any field star with $0.45 \leq B-V \leq 0.9$ and a measurement of $R'_{HK}$ and/or $0.35 \leq B-V \leq 1.9$ and measured Li EW. We compile colors, $R'_{HK}$, and Li EW from over 2630 nearby field stars from the literature and present the derived BAFFLES age posterior for each star.

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models.

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 a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.

Ground-based direct imaging surveys like the Gemini Planet Imager Exoplanet Survey (GPIES) rely on Adaptive Optics (AO) systems to image and characterize exoplanets that are up to a million times fainter than their host stars. One factor that can reduce AO performance is turbulence induced by temperature differences in the instrument's immediate surroundings (e.g.: "dome seeing" or "mirror seeing"). In this analysis we use science observations, AO telemetry, and environmental data from September 2014 to February 2017 of the GPIES campaign to quantify the effects of "mirror seeing" on the performance of the GPI instrument. We show that GPI performance is optimal when the primary mirror (M1) is in equilibrium with the outside air temperature. We then examine the characteristics of mirror seeing by calculating the power spectral densities (PSD) of spatial and temporal Fourier modes. Inside the inertial range of the PSDs, we find that the spatial PSD amplitude increases when M1 is out of equilibrium and that the integrated turbulence may exhibit deviations from Kolmogorov atmospheric turbulence models and from the 1-layer frozen flow model. We conclude with an assessment of the current temperature control and ventilation strategy at Gemini South.

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.

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 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}$).

High-contrast medium resolution spectroscopy has been used to detect molecules such as water and carbon monoxide in the atmospheres of gas giant exoplanets. In this work, we show how it can be used to derive radial velocity (RV) measurements of directly imaged exoplanets. Improving upon the traditional cross-correlation technique, we develop a new likelihood based on joint forward modelling of the planetary signal and the starlight background (i.e., speckles). After marginalizing over the starlight model, we infer the barycentric RV of HR 8799 b and c in 2010 yielding -9.2 +- 0.5 km/s and -11.6 +- 0.5 km/s respectively. These RV measurements help to constrain the 3D orientation of the orbit of the planet by resolving the degeneracy in the longitude of ascending node. Assuming coplanar orbits for HR 8799 b and c, but not including d and e, we estimate \Omega = 89 (+27,-17) deg and i = 20.8 (4.5,-3.7) deg.

We present constraints on the presence of substellar companions to the nearby ($d\sim7.6$ pc), young ($440\pm40$ Myr) K4Ve star TW Piscis Austrini, the wide ($\sim$0.3 pc) companion to the A4V star Fomalhaut. We combined absolute astrometry from Hipparcos and Gaia with literature radial velocity measurements and dedicated high-contrast imaging observations obtained with Keck/NIRC2 to achieve sensitivity to brown dwarf and planetary-mass companions ($\gtrsim2$ M$_{\rm Jup}$) over many decades of orbital period ($\lesssim 10^3$ yr). The significant astrometric acceleration measured between the Hipparcos and Gaia catalogues, reported previously in the literature, cannot be explained by the orbital motion of TW PsA around the barycenter of the Fomalhaut triple system. Instead, we find that it is consistent with the reflex motion induced by an orbiting substellar companion. The combination of astrometry, radial velocities, and a deep $L^{\prime}$ imaging dataset leads to a constraint on the companion mass of $1.2_{-0.6}^{+0.7}$ M$_{\rm Jup}$. However, the period of the companion is poorly constrained, with a highly multi-modal period posterior distribution due to aliasing with the 24.25 year baseline between Hipparcos and Gaia. If confirmed through continued astrometric or spectroscopic monitoring, or via direct detection, the companion to TW PsA would represent a choice target for detailed atmospheric characterization with high-contrast instruments on the upcoming James Webb Space Telescope and Wide Field Infrared Survey Telescope.

An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay, and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low towards the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations.

We present modeled detection limits of the Gemini Planet Imager (GPI) and the Wide-Field Infrared Space Telescope (WFIRST) to an optical and infrared laser which could be used by an extraterrestrial civilization to signal their presence. GPI and WFIRST could utilize a coronagraph to search for extraterrestrial intelligence (SETI) in the present and future. We use archival data for GPI stars and simulated WFIRST observations to find the detectable flux ratio of a laser signal to residual scattered starlight around the target star. This flux ratio is then converted to detectable power as a function of distance from the parent star. For GPI, we assume a monochromatic laser wavelength of 1.55 $\mu$m, and a wavelength of 575 nm for WFIRST. We assume the lasers are projected through a 10-m aperture, and that the intensity of the laser beam follows a Gaussian profile. Our analysis is performed on 6 stars with spectral types later than F within 20 pc (with an emphasis on solar analogs at different distances). The most notable result is the detection limit for $\tau$ Ceti, a G5V star with four known exoplanets, two of those within the habitable zone (HZ). The result shows that a 24 kW laser is detectable from $\tau$ Ceti from outside of the HZ with GPI and a 7.3 W laser is detectable from within $\tau$ Ceti's HZ by WFIRST.

Based on two decades of radial velocity (RV) observations using Keck/HIRES and McDonald/Tull, and more recent observations using the Automated Planet Finder, we found that the nearby star HR 5183 (HD 120066) hosts a 3$M_J$ minimum mass planet with an orbital period of $74^{+43}_{-22}$ years. The orbit is highly eccentric (e$\simeq$0.84), shuttling the planet from within the orbit of Jupiter to beyond the orbit of Neptune. Our careful survey design enabled high cadence observations before, during, and after the planet's periastron passage, yielding precise orbital parameter constraints. We searched for stellar or planetary companions that could have excited the planet's eccentricity, but found no candidates, potentially implying that the perturber was ejected from the system. We did identify a bound stellar companion more than 15,000 au from the primary, but reasoned that it is currently too widely separated to have an appreciable effect on HR 5183 b. Because HR 5183 b's wide orbit takes it more than 30 au (1") from its star, we also explored the potential of complimentary studies with direct imaging or stellar astrometry. We found that a Gaia detection is very likely, and that imaging at 10 $\mu$m is a promising avenue. This discovery highlights the value of long-baseline RV surveys for discovering and characterizing long-period, eccentric Jovian planets. This population may offer important insights into the dynamical evolution of planetary systems containing multiple massive planets.

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.

Young exoplanets are snapshots of the planetary evolution process. Planets that orbit stars in young associations are particularly important because the age of the planetary system is well constrained. We present the discovery of a transiting planet larger than Neptune but smaller than Saturn in the 45 Myr Tucana-Horologium young moving group. The host star is a visual binary, and our follow-up observations demonstrate that the planet orbits the G6V primary component, DS Tuc A (HD 222259A, TIC 410214986). We first identified transits using photometry from the Transiting Exoplanet Survey Satellite (TESS; alerted as TOI 200.01). We validated the planet and improved the stellar parameters using a suite of new and archival data, including spectra from SOAR/Goodman, SALT/HRS and LCO/NRES; transit photometry from Spitzer; and deep adaptive optics imaging from Gemini/GPI. No additional stellar or planetary signals are seen in the data. We measured the planetary parameters by simultaneously modeling the photometry with a transit model and a Gaussian process to account for stellar variability. We determined that the planetary radius is $5.70\pm0.17$ Earth radii and that the orbital period is 8.1 days. The inclination angles of the host star's spin axis, the planet's orbital axis, and the visual binary's orbital axis are aligned within 15 degrees to within the uncertainties of the relevant data. DS Tuc Ab is bright enough (V=8.5) for detailed characterization using radial velocities and transmission spectroscopy.

All four circumbinary (CB) protoplanetary disks orbiting short-period ($P < 20$ day) double-lined spectroscopic binaries (SB2s)---a group that includes UZ Tau E, for which we present new ALMA data---exhibit sky-plane inclinations $i_{\rm disk}$ which match, to within a few degrees, the sky-plane inclinations $i_\star$ of their stellar hosts. Although for these systems the true mutual inclinations $\theta$ between disk and binary cannot be directly measured because relative nodal angles are unknown, the near-coincidence of $i_{\rm disk}$ and $i_\star$ suggests that $\theta$ is small for these most compact of systems. We confirm this hypothesis using a hierarchical Bayesian analysis, showing that 68% of CB disks around short-period SB2s have $\theta < 3.0^\circ$. Near co-planarity of CB disks implies near co-planarity of CB planets discovered by Kepler, which in turn implies that the occurrence rate of close-in CB planets is similar to that around single stars. By contrast, at longer periods ranging from $30-10^5$ days (where the nodal degeneracy can be broken via, e.g., binary astrometry), CB disks exhibit a wide range of mutual inclinations, from co-planar to polar. Many of these long-period binaries are eccentric, as their component stars are too far separated to be tidally circularized. We discuss how theories of binary formation and disk-binary gravitational interactions can accommodate all these observations.

Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs

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.

Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

The Coronagraph Instrument (CGI) for NASA's Wide Field Infrared Survey Telescope (WFIRST) will constitute a dramatic step forward for high-contrast imaging, integral field spectroscopy, and polarimetry of exoplanets and circumstellar disks, aiming to improve upon the sensitivity of current ground-based direct imaging facilities by 2-3 orders of magnitude. Furthermore, CGI will serve as a pathfinder for future exo-Earth imaging and characterization missions by demonstrating wavefront control, coronagraphy, and spectral retrieval in a new contrast regime, and by validating instrument and telescope models at unprecedented levels of precision. To achieve this jump in performance, it is critical to draw on the experience of ground-based high-contrast facilities. We discuss several areas of relevant commonalities, including: wavefront control, post-processing of integral field unit data, and calibration and observing strategies.

Rigorously quantifying the information in high contrast imaging data is important for informing follow-up strategies to confirm the substellar nature of a point source, constraining theoretical models of planet-disk interactions, and deriving planet occurrence rates. However, within the exoplanet direct imaging community, non-detections have almost exclusively been defined using a frequentist detection threshold (i.e. contrast curve) and associated completeness. This can lead to conceptual inconsistencies when included in a Bayesian framework. A Bayesian upper limit is such that the true value of a parameter lies below this limit with a certain probability. The associated probability is the integral of the posterior distribution with the upper limit as the upper bound. In summary, a frequentist upper limit is a statement about the detectability of planets while a Bayesian upper limit is a statement about the probability of a parameter to lie in an interval given the data. The latter is therefore better suited for rejecting hypotheses or theoretical models based on their predictions. In this work we emphasize that Bayesian statistics and upper limits are more easily interpreted and typically more constraining than the frequentist approach. We illustrate the use of Bayesian analysis in two different cases: 1) with a known planet location where we also propose to use model comparison to constrain the astrophysical nature of the point source and 2) gap-carving planets in TW Hya. To finish, we also mention the problem of combining radial velocity and direct imaging observations.

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.

One of the scientific goals of the Large Synoptic Survey Telescope (LSST) is to measure the evolution of dark energy by measuring subtle distortions of galaxy shapes due to weak gravitational lensing caused by the evolving dark matter distribution. Understanding the point spread function (PSF) for LSST is a crucial step to accurate measurements of weak gravitational lensing. Atmospheric contributions dominate the LSST PSF. Simulations of Kolmogorov turbulence models are commonly used to characterize and correct for these atmospheric effects. In order to validate these simulations, we compare the predicted atmospheric behavior to empirical data. The simulations are carried out in GalSim, an open-source software package for simulating images of astronomical objects and PSFs. Atmospheric simulations are run by generating large phase screens at varying altitude and evolving them over long time scales. We compare the turbulence strength and temporal behavior of atmospheres generated from simulations to those from reconstructed telemetry data from the Gemini Planet Imager (GPI). GPI captures a range of spatial frequencies by sampling the atmosphere with 18-cm subapertures. The LSST weak lensing analysis will measure correlations of galaxy ellipticity, requiring very accurate knowledge of the magnitude and correlations of PSF shape parameters. Following from the first analysis, we use simulations and sequential short exposure observations from the Differential Speckle Survey Instrument (DSSI) to study the behavior of PSF parameters - e.g., ellipticity and size - as a function of exposure time. These studies could help inform discussions of possible variable exposure times for LSST visits for example, to provide more uniform depth of visits.

Atmospheric characterization of directly imaged planets has thus far been limited to ground-based observations of young, self-luminous, Jovian planets. Near-term space- and ground- based facilities like \emph{WFIRST} and ELTs will be able to directly image mature Jovian planets in reflected light, a critical step in support of future facilities that aim to directly image terrestrial planets in reflected light (e.g. HabEx, LUVOIR). These future facilities are considering the use of photometry to classify planets. Here, we investigate the intricacies of using colors to classify gas-giant planets by analyzing a grid of 9,120 theoretical reflected light spectra spread across different metallicities, pressure-temperature profiles, cloud properties, and phase angles. We determine how correlated these planet parameters are with the colors in the \emph{WFIRST} photometric bins and other photometric bins proposed in the literature. Then we outline under what conditions giant planet populations can be classified using several supervised multivariate classification algorithms. We find that giant planets imaged in reflected light can be classified by metallicity with an accuracy of $>$90\% if they are \emph{a prior} known to not have significant cloud coverage in the visible part of the atmosphere, and at least 3 filter observations are available. If the presence of clouds is not known \emph{a priori}, directly imaged planets can be more accurately classified by their cloud properties, as oppposed to metallicity or temperature. Furthermore, we are able to distinguish between cloudy and cloud-free populations with $>$90\% accuracy with 3 filter observations. Our statistical pipeline is available on GitHub and can be extended to optimize science yield of future mission concepts.

The Gemini Planet Imager (GPI) is the dedicated high-contrast imaging facility, located on Gemini South, designed for the direct detection and characterization of young Jupiter mass exoplanets. In 2019, Gemini is considering moving GPI from Gemini South to Gemini North. Analysis of GPI's as-built performance has highlighted several key areas of improvement to its detection capabilities while leveraging its current capabilities as a facility class instrument. We present the proposed upgrades which include a pyramid wavefront sensor, broadband low spectral resolution prisms and new apodized-pupil Lyot coronagraph designs all of which will enhance the current science capabilities while enabling new science programs.

The Gemini Planet Imager (GPI), a coronagraphic adaptive optics instrument designed for spectroscopy of extrasolar planets, had first light in 2013[13]. After five years, GPI has observed more than 500 stars, producing an extensive library of science images and associated telemetry that can be analyzed to determine performance predictors. We will present a summary of on-sky performance and lessons learned. The two most significant factors determining bright star contrast performance are atmospheric coherence time and the presence of dome seeing. With a possible move to Gemini North, we are planning potential upgrades including a pyramid-sensor based AO system with predictive control; we will summarize upgrade options and the science they would enable.

The Gemini Planet Imager (GPI) is a near-infrared instrument that uses Adaptive Optics (AO), a coronagraph, and advanced data processing techniques to achieve very high contrast images of exoplanets. The GPI Exoplanet Survey (GPIES) is a 600 stars campaign aiming at detecting and characterizing young, massive and self-luminous exoplanets at large orbital distances >5 au. Science observations are taken simultaneously with environmental data revealing information about the turbulence in the telescope environment as well as limitations of GPI's AO system. Previous work has shown that the timescale of the turbulence, Tau0, is a strong predictor of AO performance, however an analysis of the dome turbulence on AO performance has not been done before. Here, we study correlations between image contrast and residual wavefront error (WFE) with temperature measurements from multiple locations inside and outside the dome. Our analysis revealed GPI's performance is most correlated with the temperature difference between the primary mirror of the telescope and the outside air. We also assess the impact of the current temperature control and ventilation strategy at Gemini South (GS).

A semi analytic framework for simulating the effects of atmospheric seeing in Adaptive Optics systems on an 8-m telescope is developed with the intention of understanding the origin of the wind-butterfly, a characteristic two-lobed halo in the PSF of AO imaging. Simulations show that errors in the compensated phase on the aperture due to servo-lag have preferential direction orthogonal to the direction of wind propagation which, when Fourier Transformed into the image plane, appear with their characteristic lemniscate shape along the wind direction. We develop a metric to quantify the effect of this aberration with the fractional standard deviation in an annulus centered around the PSF, and use telescope pointing to correlate this effect with data from an atmospheric models, the NOAA GFS. Our results show that the jet stream at altitudes of 100-200 hPa (equivalently 10-15 km above sea level) is highly correlated (13.2$\sigma$) with the strong butterfly, while the ground wind and other layers are more or less uncorrelated.

The coronagraphic instrument (CGI) on the Wide-Field Infrared Survey Telescope (WFIRST) will demonstrate technologies and methods for high-contrast direct imaging and spectroscopy of exoplanet systems in reflected light, including polarimetry of circumstellar disks. The WFIRST management and CGI engineering and science investigation teams have developed requirements for the instrument, motivated by the objectives and technology development needs of potential future flagship exoplanet characterization missions such as the NASA Habitable Exoplanet Imaging Mission (HabEx) and the Large UV/Optical/IR Surveyor (LUVOIR). The requirements have been refined to support recommendations from the WFIRST Independent External Technical/Management/Cost Review (WIETR) that the WFIRST CGI be classified as a technology demonstration instrument instead of a science instrument. This paper provides a description of how the CGI requirements flow from the top of the overall WFIRST mission structure through the Level 2 requirements, where the focus here is on capturing the detailed context and rationales for the CGI Level 2 requirements. The WFIRST requirements flow starts with the top Program Level Requirements Appendix (PLRA), which contains both high-level mission objectives as well as the CGI-specific baseline technical and data requirements (BTR and BDR, respectively)... We also present the process and collaborative tools used in the L2 requirements development and management, including the collection and organization of science inputs, an open-source approach to managing the requirements database, and automating documentation. The tools created for the CGI L2 requirements have the potential to improve the design and planning of other projects, streamlining requirement management and maintenance. [Abstract Abbreviated]

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.

We explore KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR 8799, using PyKLIP and show algorithm stability with varying KLIP parameters. We report new and re-reduced spectrophotometry of HR 8799 c, d, and e in H & K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting simulated sources and recovering them over a range of parameters. The K1/K2 spectra for HR 8799 c and d are similar to previously published results from the same dataset. We also present a K band spectrum of HR 8799 e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We show that HR 8799 c and d show significant differences in their H & K spectra, but do not find any conclusive differences between d and e or c and e, likely due to large error bars in the recovered spectrum of e. Compared to M, L, and T-type field brown dwarfs, all three planets are most consistent with mid and late L spectral types. All objects are consistent with low gravity but a lack of standard spectra for low gravity limit the ability to fit the best spectral type. We discuss how dedicated modeling efforts can better fit HR 8799 planets' near-IR flux and discuss how differences between the properties of these planets can be further explored.

Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution ($R$). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at $R = 70$ and $R = 140$ from 0.4-1.0 $\mu$m with SNR $ = 5, 10, 15, 20$ at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48-0.6 $\mu$m and $R = 50$ spectroscopy from 0.6-0.97 $\mu$m. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least $R = 70$ / SNR$\ = 15$, or $R = 140$ / SNR$\ = 10$ for improved detections. Meaningful constraints are only achieved with $R = 140$ / SNR$\ = 20$ data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass.

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.

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 procedure for calculating expected exoplanet imaging yields, which explicitly separates the effects of instrument performance from assumptions of planet distributions. This `depth of search' approach allows for fast recalculation of yield values for variations in instrument parameters. We also describe a new target star selection metric with no dependence on an assumed planet population that can be used as a proxy for single-visit completeness. This approach allows for the recovery of the total mission completeness via convolution of the depth of search grid with an equivalent grid of assumed occurrence rates and integration over the part of the grid representing the population of interest (e.g., Earth-like planets on habitable zone orbits, etc.). In this work, we discuss the practical details of calculating the depth of search and present results of such calculations for one design iteration of the WFIRST coronagraphs.

We present a new matched filter algorithm for direct detection of point sources in the immediate vicinity of bright stars. The stellar Point Spread Function (PSF) is first subtracted using a Karhunen-Lo\'eve Image Processing (KLIP) algorithm with Angular and Spectral Differential Imaging (ADI and SDI). The KLIP-induced distortion of the astrophysical signal is included in the matched filter template by computing a forward model of the PSF at every position in the image. To optimize the performance of the algorithm, we conduct extensive planet injection and recovery tests and tune the exoplanet spectra template and KLIP reduction aggressiveness to maximize the Signal-to-Noise Ratio (SNR) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal SNR loss. We also developed a complete pipeline for the automated detection of point source candidates, the calculation of Receiver Operating Characteristics (ROC), false positives based contrast curves, and completeness contours. We process in a uniform manner more than 330 datasets from the Gemini Planet Imager Exoplanet Survey (GPIES) and assess GPI typical sensitivity as a function of the star and the hypothetical companion spectral type. This work allows for the first time a comparison of different detection algorithms at a survey scale accounting for both planet completeness and false positive rate. We show that the new forward model matched filter allows the detection of $50\%$ fainter objects than a conventional cross-correlation technique with a Gaussian PSF template for the same false positive rate.

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 describe a Bayesian rejection sampling algorithm designed to efficiently compute posterior distributions of orbital elements for data covering short fractions of long-period exoplanet orbits. Our implementation of this method, Orbits for the Impatient (OFTI), converges up to several orders of magnitude faster than two implementations of MCMC in this regime. We illustrate the efficiency of our approach by showing that OFTI calculates accurate posteriors for all existing astrometry of the exoplanet 51 Eri b up to 100 times faster than a Metropolis-Hastings MCMC. We demonstrate the accuracy of OFTI by comparing our results for several orbiting systems with those of various MCMC implementations, finding the output posteriors to be identical within shot noise. We also describe how our algorithm was used to successfully predict the location of 51 Eri b six months in the future based on less than three months of astrometry. Finally, we apply OFTI to ten long-period exoplanets and brown dwarfs, all but one of which have been monitored over less than 3% of their orbits, producing fits to their orbits from astrometric records in the literature.

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.

Using the Gemini Planet Imager (GPI) located at Gemini South, we measured the near-infrared (1.0-2.4 micron) spectrum of the planetary companion to the nearby, young star $\beta$ Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of $12.9\pm0.2$ $\mathcal{M}_\mathrm{Jup}$, an effective temperature of $1724\pm15$ K, a radius of $1.46\pm0.01$ $\mathcal{R}_\mathrm{Jup}$, and a surface gravity of $\log g = 4.18\pm0.01$ [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of $1700-1800$ K and a surface gravity of $\log g = 3.5$-$4.0$ [dex] depending upon model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of $\beta$ Pic b best matches a low-surface gravity L2$\pm$1 brown dwarf. Finally comparing the spectrum to field brown dwarfs we find the the spectrum best matches 2MASS J04062677-381210 and 2MASS J03552337+1133437.

We present new spatially resolved astrometry and photometry from the Gemini Planet Imager of the inner binary of the young multiple star system V343 Normae, which is a member of the beta Pictoris moving group. V343 Normae comprises a K0 and mid-M star in a ~4.5 year orbit (AaAb) and a wide 10" M5 companion (B). By combining these data with archival astrometry and radial velocities we fit the orbit and measure individual masses for both components of M_Aa = 1.10 +/- 0.10 M_sun and M_Ab = 0.290 +/- 0.018 M_sun. Comparing to theoretical isochrones, we find good agreement for the measured masses and JHK band magnitudes of the two components consistent with the age of the beta Pic moving group. We derive a model-dependent age for the beta Pic moving group of 26 +/- 3 Myr by combining our results for V343 Normae with literature measurements for GJ 3305, which is another group member with resolved binary components and dynamical masses.

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 Gemini Planet Imager (GPI) is a next-generation high-contrast imager built for the Gemini Observatory. The GPI exoplanet survey (GPIES) consortium is made up of 102 researchers from 28 institutions in North and South America and Europe. In November 2014, we launched a search for young Jovian planets and debris disks. In this paper, we discuss how we have coordinated the work done by this large team to improve the technical and scientific productivity of the campaign, and describe lessons we have learned that could be useful for future instrumentation-based astronomical surveys. The success of GPIES lies mostly on its decentralized structure, clear definition of policies that are signed by each member, and the heavy use of modern tools for communicating, exchanging information, and processing data.

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.

Transition disks, protoplanetary disks with inner clearings, are promising objects in which to directly image forming planets. The high contrast imaging technique of non-redundant masking is well posed to detect planetary mass companions at several to tens of AU in nearby transition disks. We present non-redundant masking observations of the T Cha and LkCa 15 transition disks, both of which host posited sub-stellar mass companions. However, due to a loss of information intrinsic to the technique, observations of extended sources (e.g. scattered light from disks) can be misinterpreted as moving companions. We discuss tests to distinguish between these two scenarios, with applications to the T Cha and LkCa 15 observations. We argue that a static, forward-scattering disk can explain the T Cha data, while LkCa 15 is best explained by multiple orbiting companions.

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.

Radial velocity and transit methods are effective for the study of short orbital period exoplanets but they hardly probe objects at large separations for which direct imaging can be used. We carried out the international deep planet survey of 292 young nearby stars to search for giant exoplanets and determine their frequency. We developed a pipeline for a uniform processing of all the data that we have recorded with NIRC2/Keck II, NIRI/Gemini North, NICI/Gemini South, and NACO/VLT for 14 years. The pipeline first applies cosmetic corrections and then reduces the speckle intensity to enhance the contrast in the images. The main result of the international deep planet survey is the discovery of the HR 8799 exoplanets. We also detected 59 visual multiple systems including 16 new binary stars and 2 new triple stellar systems, as well as 2,279 point-like sources. We used Monte Carlo simulations and the Bayesian theorem to determine that 1.05[+2.80-0.70]% of stars harbor at least one giant planet between 0.5 and 14M_J and between 20 and 300 AU. This result is obtained assuming uniform distributions of planet masses and semi-major axes. If we consider power law distributions as measured for close-in planets instead, the derived frequency is 2.30[+5.95-1.55]%, recalling the strong impact of assumptions on Monte Carlo output distributions. We also find no evidence that the derived frequency depends on the mass of the hosting star, whereas it does for close-in planets. The international deep planet survey provides a database of confirmed background sources that may be useful for other exoplanet direct imaging surveys. It also puts new constraints on the number of stars with at least one giant planet reducing by a factor of two the frequencies derived by almost all previous works.

HD 142527A is one of the most studied Herbig Ae/Be stars with a transitional disk, as it has the largest imaged gap in any protoplanetary disk: the gas is cleared from 30 to 90 AU. The HD142527 system is also unique in that it has a stellar companion with a small mass compared to the mass of the primary star. This factor of ~20 in mass ratio between the two objects makes this binary system different from any other YSO. The HD142527 system could therefore provide a valuable test bed. This low-mass stellar object may be responsible for both the gap and dust trapping observed by ALMA at longer distances. We observed this system with the NACO and GPI instruments using the aperture masking technique. Aperture masking is ideal for providing high dynamic range even at very small angular separations. We present the spectral energy distribution for HD142527A and B. Brightness of the companion is now known from the R band up to the M' band. We also followed the orbital motion of HD 142527B over a period of more than two years. The SED of the companion is compatible with a T=3000+/-100K object in addition to a 1700K blackbody environment (likely a circus-secondary disk). From evolution models, we find that it is compatible with an object of mass 0.13+/-0.03Msun, radius 0.90+/-0.15Rsun, and age $1.0^{+1.0}_{-0.75}$Myr. This age is significantly younger than the age previously estimated for HD142527A. Computations to constrain the orbital parameters found a semi major axis of $140^{+120}_{-70}$mas, an eccentricity of 0.5+/-0.2, an inclination of 125+/-15 degrees, and a position angle of the right ascending node of -5+/-40 degrees. Inclination and position angle of the ascending node are in agreement with an orbit coplanar with the inner disk, not coplanar with the outer disk. Despite its high eccentricity, it is unlikely that HD142527B is responsible for truncating the inner edge of the outer disk.

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 new astrometric measurements from our ongoing monitoring campaign of the HR 8799 directly imaged planetary system. These new data points were obtained with NIRC2 on the W.M. Keck II 10 meter telescope between 2009 and 2014. In addition, we present updated astrometry from previously published observations in 2007 and 2008. All data were reduced using the SOSIE algorithm, which accounts for systematic biases present in previously published observations. This allows us to construct a self-consistent data set derived entirely from NIRC2 data alone. From this dataset, we detect acceleration for two of the planets (HR 8799b and e) at $>$3$\sigma$. We also assess possible orbital parameters for each of the four planets independently. We find no statistically significant difference in the allowed inclinations of the planets. Fitting the astrometry while forcing coplanarity also returns $\chi^2$ consistent to within 1$\sigma$ of the best fit values, suggesting that if inclination offsets of $\lesssim$20$^{o}$ are present, they are not detectable with current data. Our orbital fits also favor low eccentricities, consistent with predictions from dynamical modeling. We also find period distributions consistent to within 1$\sigma$ with a 1:2:4:8 resonance between all planets. This analysis demonstrates the importance of minimizing astrometric systematics when fitting for solutions to highly undersampled orbits.

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.

We present H and K band imaging polarimetry for the PDS 66 circumstellar disk obtained during the commissioning of the Gemini Planet Imager (GPI). Polarization images reveal a clear detection of the disk in to the 0.12'' inner working angle (IWA) in H band, almost 3 times as close to the star as the previous HST observations with NICMOS and STIS (0.35'' effective IWA). The centro-symmetric polarization vectors confirm that the bright inner disk detection is due to circumstellar scattered light. A more diffuse disk extends to a bright outer ring centered at 80 AU. We discuss several physical mechanisms capable of producing the observed ring + gap structure. GPI data confirm enhanced scattering on the East side of the disk which is inferred to be nearer to us. We also detect a lateral asymmetry in the South possibly due to shadowing from material within the inner working angle. This likely corresponds to a temporally variable azimuthal asymmetry observed in HST/STIS coronagraphic imaging.

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.

Exoplanet detections have revolutionized astronomy, offering new insights into solar system architecture and planet demographics. While nearly 1900 exoplanets have now been discovered and confirmed, none are still in the process of formation. Transition discs, protoplanetary disks with inner clearings best explained by the influence of accreting planets, are natural laboratories for the study of planet formation. Some transition discs show evidence for the presence of young planets in the form of disc asymmetries or infrared sources detected within their clearings, as in the case of LkCa 15. Attempts to observe directly signatures of accretion onto protoplanets have hitherto proven unsuccessful. Here we report adaptive optics observations of LkCa 15 that probe within the disc clearing. With accurate source positions over multiple epochs spanning 2009 - 2015, we infer the presence of multiple companions on Keplerian orbits. We directly detect H{\alpha} emission from the innermost companion, LkCa 15 b, evincing hot (~10,000 K) gas falling deep into the potential well of an accreting protoplanet.

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.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

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.

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 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 results from a Hubble Space Telescope (HST) program characterizing the atmospheres of the outer two planets, in the HR8799 system. The images were taken over 15 orbits in three near-infrared medium-band filters - F098M, F127M and F139M - using the Wide Field Camera 3. One of the three filters is sensitive to water absorption band inaccessible from ground-based observations, providing a unique probe of the thermal emission from the atmospheres of these young giant planets. The observations were taken at 30 different spacecraft rolls to enable angular differential imaging, and the full data set was analyzed with the Karhunen-Loeve Image Projection (KLIP) routine, an advanced image processing algorithm adapted to work with HST data. To achieve the required high contrast at sub arcsecond resolution, we utilized the pointing accuracy of HST in combination with an improved pipeline designed to combine the dithered, angular differential imaging data with an algorithm designed to both improve the image resolution and accurately measure the photometry. The results include F127M ($J$) detections of the outer planets, HR8799 b and c and the first detection of HR8799 b in the water-band (F139M) filter. The F127M photometry for HR8799 c agrees well with fitted atmospheric models resolving a long standing difficulty to model the near-IR flux for the planet consistently

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We are still in the early days of exoplanet discovery. Astronomers are beginning to model the atmospheres and interiors of exoplanets and have developed a deeper understanding of processes of planet formation and evolution. However, we have yet to map out the full complexity of multi-planet architectures or to detect Earth analogues around nearby stars. Reaching these ambitious goals will require further improvements in instrumentation and new analysis tools. In this chapter, we provide an overview of five observational techniques that are currently employed in the detection of exoplanets: optical and IR Doppler measurements, transit photometry, direct imaging, microlensing, and astrometry. We provide a basic description of how each of these techniques works and discuss forefront developments that will result in new discoveries. We also highlight the observational limitations and synergies of each method and their connections to future space missions.

We report the presence of scattered light from dust grains located in the giant planet formation region of the circumbinary disk orbiting the ~20-Myr-old close (~0.045 AU separation) binary system V4046 Sgr AB based on observations with the new Gemini Planet Imager (GPI) instrument. These GPI images probe to within ~7 AU of the central binary with linear spatial resolution of ~3 AU, and are thereby capable of revealing dust disk structure within a region corresponding to the giant planets in our solar system. The GPI imaging reveals a relatively narrow (FWHM ~10 AU) ring of polarized near-infrared flux whose brightness peaks at ~14 AU. This ~14 AU radius ring is surrounded by a fainter outer halo of scattered light extending to ~45 AU, which coincides with previously detected mm-wave thermal dust emission. The presence of small grains that efficiently scatter starlight well inside the mm-wavelength disk cavity supports current models of planet formation that suggest planet-disk interactions can generate pressure traps that impose strong radial variations in the particle size distribution throughout the disk.

This report describes the 2014 study by the Science Definition Team (SDT) of the Wide-Field Infrared Survey Telescope (WFIRST) mission. It is a space observatory that will address the most compelling scientific problems in dark energy, exoplanets and general astrophysics using a 2.4-m telescope with a wide-field infrared instrument and an optical coronagraph. The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Study Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality telescope assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012) and the previous WFIRST-2.4 DRM, reported by Spergel et. (2013). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The addition of an on-axis coronagraphic instrument to the baseline design enables imaging and spectroscopic studies of planets around nearby stars.

Absorption lines from water, methane and carbon monoxide are detected in the atmosphere of exoplanet HR8799b. A medium-resolution spectrum presented here shows well-resolved and easily identified spectral features from all three molecules across the K band. The majority of the lines are produced by CO and H2O, but several lines clearly belong to CH4. Comparisons between these data and atmosphere models covering a range of temperatures and gravities yield log mole fractions of H2O between -3.09 and -3.91, CO between -3.30 and -3.72 and CH4 between -5.06 and -5.85. More precise mole fractions are obtained for each temperature and gravity studied. A reanalysis of H-band data, previously obtained at similar spectral resolution, results in a nearly identical water abundance as determined from the K-band spectrum. The methane abundance is shown to be sensitive to vertical mixing and indicates an eddy diffusion coefficient in the range of 10^6 to 10^8 cm^2 s^-1, comparable to mixing in the deep troposphere of Jupiter. The model comparisons also indicate a C/O between ~ 0.58 and 0.7, encompassing previous estimates for a second planet in the same system, HR8799c. Super-stellar C/O could indicate planet formation by core-accretion, however, the range of possible C/O for these planets (and the star) is currently too large to comment strongly on planet formation. More precise values of the bulk properties (e.g., effective temperature and surface gravity) are needed for improved abundance estimates.

We report the discovery of a wide co-moving substellar companion to the nearby ($D=67.5\pm1.1$ pc) A3V star $\zeta$ Delphini based on imaging and follow-up spectroscopic observations obtained during the course of our Volume-limited A-Star (VAST) multiplicity survey. $\zeta$ Del was observed over a five-year baseline with adaptive optics, revealing the presence of a previously-unresolved companion with a proper motion consistent with that of the A-type primary. The age of the $\zeta$ Del system was estimated as $525\pm125$ Myr based on the position of the primary on the colour-magnitude and temperature-luminosity diagrams. Using intermediate-resolution near-infrared spectroscopy, the spectrum of $\zeta$ Del B is shown to be consistent with a mid-L dwarf (L$5\pm2$), at a temperature of $1650\pm200$ K. Combining the measured near-infrared magnitude of $\zeta$ Del B with the estimated temperature leads to a model-dependent mass estimate of $50\pm15$ M$_{\rm Jup}$, corresponding to a mass ratio of $q=0.019\pm0.006$. At a projected separation of $910\pm14$ au, $\zeta$ Del B is among the most widely-separated and extreme-mass ratio substellar companions to a main-sequence star resolved to-date, providing a rare empirical constraint of the formation of low-mass ratio companions at extremely wide separations.

During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.

The Gemini Planet Imager (GPI) entered on-sky commissioning and had its first-light at the Gemini South (GS) telescope in November 2013. GPI is an extreme adaptive optics (XAO), high-contrast imager and integral-field spectrograph dedicated to the direct detection of hot exo-planets down to a Jupiter mass. The performance of the apodized pupil Lyot coronagraph depends critically upon the residual wavefront error (design goal of 60 nm RMS with 5 mas RMS tip/tilt), and therefore is most sensitive to vibration (internal or external) of Gemini's instrument suite. Excess vibration can be mitigated by a variety of methods such as passive or active dampening at the instrument or telescope structure or Kalman filtering of specific frequencies with the AO control loop. Understanding the sources, magnitudes and impact of vibration is key to mitigation. This paper gives an overview of related investigations based on instrument data (GPI AO module) as well as external data from accelerometer sensors placed at different locations on the GS telescope structure. We report the status of related mitigation efforts, and present corresponding results.

The Gemini Planet Imager (GPI) entered on-sky commissioning phase, and had its First Light at the Gemini South telescope in November 2013. Meanwhile, the fast loops for atmospheric correction of the Extreme Adaptive Optics (XAO) system have been closed on many dozen stars at different magnitudes (I=4-8), elevation angles and a variety of seeing conditions, and a stable loop performance was achieved from the beginning. Ultimate contrast performance requires a very low residual wavefront error (design goal 60 nm RMS), and optimization of the planet finding instrument on different ends has just begun to deepen and widen its dark hole region. Laboratory raw contrast benchmarks are in the order of 10^-6 or smaller. In the telescope environment and in standard operations new challenges are faced (changing gravity, temperature, vibrations) that are tackled by a variety of techniques such as Kalman filtering, open-loop models to keep alignment to within 5 mas, speckle nulling, and a calibration unit (CAL). The CAL unit was especially designed by the Jet Propulsion Laboratory to control slowly varying wavefront errors at the focal plane of the apodized Lyot coronagraph by the means of two wavefront sensors: 1) a 7x7 low order Shack-Hartmann SH wavefront sensor (LOWFS), and 2) a special Mach-Zehnder interferometer for mid-order spatial frequencies (HOWFS) - atypical in that the beam is split in the focal plane via a pinhole but recombined in the pupil plane with a beamsplitter. The original design goal aimed for sensing and correcting on a level of a few nm which is extremely challenging in a telescope environment. This paper focuses on non-common path low order wavefront correction as achieved through the CAL unit on sky. We will present the obtained results as well as explain challenges that we are facing.

Using the recently installed Gemini Planet Imager (GPI), we have taken the first H-band spectrum of the planetary companion to the nearby young star beta Pictoris. GPI is designed to image and provide low-resolution spectra of Jupiter sized, self-luminous planetary companions around young nearby stars. These observations were taken covering the H-band (1.65 microns). The spectrum has a resolving power of $\sim$ 45 and demonstrates the distinctive triangular shape of a cool substellar object with low surface gravity. Using atmospheric models, we find an effective temperature of $1650 \pm 50$ K and a surface gravity of $\log(g) = 4.0 \pm 0.25$ (cgs units). These values agree well with predictions from planetary evolution models for a gas giant with mass between 10 and 12 $M_{\rm Jup}$ and age between 10 and 20 Myrs.

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

We present the first results from the polarimetry mode of the Gemini Planet Imager (GPI), which uses a new integral field polarimetry architecture to provide high contrast linear polarimetry with minimal systematic biases between the orthogonal polarizations. We describe the design, data reduction methods, and performance of polarimetry with GPI. Point spread function subtraction via differential polarimetry suppresses unpolarized starlight by a factor of over 100, and provides sensitivity to circumstellar dust reaching the photon noise limit for these observations. In the case of the circumstellar disk around HR 4796A, GPI's advanced adaptive optics system reveals the disk clearly even prior to PSF subtraction. In polarized light, the disk is seen all the way in to its semi-minor axis for the first time. The disk exhibits surprisingly strong asymmetry in polarized intensity, with the west side >9 times brighter than the east side despite the fact that the east side is slightly brighter in total intensity. Based on a synthesis of the total and polarized intensities, we now believe that the west side is closer to us, contrary to most prior interpretations. Forward scattering by relatively large silicate dust particles leads to the strong polarized intensity on the west side, and the ring must be slightly optically thick in order to explain the lower brightness in total intensity there. These findings suggest that the ring is geometrically narrow and dynamically cold, perhaps shepherded by larger bodies in the same manner as Saturn's F ring.

To directly image exoplanets and faint circumstellar disks, the noisy stellar halo must be suppressed to a high level. To achieve this feat, the angular differential imaging observing technique and the least-squares Locally Optimized Combination of Images (LOCI) algorithm have now become the standard in single band direct imaging observations and data reduction. With the development and commissioning of new high-order high-contrast adaptive optics equipped with integral field units, the image subtraction algorithm needs to be modified to allow the optimal use of polychromatic images, field-rotated images and archival data. A new algorithm, TLOCI (for Template LOCI), is designed to achieve this task by maximizing a companion signal-to-noise ratio instead of simply minimizing the noise as in the original LOCI algorithm. The TLOCI technique uses an input spectrum and template Point Spread Functions (PSFs, generated from unocculted and unsaturated stellar images) to optimize the reference image least-squares coefficients to minimize the planet self-subtraction, thus maximizing its throughput per wavelength, while simultaneously providing a maximum suppression of the speckle noise. The new algorithm has been developed using on-sky GPI data and has achieved impressive contrast. This paper presents the TLOCI algorithm, on-sky performance, and will discuss the challenges in recovering the planet spectrum with high fidelity.

We present on-sky polarimetric observations with the Gemini Planet Imager (GPI) obtained at straight Cassegrain focus on the Gemini South 8-m telescope. Observations of polarimetric calibrator stars, ranging from nearly unpolarized to strongly polarized, enable determination of the combined telescope and instrumental polarization. We find the conversion of Stokes $I$ to linear and circular instrumental polarization in the instrument frame to be $I \rightarrow (Q_{\rm IP}, U_{\rm IP}, P_{\rm IP}, V_{\rm IP}) = (-0.037 \pm 0.010\%, +0.4338 \pm 0.0075\%, 0.4354 \pm 0.0075\%, -6.64 \pm 0.56\%)$. Such precise measurement of instrumental polarization enables $\sim 0.1\%$ absolute accuracy in measurements of linear polarization, which together with GPI's high contrast will allow GPI to explore scattered light from circumstellar disk in unprecedented detail, conduct observations of a range of other astronomical bodies, and potentially even study polarized thermal emission from young exoplanets. Observations of unpolarized standard stars also let us quantify how well GPI's differential polarimetry mode can suppress the stellar PSF halo. We show that GPI polarimetry achieves cancellation of unpolarized starlight by factors of 100-200, reaching the photon noise limit for sensitivity to circumstellar scattered light for all but the smallest separations at which the calibration for instrumental polarization currently sets the limit.

The Gemini Planet Imager (GPI) is an "extreme" adaptive optics coronagraph system that is now on the Gemini South telescope in Chile. This instrument is composed of three different systems that historically have been separate instruments. These systems are the extreme Adaptive Optics system, with deformable mirrors, including a high-order 64x64 element MEMS system; the Science Instrument, which is a near-infrared integral field spectrograph; and the Calibration system, a precision IR wavefront sensor that also holds key coronagraph components. Each system coordinates actions that require precise timing. The observatory is responsible for starting these actions and has typically done this asynchronously across independent systems. Despite this complexity we strived to provide an interface that is as close to a one-button approach as possible. This paper will describe the sequencing of these systems both internally and externally through the observatory.

The Gemini Planet Imager is an extreme AO instrument with an integral field spectrograph (IFS) operating in Y, J, H, and K bands. Both the Gemini telescope and the GPI instrument are very complex systems. Our goal is that the combined telescope and instrument system may be run by one observer operating the instrument, and one operator controlling the telescope and the acquisition of light to the instrument. This requires a smooth integration between the two systems and easily operated control interfaces. We discuss the definition of the software and hardware interfaces, their implementation and testing, and the integration of the instrument with the telescope environment.

The Gemini Planet Imager (GPI) is a next-generation, facility instrument currently being commissioned at the Gemini South observatory. GPI combines an extreme adaptive optics system and integral field spectrograph (IFS) with an apodized-pupil Lyot coronagraph (APLC) producing an unprecedented capability for directly imaging and spectroscopically characterizing extrasolar planets. GPI's operating goal of $10^{-7}$ contrast requires very precise alignments between the various elements of the coronagraph (two pupil masks and one focal plane mask) and active control of the beam path throughout the instrument. Here, we describe the techniques used to automatically align GPI and maintain the alignment throughout the course of science observations. We discuss the particular challenges of maintaining precision alignments on a Cassegrain mounted instrument and strategies that we have developed that allow GPI to achieve high contrast even in poor seeing conditions.

The Gemini Planet Imager (GPI) is a complex optical system designed to directly detect the self-emission of young planets within two arcseconds of their host stars. After suppressing the starlight with an advanced AO system and apodized coronagraph, the dominant residual contamination in the focal plane are speckles from the atmosphere and optical surfaces. Since speckles are diffractive in nature their positions in the field are strongly wavelength dependent, while an actual companion planet will remain at fixed separation. By comparing multiple images at different wavelengths taken simultaneously, we can freeze the speckle pattern and extract the planet light adding an order of magnitude of contrast. To achieve a bandpass of 20%, sufficient to perform speckle suppression, and to observe the entire two arcsecond field of view at diffraction limited sampling, we designed and built an integral field spectrograph with extremely low wavefront error and almost no chromatic aberration. The spectrograph is fully cryogenic and operates in the wavelength range 1 to 2.4 microns with five selectable filters. A prism is used to produce a spectral resolution of 45 in the primary detection band and maintain high throughput. Based on the OSIRIS spectrograph at Keck, we selected to use a lenslet-based spectrograph to achieve an rms wavefront error of approximately 25 nm. Over 36,000 spectra are taken simultaneously and reassembled into image cubes that have roughly 192x192 spatial elements and contain between 11 and 20 spectral channels. The primary dispersion prism can be replaced with a Wollaston prism for dual polarization measurements. The spectrograph also has a pupil-viewing mode for alignment and calibration.

An Atmospheric Dispersion Corrector (ADC) uses a double-prism arrangement to nullify the vertical chromatic dispersion introduced by the atmosphere at non-zero zenith distances. The ADC installed in the Gemini Planet Imager (GPI) was first tested in August 2012 while the instrument was in the laboratory. GPI was installed at the Gemini South telescope in August 2013 and first light occurred later that year on November 11th. In this paper, we give an overview of the characterizations and performance of this ADC unit obtained in the laboratory and on sky, as well as the structure of its control software.

The Gemini Planet Imager (GPI) is a new facility, extreme adaptive optics (AO), coronagraphic instrument, currently being integrated onto the 8-meter Gemini South telescope, with the ultimate goal of directly imaging extrasolar planets. To achieve the contrast required for the desired science, it is necessary to quantify and mitigate wavefront error (WFE). A large source of potential static WFE arises from the primary AO wavefront sensor (WFS) detector's use of multiple readout segments with independent signal chains including on-chip preamplifiers and external amplifiers. Temperature changes within GPI's electronics cause drifts in readout segments' bias levels, inducing an RMS WFE of 1.1 nm and 41.9 nm over 4.44 degrees Celsius, for magnitude 4 and 11 stars, respectively. With a goal of $<$2 nm of static WFE, these are significant enough to require remedial action. Simulations imply a requirement to take fresh WFS darks every 2 degrees Celsius of temperature change, for a magnitude 6 star; similarly, for a magnitude 7 star, every 1 degree Celsius of temperature change. For sufficiently dim stars, bias drifts exceed the signal, causing a large initial WFE, and the former periodic requirement practically becomes an instantaneous/continuous one, making the goal of $<$2 nm of static WFE very difficult for stars of magnitude 9 or fainter. In extreme cases, this can cause the AO loops to destabilize due to perceived nonphysical wavefronts, as some of the WFS's Shack-Hartmann quadcells are split between multiple readout segments. Presented here is GPI's AO WFS geometry, along with detailed steps in the simulation used to quantify bias drift related WFE, followed by laboratory and on sky results, and concluded with possible methods of remediation.

The Gemini Planet Imager (GPI) Extreme Adaptive Optics Coronograph contains an interferometric mode: a 10-hole non-redundant mask (NRM) in its pupil wheel. GPI operates at $Y, J, H$, and $K$ bands, using an integral field unit spectrograph (IFS) to obtain spectral data at every image pixel. NRM on GPI is capable of imaging with a half resolution element inner working angle at moderate contrast, probing the region behind the coronagraphic spot. The fine features of the NRM PSF can provide a reliable check on the plate scale, while also acting as an attenuator for spectral standard calibrators that would otherwise saturate the full pupil. NRM commissioning data provides details about wavefront error in the optics as well as operations of adaptive optics control without pointing control from the calibration system. We compare lab and on-sky results to evaluate systematic instrument properties and examine the stability data in consecutive exposures. We discuss early on-sky performance, comparing images from integration and tests with the first on-sky images, and demonstrate resolving a known binary. We discuss the status of NRM and implications for future science with this mode.

The Gemini Planet Imager (GPI) is a new facility instrument for the Gemini Observatory designed to provide direct detection and characterization of planets and debris disks around stars in the solar neighborhood. In addition to its extreme adaptive optics and corona graphic systems which give access to high angular resolution and high-contrast imaging capabilities, GPI contains an integral field spectrograph providing low resolution spectroscopy across five bands between 0.95 and 2.5 $\mu$m. This paper describes the sequence of processing steps required for the spectro-photometric calibration of GPI science data, and the necessary calibration files. Based on calibration observations of the white dwarf HD 8049B we estimate that the systematic error in spectra extracted from GPI observations is less than 5%. The flux ratio of the occulted star and fiducial satellite spots within coronagraphic GPI observations, required to estimate the magnitude difference between a target and any resolved companions, was measured in the $H$-band to be $\Delta m = 9.23\pm0.06$ in laboratory measurements and $\Delta m = 9.39\pm 0.11$ using on-sky observations. Laboratory measurements for the $Y$, $J$, $K1$ and $K2$ filters are also presented. The total throughput of GPI, Gemini South and the atmosphere of the Earth was also measured in each photometric passband, with a typical throughput in $H$-band of 18% in the non-coronagraphic mode, with some variation observed over the six-month period for which observations were available. We also report ongoing development and improvement of the data cube extraction algorithm.

We present the results of both laboratory and on sky astrometric characterization of the Gemini Planet Imager (GPI). This characterization includes measurement of the pixel scale of the integral field spectrograph (IFS), the position of the detector with respect to north, and optical distortion. Two of these three quantities (pixel scale and distortion) were measured in the laboratory using two transparent grids of spots, one with a square pattern and the other with a random pattern. The pixel scale in the laboratory was also estimate using small movements of the artificial star unit (ASU) in the GPI adaptive optics system. On sky, the pixel scale and the north angle are determined using a number of known binary or multiple systems and Solar System objects, a subsample of which had concurrent measurements at Keck Observatory. Our current estimate of the GPI pixel scale is 14.14 $\pm$ 0.01 millarcseconds/pixel, and the north angle is -1.00 $\pm$ 0.03$\deg$. Distortion is shown to be small, with an average positional residual of 0.26 pixels over the field of view, and is corrected using a 5th order polynomial. We also present results from Monte Carlo simulations of the GPI Exoplanet Survey (GPIES) assuming GPI achieves ~1 milliarcsecond relative astrometric precision. We find that with this precision, we will be able to constrain the eccentricities of all detected planets, and possibly determine the underlying eccentricity distribution of widely separated Jovians.

The Gemini Planet Imager (GPI) has as its science instrument an infrared integral field spectrograph/polarimeter (IFS). Integral field spectrographs are scientificially powerful but require sophisticated data reduction systems. For GPI to achieve its scientific goals of exoplanet and disk characterization, IFS data must be reconstructed into high quality astrometrically and photometrically accurate datacubes in both spectral and polarization modes, via flexible software that is usable by the broad Gemini community. The data reduction pipeline developed by the GPI instrument team to meet these needs is now publicly available following GPI's commissioning. This paper, the first of a series, provides a broad overview of GPI data reduction, summarizes key steps, and presents the overall software framework and implementation. Subsequent papers describe in more detail the algorithms necessary for calibrating GPI data. The GPI data reduction pipeline is open source, available from planetimager.org, and will continue to be enhanced throughout the life of the instrument. It implements an extensive suite of task primitives that can be assembled into reduction recipes to produce calibrated datasets ready for scientific analysis. Angular, spectral, and polarimetric differential imaging are supported. Graphical tools automate the production and editing of recipes, an integrated calibration database manages reference files, and an interactive data viewer customized for high contrast imaging allows for exploration and manipulation of data.

The Gemini Planet Imager instrument's adaptive optics (AO) subsystem was designed specifically to facilitate high-contrast imaging. It features several new technologies, including computationally efficient wavefront reconstruction with the Fourier transform, modal gain optimization every 8 seconds, and the spatially filtered wavefront sensor. It also uses a Linear-Quadratic-Gaussian (LQG) controller (aka Kalman filter) for both pointing and focus. We present on-sky performance results from verification and commissioning runs from December 2013 through May 2014. The efficient reconstruction and modal gain optimization are working as designed. The LQG controllers effectively notch out vibrations. The spatial filter can remove aliases, but we typically use it oversized by about 60% due to stability problems.

Searching for nearby habitable worlds with direct imaging and spectroscopy will require a telescope large enough to provide angular resolution and sensitivity to planets around a significant sample of stars. Segmented telescopes are a compelling option to obtain such large apertures. However, these telescope designs have a complex geometry (central obstruction, support structures, segmentation) that makes high-contrast imaging more challenging. We are developing a new high-contrast imaging testbed at STScI to provide an integrated solution for wavefront control and starlight suppression on complex aperture geometries. We present our approach for the testbed optical design, which defines the surface requirements for each mirror to minimize the amplitude-induced errors from the propagation of out-of-pupil surfaces. Our approach guarantees that the testbed will not be limited by these Fresnel propagation effects, but only by the aperture geometry. This approach involves iterations between classical ray-tracing optical design optimization, and end-to-end Fresnel propagation with wavefront control (e.g. Electric Field Conjugation / Stroke Minimization). The construction of the testbed is planned to start in late Fall 2013.

We have built a 45-cm long x-ray deformable mirror of super-polished single-crystal silicon that has 45 actuators along the tangential axis. After assembly the surface height error was 19 nm rms. With use of high-precision visible-light metrology and precise control algorithms, we have actuated the x-ray deformable mirror and flattened its entire surface to 0.7 nm rms controllable figure error. This is, to our knowledge, the first sub-nanometer active flattening of a substrate longer than 15 cm.

The Gemini Planet Imager (GPI) is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of GPI has been tuned for maximum sensitivity to faint planets near bright stars. During first light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-sigma contrast of $10^6$ at 0.75 arcseconds and $10^5$ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-second exposure with minimal post-processing. Beta Pictoris b is observed at a separation of $434 \pm 6$ milli-arcseconds and position angle $211.8 \pm 0.5$ deg. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of three improvement in most parameters over previous solutions. The planet orbits at a semi-major axis of $9.0^{+0.8}_{-0.4}$ AU near the 3:2 resonance with the previously-known 6 AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% posterior probability of a transit of the planet in late 2017.

With a combination of adaptive optics imaging and a multi-epoch common proper motion search, we have conducted a large volume-limited (D $\le$ 75 pc) multiplicity survey of A-type stars, sensitive to companions beyond 30 au. The sample for the Volume-limited A-STar (VAST) survey consists of 435 A-type stars: 363 stars were observed with adaptive optics, 228 stars were searched for wide common proper motion companions and 156 stars were measured with both techniques. The projected separation coverage of the VAST survey extends from 30 to 45,000 au. A total of 137 stellar companions were resolved, including 64 new detections from the VAST survey, and the companion star fraction, projected separation distribution and mass ratio distribution were measured. The separation distribution forms a log-normal distribution similar to the solar-type binary distribution, but with a peak shifted to a significantly wider value of 387 (+132,-98) au. Integrating the fit to the distribution over the 30 to 10,000 au observed range, the companion star fraction for A-type stars is estimated as 33.8%+-2.6%. The mass ratio distribution of closer (<125 au) binaries is distinct from that of wider systems, with a flat distribution for close systems and a distribution that tends towards smaller mass ratios for wider binaries. Combining this result with previous spectroscopic surveys of A-type stars gives an estimate of the total companion star fraction of 68.9%+-7.0%. The most complete assessment of higher order multiples was estimated from the 156-star subset of the VAST sample with both adaptive optics and common proper motion measurements, combined with a literature search for companions, yielding a lower limit on the frequency of single, binary, triple, quadruple and quintuple A-type star systems of 56.4 (-4.0,+3.8), 32.1 (-3.5,+3.9), 9.0 (-1.8,+2.8), 1.9 (-0.6,+1.8) and 0.6 (-0.2,+1.4) per cent, respectively.

The Gemini Planet Imager (GPI) is currently in its final Acceptance & Testing stages. GPI is an XAO system based on a tweeter & woofer architecture (43 & 9 actuators respectively across the pupil), with the tweeter being a Boston Michromachines $64^2$ MEMS device. The XAO AO system is tightly integrated with a Lyot apodizing coronagraph. Acceptance testing started in February 2013 at the University of California, Santa Cruz. A conclusive acceptance review was held in July 2013 and the instrument was found ready for shipment to the Gemini South telescope on Cerro Pachon, Chile. Commissioning at the telescope will take place by the end of 2013, matching the summer window of the southern hemisphere. According to current estimates the 3 year planet finding campaign (890 allocated hours) might discover, image, and spectroscopically analyze 20 to 40 new exo-planets. Final acceptance testing of the integrated instrument can always bring up surprises when using cold chamber and flexure rig installations. The latest developments are reported. Also, we will give an overview of GPI's lab performance, the interplay between subsystems such as the calibration unit (CAL) with the AO bench. We report on-going optimizations on the AO controller loop to filter vibrations and last but not least achieved contrast performance applying speckle nulling. Furthermore, we will give an outlook of possible but challenging future upgrades as the implementation of a predictive controller or exchanging the conventional 48x48 SH WFS with a pyramid. With the ELT era arising, GPI will proof as a versatile and path-finding testbed for AO technologies on the next generation of ground-based telescopes.

We present Herschel far-infrared and submillimeter maps of the debris disk associated with the HR 8799 planetary system. We resolve the outer disk emission at 70, 100, 160 and 250 um and detect the disk at 350 and 500 um. A smooth model explains the observed disk emission well. We observe no obvious clumps or asymmetries associated with the trapping of planetesimals that is a potential consequence of planetary migration in the system. We estimate that the disk eccentricity must be <0.1. As in previous work by Su et al. (2009), we find a disk with three components: a warm inner component and two outer components, a planetesimal belt extending from 100 - 310 AU, with some flexibility (+/- 10 AU) on the inner edge, and the external halo which extends to ~2000 AU. We measure the disk inclination to be 26 +/- 3 deg from face-on at a position angle of 64 deg E of N, establishing that the disk is coplanar with the star and planets. The SED of the disk is well fit by blackbody grains whose semi-major axes lie within the planetesimal belt, suggesting an absence of small grains. The wavelength at which the spectrum steepens from blackbody, 47 +/- 30 um, however, is short compared to other A star debris disks, suggesting that there are atypically small grains likely populating the halo. The PACS longer wavelength data yield a lower disk color temperature than do MIPS data (24 and 70 um), implying two distinct halo dust grain populations.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Project Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This document presents the final report of the SDT.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. The report of the WFIRST-AFTA Science Definition Team (SDT) presents a Design Reference Mission for WFIRST that employs one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the smaller aperture designs previously considered for WFIRST, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This short article, produced as a companion to the SDT report, summarizes the key points of the WFIRST-2.4 DRM. It highlights the remarkable opportunity that the 2.4-m telescope affords for advances in many fields of astrophysics and cosmology, including dark energy, the demographics and characterization of exoplanets, the evolution of galaxies and quasars, and the stellar populations of the Milky Way and its neighbors.

Modern coronagraphic systems require very precise alignment between optical components and can benefit greatly from automated image processing. We discuss three techniques commonly employed in the fields of computer vision and image analysis as applied to the Gemini Planet Imager, a new facility instrument for the Gemini South Observatory. We describe how feature extraction and clustering methods can be used to aid in automated system alignment tasks, and also present a search algorithm for finding regular features in science images used for calibration and data processing. Along with discussions of each technique, we present our specific implementation and show results of each one in operation.

The nature and even the existence of a putative planet-mass companion ("Fomalhaut b") to Fomalhaut has been debated since 2008. In the present paper we reanalyze the multi-epoch ACS/STIS/WFC3 Hubble Space Telescope (HST) optical/near infrared images on which the discovery and some other claims were based. We confirm that the HST images do reveal an object in orbit around Fomalhaut but the detailed results from our analysis differ in some ways from previous discussions. In particular, we do not confirm flux variability over a two-year interval at 0.6 microns wavelength and we detect Fomalhaut b for the first time at the short wavelength of 0.43 microns. We find that the HST image of Fomalhaut b at $0.8\,\mu$m may be extended beyond the PSF. We cannot determine from our astrometry if Fomalhaut b will cross or not the dust ring. The optical through mid-infrared spectral energy distribution (SED) of Fomalhaut b cannot be explained as due to direct or scattered radiation from a massive planet. We consider two models to explain the SED: (1) a large circumplanetary disk around an unseen planet and (2) the aftermath of a collision during the past 50-150 years of two Kuiper Belt-like objects of radii 50 km.

Determining the atmospheric structure and chemical composition of an exoplanet remains a formidable goal. Fortunately, advancements in the study of exoplanets and their atmospheres have come in the form of direct imaging - spatially resolving the planet from its parent star - which enables high-resolution spectroscopy of self-luminous planets in Jovian-like orbits. Here, we present a spectrum with numerous, well-resolved, molecular lines from both water and carbon monoxide from a massive planet orbiting less than 40 AU from the star HR 8799. These data reveal the planet's chemical composition, atmospheric structure, and surface gravity, confirming that it is indeed a young planet. The spectral lines suggest an atmospheric carbon-to-oxygn ratio greater than the host star's, providing hints about the planet's formation.

We discuss scientific, technical and programmatic issues related to the use of an NRO 2.4m telescope for the WFIRST initiative of the 2010 Decadal Survey. We show that this implementation of WFIRST, which we call "NEW WFIRST," would achieve the goals of the NWNH Decadal Survey for the WFIRST core programs of Dark Energy and Microlensing Planet Finding, with the crucial benefit of deeper and/or wider near-IR surveys for GO science and a potentially Hubble-like Guest Observer program. NEW WFIRST could also include a coronagraphic imager for direct detection of dust disks and planets around neighboring stars, a high-priority science and technology precursor for future ambitious programs to image Earth-like planets around neighboring stars.

To image faint substellar companions obscured by the stellar halo and speckles, scattered light from the bright primary star must be removed in hardware or software. We apply the "locally-optimized combination of images" (LOCI) algorithm to 1-minute Keck Observatory snapshots of GKM dwarfs in the Hyades using source diversity to determine the most likely PSF. We obtain a mean contrast of 10^{-2} at 0.01", 10^{-4} at <1", and 10^{-5} at 5". New brown dwarf and low-mass stellar companions to Hyades primaries are found in a third of the 84 targeted systems. This campaign shows the efficacy of LOCI on snapshot imaging as well as on bright wide binaries with off-axis LOCI, reaching contrasts sufficient for imaging 625-Myr late-L/early-T dwarfs purely in post-processing.

Small-angle coronagraphy is technically and scientifically appealing because it enables the use of smaller telescopes, allows covering wider wavelength ranges, and potentially increases the yield and completeness of circumstellar environment - exoplanets and disks - detection and characterization campaigns. However, opening up this new parameter space is challenging. Here we will review the four posts of high contrast imaging and their intricate interactions at very small angles (within the first 4 resolution elements from the star). The four posts are: choice of coronagraph, optimized wavefront control, observing strategy, and post-processing methods. After detailing each of the four foundations, we will present the lessons learned from the 10+ years of operations of zeroth and first-generation adaptive optics systems. We will then tentatively show how informative the current integration of second-generation adaptive optics system is, and which lessons can already be drawn from this fresh experience. Then, we will review the current state of the art, by presenting world record contrasts obtained in the framework of technological demonstrations for space-based exoplanet imaging and characterization mission concepts. Finally, we will conclude by emphasizing the importance of the cross-breeding between techniques developed for both ground-based and space-based projects, which is relevant for future high contrast imaging instruments and facilities in space or on the ground.

Breakthrough direct detections of planetary companions orbiting A-type stars confirm the existence of massive planets at relatively large separations, but dedicated surveys are required to estimate the frequency of similar planetary systems. To measure the first estimation of the giant exoplanetary systems frequency at large orbital separation around A-stars, we have conducted a deep-imaging survey of young (8-400 Myr), nearby (19-84 pc) A- and F-stars to search for substellar companions in the 10-300 AU range. The sample of 42 stars combines all A-stars observed in previous AO planet search surveys reported in the literature with new AO observations from VLT/NaCo and Gemini/NIRI. It represents an initial subset of the International Deep Planet Survey (IDPS) sample of stars covering M- to B-stars. The data were obtained with diffraction-limited observations in H- and Ks-band combined with angular differential imaging to suppress the speckle noise of the central stars, resulting in typical 5-sigma detection limits in magnitude difference of 12 mag at 1", 14 mag at 2" and 16 mag at 5" which is sufficient to detect massive planets. A detailed statistical analysis of the survey results is performed using Monte Carlo simulations. Considering the planet detections, we estimate the fraction of A-stars having at least one massive planet (3-14 MJup) in the range 5-320 AU to be inside 5.9-18.8% at 68% confidence, assuming a flat distribution for the mass of the planets. By comparison, the brown dwarf (15-75 MJup) frequency for the sample is 2.0-8.9% at 68% confidence in the range 5-320 AU. Assuming power law distributions for the mass and semimajor axis of the planet population, the AO data are consistent with a declining number of massive planets with increasing orbital radius which is distinct from the rising slope inferred from radial velocity (RV) surveys around evolved A-stars.

We present rotational velocities for individual components of eleven very low mass (VLM) binaries with spectral types between M7 and L7.5. These results are based on observations taken with the near-infrared spectrograph, NIRSPEC, and the Keck II laser guide star adaptive optics (LGS AO) system. We find that the observed sources tend to be rapid rotators (vsini > 10 km/s), consistent with previous seeing-limited measurements of VLM objects. The two sources with the largest vsini, LP 349-25B and HD 130948C, are rotating at ~30% of their break up speed, and are among the most rapidly rotating VLM objects known. Furthermore, five binary systems, all with orbital semi-major axes <3.5 AU, have component vsini values that differ by greater than 3sigma. To bring the binary components with discrepant rotational velocities into agreement would require the rotational axes to be inclined with respect to each other, and that at least one component is inclined with respect to the orbital plane. Alternatively, each component could be rotating at a different rate, even though they have similar spectral types. Both differing rotational velocities and inclinations have implications for binary star formation and evolution. We also investigate possible dynamical evolution in the triple system HD 130948A-BC. The close binary brown dwarfs B and C have significantly different vsini values. We demonstrate that components B and C could have been torqued into misalignment by the primary star, A, via orbital precession. Such a scenario can also be applied to another triple system in our sample, GJ 569A-Bab. Interactions such as these may play an important role in the dynamical evolution of very low mass binaries. Finally, we note that two of the binaries with large differences in component vsini, LP 349-25AB and 2MASS 0746+20AB, are also known radio sources.

Micro-electro-mechanical systems (MEMS) technology can provide for deformable mirrors (DMs) with excellent performance within a favorable economy of scale. Large MEMS-based astronomical adaptive optics (AO) systems such as the Gemini Planet Imager are coming on-line soon. As MEMS DM end-users, we discuss our decade of practice with the micromirrors, from inspecting and characterizing devices to evaluating their performance in the lab. We also show MEMS wavefront correction on-sky with the "Villages" AO system on a 1-m telescope, including open-loop control and visible-light imaging. Our work demonstrates the maturity of MEMS technology for astronomical adaptive optics.

As a part of our ongoing Volume-limited A-Star (VAST) adaptive optics survey, we have obtained observations of 26 binary systems with projected separations <100 AU, 13 of which have sufficient historical measurements to allow for refinement of their orbital elements. For each system with an estimated orbit, the dynamical system mass obtained was compared with the system mass estimated from mass-magnitude relations. Discrepancies between the dynamical and theoretical system mass can be explained by the presence of a previously unresolved spectroscopic component, or by a non-solar metallicity of the system. Using this approach to infer the presence of additional companions, a lower limit to the fraction of binaries, triples, and quadruples can be estimated as 39, 46, and 15 per cent, for systems with at least one companion within 100 AU. The fraction of multiple systems with three or more components shows a relative increase compared to the fraction for Solar-type primaries resolved in previous volume-limited surveys. The observations have also revealed a pair of potentially young ($<$100 Myr) M-dwarf companions, which would make an ideal benchmark for the theoretical models during the pre-Main Sequence contraction phase for M-dwarfs. In addition to those systems with orbit fits, we report 13 systems for which further orbital monitoring observations are required, 11 of which are newly resolved as a part of the VAST survey.

With adaptive optics imaging at Keck observatory, we have discovered a substellar companion to the F6 Pleiades star HD 23514, one of the dustiest main-sequence stars known to date (L_{IR}/L_{*}~2%). This is one of the first brown dwarfs discovered as a companion to a star in the Pleiades. The 0.06 M$_\odot$ late-M secondary has a projected separation of ~360 AU. The scarcity of substellar companions to stellar primaries in the Pleiades combined with the extremely dusty environment make this a unique system to study.

Multi-wavelength observations/spectroscopy of exoplanetary atmospheres are the basis of the emerging exciting field of comparative exoplanetology. The HR 8799 planetary system is an ideal laboratory to study our current knowledge gap between massive field brown dwarfs and the cold 5-Gyr old Solar system planets. The HR 8799 planets have so far been imaged at J- to L-band, with only upper limits available at M-band. We present here deep high-contrast Keck II adaptive optics M-band observations that show the imaging detection of 3 of the 4 currently known HR 8799 planets. Such detections were made possible due to the development of an innovative LOCI-based background subtraction scheme that is 3 times more efficient than a classical median background subtraction for Keck II AO data, representing a gain in telescope time of up to a factor of 9. These M-band detections extend the broad band photometric coverage out to about 5 microns and provide access to the strong CO fundamental absorption band at 4.5microns. The new M-band photometry shows that the HR 8799 planets are located near the L/T-type dwarf transition, similar to what was found by other studies. We also confirm that the best atmospheric fits are consistent with low surface gravity, dusty and non-equilibrium CO/CH4 chemistry models.

The properties of 2M1207b, a young (~ 8 Myr) planet-mass companion, have lacked a satisfactory explanation for some time. The combination of low luminosity, red near-IR colors, and L-type near-IR spectrum (previously consistent with Teff ~ 1600K) implies an abnormally small radius. Early explanations for the apparent underluminosity of 2M1207b invoked an edge-on disk or the remnant of a recent protoplanetary collision. The discovery of a second planet-mass object (HR8799b) with similar luminosity and colors as 2M1207b indicates that a third explanation, one of a purely atmospheric nature, is more likely. By including clouds, non-equilibrium chemistry, and low-gravity, an atmosphere with effective temperature consistent with evolution cooling-track predictions is revealed. Consequently, 2M1207b, and others like it, requires no new physics to explain nor do they belong to a new class of objects. Instead they most likely represent the natural extension of cloudy substellar atmospheres down to low Teff and log(g). If this atmosphere only explanation for 2M1207b is correct, then very young planet-mass objects with near-IR spectra similar to field T dwarfs may be rare.

The Gemini Planet Imager (GPI) is a high performance adaptive optics system being designed and built for the Gemini Observatory. GPI is optimized for high contrast imaging, combining precise and accurate wavefront control, diffraction suppression, and a speckle-suppressing science camera with integral field and polarimetry capabilities. The primary science goal for GPI is the direct detection and characterization of young, Jovian-mass exoplanets. For plausible assumptions about the distribution of gas giant properties at large semi-major axes, GPI will be capable of detecting more than 10% of gas giants more massive than 0.5 M_J around stars younger than 100 Myr and nearer than 75 parsecs. For systems younger than 1 Gyr, gas giants more massive than 8 M_J and with semi-major axes greater than 15 AU are detected with completeness greater than 50%. A survey targeting young stars in the solar neighborhood will help determine the formation mechanism of gas giant planets by studying them at ages where planet brightness depends upon formation mechanism. Such a survey will also be sensitive to planets at semi-major axes comparable to the gas giants in our own solar system. In the simple, and idealized, situation in which planets formed by either the "hot-start" model of Burrows et al. (2003) or the core accretion model of Marley et al. (2007), a few tens of detected planets are sufficient to distinguish how planets form.

With an adaptive optics imaging survey of 148 B6-A7 stars, we have tested the hypothesis that unresolved lower-mass companions are the source of the unexpected X-ray detections of stars in this spectral type range. The sample is composed of 63 stars detected in X-rays within the ROSAT All-Sky Survey and 85 stars that form a control sample; both subsets have the same restricted distribution of spectral type, age, X-ray sensitivity and separation coverage. A total of 68 companion candidates are resolved with separations ranging from 0.3" to 26.2", with 23 new detections. The multiple star frequency of the X-ray sample based on companions resolved within the ROSAT error ellipse is found to be 43 (+6,-6)%. The corresponding control sample multiple star frequency is three times lower at 12 (+4,-3)% -- a difference of 31\pm7%. These results are presented in the first of a series of papers based on our Volume-limited A-Star (VAST) survey -- a comprehensive study of the multiplicity of A-type stars.

Using the integral field spectrograph OSIRIS, on the Keck II telescope, broad near-infrared H and K-band spectra of the young exoplanet HR8799b have been obtained. In addition, six new narrow-band photometric measurements have been taken across the H and K bands. These data are combined with previously published photometry for an analysis of the planet's atmospheric properties. Thick photospheric dust cloud opacity is invoked to explain the planet's red near-IR colors and relatively smooth near-IR spectrum. Strong water absorption is detected, indicating a Hydrogen-rich atmosphere. Only weak CH4 absorption is detected at K band, indicating efficient vertical mixing and a disequilibrium CO/CH4 ratio at photospheric depths. The H-band spectrum has a distinct triangular shape consistent with low surface gravity. New giant planet atmosphere models are compared to these data with best fitting bulk parameters, Teff = 1100K +- 100 and log(g) = 3.5 +-0.5 (for solar composition). Given the observed luminosity, log(Lobs/Lsun) ~ -5.1, these values correspond to a radius of 0.75 Rjup (+0.17, -0.12) and mass ~ 0.72 Mjup (+2.6,-0.6) -- strikingly inconsistent with interior/evolution models. Enhanced metallicity (up to ~ 10 times that of the Sun) along with thick clouds and non-equilibrium chemistry are likely required to reproduce the complete ensemble of spectroscopic and photometric data and the low effective temperatures (< 1000K) required by the evolution models.

Direct imaging of the HR8799 system was a major achievement in the study of exoplanets. HR8799 is a $\gamma$\,Doradus variable and asteroseismology can provide an independent constraint on the inclination. Using 650 high signal-to-noise, high resolution, full visual wavelength spectroscopic observations obtained over two weeks at Observatoire de Haute Provence (OHP) with the SOPHIE spectrograph we find that the main frequency in the radial velocity data is 1.9875 d$^{-1}$. This frequency corresponds to the main frequency as found in previous photometric observations. Using the FAMIAS software to identify the pulsation modes, we find this frequency is a prograde $\ell$=1 sectoral mode and obtain the constraint that inclination $i\gtrsim$40$^{\circ}$.

High-contrast near-infrared imaging of the nearby star HR 8799 has shown three giant planets. Such images were possible due to the wide orbits (> 25 AU) and youth (< 100 Myr) of the imaged planets, which are still hot and bright as they radiate away gravitational energy acquired during their formation. A major area of contention in the extrasolar planet community is whether outer planets (> 10 AU) more massive than Jupiter form via one-step gravitational instabilities or, rather, via a two-step process involving accretion of a core followed by accumulation of a massive outer envelope composed primarily of hydrogen and helium. Here we report the presence of a fourth planet, interior to and about the same mass as the other three. The system, with this additional planet, represents a challenge for current planet formation models as none of them can explain the in situ formation of all four planets. With its four young giant planets and known cold/warm debris belts, the HR 8799 planetary system is a unique laboratory to study the formation and evolution of giant planets at wide > 10 AU separations.

At the University of California's Lick Observatory, we have implemented an on-sky testbed for next-generation adaptive optics (AO) technologies. The Visible-Light Laser Guidestar Experiments instrument (ViLLaGEs) includes visible-light AO, a micro-electro-mechanical-systems (MEMS) deformable mirror, and open-loop control of said MEMS on the 1-meter Nickel telescope at Mt. Hamilton. In this paper we evaluate the performance of ViLLaGEs in open- and closed-loop control, finding that both control methods give equivalent Strehl ratios of up to ~ 7% in I-band and similar rejection of temporal power. Therefore, we find that open-loop control of MEMS on-sky is as effective as closed-loop control. Furthermore, after operating the system for three years, we find MEMS technology to function well in the observatory environment. We construct an error budget for the system, accounting for 130 nm of wavefront error out of 190 nm error in the science-camera PSFs. We find that the dominant known term is internal static error, and that the known contributions to the error budget from open-loop control (MEMS model, position repeatability, hysteresis, and WFS linearity) are negligible.

The InfraRed Imaging Spectrograph (IRIS) is a first-light instrument being designed for the Thirty Meter Telescope (TMT). IRIS is a combination of an imager that will cover a 16.4" field of view at the diffraction limit of TMT (4 mas sampling), and an integral field unit spectrograph that will sample objects at 4-50 mas scales. IRIS will open up new areas of observational parameter space, allowing major progress in diverse fields of astronomy. We present the science case and resulting requirements for the performance of IRIS. Ultimately, the spectrograph will enable very well-resolved and sensitive studies of the kinematics and internal chemical abundances of high-redshift galaxies, shedding light on many scenarios for the evolution of galaxies at early times. With unprecedented imaging and spectroscopy of exoplanets, IRIS will allow detailed exploration of a range of planetary systems that are inaccessible with current technology. By revealing details about resolved stellar populations in nearby galaxies, it will directly probe the formation of systems like our own Milky Way. Because it will be possible to directly characterize the stellar initial mass function in many environments and in galaxies outside of the the Milky Way, IRIS will enable a greater understanding of whether stars form differently in diverse conditions. IRIS will reveal detailed kinematics in the centers of low-mass galaxies, allowing a test of black hole formation scenarios. Finally, it will revolutionize the characterization of reionization and the first galaxies to form in the universe.

We present a new analysis of the Jupiter+Saturn analog system, OGLE-2006-BLG-109Lb,c, which was the first double planet system discovered with the gravitational microlensing method. This is the only multi-planet system discovered by any method with measured masses for the star and both planets. In addition to the signatures of two planets, this event also exhibits a microlensing parallax signature and finite source effects that provide a direct measure of the masses of the star and planets, and the expected brightness of the host star is confirmed by Keck AO imaging, yielding masses of M_* = 0.51(+0.05-0.04) M_sun, M_b = 231+-19 M_earth, M_c = 86+-7 M_earth. The Saturn-analog planet in this system had a planetary light curve deviation that lasted for 11 days, and as a result, the effects of the orbital motion are visible in the microlensing light curve. We find that four of the six orbital parameters are tightly constrained and that a fifth parameter, the orbital acceleration, is weakly constrained. No orbital information is available for the Jupiter-analog planet, but its presence helps to constrain the orbital motion of the Saturn-analog planet. Assuming co-planar orbits, we find an orbital eccentricity of eccentricity = 0.15 (+0.17-0.10) and an orbital inclination of i = 64 (+4-7) deg. The 95% confidence level lower limit on the inclination of i > 49 deg. implies that this planetary system can be detected and studied via radial velocity measurements using a telescope of >30m aperture.

This is a scientific strategy for the detection and characterization of extrasolar planets; that is, planets orbiting other stars. As such, it maps out over a 15-year horizon the techniques and capabilities required to detect and measure the properties of planets as small as Earth around stars as large as our own Sun. It shows how the technology pieces and their development fit together to achieve the primary goal of the strategy: if planets like Earth exist around stars within some tens of light years of our own Solar System, those planets will be found and their basic properties characterized. Essential to this strategy is not only the search for and examination of individual planets, but also a knowledge of the arrangement, or architecture, of planetary systems around as large a number of stars as possible; this is the second goal of the strategy. The final goal of the strategy is the study of disks around stars, important both to understand the implications of the variety of exoplanet systems for planet formation, and to determine how many nearby stars have environments around them clean enough of debris that planets may be sought and, if found, characterized.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

The Gemini Planet Imager (GPI) is an instrument designed for the Gemini South telescope to image young Jupiter-mass planets in the infrared. To achieve the high contrast needed for this, it employs an apodized pupil Lyot coronagraph (APLC) to remove most of the starlight. Current designs use a partially-transmitting apodizer in the pupil; we examine the use of binary apodizations in the form of starshaped shaped pupils, and present a design that could achieve comparable performance, along with a series of design guidelines for creating shaped pupil versions of APLCs in other systems.

High-contrast imaging of extrasolar planet candidates around a main-sequence star has recently been realized from the ground using current adaptive optics (AO) systems. Advancing such observations will be a task for the Gemini Planet Imager, an upcoming "extreme" AO instrument. High-order "tweeter" and low-order "woofer" deformable mirrors (DMs) will supply a >90%-Strehl correction, a specialized coronagraph will suppress the stellar flux, and any planets can then be imaged in the "dark hole" region. Residual wavefront error scatters light into the DM-controlled dark hole, making planets difficult to image above the noise. It is crucial in this regard that the high-density tweeter, a micro-electrical mechanical systems (MEMS) DM, have sufficient stroke to deform to the shapes required by atmospheric turbulence. Laboratory experiments were conducted to determine the rate and circumstance of saturation, i.e. stroke insufficiency. A 1024-actuator 1.5-um-stroke MEMS device was empirically tested with software Kolmogorov-turbulence screens of r_0=10-15cm. The MEMS when solitary suffered saturation ~4% of the time. Simulating a woofer DM with ~5-10 actuators across a 5-m primary mitigated MEMS saturation occurrence to a fraction of a percent. While no adjacent actuators were saturated at opposing positions, mid-to-high-spatial-frequency stroke did saturate more frequently than expected, implying that correlations through the influence functions are important. Analytical models underpredict the stroke requirements, so empirical studies are important.

Direct imaging of exoplanetary systems is a powerful technique that can reveal Jupiter-like planets in wide orbits, can enable detailed characterization of planetary atmospheres, and is a key step towards imaging Earth-like planets. Imaging detections are challenging due to the combined effect of small angular separation and large luminosity contrast between a planet and its host star. High-contrast observations with the Keck and Gemini telescopes have revealed three planets orbiting the star HR 8799, with projected separations of 24, 38, and 68 astronomical units. Multi-epoch data show counter-clockwise orbital motion for all three imaged planets. The low luminosity of the companions and the estimated age of the system imply planetary masses between 5 and 13 times that of Jupiter. This system resembles a scaled-up version of the outer portion of our Solar System.

Searches for extrasolar planets have uncovered an astonishing diversity of planetary systems, yet the frequency of solar system analogs remains unknown. The gravitational microlensing planet search method is potentially sensitive to multiple-planet systems containing analogs of all the solar system planets except Mercury. We report the detection of a multiple-planet system with microlensing. We identify two planets with masses of ~0.71 and ~0.27 times the mass of Jupiter and orbital separations of ~2.3 and ~4.6 astronomical units orbiting a primary star of mass ~0.50 solar masses at a distance of ~1.5 kiloparsecs. This system resembles a scaled version of our solar system in that the mass ratio, separation ratio, and equilibrium temperatures of the planets are similar to those of Jupiter and Saturn. These planets could not have been detected with other techniques; their discovery from only six confirmed microlensing planet detections suggests that solar system analogs may be common.

We have carried out a multiwavelength observational campaign demonstrating some of the remarkable properties of the infrared-bright variable star BP Psc. Surrounded by a compact dusty, gaseous disk, this little-studied late-G (or early-K) type star emits about 75% of its detected energy flux at infrared wavelengths. Evidence for accretion of gas in conjunction with narrow bi-polar jets and Herbig-Haro objects is apparently consistent with classification of BP Psc as a pre-main sequence star, as postulated in most previous studies. If young, then BP Psc would be one of the nearest and oldest known classical T Tauri stars. However, such an evolutionary classification encounters various problems that are absent or much less severe if BP Psc is instead a luminosity class III post-main sequence star. In this case, it would be the first known example of a first ascent giant surrounded by a massive molecular disk with accompanying rapid gas accretion and prominent jets and HH objects. In this model, the genesis of the massive dusty gaseous disk could be a consequence of the envelopment of a low mass companion star. Properties in the disk may be conducive to the current formation of planets, a gigayear or more after the formation of BP Psc itself.

In long adaptive optics corrected exposures, exoplanet detections are currently limited by speckle noise originating from the telescope and instrument optics, and it is expected that such noise will also limit future high-contrast imaging instruments for both ground and space-based telescopes. Previous theoretical analysis have shown that the time intensity variations of a single speckle follows a modified Rician. It is first demonstrated here that for a circular pupil this temporal intensity distribution also represents the speckle spatial intensity distribution at a fix separation from the point spread function center; this fact is demonstrated using numerical simulations for coronagraphic and non-coronagraphic data. The real statistical distribution of the noise needs to be taken into account explicitly when selecting a detection threshold appropriate for some desired confidence level. In this paper, a technique is described to obtain the pixel intensity distribution of an image and its corresponding confidence level as a function of the detection threshold. Using numerical simulations, it is shown that in the presence of speckles noise, a detection threshold up to three times higher is required to obtain a confidence level equivalent to that at 5sigma for Gaussian noise. The technique is then tested using TRIDENT CFHT and angular differential imaging NIRI Gemini adaptive optics data. It is found that the angular differential imaging technique produces quasi-Gaussian residuals, a remarkable result compared to classical adaptive optic imaging. A power-law is finally derived to predict the 1-3*10^-7 confidence level detection threshold when averaging a partially correlated non-Gaussian noise.

We present the results of the Gemini Deep Planet Survey, a near-infrared adaptive optics search for giant planets and brown dwarfs around nearby young stars. The observations were obtained with the Altair adaptive optics system at the Gemini North telescope and angular differential imaging was used to suppress the speckle noise of the central star. Detection limits for the 85 stars observed are presented, along with a list of all faint point sources detected around them. Typically, the observations are sensitive to angular separations beyond 0.5" with 5-sigma contrast sensitivities in magnitude difference at 1.6 micron of 9.5 at 0.5", 12.9 at 1", 15.0 at 2", and 16.5 at 5". For the typical target of the survey, a 100 Myr old K0 star located 22 pc from the Sun, the observations are sensitive enough to detect planets more massive than 2 Mjup with a projected separation in the range 40-200 AU. Second epoch observations of 48 stars with candidates (out of 54) have confirmed that all candidates are unrelated background stars. A detailed statistical analysis of the survey results, yielding upper limits on the fractions of stars with giant planet or low mass brown dwarf companions, is presented. Assuming a planet mass distribution dN/dm m^{-1.2} and a semi-major axis distribution dN/da a^{-1}, the 95% credible upper limits on the fraction of stars with at least one planet of mass 0.5-13 Mjup are 0.28 for the range 10-25 AU, 0.13 for 25-50 AU, and 0.093 for 50-250 AU. The 95% credible interval for the fraction of stars with at least one brown dwarf companion having a semi-major axis in the range 25-250 AU is 0.019 (-0.015/+0.083), irrespective of any assumption on the mass and semi-major axis distributions. The stars HD 14802, HD 166181, and HD 213845 have been resolved into binaries for the first time.

The Gemini Planet (GPI) imager is an "extreme" adaptive optics system being designed and built for the Gemini Observatory. GPI combines precise and accurate wavefront control, diffraction suppression, and a speckle-suppressing science camera with integral field and polarimetry capabilities. GPI's primary science goal is the direct detection and characterization of young, Jovian-mass exoplanets. For systems younger than 2 Gyr exoplanets more massive than 6 MJ and semimajor axes beyond 10 AU are detected with completeness greater than 50%. GPI will also discover faint debris disks, explore icy moons and minor planets in the solar system, reveal high dynamic range main-sequence binaries, and study mass loss from evolved stars. This white paper explains the role of GPI in exoplanet discovery and characterization and summarizes our recommendations to the NSF-NASA-DOE Astronomy and Astrophysics Advisory Committee ExoPlanet Task Force.

We imaged five objects near the star forming clouds of Ophiuchus with the Keck Laser Guide Star AO system. We resolved Allers et al. (2006)'s #11 (Oph 16222-2405) and #16 (Oph 16233-2402) into binary systems. The #11 object is resolved into a 243 AU binary, the widest known for a very low mass (VLM) binary. The binary nature of #11 was discovered first by Allers (2005) and independently here during which we obtained the first spatially resolved R~2000 near-infrared (J & K) spectra, mid-IR photometry, and orbital motion estimates. We estimate for 11A and 11B gravities (log(g)>3.75), ages (5+/-2 Myr), luminosities (log(L/Lsun)=-2.77+/-0.10 and -2.96+/-0.10), and temperatures (Teff=2375+/-175 and 2175+/-175 K). We find self-consistent DUSTY evolutionary model (Chabrier et al. 2000) masses of 17+4-5 MJup and 14+6-5 MJup, for 11A and 11B respectively. Our masses are higher than those previously reported (13-15 MJup and 7-8 MJup) by Jayawardhana & Ivanov (2006b). Hence, we find the system is unlikely a ``planetary mass binary'', (in agreement with Luhman et al. 2007) but it has the second lowest mass and lowest binding energy of any known binary. Oph #11 and Oph #16 belong to a newly recognized population of wide (>100 AU), young (<10 Myr), roughly equal mass, VLM stellar and brown dwarf binaries. We deduce that ~6+/-3% of young (<10 Myr) VLM objects are in such wide systems. However, only 0.3+/-0.1% of old field VLM objects are found in such wide systems. Thus, young, wide, VLM binary populations may be evaporating, due to stellar encounters in their natal clusters, leading to a field population depleted in wide VLM systems.

We present the discovery of 2MASS J21321145+1341584AB as a closely separated (0.066") very low-mass field dwarf binary resolved in the near-infrared by the Keck II Telescope using laser guide star adaptive optics. Physical association is deduced from the angular proximity of the components and constraints on their common proper motion. We have obtained a near-infrared spectrum of the binary and find that it is best described by an L5+/-0.5 primary and an L7.5+/-0.5 secondary. Model-dependent masses predict that the two components straddle the hydrogen burning limit threshold with the primary likely stellar and the secondary likely substellar. The properties of this sytem - close projected separation (1.8+/-0.3 AU) and near unity mass ratio - are consistent with previous results for very low-mass field binaries. The relatively short estimated orbital period of this system (~7-12 yr) makes it a good target for dynamical mass measurements. Interestingly, the system's angular separation is the tightest yet for any very low-mass binary published from a ground-based telescope and is the tightest binary discovered with laser guide star adaptive optics to date.

We present the results of a five year monitoring campaign of the close binary TWA 5Aab in the TW Hydrae association, using speckle and adaptive optics on the W.M. Keck 10 m telescopes. These measurements were taken as part of our ongoing monitoring of pre-main sequence (PMS) binaries in an effort to increase the number of dynamically determined PMS masses and thereby calibrate the theoretical PMS evolutionary tracks. Our observations have allowed us to obtain the first determination of this system's astrometric orbit. We find an orbital period of 5.94 +- 0.09 years and a semi-major axis of 0.066" +- 0.005". Combining these results with a kinematic distance, we calculate a total mass of 0.71 +- 0.14 M_sun (D/44 pc)^3. for this system. This mass measurement, as well as the estimated age of this system, are consistent to within 2$\sigma$ of all theoretical models considered. In this analysis, we properly account for correlated uncertainties, and show that while these correlations are generally ignored, they increase the formal uncertainties by up to a factor of five and therefore are important to incorporate. With only a few more years of observation, this type of measurement will allow the theoretical models to be distinguished.

We have re-analyzed archival HST and Subaru data of the recently discovered planetary mass companion (PMC) GQ Lup B. With these we produce the first R and I band photometry of the companion and fit a radius and effective temperature using detailed model atmospheres. We find an effective temperature of 2335 +/- 100K, and a radius of 0.38 +/- 0.05 Rsol and luminosity of log(L/Lsol) = -2.42 +/- 0.07 (at 140pc). Since we fit wavelengths that span most of the emitted radiation from GQ Lup this luminosity estimate is robust, with uncertainty dominated by the distance uncertainty (+/- 50 pc). The radius obtained for 140pc (0.38Rsol) is significantly larger than the one originally derived and larger than model predictions. The mass of the object is much more model-dependent than the radiative properties, but for the GAIA dusty models we find a mass between 10-20 MJup, in the range of the brown dwarf and PMC deuterium burning boundary. Assuming a distance of 140pc, observations fit to 1sigma the Baraffe evolution model for a ~15 MJup brown dwarf. Additionally, the F606W photometric band is significantly overluminous compared to model predictions and other brown dwarfs. Such overluminosity could be explained by a bright Halpha emission from chromospheric activity, interaction with another undetected companion, or accretion. Assuming that GQ Lup B has a bright Halpha emission line, its Halpha emission strength is 10^(-1.71 +/- 0.10) Lbol, significantly larger than field late-type dwarfs. GQ Lup B might be strongly accreting and still be in its formation phase.

Multi-band (0.9 to 1.6 um) images of the TW Hydrae Association (TWA) brown dwarf, 2MASS J1207334-393254 (also known as 2M1207), and its candidate planetary mass companion (2M1207b) were obtained on 2004 Aug 28 and 2005 Apr 26 with HST/NICMOS. The images from these two epochs unequivocally confirm the two objects as a common proper motion pair (16.0 sigma confidence). A new measurement of the proper motion of 2M1207 implies a distance to the system of 59+-7 pc and a projected separation of 46+-5 AU. The NICMOS and previously published VLT photometry of 2M1207b, extending overall from 0.9 to 3.8 um, are fully consistent with an object of a few Jupiter masses at the canonical age of a TWA member (~8 Myr) based on evolutionary models of young giant planets. These observations provide information on the physical nature of 2M1207b and unambiguously establish that the first direct image of a planetary mass companion in orbit around a self-luminous body, other than our Sun, has been secured.

Imaging faint companions (exoplanets and brown dwarfs) around nearby stars is currently limited by speckle noise. To efficiently attenuate this noise, a technique called simultaneous spectral differential imaging (SSDI) can be used. This technique consists of acquiring simultaneously images of the field of view in several adjacent narrow bands and in combining these images to suppress speckles. Simulations predict that SSDI can achieve, with the acquisition of three wavelengths, speckle noise attenuation of several thousands. These simulations are usually performed using the Fraunhofer approximation, i.e. considering that all aberrations are located in the pupil plane. We have performed wavefront propagation simulations to evaluate how out-of-pupil-plane aberrations affect SSDI speckle noise attenuation performance. The Talbot formalism is used to give a physical insight of the problem; results are confirmed using a proper wavefront propagation algorithm. We will show that near-focal-plane aberrations can significantly reduce SSDI speckle noise attenuation performance at several lambda/D separation. It is also shown that the Talbot effect correctly predicts the PSF chromaticity. Both differential atmospheric refraction effects and the use of a coronagraph will be discussed.

Accurate astrometry and photometry of saturated and coronagraphic point spread functions (PSFs) are fundamental to both ground- and space-based high contrast imaging projects. For ground-based adaptive optics imaging, differential atmospheric refraction and flexure introduce a small drift of the PSF with time, and seeing and sky transmission variations modify the PSF flux distribution. For space-based imaging, vibrations, thermal fluctuations and pointing jitters can modify the PSF core position and flux. These effects need to be corrected to properly combine the images and obtain optimal signal-to-noise ratios, accurate relative astrometry and photometry of detected objects as well as precise detection limits. Usually, one can easily correct for these effects by using the PSF core, but this is impossible when high dynamic range observing techniques are used, like coronagrahy with a non-transmissive occulting mask, or if the stellar PSF core is saturated. We present a new technique that can solve these issues by using off-axis satellite PSFs produced by a periodic amplitude or phase mask conjugated to a pupil plane. It will be shown that these satellite PSFs track precisely the PSF position, its Strehl ratio and its intensity and can thus be used to register and to flux normalize the PSF. A laboratory experiment is also presented to validate the theory. This approach can be easily implemented in existing adaptive optics instruments and should be considered for future extreme adaptive optics coronagraph instruments and in high-contrast imaging space observatories.

Angular differential imaging is a high-contrast imaging technique that reduces quasi-static speckle noise and facilitates the detection of nearby companions. A sequence of images is acquired with an altitude/azimuth telescope while the instrument field derotator is switched off. This keeps the instrument and telescope optics aligned and allows the field of view to rotate with respect to the instrument. For each image, a reference PSF is constructed from other appropriately-selected images of the same sequence and subtracted to remove quasi-static PSF structure. All residual images are then rotated to align the field and are combined. Observed performances are reported for Gemini North data. It is shown that quasi-static PSF noise can be reduced by a factor \~5 for each image subtraction. The combination of all residuals then provides an additional gain of the order of the square root of the total number of acquired images. A total speckle noise attenuation of 20-50 is obtained for one-hour long observing sequences compared to a single 30s exposure. A PSF noise attenuation of 100 was achieved for two-hour long sequences of images of Vega, reaching a 5-sigma contrast of 20 magnitudes for separations greater than 8". For a 30-minute long sequence, ADI achieves 30 times better signal-to-noise than a classical observation technique. The ADI technique can be used with currently available instruments to search for ~1MJup exoplanets with orbits of radii between 50 and 300 AU around nearby young stars. The possibility of combining the technique with other high-contrast imaging methods is briefly discussed.

It is now well established that a large fraction of the low-mass stars are binaries or higher order multiples. Similarly a sizable fraction have giant planets. In contrast to these, the situation for brown dwarf companions is complicated: While close systems seem to be extremely rare, wide systems are possibly more common. In this paper, we present new results on a survey for low-mass companions in the Hyades. After measuring precisely the radial velocity of 98 Hyades dwarf stars for 5 years, we have selected all stars that show low-amplitude long-period trends. With AO-observations of these 14 stars we found companion candidates around nine of them, where one star has two companions. The two companions of HIP 16548 have masses between 0.07 to 0.08 Mo, and are thus either brown dwarfs or very low mass stars. In the case of HAN 172 we found a companion with a mass between 0.08 to 0.10 Mo, which is again between a star and a brown dwarf. The other seven stars all have stellar companions. In two additional cases, the RV-variations are presumably caused by stellar activity, and in another case the companion could be a short-period binary. The images of the remaining two stars are slightly elongated, which might imply that even these are binaries. Because at least 12 of the 14 stars showing low-amplitude RV trends turn out to have companions with a mass greater than 70 MJupiter, or are just active, we finally estimate the number of companions with masses between 10 MJupiter and 70 MJupiter within 8 AU of the host stars in the Hyades as less equal 2%.

Most current astronomical adaptive optics (AO) systems rely on the availability of a bright star to measure the distortion of the incoming wavefront. Replacing the guide star with an artificial laser beacon alleviates this dependency on bright stars and therefore increases sky coverage, but it does not eliminate another serious problem for AO observations. This is the issue of PSF variation with time and field position near the guide star. In fact, because a natural guide star is still necessary for correction of the low-order phase error, characterization of laser guide star (LGS) AO PSF spatial variation is more complicated than for a natural guide star alone. We discuss six methods for characterizing LGS AO PSF variation that can potentially improve the determination of the PSF away from the laser spot, that is, off-axis. Calibration images of dense star fields are used to determine the change in PSF variation with field position. This is augmented by AO system telemetry and simple computer simulations to determine a more accurate off-axis PSF. We report on tests of the methods using the laser AO system on the Lick Observatory Shane Telescope. [Abstract truncated.]

AB Dor C, a low-mass companion to the young star AB Dor A has recently been imaged. Direct detection and astrometric observations have been used to solve the AB Dor C orbit and obtain a dynamical mass estimate. Near-infrared magnitudes and effective temperatures from models have then been compared with observed values. Models and observations disagree at more than the 2 sigma level. This result prompts Close et al. 2005 to claim that models are overestimating young low-mass object luminosities by roughly a factor of two. This claim is based on the hypothesis that the detected source, AB Dor C, is a single object. Another possible interpretation of the data that was not considered is that AB Dor C is an unresolved binary brown dwarf. Considering that 21 brown dwarfs companions have been found to date around stars and that three of them are confirmed or potential brown dwarf binaries (Gl569Ba & Bb, GL564B & C and Indi Ba & Bb), this suggests that brown dwarf binaries may constitute a significant fraction (14% +/- 8%) of the brown dwarf population around stars. For these reasons, the possibility that AB Dor C is an unresolved binary brown dwarf deserves attention.

We present results of near infrared imaging of the disk-galaxy-merger NGC 6240 using adaptive optics on the Keck II Telescope and reprocessed archival data from NICMOS on the Hubble Space Telescope. Both the North and South nuclei of NGC 6240 are clearly elongated, with considerable sub-structure within each nucleus. In K' band there are at least two point-sources within the North nucleus; we tentatively identify the south-western point-source within the North nucleus as the position of one of the two AGNs. Within the South nucleus, the northern sub-nucleus is more highly reddened. Based upon the nuclear separation measured at 5 GHz, we suggest that the AGN in the South nucleus is still enshrouded in dust at K' band, and is located slightly to the north of the brightest point in K' band. Within the South nucleus there is strong H2 1-0 S(1) line emission from the northern sub-nucleus, contrary to the conclusions of previous seeing-limited observations. Narrowband H2 emission-line images show that a streamer or ribbon of excited molecular hydrogen connects the North and South nuclei. We suggest that this linear feature corresponds to a bridge of gas connecting the two nuclei, as seen in computer simulations of mergers. Many point-like regions are seen around the two nuclei. These are most prominent at 1.1 microns with NICMOS, and in K'-band with Keck adaptive optics. We suggest that these point-sources represent young star clusters formed in the course of the merger.

The results of a comprehensive infrared imaging search for the putative 0.06 solar mass astrometric companion to the 4.4 pc white dwarf van Mannen 2 are reported. Adaptive optics images acquired at 3.8 microns reveal a diffraction limited core of 0.09" and no direct evidence of a secondary. Models predict that at 5 Gyr, a 50 jupiter mass brown dwarf would be only 1 magnitude fainter than van Maanen 2 at this wavelength and the astrometric analysis suggested a separation of 0.2". In the case of a chance alignment along the line of sight, a 0.4 mag excess should be measured. An independent photometric observation at the same wavelength reveals no excess. In addition, there exist published ISO observations of van Maanen 2 at 6.8 and 15.0 microns which are consistent with photospheric flux of a 6750 K white dwarf. If recent brown dwarf models are correct, there is no substellar companion with T > 500 K.

We present the first 3.8 micron image of the dusty ring surrounding the young binary system GG Tau, obtained with the W. M. Keck II 10m telescope's adaptive optics system. This is the longest wavelength at which the ring has been detected in scattered light so far, allowing a multi-wavelength analysis of the scattering properties of the dust grains present in this protoplanetary disk in combination with previous, shorter wavelengths, HST images. We find that the scattering phase function of the dust grains in the disk is only weakly dependent on the wavelength. This is inconsistent with dust models inferred from observations of the interstellar medium or dense molecular clouds. In particular, the strongly forward-throwing scattering phase function observed at 3.8 micron implies a significant increase in the population of large (>~1 micron) grains, which provides direct evidence for grain growth in the ring. However, the grain size distribution required to match the 3.8 micron image of the ring is incompatible with its published 1 micron polarization map, implying that the dust population is not uniform throughout the ring. We also show that our 3.8 micron scattered light image probes a deeper layer of the ring than previous shorter wavelength images, as demonstrated by a shift in the location of the inner edge of the disk's scattered light distribution between 1 and 3.8 micron. We therefore propose a stratified structure for the ring in which the surface layers, located ~50 AU above the ring midplane, contain dust grains that are very similar to those found in dense molecular clouds, while the region of the ring located ~25 AU from the midplane contains significantly larger grains. This stratified structure is likely the result of vertical dust settling and/or preferred grain growth in the densest parts of the ring.

HD 199143 and HD 358623 (BD-17 6128) are two sets of binary stars which are physically associated and 48 pc from Earth. We present heliocentric radial velocities and high lithium abundances which establish these stars as members of the 12 Myr-old Beta Pictoris Moving Group. We also present mid-IR photometric measurements which show no firm evidence for warm dust around all four stars.

A significant population of nearby stars have strong far-infrared excesses, now known to be due to circumstellar dust in regions analogous to the Kuiper Belt of our solar system, though orders of magnitude more dense. Recent sub-mm and mm imaging of these systems resolves the circumstellar dust and reveals complex structures, often in the form of rings with azimuthal non-axisymmetric variations. This structure might well be due to the presence of embedded brown dwarfs or planets. We have carried out deep adaptive optics imaging of two nearby stars with such asymmetric dust: Epsilon Eridani and Vega. Ten and seven candidate companions were seen in and near the dust rings of Epsilon Eridani and Vega respectively, but second-epoch proper motion measurements indicate that all are background objects. Around these two stars we can thus exclude planetary companions at spatial scales comparable to the radius of the dust structures to a level of MK=24, corresponding to 5 Jupiter masses, for Epsilon Eridani, and MK=19-21, corresponding to 6-8 Jupiter masses, for Vega.

The useful dynamic range of an image in the diffraction limited regime is usually limited by speckles caused by residual phase errors in the optical system forming the image. The technique of speckle decorrelation involves introducing many independent realizations of additional phase error into a wavefront during one speckle lifetime, changing the instantaneous speckle pattern. A commonly held assumption is that this results in the speckles being `moved around' at the rate at which the additional phase screens are applied. The intention of this exercise is to smooth the speckles out into a more uniform background distribution during their persistence time, thereby enabling companion detection around bright stars to be photon noise limited rather than speckle-limited. We demonstrate analytically why this does not occur, and confirm this result with numerical simulations. We show that the original speckles must persist, and that the technique of speckle decorrelation merely adds more noise to the original speckle noise, thereby degrading the dynamic range of the image.

Adaptive optics, large primary mirrors, and careful selection of target stars are the keys to ground-based imaging of extrasolar planets. Our near-IR survey is capable of identifying exoplanets of 1-10 M${_J}$ within 100 AU of young (t$<$60 Myr), nearby (d$<$60 pc) stars. For very young and proximate targets such as GJ 803 (12 Myr, 10 pc) we are able to detect exoplanets with parameters approaching those of planets our solar system (1 M${_J}$ at 19 AU, 2 M${_J}$ at 9 AU; 5 M${_J}$ at 5 AU). We have thus far imaged over 100 stars with Keck AO. Here we report on our progress, discuss specific observing strategies, and present detailed sensitivity limits.

A combination of high-resolution and wide-field imaging reveals two binary stars and one triple star system among the sample of the first 11 stars with planets detected by radial velocity variations. High resolution speckle or adaptive optics (AO) data probe subarcsecond scales down to the diffraction limit of the Keck 10-m or Lick 3-m, and direct images or AO images are sensitive to a wider field, extending to 10" or 38", depending upon the camera. One of the binary systems -- HD 114762 -- was not previously known to be a spatially resolved multiple system; additional data taken with the combination of Keck adaptive optics and NIRSPEC are used to characterize the new companion. The second binary system -- Tau Boo -- was a known multiple with two conflicting orbital solutions; the current data will help constrain the discrepant estimates of periastron time and separation. Another target -- 16 Cyg B -- was a known common proper motion binary, but the current data resolve a new third component, close to the wide companion 16 Cyg A. Both the HD 114762 and 16 Cyg B systems harbor planets in eccentric orbits, while the Tau Boo binary contains an extremely close planet in a tidally-circularized orbit. Although the sample is currently small, the proportion of binary systems is comparable to that measured in the field over a similar separation range. Incorporating the null result from another companion search project lowers the overall fraction of planets in binary systems, but the detections in our survey reveal that planets can form in binaries separated by less than 50 AU.

Even though the technology of adaptive optics (AO) is rapidly maturing, calibration of the resulting images remains a major challenge. The AO point-spread function (PSF) changes quickly both in time and position on the sky. In a typical observation the star used for guiding will be separated from the scientific target by 10" to 30". This is sufficient separation to render images of the guide star by themselves nearly useless in characterizing the PSF at the off-axis target position. A semi-empirical technique is described that improves the determination of the AO off-axis PSF. The method uses calibration images of dense star fields to determine the change in PSF with field position. It then uses this information to correct contemporaneous images of the guide star to produce a PSF that is more accurate for both the target position and the time of a scientific observation. We report on tests of the method using natural-guide-star AO systems on the Canada-France-Hawaii Telescope and Lick Observatory Shane Telescope, augmented by simple atmospheric computer simulations. At 25" off-axis, predicting the PSF full width at half maximum using only information about the guide star results in an error of 60%. Using an image of a dense star field lowers this error to 33%, and our method, which also folds in information about the on-axis PSF, further decreases the error to 19%.

Using adaptive optics on the W.M. Keck II telescope we imaged Titan several times during 1999 to 2001 in narrowband near-infrared filters selected to probe Titan's stratosphere and upper troposphere. We observed a bright feature around the south pole, possibly a collar of haze or clouds. Further, we find that solar phase angle explains most of the observed east-west brightness asymmetry of Titan's atmosphere, although the data do not preclude the presence of a `morning fog' effect at small solar phase angle.

Adaptive optics (AO) on 8-10 m telescopes is an enormously powerful tool for studying young nearby stars. It is especially useful for searching for companions. Using AO on the 10-m W.M. Keck II telescope we have measured the position of the brown dwarf companion to TWA5 and resolved the primary into an 0.055 arcsecond double. Over the next several years follow-up astrometry should permit an accurate determination of the masses of these young stars. We have also re-observed the candidate extrasolar planet TWA6B, but measurements of its motion relative to TWA6A are inconclusive. We are carrying out a search for new planetary or brown dwarf companions to TWA stars and, if current giant planet models are correct, are currently capable of detecting a 1 Jupiter-mass companion at ~1 arcsecond and a 5 Jupiter-mass companion at ~0.5 arcsecon around a typical TWA member.

We have undertaken an adaptive optics (AO) imaging survey of extrasolar planetary systems and stars showing interesting radial velocity trends from high precision radial velocity searches. AO increases the resolution and dynamic range of an image, substantially improving the detectability of faint close companions. This survey is sensitive to objects less luminous than the bottom of the main sequence at separations as close as 1 arcsec. We have detected stellar companions to the planet bearing stars HD 114762 and Tau Boo. We have also detected a companion to the non-planet bearing star 16 Cyg A.