Fast and scalable Gaussian process modeling with applications to astronomical time series
read the original abstract
The growing field of large-scale time domain astronomy requires methods for probabilistic data analysis that are computationally tractable, even with large datasets. Gaussian Processes are a popular class of models used for this purpose but, since the computational cost scales, in general, as the cube of the number of data points, their application has been limited to small datasets. In this paper, we present a novel method for Gaussian Process modeling in one-dimension where the computational requirements scale linearly with the size of the dataset. We demonstrate the method by applying it to simulated and real astronomical time series datasets. These demonstrations are examples of probabilistic inference of stellar rotation periods, asteroseismic oscillation spectra, and transiting planet parameters. The method exploits structure in the problem when the covariance function is expressed as a mixture of complex exponentials, without requiring evenly spaced observations or uniform noise. This form of covariance arises naturally when the process is a mixture of stochastically-driven damped harmonic oscillators -- providing a physical motivation for and interpretation of this choice -- but we also demonstrate that it can be a useful effective model in some other cases. We present a mathematical description of the method and compare it to existing scalable Gaussian Process methods. The method is fast and interpretable, with a range of potential applications within astronomical data analysis and beyond. We provide well-tested and documented open-source implementations of this method in C++, Python, and Julia.
This paper has not been read by Pith yet.
Forward citations
Cited by 37 Pith papers
-
The 35-Myr old infant planet TOI-837 b has a mildly misaligned orbit
TOI-837 b has a true obliquity of 25.9+7.5-6.3 deg, the first planet younger than 100 Myr with accessible ψ incompatible with an aligned orbit, favoring primordial disc torque followed by disc-driven migration.
-
A hierarchical Bayesian framework for cosmology using Type 1 AGN variability
A hierarchical Bayesian framework that uses the empirical anti-correlation between AGN variability amplitude and luminosity to infer cosmological parameters from moderate-baseline light curves via importance reweighting.
-
Glossy Silicate Clouds on the Scorched Dayside of LTT9779b
First statistically significant detection of dayside silicate clouds on a Neptunian-mass exoplanet, with CO and CO2 detections and C/O ratio of 0.984 from JWST spectra.
-
KRONOS I: The $1{-}2.8\mu$m JWST Transmission Spectrum of the 23 Myr V1298 Tau c
JWST NIRISS/SOSS data yield a transmission spectrum of the 23 Myr V1298 Tau c showing H2O absorption and an atmospheric metallicity of 14.8^{+56}_{-12} times solar, lower than mature planets of similar mass.
-
MAGAZ3NE: Spatially Resolved Ages and Chemical Abundances of Ultra-Massive Quiescent Galaxies at z $\sim$ 3.5 using JWST/NIRSpec IFU
JWST IFU spectroscopy of three ultra-massive quiescent galaxies at z~3.5 shows young central ages with elevated [α/Fe] and mixed age/abundance gradients, indicating rapid star formation before quenching and the need f...
-
An Aligned Very-Low-Mass Star Orbiting an M dwarf and Obliquity Patterns Across Giant Planets, Brown Dwarfs, and Binary Stars
First obliquity measurement in an M dwarf binary shows alignment, with tentative evidence that aligned orbits around cool stars and wide separations also hold for brown dwarfs and binaries.
-
Exposure-averaged Gaussian Processes for Combining Overlapping Datasets
Exposure-integrated Gaussian processes allow prediction of both latent stellar signals and instrument-specific binned versions, supporting combination of overlapping EPRV datasets with varying exposure times.
-
The GAPS Programme at TNG LXXIV. A reanalysis of the planetary systems TOI-1272 and TOI-1694 with HARPS-N and retraction of the planetary interpretation of TOI-1272 c
Reanalysis retracts TOI-1272c as a planet due to stellar activity, finds slightly eccentric orbits for TOI-1694b and c, and provides updated parameters with smaller uncertainties for both systems.
-
The X-ray emission of the long-period transient and accreting cataclysmic variable ASKAP J174508.9-505149
X-ray timing and spectral analysis of ASKAP J174508.9-505149 detects matching periodicity and features consistent with an accreting magnetic CV.
-
Probabilistic Data-Driven Modelling of Astrophysical Transients: The Neural Process Family for Ultrafast and Class-Agnostic Light Curve Reconstruction with NightLANP
Attentive Neural Processes outperform Gaussian Processes and neural networks on light curve interpolation quality, feature recovery, calibration, and speed for 15 transient classes under realistic Rubin cadences.
-
An eclipsing CEMP candidate discovered in a search for dwarf carbon stars in post-common envelope binaries
The first eclipsing CEMP candidate is identified in a 1.224-day binary, with 3.5% of dwarf carbon stars showing short-period variability consistent with tidally locked post-common envelope systems and supporting wind ...
-
The 2MIG isolated AGNs. 3. Optical--IR variability and dust reverberation in the NLSy1 galaxies Mrk~42 and Mrk~493
Optical-MIR lags of 39.1 days and 79.4 days are measured in Mrk 42 and Mrk 493, yielding dust reverberation radii of 0.032 pc and 0.065 pc with R_dust/R_BLR ratios of approximately 6-7.
-
The NUV transit of XO-3 b
XO-3 b exhibits a 30-70% deeper NUV transit depth of 0.1371 and a 22-minute late center relative to the optical ephemeris, with an extremely low estimated mass-loss rate and a bow-shock model that predicts the opposit...
-
A Statistical-AI Framework for Detecting Transient Flares in SDSS Stripe 82 Quasar Light Curves
A modular framework combining physics-informed neural networks, Ornstein-Uhlenbeck fitting, extreme value theory, and vision-language models detects 51 transient flares in 9,258 SDSS Stripe 82 quasar light curves.
-
UMI: GPU-Accelerated Asymmetric Robust Estimator for Photometric Detrending in Exoplanet Transit Searches
UMI is an asymmetric version of the Tukey bisquare estimator with upper-RMS scaling that improves transit signal recovery while running much faster on GPUs.
-
Nested Sampling for ARIMA Model Selection in Astronomical Time-Series Analysis
Nested sampling applied to ARIMA models enables Bayesian order selection and parameter inference that recovers ground truth in simulations and fits stochastic variability in sunspot, Kepler, and TESS light curves.
-
A Gaussian process framework for testing general relativity with gravitational waves
A Gaussian process framework with a localized-deviation kernel is used to test general relativity on simulated and GWTC-3 binary black hole signals, finding no evidence for deviations and limiting strain deviations to...
-
Star-formation variability on the star-forming main sequence during the Epoch of Reionization
PSD modeling of SFR scatter at six timescales shows dominant variability on 10-30 Myr scales, stronger in lower-mass galaxies at z=3-8.
-
Stellar Obliquities of Young Systems, Atmospheres Undergoing Contraction and Escape (SOYSAUCE) II: a 135 Myr planet on an aligned orbit with transit timing variations
Validation of a 135 Myr, 3.6 R_E transiting planet with aligned obliquity and TTV evidence for a near-resonant companion.
-
An Ultra-Short Period Super-Earth and a Sub-Neptune Orbiting the K dwarf TOI-4311
TOI-4311 hosts a 0.99-day super-Earth (1.38 R_earth, 4.5 M_earth) and 15-day sub-Neptune (2.47 R_earth), plus a candidate 38-day planet, with the dense inner planet potentially challenging formation theories given the...
-
Uncovering the Rapidly Evolving Orbits of the Dynamic TOI-201 System
TOI-201 has three planets whose co-transiting configuration will end in 200 years due to Kozai-Lidov oscillations driven by mutual inclinations.
-
A tidally detached super Neptune on a strongly misaligned retrograde orbit
TOI-1710 b has a true obliquity of 149 degrees indicating retrograde motion, favoring high-eccentricity migration via planet-planet scattering and Kozai-Lidov cycles for this tidally detached super-Neptune.
-
A universal relationship between the variability timescale and black hole mass in black hole jetted and non-jetted accreting systems
Mass-scaled DRW damping timescales in AGNs follow a linear relation with black hole mass (slope 0.35-0.50) for both jetted and non-jetted sources, supporting universal accretion physics.
-
Fundamental effective temperature measurements for eclipsing binary stars -- VIII. NIRPS spectroscopy of CD-27 2812
Model-independent masses and radii plus effective temperatures are derived for the F9 V primary (M1=1.36 Msun, R1=1.72 Rsun, Teff=6197 K) and M-dwarf secondary (M2=0.56 Msun, R2=0.53 Rsun, Teff=3770 K) in CD-27 2812.
-
Stochastic problems in pulsar timing
Analytical solutions to Langevin equations for red noise and GWB in pulsars show that an Ornstein-Uhlenbeck spin frequency model is inconsistent with stationary signals, while an overdamped oscillator model and a two-...
-
A Systematic Search for Active Galactic Nucleus Flares in ZTF Data Release 23
Systematic search of ZTF DR23 data yields catalogs of 28,504 coarse and 1,984 refined AGN flares with public release.
-
Understanding eccentric temperate giants: an in-depth study of the architecture and stellar obliquity of the TOI-2134 system
Updated analysis of TOI-2134 with new TESS sectors and spectra confirms an inner mini-Neptune and outer eccentric sub-Saturn, measures their masses and radii, and reports a 59 degree obliquity for the outer planet via...
-
The GAPS programme at TNG ?. TOI-1533: a compact system hosting a super-Neptune-mass pair with disparate radii
TOI-1533 hosts an inner sub-Neptune (P=3.63 d, R=3.15 R⊕) and outer super-Neptune-mass hot giant (P=8.06 d, R>7.5 R⊕, M≈40 M⊕, ρ<0.48 g cm⁻³) both transiting an active K-dwarf.
-
Revisiting the Exo-Mercury Candidate GJ 367 b with ESPRESSO and a Self-Consistent Tidal Distortion Model
Revised mass of 0.503 M_Earth and radius of 0.736 R_Earth for GJ 367 b give a density of 6.9 g cm^{-3} and an iron fraction of 50-70% via new tidal and composition modeling.
-
The Pan-Pacific Planet Search -- IX. A menagerie of companions orbiting evolved stars
Resolves six speculative companions into one giant planet, one eccentric brown dwarf, two low-mass stars, and two stars with no detectable companions.
-
The NUV transit of XO-3 b
NUV transit depth of XO-3b measured at 0.1371 with 22-minute late center; X-ray data yield mass-loss rate ~10^4 g/s; bow-shock model predicts early rather than late transit.
-
TOI-159 b: an eccentric hot-Jupiter planet around a young, pulsating $\gamma$ Doradus star
TOI-159 b is confirmed as the hottest known eccentric hot Jupiter (e = 0.24) with a 13-sigma Keplerian detection around a young gamma Doradus star, including a preliminary low-resolution transmission spectrum.
-
POSEIDON I: The Dynamical Origins of Transiting Neptunes
New obliquity measurements for two Neptunes update the sample distribution to favor aligned systems plus a random component, resembling that of more massive planets and implying shared dynamical origins.
-
Understanding eccentric temperate giants: an in-depth study of the architecture and stellar obliquity of the TOI-2134 system
New observations confirm a mini-Neptune (9.23 d, circular) and eccentric sub-Saturn (95.85 d, e=0.31) around TOI-2134 with 59 deg obliquity for the outer planet.
-
The GAPS Programme at TNG LXXIV. A reanalysis of the planetary systems TOI-1272 and TOI-1694 with HARPS-N and retraction of the planetary interpretation of TOI-1272 c
Reanalysis of TOI-1272 and TOI-1694 retracts TOI-1272 c as planetary, attributes the signal to stellar activity via Gaussian process modeling, and refines orbital parameters for the systems.
-
TOI-2147 b and TOI-6019 b: Two eccentric warm Jupiters detected and characterized with TESS and MaHPS
Detection and characterization of two eccentric warm Jupiters TOI-2147 b (P=26.2 d, e=0.29, M=116 M⊕) and TOI-6019 b (P=14.5 d, e=0.48, M=149 M⊕) with TESS and MaHPS data, showing mildly inflated radii consistent with...
-
Star-Planet Interactions: Observational Techniques and Methods
This review chapter summarizes observational techniques and challenges for studying star-planet interactions.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.