Для удобства в будущем обозначим
ф^†н"="∑ДжАдж , нс^†Дж
х^†н"="∑ДжАдж , нс^†j + 1
Среднее значение, которое вы хотите рассчитать, читает тогда
⟨ Ψ |ЧАС^0| Ψ⟩=-т⟨ 0 |∏нф^н(∑Джс^†j + 1с^Дж)∏нф^†н| 0⟩
Начиная с обычных CCR-ов,
[с^Дж,с^†к]∓"="дельтаj k,[с^Дж,с^к]∓"="[с^†Дж,с^†к]∓= 0
получать
[∑Джс^†j + 1с^Дж,ф^†н] =х^†н
и
ф^†нх^†м= ±х^†мф^†н
ф^нх^†м= ±х^†мф^н+∑ДжА*j + 1 , пАДж , м
Теперь в среднем вам нужно рассчитать, использовать вышеизложенное для последовательного перемещения
(∑Джс^†j + 1с^Дж)
мимо орбитальных операторов справа:
(∑Джс^†j + 1с^Дж)∏нф^†н"="ф^†1(∑Джс^†j + 1с^Дж)∏п > 1ф^†н+х^†1∏п ≠ 1ф^†н"="
"="∏м = 1м = 2ф^†м(∑Джс^†j + 1с^Дж)∏п > 2ф^†н+х^†1∏п ≠ 1ф^†н+ф^†1х^†2∏п > 2ф^†н"="
"="∏м = 1м = 2ф^†м(∑Джс^†j + 1с^Дж)∏п > 2ф^†н+х^†1∏п ≠ 1ф^†н±х^†2∏п ≠ 2ф^†н= ⋯ =
"="∏мф^†м(∑Джс^†j + 1с^Дж) +∑м( ± 1)м - 1х^†м∏п ≠ мф^†н
Первый член аннулирует правый вакуум, поэтому остается только сумма. Теперь введите орбитальные операторы слева и переверните каждый
х^†м
мимо них:
(∏нф^н)х^†м= ± (∏п > 1ф^н)х^†мф^1+ (∏п ≠ 1ф^н)∑ДжА*j + 1 , 1АДж , м"="
(∏п > 1ф^н)х^†м∏л = 1л = 2ф^л+ (∏п ≠ 1ф^н)∑ДжА*j + 1 , 1АДж , м± (∏п ≠ 2ф^н)∑ДжА*j + 1 , 2АДж , м= ⋯ =
"="х^†м(∏нф^н) +∑л( ± 1)л - 1(∏п ≠ лф^н)∑ДжА*у + 1 , лАДж , м
Первый член теперь аннулирует левый вакуум, и после подстановки всего искомое среднее становится
⟨ Ψ |ЧАС^0| Ψ⟩=-т∑дж , л , м( ± 1)м - 1( ± 1)л - 1А*у + 1 , лАДж , м⟨ 0 | (∏п ≠ лф^н)⎛⎝∏н′≠ мф^†н′⎞⎠| 0⟩=
⟨ Ψ |ЧАС^0| Ψ⟩=-т∑дж , л , м( ± 1)м - 1( ± 1)л - 1А*у + 1 , лАДж , мдельтал , м
и наконец
⟨ Ψ |ЧАС^0| Ψ⟩=-т∑Дж , мА*j + 1 , мАДж , м
Могут быть некоторые ошибки, которые я пропустил, но это общая идея.
удрв
Мерлин1896