Oscillateur harmonique quantique

L'oscillateur harmonique quantique correspond au traitement par les outils de la mécanique quantique de l'oscillateur harmonique classique.

De façon générale, un oscillateur est un système dont l'évolution dans le temps est périodique. Il est dit de plus harmonique si les oscillations effectuées sont sinusoïdales, avec une amplitude et une fréquence qui ne dépendent que des caractéristiques intrinsèques du système et des conditions initiales. Cela est le cas en mécanique classique pour une particule évoluant à une dimension dans un potentiel quadratique, de forme générale , k étant une constante positive.

Loin d'être un cas particulier purement académique, cette forme de potentiel est obtenue notamment dans le cas d'oscillations de faible amplitude autour d'une position d'équilibre stable dans un potentiel quelconque, à l'exception notoire du potentiel en 1/r rencontré en gravitation ou pour des charges ponctuelles comme celle des protons et électrons (cf. article sur l'oscillateur harmonique classique), car au voisinage de cette position d'équilibre le potentiel prend cette forme. Pour cette raison, le concept d'oscillateur harmonique joue un rôle majeur dans de nombreuses applications de la physique.

La mécanique quantique a révolutionné un grand nombre de concepts fondamentaux. L'oscillateur harmonique a aussi subi une reformulation dans ce cadre quantique, ce qui a permis d'élucider plusieurs résultats expérimentaux, notamment en physique de la matière condensée. Son étude amène à introduire des outils mathématiques d'un intérêt considérable en physique, notamment en théorie des champs : les opérateurs de création et d'annihilation de quanta.

Evolution temporelle de l'oscillateur harmonique d'après les lois de Newton de la mécanique classique (A–B), et d'après la mécanique quantique (C–H). En A et B, la particule (représentée par une balle attachée à un ressort) est animée d'un mouvement d'oscillation autour de sa position d'équilibre. Sur les figures C-H, l'axe horizontal indique la position, et l'axe vertical indique les parties réelle (bleu) et imaginaire (rouge) de la fonction d'onde. Les états C, D, E, F sont des états stationnaires, d'énergie fixée, pour lesquels la probabilité de trouver la particule dans une région donnée n'évolue pas au cours du temps. G est une superposition aléatoire de ces quatre états (C-F). H est un état cohérent dont le comportement est proche d'un oscillateur classique.

L'oscillateur harmonique classique à une dimension

modifier

Avant d'aborder le traitement quantique de l'oscillateur harmonique, il est intéressant de résumer les principaux résultats du traitement classique de cet oscillateur.

Définition - Importance physique

modifier

Un oscillateur harmonique classique à une dimension est modélisé par un potentiel quadratique, typiquement (m étant la masse du système) :

 
Comparaison pour une molécule diatomique entre la courbe de potentiel « réelle », représentée par le potentiel de Morse et celle d'un oscillateur harmonique. Le caractère non-harmonique du potentiel réel conduit à un resserrement des niveaux d'énergie, qui sont également espacés pour un oscillateur harmonique « pur », cf. plus bas dans l'article.

 , où   est une grandeur homogène à une fréquence (angulaire) appelée pulsation propre de l'oscillateur. De façon générale celle-ci est liée à la « force » ou « raideur » du potentiel[1].

L'énergie mécanique d'un oscillateur harmonique unidimensionnel est donc  .

L'oscillateur harmonique a un intérêt considérable en physique, car tout système évoluant dans un potentiel au voisinage d'une position d'équilibre stable, donc un minimum de potentiel, peut être modélisé par un oscillateur harmonique pour les petites oscillations au voisinage de cette position d'équilibre.

En effet, dans le cas d'un potentiel quelconque unidimensionnel  , avec un minimum en  , il vient au voisinage de celui-ci :

 , avec nécessairement pour un minimum  ,

ce qui correspond bien à un oscillateur harmonique tel que  , l'expression précédente du potentiel correspondant au choix de   et   comme origines respectivement des coordonnées et des énergies.

L'approximation « harmonique » est valable au voisinage de la position d'équilibre stable du système : la figure ci-contre illustre la modélisation par un oscillateur harmonique de la courbe d'énergie potentielle d'une molécule diatomique.

Équations classiques du mouvement

