Микроэлектроника, 2021, T. 50, № 2, стр. 127-134

Влияние диэлектриков затворов полевых графеновых транзисторов на вольт-амперные характеристики

И. И. Абрамов a*, Н. В. Коломейцева a, В. А. Лабунов a, И. А. Романова a, И. Ю. Щербакова a

a Учреждение образования “Белорусский государственный университет информатики и радиоэлектроники”
220013 Минск, ул. П. Бровки, 6, Республика Беларусь

* E-mail: nanodev@bsuir.edu.by

Поступила в редакцию 05.08.2020
После доработки 08.08.2020
Принята к публикации 20.08.2020

Полный текст (PDF)

Аннотация

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

1. ВВЕДЕНИЕ

Значительное внимание в настоящее время в области микро- и наноэлектроники уделяется разработке приборных структур на основе 2D материалов, в частности, на графене [1, 2]. Графен обладает рядом уникальных свойств для создания новой элементной базы микро- и наноэлектроники, а именно: высокой подвижностью носителей заряда, амбиполярностью, высокой проводимостью. Так, полевые графеновые транзисторы (ПГТ) могут найти применение в высокоскоростной электронике, особенно в аналоговых и радиочастотных устройствах [3].

Существенное различие физики функционирования полевых транзисторов на основе полупроводниковых материалов и полевых графеновых транзисторов обусловило необходимость создания моделей именно для ПГТ. Принципиальное отличие наглядно проявляется на передаточных характеристиках ПГТ, имеющих области дырочной проводимости слева от точки Дирака и электронной проводимости справа от точки Дирака. В последние годы моделям графеновых полевых транзисторов, а также получению экспериментальных образцов этих приборов посвящено большое число работ, что, в частности, отражено в обзоре по компактным моделям ПГТ [4].

Такие модели могут разрабатываться на основе: диффузионно-дрейфовой модели с различного рода модификациями [57]; кинетического уравнения Больцмана с модификациями [8, 9]; функций Грина [10]; формулы Ландауэра [11]. Модели также могут предназначаться для однослойных [7, 11], двухслойных [12, 13] и многослойных [14] ПГТ. В работе [15] описана модель для предложенной конструкции туннельного транзистора на графене.

В данной работе приведено краткое описание предложенной упрощенной комбинированной самосогласованной модели полевого транзистора на основе однослойного графена (рис. 1), а также результаты моделирования, полученные с ее использованием. Исследовано влияние диэлектриков затворов ПГТ на их вольт-амперные характеристики (ВАХ).

Рис. 1.

Поперечное сечение двухзатворного ПГТ на основе однослойного графена: S – исток; G – верхний затвор; BG – нижний затвор; D – сток; 1 – подложка; 2 – диэлектрик нижнего затвора; 3 – диэлектрик верхнего затвора; 4 – металл верхнего затвора; 5 – металл нижнего затвора; 6 – металл истока и стока; 7 – графен.

2. МОДЕЛЬ

Предложенная модель построена на основе уравнений квантовой диффузионно-дрейфовой модели [16, 17] и относится к классу комбинированных согласно классификации [17, 18]. В модели электростатический потенциал канала и квантовая емкость рассчитываются самосогласованно с помощью итерационного метода. При реализации алгоритма передаточных характеристик использован метод дихотомии для учета сопротивлений стока Rd и истока Rs. Детальное описание модели приведено в статьях [1921]. Поэтому здесь приведем основные соотношения модели, являющейся комбинацией физико-топологической и электрической моделей.

На первом этапе в модели вычисляются емкости верхнего ${{С}_{{gt}}}$ и нижнего ${{С}_{{gb}}}$ затворов. Далее самосогласованно рассчитывается электростатический потенциал $V$ с применением итерационного метода по формулам [22, 23]:

