Молекулярный гамильтониан

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

Элементарные части молекулы - это ядра, характеризуемые их атомными номерами , Z , и электроны, которые имеют отрицательный элементарный заряд , - e . Их взаимодействие дает заряд ядра Z  +  q , где q  = - eN , где N равно количеству электронов. Электроны и ядра в очень хорошем приближении являются точечными зарядами и точечными массами. Молекулярный гамильтониан представляет собой сумму нескольких членов: его главными членами являются кинетические энергии электронов и кулоновские (электростатические) взаимодействия.между двумя видами заряженных частиц. Гамильтониан, содержащий только кинетические энергии электронов и ядер и кулоновские взаимодействия между ними, известен как кулоновский гамильтониан . В нем отсутствует ряд мелких членов, большая часть которых связана с электронным и ядерным спином .

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

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

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

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

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

Кулоновский гамильтониан [ править ]

Алгебраическая форма многих наблюдаемых, т. Е. Эрмитовых операторов, представляющих наблюдаемые величины, получается с помощью следующих правил квантования :

• Запишите классическую форму наблюдаемой в форме Гамильтона (как функцию импульсов p и положений q ). Оба вектора выражаются относительно произвольной инерциальной системы отсчета , обычно называемой лабораторной системой отсчета или системой отсчета с фиксированным пространством .
• Замените p на и интерпретируйте q как мультипликативный оператор. Вот это набли оператор, векторный оператор , состоящий из первых производных. Известные коммутационные соотношения для операторов p и q непосредственно следуют из правил дифференцирования.${\ displaystyle -i \ hbar {\ boldsymbol {\ nabla}}}$${\ displaystyle {\ boldsymbol {\ nabla}}}$

Классически электроны и ядра в молекуле имеют кинетическую энергию формы p ² / (2 m) и взаимодействуют посредством кулоновских взаимодействий , которые обратно пропорциональны расстоянию r _ij между частицей i и j .

${\displaystyle r_{ij}\equiv |\mathbf {r} _{i}-\mathbf {r} _{j}|={\sqrt {(\mathbf {r} _{i}-\mathbf {r} _{j})\cdot (\mathbf {r} _{i}-\mathbf {r} _{j})}}={\sqrt {(x_{i}-x_{j})^{2}+(y_{i}-y_{j})^{2}+(z_{i}-z_{j})^{2}}}.}$

В этом выражении r _i обозначает вектор координат любой частицы (электрона или ядра), но с этого момента мы сохраним заглавную букву R для обозначения ядерной координаты и строчную букву r для электронов системы. Координаты могут быть выражены относительно любой декартовой системы отсчета с центром в любом месте пространства, потому что расстояние, являющееся внутренним продуктом, инвариантно при вращении кадра и, будучи нормой вектора разности, расстояние инвариантно при перемещении рама тоже.

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

1. Операторы кинетической энергии для каждого ядра в системе;
2. Операторы кинетической энергии для каждого электрона в системе;
3. Потенциальная энергия между электронами и ядрами - полное электрон-ядерное кулоновское притяжение в системе;
4. Потенциальная энергия, возникающая в результате кулоновского электрон-электронного отталкивания
5. Потенциальная энергия, возникающая в результате отталкивания кулоновских ядер и ядер, также известная как энергия ядерного отталкивания. См. Электрический потенциал для получения более подробной информации.

1. ${\displaystyle {\hat {T}}_{n}=-\sum _{i}{\frac {\hbar ^{2}}{2M_{i}}}\nabla _{\mathbf {R} _{i}}^{2}}$
2. ${\displaystyle {\hat {T}}_{e}=-\sum _{i}{\frac {\hbar ^{2}}{2m_{e}}}\nabla _{\mathbf {r} _{i}}^{2}}$
3. ${\displaystyle {\hat {U}}_{en}=-\sum _{i}\sum _{j}{\frac {Z_{i}e^{2}}{4\pi \epsilon _{0}\left|\mathbf {R} _{i}-\mathbf {r} _{j}\right|}}}$
4. ${\displaystyle {\hat {U}}_{ee}={1 \over 2}\sum _{i}\sum _{j\neq i}{\frac {e^{2}}{4\pi \epsilon _{0}\left|\mathbf {r} _{i}-\mathbf {r} _{j}\right|}}=\sum _{i}\sum _{j>i}{\frac {e^{2}}{4\pi \epsilon _{0}\left|\mathbf {r} _{i}-\mathbf {r} _{j}\right|}}}$
5. ${\displaystyle {\hat {U}}_{nn}={1 \over 2}\sum _{i}\sum _{j\neq i}{\frac {Z_{i}Z_{j}e^{2}}{4\pi \epsilon _{0}\left|\mathbf {R} _{i}-\mathbf {R} _{j}\right|}}=\sum _{i}\sum _{j>i}{\frac {Z_{i}Z_{j}e^{2}}{4\pi \epsilon _{0}\left|\mathbf {R} _{i}-\mathbf {R} _{j}\right|}}.}$

