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

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

Radial Velocity Measurements of HR 8799 b and c with Medium Resolution Spectroscopy

High-contrast medium resolution spectroscopy has been used to detect molecules such as water and carbon monoxide in the atmospheres of gas giant exoplanets. In this work, we show how it can be used to derive radial velocity (RV) measurements of directly imaged exoplanets. Improving upon the traditional cross-correlation technique, we develop a new likelihood based on joint forward modeling of the planetary signal and the starlight background (i.e., speckles). After marginalizing over the starlight model, we infer the barycentric RV of HR 8799 b and c in 2010 yielding −9.2 ± 0.5 km s⁻¹ and −11.6 ± 0.5 km s⁻¹, respectively. These RV measurements help to constrain the 3D orientation of the orbit of the planet by resolving the degeneracy in the longitude of an ascending node. Assuming coplanar orbits for HR 8799 b and c, but not including d and e, we estimate Ω = 89°⁺²⁷₋₁₇ and i = 20°.8^(4.5)_(-3.7)

Gemini Planet Imager Observational Calibrations VI: Photometric and Spectroscopic Calibration for the Integral Field Spectrograph

The Gemini Planet Imager (GPI) is a new facility instrument for the Gemini Observatory designed to provide direct detection and characterization of planets and debris disks around stars in the solar neighborhood. In addition to its extreme adaptive optics and corona graphic systems which give access to high angular resolution and high-contrast imaging capabilities, GPI contains an integral field spectrograph providing low resolution spectroscopy across five bands between 0.95 and 2.5 $\mu$m. This paper describes the sequence of processing steps required for the spectro-photometric calibration of GPI science data, and the necessary calibration files. Based on calibration observations of the white dwarf HD 8049B we estimate that the systematic error in spectra extracted from GPI observations is less than 5%. The flux ratio of the occulted star and fiducial satellite spots within coronagraphic GPI observations, required to estimate the magnitude difference between a target and any resolved companions, was measured in the $H$-band to be $\Delta m = 9.23\pm0.06$ in laboratory measurements and $\Delta m = 9.39\pm 0.11$ using on-sky observations. Laboratory measurements for the $Y$, $J$, $K1$ and $K2$ filters are also presented. The total throughput of GPI, Gemini South and the atmosphere of the Earth was also measured in each photometric passband, with a typical throughput in $H$-band of 18% in the non-coronagraphic mode, with some variation observed over the six-month period for which observations were available. We also report ongoing development and improvement of the data cube extraction algorithm.Comment: 15 pages, 6 figures. Proceedings of the SPIE, 9147-30

Deep exploration of the planets HR 8799 b, c, and d with moderate-resolution spectroscopy

Funding: J.-B.R. acknowledges support from the David and Ellen Lee Prize Postdoctoral Fellowship. The research was supported by grants from NSF, including AST-1411868 (J.-B.R., B.M.) and 1614492 (T.S.B.). Material presented in this work is supported by the National Aeronautics and Space Administration under Grants/Contracts/Agreements No. NNX17AB63G (Q.M.K., T.S.B., and K.K.W.) issued through the Astrophysics Division of the Science Mission Directorate and NNX15AD95G (J.-B.R., R.J.D.R.).The four directly imaged planets orbiting the star HR 8799 are an ideal laboratory to probe atmospheric physics and formation models. We present more than a decade's worth of Keck/OSIRIS observations of these planets, which represent the most detailed look at their atmospheres to date by its resolution and signal-to-noise ratio. We present the first direct detection of HR 8799 d, the second-closest known planet to the star, at moderate spectral resolution with Keck/OSIRIS (K band; R ≈ 4000). Additionally, we uniformly analyze new and archival OSIRIS data (H and K band) of HR 8799 b, c, and d. First, we show detections of water (H2O) and carbon monoxide (CO) in the three planets and discuss the ambiguous case of methane (CH4) in the atmosphere of HR 8799 b. Then, we report radial-velocity (RV) measurements for each of the three planets. The RV measurement of HR 8799 d is consistent with predictions made assuming coplanarity and orbital stability of the HR 8799 planetary system. Finally, we perform a uniform atmospheric analysis on the OSIRIS data, published photometric points, and low-resolution spectra. We do not infer any significant deviation from the stellar value of the carbon-to-oxygen ratio (C/O) of the three planets, which therefore does not yet yield definitive information about the location or method of formation. However, constraining the C/O for all the HR 8799 planets is a milestone for any multiplanet system, and particularly important for large, widely separated gas giants with uncertain formation processes.Publisher PDFPeer reviewe

