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

книги / Математические модели элементов интегральной электроники

..pdf
Скачиваний:
21
Добавлен:
12.11.2023
Размер:
14.39 Mб
Скачать

рые с помощью численных методов решаются нелегко. Их применение ограничивается быстродействием и объе­ мом памяти вычислительной машины. Кроме этого, статистические разбросы значений исходных параметров, свойственные любому реальному технологическому про­ цессу, могут свести на нет высокую точность таких мо­ делей. С другой стороны, использование упрощенных ММ также приводит к снижению точности моделирова­ ния.

Рис. 1.1. Классификация математических моделей компонентов.

В зависимости от вида проектируемой схемы (мели-, нейная или линейная, низкочастотная или высокочастот­ ная и т. д.) разработчик должен применять те или иные ММ. При анализе нелинейных схем необходимо исполь­ зовать модели, которые описывают характеристики при­ боров в широком диапазоне значений напряжений и токов; для малосигнальных линейных схем важно повы­ сить точность модели только вблизи рабочей точки; в высокочастотных схемах возрастают роль собственной инерционности приборов и влияние паразитных элемен­ тов. При анализе схем, работающих при изменяющихся внешних условиях (температура, облучение, электромаг­ нитные поля и т. д.) в системе параметров и в ММ при­ боров следует учесть влияние всех этих воздействующих факторов.

Существуют две наиболее общие разновидности мо­ делей: формальные и физические (рис. 1.1).

Формальные модели используются для аппроксима­ ции характеристик приборов в тех случаях, когда физи*

и

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

 

 

 

н

О

 

 

при U < U noр,

 

 

( 1. 1)

 

 

 

 

 

 

 

 

 

 

 

k ( U — С/пор)

При U^>U nop*

 

 

 

 

 

 

 

 

Здесь / — ток диода;

V — приложенное

напряжение;

^пор — пороговое напряжение.

 

 

 

 

 

 

 

Как правило, достаточно точные формальные модели

требуют

большого

количества

измерений

и

аппрокси­

 

 

 

 

мируют

характеристики

приборов

 

 

 

 

только в ограниченной области ра­

 

 

 

 

бочего

диапазона.

Обычно они

ис­

 

 

 

 

пользуются для

 

предварительных

 

 

 

 

оценочных

расчетов.

 

в той

или

 

 

 

 

Физическая

модель

 

 

 

 

лной мере отражает процессы, про­

 

 

 

 

текающие в полупроводниковом при­

 

 

 

 

боре. В отличие от формальной мо­

Рис.

1.2. Аппроксима­

дели уравнения

физической модели

выводятся на основании теории ра­

ция

(-----------)

вольт-

боты прибора. Основным признаком

амперной

характери­

стики диода кусочно­

классификации

моделей

компонен­

линейной

функцией.

тов (с точки зрения разработчика

ра

исходных

 

ИС)

является характеристика набо­

параметров.

Именно этим

определяется

область

возможного использования модели:

электриче­

ское, конструкторско-топологическое

и

технологическое

проектирование. Эта классификация предполагает по­ строение математических моделей трех типов.

Электрическая модель — модель, в которой исходны­ ми являются электрические параметры компонентов (коэффициент усиления транзистора, обратный ток, емкость р—я-перехода и т. п.), определяемые из внеш­ них электрических измерений компонентов. В большин­

стве случаев эти параметры должны

входить и входят

в технические условия на прибор

(для дискретного

исполнения).

 

Физико-топологическая *) модель — модель, в которой исходными параметрами являются геометрические раз­ меры областей компонентов и физические характери­ стики полупроводника (концентрация примесных атомов и ее зависимость от координаты, время жизни неравно­ весных носителей заряда, плотность поверхностных состояний, подвижность носителей и т. п.). Геометриче­ ские размеры определяются рисунком фотошаблона (или топологическим чертежом), а также толщинами слоев диффузионно-эпитаксиального профиля; физи­ ческие параметры определяются из соответствующих электрофизических измерений. В большинстве своем параметры этой группы приводятся в конструкторских чертежах на ИС.

Технологическая модель —модель, в которой исход­ ными являются параметры технологических режимов (температура и время диффузии или эпитаксии, ско­ рость потока газачюсителя, концентрация диффузанта, состав и температура травителя и т. п.). Эти параметры определяются с помощью соответствующей контрольнотехнологической аппаратуры; допуски на их значения содержатся в технологической документации.

