29 «магических» кластеров и рекордно малое кольцо: как рассчитали устойчивость наноструктур сульфидов молибдена
Исследователи лаборатории компьютерного дизайна материалов Центра фотоники и двумерных материалов МФТИ провели масштабное моделирование соединений молибдена и серы и составили практически полную «карту» устойчивых нанокластеров Mo–S. В ходе вычислительного перебора более четырех тысяч возможных структур они выявили 29 особо стабильных конфигураций — так называемых магических кластеров. Лишь семь из них ранее удавалось получить и зафиксировать экспериментально.
Помимо этого, ученые описали, вероятно, одну из предельно малых плоских наноструктур для системы Mo–S: метастабильное кольцо, состоящее из 10 атомов молибдена и 20 атомов серы (Mo10S20). Такая архитектура представляет собой двумерное кольцевое образование, которое может рассматриваться как зачаточный фрагмент ультратонких пленок или нанопроволок. Работа открывает новые перспективы для проектирования катализаторов и материалов для литий‑ионных аккумуляторов на основе точно просчитанных нанокластеров.
Почему молибденовые сульфиды так важны
Интерес к системе молибден–сера связан прежде всего с дисульфидом молибдена (MoS₂) — слоистым материалом, который в виде одного или нескольких атомных слоев проявляет свойства, позволяющие рассматривать его как альтернативу кремнию в электронике. В двумерной форме MoS₂ обладает подходящей шириной запрещенной зоны, высокой механической прочностью и хорошими электрооптическими характеристиками.
Но путь к реальным устройствам лежит не только через сплошные пленки. Переход к наномасштабу выводит на первый план изучение отдельных кластеров Mo–S — небольших совокупностей атомов, которые могут служить «кирпичиками» для создания более сложных наноструктур. Такие кластеры уже сегодня рассматриваются как активные центры в катализаторах, чувствительные элементы сенсоров и компоненты электродных материалов.
До сих пор экспериментально были достоверно идентифицированы всего семь устойчивых нанокластеров молибденовых сульфидов:
- Mo₃S₁₃
- Mo₄S₆
- Mo₆S₈
- Mo₉S₁₁
- Mo₁₂S₁₄
- Mo₃S₄
- Mo₆S₄
Все они уже нашли применение в разработке новых катализаторов и перспективных литий‑ионных батарей. Однако их кристаллическая структура и возможные аналоги с иным числом атомов остаются малоизученными. Нанокластеры заметно отличаются по строению от объемных кристаллов: меняется координация атомов, тип связей, распределение зарядов, а вместе с этим — и физико‑химические свойства. Экспериментальным методам часто попросту не хватает разрешения или чувствительности, чтобы однозначно определить структуру столь малых систем, особенно если они существуют в виде смеси нескольких разновидностей.
Как вычисляют «невидимую» структуру
Чтобы обойти ограничения эксперимента, команда МФТИ опиралась на мощный арсенал вычислительных инструментов. В работе был использован связанный подход, включающий:
1. Эволюционный алгоритм USPEX — программный комплекс, который перебирает огромное количество возможных конфигураций атомов и отбирает более стабильные варианты по принципу «естественного отбора».
2. Метод теории функционала плотности (DFT) — квантово‑механический подход, позволяющий с высокой точностью рассчитывать энергию и электронную структуру систем на атомном уровне.
3. Машинно‑обучаемые межатомные потенциалы MTP (Moment Tensor Potentials) — модель, которая после обучения на результатах DFT‑расчетов способна быстро и достаточно точно оценивать энергию и динамику больших систем, где прямые квантовые расчеты были бы слишком затратны.
MTP‑потенциалы использовались в частности для анализа динамической стабильности нанокластеров и монослоев, содержащих более 30 атомов. Это позволило существенно расширить диапазон исследуемых структур без необходимости проводить для каждой из них дорогие по ресурсам DFT‑вычисления.
По словам заведующего лабораторией компьютерного дизайна материалов Ивана Круглова, для нанокластеров большого размера квантово‑механические методы требуют колоссального объема вычислений, поэтому ключевым этапом стала подготовка и обучение межатомных потенциалов MTP. С их помощью в том числе оценивалась стабильность необычных плоских колец Mo10S20, образующих монослойные структуры.
Два критерия стабильности и понятие «магических» кластеров
Физики предложили строгую систему критериев, чтобы отобрать действительно устойчивые образования. Для каждого кластера рассчитывали:
- Минимальную энергию перестройки (Δmin) — ту энергию, которую система должна «заплатить», чтобы изменить свою конфигурацию при добавлении или удалении одного атома молибдена или серы. Чем выше эта величина, тем менее выгодно кластеру трансформироваться в другой.
- Энергию диссоциации (Ediss) — количество энергии, необходимое для полного распада кластера на простые фрагменты или отдельные атомы. Положительная энергия диссоциации означает, что самопроизвольный распад невозможен: для разрушения кластера нужно внешнее воздействие.
Кластеры, которые одновременно обладают высокой минимальной энергией перестройки и положительной энергией диссоциации, исследователи отнесли к разряду магических. Такое название подчеркивает, что у них особенно выгодная геометрия и электронное строение: они заметно стабильнее соседей с близким числом атомов.
На первом этапе этих условий удовлетворял 31 кандидат. Затем ученые перешли к более тонкому анализу — изучили колебательные спектры структур, чтобы проверить их динамическую устойчивость. В рамках такой проверки ищут «мнимые» (воображаемые) частоты колебаний: их наличие свидетельствует о том, что структура легко «сползает» в соседнее состояние и, по сути, не является устойчивым минимумом энергии. После этого шага два кластера исключили, и в итоге список магических сузился до 29 устойчивых нанокластеров.
Пять семейств устойчивых наноструктур
Все найденные структуры исследователи разбили на пять групп в зависимости от соотношения атомов молибдена и серы:
1. Кластеры с явным преобладанием молибдена (обедненные серой).
2. Системы с примерно равным числом атомов Mo и S.
3. Соединения с отношением Mo:S ≈ 1:1,5.
4. Кластеры серия 1:2 (типичные дисульфидные составы).
5. Структуры с значительным избытком серы, где она формирует насыщенные оболочки вокруг ядер из атомов молибдена.
Такое разбиение помогает понять, как меняются свойства кластеров при изменении химического состава. Например, избыток серы часто стабилизирует внешнюю оболочку и может усиливать каталитическую активность за счет большего числа реакционно‑способных центров на поверхности. В то же время молибдено‑богатые кластеры ближе по поведению к металлическим наночастицам и интересны с точки зрения проводимости и структурной жесткости.
Рекордное плоское кольцо Mo10S20
Особое внимание в работе уделено необычной наноструктуре — плоскому кольцу Mo10S20, которое исследователи относят к метастабильным состояниям. Такое кольцо можно представить как фрагмент двумерного слоя, замкнутый в округлую форму.
Моделирование показало, что при определенных условиях это кольцо сохраняет свою форму и не распадается, а его динамическая устойчивость подтверждается расчетами с использованием MTP‑потенциалов. Авторы считают, что это одна из минимально возможных плоских структур в системе Mo–S, которая все еще удовлетворяет критериям устойчивости.
Интерес к таким кольцевым образованиям не только академический. Кольца и цепочки могут служить строительными блоками для создания нанопроволок, нанокольцевых сеток и других одномерных и квазидвумерных систем. Они потенциально пригодны в наноэлектронике, фотонике или в качестве специфических активных центров в катализе, поскольку локальная геометрия и краевые атомы в таких структурах заметно отличаются от классических плоских слоев.
Значение работы для катализаторов
Многие промышленные катализаторы уже основаны на сульфидах молибдена, однако до сих пор точное понимание, какие именно нанофрагменты отвечают за максимальную активность, оставалось неполным. Наличие карты из 29 магических кластеров позволяет целенаправленно проектировать катализаторы, где преобладают именно наиболее устойчивые и функционально активные конфигурации.
Например, зная, что определенный кластер обладает высокой энергией диссоциации, можно ожидать, что он сохранит свою структуру в жестких условиях реакции — при высоких температурах и давлениях или в агрессивных средах. В то же время свойства активных центров (краевых атомов, дефектов, открытых площадок) можно настраивать за счет выбора состава (богатый серой или молибденом кластер) и его окружения.
Также перспективным направлением выглядит модульный подход: катализатор может быть спроектирован как ансамбль разных магических кластеров, где каждая разновидность отвечает за свой этап сложной реакции — адсорбцию, активацию молекулы, перенос заряда и последующую десорбцию продуктов.
Нанокластеры и литий‑ионные аккумуляторы
Сульфиды молибдена активно рассматриваются как материалы для электродов литий‑ионных и натрий‑ионных батарей. При заряде и разряде ионы лития внедряются в структуру материала, изменяя ее объем, локальную химию и электронные свойства. От того, насколько устойчивы нанофрагменты внутри электрода, зависит долговечность и емкость аккумулятора.
Картина из 29 устойчивых кластеров может служить своеобразным «каталогом» возможных структурных мотивов в катодах и анодах на основе Mo–S. Зная, какие кластеры сохраняются при многократном циклировании, а какие распадаются или переходят в другие формы, исследователи могут прогнозировать деградацию материала и подбирать оптимальные условия синтеза, чтобы максимизировать долю устойчивых нанофаз.
Кроме того, модельные расчеты взаимодействия лития с конкретными магическими кластерами позволят заранее оценить, какие структуры обеспечат наилучший баланс между емкостью, скоростью диффузии и стабильностью при циклировании. Это особенно важно для создания высокоэнергетических и одновременно долговечных аккумуляторов нового поколения.
Роль машинного обучения в материаловедении
Использование MTP‑потенциалов демонстрирует, как машинное обучение меняет подход к исследованию материалов. Раньше каждая новая структура требовала отдельного тяжелого квантового расчета, и перебор тысяч вариантов был почти нереализуем. Теперь обученная модель способна за доли секунды оценивать энергию и динамику новой конфигурации, оставаясь при этом близкой по точности к DFT.
Для области нанокластеров это критично: число возможных комбинаций и геометрий растет лавинообразно с увеличением количества атомов. Эволюционные алгоритмы, работающие в связке с машинно‑обучаемыми потенциалами, превращают эту комбинаторную проблему в решаемую задачу и открывают путь к полноценному компьютерному дизайну материалов, когда нужная структура подбирается или даже «выводится» математически до эксперимента.
Что дальше: от теории к эксперименту
Следующий очевидный шаг — экспериментальное подтверждение предсказанных кластеров. Это задача не из простых: нужно научиться стабилизировать конкретные структуры, отделять их от соседних состояний и надежно регистрировать с помощью спектроскопии, дифракционных методов или микроскопии с атомным разрешением.
Однако наличие четко рассчитанных моделей заметно облегчает работу экспериментаторов. Они получают конкретные «мишени» — набор нанокластеров с описанной геометрией, энергией и предполагаемыми спектральными характеристиками. Это значит, что можно целенаправленно подстраивать условия синтеза: температуру, давление, состав газовой фазы, тип подложки, чтобы увеличить выход тех или иных магических кластеров.
В долгосрочной перспективе подобные исследования создают базу для рационального проектирования наноструктур: вместо случайного подбора параметров синтеза появится возможность осознанно собирать материалы из заранее просчитанных устойчивых блоков.
Перспективы для наноэлектроники и фотоники
Кластеры Mo–S интересны не только как катализаторы и электродные материалы. Их электронные и оптические свойства могут оказаться полезными для наноэлектроники, фотоники и сенсорики. Например, изменяя размер и состав кластера, можно тонко настраивать ширину энергетической щели и характер переходов, а значит — управлять поглощением и испусканием света в нужном диапазоне.
Возможность формировать из магических кластеров упорядоченные ансамбли — массивы, цепочки, кольцевые структуры — открывает перспективы создания миниатюрных квантовых точек, светочувствительных элементов, нелинейных оптических компонентов. В сочетании с уже известными свойствами двумерного MoS₂ это может привести к появлению гибридных устройств, где сплошные слои и отдельные кластеры работают совместно.
***
Работа коллектива МФТИ фактически создает фундаментальный справочник по устойчивым наноструктурам в системе молибден–сера: от привычных объемных аналогов исследователи переходят к уровню отдельных нанокластеров и минимальных плоских колец. Выявленные 29 магических кластеров и описанное кольцо Mo10S20 задают ориентиры для дальнейшего развития катализаторов, аккумуляторов и наноэлектронных устройств, в основе которых лежит тонкая инженерия структуры на уровне нескольких десятков атомов.
