Метод МО ЛКАО

Матеріал з Вікіпедії — вільної енциклопедії.
Перейти до: навігація, пошук

Ме́тод МО ЛКАО (англ. LCAO approximation) — квантово-хімічний метод розрахунку електронної структури молекул, заснований на наближенні, в якому молекулярні орбіталі (МО) представлені у вигляді лінійної комбінації атомних орбіталей (ЛКАО).

ЛКАО — найпростіший метод визначення хвильових функцій молекулярних орбіталей. Він розглядає хвильові функції молекулярних орбіталей як лінійні комбінації хвильових функцій атомних орбіталей. Для точного визначення хвильової функції молекулярної орбіталі необхідно розв'язати складну навіть для найпростіших молекул задачу про рух одного електрона в самоузгодженому полі, створеному атомними ядрами і рештою електронами всіх атомів, що входять в молекулу. Тому в методі МО ЛКАО використовуються припущення, які спрощують вихідну задачу.

У методі молекулярних орбіталей повна хвильова функція молекули складається із хвильових функцій, що описують поведінку окремих електронів в полі, що створюється іншими електронами і всіма атомними ядрами. Тим самим, концепція МО близька до концепції АО з тією тільки різницею, що МО є багатоцентровими орбіталями.

В принципі для знаходження одноелектронних хвильових функцій МО можна використувати метод Хартрі-Фока і отримати таблиці їх числових значень подібно до того, як це робиться для атома. Недоліком такого підходу крім відсутності рішення в аналітичному вигляді є великі математичні труднощі розрахунку хвильової функції. Бажано одержати аналітичне уявлення молекулярної орбіталі, при чому таке, щоб зберігалася можливість її модифікації в рамках відпрацьованих процедур. Найпоширенішим наближенням до істинної МО є уявлення молекулярної орбіталі у вигляді лінійної комбінації атомних орбіталей, що утворюють молекулу, тобто:

Форма данного рівняння має дві суттєві переваги. По-перше, коли електрон знаходиться на деякій МО поряд з ядром атомa v, його поведінка і хвильова функція повинні бути близькими до відповідних атомних характеристик. Ця вимога добре запезпечується розкладом МО ЛКАО. По-друге, вигляд атомних орбіталей залишеється незмінним, а покращення хвильової функції МО можливе шляхом підбору оптимальних коефіцієнтів розкладу для знаходження яких використовується варіаційний метод Рітца.

Наближення ЛКАО до пошуку вигляду МО і уявлення повної хвильової фукнції молекули у вигляді слетерівського визначника веде, в рамках методу Хартрі-Фока з використанням електронного гамільтоніана, до рівнянь, вперше отриманих Рутаном (1951). Ці рівняння є наближенням до рівнянь Хартрі-Фока і лежать в основі майже всіх сучасних неемпіричних методів розрахунку складни молекулярних систем. Вони служать також вихідними для розвитку всіх основних напівемпіричних теорій методу МО.

Рівняння Рутаана[ред.ред. код]

Розглянемо систему із закритою оболонкою, що складається з 2n електронів на n молекулярних орбіталях, локалізованих на N ядрах атомів. Повна енергія молекули визначається співвідношенням

Підставимо в цей вираз оператор Гамільтона в наближенні Борна-Оппенгеймера, а також хвильову функцію у вигляді детермінанта Слетера, елементами якого є молекулярні спін-орбіталі р.

Гамільтоніан молекули з N ядрами і n електронами має вигляд:

Тут грецькі символи належать ядрам, а римські - електронам. Оператор доповнений членом, що враховує кінетичну енергію ядер, а оператор потенціальної енергії тепер містить три складових: енергію електрон-електронного відштовхування, енергію притягання кожного електрону до ядра молекули і потенціальну енергію відштовхування позитивно заряджених ядер. Неважко зауважити, що гамільтоніан молекули залежить не тільки від координат електронів, але і від ядерних координат.

Так само і для багатоелектронного атома, умова ортонормованості хвильових фунцій дозволяє записати вираз для повної енергії молекули у вигляді:

Врахувавши, що хвильові функції МО, які входять до виразу для кулонівського і обмінного інтегралів, є лінійними комбінаціями АО, кожна з яких помножена на невідомий коефіцієнт розкладу . Коефіцієнти знаходять за допомогою варіаційного методу Рітца із умови мінімуму повної енергії молекули:

Повторюючи повністю всі міркування, зроблені при виведенні рівняння Хартрі-Фока, використовуючи варіаційний принцип Рітца, прийдемо до рівняння Рутана:

Де - одноелектронна енергія молекулярної орбіталі , - елемент N xN матриця, що називається матрицею перекривання

Накінець, - елемент другої NxN матриці, що називається матрицею Фока

В цьому виразі є матричним уявленням енергії електрону в полі "голих" ядер:

де

Величина є двохелектронний інтеграл відштовхування:

Вони помножені на одноелектронну матрицю рядків зв'язку або матрицю густини :

Сумування ведеться за всіма занятими електронами МО. Коефіцієнт 2 відбиває, що кожну молекулярну орбіталь займають два електрони. Розглянемо параметри, що входять у рівняння Рутана, можемо записати вираз для електронної енергії молекули у наступному вигляді:

Повна енергія молекули отримується доданням до виразу, що враховує міжядерне відштовхування:

Рівняння Рутана є нелінійним однорідним рівнянням відносно невідомих величин , так як в ці рівняння входить параметр, що залежний від цих невідомих коефіцієнтів. Систему рівнянь Рутана розв'язують за допомогою процедури самоузгодження. Спочатку здаються початковими коефіцієнтами, з допомогою яких розраховують початкове наближення матриці Фока. Враховуючи на цьому етапі, що вона не залежить від шуканих коефіцієнтів, розклавши її отримають систему однорідних лінійних рівнянь:

Ця система маэ нетривыальний розвязок, коли її детермінант рівний нулю

Звідки набір , підставляючи його в попереднє рівняння, знаходять . З допомогою цих коефіцієнтів знову розраховують елементи матриці Фока, потім . Таким чином, проводять цикл розрахунків:

до тих пір, поки повна енергія і матриця порядків звязків для двох послідовних розрахунків (ітерацій) не будуть співпадати із заданою точністю:

де i - малі параметри, що визначають умови збігу розвязків за повною енергією і матрицею густини відповідно.

Є також метод ППП (який запропонували Р.Паризер, Р.Г.Пар та Дж. Попл), який є напівемпіричним варіантом методу Рутана. В свою чергу, метод Рутана - це ЛКАО наближення до методу самоузгодженого поля Хартрі-Фока. У методі ППП використовується наближення нульового деференціального перекривання (НДП), тобто нехтують не тільки інтегралами перекривання, а й усіма інтегралами, які містять добутки типу , де і - атомні орбіталі1.

Приклад[ред.ред. код]

Приклад здійснення процедури самоузгодження для знаходження повної енергії молекули води методом Хартрі-Фока-Рутаана.

Iteration Ex Total energy E Change Density change
1 0 -75.368188477 -75.368188477 0.539575981
2 1 -75.571980468 -0.203791991 0.137147099
3 2 -75.582758173 -0.010777705 0.051588713

Збіг розвязку хвильового рівняння досягається за 8 ітерацією при критерію за матрицею . Після досягнення збігу проводиться одна додаткова ітерація

Iteration Ex Total energy E Change Density change
4 3 -75.583872052 -0.001113878 0.002713332
5 4 -75.583913740 -0.000041689 0.001781922
6 5 -75.583923021 -0.000009280 0.000183589
7 6 -75.583923056 -0.000000035 0.000031273
8 7 -75.583923062 -0.000000007 0.000004804
9 8 -75.583923062 0.000000000 0.000001855

Розраховане значення повної енергії молекули води для даного прикладу рівне -75.583923 Хартрі

(1 Хартрі=(/a)=627,5096 ккал/моль=2625,500 кДж/моль). Шість значущих цифр після крапки відповідають точності розрахунку не гірше 0.01 кДж/моль. Як видно із прикладу, необхідна точність була досягнута уже на шостому циклі ітераційного процесу, тобто збіг за енергією зазвичай досягається раніше.

Розв'язок електронного хвильового рівняння методом молекулярних орбіталей в наближенні лінійних комбінацій атмоних орбіталей вимагає використання процедури самоузгодження, коли шукані хвильові функції МО повинні бути узгоджені з ефективним полем молекули, що називається методом самоузгодженого поля.

Кількість власних значень і відповідних їм власних значень функцій ( рівне числу базисних функцій N в розкладі ЛКАО). Визначник Слетера для повної хвильової функції молекули будуються із n=N/2 зайнятих електронами МО. В мінімізації повної енергії беруть участь тільки зайняті МО і, так як матричні елементи залежать тільки від , а порядок зв'язку розраховується із хвильових функцій тільки зайнятих орбіталей, лише вони можуть розглядатися як фізично визначені. Не зайняті МО, отримані із рівнянь Рутаана, не беруть участі в мінімізації повної енергії системи, тому їх відповідність істинним енергетичним рівням молекули не зовсім визначено. Такі рівняння називають віртуальними.

1Віталій Яцимирський - Фізична хімія, с.257

Див. також[ред.ред. код]

Посилання[ред.ред. код]

Джерела[ред.ред. код]

  • Глосарій термінів з хімії // Й. Опейда, О. Швайка. Ін-т фізико-органічної хімії та вуглехімії ім. Л. М. Литвиненка НАН України, Донецький національний університет. — Донецьк: Вебер, 2008. — 758 с. — ISBN 978-966-335-206-0.
  • Базилевский М. В. Метод молекулярных орбит и реакционная способность органических молекул. — М. : Химия, 1969. — 304 с. (рос.)


CHEM: Це незавершена стаття з хімії.
Ви можете допомогти проекту, виправивши або дописавши її.