В электрической модели искомыми выходными ха­ рактеристиками являются токи и напряжения на выво­ дах компонента. В физико-топологической модели в упрощенном варианте выходные характеристики —это электрические параметры компонентов I(T . е. исходные параметры электрической модели), в общем случае — токи и напряжения на выводах. Технологическая модель может быть трех степеней сложности, обеспечивая воз­ можность расчета либо физико-топологических параме­ тров, либо минуя эту стадию возможность расчета параметров приборов, либо, наконец, определения токов

и напряжений на

выводах

прибора.

Как

правило,

попользуются модели первой и второй

степени слож­

ности.

модели

получаются

из

упрощенной

Электрические

физической теории и учитывают только основные эффек­ ты. Параметры этих ММ определяются из электриче­ ских измерений. Отличие электрической модели от фор-

В литературе иногда используется (в том числе в публика­ циях авторов этой книги) термин «структурная» модель. По-нашему мнению, этот термин юе совсем правильно отражает существо дела, а в ряде случаев может привести к путанице,

мальной состоит в том, что вид уравнений модели, или вид эквивалентной схемы, получается иа основании физической теории прибора, поэтому область действия такой модели значительно шире, чем формальной. Фи­ зические эффекты, которые не были учтены в упрощен­ ной теории, могут быть отражены в модели введением формальных коэффициентов и аппроксимаций. Напри­ мер, из диффузионной теории р—/г-перехода следует, что статическая характеристика диода описывается экспоненциальной функцией вида

Ujmtslf

( 1.2)

/ = / о(е

I - 1),

где фT = kTlq — температурный

потенциал.

Здесь коэф­

фициент т отражает отклонение характеристик реаль­ ного прибора от теоретических.

Выражение (1.2) содержит два параметра /0 н /н, которые определяются из измерений вольт-амперных характеристик реального диода. Измеренные значения /0 и т отражают наличие не только диффузионного тока, но и влияние тока генерации-рекомбинации в рп- переходе, токи утечки и т. д. (т. е. эффекты, которые не учитываются диффузионной теорией).

Электрические модели широко применяются для проектирования схем на дискретных элементах. Естест­ венно, готовые элементы удобно характеризовать систе­ мой электрических параметров (например, для транзи­ стора: нормальным и инверсным коэффициентами усиле­ ния, обратными токами и т. д.). Такие модели также используются для расчета ИС. В этом случае их пара­ метры определяются из измерений на специальных тестовых образцах или рассчитываются с помощью фи­ зико-топологической модели.

Ко второму типу относятся физико-топологические модели. Эти модели должны учитывать все основные эффекты, влияющие на работу прибора, так как они используются для проектирования самих полупроводни­ ковых компонентов, а не для аппроксимации характе­ ристик уже изготовленных приборов. Поэтому, как пра­ вило, физико-топологические модели значительно слож­ нее электрических. Например, физико-топологическая модель диода (одномерная, статическая, для низких

Н

уровней инжекции) имеет вид

/ =

Ь

ЧтА

р=б

lcU’, -a''ir n

I

1&+1)г U

 

V

 

 

JLqnidA

1р иР-п№т _

(1.3а)

 

 

хл + ~Р

 

 

 

 

U =

Up_n -\-IR б,

Яб = РбОУб/А

 

 

 

b =

м«/ р-р.

(1.36)

Здесь цп, тл, [Лр, Тр — подвижность и время жизни элек­ тронов и дырок соответственно; Л — площадь рп-\пере­ хода; ре, Шб — удельное сопротивление и толщина базы; pt —удельное сопротивление полупроводника с собст­ венной концентрацией /г*; d - - ширина рп-перехода; LK—диффузионная длина; £/р_п— напряжение на рп- переходе.

Выражения (1.3а) и (1.36) более точно по сравне­ нию с (1.2) описывают характеристику диода: диффу­ зионный ток (экспонента с показателем <рт), ток гене­ рации-рекомбинации (показатель 2срт)> сопротивление базы (1.36). Основным признаком, характеризующим эту модель как физико-топологическую, является то, что в качестве параметров она содержит геометрические размеры (Л, WQ, d) и электрофизические параметры (М-п, М-р , тп, тР) Lfl, рб, pi). ,При необходимости можно использовать более сложную модель диода, учитываю­ щую влияние профиля распределения примесей в рп- переходе, поверхностную генерацию-рекомбинацию но­ сителей и т. д. Физико-топологические (структурные) ММ используются при разработке дискретных приборов и главным образом при проектировании монолитных ИС, так как этот процесс включает в себя и разработку всех

еекомпонентов.

Впрактике проектирования ИС применяются элек­

