Issue |
A&A
Volume 512, March-April 2010
|
|
---|---|---|
Article Number | A46 | |
Number of page(s) | 10 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/200913525 | |
Published online | 30 March 2010 |
The HARPS search for southern extra-solar planets
XX. Planets around the active star BD
,![[*]](/icons/foot_motif.png)
G. Hébrard1 - S. Udry2 - G. Lo Curto3 - N. Robichon4 - D. Naef2,5 - D. Ehrenreich1,6 - W. Benz7 - F. Bouchy1,8 - A. Lecavelier des Etangs1 - C. Lovis2 - M. Mayor2 - C. Moutou9 - F. Pepe2 - D. Queloz2 - N. C. Santos10 - D. Ségransan3
1 - Institut d'Astrophysique de Paris, UMR7095 CNRS, Université Pierre & Marie Curie,
98bis boulevard Arago, 75014 Paris, France
2 -
Observatoire de Genève, Université de Genève, 51 Chemin des Maillettes, 1290 Sauverny, Switzerland
3 -
ESO, Karl-Schwarzschild-Strasse 2, 85748 Garching bei München, Germany
4 -
Observatoire de Paris, GEPI, 5 Place Jules Janssen, 92195 Meudon, France
5 -
ESO, Alonso de Cordoba 3107, Vitacura Casilla 19001, Santiago, Chile
6 -
Laboratoire d'Astrophysique de Grenoble,
CNRS (UMR 5571), Université J. Fourier, BP 53, 38041 Grenoble Cedex 9, France
7 -
Physikalisches Institut Universität Bern, Sidlerstrasse 5, 3012 Bern, Switzerland
8 -
Observatoire de Haute-Provence, CNRS/OAMP, 04870 Saint-Michel-l'Observatoire, France
9 -
Laboratoire d'Astrophysique de Marseille, Université de Provence, CNRS (UMR 6110),
BP 8, 13376 Marseille Cedex 12, France
10 -
Centro de Astrofísica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal
Received 22 October 2009 / Accepted 8 December 2009
Abstract
We report the detection of a planetary system around BD
that includes at least one
Uranus-mass planet and one Saturn-mass planet. This discovery serendipitously originates from
a search for planetary transits in the Hipparcos photometry database. This program preferentially selected
active stars and did not allow the detection of new transiting planets.
It allowed however the identification of the K3V star BD
as a target
harboring a multiplanet system, which we secured and
characterized thanks to an intensive monitoring with the HARPS spectrograph at the 3.6-m
ESO telescope in La Silla. The stellar activity level of BD
complicates the analysis but does
not prohibit the detection of two planets around this star.
BD
b has a minimum mass of
and an orbital period of
5.60 days, whereas BD
c has a minimum mass of
and an orbital
period of 237.6 days.
This new system strengthens the observation that low-mass planets are preferentially found in
multiplanetary systems, but not around high-metallicity stars as is the case for
massive planets. It also supports the belief
that active stars should not be neglected in exoplanet searches even when searching
for low-mass planets.
Key words: planetary systems - techniques: radial velocities - stars: individual: BD
- techniques: photometric
1 Introduction
The HARPS spectrograph (Mayor et al. 2003) is operating since 2003 at the 3.6-m ESO telescope in La Silla, Chile. This is a fiber-fed, environmentally stabilized high-resolution echelle spectrograph dedicated to high-precision radial velocity measurements. Thanks to the exoplanetology programs that are conducted with it, it allowed numerous extra-solar planets studies and discoveries in the ranges of Jupiter-mass planets (e.g. Pepe et al. 2004; Moutou et al. 2009a), low-mass planets (e.g. Lovis et al. 2006; Mayor et al. 2009), planets around early- (e.g. Desort et al. 2008; Lagrange et al. 2009) or late-type stars (e.g. Bonfils et al. 2007; Forveille et al. 2009), or transiting planets (e.g. Bouchy et al. 2008; Queloz et al. 2009). The essential quality of HARPS (the High Accuracy Radial velocity Planet Searcher) is its high stability, which results in a sub-m s-1 accuracy of the radial velocity measurements on time scales of several years. This allows the detection of planets in the Neptune- or Super-Earth mass ranges on progressively increasing orbital periods as the time span of the monitoring increases. This instrument hitherto plays a major role in the improvement of the knowledge of exoplanets.
We report here the detection of a new multiplanetary system. It serendipitously
originates from the radial velocity follow-up that accompanied a search
for planetary transits in the Hipparcos epoch photometry annex
(Perryman et al. 1997). Despite its significant
stellar activity, the star BD
was identified during this follow-up as a possible
host of extra-solar planets. An intensive monitoring of its radial velocities was
subsequently carried out with HARPS and led to the detection of a planetary
system around this star, including at least two planets.
We briefly present the search for planetary transits in the Hipparcos photometry in
Sect. 2. The HARPS observations of BD
are presented
in Sect. 3 together with the stellar properties of this
target. The planetary system is characterized in Sect. 4,
the results are then discussed in Sect. 5.
2 Search for transiting planets in the Hipparcos database
The goal of our initial search in the Hipparcos database was to try to find new transiting planets, especially around bright stars. Transiting planets could allow numerous important studies to be performed. These studies include planetary radii, masses and density measurements, accurate determination of the inclinations and the eccentricities of the orbits, possible detections of transit timing variations due to additional companions, measurements of the spin-orbit (mis)alignment angles, detections of the emission or absorption planetary spectra, or even satellites or rings detections. The famous planets transiting in front of the bright stars HD 209458 and HD 189733 are those that allow the most accurate measurements and the deepest investigations.
The Hipparcos epoch photometry annex contains between 40 and
300 photometric measurements performed during the 1226-day duration
of the mission for each of the 118 204 stars of the catalog with the
magnitude limit of
.
With a chance of about 0.1% that a given star harbors a transiting extra-solar
planet, the Hipparcos all-sky survey must contain photometric measurements for
tens of transiting hot Jupiters. This belief is reinforced by the two cases of
transiting planets a posteriori re-discovered in the Hipparcos data:
HD 209458b (Robichon & Arenou 2000; Castellano et
al. 2000; Söderhjelm 1999) and HD 189733b
(Bouchy et al. 2005; Hébrard & Lecavelier des
Etangs 2006). Thanks to the long time-baseline, these data also
allow the orbital period to be measured with an accuracy on the order
of seconds.
Up to now, there are the only two transiting planets that have been found in the Hipparcos database. Indeed, with photometric variations of about 1% or less, such transits are difficult to identify in the Hipparcos data, which present a mean photometric accuracy on the same order of magnitude. The poor time coverage compared with dedicated surveys like SuperWASP, XO, CoRoT or Kepler, is another difficulty. Jenkins et al. (2002) concluded that because of its poor photometric quality, the Hipparcos catalog is no likely place to detect planets in the absence of other informations. It might however provide planetary transit candidates for follow-up observations. Laughlin (2000) has searched in the Hipparcos epoch photometry for transiting planets within 206 metal-rich stars. None have been confirmed thereafter.
We present in Appendix A the systematic search we managed in the Hipparcos epoch photometry annex for periods compatible with planetary transits. We constructed a ranked list of candidates for follow-up radial velocity measurements, based on the depth and significance of the planetary transits that could be detected in the Hipparcos photometry. The simulations we performed indicated that our detection rate for transiting planets is low, on the order of 2%; so the chances to find new transiting planets in the Hipparcos data are small (see Appendix A). Still, the particularly high interest for the potential discovery of a planet transiting a bright star pushed us to perform a radial velocity follow-up of our candidates.
We observed 194 of these selected, ranked targets with HARPS
to search for radial velocity variations in agreement with
the transiting candidates found in the Hipparcos photometry. The HARPS setup,
the spectra extraction and the radial velocity measurements were identical
to those described below in Sect. 3.1. Most of these
observations were performed in December 2004 as part of the
program 074.C-0364.
Radial velocity variations larger than 20 m s-1 were measured for 37 stars,
i.e. 19% of our observed sample. Those variations are mainly caused by stellar
activity. Indeed, most of the targets with large radial velocity variations show
emissions in the Ca II H & K lines (3934.8 Å and 3969.6 Å).
The level of these emissions is quantified with the activity S-index (Mount Wilson System),
which is converted to the
index (Santos et al. 2000).
The majority of the variable targets have
indexes larger
than -4.7, indicating prominent chromospheric activity; this can explain such a
level of radial velocity variations (Santos et al. 2000).
Some targets first seemed to exhibit periodic radial velocity variations, but subsequent monitoring with HARPS showed that these variations were not persistent with time, in agreement with transient activity processes on the stellar surface due to activity (flares, spots, plages...) modulated by the stellar rotation. In addition, in most of the cases an analysis of the line profiles using bisectors of the cross-correlation functions show variations in the shape of the lines with time. This indicates that the observed radial velocity variations are not due to Doppler shifts of the lines, but rather to deformations of the shape of the spectral lines. Some cases revealed a clear anti-correlation between radial velocity and line-bisector orientation, which could be understood as the signature of cool spots on the stellar photospheres (see, e.g., Queloz et al. 2001; Melo et al. 2007; Desort et al. 2007; Boisse et al. 2009).
Thus, our procedure seems to preferentially select active stars. This should be contrasted with transiting candidates obtained from the photometric surveys dedicated to transiting planet searches, which are mainly eclipsing binaries or transiting planets (see, e.g., Pont et al. 2005). None of the targets that we observed through this program have shown radial velocity variations in agreement with a reflex motion due to a hot Jupiter. Instead, most of the radial velocity variations seem to be caused by stellar activity.
However, these observations allowed the fortunate discovery of a
new planetary system without a transit detection in the Hipparcos data. This system
orbits BD
,
which is an active star, but among those with the lowest
indexes within our monitored targets. This activity
level complicates planet detection, but does not prohibit it in this case.
We concentrate below on this target and the detection of its harbored
planetary system.
3 Observations and properties of BD
3.1 HARPS observations of BD
We observed BD
using the HARPS spectrograph
at the 3.6-m ESO telescope in La Silla.
The bandpass of the spectra ranges from 3800 Å to 6900 Å, and the
resolution power is R=115 000, with a fiber diameter of 1 arcsec.
The spectra were extracted from the detector images with the HARPS pipeline,
which includes localization of the spectral
orders on the 2D-images, optimal order extraction,
cosmic-ray rejection, wavelength calibration and corrections of flat-field.
The pipeline then performs a cross-correlation of the extracted spectra
with a numerical mask (K5-type in this case), and finally measures the radial
velocities from Gaussian fits of the cross-correlation functions (CCFs),
following the method described by Baranne et al. (1996)
and Pepe et al. (2002).
The full dataset we use for BD
includes 83 spectra.
All the exposures (but the first one) were obtained without simultaneous thorium calibration.
The exposure times for BD
range between 4 and 15 minutes, allowing a
signal-to-noise ratio per pixel between 30 and 80 to be reached around 550 nm.
The total exposure time is about 11 hours.
The HARPS radial velocities of BD
as a function of time are plotted in the upper panel
of Fig. 1. They span 1841 days, corresponding to 5.0 years.
A few measurements were first obtained late in 2004 as part of the program 074.C-0364
(see Sect. 2), after the target was regularly monitored.
A few extra measurements were acquired in late 2006/early 2007 as part of the
program 078.C-0044, but most of the data were secured as part of the
Guaranteed Time Observations (GTO) survey
program 072.C-0488 (Mayor et al. 2003). This is the case in particular
for the intensive
series in 2007-2008 (BDJ
-2 400 000 between 54 150 and 54 650).
![]() |
Figure 1:
Top: radial velocity HARPS measurements of BD
|
Open with DEXTER |
The derived radial velocities are reported in Table 1. The accuracies are between 1.0 and 2.4 m s-1, typically around 1.6 m s-1. This includes photon noise but not the jitter due to stellar oscillations or activity (Sect. 3.2). Wavelength calibration and spectrograph drift uncertainties are negligible with respect to the photon noise. Four exposures with radial velocity uncertainties larger than 2.4 m s-1 were not included in the final dataset of 83 spectra.
The radial velocities show a significant 11-m s-1 dispersion (43 m s-1 peak-to-peak),
well above the expected accuracy (Fig. 1, upper panel). The CCF from which
those radial velocities were measured show parameters that also significantly vary
with time. Their full widths at half maximum fluctuates between
and
km s-1, and their contrasts from
and
%
of the continuum. The bisector span of the CCF also show a dispersion at the level
of 7 m s-1 (see Fig. 2, upper panel).
Thus, Doppler shifts of the spectral lines of BD
are not necessarily the only
explanation for the observed radial velocity variations; changes in the shape of the lines,
as those due to stellar activity, are likely to be the cause of at least a part of those variations.
Table 1:
HARPS measurements of BD
(full table available electronically)
.
3.2 Stellar properties and activity
Table 2 summarizes the stellar parameters of BD
(HIP 49067, SAO 137286). According to the SIMBAD database, this is a
K3V star of the magnitude V=9.86. Its Hipparcos parallax (
mas)
implies a distance of
pc. The Hipparcos color is
.
From the spectral analysis of the HARPS data using the method presented in Santos et al. (2004a), we derive the temperature
K, the
gravity
,
the metallicity
,
and an uncertain age of
Gyr.
The stellar mass we obtain is
,
with a formal error bar of
.
Following Fernandes & Santos (2004),
we rather adopt a conservative
% uncertainty,
corresponding to
.
We derive a projected rotational velocity
km s-1 from the parameters of the CCF using
a calibration similar to that presented by Santos et al. (2002).
Acoustic oscillations are not averaged out in the data as such oscillations are on
timescales similar or shorter than the exposure times we used. However,
for a K-type star, oscillation amplitudes are expected to be negligible by
comparison with stellar activity effects on the measured radial velocities.
Table 2:
Stellar parameters for BD
.
![]() |
Figure 2: Bisector span as a function of the radial velocity ( upper panel) and of the residuals of the fit with two planets as shown in Figs. 1 and 4 ( lower panel). On the lower panel the filled circles indicate the measurements secured on the 101-day interval, where we fitted the stellar jitter (Sect. 4.2 and Fig. 6). The ranges have the same extents in thex- and y-axes on both panels. |
Open with DEXTER |
The core of the large H & K Ca II absorption lines of BD
shows emissions, which indicates chromospheric activity.
We adopted the averaged value
,
but it significantly varies with time, between
extrema
and
.
Such stellar activity would imply a significant jitter on the observed stellar radial
velocities. For a K-type star with this level of activity, Santos et al. (2000)
predict a dispersion of up to 10 m s-1 for the stellar jitter. This is the order of magnitude
of the dispersion of our measurements of BD
(Sect. 3.1).
A classic method to identify stellar activity as the main cause of radial velocity variations
is to look for anti-correlation between the radial velocity and the bisector span (see, e.g.,
Queloz et al. 2001; Boisse et al. 2009). The upper panel of
Fig. 2 shows the bisector spans as a function of the radial velocities
for BD
.
As seen above, the bisector is varying significantly; but no clear
relations are seen between the two signals. Yet this technic is less sensitive
for stars with slow rotations as BD
.
Indeed, according to
Noyes et al. (1984) and Mamajek & Hillenbrand (2008),
the activity level for BD
implies a stellar rotation period of
37 days.
The Hipparcos data of BD
include photometric measurements for 63 epochs only;
this does not allow a rotation period measurement.
![]() |
Figure 3:
Lomb-Scargle periodograms for BD
|
Open with DEXTER |
A different estimate of the rotation period could be obtained from the HARPS data through
the stellar activity indicators, as shown in Fig. 3. In this figure are
plotted the Lomb-Scargle periodograms (Press et al. 1992) of five HARPS signals:
the radial velocities, full widths at half maximum (FWHMs), contrasts and bisectors of the CCFs,
and of the
activity indexes (Table 1).
The four last parameters could show signatures of the stellar activity and in
particular their modulations with the stellar rotation; the radial velocity variations could also show
signatures of the stellar activity and rotation, but with added Doppler reflex motion due to companion(s).
The five periodograms in Fig. 3 all clearly show a peak around 26 days.
We interpret it as the signature of modulations of the stellar surface due to activity (flares, spots,
plages...). These kinds of phenomena are expected for such an active star. As there are sporadic
events that alter the surface of the star and have limited lifetimes, they would imply quasi-periodic
variations of the shapes of the spectral lines with periods near the stellar rotation period and
unconserved phases.
The rotation period of BD
consequently appears to be around 26 days.
Most of the signals also show a small peak near 13 days, which is the first harmonic of
the main signal; it is not detected in the radial velocities, however.
In addition, the five periodograms in Fig. 3 show a signal between 65
and 95 days. It may be the signature of the typical duration of cycles of
these phenomena.
The radial velocity periodogram shows two extra peaks that are not seen for other
signals in Fig. 3, at 5.6 and 238 days. The fact that these periodic signals
are detected only in radial velocities suggests it to be caused by Doppler reflex motion due to
companions for BD
rather than jitter due to stellar activity. Moreover, those two periods
do not correspond to any harmonics nor aliases of the 26 days rotation period.
The amplitude of the variations (43 m s-1 peak-to-peak, see
upper panel of Fig. 1)
and the two orbital periods (near 5.6 and 238 days) imply projected masses
well below the mass of Jupiter.
The radial velocity variations detected by HARPS would thus originate both from
stellar activity jitter and planetary companions.
4 A planetary system around BD
4.1 Fit without stellar activity
We first fitted the radial velocities using a two-planet Keplerian model without mutual
interactions. The results are plotted in Figs. 1 and 4.
The inner planet, BD
b, produces radial velocity variations with a small
semi-amplitude K=6.5 m s-1, corresponding to a planet with a minimum mass
,
thus similar to the mass of Uranus. Its orbit has a period of
5.60 days and is circular or with low eccentricity. The outer
planet, BD
c, yields a larger semi-amplitude, namely K=13.4 m s-1; this corresponds to a planet with a minimum mass
,
slightly above the Saturn mass.
The orbital period is 237.6 days, and the orbit is non-circular (e=0.19).
The derived orbital parameters of the system
are summarized in Table 3.
![]() |
Figure 4:
Phase-folded radial velocity curves for BD
|
Open with DEXTER |
Table 3:
Fitted orbits and planetary parameters for the BD
system, with 1-
error bars.
The reduced
of the Keplerian fit is 3.2, and the standard deviation of the residuals is
m s-1. This is reduced when compared to the 11-m s-1
dispersion of the original radial velocities, but is clearly higher than the 1.6-m s-1 typical
error bars on the individual measurements. The extra dispersion on the order of 4 m s-1 is
mainly due to stellar activity jitter as seen above. By fitting the two planets without including
the stellar activity, we assume that the stellar jitter would be averaged out over the five years of
observations, as activity induces quasi-periodic effects that are not exactly duplicated
in time
with the stellar rotation. This results in a quite large dispersion around the fit with
respect to the error bars, whereas the periodic signal of the two planets remains coherent
over the five years of observations.
Figure 5 strengthens the interpretation in term of planets of the two
signals at 5.6 and 238 days. The upper panel shows the periodogram
of the radial velocities, exactly like the upper panel of Fig. 3.
On the second panel of Fig. 5 is shown the periodogram of the
radial velocity residuals, after a fit including BD
c only. The standard deviation of the
residuals of this fit is
m s-1. On this periodogram,
the peak at 238 days of course is
no longer visible. The main peak is the one at 5.6 days.
Its false-alarm probability is
.
It is clearly stronger than the peaks due to activity.
We also note the presence of two other peaks at 0.8 and 1.2 days that are the
one-day aliases of the 5.6-day signal in the frequency space (
). We adopted
5.6 days as the actual period of this signal instead of one of the two aliases,
because it is more likely according to the sampling.
In addition, its peak is higher than those of the two aliases, and Keplerian fits performed with
0.8 or 1.2 day of orbital periods produce higher residuals dispersion and eccentric orbits.
Similarly, the third panel of Fig. 5 shows the periodogram of the
radial velocity residuals, after a fit including BD
b only. The peak at 5.6 days is removed
as well as its aliases, and the 238-day signal remains, clearly above the activity peaks.
Finally, the bottom panel of Fig. 5 shows the periodogram of the
residuals after the 2-planet fit shown in Figs. 1 and 4. Only
lower peaks are remaining, most of them caused by activity.
Their false-alarm probabilities are smaller than of
.
We note that in addition to the clear peak at one day
(which is due to the aliases of all the detected signals), all the four panels in
Fig. 5 show a possible peak near 700 days; this could be the signature
of a third, outer planet. Fits with a curvature in addition to the two planets also suggest a
possible long-period signal (without significant effects on the parameters of
BD
b and BD
c).
Such long-period signals are however not strong enough to claim any detection. Further
observations of this target on a longer time-baseline are mandatory to establish or not
the presence of the hypothetic planet BD
d.
![]() |
Figure 5:
Lomb-Scargle periodograms of the HARPS radial velocities. The upper panel shows the periodogram computed on the initial radial velocities without any fit
removed. The second and third panels show the periodograms computed on the residuals
of the fits including BD
|
Open with DEXTER |
The lower panel of Fig. 2 shows the bisector spans as a function of the radial velocities
residuals after the 2-planet fit. This does not show an anti-correlation between those two signals
as it could be expected in case of a spotted stellar surface due to activity (see, e.g., Melo et
al. 2007; Boisse et al. 2009). The absence of anti-correlation whereas
activity is supposed to be the main part of the radial velocity residuals in this case
could be due to
the fact that the rotation velocity is low, or to activity processes more complex than simple
cool stellar spots. It could also be due to the presence of additional companions
(as the hypothetic BD
d), which would induce radial velocity variations uncorrelated
with the bisectors.
We note also that we did not find any anti-correlation between the bisector spans and the radial velocities with either one or the other planet removed. Such anti-correlation with just one planet removed would argue that the other signal is due to activity; this is not the case.
4.2 Fit with stellar activity
In the previous section we have shown a Keplerian fit of the two planets without attempting any fit of the stellar activity jitter signal in the radial velocities. Our comprehension of the activity is poor; for example the number, locations and sizes of potential stellar spots are far to be known. We attempt here a naive, phenomenological approach, based on the fact that activity shows signals in the periodograms, so stellar jitter has a periodical nature that is linked to the rotation of the star. Unfortunately these signals are only quasi-periodical: they do not reproduce themselves periodically in an identical form for a long time, as all these phenomena have limited lifetimes. Fitting the stellar jitter with sinusoids over the five-year time span does not provide satisfactory fits.
![]() |
Figure 6:
Fit of the two planets and the stellar activity on a 101-day interval.
The upper plot shows
the HARPS radial velocities of BD
|
Open with DEXTER |
We thus tried such sinusoid fits on a shorter time span. We chose a 101-day interval, between
BJD-2 400 000 = 54 547.5 and 54 648.5. This interval has a good sampling (27
measurements), including two nine-night sequences with almost one measurement
per night (eight and nine measurements in the two sequences) at high accuracy (uncertainties
between 1.0 and 1.7 m s-1). This should improve the
coherence of the stellar jitter in terms of the periodicity, as well as allowing the 5.60-day-period
planet to be well sampled. Using sinusoids of periods near 26 days, we fitted
the activity indicators studied in Fig. 3
on this 101-day interval,
namely the FWHMs, contrasts and bisectors of the CCFs, and the
indexes. The best solution we obtained was for a
26.6-day period. Together with the periodograms shown in Fig. 3
this allowed us to determine the rotation period:
days.
We fitted the radial velocities on the same interval,
using Keplerian models for the two planets, and an extra sinusoid to approximate the
stellar jitter. The parameters of the inner planet were free, as was the mass of the outer one.
The period, eccentricity and orientation of the orbit of BD
c
were fixed to the values obtained above, the 101-day interval being too short to
constrain this 238-day-period orbit.
The results of these fits agree with those obtained in Sect. 4.1.
The fits are plotted in Fig. 6, which shows on the upper panel a
good match of the radial velocities and the model, especially on the two high-frequency
observation sequences that show the variations due to BD
b.
The lower panels show the activity signals.
The 26.6-day-period sinusoids provide acceptable approximations for their variations,
which still remain more complex than this simple sinusoid shape. The lower panel shows
the stellar jitter only, i.e. the sinusoid that is added to the two Keplerian in the upper
panel. It has an amplitude of
m s-1.
The residuals of the radial velocities fit with two Keplerian and a sinusoid is plotted in
the middle panel of Fig. 6. Its dispersion is reduced
to 2.5 m s-1, which is an improvement compared to the 4.3 m s-1 dispersion
obtained on the full dataset in Sect. 4.1 without an attempt to
fit the stellar jitter
(or 4.1 m s-1 dispersion if computed only on the 101-day interval, again without
stellar jitter fit).
It remains larger than the error bars because of the imperfection of
the sinusoid model we used for the stellar jitter
(and perhaps also because of possible extra planets).
We restrained this fit to a single
sinusoid for the stellar jitter without including extra sinusoids at the periods of
the stellar rotation harmonics. Indeed, these harmonics are not detected in the
radial velocities (see Sect. 3.2 and Fig. 3).
Despite its simplicity, the simple sinusoids allow correlations to be seen between
the different parameters on this 101-day time span (Fig. 6, lower panels).
The FWHMs and contrasts of the
CCFs are anti-correlated, implying spectral line deformations due to activity that
let their equivalent widths nearly constant. This could be understood as the
surface of the CCF is an indicator of the stellar metallicity (Santos et al. 2002).
The velocity jitter is barely correlated with
,
or with the FWHM and the bisector of the CCF. The correlation between the
FWHM and
is different from the observations made
for example in
the case of the spotted star CoRoT-7, for which those two values are rather
anti-correlated (Queloz et al. 2009). Also, the apparent correlation
between
and the radial velocity jitter is different from
the picture seen in the cases of the active stars
HD 166435 (Queloz et al. 2001), GJ 674 (Bonfils et al. 2007)
or HD 189733 (Boisse et al. 2009).
And finally, there is an apparent correlation between the bisector and the radial
velocity jitter on this 101-day time span, possibly with a small phase offset
(see Fig. 6). This possible correlation is also shown
on the lower panel of Fig. 2.
All these relations draw a picture for BD
activity which is quite different from a
simple scenario where the radial velocity jitter is mainly due to dark spots on the
stellar surface that modulate the shape of the lines as the star is rotating. Other
phenomena are likely to occur on this star; they may include pulsations, convections,
flares, plages, hot spots...
The effect of the stellar activity on the observed radial velocities of BD
is on the
order of 4 m s-1. This is lower than the semi-amplitudes measured for the two detected
planets, but not negligible. The error bars reported in Table 3
where derived from
variations and Monte Carlo experiments with and without
stellar activity modeling as well as from
trials and errors with different kinds of sinusoids for the stellar jitter.
Our poor understanding of the stellar activity of BD
made us remain cautious
on the obtained error bars.
The reality
of the two planets in this range of mass is well established however, their
signatures being clearly detected at two periods that show no signals linked to
stellar activity.
5 Discussion
Finding transiting planets in the Hipparcos epoch photometry annex does not look promising. Our attempt for a priori detections did not succeed, and up to now, only two a posteriori detections were performed within the Hipparcos data, in the cases of HD 209458b and HD 189733b which were first revealed from ground observations. Three main limitations make it difficult. First, the error bars on individual photometric measurements are on the same order of magnitude than the expected signal for transits of giant planets, or even slightly larger. Second, the sparse time-coverage allows only a few points to be obtained in a potential transit. These two limitations put any transit identification at the limit of detection. Third, stellar activity could produce false positives, which are difficult to identify with the sparse time-coverage. The ESA mission Gaia (Perryman et al. 2001) is awaited as the successor of Hipparcos. Its potential for planetary transit discoveries in front of bright stars could be considered, as its sensitivity and accuracy would be better than Hipparcos. Yet, the time coverage will not be better than that of Hipparcos, which will prohibit well-resolved light curve studies. Stellar activity should thus also be a limitation for transit detections with Gaia. Dedicated surveys have proven to be more efficient for transit detections. Ground-based programs allow Jupiter-size planets to be detected with improving accuracies which now point toward planets with smaller radii (Bakos et al. 2010), whereas space-based programs allow planets with even smaller radii to be detected (Léger et al. 2009; Queloz et al. 2009), as well as planets with longer periods (Moutou et al. 2009b). If accepted, the ESA-proposed space mission PLATO (Catala et al. 2009) could permit the detection of this kind of transiting planets in front of brighter stars.
If this search in the Hipparcos data did not provide any detection
of new transiting planets, it allowed the fortunate discovery of a new
planetary system around BD
.
This target was first identified as presenting
promising radial velocity variations, in which we identified the signature of
two new extra-solar planets thanks to an intensive monitoring with HARPS.
The inner planet, BD
b, is a hot Neptune with a minimum mass of
14.4
.
Its 5.60-day orbit could be circular or slightly eccentric.
The outer planet, BD
c, is slightly more massive than Saturn, with a
minimum mass of 0.33
.
Its orbit is moderately but significantly
eccentric and has a period of 237.6 days. As the masses are low and the
orbits are distant and nearly circular, the mutual interactions between the
two planets are negligible.
The reflex motions that these two planets induce on their host star have semi-amplitudes of 6.5 and 13.4 m s-1, which can be distinguished from the 4-m s-1 jitter due to stellar activity. We summarize here the arguments that allow us to conclude that these two signals are due to planets and not to stellar activity:
- the signals with periods of 5.60 and 237.6 days are seen only in the radial velocities and are not seen in the shapes of the lines nor in the activity indexes;
- these two periods do not correspond to the rotation period of the star, nor to its harmonics or its aliases;
- with the available data, the signal is coherent over the five-year time span of the observations.

