太赫茲互補金屬氧化物半導體場效應管探測器理論模型中擴散效應研究?

張鏡水 孔令琴 董立泉 劉明 左劍 張存林 趙躍進?

1)(北京理工大學光電學院,北京 100081)

2)(首都師范大學物理系,北京 100048)

太赫茲互補金屬氧化物半導體場效應管探測器理論模型中擴散效應研究?

張鏡水1)孔令琴1)董立泉1)劉明1)左劍2)張存林2)趙躍進1)?

1)(北京理工大學光電學院,北京 100081)

2)(首都師范大學物理系,北京 100048)

(2017年3月1日收到;2017年4月5日收到修改稿)

針對基于經典動力學理論傳統模型中忽略擴散效應的問題,通過對基于玻爾茲曼理論的場效應管傳輸線模型的理論分析,建立了包含擴散效應的太赫茲互補金屬氧化物半導體(CMOS)場效應管探測器理論模型,研究擴散效應對場效應管電導及響應度的影響.同時,將此模型與忽略了擴散效應的傳統模型進行了對比仿真模擬,給出了兩種模型下的電流響應度隨溫度及頻率變化的差別.依據仿真結果,并結合3σ原則明確了場效應管傳輸線模型中擴散部分省略的依據和條件.研究結果表明:擴散部分引起的響應度差異大小主要由場效應管的工作溫度及工作頻率決定.其中工作頻率起主要作用,溫度變化對差異大小影響較為微弱;而對于工作頻率而言,當場效應管工作頻率小于1 THz時,模型中的擴散部分可以忽略不計;而當工作頻率大于1 THz時,擴散部分不可省略,此時場效應管模型需同時包含漂移、散射及擴散三個物理過程.本文的研究結果為太赫茲CMOS場效應管理論模型的精確建立及模擬提供了理論支持.

互補金屬氧化物半導體場效應管探測器,太赫茲,模型模擬,玻爾茲曼理論

1 引 言