трические и физико-топологические модели разной сложности, а также их сочетания. Например, при изго­ товлении высокочастотных биполярных ИС режимы диффузии зафиксированы, в результате чего все типы схем имеют одинаковый диффузионный профиль. Следо­ вательно, часть параметров модели транзистора, зави­ сящих только от диффузионного профиля (удельные токи и емкости р—гг-переходов, коэффициент усиления

it т. д.), целесообразно определять из внешних электро ческих измерений на выводах транзистора.

Несколько особняком от электрических и физикотопологических моделей стоят технологические модели, в которых исходными являются па.рамет,ры технологиче­ ского /процесса. Например, глубины эмиттерного и кол­ лекторного р—«-переходов биполярного транзистора определяются следующими технологическими параме­ трами: количеством диффузанта, временем диффузии, температурой и т. д. Поскольку ход технологических операций определяется большим числом факторов, то для построения технологических моделей должны широ­ ко применяться статистические методы, позволяющие выделять основные факторы. Технологические модели могут быть применены как для анализа, так и для опти­ мизации технологических -режимов изготовления полу­ проводниковых приборов и ИС, а также для расчета исходных параметров физико-топологических моделей.

Разработка моделей этого вида находится пока еще в самой начальной стадии. До настоящего времени не создано ни одной модели, связывающей напряжения и токи на выходах транзистора с технологическими режи­ мами его изготовления. По существу, ведется лишь моделирование тех или иных технологических операций, т. е. построение неких «промежуточных» моделей, в которых искомыми являются структурные характери­ стики (толщины эпитаксиальных слоев, вид диффузион­ ного профиля и т. п.), которые сами в дальнейшем становятся исходными параметрами физико-топологиче­ ской модели.

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

прибора. Кроме того, -различают модели для малого И для большого сигналов. Критерий «малости» сигнала может быть своим специфическим для каждого модели­ руемого объекта; как правило, для полупроводниковых приборов это условие выполняется, если амплитуды воздействующих сигналов не превышают значения срт- В моделях для большого сигнала часть параметров, та­ ких, как коэффициенты усиления по току, емкости р—п-переходов и т. д., зависят от электрического режи­ ма. В малосигнальных ММ эти параметры модели могут

считаться

постоянными, равными своим

значениям

в рабочей

точке, относительно

которой воздействуют

малые по амплитуде управляющие сигналы.

координатам

В зависимости от того, по

скольким

учитываются происходящие в приборе процессы, модели делятся на одно-, двух- и трехмерные. Например, в со­ временных планарных транзисторах физические процес­ сы, происходящие в базовой области, носят ярко выра­ женный двухмерный характер вследствие эффекта вы­ теснения тока эмиттера в периферийные области.

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

Заметим, что вопрос о .степени сложности »и детали­ зации используемой модели зависит от типа решаемой задачи. На практике существует иерархия физических моделей (это относится к электрическим, физико-топо- логическим и технологическим моделям), где каждая последующая модель является более точной и учитывает большее количество эффектов. Появление вычислитель­ ных машин с высоким быстродействием и большим объемом памяти позволяет применять более точные модели, 4TOt в свою очередь, влечет за собой повышение точности и качества проектирования полупроводниковых приборов и ИС.

1.2. Исходные уравнения

Расчет полупроводниковых приборов и компонентов ИС заключается в выборе параметров исходного полу­ проводникового материала с необходимыми электрофи­ зическими свойствами, параметров профиля примеси,

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

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

Физические процессы в любом полупроводниковом приборе могут быть описаны системой уравнений, вклю­ чающей уравнения непрерывности для электронов и ды­ рок и уравнение Пауссона [4]. Постоянные параметры исходного материала еш ед, электрофизические парамет­ ры пи |in, р,р, %пу хР и т. д. и концентрация примеси N.(x, уу г) непосредственно входят в эту систему урав­ нений в качестве параметров. Геометрические размеры и конфигурация диффузионных областей и омических контактов учитываются граничными условиями.

Для общего трехмерного случая система основных уравнений в векторной форме записывается следующим

образом:

 

 

 

 

 

dnjdt = (\/q)bJn + gn — гnr

 

(1.4)

 

in == q (р*л#Е

Dn&n)t

 

(1.5)

 

^ / ^ = -(l/<7)Ajp+ gp- r p,

 

(1.6)

 

)P = qi}>;,pb-Dpb.p),

 

(1.7)

 

j = Ы - i« "h JcM*

 

(1-8)

 

A2<P = ДЕ = (gr/e„eo) (p

n - Na).

(1.9)

Здесь n,

p — концентрация электронов и дырок соответ­

ственно;

ф — электростатический потенциал;

Е — вектор

напряженности электрического

