删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

纳米线环栅隧穿场效应晶体管的电容模型

本站小编 Free考研考试/2021-12-29

摘要:纳米线环栅隧穿场效应晶体管相比于其他多栅器件具有更强的短沟道效应抑制能力及更优异的电学特性. 器件模型能够模拟器件电学特性, 对于器件及电路的实际应用极为关键. 目前, 已有纳米线环栅隧穿场效应晶体管的电流模型报道, 但是尚没有电容模型的相关报道. 电容模型主要用于瞬态特性模拟, 对于评估电路速度转换和频率特性至关重要. 由于没有可用的电容模型, 纳米线环栅隧穿场效应晶体管电路方面的研究主要通过数值迭代的方法开展, 该方法不仅对硬件平台要求高, 且耗时长, 还容易出现收敛性问题, 只能勉强用于极小规模电路模块, 对于包含晶体管数目较多的电路无能为力. 本文针对以上问题, 从基本的器件物理出发, 建立了纳米线环栅隧穿场效应晶体管的电容模型, 该模型不涉及任何数值迭代过程. 相比于数值模型, 该模型计算速度快、过程稳定, 能够加速纳米线环栅隧穿场效应晶体管器件及电路的相关研究.
关键词: 隧穿场效应晶体管/
带带隧穿/
纳米线/
电容模型

English Abstract


--> --> -->
金属氧化物场效应晶体管(metal-oxide-semiconductor field-effect-transistor, MOSFET)的特征尺寸随着摩尔定律不断缩减, 使得集成电路集成度不断提高的同时, 器件短沟道效应也越来越严重, 导致其静态功耗已超过动态功耗成为集成电路的主要功耗来源. 隧穿场效应晶体管(tunneling field-effect-transistor, TFET)是1种利用量子隧穿机制导通电流的新型器件, 其可以突破MOSFET室温下60 mV/dec 亚阈值摆幅极限, 实现更加陡峭的开关曲线, 同时降低泄漏电流及工作电压, 进一步降低集成电路功率损耗[1-3]. 目前TFET 已成为微电子领域的研究热点, 被认为是后摩尔时代最有可能取代MOSFET推动高性能、低功耗集成电路持续发展的技术之一[4].
尽管TFET能够实现十分陡峭的开关特性, 但是TFET也面临短沟道效应(short channel effects, SCEs)和漏感应势垒降低(drain-induced barrier lowering, DIBL)等非理想效应[5,6], 这些效应导致TFET特性退化, 阻碍了TFET器件及电路的实际应用. 为了解决这些问题, 研究人员在TFET结构上做出了很多创新, 其中基于纳米线结构的环栅TFET (gate-all-around TFET, GAA-TFET)具有最高的沟道电势调控效率, 能够在很大程度上抑制SCEs和DIBL等非理想效应[7], 同时由于栅电极对隧穿结的控制增强, GAA-TFET比其他多栅结构表现出了更低的亚阈值摆幅和更高的隧穿电流[6].
鉴于GAA-TFET对非理想效应出色的抑制能力和优异的器件特性, 人们针对GAA-TFET的各个方面开展了大量研究, 其中器件模型可用于计算器件及电路特性, 对于GAA-TFET器件及电路的实际应用极为关键. 器件模型包括电流模型和电容模型, 其中电流模型可用于计算电路稳态支路电流和节点电压, 而电容模型可用于计算电路瞬态特性, 评估电路的速度转换和频率特性. 目前, 针对GAA-TFET的电流模型已有报道[7-12], 但是尚未见其电容模型的报道. 由于缺乏GAA-TFET电容模型, 目前研究人员主要通过数值迭代的方法开展GAA-TFET电路设计方面的研究. 但是, 数值迭代对于硬件计算平台要求高, 而且除了计算速度慢, 耗时长外, 还容易出现收敛性问题, 仅能够勉强用于只包含几个器件的极小规模电路模块, 对于包含成百上千个器件的电路, 数值迭代方法则无能为力, 这也是目前关于GAA-TFET电路研究方面报道很少的原因之一.
本文针对以上问题, 从基本的器件物理出发, 建立了GAA-TFET的电容模型, 且该模型不涉及任何数值迭代过程. 相比于数值模型, 该模型计算速度快, 计算过程稳定, 能够加速GAA-TFET器件及电路的相关研究.
2
2.1.电势模型
-->图1(a)为一个n型GAA-TFET的三维结构示意图, 其中纳米线的半径R = 7.5 nm, 采用HfO2作为栅氧化层, 介电常数εox = 22, 厚度为Tox = 2 nm, 栅金属功函数为4.2 eV. 图1(b)是GAA-TFET沿纳米线直径方向的剖面图, 本文主要关注器件核心物理机理的分析, 在模型推导过程中, 暂不考虑短沟道效应等二次非理想效应的影响, 因此器件栅长设置相对较长, 为Lg = 50 nm[7]. 源区、沟道和漏区的掺杂浓度分别为NS = 1 × 1020 cm–3(P型), NC = 1 × 1015 cm–3 (N型) 和ND = 5 × 1019 cm–3 (N型). 器件剖面图可划分为3个区域, 分别为源耗尽区I区、沟道耗尽区II区以及沟道积累区III区, 其中I区和II区的宽度分别为LILII.
图 1 (a) GAA-TFET的三维结构示意图; (b) GAA-TFET沿沟道方向的剖面示意图; (c) GAA-TFET垂直于沟道方向的剖面示意图
Figure1. (a) Three-dimensional structure of GAA TFETs. Schematic cross section of an n-type GAA TFET (b) along the channel and (b) normal to the channel direction.

