А.М. Попов, О.В. Тихонова - Лекции по атомной физике (1121316), страница 51
Текст из файла (страница 51)
на рис.15.2):Hˆ = Hˆ + Tˆ + Tˆ + e 2 R .e12Здесь Hˆ e = Tˆe + Ve - электронная часть гамильтониана, Tˆe - кинетическая энергия элекe2 e2−- энергия взаимодействия электрона с протонами, r1 и r2 - удалеr1 r2ние электрона от первого и второго протонов соответственно, Tˆ1 , Tˆ2 - кинетическаяrvэнергия ядер, R = R1 − R2 - межъядерное расстояние. Стационарные состояния системытрона, Ve = −могут быть определены из решения уравнения Шредингераr r rr r rHˆ Φ (r , R1 , R2 ) = EΦ (r , R1 , R2 ) .(15.4)Дальнейшее рассмотрение нам будет удобнопроводить в системе координат, связанной сцентром масс молекулы.
Пренебрегая массойэлектрона m по сравнению с массой ядра M ,введем координату центра масс2rrrℜ = R1 + R1 / 2и координату относительного движенияr rrR = R2 − R1 .rrТогда координаты R1 и R2 каждого из ядер определяются какrr rR1 = ℜ − R 2 ,(15.5)rr rR2 = ℜ + R 2 .Получим выражение для кинетической энергииядер в системе координат, связанной с центром масс. Учитывая, что1 ∂∂∂∂∂1 ∂r =r − r и r =r + r,∂R1 2 ∂ℜ ∂R∂R2 2 ∂ℜ ∂Rа также∂2∂21 ∂21 ∂2∂2∂2∂2∂2r2 =r 2 + r2 − r r ,r2 =r 2 + r2 + r r ,4 ∂ℜ∂ℜ∂R ∂R2 4 ∂ℜ∂ℜ∂R∂R∂R∂R1получимh2 ∂2h2 ∂2h2h2 ∂2∂2r2 −r2 = −r2 −r ,Tˆ1 + Tˆ2 = −4 M ∂ℜM ∂R2(2 M ) ∂ℜ2µ ∂R 2(2)Мы здесь предполагаем, что массы ядер атомов, образующих молекулу, равны M 1 = M 2 = M .
В слу-чае, если массы ядер являются различными, для ℜ имеем(rrrℜ = M 1 R1 + M 2 R2) (M1+ M2 ).201202где µ = M 2 - приведенная масса молекулы. Таким образом, мы получили, что полныйгамильтониан молекулы представим в виде⎛r r rr r h2 ∂2e2 ⎞∂2h2r 2 + ⎟⎟ −r .Hˆ (r , R, ℜ) = ⎜⎜ Tˆe + Ve (r , R) −2µ ∂RR ⎠ 2(2 M ) ∂ℜ 2⎝rЗдесь r - координата электрона, отсчитываемая от центра масс молекулы, при этомr r rr r rr1 = r + R 2 , r2 = r − R 2 .(15.6)rrПолученное выражение означает, что полный гамильтониан молекулы Hˆ (r , R, ℜ) можнопредставить в виде суммы гамильтонианов двух подсистем, одна из которых определяется совокупностью координат, описывающих относительное движение электронов иrядер в молекуле, а другая – координатой центра масс молекулы ℜ .
Значит, волноваяr r rфункция стационарного состояния молекулярной системы Φ (r , R, ℜ) может быть представлена в виде произведения волновых функцийrrr r rr rΦ (r , R, ℜ) = Ψ (r , R) ⋅ exp(iKℜ) ,где второй сомножитель определяет волновую функцию свободного движения молекулыrкак целого в пространстве, K - волновой вектор молекулы. Вся информация о внутренr rней структуре молекулы содержится в функции Ψ (r , R ) , которая является собственнойфункцией гамильтонианаr rr r h2 ∂2e2r2 + .Hˆ (r , R) = Tˆe + Ve (r , R ) −(15.7)2µ ∂RRРешение задачи на собственные значения и собственные функции гамильтониана(15.7) обычно проводят в адиабатическом приближении (приближении Борна - Оппенгеймера3).
Суть этого приближения заключается в следующем. Молекула состоит издвух подсистем – электронной (легкой) и ядерной (тяжелой). Движение электронов происходит с существенно больше скоростью, чем движение ядер. Поэтому приближенноможно считать, что движение электронов происходит при неподвижных ядрах. При рассмотрении ядерной подсистемы мы будем полагать, что электронная подсистема успевает подстроиться под мгновенное положение ядер.Таким образом, в рамках приближения Борна – Оппенгеймера волновая функциямолекулярной системы представляется в видеrr rr rΨ ( r , R ) = ψ ( e ) ( r , R )φ ( N ) ( R ) ,(15.8)r(e) rψ (r , R) - электронная волновая функция молекулы, зависящая, в том числе, от ядерrной координаты, как от параметра, φ ( N ) ( R ) - волновая функция ядерной подсистемы молекулы.
При этом волновая функция электронной подсистемы является решением задачи на собственные значения для электронного гамильтонианаr rr rr rTˆe + Ve (r , R ) ψ ( e ) (r , R ) = E ( e ) ( R )ψ ( e ) (r , R) .(15.9)()(e)Здесь E ( R) - энергия электронной подсистемы молекулы, которая тоже зависит отмежъядерного расстояния как от параметра. Решая (15.9) мы получаем полный наборэлектронных состояний молекулы, зависящих, в том числе, и от расстояния между ядрами.Будем теперь искать стационарные состояния гамильтониана (15.7) в виде разложения (15.8), где электронные волновые функции удовлетворяют (15.9).
Получим:3R.Oppenheimer (1904-1967) – американский физик – теоретик.202203(Hˆ ψe(e))()rr⎛r rr re2 ⎞(r , R) φ ( N ) ( R ) + ⋅⎜⎜ TˆN + ⎟⎟ ψ ( e ) (r , R )φ ( N ) ( R ) = Eψ ( e ) φ ( N ) .R⎠⎝(15.10)Здесьh2 ∂2h2r 2 = − ∇ 2RTˆN = −2µ ∂R2µ- оператор кинетической энергии ядерной подсистемы.Наиболее существенно для дальнейшего, что слагаемые, содержащее дифференцирование электронной волновой функции по ядерной координате (их называют неадиабатическими), малы по сравнению с другими членами, входящими в уравнение(15.10). Мы не будем останавливаться на доказательстве этого утверждения4, однако, сматематической точки зрения именно это утверждение представляет собой основу адиабатического приближения.Пренебрегая неадиабатическими членами в уравнении (15.10), разделим обе часrr rти уравнения на ψ ( e ) (r , R )φ ( N ) ( R ) .
Тогда имеем1 ˆ (e) r r1 ⎛ˆe2 ⎞ (N)⎜(15.11)H e ψ (r , R ) = E − ( N ) ⎜ TN + ⎟⎟φ ( R ) .R⎠ψ (e)φ ⎝Учитывая, что в силу (15.9) левая часть (15.11) есть энергия электронной подсистемы,получим следующее уравнение для определения волновых функций стационарных состояний ядерной подсистемыrr⎛ˆ⎞e2⎜⎜ TN ++ E ( e ) ( R ) ⎟⎟φ ( N ) ( R ) = Eφ ( N ) ( R ) .(15.12)R⎝⎠Как видно, мы получили «обычное» уравнение Шредингера с эффективным потенциалом взаимодействия ядер в молекулеe2( eff )( R) =+ E (e) ( R) .(15.13)VRЭффективный потенциал учитывает не только кулоновское отталкивание между ядрами, но иналичие электронной плотностив системе, которая обеспечиваетпритяжение в системе, и можетпривести к связыванию атомов вмолекулу.
Для этого, очевидно,необходимо, чтобы зависимостьV ( eff ) ( R) характеризовалась кривой с локальным минимумом(см. рис.15.3, кривая 1).Зависимость эффективной энергии взаимодействияядер в молекуле от межъядерного расстояния называют электронным термом молекулы. Эту(4)Можно показать, что эти слагаемые имеют, по крайней мере, порядок малости ξ =корнем стоит отношение масс электронной и ядерной подсистемы.m M , где под203204же зависимость часто называют потенциальной кривой (потенциальной поверхностью5).Может так оказаться, что потенциальная кривая характеризуется отсутствием минимума энергии (такой случай также приведен на рис. 15.3).
Тогда, очевидно, при любом расстоянии между атомами существуют силы отталкивания, то есть устойчивое связывание атомов в молекулу невозможно. В этом случае говорят, что электронный терммолекулы является разлетным.Итак, подведем некоторые итоги. При рассмотрении структуры молекулы в рамках адиабатического приближения необходимо сначала решить задачу на собственныезначения и собственные функции для гамильтониана электронной подсистемы молекулы. В качестве решения получаются набор собственных функций и собственных значений, параметрически зависящих от расстояния между ядрами в молекуле.
Далее мыможем построить набор электронных термов молекулы (каждому электронному состоянию молекулы соответствует конкретный вид зависимости V ( eff ) ( R) ) и проанализировать задачу о пространственном строении молекулы и особенностях движенияядер для каждого из термов.Рассмотрим теперь более подробно электронное строение ряда простейших молекулярных систем.Молекулярный ион водорода.Приближенную электронную волновую функцию основного состояния молекулярного иона водорода будем искать в виде суперпозиции двух невозмущенных функций водородного атома6 (обозначения см. на рис.15.4):r rrrψ (r , R) = αψ 1s (r1 ) + βψ 1s (r2 ) .(15.14)rrЗдесь r1 и r2 - удаление электрона от каждого из протонов, каждая из этих величин выrrражается через переменные r и R . Приближение (15.14), очевидно, выполнено длябольших межъядерных расстояний R >> a 0 .
Однако, для нас в дальнейшем существенно, чтопредставление волновой функции электрона в молекуле (молекулярной орбитали) в виде линейнойкомбинации атомных волновых функций (атомных орбиталей) приводит к правильным результатам, хотя бы и на качественном уровне, для произвольного соотношения параметров R и a0 .Вследствие тождественности ядерных центров естественно предположить, что вероятностиобнаружить электрон около каждого из них равныr r22между собой, то есть α = β . Отсюда ясно, что функция ψ (r , R) представима в видесимметричной, или антисимметричной комбинации атомных функций7r rrrψ ± (r , R) = C ± ( R )[ψ 1s (r1 ) ± ψ 1s (r2 )] ,(15.15)±где C ( R) - нормировочная константа, определяемая из условия5Термин «потенциальная поверхность» (potential surface) используется, прежде всего, для многоатомныхмолекул, для которых обсуждаемая зависимость определяется функцией от совокупности ядерных координат и описывается поверхностью в многомерном пространстве.6Это приближение иногда называют приближением разъединенных атомов.7rrЗдесь координаты r1 и r2 связаны с координатамиr rr и R соотношениями (15.6)204205∫ψ±r rr rr r 2(r , R ) d 3 r = ψ ± (r , R ) ψ ± (r , R ) = 1 .(15.16)Вычисление интеграла (15.16) дает(C ± ( R))2 ( ψ1s (rr1 ) ψ1s (rr1 ) + ψ1s (rr2 ) ψ1s (rr2 ) ± 2 ψ1s (rr1 ) ψ1s (rr2 ) ) = 1 ,откуда находим11,(15.17)C ± ( R) =2 1 ± S ( R)rrгде S ( R ) = ψ 1s (r1 ) ψ 1s (r2 ) - интеграл перекрытия.