Разделение Борна-Оппенгеймера в скобках Дирака

Большинство выводов приближения Борна-Оппенгеймера, которые я видел, основаны на волновых функциях. Чтобы лучше понять это, я пытался написать вывод, используя нотацию Дирака, но я застрял. Я собираюсь опубликовать то, что я сделал до сих пор, чтобы вы, ребята, могли мне помочь.

Гамильтониан молекулы можно записать в виде суммы двух частей: ЧАС моль "=" ЧАС Эль + ЧАС нук , где ЧАС нук является гамильтонианом ядер самих по себе, и ЧАС Эль – гамильтониан электронов, взаимодействующих с ядрами:

ЧАС Эль "=" Т Эль + В эль-эль + В эль-нук
ЧАС нук "=" Т нук + В нук-нук

Мы хотим найти энергетические уровни молекулы. То есть мы хотим решить ЧАС моль | Е "=" Е | Е , где Е и | Е являются собственной энергией и соответствующим собственным набором молекулы.

Пространство состояний молекулы можно разделить на электронную и ядерную части: С моль "=" С Эль С нук . Позволять | р е С нук – позиционный базис ядер, где р обозначает координаты всех ядер.

Определять ЧАС Эль ( р ) "=" р | ЧАС Эль | р , который является оператором в С Эль . Позволять Е а ( р ) и | Е а ( р ) — собственные значения и соответствующие собственные наборы ЧАС Эль ( р ) в С Эль , так что ЧАС Эль ( р ) | Е а ( р ) "=" Е а ( р ) | Е а ( р ) . Для каждого р , кеты | Е а ( р ) сделать основу для С Эль .

Набор кетов | Е а ( р ) | р е С моль для всех р и а затем сделать основу для С моль . Используя этот базис, можно записать состояние молекулы:

| ψ "=" а х а ( р ) | Е а ( р ) | р д р

где х а ( р ) является амплитудной функцией.

Обратите внимание, что:

ЧАС Эль | Е а ( р ) | р "=" [ ЧАС Эль ( р ) | Е а ( р ) ] | р "=" Е а ( р ) | Е а ( р ) | р

Таким образом, молекулярная собственная проблема ЧАС моль | ψ "=" Е | ψ можно написать:

а ( Е а + Т нук + В нук-нук Е ) х а ( р ) | Е а ( р ) | р д р "=" 0

Умножение на Е а ( р ) | слева:

( Е а + Т нук + В нук-нук Е ) х а ( р ) | р д р "=" 0

Наконец, мы определяем кет | х а е С нук такой, что х а ( р ) "=" р | х а :

| х а "=" х а ( р ) | р д р

Тогда мы можем написать:

( Т нук + В нук-нук + Е а Е ) | х а "=" 0

ЗДЕСЬ ЗАКАНЧИВАЕТСЯ МОЯ ВЫВОД.

Должно быть, я где-то сделал что-то не так, потому что окончательное уравнение, которое я получаю, является, насколько я могу судить, приближением Борна-Оппенгеймера , но здесь я получаю его как точное уравнение. Что я сделал не так?

Кроме того, если кто-нибудь знает какой-либо справочник, учебник или документ, посвященный приближению Борна-Оппенгеймера в обозначениях Дирака, пожалуйста, опубликуйте его.

Ответы (1)

Ваша ошибка состоит в том, что вы умножили уравнение

а д р [ ( Е а + Т нук + В нук-нук Е ) х а ( р ) | Е а ( р ) | р ] "=" 0
слева от Е а ( р ) | . Левая часть представляет собой вектор в С моль , и хотя вы можете спроецировать его на С нук путем умножения его слева на бюстгальтер из С Эль , переменные р и а немы в том смысле, что не имеют значения вне интеграла/суммирования.

Таким образом, вы должны умножить слева на (немного, но важно другое) бюстгальтер Е а ( р ) | :

Е а ( р ) | а д р [ ( Е а + Т нук + В нук-нук Е ) х а ( р ) | Е а ( р ) | р ] "=" 0.
Эквивалентно,
( Е а + Т нук + В нук-нук Е ) а д р [ Е а ( р ) | Е а ( р ) ] х а ( р ) | р "=" 0.

(Заметим, что здесь зависимость Е а на р устраняется, делая его оператором Е а "=" д р Е а ( р ) | р р | на С нук .)

Теперь, если векторы | Е а ( р ) и | Е а ( р ) ортогональны, т.е. если

Е а ( р ) | Е а ( р ) "=" дельта а а ,
тогда вы можете исключить как суммирование, так и интеграл, чтобы получить ядерную ТИСЭ Борна-Оппенгеймера. Однако на самом деле для этого нет никаких причин, поскольку они являются электронными собственными состояниями для разных гамильтонианов. Именно этими недиагональными членами следует пренебречь в приближении Борна-Оппенгеймера; предположение, что они равны нулю, эквивалентно постулированию структуры тензорного произведения всех собственных состояний, которая лежит в основе BOA.

На практике, конечно, они близки к нулю, поскольку молекула не сильно растягивается и не изгибается, а структура собственных электронных состояний не слишком сильно изменяется. Суть в том, чтобы помнить, что они приблизительно равны нулю.

По сравнению с моим намного лучше :-)
Вопрос: как вы извлекаете Е а ( р ) от ЧАС е л ? У вас есть (учитывая только электронную часть)
( Е а ( р ) | р | ) ЧАС е л ( а д р х а ( р ) | Е а ( р ) | р ) ,
и иметь Е а ( р ) вам понадобится
ЧАС е л ( р ) р | ЧАС е л | р ,
пока у тебя есть
р | ЧАС е л | р ,
с р р в общем. Мы также предполагаем, что
р | ЧАС е л | р "=" дельта ( р р ) ЧАС е л ( р ) ,
(что на самом деле имело бы смысл)?
@взгляд Да, это почти все.