The HARPS search for southern extra-solar planets - XVIII. An Earth-mass planet in the GJ 581 planetary system (original) (raw)

A&A 507, 487-494 (2009)

XVIII. An Earth-mass planet in the GJ 581 planetary system

M. Mayor1- X. Bonfils2,3 - T. Forveille2- X. Delfosse2 - S. Udry1- J.-L. Bertaux4 - H. Beust2- F. Bouchy5 - C. Lovis1- F. Pepe1 - C. Perrier2- D. Queloz1 - N. C. Santos1,6

1 - Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, 1290 Sauverny, Switzerland
2 - Laboratoire d'Astrophysique, Observatoire de Grenoble, Université J. Fourier, CNRS (UMR5571), BP 53, 38041 Grenoble, Cedex 9, France
3 - Centro de Astronomia e Astrofísica da Universidade de Lisboa, Observatório Astronómico de Lisboa, Tapada da Ajuda, 1349-018 Lisboa, Portugal
4 - Service d'Aéronomie du CNRS, BP 3, 91371 Verrières-le-Buisson, France
5 - Institut d'Astrophysique de Paris, CNRS, Université Pierre et Marie Curie, 98bis Bd Arago, 75014 Paris, France
6 - Centro de Astrofísica, Universidade do Porto, Rua das Estrelas, P4150-762 Porto, Portugal

Received 24 March 2009 / Accepted 14 June 2009

Abstract
The GJ 581 planetary system is already known to harbour three planets, including two super-Earth planets that straddle its habitable zone. We report the detection of an additional planet - GJ 581e - with a minimum mass of 1.9 $M_\oplus$. With a period of 3.15 days, it is the innermost planet of the system and has a $\sim$5% transit probability. We also correct our previous confusion about the orbital period of GJ 581d (the outermost planet) with a one-year alias, benefitting from an extended time span and many more measurements. The revised period is 66.8 days, and positions the semi-major axis inside the habitable zone of the low mass star.

The dynamical stability of the 4-planet system imposes an upper bound on the orbital plane inclination. The planets cannot be more massive than approximately 1.6 times their minimum mass.

Key words: planetary systems - stars: late-type - techniques: radial velocities - stars: individual:GJ 581

1 Introduction

HARPS is a vacuum spectrograph designed to measure precise radial velocities, with the specific goal of searching for exoplanets in the southern hemisphere (Mayor et al. 2003). This high-resolution Echelle spectrograph (_R_=115 000) is fiber-fed by the ESO 3.6-m telescope at La Silla Observatory. The consortium that constructed this instrument was granted 500 observing nights over five years to conduct a comprehensive search for exoplanets. This large programme addresses several key exoplanet questions, including the characterization of very low mass exoplanets. A significant fraction of the observing time was devoted to the study of planets orbiting stars at the bottom of the main sequence. The M-dwarf sub-programme includes over 100 stars, which form a volume-limited sample.

Our HARPS search for southern exoplanets orbiting M-dwarfs is ongoing, and has been expanded to a larger sample of some 300 low-mass stars. Several planets have already been detected in the M-dwarf survey: GJ 581 b (Bonfils et al. 2005b), GJ 581 c and d (Udry et al. 2007), GJ 674 (Bonfils et al. 2007), GJ 176 (Forveille et al. 2009), and here a fourth planet in theGJ 581 system. These 6 planets all have minimum masses of less than approximately 15 Earth-masses. They represent approximately half of the known inventory of planets orbiting M stars, and most of its lowest mass members. Besides the HARPS planets, this inventory includes three planets around GJ 876 (Marcy et al. 1998; Delfosse et al. 1998; Rivera et al. 2005; Marcy et al. 2001), a single Neptune-mass planet orbitingGJ 436 (Butler et al. 2004), Jovian planets bound to both GJ 849 (Butler et al. 2006) and to GJ 317 (Johnson et al. 2007), and finally a gaseous giant planet orbiting GJ 832 (Bailey et al. 2009). Altogether, 8 planetary systems centered on M dwarfs have been identified by the radial velocity technique, for a total of 12 planets. Although the statistics are still limited, and come from surveys with different sensitivities, they indicate that multi-planet systems are common: the fraction of known multiple systems around M dwarfs is already 2/8 (25%), and it reaches 5/8 (63%) if one accounts for systems where a radial velocity drift indicates an additional long-period planet.

