Voir la source de Oscillateur harmonique quantique

{{Homon|oscillateur}}
L''''oscillateur harmonique quantique''' correspond au traitement par les outils de la [[mécanique quantique]] de l'[[oscillateur harmonique|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 [[oscillateur harmonique|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 <math>V(x)=\dfrac{1}{2}kx^2</math>, ''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.

[[File:QuantumHarmonicOscillatorAnimation.gif|thumb|300px|right|Evolution temporelle de l'oscillateur harmonique d'après les [[lois du mouvement de Newton|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 [[État stationnaire (physique quantique)|é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 ==
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.
{{Article détaillé|oscillateur harmonique}}
=== Définition - Importance physique ===
Un [[oscillateur harmonique|oscillateur harmonique classique]] à une dimension est modélisé par un potentiel quadratique, typiquement (''m'' étant la masse du système) :
[[Fichier:Morse-potential-fr.png|vignette|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.]]
<math>V(x) = \frac{1}{2}m\omega^2 x^2</math>,
où <math>\omega</math> 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<ref>Plus précisément à la « pente » de la parabole correspondant à la courbe représentative de ''V(x)'', voir ci-après.</ref>.

L'énergie mécanique d'un oscillateur harmonique unidimensionnel est donc <math>H = \frac{p_{x}^2}{2m} + \frac{1}{2}m\omega^2 x^2</math>.

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 <math>V\left( x \right)</math>, avec un minimum en <math>x = x_0</math>, il vient au voisinage de celui-ci :
:<math>V\left(x\right) \approx V(x_0)+\frac{1}{2}\left(\frac{d^2 V}{dx^2}\right)_{x_0}\left(x-x_0\right)^2</math>, avec nécessairement pour un ''minimum'' <math>\left(\frac{d^2 V}{dx^2}\right)_{x_0} > 0</math>,
ce qui correspond bien à un oscillateur harmonique tel que <math>m\omega^2 \equiv \left(\frac{d^2 V}{dx^2}\right)_{x_0}</math>, l'expression précédente du potentiel correspondant au choix de <math>x_0</math> et <math>V\left(x_0\right)</math> 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 ===
Pour la discussion des états quasi-classiques, il est utile de rappeler des éléments sur la nature du mouvement classique.

Les [[Mécanique hamiltonienne|équations du mouvement de Hamilton]] pour le hamiltonien classique précédent s'écrivent sans difficultés sous la forme :
:<math>\begin{cases} \displaystyle{\frac{dx}{dt}} = \frac{p_x}{m}, \\ \\ \displaystyle{\frac{dp_x}{dt}} = -m\omega^2 x \end{cases}</math>.

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 :
:<math>P_x \equiv \frac{p_x}{\sqrt{m\omega}}</math> et <math>X \equiv \sqrt{m\omega}x</math>,
ce qui permet de mettre le hamiltonien sous la forme :
:<math>H = \frac{\omega}{2}\left(\frac{p_{x}^2}{m\omega}+m\omega x^2\right)=\frac{\omega}{2}\left(P_{x}^2+X^2\right)</math>.
Comme <math>\{ X,P_X \} =1</math> la transformation <math>\left(x,p_x\right)\rightarrow \left(X,P_X\right)</math> est [[Transformation canonique|canonique]], et donc laisse inchangée la forme des équations de Hamilton. Les équations du mouvement précédentes deviennent alors :
:<math>\begin{cases} \displaystyle{\frac{d X}{dt}}=\omega P_x \\ \\ \displaystyle{\frac{dP_x}{dt}} = - \omega X \end{cases}</math>.

L'introduction de la grandeur complexe <math>\alpha \equiv \frac{X+iP_x}{\sqrt{2}}</math>, permet de regrouper ces deux équations en une équation unidimensionnelle :
:<math>\frac{d\alpha}{dt}= - i\omega \alpha</math>,
dont la solution générale est de la forme <math>\alpha(t)=\alpha_0 e^{-i\omega t}</math>, avec <math>\alpha_0=\alpha(t=0)</math>.
Par suite une [[Particule (physique)|particule]] de masse <math>m</math> dans ce potentiel a un mouvement [[sinusoïdal]] de [[Vitesse angulaire|pulsation]] <math>\omega</math>. L'énergie de la particule est une constante du mouvement positive ou nulle <math>H=\frac{1}{2}m\left(\dot{x}^2+\omega^2 x^2\right)</math>, 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 <math>\alpha</math> sous la forme <math> H = \omega \alpha^* \alpha</math>.

== L'oscillateur harmonique quantique à une dimension ==
=== Position du problème ===

En mécanique quantique l'[[opérateur hamiltonien|hamiltonien]] d'un oscillateur harmonique unidimensionnel s'écrit sous la forme :

<math>\hat H = \frac{\hat p_{x}^2}{2m} + \frac{1}{2}m \omega^2 \hat x^2</math>

où <math>\hat p_{x}</math> est la composante sur l'axe <math>x</math> de l'opérateur [[impulsion (physique)|impulsion]] de la particule.

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

<math>\hat H \left | \psi \right \rangle = E\left | \psi \right \rangle</math>, où <math>E</math> est l'énergie associée à un [[fonction d'onde|état propre]] <math>\left | \psi \right \rangle</math> du système<ref>Conformément à une convention courante l'état propre de l'équation de Schrödinger dépendante du temps s'écrit avec Ψ majuscule : <math>\left|\Psi\right\rangle=\left|\Psi(x,t)\right\rangle</math> tandis que celui de l'équation indépendante du temps s'écrit avec ψ minuscule : <math>\left|\psi\right\rangle=\left|\psi(x)\right\rangle</math>. Il en est de même des fonctions d'onde associées dans les diverses représentations.</ref>.

Il est alors possible d'écrire cette équation en représentation position, compte tenu de <math>\hat x = x \hat 1</math> et <math>\hat p=\frac{\hbar}{i}\frac{d}{dx}</math>, ce qui donne l'équation différentielle suivante pour la [[fonction d'onde]] <math>\psi(x)=\left\langle x|\psi\right\rangle</math> :
:<math>\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2}+\left(E-\frac{m\omega^2}{2}x^2\right)\psi(x)=0</math>.

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 ===
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.

:* '''{{Souligner|Les valeurs propres ''E'' de <math>\hat H</math> sont toujours positives ou nulles}}''' : en effet pour un état <math>\left|\psi\right\rangle</math> quelconque, il est évident que <math>\|\hat x\left|\psi\right\rangle\|^2\geqslant 0</math> or <math>\|\hat x\left|\psi\right\rangle\|^2= \left\langle\psi\right|\hat x^{\dagger}\hat x\left|\psi\right\rangle=\left\langle\psi\right|\hat x^2\left|\psi\right\rangle=\left\langle\hat x^2\right\rangle \geqslant 0</math>.
:De façon évidente, il en est de même pour le terme en <math>\hat p_{x}^2</math>, donc pour tout état <math>\left|\psi\right\rangle</math>, <math>\left\langle\psi\right|\hat H\left|\psi\right\rangle \geqslant0</math>.

:* '''{{Souligner|L'énergie de l'état fondamental est nécessairement non nulle}}''' : en introduisant les variances des valeurs moyennes de <math>\hat x</math> et de <math>\hat p_x</math>, soit <math>\left(\Delta x\right)^2=\left\langle\hat x^2\right\rangle - \left\langle\hat x\right\rangle^2</math> et <math>\left(\Delta p_x\right)^2=\left\langle\hat p_{x}^2\right\rangle - \left\langle\hat p_{x}\right\rangle^2</math>, il est possible de montrer que <math>\left\langle\hat H\right\rangle\geqslant\frac{\hbar\omega}{2}</math>.
{{démonstration|
En choisissant un état « d'essai » <math>\left|\psi\right\rangle</math> tel que <math>\left\langle\hat x\right\rangle=0</math> et <math>\left\langle\hat p_{x}\right\rangle=0</math>, la valeur moyenne de <math>\hat H</math> s'exprime sous la forme : <br><math>\left\langle\hat H\right\rangle=\frac{1}{2m}\left\langle\hat p_{x}^2\right\rangle+\frac{m\omega^2}{2}\left\langle\hat x^2\right\rangle</math>.
:La [[Principe d'incertitude|relation d'incertitude]] <math>\Delta x\Delta p_x \geqslant \frac{\hbar}{2}</math> implique alors ''qu'au mieux'' <math>\left\langle\hat p_{x}^2\right\rangle = \frac{\hbar^2}{4\left\langle\hat x^2\right\rangle}</math>, la valeur moyenne du hamiltonien pour un état quelconque est donc ''a minima'' telle que : <math>\left\langle\hat H\right\rangle \geqslant \frac{\hbar^2}{8m\left\langle\hat x^2\right\rangle}+\frac{m\omega^2\left\langle\hat x^2\right\rangle}{2}</math>.
:Cette valeur peut être minimisée par rapport à <math>u\equiv\left\langle\hat x^2\right\rangle</math>, considéré comme variable<ref>Ce procédé, qui consiste à prendre une « fonction d'onde d'essai », et d'obtenir à partir de là une expression de <math>\left\langle\hat H\right\rangle</math>, dépendant de certains paramètres que l'on considère alors comme des variables, pour minimiser <math>\left\langle\hat H\right\rangle</math>, 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 <math>H_2^+</math> ou de façon numérique, par exemple dans la méthode de Monte-Carlo variationnelle).</ref>, ce qui implique <math>\frac{\partial\left\langle\hat H\right\rangle}{\partial u}=0=-\frac{\hbar^2}{8m\left\langle\hat x^2\right\rangle^2}+\frac{m\omega^2}{2}</math> soit finalement <math>\left\langle\hat x^2\right\rangle=\frac{\hbar}{2m\omega}</math> pour la valeur de <math>\left\langle\hat x^2\right\rangle</math> minimisant <math>\left\langle\hat H\right\rangle</math>.
: Il vient par substitution dans l'expression précédente : <math>\left\langle\hat H\right\rangle\geqslant\frac{\hbar\omega}{2}</math><ref>En fait, comme indiquée plus bas la valeur de l'état fondamental est bien <math>\hbar\omega/2</math>.</ref>}}
: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 <math>\hat x</math> et <math>\hat p_{x}</math>.
:* 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é (physique et chimie)|degré de liberté]], il n'y aura donc qu'un seul [[nombre quantique]].

=== Introduction d'opérateurs sans dimension ===

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é <math>\hbar\omega</math>, par suite (<math>\hbar</math> est la [[constante de Planck]] réduite) :
<math>\hat H=\frac{\hat p_{x}^2}{2m} + \frac{1}{2}m \omega^2 \hat x^2=\hbar\omega\left[\frac{\hat p_{x}^2}{2m\hbar\omega} + \frac{m\omega\hat x^2}{2\hbar}\right]= \frac{\hbar\omega}{2}\left[\left(\frac{\hat p_{x}}{\sqrt{m\hbar\omega}}\right)^2+\left(\sqrt{\frac{m\omega}{\hbar}}\hat x\right)^2\right]</math>

Il est facile de simplifier l'écriture de <math>\hat H</math>, en utilisant les opérateurs sans dimension (mais hermitiens) suivants :

<math>\hat P=\frac{\hat p_{x}}{\sqrt{m\hbar\omega}}</math>, et <math>\hat X=\sqrt{\frac{m\omega}{\hbar}} \hat x</math>.

Compte tenu de la relation de commutation <math>\left[\hat x,\hat p\right]=i\hbar</math>,la relation de commutation entre ces nouveaux opérateurs s'écrit :
:<math>\left[ \hat X,\hat P\right] = i\hat 1</math>.

En termes de ces opérateurs « réduits » le hamiltonien de l'oscillateur harmonique unidimmensionnel s'écrit sous la forme : <math>\hat H=\frac{\hbar\omega}{2}\left(\hat X^2+\hat P^2 \right)</math>.

{{Souligner|'''Remarques :'''}}
:* La quantité <math>\sqrt{\frac{\hbar}{m\omega}}</math> 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é <math>\hbar</math>.

=== Opérateurs d'échelle ===
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 :
:<math>\left(\hat X-i\hat P\right)\left(\hat X+i\hat P\right)=\hat X^2+\hat P^2+i\left(\hat X\hat P-\hat P\hat X\right) = \frac{2}{\hbar\omega} \hat H-\hat 1</math>.

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

<math>\begin{cases}
\hat a=\frac{\left( \hat X+i\hat P\right)}{\sqrt{2}}  \\ \\
\hat a^{\dagger}=\frac{ \left( \hat X-i\hat P\right)}{\sqrt{2}}
\end{cases}</math>

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 (opérateur)|commutateur]] entre ces deux opérateurs : <math>\left[ \hat a,\hat a^{\dagger} \right]=\hat 1\qquad\text{, }</math>.

