Ansatz von Landau & Lifshitz (Konturmethode) zu den WKB-Verbindungsformeln

Hintergrund der Frage (siehe S. 161, Abschnitt 47 im Lehrbuch der Quantenmechanik von Landau & Lifshitz, Bd. 3, 2. Auflage, Pergamon Press). Wir haben folgendes Potential gut

U ( X ) E für X A ,
(47,0) U ( X ) > E für X > A .

Die WKB-Lösungen rechts und links vom Wendepunkt sind

(47.1) ψ = C 2 P exp ( 1 | A X P D X | ) für X > A ,

(47.2) ψ = C 1 P exp ( ich A X P D X ) + C 2 P exp ( ich A X P D X ) für X < A ,

bzw. Die meisten Lehrbücher der Quantenmechanik bestimmen die Beziehung zwischen C Und C ich indem man die exakte Lösung in der Nähe des Wendepunkts findet. Und dann die exakte Lösung mit den WKB-Lösungen abgleichen lassen.

In den Lehrbüchern zur Quantenmechanik von Landau & Lifshitz (Bd. 3, Abschnitt 47) lassen sie jedoch zu X variieren in der komplexen Ebene und passieren den Wendepunkt A von rechts nach links durch einen großen Halbkreis in der oberen komplexen Ebene. Landau Ansprüche ab + , bei der Ankunft (links zu A ), gibt es eine Phasenverstärkung π im Nenner des Vorfaktors in Gleichung (2). Daraus können wir bestimmen

(47.4a) C 2 = C 2 exp ( ich π 4 ) .

Sie behaupten auch, dass der erste Term entlang des Halbkreises in der oberen Halbebene exponentiell zerfallen wird. Frage ist warum? Können wir zeigen

( A X P D X ) ,
Wo steht für den Imaginärteil, ist positiv?

Antworten (1)

I) L&L beziehen sich auf ein linearisiertes Modell, bei dem die TISE

(1) 2 ψ ' ' ( X ) + P 2 ( X ) ψ ( X )   =   0 , P ( X )   :=   2 M ( E v ( X ) )   =   | P ( X ) | e ich θ ( X ) ,

wird zur Airy-Differentialgleichung . Kurz vor einem Wendepunkt X 0 , können wir den quadratischen Impuls approximieren

(2) P 2 ( X )     a ( X X 0 ) , ( P 2 ) ' ( X )     a   =   2 M F 0     0 ,

mit einer affinen Funktion , wo a R { 0 } ist eine reelle Konstante ungleich Null. (Im Beispiel von L&L a < 0 ist negativ).

II) Lassen Sie uns zunächst das Argument von L&L wiedergeben. Im klassisch erlaubten Bereich X < X 0 , L&L wählen das Momentum P > > 0 positiv. (Die hochgestellten > Und < beziehen sich auf die klassisch erlaubten bzw. verbotenen Bereiche.) Daraus wird die WKB-Wellenfunktion (47.2).

(3) ψ >   = ( 47.2 )   1 P > ± C ± exp ( ± ich X 0 X P >   D X )   =   1 | P | ± C ± exp ( 2 ich | a | 3 ( X 0 X ) 3 2 ) für X < X 0 .

III) L&L setzt als nächstes die WKB-Wellenfunktion (47.1) (die im klassisch verbotenen Bereich gilt) entlang eines (oberen/unteren) Halbkreises im Komplex analytisch fort X -Ebene

(47.4a) X X 0   =   ρ e ± ich ϕ , ϕ     [ 0 , π ] ,

um die WKB-Verbindungsformeln abzuleiten . L&L wählen die Impulsvariable

(4) P <   =   | a | ( X X 0 )   =   | a | ρ 1 2 e ± ich ϕ 2 , ρ   >   0 ,

damit es positiv anfängt

(5) P < | ϕ = 0   =   | P |   >   0.

Die analytische Fortsetzung der WKB-Wellenfunktion (47.1) wird dann

(6) ψ <   = ( 47.1 )   C 2 P < exp ( 1 X 0 X P <   D X )   =   C 2 | P | exp ( 2 | a | ρ 3 2 e ± ich 3 ϕ 2 3 ich ϕ 4 ) .

