JP3930223B2 - 半導体デバイスのシミュレーション方法およびシミュレーションプログラムを記録した記録媒体 - Google Patents

Description

【０００１】
【発明の属する技術分野】
本発明は半導体デバイスのシミュレーション方法およびシミュレーションプログラムを記録した記録媒体に関する。
【０００２】
【従来の技術】
一般に、ＬＳＩ産業の発展を支える半導体デバイスの微細化に伴い、半導体デバイスの動作および性能の確認のために半導体デバイスのシミュレーションが行われている。
【０００３】
従来のシミュレーション方法においては、ドナー濃度にドナーのイオン化率を乗じたものが電子濃度であり、アクセプター濃度にアクセプターのイオン化率を乗じたものが正孔濃度であるとして計算していた。
【０００４】
【発明が解決しようとする課題】
半導体デバイスの微細化に伴って、半導体基板に形成されるｐｎ接合間の距離が急激に狭くなってきている。このため、短チャネル効果やラッチアップ現象が半導体デバイスの動作に与える影響が顕著になり、ｐｎ接合が半導体デバイスのシミュレーションに重大な影響を与える。
【０００５】
これを図５に示したｎＭＯＳＦＥＴの断面図と、図６に示したバンド図を用いて、詳しく説明する。まず、図５に斜線で示した部分１０では、基板２に注入されたアクセプターと拡散層４ａ，４ｂに分布しているドナーの割合が比較的１に近い。界面近傍では不純物濃度を低く抑えるので、少し離れたこの部分１０において量子多体効果の影響が大きくなりシミュレーション結果に影響を及ぼす。この領域１０内で、拡散層４ａ側では、ドナーが優勢な不純物であり、アクセプターが劣勢な不純物である。反対に、この領域１０内で、基板２側では、アクセプターが優勢な不純物で、ドナーが劣勢な不純物である。ここで重要なことは、ｐｎ接合付近では、劣勢不純物の影響が無視出来ないということである。以下、図６を用いて、劣勢不純物の影響を具体的に説明する。図６では、ドナーを優勢とし、アクセプターを劣勢としたｎ型の場合を示してある。この図から直ちに分かるように、ドナーとアクセプターが直ぐ近くにある場合、ドナーの電子をアクセプターが受け入れて、安定化する。このとき、ドナーもアクセプターも両方ともイオン化している。即ち、従来のシミュレーション技術では、ドナー濃度にドナーのイオン化率を乗じたものが電子濃度であると近似されていたが、アクセプターの影響を考慮に入れると、この近似が成り立たないのである。ドナーのイオン化率ｒｄは、ドナーが電子を放出する過程のイオン化率ｒｄ_１とアクセプターに渡す過程のイオン化率ｒｄ_２の和で表される。
【０００６】
このとき、電子濃度は、ｒｄ_１にドナー濃度を乗じたものになる。こうして、電子濃度の計算結果は、劣勢不純物を無視した従来の計算結果に比べ、低くなることが判る。一方、正孔濃度は、価電子が伝導帯に立ち上がる過程に加え、アクセプターに受容される過程も加わる為、増大する。結果として、キャリア濃度に与える影響は小さくない。また、微細化に伴って基板の不純物濃度（劣勢不純物の濃度）が高くなるにつれ、このような劣勢不純物の影響は大きくなる。上述した議論は、アクセプターが優勢でドナーが劣勢な場合にも同様に成り立つ。
【０００７】
そして、図７（ａ）に示すように、ＭＯＳＦＥＴが微細化してゲート長が短くなったり、図７（ｂ）に示すようにバイポーラ素子のベース幅が狭くなると、ｐｎ接合付近の領域１０では劣勢不純物の影響が大きくなる。劣勢不純物の影響を完全に無視した従来のデバイスシミュレーション技術においては、劣勢な不純物の影響は無視されており、ドナーとアクセプターの割合が比較的１に近いｐｎ接合付近の現実の物理現象を正確に取り込んでいない。このため従来のシミュレーションは不正確であった。
【０００８】
本発明は上記事情を考慮してなされたものであって、可及的に正確なシミュレーションを行うことのできる半導体デバイスのシミュレーション方法およびシミュレーションプログラムを記録した記録媒体を提供することを目的とする。
【０００９】
【課題を解決するための手段】
本発明による半導体デバイスのシミュレーション方法は、電位が無く、かつ電荷の中性条件が成立すると仮定し、量子多体効果を無視した場合のフェルミエネルギーを求めるステップと、前記フェルミエネルギーを用いて、量子多体効果を考慮したドナーおよびアクセプターのイオン化率を求めるステップと、を備えたことを特徴とする。
【００１０】
なお、前記イオン化率を求めるステップは、量子多体効果によるフェルミ面のシフトを計算するステップと、このシフトを用いてイオン化率を計算するステップと、を備えたことを特徴とする。
【００１１】
なお、量子多体効果を考慮した、電荷の中性条件が成立している場合のキャリアの準粒子エネルギーシフトを求めるステップと、ポアッソン方程式と電荷の連続式を同時に満たす、ポテンシャル、電子濃度、および正孔濃度を求めるステップと、ポテンシャルが存在することによる準粒子エネルギーシフトへの付加項と、電荷の連続式から来るキャリア数の平衡状態からのずれによる付加項とを電子および正孔の場合について求めるステップと、キャリアの準粒子エネルギーシフトおよびイオン化率を求めるステップと、を備えたことを特徴とする。
【００１２】
また本発明によるシミュレーションプログラムを記録した記録媒体は、電位が無く、かつ電荷の中性条件が成立すると仮定し、量子多体効果を無視した場合のフェルミエネルギーを求める手順と、前記フェルミエネルギーを用いて、量子多体効果を考慮したドナーおよびアクセプターのイオン化率を求める手順と、を備えたことを特徴とする。
【００１３】
【発明の実施の形態】
本発明による半導体デバイスのシミュレーション方法の一実施の形態を図１および図２を参照して説明する。この実施の形態のシミュレーション方法の処理手順を図１に示す。
【００１４】
まず、半導体デバイスに電位がなく、次の（１）式に示す電荷の中性条件
Ｎ_Ｄ ^＋−Ｎ_Ａ ⁻＋ｐ_０−ｎ_０＝０ … （１）
が成り立つ場合を考える。ここで、Ｎ_Ｄ ^＋はドナーのイオン数、Ｎ_Ａ ⁻はアクセプターのイオン数、ｐ_０は正孔濃度、ｎ_０は電子濃度を示す。上記電荷の中性条件は、劣勢不純物の影響を加味してバンド構造を決定するのに最も適した条件である。
【００１５】
まず図１のステップＦ１に示すように、半導体デバイスのシミュレーションすべき領域のドナー濃度Ｎ_Ｄ、アクセプター濃度Ｎ_Ａ、ドナーのイオン化エネルギε_Ｄ、アクセプターのイオン化エネルギーε_Ａ、温度Ｔ、および真性エネルギーギャップＥＧ_ｉｎｔを入力する。なお、定数となる有効伝導帯状態密度Ｎ_Ｃ、有効価電子帯状態密度Ｎ_Ｖ、およびボルツマン定数ｋ_Ｂはシミュレーションを行う前に予め与えておく。
【００１６】
次に上述の条件の下で、量子多体効果を無視した場合のフェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を求める（図１のステップ２参照）。ここでＥ_Ｆ００は量子多体効果を無視したフェルミ準位、Ｅ_Ｃ００は量子多体効果を無視した伝導帯端である。
【００１７】
上記フェルミエネルギーは以下のようにして求められる。
【００１８】
まずフェルミ・ディラック統計に従い、電子濃度ｎ_００および正孔濃度ｐ_００は、次の（２）式によって表される。
【００１９】
【数１】

