6.2: Орбіталі
Хартрі Опис
Вважається, що енергії та хвильові функції в межах найбільш часто використовуваних теорій атомної структури виникають як розв'язки рівняння Шредінгера, гамільтоніанhe(r) якого має три види енергій:
- Кінетична енергія, середнє значення якої обчислюється шляхом прийняття очікуваного значення оператора−ℏ22m∇2 кінетичної енергії щодо будь-якого конкретного рішенняϕj(r) рівняння Шредінгера:KE=⟨ϕj|−ℏ22m∇2|ϕj⟩
- Куломбічна енергія тяжіння з ядром зарядуZ:⟨ϕj|−Z2er|ϕj⟩
- Кулонівські енергії відштовхування з усімаN−1 іншими електронами, які, як передбачається, займають інші атомні орбіталі (AoS) позначаютьсяϕK, причому ця енергія обчислюється як
∑K⟨ϕj(r)ϕK(r′)|e2|r−r′||ϕj(r)ϕK(r′)⟩.
Позначення Дірака⟨ϕj(r)ϕK(r′)|e2|r−r′||ϕj(r)ϕK(r′)⟩ використовується для представлення шестивимірного інтеграла Кулона
JJ,K=∫|ϕj(r)|2|ϕK(r′)|2e2r−r′drdr′
що описує кулонівське відштовхування між щільністю заряду|ϕj(r)|2 для електрона вϕj і щільністю заряду|ϕK(r′)|2 для електрона вϕK. Звичайно, сума надK повинна бути обмежена, щоб виключити,K=J щоб уникнути підрахунку «самовзаємодії» електрона в орбіталіϕj з самим собою.
Загальна енергіяϵJ орбітиϕj, являє собою суму вищевказаних трьох внесків:
ϵJ=⟨ϕj|−ℏ22m∇2|ϕj⟩+⟨ϕj|−Z2er|ϕj⟩+∑K⟨ϕj(r)ϕK(r′)|e2|r−r′||ϕj(r)ϕK(r′)⟩.
Ця обробка електронів та їх орбіталів називається рівнем теорії Хартрі. Як зазначено вище, коли екрановані гідрогенні АО використовуються для наближенняϕj таϕK орбіталей, результуючіϵJ значення не дають точних прогнозів. Наприклад, негативϵJ повинен наближатися до енергії іонізації для видалення електрона з АТϕj. Такі потенціали іонізації (ІС) можуть бути виміряні, а виміряні значення погано узгоджуються з теоретичними значеннями, коли робиться наближення грубого скринінгу для АТ.
Розширення LCAO
Для вдосконалення використання екранованих гідрогенних AoS найбільш часто наближається кожен з Hartree AoS {ϕK} як лінійну комбінацію так званих базисних AoS {χμ}:
ϕJ=∑μCJ,μχμ.
використовуючи те, що називається розширенням лінійної комбінації атомних орбіталів (LCAO). У цьому рівнянні коефіцієнти розширення {CJ,μ} - це змінні, які повинні бути визначені шляхом розв'язання рівняння Шредінгера
heϕJ=ϵJϕJ.
Підставивши розширення LCAO дляϕJ цього рівняння Шредінгера, помноживши зліва на одну з базисних AoSχν, а потім інтегруючи над координатами електрона вϕJ, отримуємо
∑μ⟨χν|he|χμ⟩CJ,μ=ϵJ∑μ⟨χν|χμ⟩CJ,μ.
Це матричне рівняння власних значень, в якомуϵJ і {CJ,μ} відображаються як власні значення та власні вектори. Матриці⟨χν|he|χμ⟩ і⟨χν|χμ⟩ називаються гамільтонових і перекриваються матрицями відповідно. Явний вираз для першого виходить при введенні більш раннього визначення he:
⟨χν|he|χμ⟩=⟨χν|−ℏ22m∇2|χμ⟩+⟨χν|−Ze2r|χμ⟩
+∑η,γ∑KCK,ηCK,γ⟨χν(r)χη(r′)|e2|r−r′||χμ(r)χγ(r′)⟩.
Важливим моментом, на який слід звернути увагу на форму матричних рівнянь Хартрі, є те, що для обчислення гамільтонової матриці потрібно знати коефіцієнти LCAO {CK,γ} орбіталів, які займають електрони. З іншого боку, ці коефіцієнти LCAO повинні бути знайдені шляхом розв'язання рівнянь власних значень матриці Хартрі. Цей парадокс призводить до необхідності розв'язання цих рівнянь ітераційно в так званій техніці самоузгодженого поля (SCF). У процесі SCF вводять початкове наближення до коефіцієнтів {CK,γ}. Це дозволяє сформувати гамільтонову матрицю, визначену вище. Матричні рівняння Хартрі
∑μ⟨χν|he|χμ⟩CJ,m=ϵJ∑μ⟨χν|χμ⟩CJ,μ
потім вирішуються для нових {CK,γ} коефіцієнтів і для орбітальних енергій {ϵK}. Нові коефіцієнти LCAO тих орбіталей, які зайняті, потім використовуються для формування нової гамільтонової матриці, після чого рівняння Хартрі знову вирішуються для іншого покоління коефіцієнтів LCAO і орбітальних енергій. Цей процес продовжують до тих пір, поки орбітальні енергії і коефіцієнти LCAO, отримані в послідовних ітераціях, не відрізняються помітно. При такій збіжності можна сказати, що самоузгоджене поле було реалізовано, оскільки коефіцієнти {CK,γ} використовуються для формування потенціалу кулонівського поля, який деталізує електронно-електронні взаємодії.
Базисні набори
Орбіталі типу Слейтера та орбіталі типу Гауса
Як зазначалося вище, можна використовувати екрановані гідрогенні орбіталі як {χμ}. Однак багато зусиль було витрачено на розробку альтернативних наборів функцій для використання в якості базових орбіталів. Результатом цих зусиль стало створення двох видів функцій, які в даний час широко використовуються. Базові орбіталі, які зазвичай використовуються в процесі LCAO, діляться на два початкові класи:
- Орбіталі типу Слейтера (cn,l,m(r,θ,ϕ)=Nn,l,m,zYl,m(θ,ϕ)rn−1e−zrSTO) характеризуютьсяnl квантовими числамиm та показниками (які характеризують радіальний розмір орбіти)z. СимволNn,l,m,z позначає константу нормалізації.
- Декартові орбіталі типу Гаусса (ГТО)ca,b,c(r,θ,ϕ)=N′a,b,c,axaybzce−ar2характеризуються квантовими числамиab, іc, які деталізують кутову форму і напрямок орбіти, і експоненти,a які регулюють радіальний розмір.
Для обох типів АО координатиrθ, іϕ відносяться до положення електрона щодо набору осей, прикріплених до ядра, на якому розташована орбітальна основа. Зверніть увагу, що орбіталі типу Слейтера (STO) схожі на гідрогенні орбіталі в області, близькій до ядра. Зокрема, вони мають ненульовий нахил біля ядра
ddr(e−zr)r=0=−z.
На відміну від цього, GTO, мають нульовий нахил поблизу,r=0 оскільки
ddr(e−ar2)r=0=0.
Ми говоримо, що STO демонструють стрибокr=0, що характерно для гідрогенних розчинів, тоді як ГТО цього не роблять. Хоча STO мають належну поведінку cusp поблизу ядер, вони використовуються в основному для обчислень атомних та лінійних молекул, оскільки багатоцентрові інтеграли
⟨χμ(1)χκ(2)|e2|r1−r2||χν(1)χγ(2)⟩
які виникають при розрахунках поліатомних молекул (ці інтеграли ми обговоримо пізніше в цій главі) не можуть бути ефективно оцінені при використанні STO. На відміну від цього, такі інтеграли можуть бути зазвичай обчислені при використанні GTO. Ця фундаментальна перевага ГТО призвела до домінування цих функцій у молекулярній квантовій хімії.
Щоб подолати первинну слабкість функцій GTO (тобто їх радіальні похідні зникають у ядрі), прийнято поєднувати два, три або більше ГТО з комбінаційними коефіцієнтами, які фіксуються і не розглядаються як параметри LCAO, у нові функції, які називаються контрактними ГТО (CGTO). Як правило, ряд радіально щільних, середніх і сипучих ГТО множаться на коефіцієнти стиснення і підсумовуються для отримання CGTO, який наближає належний крок в ядерному центрі (хоча жодна така комбінація ГТО не може точно виробляти такий крок, оскільки кожен GTO має нульовий нахил наr=0.
Хоча більшість розрахунків на молекулах зараз виконуються за допомогою гаусових орбіталів, слід зазначити, що інші базові множини можуть бути використані до тих пір, поки вони охоплюють достатню кількість областей простору (радіальних і кутових), де знаходиться значна електронна щільність. Насправді можна використовувати плоскі хвильові орбіталі виду
χ(r,θ,ϕ)=Nexp[i(kxrsinθcosϕ+kyrsinθsinϕ+kzrcosθ)],
деN є константа нормалізації іkxky,, іkz є квантовими числами деталізації моментів або довжини хвилі орбіти уздовжx,y, іz декартових напрямків. Перевага використання таких простих орбіталів полягає в тому, що інтеграли, які потрібно виконувати, набагато простіше впоратися з такими функціями. Недоліком є те, що потрібно використовувати багато таких функцій для точного опису різко пікових розподілів зарядів, наприклад, внутрішньооболонкових орбіталей ядра, зберігаючи при цьому достатню гнучкість, щоб також описати набагато більш гладку електронну щільність у валентних областях. Багато зусиль було присвячено розробці та складання таблиць у широко доступних місцях наборів базових орбіталів STO або GTO для елементів основної групи та перехідних металів. Ці постійні зусилля спрямовані на забезпечення стандартних базових бібліотек, які:
- Вихід передбачуваної хімічної точності в результуючих енергіях.
- Економічно ефективні для використання в практичних розрахунках.
- Є відносно переносними, так що основа даного атома є достатньо гнучкою, щоб її можна було використовувати для цього атома в різних середовищах зв'язку (наприклад, гібридизація та ступінь іонізації).
Фундаментальне ядро і валентність основи
При побудові атомно-орбітального базису можна вибрати один з декількох класів функцій. По-перше, повинні бути вказані розмір і характер первинного ядра і валентної основи. У межах цієї категорії поширені наступні варіанти:
- Мінімальна основа, при якій кількість орбіталей CGTO дорівнює числу ядра і валентних атомних орбіталей в атомі.
- Основа подвійної дзети (DZ), в якій використовується вдвічі більше ЦГТО, ніж є основні та валентні атомні орбіталі. Використання більшої кількості базисних функцій мотивовано бажанням забезпечити додаткову варіаційну гнучкість, тому процес LCAO може генерувати молекулярні орбіталі змінної дифузності, оскільки локальна електронегативність атома змінюється. Валентна подвійна дзета (VDZ) основа має лише один CGTO для представлення орбіталів внутрішньої оболонки, але використовує два набори CGTO для опису валентних орбіталів.
- Потрійна дзета (TZ) основа, в якій використовується втричі більше ЦГТО, ніж кількість ядра і валентних атомних орбіталей (звичайно, існують також чотирикратние-дзета і вищі дзета-основи). Більш того, існують бази ВТЗ, які обробляють орбіталі внутрішньої оболонки одним CGTO і валентні орбіталі з трьома CGTO.
Оптимізація орбітальних показників (z або a) і коефіцієнтів стиснення GTO-CGTO для описаного вище виду баз зазнала значного зростання в останні роки. Теорічна група в Тихоокеанських північно-західних національних лабораторіях (PNNL) пропонує всесвітній веб-сайт, з якого можна знайти (і навіть завантажити у формі, підготовленій для введення в будь-який з декількох часто використовуваних електронних кодів структури) широкий спектр гаусових атомних базисних наборів. Цей сайт можна отримати тут. Професор Кірк Петерсон з Університету штату Вашингтон бере участь у проекті розвитку базового набору PNNL, але він також проводить свій власний базовий сайт.
Поляризаційні функції
Один зазвичай підсилює будь-який основний і валентний базовий набір з набором так званих поляризаційних функцій. Вони є функціями одного вищого моменту моменту, ніж з'являється в валентному орбітальному просторі атома (наприклад,d -функції для C, N та O таp -функцій для H), і вони мають показники (zабоa), які змушують їх радіальні розміри бути подібними до розмірів валентних орбіталів ( тобто поляризаційніp орбіталі атома Н схожі за розмірами на1s орбітальну, а не2s валентну орбіталь водню). Таким чином, вони не є орбіталями, які описують валентну орбіталь атома з одним більшим l-значенням; такі орбіталі з вищою валентністю l були б радіально більш дифузними.
Основною метою поляризаційних функцій є надання додаткової кутової гнучкості процесу LCAO при формуванні зв'язкових орбіталей між парами валентних атомних орбіталей. Це проілюстровано на малюнку 6.1.2, де поляризація dp орбіталів на C і O, як видно, сприяє утвореннюp орбітальної зв'язку карбонільної групи, дозволяючи поляризуватиpπ орбіталі атома вуглецю вправо таpπ орбіталі атома кисню вліво.

Функції поляризації мають важливе значення в напружених кільцевих сполуках, таких як циклопропан, оскільки вони забезпечують кутову гнучкість, необхідну для направлення електронної щільності в області між зв'язаними атомами, але вони також важливі в ненапружених сполуках, коли потрібна висока точність.
Дифузні функції
Маючи справу з аніонами або станами Рідберга, потрібно додатково збільшити базовий набір АТ шляхом додавання так званих дифузних базисних орбіталів. Описані вище функції валентності та поляризації не забезпечують достатньої радіальної гнучкості, щоб адекватно описати будь-який із цих випадків. База даних веб-сайту PNNL, наведена вище, пропонує гарне джерело для отримання дифузних функцій, відповідних різним атомам, як це робить сайт професора Кірка Петерсона.
Після визначення атомної орбітальної основи для кожного атома в молекулі, процедура LCAO-MO може бути використана для визначенняχμ,i коефіцієнтів, що описують зайняті та віртуальні (тобто незайняті) орбіталі. Важливо пам'ятати, що базові орбіталі самі по собі не є орбіталями SCF ізольованих атомів; навіть належні атомні орбіталі є комбінаціями (з атомними значеннями дляχμ,i коефіцієнтів) базисних функцій. Сам процес LCAO-MO-SCF визначає величини і ознакиχν,i. Зокрема, саме чергування в знаках цих коефіцієнтів дозволяють утворюватися радіальні вузли.