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 place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.

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.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

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.

The orbital eccentricities of directly imaged exoplanets and brown dwarf companions provide clues about their formation and dynamical histories. We combine new high-contrast imaging observations of substellar companions obtained primarily with Keck/NIRC2 together with astrometry from the literature to test for differences in the population-level eccentricity distributions of 27 long-period giant planets and brown dwarf companions between 5-100 AU using hierarchical Bayesian modeling. Orbit fits are performed in a uniform manner for companions with short orbital arcs; this typically results in broad constraints for individual eccentricity distributions, but together as an ensemble these systems provide valuable insight into their collective underlying orbital patterns. The shape of the eccentricity distribution function for our full sample of substellar companions is approximately flat from e=0-1. When subdivided by companion mass and mass ratio, the underlying distributions for giant planets and brown dwarfs show significant differences. Low mass ratio companions preferentially have low eccentricities, similar to the orbital properties of warm Jupiters found with radial velocities and transits. We interpret this as evidence for in situ formation on largely undisturbed orbits within massive, extended disks. Brown dwarf companions exhibit a broad peak at e $\approx$ 0.6-0.9 with evidence for a dependence on orbital period. This closely resembles the orbital properties and period-eccentricity trends of wide (1-200 AU) stellar binaries, suggesting that brown dwarfs in this separation range predominantly form in a similar fashion. We also report evidence that the "eccentricity dichotomy" observed at small separations extends to planets on wide orbits: the mean eccentricity for the multi-planet system HR 8799 is lower than for systems with single planets. (Abridged)

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

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

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

We 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.

orbitize! is an open-source, object-oriented software package for fitting the orbits of directly-imaged objects. It packages the Orbits for the Impatient (OFTI) algorithm and a parallel-tempered Markov Chain Monte Carlo (MCMC) algorithm into a consistent and intuitive Python API. orbitize! makes it easy to run standard astrometric orbit fits; in less than 10 lines of code, users can read in data, perform one fit using OFTI and another using MCMC, and make two publication-ready figures. Extensive pedagogical tutorials, intended to be navigable by both orbit-fitting novices and seasoned experts, are available on our documentation page. We have designed the orbitize! API to be flexible and easy to use/modify for unique cases. orbitize! was designed by members of the exoplanet imaging community to be a central repository for algorithms, techniques, and know-how developed by this community. We intend for it to continue to expand and change as the field progresses and new techniques are developed, and call for community involvement in this process. Complete and up-to-date documentation is available at orbitize.info.

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.

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.

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 Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein.

Gas giants are the most readily detectable exoplanets but fundamental questions about their system architectures, formation, migration, and atmospheres have been unanswerable with the current generation of ground- and space-based facilities. The dominant techniques to detect and characterize giant planets $-$ radial velocities, transits, direct imaging, microlensing, and astrometry $-$ are each isolated to a limited range of planet masses, separations, ages, and temperatures. These windows into the arrangement and physical properties of giant planets have spawned new questions about the timescale and location of their assembly; the distributions of planet mass and orbital separation at young and old ages; the composition and structure of their atmospheres; and their orbital and rotational angular momentum architectures. The ELTs will address these questions by building bridges between these islands of mass, orbital distance, and age. The angular resolution, collecting area, all-sky coverage, and novel instrumentation suite of these facilities are needed to provide a complete map of the orbits and atmospheric evolution of gas giant planets (0.3$-$10 $M_\mathrm{Jup}$) across space (0.1$-$100 AU) and time (1 Myr to 10 Gyr). This white paper highlights the scientific potential of the GMT and TMT to address these outstanding questions, with a particular focus on the role of direct imaging and spectroscopy of large samples of giant planets that will soon be made available with $Gaia$.

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.

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources.

We present SCExAO/CHARIS high-contrast imaging/$JHK$ integral field spectroscopy of $\kappa$ And b, a directly-imaged low-mass companion orbiting a nearby B9V star. We detect $\kappa$ And b at a high signal-to-noise and extract high precision spectrophotometry using a new forward-modeling algorithm for (A-)LOCI complementary to KLIP-FM developed by Pueyo (2016). $\kappa$ And b's spectrum best resembles that of a low-gravity L0--L1 dwarf (L0--L1$\gamma$). Its spectrum and luminosity are very well matched by 2MASSJ0141-4633 and several other 12.5--15 $M_{\rm J}$ free floating members of the 40 $Myr$-old Tuc-Hor Association, consistent with a system age derived from recent interferometric results for the primary, a companion mass at/near the deuterium-burning limit (13$^{+12}_{-2}$ M$_{\rm J}$), and a companion-to-primary mass ratio characteristic of other directly-imaged planets ($q$ $\sim$ 0.005$^{+0.005}_{-0.001}$). We did not unambiguously identify additional, more closely-orbiting companions brighter and more massive than $\kappa$ And b down to $\rho$ $\sim$ 0.3" (15 au). SCExAO/CHARIS and complementary Keck/NIRC2 astrometric points reveal clockwise orbital motion. Modeling points towards a likely eccentric orbit: a subset of acceptable orbits include those that are aligned with the star's rotation axis. However, $\kappa$ And b's semimajor axis is plausibly larger than 75 au and in a region where disk instability could form massive companions. Deeper $\kappa$ And high-contrast imaging and low-resolution spectroscopy from extreme AO systems like SCExAO/CHARIS and higher resolution spectroscopy from Keck/OSIRIS or, later, IRIS on the Thirty Meter Telescope could help clarify $\kappa$ And b's chemistry and whether its spectrum provides an insight into its formation environment.

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.

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 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 Exoplanet Survey (GPIES) is a direct imaging campaign designed to search for young, self-luminous, giant exoplanets. To date, GPIES has observed nearly 500 targets, and generated over 30,000 individual exposures using its integral field spectrograph (IFS) instrument. The GPIES team has developed a campaign data system with a database incorporating all of the metadata for all individual raw data products, including environmental conditions and instrument performance metrics. The same database also indexes metadata associated with multiple levels of reduced data products, including contrast measures for individual images and combined image sequences, which serve as the primary metric of performance for the final science products. The database is also used to track telemetry products from the adaptive optics subsystem, and associate these with corresponding IFS data. Here, we discuss several data exploration and visualization projects enabled by the GPIES database. Of particular interest are any correlations between instrument performance and environmental or operating conditions. We show single and multiple-parameter fits of single-image and observing sequence contrast as functions of various seeing measures, and discuss automated outlier rejection and other fitting concerns. Supervised learning techniques are employed in order to partition the space of raw (single image) to final (full sequence) contrast in order to better predict the value of the final data set from the first few completed observations. Finally, we discuss the particular features of the database design that aid in performing these analyses, and suggest potential future upgrades and refinements.

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.

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.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

