Lineprofile variations in radialvelocity measurements
Two alternative indicators for planetary searches
^{1}
Centro de Astrofísica, Universidade do Porto, Rua das Estrelas, 4150762
Porto, Portugal
email:
pedro.figueira@astro.up.pt
^{2}
Departamento de Física e Astronomia, Faculdade de Ciências,
Universidade do Porto, Portugal
^{3}
Observatoire Astronomique de l’Université de Genève,
51 Ch. des Maillettes, 1290
Sauverny, Versoix,
Suisse
^{4}
Laboratoire Lagrange, UMR 7293, UNS/CNRS/OCA,
06300
Nice,
France
Received:
22
November
2012
Accepted:
8
July
2013
Aims. We introduce two methods to identify falsepositive planetary signals in the context of radialvelocity exoplanet searches. The first is the biGaussian crosscorrelation function fitting (and monitoring of the parameters derived from it), and the second is the measurement of asymmetry in radialvelocity spectral line information content, V_{asy}. We assess the usefulness of each of these methods by comparing their results with those delivered by current indicators.
Methods. We make a systematic analysis of the most used common line profile diagnosis, Bisector Inverse Slope and Velocity Span, along with the two proposed ones. We evaluate all these diagnosis methods following a set of welldefined common criteria and using both simulated and real data. We apply them to simulated crosscorrelation functions that are created with the program SOAP and which are affected by the presence of stellar spots. We consider different spot properties on stars with different rotation profiles and simulate observations as obtained with highresolution spectrographs. We then apply our methodology to real crosscorrelation functions, which are computed from HARPS spectra, for stars with a signal originating in activity (thus spots) and for those with a signal rooted on a planet.
Results. We demonstrate that the biGaussian method allows a more precise characterization of the deformation of line profiles than the standard bisector inverse slope. The calculation of the deformation indicator is simpler and its interpretation more straightforward. More importantly, its amplitude can be up to 30% larger than that of the bisector span, allowing the detection of smalleramplitude correlations with radialvelocity variations. However, a particular parametrization of the bisector inverse slope is shown to be more efficient on highsignaltonoise data than both the standard bisector and the biGaussian. The results of the V_{asy} method show that this indicator is more effective than any of the previous ones, being correlated with the radialvelocity with more significance for signals resulting from a line deformation. Moreover, it provides a qualitative advantage over the bisector, showing significant correlations with RV for active stars for which bisector analysis is inconclusive.
Conclusions. We show that the two indicators discussed here should be considered as standard tests to check for the planetary nature of a radialvelocity signal. We encourage the usage of different diagnosis as a way of characterizing the often elusive line profile deformations.
Key words: planetary systems / techniques: radial velocities / line: profiles / methods: data analysis
© ESO, 2013
1. Introduction
The detection of extrasolar planets relies heavily on the radial velocity (RV) method, which accounts for the vast majority of the planetary discoveries up to now. Yet, it has been known since the very first announcement by Mayor & Queloz (1995) that this approach is vulnerable to false positives. The potential parasitic signals are often rooted in stellar photospheric effects that deform spectral line profiles and introduce a RV signal which is not rooted on a Doppler shift of the spectral lines.
The relatively high frequency of false positives led to the establishment of studies of line profile properties for detected RV signals as a way of proving – or better said, disproving – their planetary nature. The first indepth study was done by Queloz et al. (2001), whom established the bisector and the bisector inverse slope (BIS) as the paradigm for its genre. The authors showed that there was a clear correlation between BIS and the measured RV for the case of a signal created by a photospheric spot. Numerous studies have made use of this principle (e.g. Martínez Fiorenzano et al. 2005; Niedzielski et al. 2009; Hébrard et al. 2010; Figueira et al. 2010a, just to cite a few) to reject the planetary nature of a signal, and planethunting surveys now routinely search for a correlation between RV and BIS variations as a way of pinpointing a parasite signal.
Several alternative methodologies exist for identifying these signals. For stabilized spectrographs, the monitoring of the fullwidth at halfmaximum (FWHM) of the crosscorrelation function (CCF) can reveal profile variations (e.g. Dumusque et al. 2011), and the V_{span} was introduced by Boisse et al. (2011) as an alternative to BIS to handle cases of low signaltonoise ratio (S/N) data. Moreover, several campaigns compared RV measurements made in different wavelengths to exclude planetary candidates, since a cold spot will induce a signal whose amplitude depends on the contrast relative to the surrounding stellar disk (e.g. Martín et al. 2006; Huélamo et al. 2008; Figueira et al. 2010b). When the range of possible applications is wider, an indicator is more used, and the bisector, making very few assumptions and being easy to calculate, is by far the most popular.
In a general sense, the efficiency of a diagnosis criteria in identifying RV signals created by line deformations depends mostly on two properties. First, the indicator should be proportional to the RV signal created (if possible linearly proportional, as it eases its identification and interpretation), and second this proportionality should be characterized by a proportionality constant as large as possible (when compared with the average error bars of the indicator).
Even though BIS is undoubtedly useful in many situations, studies such as those of Saar & Donahue (1997), Santos et al. (2003b), and Desort et al. (2007) pointed to an important limitation. For many cases of line profile deformations, characterized by certain stellar and spot properties, the BIS is intrinsically less sensitive to the presence of a deformation than the RV induced by the deformation itself. In other words, the proportionality constant between BIS and RV for these configurations has an absolute value ≪1. On top of this, it is well known that the relation between BIS and RV is not expected to be monotonic (see e.g. Boisse et al. 2011). While for many practical cases the correlation can be assumed to be a straight line, the real shape is similar to that of a tilted “8”; the proportionality constant between the two variables can even change from positive to negative or viceversa for some cases. These two properties can preclude the detection of a correlation between BIS and RV and render the diagnostic ineffective.
As we expand our searches, including more varied stellar hosts and pushing the precision of our surveys further to detect loweramplitude signals, there is an increasing demand for alternative and more effective ways of identifying falsepositive signals. There have been few attempts at finding alternative methods of characterizing line profile variations in a way useful for planetary RV searches (with the work of Dall et al. 2006, being a rare exception).
In this paper, we consider two different diagnosis methods and compare their performance with those of current lineprofile indicators for a wide variety of cases. We start by presenting the definition and properties of BIS and V_{span}, our comparison references, in Sect. 2. Then we describe the data sets used for this study and the principles of our comparison in Sect. 3. We employ the standard lineprofile indicators and evaluate their results in Sect. 4. We move on to apply the first of our two methods, the BiGaussian fitting, in Sect. 5 and the second, V_{asy}, in Sect. 6. We discuss the results in Sect. 7 and conclude on the usefulness of these new indicators and how they help us understand the impact of lineprofile variation on RV in Sect. 8.
2. Current diagnosis methods
The spectral information content of a stellar spectrum can be condensed in a CCF, which is then used to measure precise RVs (Baranne et al. 1996; Pepe et al. 2002). The RV is extracted by fitting a Gaussian function to the CCF with the center of the Gaussian delivering the star’s RV; this quantity is represented throughout the paper as RV_{G, center}. We apply our methods to the CCF and its derived quantities and compare the results with those obtained for known lineprofile indicator variations. For reference lineprofile indicators, we consider the BIS, the most common diagnosis method for the presence of parasitic signals, and V_{span}. In the next sections we describe these methods in detail and how they were implemented.
2.1. The BIS
In the BIS analysis, we used the definition provided by Queloz et al. (2001). Consider a CCF and calculate its midpoint (i.e., the point at equal distance between its blue and red wing) for each flux level between the continuum level and the bottom of the line. Since the flux level will only coincidentally fall on a pixel value, the flux values are calculated by interpolation from the neighboring CCF pixels. The BIS value is then simply the difference between the velocity of the top and bottom sections of the line, with the top being the average midpoint of the flux levels located between 60 and 90% and the bottom of those between 10 and 40% of the line full depth. We refer the reader to the Fig. 5 of Queloz et al. (2001) for a helpful schematic representation. Since we are working in the velocity space, the line mid points are measured in km s^{1} (or any multiple of it), as is BIS.
The interpretation of the BIS is quite straightforward. Any RV signal associated with a variation between the top and bottom sections of the line is not rooted in the coordinated movement of it but rather on a line deformation, and a plot of BISRV will reveal a correlation.
Even though the vast majority of works follow this prescription (see references in the introduction), the top and bottom regions of the BIS can be defined differently. In particular, one can consider narrower bands located at different heights, leading to different slopes when the line is subjected to deformations. This different leverage capability will lead to different indicator amplitude variations at the expense of using only a smaller section of the line, and thus working at a lower S/N when applying it to real spectra. To explore this posibility, we considered two different parametrizations, representative of extreme cases: a maximum leverage case BIS^{+} in which the limits of top and bottom regions are 80−90% and 10−20% of the total flux level, respectively, and a minimum leverage case BIS^{−} with top at 60−70% and bottom at 30−40%, respectively.
2.2. The V_{span}
We employ the V_{span} indicator as defined in Boisse et al. (2011). This indicator represents the difference RV_{high} − RV_{low} between two Gaussian fits of the CCF which consider exclusively the upper and lower part of the CCF, respectively. For the first, RV_{high}, one takes into account only the top of the CCF for the fit by considering only the points in the range [−∞:−σ][+ σ:+ ∞] as measured relative to the center of the Gaussian fit; for the second RV_{low} one considers the points in the range [∞:−3σ]∨[σ:+ σ]∨[+ 3σ:+ ∞]. The former quantity is sensitive to variations on the top of the line and the latter on the bottom.
The construction of the V_{span} indicator was motivated by the analysis of lineprofile variations for cases of low S/N and was shown to be advantageous relative to the BIS for one of these cases. We refer the interested reader to Boisse et al. (2011) for details.
3. Choice of datasets and comparison principles
3.1. The datasets
We simulated CCFs as affected by the presence of a stellar spot using the program SOAP, and refer the interested reader to Boisse et al. (2012) for details on the program. All simulations share two basic assumptions: we acquire our spectra with a spectrograph with a resolution of 100 000 (which defines the FWHM of a nonrotating star), and we observe a star with a single dark spot (i.e., with zero emissivity). We consider different v sini for the star, different inclinations of the stellar rotation axis relative to the line of sight, different spot latitudes, and different spot filling factors. When performing our simulations, we use a grid resolution of 2000 and a spot circumference resolution of 50; the orbital phase is homogeneously sampled by 100 points. For the steps in RV for the CCF calculation, we considered a very fine grid with steps of 10 m/s from −30 km s^{1} to +30 km s^{1}. On real data the CCF is calculated on a much coarser grid, usually of ~100 m/s or higher. However, a coarse discretization can introduce an appreciable deviation in the fit and on the fit parameters. Since these synthetic data have, by definition, no noise, any fit instability steming from the small number of points can introduce a detectable difference relative to high S/N spectra with a coarser grid. Since computational power is not an issue for such a small number of simulations, we choose to calculate the CCF on a very fine grid to evaluate the application of our indicators in the most precise way possible.
We started by simulating the properties of standard wellstudied cases, with the first being that of HD 166435, already discussed both in Queloz et al. (2001) and in Desort et al. (2007). The star has a vsini = 7.0 km s^{1} and is observed at an inclination of 30°; its spot has a filling factor of 1% and is located at a latitude of 60°. This is the most extensively studied spot configuration in the literature, and we included it here for that reason^{1}.
We then moved to more general cases with the objective of studying the impact of key stellar and spot parameters on the amplitude of the measured RVs and the associated indicator’s response. First we studied the impact of equatorial spots (latitude of 0°) with different filling factor values of 1.0, 2.0, 5.0, and 10.0% on a edgeon star (I = 0°) with vsini = 7.0 km s^{1}. Then, we repeated the exercise for a 1% filling factor equatorial spot on a edgeon star with rotational velocities of 2, 5, and 10 km s^{1}. These tests were also inspired from those of Desort et al. (2007), who chose these inclination and spot latitudes to induce the highest RV amplitude for a given filling factor and vsini. The spot filling factors were chosen to depict typical spot sizes on mainsequence stars, and the vsini were chosen to illustrate several cases ranging from a broadening inferior to that introduced by the resolution of the spectrograph, to a star rotating so fast that the RV precision is significantly degradated.
We finally tested our approach on real data, and to do so we used HARPS spectra and CCFs as delivered by the last version of the DRS pipeline (Mayor et al. 2003). We selected both active stars (log or higher), for which a correlation between RV and the indicators (ΔV and BIS) is expected, and a planethost star HD 216770, for which it is not (or is not expected at the level of the measured RV signal amplitude). The planetary orbits measured on this star were discussed in Mayor et al. (2004), and the reader is referred to the paper for details. The stars are listed in Table 1, along with the log (), vsini, and the number of data points taken at the time of the publication and considered for the analysis. The v sini was derived in Santos et al. (2003a) and log () was derived in Santos et al. (2000). The only exception was HD 216770, of which log () was taken from Henry et al. (1996).
HARPS star properties.
When the activity of the star is larger, the impact on both measured RV and lineprofile variation indicators is larger. To maximize its impact on RV, one should measure it at least over one full activity cycle. In practice, this is very difficult, and recent works have shown that several longerterm activity cycles are present in a star on top of the typical rotationmodulated cycle, with a strong impact on RV (e.g. Dumusque et al. 2011; Lovis et al. 2011). While keeping in mind that the characterization of the activityinduced effects is still work in progress, we concentrate our analysis on the relative variations between RV and lineprofile indicators and how the different indicators compare with each other.
Comparison between correlations of BIS and V_{span} with RV for simulated data.
Comparison between BIS and V_{span} correlations with RV for HARPS data for the parameters that do not take into consideration the error bars.
Comparison between BIS and V_{span} correlations with RV for HARPS data for the parameters that consider the error bars, i.e. the slope fit parameters.
3.2. Comparing the datasets
We aim to compare the different indicators in the most systematic way possible. To do so, we apply them to the same datasets and evaluate them according to the same criteria.
We first compare the amplitude of the indicator’s variation (relative to its average value) with those of the RV and other indicators. We pay particular attention to how the indicators perform compared to BIS, the most commonly used lineprofile indicator. If one assumes that the error bars are similar, a larger amplitude of the indicator’s variation results in a greater capability to pinpoint a signal. Yet, the indicator’s capabilities must be assessed through a rigorous quantitative analysis. To compare the correlation of an indicator with RV, we also consider:

