L'OBSERVATION D'EXOPLANETES
PAR LA METHODE DES VITESSES RADIALES

Les cas HD189733 b, HD 195019 b, 51 Peg b et tau Bootis b
 
EXTRASOLAR PLANET OBSERVATIONS
BY RADIAL VELOCITY METHOD

The situation of HD189733 b, HD 195019 b, 51 Peg b and tau Bootis b 


 

Cette page présente un compte rendu d'observation de 4 exoplanètes déjà découvertes par la communauté professionnelle, mais dans les cas HD189733b, HD195019b et 51 Peg b, pour la première fois observée avec la technique des vitesses radiales par des astronomes non professionnels. A cette liste s'ajoute la planète extrasolaire de l'étoile Tau Bootis, découverte par  P. Butler and G. Marcy en 1997, et déja été détectée par des amateurs (Thomas Kaye & all, 1996).

L'observation de Tau Bootis b a aussi fait l'objet d'un compte rendu séparé, consultable ici. On y décrit le principe de détection employé, qui n'est que brièvement rappelé ici. Les résultats pour HD189733, HD195019 et 51 Peg sont en revanche originaux.

Le spectrographe utilisé est le modèle eShel, commercialisé par la société Shelyak. Il s'agit d'un spectrographe échelle ayant un pouvoir de résolution de R=11000, pouvant capturer l'ensemble d'un spectre stellaire de 430 à 700 nm en 20 ordres. Pour une revue de quelques-unes des possibilités de ce spectrographe, cliquer ici.

Pour cette observation de planètes extrasolaires, seule la partie du spectre comprise entre 440 et 640 nanomètres a été conservée afin de tenter de réduire les erreurs de mesures causées par les éventuelles variations temporelles des principales raies telluriques. L'ensemble des logiciels nécessaire à la réduction des données et à l'analyse (prétraitement, calcul de la vitesse radiale par corrélation, périodogramme) on été écrits spécialement pour cette observation en langage C#.

Une lampe  thorium-argon est employée pour l'étalonnage spectral. Des lampes tungstène et LED servent par ailleurs à trouver la fonction de blaze et la position des ordres dans le spectre échelle.

Le télescope et le spectrographe sont reliés par une fibre optique de 50 microns de diamètre et de 15 mètres de long. Une caméra vidéo Watek 120N sert à l'autoguidage via le logiciel Audela. La caméra principale du spectrographe est un modèle QSI532 (CCD Kodak KAF-3200ME) piloté, elle aussi par le logiciel Audela.

Il a fallu surmonter deux difficultés importantes :

Les mesures de vitesses radiales se sont déroulées entre mars 2009 et aout 2009. Les observations à Castanet-Tolosan ont été réalisée par Christian Buil, et celles au Pic du Midi par Christian Buil, Valérie Desnoux, Michel Pujol et Olivier Thizy. Un grand merci au comité des programmes du télescope de 60 cm de nous avoir octroyer du temps de télescope pour cette étude ainsi qu'à l'équipe technique de l'Association T60 qui nous a permis de disposer un télescope performant et opérationnel.

This web page is a summary of the observation of four exoplanets already discovered by professional community. In the case of HD189733-b, HD195019-b and 51 Peg-b, this is the first Radial Velocity observation by amateur astronomers using a commercial echelle spectrograph. In addition, tau Bootis-b has also been observed but was already detected by amateurs using a home built spectrograph (see Thomas Kaye & all, 1996).

Results on tau Bootis-b have been previously published on this web; check out tau Boo-b results. We describe there the general principle of the method to detect exoplanet, which are briefly described here. Results on HD189733-b, HD195019-b et 51 Peg-b are genuine.

Spectrograph used is an off-the-shelf eShel model which is commercially available from Shelyak Instruments. It is an optical fibre fed echelle spectrograph with a resolution power of R=11000, covering spectral domain 430nm-700nm in 20 orders. For a more detailed review of this spectrograph and some other applications, check our eShel review.

For this work on exoplanets, only the spectral domain 440nm-640nm has been used to minimize the impact of telluric lines temporal evolution (absorption lines from our Earth atmosphere). All software used for this specific project (preprocessing, Radial Velocity calculation, periodigram...) have been specially writen in C# language.

A Thorium-Argon lamp is provided with the eShel spectrograph and has been used for wavelength calibration. Tungsten and LED have been used for the blaze function and to get the geometry of each echelle orders.

Telescope and spectrograph are linked by a 50 microns (core diameter) optical fibre (15 meters long). A Watec 120N video camera has been used for autoguiding using Audela software. Main camera on the spectrograph to capture spectra is a QSI532 (CCD Kodak KAF-3200ME) also managed by Audela software.

We had to overcome two major difficulties :

