Ожидаемое значение при втором квантовании

Я застрял в вычислении простого ожидаемого значения для оператора, которое выражается в секундном квантовании. Я знаю результат, но не могу его доказать.

Допустим, у меня есть одночастичная волновая функция. | ф н данный | ф н "=" Дж "=" 1 К | α Дж А Дж , н , где К - количество орбиталей/сайтов в системе и А Дж , н - амплитуды вероятности орбиталей | α Дж . Индекс н помечает частицы в системе, из которых мы имеем Н .

Орбитали ортонормированы, т.е.

Дж А Дж , м * А Дж , н "=" дельта м , н .

Давайте игнорировать любые спиновые степени свободы. Волновая функция многих частиц теперь определяется выражением

| Ψ "=" ( н "=" 1 Н Дж "=" 1 К с ^ Дж А Дж , н ) | вакуум ,
где с ^ Дж обычный оператор создания на сайте Дж .

Теперь я хочу рассчитать ожидаемое значение

Ψ | ЧАС прыгать | Ψ
с
ЧАС прыгать "=" т Дж "=" 1 К с ^ Дж + 1 с ^ Дж + час . с .
обычный скачкообразный гамильтониан. У меня есть сильное ощущение (и один пример расчета подтвердил это), что результат просто
Ψ | ЧАС прыгать | Ψ "=" т н "=" 1 Н Дж "=" 1 К ( А Дж , н * А Дж + 1 , н + А Дж + 1 , н * А Дж , н )

Я думаю, что этот результат тривиально связан с правилами Слейтера-Кондона, но я не вижу связи. Кроме того, я не могу явно вычислить ожидаемое значение, которое содержит сумму и произведения операторов создания/уничтожения.

Как лучше всего подтвердить свой результат?

Подсказка: используйте CCR, связанные с с ^ Дж -с. Разве ваш примерный расчет уже не предполагает что-то подобное?
Что ж, готов поспорить, что в какой-то момент мне придется использовать коммутационные соотношения. Но я действительно не вижу, как я могу прийти к общему результату. Так как скалярное произведение является линейным, достаточно просто рассмотреть один член из ЧАС прыгать , так что, возможно, можно установить wlog Дж "=" 1 . Но тогда у меня еще много с ^ из | Ψ (детерминанты Слейтера). Как мне поступить?

Ответы (1)

Для удобства в будущем обозначим

ф ^ н "=" Дж А Дж , н с ^ Дж
х ^ н "=" Дж А Дж , н с ^ Дж + 1
Среднее значение, которое вы хотите рассчитать, читает тогда
Ψ | ЧАС ^ 0 | Ψ "=" т 0 | н ф ^ н ( Дж с ^ Дж + 1 с ^ Дж ) н ф ^ н | 0
Начиная с обычных CCR-ов,
[ с ^ Дж , с ^ к ] "=" дельта Дж к , [ с ^ Дж , с ^ к ] "=" [ с ^ Дж , с ^ к ] "=" 0
получать
[ Дж с ^ Дж + 1 с ^ Дж , ф ^ н ] "=" х ^ н
и
ф ^ н х ^ м "=" ± х ^ м ф ^ н
ф ^ н х ^ м "=" ± х ^ м ф ^ н + Дж А Дж + 1 , н * А Дж , м
Теперь в среднем вам нужно рассчитать, использовать вышеизложенное для последовательного перемещения ( Дж с ^ Дж + 1 с ^ Дж ) мимо орбитальных операторов справа:
( Дж с ^ Дж + 1 с ^ Дж ) н ф ^ н "=" ф ^ 1 ( Дж с ^ Дж + 1 с ^ Дж ) н > 1 ф ^ н + х ^ 1 н 1 ф ^ н "="
"=" м "=" 1 м "=" 2 ф ^ м ( Дж с ^ Дж + 1 с ^ Дж ) н > 2 ф ^ н + х ^ 1 н 1 ф ^ н + ф ^ 1 х ^ 2 н > 2 ф ^ н "="
"=" м "=" 1 м "=" 2 ф ^ м ( Дж с ^ Дж + 1 с ^ Дж ) н > 2 ф ^ н + х ^ 1 н 1 ф ^ н ± х ^ 2 н 2 ф ^ н "=" "="
"=" м ф ^ м ( Дж с ^ Дж + 1 с ^ Дж ) + м ( ± 1 ) м 1 х ^ м н м ф ^ н
Первый член аннулирует правый вакуум, поэтому остается только сумма. Теперь введите орбитальные операторы слева и переверните каждый х ^ м мимо них:
( н ф ^ н ) х ^ м "=" ± ( н > 1 ф ^ н ) х ^ м ф ^ 1 + ( н 1 ф ^ н ) Дж А Дж + 1 , 1 * А Дж , м "="
( н > 1 ф ^ н ) х ^ м л "=" 1 л "=" 2 ф ^ л + ( н 1 ф ^ н ) Дж А Дж + 1 , 1 * А Дж , м ± ( н 2 ф ^ н ) Дж А Дж + 1 , 2 * А Дж , м "=" "="
"=" х ^ м ( н ф ^ н ) + л ( ± 1 ) л 1 ( н л ф ^ н ) Дж А Дж + 1 , л * А Дж , м
Первый член теперь аннулирует левый вакуум, и после подстановки всего искомое среднее становится
Ψ | ЧАС ^ 0 | Ψ "=" т Дж , л , м ( ± 1 ) м 1 ( ± 1 ) л 1 А Дж + 1 , л * А Дж , м 0 | ( н л ф ^ н ) ( н м ф ^ н ) | 0 "="
Ψ | ЧАС ^ 0 | Ψ "=" т Дж , л , м ( ± 1 ) м 1 ( ± 1 ) л 1 А Дж + 1 , л * А Дж , м дельта л , м
и наконец
Ψ | ЧАС ^ 0 | Ψ "=" т Дж , м А Дж + 1 , м * А Дж , м
Могут быть некоторые ошибки, которые я пропустил, но это общая идея.