modifier

Pour la discussion des états quasi-classiques, il est utile de rappeler des éléments sur la nature du mouvement classique.

Les équations du mouvement de Hamilton pour le hamiltonien classique précédent s'écrivent sans difficultés sous la forme :

 .

Il est possible de donner une forme plus symétrique à ces équations en mettant le hamiltonien précédent sous une forme « réduite », en introduisant de nouvelles variables canoniques :

  et  ,

ce qui permet de mettre le hamiltonien sous la forme :

 .

Comme   la transformation   est canonique, et donc laisse inchangée la forme des équations de Hamilton. Les équations du mouvement précédentes deviennent alors :

 .

L'introduction de la grandeur complexe  , permet de regrouper ces deux équations en une équation unidimensionnelle :

 ,

dont la solution générale est de la forme  , avec  . Par suite une particule de masse   dans ce potentiel a un mouvement sinusoïdal de pulsation  . L'énergie de la particule est une constante du mouvement positive ou nulle  , dont la valeur dépend des conditions initiales de façon continue, sans restriction.

Le hamiltonien « réduit » précédent s'écrit en fonction de   sous la forme  .

L'oscillateur harmonique quantique à une dimension

modifier

Position du problème

modifier

En mécanique quantique l'hamiltonien d'un oscillateur harmonique unidimensionnel s'écrit sous la forme :

 

  est la composante sur l'axe   de l'opérateur impulsion de la particule.

Il faut alors résoudre l'équation de Schrödinger indépendante du temps associée à cet hamiltonien :

 , où   est l'énergie associée à un état propre   du système[2].

Il est alors possible d'écrire cette équation en représentation position, compte tenu de   et  , ce qui donne l'équation différentielle suivante pour la fonction d'onde   :

 .

La résolution mathématique de cette équation ne présente pas de difficultés majeures, bien qu'elle conduise à des calculs assez complexes. Toutefois une telle méthode est surtout peu explicite physiquement, aussi est-il préférable d'employer une autre approche, très féconde, développée par Paul Dirac, qui non seulement permet d'obtenir les valeurs propres du hamiltonien sans résoudre explicitement l'équation différentielle précédente, mais qui peut être généralisée à d'autres situations (cf. quantification du champ électromagnétique classique, ou étude des vibrations dans un cristal, par exemple).

Propriétés générales des états propres

modifier

Avant toute résolution, il est possible de déduire certaines propriétés importantes des valeurs et états propres du hamiltonien quantique précédent.

  • Les valeurs propres E de   sont toujours positives ou nulles : en effet pour un état   quelconque, il est évident que   or  .
De façon évidente, il en est de même pour le terme en  , donc pour tout état  ,  .
  • L'énergie de l'état fondamental est nécessairement non nulle : en introduisant les variances des valeurs moyennes de   et de  , soit   et  , il est possible de montrer que  .
Par suite l'état fondamental de l'oscillateur harmonique quantique a une énergie non nulle (énergie du point zéro), contrairement au cas classique, et ceci résulte directement de la relation d'incertitude quantique entre   et  .
  • Le confinement de la particule dans ce potentiel indique que le spectre de ces énergies sera discret.
  • Le système ne fait apparaître qu'un seul degré de liberté, il n'y aura donc qu'un seul nombre quantique.

Introduction d'opérateurs sans dimension

modifier

Comme dans le cas classique il est utile d'introduire un hamiltonien sans dimension, dit réduit. Ceci est aisé dans la mesure où il est aisé d'introduire une « échelle » d'énergie pour le système en utilisant la quantité  , par suite (  est la constante de Planck réduite) :  

Il est facile de simplifier l'écriture de  , en utilisant les opérateurs sans dimension (mais hermitiens) suivants :

 , et  .

Compte tenu de la relation de commutation  ,la relation de commutation entre ces nouveaux opérateurs s'écrit :

 .

En termes de ces opérateurs « réduits » le hamiltonien de l'oscillateur harmonique unidimmensionnel s'écrit sous la forme :  .

Remarques :

  • La quantité   a la dimension d'une longueur.
  • Les opérateurs « réduits » d'énergie, de position et d'impulsion quantiques ont des expressions qui rappellent celles utilisées dans le cas classique, sauf pour l'utilisation de la quantité  .