(1)
$V = \frac{{{{С}_{{gt}}}\left( {{{V}_{{tg}}} - {{V}_{{tg0}}} - {{V}_{a}}} \right) + {{С}_{{gb}}}\left( {{{V}_{{beff}}} - {{V}_{a}}} \right)}}{{{{C}_{{gt}}} + {{C}_{{gb}}} + \alpha {{С}_{q}}}},$
(2)
$\alpha = \frac{1}{2}\left( {\frac{1}{{1 + \frac{{{{q}^{2}}{{V}^{2}}}}{{{{{(kT\ln (4))}}^{2}}}}}} + 1} \right),$
(3)
${{С}_{q}} = \frac{{2{{q}^{3}}}}{{\pi {{{(\hbar {{V}_{F}})}}^{2}}}}\left| V \right|,$
где ${{С}_{{gt}}}$, ${{С}_{{gb}}}$ – емкости верхнего, нижнего затворов, $\alpha $ – коэффициент емкости ($\alpha = 0.5$ при $q\left| V \right| \gg kT$), $q$ – заряд электрона, $\hbar $ – постоянная Планка, деленная на $2\pi $, ${{V}_{F}}$ – скорость Ферми, равная ${{V}_{F}}$ = 106 м/c, ${{V}_{{tg}}}$ – напряжение на верхнем затворе, ${{V}_{{beff}}}$ – эффективное напряжение на нижнем затворе, ${{V}_{{tg0}}}$ – напряжение Дирака для верхнего затвора, $V$ – электростатический потенциал канала, ${{V}_{a}}$ изменяется в диапазоне от 0 до ${{V}_{{ch}}}$, ${{V}_{{ch}}}$ – полное падение напряжения в канале, $L$ – длина канала, $k$ – постоянная Больцмана, $T$ – температура окружающей среды.

Система уравнений (1)–(3) решается итерационно до тех пор, пока поправка электростатического потенциала в канале не достигнет заданного значения.

На втором этапе на основе найденного электростатического потенциала в канале рассчитывается общая плотность заряда листа графена $Q$ согласно [24]:

(4)
$Q = \frac{{q\pi {{{\left( {kT} \right)}}^{2}}}}{{3{{{(\hbar {{V}_{F}})}}^{2}}}} + \frac{{{{q}^{3}}{{V}^{2}} + q{{\Delta }^{2}}}}{{\pi {{{(\hbar {{V}_{F}})}}^{2}}}}.$

Далее вычисляется эффективная подвижность носителей заряда в графене $\mu $, учитывающая подвижности электронов и дырок, а именно [23]:

(5)
$\mu = \left( {\frac{{14({{\mu }_{p}} - {{\mu }_{n}})\left( {\frac{{\pi {{{\left( {kT} \right)}}^{2}}}}{{3{{{(\hbar {{V}_{F}})}}^{2}}}} + \frac{{{{q}^{2}}}}{{\pi {{{(\hbar {{V}_{F}})}}^{2}}}}{{V}^{2}}} \right)V}}{{\left( {\frac{{\pi {{{\left( {kT} \right)}}^{2}}}}{{3{{{(\hbar {{V}_{F}})}}^{2}}}} + \frac{{{{q}^{2}}{{V}^{2}} + {{\Delta }^{2}}}}{{\pi {{{(\hbar {{V}_{F}})}}^{2}}}}} \right)\sqrt {1 + \frac{{{{q}^{2}}}}{{{{{(kT\ln (4))}}^{2}}}}{{V}^{2}}} }} + h} \right)\frac{m}{{m + {{V}^{2}}}},$
где ${{\mu }_{p}}$ и ${{\mu }_{n}}$ – подвижности электронов и дырок соответственно, $\Delta $ – неоднородность электростатического потенциала, $h$ и $m$ – согласующие параметры модели.

Ток стока определяется по формуле, полученной из квантовой диффузионно-дрейфовой модели, а именно [25, 26]:

(6)
${{I}_{d}} = \frac{{W\int\limits_0^{{{V}_{{ch}}}} {\mu QdV} }}{{L + \left| {\int\limits_0^{{{V}_{{ch}}}} {\frac{\mu }{{{{V}_{{sat}}}}}dV} } \right|}},$
где $W$ – ширина канала, ${{V}_{{sat}}}$ – скорость насыщения.

Скорость насыщения вычисляется на основе модели [24] и позволяет учитывать рассеяние на оптических фононах:

