A statistical model of stellar variability. I. FENRIR: a physics-based model of stellar activity, and its fast Gaussian process approximation

Abstract

The detection of terrestrial planets by radial velocity and photometry is hindered by the presence of stellar signals. Those are often modeled as stationary Gaussian processes, whose kernels are based on qualitative considerations, which do not fully leverage the existing physical understanding of stars. Our aim is to build a formalism which allows to transfer the knowledge of stellar activity into practical data analysis methods. In particular, we aim at obtaining kernels with physical parameters. This has two purposes: better modelling signals of stellar origin to find smaller exoplanets, and extracting information about the star from the statistical properties of the data. We consider several observational channels such as photometry, radial velocity, activity indicators, and build a model called FENRIR to represent their stochastic variations due to stellar surface inhomogeneities. We compute analytically the covariance of this multi-channel stochastic process, and implement it in the S+LEAF framework to reduce the cost of likelihood evaluations from $O(N^3)$ to $O(N)$. We also compute analytically higher order cumulants of our FENRIR model, which quantify its non-Gaussianity. We obtain a fast Gaussian process framework with physical parameters, which we apply to the HARPS-N and SORCE observations of the Sun, and constrain a solar inclination compatible with the viewing geometry. We then discuss the application of our formalism to granulation. We exhibit non-Gaussianity in solar HARPS radial velocities, and argue that information is lost when stellar activity signals are assumed to be Gaussian. We finally discuss the origin of phase shifts between RVs and indicators, and how to build relevant activity indicators. We provide an open-source implementation of the FENRIR Gaussian process model with a Python interface.

arXiv commentSubmitted to Astronomy \& Astrophysics

References in corpus (43)

• Fast and scalable Gaussian process modeling with applications to astronomical time series
• 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
• A simple method to estimate radial velocity variations due to stellar activity using photometry
• A Gaussian process framework for modelling instrumental systematics: application to transmission spectroscopy
• Using the Sun to estimate Earth-like planets detection capabilities II. Impact of plages
• SOAP 2.0: A tool to estimate the photometric and radial velocity variations induced by stellar spots and plages
• The HARPS search for southern extra-solar planets. X. A m sin i = 11 Mearth planet around the nearby spotted M dwarf GJ 674
• Search for exoplanets with the radial-velocity technique: quantitative diagnostics of stellar activity
• SOAP. A tool for the fast computation of photometry and radial velocity induced by stellar spots
• Three years of Sun-as-a-star radial-velocity observations on the approach to solar minimum
• The HARPS search for southern extra-solar planets: XIV. Gl 176b, a super-Earth rather than a Neptune, and at a different period
• Three Years of HARPS-N High-Resolution Spectroscopy and Precise Radial Velocity Data for the Sun
• Are solar brightness variations faculae- or spot-dominated?
• Reconstructing the solar integrated radial velocity using MDI/SOHO
• Using the Sun to estimate Earth-like planets detection capabilities. V. Parameterizing the impact of solar activity components on radial velocities
• Rapid contraction of giant planets orbiting the 20 million-years old star V1298 Tau
• The HARPS search for southern extra-solar planets. XLIII. A compact system of four super-Earth planets orbiting HD 215152
• Unsigned magnetic flux as a proxy for radial-velocity variations in Sun-like stars
• Mapping stellar surfaces I: Degeneracies in the rotational light curve problem
• Efficient modeling of correlated noise II. A flexible noise model with fast and scalable methods
• Efficient modeling of correlated noise. III. Scalable methods for jointly modeling several observables' time series with Gaussian processes
• Toward Extremely Precise Radial Velocities: II. A Tool For Using Multivariate Gaussian Processes to Model Stellar Activity
• Mapping stellar surfaces II: An interpretable Gaussian process model for light curves
• Statistical methods for exoplanet detection with radial velocities
• The first close-up of the "flip-flop" phenomenon in a single star
• The Impact of Stellar Surface Magnetoconvection and Oscillations on the Detection of Temperate, Earth-Mass Planets Around Sun-Like Stars
• Activity time series of old stars from late F to early K. I. Simulating radial velocity, astrometry, photometry, and chromospheric emission
• Diagnosing large-scale stellar magnetic fields using PCA on spectropolarimetric data
• Mitigating flicker noise in high-precision photometry. I -- Characterization of the noise structure, impact on the inferred transit parameters, and predictions for CHEOPS observations
• The GAPS Programme with HARPS-N at TNG XVII. Line profile indicators and kernel regression as diagnostics of radial-velocity variations due to stellar activity in solar-like stars
• Latitude Distribution of Sunspots: Analysis Using Sunspot Data and A Dynamo Model
• Testing whether a signal is strictly periodic. Application to disentangling planets and stellar activity in radial velocities
• The HARPS search for southern extra-solar planets XLV. Two Neptune mass planets orbiting HD 13808: a study of stellar activity modelling's impact on planet detection
• Activity time series of old stars from late F to early K. II. Radial velocity jitter and exoplanet detectability
• A Fast, 2D Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization
• Deriving stellar inclination of slow rotators using stellar activity
• A maximum entropy approach to detect close-in giant planets around active stars
• Faculae cancel out on the surfaces of active Suns
• Long-Term Variations in the Growth and Decay Rates of Sunspot Groups
• Modelling stellar activity with Gaussian process regression networks
• Joint Modeling of Radial Velocities and Photometry with a Gaussian Process Framework
• How do wavelength correlations affect transmission spectra? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy

Cited by in corpus (3)

• Using Doppler Imaging to model stellar activity and search for planets around Sun-like stars
• TOI-2322: two transiting rocky planets close to the stellar rotation period and its first harmonic
• Exposure-averaged Gaussian Processes for Combining Overlapping Datasets