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.

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

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

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

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

The continual push to directly image exoplanets at lower masses and closer separations orbiting around bright stars remains limited by both quasi-static and residual adaptive optics (AO) aberration. In previous papers we have proposed a modification of the self-coherent camera (SCC) design to address both of these limitations, called the Fast Atmospheric SCC Technique (FAST). In this paper we introduce an additional modification to the FAST focal plane mask design, including the existing Tip/tilt and Gaussian components and adding a charge four Vortex (TGV) component. In addition to boosting SCC fringe signal-to-noise ratio (S/N) as in our previous design, we show that the FAST TGV mask is also optimized to reach high contrast at separations closer to the star. In this paper we use numerical simulations to consider the performance improvement on correcting quasi-static aberration using this new mask compared to the previously proposed Tip/tilt+Gaussian mask. Using active deformable mirror control to generate a calibrated half dark hole improves contrast by a factor of about 200 at 2 - 5 $\lambda/D$ and up to a factor of 10 at 5 - 20 $\lambda/D$. The new methodology presented in this paper, now simultaneously considering both contrast and fringe S/N, opens the door to a new ideology of coronagraph design, where the coronagraph is now considered in duality as both a diffraction attenuator and a wavefront sensor.

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

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

Direct detection and detailed characterization of exoplanets using extreme adaptive optics (ExAO) is a key science goal of future extremely large telescopes (ELTs). However, wavefront errors will limit the sensitivity of this endeavor. Limitations for ground-based telescopes arise from both quasi-static and residual AO-corrected atmospheric wavefront errors, the latter of which generates short-lived aberrations that will average into a halo over a long exposure. We have developed and tested the framework for a solution to both of these problems using the self-coherent camera (SCC), to be applied to ground-based telescopes, called the Fast Atmospheric SCC Technique (FAST). In this paper we present updates of new and ongoing work for FAST, both in numerical simulation and in the laboratory. We first present numerical simulations that illustrate the scientific potential of FAST, including, with current 10-m telescopes, the direct detection of exoplanets reflected light and exo-Jupiters in thermal emission and, with future ELTs, the detection of habitable exoplanets. In the laboratory, we present the first characterizations of our proposed, and now fabricated, coronagraphic masks.

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

We describe the current performance of the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument on the Subaru telescope on Maunakea, Hawaii and present early science results for SCExAO coupled with the CHARIS integral field spectrograph. SCExAO now delivers H band Strehl ratios up to $\sim$ 0.9 or better, extreme AO corrections for optically faint stars, and planet-to-star contrasts rivaling that of GPI and SPHERE. CHARIS yield high signal-to-noise detections and 1.1--2.4 $\mu m$ spectra of benchmark directly-imaged companions like HR 8799 cde and kappa And b that clarify their atmospheric properties. We also show how recently published as well as unpublished observations of LkCa 15 lead to a re-evaluation of its claimed protoplanets. Finally, we briefly describe plans for a SCExAO-focused direct imaging campaign to directly image and characterize young exoplanets, planet-forming disks, and (later) mature planets in reflected light.

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.

We present an analysis of instrument performance using new observations taken with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS) instrument and the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system. In a correlation analysis of our datasets (which use the broadband mode covering J through K band in a single spectrum), we find that chromaticity in the SCExAO/CHARIS system is generally worse than temporal stability. We also develop a point spread function (PSF) subtraction pipeline optimized for the CHARIS broadband mode, including a forward modelling-based exoplanet algorithmic throughput correction scheme. We then present contrast curves using this newly developed pipeline. An analogous subtraction of the same datasets using only the H band slices yields the same final contrasts as the full JHK sequences; this result is consistent with our chromaticity analysis, illustrating that PSF subtraction using spectral differential imaging (SDI) in this broadband mode is generally not more effective than SDI in the individual J, H, or K bands. In the future, the data processing framework and analysis developed in this paper will be important to consider for additional SCExAO/CHARIS broadband observations and other ExAO instruments which plan to implement a similar integral field spectrograph broadband mode.

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.

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.

Exoplanet detection and characterization through extreme adaptive optics (ExAO) is a key science goal of future extremely large telescopes. This achievement, however, will be limited in sensitivity by both quasi-static wavefront errors and residual AO-corrected atmospheric wavefront errors. A solution to both of these problems is to use the science camera of an ExAO system as a wavefront sensor to perform a fast measurement and correction method to remove these aberrations as soon as they are detected. We have developed the framework for one such method, using the self-coherent camera (SCC), to be applied to ground-based telescopes, called Fast Atmospheric SCC Technique (FAST; Gerard et al., submitted). Our FAST solution requires an optimally designed coronagraph (the SCC FPM) and post-processing algorithm and is in principle able to reach a "raw" contrast of a few times the photon noise limit, continually improving with integration time. In this paper, we present new ongoing work in exploring the manufacturing limitations of the SCC FPM as well as a general framework to implement and optimize a FAST deformable mirror control loop.