The occurrence rate of young giant planets from direct imaging surveys is a fundamental tracer of the efficiency with which planets form and migrate at wide orbital distances. These measurements have progressively converged to a value of about 1% for the most massive planets ($\approx$5-13 $M_\mathrm{Jup}$) averaged over all stellar masses at separations spanning a few tens to a few hundreds of AU. The subtler statistical properties of this population are beginning to emerge with ever-increasing sample sizes: there is tentative evidence that planets on wide orbits are more frequent around stars that possess debris disks; brown dwarf companions exist at comparable (or perhaps slightly higher) rates as their counterparts in the planetary-mass regime; and the substellar companion mass function appears to be smooth and may extend down to the opacity limit for fragmentation. Within a few years, the conclusion of second-generation direct imaging surveys will enable more definitive interpretations with the ultimate goal of identifying the dominant origin of this population and uncovering its relationship to planets at smaller separations.

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

Our Keck/NIRC2 imaging survey searches for stellar companions around 144 systems with radial velocity (RV) detected giant planets to determine whether stellar binaries influence the planets' orbital parameters. This survey, the largest of its kind to date, finds eight confirmed binary systems and three confirmed triple systems. These include three new multi-stellar systems (HD 30856, HD 86081, and HD 207832) and three multi-stellar systems with newly confirmed common proper motion (HD 43691, HD 116029, and HD 164509). We combine these systems with seven RV planet-hosting multi-stellar systems from the literature in order to test for differences in the properties of planets with semimajor axes ranging between 0.1-5 au in single vs multi-stellar systems. We find no evidence that the presence or absence of stellar companions alters the distribution of planet properties in these systems. Although the observed stellar companions might influence the orbits of more distant planetary companions in these systems, our RV observations currently provide only weak constraints on the masses and orbital properties of planets beyond 5 au. In order to aid future efforts to characterize long period RV companions in these systems, we publish our contrast curves for all 144 targets. Using four years of astrometry for six hierarchical triple star systems hosting giant planets, we fit the orbits of the stellar companions in order to characterize the orbital architecture in these systems. We find that the orbital plane of the secondary and tertiary companions are inconsistent with an edge-on orbit in four out of six cases.

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

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

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

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

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.

We have conducted an angular differential imaging survey with NIRC2 at Keck in search of close-in substellar companions to a sample of seven systems with confirmed planetary-mass companions (PMCs) on wide orbits (>50 AU). These wide-separation PMCs pose significant challenges to all three possible formation mechanisms: core accretion plus scattering, disk instability, and turbulent fragmentation. We explore the possibility that these companions formed closer in and were scattered out to their present-day locations by searching for other massive bodies at smaller separations. The typical sensitivity for this survey is \Delta K ~ 12.5 at 1". We identify eight candidate companions, whose masses would reach as low as one Jupiter mass if gravitationally bound. From our multi-epoch astrometry we determine that seven of these are conclusively background objects, while the eighth near DH Tau is ambiguous and requires additional monitoring. We rule out the presence of >7 M$_{\rm Jup}$ bodies in these systems down to 15 - 50 AU that could be responsible for scattering. This result combined with the totality of evidence suggests that dynamical scattering is unlikely to have produced this population of PMCs. We detect orbital motion from the companions ROXs 42B b and ROXs 12 b, and from this determine 95% upper limits on the companions' eccentricities of 0.58 and 0.83 respectively. Finally, we find that the 95% upper limit on the occurrence rate of additional planets with masses between 5 - 15 M$_{\rm Jup}$ outside of 40 AU in systems with PMCs is 54%.

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

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

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

We present the discovery of a new satellite of asteroid (130) Elektra - S/2014 (130) 1 - in differential imaging and in integral field spectroscopy data over multiple epochs obtained with SPHERE/VLT. This new (second) moonlet of Elektra is about 2 km across, on an eccentric orbit and about 500 km away from the primary. For a comparative study, we also observed another triple asteroid system (93) Minerva. For both systems, component-resolved reflectance spectra of the satellites and primary were obtained simultaneously. No significant spectral difference was observed between the satellites and the primary for either triple system. We find that the moonlets in both systems are more likely to have been created by sub-disruptive impacts as opposed to having been captured.

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.

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

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

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

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

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

