WikiDer > Приближение Борна – Оппенгеймера

Born–Oppenheimer approximation

В квантовая химия и молекулярная физика, то Борн-Оппенгеймер (BO) приближение является наиболее известным математическим приближением в молекулярной динамике. В частности, это предположение, что движение атомные ядра и электроны в молекуле можно рассматривать отдельно, исходя из того, что ядра намного тяжелее электронов. Подход назван в честь Макс Борн и Дж. Роберт Оппенгеймер кто предложил это в 1927 году,[1] в ранний период квантовой механики.

Приближение широко используется в квантовой химии для ускорения вычисления молекулярных волновых функций и других свойств больших молекул. Бывают случаи, когда предположение об отделимости движения больше не выполняется, что делает приближение утратившим силу (говорят, что оно «ломается»), но затем часто используется в качестве отправной точки для более точных методов.

В молекулярном спектроскопия, с использованием BO приближение означает рассмотрение молекулярной энергии как суммы независимых членов, например: Эти члены имеют разный порядок величины, а энергия ядерного спина настолько мала, что ее часто опускают. Электронные энергии состоят из кинетических энергий, межэлектронного отталкивания, межъядерного отталкивания и электрон-ядерного притяжения, которые обычно используются при вычислении электронной структуры молекул.

Пример

В бензол Молекула состоит из 12 ядер и 42 электронов. В Уравнение Шредингера, которое необходимо решить, чтобы получить уровни энергии а волновая функция этой молекулы - уравнение на собственные значения в частных производных в трехмерных координатах ядер и электронов, что дает 3 × 12 + 3 × 42 = 36 ядерных + 126 электронных = 162 переменных для волновой функции. В вычислительная сложность, то есть вычислительная мощность, необходимая для решения уравнения на собственные значения, увеличивается быстрее, чем квадрат числа координат.[2]

При применении приближения БО можно использовать два меньших последовательных шага: для данного положения ядер электронный Уравнение Шредингера решается, при этом ядра рассматриваются как стационарные (не «связанные» с динамикой электронов). Соответствующий собственное значение Тогда проблема состоит только из 126 электронных координат. Затем этот электронный расчет повторяется для других возможных положений ядер, то есть деформаций молекулы. Для бензола это можно сделать с помощью сетки из 36 возможных координат положения ядра. Электронные энергии в этой сетке затем соединяются, чтобы дать поверхность потенциальной энергии для ядер. Затем этот потенциал используется для второго уравнения Шредингера, содержащего только 36 координат ядер.

Итак, взяв наиболее оптимистичную оценку сложности, вместо большого уравнения, требующего по крайней мере гипотетические шаги расчета, серия небольших расчетов, требующих N количество точек сетки для потенциала) и очень небольшой расчет, требующий шаги могут быть выполнены. На практике масштаб проблемы больше, чем , и в вычислительная химия чтобы еще больше уменьшить количество переменных и размеров.

Наклон поверхности потенциальной энергии можно использовать для моделирования молекулярная динамика, используя его для выражения средней силы, действующей на ядра, вызываемой электронами, и тем самым пропуская расчет ядерного уравнения Шредингера.

Подробное описание

Приближение БО учитывает большую разницу между электрон масса и массы атомных ядер и, соответственно, временные масштабы их движения. При таком же количестве кинетической энергии ядра движутся намного медленнее, чем электроны. С математической точки зрения приближение БО состоит из выражения волновая функция () молекулы как продукта электронной волновой функции и ядерной (колебательный, вращающийся) волновая функция. . Это позволяет разделить Гамильтонов оператор на электронные и ядерные термины, где перекрестные члены между электронами и ядрами не учитываются, так что две меньшие и разделенные системы могут быть решены более эффективно.

На первом этапе ядерная кинетическая энергия пренебрегается,[примечание 1] то есть соответствующий оператор Тп вычитается из суммы молекулярный гамильтониан. В оставшемся электронном гамильтониане ЧАСе положения ядер больше не являются переменными, но являются постоянными параметрами (они входят в уравнение «параметрически»). Электронно-ядерные взаимодействия имеют вид нет удалены, т.е. электроны все еще "чувствуют" Кулоновский потенциал ядер, зажатых в определенных положениях в пространстве. (Этот первый шаг приближения БО поэтому часто называют зажатые ядра приближение.)

Электронный Уравнение Шредингера

