Точный подход DFT-1/2 для слабых дефектных состояний: эффективный расчёт энергий связи доноров в кремнии

Почему важны крошечные изменения в кремнии

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

Атомы, дарящие дополнительные электроны

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

Почему стандартные расчёты бывают недостаточны

Компьютерные модели на основе теории функционала плотности (DFT) — основной инструмент проектирования материалов, но они склонны занижать степень локализации электронов и неверно позиционировать края зон в полупроводниках. Для слабых доноров, чьи электронные облака простираются на десятки атомов, это означает, что DFT обычно предсказывает слишком маленькие энергии связи. Более продвинутые методы, такие как гибридные функционалы и расчёты GW, могут исправить эти недостатки, но стоят чрезвычайно дорого по вычислительным ресурсам, особенно когда требуются большие вычислительные ящики с тысячами атомов, чтобы корректно отразить протяжённое донорное состояние. Ранее применявшиеся «тандемные» подходы приходилось комбинировать на разных уровнях теории и в разных размерах ячеек, а затем стыковать результаты, что делало рабочий процесс сложным и зависимым от системы.

Простая корректировка с большим эффектом

Авторы опираются на метод DFT-1/2, который вносит приближённую коррекцию собственной энергии непосредственно в стандартный DFT. На практике это означает небольшую модификацию эффективного потенциала выбранных атомов путём мысленного удаления половины электрона из определённых атомных орбиталей. Сначала они применяют эту корректировку к объёму кремния, чтобы рассчитанная ширина запрещённой зоны лучше соответствовала эксперименту и обеспечивала надёжную привязку для зоны проводимости. Затем анализируют электронную природу донорного состояния и обнаруживают, что для всех доноров пятой группы оно доминирует s-орбиталью атома донора. После этого они накладывают целевую полууэлектронную коррекцию на эту орбиталь и подбирают единственный радиус отсечки, чтобы максимизировать разделение между донорным уровнем и ближайшим пустым состоянием зоны проводимости. Важно, что эта оптимизированная корректировка остаётся применимой при увеличении размера расчётной ячейки, поэтому её можно повторно использовать в суперячейках размером до тысяч атомов.

Насколько хорошо работает метод

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

Практичный инструмент для проектирования будущих чипов

Для неспециалистов ключевой вывод в том, что авторы разработали рецепт, который сохраняет низкую стоимость и простоту стандартных расчётов DFT, одновременно достигая точности гораздо более тяжёлых методов. Систематически корректируя и общую зонную структуру материала-хоста, и локальную среду вокруг донора, их протокол DFT-1/2 обеспечивает надёжные энергии связи доноров в очень больших расчётных ячейках. Это делает его мощным и универсальным инструментом для изучения примесей, управляющих рабочими характеристиками повседневной электроники и новых квантовых устройств, помогая инженерам проектировать материалы, где отдельные примеси ведут себя именно так, как задумано.

Цитирование: Claes, J., Partoens, B., Lamoen, D. et al. An accurate DFT-1/2 approach for shallow defect states: efficient calculation of donor binding energies in silicon. npj Comput Mater 12, 153 (2026). https://doi.org/10.1038/s41524-026-02003-2

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