Direct imaging and spectral characterization of exoplanets using extreme adaptive optics (ExAO) is a key science goal of future extremely large telescopes and space observatories. However, quasi-static wavefront errors will limit the sensitivity of this endeavor. Additional limitations for ground-based telescopes arise from residual AO-corrected atmospheric wavefront errors, generating millisecond-lifetime speckles that average into a halo over a long exposure. A solution to both of these problems is to use the science camera of an ExAO system as a wavefront sensor to perform a fast measurement and correction method to minimize these aberrations as soon as they are detected. We develop the framework for one such method based on the self-coherent camera (SCC) to be applied to ground-based telescopes, called Fast Atmospheric SCC Technique (FAST). We show that with the use of a specially designed coronagraph and coherent differential imaging algorithm, recording images every few milliseconds allows for a subtraction of atmospheric and static speckles while maintaining a close to unity algorithmic exoplanet throughput. Detailed simulations reach a contrast close to the photon noise limit after 30 seconds for a 1 % bandpass in H band on both 0$^\text{th}$ and 5$^\text{th}$ magnitude stars. For the 5th magnitude case, this is about 110 times better in raw contrast than what is currently achieved from ExAO instruments if we extrapolate for an hour of observing time, illustrating that sensitivity improvement from this method could play an essential role in the future detection and characterization of lower mass exoplanets.

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.

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.

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.

Most current high contrast imaging point spread function (PSF) subtraction algorithms use some form of a least-squares noise minimization to find exoplanets that are, before post-processing, often hidden below the instrumental speckle noise. In the current standard PSF subtraction algorithms, a set of reference images is derived from the target image sequence to subtract each target image, using Angular and/or Simultaneous Spectral Differential Imaging (ADI, SSDI, respectively). However, to avoid excessive exoplanet self-subtraction, ADI and SSDI (in the absence of a strong spectral feature) severely limit the available number of reference images at small separations. This limits the performance of the least-squares algorithm, resulting in lower sensitivity to exoplanets at small angular separations. Possible solutions are to use additional reference images by acquiring longer sequences, use SSDI if the exoplanet is expected to show strong spectral features, or use images acquired on other targets. The latter option, known as Reference Star Differential Imaging (RSDI), which relies on the use of reference images that are highly correlated to the target image, has been ineffective in previous ground-based high contrast imaging surveys. We present the results of work to optimize PSF subtraction with the GPIES reference library using a least-squares algorithm designed to minimize speckle noise and maximize planet throughput, thus maximizing the planet signal to noise ratio (SNR). Using December 2014 51 Eri GPI data in the inner 100 mas to 300 mas annulus, we find no apparent improvement in SNR when using RSDI and/or our optimization scheme. This result, while still being investigated, seems to show that current algorithms on ADI+SSDI data sets are optimized, and that limited gains can be achieved by using a PSF archive.

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.

One of the key science goals for extremely large telescopes (ELTs) is the detailed characterization of already known directly imaged exoplanets. The typical adaptive optics (AO) Nyquist control region for ELTs is ~0.4 arcseconds, placing many already known directly imaged planets outside the DM control region and not allowing any standard wavefront control scheme to remove speckles that would allow higher SNR images/spectra to be acquired. This can be fixed with super-Nyquist wavefront control (SNWFC), using a sine wave phase plate to allow for wavefront control outside the central DM Nyquist region. We demonstrate that SNWFC is feasible through a simple, deterministic, non-coronagraphic, super-Nyquist speckle nulling technique in the adaptive optics laboratory at the National Research Council of Canada. We also present results in simulation of how SNWFC using the self coherent camera (SCC) can be used for high contrast imaging. This technique could be implemented on future high contrast imaging instruments to improve contrast outside the standard central dark hole for higher SNR characterization of exoplanets.

We report the results of a Green Bank Telescope survey for water masers, ammonia (1,1) and (2,2), and the H66-alpha recombination line toward 506 luminous compact 24 micron-emitting regions in the Andromeda Galaxy (M31). We include the 206 sources observed in the Darling (2011) water maser survey for completeness. The survey was sensitive enough to detect any maser useful for ~10 microarcsecond/yr astrometry. No new water masers, ammonia lines, or H66-alpha recombination lines were detected individually or in spectral stacks reaching rms noise levels of ~3 mJy and ~0.2 mJy, respectively, in 3.1-3.3 km/s channels. The lack of detections in individual spectra and in the spectral stacks is consistent with Galactic extrapolations. Contrary to previous assertions, there do not seem to be additional bright water masers to be found in M31. The strong variability of water masers may enable new maser detections in the future, but variability may also limit the astrometric utility of known (or future) masers since flaring masers must also fade.

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 a new algorithm for space telescope high contrast imaging of close-to-face-on planetary disks called Optimized Spatially Filtered (OSFi) normalization. This algorithm is used on HR 8799 Hubble Space Telescope (HST) coronagraphic archival data, showing an over-luminosity after reference star point spread function (PSF) subtraction that may be from the inner disk and/or planetesimal belt components of this system. The PSF-subtracted radial profiles in two separate epochs from 2011 and 2012 are consistent with one another, and self-subtraction shows no residual in both epochs. We explore a number of possible false-positive scenarios that could explain this residual flux, including telescope breathing, spectral differences between HR 8799 and the reference star, imaging of the known warm inner disk component, OSFi algorithm throughput and consistency with the standard spider normalization HST PSF subtraction technique, and coronagraph misalignment from pointing accuracy. In comparison to another similar STIS dataset, we find that the over-luminosity is likely a result of telescope breathing and spectral difference between HR 8799 and the reference star. Thus, assuming a non-detection, we derive upper limits on the HR 8799 dust belt mass in small grains. In this scenario, we find that the flux of these micron-sized dust grains leaving the system due to radiation pressure is small enough to be consistent with measurements of other debris disk halos.

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

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

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