the slope resulting from the fit of a linear function by leastsquare minimization (m_{ind});

the Pearson’s correlation coefficient value for the data sets (ρ_{ind}); and

the probability of obtaining an equal or larger correlation coefficient absolute value by chance.
The slope of the correlation between an indicator and the RV is a different way of relating a variation in the indicator with that of the RV. However, it does not evaluate the strength of the correlation betwen the two data sets. To do this, we use the Pearson’s correlation coefficient, which is a measure of linear dependence between the two variables. To test if the correlation coefficient’s value is significant, we test if it can be obtained by a fortuitous pairmatching between the indicator and the RV. To do so, we bootstrap 100 000 times the data pairs and calculate the correlation coefficient ρ_{ind} for each of these synthetic uncorrelated data sets. We then calculate the average and standard deviation (σ_{ρ,ind}) of these ρ_{ind} values; by assuming a Gaussian distribution, one can calculate the probability that the ρ_{ind} of the original dataset was obtained by pure chance. While the probability value is conditioned by the hypothesis that the indicators follow a Gaussian distribution, its distance in σ is still a good indicator on how the measured value compares with that of an uncorrelated distribution.
Comparison between correlations of BIS^{+} and BIS^{−} with RV for simulated data.
Comparison between BIS^{+} and BIS^{−} correlations with RV for HARPS data for the parameters that do not consider the error bars.
Comparison between BIS^{+} and BIS^{−} correlations with RV for HARPS data for the parameters that consider the error bars, i.e. the slope fit parameters.
When using real data, we have two important additional information elements. The first of these is the RV measurement uncertainty, as estimated from the RV information content of the spectra (e.g. Bouchy et al. 2001). Moreover, we know if the RV signal introduced is of planetary or stellar origin and, as a consequence, if a line profile deformation is present or not; we also know then if the indicator is expected to be correlated with the measured RV variation or not. It follows then that we can add the following criteria when evaluating the real data sets:

the error bars on the correlation slope (m)_{ind} which result from bootstrapping the RV data and indicator with respect to the respective error bars, performed by drawing 100 000 times values around their error bars;

the relationship of the slope and correlation values to the stellar activity or presence of planets around the star.
We note that the bootstrapping is a nonparametric test, depending on very few assumptions, and thus being very robust. However, to take full advantage of the real data, one has to estimate the error bars on the indicators themselves. In spite of the abundant literature on BIS, there are virtually no studies on the error distribution of this estimator. It is commonly reported that the error bar on this difference is on the order of times the photon noise on the CCF because BIS corresponds to a difference of velocities between the upper and the lower sections of the CCF. Even though these error bars are illjustified, our objective is to compare the performance among indicators, and and so we choose to use the commonly accepted uncertainties for BIS for the other indicators. To allow a more straightforward interpretation of our results, we separate those that depend on the assumptions on the error bars from those which do not and present them in different tables.
4. Standard profile indicators: BIS and V_{span}
4.1. BIS vs. V_{span}
We start by applying the two standard lineprofile indicators to our datasets. These have been discussed in their presentation papers and an exhaustive analysis is not the objective of this paper, so we focus instead on their relative strenghts and weaknesses.
Both the analysis of synthetic (Table 2) and real data (Tables 3 and 4) show that the BIS and V_{span} yields similar results for most of the cases. Their variation is proportional to both the spot filling factor and the star rotational velocity. For most of the considered cases, the results are of the same order of magnitude for the data and the different evaluation criteria considered. BIS exhibits higheramplitude variations and more significant correlations in a consistent way for all synthetic datasets but does not prove to be advantageous when applied to real data. V_{span} yielded poorer results than those delivered by the BIS on what concerns correlation slopes and associated errors, but we must stress that the latter depends on the assumed error bars for the bootstrapping.
4.2. Different parametrizations of BIS
The first aspect coming out of Table 5 is that the variation of BIS^{+} is larger than that of BIS^{−} by a factor of 2−3 for all the cases presented. This is associated with stronger correlations and correlation slopes for the former indicator. Moreover, one concludes that the BIS^{+} leads to a variation that is larger than that of BIS by ~50% for all the configurations when compared with the corresponding values of the tables in the previous section; the correlations slopes are larger, with the correlation coefficient being very similar.
The analysis of the results on real data, presented in Tables 6 and 7, are not easy to interpret. The amplitude variations are larger for BIS^{+} than for BIS^{−}; the correlation coefficients and correlation slopes are similar or larger for BIS^{+} than for BIS^{−} with one clear exception, BD213253. For the latter, the BIS^{−} correlation coefficient even approaches the 10% probability level, hinting at the existence of a correlation. This is very interesting in itself, because neither the BIS nor the V_{span} alone could reach this level. However, the high probability of obtaining these correlations by chance prevent us from taking this discussion further.
Comparison between correlations of BIS and ΔV with RV for simulated data.
It is also interesting to compare the results of BIS with those of BIS^{+} for the real data. The amplitude variation of BIS^{+} can be larger by a factor of more than 2 than that of the BIS for some cases, but they deliver essentially the same results when it comes to the detection of correlations. Interestingly, an increase in amplitude can even be accompanied by a decrease in the correlation factor; the results on the planethost star are an interesting example of this. The same happened on BD213253, decreasing the correlation coefficient to completely negligible levels. This behavior underlines the fact that the variation in amplitude cannot be interpreted straightforwardly. Moreover, it is important to note that that the spectra considered here are high S/N spectra and that the CCF has then a very high S/N; BIS^{+} is expected to be much less advantageous when the S/N decreases.
The main conclusion from the analysis of the BIS^{+} is that it tends to be advantageous for the cases considered here, but it presents a complex behavior, ranging from clearly advantageous to slightly disadvantageous, and it did not allow to detect any new correlations in the data.
5. The biGaussian fitting
5.1. The method
The biGaussian is a simple function that considers a Gaussian with wings characterized by two different HWHM (halfwidth at halfmaximum), or, equivalently, an asymmetry parameter associated with the FWHM. It was first used in the context of RV measurements by Nardetto et al. (2006) as a way of simultaneously measuring the asymmetry and the center of spectral lines on pulsating stars. It can be represented by (1)and (2)in which D is the depth of the line, and C the continuum level, as measured in photoelectrons^{2}. Since we apply it to the CCF, we work in the RV space, and both the biGaussian center RV_{center} and FWHM are expressed in km s^{1}. The asymmetry A is given as a percentage (of the FWHM). For illustration, we plot the effect of introducing an asymmetry on a Gaussian function in Fig. 1.
Fig. 1 Comparison of biGaussian functions with different asymmetry values: 0% (Gaussian function – solid line), 10% (dashed line), 25% (dasheddotted line), and 50% (dotted line). These normalized functions have a contrast of 30%, a FWHM of 5 km s^{1}, and are all centered at 0 m/s. 

