Курсовая: Расчет распределения примесей в кремнии при кристаллизационной очистке и диффузионном легировании - текст курсовой. Скачать бесплатно.
Банк рефератов, курсовых и дипломных работ. Много и бесплатно. # | Правила оформления работ | Добавить в избранное
 
 
   
Меню Меню Меню Меню Меню
   
Napishem.com Napishem.com Napishem.com

Курсовая

Расчет распределения примесей в кремнии при кристаллизационной очистке и диффузионном легировании

Банк рефератов / Радиоэлектроника

Рубрики  Рубрики реферат банка

закрыть
Категория: Курсовая работа
Язык курсовой: Русский
Дата добавления:   
 
Скачать
Архив Zip, 146 kb, скачать бесплатно
Заказать
Узнать стоимость написания уникальной курсовой работы

Узнайте стоимость написания уникальной работы

27 Государственный комитет РФ по высшему образованию Новгородск ий Государственный университет им . Ярослава Мудрого Кафедра физики твердого тела и микроэлектроники ”РАСЧЕТ РАСПРЕДЕЛЕНИЯ ПРИМЕСЕЙ В КРЕМНИИ ПРИ КРИСТАЛЛИЗАЦИОНОЙ ОЧИСТКЕ И ДИФФУЗИОННОМ ЛЕГИРОВАННИИ” Пояснительная записка к курсовой работе по дисциплине ”Физико-химические основы технологии микроэлектроники” Выполнил Студент группы 7033у Н.Е.Родин Проверил Преподаватель к афедры ФТТМ Б.М.Шишлянников 1998 Техническое задание на курсовую работу по дисциплине “Физико-химические ос новы технологии микроэлектроники” Студенту гр . 7033 Родину Н.Е. 1. Рассчитать распределение примесей вдоль слитка полупроводникового материала при очистке зонной плавкой (один проход расплавленной зоной ). материал кремн ий примеси - Ga,P и Sb исходное содержание примесей (каждой ) 0,02% (массовых ) Для трех скоростей перемещения зоны V кр =1,5 ; 5 и 15 мм /мин . 2. Проанализировать бинарную диаграмму состояния Si -Ga и представить графически об ласть существования твердых растворов примеси , найти предельную твердую растворимость примеси и температуру предельной растворимости. Рассчитать и построить распределение указанной выше примеси (Ga) в полупроводнике после диффузионного отжига при различ ных условия диффузии : · при условии бесконечного источника примеси на поверхности пластины и при температуре , соответствующей максимальной растворимости примеси в полупроводнике ; время диффузии – 30 мин . · при температуре 950 о С ; время диффузии – 30 мин . · после перераспределения примеси , накопленной в приповерхностном слое полупроводника при температуре 950 о С и времени диффузии – 30 мин . Условия перераспределения - полностью отражающая граница , температура 1150 о С , время 2 часа . Срок сдачи законченной работы руководителю - июнь 1998 г. Преподаватель........................................Б.М . Шишлянников Студент .....................................................Н.Е . Родин Реферат. В курсовом проекте производится расчет распределения примеси вдоль слитка кремния зонной плавкой . Расчет производится для трех примесей ( Ga,P и Sb) для трех скоростей ( V кр =1,5 ; 5 и 15 мм /мин ). Кроме того расчет распределения Ga в кремнии после диффузионного отжига при различных условиях диффузии. Курсовая работа содержит графики распределения примеси как при зонной плавке , так и при диффузии . Содержание. Введение………………………………………………………………… 5 1. Расчетная часть……………………………………………………….. 6 1.1 Описание процесса зонной плавки и ее математ ическая модель…… 6 1.2 Расчет распределения примеси вдоль слитка кремния после зонной плавки (один проход расплавленной зоной )……………………….. 10 1.2.1 Расчет распределения Si-Ga…………………………………………… 10 1.2.2 Расчет распределения Si-P……………………………………………... 13 1.2.3 Расчет распределения Si-Sb……………………………………………. 14 1.3 Распределение примесей после диффузии……………………………. 18 1.3.1 Распределение примеси при диффузии из полубесконечного пространства (диффузия из концентрационного порога )……………….. 21 1.3.2 Распреде ление примеси при диффузии из постоянного источника в полубесконечное тело………………………………………………….. 22 1.3.3 Распределение примеси при диффузии из слоя конечной толщины (диффузия из ограниченного источника ) в полубесконечное тело с отражающей границей……………… …………………………………. 24 1.3.4 Распределение примеси при диффузии из бесконечно тонкого слоя в полубесконечное тело с отражающей границей…………………… 25 1.4 Расчет распределения примеси после диффузионного легирования. 28 1.4.1 Диффузия из бесконечного источник а примеси на поверхности пластины и при температуре , соответствующей максимальной растворимости примеси в полупроводнике ; время диффузии 30мин……… 28 1.4.2 Диффузия из бесконечного источника примеси на поверхности пластины при Т =950 0 С , и времени диффуз ии 30 мин………………. 29 1.4.3 Распределение примеси после перераспределения примеси накопленной в приповерхностном слое полупроводника при Т =950 О С и времени диффузии 30мин . Условие перераспределения полностью отражающая граница . Т =1150 О С , время 2 часа…………… …………. 30 Заключение……………………………………………………………... 32 Литература…………………………………………………………….... 33 Введение. Каждое вещество может находится в состоянии которое характеризуется содержанием примеси в нем ниже некоторого определенного предела . Пред ел определяется различными условиями связанными со свойствами , областью применения веществ . Для полупроводниковых материалов достижения собственных свойств или близких к ним является тем необходимым пределом до которого материалы должны очищаться . При обо с новании необходимой очистки нужно руководствоваться и экономической целесообразности очистки. Для очистки полупроводниковых материалов в технологии микроэлектронных устройств используется метод зонной плавки (перекристаллизация ). В некоторых случаях в техн ологии полупроводниковых материалов выращивают монокристаллы методом зонной плавки . Достоинством метода является совмещение процесса глубокой очистки полупроводника с последующим выращиванием его монокристалла . В технологии разлагающихся полупроводниковых соединений применение этого метода позволяет совмещать в одном технологическом цикле сразу три операции : синтез , очистку синтезированного соединения и выращивание его монокристалла. Для введения в полупроводник примеси используется процесс диффузии . Для из готовления p-n переходов используется химическая диффузия примесных (растворимых ) атомов , которые вводятся в кристаллическую решетку для изменения ее электрофизических свойств . Кроме того диффузия используется для перераспределения примеси в полупроводник е . 1.Расчетная часть. 1.1 Описание процесса зонной плавки и ее математическая модель. Очистку полупроводниковых материалов методом зонной плавки предложил в 1952 году Пфанн . В связи с различной растворимостью примесей в твердой и жидкой ф азах зонная плавка является одним из наиболее эффективных и производительных методов глубокой очистки монокристаллов . При его реализации перед началом кристаллизации расплавляется не вся твердая фаза кристалла (рис .1 ), а только узкая расплавленная зона , к оторую перемещают вдоль слитка. Рисунок 1 – Схема зонной плавки. Большинство примесей обладает хорошей растворимостью в жидкой фазе по сравнению с твердой (равновесный коэффициент сегрегации k 0 < 1), поэтому по мере пр одвижения зона плавления все больше насыщается примесями , которые скапливаются на конце слитка . Обычно процесс зонной плавки повторяют несколько раз , по окончании очистки загрязненный конец слитка отрезают . Для ускорения процесса очистки вдоль контейнера с тавят несколько индукторов для образования ряда зон плавления . Для материалов с k 0 > 1 очистка материалов зонной плавкой практически невозможна. Распределение примесей после одного прохода расплавленной зоной при зонной плавке вдоль слитка представляется уравнением (1) где N тв – концентрация примеси в закристаллизовавшейся фазе на расстоянии x от начала слитка ; N о – исходная концентрация примеси в очищаемом материале ; x – текущая координата (расстояние от начала слитка ); l – длина расплавленной зоны ; k o – равновесный коэффициент распределения. Если измерять длину слитка в длинах расплавленной зоны a = x/l , выражение (1) следует записать иначе : (2) Приведенные уравнения (1) и (2) , являющиеся математическим описанием процесса зонной плавки , в ыведены при определенных допущениях , сформулированных автором метода зонной очистки В . Пфанном при выводе этих уравнений . Эти допущения в литературе принято называть пфанновскими , их суть в следующем : 1 Процессами диффузионного перераспределения компонен тов системы в объеме слитка можно пренебречь , т.е . коэффициенты диффузии компонентов в твердой фазе принимаются равными нулю ( D тв = 0 ). 2 Диффузия компонентов системы в жидкой фазе совершенна - концентрация компонентов постоянна по объему расплава в лю бой момент процесса ; 3 Коэффициент распределения примеси – величина постоянная и не зависит от концентрации примеси в кристаллизующемся веществе (кривые солидус и ликвидус диаграммы состояния прямолинейны ); 4 Начальная концентрация компонентов в исходном материале (слитке ) одинакова по всем сечениям ; 5 Геометрия подвергаемого зонной плавке слитка (длина и по перечное сечение ) в ходе процесса остаются постоянными , плотности твердой и жидкой фаз равны ( r тв =r ж =r). 6 Расплав и твердая фаза при зонной плавке н е взаимодействуют с окружающей средой - атмосферой и контейнером . Другими словами , в системе нет летучих и диссоциирующих компонентов , отсутствует поглощение примесей расплавом из атмосферы , материал контейнера не растворяется в жидкой фазе. Уравнения (1) и (2) справедливы только на участках слитка , на которых зона имеет две границы раздела фаз (постоянный объем ). Когда в системе остается только кристаллизующаяся граница , распределение примеси представляется другим уравнением , соответствующим процессу нор м альной направленной кристаллизации . Другими словами , если длина очищаемого слитка в длинах зон равна A= L/l , то уравнения (1) и (2) справедливы на длине a = (L - l)/l = A-1. При a > A-1 , (3) где g - доля закристаллизовавшегося расплава последнего участка. Только при проведении процесса при условиях , когда удовлетворяются все требования , приведенные выше , реальное распр еделение примеси в слитке после зонной плавки будет соответствовать закону , представленному выражениями (1) и (2). Анализ показывает , что в большинстве реально протекаемых процессов зонной очистки полупроводниковых материалов пфанновские допущения не реали зуются . Вместе с тем , вывод уравнений (1) и (2) без них был бы невозможен , а менее жесткие допущения приводят к существенному усложнению получаемых выражений. Наиболее жесткими являются условия 2 и 3. Допущение 2 в данной формулировке может выполняться тол ько при бесконечно малых скоростях кристаллизации (скорости движения зоны ). В этом случае сравнительно быстрая (по сравнению с диффузией в твердой фазе ) диффузия в жидкой фазе в состоянии постоянно выравнивать концентрации компонентов системы в объеме рас п лавленной зоны. Использовании выражений (1) и (2) для представления распределения примеси при реальных скоростях кристаллизации приводит к необходимости изменить формулировку допущения 2. Выполнение условия постоянства концентрации компонентов по объему р асплава возможно в данной ситуации только при реализации полного ( идеального ) перемешивания жидкой фазы . Предполагается , что в этом случае перераспределение компонентов и выравнивание состава в жидкой фазе происходит мгновенно - т . е . эффективный коэффицие нт диффузии в жидкой фазе D ж = Ґ . Условие полного перемешивания на практике реализовать невозможно . Процессы массопереноса в расплавленной зоне при реальных скоростях кристаллизации и разумной интенсивности перемешивании всегда приводят к образованию ди ффузионного слоя на границе раздела фаз в области кристаллизации . Наличие слоя жидкости с концентрационным пиком , из которого и происходит кристаллизация , влияние его на условия разделения компонентов учитывается введением в выражения (1) и (2) эффектив н ого коэффициента распределения k эфф вместо равновесного k o . Равновесный коэффициент сегрегации связан с эффективным соотношением Бартона-Прима-Слихтера : (4) где V кр - скорость перемещения расплавленной зоны (скорость кристаллизации ); d - толщина диффузионного слоя ; D ж - коэффициент диффузии примеси в жидкой фазе. Эта замена является лишь более или менее удачным приближением к реальной ситуации , и не соответствует требованию условия постоянства концентрации. Распределение примеси после зонной плавки для реальных процессов описываетс я выражением (5) Данное выражение позволяет анализировать влияние на сегрегационные процессы скорости перемещения зоны и услови й перемешивания жидкой фазы. Условие 3 справедливо только для сильно разбавленных растворов , т.е . при малых концентрациях примеси в системе . Кроме того , условие малости концентрации должно соблюдаться на протяжении всего процесса зонной плавки . Для того , ч тобы допущение 3 оказалось состоятельным , требуется использовать при кристаллизационной очистке исходные материалы прошедшие предварительную очистку. 1.2 Расчет распределения примеси вдоль слитка кремния после зонной плавки (один проход расплавленной зон ой ) 1.2.1 Расчет распределения Si-Ga. Рассчитаем распределение галия в слитке кремния для трех скоростей перемещения зоны V кр =1,5 ; 5 и 15 мм /мин. N 0 =0.02% (массовых ). Длина зоны l составляет 10% от длины слитка L. Испарением примеси при переплавке пренебр ечь. Распределение сурьмы вдоль слитка определяется уравнением (5) на длине слитка a = (L - l)/l = A-1, т.е . при 0 Ј a Ј 9. При a > 9 распределение примеси представляется уравнением (3). Доля закристаллизовавшегося расплава g на этом участке изменяется от нуля до величины , близкой к единице . Для g = 1 уравнение (3) не имеет смысла. Прежде чем приступить к расчету переведем N 0 из % (массовых ) в % (атомные ), а затем в см -3 . Для этого воспользуемся формулой перевода. (6) где А 1 , А 2 – атомные массы компонентов ; N 2 – второй компонент смеси. Атомная масса– для галлия = 69,72 [3] – для кремния = 28,08 [3] Концентрация собственных атомов в кристаллической решетке кремния N соб =5 Ч 10 22 см -3 . Следовательно , исходная концентрация галлия в слитке : N 0 =8,06 Ч 10 -5 Ч 5 Ч 10 22 =4,03 Ч 10 18 см -3 Для расчета эффективного коэффициента сегрегации воспользуемся выражением (4). Для галлия в кремнии k 0 =8 Ч 10 -3 [1]. Отношение d /D ж =200 с /см из задания. Подставляя значения k 0 , d /D ж , V кр в (4), вычислим k эфф . Для этого V кр переведем из мм /мин в см /с , получим V кр =2,5 Ч 10 -3 ; 8,33 Ч 10 -3 ; 2,5 Ч 10 -2 см /с . Соответственно получим k эфф =1,3 Ч 10 -2 ; 4,09 Ч 10 -2 ; 0,545 · Заполняем расчетную таблицу , меняя с вы бранным шагом расстояние от начала слитка в длинах зоны a (на участке зонной плавки ). Последний участок слитка , на котором примесь распределяется в соответствии с уравнением (3), разбиваем , меняя расстояние от начала этого участка , пропорционально доле зак ристаллизовавшегося расплава g. · Полученные результаты используются для построения графика распределения примеси N тв вдоль слитка. При построении профиля , как правило , используют полулогарифмический масштаб , т.к . значения концентрации изменяются практичес ки на три порядка. · Определить распределение удельного сопротивления вдоль слитка можно либо расчетным методом , либо по кривым Ирвина. Таблица 1 - Распределение галлия и удельного сопротивления вдоль слитка кремния после зонной плавки (один проход расп лавленной зоной ). Участок зонной плавки Участок направленной кристаллизации а N тв , см -3 r , Ом Ч см (по кривым Ирвина ) g (a=10) N тв , см -3 r , Ом Ч см (по кривым Ирвина ) V кр =2,5 Ч 10 -3 см /с 0 5,24 10 16 0,42 0 4,92 10 17 0,098 1 1,04 10 17 0,28 0,2 6,13 10 17 0,085 2 1,54 10 17 0,21 0,4 8,15 10 17 0,071 3 2,04 10 17 0,18 0,6 1,22 10 18 0,06 4 2,54 10 17 0,15 0,8 2,41 10 18 0,032 5 3,03 10 17 0,14 0,9 4,77 10 18 0,02 6 3,51 10 17 0,13 0,99 4,63 10 19 0,0028 7 3,98 10 17 0,11 – – – 8 4,45 10 17 0,1 – – – 9 4,92 10 17 0,098 – – – V кр =8,33 Ч 10 -3 см /с 0 1,6 10 17 0,2 0 1,35 10 18 0,05 1 3,2 10 17 0,135 0,2 1,67 10 18 0,048 2 4,68 10 17 0,098 0,4 2,2 10 18 0,036 3 6,11 10 17 0,085 0,6 3,25 10 18 0,028 4 7,48 10 17 0,075 0,8 6,32 10 18 0,017 6 1,0 10 18 0,0 61 0,9 1,23 10 19 0,009 7 1,13 10 18 0,055 0,99 1,12 10 20 0,0011 8 1,24 10 18 0,051 – – – 9 1,35 10 18 0,05 – – – V кр =2,5 Ч 10 -2 см /с 0 2,2 10 18 0,036 0 4,02 10 18 0,0215 1 2,97 10 18 0,029 0,2 4,45 10 18 0,021 2 3,41 10 18 0,025 0,4 5,07 10 18 0,019 3 3,67 10 18 0,023 0,6 6,1 10 18 0,017 4 3,82 10 18 0,0225 0,8 8,36 10 18 0,0125 6 3,96 10 18 0,0222 0,9 1,15 10 19 0,01 7 3,98 10 18 0,022 0,99 3,27 10 19 0,0037 8 4,01 10 18 0,0215 – – – 9 4,02 10 18 0,0215 – – – 1.2.2 Ра счет распределения Si-P. Расчет распределения фосфора в кремнии будем производить аналогично расчету галлия в слитке кремния (пункт 1.2.1), при тех же условиях зонной плавки . Переведем N 0 в см -3 . Атомная масса фосфора = 30,97 N 0 =0,02 % (массовых ) = 1,81 Ч 10 -2 % (атомных ) = 9,05 Ч 10 18 см -3 . Для расчета эффективного коэффициента сегрегации k эфф воспользуемся выражением (4). Для фосфора в кремнии k 0 =3,5 Ч 10 -1 [1]. Отнош ение d /D ж =200 с /см из задания. Подставляя значения k 0 , d /D ж , V кр в (4), вычислим k эфф . Для трех скоростей кристаллизации V кр =2,5 Ч 10 -3 ; 8,33 Ч 10 -3 ; 2,5 Ч 10 -2 см /с соответственно получим k эфф =0,47; 0,74; 0,99. Заполним расчетную таблицу. Таблица 2 - Распределение фосфора и удельного сопротивления вдоль слитка кремния после зонной плавки (один проход расплавл енной зоной ). Участок зонной плавки Участок направленной кристаллизации а N тв , см -3 r , Ом Ч см (по кривым Ирвина ) g (a=10) N тв , см -3 r , Ом Ч см (по кривым Ирвина ) V кр =2,5 Ч 10 -3 см /с 0 4,25 10 18 0,012 0 9 10 18 0,0068 1 6,05 10 18 0,009 0,2 1,01 10 19 0,006 2 7,18 10 18 0,0075 0,4 1,18 10 19 0,0058 3 7,88 10 18 0,0073 0,6 1,46 10 19 0,0042 4 8,32 10 18 0,0071 0,8 2,11 10 19 0,0034 5 8,6 10 18 0,007 0,9 3,05 10 19 0,0024 6 8,76 10 18 0,0069 0,99 1,03 10 20 0,00085 7 8,87 10 18 0,0069 – – – 9 9 10 18 0,0068 – – – V кр =8,33 Ч 10 -3 см /с 0 6,69 10 18 0,0085 0 9,05 10 18 0,0066 1 7,93 10 18 0,0075 0,2 9,59 10 18 0,0063 2 8,51 10 18 0,0071 0,4 1,03 10 19 0,006 3 8,8 10 18 0,0069 0,6 1,15 10 19 0,0057 4 8,93 10 18 0,0068 0,8 1,38 10 19 0,0052 5 9 10 18 0,0068 0,9 1,65 10 19 0,0045 7 9,03 10 18 0,0066 0,99 3 10 19 0,0024 8 9,04 10 18 0,0066 – – – 9 9,05 10 18 0,0066 – – – V кр =2,5 Ч 10 -2 см /с 0 8,96 10 18 0,0068 0 9,05 10 18 0,0066 1 9,01 10 18 0,0068 0,2 9,07 10 18 0,0066 2 9,03 10 18 0,0066 0,4 9,1 10 18 0,0065 3 9,05 10 18 0,0066 0,6 9,13 10 18 0,0065 4 – – 0,8 9,2 10 18 0,0064 5 – – 0,9 9,26 10 18 0,0064 9 9,05 10 18 0,0066 0,99 9,48 10 18 0,0063 1.2.3 Расчет распределения Si-Sb. Расчет распределения сурьмы в кремнии будем производить аналогично расчету галлия в слитке кре мния (пункт 1.2.1), при тех же условиях зонной плавки . Переведем N 0 в см -3 . Атомная масса сурьмы = 121,7 N 0 =0,02 % (массовых ) = 4,62 Ч 10 -3 % (атомных ) = 2,31 Ч 10 18 см -3 . Для расчета эффективного коэффициента сегрегации k эфф воспользуемся выражением (4). Для сурьмы в кремнии k 0 =2,3 Ч 10 -3 [1]. Отношение d /D ж =200 с /см из задания. Подставляя значения k 0 , d /D ж , V кр в (4), вычислим k эфф . Для трех скоростей кристаллизации V кр =2,5 Ч 10 -3 ; 8,33 Ч 10 -3 ; 2,5 Ч 10 -2 см /с соответственно получим k эфф =3,74 Ч 10 -2 ; 0,11; 0,78. Заполним расчетную таблицу. Таблица 3 - Распределение сурьмы и удельного сопротивления вдоль слитка кремния после зонной плавки (один проход расплавленной зоной ). Участок зонной плавки Участок направленной Кристаллизации а N т в , см -3 r , Ом Ч см (по кривым Ирвина ) g (a=10) N тв , см -3 r , Ом Ч см (по кривым Ирвина ) V кр =2,5 Ч 10 -3 см /с 0 8,64 10 16 0,11 0 7,22 10 17 0,028 1 1,68 10 17 0,075 0,2 8,95 10 17 0,023 2 2,47 10 17 0,052 0,4 1,18 10 18 0,0215 3 3,22 10 17 0,047 0,6 1,74 10 18 0,0192 4 3,95 10 17 0,04 0,8 3,4 10 18 0,014 5 4,66 10 17 0,038 0,9 6,62 10 18 0,0082 6 5,33 10 17 0,033 0,99 6 10 19 0,00135 7 6 10 17 0,031 – – – 9 7,22 10 17 0,028 – – – V кр =8,33 Ч 10 -3 см /с 0 2,54 10 17 0,051 0 1,55 10 18 0,02 1 4,68 10 17 0,038 0,2 1,89 10 18 0,018 2 6,6 10 17 0,03 0,4 2,44 10 18 0,016 3 8,32 10 17 0,027 0,6 3,5 10 18 0,013 4 9,86 10 17 0,024 0,8 6,49 10 18 0,0085 5 1,12 10 18 0,022 0,9 1,2 10 19 0,0055 7 1,36 10 18 0,0205 0,99 9,3 10 19 0,00088 8 1,46 10 18 0,02 – – – 9 1,55 10 18 0,02 – – – V кр =2,5 Ч 10 -2 см /с 0 1,8 10 18 0,019 0 2,31 10 18 0,0157 1 2,08 10 18 0,017 0,2 2,42 10 18 0,0156 2 2,2 10 18 0,016 0,4 2,58 10 18 0,015 3 2,26 10 18 0,0158 0,6 2,82 10 18 0,014 5 2,3 10 18 0,0157 0,8 3,29 10 18 0,0137 7 2,31 10 18 0,0157 0,9 3,83 10 18 0,012 9 2,31 10 18 0,0157 0,99 6,36 10 18 0,0086 На основании полученных данных построим графики распределения примесей вдоль слитка кремния после зонной плавки (один проход расплавленной зоной ). 1.3. Распределение примесей после диффузии. Основой математического описания процессов диффузии являются два дифференциальных уравнения Фика (немецкий ученый A. Fick предложил их в 1855 г .). Первое уравнение (первый закон Фика ) записывается следу ющим образом : J = - D grad N (7) где J - плотность потока диффундирующего вещества , т.е . количество вещества , проходящего за единицу времени через единичную площадь поверхности , перпендикулярной направлению переноса вещества ; N - концентрация атомов примеси. D - коэффициент диффузии. Физический смысл этого уравнения — первопричиной диффузионного массопереноса веще ства является градиент его концентрации . Скорость переноса пропорциональна градиенту концентрации , а в качестве коэффициента пропорциональности вводится коэффициент диффузии . Знак минус в правой части (7) указывает на то , что диффузия происходит в направ л ении убывания концентрации . Другими словами , диффузия идет благодаря стремлению системы достичь физико-химического равновесия . Процесс будет продолжаться до тех пор , пока химические потенциалы компонентов всей системы не станут равными . Уравнение (7) опис ы вает стационарный (установившийся ) процесс - процесс , параметры которого не зависят от времени. В макроскопическом представлении коэффициент диффузии определяет плотность потока вещества при единичном градиенте концентрации и является , таким образом , меро й скорости выравнивания градиента концентрации . Размерность коэффициента диффузии - м 2 /с . В общем случае диффузия анизотропна и коэффициент диффузии - симметричный тензор второго ранга. Согласно микроскопическому определению , компонента D x коэффициента д иффузии D по координате x связана со среднеквадратичным смещением диффундирующих атомов по координате x и интервалом времени D t , в течение которого это смеще ние произошло соотношением Когда концентрация вещества изменяется только в одном направлении (одномерная диффузия ) и при диффузии в изотропной среде (коэффиц иент диффузии - скаляр ) первое уравнения Фика имеет следующий вид : (8) При простейшем анализе структур и в простейших моделях процессов легирования в технологии изготовления ИМС предполагаются именно такие условия диффузии. Второе уравнение диффузии (второй з акон Фика ) получается путем сочетания первого закона и принципа сохранения вещества , согласно которому изменение концентрации вещества в данном объеме должно быть равно разности потоков этого вещества на входе в объем и выходе из него. В общем случае второ е уравнение диффузии имеет следующий вид (9) Для одномерной диффузии в изотропной среде уравнение (9) м ожно записать (10) Второй закон Фика характеризует процесс изменения концентрации диффундирующей примеси во времени в различных точках среды и является математической моделью нестационарного (развивающегося ) состояния системы (описывает период времени от начала процесса до установления стационарного состояния ). При постоянстве коэффициента диффузии D (независимости его от концентрации примеси ) уравнение (10) упрощается (11) Допущение о по стоянстве коэффициента диффузии справедливо в большом количестве случаев , анализируемых в технологии ИМС. Уравнения диффузии являются чисто феноменологическими , т.е . они не содержат никаких сведений о механизмах диффузии - о диффузионном процессе на атомно м , уровне . Кроме того , уравнения (7) - (11) не содержат информации о зарядовом состоянии диффундирующих частиц. Процессы диффузии , используемые для изготовления интегральных структур , обычно анализируются с помощью частных решений уравнения (11) т.к ., в о тличие от (8), именно оно содержит важный параметр - время установления некоторого анализируемого состояния системы . Основная цель решения уравнения - найти распределение примеси N(x,t) в полупроводнике после диффузии в течение определенного времени t при различных условиях осуществления процесса. Общее решение уравнения (11) для бесконечного твердого тела при заданном в общем , виде начальном распределении примеси N(x,0) = f(x) может быть найдено методом разделения переменных . Оно имеет вид , (12) здесь x - текущая координата интег рирования. Бесконечным в одномерном представлении называют тело , простирающееся от x=0 до x=- Ґ и до x=+ Ґ . Часто при поиске распределения концентрации примеси в полупроводнике необходимо решение уравнения (11) для полубесконечного твердого тела . Полубеск онечным в одномерном представлении называют тело , простирающееся от x=0 до x=+ Ґ . Для этого случая выражение (12) может быть приведено к виду (13) В выражении (13) знак плюс относится к ситуации , когда граница твердого тела (x=0) является непроницаемой для диффундирующего вещества , находящегося в области x>0, ( отражающая граница ), а знак минус - к случаю , когда на границе твердого тела в любой момент времени концентрация диффундирующего вещества , также находящегос я в области x>0, равна нулю - связывающая граница . Представленные решения позволяют находить распределения примеси в твердом теле при любых начальных условиях . Решение конкретной задачи сводится к подстановке в (12) или (13) соответствующих ситуации начал ьных условий с последующими , как правило , очень громоздкими преобразованиями. 1.3.1 Распределение примеси при диффузии из полубесконечного пространства (диффузия из концентрационного порога ) Диффундирующая примесь ( диффузант ) поступает в полубесконечное т ело через плоскость x=0 из второго полубесконечного тела ( источника ) с равномерным распределением примеси . Концентрация примеси в источнике - N o . Полагается , что в принимающем диффузант теле нет рассматриваемой примеси . Начальное распределение концентраци й для этого случая задается в виде N(x,0) = N o для x<0 N(x,0) = 0 для x>0 Решением уравнения (11) для этого случая является выражение (14) Второе слагаемое в квадратных скобках называют интегралом ошибок Гаусса или , иначе , функцией ошибок - error function и сокращенно обозначают erf (z). В соответствии с сокр ащением это распределение называют erf - распределением . (15) В математике часто используют как самостоятель ную и другую функцию erfc z = 1- erf z (16) которая называется дополнением функции ошибок до единицы или дополнительной функцией ошибок - error function complement. Обе функции табулированы . Таким образом , выра жение (14) можно записать (17) Величина имеет размерность длины и носит название диффузионной длины или длины диффузии . Физический смысл этого параметра - среднее расстояние , которое преодолели диффундирующие частицы в направлении выравнивани я градиента концентрации за время t. Рассмотренное решение можно использовать как простейшую модель , представляющую распределение примеси в автоэпитаксиальной структуре . При этом , в качестве независимых источников примеси выступает как подложка , так и эпи таксиальный слой . Процессы диффузии с каждой стороны рассматриваются в этом случае как независящие друг от друга , а реальное распределение примесей на границе раздела будет представлять собой сумму отдельных решений. 1.3.2 Распределение примеси при диффу зии из постоянного источника в полубесконечное тело. Диффузант поступает в полубесконечное тело через плоскость x=0 из источника , обеспечивающего постоянную концентрацию примеси N o на поверхности раздела твердое тело - источник в течение любого времени . Та кой источник называют бесконечным или источником бесконечной мощности . Полагается , что в принимающем диффузант теле нет рассматриваемой примеси . Начальное распределение концентраций и граничные условия для этого случая задаются в виде N(x,t) = N o для x=0 N(x,0) = 0 для x>0 Решением уравнения (16) для данных условий является выражение (18) Если в объеме по лупроводникового материала до диффузии имелась примесь противоположного типа по отношению к диффундирующей , эта примесь распределена по объему равномерно и её концентрация равна N b , то в этом случае в полупроводнике образуется электронно-дырочный переход . Его положение ( глубина залегания ) x j определяется условием N(x,t)=N b , откуда (19) и (20) здесь запись erfc -1 обозначает аргумент z функции erfc . Пр и решении практических задач , связанных с анализом диффузионных процессов необходимо знать количество примеси Q , накопленной в твердом теле при диффузии в течение времени t . Эта величина определяется по формуле (21) где J (0,t) - поток диффузанта в объем через плоскость x=0 (22) отсюда (23) Следует обратить внимание на возрастающее со вр еменем значение накопленной в диффузионном слое примеси при диффузии с данными граничными условиями. Рассмотренная модель диффузионного процесса с постоянным источником описывает процесс диффузионного легирования полупроводникового материала из газовой или паровой фазы . Этот процесс используется при создании сильно легированных диффузионных слоев (например , эмиттерных ) с поверхностными концентрациями N o близкими к значениям предельной твердой растворимости примеси в данном полупроводниковом материале . Твер дое тело можно считать полубесконечным ( или бесконечным ) в том случае , если его размеры в направлении движения диффузанта много больше длины диффузии. 1.3.3 Распределение примеси при диффузии из слоя конечной толщины (диффузия из ограниченного источника ) в полубесконечное тело с отражающей границей. Диффундирующая примесь поступает в полубесконечное тело из источника , который представляет собой примыкающий к границе тела слой толщиной h , примесь в котором распределена равномерно . Такой источник называют о граниченным . Концентрация примеси в источнике - N o . Полагается , что в принимающем диффузант твердом теле нет рассматриваемой примеси . При абсолютно непроницаемой для диффузанта (отражающей ) границе поток примеси через поверхность x=0 должен обращаться в н уль при всех t і 0 для t і 0 (24) Начальное распределение концентраций для рассматриваемого случая задаётся в виде N( x,0) = N o для 0 Ј x Ј h N(x,0) = 0 для x>h Граничным условием является , определяемое условием (24), постоянство количества примеси в источнике и полупроводнике Для реализации начального распределения такого типа диффундирующая примесь должна быть введена в твердое тело до начала диффузии . Решением уравнения (16) в данной ситуации является выражение (25) Здесь следует отметить , что erfс (-z) + erfс (z) є 2. В отличие от диффузии из постоянного источника при диффузии из слоя конечной толщины количество диффузанта ограничено значение м Q=N o h . В процессе диффузии происходит только его перераспределение и , следовательно , уменьшение со временем концентрации примеси на поверхности твердого тела. Примером диффузии примеси из слоя конечной толщины в полубесконечное тело с отражающей границей является диффузия в кремниевую пластину из эпитаксиального , имплантированного или диффузионного слоя и покрытую слоем двуокиси кремния SiO 2 или нитрида кремния Si 3 N 4 . Границу пластины и пленки можно с большой долей правдоподобия принять отражающей , т.к . к оэффициенты диффузии большинства примесей в кремнии на несколько порядков больше , чем в двуокиси кремния и нитриде . Однако , равномерность распределения примеси в источнике , особенно при его создании методом диффузии или имплантации - весьма грубое и вынуж д енное приближение. 1.3.4 Распределение примеси при диффузии из бесконечно тонкого слоя в полубесконечное тело с отражающей границей Решение диффузионного уравнения при этих условиях находится из предыдущего при h ® 0 и условии , что количество диффузанта в источнике Q=N o h. (26) Приведенное выражение представляет собой Гауссово распределение. Тонкий слой на поверхности по лупроводниковой пластины является источником , который очень быстро истощается . Непрерывная диффузия в этом случае приводит к постоянному понижению поверхностной концентрации примеси в полупроводнике . Эту особенность данного процесса используют в полупрово д никовой технологии для получения контролируемых значений низкой поверхностной концентрации примеси , например , для создания базовых областей кремниевых транзисторных структур дискретных приборов или ИМС. На первом этапе процесса проводится кратковременная д иффузия (при пониженных температурах ) из постоянного источника , распределение примеси после которой описывается выражением (18). Значение N o при этом велико и определяется либо пределом растворимости данной примеси в полупроводниковом материале , либо конце нтрацией примеси в стеклообразном слое на поверхности полупроводника . Этот этап часто называют загонкой. После окончания первого этапа пластины помещают в другую печь для последующей диффузии , обычно , при более высоких температурах . В этой печи нет источни ка примеси , а если он создается на первой стадии в виде стеклообразного слоя на поверхности пластин , его предварительно удаляют . Таким образом , тонкий слой , полученный на первом этапе , является источником перераспределяемой примеси при проведении второй с т адии процесса . Для создания отражающей границы второй этап (часто называемый разгонкой ) проводят в окислительной атмосфере . При этом на поверхности растет слой SiO 2 . Существует заметное несоответствие между распределением примеси в источнике , сформирован ном при загонке , с декларируемым при выводе выражений (25) и (26) - ступенчатым . Это несоответствие должно отразиться на точности описания реального распределения примеси после второй стадии диффузии выражением (26). Не существует и объективного количес т венного критерия “ тонкости ” источника — нет каких-либо признаков , согласно которым для представления результатов данного процесса следует использовать выражение (26), а на (25) и наоборот. При моделировании двухстадийной диффузии и анализе результатов пр оцесса полагают , что выражение (26) достаточно точно соответствует реальному при условии , если величина произведения D 1 t 1 для первого этапа процесса легирования значительно меньше , чем D 2 t 2 для второго - . Это условие быстрой истощаемости источника . В этом случае , учитывая , что количество накопленной при первом этапе примеси определяется соотношением из (26) получим (27) Величины D 2 и t 2 относятся ко второй стадии диффузии. В случае , если продолжительность второй стадии не очень велика по сравнению с первой , или , иными словами , D 2 t 2 і D 1 t 1 , предположение о том , что диффузионный слой , образовавшийся в результате загонки , будет вести себя как тонкий источник неверно . В этом случае решение диффузионного уравнения будет выглядеть следующим образом (28) где и Поверхностная концентрация примеси после второй стадии диффузии выражается при данных условиях соотношением (29) Выражение (25) используется для представления распределения при условии , что D 1 t 1 >D 2 t 2 – . При этом полагают , что . 1.4 Расчет распределения примеси после диффузионного легирования . 1.4.1 Дифф узия из бесконечного источника примеси на поверхности пластины и при температуре , соответствующей максимальной растворимости примеси в полупроводнике ; время диффузии 30 мин .= 1800с. Материал – кремний ; примесь – галлий. Условия проведения диффузии соответ ствуют решению , представляемому уравнением (18). . Температуру соответствующую максимальной растворимости галлия в кремнии , а так же и саму предельную растворимость найдем из графика предельной растворимости примеси в кремнии. N пред . раств. =N 0 =6 Ч 10 19 см -3 , Т =1523 К. Коэффициент диффузии сурьмы при температуре диффузии найдем используя известное выражение в форме уравнения Аррениуса , где предэкспоненциальный мн ожитель (постоянная диффузии ) D o и энергия активации диффузии DE - справочные величины . k - постоянная Больцмана , T - температура процесса в Кельвинах. Из [5] для галлия D o =0,374 см 2 /с , DE = 3,41 эВ , при T = 1523 K D = 1,94 Ч 10 -12 см 2 /с. Заполняем расчетную таблицу , меняя расстояние от x поверхности с необходимой частотой , до значения при котором значение N(x) имеет порядок не более 10 12 . В первый столбец записываем выбранные значения x , во второй - . Затем находим значения erfc(z), воспользовавшись таблицей интеграла ошибок в [4], и вносим эти значения в третий столбец . После чего рассчитываем концентрации N(x), соответствующие каждому значению x и результаты записываем в четвертый столбец. Таблица 4 - Результаты расчета распределения галлия в кремнии x, мкм erfc(z) N(x), см -3 x, мкм erfc(z) N(x), см -3 0 0 1 6 Ч 10 19 2,5 2,12 0,002716 1,63 Ч 10 17 0,5 0,42 0,552532 3,32 Ч 10 19 3 2,54 0,000328 1,968 Ч 10 16 1 0,85 0,229332 1,376 Ч 10 19 3,5 2,96 0,000028 1,68 Ч 10 15 1,5 1,27 0,072486 4,35 Ч 10 18 4 3,38 0,000001753 1,052 Ч 10 14 2 1,69 0,016847 1,012 Ч 10 18 4,5 3,8 0,000000077 4,62 Ч 10 12 Полученные результаты используются для построени я графика N = f(x) - примесного профиля . При построении профиля , как правило , используют полулогарифмический масштаб. 1.4.2 Диффузия из бесконечного источника примеси на поверхности пластины при Т =950 0 С =1223 К , и време ни диффузии 30 мин .=1800 с. Коэффициент диффузии галлия в кремнии при Т =950 0 С , N 0 = 3 Ч 10 19 см -3 . Диффузия проходит согласно выражению (18).Дальнейший ход работы идет аналогично пункту 1.4.1. Заполняем расчетную таблицу. Т аблица 5 - Результаты расчета распределения галлия в кремнии x, мкм erfc(z) N(x), см -3 x, мкм erfc(z) N(x), см -3 0 0 1 3 Ч 10 19 0,1 2,05 0,003742 1,123 Ч 10 17 0,02 0,41 0,562031 1,6861 Ч 10 19 0,12 2,46 0,000503 1,50910 16 0,04 0,82 0,246189 7,386 Ч 10 18 0,14 2,87 0,000049 1,47 Ч 10 15 0,06 1,23 0,08195 2,4585 Ч 10 18 0,16 3,28 0,0000035 1,05 Ч 10 14 0,08 1,64 0,020378 6,1134 Ч 10 17 0,18 3,69 0,00000018 5,4 Ч 10 12 Полученные результаты используются для построения графика N = f(x) - примесного профиля . При построении профиля , как правило , используют пол улогарифмический масштаб. 1.4.3 Распределение примеси после перераспределения примеси накопленной в приповерхностном слое полупроводника при Т =950 О С =1223 К и времени диффузии 30мин =1800с . Условие перераспределения полностью отражающая граница. Т =1150 О С =1423 К , время 2 часа =7200с. Произведение D 1 t 1 для процесса загонки равно : D 1 t 1 = 3,31 Ч 10 -15 Ч1800= 5,958 Ч 10 -12 см 2 Коэффициент диффузии для процесса перераспределения примеси ( D o =0,374 см 2 /с , DE = 3,41 эВ , T = 1423 K) равен D = 3,128 Ч 10 -13 с м 2 /с . Произведение D 2 t 2 = 3,128 Ч 10 -13 Ч7200= 2,25 Ч 10 -9 см 2 . D 2 t 2 > D 1 t 1 (в 377 раз ), т.е . условия быстрой истощаемости источника , следовательно , пользуемся для расчета распределения примеси выражением (27). В первый столбец таблицы (6) заносим значения x , во второй значения exp(-x 2 /4D 2 t 2 ) , рассчитанные значения N s заносим в третий столбец. Таблица 6 - Результаты расчета распределения галлия в кремнии при диффузии из приповерхностного слоя. x, мкм N s , см -3 x, мкм N s , см -3 0 1,0 9,823 Ч 10 17 1,8 0,02742 2,693 Ч 10 16 0,36 0,866 8,507 Ч 10 17 2,16 0,005633 5,534 Ч 10 15 0,72 0,5624 5,525 Ч 10 17 2,52 0,0008681 8,527 Ч 10 14 1,08 0,274 2,69 Ч 10 17 2,88 0,0001 9,854 Ч 10 13 1,44 0,1 9,831 Ч 10 16 3,24 0,0000087 8,541 Ч 10 12 Полученные результаты используются для построения графика N = f(x) - примесного профиля. Заключение. В данном курсовом проекте были рассмотрены процесс очистки полупроводникового вещества – зонная пла вка и способ введения примеси в полупроводник – диффузия примеси. Для процесса зонной плавки произведен расчет для трех очищаемых примесей : фосфор , галлий , сурьма . Результаты расчета представлены в виде таблиц и графиков : распределение удельного сопротивл ения и распределения каждой примеси вдоль слитка кремния после очистки зонной плавкой (один проход расплавленной зоной ). Эффективность очистки зависит от скорости кристаллизации : чем меньше скорость кристаллизации в донной примеси , тем лучше она очищает ся , таким образом при V кр ® 0 k эфф ® k 0 ; V кр ® Ґ k эфф ® 1. Но это не означает , что если мы уменьшим скорость кристаллизации до нуля , то получим исходное вещество в чистом виде – это лишь одно из условий очистки вещества . Определяющим является также равновесный коэффициент сегрегации (К 0 ) , который отражает эффективность перераспределения между жидкой и твердой фазой , он должен отличаться от еденицы в большую или меньшую сторону . В нашем случае k 0 Sb
1Архитектура и строительство
2Астрономия, авиация, космонавтика
 
