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.

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

References in corpus (56)

• Astropy: A Community Python Package for Astronomy
• The generalised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms
• Catalog of Nearby Exoplanets
• No planet for HD 166435
• The LCES HIRES/Keck Precision Radial Velocity Exoplanet Survey
• The HARPS search for southern extra-solar planets XVIII. An Earth-mass planet in the GJ 581 planetary system
• Assessing statistical significance of periodogram peaks
• The HARPS-TERRA project I. Description of the algorithms, performance and new measurements on a few remarkable stars observed by HARPS
• Rotation and differential rotation of active Kepler stars
• 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
• Detectability of Extrasolar Planets in Radial Velocity Surveys
• Gliese 581d is the first discovered terrestrial-mass exoplanet in the habitable zone
• Long-term magnetic activity of a sample of M-dwarf stars from the HARPS program I. Comparison of activity indices
• The HARPS search for southern extra-solar planets. VI. A Neptune-mass planet around the nearby M dwarf Gl 581
• Search for exoplanets with the radial-velocity technique: quantitative diagnostics of stellar activity
• The Lick-Carnegie Exoplanet Survey: A 3.1 M_Earth Planet in the Habitable Zone of the Nearby M3V Star Gliese 581
• Fluxes in Hαand Ca II H and K for a sample of Southern stars
• Ghost in the time series: no planet for Alpha Cen B
• EXPRES: A Next Generation RV Spectrograph in the Search for Earth-like Worlds
• BGLS: A Bayesian formalism for the generalised Lomb-Scargle periodogram
• The 55 Cnc system reassessed
• H-alpha Activity of Old M Dwarfs: Stellar Cycles and Mean Activity Levels For 93 Low-Mass Stars in the Solar Neighborhood
• 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
• Bayesian Re-analysis of the Gliese 581 Exoplanet System
• The CARMENES search for exoplanets around M dwarfs: A low-mass planet in the temperate zone of the nearby K2-18
• Gliese 581g as a scaled-up version of Earth: atmospheric circulation simulations
• Radial Velocity Planet Detection Biases at the Stellar Rotational Period
• Disentangling Planets and Stellar Activity for Gliese 667C
• Evolution of coexisting long and short period stellar activity cycles
• Characterisation of the radial velocity signal induced by rotation in late-type dwarfs
• The Impact of Stellar Rotation on the Detectability of Habitable Planets Around M Dwarfs
• The CARMENES search for exoplanets around M dwarfs: Radial-velocity variations of active stars in visual-channel spectra
• New Methods for Finding Activity-Sensitive Spectral Lines: Combined Visual Identification and an Automated Pipeline Find a Set of 40 Activity Indicators
• X-Ray, FUV, and UV Observations of alpha Centauri B: Determination of Long-term Magnetic Activity Cycle and Rotation Period
• 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
• A Goldilocks principle for modeling radial velocity noise
• The Changing Face of $\alpha$ Centauri B: Probing plage and stellar activity in K-dwarfs
• Two New Long-Period Giant Planets from the McDonald Observatory Planet Search and Two Stars with Long-Period Radial Velocity Signals Related to Stellar Activity Cycles
• Vertical Coherence of Turbulence in the Atmospheric Surface Layer: Connecting the Hypotheses of Townsend and Davenport
• Stellar Activity and Exclusion of the Outer Planet in the HD 99492 System
• 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
• GJ 581 update: Additional Evidence for a Super-Earth in the Habitable Zone
• Stellar activity analysis of Barnard's Star: Very slow rotation and evidence for long-term activity cycle
• A HARPS-N mass for the elusive Kepler-37d: a case study in disentangling stellar activity and planetary signals
• A Gaussian Process Regression Reveals No Evidence for Planets Orbiting Kapteyn's Star
• Robust detection of quasi-periodic variability: A HAWKI mini survey of late T dwarfs
• Giant Outer Transiting Exoplanet Mass (GOT 'EM) Survey. II. Discovery of a Failed Hot Jupiter on a 2.7 Year, Highly Eccentric Orbit
• HADES RV Programme with HARPS-N at TNG XIII. A sub-Neptune around the M dwarf GJ 720 A
• Habitability of the Goldilocks Planet Gliese 581g: Results from Geodynamic Models
• 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

Cited by in corpus (4)

• Gaia DR3 astrometric orbit determination with Markov Chain Monte Carlo and Genetic Algorithms. Systems with stellar, substellar, and planetary mass companions
• Optimal frequency-domain analysis for spacecraft time series: Introducing the missing-data multitaper power spectrum estimator
• HADES RV Programme with HARPS-N at TNG. XVI. A super-Earth in the habitable zone of the GJ 3998 multi-planet system
• Dual-Backend Multibeam Position Switching Targeted SETI Observations toward Nearby Active Planet-Hosting Systems with FAST