Теория квантовой химии

ТЕОРИЯ КВАНТОВОЙ 

ХИМИИ

В.В. КОСТЮКОВ

Рекомендовано 

Межрегиональным учебно-методическим советом 

профессионального образования в качестве учебного пособия 

для студентов высших учебных заведений, обучающихся 

по химическим и биологическим направлениям подготовки 

(квалификация (степень) «бакалавр») 

(протокол № 6 от 16.06.2021)

Москва
ИНФРА-М

202УЧЕБНОЕ ПОСОБИЕ

УДК 544.18(075.8)
ББК 24.5я73
 
К72

Р е ц е н з е н т ы:

Л.А. Яковишин, доктор химических наук, доцент, профессор ка-

федры химии и химических технологий Севастопольского государ-
ственного университета;

Э.В. Ткаченко, кандидат химических наук, доцент кафедры физики 

и общетехнических дисциплин Черноморского высшего военно-мор-
ского училища имени П.С. Нахимова

ISBN 978-5-16-016258-4 (print)
ISBN 978-5-16-108567-7 (online)
© Костюков В.В., 2021

Костюков В.В.

К72  
Теория квантовой химии : учебное пособие / В.В. Костюков. — 

Москва : ИНФРА-М, 2023. — 236 с. — (Высшее образование: Бакалав-
риат). — DOI 10.12737/1090584.

ISBN 978-5-16-016258-4 (print)
ISBN 978-5-16-108567-7 (online)
В учебном пособии кратко изложены основные теории квантовой хи-

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

Соответствует требованиям федеральных государственных образова-

тельных стандартов высшего образования последнего поколения.

Предназначено для студентов бакалавриата высших учебных заведе-

ний; может быть использовано аспирантами, исследующими вопросы ма-
териаловедения, структурной, органической и физической химии, молеку-
лярной биологии и биофизики, биотехнологии.

УДК 544.18(075.8)

ББК 24.5я73

Список сокращений и условных обозначений

CASSCF — метод полного самосогласованного поля в активном 
пространстве (Complete Active Space Self-Consistent Field)
CC — теория сопряженного кластера (Coupled Cluster)
CI — теория конфигурационного взаимодействия (Configurational 
Interaction)
CNDO — полное пренебрежение дифференциальным перекры-
ванием (Complete Neglect of Differential Overlap)
CSF — функция конфигурационного состояния (Configurational 
State Functions)
DBOC — диаго нальная поправка Борна — Оппенгеймера 
(Diagonal Born-Oppenheimer Correction)
DFT — теория функцио нала плотности (Density Functional 
Theory)
EHT — расширенная теория Хюккеля (Extended Huckel Theory)
FFT — быстрое преобразование Фурье (Fast Fourier Transform)
FMM — метод быстрых мультиполей (Fast Multipole Method)
GGA — приближение обобщенной градиентной поправки 
(Generalized Gradient Approximation)
GVB — теория обобщенной валентной связи (Generalized 
Valence Bond)
HF — теория Хартри — Фока (Hartree-Fock)
HOMO — наивысшая занятая молекулярная орбиталь (Highest 
Occupied Molecular Orbital)
INDO — промежуточное пренебрежение дифференциальным 
перекрыванием (Intermediate Neglect of Differential Overlap)
KS — теория Кона — Шама (Kohn-Sham)
LDA — приближение локальной плотности (Local Density 
Approximation)
LSDA — приближение локальной спиновой плотности (Local 
Spin Density Approximation)
LUMO — наинизшая незанятая (виртуальная) молекулярная 
орбиталь (Lowest Unoccupied Molecular Orbital)
MBPT — теория возмущений многих тел (Many-Body 
Perturbation Theory)
MCSCF — теория мультиконфигурационного самосогласован-
ного поля (Multi-Configuration Self-Consistent Field)
MINDO — модифицированное промежуточное пренебрежение 
дифференциальным перекрыванием (Modified Intermediate Neglect 
of Differential Overlap)