De même, les opérateurs <math>\hat X</math> et <math>\hat P</math> s'expriment sous la forme :
<math>\begin{cases}
\hat X=\frac{1}{\sqrt{2}}\left(\hat a^\dagger+\hat a \right)  \\
\hat P=i\frac{1}{\sqrt{2}}\left(\hat a^{\dagger}-\hat a\right)
\end{cases}</math>

Le hamiltonien du système peut alors s'écrire sous la forme :
:<math>\hat H=\hbar\omega\left(\hat a^{\dagger}\hat a+\frac{1}{2}\hat 1\right)=\hbar\omega\left(\hat N+\frac{1}{2}\hat 1\right)</math>,
où il a été introduit l'opérateur :
:<math>\hat N \equiv\hat a^{\dagger}\hat a</math>, 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 <math>\hat N</math>, qui est évidemment hermitien puisque <math>\hat N^{\dagger}=\left(\hat a^{\dagger}\hat a\right)^{\dagger}=\hat a^{\dagger}\left(\hat a^{\dagger}\right)^{\dagger}=\hat N</math>. Par suite, ses valeurs propres sont ''réelles''.

{{Souligner|'''Remarques :'''}}
:* l'opérateur d'échelle <math>\hat a</math> peut être considéré comme le « pendant quantique » de la quantité complexe <math>\alpha</math> 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 <math>\hbar</math>, 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 <math>\hat x</math> et <math>\hat p_x</math>), qui donnera lieu à l'apparition de l'énergie dite de « point zéro ».

