Neutrales Quantenteilchen im inhomogenen Magnetfeld

Ich versuche, das Stern-Gerlach-Experiment auf einer Computerebene zu verstehen. Angenommen, wir haben ein neutrales Teilchen mit magnetischem Moment (z. B. ein Neutron) und wenden ein inhomogenes Magnetfeld an (es soll sich linear mit der Koordinate ändern). Soweit ich weiß, würde der Hamiltonian so aussehen:

H ^ = 2 2 M 2 + ( e M C ) S ^ B

Jetzt ist der Spin-Operator

S ^ ich = 2 σ ich ,

Wo σ ich Ist ich te Pauli-Matrix .

Also für Magnetfelder B = e X B 0 X Wir hätten die Schrödinger 1D-Gleichung (Y- und Z-Richtung können aufgrund der Translationssymmetrie getrennt werden):

2 2 M 2 ψ X 2 + ( e 2 M C ) σ X B 0 X ψ = ich ψ T .

Ich versuche nun, diese Gleichung numerisch zu lösen, indem ich die anfängliche Wellenfunktion in der folgenden Form annehme:

ψ ( X , T = 0 ) = ( ψ 0 ( X ) ψ 0 ( X ) ) ,

Wo ψ 0 ( X ) ist ein Gaußsches Wellenpaket mit einem durchschnittlichen Impuls von Null.

Die Probleme beginnen, wenn ich wähle σ X wie gewöhnlich angegeben:

σ X = ( 0 1 1 0 ) .

Die Lösung scheint wie unten gezeigt auszusehen. Dh beide Wellenfunktionskomponenten beschleunigen nach links!

Geben Sie hier die Bildbeschreibung ein

Ich dachte, was wäre, wenn ich eine andere Achse als wählen würde X , also habe ich versucht, dasselbe mit zu tun σ j :

σ j = ( 0 ich ich 0 ) .

Das Ergebnis in der Animation unten. Jetzt ist es etwas besser: Die Wellenfunktion teilt sich zumindest in zwei Teile auf, einen nach links, einen nach rechts. Trotzdem bestehen beide Teile aus einer Mischung von Spin-up- und Spin-down-Zuständen, also nicht wirklich das, was man von einem Stern-Gerlach-Experiment erwarten würde.

Geben Sie hier die Bildbeschreibung ein

Schließlich habe ich die letzte Option ausprobiert – using σ z :

σ z = ( 1 0 0 1 ) .

Das Ergebnis ist unten noch einmal dargestellt. Schließlich bekomme ich die Aufteilung in "unabhängige" Spin-Teile, dh ein Spin-Teil geht nach links, ein anderer nach rechts.

Geben Sie hier die Bildbeschreibung ein

Nun die Frage : Wie sind diese Ergebnisse zu interpretieren? Warum führt die Wahl der aktiven Achse zu solch drastischen Unterschieden in den Ergebnissen? Wie hätte ich stattdessen vorgehen sollen, um aussagekräftige Ergebnisse zu erhalten? Sollte die Permutation von Pauli-Matrizen die Ergebnisse nicht beeinflussen?

Antworten (1)

Ich werde einige kosmetische Änderungen an Ihrer Gleichung vornehmen:

  1. Ihr Magnetfeld sollte sein B = e X B 0 X / X 0 , so dass B Und B 0 gleiche Maße haben. Dieses Feld gehorcht nicht B = 0 , aber das lassen wir jetzt mal.

  2. ich werde schreiben μ σ B für die magnetische Energie. Beachten Sie, dass das magnetische Moment des Neutrons ziemlich klein ist, etwa 50 neV/T. Sogenannte "ultrakalte" Neutronen mit Energien unter etwa 100 neV können durch ein Magnetfeldminimum eingeschlossen werden, aber für kalte oder thermische ~meV-Neutronen ist die Stern-Gerlach-Steuerung im Allgemeinen vernachlässigbar. (Wenn es nicht vernachlässigbar wäre, könnten Sie die Stern-Gerlach-Trennung verwenden, um einen Neutronenstrahl in zwei polarisierte Strahlen umzuwandeln, die in verschiedene Richtungen gehen; leider absorbieren echte Neutronenpolarisatoren alle einen Spinzustand.)

  3. Ihre Bewegungsgleichung ist auch zeitlich trennbar, daher werde ich nur den räumlichen Teil davon betrachten; wir können einen Faktor von anhängen e ich ω T später.

  4. Ich werde verwenden ϕ Und χ um den beiden Komponenten Ihres Spinors getrennte Namen zu geben.

