Erwartungswert in zweiter Quantisierung

Ich stecke fest bei der Berechnung eines einfachen Erwartungswerts für einen Operator, der sich in einer zweiten Quantisierung ausdrückt. Ich kenne das Ergebnis, aber ich kann es nicht beweisen.

Nehmen wir an, ich habe eine Ein-Teilchen-Wellenfunktion | ϕ N gegeben von | ϕ N = J = 1 K | a J A J , N , Wo K ist die Anzahl der Orbitale/Stellen im System und die A J , N sind die Wahrscheinlichkeitsamplituden der Orbitale | a J . Der Index N kennzeichnet die Teilchen im System, von denen wir haben N .

Die Orbitale sind orthonormal, dh

J A J , M A J , N = δ M , N .

Lassen Sie uns alle Spin-Freiheitsgrade ignorieren. Die Vielteilchenwellenfunktion ist nun gegeben durch

| Ψ = ( N = 1 N J = 1 K C ^ J A J , N ) | vac ,
bei dem die C ^ J ist der übliche Erstellungsoperator vor Ort J .

Was ich jetzt berechnen möchte, ist der Erwartungswert

Ψ | H hüpfen | Ψ
mit
H hüpfen = T J = 1 K C ^ J + 1 C ^ J + H . C .
der übliche hüpfende Hamiltonian. Ich habe das starke Gefühl (und eine Beispielrechnung hat dies unterstützt), dass das Ergebnis gerecht ist
Ψ | H hüpfen | Ψ = T N = 1 N J = 1 K ( A J , N A J + 1 , N + A J + 1 , N A J , N )

Ich denke, dieses Ergebnis hängt trivial mit den Slater-Condon-Regeln zusammen, aber ich sehe die Verbindung nicht. Außerdem verzichte ich darauf, den Erwartungswert explizit zu berechnen, der die Summe und die Produkte der Erzeugungs-/Vernichtungsoperatoren enthält.

Was ist ein guter Weg, um mein Ergebnis zu beweisen?

Hinweis: Verwenden Sie die CCR-s, ​​die mit dem verknüpft sind C ^ J -S. Lässt Ihr Rechenbeispiel nicht schon so etwas vermuten?
Nun, ich wette, ich muss irgendwann die Kommutierungsrelationen verwenden. Aber ich sehe wirklich nicht ein, wie ich zu dem allgemeinen Ergebnis kommen soll. Da das Skalarprodukt linear ist, reicht es sicherlich aus, nur einen Term von zu betrachten H hüpfen , also kann man vielleicht wlog setzen J = 1 . Aber dann habe ich noch jede Menge C ^ von dem | Ψ (Slater-Determinanten). Wie gehe ich vor?

Antworten (1)

Zur späteren Bequemlichkeit bezeichnen

ϕ ^ N = J A J , N C ^ J
χ ^ N = J A J , N C ^ J + 1
Der Durchschnitt, den Sie berechnen möchten, lautet dann
Ψ | H ^ 0 | Ψ = T 0 | N ϕ ^ N ( J C ^ J + 1 C ^ J ) N ϕ ^ N | 0
Beginnend mit den regulären CCRs,
[ C ^ J , C ^ k ] = δ J k , [ C ^ J , C ^ k ] = [ C ^ J , C ^ k ] = 0
erhalten
[ J C ^ J + 1 C ^ J , ϕ ^ N ] = χ ^ N
Und
ϕ ^ N χ ^ M = ± χ ^ M ϕ ^ N
ϕ ^ N χ ^ M = ± χ ^ M ϕ ^ N + J A J + 1 , N A J , M
Verwenden Sie nun im Durchschnitt, den Sie berechnen müssen, das obige, um sich sukzessive zu bewegen ( J C ^ J + 1 C ^ J ) an den Orbitaloperatoren rechts vorbei:
( J C ^ J + 1 C ^ J ) N ϕ ^ N = ϕ ^ 1 ( J C ^ J + 1 C ^ J ) N > 1 ϕ ^ N + χ ^ 1 N 1 ϕ ^ N =
= M = 1 M = 2 ϕ ^ M ( J C ^ J + 1 C ^ J ) N > 2 ϕ ^ N + χ ^ 1 N 1 ϕ ^ N + ϕ ^ 1 χ ^ 2 N > 2 ϕ ^ N =
= M = 1 M = 2 ϕ ^ M ( J C ^ J + 1 C ^ J ) N > 2 ϕ ^ N + χ ^ 1 N 1 ϕ ^ N ± χ ^ 2 N 2 ϕ ^ N = =
= M ϕ ^ M ( J C ^ J + 1 C ^ J ) + M ( ± 1 ) M 1 χ ^ M N M ϕ ^ N
Der erste Term vernichtet das RHS-Vakuum, sodass nur die Summe übrig bleibt. Bringen Sie nun die Orbitaloperatoren auf der linken Seite ein und drehen Sie sie um χ ^ M an ihnen vorbei:
( N ϕ ^ N ) χ ^ M = ± ( N > 1 ϕ ^ N ) χ ^ M ϕ ^ 1 + ( N 1 ϕ ^ N ) J A J + 1 , 1 A J , M =
( N > 1 ϕ ^ N ) χ ^ M l = 1 l = 2 ϕ ^ l + ( N 1 ϕ ^ N ) J A J + 1 , 1 A J , M ± ( N 2 ϕ ^ N ) J A J + 1 , 2 A J , M = =
= χ ^ M ( N ϕ ^ N ) + l ( ± 1 ) l 1 ( N l ϕ ^ N ) J A J + 1 , l A J , M
Der erste Term vernichtet nun das linke Vakuum, und nachdem alles ersetzt wurde, wird der gewünschte Durchschnitt
Ψ | H ^ 0 | Ψ = T J , l , M ( ± 1 ) M 1 ( ± 1 ) l 1 A J + 1 , l A J , M 0 | ( N l ϕ ^ N ) ( N ' M ϕ ^ N ' ) | 0 =
Ψ | H ^ 0 | Ψ = T J , l , M ( ± 1 ) M 1 ( ± 1 ) l 1 A J + 1 , l A J , M δ l , M
und schlussendlich
Ψ | H ^ 0 | Ψ = T J , M A J + 1 , M A J , M
Es mag einige Fehler geben, die ich übersehen habe, aber das ist die allgemeine Idee.