Opérateurs d'échelle

modifier

En utilisant la relation de commutation entre les opérateurs « réduits » de position et d'impulsion il est facile de vérifier l'identité suivante :

 .

Ceci suggère l'introduction des opérateurs non-hermitiques, appelés « opérateurs d'échelle », adjoints l'un de l'autre :

 

Ces opérateurs sont dits respectivement d’annihilation et de création de quantum d'énergie, pour des raisons qui apparaîtront plus bas. Il est facile de déterminer le commutateur entre ces deux opérateurs :  .

De même, les opérateurs   et   s'expriment sous la forme :  

Le hamiltonien du système peut alors s'écrire sous la forme :

 ,

où il a été introduit l'opérateur :

 , appelé parfois « opérateur nombre de quanta de vibration ».

Par suite, la détermination des valeurs propres du hamiltonien se ramène à celles de l'opérateur sans dimensions  , qui est évidemment hermitien puisque  . Par suite, ses valeurs propres sont réelles.

Remarques :

  • l'opérateur d'échelle   peut être considéré comme le « pendant quantique » de la quantité complexe   définie dans le cas classique, partie 1.
  • il est intéressant alors de comparer l'expression du hamiltonien classique « réduit » avec celle obtenu dans le cas quantique. En dehors du facteur  , la principale différence vient de l'apparition du terme supplémentaire en 1/2 dans le cas quantique, lié au caractère non commutatif des opérateurs d'échelle (en fait, de   et  ), qui donnera lieu à l'apparition de l'énergie dite de « point zéro ».

Calcul des valeurs propres

modifier

Les états propres commun de   et   sont notés  , avec pour valeurs propres :  , la valeur propre ν étant a priori un réel. Les valeurs propres de   sont alors de la forme :  . Dans la mesure où cela ne change pas les résultats pour les valeurs propres il n'est pas nécessaire dans un premier temps de prendre en compte une éventuelle dégénérescence de ces états propres.

Relations de commutation utiles

modifier

À partir de la relation de commutation précédente entre les opérateurs d'échelle, il vient facilement :

  •  
  •  

Ces diverses relations rendent possible la détermination des valeurs propres de  .

Action des opérateurs d'échelle sur les états propres  

modifier

En utilisant les notations et relations précédentes, il vient facilement :

 , ce qui donne :
 ,
autrement dit   est un vecteur propre de   de valeur propre (ν-1), il est possible de poser :
 , avec α nombre complexe qui peut être déterminé à une phase près en utilisant le fait que :
 ,

et comme nécessairement   il est possible d'en déduire successivement :

  1. que les valeurs propres ν de   sont positives ou nulles;
  2. que si en particulier ν = 0, alors   et inversement si   alors   et donc tout état respectant cette condition est état propre de   de valeur propre ν = 0;
  3. qu'en effectuant un choix de phase tel que α soit réel, il vient :  , soit encore :  ,

par suite l'opérateur   abaisse la valeur propre ν d'une unité, d'où le nom d'opérateur d'annihilation de quantum (de vibration) donné plus haut.

En procédant de même avec   il est possible de vérifier avec la même convention de phase que  , par suite l'opérateur   accroît la valeur propre ν d'une unité, d'où le nom d'opérateur de création de quantum (de vibration) donné plus haut.

Valeurs propres de l'opérateur N

modifier

Soit un état propre   avec ν non entier, positif, N étant le premier entier tel que N > ν, il résulte de la relation   que l'application à N reprises de l'opérateur d'annihilation   permettait d'obtenir un état propre de   avec une valeur propre négative, ce qui n'est pas possible puisque les valeurs propres de   sont positives ou nulles.

En revanche si ν = n, n entier strictement positif, n applications successives de   donnent :  , et une nouvelle application de   donne alors le vecteur nul.

Par suite les seules valeurs propres ν de l'opérateur   sont les entiers n positif ou nul. L'appellation d'opérateur « nombre de quanta de vibration » donnée plus haut est donc justifiée.

Les résultats précédents donnent alors :

Les énergies accessibles par l'oscillateur sont :   avec n entier positif ou nul

