Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

3509

.pdf
Скачиваний:
2
Добавлен:
08.01.2021
Размер:
575.24 Кб
Скачать

Федеральное агентство по образованию Государственное образовательное учреждение высшего профессионального образования Воронежская государственная лесотехническая академия

Н.Н. Панюшкин, А.Н. Панюшкин

ФИЗИЧЕСКИЕ ОСНОВЫ ПРОМЫШЛЕННОЙ ЭЛЕКТРОНИКИ

Лабораторный практикум

ВОРОНЕЖ 2007

Федеральное агентство по образованию Государственное образовательное учреждение высшего профессионального образования Воронежская государственная лесотехническая академия

Н.Н. Панюшкин, А.Н. Панюшкин

ФИЗИЧЕСКИЕ ОСНОВЫ ПРОМЫШЛЕННОЙ ЭЛЕКТРОНИКИ

Лабораторный практикум

ВОРОНЕЖ 2007

УДК 621.38

Панюшкин, Н.Н. Физические основы промышленной электроники [Текст]: лаб. практикум / Н.Н. Панюшкин, А.Н. Панюшкин; Фед. агентство по образованию, ГОУ ВПО “ВГЛТА”. – Воронеж, 2007. – 48с.–ISBN 5-7994-0167-0.

Лабораторный практикум состоит из четырех лабораторных работ, посвященных анализу основных типов полупроводниковых приборов: полупроводниковому диоду на основе р-п перехода, биполярному транзистору, тиристору и МДП-транзистору. Лабораторные работы выполняются методами математического моделирования в системе алгебраических вычислений MathCAD. Они позволяют студентам изучить основные физические процессы в полупроводнике, определяющие характеристики полупроводниковых приборов, провести анализ зависимости электропараметров функциональных ячеек биполярных и МДП интегральных схем от характеристик входящих в них транзисторов, исследовать влияние температуры полупроводника на параметры приборов в диапазоне от –60 0С до +125 0С.

Предназначен для студентов специальности 240400 (190702)- Организация и безопасность движения (Автомобильный транспорт).

Табл. 2. Ил. 21. Библиогр.: 5 наим.

Печатается по решению редакционно-издательского совета ГОУ ВПО “ВГЛТА”

Рецензент: кафедра естественнонаучных дисциплин Российского государственного социального университета Воронежский филиал;

проф. кафедры физики полупроводников и микроэлектроники ВГУ, д -р техн. наук Б.К. Петров

©Панюшкин Н.Н., Панюшкин А.Н., 2007 ©ГОУ ВПО “Воронежская государственная

лесотехническая академия”, 2007

Порядок выполнения лабораторных работ

Правила поведения в компьютерном классе

1.К выполнению лабораторных работ допускаются студенты, получившие допуск у преподавателя. Для получения допуска студент должен подготовить конспект, включающий номер и название лабораторной работы, цель и краткое описание теоретического материала с основными расчетными формулами. При получении допуска студент должен продемонстрировать знания в объёме, необходимом для её выполнения. Получение допуска возможно только после получения зачёта по предыдущей работе.

2.В процессе выполнения работы не разрешается самостоятельно включать/выключать компьютер или дополнительные устройства. Не разрешается также самовольное подключение или отключение от электросети любых потребителей электроэнергии. При необходимости перезапуска компьютера или включения принтера студент должен пригласить преподавателя.

3.При выполнении работы студенту категорически запрещается выходить за пределы рабочей программы и пользоваться принесенными магнитными дисками, CD –дисками и наушниками.

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

Оформление отчета по лабораторной работе

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

Отчет должен содержать разделы:

формулировку цели работы;

краткие теоретические сведения;

результаты исследований на компьютере с выводами по каждому заданию;

заключение (выводы), где указывается

какие исследования проводились в работе;

анализ полученных результатов.

Лабораторная работа № 1

Математическое моделирование вольтамперной характеристики p-n перехода

Цель работы: исследовать вольтамперную характеристику (ВАХ) p-n перехода.

Описание модели, постановка задачи

Электронно-дырочным переходом называется контакт полупроводников различного типа проводимости. После приведения полупроводников в контакт начинается процесс перераспределения носителей заряда (НЗ) под действием градиентов концентраций. Электроны начинают переходить из полупроводника n-типа в полупроводник р-типа, а дырки в противоположном направлении. Движение НЗ создаёт диффузионные токи, направленные навстречу друг другу. Величины этих токов прямо пропорциональны градиентам концентраций НЗ между контактирующими слоями:

jpдиф

= −qe Dp

dp

;

(1.1а)

 

 

 

