CN110008489B - 一种THz频段InP DHBT器件在片测试结构建模方法 - Google Patents

一种THz频段InP DHBT器件在片测试结构建模方法 Download PDF

Info

Publication number
CN110008489B
CN110008489B CN201811405567.1A CN201811405567A CN110008489B CN 110008489 B CN110008489 B CN 110008489B CN 201811405567 A CN201811405567 A CN 201811405567A CN 110008489 B CN110008489 B CN 110008489B
Authority
CN
China
Prior art keywords
formula
short
frequency
representing
resistance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811405567.1A
Other languages
English (en)
Other versions
CN110008489A (zh
Inventor
刘军
刘杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
CETC Information Science Research Institute
Original Assignee
Hangzhou Dianzi University
CETC Information Science Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University, CETC Information Science Research Institute filed Critical Hangzhou Dianzi University
Priority to CN201811405567.1A priority Critical patent/CN110008489B/zh
Publication of CN110008489A publication Critical patent/CN110008489A/zh
Application granted granted Critical
Publication of CN110008489B publication Critical patent/CN110008489B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/39Circuit design at the physical level
    • G06F30/398Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种THz频段InP DHBT器件在片测试结构建模方法,模型拓扑结构基于测试物理结构建立,并对其在亚毫米波段的高频寄生进行相对完整的考虑。模型的容性和阻性寄生采用解析提取技术,从开路结构低频测试数据中获取。模型的高频趋肤效应采用传统物理公式计算初值,并结合短路测试结构的低频解析提取结果对计算公式进行修正,使其适用于实际测试结构建模。

Description

一种THz频段InP DHBT器件在片测试结构建模方法
技术领域
本发明涉及集成电路领域,更具体的说,它涉及一种THz频段InP DHBT器件在片测试结构建模方法。
背景技术
与GaAs系HBT相比,InP系HBT由于材料系统固有的特性,使其在直流增益、高频性能、功耗和1/f噪声等诸多方面体现出优势。经过几十年的发展,作为InP HBT器件的典型,II型InP DHBT在应用频段上逐渐逼近THz频段,成为高速、低功耗毫米波、亚毫米波和超高速器件及相关集成电路技术领域的研究热点。
基于G-S-G探针接触的在片测试是获取晶体管高频行为的有效手段。为去除测试结构引入的高频寄生的影响,通常需要设计出包括Open、Short和TRL等不同的测试结构,支持不同频段、不同去嵌入算法的执行。测试结构的精确建模,也是准确评估去嵌入算法的正确性、算法执行的精度以及新的去嵌入算法推导的基础。以等效电路形式建立G-S-GPAD以及Open、Short等结构模型的工作,已有不少文献见于报道。但无论是RF CMOS工艺还是化合物工艺,已开展工作多集中在110GHz频段以下,少量工作达到220GHz频段。对建立的模型,又多采用解析提取技术获取模型参数值。当模型拓扑结构愈来愈复杂,解析提取算法的推导工作变得极其困难,且难以验证、难以随测试结构尺寸变化进行可缩放应用,需要提出新的参数提取算法以适用于不同模型的研究。
发明内容
本发明克服了现有技术的不足,提出一种THz频段InP DHBT器件在片测试结构建模方法,用于解决目前模型拓扑结构不准确及其提参方法不精准的问题。
本发明的技术方案如下:
一种THz频段InP DHBT器件在片测试结构建模方法,具体包括如下步骤:
101)建立模型步骤:建立短路结构和开路结构,其中短路结构相对于开路结构增加了短路连接线,短路结构和开路结构都包括测试地平面的寄生电感和测试地平面的寄生电阻;
102)低频下参数获取步骤:低频状况下,地平面的寄生电感、地平面的寄生电阻不起作用,金属连接线理想接地,其等效的拓扑结构包括Cmx、Cmi、Csl、Ctl、Ctr、Csr、Rsl、Rsr表示G-S-G PAD和金属传输线与地平面之间的总寄生电容、寄生电阻;Rsl为G-S-G PAD和金属传输线之间的寄生电阻,Rsr为G-S-G PAD和地平面之间的寄生电阻;其中Ctl=Cstl//Csdl,Ctr=Cstr//Csdr
将等效的拓扑结构转换为π型网络,该π型网络的各部分参数Y关系如下:Yl表示Rsl、Csl和Ctl网络,Yr表示Rsr、Csr和Ctr网络,Ym表示Rmx、Cmx和Cmi网络,可依次推导为Ym=-Y21=-Y12,Yl=Y11+Y12,Yr=Y22+Y12
其中Ym还可根据等效拓扑结构表示为如下公式:
Ym=jωCmi+[Rmx+(jωCmx)-1]-1 公式(1)
其中,Ym的实部和虚部表达式如下:
[Re(Ym)]-1=Rmx+(ω2Cmx 2Rmx)-1 (2)
Im(Ym)=[Cmi+Cmx(1+ω2Cmx 2Rmx 2)-1]×ω (3)
其中根据实际检测数据,可得[Re(Ym)]-1关于ω-2的截距和斜率,从而依次得到Rmx和Cmx,将提取所得Rmx和Cmx代入公式(3),可得到Cmi的值;
Cmi=ω-1Im(Ym)-Cmx(1+ω2Rmx 2Cmx 2)-1 公式(4)
同理Yl和Yr与Ym有相同的拓扑结构,提取方法同上;从而同理得到Ctl和Ctr,Ctl和Ctr参数提取后,分别由Ctl和Ctr计算Cstl、Csdl以及Cstr、Csdr的值时,引入了比例因子Kl和Kr,计算公式分别为Csdl=KlCtl、Cstl=(1-Kl)Ctl以及Csdr=KrCtr、Cstr=(1-Kr)Cr,其中Kl和Kr的值在实际检测中确定,其位常数;
根据已知Signal-PAD的半径和Signal-PAD到待测器件连接线的尺寸以及地平面的长宽,得到其有效面积;由于寄生电容与金属结构面积成正比,寄生电阻与面积成反比,从低频提取所得总的电阻和电容分量中,根据面积比例可分离算出Cstfl、Csfl、Rsfl和Cstgl、Csgl、Rsgl,(Cstfl、Csfl、Rsfl和Cstgl、Csgl、Rsgl表示G-S-G PAD到地平面的寄生电容和电阻)并得到其值;根据电阻和电容的单位面积因子,得到在高频下考虑的地平面的寄生电容和电阻Czi、Czx和Rzx的值,其中Czx、Czi和Rzx表示地平面的高频寄生电容和电阻,对于Czl和Czr采用以下方法计算Czl=2×Cmi,Czl=2×Cmi;其中Czl和Czr是信号端到地平面的寄生电容;
103)高频下参数获取步骤:高频状况下,金属感性寄生电阻和金属感性寄生电容,其与参数Z之间的关系如下:
ZShort′=(YShort-YOpen)-1 公式(5)
YShort表示采用Short结构的容性寄生,YOpen表示Open结构的容性寄生,ZShort表示Short结构的参数Z;
高频状况下,金属连接线和地平面的感性寄生电阻和感性寄生电容变得显著,在不考虑各部分的趋肤效应的情况下,其等效电路模型包括Za,Zb和Zm,其中通过参数Z可分别确定如下:
Za=[ZShort′]11-[ZShort′]12 公式(6)
Zb=[ZShort′]22-[ZShort′]12 公式(7)
Zm=[ZShort′]12=[ZShort′]21 公式(8)
Za,Zb和Zm还可进一步表示为:
Za=(Rl+Rstubl)+jω(Ll+Lstubl) 公式(9)
Zb=(Rr+Rstubr)+jω(Lr+Lstubr) 公式(10)
Zm=(Rm+Rstubm)+jω(Lm+Lstubm) 公式(11)
Rl、Ll表示金属连接线的感性寄生电阻和电感,Rm、Lm表示地平面的感性寄生电阻和电感,Rstubl、Rstubr、Rstubm和Lstubl、Lstubr、Lstubm表示短路连接线的寄生电阻和电感;
其中Rstubl和Lstubl的提取公式如下:
Figure BDA0001877258740000041
Rstubl=Rshrink1×MRS×Ll/Wl 公式(13)
其中Lshrink1和Rshrink1为缩放因子,MT表示金属的厚度,MRS表示金属的面电阻,Ll为金属短路线的左侧长度,Wl为金属短路线的左侧宽度;Rstubr、Rstubm和Lstubr、Lstubm处理公式相同;从而由公式(9)至公式(11)的实虚部,可直接得到Rl、Rr、Rm和Ll、Lr、Lm的值;
当频率达到一定的数值,趋肤效应变得显著,趋肤效应部分参数的提取中,当Rstubl、Rstubr、Rstubm和Lstubl、Lstubr、Lstubm以及Rl、Rr、Rm和Ll、Lr、Lm的值被确定时,采用高频下的阻抗网络的参数Z减去低频下的参数Z,可得到趋肤效应结构对应的高频等效阻抗,来进行修正,具体关于高频趋肤效应部分的参数的初值过程如下:
金属连接线的电导率为:
condtls=1/(MRS·MT) 公式(14)
金属连接线在最高频率fmax下的趋肤深度为:
Figure BDA0001877258740000051
金属连接线每一部分的趋肤厚度分别为:
h1s=δmax_m/3 公式(16)
h2s=2·δmax_m/3 公式(17)
h3s=MT-2·δmax_m 公式(18)
金属连接线每一部分对应的趋肤宽度分别为:
w1s=2·(WE+MT-2·h1s) 公式(19)
w2s=2·(WE+MT-4·h1s-2·h2s) 公式(20)
w3s=WE-2·δmax_m 公式(21)
根据公式可得,金属连接线高频趋肤各部分的电阻值分别为:
Rli=Rshrink2·Ltotal/(condtls·his·wis)(i=1,2,3) 公式(22)
根据公式可得,金属连接线高频趋肤各部分的电感值分别为:
Figure BDA0001877258740000052
其中,Lshrink2和Rshrink2为缩放因子,MT为金属的厚度,MRS为金属的面电阻,MUO为磁介电常数,Ltotal为左侧金属连接线总长度,WE为金属连接线的有效宽度;
从而可得Rl1,Rl2,Rl3和Ll1,Ll2,Ll3的初始值,并计算得到高频趋肤部分对应的总阻抗参数。
进一步的,开路短路模型拓扑结构的PAD的长度和宽度分别为380um、195um,圆形Signal-PAD的半径是30um,Signal-PAD和Ground-PAD之间的间距是22um。
本发明相比现有技术优点在于:本发明从对Open和Short结构的高频寄生分析出发,尝试由测试物理结构,建立至亚毫米波频段精确的测试结构模型。对建立的模型,其容性和阻性寄生采用解析提取技术,从Open结构低频测试数据中获取。模型的高频趋肤效应采用物理公式计算初值,结合短路测试结构的低频解析提取结果对理论公式进行修正,使之对实际测试结构适用。提出的模型对测试结构的各个子连接结构所产生的高频寄生,进行了准确的区分。模型和模型参数提取方法,采用0.5um InP DHBT工艺上设计所得开路、短路测试结构进行验证,模型仿真和测试所得S参数在0.1-325GHz频段内有很好吻合。
附图说明
图1为Open结构的顶视图和等效电路
图2为Short结构的顶视图和等效电路。
图3为open结构模型等效电路拓扑。
图4为short结构模型等效电路拓扑。
图5为Zl的等效电路模型。
图6为Zr的等效电路模型。
图7为Zgnd的等效电路模型。
图8为简化后的开路等效拓扑结构。
图9为Rmx、Cmx提取示意图。
图10为阻抗网络拓扑结构。
图11为图10的T型等效。
图12为图10等效电路模型。
图13为模型仿真和测试所得开路结构S参数对比。
图14为模型仿真和测试所得短路结构S参数对比。
具体实施方式
下面结合附图和具体实施方式对本发明进一步说明。
如图1至图14所示,一种THz频段InP DHBT器件在片测试结构建模方法,具体包括如下步骤:
101)建立模型步骤:建立短路结构和开路结构,其中短路结构相对于开路结构增加了短路连接线,短路结构和开路结构都包括测试地平面的寄生电感和测试地平面的寄生电阻;开路短路模型拓扑结构的PAD的长度和宽度分别为380um、195um,圆形Signal-PAD的半径是30um,Signal-PAD和Ground-PAD之间的间距是22um。Signal-PAD和待测器件之间连接线的长度为75um,宽度为9um。短路结构和开路结构相似,区别在于短路结构增加了短路连接线(Stub Line),图2中以红色十字结表示,其横向长宽为15um×9um,纵向长宽为21um×8um。图1所示建立的结构开路、短路连接状态同时给出了两个结构可能存在的寄生。和传统建立测试结构模型不同之处在于,本方案着重考虑了测试结构地平面的寄生(图中对其采用4Zgnd表示),实际分析结果表明,地平面阻抗或是引起高频下Open+Short去嵌入算法不精确的主要原因。
图1中Rsl和Csl串联支路表示左侧G-S-G PAD(地-信号端-地焊盘)和金属传输线与地平面之间的寄生电容、寄生电阻,Cstl表示左侧G-S-G PAD和金属传输线与地平面之间的寄生电容。Rsr和Csr串联支路表示右侧G-S-G PAD和金属传输线与地平面之间的寄生电容、寄生电阻,Cstr表示右侧G-S-G PAD和金属传输线与地平面之间的寄生电容。Cdl表示左侧S-PAD和G-PAD之间的寄生电容,Cdr表示右侧S-PAD和G-PAD之间的寄生电容。Zl表示左侧金属连接线的感性寄生网络,Zr表示右侧金属连接线的感性寄生网络,4Zgnd表示地平面的感性寄生网络,Zstub表示金属短路线的寄生网络。考虑版图结构的对称性以及对相关寄生参数的合并化简,图1拓扑结构可简化为如图3所示拓扑结构。Zl、Zr和Zgnd的等效模型如图5至图7所示。附图3中由等效关系可得Csdl=2×Cdl,Csdr=2×Cdr,另外Csfl和Rsfl串联支路表示左侧G-S-G PAD到地平面的寄生电容和电阻,Cstfl表示左侧G-S-G PAD到地平面的寄生电容,Csgl和Rsgl串联支路表示左侧金属连接线到地平面的寄生电容和电阻,Cstgl表示左侧金属连接线到地平面的寄生电容。Csfr和Rsfr串联支路表示右侧G-S-G PAD到地平面的寄生电容和电阻,Cstfr表示右侧G-S-G PAD到地平面的寄生电容,Csgr和Rsgr串联支路表示右侧金属连接线到地平面的寄生电容和电阻,Cstgr表示右侧金属连接线到地平面的寄生电容,其中Rsl=Rsfl//Rsgl,Csl=Csfl//Csgl,Cstl=Cstfl//Cstgl以及Rsr=Rsfr//Rsgr,Csr=Csfr//Csgr,Cstr=Cstfr//Cstgr。Cmx和Rmx串联之路表示输入和输出金属连接线之间的寄生电容和电阻,Cmi表示输入和输出金属连接线之间的寄生电容。Czl和Czr是信号端到地平面的寄生电容。Czx、Czi和Rzx表示地平面的高频寄生电容和电阻。
图5所示Rl、Ll表示金属连接线的感性寄生电阻和电感,Rl1和Ll1串联支路表示连接线的一阶高频趋肤电感和电阻,Rl2和Ll2串联支路表示连接线的二阶高频趋肤电感和电阻,Rl3和Ll3串联支路表示连接线的三阶高频趋肤电感和电阻,由于频率达到THz频段,本方案考虑了三阶高频趋肤,图6中的参数有相似的含义。图7中,Rm、Lm表示地平面的感性寄生电阻和电感,Rm1和Lm1串联支路表示地平面的一阶高频趋肤电感和电阻,Rm2和Lm2串联支路表示地平面的二阶高频趋肤电感和电阻,Rm3和Lm3串联支路表示地平面的三阶高频趋肤电感和电阻。
102)低频下参数获取步骤:低频状况下,地平面的寄生电感、地平面的寄生电阻不起作用,金属连接线理想接地,其等效的拓扑结构包括Cmx、Cmi、Csl、Ctl、Ctr、Csr、Rsl、Rsr、Rmx。其中Cmx和Rmx串联支路表示信号端口之间的寄生电容和电阻;Cmi表示信号端口之间的寄生电容;Csl和Rsl串联支路表示左侧G-S-G PAD和金属传输线与地平面之间的总串联寄生电容和电阻;Ctl表示左侧G-S-G PAD和金属传输线与地平面之间的总并联寄生电容;Csr和Rsr串联支路表示右侧G-S-G PAD和金属传输线与地平面之间的总串联寄生电容和电阻;Ctr表示右侧G-S-G PAD和金属传输线与地平面之间的总并联寄生电容;其中Ctl=Cstl//Csdl,Ctr=Cstr//Csdr
将等效的拓扑结构转换为π型网络,该π型网络的各部分参数Y关系如下:Yl表示Rsl、Csl和Ctl网络,Yr表示Rsr、Csr和Ctr网络,Ym表示Rmx、Cmx和Cmi网络,可依次推导为Ym=-Y21=-Y12,Yl=Y11+Y12,Yr=Y22+Y12,其中Y11、Y12、Y21和Y22分别表示二端口网络中的端口间的导纳。
其中Ym还可根据等效拓扑结构表示为如下公式:
Ym=jωCmi+[Rmx+(jωCmx)-1]-1 公式(1)
其中,j为参数;ω表示角频率,Ym的实部和虚部表达式如下:
[Re(Ym)]-1=Rmx+(ω2Cmx 2Rmx)-1 (2)
Im(Ym)=[Cmi+Cmx(1+ω2Cmx 2Rmx 2)-1]×ω (3)
其中根据实际检测数据,可得[Re(Ym)]-1关于ω-2的截距和斜率,从而依次得到Rmx和Cmx,将提取所得Rmx和Cmx代入公式(3),可得到Cmi的值;
Cmi=ω-1Im(Ym)-Cmx(1+ω2Rmx 2Cmx 2)-1 公式(4)
同理Yl和Yr与Ym有相同的拓扑结构,提取方法同上;从而同理得到Ctl和Ctr,Ctl和Ctr参数提取后,分别由Ctl和Ctr计算Cstl、Csdl以及Cstr、Csdr的值时,引入了比例因子Kl和Kr,计算公式分别为Csdl=KlCtl、Cstl=(1-Kl)Ctl以及Csdr=KrCtr、Cstr=(1-Kr)Cr,其中Kl和Kr的值在实际检测中确定,其位常数;在实验过程中确定最后取值分别为0.155、0.11。
根据已知Signal-PAD的半径和Signal-PAD到待测器件连接线的尺寸以及地平面的长宽,得到其有效面积;由于寄生电容与金属结构面积成正比,寄生电阻与面积成反比,从低频提取所得总的电阻和电容分量中,根据面积比例可分离算出Cstfl、Csfl、Rsfl和Cstgl、Csgl、Rsgl,并得到其值;根据电阻和电容的单位面积因子,得到在高频下考虑的地平面的寄生电容和电阻Czi、Czx和Rzx的值,其中Czx和Rzx串联支路表示地平面的高频寄生电容和电阻,Czi表示地平面的高频寄生电容。对于Czl和Czr采用以下方法计算Czl=2×Cmi,Czl=2×Cmi;其中Czl表示左侧信号端到地平面的寄生电容,Czr表示右侧信号端到地平面的寄生电容。
103)高频下参数获取步骤:高频状况下,金属感性寄生电阻和金属感性寄生电容,其与参数Z之间的关系如下:即针对Short结构,采用Short结构去除Open结构的容性寄生,而后转换为Z参数。
ZShort′=(YShort-YOpen)-1 公式(5)
YShort表示采用Short结构的容性寄生,YOpen表示Open结构的容性寄生,ZShort表示Short结构的参数Z;
高频状况下,金属连接线和地平面的感性寄生电阻和感性寄生电容变得显著,在不考虑各部分的趋肤效应的情况下,其等效电路模型包括Za,Zb和Zm分别代表T形二端口网络中左侧部分、右侧部分和中间部分的阻抗。其中通过参数Z可分别确定如下:
Za=[ZShort′]11-[ZShort′]12 公式(6)
Zb=[ZShort′]22-[ZShort′]12 公式(7)
Zm=[ZShort′]12=[ZShort′]21 公式(8)
Za,Zb和Zm还可进一步表示为:
Za=(Rl+Rstubl)+jω(Ll+Lstubl) 公式(9)
Zb=(Rr+Rstubr)+jω(Lr+Lstubr) 公式(10)
Zm=(Rm+Rstubm)+jω(Lm+Lstubm) 公式(11)
Rl、Ll表示左侧金属连接线的感性寄生电阻和电感,Rr、Lr表示右侧金属连接线的感性寄生电阻和电感,Rm、Lm表示地平面的感性寄生电阻和电感;Rstubl和Lstubl串联支路表示左侧短路连接线的寄生电阻和电感,Rstubr和Lstubr串联支路表示右侧短路连接线的寄生电阻和电感,Rstubm和Lstubm串联支路表示中间短路连接线的寄生电阻和电感。
其中Rstubl和Lstubl的提取公式如下:
Figure BDA0001877258740000111
Rstubl=Rshrink1×MRS×Ll/Wl 公式(13)
其中Lshrink1和Rshrink1为缩放因子,MT表示金属的厚度,MRS表示金属的面电阻,Ll为金属短路线的左侧长度,Wl为金属短路线的左侧宽度;具体初始值如表1。Rstubr、Rstubm和Lstubr、Lstubm处理公式相同;从而由公式(9)至公式(11)的实虚部,可直接得到Rl、Rr、Rm和Ll、Lr、Lm的值;
当频率达到一定的数值,趋肤效应变得显著,趋肤效应部分参数的提取中,当Rstubl、Rstubr、Rstubm和Lstubl、Lstubr、Lstubm以及Rl、Rr、Rm和Ll、Lr、Lm的值被确定时,采用高频下的阻抗网络的参数Z减去低频下的参数Z,可得到趋肤效应结构对应的高频等效阻抗,来进行修正,具体关于高频趋肤效应部分的参数的初值过程如下:
金属连接线的电导率为:
condtls=1/(MRS·MT) 公式(14)
金属连接线在最高频率fmax下的趋肤深度为:
Figure BDA0001877258740000121
金属连接线每一部分的趋肤厚度分别为:
h1s=δmax_m/3 公式(16)
h2s=2·δmax_m/3 公式(17)
h3s=MT-2·δmax_m 公式(18)
其中h1s表示一阶高频趋肤效应的厚度,h2s表示二阶高频趋肤效应的厚度,h3s表示三阶高频趋肤效应的厚度。
金属连接线每一部分对应的趋肤宽度分别为:
w1s=2·(WE+MT-2·h1s) 公式(19)
w2s=2·(WE+MT-4·h1s-2·h2s) 公式(20)
w3s=WE-2·δmax_m 公式(21)
其中w1s表示一阶高频趋肤效应的宽度,w2s表示二阶高频趋肤效应的宽度,w3s表示三阶高频趋肤效应的宽度。
根据公式可得,金属连接线高频趋肤各部分的电阻值分别为:
Rli=Rshrink2·Ltotal/(condtls·his·wis)(i=1,2,3) 公式(22)
根据公式可得,金属连接线高频趋肤各部分的电感值分别为:
Figure BDA0001877258740000122
其中,Lshrink2和Rshrink2为缩放因子,MT为金属的厚度,MRS为金属的面电阻,MUO为磁介电常数,Ltotal为左侧金属连接线总长度,WE为金属连接线的有效宽度;具体初始值如表1:
符号 单位 初始值 物理意义
MT um 2 金属的厚度
MRS 24 金属的面电阻
MUO u 1.257 磁介电常数
Ll um 7.5 金属短路线的左侧长度
Wl um 9 金属短路线的左侧宽度
Ltotal um 75 左侧金属连接线总长度
WE um 9 金属连接线的有效宽度
从而可得Rl1,Rl2,Rl3和Ll1,Ll2,Ll3的初始值,并计算得到高频趋肤部分对应的总阻抗参数。
104)应用验证步骤:为验证,对0.5um InP HBT工艺上制作的测试结构开、短路结构进行测试和建模。模型参数提取以及对参数的优化在Agilent IC-CAP软件中完成,提取所得模型参数值给在表2中,表中的优化值是在初值基础上、对比测试数据微调得到。
表2提取所得测试结构模型参数值以及优化值
Figure BDA0001877258740000131
Figure BDA0001877258740000141
将提取所得模型,在ADS环境下仿真,图13和图14分别给出了模型仿真和测试所得开路、短路结构S参数的对比,在0.2-67GHz,75-110GHz,140-220GHz和220-325GHz的频率范围内,测试和仿真的幅度和相位值有很好的吻合,表明提出的测试结构模型和参数提取方法的有效性。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员,在不脱离本发明构思的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明保护范围内。

