Unendlicher quadratischer Brunnen, der plötzlich an Größe abnimmt

Eine bekannte Übung in der grundlegenden Quantenmechanik ist die plötzliche (diabatische) Zunahme der Länge eines unendlichen quadratischen Brunnens .

Betrachten Sie nun ein Teilchen in einem Eigenzustand eines unendlichen Brunnens, dessen Länge plötzlich abnimmt . Auf den ersten Blick erscheint dies problematisch, da der Eigenzustand des Anfangstopfes nicht in die Eigenzustände des kleineren Topfes entwickelt werden kann.

Also dachte ich darüber nach, den anfänglichen Eigenzustand innerhalb des kleineren Brunnens wieder zu normalisieren, bevor ich die Überlappung berechnete. Das ist, was ich bisher bekommen habe:

Befindet sich das Teilchen zunächst im Grundzustand eines unendlichen Brunnens aus 0 Knirps 2 a , die renormierte Wellenfunktion in einem kleineren Brunnen aus 0 zu a ist

ϕ 1 = 2 a Sünde π x 2 a .
Die Eigenzustände des kleineren Wells sind
ψ n = 2 a Sünde n π x a ,
und so ist die Überlappung
c n = 2 a 0 a d x Sünde π x 2 a Sünde n π x a = ( 1 ) n 8 n π ( 1 4 n 2 )
Das Problem ist nun, dass der Erwartungswert der Energie divergiert,
E = n | c n | 2 E n = 32 2 m a 2 n n 4 ( 1 4 n 2 ) 2 .

Dies liegt daran, dass die Wellenfunktion bei diskontinuierlich ist x = a die erste Ableitung enthält also einen Term proportional zu δ ( x a ) .

Daher kam ich zu dem Schluss, dass diese Methode unbefriedigend ist, da sie nicht-physikalische Ergebnisse liefert.

Gibt es eine Lösung oder ist das Problem schlecht definiert? Die adiabatische Grenze ist sicherlich wohldefiniert.

BEARBEITEN

Dem Vorschlag von Kevin Zhou folgend, zog ich stattdessen dieses System in Betracht

                        

Zum E < v 0 , die Lösungen sind gegeben durch

ψ 1 ( x ) = Sünde k x ψ 2 ( x ) = Sünde k a Sünde λ a Sünde λ ( 2 a x ) ,
mit k = 2 m E / und λ = 2 m ( v 0 E ) / . Und das Spektrum wird aus gefunden
λ bräunen k a = k Tanh λ a ,
was numerisch gelöst werden muss. Zum E > v 0 wir lassen einfach λ ich q mit q = 2 m ( E v 0 ) / .

Huh, das ist eine interessante Frage. Auf den ersten Blick ist es eigentlich nicht verwunderlich, dass Sie unendlich sind, denn dies ist, als würde man das Teilchen unendlich hart treffen.
Sie könnten jedoch wahrscheinlich eine nicht unendliche Antwort erhalten, indem Sie das neue Potenzial als eine hohe Konstante annehmen v 0 statt unendlich. In diesem Fall sehen die Eigenzustände wie Eigenzustände des kleineren Brunnens aus, aber sie fließen in den über v 0 Region ein wenig. Dadurch können Sie ein vollkommen gutes inneres Produkt berechnen; Ich vermute, das Endergebnis wäre E v 0 , obwohl dies möglicherweise nicht analytisch erhältlich ist.
@KevinZhou Ich mag deine Idee. Die Skalierung des endlichen Brunnens ist, ähm, gut verstanden.
@KevinZhou Danke für deinen Vorschlag. Ich habe es numerisch gelöst und der Erwartungswert skaliert mit der Höhe v 0 wie erwartet.
Meine Intuition sagt mir, dass es keine zuverlässige Lösung für die Frage geben wird, weil die Situation physikalisch wirklich unmöglich ist. Zum einen gibt es keinen unendlichen quadratischen Brunnen. Außerdem kann man im Falle einer plötzlichen Expansion an eine sehr schnelle Expansion denken. Aber nicht für die Kontraktion. Selbst eine sehr schnelle Kontraktion wäre nicht adiabat. Wie Sie festgestellt haben, ist eine mögliche Analyse ( v 0 ) gibt eine unsinnige Antwort
Der Koeffizient für ϕ 1 ist falsch
@Korra Beachten Sie, dass diese Wellenfunktion in Bezug auf die kleinere Breite des Brunnens normalisiert ist a . 0 a d x 2 a Sünde 2 π x 2 a = 1
@Praan, aber wenn die Wellenfunktion für einen Brunnen mit einer Breite von 0 bis 2a gilt, warum sollte sie für die kleinere Brunnenbreite normalisiert werden?
@Korra Ich gehe davon aus, dass sich das Partikel anfangs im größeren Brunnen befindet und schrumpft den Brunnen dann schnell genug, sodass die Wellenfunktion gleich bleibt (bis zu einer Konstante). Die Normalisierung muss sich ändern, da das Teilchen nicht austreten kann. Die Wellenfunktion muss vor und nach der plötzlichen Änderung der Breite des Brunnens normalisiert bleiben.

Antworten (3)