dx

 

 

dn

(1.1б)

 

 

 

jnдиф = qe Dn dx ,

 

где qe=1,6·10-19 Кл – элементарный заряд; Dn , Dp –коэффициенты диффузии для электронов и дырок соответственно; p, n – концентрации дырок и электронов.

Dn = ϕT µn ,

Dp = ϕT µ p ,

(1.2)

где ϕT -температурный потенциал; µn ,µ p -подвижности для электронов и ды-

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

φ =

kT

,

(1.3)

Т

qe

где k=1,38·10-23 Дж·К-1 – постоянная Больцмана; Т- абсолютная температура. Знак “–” в выражении (1.1а) отражает тот факт, что движение положи-

тельных зарядов происходит в направлении убывания концентрации (производная dp/dx отрицательная).

После ухода электронов из полупроводника n-типа в нём обнажается некомпенсированный положительный заряд ионов донорной примеси. Соответственно в полупроводнике p-типа обнажается отрицательный заряд ионов акцепторной примеси. Внутри p-n перехода возникает электрическое поле, напряженность Е которого в соответствии с уравнением Пуассона будет расти по мере увеличения количества НЗ, перешедших из одной области в другую.

Вектор Е направлен от n-слоя к р- слою, его появление приводит к возникновению дрейфовых токов НЗ, противоположных диффузионным (поле стремится вернуть обратно перешедшие НЗ). Величины диффузионных токов прямо пропорциональны напряженности электрического поля p-n перехода:

jnдр = qe µn n E; jpдр = qe µ p p E . (1.4)

Таким образом, параллельно процессу перехода НЗ под действием градиента концентрации будет постепенно возрастать интенсивность процесса возврата НЗ под действием электрического поля p-n перехода. Полный ток через переход будет равен сумме дрейфовых и диффузионных токов НЗ

j = jnдиф + jpдиф + jnдр + jpдр .

(1.5)

Процесс перераспределения НЗ прекратится, когда дрейфовые токи уравновесят диффузионные. Полный ток перехода при этом будет равен нулю.

Если теперь уменьшить напряжённость поля перехода с помощью внешнего источника, то дрейфовые составляющие тока уменьшатся, и полный ток уже не будет равен нулю. Чем больше будет разность потенциалов, создаваемых внешним источником на p-n переходе, тем меньше будет напряжённость поля внутри перехода, тем меньше будут дрейфовые токи, и соответственно больше будет полный ток перехода. Такое включение p-n перехода называется прямым, оно соответствует подключению положительного электрода внешнего источника к р- области, а отрицательного к n – области.

Если увеличить напряжённость поля перехода путём обратного подключения внешнего источника (положительного электрода к n - области, а отрицательного к р – области), то дрейфовые составляющие тока не могут увеличиться, так как диффузионные токи не изменяются и, следовательно, не может измениться количество возвращаемых НЗ. Незначительное увеличение полного

тока при обратном включении происхо-

I

 

дит за счёт неосновных носителей заря-

 

 

 

 

 

 

 

 

 

да, которые приближаются к p-n перехо-

 

 

ду в процессе теплового движения. По-

 

 

этому ток p-n перехода при обратном

 

 

включении

называют также

тепловым.

 

 

Вольтамперная характеристика p-n пе-

 

 

рехода показана на рис. 1.1.

 

 

 

Для численного определения тока

 

 

через переход необходимо знать про-

0

U

странственно-временное

распределение

концентраций

НЗ,

которое

находится

 

 

Рис. 1.1 Вольтамперная характери-

решением

уравнений

непрерывности.

стика p-n перехода

 

Если пренебречь влиянием

электриче-

 

 

ского поля и процессами ионизации по-

 

 

лупроводника

под

действием внешних

 

 

факторов, уравнения непрерывности упрощаются и называются уравнениями диффузии:

дp

= −

p p

0

+ D

2 p

;

дt

τ

p

 

 

 

x

2

 

 

 

 

 

 

 