(7)
${{V}_{{sat}}} = \left\{ \begin{gathered} \frac{{2{{V}_{F}}}}{\pi },\,\,{\text{если}}\,\,\left| {{{Q}_{{net}}}} \right| \leqslant q\left| {{{\rho }_{{crit}}}} \right| \hfill \\ \frac{{2q\Omega }}{{{{\pi }^{2}}\hbar {{V}_{F}}\left| {{{Q}_{{net}}}} \right|}}\sqrt {\frac{{\pi {{{(\hbar {{V}_{F}})}}^{2}}\left| {{{Q}_{{net}}}} \right|}}{q} - {{{\left( {\frac{{\hbar \Omega }}{2}} \right)}}^{2}}} ,\,\,\,\,{\text{если}}\,\,\,\left| {{{Q}_{{net}}}} \right| > q\left| {{{\rho }_{{crit}}}} \right|, \hfill \\ \end{gathered} \right.$
${{\rho }_{{crit}}} = \frac{1}{{2\pi }}\left( {\frac{\Omega }{{{{V}_{F}}}}} \right),$
где ${{\rho }_{{crit}}}$ – критическая плотность носителей заряда, ${{Q}_{{net}}}$ – плотность заряда носителей в канале (${{Q}_{{net}}} = \left| {\alpha {{С}_{q}}V} \right|$); $\hbar \Omega $ – эффективная энергия оптических фононов.

Отметим, что, если не учитывать квантовые явления и падения напряжений на областях стока и истока, то в этом случае ${{V}_{{ch}}} = {{V}_{d}}$ (${{V}_{d}}$ – напряжение, прикладываемое к стоку относительно истока), а соотношение (6) приводит к традиционной диффузионно-дрейфовой модели, используемой многими авторами (иногда с модификациями и коррекциями [57, 22, 27]) при расчете тока стока ПГТ.

На заключительном этапе получения передаточных ВАХ ПГТ для учета падений напряжений на сопротивлениях областей стока ${{R}_{d}}$ и истока ${{R}_{s}}$, необходимо перерасчитать ток стока ${{I}_{d}}$ при фиксированном напряжении стока ${{V}_{d}}$, применяя оптимизационный метод. Нами использовался метод дихотомии [1921]. Этот метод с успехом применялся для решения гораздо более сложных задач многомерного численного моделирования фрагментов кремниевых СБИС [28, 29]. Только в частном случае, когда не учитываются сопротивления областей стока ${{R}_{d}}$ и истока ${{R}_{s}}$, т.е. ${{R}_{d}} = {{R}_{s}} = 0$, перерасчет тока не требуется [26, 30].

3. РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Тестирование разработанных модели и соответствующего программного обеспечения проводилось на примере расчета передаточных и выходных ВАХ двухзатворного ПГТ [31] с Al2O3 в качестве верхнего подзатворного диэлектрика при комнатной температуре. Как следует из рис. 2 и 3, было достигнуто хорошее согласование результатов моделирования с экспериментальными данными как для передаточных характеристик11, так и для выходных характеристик ПГТ. Это позволяет с большим доверием относится к приводимым далее теоретическим результатам.

Рис. 2.

Передаточные ВАХ двухзатворного ПГТ со следующими диэлектриками верхнего затвора: 1 – HfO2; 2 – SiC; 3 – Al2O3; 4 – h-BN; 5 – SiO2; 6 – экспериментальные данные для ПГТ с диэлектриком Al2O3.

Рис. 3.

Выходные ВАХ двухзатворного ПГТ со следующими диэлектриками верхнего затвора: 1 – HfO2; 2 – SiC; 3 – Al2O3; 4 – h-BN; 5 – SiO2; 6 – экспериментальные данные для ПГТ с диэлектриком Al2O3.

На рис. 2 и 3 представлены расчеты для пяти различных ПГТ с одинаковыми геометрическими размерами, но различающиеся материалом диэлектрика верхнего затвора. Помимо Al2O3 [32] прочие диэлектрики выбирались по принципу используемых в приводимых экспериментальных научно-исследовательских работах, посвященных графеновым транзисторам, а именно: диоксид кремния SiO2 [33], карбид кремния SiC [34], диоксид гафния HfO2 [11, 3537], гексагональный нитрид бора h-BN [38]. Материалом подложки для всех исследуемых ПГТ являлся легированный кремний, а материалом диэлектрика нижнего затвора – диоксид кремния.

На рис. 2 приведены результаты расчетов передаточных ВАХ исследуемых приборов при комнатной температуре и напряжении на стоке Vd = 0.6 В, а на рис. 3 – выходные характеристики ПГТ при фиксированном напряжении на верхнем затворе равном Vtgt = –1 В. Из рис. 3 видно насколько существенно влияние на ВАХ материала диэлектрика верхнего затвора ПГТ. Кривая 1 соответствует ПГТ с диэлектриком верхнего затвора HfO2, кривая 2 – SiC, кривая 3 – Al2O3, кривая 4 – h-BN, кривая 5 – SiO2.

В табл. 1 приведены результаты расчетов размаха тока для передаточных характеристик всех пяти ПГТ. Если следовать принципам работы полевых транзисторов, то наилучшими являются приборы с наибольшим размахом тока, т.е. отношением максимального и минимального токов стока прибора.

Таблица 1.  

Параметры передаточных характеристик ПГТ

ПГТ Idmax/Idmin
(в точке –1 В от точки Дирака)
Idmax/Idmin
(в точке –2 В от точки Дирака)
Idmax/Idmin
(в точке +1 В от точки Дирака)
Кривая 1 (HfO2) 2.45 3.00 2.43
Кривая 2 (SiC) 2.32 3.03 2.31
Кривая 3 (Al2O3) 2.25 2.97 2.24
Кривая 4 (h-BN) 2.00 2.73 2.01
Кривая 5 (SiO2) 1.83 2.51 1.76

Рассмотрим область дырочной проводимости, находящуюся левее точки Дирака. При меньшем размахе в точке –1 В и при большем размахе в точке –2 В относительно точки Дирака лучшими являются ПГТ с HfO2 (табл. 1), и далее по нисходящей SiC, Al2O3, h-BN, SiO2..

Влияние диэлектриков на ВАХ ПГТ обусловлено их различной диэлектрической проницаемостью, которая влияет на емкости затворов, а также различной энергией оптических фононов для данных материалов с расположенным на них слоем графена. Соответствующие данные для расчетов были взяты из работы [39] для диэлектриков SiO2, HfO2, SiC, h-BN и из работ [23, 40] для Al2O3, (табл. 2). Интересно заметить, что в области электронной проводимости правее точки Дирака значения размахов Idmax/Idmin близки к аналогичным значениям в области дырочной проводимости, что связано с определенной симметрией передаточных ВАХ относительно точки Дирака.

Таблица 2.  

Диэлектрические проницаемости и параметры рассеяния для графена на различных подложках

Диэлектрик HfO2 SiC Al2O3 h-BN SiO2
Диэлектрическая проницаемость $\varepsilon $ 22.0 9.7 8.9 5.09 3.9
$\hbar {{\omega }_{{opt}}}$ , мэВ 21.6 116.0 55 101.7 58.9

Асимметричность поведения передаточных характеристик в областях электронной и дырочной проводимости обусловлена главным образом контактным сопротивлением. На рис. 2 видно смещение точки Дирака вправо, если проследить поведение кривых 15 (рис. 2), для соответствующих ПГТ с различными диэлектриками верхнего затвора. Это обусловлено уменьшением диэлектрической проницаемости диэлектриков верхнего затвора, от наибольшей для HfO2$\varepsilon = $ 22.0 до наименьшей для SiO2$\varepsilon = $ 3.9 (табл. 2).

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

На рис. 4 иллюстрируется влияние $\varepsilon $ диэлектрика нижнего затвора на передаточные ВАХ двухзатворного ПГТ. На рис. 4 наглядно проявляется смещение точки Дирака в область более высоких напряжений по мере увеличения $\varepsilon $ диэлектрика нижнего затвора.

Рис. 4.

Передаточные ВАХ двухзатворного ПГТ при различных ${{\varepsilon }_{{}}}$ диэлектрика нижнего затвора: 1 – SiO2 (${{\varepsilon }_{b}}$ = 3.9); 2 – h-BN (${{\varepsilon }_{b}}$= 5.09); 3 – Al2O3 (${{\varepsilon }_{b}}$ = 8.9); 4 – SiC (${{\varepsilon }_{b}}$ = 9.7); 5 – HfO2 (${{\varepsilon }_{b}}$ = 22.0).

Для сравнения рассмотрим изменение диэлектрической проницаемости верхнего затвора (рис. 5). Здесь наблюдается обратная зависимость, т.е. чем больше диэлектрическая проницаемость, тем сильнее смещается точка Дирака в область более низких напряжений, а ток стока и соотношение максимального тока к минимальному растут. Изменение по току гораздо более заметно.

Рис. 5.

Передаточные ВАХ двухзатворного ПГТ при различных ${{\varepsilon }_{{}}}$ диэлектрика верхнего затвора: 1 – SiO2 (${{\varepsilon }_{t}}$ = 3.9); 2 – h-BN (${{\varepsilon }_{t}}$ = 5.09); 3 – Al2O3 (${{\varepsilon }_{t}}$ = 8.9); 4 – SiC (${{\varepsilon }_{t}}$ = 9.7); 5 – HfO2 (${{\varepsilon }_{t}}$ = 22.0).

Таким образом, видно, что влияние диэлектрика верхнего затвора на ток стока сильнее, чем диэлектрика нижнего затвора, что объясняется значительным различием в емкостях затворов, которые вычисляются на основе $\varepsilon $ диэлектрика соответствующего затвора.

Рассмотрим пример ПГТ с верхним подзатворным диэлектриком HfO2 и нижним – SiO2. В этом случае емкость на единицу площади верхнего затвора составляет 1.2 × 10–2 Ф/м2, а нижнего – 1.5 × 10–4 Ф/м2. Различие этих емкостей на два порядка влияет на вычисляемый в дальнейшем электростатический потенциал в канале (см. формулу (1)), а следовательно, и на ток стока, т.е. ВАХ прибора. Смещение графиков влево или вправо также определятся значением электростатического потенциала, и, следовательно, соответствующих параметров (кроме диэлектрических проницаемостей, важно значение прикладываемого смещения к верхнему и нижнему затворам, а также напряжение Дирака для верхнего затвора). Как показали исследования, при неизменных прочих параметрах, если $\varepsilon $ диэлектрика верхнего затвора больше, чем $\varepsilon $ диэлектрика нижнего затвора, смещение точки Дирака на передаточных ВАХ будет вправо, и чем больше это отличие, тем сильнее смещение (рис. 5).

Из рис. 6 видно, что изменение рассеяния на оптических фононах оказывает незначительное влияние на характеристики приборов. Значения рассеяния на оптических фононах для данных расчетов взяты из различных источников: кривая 1 соответствует $\hbar {{\omega }_{{opt}}}$ = 55.0 мэВ [40], кривая 2 $\hbar {{\omega }_{{opt}}}$ = = 75.0 мэВ [24]. Значение $\hbar {{\omega }_{{opt}}}$ = 85.0 мэВ (кривая 3) приведено для сравнения. Как следует из проведенных расчетов, влияние этого параметра следует рассматривать совместно с диэлектрической проницаемостью материала.

Рис. 6.

Передаточные ВАХ двухзатворного ПГТ с различными значениями рассеяния на оптических фононах: 1$\hbar {{\omega }_{{opt}}}$ = 55 мэВ [41]; 2$\hbar {{\omega }_{{opt}}}$ = 75 мэВ [24]; 3$\hbar {{\omega }_{{opt}}}$ = 85 мэВ.

Влияние оптических фононов для различных материалов нижнего затвора проиллюстрировано на рис. 7 и 8. В данных примерах использован верхний затвор из Al2O3, (${{\varepsilon }_{t}}$ = 8.9), а диэлектрики нижних затворов изменяются. Было оценено влияние изменения, вызванное рассеянием на оптических фононах для различных материалов относительно соответствующего Al2O3, т.е. $\hbar {{\omega }_{{opt}}}$ = 55.0 мэВ. Для каждого из материалов это влияние незначительно (например, для кривых 3 и 4 рис. 7, соответствующих нижнему диэлектрику SiO2). Это объясняется малым отличием рассеяния на оптических фононах для SiO2 ($\hbar {{\omega }_{{opt}}}$ = 58,9 мэВ) и для Al2O3 ($\hbar {{\omega }_{{opt}}}$= 55.0 мэВ). Наибольшее отличие имеет место для кривых 1 и 2 рис. 7, соответствующих нижнему диэлектрику HfО2. Отметим, что так как $\hbar {{\omega }_{{opt}}}$= 21.6 мэВ для HfО2 отличается от Al2O3 с $\hbar {{\omega }_{{opt}}}$ = 55.0 мэВ в меньшую сторону, то и токи стока соответственно уменьшаются, в отличие от всех других случаев (кривые 1 и 2 рис. 8 для SiC, кривые 3 и 4 рис. 8 для h-BN, кривые 2 и 3 рис. 7 для SiO2), для которых ток стока для соответствующего материала увеличивается.

Рис. 7.

Передаточные ВАХ двухзатворного ПГТ с рассеянием на оптических фононах соответствующие различным диэлектрикам нижнего затвора: 1${{\varepsilon }_{b}}$ = 22.0 (HfО2), $\hbar {{\omega }_{{opt}}}$ = 55 мэВ; 2${{\varepsilon }_{b}}$ = 22.0 (HfО2), $\hbar {{\omega }_{{opt}}}$ = 21.6 мэВ; 3${{\varepsilon }_{b}}$ = 3.9 (SiO2), $\hbar {{\omega }_{{opt}}}$= 58.9 мэВ; 4${{\varepsilon }_{b}}$ = 3.9 (SiO2), $\hbar {{\omega }_{{opt}}}$= 55.0 мэВ.

Рис. 8.

Передаточные ВАХ двухзатворного ПГТ с рассеянием на оптических фононах соответствующие различным диэлектрикам нижнего затвора: 1${{\varepsilon }_{b}}$ = = 9,7 (SiC), $\hbar {{\omega }_{{opt}}}$ = 116.0 мэВ; 2${{\varepsilon }_{b}}$ = 9,7 (SiC), $\hbar {{\omega }_{{opt}}}$ = 55.0 мэВ; 3${{\varepsilon }_{b}}$ = 5.09 (h-BN), $\hbar {{\omega }_{{opt}}}$ = 101.7 мэВ; 4${{\varepsilon }_{b}}$ = 5.09 (h-BN), $\hbar {{\omega }_{{opt}}}$ = 55.0 мэВ.

Было также проведено исследование влияния толщины верхнего подзатворного диэлектрика на передаточные характеристики ПГТ с диэлектриком Al2O3. Его толщина варьировалась от 15 до 30 нм с шагом 5 нм, что соответствует кривым с 1 по 4 на рис. 9. Как следует из рис. 9, точка Дирака с увеличением толщины диэлектрика смещается в область более высоких значений напряжений на верхнем затворе, а размах токов снижается. Это связано с тем, что увеличение толщины диэлектрика приводит к уменьшению емкости верхнего затвора, и как следствие к изменению электростатического потенциала в канале ПГТ, который в свою очередь влияет на ВАХ прибора.

Рис. 9.

Передаточные ВАХ двухзатворного ПГТ с различными толщинами диэлектрика верхнего затвора: 1dtg = 15 нм; 2dtg = 20 нм; 3dtg = 25 нм; 4dtg = 30 нм.

На рис. 10 приведены передаточные ВАХ двухзатворного ПГТ с различными толщинами диэлектрика нижнего затвора. Толщина диэлектрика варьировалась от 300 до 420 нм с шагом 60 нм, что соответствует кривым с 1 по 3. Установлено, что влияние на ВАХ толщины нижнего затвора незначительно. Видно, что с увеличением толщины диэлектрика точка Дирака немного смещается в область более низких значений напряжения на верхнем затворе, а размах токов на характеристиках остается фактически неизменным. Это связано с тем, что емкость на единицу площади нижнего затвора значительно меньше соответствующей емкости верхнего затвора, а следовательно, оказывает меньшее влияние на изменение электростатического потенциала в канале, и в результате на передаточные ВАХ ПГТ.

Рис. 10.

Передаточные ВАХ двухзатворного ПГТ с различными толщинами диэлектрика нижнего затвора: 1dbg = 300 нм; 2dbg = 360 нм; 3dbg = 420 нм.

4. ЗАКЛЮЧЕНИЕ

Таким образом, с помощью предложенной упрощенной комбинированной самосогласованной модели ПГТ показано существенное влияние материала диэлектрика верхнего затвора на выходные и передаточные вольт-амперные характеристики полевых транзисторов на основе однослойного графена. Проведен анализ рассчитанных характеристик пяти ПГТ, различающихся диэлектриком верхнего затвора. Установлено, что наибольшим соотношением максимального тока стока к минимальному на ВАХ обладает ПГТ с HfO2, а наименьшим – ПГТ с SiO2. В целом, результаты для всех пяти ПГТ приводят к выводу о приемлемых рабочих характеристиках исследованных транзисторов.

Установлено, что изменение толщины диэлектрика нижнего затвора (в пределах десятков нанометров) незначительно влияет на ВАХ ПГТ, но даже небольшое изменение толщины диэлектрика верхнего затвора (несколько единиц нанометров) приводит к видимым изменениям характеристик прибора, что очевидно будет важным при изготовлении приборов.

Программы, реализующие модель, включены в систему моделирования наноэлектронных устройств NANODEV [41, 42], разрабатываемую в БГУИР с 1995 года.

Работа выполнена при финансовой поддержке Государственной программы научных исследований “Конвергенция” Республики Беларусь и была частично доложена на 29-ой Международной Крымской конференции “СВЧ-техника и телекоммуникационные технологии”, Севастополь, Россия, 8–14 сентября 2019 г.

Список литературы

  1. Ferrari A.C., Bonaccorso F., Fal’ko V., Novoselov K.S. et al. Science and technology roadmap for graphene, related two-dimensional crystals, and hybrid systems // Nanoscale. 2015. V. 7. P. 4598–4810.

  2. Ратников П.В., Силин А.П. Двумерная графеновая электроника: современное состояние и перспективы // Успехи физических наук. 2018. Т. 188. № 12. С. 1249–1287.

  3. Beyond CMOS. International Roadmap for Devices and Systems. 2018.

  4. Lu N., Wang L., Li L., Liu M. A review for compact model of graphene field-effect transistors // Chin. Phys. B. 2017. V. 26. № 3. P. 036804-1–18.

  5. Zhuang F., Deng W., Ma X., Huang J. A model of carrier density and drain current for monolayer graphene field-effect transistors // AIP Advances. 2019. V. 9. P. 025222-1–7.

  6. Aguirre-Morales J.-D., Frégonèse S., Mukherjee C., Wei W., Happy H., Maneux C., Zimmer T. A Large-signal monolayer graphene field-effect transistor compact model for RF-circuit applications // IEEE Trans. Electron Dev. 2017. V. 64. № 10. P. 4302–4309.

  7. Meric I., Han M.Y., Young A.F., Ozyilmaz B., Kim P., Shepard K.L. Current saturation in zero-bandgap, top-gated graphene field-effect transistors // Nature Nanotech. 2008. V. 3. № 11. P. 654–659.

  8. Scott B., Leburton J. Modeling of the output and transfer characteristics of graphene field-effect transistors // Nanotechnology. IEEE Trans. 2011. V. 10. № 5. P. 1113–1119.

  9. Wang L., Peng S., Wang W., Xu G., Ji Z., Lu N., Li L., Jin Z., Liu M. Surface-potential-based physical compact model for graphene field effect transistor // J. Appl. Phys. 2016. V. 120. № 8. P. 084509-1–6.

  10. Nguyen V.H., Bournel A., Chassat C., Dollfus P. Quantum transport of Dirac fermions in graphene field effect transistors // International Conference on Simulation of semiconductor processes and devices (SISPAD). 2010. P. 9–12.

  11. Pugnaghi C., Grassi R., Gnudi A., Lecce V.D., Gnani E., Reggiani S., Baccarani G. Semianalytical quantum model for graphene field-effect transistors // J. Appl. Phys. 2014. V. 116. № 11. P. 114505-1–5.

  12. Szafranek B.N., Fiori G., Schall D., Neumaier D., Kurz H. Current saturation and voltage gain in bilayer graphene field effect transistors // Nano lett. 2012. V. 12. № 3. P. 1324–1328.

  13. Pasadas F., Jiménez D. Large-signal model of the bilayer graphene field-effect transistor targeting radio-frequency applications: Theory versus experiment // J. Appl. Phys. 2015. V. 118. № 24. P. 244501-1–7.

  14. Umoh I.J., Kazmierski T.J., Al-Hashimi B. Multilayer graphene FET compact circuit-level model with temperature effects // IEEE Trans. Nanotechnology. 2014. V. 13. № 4. P. 805–813.

  15. Свинцов Д.А., Вьюрков В.В., Лукичёв В.Ф., Орликовский А.А., Буренков А., Охснер Р. Туннельные полевые транзисторы на основе графена // Физика и техника полупроводников. 2013. Т. 47. Вып. 2. С. 244–250.

  16. Абрамов И.И. Проблемы и принципы физики и моделирования приборных структур микро- и наноэлектроники. VIII. Нанотранзисторы с МДП-структурой // Нано- и микросистемная техника. 2010. № 9. С. 27–37. № 10. С. 28–41. № 11. С. 29–42.

  17. Абрамов И.И. Основы моделирования элементов микро- и наноэлектроники. LAP LAMBERT Academic Publishing, Saarbrücken, Germany, 2016. 444 с.

  18. Абрамов И.И. Лекции по моделированию элементов интегральных схем. М.–Ижевск: НИЦ РХД, 2005. 152 с.

  19. Абрамов И.И., Коломейцева Н.В., Лабунов В.А., Романова И.А., Щербакова И.Ю. Моделирование полевых графеновых транзисторов с одним и двумя затворами в различных режимах функционирования // Нанотехнологии, разработка, применение: XXI век. 2018. № 3. С. 16–24.

  20. Абрамов И.И., Коломейцева Н.В., Лабунов В.А., Романова И.А., Щербакова И.Ю. Моделирование передаточных характеристик двухзатворных полевых графеновых транзисторов // Нано- и микросистемная техника. 2018. № 11. С. 643–650.

  21. Abramov I., Labunov V., Kolomejtseva N., Romanova I., Shcherbakova I. Simulation of graphene field-effect transistors and resonant tunneling diodes based on carbon nanomaterials // Proc. SPIE. Editors Vladimir F. Lukichev, Konstantin V. Rudenko. 2019. V. 11022. P. 110220F-1–11.

  22. Thiele S.A., Schaefer J.A., Schwierz F. Modeling of graphene metal-oxide-semiconductor field-effect transistors with gapless large-area graphene channels // J. Appl. Phys. 2010. V. 107. № 9. P. 094505-1–8.

  23. Tian J., Katsounaros A., Smith D., Hao Y. Graphene field-effect transistor model with improved carrier mobility analysis // IEEE Trans. Electron. Dev. 2015. V. 62. № 10. P. 3433–3440.

  24. Landauer G.M., Jimenez D., Gonzalez J.L. An accurate Verilog-A compatible compact model for graphene field-effect transistors // IEEE Trans. Nanotechnol. 2014. V. 13. № 5. P. 895–904.

  25. Абрамов И.И., Коломейцева Н.В., Лабунов В.А., Романова И.А. Моделирование полевых графеновых транзисторов с одним и двумя затворами // Нано- и микросистемная техника. 2017. Т. 19. № 12. P. 714–721.

  26. Abramov I.I., Labunov V.A., Kolomejtseva N.V., Romanova I.A. Simulation of field-effect transistors and resonant tunneling diodes based on graphene // Proc. SPIE. Editors Vladimir F. Lukichev, Konstantin V. Rudenko. 2016. V. 10224. P. 102240V-1–10.

  27. Zhu W., Linghu C., Zhang J., Zhang L., Yu Z. Compact model of graphene field effect transistors and its application in circuit simulation of RF mixer consisting of GFETs and CMOS // International Conference on Simulation of semiconductor processes and devices (SISPAD). 2012. P. 79–82.

  28. Абрамов И.И., Харитонов В.В. Методы и алгоритмы трехмерного численного моделирования полупроводниковых приборов и структур // Электронное моделирование. 1990. Т. 12. № 5. С. 39–45.

  29. Абрамов И.И., Харитонов В.В. Численное моделирование физических процессов в элементах и фрагментах кремниевых БИС и СБИС // Электронное моделирование. 1992. Т. 14. № 4. С. 49–56.

  30. Абрамов И.И., Коломейцева Н.В., Лабунов В.А., Романова И.А. Моделирование электрических характеристик полевых графеновых транзисторов с использованием упрощенной комбинированной модели // Материалы 26-й Международной Крымской конференции “СВЧ-техника и телекоммуникационные технологии”, 2016, Россия, Севастополь. Т. 7. С. 1534–1540.

  31. Wang H., Hsu A., Antoniadis D.A., Palacios T. Compact virtual-source current-voltage model for top- and back-gated graphene field-effect transistors // IEEE Trans. Electron. Dev. 2011. V. 58. № 5. P. 1523–1533.

  32. ang H., Qin S., Zheng X., Wang G., Tan Y., Peng G., Zhang X. An Al2O3 gating substrate for the greater performance of field effect transistors based on two-dimensional materials // Nanomaterials. 2017. V. 7. № 10. P. 286–294

  33. Mackin C.E. Electrolyte-gated graphene field-effect transistors: modeling and applications. Thesis. Massachusetts Institute of Technology, Department of Electrical Engineering and Computer Science, 2015. 55 p.

  34. Moon J., Curtis D., Bui S., Hu M., Gaskill D., Tedesco J., Asbeck P., Jernigan G., VanMil B., Myers-Ward R., Eddy C., Campbell P., Weng X. Top-gated epitaxial graphene FETs on Si-face SiC wafers with a peak transconductance of 600 ms/mm // IEEE Electron Dev. Lett. 2010. V. 31. № 4. P. 260–262.

  35. Kedzierski J., Hsu P.-L., Reina A., Kong J., Healey P., Wyatt P., Keast C. Graphene-on-insulator transistors made using c on Ni chemical-vapor deposition // IEEE Electron Dev. Lett. 2009. V. 30. № 7. P. 745–747.

  36. Petrone N., Meric I., Chari T., Shepard K.L., Hone J. Graphene field-effect transistors for radio-frequency flexible electronics // J. Electron. Dev. Society. 2015. V. 3. № 1. P. 44–48.

  37. Wei W., Belhaj M., Deokar G., Mele D. Fabrication and characterization of CVD-grown graphene based field-effect transistor // Proc. 44th Eur. Microw. Conf. (EuMC). 2014. P. 367–370.

  38. Meric I., Dean C., Young A., Hone J., Kim P., Shepard K.L. Graphene field-effect transistors based on boron nitride gate dielectrics // Proc. IEEE Electron Dev. Meeting. 2010. P. 23.2.1–23.2.4.

  39. Perebeinos V., Avouris P. Inelastic scattering and current saturation in graphene // Phys. Rev. B. 2010. V. 81. № 19. P. 195442-1–8.

  40. Konar A., Fang T., Jena D. Effect of high-k gate dielectrics on charge transport in graphene-based field-effect transistors // Phys. Rev. B. Condens. Matter Mater. Phys. 2010. V. 82. P. 115452-1–3.

  41. Абрамов И.И., Гончаренко И.А., Игнатенко С.А., Королев А.В., Новик Е.Г., Рогачев А.И. Система моделирования наноэлектронных приборов – NANODEV // Микроэлектроника. 2003. Т. 32. № 2. С. 124–133.

  42. Abramov I.I., Baranoff A.L., Goncharenko I.A., Kolomejtseva N.V., Bely Y.L., Shcherbakova I.Y. A nanoelectronic device simulation software system NANODEV: New opportunities // Proc. of SPIE. Editors Kamil A. Valiev, Alexander A. Orlikovsky. 2010. V. 7521. P. 75211E1-1–11.

Дополнительные материалы отсутствуют.