SOAP 2.0: A tool to estimate the photometric and radial velocity variations induced by stellar spots and plages

This paper presents SOAP 2.0, a new version of the SOAP code that estimates in a simple way the photometric and radial velocity variations induced by active regions. The inhibition of the convective blueshift inside active regions is considered, as well as the limb brightening effect of plages, a quadratic limb darkening law, and a realistic spot and plage contrast ratio. SOAP 2.0 shows that the activity-induced variation of plages is dominated by the inhibition of the convective blueshift effect. For spots, this effect becomes significant only for slow rotators. In addition, in the case of a major active region dominating the activity-induced signal, the ratio between the full width at half maximum (FWHM) and the RV peak-to-peak amplitudes of the cross correlation function can be used to infer the type of active region responsible for the signal for stars with \vsini$\le8$\kms. A ratio smaller than three implies a spot, while a larger ratio implies a plage. Using the observation of HD189733, we show that SOAP 2.0 manages to reproduce the activity variation as well as previous simulations when a spot is dominating the activity-induced variation. In addition, SOAP 2.0 also reproduces the activity variation induced by a plage on the slowly rotating star $\alpha$ Cen B, which is not possible using previous simulations. Following these results, SOAP 2.0 can be used to estimate the signal induced by spots and plages, but also to correct for it when a major active region is dominating the RV variation.Comment: 28 pages, 15 figures, accepted for publication in ApJ after minor revisions (taken into account in this version

On the long-term correlation between the flux in the Ca II H & K and Halpha lines for FGK stars

The re-emission in the cores of the Ca II H & K and H$\alpha$ lines, are well known proxies of stellar activity. However, these activity indices probe different activity phenomena, the first being more sensitive to plage variation, while the other one being more sensitive to filaments. In this paper we study the long-term correlation between $\log R'_{HK}$ and $\log I_{H\alpha}$, two indices based on the Ca II H & K and H$\alpha$ lines respectively, for a sample of 271 FGK stars using measurements obtained over a $\sim$9 year time span. Because stellar activity is one of the main obstacles to the detection of low-mass and long-period planets, understanding further this activity index correlation can give us some hints about the optimal target to focus on, and ways to correct for these activity effects. We found a great variety of long-term correlations between $\log R'_{HK}$ and $\log I_{H\alpha}$. Around 20% of our sample has strong positive correlation between the indices while about 3% show strong negative correlation. These fractions are compatible with those found for the case of early-M dwarfs. Stars exhibiting a positive correlation have a tendency to be more active when compared to the median of the sample, while stars showing a negative correlation are more present among higher metallicity stars. There is also a tendency for the positively correlated stars to be more present among the coolest stars, a result which is probably due to the activity level effect on the correlation. Activity level and metallicity seem therefore to be playing a role on the correlation between $\log R'_{HK}$ and $\log I_{H\alpha}$. Possible explanations based on the influence of filaments for the diversity in the correlations between these indices are discussed in this paper.Comment: 18 pages, 13 figures, 4 tables, accepted for publication in Astronomy and Astrophysic

Measuring precise radial velocities and cross-correlation function line-profile variations using a Skew Normal density

Context. Stellar activity is one of the primary limitations to the detection of low-mass exoplanets using the radial-velocity (RV) technique. Stellar activity can be probed by measuring time-dependent variations in the shape of the cross-correlation function (CCF). It is therefore critical to measure with high-precision these shape variations to decorrelate the signal of an exoplanet from spurious RV signals caused by stellar activity. Aims. We propose to estimate the variations in shape of the CCF by fitting a Skew Normal (SN) density which, unlike the commonly employed Normal density, includes a Skewness parameter to capture the asymmetry of the CCF induced by stellar activity and the convective blueshift. Methods. We compared the performances of the proposed method to the commonly employed Normal density using both simulations and real observations with different levels of activity and signal-to-noise ratios. Results. When considering real observations, the correlation between the RV and the asymmetry of the CCF and between the RV and the width of the CCF are stronger when using the parameters estimated with the SN density rather than those obtained with the commonly employed Normal density. In particular, the strongest correlations have been obtained when using the mean of the SN as an estimate for the RV. This suggests that the CCF parameters estimated using a SN density are more sensitive to stellar activity, which can be helpful when estimating stellar rotational periods and when characterizing stellar activity signals. Using the proposed SN approach, the uncertainties estimated on the RV defined as the median of the SN are on average 10% smaller than the uncertainties calculated on the mean of the Normal. The uncertainties estimated on the asymmetry parameter of the SN are on average 15% smaller than the uncertainties measured on the Bisector Inverse Slope Span (BIS SPAN), which is the commonly used parameter to evaluate the asymmetry of the CCF. We also propose a new model to account for stellar activity when fitting a planetary signal to RV data. Based on simple simulations, we were able to demonstrate that this new model improves the planetary detection limits by 12% compared to the model commonly used to account for stellar activity. Conclusions. The SN density is a better model than the Normal density for characterizing the CCF since the correlations used to probe stellar activity are stronger and the uncertainties of the RV estimate and the asymmetry of the CCF are both smaller.Peer reviewe

Noise Sources in Photometry and Radial Velocities

The quest for Earth-like, extrasolar planets (exoplanets), especially those located inside the habitable zone of their host stars, requires techniques sensitive enough to detect the faint signals produced by those planets. The radial velocity (RV) and photometric transit methods are the most widely used and also the most efficient methods for detecting and characterizing exoplanets. However, presence of astrophysical "noise" makes it difficult to detect and accurately characterize exoplanets. It is important to note that the amplitude of such astrophysical noise is larger than both the signal of Earth-like exoplanets and state-of-the-art instrumentation limit precision, making this a pressing topic that needs to be addressed. In this chapter, I present a general review of the main sources of noise in photometric and RV observations, namely, stellar oscillations, granulation, and magnetic activity. Moreover, for each noise source I discuss the techniques and observational strategies which allow us to mitigate their impact.Comment: 11 pages, 2 tables, Lecture presented at the IVth Azores International Advanced School in Space Sciences on "Asteroseismology and Exoplanets: Listening to the Stars and Searching for New Worlds" (arXiv:1709.00645), which took place in Horta, Azores Islands, Portugal in July 201

YARARA V2: reaching sub-m s−1 precision over a decade using PCA on line-by-line radial velocities

Context. The detection of Earth-like planets with the radial velocity (RV) method is extremely challenging today due to the presence of non-Doppler signatures such as stellar activity and instrumental signals that mimic and hide the signals of exoplanets. In a previous paper, we presented the YARARA pipeline, which implements corrections for telluric absorption, stellar activity, and instrumental systematics at the spectral level, and then it extracts line-by-line (LBL) RVs with a significantly better precision than standard pipelines. Aims. In this paper, we demonstrate that further gains in RV precision can be achieved by performing principal component analysis (PCA) decomposition on the LBL RVs. Methods. The mean-insensitive nature of PCA means that it is unaffected by true Doppler shifts, and thus can be used to isolate and correct nuisance signals other than planets. Results. We analysed the data of 20 intensively observed HARPS targets by applying our PCA approach on the LBL RVs obtained by YARARA. The first principal components show similarities across most of the stars and correspond to newly identified instrumental systematics for which we can now correct. For several targets, this results in an unprecedented RV root-mean-square of around 90 cm s−1 over the full lifetime of HARPS. We used the corrected RVs to confirm a previously published 120-day signal around 61 Vir, and to detect a super-Earth candidate (K ~ 60 ± 6 cm s−1, m sin i = 6.6 ± 0.7 M⊕) around the G6V star HD 20794, which spends part of its 600-day orbit within the habitable zone of the host star. Conclusions. This study highlights the potential of LBL PCA to identify and correct hitherto unknown, long-term instrumental effects and thereby extend the sensitivity of existing and future instruments towards the Earth analogue regime

Long-term magnetic activity of a sample of M-dwarf stars from the HARPS program II. Activity and radial velocity

Due to their low mass and luminosity, M dwarfs are ideal targets if one hopes to find low-mass planets similar to Earth by using the radial velocity (RV) method. However, stellar magnetic cycles could add noise or even mimic the RV signal of a long-period companion. Following our previous work that studied the correlation between activity cycles and long-term RV variations for K dwarfs we now expand that research to the lower-end of the main sequence. Our objective is to detect any correlations between long-term activity variations and the observed RV of a sample of M dwarfs. We used a sample of 27 M-dwarfs with a median observational timespan of 5.9 years. The cross-correlation function (CCF) with its parameters RV, bisector inverse slope (BIS), full-width-at-half- maximum (FWHM) and contrast have been computed from the HARPS spectrum. The activity index have been derived using the Na I D doublet. These parameters were compared with the activity level of the stars to search for correlations. We detected RV variations up to ~5 m/s that we can attribute to activity cycle effects. However, only 36% of the stars with long-term activity variability appear to have their RV affected by magnetic cycles, on the typical timescale of ~6 years. Therefore, we suggest a careful analysis of activity data when searching for extrasolar planets using long-timespan RV data.Comment: 20 pages, 12 figures, 3 tables, accepted for publication in Astronomy and Astophysic

YARARA: Significant improvement of RV precision through post-processing of spectral time-series

Aims: Even the most-precise radial-velocity instruments gather high-resolution spectra that present systematic errors that a data reduction pipeline cannot identify and correct for efficiently. In this paper, we aim at improving the radial-velocity precision of HARPS measurements by cleaning individual extracted spectra using the wealth of information contained in spectra time-series. Methods: We developed YARARA, a post-processing pipeline designed to clean high-resolution spectra from instrumental systematics and atmospheric contamination. Spectra are corrected for: tellurics, interference pattern, detector stitching, ghosts and fiber B contaminations as well as more advanced spectral line-by-line corrections. YARARA uses Principal Component Analysis on spectra time-series with prior information to disentangle contaminations from real Doppler shifts. We applied YARARA on three systems: HD10700, HD215152 and HD10180 and compared our results to the HARPS standard Data Reduction Software and the SERVAL post-processing pipeline. Results: On HD10700, we obtain radial-velocity measurements that present a rms smaller than 1 m/s over the 13 years of the HARPS observations, which is 20 and 10 % better than the HARPS Data Reduction Software and the SERVAL post-processing pipeline, respectively. We also injected simulated planets on the data of HD10700 and demonstrated that YARARA does not alter pure Doppler shifted signals. On HD215152, we demonstrated that the 1-year signal visible in the periodogram becomes marginal after processing with YARARA and that the signals of the known planets become more significant. Finally, on HD10180, the known six exoplanets are well recovered although different orbitals parameters and planetary masses are provided by the new reduced spectra.Comment: 23 pages, 19 figure

The HARPS search for southern extrasolar planets XXV. Results from the metal-poor sample

Searching for extrasolar planets around stars of different metallicity may provide strong constraints to the models of planet formation and evolution. In this paper we present the overall results of a HARPS (a high-precision spectrograph mostly dedicated to deriving precise radial velocities) program to search for planets orbiting a sample of 104 metal-poor stars (selected [Fe/H] below -0.5). Radial velocity time series of each star are presented and searched for signals using several statistical diagnostics. Stars with detected signals are presented, including 3 attributed to the presence of previously announced giant planets orbiting the stars HD171028, HD181720, and HD190984. Several binary stars and at least one case of a coherent signal caused by activity-related phenomena are presented. One very promising new, possible giant planet orbiting the star HD107094 is discussed, and the results are analyzed in light of the metallicity-giant planet correlation. We conclude that the frequency of giant planets orbiting metal-poor stars may be higher than previously thought, probably reflecting the higher precision of the HARPS survey. In the metallicity domain of our sample, we also find evidence that the frequency of planets is a steeply rising function of the stellar metal content, as found for higher metallicity stars.Comment: Accepted for publication in A&

The HARPS search for Earth-like planets in the habitable zone: I -- Very low-mass planets around HD20794, HD85512 and HD192310

In 2009 we started an intense radial-velocity monitoring of a few nearby, slowly-rotating and quiet solar-type stars within the dedicated HARPS-Upgrade GTO program. The goal of this campaign is to gather very-precise radial-velocity data with high cadence and continuity to detect tiny signatures of very-low-mass stars that are potentially present in the habitable zone of their parent stars. Ten stars were selected among the most stable stars of the original HARPS high-precision program that are uniformly spread in hour angle, such that three to four of them are observable at any time of the year. For each star we recorded 50 data points spread over the observing season. The data points consist of three nightly observations with a total integration time of 10 minutes each and are separated by two hours. This is an observational strategy adopted to minimize stellar pulsation and granulation noise. We present the first results of this ambitious program. The radial-velocity data and the orbital parameters of five new and one confirmed low-mass planets around the stars HD20794, HD85512, and HD192310 are reported and discussed, among which is a system of three super-Earths and one that harbors a 3.6 Earth-mass planet at the inner edge of the habitable zone. This result already confirms previous indications that low-mass planets seem to be very frequent around solar-type stars and that this may occur with a frequency higher than 30%Comment: 18 pages, 22 figures, accepted by A&A on 15/08/2011 with reference AA/2011/17055. Radial velocity data will be available through CD