在柱状坐标系下, I区和II区的二维泊松方程如下:
$ \frac{1}{r}\frac{\partial }{{\partial r}}\left( {r\frac{{\partial {\psi _j}\left( {r,z} \right)}}{{\partial r}}} \right) + \frac{{{\partial ^2}{\psi _j}\left( {r,z} \right)}}{{\partial {z^2}}} = - \frac{{q{N_j}}}{{{\varepsilon _{{\text{Si}}}}}} , $
其中下标j = I, II; q为电子电荷; ψj (r, z)为二维电势分布; εSi为沟道介电常数; Nj为掺杂浓度. 通过抛物近似[13], 沿纳米线直径方向的电势可以写为
$ {\psi _j}\left( {r,z} \right) = {\varphi _{m,j}}\left( z \right) + {p_{1,j}}\left( z \right)r + {p_{2,j}}\left( z \right){r^2} , $
式中φm, j (z)为r = 0处沿着沟道中心线方向的电势分布; φm,j (z), p1, j (z)和p2, j (z) 是未知参量, 可通过边界条件求得. 由于器件的对称结构, 可得沿直径方向的电场在r = 0处等于0, 由此可得
$ \left\{\begin{aligned}&{\frac{\partial {\psi }_{j}\left(r,z\right)}{\partial r}\bigg|}_{r=0}=0\text{, }\\&{{\psi }_{j}\left(r,z\right)|}_{r=R}={\varphi }_{\text{S},j}\left(z\right)\text{, }\end{aligned}\right. $
其中φS, j (z)为r = R处沿沟道方向的电势分布, 即表面电势.
此外, 在垂直氧化层与硅界面的方向上电位移矢量保持连续, 为求解氧化层中靠近沟道界面处的电场, 图1(c)给出了垂直沟道方向的截面图, 对于任意半径为r的闭环, 假设电场为E(r), 沟道电荷为Q, 根据高斯定理, 得2πrLgE(r) εox = Q. 对于某一固定栅电压, Q是1个常数, 因此对等式两边进行微分, 可得
$ E\left(r\right)+r\frac{\text{d}E\left(r\right)}{\text{d}r}=0\text{, } $
其通解为
$ E\left( r \right) = {F}/{r} , $
其中F为常数. 沟道表面电势与栅电极之间的电势差为
$ {V_{\text{G}}} - {V_{{\rm{FB}}}} - {\varphi _{\text{S}}}\left( z \right) = \int_R^{R + {T_{{\text{ox}}}}} {\frac{F}{r}{\text{d}}r} , $
式中VFB为平带电压, VG为栅电压. 由(6)式求得F并代入(5)式可得
$ E\left( r \right) = \frac{{{V_{\text{G}}} - {V_{{\rm{FB}}}} - {\varphi _{\text{S}}}\left( z \right)}}{{r\ln \left( 1 +T_{\rm{ox}}/R \right)}} . $
由此可以得到栅氧化层中沟道附近的电场为$ E(R)?= (V_{\rm G} - V_{\rm FB} - \varphi_{\rm S} (z))/t_{\rm ox} $, 其中$ t_{\rm ox} = R?\ln(1+ $$ T_{\rm ox}/R) $相当于平面器件中的栅氧化层厚度. 需要注意的是, 在平面器件中, 考虑栅氧化层边界泄漏电场对I区的影响, I区可被等效为环绕了一层厚度为πtox/2的氧化层[14]. 由此根据电位移矢量连续可得
$ {\left. {\frac{{\partial {\psi _j}\left( {r,z} \right)}}{{\partial r}}} \right|_{r = R}} = \frac{{{\varepsilon _{{\text{ox}}}}}}{{{\varepsilon _{{\text{Si}}}}}}\left( {\frac{{{V_{\text{G}}} - {V_{{\rm{FB}},j}} - {\varphi _{{\text{S}},j}}\left( z \right)}}{{{t_{{\text{ox}},j}}}}} \right) , $
其中对于I区, tox,j = πtox/2, 对于II区, tox,j = tox. 联立求解(3)式和(8)式可得
$ \left\{ \begin{aligned} & {p_{1,j}}\left( z \right) = 0, \\ & {p_{2,j}}\left( z \right) = \frac{{{\varepsilon _{{\text{ox}}}}}}{{2R{\varepsilon _{{\text{Si}}}}}}\left( {\frac{{{V_{\text{G}}} - {V_{{\rm{FB}},j}} - {\varphi _{{\text{S}},j}}\left( z \right)}}{{{t_{{\text{ox}},j}}}}} \right), \\ & {\varphi _{m,j}}\left( z \right) = {\varphi _{{\text{S}},j}}\left( z \right) - {p_{2,j}}\left( z \right){R^2}.\end{aligned} \right. $
把(2)式与(9)式代入(1)式, 即可得到简化的一维泊松方程:
$ \frac{{{{\text{d}}^2}{\varphi _{{\text{S}},j}}\left( z \right)}}{{{\text{d}}{z^2}}} - \frac{{{\varphi _{{\text{S}},j}}\left( z \right) - \left( {{V_{\text{G}}} - {V_{{\text{FB}},j}}} \right)}}{{\lambda _j^2}} = - \frac{{q{N_j}}}{{{\varepsilon _{{\text{Si}}}}}} , $
$\lambda_{\rm I} = [R\varepsilon_{\rm Si} \pi t_{\rm ox}/(4\varepsilon_{\rm ox})]^{1/2}, ~~ \lambda_{\rm II} = [R\varepsilon_{\rm Si} t_{\rm ox}/(2\varepsilon_{\rm ox})]^{1/2},$ (10)式的通解为
$\begin{split}{\varphi _{{\text{S}},j}}\left( z \right) =\;& {A_j}\exp \left( { - \frac{z}{{{\lambda _j}}}} \right) + {B_j}\exp \left( { + \frac{z}{{{\lambda _j}}}} \right) \\&+ \left( {{V_{\text{G}}} - {V_{{\text{FB}},j}} - \frac{{q{N_j}}}{{{\varepsilon _{{\text{Si}}}}}}\lambda _j^2} \right) , \end{split}$
其中AjBj为待求解系数. 对于I区, 将(11)式在z = –LI附近展开并仅保留常数项与二次项, 可以得到I区的表面电势为
$\begin{split}{\varphi }_{\text{S,I}}\left(z\right)=\;&C\left[1+\frac{1}{2{\lambda }_{\text{I}}^{2}}{\left(z+{L}_{\text{I}}\right)}^{2}\right]\\&+\left({V}_{\text{G}}-{V}_{\text{FB,I}}-\frac{q{N}_{\text{I}}}{\varepsilon {}_{\text{Si}}}{\lambda }_{\text{I}}^{2}\right) , \end{split}$
其中C未知, 可以通过边界条件$\varphi_{\rm S,I}( - L_{\rm I})= V_{\rm BS}$求得, $V_{\rm BS} = V_{\rm S} - V_{\rm t} \log(N_{\rm I}/n_{i,\rm I})$为中性源区电势, VS为源电极电压, Vt是热电压, ni, I是本征载流子浓度, 求得C
$ C = {V_{{\text{BS}}}} - \left( {{V_{\text{G}}} - {V_{{\text{FB,I}}}} - \frac{{q{N_{\text{I}}}}}{{\varepsilon {}_{{\text{Si}}}}}\lambda _{\text{I}}^2} \right). $
将(13)式代入(12)式求得I区表面电势分布为
$ \left\{ \begin{aligned} & {\varphi _{{\text{S,I}}}}\left( z \right) = \frac{{q{N_{{\text{eff}}}}}}{{2{\varepsilon _{{\text{Si}}}}}}{\left( {z + {L_{\text{I}}}} \right)^2} + {V_{{\text{BS}}}},\\& {N_{{\text{eff}}}} = {N_{\text{I}}} - \frac{{{\varepsilon _{{\text{Si}}}}}}{{q\lambda _{\text{I}}^{\text{2}}}}\left( {{V_{\text{G}}} - {V_{{\text{FB,I}}}} - {V_{{\text{BS}}}}} \right). \end{aligned} \right. $
对于II区, II区与III区边界面上的电场和电势保持连续, 然而需要注意的是III区为积累区, 在III区中存在高浓度的电子, 这些电子在很大程度上屏蔽了栅、漏电场的作用, 因此III区中的电场很小, 为了简化模型推导, 假设III区中电场为零, 即可求得系数AIIBII:
$ \left\{\begin{aligned}&{A}_{\rm{II}}=\frac{1}{2}\left[{\varphi }_{\text{ch}}-\left({V}_{\text{G}}-{V}_{\text{FB,}\rm{II}}-\frac{{q}{N}_{\rm {II}}{\lambda }_{\rm {II}}^{2}}{{\varepsilon }_{\text{Si}}}\right)\right]\\ & ~~~~~~\times \text{exp}\left(+\frac{{L}_{\rm{II}}}{{\lambda }_{\rm{II}}}\right),\\ &{B}_{\rm{II}}=\frac{1}{2}\left[{\varphi }_{\text{ch}}-\left({V}_{\text{G}}-{V}_{\text{FB,}\rm{II}}-\frac{{q}{N}_{\rm{II}}{\lambda }_{\rm{II}}^{2}}{{\varepsilon }_{\text{Si}}}\right)\right]\\ &~~~~~~\times \text{exp}\left(-\frac{{L}_{\rm{II}}}{{\lambda }_{\rm{II}}}\right),\end{aligned}\right. $
式中φch是III区表面电势, 将(15)式代入(11)式, 求得II区表面电势分布为
$\begin{split}&{\varphi }_{\text{S,}\text{II}}\left(z\right)=\left({V}_{\text{G}}-{V}_{\text{FB,}\text{II}}-\frac{{q}{N}_{\text{II}}{\lambda }_{\text{II}}^{2}}{{\varepsilon }_{\text{Si}}}\right)\\&+\left[{\varphi }_{\text{ch}}-\left({V}_{\text{G}}-{V}_{\text{FB,}\text{II}}-\frac{{q}{N}_{\text{II}}{\lambda }_{\text{II}}^{2}}{{\varepsilon }_{\text{Si}}}\right)\right]\text{cosh}\left(\frac{z-{L}_{\text{II}}}{{\lambda }_{\text{II}}}\right) . \end{split}$
对于III区, 其表面电势φch可写为[15]
$ \begin{split}{\varphi }_{\text{ch}}=\;&F+{V}_{\text{t}}\text{ln}\Bigg(\frac{1}{{V}_{\text{t}}}\Bigg\{{V}_{\text{t}}+\frac{\sqrt{F}}{\sqrt{F}+\gamma }\left({V}_{\text{G}}-{V}_{\text{FB,}\text{II}}-F\right)\\ &+\frac{1}{2}\left[\frac{\sqrt{F}}{{\left(\sqrt{F}+\gamma \right)}^{2}}-\frac{\gamma \left(F-2\right)}{2{\left(\sqrt{F}+\gamma \right)}^{3}}\right]\\&\times{\left({V}_{\text{G}}-{V}_{\text{FB}\text{, }\text{II}}-F\right)}^{2}\Bigg\}\Bigg),\\[-18pt]\end{split} $
其中Fφch,dep
$ F = \frac{1}{2}\Big[ {{V_{\text{D}}} + \phi + {\varphi _{{\text{ch,dep}}}} - \sqrt {{{\left( {{\varphi _{{\text{ch,dep}}}} - {V_{\text{D}}} - \phi } \right)}^2} + {\delta ^2}} }\, \Big], $
$ {\varphi }_{\text{ch,dep}}={\left[\sqrt{{V}_{\text{G}}-{V}_{\text{FB,}\text{II}}+\frac{{\gamma }^{2}}{4}}-\frac{\gamma }{2}\right]}^{2} , $
式中$\gamma = (2\varepsilon_{\rm Si} qN_{\rm II})^{1/2}/(\varepsilon_{\rm ox}/t_{\rm ox}),~\delta = 0.04$为平滑因子, $\phi = {V_{\text{t}}}\ln \left( {{N_{{\text{inv}}}}/{n_{i, {\text{II}}}}} \right)$是屏蔽电势, Ninv = 2.5 × 1019 cm–3是1个从TCAD仿真中提取出来的经验参数.
至此未知量只剩I区和II区的宽度LILII. 假设z = 0处, 即I/II界面处表面电势为φS(0), 由(14)式和(16)式可以得到LILII如下:
$ \left\{\begin{aligned}&{L}_{\text{I}}=\sqrt{\frac{2{\varepsilon }_{\text{Si}}}{q{N}_{\text{eff}}}}\sqrt{{\varphi }_{\text{S}}\left(0\right)-{V}_{\text{BS}}},\\&{L}_{\text{II}}={\lambda }_{\text{II}}{\text{cosh}}^{-1}\left[\frac{{\varphi }_{\text{S}}\left(0\right)-\left({V}_{\text{G}}-{V}_{\text{FB,}\text{II}}-\dfrac{{q}{N}_{\text{II}}{\lambda }_{\text{II}}^{2}}{{\varepsilon }_{\text{Si}}}\right)}{{\varphi }_{\text{ch}}-\left({V}_{\text{G}}-{V}_{\text{FB,}\text{II}}-\dfrac{{q}{N}_{\text{II}}{\lambda }_{\text{II}}^{2}}{{\varepsilon }_{\text{Si}}}\right)}\right].\end{aligned}\right. $
φS(0)是目前唯一的未知量, 为了求解φS(0), 可将LILII分别代入(14)式和(16)式, 利用z = 0处电势连续, 可以得到φS(0)满足如下方程:
$ \left\{ \begin{aligned} &0 = \frac{{2q{N_{{\text{eff}}}}\lambda _{\text{II}}^2}}{{{\varepsilon _{{\text{Si}}}}}}\left[ {{\varphi _{\text{s}}}\left( 0 \right) - {V_{{\text{BS}}}}} \right] - \frac{{{\varepsilon _{{\text{Si}}}}^2}}{{{\lambda _{{\text{II}}}^2}}} \bigg[{\varphi _{{\text{ch}}}} \\ & \;\;\quad - \left( {{V_{\text{G}}} - {V_{{\text{FB,II}}}} - \frac{{q{N_{{\text{II}}}}{\lambda _{{\text{II}}}}^2}}{{{\varepsilon _{{\text{Si}}}}}}} \right) \bigg]^2 {R^2}, \\ &H = \frac{{{\varphi _{\text{s}}}\left( 0 \right) - \left( {{V_{\text{G}}} - {V_{{\text{FB,II}}}} - \dfrac{{q{N_{{\text{II}}}}{\lambda ^2_{{\text{II}}}}}}{{{\varepsilon _{{\text{Si}}}}}}} \right)}}{{{\varphi _{{\text{ch}}}} - \left( {{V_{\text{G}}} - {V_{{\text{FB,II}}}} - \dfrac{{q{N_{{\text{II}}}}{\lambda^2 _{{\text{II}}}}}}{{{\varepsilon _{{\text{Si}}}}}}} \right)}}, \\ & R = \sinh \left( {{{\cosh }^{ - 1}}H} \right). \end{aligned} \right. $
借助公式变换$\sinh(\cosh^{-1} H)=(H^2 - 1)^{1/2}$, 求解(21)式, 可得φS(0)为
$ \left\{\begin{aligned}&{\varphi }_{\text{S}}\left(0\right)=\frac{-b-\sqrt{{b}^{2}-4c}}{2}\text{, }\\ &b=2\left({V}_{\text{G}}-{V}_{\text{FB,II}}-\frac{q{N}_{\text{II}}{\lambda }_{\text{II}}^2}{{\varepsilon }_{\text{Si}}}+\frac{q{N}_{\text{eff}}{\lambda }_{\text{II}}^2}{{\varepsilon }_{\text{Si}}}\right),\\&c=2{\varphi }_{\text{ch}}\left({V}_{\text{G}}-{V}_{\text{FB,II}}-\frac{q{N}_{\text{II}}{\lambda }_{\text{II}}^2}{{\varepsilon }_{\text{Si}}}\right)\\&\quad\;\;-{\varphi }_{\text{ch}}^{2}+\frac{2q{N}_{\text{eff}}{\lambda }_{\text{II}}^{2}{V}_{\text{BS}}}{{\varepsilon }_{\text{Si}}}.\end{aligned}\right. $

2
2.2.电容模型
-->在TFET中, 源区与沟道区之间的隧穿结本质上是个高势垒的反偏结, 势垒两侧载流子很难通过热发射相互流通, 因此, 可以假设源电荷QS 主要为势垒左侧的源区耗尽区电荷, 而漏电荷QD 主要为势垒右侧的电荷, 这样器件源端电荷QS 可以表示为
$ {Q}_{\text{S}}=-q{N}_{\text{eff}}\text{π}{R}^{2}{L}_{\text{Ⅰ}} . $
漏端电荷QD主要包括沟道中III区中的可动载流子以及漏区边缘电荷:
$ \begin{split}{Q}_{\text{D}}=\;&-2\text{π}R\left({L}_{\text{g}}-{L}_{\text{II}}\right)(\varepsilon_{\rm{ox}}/\varepsilon_{\rm{ox}})({V}_{\text{G}}-{V}_{\text{FB,}\text{III}}\\&-{\phi }_{\text{ch,}\text{III}} ) +\text{π}{R}^{2}{\varepsilon }_{\text{Si}}{E}_{{\max}} , \end{split}$
其中Emax = (VBD φch)/λII, VBD = VD + Vt × log(ND/ni)为中性漏区电势, VD为漏电极电压. 栅电荷QG、源电荷QS和漏电荷QD满足电荷守恒关系:
$ {Q_{\text{G}}}{\text{ + }}{Q_{\text{S}}} + {Q_{\text{D}}}{\text{ = 0}}. $
将端电荷分别对相应的端电压进行求导便可得到端电容, 然而由于上述电势的表达式十分复杂, 难以推导出解析的导数表达式. 采用差分的方式, 求解端电容的表达式如下:
$ \left\{ \begin{aligned} &{C_{{\text{gg}}}}({V_{\text{G}}}) = \frac{{{\text{d}}{Q_{\text{G}}}}}{{{\text{d}}{V_{\text{G}}}}}\bigg| _{{V_{\text{G}}}} \\ =\,& \frac{{{Q_{\text{G}}}({V_{\text{G}}} + \Delta V) - {Q_{\text{G}}}({V_{\text{G}}} - \Delta V)}}{{2\Delta V}}, \\ & {C_{{\text{gs}}}}({V_{\text{G}}}) = \frac{{{\text{d}}{Q_{\text{S}}}}}{{{\text{d}}{V_{\text{G}}}}}\bigg| _{{V_{\text{G}}}} \\ =\,& \frac{{{Q_{\text{S}}}({V_{\text{G}}} + \Delta V) - {Q_{\text{S}}}({V_{\text{G}}} - \Delta V)}}{{2\Delta V}}, \\ & {C_{{\text{gd}}}}({V_{\text{G}}}) = \frac{{{\text{d}}{Q_{\text{D}}}}}{{{\text{d}}{V_{\text{G}}}}}\bigg| _{{V_{\text{G}}}} \\ =\,& \frac{{{Q_{\text{D}}}({V_{\text{G}}} + \Delta V) - {Q_{\text{D}}}({V_{\text{G}}} - \Delta V)}}{{2\Delta V}},\end{aligned} \right. $
式中Cgg, CgsCgd分别为栅电容、栅源电容和栅漏电容, ΔV为差分步长.
通过TCAD二维数值仿真验证所建立模型的准确性, 为了能够考虑非均匀电场对隧穿过程的影响, TCAD仿真采用非局域带带隧穿模型计算载流子隧穿路径和概率. 此外, 还考虑了掺杂相关的迁移率模型、高电场速度饱和模型以及Shockley-Read-Hall复合等模型.
图2给出了所计算的表面电势分布, 模型计算结果与TCAD结果符合一致, 表明所建立的电势模型的准确性. 图2(a)VD = 1.0 V时不同VG下的电势分布, 可以看到, 随着VG的增大, 沟道电势逐渐升高, 隧穿结附近电势逐渐变得陡峭, 意味着隧穿电场也在逐渐变大, 因此隧穿概率和电流也会随着VG增大而逐步增大. 然而, 当VG逐步增大到1.2 V后, 表面电势的升高开始放缓, 这是由于随着电势的升高, 沟道中感应出了大量可动载流子, 屏蔽了VG对表面电势的影响.
图 2 表面电势随不同的 (a) VG和(b) VD的变化
Figure2. Variation of surface potential with different (a) VG and (b) VD.

虽然沟道中感应出的大量载流子弱化了VG对沟道电势的控制效率, 但是也减小了沟道电阻, 使得漏电场很容易穿透进入沟道, 甚至到达隧穿结附近, 调控隧穿电场及隧穿电流, 如图2(b)所示. 在VG = 1.1 V时, 沟道电势随VG增大而增大, 然而当VD增大到1.0 V以后, 沟道电势几乎不随VD变化而变化, 这是因为随着VD的增大, 沟道中大量载流子被吸引回到漏区, 沟道逐步变为耗尽状态, 沟道电阻增大, 漏电场无法穿透进入沟道, 而是被阻挡在沟道与漏区界面附近, 因此VD失去了对沟道电势的调控能力, 这时沟道电势又回到受VG控制的状态.
由以上分析可以发现, 沟道电势是受到VGVD的交替调控的, 当沟道处于耗尽状态时, 主要受到VG控制; 当沟道处于积累状态时, 主要受到VD控制.
图3是各个端电荷随偏置电压的变化情况.可以看出, QG主要由QD组成, 而QS的影响很小, 几乎可以忽略. 此外, QG随着VG的增大而增大, 随着VD增大而减小, 这是因为栅电场的影响是把电子从漏区吸引到沟道中, 倾向于使得沟道变为积累状态, 而漏电场的影响则是把沟道中的电子吸回到漏区, 倾向于使沟道变为耗尽状态, 这与VGVD对沟道电势的交替调控作用是一致的.
图 3 端电荷随不同的 (a) VG和(b) VD的变化
Figure3. Variation of terminal charges with different (a) VG and (b) VD.

图3(a)中的插图给出了QSVG的变化情况. 可以看出, QS随着VG增大而增大, QS主要取决于源区耗尽宽度LI, 随着VG增大, 隧穿结附近电势升高, 源耗尽宽度展宽, 因此QS变大. 然而, 图3(b)插图中QSVD的变化情况与VG有所不同, 随着VD增大, QS先是逐渐增大, 而后保持不变, 这是由于沟道中载流子浓度较高时, 沟道电阻较小, 漏电场可以穿透沟道, 到达隧穿结附近, 调控隧穿结附近电势及源耗尽宽度LI, 但是随着VD增大, 沟道载流子浓度降低, 电阻变大, 漏电场失去了对隧穿结附近电势的控制, 因此LIQS也不受VD调控.
模型计算的端电荷与TCAD结果取得了较好的一致性, QS的计算误差较大是因为在计算I区内耗尽电荷时, 简单地把I区看作一个圆柱体, 也就是在图1(b)中将I区作为一个矩形处理, 实际上φS(0)在整个源/沟道界面上是变化的, 因此I区的耗尽宽度在r方向上并非常数, 也就是图1(b)中I区的形状并不是严格的矩形. 矩形假设以及在I区电势推导中采用的泰勒级数展开所引入的误差, 最终导致QS计算偏差较大, 但是所建立的模型从基本的器件物理出发, 仍然能够正确地解释QS随偏置电压的变化趋势和规律.
图4给出了模型和TCAD所计算的电容特性. 可以看出, 模型所计算的电容与TCAD结果取得了很好的一致性, 能够准确描述电容变化的基本规律. CggVG增大而增大, 随着VD增大而减小, 与图3中端电荷的变化规律一致. 此外, 还能够看到, 栅电容Cgg主要由栅漏电容Cgd组成, 栅源电容Cgs的贡献几乎可以忽略, 这就导致GAA-TFET具有较大的米勒电容. 在数字电路中, 米勒电容效应将引起大的信号过冲以及延迟, 不利于电路特性的提升, 这一效应的影响可以通过采用异质栅结构进行缓解.
图 4 端电容随不同的 (a) VG和(b) VD的变化
Figure4. Variation of terminal capacitances with different (a) VG and (b) VD.

另外, 需要注意Si 材料禁带宽度大, 且属于间接带隙半导体, 载流子隧穿需要额外声子辅助, 载流子隧穿概率和隧穿电流十分低, 难以满足器件实际应用需求[3]. 相比之下, InAs材料禁带宽度窄, 属于直接带隙, 利于提高载流子隧穿概率, 提高TFET器件特性[16]. 从以上分析可以看出, 虽然本文主要研究了Si基GAA-TFET器件, 但是模型推导并非针对某一特定材料, 因此, 本文所提出的模型具有一定程度的材料普适性, 能够适用于更适合TFET设计的InAs等III-V族材料. 图5给出了模型所计算的InAs GAA-TFET器件表面电势及端电容特性, 可以看到模型能够精准预测器件表面电势变化, 准确描述器件电容基本变化规律.
图 5 InAs基GAA-TFET器件(a)表面电势和(b)端电容随VG的变化
Figure5. Variation of the (a) surface potential and (b) terminal capacitance with different VG for an InAs based GAA-TFET.

本文从GAA-TFET的基本器件物理出发, 建立了其表面电势和电容模型, 所建立模型计算结果与TCAD数值计算结果取得了很好的一致性, 能够准确描述GAA-TFET中栅、漏电压对表面电势的交替调控原理以及端电容的变化规律, 表明所建立模型的准确性. 此外, 所建立的模型不涉及任何数值迭代过程, 相比于数值模型, 该模型计算速度快, 计算过程稳定, 能够用于GAA-TFET器件及电路的相关研究.
相关话题/计算 电路 纳米 电压 结构

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于MXene涂层保护Cs<sub>3</sub>Sb异质结光阴极材料的计算筛选
    摘要:以锑化铯(Cs3Sb)为代表的碱金属型半导体光阴极具有高量子效率、低电子发射度、光谱响应快等特点,可作为理想的新型电子发射源.然而Cs3Sb中碱金属敏感于含氧气体,从而导致其结构不稳定,工作寿命低,影响电子发射效率.利用超薄层状的二维材料进行涂层保护Cs3Sb基底,有望构建新型高性能光阴极材料 ...
    本站小编 Free考研考试 2021-12-29
  • 渔网超结构的等离激元模式及其对薄膜电池的陷光调控
    摘要:渔网超结构具有平面、近光学无损、特定光场中可以激发表面等离激元等特点,在增强光子器件的响应效率方面极具潜力.本文基于时域有限差分方法和严格耦合波分析,系统研究了渔网超结构的等离共振模式及其对晶硅薄膜电池的光波调控性能.研究结果表明,渔网结构对光波的吸收、散射和消光特性强烈依赖金属层的厚度、线宽 ...
    本站小编 Free考研考试 2021-12-29
  • 单缺陷对Sc, Ti, V修饰石墨烯的结构及储氢性能的影响
    摘要:寻找稳定高效的储氢材料是实现氢经济的关键.过渡金属修饰石墨烯储氢材料在理论上被广泛研究,但存在H2解离和金属团聚的问题.本文基于密度泛函理论对Sc,Ti,V修饰单缺陷石墨烯的结构及储氢性能进行计算.结果表明:单缺陷使Sc,Ti,V与石墨烯的结合能提高4—5倍;Sc,Ti,V离子特性增强,可以通 ...
    本站小编 Free考研考试 2021-12-29
  • 二维共价键子结构Zintl相热电材料研究及进展
    摘要:热电材料可以实现热能和电能间的直接相互转换,在半导体制冷和热能回收方面有着重要应用.Zintl相热电材料由电负性差异较大的阴阳离子组成,其输运特征符合“声子玻璃,电子晶体”的概念,因此受到了广泛的研究,特别是具有二维共价键子结构Zintl相热电材料凭借优异的电性能更是被寄予厚望.本文综述了具有 ...
    本站小编 Free考研考试 2021-12-29
  • CuBiI三元化合物晶体结构预测及光电性能第一性原理研究
    摘要:作为潜在的新型光电材料,三元金属卤化物一直以来广受关注.本文通过基于遗传算法的晶体结构预测软件USPEX,对三元CuBiI化合物(CuBi2I7,Cu2BiI5,Cu2BiI7,Cu3BiI6,Cu3Bi2I9,CuBi3I10,Cu4BiI7)在常压、绝对零度下的稳定晶体结构进行了全局搜索. ...
    本站小编 Free考研考试 2021-12-29
  • 二维双金属铁磁半导体CrMoI<sub>6</sub>的电子结构与稳定性
    摘要:二维磁性半导体由于兼具磁性、半导体性和特殊的二维结构而受到人们的广泛关注,为纳米级自旋电子和光电子器件的研发应用和相关的基础理论研究提供了新的思路和平台.基于第一性原理计算,在对一系列二维双金属碘化物CrTMI6的交换能进行初步筛选的基础上,选出了具有铁磁性的CrMoI6单层结构.进一步计算表 ...
    本站小编 Free考研考试 2021-12-29
  • 混合时钟驱动的自旋神经元器件激活特性和计算性能
    摘要:自旋神经元是一种新兴的人工神经形态器件,其具有超低功耗、强非线性、高集成度和存算一体等优点,是构建新一代神经网络的强有力候选者.本文提出了一种磁场辅助磁弹应变驱动的混合时钟自旋神经元,利用OOMMF微磁学仿真软件建立了该神经元器件的微磁学模型,基于LLG方程建立了其数值仿真模型,利用所设计的自 ...
    本站小编 Free考研考试 2021-12-29
  • 磁制冷材料LaFe<sub>11.5</sub>Si<sub>1.5</sub>基合金成分与磁相变温度关系的高通量计算
    摘要:获得具有不同磁相变温度的La(Fe,Si)13基合金对拓宽磁制冷工作温区具有重要意义.借助第一性原理模拟软件AMS-BAND模块并结合平均场理论对LaFe11.5Si1.5基磁制冷合金的磁相变温度进行了高通量计算.研究了Mn,Co,Ni,Al和Fe缺位掺杂对LaFe11.5Si1.5基合金体系 ...
    本站小编 Free考研考试 2021-12-29
  • 基于时域有限差分法的核壳双金属纳米颗粒光吸收率反转行为
    摘要:双金属纳米颗粒能够有效整合两种金属的物理和化学性质并同时表达每种金属的独特性质,是提高光散射、光热转换、等离激元共振衰变和光子激发的重要材料.基于单独纳米颗粒的研究可以避免实验研究过程中纳米颗粒之间的相互影响,更能够有效分析入射光与纳米颗粒之间的相互作用.本文采用时域有限差分法计算了等离激元双 ...
    本站小编 Free考研考试 2021-12-29
  • 金纳米粒子修饰氧化铟锡阳极的高效率红光钙钛矿发光二极管
    摘要:金纳米颗粒对提升钙钛矿发光二极管的外量子效率有重要作用.为了避免金纳米颗粒与发光层直接接触,先前工作合成的AuNPs@SiO2影响器件电荷传输且不易合成;而将金纳米颗粒共混在聚(3,4-乙烯二氧噻吩):聚(苯乙烯磺酸酯)中时,金纳米密度又不易控制,不适合做理想的空穴传输层.于是,本文采用静电吸 ...
    本站小编 Free考研考试 2021-12-29