The competing teams monitoring M dwarfs with Doppler spectroscopy have together observed over 300 stars for several years. Planetary systems with at least one identified gaseous giant planet have been found around $\sim$1.5% of these low mass stars, and that proportion is clearly lower than for K and G dwarfs. Up to periods of several hundred days, this comparative deficit of Jupiter-mass planets for M dwarfs is statistically robust (Endl et al. 2006; Johnson et al. 2007; Bonfils et al. 2006). In contrast, planets less massive than $\sim$25 $M_\oplus$are significantly more frequent around M dwarfs (Bonfils et al. 2007, and 2009 in preparation).

The GJ 581 system is of particular interest, since two of its previously detected planets are located on the two edges of the habitable zone (HZ) of the M3V host star (von Bloh et al. 2007; Selsis et al. 2007). Despite the uncertainties in the exact location of the ``liquid water'' zone for planets more massive than the Earth, these detections demonstrate that super-Earths in the HZ can be detected by Doppler spectroscopy, for planets orbiting M dwarfs.

In the present paper, we report the detection of a fourth planet in theGJ 581 system, ``e'', with a minimum mass of 1.94 $M_\oplus$. We also correct the period of the outer planet GJ 581d, and find that the revised period implies that GJ 581d is in the habitable zone. In Sect. 2, we briefly recall the main characteristics of the host star GJ 581, and of its first three planets. An enlarged set of measurements (Sect. 3) allows a re-examination of the structure of the GJ581 planetary system, and the discovery of an additional very-low-mass planet (Sect. 4). Since Doppler spectroscopy, in particular for M dwarfs, can be confused by stellar surface inhomogeneities, we dedicate Sect. 5to discussing the magnetic activity ofGJ 581. In Sect. 6, we use dynamical stability considerations to derive upper limits to the planetary masses. To conclude, we briefly discuss how this system compares with theoretical predictions, and consider the prospects for finding even lower mass planets in the habitable zone of M dwarfs.

2 Stellar characteristics and planetary system

Our previous report of the first Neptune-mass planet on a 5.36-d orbit around GJ 581 (Bonfils et al. 2005b) describes extensively the properties of the host star. Here we therefore summarize those characteristics.

GJ 581 (HIP 74995, LHS 394) is an M 3 dwarf at a distance from the Sun of 6.3 pc. Its estimated luminosity is only 0.013 $L_\odot$.Bonfils et al. (2005a) find that GJ 581 is slightly metal-poor ([Fe/H] = -0.25), in contrast to the supersolar metallicities (Santos et al. 2004,2001) of most solar-type stars hosting giant planets. All indirect tracers (kinematics, stellar rotation, X-ray luminosity, and chromospheric activity) suggest a minimum age of 2 Gyr. The HARPS spectra exhibit weak Ca II H and K emission, in the lower quartile of stars with similar spectral types. This very weak chromospheric emission is indicative of a long stellar rotation period.

The first planet was easily detected after only 20 HARPS observations (Bonfils et al. 2005b), and the periodogram of the residuals of that first solution inferred a power excess around 13 days. As 30 additional observations were accumulated for a total time span of 1050 days, this coherent signal strengthened to strong significance, and a third planet appeared, although still with a significant false alarm probability (FAP): GJ 581c had $m\sin i = 5~M_\oplus$,_P_=12.9 days, and a $\sim$0.28% FAP, while the more distant GJ 581d had $m \sin i= 7.7~M_\oplus$,P = 83 days, and a $\sim$3% FAP. (Udry et al. 2007). Udry et al. (2007)therefore ended on a call for confirmation by additional radial velocity measurements, and based on the very-low magnetic activity ofGJ 581, suggested that the significant residuals might reflect additional planets in the system.

3 A new data set

