Моделирование кремниевого транзистора со статической индукцией Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

2005

Доклады БГУИР

апрель- июнь

№ 2

УДК[621.382.3: 546.28]: 537.852

МОДЕЛИРОВАНИЕ КРЕМНИЕВОГО ТРАНЗИСТОРА СО СТАТИЧЕСКОЙ ИНДУКЦИЕЙ

Н.Л. ДУДАР

УП "БелМикроСистемы"НПО "Интеграл" Минск, Республика Беларусь

Поступила в редакцию 25 ноября 2004

Моделирование полупроводниковых структур позволяет снизить затраты на изготовление экспериментальных образцов и оптимизировать параметры элементов интегральных схем и приборов. В данной работе представлены результаты моделирования кремниевого транзистора со статической индукцией. Транзисторы со статической индукцией (СИТ) представляют собой структуры с коротким каналом, работающие на полевом эффекте и пригодные для работы в условиях высокой мощности, высокой температуры и высокой частоты. СИТ, для которого проводился расчет, изготовлен на кремниевой пластине с эпитаксиальной пленкой. В эпитаксиальной пленке п-типа были сформированы области базы р-типа и эмиттера п+-типа. Маршрут изготовления СИТ также включает в себя выращивание локального окисла, который обеспечивает надежную изоляцию р - и п+-областей транзистора друг от друга. Была разработана одномерная модель транзистора со статической индукцией, которая использовалась для исследований его электрических характеристик. Модель описывается с помощью системы уравнений полупроводника, решение которой в одномерном случае не требует значительных затрат машинного времени. С использованием результатов моделирования технологического маршрута изготовления были рассчитаны распределения плотностей токов и вольт-амперные характеристики СИТ. Моделирование технологического маршрута было выполнено с применением программного средства 8иРРБМ3. Программа 8иРКБМ3 создает файл выходных данных, содержащий таблицы числовых значений моделируемых параметров структуры (распределения концентраций легирующих примесей и подвижных носителей заряда, потенциала и т.п.). Расчеты распределений токов по глубине структуры и расчеты вольт-амперных характеристик прибора проводились численным методом Гуммеля.

Ключевые слова: моделирование, транзистор со статической индукцией, численный метод, технологический маршрут.

Введение

Проектирование и моделирование полупроводниковых приборов является важным и необходимым этапом разработки и изготовления новых изделий на современных предприятиях электронной промышленности [1]. Моделирование полупроводниковых структур позволяет снизить затраты на проведение экспериментальных исследований и оптимизацию параметров элементов интегральных схем и приборов. Одним из средств моделирования технологии изготовления полупроводниковых структур является программный пакет 8иРКЕМ3, позволяющий получить распределения концентраций легирующих примесей, подвижных носителей заряда и потенциала, которые формируются при заданном технологическом маршруте изготовления из-

делия. Результаты моделирования технологии являются входными данными для моделирования электрических характеристик прибора.

В настоящей работе приводятся результаты исследований электрических характеристик (распределений плотностей токов по глубине структуры и вольт-амперных характеристик), полученных с использованием разработанной автором одномерной модели транзистора со статической индукцией (СИТ). СИТ используется как дискретный прибор в условиях высокой мощности, высокой температуры и высокой частоты. Модель СИТ описывается с помощью системы уравнений полупроводника, решение которой в одномерном случае не требует значительных затрат машинного времени. Полученные в программе для технологического моделирования результаты далее используются в качестве исходных данных для расчета прибора численным методом [2].

Модель транзистора со статической индукцией

Транзисторы со статической индукцией представляют собой структуры с коротким каналом, работающие на полевом эффекте. Конструкция структуры СИТ показана на рис. 1. Такая структура формируется на кремниевой пластине с эпитаксиальным слоем. Область р-базы формируется путем разгонки примеси бора при температуре 1200°С после ионного легирования. Области р-базы и р-омических контактов формируются посредством ионного легирования бором через маску нитрида кремния. п -область эмиттера формируется посредством ионного легирования мышьяком в свободную от локального окисла поверхность кремния. Локальный окисел обеспечивает надежную изоляцию р - и п+-областей друг от друга. После формирования областей р-базы и пэмиттера наносится слой поликремния. Роль коллектора в описываемой структуре выполняют эпитаксиальный слой и подложка.

Металл

Коллектор

Рис. 1. Структура транзистора со статической индукцией

В данной работе электрические характеристики СИТ рассчитываются с использованием разработанной автором одномерной модели транзистора.

В основе используемой в расчетах программы для моделирования переноса носителей заряда лежат:

уравнение Пуассона

д2ф- 4 ( - п + N. - Ма), (1)

дх2 вво ё а

где ф — электростатический потенциал; х — координата; р, п — концентрации электронов и дырок соответственно; ц — заряд электрона; в — относительная диэлектрическая проницаемость; во — диэлектрическая проницаемость вакуума; N — концентрация доноров; Ыа — концентрация акцепторов;

