More than a dozen young stars host spiral arms in their surrounding protoplanetary disks. The excitation mechanisms of such arms are under debate. The two leading hypotheses -- companion-disk interaction and gravitational instability (GI) -- predict distinct motion for spirals. By imaging the MWC 758 spiral arm system at two epochs spanning ${\sim}5$ yr using the SPHERE instrument on the Very Large Telescope (VLT), we test the two hypotheses for the first time. We find that the pattern speeds of the spirals are not consistent with the GI origin. Our measurements further evince the existence of a faint "missing planet" driving the disk arms. The average spiral pattern speed is $0.\!^\circ22\pm0.\!^\circ03$ yr$^{-1}$, pointing to a driver at $172_{-14}^{+18}$ au around a $1.9$ $M_\odot$ central star if it is on a circular orbit. In addition, we witness time varying shadowing effects on a global scale that are likely originated from an inner disk.

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.

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.

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

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

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

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

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

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

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

Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and co-rotate with giant planets, which has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10-yr timescale, we are able to provide dynamical constraints on the locations of their perturbers. We present reprocessed Hubble Space Telescope (HST)/NICMOS F110W observations of the target in 2005, and the new Keck/NIRC2 $L'$-band observations in 2017. MWC~758's two well-known spiral arms are revealed in the NICMOS archive at the earliest observational epoch. With additional Very Large Telescope (VLT)/SPHERE data, our joint analysis leads to a pattern speed of ${0.6^\circ}^{+3.3^\circ}_{-0.6^\circ}\, \mathrm{yr}^{-1}$ at $3\sigma$ for the two major spiral arms. If the two arms are induced by a perturber on a near-circular orbit, its best fit orbit is at $89$ au ($0.59"$), with a $3\sigma$ lower limit of 30 au ($0.20"$). This finding is consistent with the simulation prediction of the location of an arm-driving planet for the two major arms in the system.

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.

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.

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

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

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.

We describe a "super-Earth"-size ($2.30\pm0.15R_{\oplus}$) planet transiting an early K-type dwarf star in the Campaign 4 field observed by the K2 mission. The host star, EPIC 210363145, was identified as a member of the approximately 120-Myr-old Pleiades cluster based on its kinematics and photometric distance. It is rotationally variable and exhibits near-ultraviolet emission consistent with a Pleiades age, but its rotational period is ~20 d and its spectrum contains no H$\alpha$ emission nor the Li I absorption expected of Pleiades K dwarfs. Instead, the star is probably an interloper that is unaffiliated with the cluster, but younger (< 1 Gyr) than the typical field dwarf. We ruled out a false positive transit signal produced by confusion with a background eclipsing binary by adaptive optics imaging and a statistical calculation. Doppler radial velocity measurements limit the companion mass to <2 times that of Jupiter. Screening of the lightcurves of 1014 potential Pleiades candidate stars uncovered no additional planets. An injection-and-recovery experiment using the K2 Pleiades lightcurves with simulated planets, assuming a planet population like that in the Kepler prime field, predicts only 0.8-1.8 detections (vs. ~20 in an equivalent Kepler sample). The absence of Pleiades planet detections can be attributed to the much shorter monitoring time of K2 (80 days vs. 4 years), increased measurement noise due to spacecraft motion, and the intrinsic noisiness of the stars.

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

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

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

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

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

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

We present a new technique for forward-modeling self-subtraction of spatially extended emission in observations processed with angular differential imaging (ADI) algorithms. High-contrast direct imaging of circumstellar disks is limited by quasi-static speckle noise and ADI is commonly used to suppress those speckles. However, the application of ADI can result in self-subtraction of the disk signal due to the disk's finite spatial extent. This signal attenuation varies with radial separation and biases measurements of the disk's surface brightness, thereby compromising inferences regarding the physical processes responsible for the dust distribution. To compensate for this attenuation, we forward-model the disk structure and compute the form of the self-subtraction function at each separation. As a proof of concept, we apply our method to 1.6 and 2.2 micron Keck AO NIRC2 scattered-light observations of the HD 32297 debris disk reduced using a variant of the "locally optimized combination of images" (LOCI) algorithm. We are able to recover disk surface brightness that was otherwise lost to self-subtraction and produce simplified models of the brightness distribution as it appears with and without self-subtraction. From the latter models, we extract radial profiles for the disk's brightness, width, midplane position, and color that are unbiased by self-subtraction. Our analysis of these measurements indicates a break in the brightness profile power law at r~110 AU and a disk width that increases with separation from the star. We also verify disk curvature that displaces the midplane by up to 30 AU towards the northwest relative to a straight fiducial midplane.