Ces états ont les propriétés générales suivantes :

  1. Les énergies accessibles par l'oscillateur sont quantifiées. Ce résultat a de nombreuses répercussions en physique statistique par exemple : de fait, il est impossible de décrire correctement les propriétés thermiques des solides (capacité calorifique par exemple) sans tenir compte du caractère quantifié des états de vibrations des atomes qui les constituent.
  2. Les états d'énergie sont espacés de la même quantité  . Ceci permet de considérer que le passage de l'état n à l'état n + p correspond à « l'absorption » de p quanta d'énergie par le système. Cette façon de penser est en fait très féconde : ainsi l'étude des vibrations dans les solides conduit-elle à introduire la notion de phonons, qui se comportent comme des quasi-particules. De même, et de façon plus significative, la quantification du champ électromagnétique conduit à introduire les modes propres de ce champ, appelé « photons ».
  3. Comme cela a déjà mis en évidence qualitativement, l'énergie de l'état fondamental est non nulle, avec  , appelé souvent « énergie du point zéro » de l'oscillateur. Il a déjà été montré que cette situation, très différente du cas classique, résulte de la relation d'incertitude de Heisenberg.

Détermination des états propres

modifier

Non-dégénérescence de l'état fondamental

modifier

En supposant que l'état fondamental est k-fois dégénéré, il est possible de poser   avec λ = 1,2,..., k. Comme pour tout λ  , il est possible de revenir à la définition l'opérateur d'annihilation   et de donner l'équation différentielle correspondante en représentation position que doivent satisfaire les fonctions d'onde   pour tout k :

 .

Cette équation est aisément soluble par séparation des variables et aboutit à une même fonction, à une constante (complexe en général) de normalisation près :  .

La normalisation de ce résultat, permet d'obtenir  .

Par suite, les solutions étant proportionnelles entre elles, l'état fondamental est non dégénéré.

Il est alors facile de voir que tous les états propres du hamiltonien sont non dégénérés : en effet si c'était le cas il serait possible à partir de deux états   et   de même énergie   d'obtenir, par application à n reprises de l'opérateur d'échelle   deux états fondamentaux dégénérés   et   distincts, ce qui est impossible du fait de la non-dégénérescence de l'état fondamental[5].

États propres de l'opérateur N

modifier

Une démonstration par récurrence utilisant l'opérateur de création de quantum de vibration montre que les états propres de   s'écrivent :  .

 
Représentation des sept premiers niveaux d'énergie et des fonctions d'onde associées  ) de l'oscillateur harmonique quantique unidimensionnel.

En représentation position, il suffit de substituer l'expression de   et de   pour obtenir l'expression de la fonction d'onde   sous la forme :

 , avec   position réduite.
 
Représentation similaire à la précédente, avec cette fois-ci les allures des densités de probabilité de présence   pour les premiers états d'énergie de l'oscillateur harmonique unidimensionnel.

En introduisant les polynômes d'Hermite   sous leur forme dite « physique » définie par   il est finalement possible d'obtenir l'expression générale de la fonction d'onde   :

 .

Il est facile de vérifier que ces fonctions d'ondes sont orthogonales entre elles et normalisées à l'unité, en utilisant les relations d'orthogonalité des polynômes d'Hermite. Par ailleurs, la parité de la fonction d'onde   est celle de n.

Cette dernière relation nous permet de retrouver explicitement autant de fonctions d'ondes que nécessaire. Par exemple pour   il vient :
 

Représentation matricielle

modifier
  • La matrice représentative de l'hamiltonien   sur la base des   est, par construction, diagonale.

On a  

  • Sachant que  , en multipliant à gauche par   il vient  

La matrice représentative de   sur la base des   est donc  

Puisque  , la matrice   est construite en transposant   :  

  • Il est alors facile de construire les matrices représentatives des observables   et   puisque :

 et  

donc   et  

Il est alors possible de vérifier les propriétés suivantes pour les observables   et   :

  • Les seuls éléments de matrices non nuls de ces opérateurs sont donc ceux pris entre les états   et  .
  • Les valeurs moyennes de la position et de l'impulsion sont nulles lorsque le système est dans un état propre  [6].

Évolution des valeurs moyennes dans le temps - limite classique

