CN101630339A - 一种获得考虑版图相关应力后电路性能的方法 - Google Patents
一种获得考虑版图相关应力后电路性能的方法 Download PDFInfo
- Publication number
- CN101630339A CN101630339A CN200910091357A CN200910091357A CN101630339A CN 101630339 A CN101630339 A CN 101630339A CN 200910091357 A CN200910091357 A CN 200910091357A CN 200910091357 A CN200910091357 A CN 200910091357A CN 101630339 A CN101630339 A CN 101630339A
- Authority
- CN
- China
- Prior art keywords
- transistor
- layout
- stress
- mobility
- domain
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000009826 distribution Methods 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims abstract description 13
- 238000000605 extraction Methods 0.000 claims description 19
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 238000012821 model calculation Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 7
- 238000010168 coupling process Methods 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 239000004065 semiconductor Substances 0.000 claims description 7
- 239000013078 crystal Substances 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 229910021420 polycrystalline silicon Inorganic materials 0.000 claims description 6
- 229920005591 polysilicon Polymers 0.000 claims description 6
- 239000002210 silicon-based material Substances 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 230000002035 prolonged effect Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 8
- 239000000463 material Substances 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 6
- 238000002955 isolation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 3
- 230000003647 oxidation Effects 0.000 description 3
- 238000007254 oxidation reaction Methods 0.000 description 3
- 229910052710 silicon Inorganic materials 0.000 description 3
- 239000010703 silicon Substances 0.000 description 3
- 229910000577 Silicon-germanium Inorganic materials 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000000137 annealing Methods 0.000 description 2
- 230000003139 buffering effect Effects 0.000 description 2
- 238000005530 etching Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 229910052581 Si3N4 Inorganic materials 0.000 description 1
- INQLNSVYIFCUML-QZTLEVGFSA-N [[(2r,3s,4r,5r)-5-(6-aminopurin-9-yl)-3,4-dihydroxyoxolan-2-yl]methoxy-hydroxyphosphoryl] [(2r,3s,4r,5r)-5-(4-carbamoyl-1,3-thiazol-2-yl)-3,4-dihydroxyoxolan-2-yl]methyl hydrogen phosphate Chemical compound NC(=O)C1=CSC([C@H]2[C@@H]([C@H](O)[C@@H](COP(O)(=O)OP(O)(=O)OC[C@@H]3[C@H]([C@@H](O)[C@@H](O3)N3C4=NC=NC(N)=C4N=C3)O)O2)O)=N1 INQLNSVYIFCUML-QZTLEVGFSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- HQVNEWCFYHHQES-UHFFFAOYSA-N silicon nitride Chemical compound N12[Si]34N5[Si]62N3[Si]51N64 HQVNEWCFYHHQES-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Design And Manufacture Of Integrated Circuits (AREA)
Abstract
本发明涉及一种获得考虑版图相关应力后电路性能的方法,属于集成电路设计技术领域,该方法包括:提取出版图各个晶体管的区域;对版图中每个晶体管划分出有效区域,以将大规模电路划分成小单元;依据版图构造每个晶体管有效区域的三维结构,设定初始条件及边界条件,采用通用的有限单元方法求解,得到每个晶体管有效区域的三维应力分布;计算考虑版图相关应力影响之后各个晶体管的迁移率,使用该迁移率更新原有的晶体管模型,使用新的晶体管模型进行计算,得到考虑版图相关应力影响之后的电路特性。本发明具有网格更新过程简化,提高了模拟速度和规模,计算结果精确,可以处理复杂的版图结构等特点。
Description
技术领域
本发明属于集成电路设计技术领域,涉及考虑版图相关应力对电路特性影响的集成电路设计,尤其涉及获得STI工艺引入应力对电路特性影响的方法。
背景技术
随着集成电路技术的发展,应力作为一种影响晶体管器件性能的因素得到了广泛的研究。多种不同的工艺被采用来人为引入额外的应力以获得所需的器件性能改变。一种典型的应变沟道技术通过在SiGe衬底上外延生长一层用于制作器件的Si层。由于Si与SiGe层的晶格常数不同,在Si层中将引入沟道平面内的双轴拉应力。这一应力的引入影响Si的价带特性,使得电子的有效质量减小,迁移率增大,因而可以获得更大的驱动电流和更快的器件速度。载流子迁移率描述了半导体材料中载流子运动受外加电场影响的特性。半导体材料如硅中的机械应力会改变材料的能带特性,进而影响包括载流子迁移率在内的多方面因素。对于硅材料器件,沟道平面内的双轴拉应力会改善nMOS器件的特性,而沿着沟道方向的单轴压应力有助于提高pMOS器件的载流子迁移率。
在人为引入并得到良好控制的应力因素之外,部分与电路的版图特性直接相关的应力来源也起到了显著的影响。在集成电路制作工艺中,每个集成电路需要在不同的层次上进行描述。对于同一个集成电路,电路网单描述了电路的拓扑结构以及各晶体管的参数,版图描述了制造集成电路时的各工艺步骤形成材料的几何形状。版图使用多层不同的掩模版描述制造过程中对应的多层不同材料的几何形状。在这些层中,每个多晶硅层与有源区层的独立重合区域构成一个沟道区,每个这样的沟道区对应一个晶体管。版图提取工具用于分析集成电路的版图以及电路网单,获得版图对应的器件及其连接关系,并得到电路网单与版图中各晶体管的对应关系。集成电路制造工艺中的很多步骤都会引入额外的应力,这些工艺包括材料的生长、淀积、氧化、刻蚀以及退火等等。这些工艺引入的额外应力,或者称为无意引入的应力,会叠加在人为引入的应力上,改变晶体管中的应力分布,进而影响晶体管以及整个电路的特性。在这些无意引入的应力中,一部分不仅仅依赖于产生应力的工艺本身的参数或材料本身的特性,而是与待制作电路的版图形状直接相关,不同的版图结构会最终造成不同的应力分布,这类应力被称为版图相关应力。版图相关应力的典型例子是用于器件隔离的STI(Shallow Trench Isolation)工艺所引入的应力。在STI工艺中,首先在需要隔离的区域通过刻蚀在硅材料上形成隔离槽,然后通过氧化及淀积在隔离槽中填充SiO2,实现对不同区域的隔离。在这一工艺过程中,材料的刻蚀、氧化、淀积、退火以及不同材料的热膨胀系数差异等均会在器件中引入额外的应力。器件的位置、有源区的形状与大小以及周围器件的形状都会影响STI工艺所引入的应力的分布。由于这些过程导致的应力分布与具体的版图几何结构直接相关,最终得到的器件性能也依赖于具体的版图模式。
电子设计自动化(EDA)工具被用于在设计阶段分析预测集成电路的性能。电路仿真工具,例如HSPICE是用于电路分析的EDA工具中的一种。电路仿真工具通过采用一定的器件模型来计算电路网单所描述的集成电路的电学特性。晶体管的器件模型使用一组参数和一定的数学模型来描述特定工艺所制造出的晶体管的特性。为了在电路设计阶段获得这些与版图相关的应力对最终电路性能的影响,以实现对不同来源的应力做出优化来得到最佳的电路性能,需要有准确而易用的模型与分析方法。由于开发对应版图相关应力的模拟模型比较困难,且缺乏有效的仿真方法,尚未在目前应用的电子设计自动化(EDA)工具中包括。在最新版本的BSIM4晶体管模型中仅包括了版图相关应力的部分影响,通过扩散区长度(LOD)参数来表征不同宽度STI结构对载流子迁移率的影响。然而实际的应力分布是三维的,而且受到多种不同因素的影响,这种简化的一维或伪二维模型远不能准确描述版图相关应力的影响。多种其它未考虑的因素,如沟道宽度、多指晶体管、非矩形的有效区域、邻近器件的影响等,均会改变晶体管中的应力分布,对载流子迁移率产生显著的影响。因此,为了得到版图相关应力对电路特性的影响,需要计算晶体管中的三维应力分布状况。
采用有限单元法的应力数值计算能够处理复杂的几何结构并给出准确的结果,通过设置合适的边界条件和初始条件,有限单元法应力数值计算可以用于对晶体管中三维应力分布的计算,进而用于得到晶体管考虑应力后的性能。采用初始应力可以描述多种应力类型,包括材料自身初始应力,材料生长体积膨胀引入的初始应力,材料热膨胀系数差异引入的热失配应力,以及作为对特定工艺校准的初始应力。现有的相关商用仿真工具均采用有限元法来进行应力分布的计算。但现有的此类工具大都为工艺模拟(TCAD)而设计开发,基于预设的工艺参数以及版图,通过持续的网格更新以及方程求解的迭代过程,对集成电路制造工艺中各个工艺步骤的逐步模拟来得到最终器件的特性。但受限于可用的计算机硬件资源,现有的此类模拟工具仅能处理包含数个晶体管的电路单元,且模拟过程非常耗时,无法应用于较大规模电路中版图相关应力的计算。此外,在实际的集成电路设计中,设计者所面对的是作为一个整体实现某一功能的电路,往往包括由数万以至数百万晶体管构成的大量不同功能单元,对单个晶体管进行上述的计算并不能给出版图相关应力对电路整体性能影响的预测。
实际的应用需要能有一套完整有效的方法来准确而高效的实现从电路网单及版图出发,最终给出考虑版图相关应力后的电路性能参数。
发明内容
本发明的目的在于解决现有技术难以解决较大规模电路的版图相关应力对电路特性影响的问题,提出一种获得集成电路中版图相关应力对电路特性影响的方法,本发明通过采用压阻迁移率模型计算版图相关应力造成的载流子迁移率的变化量,并应用于晶体管模型中,最终得到考虑版图相关应力后的电路特性。具有网格更新过程简化,提高了模拟速度和规模,计算结果精确,可以处理复杂的版图结构等到特点。
本发明提出一种获得考虑版图相关应力后电路性能的方法,其特征在于,该方法包括以下步骤:
1)提取出版图中各个晶体管的区域;其具体步骤为:
11)版图划分:在整个版图范围内,检测多晶硅及有源区所对应的版图层的重合区域,得到的每一个独立的重合区域对应一个晶体管沟道区;
12)版图提取:使用版图提取工具对版图及版图所对应的电路网单进行提取处理,得到版图与电路网单中各晶体管的对应关系以及版图中各晶体管沟道区左下角坐标作为提取结果;
13)沟道区与版图提取结果匹配:将步骤11)中得到的各晶体管沟道区与步骤12)中的提取结果相比较,确定出版图中每个晶体管沟道区所对应的电路网单中的晶体管;
14)划分各晶体管有效区域:对版图中每个晶体管划分出有效区域,以将大规模电路划分成小单元;
2)依据版图构造每个晶体管有效区域的三维结构,设定初始条件及边界条件,采用通用的有限单元方法求解,得到每个晶体管有效区域的三维应力分布;具体包括以下步骤:
21) 构造晶体管有效区域的三维结构:根据每个晶体管的有效区域以及该区域的电路版图包含的平面内二维几何结构信息,以及具体的半导体制造工艺步骤包含的厚度信息,构造出晶体管有效区域的三维几何结构;
22)设定初始条件及边界条件:将每个晶体管有效区域三维结构的各个部分内的各种应力来源对应的初始应力值单独计算后线性叠加,作为该部分有限单元法计算的初始条件;在每个晶体管有效区域三维结构的四个侧面边界上采用反对称边界条件;在每个晶体管有效区域的下底面上采用固定边界条件;在有效区域的上表面采用自由边界条件;
23)根据每个晶体管三维结构对应的边界条件及各部分的初始条件,采用通用的有限单元法,计算得到每个晶体管有效区域的三维应力分布作为该晶体管的考虑版图相关应力;
3)计算考虑版图相关应力影响后各个晶体管的迁移率,使用该迁移率更新原有的晶体管模型,使用新的晶体管模型进行计算,得到考虑版图相关应力影响之后的电路特性;其具体步骤为:
31)计算各个晶体管的考虑版图相关应力作用后的迁移率:采用压阻迁移率模型计算得到考虑版图相关应力后的载流子迁移率的变化量;
32)电路性能计算:对每个晶体管,使用考虑版图相关应力后的新迁移率值替换原晶体管模型中的迁移率;根据替换迁移率后的晶体管模型,使用电路仿真工具,对待分析的电路进行电路性能计算,得到考虑版图相关应力后的电路性能。
上述步骤13)中确定出版图中每个晶体管沟道区所对应的电路网单中的晶体管的具体的方法可为:由步骤12)得到的各个晶体管沟道区左下角坐标,对每个版图提取结果中给出的晶体管,寻找步骤11)中得到的沟道区中覆盖其坐标值的区域,该区域即为该坐标信息对应晶体管的匹配沟道区。重复该步骤即实现每个晶体管与其对应沟道区的匹配。
上述步骤14)中对版图中每个晶体管划分出有效区域的具体方法可为,对每个晶体管,取其沟道区域的外接矩形,以该矩形中心为参考点,将该矩形沿x方向的边长向两侧分别延长一定距离Dx,沿y方向的边长向两侧分别延长一定距离Dy,得到一个新的矩形区域,作为该晶体管的有效区域。
上述步骤31)中采用压阻迁移率模型计算得到考虑版图相关应力后的载流子迁移率的变化量的具体步骤可为:根据每个晶体管的版图相关应力的三维分布,取该晶体管的硅材料中上表面内处于沟道区中每个格点上的版图相关应力值,采用压阻迁移率模型计算各格点上由所述版图相关应力造成的迁移率变化量,以整个沟道区中各个格点上的迁移率变化量的平均值作为该晶体管的迁移率变化量;将该晶体管的迁移率变化量叠加到原晶体管模型中的迁移率值上,作为该晶体管的考虑版图相关应力后的新迁移率值。
本发明的特点及有益效果
本发明首先采用划分的方法,在大规模集成电路中,提取出每个晶体管在版图上对应的沟道区域。对于每个晶体管,划分出影响其沟道区应力分布的小区域,仅对该区域进行有限单元法计算以得到对应晶体管的应力。保证了并行处理,提高计算速度,实现了对大规模电路直接的数值模拟,以及直接由电路版图得到考虑版图相关应力后的器件性能变化。
本发明使用初始应力来描述多种应力效应,简化了在工艺模拟工具中的网格更新过程,实现了模拟速度和规模的提高。
本发明可以给出比现有基于解析模型的应力计算精确的结果,可以处理复杂的版图结构。
附图说明
图1为本发明的一个实施例电路-带缓冲的SR触发器电路原理图;
图2为对应图1中实施例电路的版图;
图3为一个半导体制造工艺步骤实施例的厚度信息的示意简图;
图4为对应图1所示电路中一个晶体管的有效区域三维结构示意图;
图5为对应图1所示电路考虑应力后的电路性能计算结果与未考虑应力时的电路性能计算结果比较;
具体实施方式
下面结合附图及实例说明本方法的具体实施方式:
本发明提出的一种获得考虑版图相关应力后电路性能的方法,包括以下步骤:
1)提取出版图中各个晶体管的区域;其具体步骤为:
11)版图划分
在整个版图范围内,检测多晶硅及有源区所对应的版图层的重合区域,得到的每一个独立的重合区域对应一个晶体管沟道区(即为晶体管的核心区域);图1为本发明的一个实施例的电路(带缓冲的SR触发器电路),图中示出该电路包括16个晶体管M1-M16;图2为该实施例电路对应的版图,图中示出16个独立的重合区域对应16个晶体管沟道区M1-M16。图2中只显示了有源区层(图中的点阵区域)与多晶硅层(图中的水平线段区域),二者的重合区域即沟道区。
12)版图提取
使用版图提取工具对版图及版图所对应的电路网单进行提取处理,得到版图与电路网单中各晶体管的对应关系以及版图中各晶体管沟道区左下角坐标作为提取结果;
13)沟道区与版图提取结果匹配
将步骤11)中得到的各晶体管沟道区与步骤12)中的提取结果相比较,确定出版图中每个晶体管沟道区所对应的电路网单中的晶体管;
具体的方法是,由步骤12)得到的各个晶体管沟道区左下角坐标,对每个版图提取结果中给出的晶体管,寻找步骤11)中得到的沟道区中覆盖其坐标值的区域,该区域即为该坐标信息对应晶体管的匹配沟道区。重复该步骤即实现每个晶体管与其对应沟道区的匹配;
14)划分各晶体管有效区域
对版图中每个晶体管划分出有效区域,(即对该晶体管沟道区应力分布有不可忽略贡献的区域),以将大规模电路划分成小单元;
具体方法为,对每个晶体管,取其沟道区域的外接矩形,以该矩形中心为参考点,将该矩形沿x方向的边长向两侧分别延长一定距离Dx,沿y方向的边长向两侧分别延长一定距离Dy,得到一个新的矩形区域,作为该晶体管的有效区域。其中延长距离Dx与Dy是根据具体工艺及精度要求确定(本实施例中设定边界处应力源对沟道区应力影响小于5%,所采用的Dx值为2微米,Dy值为0.5微米),使得边界处应力源对沟道区应力影响小于所需的计算精度。
2)依据版图构造每个晶体管有效区域的三维结构,设定初始条件及边界条件,采用通用的有限单元方法求解,得到每个晶体管有效区域的三维应力分布;具体包括以下步骤:
21)构造晶体管有效区域的三维结构:
根据每个晶体管的有效区域以及该区域的电路版图(平面内二维几何结构信息),以及具体的半导体制造工艺步骤(厚度信息),构造出晶体管有效区域的三维几何结构;图3中给出了一个半导体制造工艺步骤实施例的厚度信息,从硅衬底11开始,制造工艺依次为氧化层12、多晶硅层13和氮化硅层14。图4中给出了对应图1所示电路中一个晶体管的一个有效区域的三维几何结构实例,图中的网格区41为沟道区;
22)设定初始条件及边界条件
将每个晶体管有效区域三维结构的各个部分内的各种应力来源(每个晶体管有效区域三维结构的不同部分包含不同的应力来源)对应的初始应力(所述的初始应力种类包括但不限于材料自身初始应力,材料生长体积膨胀引入的初始应力,材料热膨胀系数差异引入的热失配应力,以及作为对特定工艺校准的初始应力)值单独计算后线性叠加,作为该部分有限单元法计算的初始条件;
在每个晶体管有效区域三维结构的四个侧面边界上采用反对称边界条件(即限制垂直平面方向的位移同时允许平面内的位移);在每个晶体管有效区域的下底面上采用固定边界条件(即限制沿所有方向的位移);在有效区域的上表面采用自由边界条件(即允许沿任意方向的位移);
23)根据每个晶体管三维结构对应的边界条件及各部分的初始条件,采用通用的有限单元法,计算得到每个晶体管有效区域的三维应力分布作为该晶体管的考虑版图相关应力;
3)计算考虑版图相关应力影响后各个晶体管的迁移率,使用该迁移率更新原有的晶体管模型,使用新的晶体管模型进行计算,得到考虑版图相关应力影响之后的电路特性;其具体步骤为:
31)计算各个晶体管的考虑版图相关应力作用后的迁移率
采用压阻迁移率模型(采用C.S.Smith,“Piezoresistance effect ingermanium and silicon,”Physical review,vol.94,no.1,pp.42-49,1954.中公开的方法)计算得到考虑版图相关应力后的载流子迁移率的变化量;具体步骤为:根据每个晶体管的版图相关应力的三维分布,取该晶体管的硅材料中上表面内处于沟道区中每个格点上的版图相关应力值,采用压阻迁移率模型计算各格点上由所述版图相关应力造成的迁移率变化量,以整个沟道区中各个格点上的迁移率变化量的平均值作为该晶体管的迁移率变化量;将该晶体管的迁移率变化量叠加到原晶体管模型中的迁移率值上,作为该晶体管的考虑版图相关应力后的新迁移率值;
32)电路性能计算
对每个晶体管,使用考虑版图相关应力后的新迁移率值替换原晶体管模型中的迁移率;根据替换迁移率后的晶体管模型,使用电路仿真工具,对待分析的电路进行电路性能计算,得到考虑版图相关应力后的电路性能。图5中给出了对本实施例考虑应力后的电路性能计算结果与未考虑应力时的电路性能计算结果的比较。实线与点划线为输入1与输入2端口的电压,长虚线为未包括版图相关应力的输出1端口电压,短虚线为包括版图相关应力后的输出1端口电压,即得到考虑版图相关应力后电路性能。
Claims (4)
1、一种获得考虑版图相关应力后电路性能的方法,其特征在于,该方法包括以下步骤:
1)提取出版图中各个晶体管的区域;其具体步骤为:
11)版图划分:在整个版图范围内,检测多晶硅及有源区所对应的版图层的重合区域,得到的每一个独立的重合区域对应一个晶体管沟道区;
12)版图提取:使用版图提取工具对版图及版图所对应的电路网单进行提取处理,得到版图与电路网单中各晶体管的对应关系以及版图中各晶体管沟道区左下角坐标作为提取结果;
13)沟道区与版图提取结果匹配:将步骤11)中得到的各晶体管沟道区与步骤12)中的提取结果相比较,确定出版图中每个晶体管沟道区所对应的电路网单中的晶体管;
14)划分各晶体管有效区域:对版图中每个晶体管划分出有效区域,以将大规模电路划分成小单元;
2)依据版图构造每个晶体管有效区域的三维结构,设定初始条件及边界条件,采用通用的有限单元方法求解,得到每个晶体管有效区域的三维应力分布;具体包括以下步骤:
21)构造晶体管有效区域的三维结构:根据每个晶体管的有效区域以及该区域的电路版图包含的平面内二维几何结构信息,以及具体的半导体制造工艺步骤包含的厚度信息,构造出晶体管有效区域的三维几何结构;
22)设定初始条件及边界条件:将每个晶体管有效区域三维结构的各个部分内的各种应力来源对应的初始应力值单独计算后线性叠加,作为该部分有限单元法计算的初始条件;在每个晶体管有效区域三维结构的四个侧面边界上采用反对称边界条件;在每个晶体管有效区域的下底面上采用固定边界条件;在有效区域的上表面采用自由边界条件;
23)根据每个晶体管三维结构对应的边界条件及各部分的初始条件,采用通用的有限单元法,计算得到每个晶体管有效区域的三维应力分布作为该晶体管的考虑版图相关应力;
3)计算考虑版图相关应力影响后各个晶体管的迁移率,使用该迁移率更新原有的晶体管模型,使用新的晶体管模型进行计算,得到考虑版图相关应力影响之后的电路特性;其具体步骤为:
31)计算各个晶体管的考虑版图相关应力作用后的迁移率:采用压阻迁移率模型计算得到考虑版图相关应力后的载流子迁移率的变化量;
32)电路性能计算:对每个晶体管,使用考虑版图相关应力后的新迁移率值替换原晶体管模型中的迁移率;根据替换迁移率后的晶体管模型,使用电路仿真工具,对待分析的电路进行电路性能计算,得到考虑版图相关应力后的电路性能。
2、如权利要求1所述方法,其特征在于,所述步骤13)中确定出版图中每个晶体管沟道区所对应的电路网单中的晶体管的具体的方法是:由步骤12)得到的各个晶体管沟道区左下角坐标,对每个版图提取结果中给出的晶体管,寻找步骤11)中得到的沟道区中覆盖其坐标值的区域,该区域即为该坐标信息对应晶体管的匹配沟道区。重复该步骤即实现每个晶体管与其对应沟道区的匹配。
3、如权利要求1所述方法,其特征在于,所述步骤14)中对版图中每个晶体管划分出有效区域的具体方法为,对每个晶体管,取其沟道区域的外接矩形,以该矩形中心为参考点,将该矩形沿x方向的边长向两侧分别延长一定距离Dx,沿y方向的边长向两侧分别延长一定距离Dy,得到一个新的矩形区域,作为该晶体管的有效区域。
4、如权利要求1所述方法,其特征在于,所述步骤31)中采用压阻迁移率模型计算得到考虑版图相关应力后的载流子迁移率的变化量的具体步骤为:根据每个晶体管的版图相关应力的三维分布,取该晶体管的硅材料中上表面内处于沟道区中每个格点上的版图相关应力值,采用压阻迁移率模型计算各格点上由所述版图相关应力造成的迁移率变化量,以整个沟道区中各个格点上的迁移率变化量的平均值作为该晶体管的迁移率变化量;将该晶体管的迁移率变化量叠加到原晶体管模型中的迁移率值上,作为该晶体管的考虑版图相关应力后的新迁移率值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100913574A CN101630339B (zh) | 2009-08-21 | 2009-08-21 | 一种获得考虑版图相关应力后电路性能的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100913574A CN101630339B (zh) | 2009-08-21 | 2009-08-21 | 一种获得考虑版图相关应力后电路性能的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101630339A true CN101630339A (zh) | 2010-01-20 |
CN101630339B CN101630339B (zh) | 2011-08-10 |
Family
ID=41575445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100913574A Expired - Fee Related CN101630339B (zh) | 2009-08-21 | 2009-08-21 | 一种获得考虑版图相关应力后电路性能的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101630339B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101789049A (zh) * | 2010-03-04 | 2010-07-28 | 苏州锐调科技有限公司 | 集成电路器件匹配的构图方法 |
CN101840458A (zh) * | 2010-05-20 | 2010-09-22 | 清华大学 | 载流子迁移率的提取方法 |
CN102663160A (zh) * | 2012-03-16 | 2012-09-12 | 苏州芯禾电子科技有限公司 | 一种构建三维物理电路设计模型的方法 |
CN103838888B (zh) * | 2012-11-23 | 2017-03-29 | 中芯国际集成电路制造(上海)有限公司 | 集成电路spice模型中的应力影响模型的优化 |
CN114021505A (zh) * | 2022-01-06 | 2022-02-08 | 青岛展诚科技有限公司 | 集成电路FinFET复杂三维结构描述文件的生成方法和系统 |
-
2009
- 2009-08-21 CN CN2009100913574A patent/CN101630339B/zh not_active Expired - Fee Related
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101789049A (zh) * | 2010-03-04 | 2010-07-28 | 苏州锐调科技有限公司 | 集成电路器件匹配的构图方法 |
CN101789049B (zh) * | 2010-03-04 | 2013-08-14 | 何伟 | 集成电路器件匹配的构图方法 |
CN101840458A (zh) * | 2010-05-20 | 2010-09-22 | 清华大学 | 载流子迁移率的提取方法 |
CN102663160A (zh) * | 2012-03-16 | 2012-09-12 | 苏州芯禾电子科技有限公司 | 一种构建三维物理电路设计模型的方法 |
CN103838888B (zh) * | 2012-11-23 | 2017-03-29 | 中芯国际集成电路制造(上海)有限公司 | 集成电路spice模型中的应力影响模型的优化 |
CN114021505A (zh) * | 2022-01-06 | 2022-02-08 | 青岛展诚科技有限公司 | 集成电路FinFET复杂三维结构描述文件的生成方法和系统 |
CN114021505B (zh) * | 2022-01-06 | 2022-05-17 | 青岛展诚科技有限公司 | 集成电路FinFET复杂三维结构描述文件的生成方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101630339B (zh) | 2011-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101630339B (zh) | 一种获得考虑版图相关应力后电路性能的方法 | |
JP3670553B2 (ja) | 半導体集積回路解析装置、半導体集積回路解析方法及び半導体集積回路解析方法を実行するためのプログラムを記録した記録媒体 | |
US8887116B2 (en) | Flexible pattern-oriented 3D profile for advanced process nodes | |
US8615728B2 (en) | Analysis of stress impact on transistor performance | |
TWI476868B (zh) | 用於放置並安排路由系統中之設計最佳化的填充單元 | |
US9311443B2 (en) | Correcting for stress induced pattern shifts in semiconductor manufacturing | |
JP4973572B2 (ja) | 半導体回路の設計方法および製造方法 | |
KR101573130B1 (ko) | 고밀도 집적 회로 디바이스들을 제조하는 방법 | |
CN102184879A (zh) | 一种soi场效应晶体管的tcad仿真校准方法 | |
JP2009032199A (ja) | シミュレーション方法及びプログラム | |
US20100082308A1 (en) | Circuit simulation based on gate spacing from adjacent MOS transistors | |
US9378328B2 (en) | Modeling multi-patterning variability with statistical timing | |
CN104750904B (zh) | 用以改进晶体管匹配的方法 | |
US20080216041A1 (en) | Integrated circuit simulation method considering stress effects | |
CN103838888B (zh) | 集成电路spice模型中的应力影响模型的优化 | |
TWI792420B (zh) | 互補式場效電晶體中的源極和汲極材料的磊晶生長 | |
CN108875200B (zh) | 通用型wpe优化模型及其提取方法 | |
Joshi et al. | Closed-form modeling of layout-dependent mechanical stress | |
US9836569B2 (en) | Leakage reduction using stress-enhancing filler cells | |
Verghese et al. | Modeling Stress-Induced Variability Optimizes IC Timing Performance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110810 Termination date: 20150821 |
|
EXPY | Termination of patent right or utility model |