решается приблизительно[заметка 2] Количество р обозначает все электронные координаты и р для всех ядерных координат. Электронная энергия собственное значение Eе зависит от выбранных позиций р ядер. Варьируя эти позиции р маленькими шагами и многократно решая электронные Уравнение Шредингера, получается Eе как функция р. Это поверхность потенциальной энергии (PES): Eе(р). Поскольку эта процедура пересчета электронных волновых функций как функции бесконечно малой изменяющейся ядерной геометрии напоминает условия для адиабатическая теорема, этот способ получения PES часто называют адиабатическое приближение а сам PES называется адиабатическая поверхность.[заметка 3]

На втором шаге приближения БО кинетическая энергия ядра Тп (содержащие частные производные по компонентам р) вновь вводится, и уравнение Шредингера для движения ядра[примечание 4]

решено. Этот второй шаг приближения БО включает разделение колебательных, поступательных и вращательных движений. Это может быть достигнуто применением Условия Эккарта. Собственное значение E - полная энергия молекулы, включая вклады электронов, ядерных колебаний, а также общего вращения и трансляции молекулы.[требуется разъяснение] В соответствии с Теорема Геллмана – Фейнмана, ядерный потенциал принимается как среднее по электронным конфигурациям суммы электронно-ядерного и межъядерного электрических потенциалов.

Вывод

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

Будет показано, что приближению БО можно доверять, если ППЭ, полученные из решения электронного уравнения Шредингера, хорошо разделены:

.

Начнем с точный нерелятивистский, не зависящий от времени молекулярный гамильтониан:

с

Позиционные векторы электронов и векторов положения ядер относятся к декартовой инерциальная система отсчета. Расстояния между частицами записываются как (расстояние между электронами я и ядро А) и аналогичные определения верны для и .

Мы предполагаем, что молекула находится в однородном (без внешней силы) и изотропном (без внешнего крутящего момента) пространстве. Единственные взаимодействия - это двухчастичные кулоновские взаимодействия между электронами и ядрами. Гамильтониан выражается в атомные единицы, так что мы не видим в этой формуле постоянную Планка, диэлектрическую проницаемость вакуума, заряд или массу электронов. В формулу явно входят только константы: ZА и MА - атомный номер и масса ядра А.

Полезно ввести полный импульс ядра и переписать оператор кинетической энергии ядра следующим образом:

Предположим, у нас есть K собственные электронные функции из , то есть мы решили

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

Например, в молекулярно-орбитальной-линейной-комбинации атомных орбиталей (ЛКАО-МО) приближение, представляет собой молекулярную орбиталь (МО), заданную как линейное расширение атомных орбиталей (АО). АО явно зависит от координат электрона, но ядерные координаты не являются явными в МО. Однако при изменении геометрии, т.е. изменении коэффициенты ЛКАО принимают разные значения и мы видим соответствующие изменения в функциональной форме МО .

Будем предполагать, что параметрическая зависимость непрерывна и дифференцируема, поэтому имеет смысл рассмотреть

которых вообще не будет нуля.

Полная волновая функция расширяется с точки зрения :

с

и где нижний индекс указывает, что интеграция, подразумеваемая обозначение бюстгальтера, находится только над электронными координатами. По определению матрица с общим элементом

диагональный. После умножения на действительную функцию слева и интегрирование по электронным координатам полное уравнение Шредингера

превращается в набор K связанные уравнения на собственные значения, зависящие только от ядерных координат

Вектор-столбец имеет элементы . Матрица диагональна, а ядерная матрица Гамильтона недиагональна; его недиагональный (вибронная муфта) термины дополнительно обсуждаются ниже. Вибронная связь в этом подходе осуществляется посредством ядерной кинетической энергии.

Решение этих связанных уравнений дает приближение для энергии и волновой функции, которое выходит за рамки приближения Борна – Оппенгеймера. К сожалению, с недиагональными членами кинетической энергии обычно трудно справиться. Вот почему часто диабетический применяется преобразование, которое сохраняет часть членов ядерной кинетической энергии на диагонали, удаляет члены кинетической энергии из недиагонали и создает члены связи между адиабатическими ППЭ на недиагонали.

Если мы сможем пренебречь недиагональными элементами, уравнения рассоединятся и резко упростятся. Чтобы показать, когда это пренебрежение оправдано, мы опускаем координаты в обозначениях и пишем, применяя Правило Лейбница для дифференцирования матричные элементы в качестве

Диагональ () элементы матрицы оператора исчезают, потому что мы предполагаем инвариантность обращения времени, поэтому можно выбрать, чтобы он всегда был реальным. Недиагональные матричные элементы удовлетворяют