3Безопасность жизнедеятельности
4Биология
 
5Военная кафедра, гражданская оборона
 
6География, экономическая география
7Геология и геодезия
8Государственное регулирование и налоги
 
9Естествознание
 
10Журналистика
 
11Законодательство и право
12Адвокатура
13Административное право
14Арбитражное процессуальное право
15Банковское право
16Государство и право
17Гражданское право и процесс
18Жилищное право
19Законодательство зарубежных стран
20Земельное право
21Конституционное право
22Конституционное право зарубежных стран
23Международное право
24Муниципальное право
25Налоговое право
26Римское право
27Семейное право
28Таможенное право
29Трудовое право
30Уголовное право и процесс
31Финансовое право
32Хозяйственное право
33Экологическое право
34Юриспруденция
 
35Иностранные языки
36Информатика, информационные технологии
37Базы данных
38Компьютерные сети
39Программирование
40Искусство и культура
41Краеведение
42Культурология
43Музыка
44История
45Биографии
46Историческая личность
47Литература
 
48Маркетинг и реклама
49Математика
50Медицина и здоровье
51Менеджмент
52Антикризисное управление
53Делопроизводство и документооборот
54Логистика
 
55Педагогика
56Политология
57Правоохранительные органы
58Криминалистика и криминология
59Прочее
60Психология
61Юридическая психология
 
