A Fast, 2D Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization

Abstract

Gaussian processes (GPs) are commonly used as a model of stochastic variability in astrophysical time series. In particular, GPs are frequently employed to account for correlated stellar variability in planetary transit light curves. The efficient application of GPs to light curves containing thousands to tens of thousands of datapoints has been made possible by recent advances in GP methods, including the celerite method. Here we present an extension of the celerite method to two input dimensions, where, typically, the second dimension is small. This method scales linearly with the total number of datapoints when the noise in each large dimension is proportional to the same celerite kernel and only the amplitude of the correlated noise varies in the second dimension. We demonstrate the application of this method to the problem of measuring precise transit parameters from multiwavelength light curves and show that it has the potential to improve transit parameters measurements by orders of magnitude. Applications of this method include transit spectroscopy and exomoon detection, as well a broader set of astronomical problems.

arXiv commentSubmitted to AAS Journals Updated 7/22/20 to include links to code

References in corpus (42)

• Seven temperate terrestrial planets around the nearby ultracool dwarf star TRAPPIST-1
• Fast and scalable Gaussian process modeling with applications to astronomical time series
• On detecting terrestrial planets with timing of giant planet transits
• DHP Framework: Digital Health Passports Using Blockchain -- Use case on international tourism during the COVID-19 pandemic
• Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations?
• Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk
• Use and Abuse of the Fisher Information Matrix in the Assessment of Gravitational-Wave Parameter-Estimation Prospects
• GPyTorch: Blackbox Matrix-Matrix Gaussian Process Inference with GPU Acceleration
• Non-linear X-ray variability in X-ray binaries and active galaxies
• A Gaussian process framework for modelling stellar activity signals in radial velocity data
• The Flat Transmission Spectrum of the Super-Earth GJ1214b from Wide Field Camera 3 on the Hubble Space Telescope
• An Alternative Approach To Measuring Reverberation Lags in Active Galactic Nuclei
• The nature of the TRAPPIST-1 exoplanets
• Modeling reverberation mapping data II: dynamical modeling of the Lick AGN Monitoring Project 2008 dataset
• The connection between stellar granulation and oscillation as seen by the Kepler mission
• Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets
• Wotan: Comprehensive time-series de-trending in Python
• Analytic Approximations for Transit Light Curve Observables, Uncertainties, and Covariances
• Evidence for a Large Exomoon Orbiting Kepler-1625b
• Exoplanet Transit Spectroscopy Using WFC3: WASP-12 b, WASP-17 b, and WASP-19 b
• Detection of the ellipsoidal and the relativistic beaming effects in the CoRoT-3 lightcurve
• Student-t Processes as Alternatives to Gaussian Processes
• Large eccentricity, low mutual inclination: the three-dimensional architecture of a hierarchical system of giant planets
• Avocado: Photometric Classification of Astronomical Transients with Gaussian Process Augmentation
• The Hunt for Exomoons with Kepler (HEK): II. Analysis of Seven Viable Satellite-Hosting Planet Candidates
• Radial Velocity Observations and Light Curve Noise Modeling Confirm That Kepler-91b is a Giant Planet Orbiting a Giant Star
• No Evidence for Lunar Transit in New Analysis of Hubble Space Telescope Observations of the Kepler-1625 System
• A Sparse Gaussian Process Framework for Photometric Redshift Estimation
• A Class of Fast Methods for Processing Irregularly Sampled or Otherwise Inhomogeneous One-Dimensional Data
• Precise photometric transit follow-up observations of five close-in exoplanets : update on their physical properties
• Improving transit characterisation with Gaussian process modelling of stellar variability
• Standardizing Type Ia Supernova Absolute Magnitudes Using Gaussian Process Data Regression
• Gaussian Process modelling of granulation and oscillations in red-giant stars
• Robust Strong Lensing Time Delay Estimation
• Loose Ends for the Exomoon Candidate Host Kepler-1625b
• Mitigating flicker noise in high-precision photometry. I -- Characterization of the noise structure, impact on the inferred transit parameters, and predictions for CHEOPS observations
• The Stellar Variability Noise Floor for Transiting Exoplanet Photometry with PLATO
• Stellar pulsation and granulation as noise sources in exoplanet transit spectroscopy in the ARIEL space mission
• Next Generation Strong Lensing Time Delay Estimation with Gaussian Processes
• JexoSim: A time domain simulator of exoplanet transit spectroscopy with JWST
• Forward modeling the orbits of companions to pulsating stars from their light travel time variations
• The SDSS-HET Survey of Kepler Eclipsing Binaries. Description of the Survey and First Results

Cited by in corpus (18)

• exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series
• Efficient modeling of correlated noise. III. Scalable methods for jointly modeling several observables' time series with Gaussian processes
• Statistical methods for exoplanet detection with radial velocities
• Mira variables in the Milky Way's nuclear stellar disc: discovery and classification
• Modelling the Multiwavelength Variability of Mrk 335 using Gaussian Processes
• Applications of a Gaussian Process Framework for Modelling of High-Resolution Exoplanet Spectra
• Joint Modeling of Radial Velocities and Photometry with a Gaussian Process Framework
• Constraining the Milky Way's Ultraviolet to Infrared SED with Gaussian Process Regression
• A Generative Model for Quasar Spectra
• Model selection and signal extraction using Gaussian Process regression
• A Population Analysis of 20 Exoplanets Observed from the Optical to the Near-infrared Wavelengths with HST: Evidence for Widespread Stellar Contamination
• How do wavelength correlations affect transmission spectra? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy
• A statistical model of stellar variability. I. FENRIR: a physics-based model of stellar activity, and its fast Gaussian process approximation
• Sensitivity to Sub-Io-sized Exosatellite Transits in the MIRI LRS Lightcurve of the Nearest Substellar Worlds
• Scalable and Robust Multiband Modeling of AGN Light Curves in Rubin-LSST
• L 363-38 b: a planet newly discovered with ESPRESSO orbiting a nearby M dwarf star
• Scalable Gaussian Processes for Integrated and Overlapping Measurements Via Augmented State Space Models
• Modeling the Solar System I: Characterization Limits from Analytic Timing Variations