Magnitude-squared coherence: A powerful tool for disentangling Doppler planet discoveries from stellar activity

Abstract

If Doppler searches for earth-mass, habitable planets are to succeed, observers must be able to identify and model out stellar activity signals. Here we demonstrate how to diagnose activity signals by calculating the magnitude-squared coherence ${\hat{C}}_{xy}^2(f)$ between an activity indicator time series $x_t$ and the radial velocity (RV) time series $y_t$. Since planets only cause modulation in RV, not in activity indicators, a high value of ${\hat{C}}_{xy}^2(f)$ indicates that the signal at frequency $f$ has a stellar origin. We use Welch's method to measure coherence between activity indicators and RVs in archival observations of GJ 581, alpha Cen B, and GJ 3998. High RV-H$\alpha$ coherence at the frequency of GJ 3998 b, and high RV-S index coherence at the frequency of GJ 3998 c, indicate that the planets may actually be stellar signals. We also replicate previous results showing that GJ 581 d and g are rotation harmonics and demonstrate that alpha Cen B has activity signals that are not associated with rotation. Welch's power spectrum estimates have cleaner spectral windows than Lomb-Scargle periodograms, improving our ability to estimate rotation periods. We find that the rotation period of GJ 581 is 132 days, with no evidence of differential rotation. Welch's method may yield unacceptably large bias for datasets with $N<75$ observations and works best on datasets with $N>100$. Tapering the time-domain data can reduce the bias of the Welch's power spectrum estimator, but observers should not apply tapers to datasets with extremely uneven observing cadence. A software package for calculating magnitude-squared coherence and Welch's power spectrum estimates is available on github.

Submitted to AJ, responded to input from referee. 35 pages (29 in main body of article) include 21 figures

References in corpus (20)

• The generalised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms
• Assessing statistical significance of periodogram peaks
• Stellar Activity Masquerading as Planets in the Habitable Zone of the M dwarf Gliese 581
• The HARPS search for southern extra-solar planets. X. A m sin i = 11 Mearth planet around the nearby spotted M dwarf GJ 674
• TW Hydrae: evidence of stellar spots instead of a Hot Jupiter
• Search for exoplanets with the radial-velocity technique: quantitative diagnostics of stellar activity
• Fluxes in Hαand Ca II H and K for a sample of Southern stars
• BGLS: A Bayesian formalism for the generalised Lomb-Scargle periodogram
• The impact of red noise in radial velocity planet searches: Only three planets orbiting GJ581?
• Two planets around Kapteyn's star : a cold and a temperate super-Earth orbiting the nearest halo red-dwarf
• Stellar activity mimics a habitable-zone planet around Kapteyn's star
• Disentangling Planets and Stellar Activity for Gliese 667C
• The HADES RV Programme with HARPS-N@TNG - GJ 3998: An early M-dwarf hosting a system of Super-Earths
• The HADES RV Programme with HARPS-N@TNG II. Data treatment and simulations
• The Changing Face of $\alpha$ Centauri B: Probing plage and stellar activity in K-dwarfs
• Target Prioritization and Observing Strategies for the NEID Earth Twin Survey
• Auto-correlation functions of astrophysical processes, and their relation to Gaussian processes; Application to radial velocities of different starspot configurations
• A HARPS-N mass for the elusive Kepler-37d: a case study in disentangling stellar activity and planetary signals
• Robust detection of quasi-periodic variability: A HAWKI mini survey of late T dwarfs
• HADES RV programme with HARPS-N at TNG XIV. A candidate super-Earth orbiting the M-dwarf GJ 9689 with a period close to half the stellar rotation period