Comprendre

Auteur: Philippe Thébault

Datation Radio-Isotopique


Datation radio-isotopique: Principe

Certains noyaux atomiques, appelés radioisotopes, sont naturellement instables et peuvent spontanément se désintégrer en noyaux moins massifs et stables, libérant de l’énergie sous forme de rayonnement. La décroissance du nombre N d’un type de radioisotopes {}^{p}A en un élément {}^{q}B suit une décroissance exponentielle dN/dt prop -lambda*N. Où la quantité t_{1/2} = ln(2)/\lambda est le temps de demi-vie caractéristique de l’élément. Dans le cas le plus simple où aucun élément B n’est présent initialement le temps t peut être trouvé directement par le rapport B/A. Mais en réalité les choses ne sont jamais aussi simples, et ce pour au moins 2 raisons :

1) Il y a a priori toujours du {}^{q}B présent à t=0 ou, du moins, il est impossible d’exclure cette présence. Et on ne connait pas a priori cette quantité initiale de {}^{q}B

2) l’isotope {}^{q}B n’est souvent pas le seul possible pour l’élément B, qui peut également exister sous la forme {}^{q'}B, ce qui peut fortement compliquer les choses. En effet, dans un matériau à l'état gazeux ou liquide, les isotopes {}^{q}B et {}^{q'}B vont naturellement s'équilibrer entre eux à une valeur d'équilibre. En conséquence, dès qu'un matériau fond ou fusionne, toute information sur la quantité d'isotope {}^{q}B produite par désintégration de {}^{p}A va être perdue par cette mise à l'équilibre isotopique (autrement dit, toute fusion est un "reset" des isotopes de B).

Ceci va rendre la datation plus compliquée, mais elle reste néanmoins possible, du moins pour remonter jusqu'au moment de la dernière condensation du matériau. On peut grosso-modo distinguer 2 types de datation : datation absolue et datation relative.


Datation absolue

La désintégration de {}^{238}U en {}^{206}Pb a un temps de ½ vie de 4.47 109 ans, idéal pour mesurer l’âge des plus anciens corps du système solaire. Mais {}^{206}Pb n’est pas l’isotope naturel du Pb, qui est {}^{204}Pb. On obtient alors la relation suivante, liant les abondances de {}^{238}U, {}^{206}Pb et {}^{204}Pb:

\left (\frac{{}^{206}Pb}{{}^{204}Pb} \right )_P = \left (\frac{{}^{206}Pb}{{}^{204}Pb} \right )_I + \left (\frac{{}^{238}U}{{}^{204}Pb} \right )_I  \left ( 1- e^{-\lambda_{238}t}  \right )

= \left (\frac{{}^{206}Pb}{{}^{204}Pb} \right )_I + \left (\frac{{}^{238}U}{{}^{204}Pb} \right )_{P}  \left ( e^{\lambda_{238}t} -1  \right )

où les indices P et I se réfèrent aux abondances actuelles et initiales, respectivement. Le moment "initial" correspond à l'instant où l'objet en question s'est solidifié pour la dernière fois. En effet, dès que le corps fond ou se sublime en gaz, les proportions des 2 isotopes {}^{206}Pb et {}^{204}Pb s'équilibrent rapidement à leur proportion "naturelle" et toute information sur la désintégration de {}^{238}U est perdue (voir page précédente ). A cet instant initial le rapport {}^{206}Pb/{}^{204}Pb est donc égal à la valeur d'équilibre. En revanche, une fois le corps solidifié, un excès de l'isotope {}^{206}P va petit à petit se créer à mesure que {}^{238}U se désintègre. La variable inconnue est ici la quantité initiale absolue de {}^{206}Pb (ou de {}^{204}Pb), que l'on ne connaît pas a priori. Heureusement, il existe un deuxième type de désintégration d’U en Pb, la réaction {}^{235}U \Rightarrow {}^{207}Pb , dont le temps de vie est de 704 106 ans, et qui va nous permettre de contraindre les abondances initiales. Les équations sont alors:

\left (\frac{{}^{207}Pb}{{}^{204}Pb} \right )_P = \left (\frac{{}^{207}Pb}{{}^{204}Pb} \right )_I + \left (\frac{{}^{235}U}{{}^{204}Pb} \right )_P  \left (e^{\lambda_{235}t} -1 \right )

\left (\frac{{}^{206}Pb}{{}^{204}Pb} \right )_P = \left (\frac{{}^{206}Pb}{{}^{204}Pb} \right )_I + \left (\frac{{}^{238}U}{{}^{204}Pb} \right )_P  \left (e^{\lambda_{238}t} -1 \right )

Et donc: F = \left [\frac{\left (\frac{{}^{207}Pb}{{}^{204}Pb} \right )_P -\left (\frac{{}^{207}Pb}{{}^{204}Pb} \right )_I }{ \left (\frac{{}^{206}Pb}{{}^{204}Pb} \right )_P - \left (\frac{{}^{206}Pb}{{}^{204}Pb} \right )_I} \right ]=\left ( \frac{1}{137.88} \right )\left ( \frac{e^{\lambda {}_{235}t} -1}{e^{\lambda {}_{238}t} -1} \right )

où 137.88 est la valeur présente de \frac{{}^{238}U}{{}^{235}U}, qui est une constante globale du système solaire actuel, et \lambda_{235}= \ln(2) / \tau_{235}, \lambda_{238}= \ln(2) / \tau_{238}. Cette relation est directement exploitable pour toute météorite non-homogène initialement, mais dont tous les composants se sont formés à la même époque. En effet, dans ce cas, les rapports initiaux {}^{206}Pb/{}^{204}Pb et ^{207}Pb/{}^{204}Pb sont les mêmes partout dans la météorite et sont égaux à leurs valeurs d'équilibre (indiquées par a0 et b0 sur la figure). Par conséquent, dans un graphe {}^{206}Pb/{}^{204}Pb vs. ^{207}Pb/{}^{204}Pb , toutes les mesures du rapport F doivent se situer une même droite, appelée isochrone, dont la pente va directement donner l’âge de la météorite (cf. Figure).

images/Approfondissement-1-Pb-Pb.png
Isochrone isotopique construite à partir des abondances actuelles relatives de 207Pb, 206Pb et 204Pb en différents endroits d’une même météorite
Crédit : Observatoire de Paris

Datation relative

La désintégration radioactive permet également de dater des corps même bien après la disparition des radionucléotides concernés (c'est à dire bien au delà du temps de 1/2-vie de la désintégration concernée). C’est le cas par exemple de la désintégration {}^{26}Al \Rightarrow{}^{26}Mg, dont le temps de ½ vie est de « seulement » 720 000ans. Notons que ni {}^{26}Al, ni {}^{24}Mg ne sont des isotopes naturels de leurs éléments, qui sont, respectivement, {}^{27}Al et {}^{24}Mg. A la différence de la désintégration d’{}^{238}U, il n’existe aujourd’hui plus de {}^{26}Al que l’on puisse mesurer. En principe, on a donc :

({}^{26}Mg)_P=({}^{26}Mg)_I+({}^{26}Al)_I

Cette équation n’est pas d’une grande aide en elle-même, mais, comme pour la datation absolue, on peut tirer parti de la non-homogénéité d’une météorite donnée. Si en effet deux endroits de cette météorite avaient initialement des teneurs totales en Al (tous isotopes confondus) différentes, mais que la proportion de \frac{{}^{26}Al}{{}^{27}Al} était, elle, la même, alors l’excès de {}^{26}Mg ne sera aujourd’hui pas partout le même, et cet excès sera relié à l’abondance actuelle locale de \frac{{}^{27}Al}{{}^{24}Mg} par la relation :

\left ( \frac{{}^{26}Mg}{{}^{24}Mg} \right)_P = \left ( \frac{{}^{26}Mg}{{}^{24}Mg} \right)_I + \left ( \frac{{}^{26}Al}{{}^{27}Al} \right)_I + \left ( \frac{{}^{27}Al}{{}^{24}Mg} \right)_P

Les mesures de \frac{{}^{26}Mg}{{}^{24}Mg} et \frac{{}^{27}Al}{{}^{24}Mg} en différents endroits de la météorite vont alors tracer une isochrone dont la pente donnera la teneur initiale en \frac{{}^{26}Al}{{}^{27}Al} (voir Figure).

Maintenant, si on compare les teneurs initiales de \frac{{}^{26}Al}{{}^{27}Al} obtenus pour différentes météorites, on peut obtenir une datation relative des temps de formation de ces météorites. En effet, étant donné le temps de vie très court de {}^{26}Al, sa teneur par rapport à {}^{27}Al pourra être très différente suivant l’instant où la météorite s’est formée. Si on compare les valeurs \left (\frac{{}^{26}Al}{{}^{27}Al} \right )_Idans 2 météorites différentes, on obtient ainsi la datation relative de leur formation par la formule:

\left (\frac{{}^{26}Al}{{}^{27}Al} \right)_I^1 = \left (\frac{{}^{26}Al}{{}^{27}Al} \right)_I^2 e^{-\lambda(t_1- t_2)}

t_1 et t_2 sont les instants de formation des 2 météorites, et \lambda est le taux de désintégration de la réaction {}^{26}Al \Rightarrow {}^{26}Mg

Il faut cependant ici bien faire attention à deux points très importants :

  1. ces mesures relatives ne sont possibles que si {}^{26}Al était uniformément distribué dans la nébuleuse initiale. Les dernières recherches semblent cependant montrer que tel était bien le cas. 
  2. la fraction de \left (\frac{{}^{26}Al}{{}^{27}Al} \right )_{0} estimée correspond à ce qu’elle était au moment de la dernière condensation de l’objet. En effet, quand le corps est dans un état fluide (gaz, liquide), les abondances de {}^{26}Mg et {}^{24}Mg se ré-équilibrent automatiquement à leurs proportions « naturelles » et tout excès antérieur de {}^{26}Mg est effacé. Autrement dit, le chronomètre isotopique se « reset » lors de tout épisode de très forte température.
images/Approfondissement-1-Al-Mg-Alliende.png
Isochrone obtenue à partir des abondances actuelles des différents isotopes de Mg et Al dans la météorite carbonée « Alliende »
Crédit : Observatoire de Paris

La barrière du mètre

Comme nous l’avons vu (cf. lien), une des étapes les plus délicates du scénario standard de formation planétaire est celle qui fait passer des premières poussières condensées dans la nébuleuse aux planétésimaux kilométriques. La principale difficulté étant atteinte pour des corps de 10cm-1m, pour lesquels les vitesses de collision deviennent trop élevées pour permettre l’accrétion, et qui vont de plus avoir un mouvement de dérive très rapide vers l’étoile centrale. Ces deux problèmes sont tous deux liés à l’action du gaz sur les corps solides. En effet, toute particule solide plongée dans un milieu gazeux subit la friction de ce gaz, qui va être proportionnelle à la surface de contact entre le gaz et l’objet. Cette friction peut ainsi s’exprimer sous la forme

F \propto K\sigma\Delta_V

\Sigma est la section efficace du corps ( \pi r^2) et \Delta_V la différence de vitesse entre le gaz et la particule. Si maintenant on applique le principe fondamental de la dynamique, on obtient que l’accélération a due au gaz vaut:

a=\frac{1}{M}K\sigma\Delta_V \propto \frac{1}{r}\Delta_V(r)

On pourrait a priori se dire que, dans le cas présent, \Delta_V est nul car aussi bien le gaz que les corps solides orbitent autour de l’étoile suivant les mêmes lois de Kepler, et donc en principe à la même vitesse orbitale v_{kep}=(GM/R)^{1/2}. Mais il y a en fait une différence, car le gaz est lui, en plus, soumis à une force de pression due au gradient de température et de densité dans la nébuleuse. Dans un disque proto-planétaire de type MMSN, cette force de pression s’exerce de l’intérieur vers l’extérieur et tend donc à contrebalancer la gravitation de l’étoile. Tout se passe donc comme si le gaz « percevait » une étoile de masse \beta M_* < M_* et aura donc une vitesse Képlérienne v_{kep'} inférieure au v_{kep} d’un corps solide orbitant dans la vide. On dit alors que le disque de gaz est « sub-Képlérien ».

Le comportement de corps solides plongés dans ce disque de gaz est alors compris entre 2 extrêmes : les particules les plus petites (<mm) sont piégées dans le gaz et bougent avec lui, et on a dans ce cas \Delta_V= 0 et donc a=0. A l’autre extrême, les planétésimaux très massifs sont, eux, découplés du gaz et subissent de la friction, mais ne vont pas beaucoup en pâtir car le rapport \frac{1}{r} est tout petit et ils vont donc être très difficile à bouger par le gaz. Donc la aussi a=0. Entre ces 2 extrêmes, il existe un régime intermédiaire avec des corps suffisamment gros pour être découplés du gaz ( \Delta_V > 0) mais pas suffisamment massifs pour être insensibles à la friction. Ce régime de taille intermédiaire se situe autour de 10cm-1m. Pour des corps de cette taille l’effet de la friction gazeuse est maximal (cf. image). Et comme v_{kep}(gaz) < V_{kep}, cette friction aura tendance à ralentir les corps solides et à les faire dériver vers l’étoile. Cette vitesse de dérive peut atteindre plus de 50m/s, ce qui correspond à 1UA en moins de 100ans !

images/Approfondissement-2-meter.png
Vitesse de dérive radiale (en cm/s) vers l’étoile centrale, en fonction de taille des corps solides présents dans la nébuleuse protoplanétaire de type MMSN. On voit très clairement qu’un pic est atteint pour des corps proches de la taille du mètre.
Crédit : Observatoire de Paris

Dynamique d'un disque de planétésimaux

Une population de planétésimaux orbitant le soleil ne pourra former des corps plus massifs que si les vitesses de rencontres mutuelles \Delta V_{coll} sont, pour une fraction importante de ces rencontres, inférieures à \beta V_{lib} ; où V_{lib} est la vitesse de libération à la surface de 2 planétésimaux en collision et \beta<1 un paramètre prenant en compte la dispersion d’énergie lors de l’impact. Toute la question est alors de savoir si, effectivement, ce critère \Delta V_{coll} < \beta V_{lib} va être rempli dans un disque de planétésimaux kilométriques laissé à lui même.

L’état dynamique d’un tel disque va dépendre de l’équilibre entre plusieurs mécanismes: la gravité mutuelle des planétésimaux, la force de friction du gaz primordial toujours présent à ce stade, la dissipation d’énergie lors des collisions physiques, et, bien sur, la gravité de l’étoile autour de laquelle tous les corps orbitent. Si on fait l’hypothèse simplificatrice que tous les planétésimaux ont la même taille r, alors on peut montrer que le disque va tendre vers un état stationnaire où les vitesses de rencontres vont en moyenne être de l’ordre de V_{lib}(r). Si en effet \Delta V_{coll} \ll V_{lib}, alors les déflections gravitationnelles lors de rencontres proches vont automatiquement augmenter \Delta V_{coll}. Et si, à l’inverse, \Delta V_{coll} \gg V_{eva}, alors la dissipation d’énergie cinétique lors des collisions va être très forte et fera diminuer \Delta V_{coll} . Cet équilibre autour de \Delta V_{coll} = V_{lib} est plutôt une bonne nouvelle, car il entraine qu’une fraction des collisions vont effectivement permettre l’accrétion des corps (sachant qu’il y aura toujours une dispersion des vitesses de collision autour leur valeur moyenne).

Sachant que tous ces corps sont en orbite autour d’une étoile, par exemple le Soleil, les vitesses relatives de collisions vont être directement liées à l’excentricité (et à leur inclinaison si on est en 3D) de leurs orbites : plus les orbites sont excentriques, plus \Delta V_{coll} augmente, plus elles sont circulaires, plus \Delta V_{coll} tend vers zero. Pour de petites excentricités, à l’ordre zero on peut écrire que

\Delta V_{coll} \sim eV_{kep}

Il faut réaliser que, pour des planétésimaux kilométriques, on a affaire à des vitesses très faibles, car les V_{lib} de tels corps sont de l’ordre de quelques mètres par seconde seulement. Ceci se traduit par des excentricités orbitales très faibles, de l’ordre de 0.0001 ! (cf. exercice).


Accrétion «Boule de neige»

Si tous les planétésimaux du disque avaient exactement la même taille r et grandissaient tous à la même vitesse, alors leur taux de croissance serait égal à \frac{dm}{dt} = \pi (r+r)^2 \sigma   \frac{V}{h} dt\sigma est la densité surfacique de planétésimaux, V leur vitesse relative de collision et h l’épaisseur du disque. Si on fait l’approximation raisonnable que V \sim e V_{kep} (cf. page précédente) et que h = i a (a, distance à l’étoile et e excentricité moyenne de l’orbite des planétésimaux), alors on obtient \frac{dm}{dt} = 4 \pi r^2 \sigma \frac{e}{i} \Omega_k avec \Omega_k = \frac{V_{kep}}{a}, vitesse angulaire Keplerienne. On trouve alors que \frac{dr}{dt} = cste, et que la croissance en taille est linéaire avec le temps. Pour une MMSN à 1UA on trouve qu’il faut alors quelques 106 ans pour former un corps de 1000km (cf. EXERCICE)

Mais il semble qu’en réalité l’accrétion suive un chemin beaucoup plus rapide et efficace, mais très sélectif, appelé accrétion « boule de neige ». Il est en effet plus que probable que, dans tout disque réel, toutes les tailles ne sont pas identiques et que, localement, certains planétésimaux soient, par hasard, légèrement plus grands (de taille r_1>r) que ceux qui les entourent. De ce fait, ils ont une vitesse de libération V_{lib}(r_1) supérieure à celle des corps environnants. En conséquence, ils vont légèrement infléchir la trajectoire des autres corps vers eux. On peut paramétriser cette déflection en considérant que le corps r_1 a une section efficace « effective » \Sigma plus grande que sa simple section efficace géométrique \pi r_1^2. On a alors

\Sigma = \pi (r_1+r)^2 \left [ 1+ \left (\frac{V_{lib}(r_1,r)}{\Delta V} \right )^ 2\right ]

\left (1+ \left (\frac{V_{lib}(r_1,r)}{\Delta V} \right )^ 2 \right ) est appelé le terme de « focalisation gravitationnelle ». Du fait de cette surface efficace« dilatée », le corps r_1 va croître plus vite que les autres. Le rapport \frac{r_1}{r} va donc augmenter, ce qui a pour effet d’encore augmenter la focalisation gravitationnelle, et donc le taux de croissance de r_1, et ainsi de suite. La croissance de ce corps initialement légèrement privilégié va donc rapidement s’emballer.

images/Approfondissement-4-Runaway.png
Accrétion boule de neige : un corps initialement plus gros que les planétésimaux qui l‘entourent va légèrement infléchir la trajectoire de ceux-ci vers lui (Figure A). De ce fait, il va croître plus rapidement que les corps qui l’entourent, ce qui va encore accentuer sa tendance à infléchir vers lui l’orbite de ceux-ci (Figure B). Le processus s’auto-amplifie de lui même et conduit à la formation rapide d’un embryon planétaire alors que la majorité des autres planétésimaux n’a pas accrété de matière (Figure C).
Crédit : Philippe Thebault

Photo-évaporation et dispersion du gaz

L’un des moyens les plus efficaces pour disperser le disque de gaz en moins de 10 millions d’années (voir la page de cours) est le couplage entre (l'accrétion visqueuse du disque) et la photo-évaporation du gaz. Ce dernier mécanisme est la conséquence de l’effet du rayonnement ultra-violet (UV) de l’étoile sur les molécules de gaz, essentiellement H et He. L’interaction des photons UV va en effet chauffer le gaz, c’est à dire lui donner une plus grande agitation thermique, et si cette vitesse d’agitation thermique dépasse la vitesse Képlerienne locale, alors le gaz est éjecté. Comme V_{kep} est \propto a^{-0.5} mais que l’énergie transportée par un photon UV ne diminue pas avec la distance à l’étoile, ce sont les molécules des régions extérieures qui seront le plus facilement éjectées lors d’interaction photon-gaz (mais le flux de photon, et donc le taux d’interaction avec le gaz, va, lui, diminuer avec a). On peut ainsi calculer qu’il existe un rayon critique, appelée r_g, au delà duquel l’énergie déposée par photo-évaporation dépasse l’énergie orbitale :

r_g = \frac{(\gamma -1)}{2\gamma} \frac{GM \mu}{K_B T} \approx 1.4 \frac{M/M_{\odot}}{T/10^4 K} AU

T est la température du gaz par suite de l’interaction avec un photon, \mu est le poids atomique moyen du gaz,K_B la constante de Boltzman, \gamma le ratio des chaleurs spécifiques (5/3 pour un gaz mono-atomique), G la constante gravitationnelle et M la masse de l’étoile. Dans les faits, à cause de la rotation du gaz, le rayon critique de dispersion est plutôt égal à r_{cr} \sim 0.2 r_g. Pour un rayonnement UV typique, on a T \sim 10^4 K. Pour une étoile de type solaire et un disque de type MMSN , on obtient alors r_{cr} = 2 \verb?UA?.

Cependant, tant que le disque est très dense, à la distance r_{cr} le flux de matière  \dot{M_{acc}} vers l’intérieur du disque dû à la viscosité est supérieur au flux de matière M_w éjecté par photo-évaporation (cf. IMAGE, CASE 1). Mais à mesure que la masse du disque diminue par accrétion sur l’étoile, \dot{M_{acc}} va décroître, jusqu’à ce qu’on atteigne le point où \dot{M_{acc}}(r_{cr}) = M_w(r). A partir de ce moment, la matière gazeuse au-delà de r_{cr} est éjectée du système avant de pouvoir franchir la frontière r_{cr} pour se mettre « à l’abri ». Une ouverture est alors créée dans le disque, qui se retrouve coupé en deux. La partie interne du disque est protégée de la photo-évaporation, mais est tellement petite (<2UA) qu’elle va très rapidement être accrétée sur l’étoile (alors qu’auparavant le flux de matière spirallant sur l’étoile était compensé par de la matière venant de plus loin dans le disque). La partie externe du disque va elle se disperser progressivement de l’intérieur (r_{cr}) vers l’extérieur. Les modèles numériques indiquent que la troncature du disque se fait au niveau de r_{cr} sur une échelle de quelques millions d’années, tandis que la dispersion qui s’en suit est beaucoup plus rapide, quelques 105 ans peut-être.

images/Approfondissement-5-PHOTOEVAP.png
Différentes étapes du processus de dispersion d’un disque proto-planétaire gazeux par photoevaporation et accrétion visqueuse du disque.
Crédit : Observatoire de Paris