Мир современной теоретической химии и физики долгое время напоминал попытку усидеть на двух стульях. С одной стороны – строгое квантово-механическое описание. Оно безупречно точно, но практически неприменимо для реальных систем: как только молекула становится чуть сложнее, вычислительные затраты растут экспоненциально.
С другой стороны – классическая механика. Она «дешевле» для суперкомпьютеров, но совершенно бессильна перед квантовыми эффектами: интерференцией, туннелированием и нулевой колебательной энергией.
Ученые десятилетиями строили мост над этой пропастью. И вот, сотрудники Кольского научного центра РАН совместно со швейцарскими коллегами из Лозанны предложили решение, которое может изменить правила игры. Статья об этом опубликована в The Journal of Chemical Physics.
В чем была «загвоздка» десятилетий?
Одним из главных инструментов в этой области считается пропагатор Германа-Клука. Его идея безусловно элегантна: мы берем ансамбль классических траекторий и снабжаем каждую специальной «фазой» и амплитудой, чтобы учесть квантовые свойства.
Но здесь крылась некоторая ловушка: вклад каждой траектории постоянно осциллирует (колеблется). При росте размерности системы эти колебания становятся настолько хаотичными, что статистическая сходимость метода резко ухудшается. Проще говоря, компьютер начинает выдавать «цифровой шум» вместо точного графика, и расчет становится практически невозможным.
Клеточная фильтрация: от «ручного управления» к автоматике
Чтобы побороть этот шум, раньше использовали «фильтрацию Филинова» – фазовое пространство разбивали на ячейки. Но у нее был принципиальный изъян: размер ячеек вводился как внешний параметр, не связанный с числом траекторий. Даже если вы запускали бесконечное количество расчетов, метод не гарантировал идеальной точности.
Авторы новой работы – Фабиан Кренингер, Иржи Ваничек и Сергей Антипов – предложили элегантный выход: связать размер ячеек в фазовом пространстве с числом используемых траекторий.
Суть инновации заключается в том, чтобы вместо фиксирования ширины ячеек заранее, заставить их «масштабироваться». Когда вы увеличиваете число траекторий, ячейки автоматически уменьшаются, а плотность выборки расширяется.
Почему это открытие является «фундаментальным мостом»?
Новая схема впервые в истории науки непрерывно связала два подхода, которые раньше считались концептуально разными:
Первый – Предел одной траектории. Метод автоматически превращается в знаменитое «приближение оттаявшего гауссова волнового пакета» Эдварда Хеллера (1975 г.).
Второй – Предел бесконечных траекторий. Метод строго переходит в исходный пропагатор Германа-Клука, устраняя его главный концептуальный недостаток.
Проверка открытия на прочность
Теорию проверили на практике, и результаты просто впечатляют!
● Молекула NCO: В расчетах колебательных спектров этой молекулы новый метод «практически идеально воспроизвел результаты», но при существенно меньших вычислительных затратах.
● Победа над хаосом: В системах с сильным хаосом старые методы «взрывались»: выдавали нефизические, безумные результаты из-за неустойчивости траекторий. Новая «уточненная клеточная фильтрация» эффективно подавляет эти аномалии, позволяя получить осмысленное описание динамики там, где раньше был тупик.
Где же пригодятся данные исследования? (спойлер: везде)
Улучшение стабильности расчетов без усложнения стоимости одной траектории открывает двери в мир, который раньше был заперт.
Например:
– В фундаментальной химии: для понимания того, как на самом деле эволюционируют молекулярные системы во времени.
– В спектроскопии: для расшифровки сложнейших экспериментальных спектров.
– В дизайне материалов: когда нужно предсказать поведение сложных молекул еще до того, как они будут созданы в пробирке.
Ученые фактически создали универсальный инструмент, который сам настраивает баланс между точностью и скоростью. Это и есть тот самый «золотой стандарт», к которому стремится современная вычислительная химия.