We present optical and near-infrared adaptive optics (AO) imaging and spectroscopy of 13 ultracool (>M6) companions to late-type stars (K7-M4.5), most of which have recently been identified as candidate members of nearby young moving groups (YMGs; 8-120 Myr) in the literature. The inferred masses of the companions (~10-100 Mjup) are highly sensitive to the ages of the primary stars so we critically examine the kinematic and spectroscopic properties of each system to distinguish bona fide YMG members from old field interlopers. 2MASS J02155892-0929121 C is a new M7 substellar companion (40-60 Mjup) with clear spectroscopic signs of low gravity and hence youth. The primary, possibly a member of the ~40 Myr Tuc-Hor moving group, is visually resolved into three components, making it a young low-mass quadruple system in a compact (<100 AU) configuration. In addition, Li 1 $\lambda$6708 absorption in the intermediate-gravity M7.5 companion 2MASS J15594729+4403595 B provides unambiguous evidence that it is young (<200 Myr) and resides below the hydrogen burning limit. Three new close-separation (<1") companions (2MASS J06475229-2523304 B, PYC J11519+0731 B, and GJ 4378 Ab) orbit stars previously reported as candidate YMG members, but instead are likely old (>1 Gyr) tidally-locked spectroscopic binaries without convincing kinematic associations with any known moving group. The high rate of false positives in the form of old active stars with YMG-like kinematics underscores the importance of radial velocity and parallax measurements to validate candidate young stars identified via proper motion and activity selection alone. Finally, we spectroscopically confirm the cool temperature and substellar nature of HD 23514 B, a recently discovered M8 benchmark brown dwarf orbiting the dustiest-known member of the Pleiades. [Abridged]

We report here the highest resolution near-IR imaging to date of the HD 141569A disc taken as part of the NICI Science Campaign. We recover 4 main features in the NICI images of the HD 141569 disc discovered in previous HST imaging: 1) an inner ring / spiral feature. Once deprojected, this feature does not appear circular. 2) an outer ring which is considerably brighter on the western side compared to the eastern side, but looks fairly circular in the deprojected image. 3) an additional arc-like feature between the inner and outer ring only evident on the east side. In the deprojected image, this feature appears to complete the circle of the west side inner ring and 4) an evacuated cavity from 175 AU inwards. Compared to the previous HST imaging with relatively large coronagraphic inner working angles (IWA), the NICI coronagraph allows imaging down to an IWA of 0.3". Thus, the inner edge of the inner ring/spiral feature is well resolved and we do not find any additional disc structures within 175 AU. We note some additional asymmetries in this system. Specifically, while the outer ring structure looks circular in this deprojection, the inner bright ring looks rather elliptical. This suggests that a single deprojection angle is not appropriate for this system and that there may be an offset in inclination between the two ring / spiral features. We find an offset of 4+-2 AU between the inner ring and the star center, potentially pointing to unseen inner companions.

We present a new algorithm designed to improve the signal to noise ratio (SNR) of point and extended source detections in direct imaging data. The novel part of our method is that it finds the linear combination of the science images that best match counterpart images with signal removed from suspected source regions. The algorithm, based on the Locally Optimized Combination of Images (LOCI) method, is called Matched LOCI or MLOCI. We show using data obtained with the Gemini Planet Imager (GPI) and Near-Infrared Coronagraphic Imager (NICI) that the new algorithm can improve the SNR of point source detections by 30-400% over past methods. We also find no increase in false detections rates. No prior knowledge of candidate companion locations is required to use MLOCI. While non-blind applications may yield linear combinations of science images which seem to increase the SNR of true sources by a factor > 2, they can also yield false detections at high rates. This is a potential pitfall when trying to confirm marginal detections or to re-detect point sources found in previous epochs. Our findings are relevant to any method where the coefficients of the linear combination are considered tunable, e.g. LOCI and Principal Component Analysis (PCA). Thus we recommend that false detection rates be analyzed when using these techniques.

We describe the astrometric calibration of the Gemini NICI Planet-Finding Campaign. The Campaign requires a relative astrometric accuracy of $\approx$ 20 mas across multi-year timescales in order to distinguish true companions from background stars by verifying common proper motion and parallax with their parent stars. The calibration consists of a correction for instrumental optical image distortion, plus on-sky imaging of astrometric fields to determine the pixel scale and image orientation. We achieve an accuracy of $\lesssim 7$ mas between the center and edge of the 18$''$ NICI field, meeting the 20 mas requirement. Most of the Campaign data in the Gemini Science Archive are accurate to this level but we identify a number of anomalies and present methods to correct the errors.