=== Calcul des valeurs propres ===

Les états propres commun de <math>\hat N</math> et <math>\hat H</math> sont notés <math>\left|\nu\right\rangle</math>, avec pour valeurs propres : <math>\hat N\left|\nu\right\rangle=\nu\left|\nu\right\rangle</math>, la valeur propre ν étant ''a priori'' un réel.
Les [[Valeur propre (synthèse)|valeurs propres]] de <math>\hat H</math> sont alors de la forme : <math> E_{\nu}=\hbar\omega\left(\nu+\frac{1}{2}\right)</math>. 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 ====
À partir de la relation de commutation précédente entre les opérateurs d'échelle, il vient facilement :
* <math>\left[\hat N,\hat a\right]=\hat N\hat a-\hat a\hat N=\hat a^{\dagger}\hat a\hat a-\hat a\hat a^{\dagger}\hat a=-\left[\hat a,\hat a^{\dagger}\right]\hat a=-\hat a</math>
*<math>\left[\hat N,\hat a^{\dagger}\right]=\hat N\hat a^{\dagger}-\hat a^{\dagger}\hat N=\hat a^{\dagger}\hat a\hat a^{\dagger}-\hat a^{\dagger}\hat a^{\dagger}\hat a=\hat a^{\dagger}\left[\hat a,\hat a^{\dagger}\right]=\hat a^{\dagger}</math>

Ces diverses relations rendent possible la détermination des valeurs propres de <math>\hat N</math>.