A Bayesian framework for exoplanet direct detection and non-detection

Funding: This research was supported by grants from NSF, including AST-1411868 (J.-B.R., B.M.) and AST-1518332 (R.J.D.R.). Support was provided by grants from NASA, including NNX14AJ80G (B.M., J.-B.R.), NNX15AD95G (R.J.D.R.) and NNX15AC89G (R.J.D.R.).Rigorously quantifying the information in high-contrast imaging data is important for informing follow-up strategies to confirm the substellar nature of a point source, constraining theoretical models of planet-disk interactions, and deriving planet occurrence rates. However, within the exoplanet direct imaging community, non-detections have almost exclusively been defined using a frequentist detection threshold (i.e., contrast curve) and associated completeness. This can lead to conceptual inconsistencies when included in a Bayesian framework. A Bayesian upper limit is such that the true value of a parameter lies below this limit with a certain probability. The associated probability is the integral of the posterior distribution with the upper limit as the upper bound. In summary, a frequentist upper limit is a statement about the detectability of planets while a Bayesian upper limit is a statement about the probability of a parameter to lie in an interval given the data. The latter is therefore better suited for rejecting hypotheses or theoretical models based on their predictions. In this work we emphasize that Bayesian statistics and upper limits are more easily interpreted and typically more constraining than the frequentist approach. We illustrate the use of Bayesian analysis in two different cases: (1) with a known planet location where we also propose to use model comparison to constrain the astrophysical nature of the point source and (2) gap-carving planets in TW Hya. To finish, we also mention the problem of combining radial velocity and direct imaging observations.Publisher PDFPeer reviewe

Exoplanet Imaging Data Challenge, phase II: Characterization of exoplanet signals in high-contrast images

Today, there exists a wide variety of algorithms dedicated to high-contrast imaging, especially for the detection and characterisation of exoplanet signals. These algorithms are tailored to address the very high contrast between the exoplanet signal(s), which can be more than two orders of magnitude fainter than the bright starlight residuals in coronagraphic images. The starlight residuals are inhomogeneously distributed and follow various timescales that depend on the observing conditions and on the target star brightness. Disentangling the exoplanet signals within the starlight residuals is therefore challenging, and new post-processing algorithms are striving to achieve more accurate astrophysical results. The Exoplanet Imaging Data Challenge is a community-wide effort to develop, compare and evaluate algorithms using a set of benchmark high-contrast imaging datasets. After a first phase ran in 2020 and focused on the detection capabilities of existing algorithms, the focus of this ongoing second phase is to compare the characterisation capabilities of state-of-the-art techniques. The characterisation of planetary companions is two-fold: the astrometry (estimated position with respect to the host star) and spectrophotometry (estimated contrast with respect to the host star, as a function of wavelength). The goal of this second phase is to offer a platform for the community to benchmark techniques in a fair, homogeneous and robust way, and to foster collaborations.Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2022, Adaptive Optics Systems VIII, Paper 12185-

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

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.Comment: ApJ accepte

GPI spectra of HR 8799 c, d, and e from 1.5 to 2.4μ\mum with KLIP Forward Modeling

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.Comment: Accepted to AJ, 25 pages, 16 Figure