We present new astrometry for the young (12--21 Myr) exoplanet beta Pictoris b taken with the Gemini/NICI and Magellan/MagAO instruments between 2009 and 2012. The high dynamic range of our observations allows us to measure the relative position of beta Pic b with respect to its primary star with greater accuracy than previous observations. Based on a Markov Chain Monte Carlo analysis, we find the planet has an orbital semi-major axis of 9.1 (+5.3, -0.5) AU and orbital eccentricity <0.15 at 68% confidence (with 95% confidence intervals of 8.2--48 AU and 0.00--0.82 for semi-major axis and eccentricity, respectively, due to a long narrow degenerate tail between the two). We find that the planet has reached its maximum projected elongation, enabling higher precision determination of the orbital parameters than previously possible, and that the planet's projected separation is currently decreasing. With unsaturated data of the entire beta Pic system (primary star, planet, and disk) obtained thanks to NICI's semi-transparent focal plane mask, we are able to tightly constrain the relative orientation of the circumstellar components. We find the orbital plane of the planet lies between the inner and outer disks: the position angle (PA) of nodes for the planet's orbit (211.8 +/- 0.3 degrees) is 7.4 sigma greater than the PA of the spine of the outer disk and 3.2 sigma less than the warped inner disk PA, indicating the disk is not collisionally relaxed. Finally, for the first time we are able to dynamically constrain the mass of the primary star beta Pic to 1.76 (+0.18, -0.17) solar masses.

We present J, H, CH_4 short (1.578 micron), CH_4 long (1.652 micron) and K_s-band images of the dust ring around the 10 Myr old star HR 4796 A obtained using the Near Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1 meter Telescope. Our images clearly show for the first time the position of the star relative to its circumstellar ring thanks to NICI's translucent focal plane occulting mask. We employ a Bayesian Markov Chain Monte Carlo method to constrain the offset vector between the two. The resulting probability distribution shows that the ring center is offset from the star by 16.7+/-1.3 milliarcseconds along a position angle of 26+/-3 degrees, along the PA of the ring, 26.47+/-0.04 degrees. We find that the size of this offset is not large enough to explain the brightness asymmetry of the ring. The ring is measured to have mostly red reflectivity across the JHK_s filters, which seems to indicate micron-sized grains. Just like Neptune's 3:2 and 2:1 mean-motion resonances delineate the inner and outer edges of the classical Kuiper Belt, we find that the radial extent of the HR 4796 A and Fomalhaut rings could correspond to the 3:2 and 2:1 mean-motion resonances of hypothetical planets at 54.7 AU and 97.7 AU in the two systems, respectively. A planet orbiting HR 4796 A at 54.7 AU would have to be less massive than 1.6 Mjup so as not to widen the ring too much by stirring.

We present the first ground-based CCD ($\lambda < 1\mu$m) image of an extrasolar planet. Using MagAO's VisAO camera we detected the extrasolar giant planet (EGP) $\beta$ Pictoris b in $Y$-short ($Y_S$, 0.985 $\mu$m), at a separation of $0.470 \pm 0.010''$ and a contrast of $(1.63 \pm 0.49) \times 10^{-5}$. This detection has a signal-to-noise ratio of 4.1, with an empirically estimated upper-limit on false alarm probability of 1.0%. We also present new photometry from the NICI instrument on the Gemini-South telescope, in $CH_{4S,1\%}$ ($1.58$ $\mu m$), $K_S$ ($2.18\mu m$), and $K_{cont}$ (2.27 $\mu m$). A thorough analysis of our photometry combined with previous measurements yields an estimated near-IR spectral type of L$2.5\pm1.5$, consistent with previous estimates. We estimate log$(L_{bol}/L_{Sun})$ = $-3.86 \pm 0.04$, which is consistent with prior estimates for $\beta$ Pic b and with field early-L brown dwarfs. This yields a hot-start mass estimate of $11.9 \pm 0.7$ $M_{Jup}$ for an age of $21\pm4$ Myr, with an upper limit below the deuterium burning mass. Our $L_{bol}$ based hot-start estimate for temperature is $T_{eff}=1643\pm32$ K (not including model dependent uncertainty). Due to the large corresponding model-derived radius of $R=1.43\pm0.02$ $R_{Jup}$, this $T_{eff}$ is $\sim$$250$ K cooler than would be expected for a field L2.5 brown dwarf. Other young, low-gravity (large radius), ultracool dwarfs and directly-imaged EGPs also have lower effective temperatures than are implied by their spectral types. However, such objects tend to be anomalously red in the near-IR compared to field brown dwarfs. In contrast, $\beta$ Pic b has near-IR colors more typical of an early-L dwarf despite its lower inferred temperature.

We present the high-contrast image processing techniques used by the Gemini NICI Planet-Finding Campaign to detect faint companions to bright stars. NICI (Near Infrared Coronagraphic Imager) is an adaptive optics instrument installed on the 8-m Gemini South telescope, capable of angular and spectral difference imaging and specifically designed to image exoplanets. The Campaign data pipeline achieves median contrasts of 12.6 magnitudes at 0.5" and 14.4 magnitudes at 1" separation, for a sample of 45 stars (V= 4.3-13.9 mag) from the early phase of the Campaign. We also present a novel approach to calculating contrast curves for companion detection based on 95% completeness in the recovery of artificial companions injected into the raw data, while accounting for the false-positive rate. We use this technique to select the image processing algorithms that are more successful at recovering faint simulated point sources. We compare our pipeline to the performance of the LOCI algorithm for NICI data and do not find significant improvement with LOCI.

