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

Кругляк_Физика и моделирование нанотранзисторов_2018_314 стр_обложки

.pdf
Скачиваний:
37
Добавлен:
24.02.2019
Размер:
17.69 Mб
Скачать
TLIN , а (19) – через app . Уравнивая
vsat
vinj

Таким образом, было показано, что ток насыщения (9) в модели прохождения может быть записан в традиционной форме (14) через скорость насыщения, принятой в модели VS, если традиционную диффузионно контролируемую скорость заменить на скорость впрыскивания (26).

Предельным значением скорости впрыскивания является баллистическая скорость впрыскивания vT .

Вписывая экспериментальные данные в модель виртуального истока, мы находим кажущуюся подвижность app и скорость впрыскивания vinj . Если экспериментальные данные относятся к длинноканальному транзистору, то определим и стандартную диффузионную подвижность . Позже мы покажем, что традиционную подвижность можно определить и в случае коротких каналов в нанотранзисторах. Предполагая, что нам известны эти три параметра, покажем как вычислить коэффициент прохождения в линейной области и в области насыщения.

Уравнение (9) дает линейный ток через эти два выражения, получаем

TLIN

app

v 1

 

app

,

L

2kT / q

B

 

 

 

T

 

 

 

а используя далее определение кажущейся подвижности по (20), окончательно находим, что

TLIN

 

B

 

 

 

1

 

 

 

 

 

.

(27)

B

B

B

 

 

 

 

 

 

Для вычисления TSAT воспользуемся выражением (2)

 

v

 

v

 

 

 

TSAT

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 T

 

 

 

inj

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

SAT

 

 

 

 

 

откуда для TSAT имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TSAT

 

 

 

 

2

 

 

 

.

 

 

(28)

1 vT /vinj

 

 

 

 

 

 

 

 

 

 

 

Скорость впрыскивания vinj определяется путем вписывания опытных данных в

модель VS, однако, баллистическую скорость

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

труднее. Ее можно извлечь из анализа

экспериментальных входных

250

 

VDS

характеристик [5], вместе с тем ее часто вычисляют пользуясь известными эффективными массами и гипотезой о числе занятых подзон.

Оценим B , app и TLIN для линейной области n-MOSFET с длиной канала L 22 нм при Т 300 К . Пусть в подложке Si (100) заселена только одна нижняя

подзона, а подвижность 250 см2/В с .

 

 

 

 

В этих условиях, согласно

(15/гл.6),

v 1.2 107 см / с . Баллистическую

 

 

 

 

 

T

 

 

 

подвижность находим по (18)

 

 

 

 

 

B

v L

 

(1.2 107 ) (22 10 7 )

2

/В с .

(29)

T

 

2 0.026

 

508 см

 

2kT/q

 

 

 

 

 

 

Поскольку B сопоставима

с ,

этот

транзистор

работает в

квазибаллистическом режиме.

 

 

 

 

 

Кажущуюся подвижность находим из (20), а именно:

app

B

 

250 508

191

2

 

 

 

см

/В

B

250 508

 

 

 

 

 

Как и ожидалось, кажущаяся подвижность меньше подвижностей – баллистической и диффузионной. коэффициент прохождения в линейной области по (27):

TLIN

 

 

250

0.33 .

B

508 250

 

 

 

с .

наименьшей из двух И, наконец, находим

(30)

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

251

eff

не носит фундаментальный характер, оно возникает только из-за трудностей

вычисления зависимости T (VDS ) .

 

 

 

Мы

до сих

пор встречались

с тремя

разными

подвижностями:

1) традиционная

диффузионная подвижность

, 2)

баллистическая

подвижность B

и 3) кажущаяся подвижность app . В теории MOSFET нередко

используется т.

наз. эффективная подвижность eff [6]. Термин «эффективная

подвижность»

– неудачный, он подразумевает то обстоятельство, что

электроны

вблизи

границы раздела

Si/SiO2

должны

иметь меньшую

подвижность по сравнению с электронами в глубине канала проводимости, в результате дополнительного рассеяния на шероховатостях поверхности раздела. Эффективная подвижность – это подвижность электронов в канале

проводимости, усредненная по глубине. Для Si MOSFET eff

намного меньше

диффузионной подвижности в массивной кремниевой

подложке. Для

транзисторов III-V HEMT высокая подвижность сохраняется благодаря технологической возможности создания атомарно регулярной границы раздела. В современных нанотранзисторах, однако, квантовые ограничения столь сильны, что все электроны в канале проводимости испытывают рассеяние на шероховатостях поверхности раздела. Говорить об усреднении подвижности по глубине применительно к нанотранзисторам не приходится. Обе подвижности стоит уравнять, , и рассматривать их как диффузионную подвижность.

