Gemini Planet Imager Exoplanet Survey Research Articles

A Uniform Analysis of Debris Disks with the Gemini Planet Imager. I. An Empirical Search for Perturbations from Planetary Companions in Polarized Light Images

The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey imaged 24 debris disks in polarized H-band light, while other programs observed half of these disks in polarized J and/or K1 bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, particularly those due to planet–disk interactions. The multiwavelength surface brightness, disk color, and geometry permit the identification of any asymmetries such as warps or disk offsets from the central star. We find that 19 of the disks in this sample exhibit asymmetries in surface brightness, disk color, disk geometry, or a combination of the three, suggesting that for this sample, perturbations, as seen in scattered light, are common. The relationship between these perturbations and potential planets in the system is discussed. We also explore correlations among stellar temperatures, ages, disk properties, and observed perturbations. We find significant trends between the vertical aspect ratio and the stellar temperature, disk radial extent, and the dust grain size distribution power law, q. We also confirm a trend between the disk color and stellar effective temperature, where the disk becomes increasingly red/neutral with increasing temperature. Such results have important implications for the evolution of debris disk systems around stars of various spectral types.

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

Abstract 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−5), 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.

Effects of mirror seeing on high-contrast adaptive optics instruments

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.

Giant Planets and Brown Dwarfs Form in Different Ways

Once thought to be part of the same population, planets larger than Jupiter and “failed stars” likely grow via different mechanisms, the Gemini Planet Imager Exoplanet Survey has shown.

The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics from 10 to 100 au

Abstract We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey. 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, semimajor axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M * > 1.5 M ⊙ more likely to host planets with masses between 2 and 13M Jup and semimajor axes of 3–100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semimajor axis (a) for planet populations around high-mass stars (M * > 1.5 M ⊙) of the form , finding α = −2.4 ± 0.8 and β = −2.0 ± 0.5, and an integrated occurrence rate of % between 5–13M Jup and 10–100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with % of stars hosting a brown dwarf companion between 13–80M Jup and 10–100 au. Brown dwarfs also appear to be distributed differently in mass and semimajor axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semimajor axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the radial velocity method, our results are consistent with a peak in occurrence of giant planets between ∼1 and 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.

Automated data processing architecture for the Gemini Planet Imager Exoplanet Survey

The Gemini Planet Imager Exoplanet Survey (GPIES) is a multi-year direct imaging survey of 600 stars to discover and characterize young Jovian exoplanets and their environments. We have developed an automated data architecture to process and index all data related to the survey uniformly. An automated and flexible data processing framework, which we term the Data Cruncher, combines multiple data reduction pipelines together to process all spectroscopic, polarimetric, and calibration data taken with GPIES. With no human intervention, fully reduced and calibrated data products are available less than an hour after the data are taken to expedite follow-up on potential objects of interest. The Data Cruncher can run on a supercomputer to reprocess all GPIES data in a single day as improvements are made to our data reduction pipelines. A backend MySQL database indexes all files, which are synced to the cloud, and a front-end web server allows for easy browsing of all files associated with GPIES. To help observers, quicklook displays show reduced data as they are processed in real-time, and chatbots on Slack post observing information as well as reduced data products. Together, the GPIES automated data processing architecture reduces our workload, provides real-time data reduction, optimizes our observing strategy, and maintains a homogeneously reduced dataset to study planet occurrence and instrument performance.

Characterizing 51 Eri b from 1 to 5 μm: A Partly Cloudy Exoplanet

Abstract We present spectrophotometry spanning 1–5 μm of 51 Eridani b, a 2–10 planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new K1 (1.90–2.19 μm) and K2 (2.10–2.40 μm) spectra taken with the Gemini Planet Imager as well as an updated L P (3.76 μm) and new M S (4.67 μm) photometry from the NIRC2 Narrow camera. The new data were combined with J (1.13–1.35 μm) and H (1.50–1.80 μm) spectra from the discovery epoch with the goal of better characterizing the planet properties. The 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4 and 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. The model fits suggest that 51 Eri b has an effective temperature ranging between 605 and 737 K, a solar metallicity, and a surface gravity of log(g) = 3.5–4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the spectral energy distribution (SED). From the model atmospheres, we infer a luminosity for the planet of −5.83 to −5.93 ( ), leaving 51 Eri b in the unique position of being one of the only directly imaged planets consistent with having formed via a cold-start scenario. Comparisons of the planet SED against warm-start models indicate that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 .

Improving and Assessing Planet Sensitivity of the GPI Exoplanet Survey with a Forward Model Matched Filter

Abstract 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éve 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 (S/N) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal S/N loss. We also developed a complete pipeline for the automated detection of point-source candidates, the calculation of receiver operating characteristics (ROC), contrast curves based on false positives, and completeness contours. We process in a uniform manner more than 330 data sets from the Gemini Planet Imager Exoplanet Survey 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.

Integral Field Spectroscopy of the Low-mass Companion HD 984 B with the Gemini Planet Imager

Abstract We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager (GPI) as a part of the GPI Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12–1.3 μm) and H (1.50–1.80 μm) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 au (70-year) orbit with a 68% confidence interval between 14 and 28 au, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119° with a 68% confidence interval between 114° and 125°. To address the considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare them to known spectral types using field brown dwarf and low-mass star spectra and find a best-fit match of a field gravity M6.5 ± 1.5 spectral type with a corresponding temperature of K. Photometry of the companion yields a luminosity of / dex with DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34 ± 1 to 95 ± 4 M Jup. These results are consistent with previous measurements of the object.

Performance of the Gemini Planet Imager's adaptive optics system.

The Gemini Planet Imager's adaptive optics (AO) subsystem was designed specifically to facilitate high-contrast imaging. A definitive description of the system's algorithms and technologies as built is given. 564 AO telemetry measurements from the Gemini Planet Imager Exoplanet Survey campaign are analyzed. The modal gain optimizer tracks changes in atmospheric conditions. Science observations show that image quality can be improved with the use of both the spatially filtered wavefront sensor and linear-quadratic-Gaussian control of vibration. The error budget indicates that for all targets and atmospheric conditions AO bandwidth error is the largest term.

Mapping the Distributions of Exoplanet Populations with NICI and GPI

AbstractWhile more and more long-period giant planets are discovered by direct imaging, the distribution of planets at these separations (≳5 AU) has remained largely uncertain, especially compared to planets in the inner regions of solar systems probed by RV and transit techniques. The low frequency, the detection challenges, and heterogeneous samples make determining the mass and orbit distributions of directly imaged planets at the end of a survey difficult. By utilizing Monte Carlo methods that incorporate the age, distance, and spectral type of each target, we can use all stars in the survey, not just those with detected planets, to learn about the underlying population. We have produced upper limits and direct measurements of the frequency of these planets with the most recent generation of direct imaging surveys. The Gemini NICI Planet-Finding Campaign observed 220 young, nearby stars at a median H-band contrast of 14.5 magnitudes at 1”, representing the largest, deepest search for exoplanets by the completion of the survey. The Gemini Planet Imager Exoplanet Survey is in the process of surveying 600 stars, pushing these contrasts to a few tenths of an arcsecond from the star. With the advent of large surveys (many hundreds of stars) using advanced planet-imagers we gain the ability to move beyond measuring the frequency of wide-separation giant planets and to simultaneously determine the distribution as a function of planet mass, semi-major axis, and stellar mass, and so directly test models of planet formation and evolution.