ここでＮ_Ｃは有効伝導帯状態密度、Ｎ_Ｖは有効価電子帯状態密度、Ｆ_１／２はフェルミディラック積分、Ｅ_Ｖ００は量子多体効果を無視した価電子帯端である。
【００２０】
フェルミ・ディラック積分Ｆ_１／２（ａ）は次の式で表される。
【００２１】
【数２】

また価電子帯端Ｅ_Ｖ００は、量子多体効果を無視したエネルギーギャップＥＧ_ｉｎｔを用いて次の（３）式で表される。
【００２２】
Ｅ_Ｖ００＝Ｅ_Ｃ００−ＥＧ_ｉｎｔ … （３）
また、ドナーのイオン数Ｎ_Ｄ ^＋およびアクセプターのイオン数Ｎ_Ａ ⁻は、次の（４）式で表わされる。
【００２３】
Ｎ_Ｄ ^＋＝ｒｄ_００×Ｎ_Ｄ、Ｎ_Ａ ⁻＝ｒａ_００×Ｎ_Ａ … （４）
ここで、Ｎ_Ｄはドナー濃度、Ｎ_Ａはアクセプター濃度、ｒｄ_００，ｒａ_００は量子多体効果を無視したドナー、アクセプターのイオン化率を各々示す。これらのイオン化率ｒｄ_００，ｒａ_００は次の（５）式によって表される。
【００２４】
【数３】

ここでＥ_Ｄはドナー準位、Ｅ_Ａはアクセプータ準位であり、ドナーのイオン化エネルギーε_ｄ、アクセプターのイオン化エネルギーε_Ａを用いて次の（６）式で表される。ただし、
Ｅ_Ｄ＝Ｅ_Ｃ００−ε_Ｄ、Ｅ_Ａ＝Ｅ_Ｖ００＋ε_Ａ … （６）
ここで（２）式〜（６）式を（１）式に代入すれば、フェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を変数とする方程式が得られる。この方程式を、反復法または２分法等の解法を用いて逐次近似的に解き、フェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を数値的に求める。
【００２５】
次に再び図１に戻り、この求めたエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を用いて、量子多体効果によるフェルミ面のシフトｅｆ_０を計算し、量子多体効果を考慮したドナーおよびアクセプターのイオン化率ｒｄ_０（ｅｆ_０），ｒａ_０（ｅｆ_０） を計算する（図１のステップＦ３参照）。
【００２６】
まずフェルミ面のシフトｅｆ_０は以下のようにして求められる。量子多体効果の影響を次の（７）式を用いて導入する。
【００２７】
Ｅ_Ｆ−Ｅ_Ｃ＝Ｅ_Ｆ００−Ｅ_Ｃ００−ＱＰＳ_ｅｏ（ｅｆ_０）
Ｅ_Ｖ−Ｅ_Ｆ＝Ｅ_Ｖ００−Ｅ_Ｆ００−ＱＰＳ_ｈｏ（ｅｆ_０） … （７）
ここで、ＱＰＳ_ｅｏは電子の準粒子エルネギーシフト、ＱＰＳ_ｈｏは正孔の準粒子エネルギーシフトであり、量子多体効果によるフェルミ面のシフトｅｆ_０を変数に持つ形で表現することができる。また
Ｅ_Ｆは量子多体効果を考慮したフェルミ準位、
Ｅ_Ｃは量子多体効果を考慮した伝導帯端、
Ｅ_Ｖは量子多体効果を考慮した価電子帯端である。
【００２８】
この準粒子エネルギーシフトＱＰＳ_ｅｏ，ＱＰＳ_ｈｏを用いると、量子補正された電子の濃度（準粒子密度）ｎ_０および正孔の濃度（準粒子密度）ｐ_０は次の（８）式で表わされる。
【００２９】
【数４】

この（８）式から、準粒子密度ｎ_０，ｐ_０も、フェルミ面のシフトｅｆ_０の関数であることが分かる。同様に、イオン化率ｒｄ_０，ｒａ_０も次の（９）式に示すように量子補正を受け、フェルミ面のシフトｅｆ_０の関数となる。
【００３０】
【数５】

フェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）が既知なので、（７）式〜（９）式を（１）式に代入すると、フェルミ面のシフトｅｆ_０を一変数とする方程式が得られる。この方程式を、反復法または二分法等を用いて遂次近似的に解き、フェルミ面のシフトｅｆ_０を数値的に求める。これにより、ドナーおよびアクセプターの量子補正されたイオン化率ｒｄ_０（ｅｆ_０）およびｒａ_０（ｅｆ_０）を得ることができる。
【００３１】
このようにして図１のステップＦ３が終了する。また、このとき、フェルミ面のシフトｅｆ_０が求められたので、量子多体効果を考慮した、電荷の中性条件が成立している場合のキャリアの準粒子エネルギーシフトＱＰＳ_０、すなわち電子および正孔の準粒子エネルギーシフトＱＰＳ_ｅｏ（ｅｆ_０）およびＱＰＳ_ｈｏ（ｅｆ_０）をセルフコンシステントに求めることができる（図１のステップＦ５参照）。これらの和がバンドギャップナローウィングである。
【００３２】
電荷の中性条件が成立しない場合のシミュレーションを行わないならば、図１のステップＦ３の後、またはステップ４の後でシミュレーションを終了しても良い。
【００３３】
次に量子多体効果の影響を見るために、上述の計算アルゴリズムを用いて得た、量子多体効果を主因とするバンドギャップナローウィング（以下、ＢＧＮとも言う）の計算結果を図２のグラフに示す。この図２において横軸は拡散層のドナー濃度を示し、縦軸はＢＧＮを示している。この図２のグラフから分かるように、ＢＧＮはドナーとアクセプターがほぼ等しい濃度で存在するとき消滅する。これは、半導体デバイスのシミュレーションにＢＧＮを考慮に入れる上で、重要な変調を与えることになる。取分け、ＢＧＮの影響の大きいバイポーラタイプのデバイスでは重要である。
【００３４】
次に電荷の中性条件が成立しない場合のシミュレーション方法について説明する。
【００３５】
一般にキャリアの準粒子エネルギーシフトは、電子の濃度ｎ、正孔の濃度ｐ、ドナーのイオン数Ｎ_Ｄ ^＋、アクセプタのイオン数Ｎ_Ａ ⁻の関数として表わされる。なお、電子の場合をΔ_ｅ（ｎ，ｐ，Ｎ_Ｄ ^＋，Ｎ_Ａ ⁻）と表わされ、正孔の場合をΔ_ｈ（ｎ，ｐ，Ｎ_Ｄ ^＋，Ｎ_Ａ ⁻）とする。電荷の中性条件が成り立つ場合の準粒子エネルギーシフトの上述計算においては、ｎとしてｎ_０（ｅｆ_０）を、ｐとしてｐ_０（ｅｆ_０）を、Ｎ_Ｄ ^＋としてｒｄ_０（ｅｆ_０）×Ｎ_Ｄを、Ｎ_Ａ ⁻としてｒａ_０（ｅｆ_０）×Ｎ_Ａを用いていたため、準粒子エネルギーシフトＱＰＳ_０（ｅｆ_０）はフェルミ面のシフトｅｆ_０の関数となっている。
【００３６】
実際のデバイスでは、電荷の中性条件が成り立つことは少ない。電荷の輸送があれば、電荷の連続式を満たすように半導体中の各点でｎとｐは平衡状態（ｎ_０，ｐ_０）からずれた値を持つ。また、ポテンシャルΨの存在によっても平衡状態からずれる。実際的なアルゴリズムを得るには、これらを準粒子エネルギーシフトに含めるよう、上述の理論を拡張しなければならない。簡単のため、空間一次元を仮定すると、電荷の連続式（輸送方程式）は次の（１０）式のように書ける。
【００３７】
【数６】

実際のデバイスシミュレーションでは、メッシュに切られた半導体デバイスの半導体領域の各格子点上において、上述の電荷の連続式と、次の（１１）式に示すポアッソン方程式、
【数７】

を同時に満たすように数値的に計算された電子の濃度ｎ、正孔の濃度ｐ、ポテンシャルΨが求められる（図１のステップＦ５参照）。
【００３８】
ここでＥは電界であり、ポテンシャルΨの傾きに比例している。またεは半導体デバイスの誘電率、μ_ｎおよびμ_ｐは電子および正孔の各々の移動度、Ｄ_ｎおよびＤ_ｐは電子および正孔の各々の拡散係数、Ｇ_ｎおよびＧ_ｐは各々電子および正孔の生成率、Ｕ_ｎおよびＵ_ｐは電子および正孔の各々の再結合率である。
【００３９】
次にこのようにして得られたｎ，ｐ，Ψを用いて、図１のステップＦ６に示すようにポテンシャルΨが存在することによる準粒子エネルギーシフトの付加項ＱＰＳ_１（Ψ）、すなわち電子の準粒子エネルギーシフトの付加項ＱＰＳ_１ｅ（Ψ）および正孔の準粒子エネルギーシフトの付加項ＱＰＳ_１ｈ（Ψ）を次の式
【数８】

を用いて求めるとともに、電荷の連続式から来るキャリア数のずれによる付加項ＱＰＳ_２（Ψ）を次の式
【数９】

を用いて求める。ここでｑは素電荷である。
【００４０】
次に図１に示すステップＦ７に進み、キャリアの準粒子エネルギーシフトＱＰＳ＝ＱＰＳ_０＋ＱＰＳ_１＋ＱＰＳ_２およびキャリアのイオン化率ｒｄ，ｒａを求める。すなわち、電子の準粒子エネルギーシフトＱＰＳ_ｅおよび正孔の準粒子エネルギーシフトＱＰＳ_ｈは、各々
ＱＰＳ_ｅ＝ＱＰＳ _ｅ０ ＋ＱＰＳ_１ｅ＋ＱＰＳ_２ｅ
ＱＰＳ_ｈ＝ＱＰＳ _ｈ０ ＋ＱＰＳ_１ｈ＋ＱＰＳ_２ｈ
として求められる。またドナーのイオン化率ｒｄ（Ψ）およびアクセプターのイオン化率ｒａ（Ψ）は、次の式を用いて求められる。
【００４１】
【数１０】

次に図１に示すステップＦ８に進み、キャリアの準粒子エネルギーシフトＱＰＳの値およびドナーとアクセプターのイオン化率ｒｄ，ｒａの値が収束したか否かが判定され、収束しない場合にはステップＦ５に戻り、上述のことが繰り返される。また収束した場合はシミュレーションを終了する。なお、収束したか否かの判定は、計算によって求めた現在の値と１つ前の計算で求めた値との差を上記現在の値で割ったときの値の絶対値が所定値以下になったときに収束したと判定する。なお、ｒｄ_０＝ｒａ_０＝１としても良い。
【００４２】
また、上述の電荷の中性条件が成立しない場合のシミュレーションにおいては、初期条件を、ｒｄ_０＝ｒａ_０＝１００、ＱＰＳ＝０としてシミュレーションを行う。
【００４３】
以上説明したように、本実施の形態においては、劣勢不純物の影響を考慮してイオン化率を計算しているので、可及的に正確なシミュレーションを行うことができる。
【００４４】
次に、本発明による半導体デバイスのシミュレーションプログラムを記録した記録媒体の一実施の形態について図３および図４を参照して説明する。図３、図４は、本実施の形態の半導体デバイスのシミュレーションプログラムを記録した記録媒体が用いられるコンピュータシステム３０の一例を示す斜視図およびブロック図である。図３において、コンピュータシステム３０は、ＣＰＵを含むコンピュータ本体３１と、例えばＣＲＴ等の表示装置３２と、キーボードやマウス等の入力装置３３と、印刷を実行するプリンタ３４と、備えている。
【００４５】
コンピュータ本体３１は、図４に示すように、ＲＡＭより構成される内部メモリ３５と、内蔵または外付け可能なメモリユニット３６と、を備えており、メモリユニット３６としてはフレキシブルまたはフロッピディスク（ＦＤ）ドライブ３７、ＣＤ−ＲＯＭドライブ３８、ハードディスクドライブ（ＨＤ）ユニット３９が搭載されている。図３に示すように、これらのメモリユニット３６に用いられる記録媒体４０としては、ＦＤドライブ３７のスロットに挿入されて使用されるフレキシブルディスクまたはフロッピディスク（ＦＤ）４１と、ＣＤ−ＲＯＭドライブ３８に用いられるＣＤ−ＲＯＭ４２等が用いられる。
【００４６】
図３および図５に示すように、一般的なコンピュータシステムに用いられる記録媒体４０としては、ＦＤ４１やＣＤ−ＲＯＭ４２が考えられるが、本実施の形態は特に半導体デバイスのシミュレーションプログラムに関するものであるので、例えばシミュレータに内蔵させる不揮発性メモリとしてのＲＯＭチップに本発明の半導体デバイスのシミュレーションプログラムを記録させるようにしても良い。
【００４７】
本実施の形態の記録媒体４０には、図１に示すステップＦ２〜Ｆ８のアルゴリズムが記録されている。なお、電荷の中性条件が成立する場合のみしかシミュレーションを行わないならば、記録媒体４０には、少なくとも図１に示すステップＦ２及びＦ３のアルゴリズムが記録されていれば良い。
【００４８】
【発明の効果】
以上述べたように本発明によれば、可及的に正確なシミュレーションを行うことができる。
【図面の簡単な説明】
【図１】本発明による半導体デバイスのシミュレーション方法の一実施の形態の処理手順を示すフローチャート。
【図２】電荷の中立条件が成立しているときの本発明によって求めた、量子多体効果を主因とするバンドギャップナローウィングの計算結果を示すグラフ。
【図３】本発明による半導体デバイスのシミュレーションプログラムを記録した記録媒体が用いられるコンピュータシステムの一例を示す斜視図。
【図４】本発明による半導体デバイスのシミュレーションプログラムを記録した記録媒体が用いられるコンピュータシステムの一例を示すブロック図。
【図５】ｎＭＯＳＦＥＴの断面図
【図６】ｎ型半導体における劣勢不純物の影響を説明するバンド図。
【図７】劣勢不純物の影響を説明する模式図。
【符号の説明】
２ 半導体基板
４ａ，４ｂ 拡散層
６ ゲート（ゲート電極）
１０ ｐｎ接合付近の領域

Claims (4)

1. ドナーの濃度Ｎ_Ｄ、アクセプターの濃度Ｎ_Ａ、ドナーのイオン化エネルギーε_Ｄ、アクセプターのイオン化エネルギーε_Ａ、絶対温度Ｔ、真性半導体のエネルギーギャップＥＧ_ｉｎｔを入力するステップと、
   Ｅ_Ｆ００を電子及び正孔の準粒子エネルギーシフトを無視したフェルミ準位、Ｅ_Ｃ００を電子及び正孔の準粒子エネルギーシフトを無視した伝導帯端、Ｎ_Ｃを有効伝導帯状態密度、Ｎ_Ｖを有効価電子帯状態密度、Ｆ_１／２をフェルミ・ディラック積分、Ｅ_Ｖ００（＝Ｅ_Ｃ００−ＥＧ_ｉｎｔ）を電子及び正孔の準粒子エネルギーシフトを無視した価電子帯端、Ｅ_Ｄ（＝Ｅ_Ｃ００−ε_Ｄ）をドナー準位、Ｅ_Ａ（＝Ｅ_Ｖ００−ε_Ａ）をアクセプター準位、ｒｄ_００を電子及び正孔の準粒子エネルギーシフトを無視したドナーのイオン化率、ｒａ_００を電子及び正孔の準粒子エネルギーシフトを無視したアクセプターのイオン化率、Ｎ_Ｄ ^＋（＝ｒｄ_００×Ｎ_Ｄ）を電子及び正孔の準粒子エネルギーシフトを無視したドナーのイオン濃度、Ｎ_Ａ ⁻（＝ｒａ_００×Ｎ_Ａ）を電子及び正孔の準粒子エネルギーシフトを無視したアクセプターのイオン濃度としたとき、次の（１）〜（５）式
   

   

   Ｅ_Ｖ００＝Ｅ_Ｃ００−ＥＧ_ｉｎｔ （２）
   Ｎ_Ｄ ^＋＝ｒｄ_００×Ｎ_Ｄ、Ｎ_Ａ ⁻＝ｒａ_００×Ｎ_Ａ （３）
   Ｅ_Ｄ＝Ｅ_Ｃ００−ε_Ｄ、Ｅ_Ａ＝Ｅ_Ｖ００−ε_Ａ （４）
   

   

   を電荷の中性条件
   Ｎ_Ｄ ^＋−Ｎ_Ａ ⁻＋ｐ_００−ｎ_００＝０
   に代入したときのフェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を変数とする方程式を、逐次近似的に解き、前記フェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を計算するステップと、
   Ｅ_Ｆを電子及び正孔の準粒子エネルギーシフトを考慮したフェルミ準位、Ｅ_Ｃを電子及び正孔の準粒子エネルギーシフトを考慮した伝導帯端、Ｅ_ｖを電子及び正孔の準粒子エネルギーシフトを考慮した荷電子帯端とし、ＱＰＳ_ｅ０を電子の準粒子エネルギーシフト、ＱＰＳ_ｈ０を正孔の準粒子エネルギーシフト、ｅｆ_０を電子及び正孔の準粒子エネルギーシフトを考慮したフェルミ面のシフト、ｎ_０を電子及び正孔の準粒子エネルギーシフトを考慮した電子の濃度、ｐ_０を電子及び正孔の準粒子エネルギーシフトを考慮した正孔の濃度、ｒｄ_０（ｅｆ_０）を電子及び正孔の準粒子エネルギーシフトを考慮したドナーのイオン化率、ｒａ_０（ｅｆ_０）を電子及び正孔の準粒子エネルギーシフトを考慮したアクセプターのイオン化率としたとき、前記計算されたフェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）と、次の（６）〜（８）式
   Ｅ_Ｆ−Ｅ_Ｃ＝Ｅ_Ｆ００−Ｅ_Ｃ００−ＱＰＳ_ｅｏ（ｅｆ_０）
   Ｅ_Ｖ−Ｅ_Ｆ＝Ｅ_Ｖ００−Ｅ_Ｆ００−ＱＰＳ_ｈｏ（ｅｆ_０） （６）
   

   

   

   を電荷の中性条件
   Ｎ_Ｄ ^＋−Ｎ_Ａ ⁻＋ｐ_０−ｎ_０＝０
   に代入したときのフェルミ面のシフトｅｆ_０を一変数とする方程式を、逐次近似的に解き、前記フェルミ面のシフトｅｆ_０を計算するステップと、
   前記計算された前記フェルミ面のシフトｅｆ_０に基づいて、電子及び正孔の準粒子エネルギーシフトを考慮したドナーおよびアクセプターのイオン化率を求めるステップと、
   を備えたことを特徴とする半導体デバイスのシミュレーション方法。
2. 前記計算された前記フェルミ面のシフトｅｆ_０に基づいて、電子の準粒子エネルギーシフトＱＰＳ_ｅ０および正孔の準粒子エネルギーシフトＱＰＳ_ｈ０を算出するステップと、
   電荷の連続式およびポアッソン方程式を同時に満たす電子の濃度ｎ、正孔の濃度ｐ、およびポテンシャルΨを算出するステップと、
   ｑを素電荷としたとき、前記算出した電子の濃度ｎ、正孔の濃度ｐ、およびポテンシャルΨを次の式
   

   

   

   に代入し、ＱＰＳ_１ｅ、ＱＰＳ_１ｈ、ＱＰＳ_２ｅ、ＱＰＳ_２ｈを算出し、電子の準粒子エネルギーシフトＱＰＳ_ｅおよび正孔の準粒子エネルギーシフトＱＰＳ_ｈを次の式
   ＱＰＳ_ｅ＝ＱＰＳ_ｅ０＋ＱＰＳ_１ｅ＋ＱＰＳ_２ｅ
   ＱＰＳ_ｈ＝ＱＰＳ_ｈ０＋ＱＰＳ_１ｈ＋ＱＰＳ_２ｈ
   を用いて算出するステップと、
   前記算出されたＱＰＳ_ｅおよびＱＰＳ_ｈに基づいて、次の式
   

   

   を用いて、ドナーのイオン化率ｒｄ（Ψ）およびアクセプターのイオン化率ｒａ（Ψ）を求めるステップと、
   を備えたことを特徴とする請求項１記載の半導体デバイスのシミュレーション方法。
3. ドナーの濃度Ｎ_Ｄ、アクセプターの濃度Ｎ_Ａ、ドナーのイオン化エネルギーε_Ｄ、アクセプターのイオン化エネルギーε_Ａ、絶対温度Ｔ、真性半導体のエネルギーギャップＥＧ_ｉｎｔを入力する手順と、
   Ｅ_Ｆ００を電子及び正孔の準粒子エネルギーシフトを無視したフェルミ準位、Ｅ_Ｃ００を電子及び正孔の準粒子エネルギーシフトを無視した伝導帯端、Ｎ_Ｃを有効伝導帯状態密度、Ｎ_Ｖを有効価電子帯状態密度、Ｆ_１／２をフェルミ・ディラック積分、Ｅ_Ｖ００（＝Ｅ_Ｃ００−ＥＧ_ｉｎｔ）を電子及び正孔の準粒子エネルギーシフトを無視した価電子帯端、Ｅ_Ｄ（＝Ｅ_Ｃ００−ε_Ｄ）をドナー準位、Ｅ_Ａ（＝Ｅ_Ｖ００−ε_Ａ）をアクセプター準位、ｒｄ_００を電子及び正孔の準粒子エネルギーシフトを無視したドナーのイオン化率、ｒａ_００を電子及び正孔の準粒子エネルギーシフトを無視したアクセプターのイオン化率、Ｎ_Ｄ ^＋（＝ｒｄ_００×Ｎ_Ｄ）を電子及び正孔の準粒子エネルギーシフトを無視したドナーのイオン濃度、Ｎ_Ａ ⁻（＝ｒａ_００×Ｎ_Ａ）を電子及び正孔の準粒子エネルギーシフトを無視したアクセプターのイオン濃度としたとき、次の（１）〜（５）式
   

   

   Ｅ_Ｖ００＝Ｅ_Ｃ００−ＥＧ_ｉｎｔ （２）
   Ｎ_Ｄ ^＋＝ｒｄ_００×Ｎ_Ｄ、Ｎ_Ａ ⁻＝ｒａ_００×Ｎ_Ａ （３）
   Ｅ_Ｄ＝Ｅ_Ｃ００−ε_Ｄ、Ｅ_Ａ＝Ｅ_Ｖ００−ε_Ａ （４）
   

   

   を電荷の中性条件
   Ｎ_Ｄ ^＋−Ｎ_Ａ ⁻＋ｐ_００−ｎ_００＝０
   に代入したときのフェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を変数とする方程式を、逐次近似的に解き、前記フェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）を計算する手順と、
   Ｅ_Ｆを電子及び正孔の準粒子エネルギーシフトを考慮したフェルミ準位、Ｅ_Ｃを電子及び正孔の準粒子エネルギーシフトを考慮した伝導帯端、Ｅ_ｖを電子及び正孔の準粒子エネルギーシフトを考慮した荷電子帯端とし、ＱＰＳ_ｅ０を電子の準粒子エネルギーシフト、ＱＰＳ_ｈ０を正孔の準粒子エネルギーシフト、ｅｆ_０を電子及び正孔の準粒子エネルギーシフトを考慮したフェルミ面のシフト、ｎ_０を電子及び正孔の準粒子エネルギーシフトを考慮した電子の濃度、ｐ_０を電子及び正孔の準粒子エネルギーシフトを考慮した正孔の濃度、ｒｄ_０（ｅｆ_０）を電子及び正孔の準粒子エネルギーシフトを考慮したドナーのイオン化率、ｒａ_０（ｅｆ_０）を電子及び正孔の準粒子エネルギーシフトを考慮したアクセプターのイオン化率としたとき、前記計算されたフェルミエネルギー（＝Ｅ_Ｆ００−Ｅ_Ｃ００）と、次の（６）〜（８）式
   Ｅ_Ｆ−Ｅ_Ｃ＝Ｅ_Ｆ００−Ｅ_Ｃ００−ＱＰＳ_ｅｏ（ｅｆ_０）
   Ｅ_Ｖ−Ｅ_Ｆ＝Ｅ_Ｖ００−Ｅ_Ｆ００−ＱＰＳ_ｈｏ（ｅｆ_０） （６）
   

   

   

   を電荷の中性条件
   Ｎ_Ｄ ^＋−Ｎ_Ａ ⁻＋ｐ_０−ｎ_０＝０
   に代入したときのフェルミ面のシフトｅｆ_０を一変数とする方程式を、逐次近似的に解き、前記フェルミ面のシフトｅｆ_０を計算する手順と、
   前記計算された前記フェルミ面のシフトｅｆ_０に基づいて、電子及び正孔の準粒子エネルギーシフトを考慮したドナーおよびアクセプターのイオン化率を求める手順と、
   を備えたことを特徴とするコンピュータによって半導体デバイスのシミュレーションを行うためのシミュレーションプログラムを記録した記録媒体。
4. 前記計算された前記フェルミ面のシフトｅｆ_０に基づいて、電子の準粒子エネルギーシフトＱＰＳ_ｅ０および正孔の準粒子エネルギーシフトＱＰＳ_ｈ０を算出する手順と、
   電荷の連続式およびポアッソン方程式を同時に満たす電子の濃度ｎ、正孔の濃度ｐ、およびポテンシャルΨを算出するステップと、
   ｑを素電荷としたとき、前記算出した電子の濃度ｎ、正孔の濃度ｐ、およびポテンシャルΨを次の式
   

   

   

   に代入し、ＱＰＳ_１ｅ、ＱＰＳ_１ｈ、ＱＰＳ_２ｅ、ＱＰＳ_２ｈを算出し、電子の準粒子エネルギーシフトＱＰＳ_ｅおよび正孔の準粒子エネルギーシフトＱＰＳ_ｈを次の式
   ＱＰＳ_ｅ＝ＱＰＳ_ｅ０＋ＱＰＳ_１ｅ＋ＱＰＳ_２ｅ
   ＱＰＳ_ｈ＝ＱＰＳ_ｈ０＋ＱＰＳ_１ｈ＋ＱＰＳ_２ｈ
   を用いて算出する手順と、
   前記算出されたＱＰＳ_ｅおよびＱＰＳ_ｈに基づいて、次の式
   

   

   を用いて、ドナーのイオン化率ｒｄ（Ψ）およびアクセプターのイオン化率ｒａ（Ψ）を求める手順と、
   を備えたことを特徴とする請求項３記載のコンピュータによって半導体デバイスのシミュレーションを行うためのシミュレーションプログラムを記録した記録媒体。

Images