光电探测器这个方向这几年做的人越来越多尤其是雪崩光电二极管APD和PIN光电二极管这两类结构。PIN结构简单、响应快、偏压低适合做高速接收端APD则靠雪崩倍增效应把微弱光信号放大几十倍甚至上百倍适合长距离、低光功率场景。但问题在于APD的增益和噪声是一对矛盾体——增益越高过剩噪声因子越大带宽反而可能下降。怎么在仿真阶段就把这个平衡点找出来是每个做光电探测器设计的人都绕不开的坎。我这次用Silvaco TCAD做了一套从PIN到APD的完整仿真流程覆盖结构定义、掺杂剖面、电场分布、I-V特性、光谱响应、增益带宽分析这几个核心环节。整个流程跑下来踩了不少坑也积累了一些 Atlas 和 Athena 模块配合使用的经验。下面把完整的设计思路、参数设置、仿真技巧和优化方法拆开来讲适合已经接触过半导体器件仿真、想深入做光电探测器设计的工程师和研究生参考。1. 为什么PIN和APD要放在同一套仿真框架里做1.1 从PIN到APD的演进逻辑很多人做仿真习惯把PIN和APD当成两个独立项目PIN跑完I-V和响应度就收工APD单独建一个结构重新调参。这种做法在初期看起来省事但到了后期优化阶段会非常痛苦——因为你无法判断某个性能变化到底来自材料参数调整还是结构差异。PIN和APD的本质区别在于PIN是本征层吸收光子产生电子-空穴对靠外加电场分离载流子形成光电流APD则是在PIN结构基础上引入一个高电场倍增区让光生载流子在倍增区发生碰撞电离产生雪崩倍增。换句话说APD就是PIN加了一个“信号放大器”但这个放大器自带噪声。把两者放在同一套仿真框架里的好处是你可以先用PIN结构验证材料参数和光学模型是否准确确认吸收系数、载流子迁移率、复合速率这些基础参数没问题之后再在相同材料体系上叠加倍增区结构。这样每次性能变化都能追溯到具体的结构或参数改动而不是一堆变量混在一起。1.2 Silvaco TCAD中Athena与Atlas的分工Silvaco TCAD套件里Athena负责工艺仿真Atlas负责器件仿真。做光电探测器仿真时这两者的配合方式直接决定了仿真效率和精度。我的做法是用Athena跑一遍完整的工艺步骤包括外延生长、离子注入、退火、氧化生成带有实际掺杂剖面的器件结构文件。然后把Athena的输出直接导入Atlas做电学和光学仿真。这样做的好处是掺杂剖面不是理想化的高斯分布或均匀分布而是带有扩散尾部和峰值偏移的真实剖面仿真结果更接近流片实际。但这里有个坑Athena生成的网格在光学仿真时往往不够密特别是在吸收区和倍增区界面附近。如果直接用Atlas默认的网格做光电仿真光生载流子产生率会出现明显的数值振荡。我的处理方式是在Atlas中重新定义网格在吸收区、倍增区、电荷层这三个关键区域手动加密网格间距控制在1nm以内其他区域可以放宽到10nm。1.3 仿真流程的整体架构整个仿真流程我分成四个阶段第一阶段是PIN基础结构验证。用Athena生成一个标准的InP/InGaAs PIN结构吸收层厚度1.5μm本征层掺杂控制在1e15 cm^-3以下。在Atlas中跑暗电流I-V、光电流I-V、光谱响应三条曲线确认吸收系数模型和SRH复合模型参数合理。第二阶段是APD结构叠加。在PIN结构基础上在吸收层和P型层之间插入一个电荷层charge layer和一个倍增层multiplication layer。电荷层的作用是调节电场分布让吸收区电场足够低以避免隧穿漏电同时倍增区电场足够高以触发碰撞电离。第三阶段是参数扫描与优化。对电荷层掺杂浓度、倍增层厚度、吸收层厚度这三个关键参数做扫描观察增益、过剩噪声因子、带宽的变化趋势。第四阶段是工艺容差分析。在实际流片中外延层厚度和掺杂浓度都会有波动仿真阶段需要评估这些波动对器件性能的影响程度。2. 用Athena搭建InP/InGaAs PIN结构的实操细节2.1 外延结构定义与材料参数选择InP/InGaAs PIN结构是1310nm和1550nm波段最常用的探测器结构。InGaAs的吸收边在1.65μm左右对1550nm光有很好的吸收效率InP作为衬底和窗口层对1550nm光透明不会产生额外吸收。在Athena中定义外延结构时我用的典型参数如下层名称材料厚度掺杂类型掺杂浓度衬底InP350μmN型2e18 cm^-3缓冲层InP0.5μmN型1e18 cm^-3吸收层InGaAs1.5μm本征1e15 cm^-3窗口层InP0.5μmP型2e18 cm^-3接触层InGaAs0.2μmP型1e19 cm^-3这里有个细节InGaAs和InP之间的晶格匹配非常关键。In0.53Ga0.47As与InP晶格匹配如果In组分偏离这个值外延层会产生应力导致暗电流增大。在Athena中可以通过设置材料组分来模拟这个效应但更稳妥的做法是在工艺文件中明确标注组分控制精度要求。2.2 离子注入与退火步骤的仿真设置PIN结构的P型窗口层通常用锌Zn扩散或铍Be离子注入形成。我选择Be离子注入因为注入深度容易控制重复性好。Athena中的注入步骤设置implant boron dose5e14 energy40 tilt7 anneal temp650 time30这里用的是硼作为P型掺杂剂实际InP体系中常用Be或Zn但Athena的材料库中硼的参数最完整做基础验证够用。退火温度650度、时间30秒目的是激活注入离子并修复晶格损伤同时避免过度扩散导致PN结位置偏移。退火后的掺杂剖面需要用TonyPlot查看确认结深在0.3-0.5μm之间。如果结深太浅欧姆接触会不好做如果太深窗口层对光的吸收会增加量子效率下降。2.3 网格划分策略与计算效率的平衡网格划分是Athena和Atlas仿真中最容易被忽视但影响最大的环节。网格太粗仿真结果不收敛或者出现假振荡网格太密计算时间成倍增加。我的经验是分区域设置网格吸收区垂直方向网格间距0.5-1nm水平方向5-10nm倍增区APD垂直方向0.2-0.5nm水平方向5nm电极接触区垂直方向1-2nm水平方向10-20nm衬底和缓冲层垂直方向50-100nm水平方向100nm在Atlas中可以用mesh语句手动定义mesh space.mult1.0 x.mesh loc0.0 spacing0.01 x.mesh loc5.0 spacing0.05 y.mesh loc0.0 spacing0.05 y.mesh loc1.0 spacing0.001 y.mesh loc2.5 spacing0.0005 y.mesh loc3.0 spacing0.05注意网格间距变化不要超过10倍否则数值微分误差会显著增大导致I-V曲线出现台阶或毛刺。3. Atlas中光电仿真模型的选择与参数校准3.1 光学模型射线追踪还是光束传播Atlas提供两种光学仿真方式射线追踪ray tracing和光束传播法beam propagation method。对于PIN和APD这类垂直入射结构射线追踪足够用计算速度快而且可以直接得到光生载流子产生率分布。设置射线追踪时需要指定光源参数beam num1 x.origin2.5 y.origin-0.1 angle90 wavelength1.55这里wavelength1.55表示1550nm入射光angle90表示垂直入射。如果要做斜入射或者边缘耦合仿真需要调整angle和x.origin。光功率密度用power参数设置单位是W/cm^2。做响应度仿真时我一般从1e-4 W/cm^2开始这个功率下光生载流子浓度远低于背景掺杂不会引起电场屏蔽效应。3.2 载流子输运模型漂移扩散与碰撞电离PIN结构用漂移扩散模型就够了但APD必须开启碰撞电离模型。Atlas中常用的碰撞电离模型有Selberherr模型和Okuto-Crowell模型。Selberherr模型的电离系数表达式alpha_n A_n * exp(-(B_n/E)^beta_n) alpha_p A_p * exp(-(B_p/E)^beta_p)其中A_n、B_n、beta_n是材料相关的参数。对于InP材料我用的参数是参数电子空穴A (cm^-1)5.55e61.97e6B (V/cm)3.10e62.33e6beta1.01.0这些参数直接决定了APD的增益和噪声特性。如果仿真得到的增益曲线和实验数据偏差较大首先检查的就是这组参数。3.3 复合模型与暗电流来源分析暗电流是光电探测器的核心指标之一。PIN结构的暗电流主要来自三个部分SRH复合电流、表面漏电流、隧穿电流。Atlas中开启SRH复合模型models srh auger bgn fldmob其中srh是SRH复合auger是俄歇复合bgn是带隙窄化效应fldmob是场相关迁移率模型。在InGaAs材料中SRH复合寿命典型值在1-10ns之间。如果仿真暗电流比实验值低一个数量级大概率是SRH寿命设得太长。我一般从5ns开始试然后根据暗电流水平调整。提示InGaAs材料的暗电流对温度非常敏感仿真时如果要做变温分析需要同时调整SRH寿命和本征载流子浓度。4. APD结构设计电荷层与倍增区的参数博弈4.1 电荷层掺杂浓度的窗口效应APD设计中最关键的一步是电荷层charge layer的设计。电荷层的作用是“整形”电场让吸收区电场保持在较低水平避免隧穿漏电同时让倍增区电场足够高触发碰撞电离。电荷层掺杂浓度有一个很窄的窗口。掺杂太低电场会穿通到吸收区导致暗电流增大掺杂太高倍增区电场不足增益上不去。我做过一组扫描电荷层厚度固定50nm掺杂浓度从1e17到5e17 cm^-3变化电荷层掺杂 (cm^-3)吸收区电场 (V/cm)倍增区电场 (V/cm)增益0.95Vbr1e172.5e54.2e582e171.8e55.1e5253e171.2e55.8e5454e170.8e56.3e5605e170.5e56.8e535可以看到掺杂浓度在3e17到4e17之间时增益达到较高水平同时吸收区电场控制在1e5 V/cm以下。超过4e17之后倍增区电场虽然继续升高但增益反而下降这是因为电场过高导致空间电荷效应增强有效倍增长度缩短。4.2 倍增层厚度对增益带宽积的影响倍增层厚度直接决定了APD的增益带宽积GBP。薄倍增层载流子渡越时间短带宽高但增益低厚倍增层增益高但渡越时间长带宽下降。我对比了三种倍增层厚度100nm、200nm、300nm。在相同电荷层掺杂浓度下100nm增益30带宽25GHzGBP750GHz200nm增益50带宽12GHzGBP600GHz300nm增益70带宽6GHzGBP420GHz从GBP角度看100nm倍增层最优。但实际设计中还要考虑增益的绝对值——如果前置放大器噪声较大需要APD本身提供足够增益这时候200nm可能是更好的折中。4.3 吸收层厚度与响应速度的取舍吸收层厚度决定了量子效率。厚度越大光吸收越充分响应度越高。但光生载流子需要渡越整个吸收层才能到达倍增区厚度越大渡越时间越长带宽越低。对于1550nm波长InGaAs的吸收系数约为0.7μm^-1。吸收层厚度1.5μm时量子效率约65%厚度2.0μm时量子效率约75%厚度3.0μm时量子效率约88%。但带宽的变化更剧烈1.5μm吸收层的渡越时间约8ps2.0μm约11ps3.0μm约17ps。如果目标带宽是10GHz以上吸收层厚度不宜超过2.0μm。我的建议是如果系统对灵敏度要求高、带宽要求相对宽松如2.5GHz可以用2.5-3.0μm吸收层如果要求10GHz以上带宽吸收层控制在1.5μm以内靠后续跨阻放大器补增益。5. 仿真收敛问题与常见报错的排查路径5.1 碰撞电离模型导致的收敛困难APD仿真中最常见的报错就是“convergence failure”尤其是在接近击穿电压时。原因是碰撞电离率对电场呈指数依赖电场微小变化会导致电离率剧烈变化牛顿迭代法很容易发散。我的处理方式分三步第一步先用低偏压跑通确认结构和掺杂没问题。从0V开始步长0.1V逐步增加到5V。第二步在接近击穿电压时减小步长。Atlas中可以用method语句设置method newton trap maxtrap10 solve vanode0.0 vstep0.1 vfinal5.0 nameanode solve vanode5.0 vstep0.05 vfinal10.0 nameanode solve vanode10.0 vstep0.01 vfinal12.0 nameanode第三步如果还是发散开启trap和maxtrap参数让求解器在发散时自动减小步长。注意击穿电压附近的仿真结果对网格密度极其敏感。如果发现击穿电压随网格加密而显著变化说明网格还不够密需要继续加密倍增区网格。5.2 光电流仿真中的数值振荡做光谱响应仿真时如果光生载流子产生率在空间上分布不均匀容易出现数值振荡。表现为响应度曲线在某个波长附近出现锯齿状波动。根本原因是光吸收的指数衰减和网格离散化之间的不匹配。解决办法有两个一是加密吸收区网格让每个网格内的光功率衰减不超过5%二是开启Atlas的optical模块中的平滑选项。我一般优先加密网格因为平滑选项虽然能消除振荡但会引入额外的数值误差。5.3 温度相关参数设置错误做变温仿真时如果只改了temperature参数而没有同步调整材料参数结果会完全不可信。Atlas中温度相关的参数包括带隙、本征载流子浓度、迁移率、SRH寿命、碰撞电离系数。我的做法是在Atlas输入文件中显式定义温度系数material materialInGaAs eg3000.75 egalpha4.5e-4 egbeta200 material materialInGaAs taun05e-9 taup05e-9然后通过solve temperature300到solve temperature350逐步扫描。6. 从仿真曲线到设计结论关键指标提取方法6.1 增益曲线的正确提取方式APD增益定义为光电流与初级光电流之比。在Atlas中初级光电流是未开启碰撞电离模型时的光电流倍增光电流是开启碰撞电离模型后的光电流。具体操作先跑一遍不开启碰撞电离的仿真记录光电流I_primary再跑一遍开启碰撞电离的仿真记录光电流I_total。增益M I_total / I_primary。这里有个容易犯的错误直接用暗电流做分母。暗电流不包含光生载流子的倍增信息用暗电流算出来的“增益”没有物理意义。6.2 过剩噪声因子的仿真估算过剩噪声因子F(M)是APD的核心噪声指标。Atlas本身不直接输出F(M)但可以通过碰撞电离系数比kalpha_p/alpha_n来估算F(M) k*M (1-k)*(2 - 1/M)对于InP材料k值约为0.4-0.5。这意味着即使M10F(M)也在4-5左右噪声相当可观。如果仿真中需要更精确的F(M)可以用Monte Carlo方法但计算量极大。工程上一般用上面的经验公式估算就够了。6.3 带宽的交流小信号仿真设置带宽仿真需要在Atlas中开启交流小信号分析solve ac freq1e6 fstep10 nfreq100频率从1MHz扫到10GHz观察光电流的幅度随频率的变化。带宽定义为光电流幅度下降到低频值0.707倍-3dB对应的频率。做交流仿真时要注意光源的调制频率要和ac语句中的频率扫描范围匹配。如果光源是直流光交流仿真得到的是暗电流的频响不是光电流的频响。正确做法是在beam语句中设置modulate参数让光源强度按正弦规律变化。7. 工艺波动对器件性能的影响评估7.1 外延层厚度波动的敏感度分析实际外延生长中层厚控制精度一般在±5%左右。对于1.5μm吸收层厚度波动范围是±75nm。这个波动对响应度和带宽的影响需要量化。我做了三组仿真吸收层厚度1.425μm、1.5μm、1.575μm。结果显示响应度变化约±3%带宽变化约±5%。这个敏感度在可接受范围内。但倍增层厚度波动的影响更大。200nm倍增层波动±10nm时增益变化约±15%带宽变化约±12%。所以倍增层的生长精度要求比吸收层高一个等级。7.2 掺杂浓度波动的蒙特卡洛分析电荷层掺杂浓度波动对增益的影响最显著。我用蒙特卡洛方法跑了100组样本掺杂浓度按高斯分布设置均值3.5e17 cm^-3标准差5%。结果显示增益的均值为48标准差为6.2变异系数约13%。这意味着如果掺杂控制精度只有5%增益的一致性大约在±13%左右。对于要求高一致性的阵列探测器这个波动偏大需要把掺杂控制精度提高到2%以内。7.3 温度漂移对击穿电压的影响APD的击穿电压对温度非常敏感温度系数约为0.1%/K。在25度到85度的工业级温度范围内击穿电压漂移约6%。这意味着如果偏压固定增益会随温度显著变化。仿真中可以通过变温I-V曲线提取击穿电压的温度系数然后在电路设计中加入温度补偿。常见的补偿方式是使用温度传感器加可调偏压源或者用相同工艺的参考二极管做跟踪补偿。8. 一些仿真之外的实操体会仿真跑通只是第一步真正流片时还有很多仿真无法覆盖的因素。比如欧姆接触的实际比接触电阻往往比仿真值高一个数量级这是因为仿真中用的理想接触边界条件在现实中很难实现。我的做法是在仿真中把接触电阻设为一个可调参数然后根据流片后的测试数据反推实际值再更新仿真模型。另外InGaAs材料的表面复合速度对暗电流影响很大。仿真中默认的表面复合速度是1e3 cm/s但实际工艺中如果钝化不好表面复合速度可能达到1e5 cm/s甚至更高。做暗电流仿真时我一般会把表面复合速度设为一个扫描变量观察暗电流的变化范围这样在流片后就能快速判断暗电流偏高的原因是不是表面钝化问题。还有一点Atlas的光学仿真默认假设材料吸收系数是常数但实际上InGaAs的吸收系数在吸收边附近随波长变化很剧烈。如果做宽谱响应仿真需要手动导入实测的吸收系数数据而不是用材料库中的默认值。最后说一个关于仿真收敛的小技巧。如果APD仿真在击穿电压附近反复发散可以尝试先跑一个低精度的“预仿真”用粗网格和宽松的收敛条件得到一个近似解然后把这个解作为初值再用细网格和高精度条件重新求解。Atlas支持通过save和load语句保存和加载中间状态这个技巧能节省大量调试时间。