modifier

Dans le cas d'un hamiltonien   qui ne dépend pas explicitement du temps, il résulte de l'équation de Schrödinger   que le vecteur d'état   s'obtient à partir de celui pris à un temps   quelconque par :

 .

Or comme  , les coefficients complexes   étant tels que   est normalisé à l'unité[7], dans le cas d'un oscillateur harmonique l'expression du vecteur d'état à une date t(>t0) quelconque est la suivante :

 .

Par suite, la valeur moyenne de tout observable   est donnée en fonction du temps par :

 , avec  , élément de matrice de l'observable   entre   et  .

En particulier, pour les observables   et   les seuls éléments de matrices non nuls sont ceux entre   et  , par suite, l'évolution de leur grandeurs moyennes fait uniquement intervenir des termes en  , c'est-à-dire des fonctions sinusoïdales : ceci correspond bien à un comportement classique pour un oscillateur harmonique.

Ce dernier résultat peut également s'obtenir à partir du théorème d'Ehrenfest, qui compte tenu du fait que   et   ne dépendent pas explicitement du temps s'écrit pour chaque valeur moyenne :

 ,
 ,

soit finalement le système d'équations :

 .

Pour les valeurs moyennes, ces équations sont donc identiques à celles de l'oscillateur harmonique classique indiquées en première partie. Elles donnent par intégration des évolutions sinusoïdales de pulsation ω pour les valeurs moyennes de   et  .

Interprétation

modifier
 
Densité de probabilité
n(x)|2 pour les premières valeurs de l'énergie, en commençant par le fondamental (n = 0) en bas puis en augmentant la valeur de l'énergie en montant. L'axe horizontal correspond à la coordonnée x, les couleurs plus claires indiquent une densité de probabilité plus forte.

En analysant ces fonctions d'ondes, on retrouve de nombreux résultats classiques : la particule dans le puits de potentiel a une probabilité de présence plus élargie si elle a une énergie plus haute (une bille au fond d'un puits va monter plus haut sur les bords si elle a plus d'énergie), la particule a plus de chance de se retrouver sur ces positions éloignées du centre du puits (la bille a une vitesse d'autant plus petite qu'elle est haut dans le puits : elle va donc passer beaucoup plus de temps en hauteur qu'au fond du puits).

États cohérents

modifier

Les états cohérents, notés  , sont par définition les états propres de l'opérateur d'annihilation  , de valeurs propres   :  .

Propriétés générales des états cohérents  

modifier
Expression de   sur la base des états propres   du hamiltonien
modifier

Il est possible de poser de façon générale pour tout état   :

 , avec  , qui seront choisis de telle sorte que   soit normé à l'unité :  .

D'après la définition d'un état cohérent donnée plus haut, et le fait que pour tout état propre  ,   il vient :

 , où l'on a tenu compte de  .

En procédant à un changement de variable dans le premier membre, l'équation aux valeurs propres s'écrit dans la base de   :

 ,

ce qui donne par identification membre à membre la relation de récurrence suivante entre les coefficients   :

 , qui permet alors d'obtenir l'expression de   en fonction de   :
 .

La constante   peut être obtenue en utilisant la condition de normalisation précédente, qui donne alors  , soit en prenant pour convention de phase que tous les coefficients   sont réels, il vient finalement :  .

Le développement de   sur la base des états propres   s'écrit en définitive sous la forme :  .

Comme  , il est aussi possible d'écrire :

 .

Il est facile de vérifier que si α = 0 ce développement se réduit à   : l'état fondamental de l'oscillateur harmonique est donc un état cohérent, de valeur propre α = 0.

L'état cohérent comme paquet d'onde minimum
modifier

Il est facile de démontrer que les états cohérents   sont paquets d'onde minima, c'est-à-dire qu'ils minimisent les incertitudes (écarts-types en l'espèce) sur la mesure de la position et de l'impulsion : pour tout état cohérent  .

Applications

modifier

Pour tous les calculs où des particules sont dans un puits de potentiel, l'approximation harmonique est très intéressante (en effet, on remarque que le développement limité à l'ordre 2 d'un puits nous donne une parabole). Par exemple si l'on souhaite étudier un « piège harmonique » à deux dimensions (condensation de Bose-Einstein à 2D) on pourra poser le hamiltonien suivant pour débuter l'étude :  

