Initio методы - История развития и достижений квантовой химии

В переводе с латинского Ab initio означает "из первых принципов". Действительно, к данной группе относятся методы, в соответствии с которыми вычисление проводится исключительно на теоретической базе, то есть без введения в расчетную схему каких-либо параметров, полученных экспериментальным путем. При расчете все величины имеют конкретный физический смысл. Такими методами являются: метод Хартри-Фока-Рутаана, разнообразные вариации конфигурационного взаимодействия, методы теории возмущения, а также метод объединенных кластеров. К преимуществу данного подхода следует отнести приемлемую точность расчета, относительную универсальность. Недостатком же является ресурсоемкость процедуры, поэтому группа Ab initio методов стала применяться химиками позднее методов полуэмпирических (точнее, они находили применение, но лишь к простейшим системам, не намного превосходившим по сложности молекулу водорода).

По причинам, упомянутым ранее, данная группа расчетных схем долгое время была в забвении. Одним из первых реализован метод Хартри-Фока-Рутаана ("Новые разработки в теории молекулярных орбиталей", 1951 год). Затем двумя годами позже он был доработан Джоном Поплом и Р. К. Несбетом (их работа посвящена радикалам). Хотя этим методом не учитывается коррелированное движение электронов, он явился отправной точкой для так называемых пост-хартри-фоковских методов (методов учета электронной корреляции).

Одной из первых пост-хартри-фоковских схем была теория возмущения Меллера-Плессе. Сама теория возмущения применялась и раньше физиками и математиками, однако к задачам квантовой химии она была адаптирована в 1934 году (в форме концепции), реализована позднее метода Хартри-Фока-Рутаана. Суть подхода заключается в том, что система делится на две части. Первая, упрощенная, точно решаема, тогда как вторая рассматривается в виде возмущения первой. Существенным достоинством данной теории является то, что поправки к полной энергии системы, вычисленные в приближении хартри-фока, находятся неитерационным способом, то есть отсутствует необходимость многократного повторения большого объема расчетов для достижения самосогласованного решения. Дальнейшее развитие теории продолжалось на протяжении второй половины XX-го века и шло в направлении более полного учета корреляционной энергии. Так основные работы по теории возмущения второго порядка (MPPT2), учитывающей 60 -- 90 % корреляционной энергии, относятся к 70-м годам.

Другой способ учета электронной корреляции являет собой теория объединенных кластеров. Если методы теории возмущения заключаются в добавлении к основному хартри-фоковскому решению однократных, двукратных и трехкратных возбуждений, то метод теории объединенных кластеров основан на включении всех поправок данного типа вплоть до бесконечных порядков теории. Теория совершенствовалась с 70-го по 90-й годы. Некоторые разновидности методов теории объединенных кластеров в состоянии учесть до 99 % энергии электронной корреляции при разумных вычислительных затратах. Наиболее известный из разработчиков данного подхода -- Джон Попл.

Третьим и одним из самых поздних подходов является конфигурационное взаимодействие. Его принцип был рассмотрен в предыдущем разделе. Метод "обрел жизнь" в 80-е годы и развивался вплоть до 2000-го года в трудах Джона Попла, Р. Зигера и Р. Кришнана (одна из первых работ -- "Методы вариационного конфигурационного взаимодействия -- сравнение с теорией возмущения", 1977 год). Изначальным недостатком приближения было отсутствие размерной согласованности (совокупная энергия двух молекул, разделенных очень большим расстоянием, не равна сумме энергий этих же молекул, рассчитанных по отдельности) и размерной протяженности (рост погрешности расчета с увеличением числа частиц, входящих в систему). Первый изъян был устранен в 1987 году в работах упомянутых выше ученых добавлением квадратичных членов более высокого порядка (поправки Девидсона). К сожалению, второй недочет для данной группы методов устранен не был [исключение составляет метод полного конфигурационного взаимодействия (FullCI)].

Завершая описание Ab initio методов, следует отметить, что в своих наиболее совершенных (и затратных) формах выражения для волновой функции описанных выше трех пост-хартри-фоковских приближений практически совпадают.

До последнего момента рассматривались методики решения уравнения Шредингера, но ничего не говорилось о форме волновой функции. Волновую функцию было принято рассматривать в качестве линейной комбинации достаточно простых функций (например, функции Гаусса), называемой базисом или базисным набором. Каждая функция характеризуется некоторым числом подгоночных коэффициентов, обеспечивающих гибкость базиса. В соответствии с вариационным принципом, чем больше базисных функций, тем более точные решения могут быть получены. Но не стоит забывать о временных затратах: чем больше базисных функций (чем более гибкий базисный набор), тем больше времени потребуется для решения уравнения Шредингера.

Одним из первых появился минимальный базисный набор. Он был разработан под началом Р. Ф. Стюарта и Джона Попла в 1969 году. В данном базисе осуществляется представление атомных орбиталей слетеровского типа в виде комбинации гауссовых функций. Наиболее популярным в свое время был минимальный базисный набор STO-3G. Данный базис достаточно экономичен, но обладает одним существенным недостатком -- он очень жесткий, не способен подстраивать свой размер в зависимости от окружения атома, и характеризуется практически сферическим распределением заряда. В свете последнего погрешность расчета закономерно увеличивается с усложнением электронной структуры.

Более совершенными являются валентно-расщепленные базисные наборы. Они представляют атомную орбиталь в виде двух валентных функций одинаковой симметрии. Одна из них является более сжатой, а другая -- диффузной. Наиболее распространенными базисами этого типа являются 3-21G и 6-31G. Валентно-расщепленные базисные наборы были описаны в Работе Дж. С. Бинкли, Джона Попла и В. Дж. Гера 1980 года "Базисные наборы малого валентного расщепления для элементов второго периода".

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

Среди поляризационных базисных наборов получил распространение 6-31G(d).

Описанные выше базисные наборы хорошо моделируют электронейтральные системы, то есть молекулы, у которых электроны прочно удерживаются ядрами. В анионах же электрон очень слабо связан с ядром, что проявляется в значительном удалении электронной плотности от него. По этой причине в базисный набор включают диффузные функции s - и p-типа с малыми значениями экспоненциальных коэффициентов, что обусловливает большой размер и удаленность этих функций от ядра. Такие функции обозначают символом "+". В качестве примера можно привести базисный набор, включающий как поляризационные, так и диффузные функции -- 6-31+G(d).

Похожие статьи




Initio методы - История развития и достижений квантовой химии

Предыдущая | Следующая