Mapping stellar surfaces II: An interpretable Gaussian process model for light curves
arXiv:2102.01697 · doi:10.3847/1538-3881/abfdb9
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.
arXiv comment91 pages, 30 figures. Updated references. Interactive visualization at http://starry-process.flatironinstitute.org
References in corpus (38)
- MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
- dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
- Fast and scalable Gaussian process modeling with applications to astronomical time series
- Efficient, uninformative sampling of limb darkening coefficients for two-parameter laws
- Planets and Stellar Activity: Hide and Seek in the CoRoT-7 system
- A Gaussian process framework for modelling stellar activity signals in radial velocity data
- The Transit Light Source Effect: False Spectral Features and Incorrect Densities for M-dwarf Transiting Planets
- A seven-planet resonant chain in TRAPPIST-1
- STARRY: Analytic Occultation Light Curves
- The solar magnetic field
- EVEREST: Pixel Level Decorrelation of K2 Light curves
- Inferring probabilistic stellar rotation periods using Gaussian processes
- Automatic Differentiation Variational Inference
- K2SC: Flexible systematics correction and detrending of K2 light curves using Gaussian Process regression
- Asteroseismology Can Reveal Strong Internal Magnetic Fields in Red Giant Stars
- Placing the spotted T Tauri star LkCa 4 on an HR diagram
- Towards precise stellar ages: combining isochrone fitting with empirical gyrochronology
- A Volume Limited Survey of mCP Stars Within 100pc II: Rotational and Magnetic Properties
- Radial Velocity Planet Detection Biases at the Stellar Rotational Period
- The radii and limb darkenings of Alpha Centauri A and B - Interferometric measurements with VLTI/PIONIER
- Envelope Convection, Surface Magnetism and Spots in A- and late B-type Stars
- Disk-locking and the presence of slow rotators among solar-type stars in young star clusters
- Formation of starspots in self-consistent global dynamo models: Polar spots on cool stars
- Possible Bright Starspots on TRAPPIST-1
- Mapping stellar surfaces I: Degeneracies in the rotational light curve problem
- Planet-Planet Occultations in TRAPPIST-1 and Other Exoplanet Systems
- A Relationship Between Stellar Age and Spot Coverage
- Mapping stellar surfaces II: An interpretable Gaussian process model for light curves
- Toward Extremely Precise Radial Velocities: II. A Tool For Using Multivariate Gaussian Processes to Model Stellar Activity
- Persistent starspot signals on M dwarfs: multi-wavelength Doppler observations with the Habitable-zone Planet Finder and Keck/HIRES
- Gaussian Process Modelling of Asteroseismic Data
- Auto-correlation functions of astrophysical processes, and their relation to Gaussian processes; Application to radial velocities of different starspot configurations
- Improving Exoplanet Detection Power: Multivariate Gaussian Process Models for Stellar Activity
- The Radius and Entropy of a Magnetized, Rotating Fully-convective Star: Analysis With Depth-dependent Mixing Length Theories
- Modeling the Rise of Fibril Magnetic Fields in Fully Convective Stars
- The effect of spots on the luminosity spread of the Pleiades
- Biases in retrieving planetary signals in the presence of quasi-periodic stellar activity
- starry_process: Interpretable Gaussian processes for stellar light curves
Cited by in corpus (27)
- Mapping stellar surfaces I: Degeneracies in the rotational light curve problem
- Quasi-periodic Gaussian Processes for stellar activity: from physical to kernel parameters
- HST/WFC3 transmission spectroscopy of the cold rocky planet TRAPPIST-1h
- Rejuvenated accretors have less bound envelopes: Impact of Roche lobe overflow on subsequent common envelope events
- Mapping stellar surfaces II: An interpretable Gaussian process model for light curves
- Starspots and Magnetism: Testing the Activity Paradigm in the Pleiades and M67
- Starspot mapping with adaptive parallel tempering. II. Application to TESS data for M-dwarf flare stars, AU Microscopii, YZ Canis Minoris, and EV Lacertae
- Gaussian Process Regression for foreground removal in HI intensity mapping experiments
- Three K2 Campaigns Yield Rotation Periods for 1013 Stars in Praesepe
- The CARMENES search for exoplanets around M dwarfs: Stable radial-velocity variations at the rotation period of AD~Leonis -- A test case study of current limitations to treating stellar activity
- A roadmap for the atmospheric characterization of terrestrial exoplanets with JWST
- TOI-3884 b: A rare 6-R planet that transits a low-mass star with a giant and likely polar spot
- Searching for quasi-periodic oscillations in astrophysical transients using Gaussian processes
- EXPRES. III. Revealing the Stellar Activity Radial Velocity Signature of Eridani with Photometry and Interferometry
- Analytic Light Curves in Reflected Light: Phase Curves, Occultations, and Non-Lambertian Scattering for Spherical Planets and Moons
- Characterization of starspots on a young M-dwarf K2-25: multi-band observations of stellar photometric variability and planetary transits
- A Panchromatic Characterization of the Evening and Morning Atmosphere of WASP-107 b: Composition and Cloud Variations, and Insight into the Effect of Stellar Contamination
- Hazy with a chance of star spots: constraining the atmosphere of the young planet, K2-33b
- Updated Planetary Mass Constraints of the Young V1298 Tau System Using MAROON-X
- Searching for stellar cycles on low mass stars using TESS data
- A statistical model of stellar variability. I. FENRIR: a physics-based model of stellar activity, and its fast Gaussian process approximation
- Using Doppler Imaging to model stellar activity and search for planets around Sun-like stars
- Stellar activity and transits
- Flaring Latitudes in Ensembles of Low Mass Stars
- Distinguishing between light curves of ellipsoidal variables with massive dark companions, contact binaries, and semidetached binaries using principal component analysis
- gallifrey: JAX-based Gaussian process structure learning for astronomical time series
- Characterisation of starspot structure and differential rotation of Kepler-411