Indem man es einstellt ϕ = π in Gl. (6) und im Vergleich mit der WKB-Wellenfunktion (3) folgt aus der analytischen Fortsetzung entlang des (oberen/unteren) Halbkreises, dass

(47.4b) C 2     C   =   C 2 e ich π 4 , (47.4c) C 1     C +   =   C 2 e ich π 4 ,

bzw. Das ist das Hauptargument von L&L. (Wir werden auf die Frage, was mit dem anderen Zweig geschieht, in Abschnitt VII unten zurückkommen.)

IV) Alternativ kann man, da die Airy-Funktionen Fourier-Integraldarstellungen haben (vgl. z. B. meine Math.SE-Antwort hier ), stattdessen die Methode des steilsten Abstiegs verwenden , um passende asymptotische Erweiterungen auf jeder Seite des Wendepunkts abzuleiten. Es ist interessant, diese Methode mit dem obigen Argument von L&L zu vergleichen.

Man kann zeigen, dass die Airy-ähnliche Wellenfunktion

(7) ψ ( X )   =   | a | π R ich   S G N ( a )   ε D P   exp ( ich S ( P , X ) ) ,
(8) S ( P , X )   :=   P X S ~ ( P ) , S ~ ( P )   :=   P ( X 0 + P 2 3 a ) ,
(9) S ' ( P , X )   =   X S ~ ' ( P )   =   X X 0 P 2 a   ( 2 )   P 2 ( X ) P 2 a ,
(10) S ' ' ( P , X )   =   S ~ ' ' ( P )   =   2 P a ,

erfüllt die TISE (1):

(11) 2 ψ ' ' ( X ) + P 2 ( X ) ψ ( X )   = ( 7 )   | a | π R ich   S G N ( a )   ε D P   ( P 2 ( X ) P 2 ) exp ( ich S ( P , X ) )   ( 9 )   ich a | a | π R ich   S G N ( a )   ε D P   D D P exp ( ich S ( P , X ) )   =   0

für die gewählte Kontur im Komplex P -Ebene.

V) Es gibt 2 kritische Punkte:

(12) P σ   =   σ   S G N ( a )   P , S σ   =   2 σ P 3 3 | a | , S σ ' '   =   2 σ | a | P , σ     { ± 1 } .

Der steilste Abfallbeitrag von jedem kritischen Punkt kommt von einem Gaußschen Integral :

(13) ICH σ   =   | a | π 2 π ich S σ ' ' exp ( ich S σ )   =   1 P exp ( 2 ich σ P 3 3 | a | ich σ π 4 )   =   1 | P | exp ( 2 ich σ | P | 3 e 3 ich θ 3 | a | ich ( σ π 4 + θ 2 ) )

mit winkelig steilster Abstiegsrichtung σ π 4 θ 2 .

Abb. 1. Der Komplex P -Ebene mit 3 möglichen Integrationskonturen C 1 , C 2 , C 3 für a < 0 . Die schattierten Bereiche bezeichnen exponentiell abfallende Sektoren. (Wenn a > 0 , es ist umgekehrt.) Da der Integrand eine ganze Funktion in ist P , Die P -Integrale hängen nur von der Monodromie der Kontur ab. (Abbildung aus Lit. [W].)

VI) Kehren wir nun mit zum Beispiel von L&L zurück a < 0 .

Klassische verbotene Region X > X 0 : Dann P = ± ich | P | , mit θ = ± π 2 , Wo ± entspricht zwei verschiedenen möglichen Vorzeichenwahlen. Die Kontur entlang der realen P -Achse wird durch die dem kritischen Punkt entsprechende steilste Abstiegskontur wiedergegeben P 1 :

(14) ψ < ( X )     ICH 1   =   1 | P | exp ( 2 | P | 3 3 | a | ) .

Klassisch erlaubte Region X < X 0 : Dann P = | P | > 0 mit θ = 0 . Die Kontur entlang der realen P -Achse ergibt sich aus der Summe der beiden steilsten Abstiegskonturen:

(15) ψ > ( X )     σ { ± 1 } ICH σ   =   1 | P | σ { ± 1 } exp ( ich σ ( 2 | P | 3 3 | a | π 4 ) )   =   2 | P | cos ( 2 | P | 3 3 | a | π 4 )   =   2 | P | Sünde ( 2 | P | 3 3 | a | + π 4 ) .

Durch den Vergleich von Gl. (14) & (15) haben wir eine Verbindungsformel hergeleitet. Im Detail, indem man die Beiträge zu den kritischen Punkten vergleicht P 1 , leiten wir die Gleichungen von L&L ab. (47.4c) bzw. (47.4b). Die beiden Methoden stimmen also überein.

VII) Exponentiell wachsende Lösungen sind physikalisch möglich, wenn der klassisch verbotene Bereich endliche Länge hat, zB beim Quantentunneln. Im Beispiel von L&L ist der klassisch verbotene Bereich jedoch nicht kompakt, sodass exponentiell wachsende Lösungen verworfen werden müssen.

Es bleibt zu erklären, warum nur einer der beiden Schwingungszweige auf der klassisch erlaubten Seite analytisch auf die klassisch verbotene Seite fortgesetzt wird.

Einerseits argumentiert L&L, dass dem exponentiell wachsenden WKB-Zweig während des Jahres nicht vertraut werden kann 3 π 2 Phasenverschiebung im Boltzmann-Faktor, weil er momentan exponentiell unterdrückt wird. Dies unterstreicht die Tatsache, dass die Verbindungsformeln unidirektional sind.

Andererseits kreuzt man aus Sicht der Methode des steilsten Abstiegs in den Abschnitten IV-VI eine Stokes-Linie im Komplex X -Ebene, so dass die steilsten Abstiegskonturen im Komplex Monodromie aufnehmen P -Flugzeug, vgl. Abb. 1.

Weitere Informationen finden Sie auch z. B. in meiner verwandten Phys.SE-Antwort hier .

VIII) Schließlich sollten wir erwähnen, dass die π 2 Phasenverschiebung zwischen einlaufender und auslaufender Welle am Wendepunkt in Gl. (15) entspricht einem Maslov-Index | μ | = 1 , vgl. zB dieser Phys.SE Beitrag.

Verweise:

  • [LL] LD Landau & EM Lifshitz, QM, Bd. 3, 2. und 3. Auflage, 1981; § 47 & § B .

  • [W] E. Witten, Analytische Fortsetzung der Chern-Simons-Theorie, arXiv:1001.2933 ; P. 23-29, 48-49. Ein verwandter KITP-Vortrag von Witten aus dem Jahr 2015, A New Look At The Path Integral Of Quantum Mechanics, ist auf YouTube zu finden .

L&L schlagen vor, von rückwärts zu gehen π Zu 0 im Flugzeug j > 0 . Dann C 2 P e ich P ( X ) D X wird klein. Um den Begriff zu machen C 1 P e + ich P ( X ) D X , klein muss man das gleiche tun j < 0 . Könnten Sie bitte kommentieren, wie man dies möglicherweise durch Betrachten ableiten könnte C 1 P e + ich P ( X ) D X + C 2 P e ich P ( X ) D X . Ich verstehe nicht, wie die Formel X A = ρ e ich ϕ wird hier angewendet, wenn überhaupt. Ich verstehe es nicht wirklich aus deiner Antwort. Vielen Dank im Voraus
Danke für die Rückmeldung. Ich plane, dies in naher Zukunft zu überprüfen.
Ich habe es ein wenig gelehrt. Ich denke, das ist ganz einfach, wie Landau selbst schreibt. Wenn man weitermachen würde C 1 P ( X ) e ich P ( X ) D X in der komplexen Ebene via X A = ρ e ich ϕ man bekommt e Sünde 3 ϕ 2 + ich cos 3 ϕ 2 . Dasselbe kann für getan werden C 1 P ( X ) e ich P ( X ) D X = e + Sünde 3 ϕ 2 ich cos 3 ϕ 2 . Die Summe dieser beiden Terme würde das gewünschte Ergebnis ergeben, wenn man davon ausgehen würde 0 , 2 π oder rückwärts 2 π , 0 .