Approximation BKW - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Introduction

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 système 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 donc retrouver la mécanique classique dans la limite  \hbar \rightarrow 0 comme on retrouve l'optique géométrique lorsque la longueur d'onde  \lambda \rightarrow 0 dans la théorie de l'optique ondulatoire.


L'approximation BKW (pour les francophones) est également connue sous les initiales WKB (pour les anglophones), WKBJ, BWKJ et parfois 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 en général considérée 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 }

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

Sens physique

Notons le sens physique simple :

  1. Dans la région classiquement permise plus la particule va vite, plus sa probabilité de présence diminue. En effet là où E > V(R), la probabilité de présence | ψ | 2 sera proportionnelle à  1\over p .
  2. 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à où E < V(R), on a alors  p(R) = i  \sqrt{ | 2 m (E - V(R)) | } et le terme exponentiellement croissant sera en général divergent et donc 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

  1. \left| {d \bar \lambda \over d R } \right|  \ll  1 \bar \lambda(R) = \hbar / p(R) est la longueur d'onde de de Broglie dite réduite (divisée par 2 π)
  2. Une seconde condition, qui est souvent 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 comparables à 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 difficile à interpréter mais elle indique qu'il faut être prudent si le potentiel V décroit trop lentement à l'infini.

Cas des points tournants classiques (vitesse nulle)

Les points où p(R) = 0 sont appelé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 totalement fausse et il est donc nécessaire d'effectuer un traitement spécial pour ces points.

Au voisinage d'un tel point RC 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 possible 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.

  1. 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
  2. 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 en général coexister dans la région classiquement interdite.

Approximation semi-classique uniforme

Notons aussi la formule de l'approximation semi-classique uniforme (ASU) valable dans toutes les 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]

Page générée en 0.155 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales | Partenaire: HD-Numérique
Version anglaise | Version allemande | Version espagnole | Version portugaise