CN108052727B - 一种金属栅功函数变化导致栅电容统计分布的估计方法 - Google Patents
一种金属栅功函数变化导致栅电容统计分布的估计方法 Download PDFInfo
- Publication number
- CN108052727B CN108052727B CN201711298716.4A CN201711298716A CN108052727B CN 108052727 B CN108052727 B CN 108052727B CN 201711298716 A CN201711298716 A CN 201711298716A CN 108052727 B CN108052727 B CN 108052727B
- Authority
- CN
- China
- Prior art keywords
- gate capacitance
- statistical distribution
- gate
- metal gate
- work function
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明涉及金属栅电容统计分布的估计方法。一种金属栅功函数变化导致栅电容统计分布的估计方法,包括如下步骤:通过HSPICE软件,设定相关的变异参数,确定绝对标准偏差的大小,并进行10000次的蒙特卡罗仿真。提取MOSFET器件的栅电容参数C。将数据读取到MATLAB中并应用其统计工具箱,进行概率密度函数和累积概率函数拟合并表征;用假设检验判断统计的正确性。本发明利用密度函数、累计概率函数拟合并表征,使用假设检验判断,能在CMOS器件和电路设计早期,快速精确地预测由于MG‑WFV效应导致实际纳米器件和电路的制造性能变化的统计分布,以减少实际制造集成电路芯片性能的盲目性。
Description
技术领域
本发明涉及集成电路设计领域,具体涉及栅电容统计分布的估计方法。
背景技术
随着平面CMOS集成电路技术节点进入亚32nm尺度,集成器件和集成电路芯片的制造总共涉及几千个不同的工艺步骤。这些步骤,通过氧化、光刻、沉积等技术手段的物理、化学过程将MOSFET器件集成在半导体硅片上。但是,该过程实现的工艺步骤中并不能做到对硅片内部的微小粒子进行非常精确地控制,结果是引起MOSFET器件的工艺参数发生波动,这些现象严重影响MOSFET器件和电路性能,导致实际制造的集成芯片(IC)性能严重退化和成品率损失。
引起这些现象的波动源主要包括:金属栅功函数随机变化(MG-WFV)、随机离散掺杂(RDD)、线边缘粗糙(LER)等。目前,MG-WFV已经成为影响器件和电路性能不确定性的一个重要因数。这是因为:金属栅极是由大小随机差异的几个数纳米直径大小且取向不同的晶粒通过原子层沉积(ALD)过程完成,因为不同取向的晶粒具有不同的金属-半导体功函数差(Фms),导致该功函数在不同的MOSFET之间发生随机变化。以TiN金属栅为例,它具有4.4电子伏特(111)和4.6电子伏特(100)两种不同的功函数差,其出现的概率分别为40%和60%。显然,作为阈值电压组成的一部分,这将引起阈值电压的随机变化。但是,尽管目前工作原理尚不明确,MG-WFV也同样会导致栅电容(Cgg)发生随机变化现象。而这种现象,目前主要是通过采用大量样本数据得到的柱状图进行分析判断。这种方法受人为因素影响较大比较,对于大样本数据采集耗时,不易且会存在一定人为统计误差。再者统计作图预测方法传统,精确度较低。
发明内容
本发明的目的在于,提供一种金属栅功函数变化导致栅电容统计分布的估计方法,解决以上技术问题。
本发明所解决的技术问题可以采用以下技术方案来实现:
一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于,通过HSPICE软件,设定纳米尺度的MOSFET器件的模型参数与功函数差直接相关的一个变异参数;确定晶粒面积比值,并确定一个绝对标准偏差,进行多次蒙特卡罗仿真;
通过“LX18”函数,提取MOSFET器件内部参数栅电容Cgg,描述该性能参数及其变化特性;
在MATLAB软件中安装HSPICE Toolbox工具箱,将蒙特卡罗仿真产生的大量样本数据读取到MATLAB软件的相应变量中;
将MOSFET器件性能参数及其变化数据信息,应用HSPICE Toolbox工具箱,进行Burr分布的概率密度函数(PDF)和累积概率函数(CDF)拟合并表征。
本发明提供了金属栅功函数变化导致栅电容(Cgg)统计分布的快速、精确的估计方法,这种方法可以替代现仅依赖样本数据判断其性能变化的不足,且用统计学的理论来检验得到的统计分布是否可靠,构成了一个完成的性能估计体系。
本发明所使用的HSPICE软件、MATLAB软件都是现有的,本专利的创新点在于:
1.创造性的将HSPICE软件、MATLAB软件应用到金属栅功函数变化导致栅电容变化的统计分析上,可替代人工分析;
2.将两个软件进行了结合,并优化了分析过程,得出的估计结果更为精确、分析过程更加快捷。
3.本发明能在CMOS器件和电路设计早期,快速精确地预测由于MG-WFV效应导致实际纳米器件和电路的制造性能变化的统计分布,以减少实际制造集成电路芯片性能的盲目性。
可以根据实际工艺技术水平和/或金属栅面积,确定晶粒面积比值。
可以结合金属栅的金属类型确定一个绝对标准偏差。
蒙特卡罗仿真的次数优选不少于10000次,进一步优选10000次。
可以用假设检验判断统计的正确性。例如,对拟合结果的可行性通过最大似然初步判断,并进行卡方检验(设置置信概率0.05)确认其统计拟合特性的正确性。
优选的具体步骤如下:
步骤一,设定纳米尺度的MOSFET器件的模型参数与功函数差直接相关的一个变异参数;根据金属栅面积,确定晶粒面积比值,并结合金属栅的金属类型确定一个绝对标准偏差;
步骤二,将所述模型参数,所述变异参数,所述晶粒面积比值,所述标准偏差输入仿真软件,并进行蒙特卡罗仿真模拟;得到金属栅的样本空间;
步骤三,计算出所述样本空间中每个金属栅样本的栅电容C;
步骤四,将所述样本空间读取到MATLAB中;
步骤五,将步骤三中所述栅电容C读取到MATLAB中;
步骤六,利用MATLAB中的Burr统计分布,对仿真的数据进行概率密度函数和累积概率函数拟合并表征;
步骤七,对步骤六中拟合的结果进行最大似然判断,并进行假设检验;
当假设检验置信概率大于95%,则估计过程结束;否则,从步骤二重新开始。
本发明利用数学概率论原理,进行概率密度函数和累积概率函数拟合并表征,避免直接根据大量样本数据得到的柱状图进行分析判断中的偶然性与表面性。较传统方法,提高了统计估计判断的准确性与科学性。
附图说明
图1为本发明的流程图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示进一步阐述本发明。
一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于,通过HSPICE软件,设定纳米尺度的MOSFET器件的模型参数与功函数差直接相关的一个变异参数;根据实际工艺技术水平和/或金属栅面积,确定晶粒面积比值,并结合金属栅的金属类型确定一个绝对标准偏差,进行不少于10000次的蒙特卡罗仿真;
通过“LX18”函数,提取MOSFET器件内部参数栅电容Cgg,描述该性能参数及其变化特性;
在MATLAB软件中安装HSPICE Toolbox工具箱,将蒙特卡罗仿真产生的大量样本数据读取到MATLAB软件的相应变量中;
将MOSFET器件性能参数及其变化数据信息,应用HSPICE Toolbox工具箱,进行Burr分布的概率密度函数(PDF)和累积概率函数(CDF)拟合并表征。
本发明提供了金属栅功函数变化导致栅电容(Cgg)统计分布的快速、精确的估计方法,这种方法可以替代现仅依赖样本数据判断其性能变化的不足,且用统计学的理论来检验得到的统计分布是否可靠,构成了一个完成的性能估计体系。
本发明所使用的HSPICE软件、MATLAB软件都是现有的,本专利的创新点在于:
4.创造性的将HSPICE软件、MATLAB软件应用到金属栅功函数变化导致栅电容变化的统计分析上,可替代人工分析;
5.将两个软件进行了结合,并优化了分析过程,得出的估计结果更为精确、分析过程更加快捷。
6.本发明能在CMOS器件和电路设计早期,快速精确地预测由于MG-WFV效应导致实际纳米器件和电路的制造性能变化的统计分布,以减少实际制造集成电路芯片性能的盲目性。
可以用假设检验判断统计的正确性。例如,对拟合结果的可行性通过最大似然初步判断,并进行卡方检验(设置置信概率0.05)确认其统计拟合特性的正确性。
具体实施例:
步骤一,设定纳米尺度的MOSFET器件的模型参数与功函数差直接相关的一个变异参数;根据金属栅面积,确定晶粒面积比值,并结合不同金属确定一个晶粒面积比值绝对标准偏差的大小;
步骤二,将模型参数,变异参数,晶粒面积比值,标准偏差输入仿真软件,并进行蒙特卡罗仿真模拟;得到金属栅的样本空间;
步骤三,计算出样本空间中每组金属栅样本的栅电容C;
步骤四,将样本空间读取到MATLAB中;
步骤五,将步骤三中栅电容读取到MATLAB中;
步骤六,利用MATLAB中的Burr统计分布,对仿真的数据进行概率密度函数和累积概率函数拟合并表征;
步骤七,对步骤六中拟合的结果进行最大似然判断,并进行假设检验;
当假设检验置信概率大于95%,则估计过程结束;否则,从步骤二重新开始。
本发明利用数学概率论原理,进行概率密度函数和累积概率函数拟合并表征,避免直接根据大量样本数据得到的柱状图进行分析判断中的偶然性与表面性。较传统方法,提高了统计估计判断的准确性与科学性。
进一步,步骤二中仿真软件为HSPICE;仿真模拟为蒙特卡罗频率大于10000。利用HSPICE软件进行蒙特卡罗仿真,避免了现实中大量实验数据搜集的困难。进一步,步骤三中,通过LX18函数计算出每次蒙特卡罗仿真的样本栅电容C;LX18函数是一物理公式,即C=ε*S/4*π*k*d;式中:C为电容,单位F;ε为介电常数;S为电容器两极板的正对面积,单位㎡;d为极板间距,单位M;k为静电力常量,大小为9.0*109N。进一步,步骤四,步骤五中,利用MATLAB的工具箱中loadsig和evalsig函数进行读取操作;能够快速地把每次仿真得到的样本栅电容读取进MATLAB中。进一步,步骤六中,统计栅电容的分布;统计出栅电容不同分布的频率数据。进一步,步骤六中,以栅电容为随机变量,运用MATLAB中的Burr统计分布对仿真的数据进行概率密度函数和累积概率函数拟合;并使用MATLAB将拟合结果直接图形化显示。将仿真产生的大量样本数据,读取进MATLAB中,充分发挥MATLAB软件进行数学、科学类大样本数据各种形式的计算能力,并能通过其内在的数据图形化能力,将结果直接图形化显示,使得结论更加直观。进一步,步骤七中,再次用时HSPICE软件仿真出一组数据,并计算出相应的电容;对上述拟合的结果进行最大似然判断。进一步,步骤七中,再次对拟合结果进行假设检验;假设检验为u—检验法、t检验法、卡方检验、F—检验法,秩和检验中的一种。进一步,步骤七中,假设检验设置置信概率0.05。
本发明通过最大释然判断和假设检验对拟合结果进行判断,双重检验结果是否可信。体现了此方法的严谨性,科学性。这样使得最终结论真实,可靠。
以上显示和描述了本发明的基本原理和主要特征以及本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (9)
1.一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于,包括如下步骤:
步骤一,设定纳米尺度的MOSFET器件的模型参数与功函数差直接相关的一个变异参数;确定晶粒面积比值,确定一个绝对标准偏差;
步骤二,将所述模型参数,所述变异参数,所述晶粒面积比值,所述标准偏差输入仿真软件,并进行蒙特卡罗仿真模拟;得到金属栅的样本空间;
步骤三,计算出所述样本空间中每组金属栅样本的栅电容C;
步骤四,将所述样本空间读取到MATLAB中;
步骤五,将步骤三中所述栅电容读取到MATLAB中;
步骤六,利用MATLAB中的Burr统计分布,对仿真的数据进行概率密度函数和累积概率函数拟合并表征;
步骤七,对步骤六中拟合的结果进行最大似然判断,并进行假设检验;
当假设检验置信概率大于0.95,则估计过程结束;否则,从步骤二重新开始。
2.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤二中,仿真软件为HSPICE;所述仿真模拟为蒙特卡罗频率大于10000。
3.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤三中,通过LX18函数计算出每次蒙特卡罗仿真的金属栅样本栅电容C;
所述LX18函数是一物理公式,即C=ε*S/4*π*k*d;式中:C为电容,单位F;ε为介电常数;S为电容器两极板的正对面积,单位㎡;d为极板间距,单位M;k为静电力常量,大小为9.0*109N。
4.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤四、步骤五中,利用所述MATLAB的工具箱中loadsig和evalsig函数进行读取操作。
5.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤六中,统计读取进MATLAB的金属栅样本栅电容的分布;统计出栅电容不同分布的频率数据。
6.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤六中,以金属栅样本栅电容为随机变量,运用MATLAB中的Burr统计分布对仿真的数据进行概率密度函数和累积概率函数拟合;并使用MATLAB将拟合结果直接图形化显示。
7.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤七中,使用HSPICE软件仿真出一组数据,并计算出相应金属栅样本的栅电容;对上述拟合的结果进行最大似然判断。
8.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤七中,对所述拟合结果进行假设检验;所述假设检验为u—检验法、t检验法、卡方检验、F—检验法,秩和检验中的一种。
9.根据权利要求1所述的一种金属栅功函数变化导致栅电容统计分布的估计方法,其特征在于:所述步骤七中,假设检验设置置信概率0.05。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711298716.4A CN108052727B (zh) | 2017-12-08 | 2017-12-08 | 一种金属栅功函数变化导致栅电容统计分布的估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711298716.4A CN108052727B (zh) | 2017-12-08 | 2017-12-08 | 一种金属栅功函数变化导致栅电容统计分布的估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108052727A CN108052727A (zh) | 2018-05-18 |
CN108052727B true CN108052727B (zh) | 2021-05-14 |
Family
ID=62123713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711298716.4A Active CN108052727B (zh) | 2017-12-08 | 2017-12-08 | 一种金属栅功函数变化导致栅电容统计分布的估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108052727B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109508500B (zh) * | 2018-11-16 | 2022-11-25 | 杭州电子科技大学 | 一种估计金属栅晶粒随机取向引起FinFET阈值统计分布的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201852A (zh) * | 2006-12-13 | 2008-06-18 | 上海华虹Nec电子有限公司 | 用于等比例缩小工艺的spice模型方法 |
CN102436529A (zh) * | 2011-11-14 | 2012-05-02 | 华东师范大学 | 基于Weibull分布的MOS晶体管可靠性统计模型的建模方法 |
CN102916030A (zh) * | 2011-08-03 | 2013-02-06 | 尹根千 | 有机电致发光显示元件以及其制造方法 |
CN103474455A (zh) * | 2013-08-21 | 2013-12-25 | 电子科技大学 | 一种具有复合金属栅的氮化镓基高电子迁移率晶体管 |
US9728629B1 (en) * | 2016-07-14 | 2017-08-08 | Semiconductor Components Industries, Llc | Electronic device including a polycrystalline compound semiconductor layer and a process of forming the same |
CN107292026A (zh) * | 2017-06-21 | 2017-10-24 | 杭州电子科技大学 | 一种工艺参数波动引起mosfet性能变化的估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011040541A (ja) * | 2009-08-10 | 2011-02-24 | Sony Corp | ゲート絶縁膜の絶縁破壊寿命の評価方法、ゲート絶縁膜の絶縁破壊寿命の評価装置、ゲート絶縁膜の絶縁破壊寿命の評価用のプログラム |
CN104393026A (zh) * | 2014-12-12 | 2015-03-04 | 京东方科技集团股份有限公司 | Oled显示基板及其制作方法、显示装置 |
-
2017
- 2017-12-08 CN CN201711298716.4A patent/CN108052727B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201852A (zh) * | 2006-12-13 | 2008-06-18 | 上海华虹Nec电子有限公司 | 用于等比例缩小工艺的spice模型方法 |
CN102916030A (zh) * | 2011-08-03 | 2013-02-06 | 尹根千 | 有机电致发光显示元件以及其制造方法 |
CN102436529A (zh) * | 2011-11-14 | 2012-05-02 | 华东师范大学 | 基于Weibull分布的MOS晶体管可靠性统计模型的建模方法 |
CN103474455A (zh) * | 2013-08-21 | 2013-12-25 | 电子科技大学 | 一种具有复合金属栅的氮化镓基高电子迁移率晶体管 |
US9728629B1 (en) * | 2016-07-14 | 2017-08-08 | Semiconductor Components Industries, Llc | Electronic device including a polycrystalline compound semiconductor layer and a process of forming the same |
CN107292026A (zh) * | 2017-06-21 | 2017-10-24 | 杭州电子科技大学 | 一种工艺参数波动引起mosfet性能变化的估计方法 |
Non-Patent Citations (2)
Title |
---|
"Analytical modeling of cutoff frequency variability reserving correlations due to random dopant fluctuation in nanometer MOSFETs";Wiefeng Lv等;《Solid-State Electronics》;20150315;第63-69页 * |
"Statistical Modeling of Gate Capacitance Variations Induced by Random Dopants in Nanometer MOSFETs Reserving Correlations";lv weifeng 等;《Chinese Physics Letters》;20151001;第32卷(第10期);108502-1至108502-3页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108052727A (zh) | 2018-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI796494B (zh) | 高效積體電路模擬及測試 | |
TWI451336B (zh) | 預測模型之建模樣本的篩選方法及其電腦程式產品 | |
US10169506B2 (en) | Circuit design method and system | |
US8806418B1 (en) | Scaled sigma sampling | |
TW201917613A (zh) | 半導體裝置製程 | |
TW200907733A (en) | Robust design using manufacturability models | |
JP5558758B2 (ja) | 半導体素子の漏洩電流予測方法 | |
WO2011078965A1 (en) | Methods and systems for high sigma yield estimation using reduced dimensionality | |
JP2010160787A (ja) | パラメータ情報作成システム、歩留まり算出システム、プログラム及び記録媒体 | |
US10747916B2 (en) | Parameter generation for modeling of process-induced semiconductor device variation | |
WO2011078930A1 (en) | Methods and systems for high sigma yield estimation | |
US20160342719A1 (en) | Parameter generation for semiconductor device trapped-charge modeling | |
CN108052727B (zh) | 一种金属栅功函数变化导致栅电容统计分布的估计方法 | |
TWI709902B (zh) | 半導體製程模擬方法、半導體製程模擬裝置及執行半導體製程模擬的計算裝置 | |
Pan et al. | Reliability analysis of analog circuits using quadratic lifetime worst-case distance prediction | |
US20040139419A1 (en) | Minimization of microelectronic interconnect thickness variations | |
CN113111620B (zh) | 一种半导体电路良率预测方法及装置 | |
CN107292026A (zh) | 一种工艺参数波动引起mosfet性能变化的估计方法 | |
Liu et al. | Performance bound and yield analysis for analog circuits under process variations | |
KR102445989B1 (ko) | 전자 현미경을 사용한 반도체 계측 및 결함 분류 | |
Chan et al. | Design dependent process monitoring for back-end manufacturing cost reduction | |
US9378314B2 (en) | Analytical model for predicting current mismatch in metal oxide semiconductor arrays | |
Blue et al. | Spatial variance spectrum analysis and its application to unsupervised detection of systematic wafer spatial variations | |
McAndrew et al. | Unified statistical modeling for circuit simulation | |
CN102945303A (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 |