Open with DEXTER 
This function is fitted by minimizing the , just like the Gaussian function; by imposing A = 0 one recovers the typical Gaussian fitting. We denote the RV_{center} parameter obtained by Gaussian fitting by RV_{G,center}. As stated by Nardetto et al. (2006), the indicator ΔV = RV_{G,center} − RV_{center} represents the RV shift that can be explained by the asymmetry alone, and we consider it as our proxy for the line deformation. The interpretation is once again very simple. The extra degree of freedom in the biGaussian fitting allows to accommodate a line profile variation and is thus more insensitive to the line deformations than a simple Gaussian fitting. The operational similarities between BIS and ΔV then follow very naturally. In this work we considered ΔV = RV_{center} − RV_{G,center}, so that a plot ΔVRV displays a positive correlation.
5.2. ΔV: application to simulated data
The results for the simulated data are presented in Table 8. In Fig. 2 we plot the phase dependence of the most important parameters for the case of the 1% filling factor and BIS and ΔV as a function of RV_{G,center} in Fig. 3.
Fig. 2 Most important properties for the two fitted functions: Gaussian as open circles (only two upper panels, partially superposed) and BiGaussian as filled circles, as a function of orbital phase (with an arbitrary zeropoint) for the case of an equatorial spot with a a filling factor of 1% on an edgeon star with v sini of 7 km s^{1}. Upper left: center of the two fitted functions as a function of orbital phase; Upper right: FWHM of the two functions; on the two lower panels, we have the asymmetry A (bottom left) and the ΔV indicator (bottom right) for the BiGaussian fitting. 

Open with DEXTER 
It is easy to notice that it is advantageous to use ΔV relative to the BIS as an indicator of lineprofile deformations for the simulated data. The amplitude of its variation is always higher with the improvement being of 10−30% (with the most significant improvement being obtained for the largest projected rotational velocity). The value of the correlation slope is also larger, with a variable improvement, depending on the case considered, but reaches above 50%. These two factors suggest that the ΔVRV correlation and its slope will emerge more easily in cases where the BISRV slope is small when compared with the photon noise level. The absolute value of Pearson’s correlation coefficient is very high and very similar for both correlations, and the probability that these happen by chance is smaller than 0.01%. In a nutshell, we can note that the application of the indicator is advantageous when compared with that of BIS for the whole range of spot filling factors and vsini considered in the synthetic data, since the amplitude variations are higher with all other parameters being similar. From what we discussed in Sect. 4.2 we note that BIS^{+} variations are higher than those of ΔV by ~30%, while displaying similar correlation coefficients.
Fig. 3 Two indicators considered, BIS (as open circles) and ΔV (as filled circles), as a function of RV, for the case of an equatorial spot with a filling factor of 1% on an edgeon star with vsini of 7 km s^{1}. 

Open with DEXTER 
5.3. ΔV: Application to HARPS data
At this point, it is important to recall the discussion on the uncertainties of the estimators, and, in particular, of BIS, as presented in Sect. 3.2. Yet, the error bars on ΔV are a different issue: these result from the subtraction of the RV’s obtained by fitting the same CCF with two different functions. If one assumes that the fitting error is negligible and the final error is dominated by the photonnoise contribution as for BIS when fitting a function to the data, the error bars on the difference are the same as those of either fitting. In other words, the error of both fitting procedures is expected to be completely correlated. To test the two approaches on an even footing, we consider however the error bars on ΔV as being also times the photon noise, just like for BIS.
As done before, we separate the results into two tables: The first contains the quantities that do not depend on the error values and their assumptions, as seen in Table 9. The second lists the fitted correlation slopes and uncertainties drawn from bootstrapping the indicators, which depend on the indicator’s errors, as seen in Table 10. The reader thus should keep in mind the assumptions behind the results of this last Table.
Comparison between BIS and ΔV correlations with RV for HARPS data for the parameters that do not consider the error bars.
Comparison between BIS and ΔV correlations with RV for HARPS data for the parameters that consider the error bars, i.e. the slope fit parameters.
When applying ΔV to real data, the situation is slightly different than that of the application to simulated data. The amplitude of ΔV and the absolute value of Pearsons’s correlation coefficient are higher than those of BIS, but for both planethosting and nonplanet hosting stars. Consequently, the probability that the correlation happens by chance tends to be higher with the exception on both of these properties being that of BD213253. In particular, a more pronounced correlation is obtained for the planethost star HD 216770. As expected, the statistical test however shows that one has a nonnegligible probability of attaining the measured ρ by chance, which shows that the correlation is not to be trusted, just as for the BIS. We plot the application of the indicator to the star with the most significant correlation, HD 224789, and compare it to BIS in Fig. 4. Analyzing the error bars on the slope m is more difficult, and the results are more heterogeneous. As discussed before, these results stand on additional assumptions. The first point to note is that the relative error bars (in %) for the ΔV are always smaller than those for the BIS, except for the case of the star for which no correlation was detected, BD213253. As for the others, the absolute values for m are always larger, as expected, but the results are statistically very similar. We also note that all the slopes are noncompatible with zero at several sigma, if the error bars are trusted, even those from BD213253 and our planethost, HD 216770. As a concluding remark on the topic, we note that the results from this indicator are very similar to those obtained using BIS^{+} on the same stars.
Fig. 4 Dependence between the indicator ΔV (filled circles) and BIS (open circles) with RV for HD 224789. 

Open with DEXTER 
In a nutshell, the ΔV indicator proves to have a quantitative advantage over the BIS in the sense that the same results are obtained with higher significance, but no qualitative difference in the sense that the same conclusions on each signal’s nature were reached.
The ΔV, in spite of a very different construction from that of BIS, shares important traits with this other indicator; in particular, both display an 8’like shape in their correlation with RV (see Fig. 3). This suggested that we might be in presence of a common property of profile analysis indicators that have been created so far and it encouraged us to devise a completely different approach, as described in the next section.
Comparison between correlations of BIS and V_{asy} with RV for simulated data.
6. The analysis of the asymmetry in the RV calculation: V_{asy}
6.1. The V_{asy} indicator
We consider now an indicator that evaluates the RV information content of a spectral line and compares the information on the blue wing (i.e., for wavelengths shorter than that of the line center) to that on the red wing (i.e., for wavelengths longer than the center). For every flux level, the spectral information difference is calculated, and the results for all fluxes are condensed onto a weighted average value for the line in which the weight is the weight attributed to the flux level. The average spectral information between the red and blue wing is used as a proxy for the weight of the undeformed line, i.e., as a weight for the flux level. We apply this procedure to the CCF and mathematically define this indicator, V_{asy}, as (3)in which W_{i} are the weights of the point calculated at the flux level i, as defined in Eq. (8) of Bouchy et al. (2001), where (4)with λ(i) and A_{0}(i) representing the wavelength and flux of the spectra, respectively, and the detector readout noise. Since we are working in the velocity space, the wavelength is a velocity. This is standard procedure and has been applied extensively to the CCF. Moreover we are working in a very high S/N domain with the CCF and as such, . We can then rewrite the previous equation in the form of (5)in which RV(i) is the RV value of each pixel (i.e. xaxis) of the CCF (not to be mistaken by the center of the CCF obtained by Gaussian fitting, RV_{G,center}).
For every flux level i, we consider three quantities: the blue wing information, W_{i}(blue), the red wing information, W_{i}(red), and the average between the two, , as representative of the weight of the information content at this flux level for an undeformed line.
We calculate the flux levels i by quadratic interpolation between the CCF flux levels. We considered 100 equidistant flux levels and excluded the first and last 5 from the sum, to remove the points closer to the continuum and to the bottom of the line, where the spectral information is small compared with the noise.
In a nutshell, this indicator, V_{asy}, is simply the (dimensionless) average information content difference between the left and the right wing of the line. It is much closer to the RV calculation principle than line profile indicators such as BIS or ΔV. Just as for the previously discussed indicators, our indicator’s variation is expected to be correlated with that of RV for the case of a signal rooted in a line deformation.
6.2. V_{asy}: application to simulated data
The results of the application of this indicator to the simulated data are presented in Table 11. In Fig. 5, we plot the variation of RV and V_{asy} as a function of phase for the reference case of an equatorial spot with a filling factor of 1% on a edgeon star with a v sini of 7 km s^{1}, and in Fig. 6, we plot how the correlation between BIS and RV compares to that between V_{asy} and RV.
Fig. 5 RV (top) and V_{asy} (bottom) as a function of orbital phase (with an arbitrary zero point) for the case of an equatorial spot with a a filling factor of 1% on an edgeon star with vsini of 7 km s^{1}. 

Open with DEXTER 
Fig. 6 Two indicators considered, BIS (as open circles) and V_{asy} (as filled circles), as a function of RV for the case of an equatorial spot with a a filling factor of 1% on an edgeon star with vsini of 7 km s^{1}. 

Open with DEXTER 
Comparison between correlations of BIS and V_{asy} with RV for HARPS data.
The amplitude variation of the indicator and how it compares with that of BIS are now meaningless, since these indicators are not calculated in the same unit. Of greater importance and insight is that the Pearson’s correlation coefficients ρ for V_{asy} are always larger than those for the BIS or BIS^{+}; the improvement is often in the range of 10−20%. It then follows naturally that the offset in σ relative to an uncorrelated datasets are very high and that the probabilities of obtaining these correlation coefficient values by chance are very low (even though one has to stress that they are negligibly low for all cases considered, as before).
An inspection of Fig. 6 shows that our objective was accomplished in the sense that the relation between V_{asy} and RV is much closer to a linear one than that with the BIS. There is still a departure from linearity for small RV variations, but the usefulness of the method is clear.
6.3. V_{asy}: application to HARPS data
In Sect. 5.3, we applied ΔV to real data and we now repeat the same procedure and analysis for V_{asy}. The only difference is that we refrain from attempting to calculate the error bars on the indicator. As we saw in the mentioned section, this was the weakest point of the analysis, requiring more assumptions and making the interpretation of results more difficult. As a consequence, we also chose not to calculate the error bars on the correlation slope m. The results are presented in Table 12 and in Fig. 7, we depict the correlation between V_{asy} and RV and how it compares with that of BIS for the case of HD 224789.
Fig. 7 Dependence between the indicator V_{asy} (filled circles) and BIS (open circles) with RV for HD 224789. 

Open with DEXTER 
As before, we skip the analysis of the variation of the coefficient, and we note that while the values of the correlation slopes are higher this does not mean the method is more efficient in itself. These properties are in line with what was measured for the synthetic data, but they do not allow us to reach any conclusion in themselves. However, it is interesting to note that for 4 of the 5 active stars the Pearson’s correlation coefficient was significantly improved, in one case the improvement was by a factor of two, and in the case of the planethosting star, it decreased. Consequently, the deviation in σ increased for the four active stars and the probability of the correlation arising from a fortuitous alignment of the data decreased. The opposite was obtained for the planethost star. This provides a strong argument for the usefulness of the indicator and illustrates very well how advantageous it can be when compared with the BIS.
7. Discussion
We have shown that the ΔV indicator provides a more efficient characterization of line deformations than the BIS, and very similar to that of BIS^{+}. This property can make the difference when the correlation slope is comparable to the photon noise and/or the number of observations is very small. The indicator proved to be advantageous over the standard BIS for all the simulated datasets and slightly so for the real ones for which a correlation had been found between the RV and the BIS. In this sense, we can then conclude that the ΔV provides a quantitative advantage over the BIS. Yet, it is important to note that some of the active stars show RV variations without being correlated with BIS or BIS^{+}, and our new indicator did not allow to identify new correlations (for the meager number of stars tested).
Furthermore, the error bars on ΔV and BIS can be used to assess the significance of the results in another way. If one assumes that the activity signal is the only signal present in the RV data, the residuals should follow the noise distribution of the indicator once the slope is subtracted. This is expected to be a Gaussian distribution with a center at zero and a scatter dictated by the photon noise or, to be more generic, proportional to it. We consider different multiplicative factors α for the relation between the observed scatter and the one dictated by photon noise σ_{obs} = α × σ_{ph}. We evaluate them through the KolmogorovSmirnoff analysis (as implemented in Press et al. 1992), and evaluate which of these postulated distributions G(α σ_{ph}) has the highest probability of coming from the same parent distribution as the observed one, and thus which α is more appropriate to characterize σ_{obs}. We considered α over an interval ranging from 0 to 30 with a step of 0.5. We list the highestprobability α obtained for both correlations in Table 13.
The α correction factor on photon noise used to obtain the observed noise distribution.
We can see from the Table that neither of the methods/correlations is intrinsically superior when it comes to noise distribution. In other words, the correction factors applied to the error bars on the indicators are larger for the BISRV correlation for some of the cases, and larger for the ΔVRV correlations for some other cases. They are significantly higher than unity in both situations, which should send a warning to all those who use the BIS indicator routinely. Yet, one has to bear in mind that the fit of a straight line on the data is illjustified, since we know that the real shape is closer to an 8 rather than to a straight line. While the results still hold for a relative comparison, an α larger than unity was always expected.
We stress again that nonparametric tests are particularly important when the properties of the distributions are unknown or illdefined, and we recall the usefulness of the permutation test to the reader as a way of assessing the significance of a correlation.
As discussed by Nardetto et al. (2006), the fitting makes use of the full information of the line and will still deliver precise results at low S/N, which is a level at which BIS is significantly more vulnerable. The reason for this is twofold: BIS requires interpolation to calculate the flux levels and does not make use of all the spectral information content of the line. It is calculated only over a restricted range of flux. This is particularly important for those who use CCFs constructed from a relatively low number of lines or aim at measuring RV from lowS/N spectra. Neither is common for typical RV planet searches, but the latter is frequent when using RV measurements to confirm transit candidates from space observatories, such as CoRot (Csizmadia et al. 2011) or Kepler (Endl et al. 2011; Santerne et al. 2012). Interestingly, a significant source of false positives in transit searches is the blending of the spectral lines with those of a foreground or background star. The usage of indicators, such as the BIS or those presented in this paper, can again be of help, but this is a complex issue and beyond the reach of this paper. We will discuss it in a separate paper, using meaningful false positive examples from transit searches (Santerne et al., in prep.).
The results on V_{asy} are even more promising than those of ΔV, with the indicator being constructed in a completely different way. The results point for more secure correlations between V_{asy} and the RV for most of the analyzed active stars (4 out of 5) compared to those obtained using the BIS. It provides both a qualitative and quantitative advantage over the bisector in the sense that new correlations arise for active stars, and those previously identified are detected with higher significance. Yet, it has to be noted that even though the correlation is much closer to a linear one than that of the BIS, a departure from the linear relationship is apparent for the smallest amplitude RV variations. This attests for an undersensitivity of the indicator in this regime and warrants further investigation.
In this work we made a distinction and a separated analysis of synthetic spectra, without noise, and real spectra, with noise. This was a choice prompted by the availability of the SOAP model, which delivers noisefree CCFs, but allowed us to explore and distinguish expected theoretical results and observed ones. While one could have attempted to introduce noise in the CCFs and compare it with the real data, this would have to be done before the CCF computation, i.e., at the spectral generation level, using a completely different and by far more complicated procedure, which is highly dependent on the particular assumptions on the spectra. By construction, the calculation of the CCF transforms white noise into correlated noise, simply because the correlation step for the calculation is smaller than the pixel size. Moreover, a real spectra is expected to have its own noise correlated from one pixel to the other; this is a consequence of atmospheric extinction being a function of wavelength, just like instrumental transmission and instrumental profile. As a consequence, any model of the spectra which does not fully characterize the spectrograph, would fall short in describing its correlated noise, and thus lead to an incorrect evaluation of its impact on the RV and associated indicator’s calculation.
It has been suggested recently that there is a significant difference between the results derived from a simple spot model (such as SOAP) and comprehensive RV campaigns (e.g. Crockett et al. 2012). In this work, the parametrized synthetic cases show similar orderofmagnitude variations in RV and indicators compared to the real ones. Yet, one cannot exclude that there might be physical mechanisms that are not represented and might produce an appreciable effect. Works such as those of Meunier et al. (2010) showed the importance of considering the reduction of convective blueshift at the spot’s location on the measurement of precise stellar RV. This is only one of the many points not at all adressed by SOAP, since it emulates a very simple stellar model, and does not even consider emission of any kind at the spot’s location, nor the nuances in its emitted spectra. While a coldspot model might be indeed an oversimplification, we note that the methods discussed here proved to be effective not only for simulated data but for real data too.
It is also important to stress that the RV center obtained for a deformed line depends on the RV calculation method. For the implications of this characteristic on the case of RossiterMacLaughlin measurements, the interested reader is referred to Boué et al. (2013). This is a point to keep in mind when comparing the results obtained here with those obtained for RVs acquired with the iodine cell technique (e.g. Buchhave et al. 2011).
We conclude this discussion on a different note: instead of considering these different methods as alternative options, one should look at them as independent ways of characterizing the same phenomena and use them simultaneously, while bearing in mind their relative strengths and weaknesses. The active star BD123153 illustrates this very well: it exhibits a nonsignificant correlation between BIS and RV or ΔV and RV but shows a more significant one between V_{asy} and RV or BIS^{−} and RV. The use of different tests is then very important, especially because while the presence of a correlation disproves the planetary nature of a signal, its absence does not prove its planetary nature. We tried to do a systematic assessment of the potential of the presented methods, but the work presented here should be seen merely as the first step toward a more detailed characterization of these diagnosis methods and the development of complementary approaches.
8. Conclusion
In this paper, we introduced two diagnosis methods for detecting lineprofile variations in the context of exoplanet searches. The objective is to evaluate if these variations can be pinpointed as being at the root of detected RV signals. We attempted at a systematic assessment of the efficiency and sensitivity of our presented methods when compared with that of BIS, assuming different parametrizations, and V_{span}, and following very clear criteria. Our conclusions can be outlined as follows:

The standard BIS delivers very similar results to those of V_{span} for the cases considered. Different parametrizations of BIS can be explored to increase the leverage of the diagnosis at the expense of considering a lower fraction of the line and thus working at lower S/N. Yet, statisticaly more significant results can be obtained, showing that this is an option to consider.

The biGauss/ΔV method is superior to the standard BIS in a multitude of cases. It presents no qualitative advantage, identifying the same signals than BIS as activityrooted, but provides more significant results with the characterization parameters being affected by smaller error bars. This point can make the difference for cases where the correlation slope between the indicator and RV is close to the noise level and puts the indicator in even footing with the more effective BIS^{+} parametrization.

The V_{asy} indicator evaluates the unbalance between the red and blue wings of a line in terms of spectral information, being by construction very different from the other profile indicators. It shows a significant advantage over different parametrizations of BIS and ΔV, when applied to simulated spectra and to real spectra. The indicator leads to more significant correlations, and identifies new signals rooted in activity, providing a clear qualitative advantage. It also gives rise to a correlation with RV which is much closer to a straight line than any of the other profilebased indicators.
Finally we argue for the usage of independent techniques to characterize stellar line profile variations and to identify their origins, creating a coherent view of the stellar phenomena. This is particularly important in cases where the signal and/or potential correlation slopes approach the RV uncertainty level. A more detailed knowledge of the diagnosis methods themselves is required to meet the exigent requirements of today’s scientific programs. The work presented here should be considered as a first step in this direction; a more significant amount of resources should be put into characterizing these indicators and devising new and independent ones.
We note that the spot latitude definition used in SOAP is not the same as that used by Desort et al. (2007); in the former, the latitude angle is used, while the colatitude angle is used in the latter.
It is equivalent to work with a normalized function, as Nardetto et al. (2006) chose to do.
Acknowledgments
This work was supported by the European Research Council/European Community under the FP7 through Starting Grant agreement number 239953, and by Fundação para a Ciência e a Tecnologia (FCT) in the form of grant reference PTDC/CTEAST/098528/2008 and PTDC/CTEAST/098604/2008. N.C.S. would further like to thank FCT through program Ciência 2007 funded by FCT/MCTES (Portugal) and POPH/FSE (EC). We thank the anonymous referee for insightful comments that led to a significant improvement of the overall paper’s quality. P.F. would also like to thank Alexandre Santerne and Mahmoudreza Oshagh for enlightening discussions.
References
 Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Boisse, I., Bouchy, F., Hébrard, G., et al. 2011, A&A, 528, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Boisse, I., Bonfils, X., & Santos, N. C. 2012, A&A, 545, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bouchy, F., Pepe, F., & Queloz, D. 2001, A&A, 374, 733 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Boué, G., Montalto, M., Boisse, I., Oshagh, M., & Santos, N. C. 2013, A&A, 550, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buchhave, L. A., Bakos, G. Á., Hartman, J. D., et al. 2011, ApJ, 733, 116 [NASA ADS] [CrossRef] [Google Scholar]
 Crockett, C. J., Mahmud, N. I., Prato, L., et al. 2012, ApJ, 761, 164 [NASA ADS] [CrossRef] [Google Scholar]
 Csizmadia, S., Moutou, C., Deleuil, M., et al. 2011, A&A, 531, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dall, T. H., Santos, N. C., Arentoft, T., Bedding, T. R., & Kjeldsen, H. 2006, A&A, 454, 341 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Desort, M., Lagrange, A.M., Galland, F., Udry, S., & Mayor, M. 2007, A&A, 473, 983 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dumusque, X., Lovis, C., Ségransan, D., et al. 2011, A&A, 535, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Endl, M., MacQueen, P. J., Cochran, W. D., et al. 2011, ApJS, 197, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Figueira, P., Marmier, M., Bonfils, X., et al. 2010a, A&A, 513, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Figueira, P., Pepe, F., Melo, C. H. F., et al. 2010b, A&A, 511, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hébrard, G., Udry, S., Lo Curto, G., et al. 2010, A&A, 512, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Henry, T. J., Soderblom, D. R., Donahue, R. A., & Baliunas, S. L. 1996, AJ, 111, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Huélamo, N., Figueira, P., Bonfils, X., et al. 2008, A&A, 489, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lovis, C., Dumusque, X., Santos, N. C., et al. 2011, A&A, submitted [arXiv:1107.5325] [Google Scholar]
 Martín, E. L., Guenther, E., Zapatero Osorio, M. R., Bouy, H., & Wainscoat, R. 2006, ApJ, 644, L75 [NASA ADS] [CrossRef] [Google Scholar]
 Martínez Fiorenzano, A. F., Gratton, R. G., Desidera, S., Cosentino, R., & Endl, M. 2005, A&A, 442, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
 Mayor, M., Udry, S., Naef, D., et al. 2004, A&A, 415, 391 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Meunier, N., Lagrange, A.M., & Desort, M. 2010, A&A, 519, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nardetto, N., Mourard, D., Kervella, P., et al. 2006, A&A, 453, 309 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Niedzielski, A., Nowak, G., Adamów, M., & Wolszczan, A. 2009, ApJ, 707, 768 [NASA ADS] [CrossRef] [Google Scholar]
 Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN, The art of scientific computing [Google Scholar]
 Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saar, S. H., & Donahue, R. A. 1997, ApJ, 485, 319 [NASA ADS] [CrossRef] [Google Scholar]
 Santerne, A., Díaz, R. F., Moutou, C., et al. 2012, A&A, 545, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Mayor, M., Naef, D., et al. 2000, A&A, 361, 265 [NASA ADS] [Google Scholar]
 Santos, N. C., Israelian, G., Mayor, M., Rebolo, R., & Udry, S. 2003a, A&A, 398, 363 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Udry, S., Mayor, M., et al. 2003b, A&A, 406, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Comparison between BIS and V_{span} correlations with RV for HARPS data for the parameters that do not take into consideration the error bars.