Claims (2)

1.一种THz频段InP DHBT器件在片测试结构建模方法,其特征在于,具体包括如下步骤:
101)建立模型步骤:建立短路结构和开路结构,其中短路结构相对于开路结构增加了短路连接线,短路结构和开路结构都包括测试地平面的寄生电感和测试地平面的寄生电阻;
102)低频下参数获取步骤:低频状况下,地平面的寄生电感、地平面的寄生电阻不起作用,金属连接线理想接地,其等效的拓扑结构包括Cmx、Cmi、Csl、Ctl、Ctr、Csr、Rsl、Rsr、Rmx;其中Cmx和Rmx串联支路表示信号端口之间的寄生电容和电阻;Cmi表示信号端口之间的寄生电容;Csl和Rsl串联支路表示左侧G-S-G PAD和金属传输线与地平面之间的总串联寄生电容和电阻;Ctl表示左侧G-S-G PAD和金属传输线与地平面之间的总并联寄生电容;Csr和Rsr串联支路表示右侧G-S-G PAD和金属传输线与地平面之间的总串联寄生电容和电阻;Ctr表示右侧G-S-G PAD和金属传输线与地平面之间的总并联寄生电容;其中Ctl=Cstl//Csdl,Ctr=Cstr//Csdr
将等效的拓扑结构转换为π型网络,该π型网络的各部分参数Y关系如下:Yl表示Rsl、Csl和Ctl网络,Yr表示Rsr、Csr和Ctr网络,Ym表示Rmx、Cmx和Cmi网络,可依次推导为Ym=-Y21=-Y12,Yl=Y11+Y12,Yr=Y22+Y12,其中Y11、Y12、Y21和Y22分别表示二端口网络中的端口间的导纳;
其中Ym还可根据等效拓扑结构表示为如下公式:
Ym=jωCmi+[Rmx+(jωCmx)-1]-1 公式(1)
其中,j为参数;ω表示角频率;
Ym的实部和虚部表达式如下:
[Re(Ym)]-1=Rmx+(ω2Cmx 2Rmx)-1 公式(2)
Im(Ym)=[Cmi+Cmx(1+ω2Cmx 2Rmx 2)-1]×ω 公式(3)
其中根据实际检测数据,可得[Re(Ym)]-1关于ω-2的截距和斜率,从而依次得到Rmx和Cmx,将提取所得Rmx和Cmx代入公式(3),可得到Cmi的值;
Cmi=ω-1Im(Ym)-Cmx(1+ω2Rmx 2Cmx 2)-1 公式(4)
同理Yl和Yr与Ym有相同的拓扑结构,提取方法同上;从而同理得到Ctl和Ctr,Ctl和Ctr参数提取后,分别由Ctl和Ctr计算Cstl、Csdl以及Cstr、Csdr的值时,引入了比例因子Kl和Kr,计算公式分别为Csdl=KlCtl、Cstl=(1-Kl)Ctl以及Csdr=KrCtr、Cstr=(1-Kr)Cr,其中Kl和Kr的值在实际检测中确定,其位常数;
根据已知Signal-PAD的半径和Signal-PAD到待测器件连接线的尺寸以及地平面的长宽,得到其有效面积;由于寄生电容与金属结构面积成正比,寄生电阻与面积成反比,从低频提取所得总的电阻和电容分量中,根据面积比例可分离算出Cstfl、Csfl、Rsfl和Cstgl、Csgl、Rsgl,并得到其值;根据电阻和电容的单位面积因子,得到在高频下考虑的地平面的寄生电容和电阻Czi、Czx和Rzx的值,其中Czx和Rzx串联支路表示地平面的高频寄生电容和电阻,Czi表示地平面的高频寄生电容,对于Czl和Czr采用以下方法计算Czl=2×Cmi,Czl=2×Cmi;其中Czl表示左侧信号端到地平面的寄生电容,Czr表示右侧信号端到地平面的寄生电容;
103)高频下参数获取步骤:高频状况下,金属感性寄生电阻和金属感性寄生电容,其与参数Z之间的关系如下:
ZShort′=(YShort-YOpen)-1 公式(5)
YShort表示采用Short结构的容性寄生,YOpen表示Open结构的容性寄生,ZShort表示Short结构的参数Z;
高频状况下,金属连接线和地平面的感性寄生电阻和感性寄生电容变得显著,在不考虑各部分的趋肤效应的情况下,其等效电路模型包括Za,Zb和Zm分别代表T形二端口网络中左侧部分、右侧部分和中间部分的阻抗,其中通过参数Z可分别确定如下:
Za=[ZShort′]11-[ZShort′]12 公式(6)
Zb=[ZShort′]22-[ZShort′]12 公式(7)
Zm=[ZShort′]12=[ZShort′]21 公式(8)
Za,Zb和Zm还可进一步表示为:
Za=(Rl+Rstubl)+jω(Ll+Lstubl) 公式(9)
Zb=(Rr+Rstubr)+jω(Lr+Lstubr) 公式(10)
Zm=(Rm+Rstubm)+jω(Lm+Lstubm) 公式(11)
Rl、Ll表示左侧金属连接线的感性寄生电阻和电感,Rr、Lr表示右侧金属连接线的感性寄生电阻和电感,Rm、Lm表示地平面的感性寄生电阻和电感;Rstubl和Lstubl串联支路表示左侧短路连接线的寄生电阻和电感,Rstubr和Lstubr串联支路表示右侧短路连接线的寄生电阻和电感,Rstubm和Lstubm串联支路表示中间短路连接线的寄生电阻和电感;
其中Rstubl和Lstubl的提取公式如下:
Figure QLYQS_1
Rstubl=Rshrink1×MRS×Ll/Wl 公式(13)
其中Lshrink1和Rshrink1为缩放因子,MT表示金属的厚度,MRS表示金属的面电阻,Ll为金属短路线的左侧长度,Wl为金属短路线的左侧宽度;Rstubr、Rstubm和Lstubr、Lstubm处理公式相同;从而由公式(9)至公式(11)的实虚部,可直接得到Rl、Rr、Rm和Ll、Lr、Lm的值;
当频率达到一定的数值,趋肤效应变得显著,趋肤效应部分参数的提取中,当Rstubl、Rstubr、Rstubm和Lstubl、Lstubr、Lstubm以及Rl、Rr、Rm和Ll、Lr、Lm的值被确定时,采用高频下的阻抗网络的参数Z减去低频下的参数Z,可得到趋肤效应结构对应的高频等效阻抗,来进行修正,具体关于高频趋肤效应部分的参数的初值过程如下:
金属连接线的电导率为:
condtls=1/(MRS·MT) 公式(14)
金属连接线在最高频率fmax下的趋肤深度为:
Figure QLYQS_2
金属连接线每一部分的趋肤厚度分别为:
h1s=δmax_m/3 公式(16)
h2s=2·δmax_m/3 公式(17)
h3s=MT-2·δmax_m 公式(18)
其中h1s表示一阶高频趋肤效应的厚度,h2s表示二阶高频趋肤效应的厚度,h3s表示三阶高频趋肤效应的厚度;
金属连接线每一部分对应的趋肤宽度分别为:
w1s=2·(WE+MT-2·h1s) 公式(19)
w2s=2·(WE+MT-4·h1s-2·h2s) 公式(20)
w3s=WE-2·δmax_m 公式(21)
其中w1s表示一阶高频趋肤效应的宽度,w2s表示二阶高频趋肤效应的宽度,w3s表示三阶高频趋肤效应的宽度;
根据公式可得,金属连接线高频趋肤各部分的电阻值分别为:
Rli=Rshrink2·Ltotal/(condtls·his·wis)(i=1,2,3) 公式(22)
根据公式可得,金属连接线高频趋肤各部分的电感值分别为:
Figure QLYQS_3
其中,Lshrink2和Rshrink2为缩放因子,MT为金属的厚度,MRS为金属的面电阻,MUO为磁介电常数,Ltotal为左侧金属连接线总长度,WE为金属连接线的有效宽度;
从而可得Rl1,Rl2,Rl3和Ll1,Ll2,Ll3的初始值,并计算得到高频趋肤部分对应的总阻抗参数。
2.根据权利要求1所述的一种THz频段InP DHBT器件在片测试结构建模方法,其特征在于,开路短路模型拓扑结构的PAD的长度和宽度分别为380um、195um,圆形Signal-PAD的半径是30um,Signal-PAD和Ground-PAD之间的间距是22um。
CN201811405567.1A 2018-11-23 2018-11-23 一种THz频段InP DHBT器件在片测试结构建模方法 Active CN110008489B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811405567.1A CN110008489B (zh) 2018-11-23 2018-11-23 一种THz频段InP DHBT器件在片测试结构建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811405567.1A CN110008489B (zh) 2018-11-23 2018-11-23 一种THz频段InP DHBT器件在片测试结构建模方法

