Методы теории функционала плотности - История развития и достижений квантовой химии

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

Такое положение дел, разумеется, не устраивало исследователей, и в 60-е годы начались работы по изучению электронной плотности (до этого момента использование данной характеристики для описания квантовой системы было скорее интуитивным). Но почему электронная плотность столь привлекательна? Во-первых, она физически определена и измеряема в отличие от волновой функции, не имеющей физического смысла. Во-вторых, получение полной волновой функции электронов -- задача, соответствующая 3n-измерениям (или даже 4n, если принимать во внимание спин), где n-полное число электронов, тогда как электронная плотность есть функция от трех координат независимо от числа электронов в молекуле. Здесь основная сложность заключается в том, существует ли взаимосвязь между электронной плотностью и энергией, и какова форма зависимости, если она существует?

Все началось с работ П. Хоэнберга и В. Кона, которые доказали теорему о свойствах электронной плотности ("Негомогенный электронный газ", 1964 год) [15]. Согласно теореме энергия основного состояния молекулы является функционалом электронной плотности и энергия минимальна, если она (электронная плотность) является точной для основного состояния (функционал -- соответствие числового ряда функции, которая, в свою очередь сопоставлена другому числовому ряду). Следует отметить, что теорема, являя строгое доказательство связи полной энергии с электронной плотностью, не дает никакого правила в построении этого же функционала (о его формах речь пойдет чуть позже).

Методы теории функционала плотности стали использоваться после внедрения в расчетную схему орбиталей, предложенных У. Коном и Л. Дж. Шамом в 1965 году (работа называлась "Самосогласованные уравнения, включающие обменные и корреляционные эффекты"). Основная идея теории состоит в разделении функционала кинетической энергии на две части, первая вычисляется точно с использованием формально построенных орбиталей, отвечающих системе невзаимодействующих электронов, вторая представляет собой поправку. Таким образом, описание молекулярной системы практически полностью соответствует хартри-фоковскому (что, несомненно, следует отнести к достоинству). Но не стоит забывать, что данные орбитали сконструированы лишь так, чтобы давать наилучшее описание электронной плотности, и имеют мало общего с хартри-фоковскими орбиталями.

Спустя некоторое время (конец 70-х годов) было разработано приближение локальной плотности -- простейшая форма теории. Это приближение основывается на том, в каждой точке молекулы с характерной для нее электронной плотностью свойства те же, что и в однородном электронном газе с той же плотностью. Общим случаем локальной плотности является приближение локальной спиновой плотности (LSDA -- Local Spin Density Approximation). Точность методов, основанных на LSDA, сопоставима с методами Хартри-Фока.

Другим более совершенным является приближение градиентной коррекции (GGA -- Generalized Gradient Approximation). В данном случае корреляционный и обменный потенциалы, образующие функционал, связаны не только с электронной плотностью, но и с первой и второй производной от нее. Вначале было предложено несколько функциональных зависимостей в виде поправок к LSDA (например, обменная поправка Бекке, 1988 год). Практически в то же время был предложен ряд функционалов с градиентной коррекцией в чистом виде. Один из наиболее удачных функционалов данной группы был предложен К. Ли, В. Янгом и Р. Г. Паром (LYP, 1988 год).

Третью и наиболее совершенную форму являет группа гибридных методов. Они вобрали в себя лучшее из приближений локальной плотности и градиентной коррекции и получи широкое распространение благодаря феноменально низким временным затратам в сочетании с точностью, которая сопоставима с точностью лучших Ab initio методов (MP4(SDQ), QCISD и CCSD).

В заключение хотелось бы отметить, что методы теории функционала плотности позволили получить ряд электрофизических характеристик (дипольные моменты молекул и внешние электростатические поля), прогнозы спектров в диапазоне от далекой ИК-области до УФ, изучение путей химических реакций с высокой степенью точности. Столь впечатляющие результаты не остались без внимания шведской академии наук, и в 1998 году Нобелевскую премию разделили Джон Попл "за развитие вычислительных методов в квантовой химии" и Уолтер Кон "за развитие теории функционала плотности".

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




Методы теории функционала плотности - История развития и достижений квантовой химии

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