уравнения непрерывности для электронов и дырок:

^-цЯ-4%-0, (2)

д х д ^

где Jп — электронная составляющая тока; I — время; Я — скорость рекомбинации носителей заряда;

Цх~ + цЯ + ц др = 0, (3)

д х д ^

где Jp — дырочная составляющая тока.

Предполагается, что механизм рекомбинации как в объеме, так и на поверхности полупроводника описывается уравнением Шокли-Рида-Холла [2]:

пр - и2

Я ----, (4)

т (п - п.) + т (р - п.) р - п -

где п- — собственная концентрация носителей заряда в полупроводнике; тп, тр — времена жизни электронов и дырок соответственно.

Решение системы уравнений (1)-(4) осуществляется методом Гуммеля [3]. Система уравнений (1)-(4) может быть решена только численным методом. На каждой итерации численной реализации метода Гуммеля выражения (1)-(4) в дискретной форме решаются методом Гаусса [4, 5].

Результаты моделирования технологического маршрута изготовления СИТ

Для структуры кремниевого транзистора со статической индукцией был разработан технологический маршрут, включающий: 1) изготовление эпитаксиальной пленки п-типа на подложке п-типа; 2) формирование области р-базы и р -областей омических контактов к области р-базы путем ионного легирования и последующей разгонки примеси бора; 3) формирование п-областей эмиттера; 4) формирование локального окисла; 5) нанесение слоя поликремния. Роль коллектора в сформированной структуре выполняют эпитаксиальный слой п-типа и подложка. Моделирование выполнялось в одномерном пространстве вдоль координаты Х (рис. 1). Как видно из рисунка, вдоль указанного направления Х чередуются области типа п , р, п, п .

Распределения концентраций подвижных носителей заряда показаны на рис. 2.

И я

I*

с а. к и s «

Л а>

Ё Р 55

£ О О X Ы

1.00Е+22 1.00Е+20 1.00Е+18 1.00Е+16 1.00Е+14 1.00Е+12 1.00Е+10 1,00Е+08 1.00Е+06 1.00Е+О4 1.00Е+02 1,00Е+00 1.00Е-02

Концентрация электронов

лк а л к к -к—к—

/ —♦—Концентрация дьрок

Г

ШБ)\ (К)

L Л

ШШ

10 15 20 25

Глубина структуру мкм

30

35

Рис. 2. Распределения концентраций подвижных носителей заряда по структуре СИТ, полученные при моделировании на SUPREME

Как видно из рис. 2, глубина залегания рассчитанного р-«-перехода эмиттер-база в моделируемой структуре составила ~2,5 мкм, глубина залегания перехода база-коллектор ~4,8 мкм. Максимальная концентрация примеси р-типа в области базы ~8,65х1015 см-3 находится на глубине ~2,84 мкм; максимальная концентрация примеси «-типа в области эмиттера ~1,38х1020 см-3 находится на глубине ~0,64 мкм. Концентрация легирующей примеси «-типа в области эмиттера (Э) значительно выше концентрации легирующей примеси р-типа в области базы (Б). При заданных соотношениях концентраций легирующих примесей «- и р-типов коэффициент усиления транзистора при включении его по схеме с общим эмиттером составляет более 200.

Результаты моделирования электрических характеристик СИТ

Проведены расчеты распределений токов вдоль направленияХ (рис. 3, 4) и вольт-амперной характеристики (рис. 5) транзистора со статической индукцией. Было выполнено моделирование только той части структуры, которая непосредственно работает как транзистор. Начиная примерно с 23 мкм в направлении Х (рис. 1, 2), разностная концентрация примеси |Л'а(х)-^а(х)| в моделируемой структуре не меняется и равна 4,4х1018 см-3. Толщина эпитакси-альной пленки, в которой была получена моделируемая структура СИТ, составляет ~30 мкм; удельное сопротивление эпитаксиальной пленки ру ~6 Омсм. Подложка, на которой была выращена эпитаксиальная пленка, имеет удельное сопротивление ру ~0,01 Ом-см и толщину ~460 мкм. Таким образом, в условиях низкого уровня инжекции можно пренебречь сопротивлением, которое вносит подложка в общее сопротивление СИТ. Следовательно, та часть структуры, которой пренебрегалось при расчете, может рассматриваться просто как резистор с очень низким сопротивлением, включенный последовательно с транзистором.

Рис. 3. Распределение плотностей токов по структуре СИТ вдоль направления Х при напряжении база-коллектор Кбк=3 В и при напряжении база-эмиттер Кбэ=0,1 В

1.40Е+02 1.20Е+02 ^ 1.00Е+02 § 8,00Е+01

■X о и

й 6.00Е+01 н

|| 4,00Е+01 2.00Е+01 1.00Е-05

-•—Ток электронов при Убэ=0,7 В

—■—Ток дырок при Убэ=0,7 В

—Суммарный ток

электронов и дырок при Убэ=0,7 В

10 15 20

