Кругляк_Физика и моделирование нанотранзисторов_2018_314 стр_обложки
.pdfТаким образом, было показано, что ток насыщения (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 |
|
характеристик [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
не носит фундаментальный характер, оно возникает только из-за трудностей
вычисления зависимости 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
Модель 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
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