2025-11-16T18:52:11.890181

Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods

Masteran, Gaudel, Valeev
By correlating only the 1-particle states occupied in the reference determinant the conventional design for the single-reference R12/F12 explicitly-correlated methods biases them towards the ground state description thereby making treatment of response properties of the ground state, and energies and other properties of excited states less robust. While the use of multireference methods and/or extensions of the standard SP projected geminals can achieve a more balanced description of ground and excited states, here we show that the same goals can be achieved by extending the action of F12 correlators to the occupied and valence unoccupied 1-particle states only. This design choice reflects the strong dependence of the optimal correlation lengthscale of the F12 ansatz on the orbital energies/structure, and helps to avoid the unphysical raising of the ground state energy if the F12 geminals are used to correlate pairs of all 1-particle states. The improved F12 geminal design is incorporated into the unitary transcorrelation framework to produce a unitary 2-body Hamiltonian that incorporates the short-range dynamical correlation physics for ground and low-energy excited states in a balanced manner. This explicitly-correlated effective Hamiltonian reduces the basis set requirement on the correlation-consistent basis cardinal number by 1 or more over the uncorrelated counterpart for the description of the ground state coupled-cluster singles and doubles (CCSD) energies, the vertical excitation energies and harmonic vibrational frequencies of equation-of-motion CCSD low-energy excited states.
academic

К сбалансированному описанию основного и возбуждённых состояний методами трансворрелированного F12

Основная информация

  • ID статьи: 2508.20233
  • Название: Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods
  • Авторы: Conner Masteran, Bimal Gaudel, Edward F. Valeev (Virginia Tech)
  • Классификация: physics.chem-ph
  • Дата публикации: 29 августа 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2508.20233

Аннотация

Традиционные методы явной корреляции с одним опорным состоянием R12/F12 коррелируют только занятые одночастичные орбитали в опорном определителе, что приводит к смещению в сторону описания основного состояния и делает обработку свойств отклика основного состояния, а также энергий и других свойств возбуждённых состояний недостаточно надёжной. В данной работе показано, что сбалансированное описание основного и возбуждённых состояний может быть достигнуто путём расширения действия коррелятора F12 на занятые орбитали и виртуальные одночастичные состояния валентного слоя. Такая конструкция отражает сильную зависимость оптимальной длины корреляции F12 ansatz от энергии орбиталей/структуры и помогает избежать нефизического повышения энергии основного состояния при использовании F12 geminals для корреляции всех пар одночастичных состояний. Улучшенная конструкция F12 geminal интегрирована в унитарную трансворрелированную схему, что даёт унитарный двухчастичный гамильтониан, содержащий короткодействующую динамическую корреляцию основного состояния и низколежащих возбуждённых состояний сбалансированным образом.

Научный контекст и мотивация

Проблемный фон

  1. Значимость методов явной корреляции: Методы F12 значительно ускоряют сходимость по базисному набору путём явного включения межэлектронного расстояния в волновую функцию, эффективно моделируя известное аналитическое поведение волновой функции вблизи электрон-электронной кулисы.
  2. Ограничения традиционных методов F12:
    • Подавляющее большинство разработок и применений методов F12 сосредоточено на свойствах основного состояния
    • Традиционный ansatz II (коррелирующий только пары занятых орбиталей) смещён в описании возбуждённых состояний
    • Приводит к искусственному завышению вертикальных энергий возбуждения, как впервые наблюдали Fliegl и др.
  3. Недостатки существующих решений:
    • Многоопорные методы, хотя и могут равномерно обрабатывать несколько электронных состояний, имеют высокую вычислительную стоимость
    • Расширенный SP ansatz требует дополнительной оптимизации параметров, что может привести к численной нестабильности

Научная мотивация

Основная гипотеза данной работы: однопараметрический экспоненциальный коррелятор оптимизирован в основном для валентных орбиталей и не подходит для описания пар орбиталей, включающих невалентные орбитали. Эта точка зрения подтверждается исследованиями ядерной корреляции, где Werner и др. обнаружили значительные различия в оптимальных показателях экспоненты коррелятора для ядро-ядро, ядро-валентных и валентно-валентных пар орбиталей.