Матричный элемент в числителе равен

Матричный элемент одноэлектронного оператора, появляющийся в правой части, конечен.

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

И наоборот, если все поверхности хорошо разделены, всеми недиагональными членами можно пренебречь, и, следовательно, всей матрицей фактически равен нулю. Третий член в правой части выражения для матричного элемента ТпДиагональная коррекция Борна – Оппенгеймера) приближенно можно записать в виде матрицы в квадрате и, соответственно, тоже пренебрежимо мала. Только первый (диагональный) член кинетической энергии в этом уравнении сохраняется в случае хорошо разделенных поверхностей, и в результате получается диагональная несвязанная система уравнений движения ядер:

которые являются нормальным вторым этапом описанных выше уравнений БО.

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

Приближение Борна – Оппенгеймера с правильной симметрией

Чтобы включить правильную симметрию в приближение Борна – Оппенгеймера (БО),[1][3] молекулярная система, представленная в виде (зависимых от массы) ядерных координат и образованный двумя самыми низкими поверхностями адиабатической потенциальной энергии BO (PES) и Считается. Чтобы обеспечить справедливость приближения БО, энергия E системы считается достаточно низким, чтобы превращается в замкнутую ППЭ в интересующей области, за исключением единичных бесконечно малых участков, окружающих точки вырождения, образованные и (обозначены как (1, 2) точки вырождения).

Отправной точкой является ядерное адиабатическое БО (матричное) уравнение, записанное в виде[4]

куда вектор-столбец, содержащий неизвестные ядерные волновые функции , - диагональная матрица, содержащая соответствующие адиабатические поверхности потенциальной энергии , м - приведенная масса ядер, E - полная энергия системы, это градиент оператор относительно ядерных координат , и представляет собой матрицу, содержащую элементы векторной неадиабатической связи (NACT):

Здесь являются собственными функциями электронный гамильтониан предполагается сформировать полный Гильбертово пространство в данном регионе в конфигурационное пространство.

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

куда (k = 1, 2) и является (векторным) NACT, ответственным за связь между и .

Затем вводится новая функция:[5]

и производятся соответствующие перестановки:

1. Умножая второе уравнение на я и комбинируя его с первым уравнением, получаем (комплексное) уравнение

2. Последний член в этом уравнении можно удалить по следующим причинам: В тех точках, где классически закрыто, по определению, и в тех точках, где становится классически разрешенным (что происходит в окрестности точек вырождения (1, 2)), это означает, что: , или же . Следовательно, последнее слагаемое действительно пренебрежимо мало в каждой точке интересующей области, и уравнение упрощается и становится

Чтобы это уравнение привело к решению с правильной симметрией, предлагается применить метод возмущений, основанный на упругом потенциале , что совпадает с в асимптотической области.

Уравнение с упругим потенциалом может быть решено простым способом путем подстановки. Таким образом, если является решением этого уравнения, оно представляется как

куда - произвольный контур, а экспоненциальная функция содержит соответствующую симметрию, созданную при движении по .

Функция можно показать как решение (невозмущенного / упругого) уравнения

Имея , полное решение вышеприведенного несвязанного уравнения принимает вид

куда удовлетворяет полученному неоднородному уравнению:

В этом уравнении неоднородность обеспечивает симметрию возмущенной части решения по любому контуру и, следовательно, решения в требуемой области конфигурационного пространства.

