Abstract:The use of Gaussian processes (GPs) as models for astronomical time series datasets has recently become almost ubiquitous, given their ease of use and flexibility. GPs excel in particular at marginalization over the stellar signal in cases where the variability due to starspots rotating in and out of view is treated as a nuisance, such as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest since it is not obvious how the parameters of the GP model are related to the physical properties of interest, such as the size, contrast, and latitudinal distribution of the spots. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is degenerate, ill-posed, and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many light curves. In this paper, we derive a closed-form expression for the mean and covariance of a Gaussian process model that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more stellar light curves, including the typical radii and the mean and variance of the latitude distribution of starspots. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from both light curves and radial velocity (RV) measurements, as well as for robustly modeling correlated noise in both transiting and RV exoplanet searches. Our implementation is efficient, user-friendly, and open source, available as the Python package starry-process.

State-space representation of Matérn and Damped Simple Harmonic Oscillator Gaussian processes

Markovian Gaussian Process: A Universal State-Space Representation for Stationary Temporal Gaussian Process

Identification of Gaussian Process State Space Models

Hausdorff Dimension of Set Generated by Exceptional Oscillations of a Class of N-parameter Gaussian Processes

State Space representation of non-stationary Gaussian Processes

State Space Models for Gaussian Stochastic Processes.

Infinite-Horizon Gaussian Processes

Physics-Informed Variational State-Space Gaussian Processes

Towards Efficient Modeling and Inference in Multi-Dimensional Gaussian Process State-Space Models

Data-Driven Momentum Observers With Physically Consistent Gaussian Processes

Polynomial Representation of the Gaussian Process

Short-term Prediction and Filtering of Solar Power Using State-Space Gaussian Processes

Financial Time Series Volatility Analysis Using Gaussian Process State-Space Models

Sequential Estimation of Gaussian Process-based Deep State-Space Models

Latent Gaussian Process Regression

Mixed-Stationary Gaussian Process for Flexible Non-Stationary Modeling of Spatial Outcomes

Gaussian Processes with State-Dependent Noise for Stochastic Control

High-dimensional Gaussian linear processes: Marchenko-Pastur beyond simultaneous diagonalizability

Images of Gaussian and other stochastic processes under closed, densely-defined, unbounded linear operators

Mapping stellar surfaces II: An interpretable Gaussian process model for light curves

Non-stationary Gaussian Process Surrogates