Stellar Activity Signals Research Articles

Stacked Bayesian general Lomb-Scargle periodogram: Identifying stellar activity signals

Context. Distinguishing between a signal induced by either stellar activity or a planet is currently the main challenge in radial velocity searches for low-mass exoplanets. Even when the presence of a transiting planet and hence its period are known, stellar activity can be the main barrier to measuring the correct amplitude of the radial velocity signal. Several tools are being used to help understand which signals come from stellar activity in the data. Aims. We aim to present a new tool that can be used for the purpose of identifying periodicities caused by stellar activity, and show how it can be used to track the signal-to-noise ratio (SNR) of the detection over time. The tool is based on the principle that stellar activity signals are variable and incoherent. Methods. We calculate the Bayesian general Lomb-Scargle periodogram for subsets of data and by adding one extra data point we track what happens to the presence and significance of periodicities in the data. Publicly available datasets from HARPS and HARPS-N were used for this purpose. Additionally, we analysed a synthetic dataset that we created with SOAP2.0 to simulate pure stellar activity and a mixture of stellar activity and a planet. Results. We find that this tool can easily be used to identify unstable and incoherent signals, such as those introduced by stellar activity. The SNR of the detection grows approximately as the square root of the number of data points, in the case of a stable signal. This can then be used to make decisions on whether it is useful to keep observing a specific object. The tool is relatively fast and easy to use, and thus lends itself perfectly to a quick analysis of the data.

A COMBINED SPECTROSCOPIC AND PHOTOMETRIC STELLAR ACTIVITY STUDY OF EPSILON ERIDANI

ABSTRACT We present simultaneous ground-based radial velocity (RV) measurements and space-based photometric measurements of the young and active K dwarf Epsilon Eridani. These measurements provide a data set for exploring methods of identifying and ultimately distinguishing stellar photospheric velocities from Keplerian motion. We compare three methods we have used in exploring this data set: Dalmatian, an MCMC spot modeling code that fits photometric and RV measurements simultaneously; the FF′ method, which uses photometric measurements to predict the stellar activity signal in simultaneous RV measurements; and Hα analysis. We show that our Hα measurements are strongly correlated with the Microvariability and Oscillations of STars telescope (MOST) photometry, which led to a promising new method based solely on the spectroscopic observations. This new method, which we refer to as the HH′ method, uses Hα measurements as input into the FF′ model. While the Dalmatian spot modeling analysis and the FF′ method with MOST space-based photometry are currently more robust, the HH′ method only makes use of one of the thousands of stellar lines in the visible spectrum. By leveraging additional spectral activity indicators, we believe the HH′ method may prove quite useful in disentangling stellar signals.

An apodized Kepler periodogram for separating planetary and stellar activity signals.

A new apodized Keplerian (AK) model is proposed for the analysis of precision radial velocity (RV) data to model both planetary and stellar activity (SA) induced RV signals. A symmetrical Gaussian apodization function with unknown width and centre can distinguish planetary signals from SA signals on the basis of the span of the apodization window. The general model for m AK signals includes a linear regression term between RV and the SA diagnostic log (R′hk), as well as an extra Gaussian noise term with unknown standard deviation. The model parameters are explored using a Bayesian fusion Markov chain Monte Carlo code. A differential version of the generalized Lomb–Scargle periodogram that employs a control diagnostic provides an additional way of distinguishing SA signals and helps guide the choice of new periods. Results are reported for a recent international RV blind challenge which included multiple state-of-the-art simulated data sets supported by a variety of SA diagnostics. In the current implementation, the AK method achieved a reduction in SA noise by a factor of approximately 6. Final parameter estimates for the planetary candidates are derived from fits that include AK signals to model the SA components and simple Keplerians to model the planetary candidates. Preliminary results are also reported for AK models augmented by a moving average component that allows for correlations in the residuals.

Ghost in the time series: no planet for Alpha Cen B

