• Accueil
  • À propos
  • Accrom\(\alpha\)th en PDF
  • Commanditaires
  • Contact
  • Contributions des lecteurs
  • Sites amis

Logo

Un éclairage mathématique sur la dynamique des lasers

Par Thomas Erneux et Pietro-Luciano Buono
Volume 10.1 - hiver-printemps 2015

Le premier laser a été réalisé en 1960. Il apparaît dans notre vie de tous les jours en 1974 lorsqu’il devient un moyen de lecture avec l’introduction des lecteurs de codes barres. En 1978, les disques laser sont introduits, mais les disques optiques ne deviennent d’usage courant qu’en 1982 avec le disque compact. Le laser permet alors de lire un grand volume de données.

Beaucoup de chercheurs qui étudient des lasers ont été confrontés à l’apparition de fluctuations rapides d’intensité dans le faisceau de sortie du laser. Ce type de comportement était déjà évident lors des toutes premières recherches. L’intensité de la lumière produite par le laser à rubis1 montrait des impulsions erratiques comme le montre la figure suivante.

laser-1

Ces impulsions étaient-elles le résultat des fluctuations de l’environnement mécanique ou thermique ou bien étaient-elles le résultat d’une propriété intrinsèque du laser? Après de nombreuses années de recherche, nous comprenons maintenant que ces oscillations proviennent principalement de forces qui régissent l’interaction entre le rayonnement et la matière. Les lasers que nous utilisons aujourd’hui sont parfois instables et certains sont très sensibles à la moindre perturbation extérieure.

Ces instabilités n’ont pas d’effet sur un lecteur de codes barres ou un lecteur de CD. Elles peuvent cependant être une nuisance pour la communication par des fibres optiques car le faisceau laser, devant parcourir des kilomètres, doit être très stable. Aujourd’hui les ingénieurs et les physiciens savent contrôler ces instabilités et même les utiliser pour des applications.

Nous disposons également de modèles mathématiques simples qui permettent de guider nos expériences. Ceux-ci se fondent sur les équations cinétiques du laser2.

James Bond, interprété par Sean Connery, est allongé sur une table et risque l'émasculation par un rayon qui remonte depuis ses pieds jusqu'à son entrejambe.

James Bond, interprété par Sean Connery, est allongé sur une table et risque l’émasculation par un rayon qui remonte depuis ses pieds jusqu’à son entrejambe.

Une première au Cinéma

Le film Goldfinger (1964) était le premier film à montrer un laser. Les lasers n’existaient pas en 1959 quand le livre a été écrit, pas plus que les lasers industriels de forte puissance au moment où le film a été réalisé.

Dans le roman, Goldfinger utilise une scie circulaire mais les cinéastes ont décidé de la changer par un laser. Le faisceau laser lui-même était un effet d’optique ajouté en post-production. Pour les gros plans où la flamme traverse la table métallique, un technicien tient un chalumeau sous la table3.

Modélisation de la dynamique

Une modélisation utilisant une loi de conservation est une modélisation où la variation d’une quantité X dans un intervalle de temps ∆T correspond à la différence entre le gain et la perte de X pendant cet intervalle.

Fig.Laser01

Pour obtenir une modélisation de la dynamique du laser utilisant une telle loi de conservation, considérons le cas où X = n est le nombre de photons du laser et posons N = N2 – N1 où N1 est le nombre d’atomes dans l’état fondamental et N2 le nombre d’atomes dans l’état excité.

Le gain de photons se fait par interaction matière-lumière et est égal à

GNn∆T

où G est le coefficient de gain. La perte de photons, quant à elle, est due à un ensemble de facteurs (par exemple, la dispersion) et vaut, dans notre modèle,

\[\frac{n}{T_c} \Delta T\]

où le paramètre Tc est le temps que les photons passent dans la cavité du laser.

La variation de n est alors

Fig.Laser02

En divisant cette équation par \(\Delta T\) et en prenant la limite \(\Delta T \to 0,\) nous obtenons l’équation différentielle pour n

\[\frac{dn}{dT}=GNn-\frac{n}{T_c}.\]

D’autre part, l’équation différentielle pour N est

