Inferring probabilistic stellar rotation periods using Gaussian processes

Abstract

Variability in the light curves of spotted, rotating stars is often non-sinusoidal and quasi-periodic --- spots move on the stellar surface and have finite lifetimes, causing stellar flux variations to slowly shift in phase. A strictly periodic sinusoid therefore cannot accurately model a rotationally modulated stellar light curve. Physical models of stellar surfaces have many drawbacks preventing effective inference, such as highly degenerate or high-dimensional parameter spaces. In this work, we test an appropriate effective model: a Gaussian Process with a quasi-periodic covariance kernel function. This highly flexible model allows sampling of the posterior probability density function of the periodic parameter, marginalising over the other kernel hyperparameters using a Markov Chain Monte Carlo approach. To test the effectiveness of this method, we infer rotation periods from 333 simulated stellar light curves, demonstrating that the Gaussian process method produces periods that are more accurate than both a sine-fitting periodogram and an autocorrelation function method. We also demonstrate that it works well on real data, by inferring rotation periods for 275 Kepler stars with previously measured periods. We provide a table of rotation periods for these 1132 Kepler objects of interest and their posterior probability density function samples. Because this method delivers posterior probability density functions, it will enable hierarchical studies involving stellar rotation, particularly those involving population modelling, such as inferring stellar ages, obliquities in exoplanet systems, or characterising star-planet interactions. The code used to implement this method is available online.

Submitted to MNRAS. Replaced 27/06/2017: corrections made to koi_periods.csv

References in corpus (12)

• Most 1.6 Earth-Radius Planets are not Rocky
• Planets and Stellar Activity: Hide and Seek in the CoRoT-7 system
• Measuring the rotation period distribution of field M-dwarfs with Kepler
• Exoplanet population inference and the abundance of Earth analogs from noisy, incomplete catalogs
• Inferring the eccentricity distribution
• Calibrating Gyrochronology using Kepler Asteroseismic targets
• Testing the recovery of stellar rotation signals from Kepler light curves using a blind hare-and-hounds exercise
• Large eccentricity, low mutual inclination: the three-dimensional architecture of a hierarchical system of giant planets
• Characterizing K2 Planet Discoveries: A super-Earth transiting the bright K-dwarf HIP 116454
• Radial Velocity Observations and Light Curve Noise Modeling Confirm That Kepler-91b is a Giant Planet Orbiting a Giant Star
• Gaussian Process Modelling of Asteroseismic Data
• Robust detection of quasi-periodic variability: A HAWKI mini survey of late T dwarfs