B3LYP

B3LYP — это гибридный функционал плотности, один из наиболее широко используемых и успешных приближений в рамках теории функционала плотности (DFT) в квантовой химии. Он сочетает в себе обменный функционал, вычисленный по методу Хартри — Фока, с обменным и корреляционным функционалами, полученными из теории возмущений и приближения локальной плотности. Разработанный в начале 1990-х годов, B3LYP стал стандартом де-факто для расчётов электронной структуры молекул, термохимии, кинетики и спектроскопии, обеспечивая хороший баланс между точностью и вычислительной стоимостью.

История и происхождение

Разработка функционала B3LYP связана с эволюцией теории функционала плотности. В 1988 году Аксель Бекке (Axel Becke) предложил градиентно-скорректированный обменный функционал (B88), который улучшил описание обменной энергии по сравнению с приближением локальной плотности (LDA). В 1993 году Бекке представил концепцию гибридных функционалов, смешивающих точный хартри-фоковский обмен с DFT-обменом. Его трёхпараметрическая формула (B3) была параметризована для минимизации ошибок в термохимических данных.

Одновременно с этим, Ченгте Ли, Вейцзян Янг и Роберт Парр (LYP) разработали корреляционный функционал, основанный на выражении для корреляционной энергии Колле — Сальветти. В 1988 году Ли, Янг и Парр опубликовали свой функционал (LYP), который хорошо описывает динамическую корреляцию электронов.

Комбинация обменного функционала Бекке (B3) с корреляционным функционалом LYP и привела к созданию B3LYP. Формально, функционал B3LYP не был выведен из первых принципов, а был эмпирически подогнан к набору экспериментальных данных. Несмотря на это, он оказался чрезвычайно эффективным для широкого круга задач, и его популярность стремительно возросла.

Формальное определение и параметризация

Функционал B3LYP определяется как смесь нескольких компонентов:

\[ E_{xc}^{B3LYP} = (1 - a) E_x^{LDA} + a E_x^{HF} + b \Delta E_x^{B88} + (1 - c) E_c^{LDA} + c E_c^{LYP} \]

где:

• \( E_x^{LDA} \) — обменная энергия в приближении локальной плотности (по Дираку).
• \( E_x^{HF} \) — точный обмен Хартри — Фока.
• \( \Delta E_x^{B88} \) — градиентная поправка к обмену по Бекке (B88).
• \( E_c^{LDA} \) — корреляционная энергия в приближении локальной плотности (по Восько — Вилку — Нусаиру, VWN).
• \( E_c^{LYP} \) — корреляционная энергия по Ли — Янгу — Парру.

Параметры \( a \), \( b \) и \( c \) были определены Бекке путём подгонки к экспериментальным данным. Оригинальные значения, использованные в программе Gaussian, составляют:

• \( a = 0.20 \) (доля точного обмена)
• \( b = 0.72 \) (вклад градиентной поправки к обмену)
• \( c = 0.81 \) (вклад корреляции LYP)

Важно отметить, что в стандартной реализации B3LYP (например, в Gaussian) корреляционная часть LDA заменяется не на чистый VWN, а на его модификацию (VWN5), что может приводить к небольшим различиям в результатах между разными программными пакетами.

Характеристики и свойства

Баланс точности и скорости

B3LYP обеспечивает высокую точность для термохимических расчётов (энергии реакций, энергии активации, теплоты образования) со средней абсолютной ошибкой порядка 2-5 ккал/моль для органических молекул. Это делает его значительно точнее, чем чистые функционалы LDA или GGA, и при этом намного быстрее, чем методы пост-Хартри-Фока, такие как MP2 или CCSD(T).

Применимость

Функционал хорошо работает для большинства органических и элементоорганических соединений, включая молекулы с ковалентными связями. Он успешно описывает геометрию молекул, частоты колебаний (ИК- и КР-спектры), а также барьеры внутреннего вращения. B3LYP также широко используется для расчёта химических сдвигов в ЯМР-спектроскопии.

Ограничения

Несмотря на свою популярность, B3LYP имеет ряд известных недостатков:

• Дисперсионные взаимодействия: B3LYP не учитывает дисперсионные силы (вандерваальсовы взаимодействия), что приводит к систематическим ошибкам при описании слабосвязанных комплексов, молекулярных кристаллов и биологических макромолекул. Для исправления этого используются эмпирические дисперсионные поправки, такие как D3 (Grimme) или D4.
• Переходные металлы: Для систем с переходными металлами, особенно с открытой оболочкой, B3LYP может давать неточные результаты. Часто требуются функционалы с более высокой долей точного обмена (например, B3LYP* или M06) или специальные методы.
• Длинновозбуждённые состояния: B3LYP плохо описывает состояния с переносом заряда и Rydberg-состояния. Для таких задач используются функционалы с дальнодействующей коррекцией (например, CAM-B3LYP, wB97X-D).
• Спиновые состояния: Для мультиспиновых систем (например, комплексов железа) B3LYP может неправильно предсказывать относительную энергию различных спиновых состояний.

Применение

B3LYP является одним из наиболее востребованных инструментов в вычислительной химии. Он используется в следующих областях:

• Органическая химия: расчёт механизмов реакций, поиск переходных состояний, моделирование каталитических циклов.
• Неорганическая химия: изучение структуры и свойств координационных соединений, кластеров, наночастиц.
• Биохимия: моделирование активных центров ферментов, расчёт спектров биомолекул (например, флуоресцентных белков).
• Материаловедение: расчёт электронных свойств полимеров, полупроводников, фотоактивных материалов.
• Спектроскопия: предсказание ИК-, КР-, УФ-видимой и ЯМР-спектров для идентификации веществ.

Критика и альтернативы

Основная критика в адрес B3LYP связана с его эмпирическим происхождением и отсутствием строгого теоретического обоснования. С развитием вычислительных мощностей и появлением новых функционалов (например, PBE0, M06-2X, wB97X-D, MN15) его доминирование несколько снизилось. Эти более современные функционалы часто превосходят B3LYP в конкретных задачах (например, в описании дисперсионных взаимодействий или кинетики). Тем не менее, B3LYP остаётся «рабочей лошадкой» квантовой химии благодаря своей универсальности, надёжности и огромному массиву накопленных данных.

Интересные факты

• B3LYP является одним из наиболее цитируемых функционалов в научной литературе. По состоянию на 2020-е годы, количество ссылок на оригинальные статьи, описывающие его, превышает 100 000.
• В некоторых программных пакетах (например, в Gaussian) реализация B3LYP использует не оригинальный корреляционный функционал VWN, а его модификацию VWN5, что может приводить к небольшим различиям в результатах между пакетами.
• Существует модификация B3LYP*, в которой доля точного обмена уменьшена с 20% до 15%, что улучшает описание спиновых состояний комплексов переходных металлов.

Источники

1. Becke, A. D. (1993). Density-functional thermochemistry. III. The role of exact exchange. The Journal of Chemical Physics, 98(7), 5648–5652.
2. Lee, C., Yang, W., & Parr, R. G. (1988). Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Physical Review B, 37(2), 785–789.
3. Stephens, P. J., Devlin, F. J., Chabalowski, C. F., & Frisch, M. J. (1994). Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields. The Journal of Physical Chemistry, 98(45), 11623–11627.
4. Grimme, S. (2006). Semiempirical GGA-type density functional constructed with a long-range dispersion correction. Journal of Computational Chemistry, 27(15), 1787–1799.
5. Cohen, A. J., Mori-Sánchez, P., & Yang, W. (2012). Challenges for Density Functional Theory. Chemical Reviews, 112(1), 289–320.