(57)【要約】【課題】ワイドギャップ半導体内部で起こる再結合過程を精査、検討し、精度が高く、適応範囲の広い暗電流を短時間で提供すること。【解決手段】半導体装置内部のキャリア輸送を解析するシミュレーションを行うコンピュータ装置が実行するシミュレーションの計算処理方法であって、ポアソン方程式と電流連続の式を自己無撞着に解き、静電ポテンシャル分布とキャリア密度を求めるステップと、真性キャリア密度と、キャリア寿命とから再結合率を算出するステップとを備え、印加電圧が第１のある値より低いとき、キャリア寿命をキャリア寿命の本来値より短く設定するよう変更して再結合率を算出し、印加電圧の増大に伴い、キャリア寿命を順次引き上げ、キャリア寿命の本来値まで引き上げる計算を実行するよう制御され、第１のある値は降伏電界の三分の二であることを特徴とする。【選択図】図１
2014107277.tif 000002
【特許請求の範囲】
【請求項１】
半導体装置内部のキャリア輸送を解析するシミュレーションを行うコンピュータ装置が実行するシミュレーションの計算処理方法であって、
ポアソン方程式と電流連続の式を自己無撞着に解き、静電ポテンシャル分布とキャリア密度を求めるステップと、
真性キャリア密度と、キャリア寿命とから再結合率を算出するステップと
を備え、印加電圧が第１のある値より低いとき、前記キャリア寿命を前記キャリア寿命の本来値より短く設定するよう変更して前記再結合率を算出し、前記印加電圧の増大に伴い、前記キャリア寿命を順次引き上げ、前記キャリア寿命の前記本来値まで引き上げる計算を実行するよう制御され、前記第１のある値は降伏電界の三分の二であることを特徴とする計算処理方法。
【請求項２】
請求項１に記載される計算処理方法において、
前記再結合率を算出するステップにおいて、前記キャリア寿命のかわりに、準位密度、捕獲断面積、速度のいずれか一つまたは二つの積、または全ての積を用い、前記準位密度、前記捕獲断面積、前記速度のそれぞれの本来値より大きく設定するよう変更して前記再結合率を算出し、前記印加電圧の増大に伴い、前記準位密度、前記捕獲断面積、前記速度を順次引き下げ、前記準位密度、前記捕獲断面積、前記速度の前記それぞれの本来値まで引き上げる計算を実行するよう制御されたことを特徴とする計算処理方法。
【請求項３】
請求項１または２に記載される計算処理方法において、
前記印加電圧が第２のある値より高いとき、前記静電ポテンシャル分布とキャリア密度を求めるステップ毎に、異なる捕獲準位間のエネルギー差に素電荷量を乗算し、乗算結果をボルツマン定数と絶対温度の積で除算し、除算結果の指数関数を算出し、真性キャリア密度に前記除算結果の指数関数を乗算するステップをさらに備え、前記第２のある値は降伏電界の三分の一であることを特徴とする計算処理方法。
【請求項４】
請求項１または２に記載される計算処理方法において、
前記印加電圧が前記ある値より低い領域では捕獲準位を禁止帯のほぼ中央に設定し、前記印加電圧の増大に伴い、該因子の値を所定の値に引き戻し計算を実行するよう制御されたことを特徴とする計算処理方法。
【請求項５】
請求項１乃至３のいずれか１項に記載される計算処理方法において、前記印加電圧の増大に伴い、前記キャリア寿命を順次引き上げ、または、前記準位密度、前記捕獲断面積、前記速度を順次引き下げる計算処理過程において、計算された電流値を判断因子として用いることを特徴とする計算処理方法。
【請求項６】
請求項５に記載される計算処理方法において、前記判断因子である前記電流値に対し、次に計算された電流値が前記電流値を超えた場合、順次引き上げる前記キャリア寿命の変動幅をＮ倍、または、順次引き下げる前記準位密度、前記捕獲断面積、前記速度の前記変動幅をＮ分の１とするよう制御されたことを特徴とする計算処理方法。
【請求項７】
半導体装置内部のキャリア輸送をコンピュータに求めさせるための計算処理方法であって、
少なくとも、
静電ポテンシャル分布、およびキャリア分布を求めるポアソン方程式と電流連続の式を解く手段と、
算出されたキャリアの密度と再結合捕獲プロセスにおける捕獲準位、準位密度、捕獲断面積、速度、および真性キャリア密度から再結合率を算出する手段を有し、
印加電圧が低い領域では前記準位密度、前記捕獲断面積、および前記速度のいずれか一つまたは選択された複数の因子の積の値を所定の値より大きく設定するよう変更して計算し、前記印加電圧の増大に伴い、該因子の値を順次引き下げ、所定の値まで引き下げる計算を実行するよう制御されたことを特徴とする計算処理方法。
【請求項８】
請求項７に記載される計算処理手段において前記再結合率を前記真性キャリア密度と寿命より算出する手段を有し、
前記印加電圧が低い領域では前記寿命を所定の値より短く設定し、前記印加電圧の増大に伴い、該因子の値を順次引き上げ、所定の値まで引き上げる計算を実行するよう制御されたことを特徴とする計算処理方法。
【請求項９】
請求項７または８に記載される計算処理手段において前記捕獲準位、前記準位密度、前記捕獲断面積、前記速度、および前記真性キャリア密度または前記寿命と前記真性キャリア密度から前記再結合率を算出する手段を有し、
異なる捕獲準位間のエネルギー差に素電荷量を乗算し、乗算結果をボルツマン定数と絶対温度の積で除算し、除算結果の指数関数を算出し、前記真性キャリア密度に該因子を乗算する手段を有することを特徴とする計算処理方法。
【請求項１０】
請求項９で記載される計算処理手段において、
電界強度に依存して、前記真性キャリア密度に該因子を乗算する手段を有することを特徴とする計算処理方法。
【請求項１１】
請求項７または８に記載される計算処理手段において前記捕獲準位、前記準位密度、前記捕獲断面積、前記速度、および前記真性キャリア密度から前記再結合率を算出する手段を有し、
前記印加電圧が低い領域では前記捕獲準位を禁止帯のほぼ中央に設定し、前記印加電圧の増大に伴い、該因子の値を所定の値に引き戻し計算を実行するよう制御されたことを特徴とする計算処理方法。
【請求項１２】
請求項７乃至１０のいずれかに記載される計算処理手段において前記印加電圧の増大に伴い、該因子の値を順次引き上げ、または、引き下げる計算処理過程において、計算された電流値を判断因子として用いることを特徴とする計算処理方法。
【請求項１３】
請求項１２に記載される計算処理方法において、前記判断因子である前記電流値に対し、前記計算された電流値が判断電流値を超えた場合、順次引き上げ、または、引き下げられる因子の変動幅をＮ倍またはＮ分の１とするよう制御されたことを特徴とする計算処理方法。
【請求項１４】
請求項６または１３に記載される計算処理方法において、変動幅因子Ｎを２とすることを特徴とする計算処理方法。
【請求項１５】
請求項１乃至１４のいずれかに記載される計算処理手段において、４倍精度またはそれ以上の精度によるコンパイラーを用いることを特徴とする計算処理方法。
【請求項１６】
請求項１２に記載される計算処理手段において、前記判断因子である前記電流値として、４倍精度コンパイル計算については判断電流密度を１平方センチメータ当り１０のマイナス２０乗アンペア程度とすることを特徴とする計算処理方法。
【請求項１７】
請求項５または１２に記載される計算処理方法において、前記判断因子である前記電流値として、２倍精度コンパイル計算については判断電流密度を１平方センチメータ当り１マイクロアンペア程度とすることを特徴とする計算処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は半導体装置内部の電子又は正孔の少なくとも一方、あるいは両方の輸送問題を解析する半導体デバイスシミュレーション計算処理方法に関する。
【背景技術】
【０００２】
半導体装置の応用分野は拡大の一途にあり、従来のＬａｒｇｅ−Ｓｃａｌｅ−Ｉｎｔｅｇｒａｔｉｏｎ（ＬＳＩ）のみならず、電力機器制御を可能とするハイパワー応用まで開発が精力的に進められている。ハイパワーデバイスに要求される課題は大電流化と高耐圧化である。後者については、禁制帯幅の広い「ワイドギャップ」半導体が優位である。これは、降伏電圧を決めるイオン化係数が禁制帯と逆比例して小さくなるからである。近年、これまで半導体材料として主流であったシリコン（Ｓｉ）に比べ禁制帯幅の大きい炭化ケイ素（ＳｉＣ）や窒化ガリウム（ＧａＮ）と言った新しい材料が注目されている。
【０００３】
新規材料を用いた半導体デバイスを効率的に設計するには、計算機を駆使したシミュレーション手法が有効である。これには、電位・電界分布を決めるポアソン方程式と電子・正孔分布を決める電流連続の式を自己無撞着に解く方法が適している。まずポアソン方程式を解くことにより電界強度を知ることができ、これによりイオン化率を定量的に算出可能となる。続いて、電子電流と正孔電流成分に該イオン化率を乗算する事により、雪崩増倍率を把握可能となる。
【０００４】
高耐圧化を図るデバイス構造として、基板不純物と逆導電性の電気的に浮遊なｎ型（またはｐ型）高濃度領域を設置する構造、いわゆる「ガードリング」と呼ばれる構造が採用される。電気的浮遊領域の電位はデバイス内部を流れる極微小電流による電位降下（ＩＲドロップ）により自動的に決定されるものである。
【０００５】
ｎ型領域中のｐ型領域、またはｐ型領域中のｎ型領域を浮遊電位とした所謂「ガードリング」構造の場合、ｐｎ接合を流れる電流は逆方向飽和電流と呼ばれる極微小電流（通常、暗電流と呼ばれる）で、シリコンのようなバンドギャップが１ｅＶ程度の材料では１準位の捕獲モデルから逆方向飽和電流の大きさを算出可能である。算出方法については非特許文献１や非特許文献２のなかで記載されている。
【０００６】
しかし、近年注目されているワイドギャップ半導体材料の場合、例えばＳｉＣやＧａＮのバンドギャップは約３ｅＶ程度ある。エネルギーギャップの大きさを３ｅＶとすると、真性キャリア密度は１立方センチメートル当たり１０のマイナス７乗程度（１０-7ｃｍ-3）となる。ｎｐ積一定の関係を考慮すると、キャリア密度の変化領域は５０桁を超えることになる。計算機の計算精度を表す２倍精度、４倍精度のコンパイルをもってしても５０桁の数値計算精度を確保することは困難である。８倍精度と言う方法も考えられるが、計算時間・計算機記憶容量からみて、非現実的方法論である。また真性キャリア密度が１０-7ｃｍ-3となることから、再結合率はシリコンの場合と比べて１７桁も小さな値となる。この為、実用的２倍精度コンパイルの場合、あるいは４倍精度コンパイルを用いても、ｐｎ接合逆バイアスにおける微小電流、所謂「暗電流」を従来の方法で算出することが困難である。「暗電流」の算出が困難になると、ガードリングと呼ばれる領域の浮遊電位はデバイス内部の微小電流による電位降下により決定されるので、浮遊電位そのものが求解困難となる。電位未定により数値計算は収束せず、シミュレーションそのものが事実上、実行困難となる。
【先行技術文献】
【非特許文献】
【０００７】
【非特許文献１】D. K. Schroder, “The concept of generation and recombination lifetimes in semiconductors,”(半導体内部での生成・再結合過程における寿命の概念) IEEE Trans. Electron Devices, vol. ED-29, pp. 1336-1338, Aug. 1982
【非特許文献２】山口  憲、冨澤一隆著「非平衡電子輸送論」（アドバンスソフト出版事業部）978-4-9903316-3-4
【非特許文献３】山口  憲 著「パワーデバイス用シミュレータの技術課題」（アドバンスソフト出版事業部）アドバンスシミュレーション、vol. 8, 2011
【発明の概要】
【発明が解決しようとする課題】
【０００８】
そこで、本発明は上記従来の問題点を解決すべくワイドギャップ半導体内部で起こる再結合過程を精査・検討し、その物理的性質を明らかにすることにより、合理的で精度が高く、且つ、適用範囲の広い「暗電流」を短時間で提供することを課題とする。
【課題を解決するための手段】
【０００９】
上記課題を解決して本発明の目的を達成するため、本発明による請求項１の計算処理方法は、半導体装置内部のキャリア輸送を解析するシミュレーションを行うコンピュータ装置が実行するシミュレーションの計算処理方法であって、ポアソン方程式と電流連続の式を自己無撞着に解き、静電ポテンシャル分布とキャリア密度を求めるステップと、真性キャリア密度と、キャリア寿命とから再結合率を算出するステップとを備え、印加電圧が第１のある値より低いとき、キャリア寿命をキャリア寿命の本来値より短く設定するよう変更して再結合率を算出し、印加電圧の増大に伴い、キャリア寿命を順次引き上げ、キャリア寿命の本来値まで引き上げる計算を実行するよう制御され、第１のある値は降伏電界の三分の二であることを特徴とする。
【００１０】
本発明による請求項２の計算処理方法は、請求項１に記載される計算処理方法において、再結合率を算出するステップにおいて、キャリア寿命のかわりに、準位密度、捕獲断面積、速度のいずれか一つまたは二つの積、または全ての積を用い、準位密度、捕獲断面積、速度のそれぞれの本来値より大きく設定するよう変更して再結合率を算出し、印加電圧の増大に伴い、準位密度、捕獲断面積、速度を順次引き下げ、準位密度、捕獲断面積、速度のそれぞれの本来値まで引き上げる計算を実行するよう制御されたことを特徴とする。
【００１１】
本発明による請求項３の計算処理方法は、請求項１または２に記載される計算処理方法において、印加電圧が第２のある値より高いとき、静電ポテンシャル分布とキャリア密度を求めるステップ毎に、異なる捕獲準位間のエネルギー差に素電荷量を乗算し、乗算結果をボルツマン定数と絶対温度の積で除算し、除算結果の指数関数を算出し、真性キャリア密度に除算結果の指数関数を乗算するステップをさらに備え、第２のある値は降伏電界の三分の一であることを特徴とする。
【００１２】
本発明による請求項４の計算処理方法は、請求項１または２に記載される計算処理方法において、印加電圧が前記ある値より低い領域では捕獲準位を禁止帯のほぼ中央に設定し、印加電圧の増大に伴い、該因子の値を所定の値に引き戻し計算を実行するよう制御されたことを特徴とする。
【００１３】
本発明による請求項５の計算処理方法は、請求項１乃至３のいずれか１項に記載される計算処理方法において、印加電圧の増大に伴い、キャリア寿命を順次引き上げ、または、準位密度、捕獲断面積、速度を順次引き下げる計算処理過程において、計算された電流値を判断因子として用いることを特徴とする。
【００１４】
本発明による請求項６の計算処理方法は、請求項５に記載される計算処理方法において、判断因子である電流値に対し、次に計算された電流値が判断電流値を超えた場合、順次引き上げるキャリア寿命の変動幅をＮ倍、または、順次引き下げる準位密度、捕獲断面積、速度の変動幅をＮ分の１とするよう制御されたことを特徴とする。
【００１５】
本発明による請求項７の計算処理方法は、半導体装置内部のキャリア輸送をコンピュータに求めさせるための計算処理方法であって、少なくとも、静電ポテンシャル分布、およびキャリア分布を求めるポアソン方程式と電流連続の式を解く手段と、算出されたキャリアの密度と再結合プロセスにおける捕獲準位、準位密度、捕獲断面積、速度、および真性キャリア密度から再結合率を算出する手段を有し、印加電圧が低い領域では準位密度、捕獲断面積、および速度のいずれか一つまたは選択された複数の因子の積の値を所定の値より大きく設定するよう変更して計算し、印加電圧の増大に伴い、該因子の値を順次引き下げ、所定の値まで引き下げる計算を実行するよう制御されたことを特徴とする。
【００１６】
本発明による請求項８の計算処理方法は、請求項７に記載される計算処理手段において再結合率を真性キャリア密度と寿命より算出する手段を有し、印加電圧が低い領域では寿命を所定の値より短く設定し、印加電圧の増大に伴い、該因子の値を順次引き上げ、所定の値まで引き上げる計算を実行するよう制御されたことを特徴とする計算処理方法。
【００１７】
本発明による請求項９の計算処理方法は、請求項７または８に記載される計算処理手段において捕獲準位、準位密度、捕獲断面積、速度、および真性キャリア密度または寿命と真性キャリア密度から再結合率を算出する手段を有し、異なる捕獲準位間のエネルギー差に素電荷量を乗算し、乗算結果をボルツマン定数と絶対温度の積で除算し、除算結果の指数関数を算出し、真性キャリア密度に該因子を乗算する手段を有することを特徴とする。
【００１８】
本発明による請求項１０の計算処理方法は、請求項９で記載される計算処理手段において、電界強度に依存して、真性キャリア密度に該因子を乗算する手段を有することを特徴とする。
【００１９】
本発明による請求項１１の計算処理方法は、請求項７または８に記載される計算処理手段において捕獲準位、準位密度、捕獲断面積、速度、および真性キャリア密度から再結合率を算出する手段を有し、印加電圧が低い領域では捕獲準位を禁止帯のほぼ中央に設定し、印加電圧の増大に伴い、該因子の値を所定の値に引き戻し計算を実行するよう制御されたことを特徴とする。
【００２０】
本発明による請求項１２の計算処理方法は、請求項７乃至１０のいずれかに記載される計算処理手段において印加電圧の増大に伴い、該因子の値を順次引き上げ、または、引き下げる計算処理過程において、計算された電流値を判断因子として用いることを特徴とする。
【００２１】
本発明による請求項１３の計算処理方法は、請求項１２に記載される計算処理方法において、判断因子である電流値に対し、計算された電流値が判断電流値を超えた場合、順次引き上げ、または、引き下げられる因子の変動幅をＮ倍またはＮ分の１とするよう制御されたことを特徴とする。
【００２２】
本発明による請求項１４の計算処理方法は、請求項６または１３に記載される計算処理方法において、変動幅因子Ｎを２とすることを特徴とする。
【００２３】
本発明による請求項１５の計算処理方法は、請求項１乃至１４のいずれかに記載される計算処理手段において、４倍精度またはそれ以上の精度によるコンパイラーを用いることを特徴とする。
【００２４】
本発明による請求項１６の計算処理方法は、請求項１２に記載される計算処理手段において、判断因子である電流値として、４倍精度コンパイル計算については判断電流密度を１平方センチメータ当り１０のマイナス２０乗アンペア程度とすることを特徴とする。
【００２５】
本発明による請求項１７の計算処理方法は、請求項５または１２に記載される計算処理方法において、判断因子である電流値として、２倍精度コンパイル計算については判断電流密度を１平方センチメータ当り１マイクロアンペア程度とすることを特徴とする。
【発明の効果】
【００２６】
静電ポテンシャル分布を求めるポアソン方程式とキャリア分布を求める電流連続の式を解く手段と、且つ、方程式群を自己無撞着に数値解法する手段と、再結合率を真性キャリア密度と寿命より算出する手段を具備する計算処理方法において、高耐圧デバイス設計ではｐｎ接合の暗電流を高精度、且つ、安定に解く必要がある。本発明によれば、キャリア寿命、あるいは該値の替わりに、準位密度、捕獲断面積、速度、あるいは真性キャリア密度を印加電圧と計算された電流値に応じて適宜制御することで暗電流算出を高速に、安定に、且つ高精度に実現可能とする。
【図面の簡単な説明】
【００２７】
【図１】本発明の計算処理手順を示す模式図である。
【図２】本発明の実施例１にかかる、印加電圧が小さいときの、ガードリング付きシリコンカーバイトダイオードの場合の計算処理手順を示す模式図である。
【図３】２準位型再結合を説明するバンド構造図である。
【図４】２準位型再結合に基づく再結合率の大きさを示す図である。
【図５】本発明の実施例３にかかる、再結合率の自動制御のフローチャートである。
【発明を実施するための形態】
【００２８】
ワイドギャップ半導体中のｐｎ接合が逆バイアスされた状況下におけるキャリアの生成率は真性キャリア密度に比例し、キャリア寿命に反比例することから、真性キャリア密度とキャリア寿命を計算過程で適宜制御することで暗電流の計算精度向上と計算処理速度向上の両立を達成可能とする。
【００２９】
上記の暗電流を計算するための、半導体装置内部のキャリア輸送を解析するシミュレーションを行うコンピュータ装置が実行するシミュレーションの計算処理方法を、以下に述べる。
【００３０】
具体的には、静電ポテンシャル分布を求めるポアソン方程式とキャリア分布を求める電流連続の式を解く手段と、且つ、方程式群を自己無撞着に数値解法する手段と、キャリアの生成・再結合率を真性キャリア密度とキャリア寿命より算出する手段を具備する計算処理方法において、印加電圧が低く降伏電圧に達していない場合、キャリア寿命を本来値より短く設定することでデバイス内部の電流密度を大きく見積もらせるよう計算を実行させる手段が挙げられる。キャリア寿命は準位密度、捕獲断面積、速度の積の逆数であるから該３つの物理定数の内の任意の一つ、または２つの積、または３つ全ての積を本来値より大きく設定することで電流密度と端子電流値（暗電流）を大きく見積もらせるよう計算を実行させる手段も挙げられる。ここで、本来値とは、キャリア寿命を短く設定する前の、または準位密度、捕獲断面積、速度を長く設定する前の本来の値を指す。上記の手段は、印加電圧が低い領域、具体的には印加電圧が降伏電界の三分の二以下の場合に有効である。高電圧領域、具体的には印加電圧が降伏電界の三分の一以上の場合では、上述の手段に加えて、２つの捕獲準位間隔に素電荷量を乗算し、乗算結果をボルツマン定数と絶対温度の積で除算し、除算結果の指数関数を算出し、真性キャリア密度に該除算結果の指数関数を電界強度に依存して乗算する手段が有効である。
【００３１】
まず、印加電圧を零とした時のポアソン方程式と電流連続の式を自己無撞着に解く。最初に、ポアソン方程式を解き、静電ポテンシャル分布を求める。このときのポアソン方程式は以下のとおりである。
【００３２】
【数１】
2014107277.tif 000003
【００３３】
2014107277.tif 000004
【００３４】
次に、電流連続の式を解き、キャリア密度ｎおよびｐを求める。暗電流算出時の電流連続の式は以下のとおりである。
【００３５】
【数２】
2014107277.tif 000005
【００３６】
【数３】
2014107277.tif 000006
【００３７】
ここで、Ｊは電流密度（添え字のｎとｐは電子と正孔の属性）、ｎとｐは電子濃度と正孔濃度、ｑは単位電荷、Ｒが再結合率である。逆バイアスではＲの値は負となり、生成項を表わす。これが暗電流の源となる。定常状態を対象に考えると時間微分項はゼロとなるから、  
【００３８】
【数４】
2014107277.tif 000007
【００３９】
【数５】
2014107277.tif 000008
【００４０】
を得る。真性キャリア密度をｎi、キャリア寿命をτとすると、逆バイアス動作では、
【００４１】
【数６】
2014107277.tif 000009
【００４２】
となる。端子電流（Ｉd）は全電流密度（Ｊ）
J =Jn+Jp                       （７）
を積分して求めることができる。
【００４３】
【数７】
2014107277.tif 000010
【００４４】
これにガウスの定理を用いると
【００４５】
【数８】
2014107277.tif 000011
【００４６】
となる。電流成分である電子電流と正孔電流は、ｐｎダイオードのｎ型領域、ｐ型領域で非零の値を持ち、逆領域ではほぼ零となる。従って、式（８）の積分はｎ型領域では電子電流が、ｐ型領域では正孔電流が非零となる。しかも、ｎ型領域におけるダイオード電流と、ｐ型領域におけるダイオード電流とは互いに等しい （電流不滅の法則）から、式（４）と式（５）いずれに積分を当てはめても同一の解を得ることになる。すなわち、端子電流は式（９）と式（４）、（５）から
【００４７】
【数９】
2014107277.tif 000012
【００４８】
となる。以下、端子電流Ｉdを暗電流とし、暗電流は、式（１０）より算出され、デバイス内部を流れる局所電流もまた、式（１０）の積分範囲を局所領域に限定することで評価できるものとする。
【００４９】
電流連続の式の中の電流密度ベクトルは
【００５０】
【数１０】
2014107277.tif 000013
【００５１】
【数１１】
2014107277.tif 000014
【００５２】
と書くことができ、電流連続の式を解くことで電子密度、正孔密度を求めることが可能となる。ここで、ｖn、ｖpはそれぞれ電子、正孔の移動速度、Ｄn、Ｄpはそれぞれ電子と正孔の拡散係数である。
【００５３】
このとき得られたキャリア密度を用いて、再度ポアソン方程式を解き、電位分布を求める。この繰り返しによって、電位分布とキャリア密度が自己無撞着に収束する。このときの解を初期値とする。電圧を印加するとき、キャリア寿命τを本来値より短く設定することで、電流連続の式を安定に且つ短時間で解くことを可能とする。内部電流の計算精度は向上し、電位が浮遊状態にあるガードリングが設定されている問題でも内部電流により浮遊電位が安定的に決定可能となる。
【００５４】
キャリア寿命を本来値より短く設定する替わりに、準位密度、捕獲断面積、速度の内のいずれか一つまたは２つの積、または全ての積を本来値より大きく設定する手段も有効である。
【００５５】
印加電圧を除々に上げて行くと暗電流は僅かながら増大し、電流連続の式の数値解析は更に安定化する。
【００５６】
キャリア寿命τを本来値より短く設定して電流連続の式を解いた結果、電位分布とキャリア密度が収束し、得られた暗電流がキャリア寿命を本来値としたときよりも増大したとき、キャリア寿命を本来値へ近づくように大きい値へ設定し直すことが可能となる。寿命の替わりに、準位密度、捕獲断面積、速度のいずれかのパラメータを用いた場合は、該用いたパラメータを本来値へ近づくように小さい値へ設定し直すことが可能となる。
【００５７】
印加電圧の増加に伴い増加する電流値を検知しながら、キャリア寿命を小さく再設定することを繰り返すことで本来の寿命の値での解析を実行可能とする。寿命の替わりに、準位密度、捕獲断面積、速度のいずれかのパラメータを用いた場合は、該用いたパラメータあるいは該用いたパラメータの任意の組合せの積を小さい方向へ再設定することを繰り返すことで該用いたパラメータを所定の値へ戻すことが可能となる。
【００５８】
上述の方法は印加電圧が低い領域、具体的には印加電圧が降伏電界の三分の二以下であり、通常、印加電界が１．５ＭＶ／ｃｍ以下の場合に有効である。印加電圧が増加し電界強度がある値を超えたとき、ワイドギャップ半導体では２つの捕獲準位間隔に素電荷量を乗算し、乗算結果をボルツマン定数と絶対温度の積で除算し、除算結果の指数関数を算出し、真性キャリア密度に該因子を電界強度に依存して乗算する手段を有する再結合率を電流連続式に付け加えることが有効である。ここで、上述のある値は、通常、降伏電界の三分の一であり、印加電界が１ＭＶ／ｃｍ以上の値をとる。電界強度と捕獲センターの空間距離と素電荷量を乗算し、捕獲準位間隔から該乗算した値を減算した値を電界強度依存の判断手段として用いることは有効である。判断手段として、ある領域との境界を連続的に表現可能なガウス関数を用いることは数値計算上、有利である。
【００５９】
キャリアの寿命を本来値へ戻す計算処理過程において、電流の計算精度が高くとれないので、ニュートンラプソンのような微分係数を用いる方法より、キャリア寿命をＮ倍化、Ｎ倍化….と繰り返し、電流値を概略Ｎ分の１ずつ引き下げる「挟み込み法」が有効である。
【００６０】
キャリアの寿命を本来値へ戻す計算処理過程において、計算された電流値を判断因子として用いることは有効である。
【００６１】
キャリアの寿命を本来値へ戻す計算処理過程において計算された電流値を判断因子として用いる場合、２倍精度コンパイル計算については判断電流密度を１平方センチメータ当り１マイクロアンペア程度、４倍精度コンパイル計算については判断電流密度を１平方センチメータ当り１０のマイナス２０乗アンペア程度とすることが有効である。
【００６２】
図１は、本発明の原理を模式的に示すもので、縦軸に電流、横軸に印加電圧を例に取り計算処理手順を示すものである。図１における縦軸の電流は暗電流である。キャリア寿命を印加電圧に応じて制御し、ガードリング付きシリコンダイオード構造の電流−電圧特性を計算している。所定のキャリア寿命が１マイクロ秒であることに対し、印加電圧が低い領域では寿命を１ナノ秒まで小さく設定することで暗電流を高精度に且つ安定に解析可能としている。印加電圧の増加に伴い電流値が増加するので、電流がある程度増加したところで、寿命を本来値へ向け、より大きい値へ設定し直しても暗電流を安定に且つ高精度で計算可能としている。電流値の増加を確認しながら寿命を本来値まで引き戻す計算が可能である。これに対し、所定の寿命が１マイクロ秒の状態で低い印加電圧の解析を実行すると２倍精度計算では暗電流の計算精度が不足し、安定な収束計算を実行できない。このことについては非特許文献３で既に述べられている。
【００６３】
本特許の有効性は高耐圧設計のための逆バイアス状態で起こる微小暗電流算出を高精度で且つ安定、高速に演算処理することである。
【００６４】
なお、印加電圧が小さく、電界強度が降伏電界の三分の二以下であるときは、ポアソン方程式と電流連続の式を解き、キャリアの生成・再結合率を算出する際、キャリア寿命を本来値より短く設定することでデバイス内部の電流密度を大きく見積もらせるよう計算を実行させる方法を、印加電圧が増加し、電界強度が降伏電界の三分の一以上であるときは上述の方法に加え、２つの捕獲準位間隔に素電荷量を乗算し、乗算結果をボルツマン定数と絶対温度の積で除算し、除算結果の指数関数を算出し、真性キャリア密度に該除算結果の指数関数を電界強度に依存して乗算する方法を用いることが出来るため、電界強度が降伏電界の三分の一以上三分の二以下であるときは、上述の２つの方法のいずれも用いることが出来る。ここで、前者の方法が、印加電圧が小さく、電界強度が降伏電界の三分の二以下であるとき、後者の方法が、印加電圧が増加し、電界強度が降伏電界の三分の一以上であるときに有効であることは、この印加電圧の範囲において電位分布とキャリア密度が収束したことから得られた。また、図４に示される曲線は、電界強度がある程度大きい領域で成り立つ近似であって、弱電界では、図４に示される曲線そのものが成立しない。
【００６５】
〔実施例１〕
本特許の実施例１では、印加電圧が小さく、電界強度が降伏電界の三分の二以下であるときの暗電流を計算する。本実施例を、図２を用いて説明する。図２はワイドギャップ半導体の典型的例として知られているシリコン・カーバイド（ＳｉＣ）ダイオードを用いた計算である。図２における縦軸の電流は暗電流である。ガードリングと呼ばれる電位の浮遊領域が存在する。この為キャリア寿命が本来値付近にあると、暗電流が極度に小さくなり、ＩＲドロップによる電位計算が不能となる。これに対し、本発明によれば、キャリア寿命を本来値より小さい値に設定することで暗電流を安定に計算可能とさせ、結果、浮遊電位を安定に算出可能としている。図２では印加電圧が６００Ｖと高い値に到達したところで、キャリア寿命を１０のマイナス４８乗から１０のマイナス３８乗まで１０桁に渡り、本来値へむけ変更可能としている。
【００６６】
キャリア寿命は捕獲準位密度、捕獲断面積、速度の積の逆数であるから、キャリア寿命を制御因子とする替わりに、捕獲準位密度、捕獲断面積、速度の内のいずれか一つ、または任意の２つ、または全ての積を制御因子として用いることにより極微小暗電流計算の安定計算を得ることが可能であることは明らかである。
【００６７】
〔実施例２〕
本特許の実施例２では、印加電圧が増加し、電界強度が降伏電界の三分の一以上であるときの暗電流を計算する。印加電圧が増加し、電界強度が降伏電界の三分の一を超えたとき、印加電圧が小さく電界強度が降伏電界の三分の二以下であるときの暗電流の計算方法に加えて、捕獲準位が複数個あるとしてキャリア生成率を計算する方法が、暗電流の計算の精度向上に有効である。本実施例では、捕獲準位の数を２としている。本実施例を、図３を用いて説明する。図３はバンドギャップが３ｅＶを超えるようなワイドギャップ半導体の捕獲準位を模式的に示したものである。一般的に捕獲準位は複数個あり、図３（ａ）では捕獲準位が２つの場合を、図３（ｂ）では捕獲準位が１つの場合を示している。図３（ａ）に示すように捕獲準位が２つの場合、捕獲準位間の遷移による生成・再結合過程が存在するため、１準位型再結合ではキャリアの生成・再結合率が
【００６８】
【数１２】
2014107277.tif 000015
【００６９】
ここに
【００７０】
【数１３】
2014107277.tif 000016
【００７１】
【数１４】
2014107277.tif 000017
【００７２】
であることに対し、２準位型再結合ではキャリアの生成・再結合率が
【００７３】
【数１５】
2014107277.tif 000018
【００７４】
となる生成・再結合過程が存在する。ここで、ＥTが捕獲準位、Ｅiは真性フェルミ準位で、Ｅiはエネルギー禁止帯のほぼ中央にある。また、Ｎtはトラップ密度、σは捕獲断面積、ｖは速度であり、下付きのｅとｐはそれぞれ電子と正孔を表し、ｎは電子密度、ｐは正孔密度、ｎiは真性キャリア密度、ｎ1はトラップ準位１における実効的電子の準位密度、ｐ1はトラップ準位１における実効的正孔の準位密度、Ｎt1、Ｎt2は２つの準位における密度、Δは２つの準位のエネルギー差、βはボルツマン電圧の逆数、ωは２つの準位の間の遷移レート係数、ｎ2はトラップ準位２における実効的電子の準位密度、ｐ2はトラップ準位２における実効的正孔の準位密度であり、Ｘ1Ｘ2はそれぞれ、
【００７５】
【数１６】
2014107277.tif 000019
【００７６】
である。
【００７７】
ｐｎ接合が逆バイアスの場合、ｎｐ＜＜ｎi2であるから、ｎ＝ｐ＝０として、キャリアの生成・再結合率は式（１３）に示されるように１準位過程では真性キャリア密度（ｎi）に比例する。これに対し、２準位過程におけるキャリアの生成・再結合率は真性キャリア密度の２乗に、さらに、準位間エネルギー差Δをボルツマン定数と絶対温度の積であるボルツマン電圧（βの逆数）で除算した値（Δβ）の指数関数分だけ大きな値となる。これにより暗電流の計算が安定化することは明白であり、このような手段を高電界領域で追加することで暗電流を高精度に且つ安定に解析可能とする。
【００７８】
図４は式（１６）に基づくキャリアの生成・再結合率Ｒ2-levelsを準位間の遷移確率ωを横軸にして示したものである。遷移確率ωが低い領域では真性キャリア密度（＝ｎi）がキャリアの生成・再結合率の決定要因であったことに対し、遷移確率ωの上昇と共に、キャリアの生成・再結合率の決定要因はトラップ準位における実効的電子の準位密度（＝ｎt）が支配因子となることを示している。つまり、Ｒ2-levelsの飽和値が（ｎt／２τ）となっており、ｎtの値が１準位過程の場合のＲSRHの分子にあるｎiと比べ圧倒的に大きい。仮に、Ｒ2-levelsが飽和領域に至らずとも、Ｒ2-levelsの値そのものが１準位過程の場合のＲSRHに比べはるかに大きい値となる。トラップ準位が２つある時のトラップ準位における実効的電子の準位密度は真性キャリア密度より大きな値であるから、２準位過程におけるキャリア生成を基にした計算手順が極微小電流（暗電流）の算出に有効であることは明らかである。
【００７９】
式（１４）、（１５）よりＲSRHの分母が最少となるのは捕獲準位が真性フェルミ準位と一致する位置、すなわち禁止帯のほぼ中央に位置する時となるので、捕獲準位を禁止帯のほぼ中央に設定して計算を開始する手段も有効であることは明らかである。
【００８０】
〔実施例３〕
本特許の第３の実施例を、図５を用いて説明する。図５は本特許の原理を説明する図(図１)の手順を自動化する方策を明示したものである。計算の出発点として、キャリア寿命τを本来値よりも大きく設定する（ｓ１０１）。このときのキャリア寿命の初期値をτinitとする。これにより、低印加電圧領域での極微小電流（暗電流）計算の数値精度を向上させることが安定可能となる。具体的にはキャリア寿命の本来値が１マイクロ秒であるとき、印加電圧が低い領域では、寿命の本来値より短い１ナノ秒を設定することで、暗電流を安定に且つ高精度で計算可能としている。
【００８１】
次に、印加電圧Ｖを上昇させ、印加電圧をＶ＋ΔＶとする（ｓ１０２）。続いて、暗電流Ｉを計算する（ｓ１０３）。暗電流Ｉの計算過程においては、数値計算可能な電流レベル（判断電流密度、Ｉ0）より計算された電流値が小さい時（ｓ１０３のＮＯ）は、設定寿命値をそのまま据え置く（ｓ１０２に戻る）。逆に、計算された電流値がＩ0を超えた場合（ｓ１０３のＹＥＳ）、キャリア寿命τを２倍に設定する（ｓ１０５）。次に、２倍化したキャリア寿命がキャリア寿命の本来値τ0を超えるか否かを判断し（ｓ１０６）、超える場合は本来値に設定する（ｓ１０７）。このようにしてキャリア寿命を初期値より本来値へ向け増加させながらバイアスを増加させていくことにより、安定に極微小リーク電流を計算することができる。
【００８２】
判断電流密度（Ｉ0）の大きさとしては、２倍精度コンパイル計算については１平方センチメータ当り１マイクロアンペア程度、４倍精度コンパイル計算については１平方センチメータ当り１０のマイナス２０乗アンペア程度とすることが有効である。
【００８３】
図５では、キャリア寿命を制御因子として説明したが、再結合率を決める因子である捕獲準位密度、捕獲断面積、キャリア速度もまた、制御因子としてもちいることが出来る。この場合は、キャリア寿命と逆数関係にあるので、制御因子の増大・減少を全て逆転させて考えれば良いことは明らかである。
【００８４】
本発明による手法は従来のシリコンデバイスのみならず、実施例であげたシリコンカーバイト材料や窒化ガリウム、酸化ガリウム砒素、ダイアモンド系デバイス解析にも適用可能である。
【図１】
2014107277.tif 000020
【図２】
2014107277.tif 000021
【図３】
2014107277.tif 000022
【図４】
2014107277.tif 000023
【図５】
2014107277.tif 000024