上周我把一個平面型功率二極管的Silvaco仿真結(jié)果和實(shí)測數(shù)據(jù)擺在一起反向擊穿電壓整整差了四十多伏。我一開始以為是雪崩模型參數(shù)沒選對srh、auger、impact selb換了又換收斂曲線倒是挺漂亮擊穿電壓紋絲不動。最后實(shí)在沒轍打開tonyplot把網(wǎng)格顯示勾上才終于看明白問題冶金結(jié)附近的網(wǎng)格線間距拉得太大峰值電場被數(shù)值插值直接“抹平”雪崩擊穿自然被推遲到了更高的電壓。那次之后我做Silvaco二極管仿真網(wǎng)格優(yōu)化時第一件事永遠(yuǎn)是想清楚哪里該密、哪里該稀而不是急著把x.mesh、y.mesh鋪滿整個結(jié)構(gòu)。這篇就專門聊聊二極管仿真網(wǎng)格的參數(shù)設(shè)置思路和實(shí)操技巧適合剛接觸Silvaco TCAD的學(xué)生也適合正在調(diào)器件結(jié)構(gòu)、總感覺仿真結(jié)果“差口氣”的工程師。1. 網(wǎng)格密度決定仿真精度的底層邏輯欠密區(qū)域怎么毀掉擊穿電壓1.1 一次反向擊穿電壓差40V的排查問題不在模型在網(wǎng)格先說我那個十微米寬、結(jié)深零點(diǎn)五微米的平面功率二極管。器件結(jié)構(gòu)很簡單P陽極區(qū)、N-漂移區(qū)、N襯底摻雜分布也是教科書級別的。按理說這種結(jié)構(gòu)在Silvaco里跑擊穿電壓不應(yīng)該出問題但結(jié)果就是和實(shí)測對不上。排查過程一開始完全是朝著錯誤方向走的。我懷疑是不是雪崩模型參數(shù)里某個系數(shù)設(shè)得不對于是把models語句里的impact selb參數(shù)來回調(diào)毫無變化又懷疑網(wǎng)格是不是太密導(dǎo)致收斂到偽解試著整體粗化了一倍網(wǎng)格結(jié)果擊穿電壓反而更離譜。直到某次我隨手把網(wǎng)格顯示打開才注意到一個問題在結(jié)深0.5微米附近y方向網(wǎng)格間距只有0.05微米左右而耗盡區(qū)在這個位置的電場梯度變化極其劇烈整個空間電荷區(qū)里滿打滿算只有三層網(wǎng)格。三層網(wǎng)格是什么概念在Silvaco的有限體積法框架下電場和載流子濃度都是通過網(wǎng)格節(jié)點(diǎn)之間的插值來逼近的。網(wǎng)格太少相當(dāng)于用兩點(diǎn)連直線去描述一條指數(shù)衰減曲線峰值被平均掉、陡峭被抹平。雪崩擊穿依賴的是峰值電場和碰撞電離率的空間積分峰值一旦被歪曲擊穿電壓自然跑偏。后來我把冶金結(jié)附近的網(wǎng)格間距從0.05微米一步步壓到0.002微米擊穿電壓從虛高的182V回落到實(shí)測的141V附近誤差直接縮到3%以內(nèi)。這給我留下的印象太深了——二極管仿真里網(wǎng)格永遠(yuǎn)不是背景板它本身就是計算精度的一部分。1.2 有限體積法對網(wǎng)格的一般要求物理量梯度決定網(wǎng)格步長Silvaco ATLAS求解的是泊松方程加上電子、空穴連續(xù)性方程這組偏微分方程組。數(shù)值方法上它先把連續(xù)的器件區(qū)域離散成一個個四邊形網(wǎng)格單元然后把偏微分方程在每個網(wǎng)格節(jié)點(diǎn)附近轉(zhuǎn)成代數(shù)方程再用Newton迭代或Gummel迭代來求解。這套方法的精度上限并不取決于你的物理模型有多高級而取決于你給求解器提供的物理量分辨率。說白了網(wǎng)格就是仿真世界的“像素”物理量在網(wǎng)格節(jié)點(diǎn)之間有值但在網(wǎng)格內(nèi)部是不知道的。求解器默認(rèn)用線性插值去近似網(wǎng)格點(diǎn)之間的物理量分布如果某個區(qū)域里電場或載流子濃度變化劇烈而網(wǎng)格太少插值出來的結(jié)果一定會失真。我習(xí)慣拿拍照做類比同一個場景一千像素和五千萬像素拍出來的細(xì)節(jié)完全不是一個層級。PN結(jié)耗盡區(qū)里空間電荷密度、電場強(qiáng)度、載流子濃度在亞微米尺度內(nèi)可能變化好幾個數(shù)量級你就是把整個器件的網(wǎng)格都加密到0.1微米如果結(jié)區(qū)那幾層網(wǎng)格不夠密結(jié)果照樣是糊的。更關(guān)鍵的一點(diǎn)是Silvaco的網(wǎng)格間距需要遠(yuǎn)小于物理量發(fā)生顯著變化的特征長度。對二極管來說這個特征長度由摻雜濃度梯度和耗盡區(qū)寬度共同決定實(shí)際中往往是幾十納米到幾百納米的量級。很多初學(xué)者按照“整個結(jié)構(gòu)撒一遍1微米網(wǎng)格”的思路去做反向特性勉強(qiáng)能看但正向壓降、擊穿電壓這些對梯度敏感的參數(shù)就會明顯失真。1.3 二極管里需要重點(diǎn)加密的三類區(qū)域既然網(wǎng)格分辨率不是均勻撒的那就要搞清楚二極管仿真中哪幾個區(qū)域才是決定精度的“關(guān)鍵戰(zhàn)場”。以我做過的高壓整流二極管、PiN二極管、TVS結(jié)構(gòu)為例下面這幾個位置是網(wǎng)格優(yōu)化時必須優(yōu)先照顧的。第一冶金結(jié)和它兩側(cè)的耗盡區(qū)。這里是空間電荷密度和電場梯度最大的地方也是雪崩擊穿的發(fā)生地。雪崩擊穿電壓對峰值電場的依賴是非線性的網(wǎng)格哪怕只是把峰值拉低百分之幾擊穿電壓都可能偏差十幾伏。所以冶金結(jié)附近必須用最密的網(wǎng)格且結(jié)兩側(cè)都要加密不只是結(jié)線本身那一個位置。第二陽極和陰極接觸面附近的重?fù)诫s濃度過渡區(qū)。這個區(qū)域載流子復(fù)合率、注入效率都很高直接決定正向壓降和接觸電阻的仿真精度。雖然這里的網(wǎng)格不需要像結(jié)區(qū)那樣夸張但如果完全放任不管正向IV曲線的低壓段會出現(xiàn)明顯的非線性畸變。第三反偏時耗盡區(qū)可能擴(kuò)展到的邊界位置。對高壓功率二極管來說反偏電壓升高時耗盡區(qū)會向低摻雜漂移區(qū)深處不斷擴(kuò)展。很多人以為只要結(jié)區(qū)加密就夠了結(jié)果耗盡區(qū)邊緣跑到稀疏網(wǎng)格區(qū)電場峰值在數(shù)值上被“截斷”擊穿點(diǎn)位置也跟著漂移。網(wǎng)格控制線必須提前鋪到耗盡區(qū)可能到達(dá)的范圍而不是只圍著自己的舒適區(qū)轉(zhuǎn)。2. Silvaco網(wǎng)格構(gòu)建參數(shù)逐個拆解從space.mult到ratio2.1 最小的網(wǎng)格框架x.mesh和y.mesh到底在做什么Silvaco ATLAS里構(gòu)建網(wǎng)格的核心語句就是x.mesh和y.mesh兩條。很多人第一次接觸時會以為寫一個y.mesh就夠因?yàn)槎O管結(jié)構(gòu)看起來是“一維”的——摻雜只沿深度方向變化。但實(shí)際上Silvaco的數(shù)值求解器要求網(wǎng)格必須是完整的二維四邊形網(wǎng)格哪怕x方向只有兩個位置你也得把x方向的兩條網(wǎng)格線寫出來。少了x方向deckbuild會直接報錯連求解都進(jìn)不去。先看一個最小可運(yùn)行的網(wǎng)格框架go atlas mesh space.mult1.0 x.mesh location0 spacing0.5 x.mesh location10 spacing0.5 y.mesh location0 spacing0.05 y.mesh location0.5 spacing0.002 ratio1.2 y.mesh location2.0 spacing0.05 ratio1.4 y.mesh location6.5 spacing0.3這個例子里x方向在0和10微米各放一條網(wǎng)格線間距0.5微米表示器件寬度10微米y方向從表面0微米開始到6.5微米結(jié)束中間在0.5微米處用0.002微米的間距加密。location是網(wǎng)格線的絕對位置spacing是該位置附近的網(wǎng)格間距ratio則是從這條網(wǎng)格線向外延伸時網(wǎng)格間距的放大倍數(shù)。需要理解的一點(diǎn)是x.mesh和y.mesh并不是“畫格子”而是定義每一根網(wǎng)格線的位置和該處的局部間距。求解器拿到這些控制點(diǎn)之后會自動在相鄰網(wǎng)格線之間生成過渡網(wǎng)格最終的網(wǎng)格疏密分布就是由這些控制點(diǎn)的spacing和ratio共同決定的。2.2 spacing、location、ratio的關(guān)系與配合寫法這三個參數(shù)里spacing是局部精度的直接控制旋鈕location是加密位置的錨點(diǎn)ratio是控制網(wǎng)格從密到疏過渡方式的“坡度”。舉個例子在0.5微米處的冶金結(jié)上我設(shè)置spacing0.002微米ratio1.15。那么離開這條網(wǎng)格線后下一層網(wǎng)格間距大概就是0.002乘以1.15等于0.0023微米再下一層是0.00264微米以此類推按等比數(shù)列逐步放大。大約十幾層之后間距就能從0.002微米平滑過渡到亞微米量級而不至于在某一個位置突然跳變。這里有一個經(jīng)驗(yàn)區(qū)間值得強(qiáng)調(diào)ratio一般取1.1到1.5比較穩(wěn)妥。比例太大會導(dǎo)致網(wǎng)格間距增長過快過渡區(qū)域出現(xiàn)“斷層”數(shù)值上容易產(chǎn)生非物理的震蕩比例太小則網(wǎng)格要經(jīng)過很多層才能變疏白白浪費(fèi)大量網(wǎng)格點(diǎn)求解變慢。我自己的習(xí)慣是在結(jié)區(qū)附近用1.1到1.2在離開敏感區(qū)后的過渡帶用1.3到1.5。還有一個容易忽略的細(xì)節(jié)location可以配合注釋在code里標(biāo)注“這一條是結(jié)區(qū)邊界”“這一條是空間電荷區(qū)擴(kuò)展邊界”后來回頭調(diào)試時能省不少時間。TCAD仿真一旦搞起來deck文件里網(wǎng)格語句可能十幾條沒有標(biāo)注簡直是災(zāi)難。2.3 一維PN結(jié)也要寫完整的二維網(wǎng)格框架剛才說過Silvaco不允許只有一維網(wǎng)格存在。很多初學(xué)者第一次寫二極管仿真下意識只在y方向?qū)懥嗣苊苈槁榈木W(wǎng)格線x方向一個location都沒有結(jié)果跑起來直接報“missing mesh line in x direction”一類的錯誤。正確的做法是哪怕你做的確實(shí)是深度方向的一維PN結(jié)x方向也要定義至少兩條網(wǎng)格線間距隨意但別太小。x方向兩條網(wǎng)格線之間的區(qū)域Silvaco會把它當(dāng)成一個在x方向無限延展的平板結(jié)構(gòu)來算物理上就是平行于表面的二極管家底這完全沒問題。但如果你仿真的是帶結(jié)終端保護(hù)環(huán)、場限環(huán)、或者臺面刻蝕結(jié)構(gòu)的功率二極管情況就不一樣了。這類器件的電場分布是真正的二維分布終端區(qū)域在x方向上有明顯變化這時x方向的網(wǎng)格也要像y方向一樣認(rèn)真控制。尤其是在保護(hù)環(huán)和主結(jié)之間的區(qū)域電場線彎曲得很厲害x方向網(wǎng)格不能放得太粗否則終端擊穿電壓就是仿給鬼看的。3. 二極管網(wǎng)格參數(shù)配置實(shí)操一套可直接套用的結(jié)構(gòu)方案3.1 典型PN功率二極管結(jié)構(gòu)參數(shù)與網(wǎng)格控制線布局理論扯了這么多直接上一套我能復(fù)現(xiàn)的實(shí)操配置。下面這個例子是常見的高壓PN二極管結(jié)構(gòu)P陽極區(qū)0.5微米峰值摻雜1e19每立方厘米N-漂移區(qū)5微米摻雜1e15每立方厘米N襯底區(qū)1微米摻雜1e18每立方厘米x方向?qū)挾?0微米。網(wǎng)格控制線的布局邏輯是這樣的表面接觸區(qū)和中段漂移區(qū)用寬松間距冶金結(jié)兩側(cè)用最密的網(wǎng)格離開結(jié)區(qū)后按比例逐步放大襯底與漂移區(qū)界面再補(bǔ)一層加密。完整的網(wǎng)格定義段如下# 全局網(wǎng)格間距乘數(shù)0.8表示整體加密20% mesh space.mult0.8 x.mesh location0 spacing0.5 x.mesh location10 spacing0.5 # y方向表面接觸區(qū) y.mesh location0 spacing0.02 # 結(jié)區(qū)上側(cè)漸變加密0.45微米處開始 y.mesh location0.45 spacing0.002 ratio1.15 # 冶金結(jié)精確位置0.5微米最密網(wǎng)格 y.mesh location0.5 spacing0.001 ratio1.1 # 結(jié)區(qū)下側(cè)漸疏0.6微米處開始 y.mesh location0.6 spacing0.005 ratio1.2 # 空間電荷區(qū)內(nèi)部 y.mesh location1.2 spacing0.02 ratio1.35 # 漂移區(qū)中段 y.mesh location2.5 spacing0.08 ratio1.4 # 襯底界面附近 y.mesh location5.9 spacing0.2 ratio1.5 y.mesh location6.5 spacing0.1為什么要在0.45微米和0.6微米各放一條控制線而不是只在0.5微米放一條因?yàn)镻N結(jié)兩側(cè)的電場分布都不均勻單靠結(jié)線上一個密點(diǎn)向上和向下的網(wǎng)格過渡會太快達(dá)不到需要的分辨率。兩條控制線相當(dāng)于在結(jié)區(qū)兩側(cè)各建了一個“加密緩沖區(qū)”讓網(wǎng)格間距從最密值平滑地爬升到正常水平。這套網(wǎng)格跑下來整個結(jié)構(gòu)大約會產(chǎn)生數(shù)千個網(wǎng)格點(diǎn)。對于做幾微米到幾十微米深度的二極管完全夠用單次擊穿求解在桌面工作站上也就是幾十秒的事。3.2 從工藝仿真導(dǎo)入網(wǎng)格后的補(bǔ)刀REFINEBOX局域加密很多二極管仿真并不是在ATLAS里直接搭結(jié)構(gòu)而是先用ATHENA跑完全部工藝步驟再通過結(jié)構(gòu)文件導(dǎo)到ATLAS做電學(xué)仿真。工藝仿真的網(wǎng)格會根據(jù)氧化、擴(kuò)散、注入這些步驟自動調(diào)整但它的網(wǎng)格分布是以“工藝熱歷史”為核心邏輯的和電學(xué)仿真的需求不完全一致。這時候REFINEBOX語句就能派上用場。它允許你在指定的矩形區(qū)域內(nèi)重新按步長生成網(wǎng)格非常適合在導(dǎo)入網(wǎng)格之后做定點(diǎn)補(bǔ)密。比如上面那個結(jié)構(gòu)如果導(dǎo)入后發(fā)現(xiàn)結(jié)區(qū)網(wǎng)格還是不夠密可以在ATLAS里追加一條refinebox x.min0 x.max10 y.min0 y.max0.65 x.step0.5 y.step0.002這條語句的含義是在x方向0到10微米、y方向0到0.65微米的盒形區(qū)域內(nèi)x方向步長0.5微米y方向步長0.002微米重新生成網(wǎng)格。它不會動區(qū)域外已有的網(wǎng)格因此非常適合“局部打補(bǔ)丁”。使用時唯一要注意的是step不要一上來就給太小。0.002微米在0.65微米深度內(nèi)意味著三百多層網(wǎng)格如果區(qū)域還很大網(wǎng)格點(diǎn)數(shù)量會瞬間爆炸求解時間直接從秒級跳到小時級。我一般先把step給到預(yù)估值的兩三倍在tonyplot里看效果不夠再逐步收緊。3.3 用space.mult做多檔掃描快速找敏感區(qū)域space.mult是網(wǎng)格定義段開頭那個全局參數(shù)作用是整體縮放所有網(wǎng)格間距。space.mult小于1時網(wǎng)格變密大于1時網(wǎng)格變粗等于1時保持原位。這個參數(shù)是網(wǎng)格優(yōu)化過程中性價比最高的工具沒有之一。當(dāng)你在掃描不同的摻雜濃度、漂移區(qū)厚度、結(jié)深時不需要為每個候選結(jié)構(gòu)手工調(diào)整十幾條y.mesh語句只要把space.mult從1.0換成0.7或者0.5就能在保持網(wǎng)格分布形態(tài)的前提下整體加密。我自己的操作習(xí)慣是先以space.mult1.0粗掃一批候選結(jié)構(gòu)選出兩三個性能指標(biāo)接近目標(biāo)的結(jié)構(gòu)再用space.mult0.7和0.5分別精算對比擊穿電壓、正向?qū)▔航档汝P(guān)鍵參數(shù)的變化趨勢。如果0.7和0.5的結(jié)果相差很小說明網(wǎng)格密度已經(jīng)夠用了如果差距還很大那就要回到單個結(jié)構(gòu)的網(wǎng)格布局上做定點(diǎn)加密而不是繼續(xù)依賴全局參數(shù)。4. 網(wǎng)格收斂性測試與自適應(yīng)網(wǎng)格怎么判斷網(wǎng)格真的夠了4.1 Grid Test的標(biāo)準(zhǔn)流程三組密度、兩個指標(biāo)、一次結(jié)論網(wǎng)格優(yōu)化的最終目的是在誤差可接受的前提下盡量少耗時。怎么判斷誤差可接受唯一靠譜的辦法就是做網(wǎng)格收斂性測試也就是社區(qū)里常說的Grid Test。標(biāo)準(zhǔn)操作分三步。第一步保持器件結(jié)構(gòu)和物理模型完全不變用同樣的space.mult跑三組仿真推薦1.0、0.7、0.5。第二步記錄每個方案下的關(guān)鍵電學(xué)指標(biāo)對二極管來說主要看反向擊穿電壓和正向?qū)ǖ膲航狄部梢约右粋€零偏結(jié)電容作為輔助參考。第三步比較后兩組結(jié)果。如果space.mult0.7和0.5之間的指標(biāo)差異小于1%就認(rèn)為網(wǎng)格已經(jīng)收斂space.mult0.7那一檔可以直接作為最終配置。我拿一次實(shí)際的掃參記錄舉個例子具體數(shù)值因器件而異但趨勢判斷方法是通用的網(wǎng)格方案space.mult網(wǎng)格點(diǎn)數(shù)反向擊穿電壓正向?qū)▔航祮吸c(diǎn)求解時間方案A1.0約2600158.3V0.67V約30s方案B0.7約5200162.6V0.64V約2min方案C0.5約10200163.1V0.63V約9min方案B和方案C的擊穿電壓只差0.3%正向壓降差約1.5%這時就可以判斷網(wǎng)格已經(jīng)進(jìn)入收斂平臺方案B足以作為精度可靠的配置。如果方案A到方案B的變化很大而方案B到方案C又收斂了說明默認(rèn)網(wǎng)格確實(shí)不夠解決方案就是B檔。4.2 自適應(yīng)網(wǎng)格MGOALS的正確打開方式與適用邊界Silvaco還提供了一類自適應(yīng)網(wǎng)格手段核心是METHOD語句里加ADAPT參數(shù)再用MGOALS語句定義誤差目標(biāo)。典型寫法如下method adapt newton mgoals n.points8000 curr.error1e-3它的工作邏輯是在迭代求解過程中不斷評估局部物理量的誤差如果發(fā)現(xiàn)某個區(qū)域誤差超出MGOALS設(shè)定的目標(biāo)就自動插入新的網(wǎng)格點(diǎn)然后再接著往收斂方向求解。從概念上講這東西很“智能”等于讓求解器自己決定哪里該加密。但我的實(shí)際體驗(yàn)是自適應(yīng)網(wǎng)格在二極管仿真里不能無腦用。它適合的場景是做單點(diǎn)精修比如已經(jīng)用手動網(wǎng)格跑通了一個器件想看某個偏置點(diǎn)附近峰值電場到底是多少可以用adaptive模式局部細(xì)化驗(yàn)證。它不適合的場景是全偏置IV掃描因?yàn)槊總€偏置點(diǎn)都可能觸發(fā)不同的自適應(yīng)細(xì)化網(wǎng)格點(diǎn)數(shù)量和分布隨偏置改變最終IV曲線會出現(xiàn)不少細(xì)小的臺階反而不利于結(jié)果分析。另外自適應(yīng)網(wǎng)格會顯著放大求解時間。MGOALS里的n.points設(shè)得太大求解器會為了滿足誤差目標(biāo)而瘋狂加網(wǎng)格點(diǎn)有時候一個偏置點(diǎn)的收斂要等十分鐘還不如靜下心來把手動網(wǎng)格的關(guān)鍵位置調(diào)好。我的建議是手動網(wǎng)格永遠(yuǎn)是主力MGOALS只是補(bǔ)充校驗(yàn)手段。4.3 網(wǎng)格點(diǎn)數(shù)、迭代次數(shù)與求解時間的實(shí)際賬本很多人對網(wǎng)格數(shù)量和求解時間之間的關(guān)系沒有概念以為網(wǎng)格多了就是慢一點(diǎn)而已。實(shí)際上對于ATLAS這類求解非線性偏微分方程組的程序網(wǎng)格點(diǎn)數(shù)N和求解時間不是線性關(guān)系而是接近O(N^1.5)到O(N^2)的指數(shù)放大。原因并不復(fù)雜。Newton迭代過程中每個迭代步都要組裝和求解一個規(guī)模與網(wǎng)格點(diǎn)數(shù)相關(guān)的Jacobian矩陣矩陣規(guī)模大了單次求解的時間會超線性增長。與此同時網(wǎng)格多了之后求解器為了保持?jǐn)?shù)值穩(wěn)定性往往需要更多迭代步數(shù)才能收斂。兩個因素疊加在一起就會出現(xiàn)“網(wǎng)格翻一倍時間翻五倍”的情況。所以網(wǎng)格優(yōu)化的本質(zhì)不是無限加密而是把有限的網(wǎng)格預(yù)算投放到物理量梯度最劇烈的地方。同樣是一萬個網(wǎng)格點(diǎn)勻著撒滿整個器件和集中在結(jié)區(qū)加耗盡區(qū)邊緣計算出來的擊穿電壓可能差出百分之十幾。這也是我反復(fù)強(qiáng)調(diào)要手工布置網(wǎng)格控制線、而不是依賴全局加密的原因。5. 我被Silvaco網(wǎng)格坑過的幾件事與最終建議5.1 網(wǎng)格線與冶金結(jié)錯位摻雜分布出現(xiàn)鋸齒這是我最開始玩Silvaco時踩過的一個隱蔽坑。結(jié)構(gòu)定義里我把P區(qū)和N-區(qū)的region邊界設(shè)在了y0.5微米摻雜語句也對應(yīng)做了定義但y.mesh網(wǎng)格線的控制點(diǎn)卻放在了0.45和0.55微米沒有一條正好落在0.5微米上。結(jié)果就是冶金結(jié)這個摻雜突變面沒有和網(wǎng)格線對齊而是斜穿了一排網(wǎng)格單元。在數(shù)值模型里摻雜分布被離散成了跨網(wǎng)格的插值結(jié)面變成了一個“鋸齒狀”的過渡帶。仿出來的電場分布曲線在結(jié)附近出現(xiàn)周期性抖動擊穿特性的平滑度很不對勁。解決辦法說穿了很簡單網(wǎng)格控制線的location必須和region邊界、摻雜突變面的位置嚴(yán)格對齊。建議在搭結(jié)構(gòu)時先把region用變量或注釋標(biāo)清楚再回頭寫y.mesh確保結(jié)面處一定有一條網(wǎng)格線。凡是涉及材料界面、摻雜突變、電極接觸邊界網(wǎng)格線都必須老老實(shí)實(shí)落在邊界上這是Silvaco網(wǎng)格優(yōu)化的第一紀(jì)律。5.2 反向耗盡區(qū)擴(kuò)展到稀疏網(wǎng)格區(qū)擊穿點(diǎn)漂移另一個坑出現(xiàn)在做高壓二極管的時候。結(jié)構(gòu)是P陽極然后是低摻雜的N-漂移區(qū)再往下是N陰極。我一開始按照經(jīng)驗(yàn)只在結(jié)區(qū)附近加密了網(wǎng)格N-漂移區(qū)中段和靠近N界面的地方網(wǎng)格都放得很粗。仿真時發(fā)現(xiàn)一個很奇怪的現(xiàn)象反向擊穿位置不在PN結(jié)附近反而出現(xiàn)在了N區(qū)深處而且離N襯底界面越近越嚴(yán)重。后來把耗盡區(qū)寬度一算才明白這個二極管在接近擊穿電壓時耗盡區(qū)已經(jīng)擴(kuò)展到了漂移區(qū)中后段而那一帶的網(wǎng)格間距大概有0.08微米電場梯度完全無法被準(zhǔn)確分辨峰值被強(qiáng)行拉到了網(wǎng)格相對密一點(diǎn)的襯底邊界。從那以后我做高壓二極管網(wǎng)格布局時第一件先做的事情就是估算反偏條件下耗盡區(qū)的最大擴(kuò)展范圍然后在該范圍內(nèi)部署數(shù)條漸變的網(wǎng)格控制線。哪怕估算結(jié)果有點(diǎn)偏差也要在可能結(jié)東的區(qū)域留下幾層中等密度的網(wǎng)格作為緩沖。對低摻雜漂移區(qū)100V附近的反偏耗盡區(qū)可能擴(kuò)展到十幾微米網(wǎng)格控制線必須跟著鋪過去而不是只圍著結(jié)區(qū)打轉(zhuǎn)。5.3 判斷網(wǎng)格方案是否合理的三板斧說了這么多踩坑經(jīng)歷最后總結(jié)三個我每次做完網(wǎng)格都過一遍的檢查習(xí)慣。第一板斧在tonyplot里把網(wǎng)格顯示打開和摻雜分布曲線疊在一起看。重點(diǎn)網(wǎng)格加密區(qū)是否覆蓋了摻雜梯度大的區(qū)域、結(jié)面是否被網(wǎng)格線“騎”過、有沒有明顯的網(wǎng)格間距跳變。第二板斧把電場分布和碰撞電離率分布畫出來看峰值位置是否落在加密區(qū)內(nèi)部。如果電場峰值落在網(wǎng)格稀疏區(qū)那不管IV曲線收斂成什么樣結(jié)論都不可信。第三板斧也是最快的驗(yàn)證方法跑兩組不同space.mult的關(guān)鍵指標(biāo)對比。如果1.0和0.7的結(jié)果差異巨大那沒有任何理由相信1.0的結(jié)果是對的如果0.7和0.5差異很小則說明網(wǎng)格深度已經(jīng)夠用。這三步走完一版網(wǎng)格方案基本就能讓人放心了。前陣子一個朋友問我Silvaco和另一款TCAD工具哪個好。我回他說工具都有各自的優(yōu)勢但不管你用哪款網(wǎng)格問題都是繞不過去的換工具解決不了網(wǎng)格沒布置好的問題。就說Silvaco的二極管仿真網(wǎng)格參數(shù)設(shè)置看著不起眼實(shí)際卻是決定結(jié)果可信度的那根定海神針。把結(jié)區(qū)加密、把耗盡區(qū)邊界鋪滿、用space.mult做收斂測試這套方法論用熟之后你會發(fā)現(xiàn)自己對TCAD仿真的理解又深了一層。