Adiabatische Entwicklung der Überlagerung von Zuständen

Ich soll eine bestimmte Überlagerung von Eigenzuständen eines zeitabhängigen Hamiltonoperators finden H ( T ) . Der Hamiltonian hat die Form: H ( T ) = ich ( J σ ich σ ich + 1 + H ( T ) σ ich z ) . Periodische Randbedingung.

Ich habe folgenden Zustand eingenommen T = 0 ,

| Ψ ( 0 ) = | 0 ( 0 ) + e ich ϕ | 1 ( 0 )
Nun, der Hamiltonian ändert sich mit der Zeit, das weiß ich T = T , entwickelt sich dieser Zustand zu:
| Ψ ( T ) = e ich Φ 0 | 0 ( T ) + e ich ϕ e ich Φ 1 | 1 ( T )

Wo Φ ich sind die Summe aus dynamischen und Berry-Phasen.

Das Problem, das ich habe, ist, dass die Staaten | 1 ( T ) Und | 0 ( T ) die ich durch einige Numeriken (Lanczos-Algorithmus) erhalte, haben eine willkürliche Phase, die nicht für alle Diagonalisierungen gleich ist (ich muss mehrere Diagonalisierungen für unterschiedliche Werte von durchführen T ).

Aus diesem Grund bekomme ich nicht immer genau die Entwicklung des gleichen Zustands T .

Ich möchte den Überlagerungszustand für alle Zeiten T .

Welche Struktur hat der Hamiltonoperator? Arbeiten Sie im Interaktionsbild?
Ich habe die Hamiltonsche Form hinzugefügt. Ich arbeite nicht im internationalen Bild, es ist Schrödinger Bild.
Kannst du etwas genaueres zu den Zahlen sagen? Ich gehe von Ihrer Frage aus, dass Sie den Lanczos-Algorithmus verwenden, um die 2 Spin-Konfigurationen mit der niedrigsten Energie zu erhalten (obwohl es im Allgemeinen eine geben wird 2 N Eigenvektoren) und versuchen, nach der Berry-Phase zu suchen, aber das geht aus der Frage nicht hervor.
@catalogue_number Ja, ich verwende Lanczos, um die niedrigsten zwei Eigenzustände zu erhalten. Aber sie kommen jedes Mal mit einer willkürlichen Phase. Da ich eine Überlagerung dieser Zustände nehme, bedeutet diese willkürliche Phase für verschiedene Zustände, dass ich nicht die gleiche Überlagerung von Zuständen habe, die ich ursprünglich genommen habe.

Antworten (2)

Ich werde eine Teilantwort versuchen, um den Ball ins Rollen zu bringen.

Unter der adiabatischen Näherung werden die Phasen wie folgt berechnet:

Unter einem sich langsam ändernden Hamiltonian H ( T ) mit momentanen Eigenzuständen | N ( T ) und entsprechende Energien E N ( T ) , entwickelt sich aus dem Anfangszustand ein Quantensystem

| ψ ( 0 ) = N C N ( 0 ) | N ( 0 )

zum Endzustand

| ψ ( T ) = N C N ( T ) | N ( T )

wo die Koeffizienten die Phasenänderung erfahren

C N ( T ) = C N ( 0 ) e ich θ N ( T ) e ich γ N ( T )

mit der dynamischen Phase