Abstract We re-analyse the publicly available radial velocity (RV) measurements for Alpha Cen B, a star hosting an Earth-mass planet candidate, Alpha Cen Bb, with 3.24 d orbital period. We demonstrate that the 3.24 d signal observed in the Alpha Cen B data almost certainly arises from the window function (time sampling) of the original data. We show that when stellar activity signals are removed from the RV variations, other significant peaks in the power spectrum of the window function are coincidentally suppressed, leaving behind a spurious yet apparently significant ‘ghost’ of a signal that was present in the window function's power spectrum ab initio. Even when fitting synthetic data with time sampling identical to the original data, but devoid of any genuine periodicities close to that of the planet candidate, the original model used to infer the presence of Alpha Cen Bb leads to identical conclusions: viz., the 3σ detection of a half-a-metre-per-second signal with 3.236 d period. Our analysis underscores the difficulty of detecting weak planetary signals in RV data, and the importance of understanding in detail how every component of an RV data set, including its time sampling, influences final statistical inference.

Bayesian Planet Searches for the 10 cm/s Radial Velocity Era

AbstractA new apodized Keplerian model is proposed for the analysis of precision radial velocity (RV) data to model both planetary and stellar activity (SA) induced RV signals. A symmetrical Gaussian apodization function with unknown width and center can distinguish planetary signals from SA signals on the basis of the width of the apodization function. The general model for m apodized Keplerian signals also includes a linear regression term between RV and the stellar activity diagnostic In (R'hk), as well as an extra Gaussian noise term with unknown standard deviation. The model parameters are explored using a Bayesian fusion MCMC code. A differential version of the Generalized Lomb-Scargle periodogram provides an additional way of distinguishing SA signals and helps guide the choice of new periods. Sample results are reported for a recent international RV blind challenge which included multiple state of the art simulated data sets supported by a variety of stellar activity diagnostics.

A Gaussian process framework for modelling stellar activity signals in radial velocity data

To date, the radial velocity (RV) method has been one of the most productive techniques for detecting and confirming extrasolar planetary candidates. Unfortunately, stellar activity can induce RV variations which can drown out or even mimic planetary signals - and it is notoriously difficult to model and thus mitigate the effects of these activity-induced nuisance signals. This is expected to be a major obstacle to using next-generation spectrographs to detect lower mass planets, planets with longer periods, and planets around more active stars. Enter Gaussian processes (GPs) which, we note, have a number of attractive features that make them very well suited to disentangling stellar activity signals from planetary signals. We present here a GP framework we developed to model RV time series jointly with ancillary activity indicators (e.g. bisector velocity spans, line widths, chromospheric activity indices), allowing the activity component of RV time series to be constrained and disentangled from e.g. planetary components. We discuss the mathematical details of our GP framework, and present results illustrating its encouraging performance on both synthetic and real RV datasets, including the publicly-available Alpha Centauri B dataset.

Interferometric diagnostics for magnetic spots detection

The signature of activity in general, and of stellar magnetic spots in particular, is present in every measurements, including interferometric ones. Indeed, stellar spots can be found on many stellar surfaces, their size and number varying according to their host's magnetic field and rotational velocity. To correctly determine stellar parameters, it is thus necessary to determine and extract stellar activity's signals. Interferometric observables are disturbed by activity, and this observing technique thus constitutes a good way of probing stellar surface. However, magnetic spots sometimes mimic other phenomenon, like a transiting exoplanet. In that case, the combination of several observing techniques, like photometry and interferometry, is mandatory to extract the planetary signal from the spot's one, and then characterize the exoplanet.

THE RADIAL VELOCITY DETECTION OF EARTH-MASS PLANETS IN THE PRESENCE OF ACTIVITY NOISE: THE CASE OF α CENTAURI Bb

We present an analysis of the publicly available HARPS radial velocity (RV) measurements for Alpha Cen B, a star hosting an Earth-mass planet candidate in a 3.24 day orbit. The goal is to devise robust ways of extracting low-amplitude RV signals of low mass planets in the presence of activity noise. Two approaches were used to remove the stellar activity signal which dominates the RV variations: 1) Fourier component analysis (pre-whitening), and 2) local trend filtering (LTF) of the activity using short time windows of the data. The Fourier procedure results in a signal at P = 3.236 days and K = 0.42 m/s which is consistent with the presence of an Earth-mass planet, but the false alarm probability for this signal is rather high at a few percent. The LTF results in no significant detection of the planet signal, although it is possible to detect a marginal planet signal with this method using a different choice of time windows and fitting functions. However, even in this case the significance of the 3.24-d signal depends on the details of how a time window containing only 10% of the data is filtered. Both methods should have detected the presence of Alpha Cen Bb at a higher significance than is actually seen. We also investigated the influence of random noise with a standard deviation comparable to the HARPS data and sampled in the same way. The distribution of the noise peaks in the period range 2.8 - 3.3 days have a maximum of approximately 3.2 days and amplitudes approximately one-half of the K-amplitude for the planet. The presence of the activity signal may boost the velocity amplitude of these signals to values comparable to the planet. It may be premature to attribute the 3.24 day RV variations to an Earth-mass planet. A better understanding of the noise characteristics in the RV data as well as more measurements with better sampling will be needed to confirm this exoplanet.

