Алгоритм от химиков МГУ вошел в известную программу квантовой химии

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

«В идеале необходимо использовать способ посередине между чистой теорией и чистой практикой, – добавляет Артём Митрофанов. – Это должен быть метод, использующий и известные данные о f-элементах, и подходы квантовой химии, и на балансе времени и точности надо построить рабочий инструмент для предсказания пути протекания реакций тяжёлых элементов».

Основой для такого инструмента может послужить DFT – теория функционала плотности, популярная среди химиков-вычислителей. Сейчас в мире насчитывается более 200 алгоритмов-функционалов, с помощью которой предсказываются различные свойства соединений. Их ограничения связаны с интерполяционным характером: полноценное и правильное описать соединения актиноидов с большим количеством электронов невозможно с достаточной точностью, поскольку разработчики существующих подходов не рассчитывали на особенности соединений f-элементов. Поэтому необходима разработка нового функционала с параметрами, корректно описывающими поведение этой группы элементов. Именно такую задачу поставили себе сотрудники кафедры радиохимии химического факультета МГУ, которые соединили знания из квантовой химии и машинное обучение для создания метода, обладающего как описательной, так и предсказательной силой. Исследование поддержал Российский научный фонд (проект № 19-73-20115), и работа началась.

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

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

«Новый функционал может создаваться из комбинации уже существующих, их параметров и весовых коэффициентов, приписываемых им программой. Это многопараметрическая задача, которая не решается в лоб, перебором, поэтому мы и использовали алгоритм глобальной оптимизации», – рассказывает Артём Митрофанов. Учёным удалось показать, что точность их функционалов, которые получаются в результате таких операций, не уступает существующим. Ещё одно немаловажное преимущество такого подхода к описанию химии актинидов заключается в минимальной необходимости человеческого участия в работе алгоритма. Химикам МГУ удалось создать функционал, который, не требуя математической мощности суперкомпьютера, с экспериментальной точностью может описывать химические взаимодействия соединений актинидов в газовой фазе.

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