Приближение Борна-Оппенгеймера: разделение ядер и электронов

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

Почему движение ядер и электронов вообще можно разделить

Ключевой физический факт - колоссальная разница масс: даже самое лёгкое ядро, протон, тяжелее электрона примерно в 1836 раз, а для типичных атомов это отношение измеряется десятками тысяч. При одинаковых силах, действующих со стороны электрических полей друг друга, тяжёлые ядра движутся медленно, а лёгкие электроны - быстро. Поэтому с точки зрения электрона ядра выглядят почти неподвижными: пока ядра сколько-нибудь заметно сдвинутся, электронное облако успевает много раз «подстроиться» под их текущее положение. Это и есть содержание адиабатического приближения: электронная подсистема мгновенно отслеживает медленно меняющиеся координаты ядер, оставаясь всё время в своём мгновенном стационарном состоянии.

Формально это означает, что полную волновую функцию молекулы $\Psi(r, R)$, зависящую от координат электронов $r$ и ядер $R$, ищут в виде произведения:

$\Psi(r, R) \approx \psi_e(r; R)\,\chi(R),$

где $\psi_e(r; R)$ - электронная волновая функция, зависящая от координат ядер $R$ лишь как от параметра, а $\chi(R)$ - волновая функция ядер. Такое разделение переменных было бы точным, если бы электронная часть вовсе не зависела от скорости ядер; приближение Борна-Оппенгеймера утверждает, что зависимостью от этой скорости можно пренебречь в первом порядке, и оценивает, насколько малы отброшенные поправки.

Что даёт приближение: поверхность потенциальной энергии

Рассчитать для своей задачи

Практический рецепт приближения Борна-Оппенгеймера такой: сначала при фиксированном наборе координат ядер $R$ решают электронное уравнение Шрёдингера

$\hat{H}_e\,\psi_e(r; R) = E_e(R)\,\psi_e(r; R),$

где $\hat{H}_e$ включает кинетическую энергию электронов, их взаимное отталкивание и притяжение к неподвижным ядрам. Решая эту задачу для разных значений $R$, получают функцию $E_e(R)$ - так называемую поверхность потенциальной энергии (для двухатомной молекулы это просто кривая от одной переменной, межъядерного расстояния). К ней добавляется кулоновское отталкивание самих ядер, и итоговая функция $U(R) = E_e(R) + V_{nn}(R)$ становится потенциалом, в котором на втором шаге двигаются уже ядра:

$\left[-\frac{\hbar^2}{2\mu}\nabla_R^2 + U(R)\right]\chi(R) = E\,\chi(R).$

Именно эта картина лежит в основе всей молекулярной спектроскопии: у молекулы есть не одна, а несколько адиабатических кривых $U(R)$ - по одной на каждое электронное состояние, - и переходы между ними видны в спектре как отдельные электронные полосы, внутри каждой из которых есть ещё более тонкая колебательная и вращательная структура.

Иерархия энергетических масштабов: почему спектры расслаиваются

Рассчитать для своей задачи

Разницу масс электрона и ядра удобно свести к одному безразмерному малому параметру - параметру адиабатичности:

$\kappa = \left(\frac{m_e}{\mu}\right)^{1/4},$

где $m_e$ - масса электрона, а $\mu$ - приведённая масса ядер молекулы. Оценка по порядку величины (её приводят практически во всех курсах квантовой механики при выводе приближения Борна-Оппенгеймера) показывает, что характерные масштабы энергии электронного, колебательного и вращательного движений связаны через степени этого параметра:

$\Delta E_{эл} : \Delta E_{кол} : \Delta E_{вращ} \;=\; 1 : \kappa^2 : \kappa^4.$

Иными словами, если ввести коэффициент $r = \kappa^2 = \sqrt{m_e/\mu}$, то колебательный масштаб энергии составляет примерно $\Delta E_{кол} \approx r\,\Delta E_{эл}$, а вращательный - $\Delta E_{вращ} \approx r^2\,\Delta E_{эл}$. Для лёгкой молекулы вроде H₂, где приведённая масса ядер $\mu \approx 0{,}5$ а.е.м., параметр $\kappa \approx 0{,}18$, а $r \approx 0{,}033$: при электронном переходе порядка 10 эВ колебательный масштаб получается около 0,33 эВ, а вращательный - около 0,011 эВ, то есть на два порядка меньше электронного. Именно поэтому в молекулярных спектрах отчётливо видна трёхуровневая иерархия - грубые электронные полосы, внутри них колебательная структура и внутри неё ещё более тонкая вращательная структура, - а не сплошной непрерывный набор линий.

Малый параметр и точность приближения

Чем тяжелее ядра, тем меньше параметр $\kappa$ и тем лучше выполняется приближение Борна-Оппенгеймера: для молекул с тяжёлыми ядрами (например, I₂, где приведённая масса ядер около 63,5 а.е.м.) $\kappa$ заметно меньше, чем для лёгкого H₂, и разделение электронного и ядерного движений работает точнее. Это согласуется с общей логикой адиабатического приближения: параметр $\kappa$ фактически показывает, во сколько раз период характерного электронного движения короче периода колебаний ядер, и чем сильнее это неравенство, тем надёжнее электроны успевают «подстроиться» под ядра. Поправки к энергии, которые приближение отбрасывает на первом шаге (неадиабатические поправки, связывающие разные электронные состояния через движение ядер), по порядку величины пропорциональны степеням того же $\kappa$, поэтому его малость одновременно объясняет и хорошую точность приближения, и то, почему её вообще можно контролируемо оценивать.

Где приближение работает плохо

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

Пример решения типовой задачи

Пусть нужно оценить масштабы уровней энергии молекулы CO, у которой приведённая масса ядер $\mu = 6{,}857$ а.е.м., если энергия характерного электронного перехода составляет $\Delta E_{эл} = 8$ эВ. Сначала находим коэффициент разделения масштабов:

$r = \sqrt{\frac{m_e}{\mu}} = \sqrt{\frac{5{,}486 \cdot 10^{-4}}{6{,}857}} \approx 0{,}00894.$

Отсюда колебательный масштаб энергии:

$\Delta E_{кол} \approx r\,\Delta E_{эл} = 0{,}00894 \cdot 8 \approx 0{,}072\ \text{эВ},$

а вращательный масштаб - на порядок меньше:

$\Delta E_{вращ} \approx r^2\,\Delta E_{эл} = (0{,}00894)^2 \cdot 8 \approx 0{,}00064\ \text{эВ}.$

Параметр адиабатичности при этом $\kappa = \sqrt{r} \approx 0{,}0945$ - заметно меньше, чем у лёгкого H₂, что и должно быть для молекулы с более тяжёлыми ядрами. Такой расчёт даёт лишь порядок величины, а не точное спектроскопическое значение, но именно этот порядок объясняет, почему в спектре CO электронные, колебательные и вращательные полосы так сильно разнесены друг от друга.

Частые ошибки

• Путают приближение Борна-Оппенгеймера с полным пренебрежением движением ядер. Ядра всё равно двигаются - просто по уже готовой поверхности потенциальной энергии, а не одновременно с электронами.
• Забывают, что параметр $\mu$ - это приведённая, а не полная масса ядер. Для двухатомной молекулы $\mu = m_1 m_2/(m_1+m_2)$, а не сумма масс.
• Считают, что $\kappa$ должен быть одинаков для всех молекул. На деле параметр зависит от приведённой массы ядер: чем тяжелее ядра, тем меньше $\kappa$ и тем точнее приближение.
• Путают порядок степеней в иерархии масштабов. Соотношение $1 : \kappa^2 : \kappa^4$, а не $1 : \kappa : \kappa^2$ - колебательный масштаб отличается от электронного на $\kappa^2$, а не на $\kappa$.
• Считают приближение точным всегда. Вблизи пересечений термов и в вырожденных электронных состояниях неадиабатические поправки перестают быть малыми, и однопотенциальная картина перестаёт работать.

FAQ

Что означает символ $\mu$ в формуле параметра адиабатичности? Это приведённая масса ядер молекулы (для двухатомной молекулы - $\mu = m_1 m_2/(m_1+m_2)$), а не масса электрона и не сумма масс ядер.

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

Как связаны электронные, колебательные и вращательные уровни энергии? Их характерные масштабы соотносятся как $1 : \kappa^2 : \kappa^4$, где $\kappa = (m_e/\mu)^{1/4}$ - малый параметр адиабатичности. Поэтому колебательная структура спектра в сотни раз теснее электронной, а вращательная - ещё теснее.

Коротко

Приближение Борна-Оппенгеймера разделяет задачу о молекуле на два шага: сначала электроны находят своё состояние при фиксированных ядрах, образуя поверхность потенциальной энергии, а затем ядра двигаются уже по этой поверхности. Законность разделения обеспечивается малостью параметра $\kappa = (m_e/\mu)^{1/4}$, из которого напрямую следует иерархия энергетических масштабов молекулы $1 : \kappa^2 : \kappa^4$ - электронного, колебательного и вращательного. Чем тяжелее ядра, тем меньше $\kappa$ и тем точнее приближение, а нарушается оно там, где электронные поверхности потенциальной энергии сближаются или пересекаются.