Графово‑теоретический анализ доли насыщающих отталкивающих примесей в твёрдых растворах

Почему важно, как плотно размещены атомы

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

Атомы на решётке

Авторы сосредотачиваются на замещающих твёрдых растворах, широкой категории сплавов, где каждая точка в регулярной атомной решётке (латисе) занята либо базовым атомом, либо примесным. Эксперименты показывают, что во многих системах — таких как железо‑хромовые стали, сложные сплавы с высокой энтропией и полупроводниковые сплавы IV группы, например германие‑олово — некоторые пары примесей практически никогда не оказываются соседями. Вместо этого они формируют паттерны, известные как короткодействующий порядок, где локальные конфигурации отклоняются от случайных. Эта скрытая упорядоченность сильно влияет на механические и электрические свойства, но её трудно увидеть прямо в экспериментах. Естественный, но ранее не ответленный вопрос: если примесные атомы должны избегать соседства, сколько их можно разместить в решётке, прежде чем это правило станет невыполнимым?

Простая игра на упаковку на решётке

Для ответа исследователи моделируют введение примесей как процесс случайной упаковки на решётке. Они представляют, что начинают с чистого базового материала и добавляют примесные атомы по одному. Каждый новый атом случайно ставится на узел, который ещё не занят примесью и не является соседом никакой примеси. Выбранный узел становится сайтом примеси; его соседние узлы блокируются для будущих примесей. Этот процесс продолжается, пока не исчерпаются подходящие узлы. Конечная доля узлов, занятых примесями, определяется как доля насыщения. С помощью компьютерных симуляций на 14 различных типах решёток — включая распространённые структуры, такие как объёмно-центрированная кубическая (встречается в сталях), гранецентрированная кубическая и более экзотические высокоразмерные решётки — авторы показывают, что каждая решётка имеет очень воспроизводимую долю насыщения, уникальный отпечаток того, как она вмещает отталкивающие примеси.

Графы, связи и универсальное правило

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

Когда локальные циклы меняют предел

Однако реальные кристаллы не являются полностью случайными сетями. В них присутствуют многие малые циклы связанных узлов — треугольники, квадраты, шестиугольники — которые тонко меняют емкость упаковки. Чтобы учесть это, авторы обращаются к ещё одному свойству графа, называемому girth (обхват): размер наименьшего цикла в сети. Сравнивая симуляции на реальных решётках с формулой для случайных графов, они обнаруживают систематический паттерн. Решётки, богатые трёхузловыми циклами (girth = 3), такие как гранецентрированная кубическая структура, как правило, имеют более низкие доли насыщения, чем предсказывает простая модель. Решётки, в которых преобладают четырёхузловые циклы (girth = 4), например простая кубическая и объёмно‑центрированная кубическая, могут упаковывать отталкивающие примеси плотнее, чем предполагает модель случайного графа. Структуры с большими циклами ближе к простой предсказанной зависимости. Даже одномерные цепочки и конечные кольца аккуратно вписываются в эту графово‑теоретическую картину.

От абстрактных графов к реальным материалам

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

Что это означает простыми словами

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

Цитирование: Kubo, A., Abe, Y. Graph-theoretic analyses of saturation fraction of repulsive dopants in solid solutions. Sci Rep 16, 7650 (2026). https://doi.org/10.1038/s41598-025-30829-1

Ключевые слова: отталкивающие примеси, короткодействующий порядок, случайные графы, проектация сплавов, доля насыщения