隨著太赫茲波技術在軍事和民用各個領域的不斷發展,研制高響應度、低噪聲的太赫茲探測器成為相關科研領域的研究重點之一,而互補金屬氧化物半導體(CMOS)場效應管太赫茲探測器因其成熟的技術、相對較低的成本及其他材料目前無法企及的良好性能(其最高響應度可達2.2 kV/W,等效噪聲功率可達到而受到廣泛的關注[1?5].為了減小探測器的制作成本,需要對設計好的探測器進行性能模擬測試,而該項工作的可信度依賴于太赫茲場效應管模型的精準程度,故精確的太赫茲場效應管探測器模型一直是國內外相關領域的重點研究方向[6?14].

目前常用太赫茲波段CMOS場效應模型為基于均勻電荷控制理論的線性電阻-電容-電感(resistance-capacitance-inductance,RCL)傳輸線理論模型.該模型運用1996年Dyaknov和Shur[15]提出的均勻電荷控制理論,通過經典動力學理論的一維德魯德模型描述溝道內載流子的整體運動,并將其轉變為由許多RCL單元組成的RCL傳輸線模型[16].目前針對基于該模型的理論研究只考慮載流子的漂移及散射運動,忽略了溝道內載流子濃度不均勻分布的擴散效應的影響.此外,關于擴散效應可省略的依據及條件并未見相關報道.

基于此,本文從玻爾茲曼理論出發,研究模型中擴散效應對太赫茲場效應管電導及響應度的影響.首先從玻爾茲曼電流方程推導出基于此理論的RCL模型方程,并通過與一維德魯德RCL模型的對比,明確擴散效應在模型中的表達式.其次通過理論計算與模擬分析,分別得到靜態與動態下兩種不同理論的場效應管模型的響應區別,分析擴散效應對場效應管模型響應度的影響.同時,通過分析擴散效應隨溫度與頻率的變化情況,結合3σ原則給出了擴散效應在模型中可省略的條件.研究結果為太赫茲CMOS場效應管探測器模型的精確建立提供了重要依據,并為其進一步技術及工程實現提供了關鍵理論及技術支持.

2 模型理論

現在的模型主要基于流體動力學中一維德魯德模型,其主要包含漂移和散射效應.而玻爾茲曼方程則主要從非平衡態時的分布函數出發,推導出溝道的電流情況.故本節中通過對比二者之間的差異來說明擴散效應的重要性.

2.1 經典動力學傳輸線模型的基本形式

根據經典動能理論,溝道中的載流子運動可以用德魯德模型來描述.基于該理論,對于溝道中電荷為q的電子,其動量p(x,t)為

其中E(x,t)為溝道x處的電場大小;τ為有效電子散射時間,此處考慮聲子作用與雜質碰撞.為了描述場效應管溝道中載流子的運動情況,Dyaknov和Shur[15]用經典動力學理論提出了一維德魯德模型及粒子連續邊界條件,

式中m為粒子的等效質量;n,v和Ug分別為溝道x處的載流子密度、載流子速度和柵極至溝道處的電壓幅值大小.

若只考慮溝道內等離子體的基波與一次諧波,(2a)與(2b)式可以寫成:

其中W為場效應管溝道寬度,μ為粒子遷移率,I是溝道內電流.(3a)和(3b)式與電報方程的基本形式一致,故溝道中載流子的運動實際可表示為RCL傳輸線模型,如圖1所示.其R,L和C的阻抗表達分別為[16]

(4a),(4b)及(4c)式即為基于經典動力學理論的場效應管RCL傳輸線模型方程.運用該模型方程,場效應管溝道中載流子的運動可轉化為如圖1所示的模型.

圖1 場效應管溝道RCL傳輸線理論模型,其中Ug為柵極至溝道的電壓差,Ugs為柵源極間所加直流電壓,Uac為通過天線耦合進場效應管的太赫茲波;R,L和C的計算分別如(4a),(4b)和(4c)式所示Fig.1.The transmission line model for fi eld-electrictransistor in high frequency.Ugis the gate-to-channel local voltage,Ugsis the gate-to-source DC voltage,and Uacis the terahertz signal that coupled by the antenna.R,C and L are calculated with Eq.(4a),(4b)and(4c).

2.2 基于玻爾茲曼理論傳輸線模型的基本形式

玻爾茲曼傳輸方程描述了非平衡下的熱動力學系統中流體特征的變化情況.從玻爾茲曼傳輸方程可求得溝道中一維電流密度J滿足如下公式[6]:

式中kB為玻爾茲曼常數,Tn為溝道內載流子溫度.其中在現代模型模擬中通常省略,省略后可得到只包含亞聲波流動區域的拋物線能量波矢關系曲線系統[6].同時,設定沿溝道方向的溫度變化為零,給定熱電壓VT的定義為VT=kBTn/q,則(5)式可表達為

為進一步明確玻爾茲曼理論與經典動力學理論的差別,將(3a)式轉化為電流密度形式,可得基于經典動力學理論的溝道一維電流密度為

對比(6)式和(7)式,可得玻爾茲曼理論除包含經典動力學理論中的漂移部分qnμE和散射部分外,還包含擴散部分

將(7)式轉換成傳輸線理論方程可得

結合粒子連續邊界條件,可得基于玻爾茲曼理論的RCL傳輸線模型表達式為

為方便對比兩種模型,此處用電導取代電阻,表述RCL模型單元中載流子主要的漂移運動.由(4a)和(9a)式可得,基于流體力學德魯德模型的總電導(每單位長度)如(10a)式,而基于玻爾茲曼傳輸方程的總電導(每單位長度)如(10b)式:

由(10b)式可得,基于玻爾茲曼理論的總電導由兩部分組成,第一部分為漂移電導,其形式與流體力學德魯德模型所得總電導相同;第二部分即為擴散電導,與溫度及載流子密度隨Ug的變化有關.

3 模型計算及仿真結果

為了更好地說明經典動力學傳輸線模型和玻爾茲曼理論傳輸線模型的差別,本節中根據前文理論計算,分別對兩種模型進行靜態下電導仿真和高頻電流響應度仿真,并通過二者差值對兩種模型的差異性進行闡述分析.

3.1 靜態下擴散部分對整體電導的影響

根據(10a)及(10b)式,可得漂移電導和擴散電導的計算式分別為

為探討擴散電導的影響,我們用擴散電導σDiffusion與傳統基于經典動力學理論總電導σKinetic的比值Raσ來表示擴散電導對整體電導的影響.Raσ的計算如(12)式所示:

根據(12)式可計算得到擴散部分的加入對整體電導的影響.該計算中相關取值依據如下:場效應管閾值電壓取為CMOS場效應管常用閾值電壓Vth=0.4 V,柵極電壓變化范圍為CMOS場效應管常用柵極電壓范圍0—1 V,溫度變化范圍為Si材料正常工作溫度273 K(0?C)至523 K(250?C),由此得到計算結果如圖2所示.由圖2(a)可知,擴散部分電導對整體電導的影響隨柵極電壓的增大而減小,隨溫度緩慢遞增.這是由擴散運動本質,即粒子濃度不均勻性引起的.隨著柵極電壓的增大,溝道內載流子濃度隨柵極電壓的變化逐漸降低,粒子濃度不均勻性減小,擴散效應減弱,因此擴散電導對整體電導的影響力隨之減小.當溫度增大時,載流子無規則熱運動增大,擴散效應隨之變強.

然而在圖2(a)中,溫度的影響不明顯.為更好地體現溫度變化在擴散運動中的作用,此處取柵極電壓分別為:Ugs=0.2 V,Ugs=0.4 V及Ugs=0.6 V,可得到比值Raσ在此三個柵極電壓處隨溫度的變化如圖2(b)所示.由圖2(b)可知,比值Raσ隨溫度的增大緩慢增大,且隨柵極電壓的增大變化趨勢逐漸減緩,這也與之前所得結論——擴散部分電導對整體電導的影響隨柵極電壓的增大而減小相符.

另外,一般場效應管探測器工作在室溫(T=300 K),其最佳工作點在閾值電壓附近(Ugs=Vth=0.4 V),由圖2可以看到,此時二者之間的差值約占經典動力學模型電導的25%,仍為較大比重,故此時擴散效應不可忽略.

圖2 (網刊彩色)(a)Raσ隨柵極電壓及溫度的變化情況;(b)不同柵極電壓下Raσ隨溫度的變化情況Fig.2.(color online)(a)The variation of Raσwith the gate voltage and the temperature;(b)the variation of Raσwith the temperature at di ff erent gate voltage.

3.2 擴散電導對場效應管電流響應度的影響

電流響應度是評價場效應管探測器的重要參數之一.因探測器實際工作時需配合天線完成自由空間太赫茲波耦合,故電流響應度的高低除受場效應管本身性能決定之外,還與配合天線的損耗與二者之間的匹配系數有關.此處假定天線無損耗,且天線與場效應管完全匹配,分別對基于玻爾茲曼理論模型的電流響應度?Boltzmann與基于經典動力學德魯德理論模型的電流響應度?Kinetic進行仿真分析,并為了研究二者的差別,給出其差值定義RRI為

所得到的模擬計算結果如圖3所示.由圖3可知,二者差值隨溫度升高而逐漸降低,隨頻率增大而逐漸增大.對于CMOS場效應管探測器而言,其正常工作在室溫(T=300 K)下,由圖3可知,二者差值在1 THz以下差距較小,當探測器的工作頻率大于1 THz時,二者的差距開始逐漸增大.

圖3 (網刊彩色)RRI隨柵極頻率及溫度的變化情況,其中,橫軸為溫度,縱軸為場效應管工作頻率Fig.3.(color online)The variation of RRIwith the frequency and the temperature.In the fi gure,the x axis is temperature,and the y axis is the transistor working frequency.

4 仿真結果分析

由上文可知,當頻率和溫度不同時,RRI不同.為了給出明確的擴散部分省略條件,我們依據如下兩條判斷依據給出省略原則.

1)數學計算中,認為當滿足