Following a call for more data, we recorded a total of 119 HARPS measurements of GJ 581 (Table 1, only available electronically), which together span 1570 days (4.3 years). In parallel to this extension of the data set, the precision of all measurements, including the previously published ones, has been improved. The reference wavelengths of the thorium and argon lines of the calibration lamps were revised by Lovis & Pepe (2007), and Lovis (2009, in prep.) implemented a correction for the small effect of Th-Ar calibration lamp aging. The changing internal pressure of aging lamps shifts the wavelength of their atomic transition. The effect is now corrected down to a level of just a few 0.1 m/s, using the differential pressure sensitivity of the argon and thorium lines. The HARPS pipeline also corrects the radial velocities to the Solar System barycentric reference frame (Lindegren & Dravins 2003). That correction is computed at the weighted-mean time of the exposure, which is measured by a photometer that diverts a small fraction of the stellar flux. The effective time of the exposure is measured to a higher accuracy than 1%, which at the latitude of La Silla Observatory and for 900s-long exposures translates into a $\sim$ $0.3~{\rm m~s^{-1}}$worst case error in the Earth-motion correction. The combined photon noise, calibration error, and error in the mean time of the exposure result in a typical overall uncertainty of $1~{\rm {\rm m~s^{-1}}}$, with a full range of 0.7 to 2 ${\rm m~s^{-1}}$, reflecting the weather conditions. Finally, the pipeline corrects for the $0.20~{\rm {\rm m~s^{-1}~{\rm yr}^{-1}}}$perspective acceleration from the proper motion of the star (e.g., Kürster et al. 2003).

Table 1:Radial-velocity measurements, associated uncertainties and Ca IIH+K index for Gl 581.

4 GJ 581, a system with four low-mass planets

Our preferred method for planet searches in RV data uses a heuristic algorithm, which combines standard non-linear minimizations with genetic algorithms. It efficiently explores the large parameter space of multi-planet systems, and quickly converges toward the optimal solution. We attempted to fit models consisting of between zero and five planets in Keplerian orbits, and found that the data require a 4-planet model.

To verify the robustness of that ``black-box'' solution, and gain insight into its content, we also performed a step-by-step periodogram analysis. The top panel of Fig. 1displays the window function of our RV measurements, and unsurprisingly shows that the dominant periodicities of our sampling are 1 day and 1 year. The second panel is a floating-mean periodogram (Zechmeister & Kürster 2009; Gilliland & Baliunas 1987) of the RVs. It is dominated by a $\sim$5.36-day periodicity, which corresponds to GJ 581b. To estimate the FAP of that peak, we generated virtual data sets using boostrap resampling (Press et al. 1992) of the true measurements, and examined the peak power in the periodogram of each of these virtual datasets. None of these 10 000 trials had as high a peak as detected for the true 5.36-day signal, which therefore has a FAP below 0.01%. We then adjusted the Keplerian orbit corresponding to that starting period, and examined the periodogram of the residuals of this 1-planet solution (Fig. 1, 3rd panel). This second periodogran is dominated by a peak around $\sim$12.9 day with FAP< 0.01%, corresponding to the known planet GJ 581c. We then investigated the periodogram of the residuals of the corresponding 2-planet least squares adjustment. That new periodogram has a power excess in the 50-90 day range and around 3.15 days (as well as at 1-day aliases of these periods), both with FAPs< 0.01%. The broad power excess between 50 and 90 days divides into 3 separate peaks at $\sim$59, 67, and 82 days, which are 1-year aliases of each other ($1/67{-}1/365\sim1/82, 1/67+1/365\sim1/59)$. In Udry et al. (2007), we attributed an 82-day period toGJ 581 d. The corresponding periodogram peak now has markedly less power than the 67-day periodicity, and attempts to adjust orbits with that period (as well as with a 59-day period) produced significantly larger residuals. We conclude that the true period of GJ 581 d is around 67 days and that we had previously been confused by a 1-year alias of that period. As we discuss below, the revised period locates GJ 581 d in the habitable zone of GJ 581. Subtracting the corresponding 3-planet adjustment leaves no power excess in the 50-90 day range (5th panel), and markedly enhances the 3.15 day peak. The residuals of this 4-planet solution finally display no statistically significant peak (6th panel). The strongest peak is located close to one year ($P\sim384$ day) and has a $\sim$50% FAP.