==== Action des opérateurs d'échelle sur les états propres <math>\left|n\right\rangle</math> ====
En utilisant les notations et relations précédentes, il vient facilement :
:<math>\left[\hat N,\hat a\right]\left|\nu\right\rangle=\hat N\left(\hat a\left|\nu\right\rangle\right)-\nu\left(\hat a\left|\nu\right\rangle\right)=-\hat a\left|\nu\right\rangle</math>, ce qui donne :
:<math>\hat N\left(\hat a\left|\nu\right\rangle\right)=(\nu-1)\left(\hat a\left|\nu\right\rangle\right)</math>,
:autrement dit <math>\left(\hat a\left|\nu\right\rangle\right)</math> est un ''vecteur propre'' de <math>\hat N</math> de valeur propre ''(ν-1)'', il est possible de poser :
:<math>\hat a\left|\nu\right\rangle=\alpha\left|\nu-1\right\rangle</math>, avec α [[nombre complexe]] qui peut être déterminé à une phase près en utilisant le fait que :
:<math>\left|\alpha\right|^2=\|\hat a\left|\nu\right\rangle \|^2=\left\langle \nu\right|\hat a^{\dagger}\hat a\left|\nu\right\rangle=\nu</math>,
et comme nécessairement <math>\|\hat a\left|\nu\right\rangle \|^2\geqslant 0</math> il est possible d'en déduire successivement :
# que les valeurs propres ''ν'' de <math>\hat N</math> sont {{Souligner|'''positives ou nulles'''}};
# que si en particulier ''ν = 0'', alors <math>\hat a\left|0\right\rangle = \vec{0}</math> et inversement si <math>\hat a\left|\nu\right\rangle=\vec{0}</math> alors <math>\hat a^{\dagger}\hat a\left|\nu\right\rangle=\vec{0}=\hat N\left|\nu\right\rangle</math> et donc tout état respectant cette condition est état propre de <math>\hat N</math> de valeur propre ''ν = 0'';
# qu'en effectuant un choix de phase tel que α soit réel, il vient : <math>\alpha=\sqrt{\nu}</math>, soit encore : <math> \hat a\left|\nu\right\rangle =\sqrt{\nu}\left|\nu-1\right\rangle</math>,
par suite l'opérateur <math>\hat a</math> ''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 <math>\hat a^{\dagger}</math> il est possible de vérifier avec la même convention de phase que <math> \hat a^{\dagger}\left|\nu\right\rangle=\sqrt{\nu+1}\left|\nu+1\right\rangle</math>,
par suite l'opérateur <math>\hat a^{\dagger}</math> ''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'' ====

Soit un état propre <math>\left|\nu\right\rangle</math> avec ''ν non entier'', positif, ''N'' étant le premier entier tel que ''N > ν'', il résulte de la relation <math> \hat a\left|\nu\right\rangle =\sqrt{\nu}\left|\nu-1\right\rangle</math> que l'application à ''N'' reprises de l'opérateur d'annihilation <math>\hat a</math> permettait d'obtenir un état propre de <math>\hat N</math> avec une valeur propre ''négative'', ce qui n'est pas possible puisque les valeurs propres de <math>\hat N</math> sont positives ou nulles.

En revanche si ''ν = n'', ''n'' [[Entier naturel|entier strictement positif]], ''n'' applications successives de <math>\hat a</math> donnent : <math>\left(\hat a\right)^n\left|\nu=n\right\rangle=\sqrt{n!}\left|0\right\rangle</math>, et une nouvelle application de <math>\hat a</math> donne alors le [[vecteur nul]].