Publications (2)

Publication Number Publication Date
CN110008489A CN110008489A (zh) 2019-07-12
CN110008489B true CN110008489B (zh) 2023-07-14

Family

ID=67165004

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811405567.1A Active CN110008489B (zh) 2018-11-23 2018-11-23 一种THz频段InP DHBT器件在片测试结构建模方法

Country Status (1)

Country Link
CN (1) CN110008489B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111781479B (zh) * 2020-02-13 2022-11-08 中国电子科技集团公司第十三研究所 一种在片校准件模型建立方法
CN112069755B (zh) * 2020-09-10 2023-12-29 浪潮商用机器有限公司 Pcb电源的阻抗检测方法、系统、设备及存储介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102375900A (zh) * 2010-08-19 2012-03-14 上海华虹Nec电子有限公司 考虑趋肤效应的模拟射频mim电容的电路结构
CN106202835A (zh) * 2016-08-16 2016-12-07 南京展芯通讯科技有限公司 包含沟道高级寄生元件的场效应晶体管小信号等效电路模型
CN106777483B (zh) * 2016-11-18 2019-10-11 东南大学 用于集成电路的片上电感等效电路模型及参数提取方法

Also Published As

Publication number Publication date
CN110008489A (zh) 2019-07-12

Similar Documents

