Die semiklassische WKB-Näherung aus der Quantenmechanik (benannt nach Gregor Wentzel, Hendrik Anthony Kramers und Léon Brillouin) liefert eine Näherung der Lösung der eindimensionalen, stationären Schrödingergleichung. Die Näherung basiert auf der Annahme, dass sich das Potential $ V(x) $ nur 'langsam' mit der Position, d. h. über die Ausdehnung einer Wellenlänge, ändert und sich daher eine Lösung aus dem konstanten Potential $ V(x)=V_{0} $ finden lässt.
Unter dieser Voraussetzung lautet die genäherte Lösung der Schrödingergleichung:
Die beiden Vorzeichen stehen für zwei unabhängige Lösungen.
Die Näherung wurde 1926 fast gleichzeitig und unabhängig voneinander von den Physikern Gregor Wentzel, Hendrik Anthony Kramers und Léon Brillouin im Rahmen der Quantenmechanik publiziert, deren Initialen ihr den Namen gaben. Sie findet sich aber auch schon vorher in den Arbeiten verschiedener Mathematiker und Physiker wie Francesco Carlini (1817, in der Himmelsmechanik), George Green (1837), Joseph Liouville (1837), John William Strutt, 3. Baron Rayleigh (1912), Richard Gans (1915), Harold Jeffreys[1] (1923)[2][3] verwendet. Sie wird deshalb auch manchmal WKBJ (zusätzlich nach Jeffreys) oder Liouville-Green-Methode genannt. Auch Werner Heisenberg benutzte das Verfahren 1924 in seiner Dissertation über Hydrodynamik.[4]
Aus der eindimensionalen stationären Schrödinger-Gleichung
ergibt sich bei konstantem Potential $ V(x)=V_{0} $ als Lösung die ebene Welle
mit $ p_{0}={\sqrt {2m(E-V_{0})}} $. Bei langsamer Änderung des Potentials, also einem Potential, das in der Größenordnung der deBroglie-Wellenlänge als konstant angesehen werden kann, kann man $ p(x)={\sqrt {2m(E-V(x))}} $ annehmen und daraus einen zum Problem mit konstantem Potential analogen Lösungsansatz folgendermaßen wählen.
Eingesetzt in die Schrödinger-Gleichung erhält man
Soweit wurde keine Näherung gemacht. Wir können nun $ S(x) $ folgendermaßen in Potenzen von $ \hbar $ entwickeln
Das setzt man in die Schrödingergleichung ein:
Nun kann man diese Terme bis zur gewünschten Ordnung berechnen und nach der Potenz von $ \hbar $ sammeln.
Jeder zu einer Potenz von $ \hbar $ zugehörige Term muss dann einzeln verschwinden.
Für die zweite Ordnung lautet die Schrödingergleichung:
Für die Differentialgleichung im Glied nullter Ordnung in $ \hbar $
findet man eine Lösung durch
und es folgt
Dieses Ergebnis beschreibt Lösungen einer eindimensionale Schrödingergleichung im Grenzfall $ \hbar \rightarrow 0 $, welche in der Punktmechanik mit der Hamilton-Jacobi-Gleichung ein gleichwertiges Gegenstück besitzt.
Die Taylorreihenentwicklung der Exponentialfunktion zeigt jedoch die mathematische Inkonsistenz dieser einfachsten Näherung im Rahmen der klassischen Konvergenzkriterien, da jeder auftretende Summand aufgrund der Division durch $ \hbar =0 $ divergiert und damit auch die Summe ohne passende Regularisierung nicht wohldefiniert ist. Außerdem ist in dieser Näherung die Beschreibung des Tunnelns $ (V(x)>E) $ problematisch, da einerseits zwar eine Näherungslösung für eine Schrödingergleichung konstruiert werden soll, für deren Lösung die Gültigkeit des Bornschen Wahrscheinlichkeitspostulats angenommen wird, andererseits aber der Limes $ \hbar \rightarrow 0 $ ein klassisches Hamiltonsches Wirkungsprinzip unterstellt, welches unvereinbar mit intrinsisch quantenmechanischen Tunnelprozessen ist.
Eine zusätzliche Betrachtung erster Ordnung der Wirkungsfunktion in $ \hbar $ fixiert die Konstante A und komplettiert die semiklassische WKB-Näherung. Eine genaue Rechnung zeigt, dass diese Näherung nur dann gut ist, wenn der Impuls $ p_{0}=\hbar /\lambda $ der Wellenfunktion wesentlich größer ist, als die örtliche Variation des Potenzials (siehe oben). Dies zeigt trotz der vollzogenen Wirkungsdiskretisierung, aufgelöst bis zur ersten Ordnung, die Nähe des Ansatzes zur klassischen strahlenoptischen Näherung, die in der geometrischen Optik durch das Eikonal und in der Hamilton-Jacobi-Theorie Anwendung findet.
Die WKB-Approximation wird benutzt, um nichtrechteckige Barrieren zu nähern. Dazu wird die Barriere in viele dünne rechteckige Teilbarrieren zerlegt.
Für die Tunnelwahrscheinlichkeit $ \left|T\right|^{2} $ durch diese Potentialbarriere werden die einzelnen Tunnelwahrscheinlichkeiten für jedes Segment multipliziert. Damit ergibt sich
wobei
Wie im vorherigen Abschnitt erläutert, muss die Barriere im Vergleich zur Ausdehnung der Wellenlänge der Materiewelle stückweise konstant sein, um die Näherung zu rechtfertigen. Für hinreichend große Wellenlängen, also z. B. nahe der Umkehrpunkte klassischer Teilchenbewegungen mit $ p\propto 1/\lambda \rightarrow 0 $, kann dies nicht mehr der Fall sein und es muss in diesen Regionen ein stetiger Anschluss durch exakte Lösungen, sogenannte Airy-Funktionen, erfolgen.