Approximation BKW

En physique, l'approximation BKW Brillouin - Kramers - Wentzel est une méthode développée en 1926 qui permet d'étudier le régime semi-classique d'un dispositif quantique.

Catégories :

Physique quantique - Physique mathématique

Source image : fr.wikipedia.org
Cette image est un résultat de recherche de Google Image. Elle est peut-être réduite par rapport à l'originale et/ou protégée par des droits d'auteur.

Page(s) en rapport avec ce sujet :

... et avec l'approximation de la phase stationnaire en utilisant la formule BKW... Elément matriciel ; Théorie quantique ; Fonction onde ; Espace phase... (source : cat.inist)

En physique, l'approximation BKW Brillouin^[1] - Kramers ^[2] - Wentzel ^[3] est une méthode développée en 1926 qui permet d'étudier le régime semi-classique d'un dispositif quantique. La fonction d'onde est développée asymptotiquement au premier ordre de la puissance du quantum d'action $\hbar$ .

L'idée de base de la méthode BKW est que l'équation de Schrödinger se dérive de l'équation de propagation des ondes. On doit par conséquent retrouver la mécanique classique dans la limite $\hbar \rightarrow 0$ comme on retrouve l'optique géométrique quand la longueur d'onde $\lambda \rightarrow 0$ dans la théorie de l'optique ondulatoire.

L'approximation BKW (pour les francophones) est aussi connue sous les initiales WKB (pour les anglophones), WKBJ, BWKJ et quelquefois WBK, BWK. Le J supplémentaire est pour le mathématicien Harold Jeffreys, qui a développé en 1923 une méthode générale d'approximation pour des équations différentielles linéaires du second ordre, qui inclut l'équation de Schrödinger à une dimension. Les trois physiciens BKW n'avaient apparemment pas eu connaissance de ce travail.

Formule à une dimension d'espace

De façon générale la fonction d'onde est mise sous la forme ansatz :

$\psi(\vec{r},t) \ = \ A(\vec{r},t) \ \exp \left[ \ \frac{i}{\hbar} \ S(\vec{r},t) \ \right]$

Les deux fonctions inconnues sont l'amplitude A et l'action S, l'une de ces deux fonctions est généralement reconnue comme «lentement variable». En fait seul le cas unidimentionnel où $\vec{r}=R$ est utilisé, c'est ce cas que nous allons développer ici.

Formule BKW

Notons $ψ$ la fonction d'onde, solution stationnaire de l'équation de Schrödinger, d'une particule de masse $m$ se déplaçant dans le potentiel $V (R)$ ;

$\left[ -{ \hbarˆ2 \over 2 m } {dˆ2 \over d Rˆ2} + V(R) \right] \psi(R) = E \psi(R)$

L'approximation BKW consiste à écrire la fonction d'onde sous la forme

$\psiˆ{\rm BKW} (R) = {C_+ \over \sqrt{ |p(R) | } } eˆ{ {i \over \hbar} \int p } + {C_- \over \sqrt{ |p(R) | } } eˆ{- {i \over \hbar} \int p }$

où $p(R) = \sqrt{ 2 m (E - V(R))}$ est l'impulsion locale de la particule.

Sens physique

Notons le sens physique simple :

Dans la région classiquement permise plus la particule va vite, plus sa probabilité de présence diminue. En effet à l'endroit où $E > V (R)$ , la probabilité de présence $| ψ | 2$ sera proportionnelle à $1\over p$ .
Dans la région classiquement interdite la probabilité de présence $| ψ | 2$ sera exponentiellement décroissante en $eˆ{- \frac{1}{\hbar} \int |p | }$ . En effet à l'endroit où $E < V (R)$ , on a alors $p(R) = i \sqrt{ | 2 m (E - V(R)) | }$ et le terme exponentiellement croissant sera généralement divergent et par conséquent non physique, la normalisation de la fonction d'onde impose alors $C + = 0$ )

Domaine de validité

Le domaine de validité de l'approximation est le suivant

$\left| {d \bar \lambda \over d R } \right| \ll 1$ où $\bar \lambda(R) = \hbar / p(R)$ est la longueur d'onde de de Broglie dite réduite (divisée par 2 π)
Une seconde condition, qui est fréquemment vérifiée, vient s'ajouter à celle-ci mais elle est rarement utilisée $\left| \int { 1 \over \bar \lambda (R) } \left( {d \bar \lambda \over d R } \right)ˆ2 \right|ˆ2 \ll 1$

La première condition peut s'interpréter, en faisant apparaitre le produit $V'(R) \bar \lambda (R)$ , comme une condition adiabatique, i. e. comme le fait que le potentiel $V$ doit changer lentement sur des distances identiques à longueur d'onde $\bar \lambda$ de la particule pour que celle-ci ait le temps de s'adapter au nouveau potentiel lors du mouvement.

La seconde condition est plus complexe à interpréter mais elle indique qu'il faut être prudent si le potentiel V décroit trop lentement à l'infini.

