Abstract:The use of Gaussian processes (GPs) is a common approach to account for correlated noise in exoplanet time series, particularly for transmission and emission spectroscopy. This analysis has typically been performed for each wavelength channel separately, with the retrieved uncertainties in the transmission spectrum assumed to be independent. However, the presence of noise correlated in wavelength could cause these uncertainties to be correlated, which could significantly affect the results of atmospheric retrievals. We present a method that uses a GP to model noise correlated in both wavelength and time simultaneously for the full spectroscopic dataset. To make this analysis computationally tractable, we introduce a new fast and flexible GP method that can analyse 2D datasets when the input points lie on a (potentially non-uniform) 2D grid - in our case a time by wavelength grid - and the kernel function has a Kronecker product structure. This simultaneously fits all light curves and enables the retrieval of the full covariance matrix of the transmission spectrum. Our new method can avoid the use of a `common-mode' correction, which is known to produce an offset to the transmission spectrum. Through testing on synthetic datasets, we demonstrate that our new approach can reliably recover atmospheric features contaminated by noise correlated in time and wavelength. In contrast, fitting each spectroscopic light curve separately performed poorly when wavelength-correlated noise was present. It frequently underestimated the uncertainty of the scattering slope and overestimated the uncertainty in the strength of sharp absorption peaks in transmission spectra. Two archival VLT/FORS2 transit observations of WASP-31b were used to compare these approaches on real observations. Our method strongly constrained the presence of wavelength-correlated noise in both datasets, and significantly different constraints on atmospheric features such as the scattering slope and strength of sodium and potassium features were recovered.

A Gaussian process cross-correlation approach to time delay estimation in active galactic nuclei

gaus SN: Bayesian time-delay estimation for strongly lensed supernovae

GausSN: Bayesian Time-Delay Estimation for Strongly Lensed Supernovae

Using Gaussian Processes to detect AGN flares

Disentangling the optical AGN and Host-galaxy luminosity with a probabilistic Flux Variation Gradient

How do wavelength correlations affect transmission spectra? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy

Expectations for time-delay measurements in active galactic nuclei with the Vera Rubin Observatory

Conditional Neural Process for Nonparametric Modeling of Active Galactic Nuclei Light Curves

TD-CARMA: Painless, accurate, and scalable estimates of gravitational-lens time delays with flexible CARMA processes

How do wavelength correlations affect your transmission spectrum? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy

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

Statistical recovery of 21cm visibilities and their power spectra with Gaussian constrained realisations and Gibbs sampling

First Combined Study on LIV from Observations of Energy-dependent Time Delays from Multiple-type Gamma-ray Sources -- Part I. Motivation, Method Description and Validation through Simulations of H.E.S.S., MAGIC and VERITAS Datasets

Actively Learning Numerical Relativity

Gaussian Process Foreground Subtraction and Power Spectrum Estimation for 21 cm Cosmology

A Gaussian process framework for modelling stellar activity signals in radial velocity data

On the accuracy of time-delay cosmography in the Frontier Fields Cluster MACS J1149.5+2223 with supernova Refsdal

A Gaussian process model for stellar activity in 2-D line profile time-series

Photometric AGN reverberation mapping - an efficient tool for BLR sizes, black hole masses and host-subtracted AGN luminosities

Searching for Quasi-periodic Oscillations in Astrophysical Transients Using Gaussian Processes