简介面向半导体光电器件与Micro-LED显示技术领域的科研人员和工程师这份资源是一份高光效GaN基Micro-LED仿真模型研究的完整文档核心针对芯片尺寸微缩带来的侧壁效应详细阐述利用有限差分时域法FDTD构建与优化模型的思路。文档为单个docx文件大小约2.31MB内容包含初始垂直侧壁模型建立、多量子阱位置影响分析、底部反射材料对比以及顶部透射光栅的周期、高度、占空比参数优化等关键步骤最终给出LEE较初始结构提升2.42倍的结论。已有303人学习该资源适合正在从事Micro-LED光提取效率提升、微纳光学仿真或相关课题研究的读者作为模型设计参考。1. 为什么我非要去建一个Micro-LED仿真模型先交代一下背景。我是在做第三代半导体器件研究时接触到GaN基Micro-LED的最早接触它是因为一个很现实的需求Micro-LED的样品流片周期太长、成本太高一次实验下来光代工费就够呛而且参数调试完全没有灵活性。当时团队里有个做工艺的兄弟提了个需求——能不能在流片之前先通过仿真筛一遍外延结构和量子阱参数把大概率能跑出高外量子效率的结构定下来再上产线。这个需求听着很合理但真正落地的时候才发现GaN基Micro-LED的仿真模型远比想象中难搞。先给刚开始接触这个方向的朋友简单说一下Micro-LED是什么。Micro-LED就是芯片尺寸小于100微米甚至小于50微米的发光二极管和传统的大尺寸LED相比它在亮度、响应速度、寿命上都有优势是下一代显示技术的热门候选。但Micro-LED有一个非常典型的物理困境芯片尺寸越小侧壁缺陷占比越高非辐射复合越严重效率急剧下降。这就是所谓的“尺寸效应”。如果你想在芯片做小之后还能保持高光效就必须在结构设计上做文章——比如优化量子阱层数、调整势垒厚度、设计侧壁钝化方案——而这些工作如果全靠在实验台上试错时间和成本都承受不住。仿真模型的价值就在这里它能在不流片的情况下把量子阱区的载流子分布、辐射复合速率、俄歇复合损失这些关键物理量先跑一遍告诉你哪种结构组合更有可能出高光效。但要建一个靠谱的仿真模型远不是打开软件、选几个参数、点击运行那么简单。从物理模型的选取到边界条件的设定从材料参数的标定到量子效率曲线的拟合每一步都有坑。这篇博文就是我把自己在GaN基Micro-LED仿真建模这条路上踩过的坑、验证过的方法、总结下来的经验梳理一遍给同样在搞Micro-LED仿真的朋友一个参考。这篇内容的目标读者很明确——你可能是研究生、工程师正在用Silvaco、Crosslight或者COMSOL这类工具做Micro-LED或LED器件的仿真也可能只是想搞清楚“仿真模型到底怎么搭才能贴近真实器件”。我会把从器件结构设定、物理模型选择到参数校准、结果验证的完整链路讲清楚重点放在那些仿真软件里不会明说、但实际影响结果判断的关键细节上。2. 建模仿真的第一个分水岭你想要的到底是电学模型还是光学模型很多第一次接触Micro-LED仿真的人上来就问“用哪个软件好”。这个问题的前提其实是错的。真正该先问的是你仿真要回答什么问题不同的物理问题对应的模型维度和求解方程完全不一样选错了方向后面全白搭。2.1 两类模型的本质差异Micro-LED仿真模型大体可以分成两大类。第一类是电学模型核心是求解半导体器件物理方程组——泊松方程、载流子连续性方程、漂移扩散方程通过这些方程算出量子阱区的载流子浓度分布、复合速率、电流-电压特性。大部分用Silvaco Atlas或者Crosslight APSYS做的工作都属于这一类。第二类是光学模型核心是用电磁场仿真——比如FDTD时域有限差分或者RCWA严格耦合波分析——去计算Micro-LED的光提取效率、远场分布、光谱特性。这两类模型面对的是完全不同的物理过程对网格划分、材料参数、边界条件的处理方式也截然不同。更麻烦的是Micro-LED的光提取效率问题和尺寸、侧壁角度、电极布局密切相关而这些因素对电学特性的影响往往可以忽略反之亦然。所以你经常会看到有人用两种软件各建一个模型一个算电学一个算光学然后想办法把两个结果耦合起来。这就是所谓的“多物理场协同仿真”。如果你刚起步我建议先聚焦在一个问题上先把电学仿真做扎实因为你最想优化的外量子效率首先受限于内量子效率而内量子效率是由电学仿真算出来的。2.2 仿真维度怎么选1D、2D还是3D这是个特别容易忽略、但影响巨大的决定。Micro-LED和传统大尺寸LED在仿真上最大的不同就是尺寸效应带来的三维特性。如果是做大面积功率LED毫米级别传统的做法是建一个一维或者二维模型沿着外延生长方向c轴方向做切片分析把横向效应忽略掉计算量小结果也基本够用。但Micro-LED不行。当芯片尺寸缩小到几十微米甚至几微米的量级后侧壁区域的非辐射复合对整个器件的贡献变得非常显著你必须把侧壁效应包含进来。这时候至少得用二维模型——取芯片的一个横截面同时考虑垂直方向的载流子输运和水平方向的侧壁效应。如果你做的是像素尺寸在10微米以下的Micro-LED而且涉及电流拥堵、侧壁损伤层或钝化层设计这类问题我建议直接上三维模型。三维模型的计算量会成倍增长尤其当你还需要求解量子阱区的耦合量子力学问题时对硬件和算法都是考验。但没办法Micro-LED的科学问题本身就长在三维空间里简化模型的代价是丢失你最关心的物理效应。3. 核心物理模型的选取原则ABC模型和它的“坑”确定了仿真维度之后下一个关键环节就是选物理模型。这个环节的重要性怎么强调都不过分——模型选错了参数调得再漂亮结果都是自欺欺人。3.1 ABC复合模型到底是什么在III-V族氮化物LED仿真里最核心的复合机制描述就是ABC模型。名字听起来吓人其实逻辑很简单量子阱里的载流子复合总速率由三部分构成——SRH非辐射复合A项、双分子辐射复合B项、俄歇复合C项复合速率R可以写成一个简单的三项式R A·n B·n² C·n³其中n是量子阱区载流子浓度。A项是缺陷辅助的非辐射复合只要载流子被缺陷能级捕获就会复合掉不发光和载流子浓度成线性关系B项是导带电子和价带空穴的辐射复合发光的主力和浓度的平方成正比C项是俄歇复合——两个电子碰撞把能量交给第三个载流子不发光和浓度的立方成正比。仿真软件内部其实是联立求解漂移扩散方程来间接得到这些复合速率的不是直接套这个多项式。但ABC模型的思想贯穿在仿真结果分析里你最终拟合出来的IQE内量子效率曲线、效率droop现象的解释都离不开这三项的博弈。3.2 效率droop的物理根源仿真最容易扯皮的地方GaN基LED效率droop是经典难题——电流密度增大后外量子效率先升后降。关于droop的根源学术圈争论了很久有人说是俄歇复合主导有人说是载流子泄漏主导还有人提出是极化场导致的载流子局域化。做仿真时你必须选一个立场因为你的模型参数会直接影响droop的仿真结果。我个人的经验和建议是在做Micro-LED仿真时应该默认同时开启俄歇复合和载流子泄漏通道然后通过实验数据去反推哪个因素占主导。因为Micro-LED的工作电流密度通常比传统LED大为了补偿尺寸变小后的亮度损失Micro-LED经常在高电流密度下工作在几十到几百A/cm²的电流密度范围内俄歇复合的C项三次方项会迅速增长对droop的贡献非常显著。如果不把C值设对你的内量子效率曲线在高注入下会偏得离谱。3.3 极化效应的取舍GaN材料有一个非常特殊的性质——强自发极化和压电极化。因为GaN是纤锌矿结构不是对称的闪锌矿结构所以沿着c轴方向天然存在一个内建电场。在InGaN/GaN量子阱中这个极化电场会导致量子阱的能带倾斜电子和空穴波函数在空间上分离降低辐射复合效率——这就是所谓的量子限制斯塔克效应QCSE。仿真时极化效应可以在模型中显式加入极化电荷也可以忽略。我的建议是如果你的研究重点是Micro-LED的量子阱结构优化必须加极化模型否则你的能带结构和载流子分布算出来都是错的。特别要注意的是Micro-LED的尺寸小应力状态和传统LED不一样——比如侧壁刻蚀之后应力释放的程度会影响压电极化的大小。这一块在仿真里需要结合实际工艺条件去设定没有统一的答案。4. 从建结构到跑通的完整实操流程讲完模型选取的思路我把跑通一个GaN基Micro-LED电学仿真的完整流程拆开来讲。下面这些步骤是我在Silvaco Atlas里反复验证过的用Crosslight或COMSOL的话思路类似只是界面和语法不同。4.1 器件结构的定义第一步是定义器件几何结构和掺杂。一个典型的InGaN/GaN蓝光Micro-LED结构从下往上大概是这样的衬底蓝宝石Al2O3仿真中通常简化为热边界条件 n-GaN层厚度2~4微米掺杂浓度1e18~3e18 cm^-3 多量子阱MQW层3~8个周期的InGaN/GaN量子阱 量子阱In组分15%~25%厚度2~4纳米 势垒GaN厚度8~15纳米 p-AlGaN电子阻挡层EBL厚度10~30纳米Al组分10%~20% p-GaN层厚度100~300纳米掺杂浓度1e17~5e19 cm^-3受限于Mg受主激活效率在Silvaco里定义结构的核心思路就是把每一层用region语句写出来。举个例子# 定义n型GaN层 region number1 materialGaN x.min0 x.max10 y.min0 y.max3.0 \ doping2e18 n.type # 定义量子阱层 region number2 materialInGaN x.min0 x.max10 \ y.min3.0 y.max3.004 \ # 阱厚4nm doping1e17 n.type \ # 本征或者轻掺杂 composition0.18 # In组分设定为18% # 定义势垒层 region number3 materialInGaN x.min0 x.max10 \ y.min3.004 y.max3.014 \ # 垒厚10nm doping1e17 n.type \ composition0.0 # In组分为0即纯GaN势垒这里唯一要注意的是很多初学者会直接写materialInN再调组分这是错的。正确做法是材料仍然定义为InGaN通过composition参数控制In的占比。仿真软件会自动按组分插值材料的禁带宽度、介电常数、有效质量等参数。4.2 网格划分的策略网格划分是最能体现经验的地方。量子阱只有几纳米厚但你上方有几百纳米的p型层、下方有数微米的n型层网格尺度跨越三个数量级。很多第一次做的人把整个结构划成均匀网格结果要么算不准要么算不动。我的做法是在量子阱和界面附近垂直方向网格间距设到0.2~0.5纳米确保每个量子阱里至少有3~5个网格点在p型层和n型层里网格可以放宽到5~10纳米在材料界面上必须做网格加密过渡不要突然从0.5纳米跳到10纳米否则泊松方程容易不收敛。水平方向的网格根据你选的维度来定。做二维模型时Micro-LED的电流扩展和侧壁效应都在这个方向上体现建议在电极边缘和侧壁附近加密网格。做三维模型的话网格数量会非常可观建议先用粗网格跑通再逐步加密看结果的收敛性不要一上来就追求最细网格。4.3 物理模型的参数配置物理模型参数的设定是仿真成败的关键。这里把我常用的参数和设置列出来方便参考但请注意实际使用时要结合你自己的实验数据校准。参数项常用值说明SRH寿命阱内5~50 nsMicro-LED尺寸越小取值应越低SRH寿命p-GaN1~10 nsp型层缺陷多寿命短俄歇复合系数C1e-30 ~ 1e-29 cm^6/s蓝光LED典型量级辐射复合系数B1e-11 ~ 1e-10 cm^3/s与In组分、温度相关极化电荷比例30%~80%实际器件中因屏蔽效应不会100%Mg受主激活能160~220 meVp-GaN激活率很低在Silvaco里辐射复合系数B一般是自动根据带间复合模型计算的不需要直接指定需要你手动设置的通常是SRH寿命和俄歇系数。注意SRH寿命不要设成一个固定值完事它和材料质量、缺陷密度甚至量子阱宽度都有关系做参数扫描时可以把寿命作为变量跑一遍看不同寿命下结果的敏感性。4.4 电极与边界条件电学仿真的边界条件相对直接——底部n电极设定为欧姆接触、顶部p电极设定为欧姆接触或肖特基接触取决于具体工艺四周侧壁通常设定为反射边界或表面复合边界。对于Micro-LED侧壁边界条件的设定需要特别注意。真实工艺中侧壁刻蚀损伤会导致表面态密度很高这些表面态会成为非辐射复合中心。做二维仿真时你可以通过设置侧壁表面的表面复合速度来模拟这种效应。Silvaco里通过interface语句设定interface x.min0 x.max0 y.min0 y.max5 \ \n.surf1e12 qf0.4 # 侧壁表面电荷和表面态表面复合速度的经验值在1e4~1e6 cm/s之间。如果没做过专门的侧壁钝化工艺建议取偏大的值如果侧壁做了钝化处理可以适当降低。4.5 求解与收敛最后是求解器设置。GaN基LED仿真最常见的头疼问题就是不收敛尤其在高电流密度或强极化条件下。我自己总结的几个有效策略先用小电流比如1e-6 A起步让求解器在一个相对容易收敛的点完成初始化然后逐步增加偏压把上一次的收敛解作为下一步的初值Silvaco的solve语句就是这个思路极化效应开启时建议先关闭极化跑通再加极化电荷否则初场太强极易发散GaN材料的禁带宽度大内建电压高偏压扫描的步长不要设太大0.1~0.2V一个步长比较稳。跑通之后你会得到一系列结果文件能带图、载流子浓度分布、复合速率分布、I-V特性、IQE曲线等。到这里仿真模型才算是“建起来了”但距离“可信”还有一大步——参数校准。5. 参数校准怎么做别让仿真变成“自嗨”仿真界有一句老话所有的仿真结果都是错的只有误差大小之分。一个没人校准过的仿真模型本质上就是一个高级玩具。校准的意思就是拿模型算出来的结果去对比真实实验数据反过来调整参数直到两条曲线吻合到足够好。5.1 从I-V特性校准电学参数I-V特性是最容易获得的实验数据也是校准的第一步。把测得的正向I-V曲线和仿真结果叠在一起先看开启电压和串联电阻是否吻合。开启电压不匹配通常是p型层掺杂浓度、Mg激活能或者极化电荷量设置不对。极化电荷会改变量子阱的内建电场直接影响开启电压串联电阻太高/太低检查p-GaN的载流子迁移率和激活率。Mg受主是好几百毫电子伏的深能级室温下激活率低得可怜这个参数设置不准串联电阻肯定不对低电流区域的曲线形状反映SRH复合的影响。如果实验曲线在低电流段比仿真掉得快说明SRH寿命设大了横向调小就行。5.2 用IQE曲线拟合复合系数IQE内量子效率曲线是你校准复合模型的核心参照。测试方法一般是变温PL或者变电流EL加积分球测光功率拿到数据后和仿真算出的IQE曲线对比。校准的核心思路是分开拟合IQE曲线峰值的位置和高度主要由A/B比决定缺陷复合vs辐射复合而高电流端的下降斜率主要由C项和载流子泄漏决定。所以调参数时可以分步走先调A/B比让曲线峰值高度和位置覆盖实验值。注意峰值位置对应的电流密度——如果实验峰值在1 A/cm²而仿真在10 A/cm²说明A值相对B值偏小也就是缺陷复合占比不够再调C值让高电流端的droop斜率匹配实验曲线。如果仿真曲线在高电流端下降得太急说明C值偏大需要调小如果高电流端怎么调C都拟合不上就要考虑开启载流子泄漏通道——在仿真模型里加一个p型层的泄漏路径。这里给一个很实用的建议把参数扫描做成一组平行仿真别一个一个单独调。因为A、B、C三个参数之间有耦合性——调大B值相当于相对减小了A和C的影响。我用脚本跑一个A值扫描、一个C值扫描然后把结果叠在实验数据上观察哪组组合整体拟合度最好。5.3 校准过程中容易忽略的坑有几个坑是特别容易踩的提醒一下温度参数没设对。很多人在仿真里默认300K但实际EL测试中芯片是自热的在大电流下结温可能升到350K以上。温度会影响B系数随温度升高而下降和载流子泄漏不修正温度高电流端拟合永远有偏差极化电荷设成100%。理论上的极化电荷值在实际器件中会因为缺陷、位错和自由载流子屏蔽而部分抵消设成100%会让开启电压和IQE曲线同时失真。经验值是设30%~70%具体要看你的实验数据只校准了IQE曲线一个维度。IQE曲线的拟合可以有无数个参数组合同时成立这就是所谓“参数非唯一性”。想减小非唯一性最好多点验证同时拟合I-V曲线、IQE曲线、光谱特性如果有的话三个维度都对了模型的可信度才高。6. 从仿真回到器件设计量子阱结构优化的实战思路模型校准完之后就可以开始干正事了——用仿真来指导器件设计。这里分享一个我用校准好的模型做量子阱结构优化的完整思路也算是整个建模工作的最终价值所在。6.1 量子阱数量与厚度的权衡Micro-LED的量子阱设计有一个经典的权衡问题阱的数量多了有源区体积大可以在低电流下储存更多载流子但注入均匀性变差而且每个阱分到的电流密度下降阱的数量少了有源区体积小高电流下容易饱和。我的仿真实践结论是对于目标工作电流密度在10~100 A/cm²的Micro-LED5~8个周期的多量子阱结构相对合理。阱宽方面3纳米左右比4纳米更有优势因为更薄的阱可以让电子空穴波函数重叠更好极化场导致的能带倾斜影响变小辐射复合效率更高。但阱太薄又会带来In组分波动增大的问题实际外延生长也很难控制所以3纳米是个不错的折中。这个过程特别适合用仿真去做参数扫描固定其他条件不变把阱数从3个扫描到10个把阱宽从2.5纳米扫描到5纳米看哪个组合在目标电流密度下IQE最高。每个结构跑一条IQE曲线半天时间能跑完这要放在实验上一组MOCVD生长加器件制备就是一两个月的事情。6.2 电子阻挡层EBL的优化EBL是p型层和量子阱之间的一层AlGaN作用是阻止电子溢流到p型层参与非辐射复合。在Micro-LED的高电流密度工作条件下EBL的作用尤为关键。仿真中EBL的核心优化参数有三个Al组分、厚度、掺杂浓度。Al组分越高对电子的势垒越高但空穴的注入也会更困难——这是根本矛盾。厚度增厚可以增强阻挡效果但会增加串联电阻。我的经验做法是先扫Al组分10%~25%确定一个能有效抑制电子溢流同时不显著恶化空穴注入的区间再看厚度10~30纳米。仿真结果中有两个判据特别重要一是p型层里的电子浓度如果高电流下p型层里积累了显著浓度的电子说明电子溢流很严重二是EBL/p-GaN界面处的价带形状如果出现明显的空穴势垒需要调整掺杂策略来改善空穴注入。6.3 侧壁钝化结构的仿真验证最后说一个Micro-LED独有、也是最影响光效的设计点——侧壁钝化。前文提到侧壁表面态导致非辐射复合实际工程中一般通过沉积SiO2、Al2O3等钝化层来减少表面态。仿真在这个问题上能做的事情是定量评估不同表面复合速度下的效率损失帮你判断钝化工艺做到什么程度才够用。举个例子我在一个20微米尺寸的Micro-LED模型中把侧壁表面复合速度从1e3设到1e6 cm/s做扫描结果发现当表面复合速度超过1e5 cm/s后内量子效率大幅下降尤其在低电流密度下更明显。这意味着如果你的钝化工艺只能把表面复合速度压到1e5 cm/s以上那效率损失会很严重工艺改进方向就明确了。另外仿真还可以帮你评估钝化层的厚度和介电常数对光提取效率的影响——不过这部分就属于光学仿真的范畴了需要和FDTD工具联用。7. 实测复现中的教训与工具链建议整个建模过程走下来除了物理和参数方面的经验我在实际工具链和复现流程上也有一些体会整理在这里供参考。7.1 软件选型的个人建议目前主流的Micro-LED仿真工具有Silvaco TCAD、Crosslight APSYS和COMSOL Multiphysics各有侧重Silvaco Atlas上手相对平缓语法清晰在LED电学仿真方面有大量前例可循适合初学者和快速迭代结构扫描Crosslight APSYS在III-V族光电器件领域地位稳固对极化模型、量子阱自发辐射等物理细节的处理更精细但上手门槛高License费用也更高COMSOL强在多物理场耦合如果你要同时做电、热、力、光分析COMSOL更灵活但它并不是特别为半导体器件优化过的求解漂移扩散问题时精度和效率不如前两者。我的建议是如果主要做器件电学仿真和结构优化从Silvaco或APSYS起步更合适如果要做光电热力多物理场耦合可以关注COMSOL。不用同时学三个工具把其中一个吃透建模思路是相通的。7.2 参数库的积累仿真的一个重要习惯是要建立自己的参数库。不同文献中报告的SRH寿命、俄歇系数、极化电荷比例差异极大每次建模都从零开始查文献效率极低而且容易出错。我在建模初期就把常用的InGaN/GaN材料参数整理成了一份表格包括不同In组分下的禁带宽度用带隙弯曲系数修正、折射率、有效质量、迁移率、热导率等每篇文献的数据都标注来源和适用条件。后续所有仿真都从这套参数库出发再根据目标器件的具体工艺做修正效率提升非常明显。7.3 跑仿真前先想清楚要什么结果最后一条建议可能听起来像废话但我见过太多人栽在这上面开始仿真之前先写清楚你想回答什么问题、需要输出哪些结果、判断标准是什么。是做结构对比那就要设计好对照组。是拟合实验那就要先整理好实验数据的格式和单位。是给工艺做趋势预测那参数扫描的范围和步长要先规划好。没有目标的仿真最容易时间黑洞——调了一周参数最后发现算出来的结果根本没法回答最初的问题。每次跑仿真之前拿十分钟把目标、输入、输出写下来比什么都管用。我自己现在是这么操作的先在文档里列一张小表格写着“问题是什么、用什么指标回答、哪些参数要变、哪些参数固定、最终交付什么图”想清楚了再动手。这一套流程走完之后做Micro-LED仿真对我来说才真正变成一件可控的事。希望在读的你也能沿着这条路径少走几个我已经替你踩过的坑。本文还有配套的精品资源点击获取