Joint Modeling of Radial Velocities and Photometry with a Gaussian Process Framework

Abstract

Developments in the stability of modern spectrographs have led to extremely precise instrumental radial velocity (RV) measurements. For most stars, the detection limit of planetary companions with these instruments is expected to be dominated by astrophysical noise sources such as starspots. Correlated signals caused by rotationally-modulated starspots can obscure or mimic the Doppler shifts induced by even the closest, most massive planets. This is especially true for young, magnetically active stars where stellar activity can cause fluctuation amplitudes of $\gtrsim$0.1 mag in brightness and $\gtrsim$100 m s${}^{-1}$ in RV semi-amplitudes. Techniques that can mitigate these effects and increase our sensitivity to young planets are critical to improving our understanding of the evolution of planetary systems. Gaussian processes (GPs) have been successfully employed to model and constrain activity signals in individual cases. However, a principled approach of this technique, specifically for the joint modeling of photometry and RVs, has not yet been developed. In this work, we present a GP framework to simultaneously model stellar activity signals in photometry and RVs that can be used to investigate the relationship between both time series. Our method, inspired by the ${\text{FF}}^{\prime }$ framework of (Aigrain et al. 2012), models spot-driven activity signals as the linear combinations of two independent latent GPs and their time derivatives. We also simulate time series affected by starspots by extending the $\text{starry}$ software (Luger et al. 2019) to incorporate time evolution of stellar features. Using these synthetic datasets, we show that our method can predict spot-driven RV variations with greater accuracy than other GP approaches.

19 pages, 10 figures

References in corpus (35)

• The Transiting Exoplanet Survey Satellite
• Planets and Stellar Activity: Hide and Seek in the CoRoT-7 system
• emcee v3: A Python ensemble sampling toolkit for affine-invariant MCMC
• SOAP 2.0: A tool to estimate the photometric and radial velocity variations induced by stellar spots and plages
• A planet within the debris disk around the pre-main-sequence star AU Microscopii
• Inferring probabilistic stellar rotation periods using Gaussian processes
• The Gaia-Kepler Stellar Properties Catalog. II. Planet Radius Demographics as a Function of Stellar Mass and Age
• Stellar Spectroscopy in the Near-infrared with a Laser Frequency Comb
• Four newborn planets transiting the young solar analog V1298 Tau
• The impact of red noise in radial velocity planet searches: Only three planets orbiting GJ581?
• Photoevaporation vs. core-powered mass-loss: model comparison with the 3D radius gap
• Radial-Velocity Fitting Challenge. II. First results of the analysis of the data set
• TESS Hunt for Young and Maturing Exoplanets (THYME) II: A 17 Myr Old Transiting Hot Jupiter in the Sco-Cen Association
• pyaneti II: A multidimensional Gaussian process approach to analysing spectroscopic time-series
• A hot Jupiter around the very active weak-line T Tauri star TAP 26
• Rotation of Low-Mass Stars in Taurus with K2
• Zodiacal Exoplanets in Time (ZEIT) VI: a three-planet system in the Hyades cluster including an Earth-sized planet
• Radial Velocity Observations and Light Curve Noise Modeling Confirm That Kepler-91b is a Giant Planet Orbiting a Giant Star
• K2 Discovers a Busy Bee: An Unusual Transiting Neptune Found in the Beehive Cluster
• Rapid contraction of giant planets orbiting the 20 million-years old star V1298 Tau
• The discovery and mass measurement of a new ultra-short-period planet: EPIC~228732031b
• The EXPRES Stellar Signals Project II. State of the Field in Disentangling Photospheric Velocities
• The HADES RV Programme with HARPS-N@TNG - GJ 3998: An early M-dwarf hosting a system of Super-Earths
• Efficient modeling of correlated noise II. A flexible noise model with fast and scalable methods
• Efficient modeling of correlated noise. III. Scalable methods for jointly modeling several observables' time series with Gaussian processes
• Quasi-periodic Gaussian Processes for stellar activity: from physical to kernel parameters
• Statistical methods for exoplanet detection with radial velocities
• Caught in the Act: Core-powered Mass-loss Predictions for Observing Atmospheric Escape
• Decoding the radial velocity variations of HD41248 with ESPRESSO
• Photometry as a proxy for stellar activity in radial velocity analyses
• Impact of Correlated Noise on the Mass Precision of Earth-analog Planets in Radial Velocity Surveys
• EXPRES. III. Revealing the Stellar Activity Radial Velocity Signature of $\epsilon$ Eridani with Photometry and Interferometry
• The Epoch of Giant Planet Migration Planet Search Program. I. Near-Infrared Radial Velocity Jitter of Young Sun-like Stars
• RVSPY -- Radial Velocity Survey for Planets around Young Stars. Target characterization and high-cadence survey
• A HARPS radial velocity search for planets in the Scorpius-Centaurus association. And its combination with the HARPS and SOPHIE young nearby stars (YNS) surveys

Cited by in corpus (4)

• A statistical model of stellar variability. I. FENRIR: a physics-based model of stellar activity, and its fast Gaussian process approximation
• gaspery: Optimized Scheduling of Radial Velocity Follow-Up Observations for Active Host Stars
• Scalable Gaussian Processes for Integrated and Overlapping Measurements Via Augmented State Space Models
• NIRPS tightens the mass estimate of GJ 3090 b and detects a planet near the stellar rotation period