Publication Publication Date Title
US6560567B1 (en) Method and apparatus for measuring on-wafer lumped capacitances in integrated circuits
EP2499739B1 (de) Verfahren zur bestimmung eines anpassnetzwerks basierend auf poincaré-kreisen von reflexionsfaktoren
CN110008489B (zh) 一种THz频段InP DHBT器件在片测试结构建模方法
CN102073757B (zh) 一种用于集成电路中电感模型的分析方法
DE10338072A1 (de) Verfahren und Vorrichtung zum Kalibrieren eines Meßvorrichtungswegs und zum Messen einer Testvorrichtung in dem kalibrierten Meßvorrichtungsweg
Groteluschen et al. Quasi-analytical analysis of the broadband properties of multiconductor transmission lines on semiconducting substrates
CN106777483B (zh) 用于集成电路的片上电感等效电路模型及参数提取方法
CN105552493A (zh) 平行耦合微带线带通滤波器的直接设计方法
CN109522653A (zh) 一种三维集成电路的pdn交流噪声分析方法
CN105022878A (zh) 射频soi-mos变容管衬底模型及其参数提取方法
CN107167724A (zh) 一种小信号量测的去嵌方法
CN111679171A (zh) 基于互连线单元的电路拓扑结构及互连线单元级联的去嵌方法
CN107480397A (zh) 考虑双根金丝键合线互感的微波器件路耦合传输性能预测方法
CN113655360A (zh) Rf mos器件的在片测试结构的去嵌方法
CN104933270A (zh) 带电容电源地平面建模及电容去耦半径仿真方法
CN102375101B (zh) 采用不同层金属引线连出的无源器件测试去嵌方法
US8046208B1 (en) Method and apparatus for representing high speed interconnects in closed forms
CN106021670B (zh) 一种毫米波fet的建模方法
CN104635135A (zh) 一种化合物半导体器件的去嵌入方法
CN103943922A (zh) 谐波抑制带通滤波器及制造方法
CN209692710U (zh) 一种基于多级二次谐波控制的宽带连续型功率放大器
JP2010072684A (ja) 回路特性解析方法、装置、およびプログラム
CN109581071A (zh) 一种直接变频结构的复介电常数测量系统
Abd El-Raouf et al. Application of the Z-transform technique to modelling linear lumped loads in the FDTD
CN104578755B (zh) 电源式电磁干扰抑制滤波器

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant