Sentaurus非局域隧穿模型：从理论到FTJ仿真的关键配置解析

1. 非局域隧穿模型为何成为FTJ仿真的关键选择在模拟铁电隧穿结FTJ这类新型存储器时传统隧穿模型常常力不从心。我刚开始接触FTJ仿真时发现用常规方法模拟出的电流-电压曲线总是和实验数据对不上直到尝试了Sentaurus中的非局域隧穿模型Nonlocal Tunneling Model问题才迎刃而解。这个模型的独特之处在于它能同时处理三种关键载流子行为来自导带的电子ECB、来自价带的空穴HVB以及价带电子直接隧穿EVB。举个例子当我们研究HfO₂基FTJ器件时传统局域隧穿模型只能简单考虑导带电子隧穿而实际上在铁电材料中价带电子通过缺陷态参与的隧穿过程往往占总电流的30%以上。非局域模型通过Band2Band参数的不同配置模式Simple/Full/UpsideDown可以灵活组合这些隧穿分量。实测数据显示采用Full模式仿真Zr掺杂HfO₂器件时仿真结果与实测I-V曲线的误差能从原来的50%降低到8%以内。2. Physics模块配置激活隧穿机制的神经中枢2.1 基础配置模式对比在sdevice文件中Physics模块就像控制隧穿机制的开关面板。我常用的五种配置方式各有适用场景# 最简配置仅ECB和HVB Physics { eBarrierTunneling NLM hBarrierTunneling NLM } # 带势垒降低效应的进阶配置 Physics { eBarrierTunneling NLM(BarrierLowering) } # 包含EVB分量的典型配置 Physics { eBarrierTunneling NLM(Band2BandSimple) hBarrierTunneling NLM(Band2BandSimple) }对于10nm以下厚度的FTJ我强烈推荐使用Full模式。去年在仿真Al₂O₃/Si₃N₄复合势垒结构时Simple模式漏掉了关键的缺陷辅助隧穿路径导致阈值电压预测偏差达0.3V。改用Full模式后不仅重现了实验观察到的双峰隧穿特征还准确预测了温度系数。2.2 Band2Band模式的物理内涵Simple模式相当于在ECB基础上增加镜像的EVB分量适合单晶势垒材料Full模式会计算所有可能的k空间路径对多晶铁电薄膜特别重要UpsideDown模式专门处理反常能带结构在钙钛矿型FTJ中效果显著有个容易踩的坑是忘记匹配电子和空穴的Band2Band模式。有次仿真中只设置了电子的Full模式而空穴用默认值结果空穴电流被严重低估。正确的做法是保持对称配置Physics { eBarrierTunneling NLM(Band2BandFull) hBarrierTunneling NLM(Band2BandFull) }3. Math模块构建隧穿路径的精密地图3.1 非局域网格的两种定义方式Math模块中的NonLocal定义就像绘制隧穿载流子的高速公路网。经过多次测试我发现对于FTJ这类具有明确界面结构的器件MaterialInterface比RegionInterface更可靠Math { NonLocal NLM_HfO2( MaterialInterfaceHfO2/Si Length 3e-7 # 关键参数需大于势垒厚度 Digits 6 # 高精度需求时建议≥5 EnergyResolution 0.0005 ) }最近在仿真Zr₀.₅Hf₀.₅O₂器件时发现当铁电层厚度小于5nm时必须将Length参数设置为实际厚度的1.5-2倍否则会漏掉重要的边缘隧穿效应。这个经验也写进了我们团队的内部设计手册。3.2 精度参数的设置艺术EnergyResolution参数对计算效率影响巨大。通过对比测试发现设置值计算时间电流误差0.01eV1x15%0.001eV3x5%0.0001eV8x1%对于初期参数扫描建议先用0.01eV快速筛选最终仿真时再用0.001eV。除非研究量子相干效应否则不需要用到0.0001eV。4. Parameter模块微调隧穿行为的秘密武器4.1 全局与材料局部的参数覆盖Parameter文件中可以定义多级参数这在处理复合势垒时特别有用。比如下面这个配置BarrierTunneling NLM{ g 1.0, 1.0 # 通用校正因子 mt 0.5, 0.8 # 电子/空穴有效质量比 } MaterialHfO2 { BarrierTunneling NLM{ mt 0.3, 0.6 # 覆盖全局设置 A 2.0 # 材料特定拟合系数 } }有个实用技巧先用全局参数进行初步仿真再针对特定材料微调。去年优化La:ZrO₂器件模型时通过单独调整HfO2层的mt值从0.46降到0.38使亚阈值摆幅的仿真误差从12%降到3%。4.2 关键参数的经验取值根据多个项目积累的数据总结出这些典型取值范围g因子通常1.0-1.2铁电相时建议1.05mt(电子)HfO₂系0.3-0.5BST系0.4-0.6A系数与缺陷密度相关1.0-3.0不等特别注意当仿真出现收敛困难时可以尝试先将所有g值设为1.0mt取中间值0.5等收敛后再逐步调整到真实值。这个方法帮我解决了至少三次棘手的收敛问题。