We report the results of a ∼4 yr direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager (GPI) Exoplanet Survey. We targeted nearby (≲150 pc), young (≲500 Myr) stars with high infrared (IR) excesses (L _IR/L _⋆ > 10⁻⁵), including 38 with previously resolved disks. Observations were made using the GPI high-contrast integral field spectrograph in H-band (1.6 μm) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and 3 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 except 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 nondetections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low-inclination disks (≲70°). Based on postsurvey statistics, we improved upon our presurvey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing H-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

We present H-band observations of β Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~03 (6 AU) and ~17 (33 AU), while simultaneously detecting β 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 () with a position angle (PA) (slightly offset from the main outer disk, ), that extends from an inner disk radius of to well outside GPI's field of view. In addition, we present an updated orbit for β Pic b based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semimajor axis of , with an eccentricity The PA of the ascending node is offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of β Pic to Dynamical sculpting by β 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 β Pic b; (2) the mutual inclination of the inner disk and β Pic b is when it is expected to be closer to zero; and (3) the aspect ratio of the disk () is larger than expected from interactions with β Pic b or self-stirring by the disk's parent bodies.

Using the Gemini Planet Imager located at Gemini South, we measured the near-infrared (1.0–2.4 μm) spectrum of the planetary companion to the nearby, young star β 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 ± 0.2 , an effective temperature of 1724 ± 15 K, a radius of 1.46 ± 0.01 , and a surface gravity of [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 –4.0 [dex] depending upon the model. These values agree well with other publications and with “hot-start” predictions from planetary evolution models. Further, we find that the spectrum of β Pic b best matches a low surface gravity L2 ± 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 have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.