A&A 443, L15-L18 (2005)
DOI: 10.1051/0004-6361:200500193
X. Bonfils 1,2 - T. Forveille 3,1 - X. Delfosse 1 - S. Udry 2 - M. Mayor 2 - C. Perrier 1 - F. Bouchy 4 - F. Pepe 2 - D. Queloz 2 - J.-L. Bertaux 5
1 - Laboratoire d'Astrophysique, Observatoire de Grenoble, BP 53,
38041 Grenoble Cedex 9, France
2 -
Observatoire de Genève, 51 Ch. des Maillettes, 1290 Sauverny,
Switzerland
3 -
Canada-France-Hawaii Telescope Corporation, 65-1238 Mamalahoa Highway,
Kamuela, HI96743, Hawaii, USA
4 -
Laboratoire d'Astrophysique de Marseille, Traverse du Siphon,
13013 Marseille, France
5 -
Service d'Aéronomie du CNRS, BP 3, 91371 Verrières-le-Buisson, France
Received 22 August 2005 / Accepted 24 September 2005
Abstract
We report the discovery of a Neptune-mass planet around Gl 581
(M3V, M = 0.31 M
), based on precise Doppler measurements with the HARPS
spectrograph at La Silla Observatory. The radial velocities reveal a
circular orbit of period P = 5.366 days and semi-amplitude
K1 = 13.2 m s-1. The resulting minimum mass of the planet
(
)
is only 0.052
= 0.97
= 16.6
making Gl 581b one of the lightest extra-solar planet known to date.
The Gl 581 planetary system is only the third centered on an M dwarf,
joining the Gl 876 three-planet system and the lone planet around Gl 436. Its
discovery reinforces the emerging tendency of such planets to be of low
mass, and found at short orbital periods. The statistical properties of
the planets orbiting M dwarfs do not seem to match a simple mass scaling
of their counterparts around solar-type stars.
Key words: stars: individual: Gl 581 - stars: planetary systems - stars: late-type - techniques: radial-velocity
Over 150 planets have been found orbiting main sequence stars other than the Sun, in about 140 planetary systems of which 18 have multiple planets (http://vo.obspm.fr/exoplanetes/encyclo/). These extra-solar planets are a very diverse class: their mass ranges between half the mass of Neptune and 15 times the mass of Jupiter, some have large eccentricities when others have nearly circular orbits, their periods range from slightly over a day (Konacki et al. 2003, OGLE-TR-56) to over a decade (Marcy et al. 2002, 55 Cnc). The multiple systems range from strongly resonant to fully hierarchical (Marcy et al. 2002; Rivera et al. 2005). This diversity demonstrates that our own solar system represents but one possible outcome of the planetary formation and evolution processes, and apparently not even a very common one.
The statistical properties of these exoplanets provide crucial clues to their formation mechanism. As perhaps the most dramatic example, the seminal detection of the 51 Peg planet in a 4-days orbit (Mayor & Queloz 1995) immediately forced theoreticians to recognize the critical importance of orbital migration (Ward 1997; Lin et al. 1996). The correlation between the occurence of Jupiter-mass planet and the high metallicity of the host stars (Gonzalez 1997; Santos et al. 2004a,2001) is another example. It is thought to reflect the controlling role of the condensate mass in the protoplanetary disk, but it has taken longer to converge towards that consensus.
To date, all but 2 of these 140 planetary systems orbit solar-type stars.
In part, this no doubts reflects a bias of most planet-search programmes
towards the relatively bright F to K main sequence stars, and away from their
fainter M-type counterparts (M < 0.6 M
). Nonetheless, several teams
(Wright et al. 2004; Delfosse et al. 1998b; Endl et al. 2003)
collectively monitor over 200 M dwarfs with
sufficient precision to detect a Jupiter-mass planet out to at least 2 AU.
These efforts have up to now identified the 3-planet system around Gl 876
(Marcy et al. 1998,2001; Rivera et al. 2005; Delfosse et al. 1998b),
and the single-planet Gl 436 system (Butler et al. 2004). Of
these 4 planets, 2 are in the Neptune-mass class, leaving only two of
the Gl 876 planets with approximately Jupiter-mass. By constrast
5% of solar-type stars have Jupiter-mass planets
(Marcy et al. 2000), and the comparative deficit for the
M dwarfs is therefore statistically robust (Naef et al. 2005; Butler et al. 2004).
An open question, though, is whether M dwarfs genuinely have fewer planets, or whether their planets are just as abundant, but not quite as massive. Addressing this question needs higher precision, or more measurements, than the radial-velocity surveys have achieved to date. To help answering this question, we are using the HARPS spectrograph for a high-precision survey of more than 100 nearby M dwarfs. We present here its first detection, a Neptune-mass planet around Gl 581.
Gl 581 (HIP 74995, LHS 394) is an M3 dwarf (Hawley et al. 1997)
with a distance to the Sun of 6.3 pc (
,
ESA 1997).
Its photometry (
,
B-V=1.60, Mermilliod et al. 1997;
,
Leggett 1992) and the parallax together result
in absolute magnitudes of
,
and
.
From its absolute V magnitude and the 2.08 V-band bolometric correction of
Delfosse et al. (1998a), the luminosity of Gl 581 is 0.013 L
.
The Delfosse et al. (2000) K-band mass-luminosity relation, which
has much lower intrinsic dispersion than the equivalent V-band relation,
gives a
M
mass,
and the Chabrier & Baraffe (2000) theoretical Mass-Radius relation
then a radius of 0.29 R
.
Interestingly, Bonfils et al. (2005) find Gl 581
slightly metal-poor (
), in contrast to most planet-host stars
having supersolar metallicities.
The age of Gl 581 can be estimated from its kinematic characteristics, its
magnetic activity, and its metallicity, all of which point towards a
moderate to older age. Leggett (1992) find that its UVW galactic
velocities are intermediate between those typical of the young and old
galactic disk, and Delfosse et al. (1998a) find very low X-ray emission
(
)
and a 2.1 km s-1 upper limit on the projected
rotation velocity,
.
Gl 581 has been classified as a variable star (HO Lib), but its variability
(Weis 1994) is only marginally significant. If real it would be on a time
scale of several year, with short-term variability being at most
0.006 mag.
The HARPS spectra show weak Ca II H and K emission, in
the lower quartile of stars with similar spectral types (Fig. 1).
As mentioned above, Gl 581 also has a subsolar metallicity.
Altogether, these properties suggest that it is at least 2 Gyr old, and
they ensure that the radial velocity "jitter'' from magnetic activity
must be minimal.
![]() |
Figure 1:
HARPS spectra of the Ca II H
(
![]() ![]() |
Open with DEXTER |
For the V=10.5 Gl 581 we use 15 mn exposures, and the median S/N ratio
of our 20 spectra is 40 per pixel at 550 nm. The radial
velocities (Table 2, only available electronically) were
obtained with the standard HARPS reduction pipeline, based on the
cross-correlation with a stellar template and the precise nightly
wavelength calibration with ThAr spectra (Baranne et al. 1996).
They have a median internal error of only 1.3 m s-1, which includes both
the nightly zero-point calibration uncertainty (0.8 m s-1)
and the photon noise, computed from the full Doppler
information content of the spectra (Bouchy et al. 2001).
The computed velocities exhibit an rms dispersion of 10 m s-1, much above their internal errors and also considerably more than we observe for stars with higher chromospheric activity. Of the three comparison stars with stronger chromospheric emission in Fig. 1, Gl 752 and Gl 832 have enough HARPS radial velocities to measure rms dispersions of 1.2 m s-1 (from 10 measurements) and 1.8 m s-1 (from 19 measurements). Rivera et al. (2005) report an rms dispersion of 3.42 m s-1 for the the third, Gl 908, dominated by their measurement noise.
![]() |
Figure 2: Upper panel: phased radial velocities for Gl 581. Lower panel: residuals around the fitted solution versus time. |
Open with DEXTER |
As demonstrated by Fig. 2, a circular orbit of period
P = 5.366 d and semi-amplitude K = 13.2 m s-1 is an excellent
fit to these velocities. Attempts at adjusting elliptical orbits resulted in
non-significant eccentricities. We therefore adopt a circular orbit. Its
parameters are listed in Table 1 and lead to a minimum mass
()
for the planet of only
0.052
= 0.97
= 16.6
.
The weighted rms of the residuals around the fit is 2.5 m s-1, and
twice the internal errors of the measurements. More data points are needed to
establish whether this extra dispersion is intrinsic to the star, or
whether it could denote the presence of a third body in the system. The
very stable shape of the bisector (Fig 3) gives some weight
to the later hypothesis. An overall translation dominates the evolution
of the spectral profile, leaving no doubt on the keplerian origin of the
radial-velocity variations.
![]() |
Figure 3: Bissectors of the HARPS correlation profiles for Gl 581. The shape of the bissector curve is independent of its position, eliminating the possibility that stellar spot cause the radial-velocity variations. |
Open with DEXTER |
We also have available 11 ELODIE spectra, which extend the measurement
baseline to 9 years, albeit with a 5 years gap between mid-2000 and
mid-2005. Their 17 m s-1 median error bars are too
large to reveal the K=13 m s-1 planetary signal.
They exhibit an rms dispersion of 23 m s-1 and lack
any obvious long term trend, demonstrating that the Gl 581 system
contains no Jupiter-mass planet with any period shorter than
10 years.
Table 1: Orbital and physical parameters.
The semi-major axis of the planetary orbit is only 0.042 AU or 9 solar
radii, similar to most close-in planets around solar-type stars. In the
natural length unit of its central star however, this amounts to 31 radii
of Gl 581. The geometric transit probability is thus only 3%, and
significantly less than the 10% typical of close-in planets
around solar-type stars. If transits do occur on the other hand, the
planet will cover a larger fraction of its smaller star. For a constant
planetary radius, transits would thus be correspondingly deeper and more
easily detected. At this
radius and given the 0.013 L
luminosity of the star, the expected
temperature of the planetary surface is
420 K, with large
uncertainties from the unknown albedo and energy transport. Even with
conservative error bars though, this temperature is compatible with either
a rocky planet or a gas giant, and evaporation will be negligibly small
in either configuration.
The detection of Gl 581b brings the inventory of M-dwarf which harbor
planetary systems to 3, and the number of their planets to 5. While
admittedly still very small, these samples allow us an initial peek at their
properties as a population, compared to the much more numerous planets
known around solar-type stars. One immediate observation is that none of
the three stars is metal-rich, with Gl 876, Gl 436 and Gl 581 having
metallicities of respectively
,
+0.02 and -0.25 dex
(Bonfils et al. 2005). This contrasts with the median metallicity for
solar-type stars surrounded by planets,
(Santos et al. 2005), though the significance of the difference
is obviously still modest.
As discussed in the introduction, various groups monitor over 200 M dwarfs
with 3-15 m s-1 precision,
sufficient to easily detect the >40 m s-1 reflex motion of a
0.3
star orbited by a Jupiter-mass planet out to 2 AU.
That these efforts have to date found only 5 planets, of which only
Gl 876b and c have approximately Jovian masses, demonstrates that there
are much fewer
planets around M dwarfs than the
5% (Naef et al. 2005; Marcy et al. 2000) found around solar-type stars.
The 5 planets include no hot-Jupiter, but with only
1%
solar-type star orbited by such a planet the significance of that fact is
still modest. Three of the 5 on the other hand are hot-Neptunes (Gl 436b,
;
Gl 876d, 0.44
;
Gl 581b, 0.99
),
as many as currently known around all solar-type stars. This matches the
theoretical model of Ida & Lin (2005):
the mass-distribution of close-in planets has two peaks centered
at about the masses of Jupiter and Neptune, with the former preferentially
populated around G-dwarfs and the latter around M dwarfs, reflecting
how much matter remains available in the disk for accretion during
the inward migration of the planet. Other theoreticians however take the
view that many hot-Neptunes are actually evaporated hot-Jupiters
(Baraffe et al. 2005). Better statistics on M-dwarf planets will
help determine which of these mechanisms dominate.
A final striking characteristic of the current M-dwarf planets is
that none has a period longer than the 2 months of Gl 876b.
By contrast, 66% of the 164 planets known around solar-type stars have
orbital periods above 100 days and their distribution is even observed to increase with period (Udry et al. 2003).
With 5 planets known around M dwarfs, the
probability that the long-period deficit amongst M-dwarf occurs by chance
is thus less than (0.34)
.
The sensitivity of
Doppler searches does degrade for longer periods however, and planets of
M dwarfs have often been found close to the sensitivity floor of their
respective discovery surveys, making detection biases potentially important.
A full account is not currently possible from published information (our
own high precision M-dwarf survey is still too recent to be very useful in
this respect), but for periods shorter than the observing interval the
sensitivity degrades only slowly (
P-1/3). Some of the
M-dwarf Doppler surveys have been observing for long enough
(Bonfils et al. 2004; Butler et al. 2004)
that detection biases make an unlikely full explanation for the
lack of long period planets. That deficit therefore has to be at least
in part intrinsic, perhaps reflecting smaller protoplanetary disks around
the lower mass M dwarfs.
Acknowledgements
We thank our technical and scientific collaborators of the HARPS Consortium, ESO Head Quarter, and ESO La Silla, who have contributed with passion and competence to the success of the HARPS project. This study benefited from the support of the HPRN-CT-2002-00308 European programme. We are also grateful to Damien Ségransan who contributed additional ELODIE observations, and to Jean-Christophe Leyder and collaborators for using some of their own observing time to obtain critical confirmation observations.
Table 2: Radial-velocity measurements and error bars for Gl 581. All values are relative to the solar system barycenter.