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

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

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

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

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

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

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

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

The 32 Orionis group is a co-moving group of roughly 20 young (24 Myr) M3-B5 stars 100 pc from the Sun. Here we report the discovery of its first substellar member, WISE J052857.69+090104.2. This source was previously reported to be an M giant star based on its unusual near-infrared spectrum and lack of measurable proper motion. We re-analyze previous data and new moderate-resolution spectroscopy from Magellan/FIRE to demonstrate that this source is a young near-infrared L1 brown dwarf with very low surface gravity features. Spectral model fits indicate T$_{eff}$ = 1880$^{+150}_{-70}$ K and $\log{g}$ = 3.8$^{+0.2}_{-0.2}$ (cgs), consistent with a 15-22 Myr object with a mass near the deuterium-burning limit. Its sky position, estimated distance, kinematics (both proper motion and radial velocity), and spectral characteristics are all consistent with membership in 32 Orionis, and its temperature and age imply a mass (M = 14$^{+4}_{-3}$ M$_{Jup}$) that straddles the brown dwarf/planetary-mass object boundary. The source has a somewhat red $J-W2$ color compared to other L1 dwarfs, but this is likely a low-gravity-related temperature offset; we find no evidence of significant excess reddening from a disk or cool companion in the 3-5 $\mu$m waveband.

We present spectroscopic analysis of the low mass binary star system GJ 660.1AB, a pair of nearby M dwarfs for which we have obtained separated near-infrared spectra (0.9-2.5 $\mu$m) with the SpeX spectrograph. The spectrum of GJ 660.1B is distinctly peculiar, with a triangular-shaped 1.7 $\mu$m peak that initially suggests it to be a low surface gravity, young brown dwarf. However, we rule out this hypothesis and determine instead that this companion is a mild subdwarf (d/sdM7) based on the subsolar metallicity of the primary, [Fe/H] = $-$0.63$\pm$0.06. Comparison of the near-infrared spectrum of GJ 660.1B to two sets of spectral models yields conflicting results, with a common effective temperature T$_{eff}$ = 2550-2650 K, but alternately low surface gravity ($\log{g}$ = 4.4$\pm$0.5 cgs) and very low metallicity ([M/H] = $-$0.96$^{+0.19}_{-0.24}$), or high surface gravity ($\log{g}$ = 5.0-5.5) and slightly subsolar metallicity ([M/H] =$-$0.20$^{+0.13}_{-0.19}$). We conjecture that insufficient condensate opacity and excessive collision induced H$_2$ absorption in the models bias them toward low surface gravities and a metallicity inconsistent with the primary, and points toward improvements needed in the spectral modeling of metal-poor, very-low mass dwarfs. The peculiar spectral characteristics of GJ 660.1B emphasize that care is needed when interpreting surface gravity features in the spectra of ultracool dwarfs.