62Радиоэлектроника
63Религия
 
64Сельское хозяйство и землепользование
65Социология
66Страхование
 
67Технологии
68Материаловедение
69Машиностроение
70Металлургия
71Транспорт
72Туризм
 
73Физика
74Физкультура и спорт
75Философия
 
76Химия
 
77Экология, охрана природы
78Экономика и финансы
79Анализ хозяйственной деятельности
80Банковское дело и кредитование
81Биржевое дело
82Бухгалтерский учет и аудит
83История экономических учений
84Международные отношения
85Предпринимательство, бизнес, микроэкономика
86Финансы
87Ценные бумаги и фондовый рынок
88Экономика предприятия
89Экономико-математическое моделирование
90Экономическая теория

 Анекдоты - это почти как рефераты, только короткие и смешные Следующий
В Москве начали закрывать Макдональдсы за экстремистский лозунг: "Свободная касса".
Anekdot.ru

Узнайте стоимость курсовой, диплома, реферата на заказ.

Обратите внимание, курсовая по радиоэлектронике "Расчет распределения примесей в кремнии при кристаллизационной очистке и диффузионном легировании", также как и все другие рефераты, курсовые, дипломные и другие работы вы можете скачать бесплатно.

Смотрите также:


Банк рефератов - РефератБанк.ру
© РефератБанк, 2002 - 2016
Рейтинг@Mail.ru