(1.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n n

 

 

 

 

2n

дn

= −

 

0

 

+ D

 

 

 

 

 

 

 

 

,

 

дt

 

τ

n

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где τn, τp – время жизни неравновесных носителей заряда – электронов (n) и дырок (p) соответственно; n0, р0- равновесные концентрации электронов и дырок; D –коэффициент амбиполярной диффузии, определяемый по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ p

 

 

 

 

 

 

D = ϕ

 

n

0

 

 

 

 

 

 

 

0

 

 

 

 

.

(1.7)

 

 

n

 

p

 

T

+

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Амбиполярный коэффициент диффузии является усреднённым коэффициентом диффузии при совместном движении электронов и дырок. Более подвижные электроны движутся быстрее дырок. Возникающее вследствие этого пространственное разделение зарядов создаёт электрическое поле, ускоряющее дырки и замедляющее электроны (эффект Дембера). В результате электроны и дырки диффундируют совместно с коэффициентом D. Уравнения (1.6) уже не связаны между собой уравнением Пуассона, и их можно решать отдельно, распространяя результаты решения одного из них на другое.

Плотность тока р-п перехода будет определяться суммой плотностей диффузионных токов электронов и дырок, определяемых по формулам (1.1а) и (1.1б).

ЗАДАЧА. Исследовать пространственно-временное распределение концентраций неравновесных носителей заряда при прямом и обратном смещении p-n перехода. Смоделировать ВАХ идеального p-n перехода.

Численное решение диффузионного уравнения

Для решения (1.5) воспользуемся явным методом Эйлера, согласно которому производные заменяются конечно-разностными уравнениями:

 

 

 

p

 

pi, j pi1, j

;

 

 

 

 

 

 

 

t

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 p

(pi, j+1 pi, j )(pi, j

 

pi, j1 )

=

 

 

pi, j+1 2pi, j

+ pi, j1

;

x2

 

x2

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

ni, j ni1, j

 

;

 

 

 

(1.8)

 

 

 

t

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2n

 

(ni, j+1 ni, j )(ni, j

 

ni, j1 )

=

 

ni, j+1 2ni, j

+ ni, j1

,

 

x2

 

x2

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где t, x – шаг интегрирования по времени и по координате соответственно, рi,j-значения концентрации НЗ при i-ом значении времени и j-ом значении координаты. Подставляя (1.8) в (1.6), выразим значения концентраций НЗ в j-ом узле для i-го момента времени:

 

 

 

 

 

 

t

 

 

t

 

 

pi, j+1

+ pi, j1

2pi, j

 

 

p

i+1, j

= p

1

 

 

+ p0

 

+ D

t

 

 

 

 

 

 

.

(1.9)

 

 

 

 

 

2

 

 

 

 

i, j

 

 

 

 

 

τp

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

τp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

t

 

 

ni, j+1

+ ni, j1 2ni, j

 

 

 

ni+1, j = ni, j 1

 

 

+ n0

τn

+ D t

 

 

 

 

 

 

.

 

(1.10)

 

 

 

 

x2

 

 

 

 

 

 

 

 

τn

 

 

 

 

 

 

 

 

 

 

 

 

В этих выражениях коэффициент амбиполярной диффузии определяется по формуле (1.7).

Для стационарного случая, когда

∂p = 0 и

∂n = 0 , уравнения диффузии

 

 

 

 

 

∂t

 

 

 

∂t

приобретают вид

 

 

 

 

 

 

 

 

 

 

p p0

+ D

2 p

= 0;

 

τp

 

x

2

 

(1.11)

 

 

 

 

 

 

 

 

n n0

 

 

 

2n

 

 

+ D

= 0.

 

τn

x2

 

 

 

 

 

 

 

 

 

 

Для этого случая значения концентраций НЗ в j-ом узле будут определяться по формулам:

p j

=

p j+1 + p j1

+

 

 

 

p0

;

(1.12а)

2

+ del _ p

1+

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

del _ p

 

 

 

 

 

 

 

 

 

 

 

 

 

nj

=

nj+1 + nj1

 

+

 

 

 

n0

 

,

(1.12б)

2

+ del _ n

1+

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

del _ n

где

 

 

x2

 

 

 

del _ p =

 

 

 

 

;

 

 

 

 

 

 

 

 

 

τp D

 

 

 

x2

 

 

 

del _ n =

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

τn D

 

(1.13)

(1.14)

Начальные и граничные значения для уравнений диффузии в общем случае имеют вид

p0, j

p0,N

pi,0

= p0;

 

n0, j

= p0;

 

n0,N

 

U

 

 

 

 

 

ϕ

 

 

= p0e

 

T ;

= n0;

 

 

 

= n0;

 

(1.15)

 

 

U

 

 

 

 

 

n

 

ϕ

 

 

= n0e

 

T .

i,0

 

 

 

 

Здесь U-напряжение на p-n переходе, для прямого смещения U>0, для обратного - U<0, L-глубина слоя полупроводника, x=0 соответствует границе p-n перехода.

На рис. 1.2 показано распределение НЗ по координате для нескольких значений времени в полупроводнике n-типа для прямого смещения на p-n переходе. Аналогичная зависимость для обратного смещения на p-n переходе показана на рис. 1.3. На этих рисунках (1.2 и 1.3) значение x=0 соответствует границе p-n перехода, а x=L глубине n- области.

Зная распределение НЗ, можно найти плотности диффузионных токов, используя выражения (1.1). Конечно-разностные аппроксимации этих выражений имеют вид

n(x,t),

p(x,t)

n(x,t1)

n(x,t2)

 

 

 

 

n(x,t3)

n0

 

 

 

 

p(x,t1)

p(x,t2)

p(x,t3)

p0

0

L x

Рис. 1.2 Распределение концентраций НЗ в полупроводнике n–типа при прямом смещении на p-n переходе (t3>t2>t1)

n(x,t),

p(x,t)

n0

 

 

 

n(x,t1)

n(x,t2)

 

n(x,t3)

p0

 

 

 

0

 

L

x

p(x,t1) p(x,t2)

p(x,t3)

 

 

Рис. 1.3 Распределение концентраций НЗ в полупроводнике n–типа при обратном смещении на p-n переходе (t3>t2>t1)

jn дифi, j = qe Dx (ni, j+1 ni, j );

(1.16)

jp дифi, j = −qe Dx (pi, j+1 pi, j ) .

Величина диффузионного тока прямо пропорциональна производной от концентраций НЗ по координате (согласно выражениям 1.1). Как видно из рис. 1.2, производная от концентраций НЗ максимальна (по модулю) на границе p-n перехода и стремится к нулю при приближении x к L. Поэтому величина диффузионного тока также будет максимальна при x=0 и равна 0 при x=L. В соответствии с выражением (1.4) плотность полного тока через p-n переход равна сумме дрейфовых и диффузионных составляющих. В теории полупроводниковых приборов доказывается, что плотность полного тока не зависит от координаты (что соответствует закону сохранения заряда) и равна сумме максимальных значений диффузионных токов неосновных носителей. Максималь-

ное значение плотности диффузионного тока неосновных носителей будет иметь место при x=0. Для n– слоя

jn (t) = jp диф (t)|x=0 ;

(1.17)

для p– слоя

 

jp (t) = jnдиф (t)|x=0 .

(1.18)

Здесь jn (t) и jp (t) плотности токов в n и p слоях соответственно. Конеч- но-разностные аппроксимации для них будут иметь вид

jn

 

= qe

D

(n

 

n

 

);

 

(1.19)

i

 

 

i,1

i,0

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jp

 

= −qe

D

(p

 

p

 

).

(1.20)

i

 

i,1

i,0

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Ii

= S(jni

+ jpi

)= −Sqe

D

(ni,1 ni,0 )+

D

(pi,1

pi,0

) .

(1.21)

 

 

 

 

 

 

x

x

 

 

 

Составление программы, проведение численных экспериментов

1. Из таблицы с заданиями введите исходные данные для расчёта. Задайте

L=5 10-4м; Um=0,7 В; Т=300 К; k=1,38 10-23Дж /К; qe=1,6 10-19 Кл. При вводе определений и формул учитывайте, что MathСAD выполняет операции в последовательности сверху вниз и слева направо (первыми выполняются операции, расположенные на экране выше и левее). Учитывайте, что MathСAD воспринимает только символы английского алфавита, знаком, отделяющим целую часть от дробной вместо запятой является точка, размерность не указывается. Оператором присваивания вместо = является :=. Например, для задания L=5 10-4м необходимо ввести L, затем нажать одновременно клавиши Shift и :, после этого ввести числовое значение 5 10-4. Знак умножения вводится одновременным нажатием клавиш Shift и *. Для ввода показателя степени нажмите одновременно Shift и ^ . На экране появится L:=5 10-4.

2. По формуле (1.7) рассчитайте коэффициент амбиполярной диффузии D, а по формуле (1.3) величину температурного потенциала ϕТ. Необходимо помнить, что в системе MathСAD нижний индекс является признаком массива, поэтому вместо ϕТ необходимо вводить ϕТ.

3. Задайте время моделирования ТМ=10-6 с, количество шагов (узлов) по времени М=25 и по координате N=25. Рассчитайте значения шагов по времени и

по координате t = TM ; x = L . Вместо символа можно использовать del.

MN

4.Задайте интервалы изменения индексов узлов моделирования по времени и по координате: i=1,..,M-1 и j=1..N-1. Для этого наберите i (имя переменной), одновременно Shift и : , 1 (нижний предел), М-1 (верхний предел). На экране появится i:=1..M-1 . Это означает, что величина i будет принимать целые зна-

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]