Disentangling between stellar activity and planetary signals

Photospheric stellar activity (i.e. dark spots or bright plages) might be an important source of noise and confusion in stellar radial-velocity (RV) measurements. Radial-velocimetry planet search surveys as well as follow-up of photometric transit surveys require a deeper understanding and characterization of the effects of stellar activities to differentiate them from planetary signals. We simulate dark spots on a rotating stellar photosphere. The variations in the photometry, RV, and spectral line shapes are characterized and analyzed according to the stellar inclination, the latitude, and the number of spots. We show that the anti-correlation between RV and bisector span, known to be a signature of activity, requires a good sampling to be resolved when there are several spots on the photosphere. The Lomb-Scargle periodograms of the RV variations induced by activity present power at the rotational period P rot of the star and its two first harmonics P rot /2 and P rot /3. Three adjusted sinusoids fixed at the fundamental period and its two-first harmonics allow us to remove about 90% of the RV jitter amplitude. We apply and validate our approach on four known active planet-host stars: HD 189733, GJ 674, CoRoT-7, and ι Hor. We succeed in fitting simultaneously activity and planetary signals on GJ674 and CoRoT-7. This simultaneous modeling of the activity and planetary parameters leads to slightly higher masses of CoRoT-7b and c of respectively, 5.7 ± 2.5 M Earth and 13.2 ± 4.1 M Earth . The larger uncertainties properly take into account the stellar active jitter. We exclude short-period low-mass exoplanets around ι Hor. For data with realistic time-sampling and white Gaussian noise, we use simulations to show that our approach is effective in distinguishing reflex-motion due to a planetary companion and stellar-activity-induced RV variations provided that 1) the planetary orbital period is not close to that of the stellar rotation or one of its two first harmonics; 2) the semi-amplitude of the planet exceeds ~ 30% of the semi-amplitude of the active signal; 3) the rotational period of the star is accurately known, and 4) the data cover more than one stellar rotational period.

Disentangling stellar activity and planetary signals

Photospheric stellar activity (i.e. dark spots or bright plages) might be an important source of noise and confusion in the radial-velocity (RV) measurements. Radial-velocimetry planet search surveys as well as follow-up of photometric transit surveys require a deep understanding and precise characterization of the effects of stellar activity, in order to disentangle it from planetary signals. We simulate dark spots on a rotating stellar photosphere. The variations of the RV are characterized and analyzed according to the stellar inclination, the latitude and the number of spots. The Lomb-Scargle periodograms of the RV variations induced by activity present power at the rotational period Prot of the star and its two-first harmonics Prot /2 and Prot /3. Three adjusted sinusoids fixed at the fundamental period and its two-first harmonics allow to remove about 90% of the RV jitter amplitude. We apply and validate our approach on four known active planet-host stars: HD 189733, GJ 674, CoRoT-7 and ι Hor.

Disentangling stellar activity and planetary signals

Disentangling stellar activity and planetary signals

Disentangling stellar activity and planetary signals