Здесь M _i - масса ядра i , Z _i - атомный номер ядра i , а m _e - масса электрона. Оператор Лапласа частицы я это: . Поскольку оператор кинетической энергии является внутренним произведением, он инвариантен относительно вращения декартовой системы отсчета, относительно которой выражаются x _i , y _i и z _i .${\displaystyle \nabla _{\mathbf {r} _{i}}^{2}\equiv {\boldsymbol {\nabla }}_{\mathbf {r} _{i}}\cdot {\boldsymbol {\nabla }}_{\mathbf {r} _{i}}={\frac {\partial ^{2}}{\partial x_{i}^{2}}}+{\frac {\partial ^{2}}{\partial y_{i}^{2}}}+{\frac {\partial ^{2}}{\partial z_{i}^{2}}}}$

Небольшие условия [ править ]

В 1920-х годах многие спектроскопические данные ясно показали, что в кулоновском гамильтониане отсутствуют некоторые члены. Эти члены, особенно для молекул, содержащих более тяжелые атомы, хотя и намного меньше, чем кинетическая и кулоновская энергии, не пренебрежимо малы. Эти спектроскопические наблюдения привели к введению новой степени свободы для электронов и ядер, а именно спина . Эта эмпирическая концепция была теоретически обоснована Полем Дираком, когда он ввел релятивистски правильную ( лоренц-ковариантную ) форму одночастичного уравнения Шредингера. Уравнение Дирака предсказывает, что спин и пространственное движение частицы взаимодействуют посредством спин-орбитальной связи . По аналогии спин-другая-орбитальная связьбыл представлен. Тот факт, что спин частицы обладает некоторыми характеристиками магнитного диполя, привел к спин-спиновой связи . Другими членами, не имеющими классического аналога, являются член ферми-контакта (взаимодействие электронной плотности на ядре конечного размера с ядром) и ядерная квадрупольная связь (взаимодействие ядерного квадруполя с градиентом электрического поля, обусловленного электронами). Наконец, необходимо упомянуть условие нарушения четности, предсказываемое Стандартной моделью . Хотя это чрезвычайно слабое взаимодействие, оно привлекло изрядное внимание в научной литературе, поскольку дает разные энергии энантиомерам вхиральные молекулы .

В оставшейся части этой статьи спиновые члены будут игнорироваться и будет рассмотрено решение уравнения на собственные значения (не зависящее от времени) уравнения Шредингера кулоновского гамильтониана.

Уравнение Шредингера кулоновского гамильтониана [ править ]

Кулоновский гамильтониан имеет непрерывный спектр из-за движения центра масс (ЦМ) молекулы в однородном пространстве. В классической механике легко отделить СОМ движение от системы точечных масс. Обычно движение COM не связано с другими движениями. COM движется равномерно (то есть с постоянной скоростью) в пространстве, как если бы это была точечная частица с массой, равной сумме M _tot масс всех частиц.

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

Вводя координатный вектор X центра масс как три степени свободы системы и исключая координатный вектор одной (произвольной) частицы, так что количество степеней свободы остается неизменным, можно получить линейным преобразование нового набора координат t _i . Эти координаты представляют собой линейные комбинации старых координат всех частиц (ядер и электронов). Применяя цепное правило, можно показать, что

${\displaystyle H=-{\frac {\hbar ^{2}}{2M_{\textrm {tot}}}}\nabla _{\mathbf {X} }^{2}+H'\quad {\text{with }}\quad H'=-{\frac {\hbar ^{2}}{2}}\sum _{i=1}^{N_{\textrm {tot}}-1}{\frac {1}{m_{i}}}\nabla _{i}^{2}+{\frac {\hbar ^{2}}{2M_{\textrm {tot}}}}\sum _{i,j=1}^{N_{\textrm {tot}}-1}\nabla _{i}\cdot \nabla _{j}+V(\mathbf {t} ).}$