Глубина структуры X, мкм

-щ—

25

30

35

Рис. 4. Распределение плотностей токов по глубине структуры СИТ при напряжении база-коллектор Кбк=3 В и при напряжении база-эмиттер Кбэ=0,7 В

В расчетах предполагалось, что моделируемый транзистор включен по схеме с общей базой. Токи, протекающие через транзистор, при малом смещении перехода эмиттер-база (~0-0,5 В) будут очень малы (рис. 5). По мере увеличения напряжения эмиттер-база Уэб потенциал на эмиттере понижается по отношению к потенциалу базы. В результате потенциальный барьер для подвижных носителей заряда на переходе эмиттер-база понижается. Это приводит к росту концентрации дырок в эмиттере п-типа и к росту концентрации электронов в базе р-типа. В целом с ростом Уэб при обратно смещенном переходе база-коллектор (Р"6к=3 В) транзистор постепенно открывается (рис. 3-5), и ток эмиттера растет по экспоненциальному закону (рис. 5).

Зависимость тока эмиттера от напряжения база-эмиттер при Убк=3 В

3.50Е-Ю2

^ 3.00Е-Ю2 ей"

Сч

| 2.50Е-Ю2

(П

И 2.00Е-Ю2 □

н Ё

С 1.50Е-Ю2 о

¡5 □

§ 1.00Е-Ю2

£ 5.00Е-Ю1 £ К т

0.00Е-Ю0

0,0

Рис. 5. Зависимость тока эмиттера СИТ от напряжения база-эмиттер при напряжении

база-коллектор Кбк=3 В

Распределения плотностей токов в направлении Х при различных напряжениях на переходе база-эмиттер и при Кбк=3 В показаны на рис. 3, 4. На рис. 3 изображены распределения электронной и дырочной составляющих, а также суммарной плотности тока при Кбэ=0,1 В. Как видно из рис. 3, суммарный ток является неизменной величиной в направлении, в котором проводился расчет. В области базы (Б) электронная составляющая плотности тока, протекающего через транзистор, принимает минимальные значения, а дырочная составляющая — максимальные значения. Аналогичные распределения плотностей токов получены при разных значениях напряжения база-эмиттер Кбэ=0,1, 0,3, 0,5 В, т.е в условиях низкого уровня инжекции. В условиях высокого уровня инжекции (рис. 4) электронная и дырочная составляющие плотности тока постоянны. Однако значения плотности тока электронов на несколько порядков превосходят значения плотности дырок. Поскольку уровень легирования эмиттера значительно выше уровня легирования базы (максимальное значение концентрации

20 3 15 3

примеси в области эмиттера 1,33x10 см , в области базы — 8,16x10 см ) в случае низкого уровня инжекции (рис. 3) значения плотности тока электронов в областях эмиттера и коллектора значительно превышают значения плотности тока дырок в этих же областях структуры.

Заключение

Представленные результаты демонстрируют возможности управления вольт-амперными характеристиками кремниевого транзистора со статической индукцией в различных режимах функционирования посредством целенаправленного изменения технологических параметров его изготовления.

Автор статьи выражает благодарность проф. кафедры "Микроэлектроника" БГУИР, д-ру физ.-мат. наук В.В. Нелаеву за плодотворные дискуссии и ценные замечания.

SIMULATION OF A SI BASED STATIC INDUCTION TRANSISTOR

N.L. DUDAR Abstract

The semiconductor structures simulation provides a possibility to cut down expenses for manufacturing and optimization the integrated circuit elements and devices parameters. In this paper results of a Si based static induction transistor simulation are presented. Static induction transistors (SITs) are short channel FET structures that are suitable for high power, high temperature and high frequency applications. The static induction transistor to be calculated was formed on the silicon substrate with epitaxial film. P-base and «+-emitter regions were formed on the «-epitaxial film. The routine of manufacturing SIT process involves also growing of the local oxide which provides the good isolation of separated p+- and «+-regions. The one-dimensional SIT model was developed. This model was used for the investigation of SIT electric characteristics. The base of the model is the fundamental system of semiconductor equations that can be solved without the large computing time. The current densities distributions and I-V features were calculated using results of SIT technology simulation (impurities concentrations, potential etc.) as an input data. SIT technology simulation was performed by SUPREM3 program package. The program SUPREM3 creates the output data file involving tables of the simulated structure parameters numerical values. The current densities distributions and I-V features calculation was performed by Gummel numerical method.

Литература

1. Моделирование элементов и технологических процессов / Под ред. Антонетти и др., Пер. с англ. М., 1988. 490 с.

2. Dudar N., Shahov V. // 20th International Scientific Symposium of Students and Young Researchers. 1998. Vol. II. P. 165-169.

3. Gummel H.K. // IEEE Trans. Electron. Dev. 1964. Vol. ED-11, N 10. P. 455-465.

4. Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. М., 1978.

5. Польский Б.С. Численное моделирование полупроводниковых приборов. Рига, 1986.