時,δ→0.

此處我們運用3σ原則對“→”運算進行判斷,認為當滿足

時,(1+δ)(1?δ)→1,δ→0.

2)給定命題真假評判系數α,并設定需評判的命題為模型中擴散部分可以省略.可得當擴散部分可省略成立時,命題為真,α=1;當模型中擴散部分不可省略時,命題為假,α=0.

由以上兩條原則,結合(13)式與(15)式,可得當

由(17)式得到α隨溫度及頻率的變化如圖4所示.其中藍色部分為α=1,即擴散部分可省略區域;紅色部分為α=0,即擴散部分不可省略區域.圖4中紅色占據主要部分,故可知在大部分情況下,模型中的擴散部分是不可以忽略的.只有當頻率足夠低(頻率<1 THz)、部分頻率下溫度足夠高時,場效應管模型可簡化為只包含漂移和散射部分,忽略擴散效應對模型的影響.

圖4 (網刊彩色)α隨頻率及溫度的變化情況,圖中橫軸為溫度,縱軸為場效應管工作頻率;圖中藍色部分為擴散部分可以省略區域,紅色部分為不可省略區域Fig.4.(color online)The variation of α with the frequency and the temperature.In the fi gure,the x axis is temperature,and the y axis is the transistor working frequency.The blue part is the area that the di ff usion part can be neglected in the model,and the red one is for not neglected part.