[\begin{figure} \par\includegraphics[width=8cm,clip]{12172fg1.eps} \end{figure}](/articles/aa/full%5Fhtml/2009/43/aa12172-09/img12.png) Figure 1:Window function and periodograms of GJ581 radial velocity measurements (see text).
Open with DEXTER

Both the heuristic algorithm and the spectral analysis method identify the same 4 coherent signals in our data. Once this optimal approximate 4-planet solution had been identified, we performed a final global least squares adjustment, using the iterative Levenberg-Marquardt algorithm. Since we find insignificant eccentricities for planets ``b'' and ``c'', we fixed the eccentricities and longitude of periastron to zero for both planets. Table 2 lists the resulting orbital elements, and Fig. 2 displays the corresponding Keplerian orbits.

Compared to our Udry et al. (2007) study, we thus add a fourth planet and very significantly revise the characteristics of the outer planet (GJ 581d), and in particular its period. Planets _b_to e have respective minimum masses of $m \sin i \sim 15.7$, 5.4, 7.1, and 1.9 $M_\oplus$.GJ 581e is the lowest mass exoplanet detected so far around a main-sequence star. None of the period ratios of the system comes close to any simple ratio of integers, with e.g., $P_{\rm b}$/ $P_{\rm e} = 1.70$, $P_{\rm c}/P_{\rm b} = 2.4$, and$P_{\rm d}$/$P_{\rm c} = 5.17$. Non-resonant hierarchies have similarly been found for all multi-planet systems containing only low mass planets, such as HD 40307 (Mayor et al. 2009) and HD 69830 (Lovis et al. 2006).

Table 2:Fitted orbital solution for the GJ 581 planetary system: 4 Keplerians.

5 Planets rather than spots

Unfortunately, coherent Doppler shifts do not always correspond to planets. Inhomogeneities of the stellar surface such as spots, plages, flares, or convective patterns can break the even distribution between the blue- and red-shifted halves of a rotating star and induce apparent radial velocity shifts. This can introduce noise into the RV measurements, which is usually referred to as ``jitter''. In some cases, the surface pattern can be stable over long time spans (from days to years), and the ``jitter'' noise is then correlated. In these cases, it does not average out as white noise would, and instead builds up a coherent signal. This apparent Doppler shift can easily mimic a Keplerian orbit (e.g., Queloz et al. 2001; Bonfils et al. 2007).

[\begin{figure} \par\includegraphics[width=7.75cm,clip]{12172fg2.eps} \end{figure}](/articles/aa/full%5Fhtml/2009/43/aa12172-09/img26.png) Figure 2:Radial velocity curves for planets e, b, c and d, from top to bottom. The lowest panel displays the residual to the four-planets keplerian fit.
Open with DEXTER

M dwarfs tend to have stronger magnetic activity than solar-like stars, and possible false detections caused by coherent jitter therefore need even closer scrutiny. We indeed encountered several such cases in our M star survey, including GJ 674 (Bonfils et al. 2007) and GJ 176 (Forveille et al. 2009). In both cases, photometric monitoring and spectroscopic activity diagnostics (measuring the Ca IIH&K and H$_\alpha$lines) demonstrated that an apparently Keplerian signal was caused by a spot modulation of the absorption line profiles. These two examples serve as reminders that a periodic Doppler signal does not always reflect the presence of a planet, and that additional checks are needed before robust conclusions can be made.

In our previous papers about the GJ 581 system, we could exclude any such confusion for the first three planets. The $\sim$ $13~{\rm m~s^{-1}}$radial velocity amplitude related toGJ 581 b (_P_=5.36 day - Bonfils et al. 2005b) is sufficiently large for a spot signal of that magnitude to have induced measurable distortions in the rotational profile. The lack of any correlation between the bisector span (a first order measurement of spectral line asymmetry) and the measured RVs, excludes that scenario. The strength of the Ca IIH&K emission lines in addition classifies GJ 581 in the lowest $\sim$5% magnetic activity bin of our M dwarf sample. This implies a much longer rotation period than the observed 5.4 days, as confirmed by the very low projected rotational velocity $v \sin i < 1~{\rm km~s^{-1}}$). This velocity would need to correspond to an unprobably low inclination i to be compatible with a 5 day rotation period. Additional measurements derived since Bonfils et al. (2005b) have strengthened these conclusions even further.

