Halo注入MOSFET闪烁噪声建模：非均匀沟道下的1/f噪声分析与BSIM6集成

1. 项目概述与核心挑战在模拟和射频集成电路设计中闪烁噪声也就是我们常说的1/f噪声一直是个让人头疼的“老朋友”。它就像电路背景里挥之不去的底噪频率越低能量越强直接影响着运算放大器、压控振荡器、锁相环等关键电路的性能尤其是低频下的信噪比和相位噪声。对于从事模拟IC设计或者器件建模的工程师来说准确预测和建模这种噪声是确保芯片一次流片成功、性能达标的基础。传统的闪烁噪声模型比如经典的统一噪声模型Unified Noise Model其物理基础是假设MOSFET的沟道是均匀掺杂的。这个假设在过去工艺节点较大时基本成立。然而随着工艺尺寸不断微缩到45纳米及以下为了抑制恼人的短沟道效应比如漏致势垒降低DIBLHalo注入或称口袋注入工艺成为了标准配置。简单说就是在源漏结的侧向注入高浓度的杂质形成局部的高掺杂“口袋”像给沟道两端加了个“护城河”有效控制了关态漏电。但凡事有利有弊。Halo注入在提升数字电路性能的同时却给模拟性能埋下了“地雷”。这个高掺杂区域会引入额外的晶格损伤和界面态导致局部陷阱密度激增。更重要的是它造成了沟道内阈值电压Vth的非均匀分布Halo区域Vth高中心沟道区域Vth低。这种“地形”上的剧烈变化使得载流子在沟道中的分布和输运变得异常复杂传统的均匀沟道噪声模型立刻“失灵”了。我遇到过最典型的场景就是一个用于传感器信号调理的低噪声放大器在长沟道器件上仿真结果很漂亮但实际芯片测试时在亚阈值区或近阈值电压区的噪声性能远差于预期。排查半天最后问题就出在模型没有准确反映Halo注入带来的噪声特性畸变。具体表现就是在弱反型区无论沟道长短噪声几乎不再随长度缩放而变化而对于长沟道器件其噪声功率谱密度随栅压的变化曲线会出现反常的“拐点”。这些现象用老模型完全无法解释和拟合。因此开发一个能够兼容主流紧凑模型如BSIM6、且能精确刻画Halo注入MOSFET闪烁噪声特性的新模型就成了一个既紧迫又具有高度实用价值的课题。这不仅是一个学术问题更是支撑先进工艺节点下高性能模拟电路设计的基石。2. 模型核心思路化整为零的串联晶体管法面对Halo注入带来的非均匀沟道最直观的建模思路就是“分而治之”。既然整个沟道在物理和电学特性上可以清晰地分为两个区域——两端的高掺杂Halo区域和中间的低掺杂本征沟道区域那么很自然地我们可以将整个MOSFET想象成两个虚拟的MOSFET串联而成。2.1 物理图像的建立我们可以把实际器件总长L拆解为Halo晶体管对应源、漏两端的Halo注入区域。假设每个Halo区域的长度为Lh通常对称总Halo区域长度为2Lh但在许多简化模型中特别是对于噪声建模常将其等效为一个长度为Lh的晶体管来代表总的Halo效应其掺杂浓度高Nh阈值电压高Vth,h且由于注入损伤其氧化层界面处的等效陷阱密度也更高。沟道晶体管对应中间的本征沟道区域。长度为 L - Lh或 L - 2Lh掺杂浓度低Nch阈值电压低Vth,ch。这两个虚拟晶体管串联流过相同的漏极电流Ids。这个等效电路是理解后续所有噪声推导的基石。它的巧妙之处在于将复杂的二维非均匀问题简化成了两个相对均匀、可独立分析的“子器件”的互联问题。2.2 噪声贡献因子谁是噪声输出的“放大器”关键问题来了Halo区域和沟道区域自身都会产生闪烁噪声但它们的噪声“传递”到外部测量端漏极的效率一样吗答案是否定的。这就引入了本模型最核心的概念之一噪声贡献因子。想象一下Halo区域和沟道区域就像两个串联的电阻它们内部都有自己的噪声源类似一个噪声电压源。这个内部噪声要想在总电流中体现出来必须经过所在区域以及另一个区域的“阻抗分压”。贡献因子CF, Contribution Factor本质上就是这个“分压比”或“传递函数”。通过小信号电路分析具体推导过程基于叠加原理详见原论文图3我们可以得到总漏极电流噪声功率谱密度 S_ID 的表达式S_ID S_ID,h * CF_h S_ID,ch * CF_ch其中S_ID,h和S_ID,ch分别是假设Halo晶体管和沟道晶体管独立工作时它们自身产生的漏极电流噪声功率谱密度。这部分由各区域的局部物理参数陷阱密度、迁移率、反型电荷密度等决定。CF_h和CF_ch就是贡献因子它们的表达式为CF_h (gm_ch gd_ch) / (gm_ch gd_ch gd_h)CF_ch gd_h / (gm_ch gd_ch gd_h)这里gm_ch是沟道晶体管的跨导gd_ch和gd_h分别是沟道和Halo晶体管的输出电导。注意这个贡献因子的形式非常关键。它不同于早期一些模型简单地将两个区域的噪声功率直接相加即假设CF_h CF_ch 1。那种做法会严重高估强反型区的噪声因为此时Halo区域虽然自身噪声可能不小但其噪声传递到外部的效率CF_h已经变得非常低。2.3 贡献因子的物理意义与偏置依赖性贡献因子的行为深刻揭示了噪声的偏置依赖特性在弱反型区此时沟道晶体管处于弱反型或耗尽其跨导gm_ch很小。而Halo晶体管由于阈值电压更高几乎完全关闭其输出电导gd_h极小远小于gm_ch。根据公式CF_h ≈ 1CF_ch ≈ 0。这意味着总噪声几乎完全由Halo区域的噪声S_ID,h主导。这也完美解释了图1中的现象在弱反型区噪声与沟道长度L无关因为决定噪声的是Halo区域的特性Lh而Lh通常不随设计沟道长度L等比例变化。在强反型区沟道晶体管完全开启gm_ch很大。同时随着栅压升高Halo区域也逐渐反型gd_h增大。此时CF_h迅速减小CF_ch增大。在长沟道器件中L远大于Lh沟道区域的电阻占主导导致CF_ch CF_h因此总噪声转而由沟道区域的噪声S_ID,ch主导。短沟道器件当L很小时Lh与(L-Lh)可比两个区域的“权重”相差不大。因此从弱反型到强反型的过渡非常平滑不会出现长沟道器件那样明显的噪声主导权“交接”拐点。这个基于贡献因子的框架是模型能够精确捕捉复杂偏置依赖性的灵魂。它告诉我们评估一个区域的噪声影响不仅要看它“喊”得多大声S_ID, local更要看它的声音有多少能“传出来”CF。3. 模型实现细节与BSIM6兼容性理论框架搭建好了下一步就是如何将其落地变成一个可供电路仿真器如SPICE调用的、与工业标准BSIM6模型无缝兼容的紧凑模型。这是从论文走向实践的关键一步。3.1 关键参数的计算要实现模型需要计算几个关键的中介变量主要是两个虚拟晶体管的源端和漏端归一化反型电荷密度qs, qd。从全局到局部首先通过标准的BSIM6模型或任何其他精确的DC模型我们可以得到整个器件在给定Vgs和Vds下的漏极电流Ids、表面势等全局信息。求解电荷密度对于串联的两个晶体管它们流过相同的Ids。利用电荷控制模型例如基于EKV或BSIM6的核心方程我们可以为每个虚拟晶体管建立电流方程对于Halo晶体管Ids f(qs_h, qd_h, Vth_h, Lh, ...)对于沟道晶体管Ids f(qs_ch, qd_ch, Vth_ch, L-Lh, ...) 其中整个器件的源端电荷qs_h是已知的由全局Vbs决定整个器件的漏端电荷qd_ch也是已知的由全局Vds决定。而两个晶体管连接处的电荷qd_h和qs_ch是未知的内部节点。连接点求解利用电流连续性原理串联电流相等我们可以建立方程来求解这两个内部节点的电荷qd_h和qs_ch。这通常需要求解一个关于归一化电荷的二次方程如原论文中的公式(16)和(17)。这个过程确保了在两个虚拟晶体管的连接处电流和电荷都是连续的这是物理正确性的基础。计算电导参数得到所有四个电荷qs_h, qd_h, qs_ch, qd_ch后就可以根据电荷计算各虚拟晶体管的跨导gm_ch和输出电导gd_ch, gd_h。这些参数直接用于计算前述的贡献因子CF_h和CF_ch。3.2 噪声源模型的选择与集成对于每个虚拟晶体管自身的噪声功率S_ID,h和S_ID,ch本模型采用了业界广泛接受的统一闪烁噪声模型。该模型同时考虑了载流子数波动和迁移率波动两种机制其表达式一般形式为S_ID,local (Ids^2 / (f * W * L_local^2)) * ∫(N_t(x) / N_inv(x)^2) dx*其中f是频率。W是栅宽。L_local是虚拟晶体管的长度Lh或L-Lh。N_t*(x)是等效的氧化层陷阱面密度它通常是局部反型层载流子密度N_inv(x)的函数可以建模为N_t* A BN_inv CN_inv^2。A、B、C是拟合参数分别代表与载流子密度无关的陷阱项、线性相关项和平方相关项。N_inv(x)是反型层载流子面密度与归一化电荷q(x)直接相关。这个积分可以通过将变量从位置x转换为归一化电荷q来解析求解最终得到S_ID,local关于qs和qd的闭合表达式。这里的关键在于Halo晶体管和沟道晶体管使用不同的参数组A_h, B_h, C_h和A_ch, B_ch, C_ch。这直接反映了Halo区域因离子注入导致陷阱密度增加这一物理事实。通常A_h A_ch。3.3 与BSIM6的集成方式BSIM6作为最新的体硅MOSFET工业标准紧凑模型其架构允许以“模块化”的方式集成新的物理效应模型。本噪声模型的集成通常通过以下方式实现参数扩展在BSIM6模型的参数集中新增用于描述Halo噪声的参数。例如LHALO等效Halo区域长度。NOIA_H,NOIB_H,NOIC_HHalo区域的闪烁噪声参数A_h, B_h, C_h。NOIA_CH,NOIB_CH,NOIC_CH沟道区域的闪烁噪声参数可能与标准BSIM6的NOIA, NOIB, NOIC不同或复用。可能还需要一个标志参数HALO_NOISE来控制该模型的开启与关闭。代码实现在BSIM6的模型C代码中噪声计算函数如evalNoise需要被修改或扩展。在计算闪烁噪声时首先判断是否启用了Halo噪声模型。如果启用则执行以下流程 a. 调用标准BSIM6核心函数计算全局电流、电荷、电势。 b. 根据LHALO等参数将沟道分割为Halo和沟道区域。 c. 调用内部子函数计算两个虚拟晶体管的端点电荷(qs_h, qd_h, qs_ch, qd_ch)。 d. 分别计算S_ID,h和S_ID,ch。 e. 计算gm_ch, gd_ch, gd_h进而得到CF_h和CF_ch。 f. 根据公式 S_ID S_ID,h * CF_h S_ID,ch * CF_ch 计算总噪声功率谱密度。模型提取流程对于工艺工程师提取这套新参数需要额外的步骤DC和CV参数首先用标准流程提取BSIM6的DC和CV参数。Halo区域参数估算LHALO通常可以从工艺TCAD仿真或特殊的测试结构如不同偏置下的栅电容曲线中初步估算。噪声参数分步提取 i. 测量长沟道器件在强反型区的噪声。此时噪声由沟道区域主导CF_ch ≈1可以提取出沟道区域的噪声参数NOIA_CH, NOIB_CH, NOIC_CH。这些参数可能与均匀掺杂器件提取的标准噪声参数接近。 ii. 测量长沟道器件在深弱反型区的噪声。此时噪声由Halo区域主导CF_h ≈1结合上一步得到的沟道参数可以拟合提取出Halo区域的噪声参数NOIA_H, NOIB_H, NOIC_H。通常NOIA_H会显著大于NOIA_CH。 iii. 用所有参数同时拟合不同长度长、短沟道和不同偏置从弱反型到强反型的噪声数据进行最终优化和验证。实操心得在模型提取时最容易出错的是低估了Halo区域长度LHALO的影响。它不是一个可随意调节的拟合参数而应具有物理意义。建议通过TCAD仿真或分析不同栅压下的C-V曲线拐点先对其有一个物理预估并将其作为提取过程中的一个软约束这样可以提高参数组的物理可信度和外推性。4. 模型验证、结果分析与设计启示任何模型的价值都需要通过实验数据来验证。原论文基于45纳米低功耗CMOS工艺的测试数据充分展示了本模型相较于传统模型的优越性。4.1 模型与实测数据的对比长沟道器件的复杂偏置依赖性如图1和图4所示对于长沟道器件例如L1um其归一化噪声功率S_ID/W 随漏电流Ids的变化曲线并非单调。在弱反型区低Ids噪声水平很高且曲线斜率较大随着Ids增加进入中反型区曲线会出现一个明显的“拐点”斜率突然下降进入强反型区后斜率再次趋于平缓。传统的均匀模型如BSIM4的默认噪声模型完全无法拟合这种“S”形的变化趋势。而本模型通过Halo和沟道区域噪声主导权的动态转换精确地复现了整个曲线。沟道长度缩放效应的缺失图1另一个关键结果是在弱反型区短沟道如L0.1um和长沟道器件的噪声几乎重合。这与传统认知“噪声功率反比于沟道面积W*L”相悖。本模型给出了完美解释在弱反型区噪声由Halo区域主导。而Halo区域的长度Lh主要由工艺注入条件决定与设计沟道长度L关系不大。因此当L变化时只要Lh不变弱反型噪声就基本不变。这是Halo注入工艺下噪声特性的一个标志性现象。短沟道器件的平滑过渡短沟道器件由于L与Lh可比两个区域的贡献因子变化平缓因此其噪声曲线从弱反型到强反型是平滑过渡的没有长沟道那样的明显拐点。模型同样准确地捕捉到了这一差异。漏压依赖性如图6所示模型在不同漏源电压Vds下也表现稳健。在高栅压强反型下总噪声主要由沟道区域贡献对Vds变化相对不敏感在低栅压弱反型下总噪声由Halo区域贡献其特性也得到了正确反映。4.2 对电路设计的核心启示这个模型不仅是一个更精确的仿真工具其物理内涵更能给模拟电路设计师带来直接的设计指导谨慎使用长沟道器件于低偏置电流应用如果您的电路如超低功耗传感器前端工作在亚阈值或近阈值区传统观念认为使用长沟道器件可以降低1/f噪声。但在Halo注入工艺下这个优势可能不复存在因为弱反型噪声由Halo区域决定与沟道长度无关。此时选择长沟道可能只会增加面积和电容而无法换来噪声性能的提升。需要通过仿真使用正确的模型或实测来验证。优化偏置点对于给定的工艺和器件尺寸噪声特性随偏置剧烈变化。通过仿真可以找到一个噪声相对较低的“甜蜜点”。例如对于某些长沟道器件避开中反型区的那个“噪声拐点”区域可能对优化电路的整体噪声性能有益。理解匹配性能的限制晶体管的匹配特性也与1/f噪声密切相关。Halo注入导致的阈值电压非均匀性以及由此带来的噪声非均匀性可能会影响精密匹配电路如电流镜、差分对的匹配精度尤其是在低电流下。本模型有助于更准确地预测匹配性能。工艺选择的考量对于极度追求噪声性能的模拟模块如VCO核心、高精度ADC的输入级在与工艺厂沟通时可以关注Halo注入的剂量和能量。更轻的Halo注入在满足数字电路泄漏要求的前提下可能有助于改善模拟器件的噪声特性。本模型可以为这种工艺-设计协同优化提供量化分析工具。5. 常见问题、模型局限与扩展讨论在实际应用这个模型或理解相关现象时可能会遇到一些疑问和挑战。5.1 常见问题与误区Q这个模型如此复杂是否值得在所有的仿真中都启用A并非必要。对于纯数字电路或对噪声不敏感的模拟模块使用标准均匀噪声模型可能已足够且仿真速度更快。但对于高性能模拟模块尤其是工作频率较低、或偏置点在亚阈值/近阈值区的电路启用Halo噪声模型对于准确预测相位噪声、信噪比、失调电压等至关重要是避免设计失误的关键。Q模型中假设只有两个区域一个Halo一个沟道但实际的Halo注入在源漏两端各有一个是否不够精确A论文中指出使用两个晶体管一个代表所有Halo区域一个代表沟道的等效已经足以捕捉主要的偏置和长度依赖特性。使用三个晶体管源端Halo、沟道、漏端Halo的模型虽然物理上更精确但会极大增加模型复杂度和参数提取难度而精度提升有限。在工程上两晶体管模型在精度和复杂度之间取得了很好的平衡。Q如何获取模型中的关键参数如等效Halo长度LhALh是一个关键的物理参数而非纯粹的拟合参数。最佳获取途径是通过工艺的TCAD仿真。如果无法获得可以通过分析不同沟道长度器件在特定偏置下的电学参数如阈值电压、DIBL随长度的变化来间接提取。在参数提取流程中应优先从物理层面约束Lh的范围。Q这个模型只适用于45nm工艺吗A不是。其物理机制具有普适性。只要工艺中使用了Halo或类似的非均匀掺杂技术如逆向掺杂Retrograde Well就会引入类似的效应。该模型已在更先进的工艺节点如28nm, 16nm FinFET工艺前的平面器件中得到应用和验证。对于FinFET由于其三维结构和不同的掺杂剖面噪声机制有所不同需要专门的模型但分析非均匀性影响的思路是相通的。5.2 模型的局限性与未来方向温度依赖性当前模型主要关注室温下的噪声特性。陷阱的激活能、载流子迁移率等都是温度的函数因此模型的噪声参数A, B, C可能具有温度依赖性。一个更完善的模型需要包含温度缩放公式。栅压依赖性陷阱密度模型中的等效陷阱密度N_t*是反型电荷密度即栅压的函数。但实际的陷阱能级在禁带中的分布是固定的其与费米能级的相对位置随栅压变化导致被占据的陷阱数量发生变化。更物理的模型会引入基于能级的陷阱分布函数。几何尺寸的精细缩放模型中的Lh被当作一个固定值。但在实际工艺中当沟道长度非常短时源漏的Halo区域可能会发生重叠Lh与L的关系不再是简单的相加。对于纳米尺度器件可能需要考虑这种几何上的非线性效应。与随机电报噪声的关联1/f噪声在微观上被认为是多个随机电报噪声RTN的叠加。在先进工艺节点单个陷阱引起的RTN效应显著特别是在小面积器件中。未来的模型可能需要在一个框架下统一描述1/f噪声和RTN。5.3 给模拟设计工程师的实操建议仿真设置在启动任何关键的噪声仿真前先与负责PDK的团队确认所用的工艺模型是否已经包含了类似的Halo噪声模型可能是BSIM6的一个特定版本或厂商的私有模型。如果包含确保在仿真器如Spectre, HSPICE中正确启用了相应的模型选项和参数。模型卡检查仔细阅读模型文档查看是否有NOIA_H、LHALO等新增参数。对比标准模型卡和包含Halo噪声的模型卡理解参数差异。基准测试针对你的设计中最关键的器件比如差分对输入管、电流镜负载管单独搭建测试电路扫描其从亚阈值到强反型的噪声曲线。将仿真结果与工艺厂提供的测试数据如果有或文献中的典型趋势进行对比验证模型在你关心的偏置区间是否表现合理。敏感性分析由于Halo参数特别是Lh和NOIA_H可能存在工艺波动进行蒙特卡洛分析时应考虑这些参数的统计变化对电路整体噪声性能的影响。这能帮助你评估设计的鲁棒性。这个针对Halo注入MOSFET的闪烁噪声模型从一个具体的工艺挑战出发通过清晰的物理图像和严谨的数学推导提供了一个既精确又实用的解决方案。它深刻地揭示了现代非均匀掺杂器件中噪声行为的复杂性将“沟道非均匀性”这个因素从幕后推到了台前。对于深耕于模拟电路设计的工程师而言理解并善用这类模型意味着能更精准地驾驭先进工艺在芯片设计这场微观世界的博弈中更好地控制那看似随机的“底噪”从而创造出性能更卓越的电路。