Démonstration

En faisant apparaître les différents ordres du développement en puissance de $\hbar$ on pose

$\psi(R) = eˆ{\frac{i}{ \hbar} \left[ \sigma_0 (R) + {\hbar \over i} \sigma_1 (R) + \left(\frac{\hbar}{ i}\right)ˆ2 \sigma_2 (R) + \cdots \right]}$

Ordre 0

En n'utilisant que $σ 0$ dans $ψ$ on obtient immédiatement

$- {\sigma'}ˆ2_0 (R) + i \hbar \sigma''_0 (R) + pˆ2(R) = 0$

L'ordre 0, qui se nomme l'approximation classique, consiste à ne conserver aucun terme en $\hbar$ . On obtient $\sigma_0(R) = \pm \int p$ et donc

$\psi (R) \approx eˆ{ \frac{i}{ \hbar} \pm \int p }$

Ordre 1

L'ordre suivant est l'approximation B. K. W. elle-même.

En reprenant la formule précédente, avec $\sigma_0 + \frac{\hbar}{i} \sigma_1$ au lieu de $σ 0$ , et en ne gardant que les termes en $\hbar$ on obtient immédiatement $2σ' 0 σ' 1 + σ'' 0 = 0$

En utilisant la valeur de $\sigma_0 (R) =\pm \int p(R)$ , on en déduit $\sigma_1 (R) = cte + \frac{1}{2} \ln |p(R)|$ et finalement avec les deux signes envisageables de $σ 0$

$\psiˆ{\rm BKW} (R) = {C_1 \over \sqrt{ |p(R) | } } eˆ{ {i \over \hbar} \int p } + {C_2 \over \sqrt{ |p(R) | } } eˆ{- {i \over \hbar} \int p }$

Ordre 2

Le calcul à l'ordre 2 apporte

$\psi (R) = {C'_1 \over \sqrt{ |p(R) | } } \left[ 1 - {i m \hbar \over 4} { F(R) \over pˆ3(R) } - {i mˆ2 \hbar \over 8} \int { Fˆ2 \over pˆ5 } \right] eˆ{ {i \over \hbar} \int p } + {C'_2 \over \sqrt{ |p(R) | } } \left[ 1 + {i m \hbar \over 4} { F(R) \over pˆ3(R) } + {i mˆ2 \hbar \over 8} \int { Fˆ2 \over pˆ5 } \right] eˆ{- {i \over \hbar} \int p }$

où $F = - {d V \over d R}$ sert à désigner la force à laquelle est soumise la particule. Cette formule est rarement utilisée, mais en comparant avec la formule BKW, on voit que l'approximation BKW sera valide dans le cas où $\left| { m \hbar F(R) \over pˆ3(R) }\right| \ll 4$ et $\left| mˆ2 \hbar \int { Fˆ2(R) \over pˆ5(R) } \right| \ll 8$

On préfère fréquemment réécrire ces conditions en utilisant $F = p p'/ m$ ce qui amène aux conditions données auparavant qui sont :

$\left| {d \bar \lambda \over d R } \right| \ll 4 \hbox{ et } \left| \int { 1 \over \bar \lambda (R) } \left( {d \bar \lambda \over d R } \right)ˆ2 \right|ˆ2 \ll 4/\pi$

Cas des points tournants classiques (vitesse nulle)

Les points où $p (R) = 0$ sont nommés les points de retournements classiques, en effet la vitesse v=p/m y est nulle, le mobile (ou la particule) fera demi-tour. En ces points la première condition n'est plus valable et l'approximation BKW est complètement fausse et il est par conséquent indispensable d'effectuer un traitement spécial pour ces points.

Au voisinage d'un tel point $R C$ on peut écrire un développement limité du potentiel. En s'arrêtant au premier ordre où $V(R) \approx E - (R -R_C) F_C$ l'équation de Schrödinger devient une équation fonction d'Airy dont la solution est donnée par $\psi (R) = C_A {\rm Ai} \left[ \left( R_C- R \right) \left( { 2 m F_C \over \hbarˆ2} \right)ˆ{1/3} \right]$

Fonctions de connections

En utilisant les développement asymptotiques de la fonction d'Airy il est envisageable de les raccorder aux fonctions BKW de part et d'autre d'un point tournant.

Les raccordements de deux fonctions BKW s'en suivent et sont donnés par les lois suivantes.

La fonction ${C \over 2 \sqrt{ |p(R) | } } eˆ{ - {1 \over \hbar} \left| \int_{R_C}ˆR p(R') d R' \right| }$ dans la région $V (R) > E$ devient ${C \over \sqrt{ p(R) } } \cos\left( {1 \over \hbar} \left| \int_{R_C}ˆR p(R') d R' \right| -{\pi\over 4} \right)$ dans la région $V (R) < E$
La fonction ${C \over \sqrt{ p(R) } } eˆ{ {i \over \hbar} \left| \int_{R_C}ˆR p(R') d R' \right| + { i \pi \over 4} }$ dans la région $V (R) < E$ devient ${C \over \sqrt{ |p(R) | } } eˆ{ {1 \over \hbar} \left| \int_{R_C}ˆR p(R') d R' \right| }$ dans la région $V (R) > E$

Signalons qu'il est préférable de ne pas extrapoler d'autres formules car les termes exponentiellement croissant et exponentiellement décroissant ne peuvent généralement cœxister dans la région classiquement interdite.

Approximation semi-classique uniforme

Notons aussi la formule de l'approximation semi-classique uniforme (ASU) valable dans l'ensemble des régions, donnée par :

$\psiˆ{\rm ASU} (R) = {C_4 \over \pi} \left[ {3 \hbarˆ2 \int_{R_C}ˆR p \over 2 pˆ3} \right]ˆ{1/6} {\rm Ai}\left[ \left( {3 \int_{R_C}ˆR p \over 2 \hbar } \right)ˆ{2/3} \right]$

Etats liés

L'une des applications principales de la théorie BKW concerne le calcul des fonctions d'onde dans un puits de potentiel.

En notant $R int$ le point tournant classique interne et $R ext$ le point externe et en utilisant les formules de connections en ces deux points on s'aperçoit aisément que la somme des phases des cosinus doit être un multiple de π. On en déduit la condition de quantification, qui est en fait celle trouvée par Niels Bohr et Arnold Sommerfeld en 1913 dans l'ancienne théorie des quanta mais avec le 1/2 en plus ${ 1 \over \pi \hbar } \int_{R_{\rm int}}ˆ{R_{\rm ext} } p(R) d R = v + {1 \over 2}$ où v est le nombre de zéro de la fonction d'onde ψ du v+1 ème niveau lié du potentiel (théorème d'oscillation). L'approximation BKW s'écrit

$\psi (R) = \sqrt{ 2 \omega m \over \pi p(R) } \cos\left( {1 \over \hbar} \int_{R_{\rm int}}ˆR p(R') d R' -{\pi\over 4}\right)$

où on a normalisé à la fonction d'onde en négligeant la partie classiquement interdite et utilisant l'approximation de l'oscillation rapide du cosinus ( $\cosˆ2 \approx 1/2$ ).

ω sert à désigner la pulsation du mouvement classique et T est la période d'oscillation définies par $\omega = {2 \pi \over T} = {2 \pi \over 2 m \int_{R_{\rm int}}ˆ{R_{\rm ext} } pˆ{-1} (R) d R}$

Plus v est grand, plus p l'est , et par conséquent plus l'approximation BKW sera valable (voir la première condition de validité). Il convient tout de même d'être soigneux pour les tous derniers niveaux du potentiel, car l'approximation BKW n'est plus valable (voir la seconde condition de validité).

Notes

↑ Léon Brillouin La mécanique ondulatoire de Schrödinger; une méthode générale de resolution par approximations successives, Comptes rendus (Paris) 183 24-26 (1926)
↑ H. A. Kramers Wellenmechanik und halbzahlige Quantisierung, Z. Physik. 39 828-840 (1926)
↑ Gregor Wentzel Eine Verallgemeinerung der Quantenbedingungen für die Zwecke der Wellenmechanik, Z. Physik. 38 518-529 (1926).

Aspects récents

Harald Siegfried Friedrich ; Theoretical Atomic Physics, Springer; 3 edition (October 6, 2005), ISBN-13 : 978-3540256441

André Voros ; Spectre de l'équation de Schrödinger et méthode BKW, notes de cours (Orsay - 1980/1981), Publications Mathématiques de l'Université Paris-Sud Orsay 81-09 (1981). pdf.

André Voros ; De la théorie BKW exacte vers la théorie des perturbations singulières (II) , AOKI T., KOIKE T., MAJIMA H., TAKEI Y., TOSE N., eds. (Springer-Verlag, Berlin) Mondial Conference on Algebraic Analysis of Differential Equations (from Microlocal Analysis to Exponential Asymptotics), in honor of Prof. Takahiro Kawai on the occasion of his sixtieth birthday (invitation), Kyoto University, Kyoto, Japan, July 7-14 2005. ArXiv : math-ph/0603043.

Recherche sur Amazone (livres) :

Ce texte est issu de l'encyclopédie Wikipedia. Vous pouvez consulter sa version originale dans cette encyclopédie à l'adresse http://fr.wikipedia.org/wiki/Approximation_BKW.
Voir la liste des contributeurs.
La version présentée ici à été extraite depuis cette source le 12/04/2009.
Ce texte est disponible sous les termes de la licence de documentation libre GNU (GFDL).
La liste des définitions proposées en tête de page est une sélection parmi les résultats obtenus à l'aide de la commande "define:" de Google.
Cette page fait partie du projet Wikibis.