We report results of a direct imaging survey for giant planets around 80 members of the Beta Pic, TW Hya, Tucana-Horologium, AB Dor, and Hercules-Lyra moving groups, observed as part of the Gemini NICI Planet-Finding Campaign. For this sample, we obtained median contrasts of \Delta H=13.9 mag at 1" in combined CH4 narrowband ADI+SDI mode and median contrasts of \Delta H=15.1 mag at 2" in H-band ADI mode. We found numerous (>70) candidate companions in our survey images. Some of these candidates were rejected as common-proper motion companions using archival data; we reobserved with NICI all other candidates that lay within 400 AU of the star and were not in dense stellar fields. The vast majority of candidate companions were confirmed as background objects from archival observations and/or dedicated NICI campaign followup. Four co-moving companions of brown dwarf or stellar mass were discovered in this moving group sample: PZ Tel B (36+-6 MJup, 16.4+-1.0 AU, Biller et al. 2010), CD -35 2722B (31+-8 MJup, 67+-4 AU, Wahhaj et al. 2011), HD 12894B (0.46+-0.08 MSun, 15.7+-1.0 AU), and BD+07 1919C (0.20+-0.03 MSun, 12.5+-1.4 AU). From a Bayesian analysis of the achieved H band ADI and ASDI contrasts, using power-law models of planet distributions and hot-start evolutionary models, we restrict the frequency of 1--20 MJup companions at semi-major axes from 10--150 AU to <18% at a 95.4% confidence level using DUSTY models and to <6% at a 95.4% using COND models.

We have completed a high-contrast direct imaging survey for giant planets around 57 debris disk stars as part of the Gemini NICI Planet-Finding Campaign. We achieved median H-band contrasts of 12.4 mag at 0.5" and 14.1 mag at 1" separation. Follow-up observations of the 66 candidates with projected separation < 500 AU show that all of them are background objects. To establish statistical constraints on the underlying giant planet population based on our imaging data, we have developed a new Bayesian formalism that incorporates (1) non-detections, (2) single-epoch candidates, (3) astrometric and (4) photometric information, and (5) the possibility of multiple planets per star to constrain the planet population. Our formalism allows us to include in our analysis the previously known Beta Pictoris and the HR 8799 planets. Our results show at 95% confidence that <13% of debris disk stars have a >5MJup planet beyond 80 AU, and <21% of debris disk stars have a >3MJup planet outside of 40 AU, based on hot-start evolutionary models. We model the population of directly-imaged planets as d^2N/dMda ~ m^alpha a^beta, where m is planet mass and a is orbital semi-major axis (with a maximum value of amax). We find that beta < -0.8 and/or alpha > 1.7. Likewise, we find that beta < -0.8 and/or amax < 200 AU. If we ignore the Beta Pic and HR 8799 planets (should they belong to a rare and distinct group), we find that < 20% of debris disk stars have a > 3MJup planet beyond 10 AU, and beta < -0.8 and/or alpha < -1.5. Our Bayesian constraints are not strong enough to reveal any dependence of the planet frequency on stellar host mass. Studies of transition disks have suggested that about 20% of stars are undergoing planet formation; our non-detections at large separations show that planets with orbital separation > 40 AU and planet masses > 3 MJup do not carve the central holes in these disks.

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (~1.5-2.5 M_sun) conducted to date and includes the planet hosts beta Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58 (+21, -20) M_Jup and 55 (+20, -19) M_Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M_sun stars can have giant planets greater than 4 M_Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M_Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M_Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

We report the discovery by the Gemini NICI Planet-Finding Campaign of two low-mass companions to the young A0V star HD 1160 at projected separations of 81 +/- 5 AU (HD 1160 B) and 533 +/- 25 AU (HD 1160 C). VLT images of the system taken over a decade for the purpose of using HD 1160 A as a photometric calibrator confirm that both companions are physically associated. By comparing the system to members of young moving groups and open clusters with well-established ages, we estimate an age of 50 (+50,-40) Myr for HD 1160 ABC. While the UVW motion of the system does not match any known moving group, the small magnitude of the space velocity is consistent with youth. Near-IR spectroscopy shows HD 1160 C to be an M3.5 +/- 0.5 star with an estimated mass of 0.22 (+0.03,-0.04) M_Sun, while NIR photometry of HD 1160 B suggests a brown dwarf with a mass of 33 (+12,-9) M_Jup. The very small mass ratio (0.014) between the A and B components of the system is rare for A star binaries, and would represent a planetary-mass companion were HD 1160 A to be slightly less massive than the Sun.

We present the discovery of a wide (67 AU) substellar companion to the nearby (21 pc) young solar-metallicity M1 dwarf CD-35 2722, a member of the ~100 Myr AB Doradus association. Two epochs of astrometry from the NICI Planet-Finding Campaign confirm that CD-35 2722 B is physically associated with the primary star. Near-IR spectra indicate a spectral type of L4\pm1 with a moderately low surface gravity, making it one of the coolest young companions found to date. The absorption lines and near-IR continuum shape of CD-35 2722 B agree especially well the dusty field L4.5 dwarf 2MASS J22244381-0158521, while the near-IR colors and absolute magnitudes match those of the 5 Myr old L4 planetary-mass companion, 1RXS J160929.1-210524 b. Overall, CD-35 2722 B appears to be an intermediate-age benchmark for L-dwarfs, with a less peaked H-band continuum than the youngest objects and near-IR absorption lines comparable to field objects. We fit Ames-Dusty model atmospheres to the near-IR spectra and find T=1700-1900 K and log(g) =4.5\pm0.5. The spectra also show that the radial velocities of components A and B agree to within \pm10 km/s, further confirming their physical association. Using the age and bolometric luminosity of CD-35 2722 B, we derive a mass of 31\pm8 Mjup from the Lyon/Dusty evolutionary models. Altogether, young late-M to mid-L type companions appear to be over-luminous for their near-IR spectral type compared to field objects, in contrast to the under-luminosity of young late-L and early-T dwarfs.

