- •Ю.П. Головатый, в.Г. Косушкин
- •Глава 6. Численное решение краевых задач для обыкновенных дифференциальных уравнений
- •Глава 7. Численное решение задач на собственные значения
- •Глава 8. Численное решение дифференциальных уравнений в частных производных
- •Глава 1. Формулировка математической модели
- •1.1. Тепловой баланс резистивного элемента
- •1.2. Радиационно-стимулированная диффузия
- •Глава 2.Решение модельной задачи
- •Глава 3. Метод конечных разностей
- •3.1. Сетки и сеточные функции
- •3.2. Аппроксимация производных
- •Конечноразностная аппроксимация обыкновенных дифференциальных уравнений
- •Глава 4. Решение задачи коши методом конечных разностей
- •Методы эйлера и тейлора
- •4.2 Явные методы рунге - кутты
- •4.3 Явные многошаговые методы решения задачи коши.
- •Лабораторная работа №3
- •4.4 Подвохи при применении методов рунге-кутты
- •«Жёстские» обыкновенные дифференциальные уравнения
1.2. Радиационно-стимулированная диффузия
Одной из важнейших электронных технологий является ионная имплантация – внедрение в полупроводник–мишень ионов легирующей примеси, ускоренных до энергий в десятки и сотни килоэлектронвольт. В результате столкновений с атомными ядрами и электронами мишени внедряемый ион теряет энергию и в конечном итоге останавливается на некотором расстоянии от точки падения на поверхность мишени. Процесс торможения является статистическим, поэтому предсказать точку остановки конкретного иона нельзя. Но можно вывести теоретически функцию пространственного распределения ионов в полубесконечном полупроводнике. Впервые это сделали в 1963 году датские физики Линдхард, Шарф и Шиотт. В одномерном случае, что физически соответствует облучению широким ионным пучком, функция распределения ЛШШ имеет вид
(1.12)
Если доза облучения составляет , то концентрация примеси на глубине от поверхности равна
(1.13)
Параметры и называются средним проективным пробегом и разбросом пробегов соответственно. Они монотонно возрастают с ростом энергии ионов. Метод ионной имплантации позволяет гибко и точно управлять глубиной залегания p-n –перехода. Но неизбежным сопутствующим эффектом является образование огромного числа неравновесных радиационных дефектов – вакансий, междоузельных атомов и более сложных комплексов. Поэтому после имплантации приходится проводить термический отжиг полупроводника при температуре порядка для восстановления его кристаллической структуры и электрического активирования примесных атомов. При такой температуре диффузионная подвижность примесных атомов становится столь высокой, что происходит их пространственное перераспределение и деформация исходного профиля. На перераспределение примесных атомов оказывают самое существенное влияние радиационные дефекты, так как внедрённые атомы диффундируют в составе примесно-дефектных комплексов, образованных ими с вакансиями и междоузельными атомами. При диффузии по вакансионному механизму коэффициент диффузии атома можно приближённо представить в виде
(1.14)
где - равновесный коэффициент диффузии, - равновесная концентрация вакансий, - их реальная концентрация. Избыточные нравновесные вакансии сами распределены в полупроводнике неоднородно. После имплантации их распределение может быть в первом приближении описано функцией распределения , аналогичной (1.12),
(1.15)
Параметр задаёт положение максимума распределения вакансий, он всегда меньше . есть ширина распределения, она, как правило, немного больше . При температурах отжига подвижность вакансий очень высока, они быстро диффундируют вглубь полупроводника и рекомбинируют на стоках. Скорость рекомбинации характеризуется временем жизни вакансий . В результате профиль распределения вакансий является функцией координат и времени. Следовательно, зависит от координат и времени также и коэффициент диффузии примеси . Поэтому перераспределение примеси при отжиге описывается в первом приближении двумя диффузионными уравнениями:
для вакансий (1.16)
для примеси (1.17)
Начальные условия для них задаются формулами (1.13) и (1.15),
(1.18)
(1.19)
В глубине полупроводника концентрация вакансий равна равновесной, а концентрация примеси равна нулю,
(1.20)
(1.21)
Граница полупроводника является непроницаемой для примеси и стоком для вакансий, поэтому граничные условия имеют вид
(1.22)
(1.23)
Уравнения (1.16) – (1.23) составляют математическую модель перераспределения примеси с участием радиационных вакансий.