In Matrixschreibweise lautet Ihre zeitunabhängige Schrödinger-Gleichung

( 2 2 M X 2 μ B 0 X X 0 μ B 0 X X 0 2 2 M X 2 ) ( ϕ χ ) = E ( ϕ χ ) .

Das macht es ein wenig offensichtlicher, was passiert. Wenn Sie Ihr anfängliches Ensemble definiert haben ϕ = χ = ψ 0 ( X ) , versetzen Sie Ihr System in einen Eigenzustand von σ X ! Natürlich bewegt sich das gesamte Ensemble zusammen: Es ist entlang der Feldrichtung polarisiert! Wenn Ihre Anfangsbedingung ist ϕ = χ , sollte sich ein Paket in die andere Richtung bewegen.

Wenn Sie ersetzen σ X mit σ j , z , ändern Sie effektiv die Richtung des Feldes. Denken Sie daran, dass die Energie ist μ σ B ; wenn der einzige Term ungleich Null in diesem Produkt ist σ j B 0 X / X 0 , teilen Sie Ihrem Modell mit, dass die Feldstärke mit variiert X aber Punkte entlang der j Richtung. (Da diese Felder gehorchen B = 0 , Sie sollten sie sowieso bevorzugen.) In Ihre zweite und dritte Zahl setzen Sie also ein X -polarisierte Probe in einem Feld entlang j oder entlang z , und es trennt sich. Die Trennung erfolgt eher entlang der Richtung des Feldstärkegradienten als entlang der Richtung des Feldes; die UCN-Trapper sprechen von "starken Feldsuchern" und "schwachen Feldsuchern".

Im Feld entlang zeigend z , Ihre Probe trennt sich in Ihre Spinorbasis, wodurch Sie die schönen Farben erhalten.

Im Feld entlang zeigend j , Ihre Probe trennt sich immer noch. Aber anstatt sich in die zu trennen ϕ , χ Grundlage, die Sie für Farben verwenden, sind die diagonalen Zustände ϕ ± ich χ . Dies führt zu Interferenzen, wenn sich die beiden Wellenpakete überlappen.

Wow, das ist interessant. Könnten Sie erläutern, warum Einstellung ϕ = χ macht die Wellenfunktion zu einem Eigenzustand von σ X ? Ich dachte, ein Eigenzustand wäre, wenn ich zB setze χ = 0 , ϕ 0 .
Hausaufgabe: Zeigen Sie, dass die Eigenvektoren von σ X Sind ( 1 , ± 1 ) , und dass die Eigenvektoren von σ j Sind ( 1 , ± ich ) . Finde die Eigenvektoren von σ z . :-)
OK, ich scheine langsam zu verstehen. Aber warum trennt es sich in Spinorbasis für z -gerichtetes Feld, aber nicht für j -gerichtet? Sollte es nicht genauso sein, denn wir könnten das Bezugssystem umdrehen e X beim Schalten j -Feld zu z -Feld und hat sich nichts geändert?
Ah, ich scheine zu verstehen: Dieser Unterschied ist nicht messbar: Die Wahrscheinlichkeitsdichte für Elektronen ändert sich nicht, also ist das kein wirkliches Problem.
Ich denke, dass ich dir zustimme: der einzige Grund ist die Reaktion auf das Feld entlang z besonders aussieht, ist, dass Sie zufällig das gewählt haben z Achse für Ihre Darstellung.