Radioactivité : Décroissance et Datation : La modélisation mathématique du temps nucléaire

Problématique

Comment la résolution mathématique des équations différentielles de décroissance radioactive permet-elle d’estimer des durées géologiques et archéologiques ?

---

Introduction

Bonjour à toutes et à tous. Le concept de temps est l’une des notions les plus fascinantes des sciences physiques. Si l’humanité a conçu des horloges mécaniques ou solaires pour mesurer le passage des heures, la nature possède ses propres horloges internes, logées au cœur même des noyaux atomiques instables : c’est la radioactivité. Découverte à la fin du XIXe siècle par Henri Becquerel et Pierre et Marie Curie, la radioactivité est un phénomène physique spontané qui transforme la matière. Mais au-delà de ses aspects physiques, c’est sa modélisation mathématique rigoureuse qui a offert aux scientifiques un outil d’une puissance inouïe : la possibilité de dater des objets anciens avec une précision remarquable.

Passionné par la physique nucléaire et la rigueur de l’analyse mathématique, j’ai voulu étudier comment ces deux disciplines s’articulent pour quantifier le passé. Ma problématique est la suivante : Comment la résolution mathématique des équations différentielles de décroissance radioactive permet-elle d’estimer des durées géologiques et archéologiques ?

Pour y répondre, nous analyserons d’abord le passage d’une loi probabiliste individuelle à une modélisation déterministe sous forme d’équation différentielle. Puis, nous verrons comment cette loi se traduit concrètement pour la datation archéologique au Carbone 14.

---

I. Modélisation mathématique de la décroissance (environ 4 min)

À l’échelle microscopique, la désintégration d’un noyau radioactif individuel est un phénomène régi par le hasard. Elle possède trois caractéristiques fondamentales : elle est aléatoire (on ne peut prédire quand un noyau va se désintégrer), spontanée (elle se produit sans intervention extérieure) et indépendante des autres noyaux et des conditions extérieures (comme la température ou la pression).

Pour un noyau radioactif unique, la probabilité de se désintégrer pendant un court intervalle de temps $\Delta t$ est proportionnelle à cet intervalle. On note cette probabilité :

$$P(\text{désintégration}) = \lambda \cdot \Delta t$$

Où $\lambda$ est la constante radioactive du radio-isotope considéré, exprimée en $\text{s}^{-1}$.

Lorsque l’on étudie une population macroscopique contenant un nombre extrêmement grand $N(t)$ de noyaux instables (de l’ordre de $10^{20}$ noyaux), nous pouvons passer d’un modèle probabiliste à une loi continue déterministe grâce à la loi des grands nombres. La variation moyenne $\Delta N(t)$ du nombre de noyaux pendant la durée $\Delta t$ s’écrit comme le produit du nombre de noyaux par la probabilité individuelle de désintégration, affecté d’un signe moins car la population diminue :

$$\Delta N(t) = -N(t) \cdot \lambda \cdot \Delta t$$

En divisant par $\Delta t$ et en faisant tendre la durée $\Delta t$ vers zéro, on obtient la dérivée temporelle $\frac{dN(t)}{dt}$ :

$$\lim_{\Delta t \to 0} \frac{\Delta N(t)}{\Delta t} = \frac{dN(t)}{dt} = -\lambda \cdot N(t)$$

Nous obtenons ainsi une équation différentielle linéaire du premier ordre, homogène et à coefficients constants :

$$\frac{dN(t)}{dt} + \lambda \cdot N(t) = 0$$

D’après le cours de mathématiques de Terminale, l’unique solution de cette équation différentielle de la forme $y’ + ay = 0$ vérifiant la condition initiale $N(0) = N_0$ est la fonction exponentielle :

$$N(t) = N_0 \cdot e^{-\lambda t}$$

Où $N_0$ est la population de noyaux radioactifs à l’instant initial $t=0$.

Pour caractériser la rapidité de cette décroissance, on définit la demi-vie (ou période radioactive), notée $t_{1/2}$. C’est la durée nécessaire pour que la moitié des noyaux initialement présents se désintègrent. Mathématiquement, elle vérifie l’égalité $N(t_{1/2}) = \frac{N_0}{2}$, ce qui donne :

$$N_0 \cdot e^{-\lambda t_{1/2}} = \frac{N_0}{2} \implies e^{-\lambda t_{1/2}} = \frac{1}{2}$$

En appliquant la fonction logarithme népérien ($\ln$) de chaque côté :

$$-\lambda \cdot t_{1/2} = \ln\left(\frac{1}{2}\right) = -\ln(2) \implies t_{1/2} = \frac{\ln(2)}{\lambda}$$

Cette relation fondamentale montre que la demi-vie $t_{1/2}$ est inversement proportionnelle à la constante radioactive $\lambda$. Chaque isotope possède sa propre demi-vie, variant de quelques microsecondes à des milliards d’années.

---

II. Applications pratiques : La datation au Carbone 14 (environ 3 min 30s)

En pratique, mesurer directement le nombre de noyaux $N(t)$ est difficile. On mesure plutôt l’activité $A(t)$ de l’échantillon, qui est le nombre de désintégrations par seconde, mesurée en Becquerels ($\text{Bq}$). L’activité est définie par :

$$A(t) = -\frac{dN(t)}{dt} = \lambda \cdot N(t) = \lambda \cdot N_0 \cdot e^{-\lambda t} = A_0 \cdot e^{-\lambda t}$$

Où $A_0 = \lambda \cdot N_0$ est l’activité initiale à $t=0$. L’activité suit la même loi de décroissance exponentielle.

L’application la plus célèbre de cette loi est la datation au Carbone 14 ($^{14}C$). Le carbone 14 est un isotope radioactif produit en continu dans la haute atmosphère et qui possède une demi-vie $t_{1/2} \approx 5730$ ans. Tant qu’un organisme (végétal ou animal) est vivant, il échange du carbone avec son environnement (photosynthèse, alimentation). Le rapport entre le carbone 14 instable et le carbone 12 stable reste donc constant au cours de sa vie, maintenant une activité spécifique stable $A_0 \approx 13{,}6$ désintégrations par minute et par gramme de carbone.

À la mort de l’organisme, les échanges cessent. Le carbone 14 piégé dans les tissus commence à se désintégrer sans être renouvelé. L’activité de l’échantillon diminue alors selon la loi exponentielle.

Pour dater un fossile ou un vestige archéologique en bois, on mesure son activité spécifique actuelle $A(t)$ et on la compare à l’activité de référence $A_0$ d’un organisme vivant équivalent. En inversant la loi de décroissance :

$$\frac{A(t)}{A_0} = e^{-\lambda t} \implies \ln\left(\frac{A(t)}{A_0}\right) = -\lambda t \implies t = -\frac{1}{\lambda} \ln\left(\frac{A(t)}{A_0}\right)$$

En remplaçant $\lambda$ par sa valeur en fonction de $t_{1/2}$, on obtient la formule de datation :

$$t = -\frac{t_{1/2}}{\ln(2)} \cdot \ln\left(\frac{A(t)}{A_0}\right)$$

Faisons une application numérique. Supposons que l’analyse d’un morceau de charbon de bois issu d’une grotte préhistorique révèle une activité résiduelle en carbone 14 égale à $30\ %$ de celle d’un arbre vivant aujourd’hui (soit $A(t)/A_0 = 0{,}30$). Son âge $t$ est estimé par :

$$t = -\frac{5730}{\ln(2)} \cdot \ln(0{,}30) \approx -8266 \times (-1{,}204) \approx 9952\text{ ans}$$

Le bois date donc d’environ $10,000$ ans avant notre ère.

---

Conclusion

En conclusion, la modélisation de la décroissance radioactive illustre la symbiose parfaite entre physique et mathématiques. Une équation différentielle linéaire simple permet de décrire le comportement temporel statistique d’une population de noyaux pourtant régie par le hasard à l’échelle individuelle. C’est cette loi de décroissance exponentielle qui sert de clé de lecture temporelle aux archéologues et géologues pour cartographier notre histoire.

Cependant, cette loi s’applique aussi à des problématiques modernes complexes, comme la gestion des déchets nucléaires issus des centrales. Certains isotopes, comme le plutonium 239 ($t_{1/2} \approx 24,000$ ans), exigent un confinement sur des échelles de temps qui dépassent l’histoire humaine écrite.

Pour mon projet d’orientation, je souhaite intégrer une classe préparatoire MPSI afin de poursuivre des études de physique et de modélisation mathématique. Les outils d’analyse différentielle que nous étudions cette année ne sont qu’une introduction à la modélisation des grands systèmes physiques.

Je vous remercie pour votre attention et je me tiens prêt pour l’échange.

---

Questions Potentielles du Jury

1. Démontrez par le calcul que la fonction $N(t) = N_0 \cdot e^{-\lambda t}$ est bien l’unique solution de l’équation différentielle $\frac{dN(t)}{dt} + \lambda \cdot N(t) = 0$ vérifiant $N(0) = N_0$.

Réponse :

1. Vérification qu’elle est solution : Calculons la dérivée temporelle de $N(t)$ : $$\frac{dN(t)}{dt} = \frac{d}{dt}(N_0 \cdot e^{-\lambda t}) = N_0 \cdot (-\lambda \cdot e^{-\lambda t}) = -\lambda \cdot (N_0 \cdot e^{-\lambda t}) = -\lambda \cdot N(t)$$ En injectant cette expression dans l’équation différentielle, on obtient bien : $$\frac{dN(t)}{dt} + \lambda \cdot N(t) = -\lambda \cdot N(t) + \lambda \cdot N(t) = 0$$ La fonction est donc bien solution.
2. Vérification de la condition initiale : Pour $t = 0$ : $$N(0) = N_0 \cdot e^{-\lambda \cdot 0} = N_0 \cdot e^0 = N_0 \cdot 1 = N_0$$ La condition initiale est respectée.
3. Preuve de l’unicité : Supposons qu’il existe une autre solution $f(t)$ de cette équation différentielle. Posons la fonction auxiliaire $g(t) = f(t) \cdot e^{\lambda t}$. Dérivons $g(t)$ : $$g’(t) = f’(t) \cdot e^{\lambda t} + f(t) \cdot \lambda e^{\lambda t} = (f’(t) + \lambda f(t)) \cdot e^{\lambda t}$$ Or, $f$ étant solution de l’équation différentielle, $f’(t) + \lambda f(t) = 0$. Donc $g’(t) = 0 \cdot e^{\lambda t} = 0$. La dérivée de $g$ étant nulle sur $\mathbb{R}$, $g$ est une fonction constante, notée $C$. On en déduit : $$f(t) \cdot e^{\lambda t} = C \implies f(t) = C \cdot e^{-\lambda t}$$ En utilisant la condition initiale $f(0) = N_0$, on trouve $C = N_0$, d’où $f(t) = N_0 \cdot e^{-\lambda t}$. La solution est donc unique.

2. Quelle est la signification de la constante de temps $\tau$ ? Comment se définit-elle mathématiquement et graphiquement ?

Réponse : La constante de temps $\tau$ (exprimée en secondes) est définie comme l’inverse de la constante radioactive $\lambda$ : $$\tau = \frac{1}{\lambda}$$ Elle caractérise la vitesse de décroissance de la population radioactive.

• Signification physique : Au bout d’une durée $t = \tau$, le nombre de noyaux restants vaut : $$N(\tau) = N_0 \cdot e^{-\lambda \tau} = N_0 \cdot e^{-1} \approx 0{,}37 \cdot N_0$$ Il reste environ $37\ %$ de noyaux initiaux, ce qui signifie que $63\ %$ se sont désintégrés.
• Détermination graphique : Graphiquement, sur la courbe représentant $N(t)$ en fonction du temps, la tangente à la courbe à l’instant initial $t=0$ a pour équation : $$y = N’(0) \cdot (t - 0) + N(0) = -\lambda N_0 t + N_0$$ Cette tangente coupe l’axe des abscisses (où $y=0$) à l’instant $t$ tel que : $$0 = -\lambda N_0 t + N_0 \implies \lambda N_0 t = N_0 \implies t = \frac{1}{\lambda} = \tau$$ Ainsi, la constante de temps $\tau$ est l’abscisse du point d’intersection de la tangente à l’origine avec l’axe des abscisses.

3. Pourquoi ne peut-on pas utiliser la méthode du carbone 14 pour dater des vestiges vieux de plus de $50,000$ ans ? Quelle autre méthode préconise-t-il pour des durées plus longues ?

Réponse : La demi-vie du carbone 14 est d’environ $5730$ ans. Au bout de $50,000$ ans, le temps écoulé correspond à environ $9$ demi-vies ($50,000 / 5730 \approx 8{,}7$). La population de carbone 14 est alors divisée par $2^9 = 512$, ce qui signifie qu’il reste moins de $0{,}2\ %$ du carbone 14 initial. L’activité résiduelle devient si faible qu’elle se confond avec le bruit de fond des appareils de mesure, rendant la datation imprécise ou impossible. Pour dater des échantillons beaucoup plus anciens (comme des roches volcaniques ou des fossiles de dinosaures âgés de plusieurs millions d’années), on utilise d’autres couples de radio-isotopes à demi-vie beaucoup plus longue :

• Le couple Uranium-Plomb ($^{238}U \to ^{206}Pb$) avec une demi-vie de $4{,}5$ milliards d’années (très utilisé en géologie pour estimer l’âge de la Terre).
• Le couple Potassium-Argon ($^{40}K \to ^{40}Ar$) avec une demi-vie de $1{,}25$ milliard d’années, idéal pour dater les roches volcaniques anciennes.

4. Quelle est la différence entre l’activité d’un échantillon (en Becquerels) et la dose absorbée (en Grays et Sieverts) ?

Réponse :

1. L’activité (en Becquerels, Bq) : Elle quantifie uniquement la source radioactive. $1\text{ Bq}$ correspond à $1$ désintégration spontanée par seconde, quel que soit le type de rayonnement émis ($\alpha$, $\beta$, $\gamma$).
2. La dose absorbée (en Grays, Gy) : Elle mesure l’énergie déposée par les rayonnements par unité de masse de matière cible (tissus vivants par exemple). $1\text{ Gy}$ correspond à un transfert d’énergie de $1\text{ Joule}$ dans $1\text{ kg}$ de matière.
3. La dose équivalente ou efficace (en Sieverts, Sv) : Elle prend en compte la dangerosité biologique réelle. Tous les rayonnements ne causent pas les mêmes dégâts à énergie égale (les particules $\alpha$ sont beaucoup plus destructrices pour les cellules que les photons $\gamma$). On multiplie la dose en Grays par un facteur de pondération biologique pour obtenir la dose en Sieverts. Le Sievert est l’unité utilisée pour la radioprotection et l’évaluation des risques sur la santé humaine.

5. Est-ce que la décroissance d’un noyau radioactif dépend de sa température ou de son état physique solide ou gazeux ?

Réponse : Non, la décroissance radioactive est un phénomène exclusivement nucléaire, qui se déroule au niveau du noyau de l’atome. Les forces impliquées dans les réactions chimiques, les changements d’état ou l’agitation thermique sont d’origine électromagnétique et lient les électrons périphériques de l’atome (dont l’énergie typique est de l’ordre de l’électronvolt, $eV$). Les énergies de liaison nucléaires (liées à l’interaction forte et l’interaction faible) sont de l’ordre du méga-électronvolt ($MeV$), soit un million de fois supérieures. Par conséquent, chauffer un échantillon, le soumettre à une pression gigantesque, ou le faire réagir chimiquement ne modifie en rien la structure interne du noyau instable, ni sa constante radioactive $\lambda$. La demi-vie d’un isotope est une constante physique universelle et absolue.