Radial Velocities have been measured between March and August 2009. Castanet-Tolosan observations have been done by christian Buil. Pic du Midi observations have been done by Christian Buil, Valérie Desnoux, Michel Pujol et Olivier Thizy. A big Thank You to the T60 (24inch telescope) program committee for allowing us time for this study as well as the technical team from Association T60 for providing a powerful and operational telescope.


 


 


 


 

En haut à gauche, le télescope de 28 cm de l'observatoire de Castanet-Tolosan (proche de Toulouse).
En haut à droite, le coupole du télescope de 60 cm de l'observatoire du Pic du Midi.
En bas à gauche, la bonnette fibre optique au foyer Newton du T60.
En bas à droite, le spectrographe eShel dans le laboratoire du T60 (cliquer sur l'image pour agrandir).

Top left: The 11inch telescope from Castanet-Tolosan (close to Toulouse city).
Top right: The dome for the 24inch telecope at Pic du Midi.
Bottom left: Fibre injection & guiding unit in f/3.5 prime focus of the 24inch telescope.
Bottom right: The eShel spectrograph inside the 24inch telescope control room (click to enlarge). 



De gauche à droite, le télescope de 60 cm, Olivier Thizy, Michel Pujol, Valérie Desnoux, Christian Buil. Photo V. Desnoux.

Left ot right: the 24inch telescope, Olivier Thizy, Michel Pujol, Valérie Desnoux, Christian Buil. Photo V. Desnoux.

L'observation de HD189733 b
Observation of HD189733 b


Le Champ de l'étoile HD189733 photographié avec un téléobjectif de 200 mm de focale et un Canon EOS 5D Mark II depuis le Pic du Midi. Le nord est en haut.
The field for HD189733 taken with a 200mm focal lens and a Canon EOS 5D Mark II from Pic du Midi site. North is Up.

HD189733 est une étoile de magnitude V=7.7, située aux coordonnées (2000), AD=20h00.7m, DEC=22°43'. La planète extrasolaire de HD189733 a été découverte par F. Bouchy & all en 2005.

HD189733 is a V=7.7 star (J2000: AD=20h00.7m, DEC=22°43'). The extrasolar planet around HD189733 has been discovered by F. Bouchy & all in 2005.

Journal des observations et mesures de HD189733 :
Observing log & results for HD189733: 

JD - 2450000

Radial velocity (km/s)

Observatory

5005.444
5006.412
5006.629
5007.402
5009.575
5010.473
5012.428
5013.430
5014.416
5018.431
5023.399
5024.398
5025.476
5028.483
5036.538
5037.477
5038.457
5039.382
5041.505
5042.457
5044.465
5047.438
5048.402

-2.539
-2.169
-2.005
-2.512
-2.380
-2.318
-2.277
-2.114
-2.464
-2.554
-2.554
-2.182
-2.586
-2.252
-2.554
-2.320
-2.629
-2.472
-2.436
-2.262
-2.180
-2.584
-2.187

Pic du Midi
Pic du Midi
Pic du Midi

Pic du Midi

Pic du Midi

Pic du Midi
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet

Pic du Midi
Pic du Midi
Pic du Midi
Pic du Midi
Castanet
Castanet
Castanet
Castanet
Castanet


En rouge l'éphéméride de HD189733 calculée avec les paramètres actuellement admis pour ce système :
P = 2.218573 jours (pédiode orbitale), K = 0,205 km/s (semi-ampitude vitesse radiale),
RV = -2.36 km/s (vitesse radiale moyenne), J0 = 2455019.88 (jour Julien amplitude max).
Les mesures de vitesse radiale de cette présente étude sont figurées avec des points noirs.

In red the ephemeris for HD189733-b with catalog parameters: P = 2.218573 days (orbital period), K = 0,205 km/s (Radial Velocity
half-amplitude), RV = -2.36 km/s (average Radial Velocity), J0 = 2455019.88 (max amplitude julian day).

 


A gauche, le périodogramme (algorithme de Lomb-Scargle) calculé à partir des points de mesure de HD189733. Deux périodes concurrentes se détachent, pratiquement avec la même intensité. L'une de 1,810 jours, l'autre de 2,223 jours. La période "officielle" de la planète est bien connue (2.219 jours) et on peut donc considérer que le pic franc à 2.223 jours est bien la signature de la planète, l'écart avec la valeur admise aujourd'hui n'étant que de 0,2%. La période parasite de 1.810 jours est très proche de l'alias à 1 jour de la période à 2.219 jours. Si P_orbite est la période orbitale de la planète et si P_obs est la période moyenne d'observation, alors la période de l'alias est donnée par :
P_alias = 1/(1/P_obs - 1/P_orbite)). Ici P_obs = 1 jour (l'étoile est quasiment toujours observée à la même heure de la nuit), donc P_alias = 1.820 jours.  Le pic parasite s'avère être la conséquence de l'observation régulière à la même heure de l'étoile, un problème classique en traitement de signal. Seul un décalage de l'heure sidérale d'observation ou mieux, des observations à d'autres longitudes, permettraient de réduire le pic parasite.

A droite, la courbe rouge montre l'éphéméride calculée à partir des paramètres orbitaux couramment admis, alors que la courbe en bleu est l'éphéméride calculée à partir du périodogramme des points de mesure. La valeur officielle de
K est de 205 m/s, alors que nous trouvons K = 199 m/s (mais avec une marge d'erreur de quelques dizaines de mètre probable). L'accord en phase est excellent. On trouve J0=2455006.53. Nous pensons que cette convergence entre la mesure et la prédiction permet d'annoncer une détection positive de la planète extrasolaire HD189733 b.

Left: periodogram (Lomb-Scargle algorythm) calculated from our data points for HD189733. Two periods pop up with almost the same intensity: 1.810 days and 2.223 days. Additional observations would be required to decide between the two but P=2.223 days is very close (0.2%) to the catalog period of 2.219 days and is the sign of the exoplanet. The 1.810 days peak period is a problable 1 days alias of the 2.219 days exact period of HD189733 exoplanet. The alias period is given by 1/P_alias = 1/P_obs - 1/P_orb. We find P_alias=1.820 days with P_obs=1 days (because the star is always observed nearly at the same time in the night) and P_orb=2.219 days.

Right: red curve is the ephemeris curve from catalog orbital parameters while the blue curve is the one calculated from our periodigram with K=199m/s (Vs K=205m/s in catalog). Phases are matching very well with J0=2455006.53. This match between catalog and our own observations allow us to announce the positive detection of extrasolar planet HD189733-b.
 


La courbe de phase finale de HD189733 regroupant l'ensemble des observations.
Les paramètres de l'orbite adoptés sont ceux admis pour ce système dans la litterature (courbe rouge).
La barre d'erreur des mesures est de 75 m/s. La variation sinusoïdale de la vitesse radiale du système est évidente et non ambigüe. Ce n'était pas un résultat acquis à l'avance compte tenu du faible éclat de l'étoile (magnitude 7.7, on le rappelle).

Final phase curve for HD189733 with all our observations. Red curve is the RV graph using catalog orbital parameters. Error bar is 75m/s. The sinusoidal RV curve is obvious and this is a great success considering the V=7.7 magnitude of the star !

En parallèle de l'observation par la méthode des vitesses radiale, la planète HD189733 b a été détectée le 19 juillet 2009 par la méthode photométrique à l'occasion de son passage devant son étoile. L'observation a été réalisée depuis une terrasse de l'observatoire du Pic du Midi avec un équipement très basique en regard de la discrète baisse d'éclat de lumière de l'étoile au moment du transit. Le télescope se résume à un téléobjectif Canon de 400 mm de distance focale ouvert à f/5.6 (le diamètre de la pupille est donc d'environ 70 mm seulement). La caméra est un simple appareil photographique numérique grand-public Canon 40D.

HD189733-b is a very interesting exoplanet as it is also transiting in front of the star and this is visible from Earth. In parallel to our observation with the Radial Velocity method, exoplanet HD189733-b has been detected on July 19, 2009 using photometry transit method. Observation has been conducted on the terrasse of Pic du Midi observatory with a very basic equipment: a 400mm focal Canon lens f/5.6 (diamter is around 70mm only) and a simple digital SLR Canon 40D.

 
L'équipement utilisé. En parallèle au téléobjectif, on trouve (en noir) une lunette de 80 mm de diamètre associé
à une caméra vidéo Watek 120N pour l'autoguidage.

Equipment used: 400mm Canon lens and 40D digital SLR, mounted on a 80mm refractor with Watec 120N for autoguiding.


Vue nocturne prise en fin de mesure du transit. La chaîne des Pyrénées se détache sur le fond. Contrairement à ce que laisse entrevoir
cette image, la nuit était très noire. Le cliché est en fait une pose 2 minutes réalisée avec un objectif très lumineux (Canon 85 mm f/1.2).

Night scene taken at the end of the transit. Pyrénees mountains are visible in the background.
Two-minutes exposure taken during the night (despite the feeling) with a very luminous 85mm f/1.2 lens.

Compte tenu du petit diamètre d'instrument, la scintillation atmosphérique est une source d'erreur potentiellement importante. D'après D. Dravins (Drawins & all, PASP, 110, 610, 1998), l'erreur RMS en magnitude liée à la scintillation est donné par :

Considering the small lens used for the photometry, atmospheric scintillation can be a large source of error. From D. Dravins (Drawins & all, PASP, 110, 610, 1998), the RMS error in magnitude due to scintillation is given by:

avec D, le diamètre du télescope en cm, Z, la distance zénithale de l'étoile, h, l'altitude de l'observatoire en mètres et T, le temps d'intégration en secondes.

with D the telescope diameter in cm, Z the zenithal distance of the star, h the altitude of the observing site and T the exposure time in seconds.

On vise une fluctuation maximale de l'ordre de 0,001 magnitude. Avec les paramètres de l'instrument et de l'observation (D=7,1 cm, Z = 20°, h = 3000 m), on trouve que le temps d'intégration caractéristique doit être d'environ 2 minutes pour lisser l'effet de la scintillation à ce niveau de précision. C'est le temps de pose adopté.

Pour éviter de saturer le détecteur au bout de ce laps de temps avec le signal de l'étoile HD189733 et avec celui des étoiles de comparaison, tout en accumulant la quantité de signal nécessaire à un objectif de précision proche du millième de magnitude, il est obligatoire de défocaliser notoirement l'image. Le résultat est peu esthétique mais efficace, ce qui est l'essentiel.

L'opération de défocalisation permet aussi de réduire les erreurs photométriques associées à l'organisation particulière des pixels dans le détecteur d'un appareil photo numérique (structure de Bayer) en liaison avec les défauts de guidage durant la pose. Cette initiative de défocalisation est fondamentale. Elle participe grandement au succès de l'opération. Soulignons que même avec le détecteur d'une caméra CCD astronomique noir et blanc il est très vivement recommandé de défocaliser systématiquement et assez fortement les images stellaires dès lors que l'on vise des mesures de haute précision. Ceci gomme les effets négatifs des non-uniformité de réponse intra-pixel (interne aux pixels) et d'échantillonnage de l'image (pixelisation), bien plus critiques qu'on ne pense habituellement. La figure suivante montre l'allure de l'étoile HD189733 avec le niveau de défocalisation adopté pour cette mesure (la tache fait environ 20 pixels de diamètre dans l'image RAW) :

For a maximal error of 0.001 magnitude and our equipement (D=7,1 cm, Z = 20°, h = 3000 m), the minimum exposure time should be 2 minutes. This is the adopted exposure time.

To avoid saturation in 2 minutes for HD189733 and comparison stars, it is mandatory to largely defocus the lens. It is not very nice aspect but it is efficient which is the most important point!

Defocusing also reduces the photometric error from the peculiar pixel arrangement of a digital SLR (Bayer matrix) and the guiding defect. This defocus is critical for the success of such photometry measures. Despite the general thinking, even with CCD camera it is highly recommended to defocus for high accuracy photometry; this reduce the impact of pixel sensibility differences and digitalization of the image. The following images show the HD189733 star with the defocusing (spot is around 20 pixels diameter in the RAW file) :

 

A gauche, l'image de l'étoile HD189733 telle qu'elle apparait dans le fichier RAW du Canon 40D (agrandissement
d'un facteur 4 par rapport à l'orignal). La structure en damier de la succession des pixels rouge, vert et bleu est bien visible.
A droite, la même image stellaire, mais après avoir isolé l'information des pixels "verts", les seuls utilisés pour cette présente étude.
On notera cependant au passage qu'un intérêt des appareils photographiques numériques pour la photométrie
est de fournir une information absolument simultanée dans trois couleurs de manière fort économique.

Left: HD189733 image from a RAW file of the Canon 40D (zoom x4). The grid structure comes from the successive red/green/blue pixels.
Right: the same using the green pixels only as we did for this study.
It is interesting to note here that digital SLR provide economical simultaneous data in three colors.


Carte du champ défocalisé (canal vert uniquement), où est identifié l'étoile HD189733 et les 4 étoiles de référence exploitées pour la réduction photométrique
(l'ensemble du traitement (flat-field, dark, ...) et de l'analyse a été réalisé avec le logiciel Iris). L'information photométrique est extraite par la technique de la
photométrie d'ouverture. Le fond de ciel est évalué en calculant le niveau médian dans l'anneau externe, ce qui réduit l'impact de la pollution
d'éventuelles étoiles parasites voisines. On note dans le flou la présence de la nébuleuse planétaire M27 (le nord est à gauche dans cette vue).

Star field (defocused, green pixels only) with HD189733 identified as well as the four reference stars used for photometry data reduction.
All processing (flat-field, dark...) and photometrical analysis have been done with IRIS software. The flux is calculated for each star by using circular
aperture photometric method. Median sky background level is evaluated in a circular area around each star to minimize the impact of close stars.
One can note the M27 planetary nebulae visible top-right. North is left. 


La courbe photométrique du transit de HD189733b du 19 juillet 2009. L'erreur RMS caractéristique estimée est de 0,0015 magnitude, ce qui très satisfaisant compte tenu du type de matériel utilisé. Les appareils photographiques numériques sont clairement aptes à réaliser des mesures de haute précision pour peu que l'on utilise des méthodes d'observation et réduction adaptées (à vrai dire, classiques en photomérie de précision - par exemple l'image donnée par le satellite photométrique Corot sur ces détecteurs CCD est défocalisée de manière fort semblable à la notre). En outre, la taille des capteurs des reflex numériques est très généreuse en regard du coût, ce qui permet d'envisager des programmes de recherche d'exoplanètes (ou autres objets variables) relativement économiques.

L'amplitude du transit de HD189733b est mesurée de 0,030 magnitude, alors que la valeur officiellement admise est de 0,028 magnitude. La date du début du transit est mesurée au jour Julien 2455031.4982 +/-0.0007. La date du début du transit est mesurée au jour Julien 2455031.5637 +/-0.0007. On en déduit la date médiane du transit, soit JD0 = 2455031.5309 +/- 0.0010, soit le 19 juillet 2009 à 0h44m38s UT. Le site ETD (http://var2.astro.cz/ETD) prévoyait 0h47m UT comme heure de milieu de transit. Bien que proche l'erreur de mesure, le présent transit semble indiquer une avance de l'ordre de 2 minutes sur la chronologie attendue. Cliquer ici pour télécharger le fichier photométrique.

Photometric curve of HD189733 b transit on July 19th, 2009. RMS error is estimated at 0.0015 magnitude which is very good in regards to the equipment used. With caution during acquisition and processing, digital SLR are excellent too for such projets. It is interesting to note that COROT satellite high accuracy photometry is also done with defocused images! Also, digital SLR have very large chip then field of view which should allow economical exoplanets and other variable stars hunting.

Amplitude for the transit is around 0.03 magnitude while catalogs give 0.028 magnitude. Transit start is measured as julian day 2455031.4982 +/-0.0007. End of transit is measured at JD=2455031.5637 +/-0.0007. Median transit time is JD0 = 2455031.5309 +/- 0.0010: July 19th, 2009 at 0h44m38s UT. ETD site (http://var2.astro.cz/ETD) scheduled at 0h47m UT middle of transit. The delta is 2 minutes. Click here to download photometric data. 

 

 



Image grand-champ de la région de HD189733 réalisé avec un Canon 5D Mark II est un objectif Canon 200 mm f/2.8 depuis l'observatoire du Pic du Midi le lendemain du transit (pose de 5 x 300 sec). Cliquer sur l'image pour agrandir.

Wide field around HD189733; Canon 5D Mark II and 200mm f/2.8 Canon lens from Pic du Midi observatory the following night after transit (5 x 300 sec exposure). Click on image to enlarge. 

L'observation de HD195019 b
Observation of HD195019 b


Le champ de l'étoile HD195019 photographié avec un téléobjectif de 200 mm de focale et un Canon EOS 5D Mark II. Le nord est en haut.
Field for HD195019 taken with a 200mm focal lens and a Canon EOS 5D Mark II from Pic du Midi site. North is Up. 

HD195019 est une étoile de magnitude V=6.9, située aux coordonnées (2000), AD=20h28.6m, DEC=18°46'. La planète extrasolaire de HD195019 a été découverte par D. Fischer & all en 1999.

HD195019 is a V=6.9 star (J2000: AD=20h28.6m, DEC=18°46'). Extrasolar planet HD195019-b has been discovered byD. Fischer & all in 1999 and has a long orbital period (above 18 days) compared to other exoplanets in this study.

Journal des observations et mesures de HD195019 :
Observing log & results for HD195019:

JD - 2450000

Radial velocity (km/s)

Observatory

5005.467
5006.439
5010.568
5013.491
5014.512
5019.473
5023.474
5024.530
5025.554
5028.5456
5033.533
5036.459
5037.559
5038.540
5039.461
5041.564
5042.511
5044.536
5047.512
5048.529

-91.249
-91.043
-90.962
-91.532
-91.344
-91.584
-91.194
-91.120
-91.034
-91.176
-91.371
-91.487
-91.469
-91.474
-91.364
-91.154
-91.094
-90.883
-91.080
-91.223

Pic du Midi
Pic du Midi
Pic du Midi
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet

Pic du Midi
Pic du Midi
Pic du Midi
Pic du Midi
Castanet
Castanet
Castanet
Castanet
Castanet


En rouge l'éphéméride de HD195019 calculée avec les paramètres actuellement admis pour ce système :
P = 18.202391 jours (pédiode orbitale), K = 0,275 km/s (semi-ampitude vitesse radiale),
RV = -91.3 km/s (vitesse radiale moyenne), J0 = 2455008.41 (jour Julien amplitude max).
Les mesures de vitesse radiale de cette présente étude sont figurées avec des points noirs.

In red the ephemeris using catalog orbital parameters for HD195019 system:
P = 18.202391 days (orbital period), K = 0,275 km/s (RV half amplitude),
RV = -91.3 km/s (average RV), J0 = 2455008.41 (max amplitude julian day).
Observed data points are in black. 


A gauche, le périodogramme (algorithme de Lomb-Scargle) calculé à partir des points de mesure de HD195019. La puissance maximale est observée pour la période de 18.106 jours, alors que la période orbitale issue des données professionnelles est de 18.202 jours. L'écart n'est que de 0,5%, ce qui est considéré comme raisonnablement faible, d'autant plus que seulement trois rotations consécutives de la planètes ont été couvertes par les observations. A droite, la courbe rouge montre l'éphéméride calculée à partir des paramètres orbitaux officiels alors que la courbe en bleu est l'éphéméride calculé à partir du périodogramme des points de mesure. La valeur de K officielle est de 275 m/s, alors que nous trouvons K = 261 m/s (mais avec une marge d'erreur de quelques dizaines de mètre probable). L'accord en phase trouvée est par ailleurs excellent. Le premier maximum de vitesse radiale de la période d'observation est trouvé à la date J0=2455008.78 alors que les éphémérides donnent J0=2455008.41, soit un décalage de 9 heures sur une période orbitale de 18 jours. Le décalage vertical entre les deux courbes est dû à une évaluation médiocre de la vitesse radiale de l'étoile, qui n'a pas d'incidence directe sur les paramètres orbitaux de l'exoplanète. Nous pensons que cette convergence entre la mesure et la prédiction permet d'annoncer une détection positive de la planète extrasoloaire HD195019 b.

Left, the periodogram (Lomb-Scargle algoritm) calculated with our observations of HD195019. Maximal power is at P=18.106 days; i.e. 0.5% from the catalog P=18.202 days, a good result specially as only 3 phases have been covered for this star. Right, red curve is the ephemeris while blue curve is our calculated RV curve (K=261m/s versus catalog K=275m/s). Phase matching is excellent (first maximum at J0=2455008.78 while catalog give J0=2455008.41; so a difference of 9 hours over a 18 days period). Vertical shift is due to poor average RV calculation which doesn't have any direct impact on orbital parameters calculation. We believe this good match allows us to announce a positive detection of exoplanet HD195019-b.



La courbe de phase finale de HD195019 regroupant l'ensemble des observations.
Les paramètres de l'orbite adoptés sont ceux admis pour ce système dans la litterature (courbe rouge).
La barre d'erreur des mesures est évaluée à 50 m/s. La variation sinusoïdale de la vitesse radiale du système est évidente et non ambigüe.

Final phase curve for HD195019 with all our observations. Red curve is calculated with the litterature orbital parameters. Error bar is 50m/s and the sinusoidale Radial Velocity is obvious.

L'observation de 51 Peg b
Observation of 51 Peg b


Le champ de l'étoile 51 Pegasus photographié avec un téléobjectif de 200 mm de focale et un Canon EOS 5D Mark II. Le nord est en haut.
Field for 51 Peg taken with a 200mm focal lens and a Canon EOS 5D Mark II from Pic du Midi site; North is Up.

51 Pegasis est une étoile de magnitude V=5.5, situé aux coordonnées (2000), AD=22h57.5m, DEC=20°46'. La planète extrasolaire de 51 Peg a été découverte par M. Mayor & D. Queloz en 1995 (c'est le premier cas de détection d'une exoplanète par la méthode des vitesses radiales).

51 Pegasis is a V=5.5 star (J2000: AD=22h57.5m, DEC=20°46'). 51 Peg b is discovered by M. Mayor & D. Queloz in 1995 with the Elodie spectrograph at the 193cm telescope of Observatory of Haute Provence. This is the historical first exoplanet discovered with Radial Velocity method.

Journal des observations et mesures de 51 Peg :
Observing log & results for 51 Peg: 

JD - 2450000

Radial velocity (km/s)

Observatory

5005.573
5006.595
5009.626
5013.585
5014.563
5016.568
5018.542
5023.590
5024.605
5025.611
5033.590
5036.605
5037.632
5038.619
5039.547
5041.622
5044.625
5047.591
5048.619

-33.125
-33.054
-33.104
-33.112
-32.994
-33.113
-33.099
-33.118
-33.195
-33.151
-33.197
-33.028
-33.168
-33.105
-33.045
-33.193
-33.083
-33.102
-32.981

Pic du Midi
Pic du Midi
Pic du Midi
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet

Pic du Midi
Pic du Midi
Pic du Midi
Pic du Midi
Castanet
Castanet
Castanet
Castanet



En rouge, l'éphéméride du système de 51 Peg calculée avec les paramètres actuellement admis pour ce système :
P = 4.230785 jours (pédiode orbitale), K = 0,056 km/s (semi-ampitude vitesse radiale),
RV = -33.11 km/s (vitesse radiale moyenne), J0 = 2455018.80 (jour Julien de l'amplitude max).
Les mesures de vitesse radiale de cette présente étude sont figurées avec des points noirs.
Dans le lot sélectionné pour cette étude, 51 Pegasus est de loin le système plus difficile.
L'amplitude en vitesse radiale est près de 10 fois plus faible que celle du système Tau Bootis, le
premier a être détecté par des astronomes amateurs. La demi-amplitude à mesurer (56 m/s) représente
moins de 0,2% de la résolution spectrale du spectrographe eShel (qui est de l'ordre de 28 km/s) !
Le moindre grain de sable dans la mesure détruit clairement tout espoir de détection.

Red: ephemeris using catalog parameters for 51 Peg-b: P = 4.230785 days (orbital period), K = 0,056 km/s (half-amplitude),
RV = -33.11 km/s (average RV), J0 = 2455018.80 (max amplitude julian day). Our observations are the black dots.
In this study, 51 Peg-b is clearly the most difficult object to detect with a RV amplitude 10 times lower than tau Boo-b,
the first exoplanet detected by amateurs using RV method! The half amplitude K=56m/s represents less than 0.2% of the
spectral resolution of the eShel spectrograph (around 28 km/s) ! The smallest sand dust in the process
would destroy all detection hope...

A gauche, le périodogramme (algorithme de Lomb-Scargle) calculé à partir des points de mesure de 51 Peg. Un pic de puissance est bien marqué pour la période orbitale de 4.186 jours, alors que la période reconnue de la planète extrasolaire de 51 Peg est de 4.231 jours. L'écart entre l'observation et la prévision est donc de 1%. On imagine mal une coïncidence, d'autant plus qu'à droite, est affichée en rouge l'éphéméride calculée à partir des paramètres orbitaux courants et en bleu l'éphéméride calculée à partir du périodogramme des points de mesure. Ce graphique montre un assez bon accord entre la semi-amplitude officielle de 56 m/s et le résultat d'observation, qui donne K = 65 m/s (la dispersion des points de mesure a tendance à accroitre l'amplitude). L'accord en phase est lui aussi satisfaisant. Nous pensons que cette convergence entre la mesure et la prédiction permet d'annoncer une détection positive de la planète extrasolaire 51 Peg b. Comme cela était prévu, cette observation est bien à la limite des possibilités de l'instrumentation utilisée, mais celle-ci a un gout particulier : la mythique première planète extrasolaire découverte par Michel Mayor et Didier Queloz par la méthode des vitesses  radiales est très probablement ici détectée pour la première fois par des astronomes non-professionnels !

Left, the periodogram (Lomb-Scargle algorithm) calculated with all our data point of 51 Peg. A peak in power is clearly visible with 4.186 days while known period is P=4.231 days... a proximity of 1%.  This is far from a random result specially with the graph in the right: red curve is the ephemeris and the blue curve our calculated curve from our data points. This graph shows a good match in half-amplitude with K=65m/s (Vs catalog K=56m/s). The phase is also matching well. We believe this represent a positive detection of 51 Peg-b exoplanet. Like we forecasted, this is the limit of this commercial instrument but this is a peculiar taste: this mythic exoplanet discovered in 1995 by Michel Mayor and Didier Queloz using Radial Velocity method is most probably detected for the first time here by non professional astronomers!


La courbe de phase finale de 51 Pegasus regroupant l'ensemble des observations.
Les paramètres de l'orbite adoptés sont ceux admis pour ce système dans la littérature (courbe rouge).
La barre d'erreur des mesures est évaluée à 50 m/s. L'erreur est donc du même ordre de grandeur que l'amplitude de vitesse radiale recherchée, et seule l'accumulation de mesures donne l'espoir de dégager une tendance faisant sortir la planète extrasolaire du bruit. Les conditions météorologiques ne nous ont pas été très favorables durant la période d'observation, cependant, le nombre de points engrangé commence à dessiner une forme sinusoïdale qui ne semble pas ambiguë.

Final phase graph for 51 Peg with all our data points.
Red curve is calculated from catalog orbital parameters. Error bar is around 50m/s, about the same as the half-amplitude; cumulating lot of data/measurements is the only way to get a trend and highlight the exoplanet impact on the star movement from the instrumental noise. The radial velocity curve shows a sinusoidal shape matching the ephemeris.

L'observation de Tau Bootis b
Observation of Tau Bootis b


Le champ de l'étoile Tau Boo photographié avec un téléobjectif de 200 mm de focale et un Canon EOS 5D Mark II. Le nord est en haut.
Field for tau Boo taken with a 200mm focal lens and a Canon EOS 5D Mark II from Pic du Midi site. North is Up. 

Tau Bootis est une étoile de magnitude V=4.5, situé aux coordonnées (2000), AD=13h47.2m, DEC=17°27'. La planète extrasolaire de Tau Boo a été découverte par P. Butler & all en 1997.

Tau Bootis is a V=4.5 star (J2000: AD=13h47.2m, DEC=17°27'). The extrasolar planet tau Boo-b has been discovered by P. Butler & all in 1997. 

Journal des observations et mesures de 51 Peg :
Observing log & results for 51 Peg: 

JD - 2450000

Radial velocity (km/s)

Observatory

4908.592
4909.528
4910.545
4911.528
4912.514
4913.476
4914.501
4917.558
4918.367
4920.466
4954.385
4967.410
4968.386
4970.400
4971.459
4979.399
4980.403
4981.389
4982.406
4983.461
4984.408
4985.392
4986.378
5005.407
5006.353
5007.362
5009.528
5010.399

-16.385
-16.124
-17.003
-16.778
-16.156
-16.636
-17.011
-17.135
-16.471
-17.053
-17.066
-17.251
-16.349
-17.012
-16.482
-16.422
-17.103
-16.626
-16.073
-16.929
-16.736
-16.113
-16.611
-16.142
-16.914
-17.056
-16.570
-17.183

Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet
Castanet

Pic du Midi
Pic du Midi
Pic du Midi
Pic du Mid
Pic du Midi


En rouge l'éphéméride de Tau Bootis calculée avec les paramètres actuellement admis pour ce système :
P = 3.312463 jours (pédiode orbitale), K = 0,461 km/s (semi-ampitude vitesse radiale),
RV = -16.63 km/s (vitesse radiale moyenne), J0 = 2454935.83 (jour Julien amplitude max).
Les mesures de vitesse radiale de cette présente étude sont figurées avec des points noirs.

Red curve: ephemeris using catalog orbital parameters : P = 3.312463 days (orbital period), K = 0,461 km/s (half-amplitude radial velocity),
RV = -16.63 km/s (average radial velocity), J0 = 2454935.83 (max velocity julian day).
Our data points are the black dots.


A gauche, le périodogramme (algorithme de Lomb-Scargle) calculé à partir des points de mesure de tau Bootis. La puissance maximale est observée pour la période de 3.317 jours, alors que la période orbitale issues des données professionnelles est de 3.312 jours. Ces deux valeurs s'accordent pratiquement à 0,1% près. A droite, la courbe rouge montre l'éphéméride calculée à partir des paramètres orbitaux courants alors que la courbe en bleu est l'éphéméride calculée à partir du périodogramme des points de mesure. Les deux courbes se superposent très bien. La valeur de K officielle est de 461 m/s, alors que nous trouvons ici K = 493 m/. L'accord en phase trouvée est par ailleurs excellent. Le premier maximum de vitesse radiale de la période d'observation est trouvé à la date J0=2454909.189 alors que les éphémérides donnent J0=2454909.330, soit un décalage temporel de 3 heures environ. Nous pensons que cette convergence entre la mesure et la prédiction permet d'annoncer une détection positive de la planète extrasolaire tau Boo b.

Left: periodogram (Lomb-Scargle algorithm) calculated from our data points of tau Boo. Maximal power is reached for P=3.317 days compared to predicted P=3.312 days... an accuracy of nearly 0.1% ! Right: red curve is drawn from catalog orbital parameters while blue curve is calculated from our data points. Official K=461m/s while we find K=493m/s. The phases are matching perfectly. First maximum is at J0=2454909.189 while ephemeris give J0=2454909.330, so a difference if about 3 hours. We believe this good match allows us to announce the detection of exoplanet tau Boo-b.


La courbe de phase finale de tau Bootis regroupant l'ensemble des observations.
Les paramètres de l'orbite adoptés sont ceux prédits pour ce système dans la litterature (courbe rouge). La barre d'erreur des mesures est évaluée à 75 m/s. La variation sinusoïdale de la vitesse radiale du système est évidente et non ambigüe.

Final tau Bootis phase diagram with all our observations (dots) and the predicted orbital curve (red). Error bar is estimated at 75m/s and the sinusoidal curve is well visible.

Syntèse
Summary

 
Cliquer pour agrandir
Click for enlarge

Back