pith. machine review for the scientific record. sign in

arxiv: 1703.09710 · v2 · submitted 2017-03-28 · 🌌 astro-ph.IM · astro-ph.EP· astro-ph.SR· physics.data-an· stat.AP

Recognition: unknown

Fast and scalable Gaussian process modeling with applications to astronomical time series

Authors on Pith no claims yet
classification 🌌 astro-ph.IM astro-ph.EPastro-ph.SRphysics.data-anstat.AP
keywords methodgaussianprocessastronomicaldatadatasetstimeanalysis
0
0 comments X
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.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 13 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. An Aligned Very-Low-Mass Star Orbiting an M dwarf and Obliquity Patterns Across Giant Planets, Brown Dwarfs, and Binary Stars

    astro-ph.EP 2026-04 unverdicted novelty 7.0

    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.

  2. An eclipsing CEMP candidate discovered in a search for dwarf carbon stars in post-common envelope binaries

    astro-ph.SR 2026-05 unverdicted novelty 6.0

    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 ...

  3. The 2MIG isolated AGNs. 3. Optical--IR variability and dust reverberation in the NLSy1 galaxies Mrk~42 and Mrk~493

    astro-ph.GA 2026-05 unverdicted novelty 6.0

    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.

  4. The NUV transit of XO-3 b

    astro-ph.EP 2026-05 unverdicted novelty 6.0

    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...

  5. A Statistical-AI Framework for Detecting Transient Flares in SDSS Stripe 82 Quasar Light Curves

    astro-ph.IM 2026-04 unverdicted novelty 6.0

    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.

  6. UMI: GPU-Accelerated Asymmetric Robust Estimator for Photometric Detrending in Exoplanet Transit Searches

    astro-ph.IM 2026-04 unverdicted novelty 6.0

    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.

  7. An Ultra-Short Period Super-Earth and a Sub-Neptune Orbiting the K dwarf TOI-4311

    astro-ph.EP 2026-05 unverdicted novelty 5.0

    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...

  8. Uncovering the Rapidly Evolving Orbits of the Dynamic TOI-201 System

    astro-ph.EP 2026-04 unverdicted novelty 5.0

    TOI-201 has three planets whose co-transiting configuration will end in 200 years due to Kozai-Lidov oscillations driven by mutual inclinations.

  9. A tidally detached super Neptune on a strongly misaligned retrograde orbit

    astro-ph.EP 2026-04 unverdicted novelty 5.0

    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.

  10. A universal relationship between the variability timescale and black hole mass in black hole jetted and non-jetted accreting systems

    astro-ph.HE 2026-04 unverdicted novelty 5.0

    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.

  11. Fundamental effective temperature measurements for eclipsing binary stars -- VIII. NIRPS spectroscopy of CD-27 2812

    astro-ph.SR 2026-04 unverdicted novelty 5.0

    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.

  12. Stochastic problems in pulsar timing

    astro-ph.HE 2026-04 unverdicted novelty 5.0

    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-...

  13. TOI-159 b: an eccentric hot-Jupiter planet around a young, pulsating $\gamma$ Doradus star

    astro-ph.EP 2026-05 accept novelty 4.0

    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.