Основные вклады

  1. Предложен метод валентно-ограниченного geminal: F12 коррелятор ограничен подмножеством виртуальных орбиталей валентного слоя, избегая нефизического повышения энергии основного состояния в традиционных методах
  2. Разработаны три новых ansatz коррелятора:
    • ansatz II: традиционный метод, коррелирующий только пары занятых орбиталей
    • ansatz IP: гибридный метод, коррелирующий пары занято-виртуальных орбиталей
    • ansatz PP: полный метод, коррелирующий все пары активных орбиталей
  3. Реализовано сбалансированное описание основного состояния-возбуждённых состояний: коррелятор VVO PP сохраняет точность основного состояния при значительном улучшении описания возбуждённых состояний
  4. Доказано значительное улучшение сходимости по базисному набору: потребность в базисном наборе для коррелирующей энергии и энергии возбуждения уменьшена на один или более порядков

Детальное описание методики

Определение задачи

Данное исследование направлено на разработку методов F12 трансворрелирования, способных сбалансированно описывать основное и возбуждённые состояния молекул. Входные данные: геометрия молекулы; выходные данные: высокоточные энергии и свойства основного и возбуждённых состояний.

Архитектура основного метода

1. Унитарная F12 трансворрелированная схема

Трансворрелированный гамильтониан генерируется подобным преобразованием:

Ĥ̄ = e^(-Â)Ĥe^Â

где Â — оператор, содержащий зависимость от межчастичного расстояния.

2. Три ansatz коррелятора

Ansatz II (традиционный метод):

Â_II = 1/2(G^(α₁α₂)_(i₁i₂) Ê^(i₁i₂)_(α₁α₂) - h.c.)

Ansatz IP (гибридный метод):

Â_IP = (G^(α₁α₂)_(i₁p'₁) Ê^(i₁p'₁)_(α₁α₂) - h.c.) - Â_II

Ansatz PP (полный метод):

Â_PP = 1/2(G^(α₁α₂)_(p'₁p'₂) Ê^(p'₁p'₂)_(α₁α₂) - h.c.)

3. Конструкция виртуальных орбиталей валентного слоя (VVO)

Виртуальные орбитали валентного слоя определяются посредством проекции на минимальный базис:

  1. Выполняется симметричная ортогонализация Лёвдина атомных орбиталей минимального базиса
  2. Вычисляется матрица перекрытия с неоккупированными орбиталями
  3. Виртуальные орбитали валентного слоя получаются через SVD-разложение

Технические инновации

  1. Стратегия валентного ограничения: впервые F12 geminals ограничены подпространством валентных орбиталей на основе физических соображений, связанных с масштабами энергии орбиталей
  2. Автоматизированная реализация: использован пакет символической тензорной алгебры SeQuant для автоматического генерирования и оптимизации тензорных выражений
  3. Универсальный по состояниям дизайн: ansatz PP обеспечивает лучшее описание низколежащих возбуждённых состояний при сохранении точности основного состояния

Экспериментальная установка

Тестовые молекулярные системы

Для детального исследования выбраны пять малых молекул: BH, CH₂O, N₂, H₂O, CH₂, охватывающих различные типы характеристик возбуждённых состояний.

Выбор базисного набора

  • aug-cc-pVXZ (aXZ)
  • doubly-augmented daug-cc-pVXZ (daXZ)
  • cc-pVXZ-F12 (XZ-F12)
  • Использованы рекомендуемые значения показателя geminal

Вычислительные методы

  • Основное состояние: полная переоптимизация орбиталей Хартри-Фока + CCSD
  • Возбуждённые состояния: метод EOM-CCSD с замороженным ядром
  • Оптимизация геометрии и вибрационные частоты получены посредством подгонки поверхности потенциальной энергии

Экспериментальные результаты

Основные результаты

1. Улучшение вертикальных энергий возбуждения

На примере молекулы BH (таблица 1):

  • Коррелятор II завышает энергию возбуждения примерно на 0,1 эВ
  • Коррелятор VVO PP на уровне aTZ снижает ошибку энергии возбуждения с 0,01 эВ до 0,004 эВ
  • Восстановление коррелирующей энергии основного состояния близко к уровню ansatz II

2. Производительность для различных молекулярных систем

  • CH₂O: коррелятор VVO PP значительно улучшает энергию возбуждения Ридберговского состояния, снижая ошибку на уровне DZ с 0,2-0,3 эВ до меньших значений
  • N₂: коррелятор II завышает энергию возбуждения на 0,5 эВ, корреляторы VVO IP/PP значительно снижают ошибку
  • H₂O: для энергий возбуждения, сходящихся снизу, трансворрелирование F12 значительно улучшает сходимость по базисному набору