In Udry et al. (2007), we used distinct lines of arguments for the 12.9 and the 83 day (now 67 day) signals. The arguments developed above for GJ 581b hold for the 12.9-day period, which is also too short to be compatible with the rotation period. They also apply for GJ 581e at 3.1 days. On the other hand, an 83 day (or 67 day) period would be perfectly consistent with the rotation period of an M dwarf with an activity level ofGJ 581. For these long periods, the bisector span also loses its diagnostic power, because the rotational profile becomes much too narrow to be resolved with HARPS. We therefore require to take a different approach, and we used a relation by Saar & Donahue (1997) to estimate the minimum filling factor of a spot that could produce the observed radial velocity signal. The relation that we use

\begin{displaymath}K_{\rm s} \sim 6.5 \times f_{\rm s}^{0.9} \times v \sin i \ \ \ \ \left[{\rm m~s^{-1}}\right] \end{displaymath} (1)

connects the semi-amplitude $K_{\rm s}$of the radial-velocity signal induced by an equatorial spot with its surface filling factor, $f_{\rm s}$%, and the projected rotation velocity, $v \sin i$

. This has been verified by numerous examples, including e.g., GJ 674 and GJ 176. Applying this relation to a $R_\star=0.31~{R_\odot}$stellar radius, a P

=83 days period, and a $K\sim2.5~{\rm m~s^{-1}}$

semi-amplitude, Eq. (1) readily shows that the GJ 581 d signal would have needed a $f_{\rm s} \sim 2.2\%$filling factor to be explained by a spot. Such a large spot is firmly excluded by the sub-% photometric stability of GJ 581.

