光电成像器件CAD仿真：从物理模型到像素阵列设计

简介以电场三大定理、连续超张弛迭代算法和拉格朗日插值算法为基础结合C语言编程的光电成像器件CAD实验报告完整实现了像管内电势计算、等势线绘制与对称轴电位求解的程序设计与验证。报告先推导有限差分方程再引入超张弛迭代加速收敛并处理非封闭电极边界内容覆盖数值求解、程序实现与结果分析全过程。报告结构包含摘要、理论基础、算法推导与程序实现说明重点展示了拉格朗日插值在对称轴电位计算中的应用。面向光电工程、微光夜视与电子光学方向的高年级本科生和初学者可作为课程设计或相关课题的参考模板。资源为1个PDF文件大小743KB已有170人学习。通过梳理其算法流程、边界处理技巧和程序框架可加深对电磁场数值计算与CAD方法结合应用的理解也为后续像管性能分析、电子轨迹模拟等实践提供可直接借鉴的编程思路。1. 光电成像器件CAD的真正入口把“画图”换成“求解”看到“光电成像器件计算机辅助设计CAD实验报告1”这个题目先想清楚一件事这个实验课要求的不是用AutoCAD画剖面图。光电成像器件从CMOS图像传感器的像元到CCD、InGaAs短波红外探测器的CAD是指Computer Aided Design里的辅助设计把坐标、掺杂和光照写成偏微分方程在计算机里解出电势分布、载流子浓度和电流密度。一个PN结画得再好看不跑出暗电流和量子效率就不算设计完成。所以这套流程最重的部分在器件结构建模、材料参数选取和求解器设置最后把J-V曲线、QE光谱和实测去对表实验报告才算有信息量。2. 光电成像器件仿真的数学框架与材料参数库2.1 从泊松方程到漂移-扩散光电成像器件CAD为什么是耦合求解凡是用仿真器做半导体器件的CAD核心都是三个偏微分方程泊松方程给出电势与净电荷的关系电子和空穴两条连续性方程给出载流子浓度随时间和位置的演化。电场方向漂移电流由Jn qμnnE qDn∇n描述扩散项在浓度梯度大的结区尤其重要。对于光电成像器件光照会被写进连续性方程的生成项G里而不是当成边界条件随便挂一个电流源。这意味着电势、电子浓度、空穴浓度三个未知量必须一起解。反向偏置下的耗尽区里电场很大载流子浓度低如果只解泊松方程看静态能带图那看到的是热平衡状态完全推不出光电流。光生载流子产生之后能不能被电极收集、收集效率多高取决于电场分布和寿命这两者反过来又受载流子浓度影响所以仿真器内部用Newton迭代做耦合求解。成像器件的工作条件通常是反向偏置因此产生机制和复合机制的选择比正向二极管仿真更讲究。暗电流主要来自耗尽区SRH复合产生表面复合速度对短波响应影响大高掺杂区要打开禁带变窄模型。仿真参数的初值可以由小面积二极管测试结构校准再复用到整个像元阵列这是我做彩色CMOS像素时最常用的标定路径。2.1.1 成像器件与普通二极管的仿真差异普通功率二极管仿真往往只关心导通压降和击穿成像器件仿真要关心的事更细暗电流密度要按nA/cm²量级评估量子效率要分波长看像素间串扰要算到百分之几。这要求在仿真器里把吸收系数、载流子寿命、表面复合速度设成波长、温度、掺杂浓度的函数而不是用一套默认值算到底。2.2 材料参数库带隙、吸收系数与InGaAs的低暗电流代价光电成像器件的材料选择基本决定仿真参数的量级。下面这张表把几种常见材料在300 K下的关键参数列到一起做实验报告时可以直接用来核对自己输入的材料文件材料带隙(eV)截止波长(μm)0.7μm吸收系数(cm⁻¹)典型应用Si1.121.10约2×10³CMOS图像传感器、CCDGe0.661.88约1×10⁴近红外探测、SiGe光电二极管In₀.₅₃Ga₀.₄₇As0.741.68约5×10³短波红外焦平面GaAs1.420.87约5×10³可见光/X射线直接探测带隙决定暗电流下限因为本征载流子浓度ni随带隙指数变化。同样面积下InGaAs探测器工作在室温时的暗电流比硅器件大好几个数量级这不是工艺差而是材料属性决定的。所以做CAD仿真时如果InGaAs探测器报了很高的暗电流先看温度模型和ni公式是否启用了窄禁带修正不要急着调寿命参数。能带随温度变化用Varshni公式拟合比较常见参数alpha和beta在材料的TCAD模型库里都有def varshni_bandgap(Eg0, alpha, beta, T): # Eg0是0K时的带隙alpha/beta是材料拟合参数 # 对硅: alpha4.73e-4 eV/K, beta636 K return Eg0 - (alpha * T * T) / (T beta)这段脚本的逻辑很简单但不代表可以跳过很多COMSOL和Silvaco案例文件默认把带隙设成300K的常数如果你只仿真室温暗电流误差不大一旦做温度区间扫描ni会算错数倍暗电流随温度的斜率就不对。实验报告里如果画了ln(J/T²)对1/kT的Arrhenius曲线这个温度模型就是前提。2.3 光照输入体产生率、面产生率与Beam设置的等价关系光进入半导体后按Beer-Lambert吸收定律被逐渐衰减。仿真器内部通常把光强换算成光生载流子产生率G(x)公式里最关键的是吸收系数α(λ)不同波长对应的α差了三个量级。红光的吸收深度在硅里能达到几十微米蓝光靠近表面0.1μm内就被吸收完所以表面复合速度对蓝光响应的影响要大得多。最常写进CAD脚本的参数是光子流密度phi0而不是光功率密度P。两者换算要看单光子能量import math def generation_rate(x_um, phi0, alpha, R0.1): # x_um: 距离表面的深度单位μm # phi0: 入射光子流密度单位 1/(cm2*s) # alpha: 该波长的吸收系数单位1/cm x_cm x_um * 1e-4 return phi0 * (1 - R) * alpha * math.exp(-alpha * x_cm)代码里R是表面反射率。光强不是直接乘进去的先换算成phi0再乘吸收系数否则光电流数值会差一个量级。对应到Silvaco、TCAD、Lumerical的Beam或Port设置里需要指定波长、光功率密度和入射方向求解器按介质折射率和吸收系数自动算G(x)。表层的反射率通常当成常数但在减反膜设计里要按薄膜干涉计算这就是另外一块工作。两种简化写法在实验里经常混用。面产生率适合超薄吸收层光在吸收层里衰减不明显直接当成表面边界通量体产生率适合厚衬底蓝紫光在表面就产生了这个细节决定了光谱响应曲线从短波端哪个位置开始下降。把两种光线条件各跑一次如果QE曲线差异在短波区很明显说明表面复合参数没有校准。2.4 网格策略耗尽区加密和界面陷阱对收敛的影响器件仿真器的网格密度直接影响结果和收敛行为这里有一个容易忽略的点光电成像器件的反偏耗尽区是电场峰值所在网格不够密碰撞电离和SRH复合的积分值会被低估计算出的暗电流偏低。像素顶层的光电二极管和转移栅交叠区域沟道电势突变网格尺寸建议从表面到耗尽区逐步过渡比设置为从0.5μm到1μm再到2μm这样跳。区域推荐网格尺寸原因表面0~0.1μm≤5 nm短波长光生载流子在此处集中耗尽区≤20 nm电场变化快SRH产生项积分精度要求高中性衬底100~500 nm影响串联电阻可适当放松网格突变比超过1.5时Newton迭代容易发散。遇到不收敛第一步不是改物理模型而是看网格过渡区相邻网格尺寸比例。成像器件的界面态陷阱与表面复合速度耦合如果这个区域网格太粗表面复合产生的暗电流可能被低估几个量级。提示Silvaco和Sentaurus默认的网格生成器对垂直方向的加密比较保守我一般会在掺杂浓度梯度和电场梯度两个判据都打开的情况下再手动加一处表面加密。单纯依赖自适应网格容易出现“看起来收敛、结果却不对”的情况。3. 光电成像器件仿真的J-V与量子效率提取流程3.1 从光照到J-V曲线量子效率要先于灵敏度换算光电成像器件的仿真结果通常分成两组输出暗态下扫描反向偏压得到暗J-V曲线光照下得到光J-V曲线。两线之差是光电流密度Jph量子效率的定义是收集到的电子数除以入射光子数落到电流上就是QE Jph / (q × phi)。很多同学直接把光J-V当成QE的曲线画但那只是电流换成A/W灵敏度还要再乘光子能量。灵敏度R(A/W)和QE的换算经常在实验报告的“数据分析”部分出错R QE × λ / 1.24其中λ用μm单位代入这个系数来自hc/q。0.55μm波长时理想的100% QE对应0.44 A/W看到仿真结果超过0.5 A/W先别兴奋大概率是光子流密度算错了或者把光功率密度当成光子流密度用了。3.2 最小可复现的器件仿真输入骨架不同CAD工具的命令语法不同结构骨架是共同的。以下用一个可读的配置范本说明每段参数的作用具体关键词以你手上工具为准# 器件仿真输入骨架结构上对应Silvaco/TCAD/COMSOL的公共设置 set materialSi temp300 srh_lifetime1e-6 set mesh depth5 refinementsurface:3nm, depletion:15nm set doping profileanalytic peak1e17 junction_depth0.3um set electrode anodecontact(P) cathodecontact(N) # 光照与偏置 set light sourcetop spectrummono wavelength550nm set power_density1e-3 # W/cm2室内成像环境的典型光强 solve bias-1.0V to0.6V step0.02V extract dark_current at bias-0.1V extract photocurrent at bias-0.1V这段配置要理解三件事。第一srh_lifetime是暗电流的主控参数默认从1μs开始调。第二junction_depth决定耗尽区的电场分布它由工艺仿真的注入能量和退火时间决定不能随意填。第三光功率密度1mW/cm²是相对弱的照明如果做高动态范围成像要按实际曝光量换算成等效连续光功率再仿真不能拿强光参数套用到弱光暗电流。3.3 从光谱扫描到QE曲线积分时别混用波长单位光谱响应仿真要逐波长扫描跑完350nm到1000nm的几十个点最后拟合出一条QE曲线。这个环节最容易在单位换算上翻车。下面是数据处理段的常见写法import numpy as np # wl: 波长 [nm], i_ph: 光电流密度 [A/cm2] wl np.array([400, 500, 600, 700, 800, 900, 1000]) i_ph np.array([0.10, 0.15, 0.17, 0.16, 0.14, 0.10, 0.05]) * 1e-3 def qe_from_photocurrent(wl_nm, i_ph): q 1.602e-19 h 6.626e-34 c 3.0e8 lam wl_nm * 1e-9 phi (1e-3 * lam) / (h * c) # 光子流密度 return i_ph / (q * phi) qe qe_from_photocurrent(wl, i_ph)函数里光功率密度写死为1mW/cm²与仿真输入保持一致。光子流密度用P×λ/(hc)计算因为单光子能量是hc/λ波长越短光子越多。这里有一个反向验证技巧把QE对波长积分估算出的总光电流如果和实验里的sensor响应差别非常大那多半是反射率R没有参与计算。3.4 暗电流仿真的三个翻车点暗电流密度通常在nA/cm²量级仿真器很容易在数值噪声里把暗电流算成负的或者振荡。第一个坑是SRH寿命参数统一用一个常数。器件最上层的表面区域和体内衬底寿命相差很大表面层通常要单独指定较短寿命。第二个坑是忘记关闭光照再跑暗态扫描导致“暗电流”里带了几十个nA的光电流。第三是扫描方向反偏暗电流要在热平衡解的基础上从零偏压往负方向扫描直接从正偏起始点扫过去瞬态数值抖动会污染最后几组暗电流数据。提示在Silvaco/TCAD里跑暗电流时把复合模型的“Auger”和“BGN”都打开但寿命参数先用经验值等校准之后再精调。暗电流本来就是对复合模型极其敏感的指标三组参数同时调的结果几乎无法解释。4. 光电成像器件工艺仿真到像素阵列串扰与开孔率分析4.1 工艺仿真到器件仿真的数据传递先核对结深光电成像器件的实际结构不是解析掺杂公式生成的而是扩散、注入、氧化、退火一步步推出来的。常见的流程是先在工艺仿真器里把离子注入能量、剂量、退火时间跑完得到掺杂分布和应力场再导出结构文件到器件仿真器做电学求解。常见的做法是Silvaco保存.str结构Sentaurus导出.tdrCOMSOL用导入网格与材料域的方式对接。数据传递里最容易出错的是单位。工艺文件里注入能量的单位是keV结深是nm而器件仿真器把坐标当cm处理。我一般会在导入后先数一下坐标范围再输出x从表面到衬底的掺杂浓度截面图核对结深是否与工艺处方一致。这个步骤如果跳过后面J-V曲线上的击穿电压偏大或偏小都找不到原因。4.2 用Python批量修改CAD几何对同一像元跑开孔率扫描成像器件设计里要反复回答一个问题光电二极管开口做多大。开口小了量子效率下降开大了相邻像素的串扰上升。逐个在GUI里画太慢直接写Python脚本批量生成几何再喂给光学或电学仿真器。GDS格式是版图与器件仿真之间的常用中间格式下面的脚本可以生成一组不同开口宽度的像元版图import gdspy lib gdspy.GdsLibrary() pitch 5.0 # 像元周期 5um for ratio in (0.4, 0.5, 0.6, 0.7): cell lib.new_cell(fPIXEL_R{int(ratio*100)}) open_w pitch * ratio # layer 1: 光敏区 gdspy.rectangle((0, 0), (open_w, open_w), layer1, datatype0) # layer 2: 金属挡光层覆盖光敏区之外的部分 gdspy.rectangle((0, 0), (pitch, pitch), layer2, datatype0) lib.write_gds(pixel_opening_scan.gds)这段代码的关键参数是pitch和ratio。ratio从0.4到0.7代表光敏区开口占整个像元的比例也就是开孔率。GDS导出后可以直接导入FDTD或蒙特卡洛光线追迹工具用来统计不同入射角下的能量分布。版图上的layer编号要和工艺厂掩膜层次对应否则导入后layer 2的金属挡光层会被当成电介质仿真结果完全失真。4.2.1 开孔率与微透镜的配合有微透镜的情况下光会被聚焦到光敏区中心开孔率的影响比裸片小。但微透镜的聚焦能力有限大入射角时焦点偏离光敏区串扰依然上升。仿真时至少扫描三种入射角0°、15°和30°不要只做垂直入射。4.3 像素串扰与开孔率的折中几何近似先行全电磁波光学的串扰仿真耗时很长初期设计可以直接用几何光束近似判断边界。串扰比约等于tan(θ)×深度/间距其中深度是光电二极管到感光面的距离间距是相邻像素中心的间距。用这个公式可以算出一个粗略的上下限超过目标值再上高精度仿真。开孔率集光率15°入射角串扰趋势适用场景40%低低高分辨率小像元60%中中安防照度中等80%高高低照度成像表中的趋势在30°入射角时会更明显。如果结构已经做了深沟槽隔离串扰主要由光学衍射贡献几何近似就不准了此时需要用FDTD计算点扩展函数。实验报告如果能放一张“不同开孔率下的串扰比例”柱状图比单列一组J-V曲线更有区分度。5. 光电成像器件仿真结果如何验收Arrhenius验证与行为模型5.1 温度扫描用激活能判断暗电流来源全结构仿真做出来的暗电流不能只看室温单点至少跑四个温度点230 K、260 K、290 K、320 K。把暗电流按Arrhenius形式整理成ln(J/T²)对1/kT的曲线斜率就是激活能Ea。import numpy as np from scipy.optimize import curve_fit # 仿真得到的不同温度暗电流密度单位 A/cm2 T_K np.array([230, 260, 290, 320]) J_dark np.array([2.1e-13, 5.4e-12, 8.9e-11, 1.2e-9]) kB 8.617e-5 def arrhenius(invT, prefactor, Ea): return prefactor - Ea / kB * invT ln_JT2 np.log(J_dark / T_K**2) fit, _ curve_fit(arrhenius, 1/T_K, ln_JT2) print(Ea , fit[1], eV)激活能接近材料禁带宽度Eg说明暗电流以中性区扩散为主接近Eg/2说明耗尽区SRH产生主导。成像器件工作在反偏看到Eg/2不必意外。如果拟合出的激活能小于Eg/2一截基本可以断定仿真网格在耗尽区太粗或者表面复合速度被设成了不合理的超大值。5.2 从单点拟合到全曲线校准参数时的调整顺序校准仿真和实测的差值时不要同时动寿命、表面复合速度、禁带变窄三个参数。我习惯的顺序是先用高温暗电流校准体寿命因为高温下SRH产生占主导再用低温暗电流校准表面复合速度低温下表面漏电相对突出最后才调节衬底掺杂浓度或禁带变窄修正。光电流的QE校准反过来先调反射率再调短波区的表面复合最后动长波区的吸收系数和衬底寿命。这个顺序可以保证每次只有一个变量在解释偏差。5.3 把J-V表落成行为模型仿真结果怎么往外“交货”器件仿真结果最终要交给像素电路设计用常见做法是把J-V表导出成查表模型而不是复杂TCAD网格。将仿真所得的电压、暗电流、光电流三列数据写入CSV在网表里以分段线性电流源的方式调用这样一次TCAD扫描可以支撑整列像素的时序仿真。给行为模型做验收时用三档光强回灌宏模型查看满阱电荷量和读出噪声是否在预期范围这一步基本能省掉一轮试错流片。把宏模型的表头列命名与SPICE子电路端口对应好后续的电路级仿真可以直接复用这份输出文件不用再回头改器件仿真脚本。本文还有配套的精品资源点击获取