8.3.Транспорт в нанотранзисторах с позиции модели MVS/прохождения

Офизике процессов в нанотранзисторах можно многое выяснить из анализа вольт-амперных характеристик наноразмерных MOSFET. Разработано немало подходов для такого анализа [7 – 12]. Перспективность объединенной модели VS/прохождения хорошо зарекомендовала себя при изучении транспортных явлений в нанотранзисторах [13 – 16]. Продемонстрируем это [17, 18] путем анализа экспериментальных результатов на примере как ETSOI

MOSFET [19, 20], так и III-V HEMT [16, 21].

Обзор модели MVS/прохождения

Прежде чем переходить к анализу экспериментальных данных приведем основные положения последней версии этой объединенной модели, получившей название модели MVS/прохождения, поскольку она была сформулирована в Массачусетском технологическом институте [22].

252

Q(VGi ,VDi ) вычисляется из

Модель MVS описывает ток стока в виде произведения электронного заряда на скорость электронов [22 – 24], а именно:

IDS W |Q(VGi ,VDi )

 

| FSAT (VDi )vinj ,

(31)

 

 

 

 

x 0

 

 

 

 

 

где FSAT (VDi ) vinj – это скорость электронов на виртуальном истоке. Напряжения VGi и VDi – это собственные (intrinsic) напряжения на затворе и на стоке. Электронный заряд взят по абсолютной величине, поскольку инверсионный заряд отрицательный в n-MOSFET.

В модели MVS заряд на виртуальном истоке полуэмпирического выражения, подобного (45/гл. 4) [22], а именно:

 

 

 

 

inv kT

 

q VGi

VT kT / q Ff

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|Q(V ,V )

 

| mC

 

 

 

ln 1

exp

 

 

 

 

. (32)

 

 

 

 

 

Gi Di

 

G

 

 

 

 

 

 

 

x 0

 

 

 

q

 

mkT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В этом выражении используется инверсионная переходная функция (inversion transition function) [22]

Ff

 

 

 

 

1

 

 

,

(33)

 

 

 

 

 

 

V

V kT / q / 2

 

 

 

 

 

 

 

Gi

 

T

 

 

 

 

 

1 exp

 

 

 

 

 

 

 

 

 

 

 

kT / q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которая ответственна за увеличение порогового напряжения на величинуkT /q / 2 по мере того как MOSFET переходит из подпорогового режима в

режим сильной инверсии: Ff 1 в подпороговом режиме и Ff 0 в режиме

сильной инверсии. Эмпирический параметр обычно выбирают равным 3.5 [17, 22]. В уравнении (32) пороговое напряжение зависит от напряжения на стоке согласно (47/гл. 4), а именно:

 

VT VT 0 VDi ,

(34)

где VT 0

– это пороговое напряжение в режиме сильной инверсии при VDD VDi 0 ,

а – это параметр DIBL. Подпороговый параметр наклона в (32)

 

 

m m0 m VDi ,

(35)

где m0

– это подпороговый параметр при VDD VDi 0 , а m dm /dVDi

описывает

изменение параметра m с изменением напряжения на стоке.

 

 

253

 

В модели MVS используется эмпирическая функция насыщения тока стока

(45/гл. 4) [22]

 

FSAT (VDi )

 

 

VDi /VDSATs

 

 

 

(36)

 

 

 

 

 

 

 

 

 

 

 

 

1 (V /V

 

) 1/

 

 

 

 

 

 

 

Di

DSATs

 

 

 

 