由圖4可知,場效應管探測器工作頻率越大,擴散效應產生的影響越明顯.產生此現象的原因與溝道內載流子濃度不均勻性隨頻率增大而提高有關.當場效應管探測器工作在太赫茲頻段并且溝道長度大于溝道內等離子體波的衰減長度時(衰減長度隨等離子體波振蕩頻率增大而降低,太赫茲波段衰減長度一般為10—30 nm,小于常用CMOS場效應管溝道長度),太赫茲波從場效應管溝道一段耦合進溝道內,并在傳輸過程中等離子體波振幅迅速衰減至零.這就導致了溝道一段載流子濃度隨太赫茲波振幅改變而變化,頻率越高,太赫茲波振幅變化越快,載流子濃度變化越大,溝道兩端的載流子濃度不均勻性也隨之提高,進而令擴散效應愈發明顯.因此,在高頻區域,模型中的擴散運動不可忽略.

5 結 論

綜上所述,本文針對現有太赫茲CMOS場效應管理論模型中忽略擴散效應的問題,通過理論計算及模擬分析給出模型中擴散部分的作用,明確了該部分可以省略的理論依據和條件.文中首先給出了基于玻爾茲曼理論的場效應管傳輸線模型,并將其與只包含漂移與散射效應的經典動力學理論場效應管傳輸線模型進行對比,從理論上得到二者之間的差別和擴散部分電導的具體表達式,并明確溝道中載流子擴散效應主要由溫度及載流子密度隨柵極至溝道的電壓引起.其次,本文通過模擬結果給出了靜態下的電導比及高頻下的電流響應度之比來討論兩種模型之間的差異.結果表明:靜態下,擴散部分的影響隨溫度的升高呈微弱增大,隨柵極電壓的增大而減小;常溫下,當柵極電壓為閾值電壓時,場效應管探測器處于最佳工作點,此時擴散電導占總電導的25%;高頻狀態下,擴散部分對探測器電流響應度的影響隨頻率增大而增大,隨溫度升高而降低.同時,基于計算及模擬結果,我們結合3σ原則提出擴散部分可省略的具體條件:當頻率足夠低(頻率<1 THz)或溫度足夠高時,模型中的擴散部分可以省略,此時兩個模型可認為無差別.當場效應管工作頻率高于1 THz時,模型中擴散部分的省略會導致場效應管模型的不準確.產生此現象的原因為擴散運動與溝道內載流子濃度不均勻性有關.柵極電壓增大時,溝道內載流子濃度隨柵極電壓的變化逐漸降低,粒子濃度不均勻性減小,擴散效應減弱,因此擴散電導對整體電導的影響力隨之減小;當溫度增大時,載流子無規則熱運動增大,擴散效應隨之變強;當頻率增高時,溝道一端等離子體波的振幅變化越快,而溝道另一端因衰減效應始終為零,故此時溝道兩端載流子濃度變化增大,載流子濃度不均勻性變化也隨之提高,進而令擴散效應愈發明顯.綜上所述,場效應管正常工作時的柵極電壓一般在閾值電壓附近,此時隨著溫度增高、工作頻率增大,溝道內的擴散效應影響增強,在模型中不可忽略.

本文的研究工作為CMOS場效應管模型的精確建立提供了理論依據,為更好地進行場效應管的設計工作和節約成本提供了更有效的保障.

[1]Pfei ff er U R,Grzyb J,Sherry H,Cathelin A,Kaiser A 2013 38th International Conference on Infrared,Millimeter,and Terahertz Waves(IRMMW-THz)Mainz,Germany,September 1–6,2013 p1

[2]Minoru F,Shuhei A 2015 IEICE Electron.Express 12 20152006

[3]Lu J Q,Shur M S,Hesler J L 1998 Electron Dev.Lett.19 373

[4]Hadira R A,Sherry H,Grzyb J,Zhao Y 2012 IEEE J.Solid-State Circuit 47 2999

[5]Ryu M W,Lee J S,Kim K S,Park K,Yang J R,Han S T,Kim K R 2016 IEEE Trans.Electron Dev.63 1742

[6]Grasser T,Tang T,Kosina H,Selberherr S 2003 Proc.IEEE 91 251

[7]Preu S,Kim S,Verma R,Burke P G,Sherwin M S,Gossard A C 2012 J.Appl.Phys.111 024502

[8]Gutin A,Nahar S,Hella M,Shur M 2013 IEEE Trans.Terahertz Sci.Technol.3 545

[9]Ibrahim N Y,Rafat N H,Elnahwy S E A 2013 J.Infrared Millim.Terahertz Waves 34 606

[10]Tan R B,Qin H,Sun J D,Zhang X Y,Zhang B S 2013 Appl.Phys.Lett.103 173507