Par suite {{Souligner|'''''les seules valeurs propres ν de l'opérateur <math>\hat N</math> 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 :

<center>
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-----
| bgcolor="#fff8ff" | Les énergies accessibles par l'oscillateur sont : <math>E_n=\hbar\omega\left(n+\frac{1}{2}\right)</math> avec ''n'' {{Souligner|'''entier positif ou nul'''}}
|}</center>

Ces états ont les propriétés générales suivantes :
:# '''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.
:# '''Les états d'énergie sont espacés de la même quantité <math>\hbar\omega</math>'''. 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 ''[[phonon]]s'', 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é « [[photon]]s ».
:# Comme cela a déjà mis en évidence qualitativement, {{Souligner|'''l'énergie de l'état fondamental est non nulle'''}}, avec <math>E_0=\frac{\hbar\omega}{2}</math>, 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 [[Principe d'incertitude|relation d'incertitude de Heisenberg]].

=== Détermination des états propres ===
==== Non-dégénérescence de l'état fondamental ====
En supposant que l'état fondamental est ''k-''fois dégénéré, il est possible de poser <math>\left|0\right\rangle=\left|0,\lambda\right\rangle</math> avec λ = 1,2,..., k. Comme pour tout λ <math>\hat a\left|0,\lambda\right\rangle=\vec{0}</math>, il est possible de revenir à la définition l'opérateur d'annihilation <math>\hat a</math> et de donner l'équation différentielle correspondante en représentation position que doivent satisfaire les fonctions d'onde <math>\psi_0^{\lambda}(x)</math> pour tout ''k'' :
:<math>\frac{d\psi_0^{\lambda}}{dx}+\frac{m\omega}{\hbar}x\psi_0^{\lambda}(x)=0</math>.
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 : <math>\psi_0(x)=C_0\exp{\left(-\frac{m\omega}{2\hbar}x^2\right)}</math>.

La normalisation de ce résultat, permet d'obtenir <math>C_{0} = \left(\frac{m \omega}{\pi \hbar}\right)^{\frac{1}{4}}</math>.

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 <math>\left|n,1\right\rangle</math> et <math>\left|n,2\right\rangle</math> de même énergie <math>E_n</math> d'obtenir, par application à ''n'' reprises de l'opérateur d'échelle <math>\hat a</math> deux états fondamentaux dégénérés <math>\left|0,1\right\rangle</math> et <math>\left|0,2\right\rangle</math> distincts, ce qui est impossible du fait de la non-dégénérescence de l'état fondamental<ref>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 <math>E_{n_1,\ldots,n_p} = \left(n_1+\ldots+n_p+p/2\right)\hbar\omega</math>. Il est clair qu'il y a plusieurs combinaisons de valeurs possibles des <math>n_i</math>, donc plusieurs états distincts, donnant une même valeur de l'énergie du système.</ref>.

==== États propres de l'opérateur ''N'' ====

Une [[Raisonnement par récurrence|démonstration par récurrence]] utilisant l'opérateur de création de quantum de vibration montre que les états propres de <math>\hat a^\dagger \hat a</math> s'écrivent : <math>|n\rangle=\frac{1}{\sqrt{n!}}\left(\hat a^\dagger \right)^n |0\rangle</math>.

[[Fichier:HarmOsziFunktionen.png|vignette|Représentation des sept premiers niveaux d'énergie et des fonctions d'onde associées <math>\psi_n(x)</math>) de l'oscillateur harmonique quantique unidimensionnel.]]

En représentation position, il suffit de substituer l'expression de <math>\hat a</math> et de <math>\psi_0(x)=\left\langle x|0\right\rangle</math> pour obtenir l'expression de la fonction d'onde <math>\psi_n(x)=\left\langle x|n\right\rangle</math> sous la forme :
:<math>\psi_n(x) = \frac{1}{\sqrt{2^{n} n!}}\left(\frac{m\omega}{\pi \hbar}\right)^{\dfrac{1}{4}} \left( X - \frac{\mathrm{d}}{\mathrm{d}X} \right) ^{n}  \exp{\left(-\frac{X^2}{2}\right)}</math>, avec <math>X=\sqrt{\frac{m\omega}{\hbar}}x</math> position réduite.
[[Fichier:Aufenthaltswahrscheinlichkeit harmonischer Oszillator.png|vignette|Représentation similaire à la précédente, avec cette fois-ci les allures des densités de probabilité de présence <math>| \psi_n(x) |^2</math> pour les premiers états d'énergie de l'oscillateur harmonique unidimensionnel.]]
En introduisant les [[polynôme d'Hermite|polynômes d'Hermite]] <math>H_n(X)</math> sous leur forme dite « physique » définie par <math>e^{-\frac{X^2}{2}}H_n(X)=\left( X - \frac{\mathrm{d}}{\mathrm{d}X} \right) ^{n}  e^{-\frac{X^2}{2}}</math> il est finalement possible d'obtenir l'expression générale de la fonction d'onde <math>\psi_n(x)</math> :
:<math>\psi_n(x)=\frac{1}{\sqrt{2^{n} n!}}\left(\frac{m\omega}{\pi \hbar}\right)^{\dfrac{1}{4}}H_n\left(\sqrt{\frac{m\omega}{\hbar}}x\right)e^{-\frac{m\omega}{2 \hbar}x^2}</math>.

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 <math>\psi_n(x)</math> est celle de ''n''.

Cette dernière relation nous permet de retrouver explicitement autant de [[Fonction d'onde|fonctions d'ondes]] que nécessaire. Par exemple pour <math>n=1\,</math> il vient : <br><math>\psi_{1}(x) = \left( \frac{4}{\pi} \left( \frac{m \omega}{\hbar} \right )^{3} \right )^{\frac{1}{4}} x \; e^{-\frac{1}{2} \frac{m \omega}{\hbar} x^{2} }</math>

=== Représentation matricielle ===

* La matrice représentative de l'hamiltonien <math>\hat H</math> sur la base des <math>\vert n \rangle</math> est, par construction, diagonale.

On a <math>\hat H= \hbar \omega \begin{pmatrix} 1/2 & 0 & 0 &\cdots \\ 0 & 1+1/2 & 0 &\cdots \\ 0 & 0 & 2+1/2 &\cdots \\ \cdots & \cdots & \cdots &\ddots \end{pmatrix}</math>
* Sachant que <math>\hat a^{\dagger} \vert n \rangle = \sqrt{n+1} \vert n+1 \rangle</math>, en multipliant à gauche par <math>\vert k \rangle</math> il vient <math>\hat a^{\dagger}_{k,n}=\sqrt{n+1} \langle k \vert n+1 \rangle = \sqrt{n+1} \delta_{k,n+1}</math>

La matrice représentative de <math>\hat a^{\dagger}</math> sur la base des <math>\vert n \rangle</math> est donc
<math>\hat a^{\dagger}= \begin{pmatrix} 0 & 0 & 0 & 0 & \cdots \\ 1 & 0 & 0 & 0 & \cdots \\ 0 & \sqrt{2} & 0 & 0 &\cdots \\ 0 & 0 & \sqrt{3} & 0 & \cdots \\ \cdots & \cdots & \cdots & \cdots &\ddots \end{pmatrix}</math>

Puisque <math>\hat a_{k,n}=\hat a^{\dagger}_{n,k}</math>, la matrice <math>\hat a</math> est construite en transposant <math>\hat a^{\dagger}</math> : <math>\hat a= \begin{pmatrix} 0 & 1 & 0 & 0 & \cdots \\ 0 & 0 & \sqrt{2} & 0 & \cdots \\ 0 & 0 & 0 & \sqrt{3} &\cdots \\ 0 & 0 & 0 & 0 & \cdots \\ \cdots & \cdots & \cdots & \cdots &\ddots \end{pmatrix}
</math>
* Il est alors facile de construire les matrices représentatives des observables <math>\hat x</math> et <math>\hat p</math> puisque :

<math> \hat x=\sqrt{\dfrac{\hbar}{m\omega}} \hat X= \sqrt{\dfrac{\hbar}{2m\omega}} \left(\hat a + \hat a^{\dagger}\right)</math>et  <math> \hat p=\sqrt{\hbar m\omega}\hat P =-i\sqrt{\dfrac{\hbar m\omega}{2}}\left(\hat a - \hat a^{\dagger}\right)</math>

donc <math>\hat x= \sqrt{\dfrac{\hbar}{2m\omega}} \begin{pmatrix} 0 & 1 & 0 & 0 & \cdots \\ 1 & 0 & \sqrt{2} & 0 & \cdots \\ 0 & \sqrt{2} & 0 & \sqrt{3} &\cdots \\ 0 & 0 & \sqrt{3} & 0 & \cdots \\ \cdots & \cdots & \cdots & \cdots &\ddots \end{pmatrix}</math>
et
<math>\hat p= \sqrt{\dfrac{\hbar m\omega}{2}} \begin{pmatrix} 0 & -i & 0 & 0 & \cdots \\ i & 0 & -i\sqrt{2} & 0 & \cdots \\ 0 & i\sqrt{2} & 0 & -i\sqrt{3} &\cdots \\ 0 & 0 & i\sqrt{3} & 0 & \cdots \\ \cdots & \cdots & \cdots & \cdots &\ddots \end{pmatrix}</math>

Il est alors possible de vérifier les propriétés suivantes pour les observables <math>\hat x</math> et <math>\hat p_x</math> :

:* Les seuls éléments de matrices non nuls de ces opérateurs sont donc ceux pris entre les états <math>\left|n\right\rangle</math> et <math>\left|n\pm 1\right\rangle</math>.
:* Les valeurs moyennes de la position et de l'impulsion sont nulles lorsque le système est dans un état propre <math>\left|n\right\rangle</math><ref>Ce qui ne signifie pas que sur un état donné du système de vecteur d'état <math>\left|\Psi(t)\right\rangle</math> ces valeurs moyennes soient nulles, puisque dans le cas général <math>\left|\Psi(t)\right\rangle</math> résulte d'une superposition d'états propres, voir paragraphe suivant.</ref>.

=== Évolution des valeurs moyennes dans le temps - limite classique ===
Dans le cas d'un hamiltonien <math>\hat H</math> qui ne dépend pas explicitement du temps, il résulte de l'équation de Schrödinger <math>{\partial\over \partial t}\left|\Psi(t)\right\rangle=\dfrac{-i}{\hbar}\hat H\left|\Psi(t)\right\rangle</math> que le vecteur d'état <math>\left|\Psi(t)\right\rangle</math> s'obtient à partir de celui pris à un temps <math>t_0</math> quelconque par :
:<math>\left|\Psi(t)\right\rangle=\exp{\left(-\dfrac{i(t-t_0)}{\hbar}\hat H\right)}\left|\Psi(t_0)\right\rangle</math>.

Or comme <math>\left|\Psi(t_0)\right\rangle=\sum_{n=0}^{+\infty} c_n(t_0)\left|n\right\rangle</math>, les coefficients complexes <math>c_n(t_0)</math> étant tels que <math>\left|\Psi(t_0)\right\rangle</math> est normalisé à l'unité<ref>Par suite <math>\left|c_n(t_0)\right|^2</math> est la probabilité que lors d'une mesure effectuée à <math>t=t_0</math> on trouve le système dans l'état <math>\left|n\right\rangle</math>.</ref>, dans le cas d'un oscillateur harmonique l'expression du vecteur d'état à une date ''t''(''>t0'') quelconque est la suivante :
:<math>\left|\Psi(t)\right\rangle=\sum_{n=0}^{+\infty}c_n(t_0)\exp{\left(-\dfrac{i(t-t_0)}{\hbar}\hat H\right)}\left|n)\right\rangle=\exp{\left(\frac{-i\omega(t-t_0)}{2}\right)}\sum_{n=0}^{+\infty} c_n(t_0)\exp{\left(-in\omega(t-t_0)\right)}\left|n\right\rangle</math>.

Par suite, la valeur moyenne de tout observable <math>\hat a</math> est donnée en fonction du temps par :
:<math>\left\langle\hat a\right\rangle =\left\langle\Psi(t)\right|\hat a\left|\Psi(t)\right\rangle= \sum_{p=0}^{+\infty} \sum_{n=0}^{+\infty} c_p^*(t_0)c_n(t_0)A_{mn}\exp{\left(i(m-n)\omega (t-t_0)\right)}</math>, avec <math>A_{mn}\equiv\left\langle m \right|\hat a\left|n\right\rangle</math>, élément de matrice de l'observable <math>\hat a</math> entre <math>\left|m\right\rangle</math> et <math>\left|n\right\rangle</math>.

En particulier, pour les observables <math>\hat x</math> et <math>\hat p_x</math> les seuls éléments de matrices non nuls sont ceux entre <math>\left|n\pm 1\right\rangle</math> et <math>\left|n\right\rangle</math>, par suite, l'évolution de leurs grandeurs moyennes fait uniquement intervenir des termes en <math>\exp{\left(\pm i \omega(t-t_0)\right)}</math>, 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 <math>\hat x</math> et <math>\hat p_x</math> ne dépendent pas explicitement du temps s'écrit pour chaque valeur moyenne :
:<math>\frac{d}{dt}\left\langle\hat x\right\rangle = \frac{1}{i\hbar}\left\langle\left[\hat x,\hat H\right]\right\rangle=\frac{1}{2m i\hbar}\left\langle\left[\hat x,\hat p_{x}^2\right]\right\rangle=\frac{1}{2m i\hbar}\left\langle\left(\left[\hat x,\hat p_x\right]\hat p_x + \hat p_x\left[\hat x,\hat p_x\right]\right)\right\rangle</math>,

:<math>\frac{d}{dt}\left\langle\hat p_x\right\rangle = \frac{1}{i\hbar}\left\langle\left[\hat p_x,\hat H\right]\right\rangle=\frac{m\omega^2}{2 i\hbar}\left\langle\left[\hat p_x,\hat x^2\right]\right\rangle=\frac{m\omega^2}{2 i\hbar}\left\langle\left(\left[\hat p_x,\hat x\right]\hat x + \hat x\left[\hat p_x,\hat x\right]\right)\right\rangle</math>,

soit finalement le système d'équations :
:<math>\begin{cases}
\frac{d}{dt}\left\langle\hat x\right\rangle =\frac{\left\langle\hat p_x\right\rangle}{m} \\
\frac{d}{dt}\left\langle\hat p_x\right\rangle =-m\omega^2\left\langle\hat x\right\rangle
\end{cases} </math>.

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 <math>\left\langle\hat x\right\rangle</math> et <math>\left\langle\hat p_x\right\rangle</math>.

=== Interprétation ===

[[Fichier:QHarmonicOscillator.png|vignette|Densité de probabilité<br>
<nowiki>|</nowiki>ψ<sub>n</sub>(x)<nowiki>|</nowiki>{{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 [[Fonction d'onde|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 ===
{{Article détaillé|État cohérent}}
Les ''états cohérents'', notés <math>\left|\alpha\right\rangle</math>, sont par définition les états propres de l'opérateur d'annihilation <math>\hat a</math>, de valeurs propres <math>\alpha\in\Complex</math> : <math>\hat a\left|\alpha\right\rangle=\alpha\left|\alpha\right\rangle</math>.

==== Propriétés générales des états cohérents <math>\left|\alpha\right\rangle</math> ====
===== Expression de <math>\left|\alpha\right\rangle</math> sur la base des états propres <math>\left|n\right\rangle</math> du hamiltonien =====
Il est possible de poser de façon générale pour tout état <math>\left|\alpha\right\rangle</math> :
:<math>\left|\alpha\right\rangle=\sum_{n=0}^{\infty} c_{n}(\alpha)\left|n\right\rangle</math>, avec <math>c_{n}(\alpha)\in\Complex</math>, qui seront choisis de telle sorte que <math>\left|\alpha\right\rangle</math> soit normé à l'unité : <math>\left\langle\alpha|\alpha\right\rangle = 1 \quad \Rightarrow \quad \sum_{n=0}^{\infty} \left|c_{n}(\alpha)\right|^2=1</math>.

D'après la définition d'un état cohérent donnée plus haut, et le fait que pour tout état propre <math>\left|n\right\rangle</math>, <math>\hat a\left|n\right\rangle=\sqrt{n}\left|n-1\right\rangle</math> il vient :
:<math>\hat a\left|\alpha\right\rangle=\sum_{n=1}^{\infty} c_{n}(\alpha)\sqrt{n}\left|n-1\right\rangle=\alpha\left|\alpha\right\rangle=\sum_{n=0}^{\infty}\alpha c_{n}(\alpha)\left|n\right\rangle</math>, où l'on a tenu compte de <math>\hat a\left|0\right\rangle = 0</math>.
En procédant à un changement de variable dans le premier membre, l'équation aux valeurs propres s'écrit dans la base de <math>\left|n\right\rangle</math> :
:<math>\sum_{n=0}^{\infty} \sqrt{n+1}c_{n+1}(\alpha)\left|n\right\rangle=\sum_{n=0}^{\infty} \alpha c_{n}(\alpha)\left|n\right\rangle</math>,
ce qui donne par identification membre à membre la [[Suite définie par récurrence|relation de récurrence]] suivante entre les coefficients <math>c_{n}(\alpha)</math> :
:<math>c_{n+1}(\alpha)=\frac{\alpha}{\sqrt{n+1}}c_{n}(\alpha)</math>, qui permet alors d'obtenir l'expression de <math>c_{n}(\alpha)</math> en fonction de <math>c_{0}(\alpha)</math> :
:<math>c_{n}(\alpha)=\frac{\alpha^n}{\sqrt{n!}}c_{0}(\alpha)</math>.
La constante <math>c_{0}(\alpha)</math> peut être obtenue en utilisant la condition de normalisation précédente, qui donne alors <math>\sum_{n=0}^{\infty} \frac{\left|\alpha\right|^{2n}}{n!}\left|c_{0}(\alpha)\right|^2 = \left|c_{0}(\alpha)\right|^2\exp{\left(\left|\alpha\right|^2\right)}=1</math>, soit en prenant pour convention de phase que tous les coefficients <math>c_{n}(\alpha)</math> sont réels, il vient finalement : <math>c_{n}(\alpha)=\frac{\alpha^n}{\sqrt{n!}}\exp{\left(-\tfrac{\left|\alpha\right|^2}{2}\right)}</math>.

Le développement de <math>\left|\alpha\right\rangle</math> sur la base des états propres <math>\left|n\right\rangle</math> s'écrit en définitive sous la forme :
<math>\left|\alpha\right\rangle=e^{\left(-\tfrac{\left|\alpha\right|^2}{2}\right)} \sum_{n=0}^{\infty} \frac{\alpha^n}{\sqrt{n!}}\left|n\right\rangle</math>.

Comme <math>|n\rangle=\frac{1}{\sqrt{n!}}\left(\hat a^\dagger \right)^n |0\rangle</math>, il est aussi possible d'écrire :
:<math>\left|\alpha\right\rangle=e^{\left(-\tfrac{\left|\alpha\right|^2}{2}\right)}\sum_{n=0}^{\infty} \frac{\alpha^n}{n!}\left(\hat a^{\dagger}\right)^n\left|0\right\rangle=e^{\left(-\tfrac{\left|\alpha\right|^2}{2}\right)} \exp{\left(\alpha\hat a^{\dagger}\right)}\left|0\right\rangle</math>.

Il est facile de vérifier que si ''α = 0'' ce développement se réduit à <math>\left|\alpha=0\right\rangle=\left|0\right\rangle</math> : ''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 =====
Il est facile de démontrer que les états cohérents <math>|\alpha\rangle</math> 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 <math>\Delta x\Delta p= \frac{\hbar}{2}</math>.

{{démonstration|La démonstration est aisée compte tenu des relations
:<math> \hat x=\sqrt{\dfrac{\hbar}{m\omega}} \hat x= \sqrt{\dfrac{\hbar}{2m\omega}} \left(\hat a + \hat a^{\dagger}\right)</math>,
:<math> \hat p=\sqrt{\hbar m\omega}\hat p=-i\sqrt{\dfrac{\hbar m\omega}{2}}\left(\hat a - \hat a^{\dagger}\right)</math>,
ainsi que du fait que :
:<math>\langle\alpha|\hat a^{\dagger}=\left(\hat a|\alpha\rangle\right)^{\dagger}</math>,
et <math>\hat a\hat a^{\dagger}=\hat{1}+\hat a^{\dagger}\hat a</math>.

Il est facile de vérifier que l'on obtient pour les valeurs moyennes de <math>\hat x</math>, <math>\hat x^2</math>, <math>\hat p</math> et <math>\hat p^2</math> sur un état cohérent quelconques les résultats suivants :
:<math>\langle\hat x\rangle=\sqrt{\frac{\hbar}{2m\omega}}\left(\alpha+\alpha^*\right)</math> et <math>\langle\hat x^2\rangle=\frac{\hbar}{2m\omega}\left(\alpha^2+\alpha^{*2}+2|\alpha|^2+1\right)</math>,
:<math>\langle\hat p\rangle=-i\sqrt{\frac{m\hbar\omega}{2}}\left(\alpha-\alpha^*\right)</math> et <math>\langle\hat p^2\rangle=\frac{m\hbar\omega}{2}\left(1+2|\alpha|^2-\alpha^2-\alpha^{*2}\right)</math>.

Dès lors il vient les écarts-types de <math>\hat x</math> et <math>\hat p</math> :
:<math>\Delta_x=\sqrt{\langle\hat x^2\rangle-\langle\hat x\rangle^2}=\sqrt{\frac{\hbar}{2m\omega}}</math>,
et <math>\Delta_p=\sqrt{\langle\hat p^2\rangle-\langle\hat p\rangle^2}=\sqrt{\frac{m\hbar\omega}{2}}</math>.

Ces valeurs ne dépendent pas de ''α'' et donnent bien l'égalité <math>\Delta x\Delta p= \frac{\hbar}{2}</math>.}}

=== Applications ===

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 : <math> \hat H = \left( \frac{\hat p_{x}^2+\hat p_{y}^2}{2m} + {1\over 2} m \omega^2 \left(\hat x^2+\hat y^2 \right) \right)</math>

== Généralisation à trois dimensions ==

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.

<center>
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-----
| bgcolor="#fff8ff" | Les énergies accessibles par l'oscillateur sont : <math>E_n=\hbar\omega\left(n+\frac{3}{2}\right)</math> <math>(n\ge 0)</math>
|}</center>

avec <math>n=n_{x}+n_{y}+n_{z}</math> : 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]] : <math> g_{n} = \frac{(n+1)(n+2)}{2} </math>

== Oscillateur anharmonique ==

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 <math>W</math> dans l'hamiltonien puis, à l'aide de la [[théorie des perturbations stationnaires]], on va calculer les nouvelles énergies du système.

=== Terme en x<sup>3</sup> ===

Si l'on pousse le développement limité du potentiel à l'ordre 3, on obtient le terme perturbatif suivant par rapport à l'oscillateur harmonique : <math>W = \epsilon \hbar \omega X^{3}</math> où <math>\epsilon</math> 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 : <math>E_n = \left( n + {1 \over 2} \right) \hbar \omega - {15 \over 4} \epsilon^2 \left( n + {1 \over 2} \right)^2 \hbar \omega - {7 \over 16} \epsilon^2 \hbar \omega + ...</math>

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 : <math>\Phi_n (x) = \Psi_n (x) + \epsilon \left( - \alpha (n) \Psi_{n+1} (x)+\alpha (n-1) \Psi_{n-1} (x) - \beta(n) \Psi_{n+3} (x) + \beta (n-3) \Psi_{n-3} (x) \right)  + ... </math>

=== Application ===

Une [[molécule diatomique]] dans le modèle harmonique peut absorber et émettre des [[onde électromagnétique|ondes électromagnétiques]] à la pulsation <math>\omega</math>. 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 <math>\Phi_{n_{1}}</math> et <math>\Phi_{n_{2}}</math> peuvent être reliés même si <math>n_1</math> et <math>n_2</math> différent de plus de 1. La constante <math>\epsilon</math> peut ainsi être calculée dans cet exemple en analysant le [[spectre électromagnétique|spectre]] d'absorption.

== Notes ==
{{Références}}

== Voir aussi ==
* [[oscillateur]]
* [[oscillateur harmonique|oscillateur harmonique classique]]
* [[oscillateur de Dunkl]]
* [[phonons]]
* [[propagateur de l'oscillateur harmonique]]

== Bibliographie ==
* {{Cohen}}.
* {{Messiah}}.
* {{Landau|tome 3}}.
* {{Landau|tome 4}}.
* {{Basdevant}}.
* R. Feynman, ''Mécanique Quantique''.
* Y. Ayant, E. Belorizky, ''Cours de mécanique quantique''.
* R. Shankar, ''Principles of quantum mechanics'', {{2e|édition}}, NY, Plenum publishing, 1994, {{chap.|7}}.
*  O. Henri-Rousseau and Paul Blaise,"Quantum Oscillators", Wiley,Hoboken, NJ, 2011, 647 pages
== Liens externes ==
* [http://www.u-bourgogne.fr/PHYSIQUE/teroscil/Applet.html Tracé des fonctions d'onde et de la densité de probabilité associée de l'oscillateur harmonique quantique].

{{Portail|physique}}

[[Catégorie:Modèle quantique]]