MNDO — модифицированное пренебрежение межатомным пе-
рекрыванием (Modified Neglect of Diatomic Overlap)
МР — теория возмущений Мёллера — Плессе
MRCI — многоопорное конфигурационное взаимодействие 
(Multi-Reference Configuration Interaction)
NDDO — пренебрежение межатомным дифференциальным пе-
рекрыванием (Neglect of Diatomic Differential Overlap)
PP — идеальное спаривание электронов (Perfect Pairing)
RASSCF — метод самосогласованного поля с ограниченным 
активным пространством (Restricted Active Space Self-Consistent 
Field)
SAC — адаптированная к спину конфигурация (Spin-Adapted 
Configurations)
SCF — самосогласованное поле (Self-Consistent Field)
SCVB — теория спин-связанных валентных связей (Spin-
Coupled Valence Bond)
SO — условие строгой ортого нальности волновой функции 
(Strong Orthogonality)
VB — теория валентных связей (Valence Bond)
ZDO — нулевое дифференциальное перекрывание (Zero 
Differential Overlap)
аi, сi — линейные коэффициенты разложения
С — матрица коэффициентов молекулярных орбиталей
D — матрица электронной плотности
ˆD — оператор Дирака
Е — полная энергия системы
F — матрица Фока
ˆF  — оператор Фока
ge — множитель Ланде (гиромагнитный множитель, g-фактор) 
для электрона
ˆg — оператор межэлектронного отталкивания
ˆH, ˆh — операторы Гамильтона (ОГ)
J — кулоновский интеграл
K — обменный интеграл
L — функция Лагранжа
Mbasis — число базисных функций в базисном наборе
Ne — число электронов в системе
Nn — число ядер в системе
ˆp — оператор импульса
r— радиус-вектор электрона
R
— радиус-вектор ядра

ˆs — оператор спина
S — интеграл перекрывания
ˆT — оператор кинетической энергии
ˆV  — оператор потенциальной энергии
α, β — спиновые переменные
εi, λij — множитель Лагранжа
λ — параметр возмущения
μВ — магнетон Бора
ρ — электронная плотность
τ — затраты машинного времени
Ψ, ϕ — волновая функция
χ — атомная орбиталь
а.е.Х. — атомная единица Хартри
АО — атомная орбиталь
БН — базисный набор
БФ — базисная функция
ВП — вариационный принцип
ВФ — волновая функция
ИСО — инерциальная система отсчета
КМК — квантовый метод Монте-Карло
КЭД — квантовая электродинамика
ММ — молекулярная механика
МО — молекулярная орбиталь
ОГ — оператор Гамильтона
ОС — определитель Слейтера
ПБО — приближение Борна — Оппенгеймера
ППЭ — поверхность потенциальной энергии
ПС — переходная структура
ПЭМ — полуэмпирический метод
САО — среднее абсолютное отклонение
СКО — среднеквадратичное отклонение
СОВ — спин-орбитальное взаимодействие
СТО — специальная теория относительности
ТХК — теорема Хоэнберга — Кона
УШ — уравнение Шредингера
ЭК — электронная корреляция

Введение

Квантовохимические расчеты сегодня являются неотъемлемой 
частью научных исследований не только в химии, но и в матери-
аловедении, молекулярной биологии, биофизике, биотехнологии 
и синтезе лекарственных препаратов. Практически любые научные 
изыскания, связанные со строением молекул и вещества в целом, 
обязательно включают квантовохимическое компьютерное моде-
лирование. Успехи в этой области были отмечены в 1998 г. прису-
ждением Нобелевской премии Вальтеру Кону и Джону Поплу — 
основателям современной компьютерной квантовой химии. Данная 
награда явилась признанием более чем 30-летней работы ученых. 
В 1964 г. В. Кон сформулировал теорию функцио нала электронной 
плотности, сделавшую реальным ранее невозможный в силу своей 
чрезвычайной сложности расчет многих молекулярных систем. 
В 1970 г. Д. Попл создал первый вариант компьютерной прог-
раммы для таких расчетов — GAUSSIAN, современная версия ко-
торого и на сегодняшний день остается лучшей из пакетов вычи-
слительной квантовой химии.
В настоящем пособии кратко изложены сложные, но инте-
ресные теории, возникшие на стыке квантовой физики и струк-
турной химии. Поскольку сегодня все квантовохимические рас-
четы выполняются на компьютере, уже на этапе создания этих 
теорий их авторам пришлось решать основную задачу — выбор 
соотношения «ресурсоемкость — точность». Поэтому современные 
теории квантовой химии всегда ориентированы на создание на их 
основе программных вычислительных пакетов. Впечатляющая эф-
фективность квантовохимических методов обусловлена не только 
многолетней упорной работой теоретиков, но и бурным развитием 
компьютерной техники (как производительности процессоров, так 
и емкости памяти и скоростей ее чтения (записи)). Именно появ-
ление суперкомпьютеров, обладающих огромной вычислительной 
мощностью, позволило выполнить расчеты, что ранее считалось 
невозможными ввиду их гигантской ресурсоемкости (например, 
моделирование ab initio крупных белковых молекул).
Ключевой проблемой квантовой химии является учет элект-
ронной корреляции — взаимного влияния движущихся электронов. 
Факторами, сильно усложняющими анализ и моделирование кор-
реляции электронов, являются их волновые свойства, а также на-
личие у них заряда и спина, приводящих к обменным, спин-ор-

битальным и различным электромагнитным взаимодействиям. 
Поэтому в каждой из описанных в настоящем пособии теорий зна-
чительное внимание уделено именно электронной корреляции.
Немаловажным фактором, влияющим на квантовохимические 
свойства атомов и молекул, являются релятивистские эффекты, 
которым в пособии посвящен отдельный параграф. Огромные ско-
рости движения электронов и, наоборот, очень малые расстояния 
между ними и ядрами в атоме приводят к существенному влиянию 
на электрон его собственного электромагнитного поля и к другим 
чисто релятивистским явлениям.
Современная квантовая химия — мощная, бурно развива-
ющаяся область знаний, владение которой химику, физику, био-
логу или фармацевту сегодня особенно необходимо. Цель данной 
книги — ознакомление читателя с ее теоретической базой и нераз-
рывно связанными с ней практическими (вычислительными) ас-
пектами.
После изучения данного пособия студенты должны:
 
• знать и понимать положения и следствия основных теорий 
квантовой химии;
 
• уметь интерпретировать и анализировать квантовохимические 
эффекты и закономерности;
 
• владеть навыками выбора и применения квантовохимических 
методов при решении вычислительных задач для различных 
молекулярных систем.

Глава 1. 
ТЕОРИЯ ХАРТРИ — ФОКА

1.1. КРАТКАЯ ХАРАКТЕРИСТИКА КВАНТОВОХИМИЧЕСКИХ 
ТЕОРИЙ

Единственным инструментом, используемым для детального 
описания распределения электронов, на сегодняшний день яв-
ляется квантовая механика. Рассмотрим решение стационарного 
(не зависящего от времени) уравнения Шредингера (УШ), которое 
в краткой операторной форме имеет следующий вид:

 
ˆH
E
Ψ =
Ψ. 
(1.1)

Существенной частью решения УШ является приближение 
Борна — Оппенгеймера (ПБО), в котором взаимосвязью между 
движениями ядра и электронов пренебрегают. Это позволяет ре-
шить электронную часть УШ с координатами ядра в качестве па-
раметров. Получаемая в результате поверхность потенциальной 
энергии (ППЭ) создает основу для решения задачи движения ядра. 
Поэтому бóльшую часть вычислений составляет решение элект-
ронного УШ при фиксированных координатах ядра. В теории кван-
товой химии можно условно выделить шесть основных разделов, 
изложенных в настоящем пособии.
1. Динамика многоэлектронной системы очень сложна и, сле-
довательно, требует сложных вычислительных методов для опи-
сания. Значительное упрощение как концептуально, так и в вычи-
слительном отношении может быть достигнуто путем введения мо-
делей независимых частиц, в которых движение данного электрона 
считается независимым от координат остальных; пренебрегают 
большинством взаимодействий между частицами, кроме самых 
важных, или учитывают все взаимодействия, но в усредненной 
форме. Оказалось, что только второй подход (теория Хартри — 
Фока (Hartree-Fock, HF)) имеет приемлемую точность. В HF 
каждый электрон описывается индивидуальной волновой функ-
цией (ВФ) — атомной орбиталью (АО), а полная ВФ Ψ задается 
как их произведение.
Так как электроны являются неразличимыми фермионами (ча-
стицами со спином 1/2), Ψ должна быть антисимметричной (из-

менять знак при перестановке любых двух электронов), что до-
стигается путем ее представления в виде определителя Слейтера 
(ОС, Slater). Оптимальный набор АО определяется согласно вари-
ационному принципу (ВП), т.е. они должны иметь самую низкую 
энергию в пределах ограничений однодетерминантной (заданной 
одним ОС) ВФ.
Молекулярная орбиталь (МО) учитывает притяжение ко всем 
ядрам и среднее отталкивание от всех остальных электронов. По-
скольку другие электроны описываются своими АО, уравнения HF 
зависят от их собственных решений и по это му должны решаться 
итеративно. МО записывается при помощи базисного набора (БН), 
а уравнения HF могут быть зафиксированы как матричная задача 
на собственные значения. Элемен ты в матрице Фока соответствуют 
интегралам одно- и двухэлектронных операторов по базисным 
функциям (БФ), умноженным на элемен ты матрицы плотности. 
Таким образом, уравнения HF могут быть получены путем по-
вторных диаго нализаций матрицы Фока. К HF можно добавлять 
дополнительные приближения (упрощения), приводящие к по-
луэмпирическим методам (ПЭМ). Также ее можно улучшить, до-
бавляя дополнительные ОС и тем самым создавая модели, которые 
сходятся к точному решению электронного УШ (см. (1.1); рис. 1.1).
2. ПЭМ вытекают из HF в результате отказа от вычисления 
при построении матрицы Фока всех интегралов, включающих 
более двух ядер. Поскольку HF обладает ограниченной точностью, 
сами по себе приближения ПЭМ могут привести к существенным 
ошибкам. Поэтому успех ПЭМ зависит от замены данных интег-
ралов параметрами, подогнанными к экспериментальным данным, 
особенно к энергии и геометрии. Такие методы в вычислительном 
отношении намного более эффективны, чем HF, но ограничены си-
стемами, для которых указанные параметры известны.

ˆH
E
Ψ =
Ψ

Уравнения 
Хартри — Фока

Полуэмпирические 
методы

Сходимость к точ-
ному решению

Добавление большего 
числа определителей
Дополнительные 
приближения

Рис. 1.1. Модель Хартри — Фока как исходная для более точных методов

3. HF учитывает только среднее межэлектронное взаимодей-
ствие и, следовательно, не учитывает корреляцию между дви-
жениями электронов. Методы, которые включают электронную 
корреляцию (ЭК), требуют многодетерминантной ВФ, поскольку 
лучшей однодетерминантной ВФ является хартри-фоковская. 
Многодетерминантные методы в вычислительном отношении на-
много сложнее HF, но могут давать результаты, систематически 
приближающиеся к точному решению УШ.
4. Теория функцио нала плотности (Density Functional Theory, 
DFT) Кона и Шама может рассматриваться как улучшение HF, 
в которой многочастичный эффект ЭК моделируется функцией 
электронной плотности. DFT так же, как HF, является моделью 
независимых частиц и сравним с HF в вычислительном отношении, 
но обеспечивает значительно лучшие результаты. Основным недо-
статком DFT является отсутствие систематического подхода к по-
лучению более точных решений.
5. Хотя в настоящее время в квантовой химии преобладает 
теория МО Малликена, исторически предшествующая ей теория 
валентных связей (Valence Bond, VB) Полинга позволяет наглядно 
интерпретировать объединение атомов в молекулы. VB основана 
на представлении о том, что каждая пара атомов в молекуле удер-
живается вместе при помощи одной или нескольких общих элект-
ронных пар. При образовании молекулы из атомов последние 
в значительной мере сохраняют свою электронную конфигурацию, 
а связывание атомов достигается в результате обмена электронов 
между ними и спаривания спинов двух электронов, находящихся 
на АО исходных атомов. Ключевыми концепциями VB являются 
гибридизация АО и резонанс канонических структур.
6. Пренебрежение релятивистскими эффектами оправдано 
для первых трех периодов таблицы Менделеева (т.е. Z < 40), если 
не требуется высокая точность. Однако данные эффекты стано-
вятся важными, начиная с элемен тов четвертого периода, а также 
для переходных металлов. Релятивистское происхождение имеют 
спин-зависимые эффекты (например, спин-орбитальное взаимо-
действие), хотя они могут быть представлены как поправки к ре-
шению электронного УШ (например, в рамках теории возмущений; 
см. параграф 3.4).
Для упрощения записи уравнений в настоящем пособии иногда 
будут использоваться обозначения Дирака:

*
*
*
*
ˆ
ˆ
;  
;  
;  
.
dr
H
dr
H
Ψ ≡ Ψ
Ψ ≡ Ψ
Ψ Ψ
≡ Ψ Ψ
Ψ
Ψ
≡ Ψ
Ψ
∫
∫

(1.2)