Метод молекулярних орбіталей

Мéтод молекулярних орбітáлей (англ. MO method, рос. метод МО) — це метод наближеного розв'язання електронного рівняння Шредингера для багатоелектронних молекулярних систем. Ґрунтується на побудові повної хвильової функції у вигляді антисиметризованого добутку молекулярних орбіталей, його зручно записувати як детермінант. Молекулярні орбіталі, у свою чергу, зазвичай представляють як лінійні комбінації атомних орбіталей (наближення МО ЛКАО).

Геометрія дескрипції д-орбіталей

Розрахунковий квантово-хімічний метод, в якому використано одноелектронне наближенняi кожному з електронiв у молекулi приписуються окремi хвильовi функцiї — молекулярнi орбiталі, з яких будується загальна хвильова функцiя всiх електронiв молекули. Лежить в основі як напівемпіричних методів розрахунку, так і метода ab initio. Для побудови моле-кулярних орбіталей використовуються одноелектронні атомні орбіталі.

Короткі історичні відомості

Підґрунтя методу було закладене Гундом [1] (ідея про середнє поле тут висловлена дещо раніше за Гартрі), Маллікеном [2] і Леннард-Джонсом [3] (аналіз електронних конфігурацій атомів та двохатомних молекул, в тому числі пояснення парамагнетизму кисню). Від двохатомних молекул Маллікен перейшов до багатоатомних[4], визначивши роль симетрії молекули та її вимоги до молекулярних орбіталей. Хюккель в 1930—1932 оприлюднив цикл статей про подвійний зв'язок та молекулу бензену; виявилось, що метод молекулярних орбіталей теж, як і метод валентних зв'язків, дозволяє пояснити електронне спряження. Пізніше Ґепперт-Майєр та Скляр[5] застосували метод молекулярних орбіталей до обчислення енергій збудження молекули бензену; моделювання УФ спектрів залишається актуальною задачею квантової хімії й досі.

Практично сучасної математичної форми (матричне формулювання в наближенні МО ЛКАО) методу надав Рутан[6]. Ця форма використовується в більшості квантовохімічних програм.

Математичні основи

Молекулярну спін-орбіталь представляють у вигляді добутку просторової частини φ (молекулярної орбіталі) та спінової частини σ (яка в цьому випадку відповідає орієнтації спіну «вгору» α чи «вниз» β):

\psi ({\overrightarrow {r_{1}}},s_{1})=\phi ({\overrightarrow {r_{1}}})\cdot \sigma (s_{1})

Тут набір 4-х змінних $\{{\overrightarrow {r_{1}}},s_{1}\}$ є координатами електрона під номером 1. Повна хвильова функція Ψ будується як антисиметризований добуток тих спін-орбіталей, які заповнені електронами. Його зазвичай записують як детермінант Слейтера. Приміром, так її записують для чотириелектронної системи із замкненою електронною оболонкою, як в атомі берилію (для скорочення запису замість ${\overrightarrow {r_{1}}}$ та $s_{1}$ позначено лише їхні номери):

\Psi ={\frac {1}{2{\sqrt {6}}}}{\begin{vmatrix}\phi _{1}(1)\alpha (1)&\phi _{1}(2)\alpha (2)&\phi _{1}(3)\alpha (3)&\phi _{1}(4)\alpha (4)\\\phi _{1}(1)\beta (1)&\phi _{1}(2)\beta (2)&\phi _{1}(3)\beta (3)&\phi _{1}(4)\beta (4)\\\phi _{2}(1)\alpha (1)&\phi _{2}(2)\alpha (2)&\phi _{2}(3)\alpha (3)&\phi _{2}(4)\alpha (4)\\\phi _{2}(1)\beta (1)&\phi _{2}(2)\beta (2)&\phi _{2}(3)\beta (3)&\phi _{2}(4)\beta (4)\end{vmatrix}}