Zu lang für einen Kommentar: Eine andere Möglichkeit besteht darin, den Fall des unendlichen quadratischen Topfpotentials mit einer sich bewegenden Wand zu betrachten. Dies kann genau für den Fall konstanter Geschwindigkeit gelöst werden, v (Wo für einen sich ausdehnenden Brunnen v ist psoitiv und für ein Contracting gut v ist negativ: Siehe Griffiths, Introduction to Quantum Mechanics (2005), Problem 10.1, der auf die Arbeit von Doescher und Rice Am J Phys {\bf 37} 1246 (1969) verweist (Griffiths berücksichtigt nur v positiv (ein sich ausdehnender Brunnen), aber Doescher und Rice machen deutlich, dass die Lösung auch für ein negatives gilt v ).

Wenn der ursprüngliche Brunnen breit ist a , ein vollständiger Satz von Eigenfunktionen sind

Φ n ( x , t ) 2 w Sünde n π w x e ich ( m v x 2 2 E n ich a t ) / ( 2 w )
mit w ( t ) a + v t ist die momentane Breite des Brunnens und E n ich = n 2 π 2 2 2 m a 2 ist der n t h erlaubte Energie des ursprünglichen Brunnens. Jede Lösung ist eine Linearkombination der Φ n ( x , t ) und die Koeffizienten können auf übliche Weise bestimmt werden.

Ich würde sagen, dass das nicht-adiabatische Problem schlecht definiert ist. Mal schauen.

Alle ersten Lösungen Ψ n ( x ) = EIN n Sünde ( n π x / 2 a ) einen Wahrscheinlichkeitsstrom haben Ψ ' Ψ Ψ Ψ ' überall gleich Null, so dass Sie sie im Prinzip an jedem Punkt schneiden können, ohne ein Wahrscheinlichkeitsleck zu haben. Das Problem besteht darin, dass die Randbedingungen für den Definitionsbereich jedes Hamiltonoperators – jedes selbstadjungierten Operators – eindeutig und mit Linearität kompatibel sein müssen.

Es hängt davon ab n , der zentrale Schnitt ( a / 2 , 3 a / 2 ) der ursprünglichen Eigenfunktion geht an 4 verschiedene Hamiltonianer:

n Mod 4 = 1 Ψ ( a / 2 ) = Ψ ' ( a / 2 ) , Ψ ( 3 a / 2 ) = Ψ ' ( 3 a / 2 )

n Mod 4 = 2 Ψ ' ( a 2 ) = 0 , Ψ ' ( 3 a 2 ) = 0

n Mod 4 = 3 Ψ ( a 2 ) = Ψ ' ( a 2 ) , Ψ ( 3 a 2 ) = Ψ ' ( 3 a 2 )

n Mod 4 = 0 Ψ ( a 2 ) = 0 , Ψ ( 3 a 2 ) = 0

Der Schnitt zu ( 0 , a ) geht an zwei verschiedene Hamiltonianer

n Mod 2 = 1 Ψ ( 0 ) = 0 , Ψ ' ( a ) = 0

n Mod 2 = 0 Ψ ( 0 ) = 0 , Ψ ( a ) = 0

Das ist einfacher, aber es ist immer noch das gleiche Problem: Die Eigenfunktionen einer Hamilton-Funktion liegen außerhalb des Geltungsbereichs der anderen. Dies sollte der Ursprung der Unendlichkeit sein, die Sie erhalten, wenn Sie alle Projektionen addieren, um einen Erwartungswert zu erhalten. Die Diskontinuität an der Grenze ist wirklich mathematisch irrelevant, weil Sie integrieren L 2 ( 0 , a ) ; Sie sollten ein endliches Ergebnis erhalten, wenn Sie den Satz von Funktionen in derselben Domäne überlappen und addieren, dh mit denselben Randbedingungen.

Es könnte sinnvoll sein zu behaupten, dass Sie die Hälfte des Spektrums zerstören können, denn das ist auch das, was wir im Ortsraum tun: Wir zerstören alle Eigenwerte | δ x > zum x größer als a . Aber das gleiche passiert mit dem zentralen Schnitt, und dort sind wir gezwungen, drei von vier Energiewerten zu zerstören. Das Argument sieht also nicht sehr stark aus.

Der Vorschlag von Kevin Zhou hat dagegen den Vorteil, dass wir versuchen könnten, alle Pole des Originalbrunnens und des halbierten zu kontrollieren. Die Polologie des quadratischen Wells wurde in Nussenzveig „The poles of the S-matrix of a rechteckig potential well or barrier“ Nuclear Physics, 11:409–521, 1959 untersucht . und in seinem Buch "Causality and Dispersion Relations" erneut aufgegriffen. Einige Details auf Spanisch können aus Kapitel 3.2 meiner Diplomarbeit erraten werden , aber der Artikel von B. Belchev, SG Neale und MA Walton bietet eine vollständigere Übersicht, und ich habe die Plots nicht in mein Dokument aufgenommen.

Eine unausgegorene Idee, die wir hier ausnutzen sollten, ist, dass wir die potenzielle Form verkleinern 2 a zu a kein Eigenzustand geht wirklich verloren. Was in der Streumatrix passiert, ist, dass sich der entsprechende Pol über die imaginäre Achse bewegt ich k bis es mit einem anderen Pol kollidiert, der aus einem reinen "antigebundenen" Zustand stammt ("Gamow-Vektoren" oder diskrete Zustände, die in den kontinuierlichen E> 0-Teil des Spektrums eingefügt sind, wenn ich mich richtig erinnere) und dann verlassen beide Pole symmetrisch das Imaginäre Achse werden Resonanzen in der komplexen Ebene sein (einige Autoren - ich - nennen auch Resonanzen zu den "antigebundenen" Zuständen). Sie können es in Abb. 5 des BNW-Artikels sehen.

Es könnte also sinnvoll sein, das Halbierungsproblem nicht nur als Funktion der Eigenwerte zu analysieren E n links, sondern auch der in die komplexe Ebene eingedrungenen Resonanzen. Und kombinieren Sie den Prozess der Halbierung des Potenzials mit der Erhöhung von v 0 .

(zu groß für einen Kommentar, zu unausgegoren für eine Antwort, sorry dafür ... aber es könnte ein guter Anfang sein)