Предложен новый метод описания сложных молекулярных систем

Ученые из Института химии и технологии редких элементов и минерального сырья им. И.В. Тананаева Кольского научного центра РАН совместно с зарубежными коллегами предложили новый полуклассический метод описания сложных молекулярных систем. Переосмысление математической структуры полуклассических приближений позволило авторам не только повысить вычислительную эффективность, но и прояснить связи между различными исторически сложившимися методами.

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

Одним из наиболее широко применяемых полуклассических методов является пропагатор Германа-Клука, предложенный Майклом Германом и Эриком Клуком в 1984 году. Он используется для описания временной эволюции квантовых волновых пакетов, особенно в молекулярной динамике и спектроскопии, и позволяет учитывать важные квантовые эффекты, используя при этом ансамбль классических траекторий.

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

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

Однако у стандартной клеточной фильтрации имеется принципиальный недостаток: размер ячеек вводится как внешний параметр, не связанный напрямую с числом используемых траекторий. В результате даже при бесконечном числе траекторий метод не гарантирует сходимости к исходному пропагатору Германа-Клука, а исследователю приходится дополнительно анализировать зависимость результата от ширины ячеек, что увеличивает вычислительные затраты и снижает концептуальную строгость метода.

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

Авторы иллюстрируют преимущества предложенной схемы на ряде модельных задач, охватывающих как интегрируемые, так и хаотические динамические режимы. В качестве примера относительно регулярной динамики рассматривается модель коллинеарной молекулы NCO, в которой волновой пакет эволюционирует на возбужденной электронной поверхности. Расчеты временных автокорреляционных функций и соответствующих колебательных спектров показывают, что уточненная клеточная фильтрация практически идеально воспроизводит результаты полностью сходившегося пропагатора Германа-Клука при существенно меньших вычислительных затратах. В то же время стандартная фильтрация Филинова демонстрирует заметные искажения спектральных линий при недостаточно малых ячейках, что подчеркивает ее чувствительность к выбору параметров.

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