Our revised 66.7 day period for GJ 581 d slightly weakens that argument, since a shorter period signal can be induced by a smaller spot. The revision however is minor, as the 67-day signal still requires a $f_{\rm s} = 1.7\%$spot, which is easily excluded by the improved photometric measurements that have become available since Udry et al. (2007). The strongest photometric stability constraint is provided by 6 weeks of continuous monitoring of GJ 581 by the MOST satellite. Preliminary results, presented by Matthews at the 2007 Michelson Summer School[[*]](/articles/aa/full%5Fhtml/2009/43/aa12172-09/aa12172-09.html#foot1179), indicate that the peak-to-peak photometric variability during those 6 weeks was less than 5 mmag. The 66.7-day HARPS radial velocity signal is therefore incompatible with rotational modulation of spots or granulation patches.

The variations in the chromospheric Ca IIH+K line (Table 1) provide another diagnostic of magnetic activity, which proved very valuable in the case of GJ 674 (Bonfils et al. 2007). We parametrize this chromospheric emission by

\begin{displaymath}{\rm Index} = \frac{H+K}{R+B}, \end{displaymath} (2)

where H (resp. K) is the flux measured in the Ca IIH (resp. K) line, and where R and B measure the flux in pseudo-continuum bands on both sides of the lines. Figure 3displays our Ca II H+K measurements as a function of time (top panel), and the corresponding periodogram (lower panel). The figure implies that there are long-term variations, over timescales of several thousands of days, which are much too long to be the rotational period. If those variation are real, they would be more likely to be an analog of the 11 year solar magnetic cycle. Magnetic cycles can potentially affect apparent stellar velocities, since they can change the balance between ascending (blue-shifted) and receding (red-shifted) convective elements in the atmosphere, but we see no evidence of long-term RV variations. The other significant structure in the chromospheric signal is a broad power excess between 80 and 120 days, which is far more likely to reflect the rotation period. This peak overlapped our incorrect initial estimate of the period of GJ 581 d, but as discussed above photometry limits the amplitude of any velocity signal from spots at that period to at most a quarter of the observed value. The revised period provides conclusive evidence of that concern, since the periodogram of the Ca II H+K indices has no power excess around 67 days. The Ca II H + K periodogram finally has a weak peak at $\sim$15 days, which does not match any of the observed radial-velocity periods.

[\begin{figure} \par\includegraphics[width=8.8cm]{12172fg3.eps} \end{figure}](/articles/aa/full%5Fhtml/2009/43/aa12172-09/img39.png) Figure 3:The Ca II H+K index as function of Julian date ( upper panel) and its periodogram (lower panel).
Open with DEXTER

To summarize, none of the 4 periodic signal detected in our RV measurements is compatible with rotational modulation of stellar surface patterns. Four low-mass exoplanets orbiting the star remain, by far, the most likely interpretation.

6 Dynamical evolution

We reinvestigate the dynamics of the revised GJ 581 system, using the same techniques as in Beust et al. (2008, B08). With the SyMBA symplectic _N_-body code (Duncan et al. 1998), we carry 5-body integrations over 0.1 Gyr, with a $2\times10^{-4}~$yr timestep. Starting from the nominal solutions (Table 2), which assumes initially circular orbits for Gl581e and Gl581b (the two innermost planets), we perform those integrations for inclinations of 90 to 10 degrees, assuming full coplanarity of all orbits.

Our nominal simulation uses the periastron parameters derived in the Keplerian fit, listed in Table 1. To probe the sensitivity of the simulations to these parameters, we performed additional integrations with different starting values for the periastron longitudes. Unsurprisingly, since the orbits are non-resonant, the results did not differ significantly.

Figure 4shows the maximum variations in the semi-major axes and eccentricities during those integrations. As in B08, we note that the semi-major variations gradually increase as the inclination decreases. This is clear indication that lower inclinations, i.e., higher planetary masses, produce a less stable system. Compared to our previous study, the addition of GJ 581e drastically tightens the minimum inclination constraint. This low mass planet has little effect on the stability of the more massive GJ 581b, GJ 581c, and GJ 581d, but for inclinations $5\le i \le 30\hbox{$^\circ$ }$,GJ 581e is ejected out of the system on timescales of at most a few Myrs.

The nominal integrations assume circular initial orbits for GJ 581e and GJ581b. Relaxing this constraint and setting the initial eccentricity of planet e to 0.1, we found the system to be even less stable: GJ 581e was then ejected for inclinations up to $i=40\hbox{$^\circ$ }$.

Unlike most multi-planet radial-velocity systems, the nominal least squares adjustment to the radial velocities is dynamically stable for sufficiently high inclinations (and co-planar orbits). Restricting the solutions to dynamically stable systems only implies that $i\ga40\hbox{$^\circ$ }$and therefore that the mass of each planet is no more than $\sim$1.6 times its minimum mass. For GJ 581e, b, c, and d, those upper limits are 3.1, 30.4, 10.4, and 13.8 $M_\oplus$.

7 Summary and discussion

The large number of precise HARPS measurements has now revealed four planets in a non resonant-configuration. The FAPs of all four planetary signals are below 0.01%. The orbital elements of planets GJ581b (_P_=5.36d - Bonfils et al. 2005a) and GJ581c (_P_=12.9d - Udry et al. 2007) are fully confirmed, and slightly refined. The 82 day period reported by Udry et al. (2007) for GJ 581d, on the other hand, was a one-year alias of the true 66.8 day period of that planet. Our extended time span and much larger number of measurements correct this confusion. Photometry, and monitoring of the Ca IIH+K activity indicator over the four years of our spectroscopic measurements, indicate that the signal at 66 days does not originate in stellar activity. At the present time, based on the robustly measured radial velocity signal and an estimate of the worst-case effect of stellar activity, a planet with a period of 66.8 days is by far the most probable interpretation.

[\begin{figure} \par\makebox[18cm]{\includegraphics[angle=270,width=8.8cm]{12172f4a.ps}\hfil \includegraphics[angle=270,width=8.8cm]{12172f4b.ps} }\end{figure}](/articles/aa/full%5Fhtml/2009/43/aa12172-09/img44.png) Figure 4:Stability of various configurations of the four-planet Gl 581 system, after 108 years integration.
Open with DEXTER

The new orbital elements for the outer planet GJ 581d correspond to a 0.22 AU semi-major axis. A planet with an eccentric orbit corresponds to a higher orbit-averaged flux than a planet on a circular orbit of the same semi-major axis (Williams & Pollard 2002), and GJ 581 d receives the same average flux as a planet on a a = 0.21 AU circular orbit. The Udry et al. (2007) orbital elements of GJ 581d, with a 0.25 AU semi-major axis, inferred that the planet is close to the outer edge of the habitable zone, but most likely just slightly outside (e.g., Selsis et al. 2007). Our revised parameters place GJ 581dfirmly in the habitable zone.

By analyzing 119 precise radial velocity measurements, we have discovered a fourth planet in the GJ 581 system. Its minimum mass, $m \sin i = 1.9~M_\oplus$, is by far the lowest amongst planets detected around main-sequence stars. This low-mass planet is almost certainly rocky, and its equilibrium temperature is too high to allow a substantial atmosphere to exist.

In our HARPS M-star sample, GJ 581 has the largest number of measurements, which evidently helped us to detect such a low mass planet in a four planet system. The detection however does not stretch the limits of Doppler spectroscopy, and for the same period a planet only half as massive would have been detectable in our dataset. Our HARPS measurements set an upper limit of 1.2 m/s to the stellar jitter of GJ 581, based on a quadratic subtraction of just the mean photon noise from the observed dispersion of the radial-velocity residuals to the 4-planet solution. Some small fraction of this value could still be instrumental, but a 1.2 m/s ``jitter'' allows the detection of a 2 $M_\oplus$planet with a 25 day orbital period. This planet would be earth-type and in the middle of the habitable zone of its M dwarf.

The stellar ``jitter'' has different physical components, of independent timescales. For more massive stars and shorter integration times, acoustic modes represent a major contribution, but the short acoustic timescales in this small star and the 15-mn integration ensure that they have been averaged out to below a few 0.1 m/s. The 1.2 m/s jitter therefore reflects mostly longer timescale phenomena, such as granulation noise and changes in stellar surface anisotropies during a magnetic cycle. We are developing optimized measurement strategies to average out the granulation noise contribution to the stellar jitter, but HARPS already demonstrates that the global jitter (including instrumental errors, but not photon noise) of reasonably chromospherically quiet M dwarfs is at most on the order of 1 m/s.

This new low-mass planet adds to the small set of Doppler-detected super-Earths, planets with a minimum mass of less than 10 $M_\oplus$. This planetary mass domain, as well as the Neptune-mass range (<30 $M_\oplus$) has been largely populated by the HARPS surveys. A few statistical properties have already emerged from these early discoveries (Mayor & Udry 2008):

-

The full distribution of planetary masses is bimodal, with distinct peaks corresponding to gaseous giants and super-Earths. Despite the observational bias against low mass planets, the distribution below $\sim$30 $M_\oplus$rises towards super-Earth planets (cf. Fig. 7 in Mayor & Udry 2008).

-

The majority of super-Earths and Neptunes are found in multiplanetary systems. Four of the 6 planetary systems with a known super-Earth, (GJ 876, HD 40307, HD 7924,GJ 176, GJ 581, HD 181433) have multiple planets. Two of these 4 multi-planet systems associate one super-Earth with one or two gaseous giant planets (GJ 876,HD 181433), and the other 2 have several low-mass planets on non-resonant orbits (HD 40307, GJ 581). GJ 176 andHD 7924have only one detected planet, although they could obviously have more which have yet to be detected. GJ 176 is too magnetically active for a low-mass planet to be easily detected, but the periodogram of HD 7924 does indeed provide evidence of additional planets (Howard et al. 2009).

-

In contrast to gaseous giants, low-mass planets do not seem to be more frequent around metal-rich host stars (Udry et al. 2006).

-

In a preliminary analysis, we have detected low mass close-in planets (P < 100 d and $m \sin i < 30~M_\oplus$) around 30% of solar-type stars in the HARPS high precision survey (Lovis et al. 2009, in prep.).

Systems with several low mass close-in planets provide interesting constraints on models of planetary formation. Three systems are particularly noteworthy: HD 69830 (3 planets),HD 40307 (3 planets) and GJ 581(4 planets).

Terquem & Papaloizou (2007), for instance, studied the migration of cores and terrestrial planets induced by their interaction with the protoplanetary disk, and suggested that ``if hot super-Earths or Neptunes form by mergers of inwardly migrating cores, then such planets are most likely not isolated. We would expect to always find at least one, more likely a few, companions on close and often near-commensurable orbits''. The high observed fraction of multi-planet systems matches that prediction, but the observed periods are quite far from near-commensurability.

Raymond et al. (2008) discuss observable consequences of planet formation models in systems with close-in terrestrial planets, and specifically address the case of GJ 581. In addition to system architecture, they consider information from planet composition, which unfortunately is unavailable for the GJ 581 system. These authors find two formation mechanisms consistent with the available data: in situ formation, and formation at a larger distance followed by type I migration. In situ formation, however, would correspond to a massive disk, 17-50 times above the minimum-mass disk. On the other hand, type I migration from a larger distance is expected to leave the close-in low mass planets in a MMR (mean motion resonance), at odds with the characteristics of the GJ 581 system.

Kennedy & Kenyon (2008) study hot super-Earth formation as a function of stellar mass and suggest that, with migration, short period low-mass planets are most likely to form around low-mass stars. Above approximately 1 solar-mass, the minimum protoplanet mass for migration to close-in orbits is above 10 Earth-masses and no hot super-Earths can form. Our searches for low mass planets with the HARPS spectrograph will test that prediction.

Zhou et al. (2009, in prep.) study the birth of multiple super-Earths by sequential accretion. Focusing on the formation of theHD 40307 triple super-Earths system, they found evidence that massive embryos can assemble in a gaseous protostellar disk. The final states of their numerical simulations favor models with a reduction factor of ten for type I migration. These simulations quite succesfully reproduce the global structure of the triple super-Earth system, including non-resonant period ratios.

Ogihara & Ida (2009) provides an illuminating discussion of formation mechanisms of low-mass multiple planets orbiting M stars. It shows, in particular, that the final orbital configuration could be quite sensitive to migration timescale. The water fraction in low mass planets close to the habitable zone of M dwarfs is probably very large, making GJ 581d a serious candidate ocean planet (Léger et al. 2004).

Our HARPS surveys search for low-mass close-in planets around stars with masses from 0.3 to over 1 solar-mass. The numerous multiplanetary systems that they will detect will provide important constraints on formation scenarios.

Acknowledgements

The authors thank the observers of the other HARPS GTO sub-programmes who helped measure GJ 581. We are grateful to the staff of La Silla Observatory for their contribution to the success of the HARPS project. We are grateful to our colleagues Ji-Lin Zhou, Douglas Lin, Su Wang and Katherine Kretke for having made their study of multiple super-Earth formation available in advance of publication. We wish to thank the Programme National de Planétologie (INSU-PNP) and the Swiss National Science Foundation for their continuous support of our planet-search programs. X.B. acknowledges support from the Fundação para a Ciência e a Tecnologia (Portugal) in the form of a fellowship (reference SFRH/BPD/21710/2005) and a program (reference PTDC/CTE-AST/72685/2006), as well from the Gulbenkian Foundation for funding through the ``Programa de Estímulo Investigação''. N.C.S. would like to thank the support from Fundação para a Ciência e a Tecnologia,Portugal, in the form of a grant (references POCI/CTE-AST/56453/2004 and PPCDT/CTE-AST/56453/2004), and through program Ciência 2007 (C2007-CAUP-FCT/136/2006).

References

Footnotes

... planets[[*]](/articles/aa/full%5Fhtml/2009/43/aa12172-09/aa12172-09.html#tex2html1)

Based on observations made with the HARPS instrument on the ESO 3.6 m telescope under the GTO and LP programs 072.C-0488 and 183.C-0437 at Cerro La Silla (Chile). Our radial-velocity, photometric and Ca II H+K index time series (Table 1) are 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/507/487

... School[[*]](/articles/aa/full%5Fhtml/2009/43/aa12172-09/aa12172-09.html#tex2html8)

http://nexsci.caltech.edu/workshop/2007/Matthews.pdf

Copyright ESO 2009

All Tables

Table 1: Radial-velocity measurements, associated uncertainties and Ca II H+K index for Gl 581.

Table 2: Fitted orbital solution for the GJ 581 planetary system: 4 Keplerians.