поля; Dn, Dv — коэффи­

циенты диффузии электронов и дырок; jnn,

— подвиж­

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

диэлек­

трическая проницаемость полупроводника; g — скорость генерации электронов и дырок; j, j„, jp — векторы плот­ ности полного тока, тока электронов и тока дырок; jCM= = —епео(д/д/)Дф вектор плотности тока смещения; NR, Na— концентрации доноров и акцепторов.

Уравнения (1.4) и (1.6) являются уравнениями не­ прерывности для электронов и дырок и связывают вели­ чину изменения концентрации носителей в данной точке (dn/dt) с изменением плотности тока [(l/?)Ajn] и ско-

18

ростями 'генерации (gn) н рекомбинации (г„). Уравне­ ния (1.5) и (1.7) являются уравнениями для плотности токов электронов jn и дырок jp. Плотность тока содер­ жит две составляющие: дрейфовую (qn\inE) и диффу­ зионную (qD„&n). Вектор плотности полного тока j (1.8) равен сумме векторов электронного jn и дырочного jp токов и тока смещения jCMУравнение Пуассона (1.9) устанавливает связь между изменением электрического поля и зарядом.

Система уравнений (1.4) — (1.9), являясь достаточно общей, справедлива для макроскопических процессов. Она не учитывает влияния магнитного поля и неодно­ родностей полупроводникового материала (возможные нарушения и дефекты кристаллической решетки и т. д.). Приведенные уравнения справедливы при анализе про­ цессов в полупроводниковых приборах с размерами диффузионных областей не менее сотых долей микрона, градиентами концентраций легирующих примесей не более 1024 см-4 (для Si и Ge), при плотностях тока через

р —л-переходы не более 10б

А/см2

и иа частотах, не

превышающих величину 1012

Гц [17]. Эти ограничения

вряд ли могут быть нарушены

в

реальных приборах,

поэтому на уравнениях (1.4)

(1.9)

базируются все ра­

боты по расчету полупроводниковых приборов и компо­ нентов ИС.

В общем виде граничные условия, необходимые для решения системы уравнений (1.4) —(1.9), формулируют­ ся следующим образом:

1) На свободной поверхности полупроводниковой структуры нормальные составляющие электронного и дырочного токов определяются процессами поверхност­ ной .генерации-рекомбинации. Действительно, если по­ верхность не соприкасается с каким-либо проводником, то полный ток через нее равен нулю. Следовательно, потоки электронов и дырок к поверхности должны быть равны числу рекомбинирующих носителей, т. е. превы­ шению скорости рекомбинации над скоростью гене­ рации:

= (L1°)

где 0 — нормаль к рассматриваемой поверхности.'

Нормальная составляющая электрического поля (при отсутствии у дайной поверхности внешних источников

поля) равна нулю

Е0 = -д<?1д6 = 0.

(1.11)

2) Если полупроводник контактирует с диэлектри­ ком, через который на полупроводник воздействует внешнее поле (на таком принципе основана работа МДП-приборов), то граничное значение поля определя­ ется из условия перехода вектора Е через границу двух сред:

£п£оЕ0п— ел£оЕвД = — Qпоп»

(1-12)

Etn = Eu ,

(1.13)

где QnoD — плотность на единицу поверхности заряда на границе раздела диэлектрик — полупроводник; £*п, Etд — тангенциальные составляющие векторов напря­ женности электрического поля в полупроводнике и ди­ электрике. Величины в выражениях (1.12), (1.13) с ин­

дексом П относятся к

полупроводнику, а с индексом

Д — к диэлектрику.

 

3) Наповерхности полупроводника, контактирующей

с металлом (контакт

металл — полупроводник или

барьер Шоттки), имеет место термодинамическое рав­ новесие и суммарный заряд равен нулю

 

np = nzi,

(1.14)

9 =>q(p-\-NA— n — N J = 0.

(1.15)

4) Ток через выводы прибора определяется интегра­

лом

 

 

 

I = \ j ( x , y,z)dA,

(1.16)

 

А

 

где А — площадь

омического контакта;

dA — элемент

площади; / (дс, у, z)

—плотность тока.

 

При анализе (полупроводниковых приборов необходи­ мо рассматривать не только область полупроводника, но и примыкающие к нему области диэлектрика и металла. В этих случаях приходится совместно с системой урав­ нений (1.4) —(1.9) рассматривать уравнения, описываю­ щие физические процессы в смежных областях, напри­ мер, уравнение Лапласа, которое описывает распреде­ ление потенциала в диэлектрике А2ср=0.

20

Соседние файлы в папке книги