Our team is carrying out a multi-year observing program to directly image and characterize young extrasolar planets using the Near-Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1-meter telescope. NICI is the first instrument on a large telescope designed from the outset for high-contrast imaging, comprising a high-performance curvature adaptive optics system with a simultaneous dual-channel coronagraphic imager. Combined with state-of-the-art observing methods and data processing, NICI typically achieves ~2 magnitudes better contrast compared to previous ground-based or space-based programs, at separations inside of ~2 arcsec. In preparation for the Campaign, we carried out efforts to identify previously unrecognized young stars, to rigorously construct our observing strategy, and to optimize the combination of angular and spectral differential imaging. The Planet-Finding Campaign is in its second year, with first-epoch imaging of 174 stars already obtained out of a total sample of 300 stars. We describe the Campaign's goals, design, implementation, performance, and preliminary results. The NICI Campaign represents the largest and most sensitive imaging survey to date for massive (~1 Mjup) planets around other stars. Upon completion, the Campaign will establish the best measurements to date on the properties of young gas-giant planets at ~5-10 AU separations. Finally, Campaign discoveries will be well-suited to long-term orbital monitoring and detailed spectrophotometric followup with next-generation planet-finding instruments.

We report the discovery of a tight substellar companion to the young solar analog PZ Tel, a member of the Beta Pictoris moving group observed with high contrast adaptive optics imaging as part of the Gemini NICI Planet-Finding Campaign. The companion was detected at a projected separation of 16.4 +/- 1.0 AU (0.33 +/- 0.01") in April 2009. Second-epoch observations in May 2010 demonstrate that the companion is physically associated and shows significant orbital motion. Monte Carlo modeling constrains the orbit of PZ Tel B to eccentricities > 0.6. The near-IR colors of PZ Tel B indicate a spectral type of M7+/-2 and thus this object will be a new benchmark companion for studies of ultracool, low-gravity photospheres. Adopting an age of 12 +8 -4 Myr for the system, we estimate a mass of 36 +/- 6 Mjup based on the Lyon/DUSTY evolutionary models. PZ Tel B is one of few young substellar companions directly imaged at orbital separations similar to those of giant planets in our own solar system. Additionally, the primary star PZ Tel A shows a 70 um emission excess, evidence for a significant quantity of circumstellar dust that has not been disrupted by the orbital motion of the companion.

We expand on the results of Nielsen et al. (2008), using the null result for giant extrasolar planets around the 118 target stars from the VLT NACO H and Ks band planet search (Masciadri et al. 2005), the VLT and MMT Simultaneous Differential Imaging (SDI) survey (Biller et al. 2007), and the Gemini Deep Planet Survey (Lafreniere et al. 2007) to set constraints on the population of giant extrasolar planets. Our analysis is extended to include the planet luminosity models of Fortney et al. (2008), as well as the correlation between stellar mass and frequency of giant planets found by Johnson et al. (2007). Doubling the sample size of FGKM stars strengthens our conclusions: a model for extrasolar giant planets with power-laws for mass and semi-major axis as giving by Cumming et al. (2008) cannot, with 95% confidence, have planets beyond 65 AU, compared to the value of 94 AU reported in Nielsen et al. (2008), using the models of Baraffe et al. (2003). When the Johnson et al. (2007) correction for stellar mass (which gives fewer Jupiter-mass companions to M stars with respect to solar-type stars) is applied, however, this limit moves out to 82 AU. For the relatively new Fortney et al. (2008) models, which predict fainter planets across most of parameter space, these upper limits, with and without a correction for stellar mass, are 182 and 234 AU, respectively.

We examine the implications for the distribution of extrasolar planets based on the null results from two of the largest direct imaging surveys published to date. Combining the measured contrast curves from 22 of the stars observed with the VLT NACO adaptive optics system by Masciadri et al. (2005), and 48 of the stars observed with the VLT NACO SDI and MMT SDI devices by Biller et al. (2007) (for a total of 60 unique stars; the median star for our survey is a 30 Myr K2 star at 25 pc), we consider what distributions of planet masses and semi-major axes can be ruled out by these data, based on Monte Carlo simulations of planet populations. We can set this upper limit with 95% confidence: the fraction of stars with planets with semi-major axis from 20 to 100 AU, and mass >4 M_Jup, is 20% or less. Also, with a distribution of planet mass of dN/dM ~ M^-1.16 between 0.5-13 M_Jup, we can rule out a power-law distribution for semi-major axis (dN/da ~ a^alpha) with index 0 and upper cut-off of 18 AU, and index -0.5 with an upper cut-off of 48 AU. For the distribution suggested by Cumming et al. (2007), a power-law of index -0.61, we can place an upper limit of 75 AU on the semi-major axis distribution. At the 68% confidence level, these upper limits state that fewer than 8% of stars have a planet of mass >4 M_Jup between 20 and 100 AU, and a power-law distribution for semi-major axis with index 0, -0.5, and -0.61 cannot have giant planets beyond 12, 23, and 29 AU, respectively. In general, we find that even null results from direct imaging surveys are very powerful in constraining the distributions of giant planets (0.5-13 M_Jup) at large separations, but more work needs to be done to close the gap between planets that can be detected by direct imaging, and those to which the radial velocity method is sensitive.

We present the results of a survey of 45 young (<250 Myr), close (<50 pc) stars with the Simultaneous Differential Imager (SDI) implemented at the VLT and the MMT for the direct detection of extrasolar planets. As part of the survey, we observed 54 objects total, of which 45 were close, young stars, 2 were somewhat more distant (<150 pc), very young (<10 Myr) stars, 3 were stars with known radial velocity planets, and 4 were older, very nearby (<20 pc) solar analogues. Our SDI devices use a double Wollaston prism and a quad filter to take images simultaneously at 3 wavelengths surrounding the 1.62 um methane absorption bandhead found in the spectrum of cool brown dwarfs and gas giant planets. In our VLT data, we achieved H band contrasts > 10 mag (5 sigma) at a separation of 0.5" from the primary star on 45% of our targets and H band contrasts of > 9 mag at a separation of 0.5'' on 80% of our targets. With this degree of attenuation, we should be able to image (5sigma detection) a 5 M_{Jup} planet 15 AU from a 70 Myr K1 star at 15 pc or a 5 M_{Jup} planet at 2 AU from a 12 Myr M star at 10 pc. We believe that our SDI images are the highest contrast astronomical images ever made from ground or space for methane rich companions <1'' from their star. For the best 20 of our survey stars, we attained 50% 5 sigma completeness for 6-10 M_Jup planets at semi-major axes of 20-40 AU. Thus, our completeness levels are sufficient to significantly test theoretical planet distributions. From our survey null result, we can rule out (at the 98% confidence/2.0sigma level) a model planet population using a planet distribution where N(a) $\propto$ constant out to a distance of 45 AU (further model assumptions discussed within).

We present an extension of the spectral deconvolution method (Sparks & Ford 2002) to achieve very high contrast at small inner working radii. We apply the method to the specific case of ground based adaptive optics fed integral field spectroscopy (without a coronagraph). Utilising the wavelength dependence of the Airy and speckle patterns, we make an accurate estimate of the PSF that can be scaled and subtracted from the data cube. The residual noise in the resulting spectra is very close to the photon noise from the starlight halo. We utilise the technique to extract a very high SNR H & K band spectrum of AB Dor C, the low mass companion to AB Dor A. By effectively eliminating all contamination from AB Dor A, the extracted spectrum retains both continuum and spectral features. The achieved 1 sigma contrast is 9 mag at 0.2", 11 mag at 0.5", in 20 mins exposure time, at an effective spectral bandwidth of 5.5 nm, proving that the method is applicable even in low Strehl regimes. The spectral deconvolution method clearly demonstrates the efficacy of image slicer based IFUs in achieving very high contrast imaging spectroscopy at small angular separations, validating their use as high contrast spectrographs/imagers for extreme adaptive optics systems.

We present new photometric and spectroscopic measurements for the unique, young, low-mass evolutionary track calibrator AB Dor C. While the new Ks photometry is similar to that previously published in Close et al. (2005) the spectral type is found to be earlier. Based on new H & K IFS spectra of AB Dor C (Thatte et al. 2007; paper 1) we adopt a spectral type of M5.5+/-1.0 for AB Dor C. This is considerably earlier than the M8+/-1 estimated in Close et al. (2005) and Nielsen et al. (2005) yet is consistent with the M6+/-1 independently derived by Luhman & Potter (2005). However, the spectrum presented in paper 1 and analyzed here is a significant improvement over any previous spectrum of AB Dor C. We also present new astrometry for the system which further supports a 0.090+/-0.005 Msun mass for the system. Once armed with an accurate spectrum and Ks flux we find L=0.0021+/-0.0005 Lsun and Teff=2925{+170}{-145}K for AB Dor C. These values are consistent with a ~75 Myr 0.090+/-0.005 Msun object like AB Dor C according to the DUSTY evolutionary tracks (Chabrier et al. 2000). Hence masses can be estimated from the HR diagram with the DUSTY tracks for young low-mass objects like AB Dor C. However, we cautiously note that underestimates of the mass from the tracks can occur if one lacks a proper (continuum preserved) spectra or is relying on NIR fluxes alone.

The active K2V star $\epsilon$ Eri hosts the most nearby known extrasolar planet. With an angular separation of about 1" on average, and an age of a few to several hundred Myrs, $\epsilon$ Eri b is one of the prime candidates for becoming the first definitive extrasolar planet imaged directly. We present a multi-epoch deep differential imaging survey performed with NACO-SDI at the VLT with the aim of finding the planet. The results are combined with recent astrometry in an attempt to further constrain the detection limits. No convincing candidate is found among the many coherent structures that constitute the residual speckle noise, which is the dominant noise at small angular scales. We present our detection limits, compare them with the estimated brightness of $\epsilon$ Eri b, and analyze how the limits can be improved further. It is found that integration time remains a very important parameter for achieving good results, even in the speckle-dominated regimes. The results yield new, improved upper 3$\sigma$ limits on the absolute H-band (1.6 $\mu$m) brightness of the 1.55 $M_{\rm jup}$ companion of 19.1 to 19.5 mag, depending on the specific age of the system.

We present resolved NACO photometry of the close binary AB Dor B in H- and Ks-band. AB Dor B is itself known to be a wide binary companion to AB Dor A, which in turn has a very low-mass close companion named AB Dor C. These four known components make up the young and dynamically interesting system AB Dor, which will likely become a benchmark system for calibrating theoretical pre-main sequence evolutionary mass tracks for low-mass stars. However, for this purpose the actual age has to be known, and this subject has been a matter of discussion in the recent scientific literature. We compare our resolved photometry of AB Dor Ba and Bb with theoretical and empirical isochrones in order to constrain the age of the system. This leads to an age estimate of about 50 to 100 Myr. We discuss the implications of such an age range for the case of AB Dor C, and compare with other results in the literature.

We report laboratory results of a coronagraphic test bench to assess the intensity reduction differences between a "Gaussian" tapered focal plane coronagraphic mask and a classical hard-edged "Top Hat" function mask at Extreme Adaptive Optics (ExAO) Strehl ratios of ~94%. However, unlike a traditional coronagraph design, we insert a reflective focal plane mask at 45 degree to the optical axis. We also used an intermediate secondary mask ("Mask_2") before a final image in order to block additional mask-edge diffracted light. The test bench simulates the 8.1m Gemini North telescope. It includes one spider vane, different mask radii (r= 1.9, 3.7, 7.4 lambda/D) and two types of reflective focal plane masks (hard-edged "Top Hat" and "Gaussian" tapered profiles). In order to investigate the relative performance of these competing coronagraphic designs with regard to extra-solar planet detection sensitivity, we utilize the simulation of realistic extra-solar planet populations (Nielson et al. 2006). With an appropriate translation of our laboratory results to expected telescope performance, a "Gaussian" tapered mask radius of 3.7 lambda/D with an additional mask ("Mask_2") performs best (highest planet detection sensitivity). For a full survey with this optimal design, the simulation predicts ~30% more planets detected compared to a similar sized "Top Hat" function mask with "Mask_2." Using the best design, the point contrast ratio between the stellar PSF peak and the coronagraphic PSF at 10 lambda/D (0.4" in H band if D = 8.1m) is ~10 times higher than a classical Lyot "Top Hat" coronagraph. Hence, we find a Gaussian apodized mask with an additional blocking mask is a superior (~10x higher contrast) than use of a classical Lyot coronagraph for ExAO-like Strehls.

We discuss the preliminary results of a survey of young (<300 Myr), close (<50 pc) stars with the Simultaneous Differential Extrasolar Planet Imager (SDI) implemented at the VLT and the MMT. SDI uses a double Wollaston prism and a quad filter to take 4 identical images simultaneously at 3 wavelengths surrounding the 1.62 um methane bandhead found in the spectrum of cool brown dwarfs and gas giants. By performing a difference of images in these filters, speckle noise from the primary can be significantly attenuated, resulting in photon noise limited data. In our survey data, we achieved H band contrasts >25000 (5 sigma Delta F1(1.575 um) > 10 mag, Delta H > 11.5 mag for a T6 spectral type) at a separation of 0.5" from the primary star. With this degree of attenuation, we should be able to image (5 sigma detection) a 2-4 Jupiter mass planet at 5 AU around a 30 Myr star at 10 pc. We are currently completing our survey of young, nearby stars. We have obtained complete datasets for 35 stars in the southern sky (VLT) and 7 stars in the northern sky (MMT). We believe that our SDI images are the highest contrast astronomical images ever made from ground or space for methane rich companions.

We discuss the instrumental and data reduction techniques used to suppress speckle noise with the Simultaneous Differential Imager (SDI) implemented at the VLT and the MMT. SDI uses a double Wollaston prism and a quad filter to take 4 identical images simultaneously at 3 wavelengths surrounding the 1.62 um methane bandhead found in the spectrum of cool brown dwarfs and gas giants. By performing a difference of images in these filters, speckle noise from the primary can be significantly attenuated, resulting in photon noise limited data past 0.5''. Non-trivial data reduction tools are necessary to pipeline the simultaneous differential imaging. Here we discuss a custom algorithm implemented in IDL to perform this reduction. The script performs basic data reduction tasks but also precisely aligns images taken in each of the filters using a custom shift and subtract routine. In our survey of nearby young stars at the VLT and MMT (see Biller et al., this conference), we achieved H band contrasts >25000 (5 sigma Delta F1(1.575 um) > 10.0 mag, Delta H > 11.5 mag for a T6 spectral type object) at a separation of 0.5" from the primary star. We believe that our SDI images are among the highest contrast astronomical images ever made from ground or space for methane rich companions.

We expand upon the results of Close et al. 2005 regarding the young, low-mass object AB Dor C and its role as a calibration point for theoretical tracks. We present an improved spectral reduction and a new orbital solution with two additional epochs. Our improved analysis confirms our spectral type of M8 (+/- 1) and mass of 0.090+/-0.003 solar masses for AB Dor C. Comparing the results for AB Dor C with other young, low-mass objects with dynamical masses we find a general trend where current evolutionary models tend to over-predict the temperature (or under-predict the mass) for low mass stars and brown dwarfs. Given our precision, there is a ~99% chance that the mass of AB Dor C is underestimated by the DUSTY tracks in the HR diagram.