AbstractHigh-precision radial-velocimetry (RV) is until now the more efficient way to discover planetary systems. Moreover, photometric transit search missions like CoRoT and Kepler, need spectroscopic RV measurements to establish the planetary nature of a transit candidate and to measure the true mass. An active star has on its photosphere dark spots and bright plages rotating with the star. These inhomogeneities of the stellar surface can induce a variation of the measurement of the RV, due to changes in lines shapes and not to a Doppler motion of the star (e.g. Queloz et al. 2001; Desort et al. 2007; Boisse et al. 2009). We study how the Keplerian fit used to search for planets in RV data is confused by spots and we test an approach to subtract RV jitter based on harmonic decomposition of the star rotation. We use simulations of spectroscopic measurements of rotating spotted stars and validate our approach on active stars monitored by high-precision spectrograph HARPS: CoRoT-7 and ι Hor.

Disentangling between stellar activity and planetary signals

AbstractPhotospheric stellar activity (i.e. dark spots or bright plages) might be an important source of noise and confusion in the radial-velocity (RV) measurements. Radial-velocimetry planet search surveys as well as follow-up of photometric transit surveys require a deeper understanding and characterization of the effects of stellar activities to disentangle it from planetary signals.We simulate dark spots on a rotating stellar photosphere. The variations of the RV are characterized and analyzed according to the stellar inclination, the latitude and the number of spots. The Lomb-Scargle periodograms of the RV variations induced by activity present power at the rotational period Prot of the star and its two-first harmonics Prot/2 and Prot/3. Three adjusted sinusoids fixed at the fundamental period and its two-first harmonics allow to remove about 90% of the RV jitter amplitude. We apply and validate our approach on four known active planet-host stars: HD 189733, GJ 674, CoRoT-7 and ι Hor. We succeed in fitting simultaneously activity and planetary signals on GJ674 and CoRoT-7. We excluded short-period low-mass exoplanets around ι Hor. Our approach is efficient to disentangle reflex-motion due to a planetary companion and stellar-activity induced-RV variations provided that 1) the planetary orbital period is not close to that of the stellar rotation or one of its two-first harmonics, 2) the rotational period of the star is accurately known, 3) the data cover more than one stellar rotational period.

The CoRoT-7 planetary system: two orbiting super-Earths

We report on an intensive observational campaign carried out with HARPS at the 3.6 m telescope at La Silla on the star CoRoT-7. Additional simultaneous photometric measurements carried out with the Euler Swiss telescope have demonstrated that the observed radial velocity variations are dominated by rotational modulation from cool spots on the stellar surface. Several approaches were used to extract the radial velocity signal of the planet(s) from the stellar activity signal. First, a simple pre-whitening procedure was employed to find and subsequently remove periodic signals from the complex frequency structure of the radial velocity data. The dominant frequency in the power spectrum was found at 23 days, which corresponds to the rotation period of CoRoT-7. The 0.8535 day period of CoRoT-7b planetary candidate was detected with an amplitude of 3.3 m s −1 . Most other frequencies, some with amplitudes larger than the CoRoT-7b signal, are most likely associated with activity. A second approach used harmonic decomposition of the rotational period and up to the first three harmonics to filter out the activity signal from radial velocity variations caused by orbiting planets. After correcting the radial velocity data for activity, two periodic signals are detected: the CoRoT-7b transit period and a second one with a period of 3.69 days and an amplitude of 4 m s −1 . This second signal was also found in the pre-whitening analysis. We attribute the second signal to a second, more remote planet CoRoT-7c . The orbital solution of both planets is compatible with circular orbits. The mass of CoRoT-7b is 4.8 ± 0. 8( M⊕) and that of CoRoT-7c is 8.4 ± 0. 9( M⊕), assuming both planets are on coplanar orbits. We also investigated the false positive scenario of a blend by a faint stellar binary, and this may be rejected by the stability of the bisector on a nightly scale. According to their masses both planets belong to the super-Earth planet category. The average density of CoRoT-7b is ρ = 5.6 ± 1. 3gc m −3 , similar to the Earth. The CoRoT-7 planetary system provides us with the first insight into the physical nature of short period super-Earth planets recently detected by radial velocity surveys. These planets may be denser than Neptune and therefore likely made of rocks like the Earth, or a mix of water ice and rocks.