θ M ( T ) = 1 0 T E M ( T ' ) D T '

und geometrische Phase

γ M ( T ) = ich 0 T M ( T ' ) | T ' | M ( T ' ) D T ' .

Die Probleme, die OP mit der zufälligen Phase der numerischen Diagonalisierung identifiziert hat, ruinieren den letzten Ausdruck. Insbesondere wenn eine zufällige Phase eingeführt wird, T | M ( T ) ist nicht mehr definiert. Dies ist nicht nur ein numerisches Problem, es ist grundlegend – es ist immer möglich (in der adiabatischen Grenze), eine U(1)-Eichtransformation einzuführen | N ( T ) e ich β ( T ) | N ( T ) das bezieht sich auf den gleichen physikalischen Zustand. β ( T ) ist eine langsam variierende glatte Funktion. Mit anderen Worten, die Phase ist an einem gegebenen "Langzeit"-Punkt nicht eindeutig.

Um mit Berrys Argument fortzufahren, müssen wir uns auf den Fall spezialisieren, dass der Hamilton-Operator von einer Reihe von Parametern abhängt R . In unserem Fall werden wir den Hamiltonoperator leicht verallgemeinern:

H ( H ) = ich J σ ich σ ich + 1 + H σ ich

Das Spektrum und die Eigenzustände von H sind nur vom 3-Vektor abhängig H , werden also umbenannt E ( H ) , | N ( H ) bzw.

Die Beerenphase ist γ N = 0 T ich N | T | N ( T ) D T = C ich N | H ich | N ( T ) D l ich , wobei der rechte Ausdruck ein Linienintegral über in ist H - Zwischenraum H ( T = 0 ) Und H ( T = T ) .

Nur im Sonderfall das H zu einem späteren Zeitpunkt auf seinen ursprünglichen Wert zurückkehrt, wird die Berry-Phase eichinvariant, also experimentell relevant/messbar.

Ich fasse nun das hier angeführte Argument zusammen . Vermuten C ist eine geschlossene Kurve, die eine Fläche begrenzt S . Definieren Sie das Beerenpotential A ich = ich N ( H ) | H ich | N ( H ) , so dass γ N = C A D H . Wenden Sie den Satz von Stokes an:

C A D H = S ich D [ N ( H ) | ich N ( H ) D H ich ]
= S ich D [ N ( H ) | ] ich N ( H ) D H ich + ich N ( H ) | D [ ich N ( H ) ] D H ich
= M S ich ich N ( H ) | M ( H ) M ( H ) | J N ( H ) D H ich D H J = S Ω ich J D H ich D H J

Dann mit

N | J H | M + N | H | J M = δ M N J E M + E M N | J M

Abhängigkeit zu beseitigen J M , kann man die Berry-Krümmung umschreiben als

Ω ich J = ich M N ich N ( H ) | M ( H ) M ( H ) | J N ( H )
= ich M N N | ich H | M M | J H | N ( J ich ) ( E N E M ) 2

was eichinvariant ist! bis zu einem insgesamt irrelevanten Faktor von  2 π

Kehren Sie zu Ihrer ursprünglichen Situation zurück, die dem Fall entspricht, dass die 3D H ( T ) Kurve ist auf die beschränkt z Achse finden wir folgendes -

  1. Das Gebiet von S ist Null.
  2. Die Berry-Krümmung ist gegeben durch Ω μ v = ich M N ich N | σ ich μ | M M | σ ich v | N ( μ v ) ( E N E M ) 2 .
  3. Sofern die Berry-Krümmung nicht irgendwo eine Singularität aufweist, ist die Berry-Phase Null. Es ist jedoch wahrscheinlich, dass es irgendwo eine Singularität geben wird H Raum.

Bitte nehmen Sie diese Berechnungen nicht als Gebot der Stunde, es können Faktoren von 2 oder fehlen ich .

Ihre Zustände haben aufgrund einer willkürlichen globalen Phase U ( 1 ) Symmetrie messen. Um es zu entfernen, müssen Sie ein gemeinsames Messgerät für alle Ihre Staaten reparieren. Ein bequemer Weg, dies zu tun, besteht darin, iterativ einen parallelen Transport aufzuerlegen, dh die Phase Ihrer zu modifizieren | N ich angeben, dass die Bedingung Ich bin N ich 1 | N ich = 0 ist erfüllt. Beachten Sie, dass in diesem Fall die (diskretisierte) Berry-Phase nicht als Summe der Phasendifferenzen zwischen Ihren nachfolgenden Zuständen entlang Ihrer geschlossenen Schleife im Hilbert-Raum angegeben wird, sondern einfach als Phasendifferenz zwischen Ihrem Anfangs- und Endzustand, ϕ = Ich bin ln N 0 | N N .