с (56/гл. 4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

vinj Leff

,

 

 

 

 

(37)

 

 

 

app

 

 

 

 

 

 

DSATs

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Leff – это

эффективная

длина

канала

проводимости

по [25].

Обращаем

внимание, что в индекс SAT

в VDSATs

добавлена буква s , чтобы учесть тот факт,

что FSAT описывает насыщение тока

стока в режиме сильной инверсии. В

условиях подпорогового режима VDSAT

kT / q [25]. Модель MVS рассматривает

этот переход

от VDSAT в подпороговом режиме

в режим

сильной

инверсии

эвристически с использованием инверсионной переходной функции (33) [22], а именно:

 

VDSAT VDSATs (1 Ff ) (kT/q)Ff .

(38)

Собственные

терминальные

напряжения VDi , VGi связанны

с внешними

терминальными напряжениями VDS , VGS следующими соотношениями:

 

VDi VDS IDS RSD0,

(39)

 

VGi VGS IDS RSD0/2,

 

 

где суммарное

терминальное

сопротивление RSD0 RS 0 RD0

есть сумма

сопротивления истокового терминала RS 0 и сопротивления стокового терминала RD0 , которые предполагаются одинаковыми и независящими от напряжений на затворе и на стоке.

Модель MVS может быть вписана в экспериментально измеренные

передаточные характеристики

IDS VGS

и

выходные характеристики IDS VDS ,

что позволяет определить несколько важных параметров MOSFET.

Дальнейший анализ сосредоточим на линейной области при низких VDS и на

области насыщения при высоких VDS .

 

 

 

 

 

 

 

 

 

При небольших напряжениях на стоке

 

 

 

 

 

 

F

VDS

 

, v

(V

 

,V

 

)

 

 

appVDS

.

 

 

 

 

 

 

 

 

 

SAT

 

 

 

x

GS

DS

 

 

x 0

 

 

 

VDSATs

 

 

 

 

 

 

Leff

 

 

 

 

 

 

 

 

 

 

 

 

 

254

 

 

 

 

 

 

Уравнение для тока (31) в линейной области становится таким:

I

DLIN

 

W

|Q(V

) | V

 

VDS

,

(40)

 

 

 

 

Leff

GS

app DS

 

Rch

 

 

 

 

 

 

 

 

где Rch – сопротивление канала проводимости.

Для больших напряжений на стоке FSAT 1 и уравнение (31) сводится к традиционному выражению для тока насыщения

IDSAT W |Q(VGS ,VDS ) | vT ,

(41)

где

 

 

 

 

 

 

 

 

v

2kT

 

vball

(42)

m *

T

inj

 

 

 

 

есть баллистическая скорость впрыскивания в случае использования статистики Максвелла – Больцмана. Обратим внимание на то, что баллистическую скорость впрыскивания на практике вычислить трудно. Физические напряжения и квантовые ограничения в нанотранзисторах могут сказаться на величине m *, к тому же уравнение (42) предполагает, что только одна нижняя подзона заселена, что не всегда верно.

Кажущаяся подвижность (20) в

модели

MVS зависит от Leff

и дается

выражением

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

1

 

,

(43)

 

 

(L

 

 

)

 

 

 

B

(L )

 

 

 

 

 

 

 

 

 

 

 

app

 

 

 

eff

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eff

 

 

где диффузионная подвижность (17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

vT 0

,

 

(44)

 

kT/q

 

2kT/q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а баллистическая подвижность (18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(L

 

)

vT Leff

 

.

 

(45)

 

B

 

 

 

 

 

 

 

 

 

 

 

 

eff

 

 

 

 

 

 

2kT/q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Скорость впрыскивания при высоком напряжении на стоке определяется

из (26), а именно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

1

 

,

 

 

 

(46)

 

 

 

v

 

 

 

v

 

D/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

inj

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

255

где

Leff и

D

vT 0

.

(47)

 

2

 

 

Мы предположили, что в линейной области средняя длина свободного пробегаLIN совпадает со средней длиной свободного пробега в области насыщенияSAT . Хотя равенство LIN SAT 0 , строго говоря, не совсем корректно, однако, физически оно вполне приемлемо [26] и подтверждается в экспериментах [17]. Наконец, полезно вспомнить как параметры в модели MVS связаны с коэффициентом прохождения. Из уравнения (27) для линейной области имеем

TLIN

 

 

 

 

app

 

 

 

,

(48)

 

 

0

 

 

 

 

 

 

 

 

 

0 Leff

 

 

 

B

 

 

B

 

 

а из (28) для области насыщения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TSAT

 

0

 

 

 

 

 

2

 

 

.

 

(49)

0

 

 

1

vT

 

 

 

 

 

 

 

 

/vinj

 

 

Измеряемая скорость впрыскивания связана с коэффициентом прохождения как

 

v

 

TSAT

 

 

v

 

 

.

(50)

2 T

inj

T

 

 

 

 

 

SAT

 

 

Мы привели основные

результаты

новейшей

версии модели

MVS/прохождения, в основе которой лежит объединенная модель прохождения и модель виртуального истока. Как только экспериментально измеренные вольт-амперные характеристики вписаны в модель MVS, результатами

измерений будем считать

суммарное терминальное сопротивление RSD0 ,

кажущуюся подвижность app

и скорость впрыскивания vinj . Мы также увидим,

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

256

TSOI

8.4. ETSOI MOSFET и III-V HEMT

Проанализируем Si MOSFET, имеющий простую и тщательно детализированную физическую структуру, что весьма благоприятствует поставленной задаче. Как показано на рис. 3, это Si SOI (silicon-on-insulator) структура c исключительно тонким слоем SOI толщиной 6.1 0.4 нм [17].

Рис. 3. Сечение ETSOI MOSFET и зависимость nS VGS , полученная

вписыванием экспериментальных данных (кружки) в модель VS (линии). Показаны SD область расширения (SD extension region/SDE), расширенный SD (raised SD/RSD), слой углубленного оксида (buried oxide/BOX) толщиной

145 нм и длина канала проводимости Leff [17].

Плоскость канала – (100) , направление транспорта электронов – 110 . Электрод затвора – poly Si и оксид затвора – SiON с толщиной эквивалентной

емкости CET 1.1 нм . Емкость затвора в режиме сильной инверсии (21/гл. 3) Сinv

 

 

 

 

 

G

определена из вольт-емкостных (C-V) измерений на длинноканальном

транзисторе

[17].

Для

рассматриваемого

конкретного

n-FET

Сinv 1.98 Ф см 2

[17].

Измеренная квазиравновесная подвижность

для

G

 

 

 

 

 

длинноканального транзистора оказалась равной 350 см2/В с , что соответствует средней длине свободного пробега 15.8 нм .

В этом транзисторе использовались прокладки для снятия механического напряжения (neutral stress liners), так что кремниевый канал номинально не

напряжен, что упрощает вычисление v . Предполагая, что m* 0.22 m , находим

T

 

 

0

vT 1.14 см / с . Технологический процесс

создания SD

областей

расширения

привел к очень резким p-n-переходам

с низким

значением

суммарного

257

 

 

 

терминального сопротивления [19]. Физическая длина электрода затвора определялась путем C-V-измерений [18]. Детальное моделирование процесса изготовления транзистора показало, что имеется перекрывание в 1–2 нм между электродом затвора и S - и D -расширениями для n-MOSFET и p-MOSFET, соответственно, так что Leff LG 2 нм для n-FET и Leff LG 4 нм для p-FET, где LG это физическая длина электрода затвора. Перечисленные длины были подтверждены путем тщательного анализа 2D электростатики [19, 20].

Транзистор с высокой подвижностью электронов HEMT это полевой транзистор, в котором для создания проводимости канала используется контакт двух полупроводниковых материалов III-V групп с различной шириной запрещенной зоны (вместо легированной области как у обычных MOSFET): широкая запрещенная зона работает как изолятор, а узкая запрещенная зона работает как канал проводимости. Моделировался III-V HEMT (рис. 4), подробно описанный в [16, 21]. Транзистор собран на подложке InP. Буфер наращивался на подложке, далее шли слои 2 нм In0.53Ga0.47As, 5 нм InAs и 3 нм In0.53Ga0.47As. Решетки слоя 2 нм In0.53Ga0.47As и подложки InP совместимы, однако, имеет место несоответствие между решетками In0.53Ga0.47As и InAs, так что слой InAs – псевдоморфный, напряжен, но он все же достаточно тонкий и не создает дефекты. Поверх этого 10 нм канала создан барьерный слой In0.52Al0.48As, играющий роль изолятора. Т-образная структура затвора понижает сопротивление затвора, что важно для радиочастотных приложений. Сильно допированные подконтактные SD слои (показаны на рисунке белым цветом) обеспечивают низкие терминальные сопротивления.

Рис. 4.III-V HEMT [21],взятый для анализа в рамках модели MVS/прохождения.

258

 

Измеренная подвижность длинноканального HEMT оказалась равной

12500 см2/В с ,

что

дает

среднюю

длину свободного

пробега 153 нм [18].

Эффективная

масса

m* 0.022 m ,

откуда v 3.62 107 см / с [18].

Слой

 

 

 

 

0

T

 

 

4нм

In0.52Al0.48As

поверх канала

проводимости

показывает

емкость

Сinv

1.08 Ф см 2 [18].

 

 

 

 

G

 

 

 

 

 

 

 

Вписывание модели MVS/прохождения в экспериментальные данные

Анализ экспериментальных характеристик обычно проводится как для выходных, так и для передаточных характеристик. Предполагается, что физические и эффективные длины затвора измерены независимо друг от друга, так же как и емкость затвора в режиме сильной инверсии. Параметр , который контролирует переход от слабой до сильной инверсии, взят равным 3.5 [17, 22]. Параметр в FSAT (36) подобран таким, чтобы вписаться в ток насыщения. Обычно он оказывается в узкой области 1.6 2.0 [22]. Чтобы вписаться в экспериментальные данные, показанные на рис. 5, в качестве подгоночных использовались четыре параметра.

Рис. 5. Экспериментальные характеристики ETSOI MOSFET с Leff 30 нм

(точки) и вписанные в них уравнения модели MVS/прохождения (линии). Данные приведены для восьми значений VGS 0.2(0.1) 0.5 В [27].

259