Нормувальний множник перед детермінантом такий, що його квадрат дорівнює числу доданків в розкритому детермінанті. Конструкція такого детермінанта також визначає загальний спін системи електронів. Наведений приклад відповідає двом молекулярним орбіталям φ₁ та φ₂, кожна з яких заповнена двома електронами, тобто загальний спін системи нульовий.

Далі, в наближенні МО ЛКАО, молекулярну орбіталь представляють як лінійну комбінацію атомних орбіталей:

\phi _{i}({\overrightarrow {r}})=\sum _{j}C_{ij}\chi _{j}({\overrightarrow {r}})

Тут χ_j — атомна орбіталь, а C_ij — числові коефіцієнти (в загальному випадку комплексні). Відповідним добором цих коефіцієнтів забезпечуються нормування та ортогональність молекулярних орбіталей; ці корисні властивості виражаються формулою

\int \phi _{i}^{*}({\overrightarrow {r}})\phi _{j}({\overrightarrow {r}})\,d{\overrightarrow {r}}={\begin{cases}0,&{\mbox{if }}i\neq j\\1,&{\mbox{if }}i=j\end{cases}}

Електронні формули атомів є фактично переліком заповнених атомних орбіталей і безпосередньо відповідають слейтерівським детермінантам (C_ij тут буде одиничною матрицею).

Коефіцієнти C_ij визначають за допомогою варіаційного принципу із умови мінімуму повної енергії. Для спрощення використовують наближення середнього поля, в якому кожен із електронів рухається в електростатичному полі, створеному ядрами та іншими електронами при усередненні їхніх позицій. Багатоелектронне рівняння розпадається в такому наближенні в набір одноелектронних із ефективним оператором енергії (оператором Фока, або фокіаном), і математично задача зводиться до знаходження власних векторів та власних значень цього оператора. Ці вектори й є (канонічними) молекулярними орбіталями, а відповідні їм власні значення — їхніми енергіями. Теорема Купманса стверджує, що у певному наближенні ці орбітальні енергії дають оцінку для потенціалів іонізації молекули.

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

Кореляційні діаграми і орбітальна симетрія

Гунд, Маллікен вивчали двохатомні молекули при зміні відстані між ядрами в діапазоні від нескінченної (ізольовані атоми) до нульової (т.зв. об'єднаний атом). Для двох чи навіть трьох точок цього діапазону (ізольовані атоми, двоатомна молекула з рівноважною довжиною зв'язку, об'єднаний атом) були відомі електронні спектри, саме це дозволило виконати аналіз та узагальнення. Кореляційні діаграми показують, як змінюється енергія орбіталі при зміні міжатомної відстані, які орбіталі переходять одна в одну. У сукупності із правилом збереження орбітальної симетрії (збереження симетрії молекули при її деформації веде до збереження типів симетрії молекулярних орбіталей) та правилом неперетинання (криві орбітальних енергій для орбіталей однакової симетрії на кореляційних діаграмах не перетинаються) такий підхід дозволив пояснити певні закономірності органічних реакцій.

Саме орбітальна симетрія та принцип її збереження дозволяє класифікувати ковалентні зв'язки на σ-тип та π-тип. Останній відіграє велику роль в поясненні реакційної здатності органічних сполук, див. концепції спряження, ароматичності тощо. Також на теорії симетрії (теорії груп) ґрунтуються такі моделі будови комплексних сполук, як теорія кристалічного поля та теорія поля лігандів.

Локалізація молекулярних орбіталей

Канонічні молекулярні орбіталі делокалізовані по всій молекулі. Це діаметрально протилежно до локалізованого двохелектронного зв'язку, яким оперує метод валентних зв'язків. Процедура локалізації молекулярних орбіталей (існує багато її варіантів[7][8][9]) не є частиною методу молекулярних орбіталей. Вона спирається на таку властивість детермінанта, як сталість його значення при утворенні будь-яких лінійних комбінацій його рядків. Локалізовані орбіталі доволі точно відтворюють ковалентні зв'язки та неподілені електронні пари, але вони, на відміну від канонічних, не характеризуються орбітальною енергією.

Граничні молекулярні орбіталі

Ідея про те, що при хімічній взаємодії вирішальну роль відіграє взаємодія вищої зайнятої молекулярної орбіталі одної молекули (нуклеофіла) із нижчою вільною іншої молекули (електрофіла), привела до створення нових способів оцінки реакційної здатності, серед яких теорія граничних орбіталей Фукуї[10] та теорія збурень молекулярних орбіталей Дьюара[11].

Шляхи уточнення методу

У деяких системах внаслідок виродження орбіталей не вдається однозначно обрати орбіталі для слетерівського детермінанту. Так, в атомі бору атомні орбіталі 2p_x, 2p_y, 2p_z характеризуються однаковою енергією, а на них розташований лише один електрон, тож в детермінант можна включити лише одну з них. У таких випадках слід використовувати лінійні комбінації еквівалентних (за енергією) детермінантів (переходити до т.зв. багатодетермінантних методів). Це врахування статичної електронної кореляції.

У загальному випадку відсутнього виродження побудова лінійних комбінацій кількох детермінантів (яким цього разу відповідають різні значення енергії) додає гнучкості хвильовій функції, тобто покращує її. Це врахування динамічної електронної кореляції.

Переваги й вади

Переваги:

На основі методу молекулярних орбіталей було створено багато інтерпретативних теорій, тобто таких, які екстрагують із повної хвильової функції хімічно значущу інформацію, безпосередньо пов'язану із структурою молекул чи реакційної здатністю.
Напівемпіричний метод Гюккеля, що використовує математичний апарат методу молекулярних орбіталей, є одним із небагатьох наближених методів, в межах якого можливі аналітичні рішення та доведення відносно широкого набору теорем.
Пристосоване до комп'ютерів матричне формулювання зумовило широке використання методу молекулярних орбіталей в сучасних квантовохімічних програмах, де він відомий під назвою методу Гартрі-Фока (англ. HF) чи методу самоузгодженого поля (англ. SCF).

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

Див. також

Джерела

Молекулярных орбиталей методы // Химическая энциклопедия : в 5 т. / гл. ред. И. Л. Кнунянц. — М. : Большая Рос. энцикл., 1992. — Т. 3 : Меди сульфиды — Полимерные красители. — Стб. 231. — Библиогр. в конце ст. — ISBN 5-85270-039-8.(рос.)
Слэтер Дж. Электронная структура молекул. — М. : Мир, 1965.(рос.)
Глосарій термінів з хімії // Й. Опейда, О. Швайка. Ін-т фізико-органічної хімії та вуглехімії ім. Л. М. Литвиненка НАН України, Донецький національний університет. — Донецьк: Вебер, 2008. — 758 с. — ISBN 978-966-335-206-0

Примітки

Hund F. Zur Deutung der Molekelspektrum. I. // Zs. Phys.. — 1927. — Т. 40. — С. 742-764. — DOI:10.1007/BF01400234.(нім.)
Mulliken R.S. The assignment of quantum number for electrons in molecules. I. // Phys. Rev.. — 1928. — Т. 32. — С. 186-222. — DOI:10.1103/physrev.32.186.(англ.)
Lennard-Jones J.E. The electronic structure of some diatomic molecules // Trans. Farad. Soc.. — 1929. — Т. 25. — С. 668-686. — DOI:10.1039/TF9292500668.(англ.)
Mulliken R.S. Electronic Structures of Polyatomic Molecules and Valence. IV. Electronic States, Quantum Theory of the Double Bond // Phys. Rev.. — 1933. — Т. 43. — С. 279-302. — DOI:10.1103/physrev.43.279.(англ.)
Goeppert-Mayer M., Sklar A.L. Calculations of the Lower Excited Levels of Benzene // J. Chem. Phys.. — 1938. — Т. 6. — С. 645-652. — DOI:10.1063/1.1750138.(англ.)
Roothaan C.C.J. New Developments in Molecular Orbital Theory // Rev. Mod. Phys.. — 1951. — Т. 23. — С. 69-89. — DOI:10.1103/RevModPhys.23.69.(англ.)
Foster J.M., Boys S.F. Canonical Configuration Interaction Procedure // Rev. Mod. Phys.. — 1960. — Т. 32. — С. 300-302. — DOI:10.1103/RevModPhys.32.300.(англ.)
Edmiston C., Ruedenberg K. Localized Atomic and Molecular Orbitals // Rev. Mod. Phys.. — 1963. — Т. 35. — С. 457-465. — DOI:10.1103/RevModPhys.35.457.(англ.)
Pipek J., Mezey P.G. A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions // J. Chem. Phys.. — 1989. — Т. 90. — С. 4916-4926. — DOI:10.1063/1.456588.(англ.)
Fukui K., Yonezawa T., Nagata C., Shingu H. Molecular Orbital Theory of Orientation in Aromatic, Heteroaromatic, and Other Conjugated Molecules // J. Chem. Phys.. — 1954. — Т. 22. — С. 1433-1442. — DOI:10.1063/1.1740412.
Dewar M.J.S. A Molecular Orbital Theory of Organic Chemistry. I. General Primciples // J. Am. Chem. Soc.. — 1952. — Т. 74. — С. 3341-3345. — DOI:10.1021/ja01133a038.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[1] Hund F. Zur Deutung der Molekelspektrum. I. // Zs. Phys.. — 1927. — Т. 40. — С. 742-764. — DOI:10.1007/BF01400234.(нім.)

[2] Mulliken R.S. The assignment of quantum number for electrons in molecules. I. // Phys. Rev.. — 1928. — Т. 32. — С. 186-222. — DOI:10.1103/physrev.32.186.(англ.)

[3] Lennard-Jones J.E. The electronic structure of some diatomic molecules // Trans. Farad. Soc.. — 1929. — Т. 25. — С. 668-686. — DOI:10.1039/TF9292500668.(англ.)

[4] Mulliken R.S. Electronic Structures of Polyatomic Molecules and Valence. IV. Electronic States, Quantum Theory of the Double Bond // Phys. Rev.. — 1933. — Т. 43. — С. 279-302. — DOI:10.1103/physrev.43.279.(англ.)

[5] Goeppert-Mayer M., Sklar A.L. Calculations of the Lower Excited Levels of Benzene // J. Chem. Phys.. — 1938. — Т. 6. — С. 645-652. — DOI:10.1063/1.1750138.(англ.)

[6] Roothaan C.C.J. New Developments in Molecular Orbital Theory // Rev. Mod. Phys.. — 1951. — Т. 23. — С. 69-89. — DOI:10.1103/RevModPhys.23.69.(англ.)

[7] Foster J.M., Boys S.F. Canonical Configuration Interaction Procedure // Rev. Mod. Phys.. — 1960. — Т. 32. — С. 300-302. — DOI:10.1103/RevModPhys.32.300.(англ.)

[8] Edmiston C., Ruedenberg K. Localized Atomic and Molecular Orbitals // Rev. Mod. Phys.. — 1963. — Т. 35. — С. 457-465. — DOI:10.1103/RevModPhys.35.457.(англ.)

[9] Pipek J., Mezey P.G. A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions // J. Chem. Phys.. — 1989. — Т. 90. — С. 4916-4926. — DOI:10.1063/1.456588.(англ.)

[10] Fukui K., Yonezawa T., Nagata C., Shingu H. Molecular Orbital Theory of Orientation in Aromatic, Heteroaromatic, and Other Conjugated Molecules // J. Chem. Phys.. — 1954. — Т. 22. — С. 1433-1442. — DOI:10.1063/1.1740412.

[11] Dewar M.J.S. A Molecular Orbital Theory of Organic Chemistry. I. General Primciples // J. Am. Chem. Soc.. — 1952. — Т. 74. — С. 3341-3345. — DOI:10.1021/ja01133a038.