Most of the Neptune-mass and Super-Earth planets found by the HARPS
GTO program are detected from a survey of about 200 non-active solar-type
stars. They are monitored with numerous high-precision measurements, using
simultaneous thorium calibration and exposures of 15-min duration to
average out the stellar oscillations (Santos et al. 2004b;
Udry et al. 2006; Lovis et al. 2006; Bouchy et al. 2009;
Mayor et al. 2009). This is not the case for BD
b, which was
discovered from lower signal-to-noise exposures without simultaneous
calibration and as a serendipitous result of a program that preferentially
selected active stars. BD
is one of the most active stars around which a
low-mass planet has been found, together with CoRoT-7
(Queloz et al. 2009) and a few M-dwarfs (Bonfils et al. 2007;
Forveille et al. 2009).
Those cases support the fact that active stars should not be neglected in exoplanet
hunts, even when searching for low-mass planets.
![]() |
Figure 7:
Stellar metallicities as a function of the planetary masses for 304 known extra-solar
planets (top), and for the 96 ones among them with orbital periods shorter than
20 days (bottom). BD
|
Open with DEXTER |
BD
has a metallicity which is slightly lower than that of the Sun. This strengthens
the observation that whereas giant planets are more frequently found around stars with
supersolar metallicities, this tendency is not seen for
Neptune-mass or Super-Earth planets.
This is shown in Fig. 7, which presents the minimum masses of known
exoplanets as a function of the metallicity of their hosting stars. The data are taken from
the compilation of the extrasolar planets encyclopedia
.
While planets with masses larger than 1
are more numerous around
overmetallic stars, they are clearly more numerous around undermetallic stars for
planetary masses below 20
.
The same effect is seen if we select only the
planets with orbital periods of 20 days or shorter, i.e. in the period regime where most
of the low-mass planets are detected up to now (lower panel of Fig. 7).
Indeed, among the
30 planets with a projected mass lower than
0.1
,
only four have orbital periods longer than 20 days (namely
Gl 581d, HD 40307d, HD 69830c and d).
Planet formation models based on the core
accretion scenario (Mordasini et al. 2009a,b) have shown that the sharply rising probability of detecting
giant planets with stellar metallicity can be at least partially accounted for because
metal rich systems favor the formation of massive planets and the radial
velocity technique is most sensitive to massive bodies. Observations also seem to
indicate that this correlation vanishes for Neptune-mass planets, a trend also found
in population synthesis calculations (Mordasini et al, in preparation). These calculations
even show that the correlation reverses for very small planetary masses
(<3-5
), a prediction that will have to be confirmed by observations.
In fact, a careful analysis of the theoretical models indicates that the critical parameter
is the overall mass of heavy elements rather than the metallicity. Since this mass is
determined by the metallicity and the mass of the initial proto-planetary disk (which is
unknown), it makes a straight interpretation more difficult. For example, models predict
that low mass objects orbiting metal rich stars or relatively massive planets orbiting
metal poor stars are also possible, but they should be rare.
The BD
system also agrees with the tendency of low-mass planets to be preferentially
found in multiple systems. Almost 70% of the
30 planets with projected masses
lower than 0.1
are found in multiplanet systems, whereas this ratio
is only
25% for all the
350 known planets. Again considering only the known
planets with orbital periods shorter than 20 days, less than
20% of them are found
to be in multiplanet systems: those
20% are mainly low-mass planets.
No photometric search for transits have been managed for BD
from
follow-up ground-based observations; depending on the unknown inclination
i of the orbit, the transit probability for BD
b is about 9%. Its expected depth
is out of reach of the Hipparcos photometric accuracy.
We would like to thank F. Pont, F. Arenou, I. Boisse, X. Bonfils, J.-M. Désert, R. Ferlet, J. Laskar, D. Sosnowska and A. Vidal-Madjar for help and discussions, as well as the different observers from other HARPS programs who have also measured BD. We are grateful to the ESO staff for their support on the HARPS instrument. G.H., F.B. and A.L.D.E. acknowledge support from the French National Research Agency (ANR-08-JCJC-0102-01). D.E. acknowledges financial support from the CNES. N.C.S. would like to thank the support by the European Research Council/European Community under the FP7 through a Starting Grant, as well from Fundação para a Ciência e a Tecnologia (FCT), Portugal, through programme Ciência 2007, and in the form of grants reference PTDC/CTE-AST/098528/2008 and PTDC/CTE-AST/098604/2008.
Appendix A: Selection of targets for transiting planets search in the Hipparcos database
We present in this Appendix the systematic search we managed in the Hipparcos epoch photometry for planetary transit candidates (Sect. 2).
First, we selected 23 304 stars among the 118 204 of the Hipparcos epoch photometry catalog, according to the following criteria:
- B-V>0.4 (F2 type and later);
-
(defined parallax);
- empty H48 field (reference flag for photometry);
- H52 field different from D, P, or R (variability types);
-
(small stellar radius).
Then we kept only targets with at most one epoch brighter than
3
from the average magnitude (to remove remaining
variables) and at least 40 different available epochs. We performed
the periods' research on the 17 800 remaining stars.
For each target, we took the two faintest epochs and scanned all
possible periods to have those two points in a transit. Then
we adopted the period producing the lowest flux in the drop and
quantified this solution with the parameter ,
defined as:

![]() |
Figure A.1:
Threshold tuning simulations.
The parameters |
Open with DEXTER |
where
and
are the averaged magnitudes
of the epochs in and out the drop respectively, and
is the
standard deviations. The higher
is, the deeper and more
significant is the drop.
We also computed the
for the fit of the epochs in the drop
with a transit curve for a planet with a radius
,
implying a
drop in the magnitude:

We finally computed the goodness of fit (gof) for this


We performed two simulations in order to tune the threshold: for the
17 800 stars selected above, we assumed in the first simulation that
all the stars are constant, whereas in the second simulation, we
assumed that all the stars host a transiting planet with a radius
chosen randomly within
.
We kept the
individual error on the photometry on each epoch.
The upper panel of Fig. A.1 shows the goodness of fit as
a function of
for the first simulation. The distribution is skew
as
and the goodness of fit are correlated.
In order to simplify the selection,
we fitted the distribution with a line and computed
,
which is
corrected (with a first-order polynomial) from this skew. The middle panel
of Fig. A.1 shows the goodness of fit as a function of
for the two simulations.
appears to be an appropriate
threshold for the transit candidates selection: indeed, 127 and 489 stars
are above this limit in the first and second simulation,
respectively. We can thus estimate the detection rate as
(489-127)/17800=2%.
With 0.06% chance that a given star harbors a
transiting extra-solar planet, the 17 800 targets of our latest sample
should include about 11 transiting planets. The expected transiting
planet detection number with our method should thus be slightly less
than unity (
).
The actual data for the 17800 selected targets are plotted in the lower
panel of Fig. A.1. The stars with high and a high goodness of fit are eclipsing binaries that were not identified
in the Hipparcos catalog. To remove those binaries we removed targets
with a goodness of fit larger than 2.5. We also removed targets with
,
and obtained a list of
candidates sorted by decreasing
.
Targets referenced in SIMBAD as binaries, active, or variable stars were removed from the obtained list, as well as targets with already known planets at the time of the observing run (HD 70642, HD 39091, HD 10647, HD 4208, HD 17051, HD 13445, HD 75289, and HD 83443 - we checked that the Hipparcos photometric data did not include periodic variations at their periods).
Using HARPS, we performed follow-up observations of 194 of these selected targets in December 2004, as part of the program 074.C-0364. The apparent magnitudes of the observed stars range from 4.9 to 11.5. Exposures of typically a few minutes' duration were obtained, allowing 70 to 90 targets to be observed each night. Errors on the measured radial velocities are typically on the order of 2 m/s. To identify hot Jupiters, targets were removed from the candidate list after one, two, or three observations, according to the following observational strategy:
- observations were stopped after one HARPS measurement if: 1) there were two peaks in the CCF; only one SB2 was found during our observations (HD 23919), as binaries were previously removed from the candidate list from Hipparcos flags and SIMBAD checks; 2) the full width at half maximum of the CCF was larger than 15 km s-1, which prohibits accurate radial velocity measurements (17 stars);
-
observations were stopped after two HARPS measurements if:
1) the radial velocity difference
between the two measurements was larger than 2.5 km s-1 (no SB1 were found); 2)
m/s (113 stars);
-
observations were stopped after three HARPS measurements if
m/s (26 stars).
Radial velocity variations larger than 20 m s-1 were measured for 37 stars, i.e. 19% of our observed sample. However, those variations are mainly caused by stellar activity (see Sect. 2).
References
- Bakos, G. A., Torres, G., Pál, A., et al. 2010, ApJ, 710, 1724 [NASA ADS] [CrossRef] [Google Scholar]
- Baranne, A., Queloz, D., Mayor, M., et al. 1994, A&AS, 119, 373 [Google Scholar]
- Boisse, I., Moutou, C., Vidal-Madjar, A., et al. 2008, A&A, 495, 959 [Google Scholar]
- Bonfils, X., Mayor, M., Delfosse, X., et al., A&A, 474, 293 [Google Scholar]
- Bouchy, F., Udry, S., Mayor, M., et al. 2005, A&A, 444, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bouchy, F., Queloz, D., Deleuil, M., et al. 2008, A&A, 482, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bouchy, F., Mayor, M., Lovis, C., et al. 2009, A&A, 496, 527 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Castellano, T., Jenkins, J., Trilling, D. E., Doyle, L., & Koch, D. 2000, ApJ, 532, L51 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Catala, C. 2009, ExA, 23, 329 [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]
- Desort, M., Lagrange, A.-M., Galland, F., et al. 2008, A&A, 491, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fernandes, J., & Santos, N. C. 2004, A&A, 427, 607 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Forveille, T., Bonfils, X., Delfosse, X., et al. 2009, A&A, 493, 645 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hébrard, G., & Lecavelier des Etangs, A. 2006, A&A, 445, 341 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jenkins, J. M., Caldwell, D. A., & Borucki, W. J. 2002, ApJ, 564, 495 [NASA ADS] [CrossRef] [Google Scholar]
- Lagrange, A.-M., Desort, M., Galland, F., Udry, S., & Mayor, M. 2009, A&A, 495, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Laughlin, G. 2000, ApJ, 545, 1064 [NASA ADS] [CrossRef] [Google Scholar]
- Léger, A., Rouan, D., Schneider, J., et al. 2009, A&A, 506, 287 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
- Lovis, C., Mayor, M., Pepe, F., et al. 2006, Nature, 441, 305 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Mamajek, E. E., & Hillenbrand, L. A. 2008, ApJ, 687, 1264 [NASA ADS] [CrossRef] [Google Scholar]
- Mayor, M., Pepe, F., Queloz, D., et al. 2003, Msngr., 114, 20 [Google Scholar]
- Mayor, M., Udry, S., Lovis, C., et al. 2009, A&A, 493, 639 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Melo, C., Santos, N. C., Gieren, W., et al. 2007, A&A, 467, 721 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mordasini, C., Alibert, Y., & Benz, W. 2009a, A&A, 501, 1139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mordasini, C., Alibert, Y., Benz, W., & Naef, D. 2009b, A&A, 501, 1161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Moutou, C., Mayor, M., Lo Curto, G., et al. 2009a, A&A, 496, 513 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Moutou, C., Bruntt, H., Guillot, T., et al. 2009b, A&A, 488, L47 [Google Scholar]
- Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763 [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]
- Pepe, F., Mayor, M., Queloz, D., et al. 2004, A&A, 423, 385 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 [NASA ADS] [Google Scholar]
- Perryman, M. A. C., de Boer, K. S., Gilmore, G., et al. 2001, A&A, 369, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pont, F., Bouchy, F., Melo, C., et al. 2005, A&A, 438, 1123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in C. The art of scientiÞc computing, 2nd edn. (Cambridge: University Press) [Google Scholar]
- Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Queloz, D., Bouchy, F., Moutou, C., et al. 2009, A&A, 506, 303 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Robichon, N., & Arenou, F. 2000, A&A, 355, 295 [NASA ADS] [Google Scholar]
- Santos, N. C., Mayor, M., Naef, D., et al. 2000, A&A, 361, 265 [NASA ADS] [Google Scholar]
- Santos, N. C., Mayor, M., Naef, D., et al. 2002, A&A, 392, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Santos, N. C., Israelian, G., & Mayor, M. 2004a, A&A, 415, 1153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Santos, N. C., Bouchy, F., Mayor, M., et al. 2004b, A&A, 426, L19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Söderhjelm, S. 1999, Informational Bulletin on Variable Stars, 4816, 1 [Google Scholar]
- Udry, S., Mayor, M., Benz, W., et al. 2006, A&A, 447, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Footnotes
- ...
- Based on observations made with HARPS spectrograph on the 3.6-m ESO telescope at La Silla Observatory, Chile, under the programs ID 072.C-0488, 074.C-0364 and 078.C-0044.
- ...
- The full version of Table 1 (HARPS
measurements of BD
) is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/512/A46
- ... encyclopedia
- http://exoplanet.eu
All Tables
Table 1:
HARPS measurements of BD
(full table available electronically)
.
Table 2:
Stellar parameters for BD
.
Table 3:
Fitted orbits and planetary parameters for the BD
system, with 1-
error bars.
All Figures
![]() |
Figure 1:
Top: radial velocity HARPS measurements of BD
|
Open with DEXTER | |
In the text |
![]() |
Figure 2: Bisector span as a function of the radial velocity ( upper panel) and of the residuals of the fit with two planets as shown in Figs. 1 and 4 ( lower panel). On the lower panel the filled circles indicate the measurements secured on the 101-day interval, where we fitted the stellar jitter (Sect. 4.2 and Fig. 6). The ranges have the same extents in thex- and y-axes on both panels. |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Lomb-Scargle periodograms for BD
|
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Phase-folded radial velocity curves for BD
|
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Lomb-Scargle periodograms of the HARPS radial velocities. The upper panel shows the periodogram computed on the initial radial velocities without any fit
removed. The second and third panels show the periodograms computed on the residuals
of the fits including BD
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Fit of the two planets and the stellar activity on a 101-day interval.
The upper plot shows
the HARPS radial velocities of BD
|
Open with DEXTER | |
In the text |
![]() |
Figure 7:
Stellar metallicities as a function of the planetary masses for 304 known extra-solar
planets (top), and for the 96 ones among them with orbital periods shorter than
20 days (bottom). BD
|
Open with DEXTER | |
In the text |
![]() |
Figure A.1:
Threshold tuning simulations.
The parameters |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.