Зависимость в полулогарифмическом масштабе от обратной температуры приведена на рис. 8.
Увеличивая масштаб:
Т, К | 1123 | 1223 | 1323 | 1423 | 1473 | 1573 |
Nd, см-3 |
Т.о., в этих диапазонах можно решать уравнение Фика в приближении отсутствия зависимости коэффициента диффузии от температуры.
Видно, что с ростом температуры поверхностная концентрация, необходимая для начала нелинейной зависимости коэффициента диффузии, немного растет. Это означает, что при более высоких температурах, коэффициент диффузии становится менее чувствительным к примеси. Это можно объяснить тем, что при повышении температуры увеличивается количество и вакансий, и междоузлий, но количество междоузлий растет быстрее, а примесь по ним перемещается медленнее.
Возможно также, что отношение созданной концентрации к собственной при более низкой температуре больше, чем при более высокой, т.е.
, ввиду этого точка начала нелинейности смещается. Т.е. эта точка определяется собственной концентрацией при данной температуре. Это объяснение хорошо согласуется с приведенным выше, касающимся нелинейности зависимости коэффициента диффузии от концентрации доноров.Предельная растворимость сурьмы в кремнии составляет
см-3 при температуре 1573 К [7] и снижается до 2*1019 при 1123 К, таким образом можно утверждать, что коэффициент диффузии будет постоянен при любой концентрации доноров, то есть можно использовать решение уравнения Фика только для случая постоянного коэффициента диффузии.Максимальная погрешность такого решения по коэффициенту диффузии составит (для 1473 К):
,(82)то есть все найденные решения будут удовлетворять условию, чтобы погрешность не превышала 5%. Далее покажем погрешность глубины залегания примеси, при расчете методом постоянного коэффициента диффузии.
Отметим сразу, что собственная концентрация при выбранном диапазоне температур превышает
, в то время как максимальная предельная растворимость < см-3, то есть концентрация сурьмы в кремнии никогда не превысит собственную концентрацию при равновесном процессе.Пусть диффузия ведется на глубину 10 мкм. Зададим четыре температуры процесса (1173 К, 1273 К, 1373 К и 1473 К) и сравним полученные профили концентрации в этом случае на пятом временном слое. Выбранное t составляет полчаса.
Для этого решим уравнение диффузии:
(83)где
- безразмерные переменные.Полный текст программы приведен в приложении 1.
Для расчетов будем использовать программу, учитывающую непостоянство коэффициента диффузии.