Engineering Papers⌕ Search

Engineering topics

Ruffio, Jean-Baptiste

Publications and source records attributed to Ruffio, Jean-Baptiste.

Octofitter: Fast, Flexible, and Accurate Orbit Modeling to Detect Exoplanets

As next-generation imaging instruments and interferometers search for planets closer to their stars, they must contend with increasing orbital motion and longer integration times. These compounding effects make it difficult to detect faint planets but also present an opportunity. Increased orbital motion makes it possible to move the search for planets into the orbital domain, where direct images can be freely combined with the radial velocity and proper motion anomaly, even without a confirmed detection in any single epoch. In this paper, we present a fast and differentiable multimethod orbit-modeling and planet detection code called Octofitter. This code is designed to be highly modular and allows users to easily adjust priors, change parameterizations, and specify arbitrary function relations between the parameters of one or more planets. Octofitter further supplies tools for examining model outputs including prior and posterior predictive checks and simulation-based calibration. We demonstrate the capabilities of Octofitter on real and simulated data from different instruments and methods, including HD 91312, simulated JWST/NIRISS aperture masking interferometry observations, radial velocity curves, and grids of images from the Gemini Planet Imager. We show that Octofitter can reliably recover faint planets in long sequences of images with arbitrary orbital motion. This publicly available tool will enable the broad application of multiepoch and multimethod exoplanet detection, which could improve how future targeted ground- and space-based surveys are performed. Finally, its rapid convergence makes it a useful addition to the existing ecosystem of tools for modeling the orbits of directly imaged planets.

79 ASTRONOMY AND ASTROPHYSICS↗

Detection and Bulk Properties of the HR 8799 Planets with High-resolution Spectroscopy

Using the Keck Planet Imager and Characterizer, we obtained high-resolution (R ~ 35,000) K-band spectra of the four planets orbiting HR 8799. We clearly detected H{sub 2}O and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and H{sub 2}O in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward-modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured vsin (i) values of 10.1{sub -2.7}{sup +2.8} km s{sup -1} for HR 8799 d and 15.0{sub -2.6}{sup +2.3} km s{sup -1} for HR 8799 e, and placed an upper limit of <14 km s{sup -1} of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anticorrelated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower-mass planets.

74 ATOMIC AND MOLECULAR PHYSICS↗

Early High-contrast Imaging Results with Keck/NIRC2-PWFS: The SR 21 Disk

High-contrast imaging of exoplanets and protoplanetary disks depends on wave front sensing and correction made by adaptive optics instruments. Classically, wave front sensing has been conducted at optical wavelengths, which made high-contrast imaging of red targets such as M-type stars or extincted T Tauri stars challenging. Keck/NIRC2 has combined near-infrared (NIR) detector technology with the pyramid wave front sensor (PWFS). With this new module we observed SR 21, a young star that is brighter at NIR wavelengths than at optical wavelengths. Compared with the archival data of SR 21 taken with the optical wave front sensing we achieved ~20% better Strehl ratio in similar natural seeing conditions. Further post-processing utilizing angular differential imaging and reference-star differential imaging confirmed the spiral feature reported by the Very Large Telescope/Spectro-Polarimetric High-contrast Exoplanet REsearch instrument polarimetric observation, which is the first detection of the SR 21 spiral in total intensity at L' band. We also compared the contrast limit of our result (10 -4 at 0.'4 and 2 × 10 -5 at 1.'0) with the archival data that were taken with optical wave front sensing and confirmed the improvement, particularly at ≤0.'5. Our observation demonstrates that the NIR PWFS improves AO performance and will provide more opportunities for red targets in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Moderate-resolution K-band Spectroscopy of Substellar Companion κ Andromedae b

We present moderate-resolution (R ∼ 4000) K-band spectra of the “super-Jupiter,” κ Andromedae b. The data were taken with the OSIRIS integral field spectrograph at Keck Observatory. The spectra reveal resolved molecular lines from H{sub 2}O and CO, and are compared to a custom PHOENIX atmosphere model grid appropriate for young planetary-mass objects. We fit the data using a Markov chain Monte Carlo forward-modeling method. Using a combination of our moderate-resolution spectrum and low-resolution, broadband data from the literature, we derive an effective temperature of T {sub eff} = 1950–2150 K, a surface gravity of logg=3.5--4.5, and a metallicity of [M/H] = −0.2–0.0. These values are consistent with previous estimates from atmospheric modeling and the currently favored young age of the system (<50 Myr). We derive a C/O ratio of 0.70{sub −0.24}{sup +0.09} for the source, broadly consistent with the solar C/O ratio. This, coupled with the slightly subsolar metallicity, implies a composition consistent with that of the host star, and is suggestive of formation by a rapid process. The subsolar metallicity of κ Andromedae b is also consistent with predictions of formation via gravitational instability. Further constraints on formation of the companion will require measurement of the C/O ratio of κ Andromedae A. We also measure the radial velocity of κ Andromedae b for the first time, with a value of −1.4 ± 0.9 km s{sup −1} relative to the host star. We find that the derived radial velocity is consistent with the estimated high eccentricity of κ Andromedae b.

79 ASTRONOMY AND ASTROPHYSICS↗

Revised astrometric calibration of the Gemini Planet Imager

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

79 ASTRONOMY AND ASTROPHYSICS↗

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet β Pictoris b from Combined Direct Imaging and Astrometry

Herein, we present new observations of the planet β Pictoris b from 2018 with the Gemini Planet Imager (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.03}^{+0.04}$ disfavors circular orbits. We consider orbit fits to several different imaging data sets, 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, β 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 (8.0 ± 2.6 M Jup) for the imaged planet β Pic b, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

Asymmetries in adaptive optics point spread functions

An explanation for the origin of asymmetry along the preferential axis of the point spread function (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 toward 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.

79 ASTRONOMY AND ASTROPHYSICS↗

An Exo–Kuiper Belt with an Extended Halo around HD 191089 in Scattered Light

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 the Kuiper Belt in radial extent (FWHM ~ 25 au, centered at ~46 au) and a halo extending to ~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°–150° scattering-angle range and find that the halo dust is more forward- and backward-scattering than the ring dust. We measure a surface density power-law index of -0.68 ± 0.04 for the halo, which indicates the slowdown 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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

Here, we introduce 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 13$M$ 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 ${d}^{2}N/({dm}\,{da})\propto {m}^{\alpha }\,{a}^{\beta }$, finding α = –2.4 ± 0.8 and β = –2.0 ± 0.5, and an integrated occurrence rate of ${9}_{-4}^{+5}$% between 5–13M Jup and 10–100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with ${0.8}_{-0.5}^{+0.8}$% of stars hosting a brown dwarf companion between 13–80$M$ 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. 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, according to this work.

79 ASTRONOMY AND ASTROPHYSICS↗

GPI Spectra of HR8799 C, D, and E in H-K Bands with KLIP Forward Modeling

We demonstrate KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR8799, using PyKLIP. We report new and re-reduced spectrophotometry of HR8799 c, d, and e from H-K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting fake sources and recovering them over a range of parameters. The K1/K2 spectra for planets c and d are similar to previously published results from the same dataset. We also present a K band spectrum of HR8799e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We compare planets c, d, and e with M, L, and T-type field objects. All objects are consistent with low gravity mid-to-late L dwarfs, however, a lack of standard spectra for low gravity late L-type objects lead to poor fit for gravity. We place our results in context of atmospheric models presented in previous publications and discuss differences in the spectra of the three planets.

planets and satellites: gaseous planets { stars: i↗