硅基SPAD仿真避坑指南ATLAS中那些容易忽略的物理模型参数设置细节在半导体器件仿真领域硅基单光子雪崩二极管SPAD因其超高灵敏度而备受关注。然而许多工程师和研究人员在使用ATLAS进行SPAD仿真时常常在物理模型参数设置环节遭遇意想不到的障碍。本文将从一个过来人的视角揭示那些容易被忽视却至关重要的参数设置细节帮助您避开常见的仿真陷阱。1. 迁移率模型的选择与参数陷阱迁移率模型的选择直接影响载流子输运行为的仿真精度。新手常犯的错误是直接采用默认设置而忽略了不同模型适用的场景差异。1.1 低场迁移率模型对比ATLAS提供了五种低场迁移率模型每种都有其独特的适用场景模型类型关键特征适用场景典型误差来源CONSTANT固定迁移率值快速验证场景忽略掺杂浓度依赖性CONMOB基于查找表的浓度依赖模型常规硅器件未考虑温度变化影响KLAASSEN统一散射模型高精度需求计算资源消耗大ANALYTIC解析表达式中等精度需求高掺杂区域精度不足CCSMOB载流子-载流子散射考虑高注入条件低估强电场下迁移率提示在SPAD仿真中雪崩区域的载流子行为尤为关键。建议至少使用KLAASSEN模型或CONMOBANALYTIC组合以获得更准确的迁移率描述。1.2 高场迁移率与速度饱和当电场强度超过1e4 V/cm时载流子速度饱和效应变得显著。ATLAS中相关参数设置常被忽视mobility conmob analytic fldmob mobility mun1500 mup450 # 基础迁移率值 mobility vsatn1.0e7 vsatp8.0e6 # 饱和速度单位cm/s注意vsatn和vsatp的默认值可能不适用于所有工艺节点。对于先进工艺的SPAD仿真建议通过实验数据校准这些参数。2. 复合模型的关键参数设置复合过程直接影响SPAD的暗计数率和光子探测效率。以下是三个最易出错的复合模型参数2.1 SRH复合的浓度依赖性标准SRH模型假设恒定的载流子寿命而实际上寿命与掺杂浓度密切相关models srh consrh material silicon taun01e-6 taup01e-6 # 基础寿命值(s) material silicon nsrhn1e16 nsrhp1e16 # 参考浓度(cm^-3)常见错误包括未激活CONSRH参数导致忽略浓度依赖性使用默认的taun0/taup0值与实际材料不符未正确设置nsrhn/nsrhp导致浓度缩放不准确2.2 Auger复合的陷阱辅助效应在高注入条件下陷阱辅助Auger复合变得重要models trap.auger material silicon augn2.8e-31 augp9.9e-32 # Auger系数(cm^6/s)实际经验在仿真SPAD的雪崩过程时Auger系数每偏差10%可能导致暗电流仿真结果偏差达30%。2.3 表面复合速度的设置表面复合对SPAD的边缘击穿特性影响显著interface s.n1e3 s.p1e3 # 表面复合速度(cm/s)警告粗糙的表面处理可能导致实际复合速度比仿真设置高1-2个数量级。建议通过TCAD实验确定最佳值。3. 碰撞电离模型的精细调整碰撞电离模型是SPAD仿真的核心也是错误最多的部分。3.1 Selberherr与Crowell-Sze模型对比两种主流模型的适用场景对比特征Selberherr模型Crowell-Sze模型理论基础经验公式物理机制驱动计算效率较高较低高场精度3e5 V/cm时可能偏差全电场范围较准确温度依赖性需手动设置多组参数内置更完善的温度模型适用场景快速仿真、初步设计高精度需求、科学研究3.2 关键参数校准技巧以Selberherr模型为例关键参数设置impact selb an7.03e5 bn1.23e6 # 电子电离参数 impact selb ap1.58e6 bp2.04e6 # 空穴电离参数校准建议从文献获取初始值在1e5-5e5 V/cm电场范围内与实验数据对比采用二分法逐步调整参数注意不同温度下的参数变化踩坑记录曾有一个项目因未校准bn参数导致击穿电压仿真结果偏差达15V耗费两周排查。4. 收敛性问题诊断与解决SPAD仿真常遭遇收敛困难以下是常见原因及对策4.1 物理模型导致的收敛问题问题现象可能原因解决方案雪崩区振荡电离系数设置过高逐步降低an/ap值低电场区发散迁移率模型不连续改用KLAASSEN统一模型温度计算不稳定热模型与电模型耦合过强启用松弛因子(relax0.3)网格点异常值复合模型参数突变检查SRH寿命的温度依赖性4.2 数值方法调整策略method newton trap # 使用Newton-Raphson方法 method itlimit50 # 增加迭代次数 method climit1e12 # 调整载流子浓度限制调试步骤先简化模型确保基础结构收敛逐步添加物理模型监控关键节点的电场和载流子浓度使用log输出详细迭代信息4.3 网格优化实践SPAD仿真的网格设置特别关键mesh width0.1 # 雪崩区网格细化(μm) mesh y.mesh0.05 # 结区垂直方向网格优化原则雪崩区网格尺寸0.1μm渐变掺杂区需要过渡网格衬底区可适当放宽网格边界处避免突变网格密度在最近一次130nm工艺SPAD仿真中通过优化网格将计算时间从8小时缩短至1.5小时同时保持了结果精度。