3. Улучшение свойств возбуждённых состояний

Сходимость по базисному набору длин связей и вибрационных частот улучшена в среднем на один порядок базиса:

  • Равновесные длины связей: результаты метода F12 практически совпадают с традиционными результатами EOM-CCSD(F12)
  • Вибрационные частоты: {средняя, максимальная} абсолютная ошибка коррелятора VVO PP составляет {2,1, 5} см⁻¹

Абляционные эксперименты

Сравнение FVO и VVO доказывает важность валентного ограничения:

  • Вариант FVO систематически производит меньшую коррелирующую энергию основного состояния
  • Вариант VVO достигает лучшего баланса между основным и возбуждёнными состояниями

Связанные работы

Данная работа опирается на следующие важные исследования:

  1. Теоретические основы F12: SP ansatz Ten-no и стандартные проекционные geminals
  2. Методы F12 для возбуждённых состояний: XSP ansatz Кёна, EOM-CC-F12 Бохана и Ten-no
  3. Трансворрелированные методы: унитарная трансворрелированная схема Янаи-Сиозаки
  4. Теория валентных орбиталей: определение виртуальных орбиталей валентного слоя Лу и др., Шмидта и др.

Заключение и обсуждение

Основные выводы

  1. Оптимальность коррелятора VVO PP: достигает наилучшего баланса между восстановлением коррелирующей энергии основного состояния и описанием возбуждённых состояний
  2. Значительное улучшение сходимости по базисному набору: сходимость коррелирующей энергии по базисному набору улучшена на два порядка базиса
  3. Универсальность метода: применим к различным типам возбуждённых состояний (валентные состояния и состояния Ридберга)

Ограничения

  1. Неполнота одночастичного базиса: для некоторых возбуждённых состояний остаётся значительная ошибка одночастичного базиса
  2. Оптимизация параметров geminal: в настоящее время используются фиксированные параметры, которые могут быть неоптимальны для всех систем
  3. Ограничение по размеру системы: в настоящее время проверено только на малых молекулах

Направления будущих исследований

  1. Обработка неполноты одночастичного базиса: разработка универсальных по состояниям методов коррекции одночастичного базиса
  2. Расширение на релятивистские эффекты: расширение метода на включение специальных релятивистских эффектов
  3. Оценка свойств отклика: оценка производительности метода для молекулярных свойств отклика

Глубокий анализ

Преимущества

  1. Высокая теоретическая новизна: впервые систематически расширены методы F12 на сбалансированное описание основного и возбуждённых состояний
  2. Ясная физическая интуиция: стратегия валентного ограничения, основанная на масштабах энергии орбиталей, имеет чёткую физическую мотивацию
  3. Достаточное экспериментальное обоснование: охватывает различные молекулярные системы и типы возбуждённых состояний, результаты убедительны
  4. Полная реализация: предоставляет полную автоматизированную схему реализации

Недостатки

  1. Ограниченность тестовых систем: проверено только на малых молекулах, производительность на больших молекулах неизвестна
  2. Недостаточный теоретический анализ: отсутствует глубокий анализ причин эффективности ограничения VVO
  3. Отсутствие анализа вычислительной стоимости: не обсуждается подробно вычислительная нагрузка относительно традиционных методов

Влияние

  1. Методологический вклад: предоставляет новую парадигму проектирования методов F12 для возбуждённых состояний
  2. Практическая ценность: может быть непосредственно применён для улучшения точности расчётов возбуждённых состояний
  3. Воспроизводимость: предоставляет детальные технические детали реализации

Области применения

  1. Фотохимические исследования: исследования, требующие высокоточных энергий и свойств возбуждённых состояний
  2. Расчёты спектроскопии: точное предсказание вертикальных энергий возбуждения и вибрационных частот
  3. Методологическое развитие: основа для дальнейшей разработки методов сбалансированного описания основного и возбуждённых состояний

Библиография

Данная работа цитирует 91 важную публикацию, охватывающую классические работы в области теории F12, трансворрелированных методов, расчётов возбуждённых состояний и других смежных областей, обеспечивая прочную теоретическую основу для исследования.

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