Généralisation à trois dimensions

modifier

Ces calculs pour une seule dimension se généralisent très bien à 3 dimensions. L'hamiltonien est alors simplement sous la forme d'une somme de trois hamiltoniens indépendants, qui sont donc étudiables séparément exactement comme nous l'avons fait précédemment.

Les énergies accessibles par l'oscillateur sont :    

avec   : l'énergie dépend de trois nombres quantiques indépendants. Pour une même énergie, il va donc être possible d'imaginer des configurations différentes : les niveaux d'énergie sont dégénérés.

On calcule le nombre de dégénérescences pour le n-ième niveau d'énergie :  

Oscillateur anharmonique

modifier

Le modèle harmonique d'un puits de potentiel est certes une approche intéressante et finalement assez simple, mais pour étudier les écarts à l'harmonicité d'un système réel, il va falloir pousser le développement limité du potentiel à un ordre supérieur. On introduit alors un terme perturbatif   dans l'hamiltonien puis, à l'aide de la théorie des perturbations stationnaires, on va calculer les nouvelles énergies du système.

Terme en x3

modifier

Si l'on pousse le développement limité du potentiel à l'ordre 3, on obtient le terme perturbatif suivant par rapport à l'oscillateur harmonique :    est très petit par rapport à 1. Le mouvement de la particule dans un tel potentiel n'est plus tout à fait symétrique. Le calcul théorique des énergies donne alors :  

L'utilisation de la méthode des perturbations pour ce calcul impose un résultat approché (correction du deuxième ordre ici). On remarque alors que les niveaux d'énergie ont été abaissés par cette correction. Le calcul des nouvelles fonctions d'onde montre qu'elles sont formées par des couplages entre les anciennes :  

Application

modifier

Une molécule diatomique dans le modèle harmonique peut absorber et émettre des ondes électromagnétiques à la pulsation  . Or dans la pratique, on observe d'autres raies absorbées ou émises : le modèle anharmonique permet de bien en rendre compte, puisque d'après la nouvelle expression des fonctions d'onde, deux états   et   peuvent être reliés même si   et   différent de plus de 1. La constante   peut ainsi être calculée dans cet exemple en analysant le spectre d'absorption.

  1. Plus précisément à la « pente » de la parabole correspondant à la courbe représentative de V(x), voir ci-après.
  2. Conformément à une convention courante l'état propre de l'équation de Schrödinger dépendante du temps s'écrit avec Ψ majuscule :   tandis que celui de l'équation indépendante du temps s'écrit avec ψ minuscule :  . Il en est de même des fonctions d'onde associées dans les diverses représentations.
  3. Ce procédé, qui consiste à prendre une « fonction d'onde d'essai », et d'obtenir à partir de là une expression de  , dépendant de certains paramètres que l'on considère alors comme des variables, pour minimiser  , et donc estimer la valeur de l'état fondamental, est très fécond. Il est à la base de la méthode variationnelle, utilisée par exemple dans la théorie élémentaire des orbitales moléculaires (cf. estimation de l'état fondamental de l'ion moléculaire   ou de façon numérique, par exemple dans la méthode de Monte-Carlo variationnelle).
  4. En fait, comme indiquée plus bas la valeur de l'état fondamental est bien  .
  5. Toutefois, à plusieurs dimensions p, la dégénérescence sera possible. En effet dans ce cas il est possible de montrer que l'on se ramène à p oscillateurs harmoniques indépendants (modes normaux), l'énergie totale étant la somme des énergies associées à chacun de ces oscillateurs soit  . Il est clair qu'il y a plusieurs combinaisons de valeurs possibles des  , donc plusieurs états distincts, donnant une même valeur de l'énergie du système.
  6. Ce qui ne signifie pas que sur un état donné du système de vecteur d'état   ces valeurs moyennes soient nulles, puisque dans le cas général   résulte d'une superposition d'états propres, voir paragraphe suivant.
  7. Par suite   est la probabilité que lors d'une mesure effectuée à   on trouve le système dans l'état  .

Voir aussi

modifier

Bibliographie

modifier

Liens externes

modifier