Как предсказать новый материал: союз машинного обучения и квантовой механики ✛ ✦

Оригинал: "Data-Driven Prediction of NaCl-Type Entropy-Stabilized Oxide Compositions from First-Principles and Supervised Learning"
arXiv:2607.04502 · 2026-07-05 · CC BY 4.0 · 4 мин
Высокопроизводительный скрининг энтропийно-стабилизированных оксидов с помощью DFT и нейросети открывает путь к ускоренному дизайну функциональных материалов.

Контекст

Поиск новых неорганических материалов с уникальными свойствами — одна из центральных задач современной химии твердого тела. Традиционный подход «синтез–тестирование» крайне трудоемок, особенно когда речь идет о многокомпонентных системах: пространство возможных комбинаций элементов огромно и напоминает поиск иголки в стоге сена. Прорыв произошел с открытием высокоэнтропийных оксидов (ВЭО) — твердых растворов из пяти и более металлов в эквимолярном соотношении, стабилизированных именно высокой конфигурационной энтропией смешения. Как предсказывал еще Людвиг Больцман, беспорядок может быть движущей силой стабильности — вступает в силу второй закон термодинамики, согласно которому высокоэнтропийное состояние выгодно при нагреве. Однако до сих пор не существовало надежных методов быстрого скрининга тысяч возможных составов на предмет того, образуют ли они однофазный ЭСО или распадаются на смесь более простых оксидов. Это ограничивало практическое применение концепции энтропийной стабилизации.

Методы

Авторы построили многоуровневый виртуальный конвейер. На первом этапе с помощью численного моделирования в рамках теории функционала плотности (DFT) они рассчитали энтальпии образования всех возможных упорядоченных оксидов — бинарных и тройных — для 16 отобранных катионов: от щелочноземельных до переходных металлов (включая железо, кобальт, никель, медь, цинк) и постпереходных элементов. Особое внимание уделили элементам с высокой металличностью — то есть металлам, способным образовывать оксиды с разной валентностью. Беспорядок в пятивкомпонентных твердых растворах моделировали с помощью специальных квазислучайных суперячеек (SQS) размером до 60 атомов, что позволяет воспроизвести статистически случайное распределение катионов по подрешетке на фоне жесткой анионной решетки из кислорода. Для обучения моделей были проведены полноценные DFT-расчеты 416 пятерных конфигураций (около миллиона ядро-часов). Затем к делу подключилось машинное обучение: три модели (линейная регрессия, случайный лес и нейросеть-перцептрон) обучались предсказывать ключевой параметр — расстояние до выпуклой оболочки стабильности в многомерном пространстве составов и вычисляемую из него температуру стабилизации Tstab. В роли дескрипторов выступали атомный номер, радиус, электроотрицательность и число валентных электронов. Термодинамический баланс описывался через квантовую термодинамику: свободная энергия смеси включала DFT-энтальпию и идеальную конфигурационную энтропию, подчиняющуюся распределению Больцмана.

Результаты

Наилучшие результаты показал многослойный перцептрон (MLP): на тестовой выборке из 104 систем ошибка составила всего 4.2 кДж/моль (против 6.1 у линейной модели). Это означает, что нейросеть улавливает тонкие нелинейные взаимосвязи между химической природой катионов и стабильностью оксида. Применив MLP ко всем 4368 возможным эквимолярным комбинациям, исследователи получили ранжированный список кандидатов. Уже известные ЭСО, такие как (Co,Cu,Mg,Ni,Zn)O, оказались среди лидеров: после коррекции более точным функционалом r2SCAN его Tstab снизилась с 2745 К до 1249 К, приблизившись к экспериментальному значению ~1100 К. Среди новых перспективных составов выделились (Ca,Cu,Fe,Mg,Zn)O с Tstab около 1676 К и (Ca,Cu,In,Sr,Zn)O с 1750 К. Экспериментальная проверка трёх отобранных композиций подтвердила основные тренды: в одних случаях возникала фаза каменной соли (хотя с частичной сегрегацией кальция), в других плавление наступало раньше расчётной стабилизации, что указало на необходимость учёта жидкой фазы. Доминирующие конкурирующие фазы модель предсказала верно.

Значение

Разработанный подход меняет парадигму поиска новых энтропийно-стабилизированных оксидов. Теперь вместо слепого перебора экспериментатор получает карту перспективных составов, основанную на строгих физических принципах и подкреплённую квантовыми расчётами. Это не только экономит годы лабораторной работы, но и позволяет рационально интерпретировать неудачи синтеза — например, выявить, какие именно конкурирующие фазы мешают образованию целевого материала. С фундаментальной точки зрения, работа подтверждает, что стабильность высокоэнтропийных оксидов не сводится к одной лишь конфигурационной энтропии: энтальпийные факторы (перенос заряда, совместимость степеней окисления, локальные искажения) играют не менее важную роль, и именно их нелинейные комбинации успешно выучивает нейросеть.

