从零开始掌握Silvaco ATLAS SPAD仿真物理模型配置实战指南1. 环境准备与基础认知在开始SPAD器件仿真之前我们需要先理解几个关键概念。SPADSingle Photon Avalanche Diode作为单光子探测的核心元件其仿真精度直接依赖于物理模型的准确配置。不同于常规二极管仿真SPAD需要特别关注雪崩击穿区的载流子行为模拟。必备软件环境Silvaco ATLAS 5.28.2.R或更新版本DeckBuild图形界面推荐新手使用TonyPlot可视化工具提示首次启动ATLAS时建议创建独立项目目录所有仿真文件保持路径无中文和空格安装完成后建议优先运行内置的example/ex18_apd.in示例文件观察标准APD的仿真流程。这个示例虽然不针对SPAD但包含了碰撞电离等关键模型的配置方法。关键参数认知误区低场迁移率μn0/μp0并非固定值实际与掺杂浓度强相关SRH复合寿命TAUN0/TAUP0默认值针对硅材料其他材料需手动调整碰撞电离模型选择直接影响击穿电压的仿真精度2. 物理模型配置全流程2.1 迁移率模型设置在DeckBuild界面中迁移率模型通过MODELS和MOBILITY语句组合配置。对于SPAD仿真推荐采用Masetti模型掺杂相关结合高场饱和效应models conmob fldmob mobility materialsilicon masetti参数对照表参数名物理意义硅默认值适用场景MUN电子迁移率1450 cm²/Vs低场区域MUP空穴迁移率500 cm²/Vs低场区域VSATN电子饱和速度1.0e7 cm/s高场区域VSATP空穴饱和速度0.8e7 cm/s高场区域注意当掺杂浓度超过1e18 cm⁻³时必须启用CONMOB参数以避免迁移率高估2.2 复合模型配置SPAD仿真的核心在于准确描述载流子的产生-复合过程。在DeckBuild中输入以下关键配置models srh auger optr material silicon taun01e-7 taup01e-7复合模型选择逻辑SRH复合必须启用描述缺陷辅助复合Auger复合高注入条件下必需光学复合光子探测场景必选典型报错排查出现SRH parameter not defined检查material语句是否正确定义taun0/taup0出现Auger coefficient missing添加augn2.8e-31 augp9.9e-32参数2.3 碰撞电离模型实战雪崩效应是SPAD工作的物理基础ATLAS提供两种主流模型impact selb # 或 impact crowell模型对比实验数据电场强度 (V/cm)Selberherr误差Crowell误差2.0e512%8%3.5e518%5%5.0e525%15%关键发现在3-4e5 V/cm典型工作区间Crowell模型表现更优3. 参数优化技巧3.1 网格划分策略SPAD的雪崩区需要特殊网格加密mesh width0.1 x.mesh loc0.00 spac0.01 x.mesh loc0.05 spac0.002 # 雪崩区加密 y.mesh loc0.00 spac0.05 y.mesh loc0.10 spac0.01 # 结区加密网格质量检查标准雪崩区节点间距 ≤ 0.1 μm纵向网格长宽比 5:1过渡区网格渐变比 ≤ 1.23.2 收敛性调优针对SPAD高场强仿真的收敛技巧采用渐进式偏压加载solve vstep0.1 vfinal30启用载流子温度模型models latt.temp carr.temp调整迭代参数method newton trap maxtraps5常见发散场景处理振荡发散减小vstep或增加maxtraps收敛停滞尝试method gummel与newton交替使用4. 结果验证与案例解析4.1 关键指标提取方法在TonyPlot中分析击穿特性时推荐使用以下TCL命令脚本set Vbd [lindex [extract curve.iv curveAnode Current] 0] puts 击穿电压 $Vbd V标准验证流程检查暗电流曲线是否符合exp(V)关系验证光子响应率在击穿区呈现显著增益确认电场分布峰值位于设计雪崩区4.2 典型SPAD结构仿真示例以p/n阱结构为例的完整输入文件框架go atlas mesh # 网格定义 ... region # 区域定义 ... doping # 掺杂分布 ... models # 物理模型选择 srh auger optr bbt impact selb mobility masetti material # 材料参数 silicon taun01e-7 taup01e-7 solve # 偏压扫描 ... save # 结果输出 ... tonyplot # 可视化设置 ...参数敏感度分析雪崩区宽度变化±20% → 击穿电压偏移8-15%掺杂浓度变化1个数量级 → 暗电流变化3-5倍复合寿命变化50% → 光子探测效率变化10-30%