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.
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.
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.
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
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.
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 .\]
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)}.\]
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.
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)
- 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). ↩
- T. Erneux et P. Glorieux. « Laser Dynamics », Camb. Univ. Press, Cambridge, UK (2010). ↩
- http://en.wikipedia.org/wiki/Goldfinger_(film) ↩