Актуальность настоящего подхода была продемонстрирована при исследовании модели с двумя каналами (содержащей один неупругий канал и один реактивный канал), для которой два адиабатических состояния были связаны Ян-Теллер коническое пересечение.[6][7][8] Было получено хорошее совпадение между лечением в одном состоянии с сохранением симметрии и соответствующим лечением в двух состояниях. Это, в частности, относится к реактивным вероятностям между состояниями (см. Таблицу III в [5a] и таблицу III [[5b]), для которых обычное приближение БО привело к ошибочным результатам, тогда как приближение БО, сохраняющее симметрию, произвело точные результаты, как они следовали из решения двух связанных уравнений.

Смотрите также

Примечания

  1. ^ Авторы часто оправдывают этот шаг тем, что «тяжелые ядра движутся медленнее, чем легкие. электроны". Обычно это утверждение имеет смысл, только если импульс п электронов и ядер одного порядка. В таком случае мпме подразумевает п2/(2мп) ≪ п2/(2ме). Легко показать, что для двух тел, вращающихся по круговым орбитам вокруг своего центра масс (независимо от индивидуальных масс), импульсы двух тел равны и противоположны, и что для любого набора частиц в системе центра масс , чистый импульс равен нулю. Учитывая, что система координат центра масс - это лабораторная система координат (где молекула неподвижна), импульс ядер должен быть равен импульсу электронов и противоположен ему. Обоснование махания рукой можно также извлечь из квантовой механики. Соответствующие операторы не содержат массы, и молекулу можно рассматривать как ящик, содержащий электроны и ядра. Поскольку кинетическая энергия равна п2/(2м), следует, что действительно кинетическая энергия ядер в молекуле обычно намного меньше кинетической энергии электронов, причем отношение масс составляет порядка 104).[нужна цитата]
  2. ^ Обычно уравнение Шредингера для молекул не может быть решено точно. Методы аппроксимации включают Метод Хартри-Фока
  3. ^ Предполагается, что в соответствии с адиабатическая теорема, что такое же электронное состояние (например, основное электронное состояние) достигается при небольших изменениях геометрии ядра. Этот метод приведет к разрыву (скачку) в PES, если произойдет переключение электронного состояния.[нужна цитата]
  4. ^ Это уравнение не зависит от времени, и для ядер получены стационарные волновые функции; тем не менее, в этом контексте традиционно используется слово «движение», хотя классически движение подразумевает зависимость от времени.[нужна цитата]

Рекомендации

  1. ^ а б Макс Борн; Дж. Роберт Оппенгеймер (1927). "Zur Quantentheorie der Molekeln" [О квантовой теории молекул]. Annalen der Physik (на немецком). 389 (20): 457–484. Bibcode:1927АнП ... 389..457Б. Дои:10.1002 / andp.19273892002.
  2. ^ Т. Х. Кормен, К. Э. Лейзерсон, Р. Л. Ривест, К. Штайн, Введение в алгоритмы, 3-е изд., MIT Press, Кембридж, Массачусетс, 2009, § 28.2.
  3. ^ Родился М.; Хуанг, К. (1954). «IV». Динамическая теория кристаллических решеток.. Нью-Йорк: Издательство Оксфордского университета.
  4. ^ «Подход Борна-Оппенгеймера: диабатизация и топологическая матрица». За пределами Борна-Оппенгеймера: условия электронной неадиабатической связи и конические пересечения. Хобокен, Нью-Джерси, США: John Wiley & Sons, Inc., 28 марта 2006 г., стр. 26–57. Дои:10.1002 / 0471780081.ch2. ISBN 978-0-471-78008-3.
  5. ^ Баер, Майкл; Энглман, Роберт (1997). «Модифицированное уравнение Борна-Оппенгеймера: приложение к коническим пересечениям и другим типам особенностей». Письма по химической физике. Elsevier BV. 265 (1–2): 105–108. Bibcode:1997CPL ... 265..105B. Дои:10.1016 / с0009-2614 (96) 01411-х. ISSN 0009-2614.
  6. ^ Баер, Рой; Чаруц, Дэвид М .; Кослофф, Ронни; Баер, Майкл (22 ноября 1996 г.). «Исследование эффектов конического пересечения на процессы рассеяния: обоснованность адиабатических приближений одной поверхности в квази-модели Яна – Теллера». Журнал химической физики. Издательство AIP. 105 (20): 9141–9152. Bibcode:1996ЖЧФ.105.9141Б. Дои:10.1063/1.472748. ISSN 0021-9606.
  7. ^ Адхикари, Сатраджит; Биллинг, Герт Д. (1999). «Эффекты конического пересечения и адиабатические одноповерхностные приближения процессов рассеяния: подход с временными волновыми пакетами». Журнал химической физики. Издательство AIP. 111 (1): 40–47. Bibcode:1999ЖЧФ.111 ... 40А. Дои:10.1063/1.479360. ISSN 0021-9606.
  8. ^ Чаруц, Дэвид М .; Баер, Рой; Баер, Майкл (1997). «Изучение эффектов вырожденной вибронной связи на процессы рассеяния: влияет ли вырожденная вибронная связь на резонансы?». Письма по химической физике. Elsevier BV. 265 (6): 629–637. Bibcode:1997CPL ... 265..629C. Дои:10.1016 / с0009-2614 (96) 01494-7. ISSN 0009-2614.

внешняя ссылка

Ресурсы, связанные с приближением Борна – Оппенгеймера: