---+ 3日目(デバイスシミュレーション) %TOC% ---++ 前回のSimple2Dの例題を走らせてみる。 まずは、Simple2Dのプロジェクトを開いてください。 <verbatim>cd ~/work/Silicon/TCAD/Sentaurus/tutorials/Simple2D swb . &</verbatim> node 27だけ実行してあったと思います。 その構造に相当するIVシミュレーション(node 45)も実行して前回同様、Inspectorを使ってIV_n45_des.plt を開いてみてください。 anodeの !OuterVoltage をx軸に、anodeのTotalCurrentをy軸に指定してプロットを描くとしたの図(左)のようになります。 最初の例題で見たグラフとずいぶん違う形になっていますが、これは130V付近で電場の集中による雪崩増幅が起こったことを示唆しています。 一般的に、p-n 接合に大きな逆方向電圧がかけられると,接合は降伏現象(Break down)を起こし,大電流が流れます。 Y軸をログスケールにするともう少しわかりやすいかもしれません。上の方真ん中付近にある logYというボタンを押してみてください。 下の図(右)のようになります。(対数グラフになっているので縦軸は絶対値になっているので注意! (0,0)は右下です。) <img alt="iv_break.png" height="400" src="%ATTACHURLPATH%/iv_break.png" /> <img alt="iv_break.png" height="400" src="%ATTACHURLPATH%/iv_break_ly.png" /> 130V付近までは依然見た様な電流電圧特性を示しますが、その後急に増えていることがわかります。 この降伏(Break down)が起こる電圧をBreak down voltageと呼びます。 Breakを起こしたときの電場分布と電子の密度分布を見てみましょう。 node 45で右クリック、 svisual で n45_des.tdrを開く。 (scaling : X to Y ratioを10くらいにして表面領域を拡大してみてください。) 電場分布 (Electric Field)や、電子密度(Electron Density)、衝突イオン化(Impact Ionization)などを見てみるとガードリングの付近で非常に電場が集中していることがわかります。 一般的に電子の持つ運動エネルギーがシリコンのバンドギャップを超えたとき、雪崩増幅が起こる。この電圧はシリコンであれば約300kV/cm以上である。(出典: SM.Sze. 『Semiconductor Device』3.7章「高電界効果」より) 電場分布をy軸に沿った軸でカットしてみると3e5V/cmを超えていることがわかります。 <img alt="field.png" height="863" src="%ATTACHURLPATH%/field.png" width="1480" /> ---++ !Simple2D のIVシミュレーションのコマンド まず、デバイスシミュレーションに関しては膨大なオプションが準備されています。マニュアルは全1574ページ!! もちろん今日だけで全部は説明できませんし、おそらく必要ないものがほとんどです。辞書的にマニュアルを使いながら理解してみてください。 IV_des.cmdという名前のファイルを見てみましょう。以下のコンポーネントからできています。 <verbatim>File {...} Electrode {...} Physics {...} CurrentPlot {...} Solve {...} Plot {...} </verbatim> ---+++ File section まず初めに入力ファイルや出力ファイルを指定する。 入力ファイルは Gridとparameter。 Gridでは sdeで作った構造体を読み込みます。Work Benchでdevice simulationを行おうとしているnodeの左にあるSimple2Dというsdeを実行したノードの番号は @node|Simple2D@で取得できます。 parameterはシミュレーションに使う細かいパラメータを設定しています。通常はsdevice.parというファイルが使われます。 出力ファイルはCurrentとPlot が結果になります。 outputでlogを出力します。 Current はすべてのシミュレーション(下のSolveで指定)の電気的な結果の値が詰まったファイルを出力します。 @plot@でnode番号 + _des.pltという名前のファイルができます。 Plotはシミュレーション終了時の電気的な変数の情報を構造体の位置情報を込みで出力します。 <verbatim>File { Grid = "n@node|Simple2D@_msh.tdr" Current = "@plot@" Plot = "@tdrdat@" output = "@log@" parameter = "@parameter@" }</verbatim> ---+++ Electrode section ここでは「電極」を指定します。IVなどのシミュレーションをする際にどこをGNDにどこを電位Vにするなどを指定する必要があります。 ただし、構造体をsdeで作った際に電極になる場所を作っておく必要があります。前回の sdegeo:define-contact-set を参照。 !Simple2D_des.cmdのなかでストリップ電極数(cathode1,2,3)+裏面(anode)の4つの電極を作って、それぞれを 0.0 Vに設定します。 eRecVelocity や hRecVelocityは再結合モビリティー(単位は cm/s)。 そのほかに直列に抵抗を入れたり、電流や電荷の境界条件や初期条件を入れることができる。 <verbatim> Electrode { !( for {set i 1} {$i<=@NumStripes@} {incr i} { puts "\{name = \"cathode${i}\" voltage = 0.0 eRecVelocity=1e7 hRecVelocity=1e7 \}" } )! {name = "anode" voltage = 0.0 eRecVelocity=1e7 hRecVelocity=1e7 } }</verbatim> ---+++ Physics section 物理モデルを指定する。この部分は来週、再来週にまとめて話します。 <verbatim> Physics { ##AreaFactor=1e8 ##Fermi } Physics (material="Silicon") { Mobility( DopingDep (Unibo) HighFieldsat (GradQuasiFermi) ) Recombination( SRH(DopingDep TempDep) Auger ##hAvalanche(UniBo) eAvalanche(UniBo) hAvalanche eAvalanche ) EffectiveIntrinsicDensity (OldSlotboom) Traps( ( name="state1" acceptor conc=@<fluence*1.613>@ Level FromConductionBand EnergyMid =0.42 eXsection=2E-15 hXsection=2E-14 ##eJfactor=1.0 hJfactor=1.0 ) ( name="state2" acceptor conc=@<fluence*100.0>@ Level FromConductionBand EnergyMid =0.46 eXsection=5E-15 hXsection=5E-14 ##eJfactor=1.0 hJfactor=1.0 ) ( name="state3" donor conc=@<fluence*0.9>@ Level FromValenceBand EnergyMid =0.36 eXsection=2.5E-14 hXsection=2.5E-15 ##eJfactor=1.0 hJfactor=1.0 ) ) }</verbatim> ---+++ !CurrentPlot section <br />CurrentPlot {<br />eLifeTime(Maximum(material="Silicon"))<br />hLifeTime(Maximum(material="Silicon"))<br />eAvalanche(Maximum(material="Silicon"))<br />hAvalanche(Maximum(material="Silicon"))<br />} ---+++ Math section <br />Math {<br /><br /> ##CDensityMin=1e-100<br /> Extrapolate<br /> Derivatives<br /> Avalderivatives<br /> Digits=5<br /> Notdamped=1000<br /> Iterations=20<br /> RelerrControl <br /> ErrRef (electron)=1e6<br /> ErrRef (hole)=1e6<br /> RhsMin =1e-15<br /> <br /> eMobilityAveraging=ElementEdge <br /> hMobilityAveraging=ElementEdge <br /> ParallelToInterfaceInBoundaryLayer (-ExternalBoundary)<br /> <br /> ##Method=Blocked<br /> ##SubMethod=ILS (set=1)<br /> ##Method=ILS (set=1)<br /> Method=Pardiso<br /> ILSrc = "<br /> set(1) { <br /> iterative( gmres(100), tolrel=1e-8, tolunprec=1e-4, tolabs=0, maxit=200 ); <br /> preconditioning( ilut(0.0001,-1), left ); <br /> ordering( symmetric=nd, nonsymmetric=mpsilst ); <br /> options( compact=yes, verbose=0, refineresidual=0 ); <br /> };<br /> "<br /><br /> <br />number_of_threads=4<br /><br />} ---+++ Solve Section <br /><br />Solve {<br /><br /> Poisson<br /> plugin { Poisson Electron Hole }<br /> Coupled { Poisson Electron Hole }<br /><br /> NewCurrent ="IV_" <br /> <br /> Quasistationary ( DoZero <br /> MaxStep =1 MinStep =1e-8 InitialStep =1e-5<br /> Increment=1.6 Decrement=4.0<br /> Goal { Name="anode" Voltage=-120} <br /> BreakCriteria {Current (Contact = "anode" minval = -1e-7 )}<br /> )<br /> { Coupled { Poisson Electron Hole } } <br /><br /> Quasistationary ( DoZero <br /> MaxStep =1 MinStep =1e-8 InitialStep =1e-5<br /> Increment=1.6 Decrement=4.0<br /> Goal { Name="anode" Voltage=@Vop@} <br /> BreakCriteria {Current (Contact = "anode" minval = -1e-7 )}<br /> )<br /> { Coupled { Poisson Electron Hole } } <br />system("")<br /><br />}<br /> ---+++ Plot Section <br />Plot {<br /> <br /> Current/Vector <br /> eCurrent/Vector<br /> hCurrent/Vector<br /> eDensity<br /> hDensity<br /> ElectricField /Vector<br /> Potential<br /> CurrentPotential <br /> DopingConcentration <br /> eMobility<br /> hMobility<br /> DonorConcentration <br /> AcceptorConcentration <br /> AvalancheGeneration <br /> <br /> eAvalanche hAvalanche<br /> eLifeTime hLifeTime<br /><br /><br />}<br /><br /><br /><br /> <br /><br /><br /><br /> -- %USERSIG{AtlasjSilicon - 2020-05-21}%
Attachments
Attachments
Topic attachments
I
Attachment
History
Action
Size
Date
Who
Comment
png
field.png
r1
manage
102.4 K
2020-05-21 - 08:53
AtlasjSilicon
png
iv_break.png
r1
manage
8.7 K
2020-05-21 - 03:11
AtlasjSilicon
png
iv_break_ly.png
r1
manage
8.5 K
2020-05-21 - 03:15
AtlasjSilicon
This topic: Main
>
TCADSimulationInformation
>
Tutorial4thTCADWorkshop
>
Tutorial6thTCADDay3
Topic revision: r4 - 2020-05-21 - AtlasjSilicon
Copyright © 2008-2025 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki?
Send feedback