[11]Zhao X H,Li C,Zhang P 2013 Acta Phys.Sin.62 130506(in Chinese)[趙曉輝,蔡理,張鵬 2013物理學報62 130506]

[12]Gutin A,Ytterdal T,Muraviev A,Shur M 2015 Solid-State Electron.104 75

[13]Kim K S 2016 M.S.Thesis(Ulsan:Ulsan National Institute of Science and Technology)

[14]Liu Y,He J,Chan M S,Du C X,Ye Y,Zhao W,Wu W,Deng W L,Wang W P 2014 Chin.Phys.B 23 097102

[15]Dyakonov M I,Shur M S 1996 IEEE Trans.Electron Dev.43 1640

[16]Khmyrova I,Seijyou Y 2007 Appl.Phys.Lett.91 143515

PACS:73.43.Cd,73.50.Mx,85.30.Tv,85.60.GzDOI:10.7498/aps.66.127302

Di ff usion part in terahertz complementary metal oxide semiconductor transistor detector model?

Zhang Jing-Shui1)Kong Ling-Qin1)Dong Li-Quan1)Liu Ming1)Zuo Jian2)Zhang Cun-Lin2)Zhao Yue-Jin1)?

1)(Institute of Opto-Electronics,Beijing Institute of Technology,Beijing 100081,China)
2)(Department of Physics,Capital Normal University,Beijing 100048,China)

1 March 2017;revised manuscript

5 April 2017)

In this paper,we discuss the di ff usion motion of carriers in the transistor channel in a terahertz frequency range,and propose an resistance-capacitance-inductance(RCL)model based on Boltzmann transport theory,and then put forward the rules to determine whether the di ff usion part in the RCL model can be neglected for terahertz fi eld-e ff ect-transistor(FET)detectors.The traditional RCL model for FET detectors is based on classic kinetic theory.In this model only the drift and the scattering motion of the carrier density in transistor channel are considered,and the di ff usion part is neglected without giving any explanation.To solve this problem,in this paper we adopt three steps: fi rst,instead of classic kinetic theory,the equations of RCL transistor model including di ff usion part are derived from Boltzmann transport equation,and by comparing the two models,the speci fi c expression for the di ff usion part is given.Second,the di ff erences between the two models are calculated and simulated,including the conductivity in quasi-static mode and the current response in high frequency mode,with di ff erent gate voltages,temperatures and working frequencies.Third,combined with the 3σ rules,the conditions to neglect the di ff usion motion in the model are put forward.The results show that the di ff usion motion of the carriers is caused by the inhomogeneity of the carrier density,a ff ected by the gate voltage,the temperature and the changing speed of the carriers with respect to the local voltage.In quasi-static mode,the role of di ff usion part will change with the gate voltage,and when the gate voltage equals threshold voltage(which is the best working point for transistor detector),the di ff usion part cannot be neglected,for which the reason is that a larger gate voltage will lead to a smaller inhomogeneity of channel carrier density and then a weaker di ff usion e ff ect,thus the e ff ect of di ff usion conductance on the whole transistor conductance becomes smaller.For the terahertz-frequency working mode,the di ff usion part will depend on temperature and frequency.With temperature increasing,the current responsivity di ff erence caused by the di ff usion part in the model slightly decreases;when the working frequency increases but below 1 THz,the di ff usion part can be neglected;however,when the working frequency is above 1 THz,the transistor model should contain drift,scattering and di ff usion part at the same time,for which the explanation is that when the temperature increases,the random thermal motion of the carrier becomes larger,thus the di ff usion e ff ect will be stronger;and if the frequency increases,the number of the carriers in one terminal of the channel will change faster,but due to the channel damping,the number of the carriers in another terminal will always be zero,thus the changing speed of the carrier density between the two terminals will be faster,then a larger inhomogeneity of carrier density and a stronger di ff usion e ff ect will appear.In conclusion,normally the transisitor works at the threshold gate voltage,and at this point,the di ff usion e ff ect in the channel will increase with working temperature and frequency increasing,thus the di ff usion part in the model cannot be neglected.The results in this paper make a signi fi cant contribution to a more accurate terahertz transistor detector model.

complementary metal oxide semiconductor transistor detector,terahertz,model simulation,Boltzmann theory

10.7498/aps.66.127302

?國家自然科學基金(批準號:61377109)資助的課題.

?通信作者.E-mail:yjzhao@bit.edu.cn

?2017中國物理學會Chinese Physical Society

http://wulixb.iphy.ac.cn

*Project supported by the National Natural Science Foundation of China(Grant No.61377109).

?Corresponding author.E-mail:yjzhao@bit.edu.cn