Развитие

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

Влияние

Результаты окажут непосредственное влияние на дизайн материалов в нескольких передовых областях: в энергетике — поиск новых катодов для литий-ионных аккумуляторов с высокой ионной проводимостью; в катализе — создание стабильных при высоких температурах катализаторов; в электронике — многофункциональные оксиды с управляемыми диэлектрическими и магнитными свойствами.

Следующие шаги

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

Ключевые проблемы

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

Взгляд машины знаний

Этот раздел адресован автору работы и его коллегам. Соседние работы нашёл смысловой поиск по нашему архиву — не по совпадению слов, а по смыслу. Мы не рецензируем: мы показываем поле вокруг.

Что сделано в работе

Авторы построили виртуальный конвейер для предсказания стабильности высокоэнтропийных оксидов: DFT-расчёты энтальпий для 16 катионов, моделирование беспорядка через SQS, обучение трёх ML-моделей на 416 пятерных конфигурациях. Лучшая модель (MLP) достигла ошибки 4.2 кДж/моль при предсказании расстояния до выпуклой оболочки и температуры стабилизации.

В чём сила

Сочетание DFT-расчётов с машинным обучением на дескрипторах из фундаментальных свойств элементов позволило быстро оценивать стабильность огромного пространства составов. Это даёт практичный инструмент для скрининга новых керамик без полного перебора.

Куда можно двигаться

Что можно попробовать

  • Попробовать использовать активное обучение: итеративно выбирать составы с наибольшей неопределённостью модели для дополнительных DFT-расчётов.
  • Расширить набор дескрипторов, добавив, например, энергии ионизации или параметры кристаллического поля, как в соседних работах.
  • Проверить переносимость обученной модели на другие классы высокоэнтропийных материалов, например, нитриды или карбиды.

Почему это может быть важно

Быстрый и точный скрининг высокоэнтропийных оксидов может ускорить открытие материалов для энергетики и высокотемпературных приложений, сократив цикл от идеи до синтеза.

Оригинальная аннотация

как её написали авторы, на английском arXiv:2607.04502

Entropy-stabilized oxides (ESOs) open access to vast multicomponent compositional spaces, but identifying promising candidates remains challenging because of the large number of possible mixtures and the need to assess their stability against competing phases. In this work, we develop a high-throughput computational framework to screen equimolar quinary ESOs in the NaCl structure type by combining density functional theory (DFT), special quasirandom structures (SQS), convex-hull thermodynamics, and supervised machine learning. A consistent reference database of binary and ternary ordered oxides, including disordered phases such as all binary cation combinations in the NaCl-type oxide, is first constructed using GGA and meta-GGA calculations. Quinary disordered phases are then described by SQS supercells and used to train machine-learning models that predict the distance to the convex hull and the corresponding stabilization temperature over the full set of 4368 possible equimolar quinary compositions generated from 16 cation species. Among the tested models, an optimized multilayer perceptron provides the best predictive performance, with a test error of about 4 kJ/mol, while requiring explicit DFT calculations for only about 10% of the quinary systems. Comparison with experimental synthesis tests and computed decomposition paths further shows that the approach captures the main stability trends and the dominant competing phases, although absolute stabilization temperatures remain affected by systematic thermodynamic approximations. These results establish an efficient route for the data-driven exploration of multicomponent oxides and provide practical guidance for the experimental search for new ESOs.

Если бы вы решили проверить все 4368 пятивкомпонентных комбинаций в лаборатории, тратя на синтез и анализ каждой всего неделю, вам потребовалось бы около 84 лет непрерывной работы. Нейросеть же справляется с задачей за доли секунды после обучения.

\Delta S_{\text{conf}} = -R \sum_{i=1}^{n} y_i \ln y_i
Конфигурационная энтропия смешения n компонентов с мольными долями y_i и универсальной газовой постоянной R.
\Delta_{\text{hull}} H(\mathbf{x}_\varphi) = \Delta_f H(\mathbf{x}_\varphi) - E_{\text{GS}}(\mathbf{x}_\varphi)
Разность между энтальпией образования соединения и энергией основного состояния (выпуклой оболочки) при том же составе; если >0 — соединение нестабильно при 0 K.

Ключевые числа

  • количество изученных пятерных систем: 4368
  • среднеквадратичная ошибка MLP-модели: 4.2 кДж/моль
  • число катионов в библиотеке: 16
  • размер суперячейки для моделирования беспорядка: 60 атомов
  • стабилизационная температура (Co,Cu,Mg,Ni,Zn)O (после коррекции): 1249 K
Учёные
Ludwig BoltzmannRudolf ClausiusJacob BekensteinStephen HawkingJames Clerk MaxwellEdward Witten
Понятия
энсо многомерный индекс ЭНСО