\[\frac{dn}{dT}=R_p – \frac{1}{T_1}N-GNn\]

où T1 est le temps caractéristique des atomes excités pour revenir à l’équilibre et Rp est le paramètre de pompe qui correspond à la source d’énergie extérieure.

Par des changements d’échelles des variables d’état et de temps, nous transformons ces deux équations sous leur forme sans dimension:

\[\left\{ \begin{array} {r c l}\displaystyle \frac{dl}{dt}&=&(D-1)I \\ \displaystyle \frac{dD}{dt}&=& \gamma(A-D-DI) \end{array} \right .\]

où \(I = 2GT_1n\) et \(D = GTN\) sont les nouvelles variables d’état, \(t = T /T_c\) est la nouvelle variable de temps et

\[A = GT_cR_pT_1\]

ainsi que

\[\gamma = T_c /T_1\]

sont les nouveaux paramètres.

Ce sont deux équations qui décrivent l’évolution de la densité des photons dans la cavité laser (intensité I) et celle des atomes excités (D) qui produisent ces photons.

Ces équations différentielles ordinaires sont non linéaires à cause des termes \(DI.\) Deux paramètres apparaissent dans ces équations, \(A\) est directement proportionnel à \(R_p\) le paramètre de pompe et est le paramètre de contrôle. Il est, par exemple, proportionnel au courant électrique appliqué à un laser à semi-conducteur. Le laser est OFF (ON) lorsque \(A < 1 (A > 1).\) Le second paramètre \(\gamma=T_c/T_1\) est le rapport entre les deux échelles de temps. Pour la plupart des lasers utilisés aujourd’hui (CO2, état solide, ou semi-conducteur), \(\gamma\) est typiquement de l’ordre de \(10^{-3}.\)

État stationnaire

Un état stationnaire d’un système d’équations différentielles

\[\left\{ \begin{array} {r c l}\displaystyle \frac{dx}{dt}&=&f(x,y) \\ \displaystyle \frac{dy}{dt}&=& g(x,y) \end{array} \right .\]

est une solution où chaque composante est constante: \((x(t); y(t)) = (x_0; y_0).\) Ceci implique

\[\left\{ \begin{array} {r c l}\displaystyle \frac{dx(t)}{dt} &=&\displaystyle \frac{dx_0}{dt}=0\\ \displaystyle \frac{dy(t)}{dt}&=& \displaystyle \frac{dy_0}{dt}=0 \end{array} \right .\]

On peut donc trouver les états stationnaires par la résolution du système d’équations

\[\left\{ \begin{array} {r c l}f(x,y)&=& 0 \\ g(x,y)&=& 0 \end{array} \right .\]

Nous pouvons déterminer les états stationnaires (voir l’encadré « État stationnaire »). Ceux-ci vérifient les équations

\[\left\{ \begin{array} {r c l}\displaystyle \frac{dl}{dt}&=&(D-1)I \\ \displaystyle \frac{dD}{dt}&=& \gamma(A-D-DI) \end{array} \right .\]

avec dI /dt = 0 et dD/dt = 0. Nous avons ainsi deux équations algébriques pour I et D

\[\left\{ \begin{array} {r c l}(D -1)I &=&0\\A – D – D I &= &0\end{array} \right .\]

La première équation est satisfaite si (a) I = 0 ou si (b) D = 1. En utilisant la deuxième équation, nous trouvons l’expression de la variable manquante. Les états stationnaires sont donc:

a) I = 0 et D = A;
b) I = A – 1 ≥ 0 et D = 1.

L’inégalité dans b) est imposée pour des raisons physiques: une intensité (nombre de photons) est toujours positive.

Nous pouvons ensuite analyser la stabilité linéaire de ces états stationnaires en déterminant comment une toute petite perturbation va évoluer dans le temps. Si la perturbation évolue vers zéro, nous disons que l’état est stable.

Si au contraire, la perturbation croît dans le temps, nous parlons d’un état instable. Cette analyse se fait à partir des équations originales

\[\left\{ \begin{array} {r c l}\displaystyle \frac{dl}{dt}&=&(D-1)I \\ \displaystyle \frac{dD}{dt}&=& \gamma(A-D-DI) \end{array} \right .\]

(voir l’encadré « Calcul de la stabilité linéaire ») et les résultats sont résumés dans la figure suivante.

laser-2

Calcul de la stabilité linéaire

Calculons la stabilité linéaire de l’état stationnaire I = 0 et D = A. Posons

Ĩ (t)=0+α(t) et \(\tilde{D}\) (t)=A+β(t),

où α(t) et β(t) sont les petites perturbations dépendant du temps. Nous allons remplacer I par Ĩ et D par \(\tilde{D}\) dans les deux équations:

\[\left\{ \begin{array} {r c l}\displaystyle \frac{dl}{dt}&=&(D-1)I \\ \displaystyle \frac{dD}{dt}&=& \gamma(A-D-DI) \end{array} \right .\]

On obtient alors:

\[\begin{array} {r c l}\displaystyle \frac{d\tilde{I}(t)}{dt} &=& \displaystyle \frac{d \alpha (t)}{dt} \\ &=& (A + \beta (t)-1) \alpha (t), \\ \displaystyle \frac{d\tilde{D}(t)}{dt} &=& \displaystyle \frac{d \beta (t)}{dt} \\ &=& – \gamma (A \alpha (t)+ \beta (t)+ \alpha (t) \beta (t)).\end{array} \]

En gardant uniquement les termes linéaires (c.à.d. de degré 1) nous avons le système

\[\begin{array} {r c l}\displaystyle \frac{d\alpha(t)}{dt} &=& (A-1) \alpha (t), \\ \displaystyle \frac{d\beta (t)}{dt} &=& – \gamma (A \alpha(t)+\beta(t)). \end{array} \]

L’équation pour \(\alpha (t)\) a comme solution \(\alpha (t) = \alpha_{0}e^{(A-1)t}\) (où \(\alpha_0 ≠ 0\) est la valeur initiale en t = 0 de \(\alpha (t)\)) et nous voyons donc que si A < 1 la perturbation décroit vers 0 lorsque t augmente.

Par contre, si A > 1, alors la perturbation augmente avec le temps et on dit que la solution stationnaire (0, A) est instable. L’équation en \(\beta\) a comme solution

\[\beta(t)= \left ( \beta_0 + \frac{\alpha_0 A \gamma}{A-1 + \gamma} \right ) e^{-\gamma t} – \frac{\alpha_0 A \gamma}{A-1 + \gamma} e^{A-1)t},\]

(où \(\beta (0) = \beta 0).\) Sa croissance/décroissance dépend aussi de la valeur de A de la même manière que \(\alpha(t).\)

Pour les physiciens, le point A = 1 de la figure ci-contre est appelé seuil du laser: Un laser n’est actif que si le paramètre de pompe dépasse la valeur critique A = 1. Le courant électrique d’un laser à semi-conducteur doit être suffisamment élevé pour que les atomes soient dans un état excité. Pour les mathématiciens, le point A = 1 est appelé point de bifurcation parce que la solution d’intensité nulle échange sa stabilité avec la solution non nulle. Le diagramme de cette figure est appelé diagramme des bifurcations. Nous apprenons de ce diagramme que l’état « ON » est toujours stable. Mais cette observation est trompeuse. La figure « Oscillations décroissantes pour I(t) » montre une simulation numérique des équations

\[\left\{ \begin{array} {r c l}\displaystyle \frac{dl}{dt}&=&(D-1)I \\ \displaystyle \frac{dD}{dt}&=& \gamma(A-D-DI) \end{array} \right .\]

laser-3

Suite à une faible perturbation de l’état stationnaire « ON », nous observons des oscillations de grande amplitude qui décroissent lentement avec le temps. La raison de ce comportement dramatique est la faible valeur de \(\gamma\) qui force nos deux variables à opérer sur deux échelles de temps distinctes. L’état stationnaire « ON » est bien stable mais une petite perturbation conduit à une évolution violente qui explique la haute sensibilité des lasers à la moindre perturbation. Ces oscillations sont appelées « oscillations de relaxation » par le physicien, voir figure ci-dessous.

Une quantité précieuse à mesurer est la fréquence de ces oscillations parce que l’analyse de stabilité linéaire nous propose une bonne approximation donnée par

\[\omega=\sqrt{\gamma(A-1)}.\]

laser-4

En mesurant la fréquence pour différentes valeurs de A; nous pouvons déterminer \(\gamma\) comme la pente de la droite \(\omega^2=\gamma(A–1)\) de la figure suivante.

laser-5

L’étude mathématique de la dynamique des lasers est un sujet très actif. Elle utilise et stimule des découvertes en théorie des systèmes dynamiques, entre autres pour les équations à retard. De plus, c’est un champ d’étude où les modélisateurs travaillent en proche collaboration avec les chercheurs en laboratoire, créant ainsi des interactions bénéfiques pour les deux approches scientifiques. Les multiples applications du laser dans la vie de tous les jours et dans des domaines de technologie de pointe vont continuer de générer des défis pour les chercheurs en dynamique des lasers.

Pour en s\(\alpha\)voir plus!

COLLINS R.J., NELSON D.F., SCHAWLOW A.L., BOND W., GARETT C.G.B., and KAISER W.,
Coherence, narrowing, directionality, and relaxation oscillations in the light emission from ruby, Phys. Rev. Lett. 5, 303-305 (1960)
NELSON D.F. and BOYLE W.S., A continuously operating ruby optical laser, Appl. Opt. 1, 181-183 (1962) ERNEUX T. and GLORIEUX P., « Laser Dynamics », Camb. Univ. Press, Cambridge UK (2010)

PDF

  1. Un laser à rubis est un laser qui utilise un rubis synthétique comme milieu amplificateur solide. Le tout premier laser opérationnel était un laser à rubis fabriqué le 16 mai 1960 par Theodore Maiman au Hughes Research Laboratories à Malibu (Californie). ↩
  2. T. Erneux et P. Glorieux. « Laser Dynamics », Camb. Univ. Press, Cambridge, UK (2010). ↩
  3. http://en.wikipedia.org/wiki/Goldfinger_(film) ↩
  • ● Version PDF
Partagez
  • tweet

Tags: Lumière

Articles récents

  • Le mouvement brownien : Du pollen de Brown à l’origine de la finance moderne

    Michel Adès, Matthieu Dufour, Steven Lu et Serge Provost
  • Le problème des \(N\) corps

    Christiane Rousseau
  • Comprendre la structure des nombres premiers

    Andrew Granville

Sur le même sujet

  • La lumière: un éclairage moderne

    André Ross
  • Construire un cadran solaire

    Christiane Rousseau
  • Pleins feux sur la lumière

    André Ross

Volumes

  • Volume 18.1 – hiver-printemps 2023
  • Volume 17.2 – été-automne 2022
  • Volume 17.1 – hiver-printemps 2022
  • Journée internationale des mathématiques: Accromath multilingue
  • Volume 16.2 – été-automne 2021
  • Volume 16.1 – hiver-printemps 2021
  • Volume 15.2 – été-automne 2020
  • Thème spécial: Les mathématiques sont partout
  • Volume 15.1 – hiver-printemps 2020
  • Volume 14.2 – été-automne 2019
  • Volume 14.1 – hiver-printemps 2019
  • Volume 13.2 – été-automne 2018
  • Volume 13.1 – hiver-printemps 2018
  • Volume 12.2 – été-automne 2017
  • Volume 12.1 – hiver-printemps 2017
  • Volume 11.2 – été-automne 2016
  • Volume 11.1 – hiver-printemps 2016
  • Volume 10.2 – été-automne 2015
  • Volume 10.1 – hiver-printemps 2015
  • Volume 9.2 – été-automne 2014
  • Volume 9.1 – hiver-printemps 2014
  • Volume 8.2 – été-automne 2013
  • Volume 8.1 – hiver-printemps 2013
  • Volume 7.2 – été-automne 2012
  • Volume 7.1 – hiver-printemps 2012
  • Volume 6.2 – été-automne 2011
  • Volume 6.1 – hiver-printemps 2011
  • Volume 5.2 – été-automne 2010
  • Volume 5.1 – hiver-printemps 2010
  • Volume 4.2 – été-automne 2009
  • Volume 4.1 – hiver-printemps 2009
  • Volume 3.2 – été-automne 2008
  • Volume 3.1 – hiver-printemps 2008
  • Volume 2.2 – été-automne 2007
  • Volume 2.1 – hiver-printemps 2007
  • Volume 1 – été-automne 2006
  • Article vedette

    Auteurs

    • Michel Adès
    • Antoine Allard
    • Jean Aubin
    • Marie Beaulieu
    • Rosalie Bélanger-Rioux
    • Claude Bélisle
    • Marc Bergeron
    • Pierre Bernier
    • André Boileau
    • Véronique Boutet
    • Pietro-Luciano Buono
    • Massimo Caccia
    • Jérôme Camiré-Bernier
    • France Caron
    • Philippe Carphin
    • Kévin Cazelles
    • Laurent Charlin
    • Pierre Chastenay
    • Noémie Chenail
    • Jocelyn Dagenais
    • Marie-France Dallaire
    • Jean-Lou de Carufel
    • Jean-Marie De Koninck
    • Lambert De Monte
    • Jean-Paul Delahaye
    • Marc-André Desautels
    • Florin Diacu
    • Jimmy Dillies
    • Nicolas Doyon
    • Philippe Drobinski
    • Hugo Drouin-Vaillancourt
    • Louis J. Dubé
    • Thierry Duchesne
    • Matthieu Dufour
    • Stéphane Durand
    • Thomas Erneux
    • Philippe Etchécopar
    • Julien Fageot
    • Charles Fleurent
    • Jérôme Fortier
    • Marlène Frigon
    • Jean-François Gagnon
    • André Garon
    • Christian Genest
    • Denis Gilbert
    • Jonathan Godin
    • Frédéric Gourdeau
    • Samuel Goyette
    • Andrew Granville
    • Jean Guérin
    • Hervé Guillard
    • Abba B. Gumel
    • James A. Hanley
    • Alain Hertz
    • Bernard R. Hodgson
    • Isabelle Jalliffier-Verne
    • Guillaume Jouvet
    • Tomasz Kaczynski
    • Patrick Labelle
    • Marc Laforest
    • Nadia Lafrenière
    • Josiane Lajoie
    • Alexis Langlois-Rémillard
    • Simon-Olivier Laperrière
    • René Laprise
    • Steffen Lauritzen
    • Denis Lavigne
    • Adrien Lessard
    • Steven Lu
    • Jean Meunier
    • Erica Moodie
    • Normand Mousseau
    • Johanna G. Nešlehová
    • Pierre-André Noël
    • Dmitry Novikov
    • Ostap Okhrin
    • Laurent Pelletier
    • Jean-François Plante
    • Serge B. Provost
    • Annie Claude Prud'Homme
    • Benoît Rittaud
    • Louis-Paul Rivest
    • Serge Robert
    • André Ross
    • Christiane Rousseau
    • Guillaume Roy-Fortin
    • Yvan Saint-Aubin
    • Maria Vittoria Salvetti
    • Charles Senécal
    • Vasilisa Shramchenko
    • Robert Smith?
    • Anik Trahan
    • Shophika Vaithyanathasarma
    • William Verreault
    • Redouane Zazoun

Sujets

Algèbre Applications Applications des mathématiques Changements climatiques Climat Construction des mathématiques COVID-19 Cristallographie cryptographie GPS Gravité Géométrie Histoire des mathématiques Imagerie Infini Informatique Informatique théorique intelligence artificielle Jeux mathématiques Logique mathématique Lumière Mathématiques de la planète Terre Mathématiques et architecture mathématiques et art Mathématiques et arts Mathématiques et astronomie Mathématiques et biologie Mathématiques et développement durable Mathématiques et littérature Mathématiques et musique Mathématiques et médecine Mathématiques et physique Mathématiques et transport Modélisation Nombres Portrait d'un mathématicien Portrait d'un physicien Probabilités Probabilités et statistique Racines Rubrique des Paradoxes Section problèmes Théorie des groupes Éditorial Épidémiologie

© 2023 Accromath