Comparison between BIS and V_{span} correlations with RV for HARPS data for the parameters that consider the error bars, i.e. the slope fit parameters.
Comparison between correlations of BIS^{+} and BIS^{−} with RV for simulated data.
Comparison between BIS^{+} and BIS^{−} correlations with RV for HARPS data for the parameters that do not consider the error bars.
Comparison between BIS^{+} and BIS^{−} correlations with RV for HARPS data for the parameters that consider the error bars, i.e. the slope fit parameters.
Comparison between BIS and ΔV correlations with RV for HARPS data for the parameters that do not consider the error bars.
Comparison between BIS and ΔV correlations with RV for HARPS data for the parameters that consider the error bars, i.e. the slope fit parameters.
The α correction factor on photon noise used to obtain the observed noise distribution.
All Figures
Fig. 1 Comparison of biGaussian functions with different asymmetry values: 0% (Gaussian function – solid line), 10% (dashed line), 25% (dasheddotted line), and 50% (dotted line). These normalized functions have a contrast of 30%, a FWHM of 5 km s^{1}, and are all centered at 0 m/s. 

Open with DEXTER  
In the text 
Fig. 2 Most important properties for the two fitted functions: Gaussian as open circles (only two upper panels, partially superposed) and BiGaussian as filled circles, as a function of orbital phase (with an arbitrary zeropoint) for the case of an equatorial spot with a a filling factor of 1% on an edgeon star with v sini of 7 km s^{1}. Upper left: center of the two fitted functions as a function of orbital phase; Upper right: FWHM of the two functions; on the two lower panels, we have the asymmetry A (bottom left) and the ΔV indicator (bottom right) for the BiGaussian fitting. 

Open with DEXTER  
In the text 
Fig. 3 Two indicators considered, BIS (as open circles) and ΔV (as filled circles), as a function of RV, for the case of an equatorial spot with a filling factor of 1% on an edgeon star with vsini of 7 km s^{1}. 

Open with DEXTER  
In the text 
Fig. 4 Dependence between the indicator ΔV (filled circles) and BIS (open circles) with RV for HD 224789. 

Open with DEXTER  
In the text 
Fig. 5 RV (top) and V_{asy} (bottom) as a function of orbital phase (with an arbitrary zero point) for the case of an equatorial spot with a a filling factor of 1% on an edgeon star with vsini of 7 km s^{1}. 

Open with DEXTER  
In the text 
Fig. 6 Two indicators considered, BIS (as open circles) and V_{asy} (as filled circles), as a function of RV for the case of an equatorial spot with a a filling factor of 1% on an edgeon star with vsini of 7 km s^{1}. 

Open with DEXTER  
In the text 
Fig. 7 Dependence between the indicator V_{asy} (filled circles) and BIS (open circles) with RV for HD 224789. 

Open with DEXTER  
In the text 