Первый член есть кинетическая энергия движения COM, который может быть обработан по отдельности , так как не зависит от X . Как только что сказано, его собственные состояния - плоские волны. Потенциал V ( t ) состоит из кулоновских членов, выраженных в новых координатах. Первый член имеет обычный вид оператора кинетической энергии. Второй член известен как член массовой поляризации . Можно показать, что трансляционно-инвариантный гамильтониан самосопряжен и ограничен снизу. То есть его наименьшее собственное значение действительно и конечно. Хотя обязательно инвариантен относительно перестановок одинаковых частиц (поскольку${\displaystyle H}$${\displaystyle H'}$${\displaystyle H'}$${\displaystyle H'}$${\displaystyle H'}$${\displaystyle H}$ и кинетическая энергия СОМ инвариантны), его инвариантность не проявляется.

Существует не так много реальных молекулярных приложений ; см., однако, основополагающую работу ^[1] о молекуле водорода для более раннего применения. В подавляющем большинстве вычислений молекулярных волновых функций электронная задача решается с помощью гамильтониана зажатого ядра, возникающего на первом этапе приближения Борна – Оппенгеймера .${\displaystyle H'}$

См. Ссылку. ^[2] для подробного обсуждения математических свойств кулоновского гамильтониана. Также в этой статье обсуждается, можно ли априори прийти к концепции молекулы (как стабильной системы электронов и ядер с четко определенной геометрией) только на основе свойств кулоновского гамильтониана.

Гамильтониан зажатого ядра [ править ]

Гамильтониан зажатого ядра описывает энергию электронов в электростатическом поле ядер, где ядра считаются стационарными по отношению к инерциальной системе отсчета. Электронный гамильтониан имеет вид

${\displaystyle {\hat {H}}_{\mathrm {el} }={\hat {T}}_{e}+{\hat {U}}_{en}+{\hat {U}}_{ee}+{\hat {U}}_{nn}.}$

Координаты электронов и ядер выражаются относительно системы отсчета, которая движется вместе с ядрами, так что ядра находятся в состоянии покоя относительно этой системы отсчета. Рама остается параллельной рамке с фиксированным пространством. Это инерциальная система отсчета, поскольку предполагается, что ядра не ускоряются внешними силами или моментами. Начало рамки произвольно, обычно она располагается на центральном ядре или в центре масс ядра. Иногда говорят, что ядра «покоятся в фиксированной пространственной системе отсчета». Это утверждение подразумевает, что ядра рассматриваются как классические частицы, потому что квантово-механическая частица не может находиться в состоянии покоя. (Это означало бы, что он имел одновременно нулевой импульс и четко определенное положение, что противоречит принципу неопределенности Гейзенберга).

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

Как объясняется в статье о приближении Борна – Оппенгеймера , достаточное количество решений уравнения Шредингера приводит к поверхности потенциальной энергии (ППЭ) . Предполагается, что функциональная зависимость V от его координат такова, что${\displaystyle H_{\text{el}}}$${\displaystyle V(\mathbf {R} _{1},\mathbf {R} _{2},\ldots ,\mathbf {R} _{N})}$

${\displaystyle V(\mathbf {R} _{1},\mathbf {R} _{2},\ldots ,\mathbf {R} _{N})=V(\mathbf {R} '_{1},\mathbf {R} '_{2},\ldots ,\mathbf {R} '_{N})}$

за

${\displaystyle \mathbf {R} '_{i}=\mathbf {R} _{i}+\mathbf {t} \;\;{\text{(translation) and}}\;\;\mathbf {R} '_{i}=\mathbf {R} _{i}+{\frac {\Delta \phi }{|\mathbf {s} |}}\;(\mathbf {s} \times \mathbf {R} _{i})\;\;{\text{(infinitesimal rotation)}},}$

где t и s - произвольные векторы, а Δφ - бесконечно малый угол, Δφ >> Δφ ² . Это условие инвариантности для PES автоматически выполняется, когда PES выражается в терминах разностей и углов между R _i , что обычно имеет место.

Гамильтониан гармонического движения ядра [ править ]

В оставшейся части статьи мы предполагаем, что молекула полужесткая . На втором этапе приближения БО снова вводится ядерная кинетическая энергия T _n и уравнение Шредингера с гамильтонианом

${\displaystyle {\hat {H}}_{\mathrm {nuc} }=-{\frac {\hbar ^{2}}{2}}\sum _{i=1}^{N}\sum _{\alpha =1}^{3}{\frac {1}{M_{i}}}{\frac {\partial ^{2}}{\partial R_{i\alpha }^{2}}}+V(\mathbf {R} _{1},\ldots ,\mathbf {R} _{N})}$

Считается. В ее решении хотелось бы распознать: движение центра масс ядра (3 степени свободы), общее вращение молекулы (3 степени свободы) и ядерные колебания. В общем, это невозможно с данной ядерной кинетической энергией, потому что она не отделяет явно 6 внешних степеней свободы (общий перенос и вращение) от 3 N  - 6 внутренних степеней свободы. Фактически, оператор кинетической энергии здесь определен относительно фиксированной в пространстве (SF) системы отсчета. Если бы мы переместили начало системы отсчета SF в центр масс ядра, то, применив правило цепочки , появились бы члены поляризации ядерной массы. Принято полностью игнорировать эти термины, и мы будем следовать этому обычаю.

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

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

${\displaystyle {\boldsymbol {\rho }}_{i}\equiv {\sqrt {M_{i}}}(\mathbf {R} _{i}-\mathbf {R} _{i}^{0})}$.

С

${\displaystyle {\frac {\partial }{\partial \rho _{i\alpha }}}={\frac {\partial }{{\sqrt {M_{i}}}(\partial R_{i\alpha }-\partial R_{i\alpha }^{0})}}={\frac {1}{\sqrt {M_{i}}}}{\frac {\partial }{\partial R_{i\alpha }}},}$

оператор кинетической энергии принимает вид

${\displaystyle T=-{\frac {\hbar ^{2}}{2}}\sum _{i=1}^{N}\sum _{\alpha =1}^{3}{\frac {\partial ^{2}}{\partial \rho _{i\alpha }^{2}}}.}$

Если мы сделаем разложение V по Тейлору вокруг равновесной геометрии,

${\displaystyle V=V_{0}+\sum _{i=1}^{N}\sum _{\alpha =1}^{3}{\Big (}{\frac {\partial V}{\partial \rho _{i\alpha }}}{\Big )}_{0}\;\rho _{i\alpha }+{\frac {1}{2}}\sum _{i,j=1}^{N}\sum _{\alpha ,\beta =1}^{3}{\Big (}{\frac {\partial ^{2}V}{\partial \rho _{i\alpha }\partial \rho _{j\beta }}}{\Big )}_{0}\;\rho _{i\alpha }\rho _{j\beta }+\cdots ,}$

и усекая после трех членов (так называемое гармоническое приближение), мы можем описать V только третьим членом. Член V ₀ может поглощаться энергией (дает новый ноль энергии). Второй член исчезает из-за условия равновесия. Оставшийся срок содержит матрицу Гессе F из V , который является симметричным и может быть диагонализуются с ортогональным- N  × 3 Н матрицей с постоянными элементами:

${\displaystyle \mathbf {Q} \mathbf {F} \mathbf {Q} ^{\mathrm {T} }={\boldsymbol {\Phi }}\quad \mathrm {with} \quad {\boldsymbol {\Phi }}=\operatorname {diag} (f_{1},\dots ,f_{3N-6},0,\ldots ,0).}$

Из инвариантности V относительно вращения и сдвига можно показать, что шесть собственных векторов F (последние шесть строк Q ) имеют нулевое собственное значение (являются модами с нулевой частотой). Они охватывают внешнее пространство . Первые 3 N  - 6 строк Q являются - для молекул в их основном состоянии - собственными векторами с ненулевым собственным значением; они являются внутренними координатами и образуют ортонормированный базис для (3 N  - 6) -мерного подпространства ядерного конфигурационного пространства R ^{3 N} , внутреннего пространства. Собственные векторы с нулевой частотой ортогональны собственным векторам с ненулевой частотой. Можно показать, что эти ортогональности на самом деле являются условиями Эккарта . Кинетическая энергия, выраженная во внутренних координатах, является внутренней (колебательной) кинетической энергией.

С введением нормальных координат

${\displaystyle q_{t}\equiv \sum _{i=1}^{N}\sum _{\alpha =1}^{3}\;Q_{t,i\alpha }\rho _{i\alpha },}$

колебательная (внутренняя) часть гамильтониана движения ядра принимает вид в гармоническом приближении

${\displaystyle {\hat {H}}_{\mathrm {nuc} }\approx {\frac {1}{2}}\sum _{t=1}^{3N-6}\left[-\hbar ^{2}{\frac {\partial ^{2}}{\partial q_{t}^{2}}}+f_{t}q_{t}^{2}\right].}$

Соответствующее уравнение Шредингера легко решается, оно разлагается на 3 N  - 6 уравнений для одномерных гармонических осцилляторов . Основное усилие в этом приближенном решении уравнения Шредингера движения ядра - это вычисление гессиана F для V и его диагонализация.

Это приближение к проблеме движения ядра, описанное в 3 N декартовых координатах, взвешенных по массе, стало стандартом в квантовой химии с тех пор (1980–1990-е годы), когда стали доступны алгоритмы для точных вычислений гессиана F. Помимо гармонического приближения, он имеет еще один недостаток, заключающийся в том, что внешние (вращательные и поступательные) движения молекулы не учитываются. Они учитываются в колебательном гамильтониане, который иногда называют гамильтонианом Ватсона .

Гамильтониан ядерного движения Ватсона [ править ]

Чтобы получить гамильтониан для внешних (поступательных и вращательных) движений, связанных с внутренними (колебательными) движениями, обычно возвращаются в этот момент к классической механике и формулируют классическую кинетическую энергию, соответствующую этим движениям ядер. Классически легко отделить поступательное движение центра масс от других движений. Однако отделить вращательное движение от колебательного сложнее и не вполне возможно. Это вращательно-колебательное разделение было впервые достигнуто Эккартом ^[3] в 1935 году путем наложения того, что сейчас известно как условия Эккарта . Поскольку проблема описана в системе отсчета («Эккарт»), которая вращается вместе с молекулой и, следовательно, является неинерциальной системой отсчета., энергии, связанные с фиктивными силами : центробежная сила и сила Кориолиса входят в кинетическую энергию.

В общем, классическая кинетическая энергия T определяет метрический тензор g = ( g _ij ), связанный с криволинейными координатами s = ( s _i ) через

${\displaystyle 2T=\sum _{ij}g_{ij}{\dot {s}}_{i}{\dot {s}}_{j}.}$

Шаг квантования - это преобразование этой классической кинетической энергии в квантово-механический оператор. Обычно вслед за Подольским ^[4] оператор Лапласа – Бельтрами записывается в тех же (обобщенных, криволинейных) координатах s, что и в классической форме. Уравнение для этого оператора требует обращения к метрическому тензору g и его определителю. Умножение оператора Лапласа – Бельтрами на дает требуемый квантовомеханический оператор кинетической энергии. Когда мы применяем этот рецепт к декартовым координатам, которые имеют единичную метрику, получается та же кинетическая энергия, что и при применении правил квантования .${\displaystyle -\hbar ^{2}}$

Гамильтониан ядерного движения был получен Уилсоном и Ховардом в 1936 году ^[5], которые следовали этой процедуре, и дополнительно уточнен Дарлингом и Деннисоном в 1940 году. ^[6] Он оставался стандартом до 1968 года, когда Уотсон ^[7] смог упростить его. это радикально, коммутируя через производные определитель метрического тензора. Мы дадим вращательно-колебательный гамильтониан, полученный Ватсоном, который часто называют гамильтонианом Ватсона . Прежде чем мы это сделаем, мы должны упомянуть, что вывод этого гамильтониана также возможен, исходя из оператора Лапласа в декартовой форме, применения преобразований координат и использования цепного правила . ^[8]Гамильтониан Ватсона, описывающий все движения ядер N , имеет вид

${\displaystyle {\hat {H}}=-{\frac {\hbar ^{2}}{2M_{\mathrm {tot} }}}\sum _{\alpha =1}^{3}{\frac {\partial ^{2}}{\partial X_{\alpha }^{2}}}+{\frac {1}{2}}\sum _{\alpha ,\beta =1}^{3}\mu _{\alpha \beta }({\mathcal {P}}_{\alpha }-\Pi _{\alpha })({\mathcal {P}}_{\beta }-\Pi _{\beta })+U-{\frac {\hbar ^{2}}{2}}\sum _{s=1}^{3N-6}{\frac {\partial ^{2}}{\partial q_{s}^{2}}}+V.}$

Первый член - это член центра масс.

${\displaystyle \mathbf {X} \equiv {\frac {1}{M_{\mathrm {tot} }}}\sum _{i=1}^{N}M_{i}\mathbf {R} _{i}\quad \mathrm {with} \quad M_{\mathrm {tot} }\equiv \sum _{i=1}^{N}M_{i}.}$

Второй член - это вращательный член, родственный кинетической энергии жесткого ротора . Здесь представлена ​​α-составляющая оператора углового момента жесткого ротора, закрепленного на теле, ее выражение в терминах углов Эйлера см. В этой статье . Оператор является компонентом оператора, известного как оператор колебательного углового момента (хотя он не удовлетворяет соотношениям коммутации углового момента),${\displaystyle {\mathcal {P}}_{\alpha }}$${\displaystyle \Pi _{\alpha }\,}$

${\displaystyle \Pi _{\alpha }=-i\hbar \sum _{s,t=1}^{3N-6}\zeta _{st}^{\alpha }\;q_{s}{\frac {\partial }{\partial q_{t}}}}$

с константой связи Кориолиса :

${\displaystyle \zeta _{st}^{\alpha }=\sum _{i=1}^{N}\sum _{\beta ,\gamma =1}^{3}\epsilon _{\alpha \beta \gamma }Q_{s,i\beta }\,Q_{t,i\gamma }\;\;\mathrm {and} \quad \alpha =1,2,3.}$

Здесь ε _αβγ - символ Леви-Чивиты . Члены, квадратичные по отношению к центробежным членам, билинейные по и являются членами Кориолиса. Величины Q _{s, iγ} являются компонентами введенных выше нормальных координат. В качестве альтернативы нормальные координаты могут быть получены с помощью метода GF Вильсона . Симметричная матрица 3 × 3 называется эффективным тензором обратной инерции . Если все Q _s были равны нулю (жесткая молекула) рама Эккарт будет совпадать с основной рамой осей (см жесткий ротор ) и${\displaystyle {\mathcal {P}}_{\alpha }}$${\displaystyle {\mathcal {P}}_{\alpha }}$${\displaystyle \Pi _{\beta }\,}$${\displaystyle {\boldsymbol {\mu }}}$${\displaystyle {\boldsymbol {\mu }}}$будет диагональным, с равновесными взаимными моментами инерции на диагонали. Если все Q _s будет равен нуль, только кинетическая энергия перевода и вращения твердого выживет.

Подобный потенциалу член U - это член Ватсона :

${\displaystyle U=-{\frac {1}{8}}\sum _{\alpha =1}^{3}\mu _{\alpha \alpha }}$

пропорциональна следу эффективного тензора обратной инерции.

Четвертый член в гамильтониане Ватсона - это кинетическая энергия, связанная с колебаниями атомов (ядер), выраженная в нормальных координатах q _s , которые, как указано выше, выражаются в терминах ядерных смещений ρ _iα формулой

${\displaystyle q_{s}=\sum _{i=1}^{N}\sum _{\alpha =1}^{3}Q_{s,i\alpha }\rho _{i\alpha }\quad \mathrm {for} \quad s=1,\ldots ,3N-6.}$

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

${\displaystyle V\approx {\frac {1}{2}}\sum _{s=1}^{3N-6}f_{s}q_{s}^{2}.}$

См. Также [ править ]

• Компьютерные программы по квантовой химии
• Адиабатический процесс (квантовая механика)
• Принцип Франка – Кондона
• Приближение Борна – Оппенгеймера.
• Метод GF
• Условия Эккарта
• Жесткий ротор

Ссылки [ править ]

Дальнейшее чтение [ править ]

• Родился, Макс ; Оппенгеймер, Роберт (25 августа 1927 г.). "Zur Quantentheorie der Molekeln". Annalen der Physik . 389 (20): 457–484. Bibcode : 1927AnP ... 389..457B . DOI : 10.1002 / andp.19273892002 .
• Мосс, RE (1973). Продвинутая молекулярная квантовая механика . Чепмен и Холл . ISBN 978-0-412-10490-9.
• Тинкхэм, Майкл (2003). Теория групп и квантовая механика . Dover Publications . ISBN 978-0-486-43247-2.
• Читаемое и подробное обсуждение членов спина в молекулярном гамильтониане находится в: McWeeny, R. (1989). Методы молекулярной квантовой механики (2-е изд.). Лондон: Академ. ISBN 978-0-12-486550-1.