JP5235124B2 - Semiconductor device simulation program - Google Patents

Semiconductor device simulation program Download PDF

Info

Publication number
JP5235124B2
JP5235124B2 JP2008282458A JP2008282458A JP5235124B2 JP 5235124 B2 JP5235124 B2 JP 5235124B2 JP 2008282458 A JP2008282458 A JP 2008282458A JP 2008282458 A JP2008282458 A JP 2008282458A JP 5235124 B2 JP5235124 B2 JP 5235124B2
Authority
JP
Japan
Prior art keywords
energy
carrier
electric field
relaxation coefficient
momentum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2008282458A
Other languages
Japanese (ja)
Other versions
JP2010092445A (en
Inventor
憲 山口
一隆 冨澤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
MEIJI UNIVERSITY LEGAL PERSON
Original Assignee
MEIJI UNIVERSITY LEGAL PERSON
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by MEIJI UNIVERSITY LEGAL PERSON filed Critical MEIJI UNIVERSITY LEGAL PERSON
Priority to JP2008282458A priority Critical patent/JP5235124B2/en
Publication of JP2010092445A publication Critical patent/JP2010092445A/en
Application granted granted Critical
Publication of JP5235124B2 publication Critical patent/JP5235124B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Insulated Gate Type Field-Effect Transistor (AREA)

Description

本発明は半導体装置内部の電子又は正孔の少なくとも一方、あるいは両方の輸送問題を解析することに関するもので、半導体装置の設計、最適化、不良解析等への応用が可能なものである。  The present invention relates to analyzing the transport problem of at least one or both of electrons and holes in a semiconductor device, and can be applied to design, optimization, failure analysis, and the like of a semiconductor device.

半導体装置の微細化は留まるところを知らず、物理的ゲート長は数十ナノメートルの領域に達している。国際的標準化団体による将来予測であるInternational Technology Roadmap for Semiconductors(ITRS)によれば、2015年にはゲート長が10ナノメートルに達し、それ以降、数ナノメートルに達すると予測している。このような極微細半導体装置内部の電子または正孔の輸送は、従来のドリフト・拡散型に対し弾道輸送(Ballistic transport)が顕著になってくる。弾道輸送が顕著になると、電流値や増幅率にプラスの効果が期待される。ITRSによる将来予測でも弾道輸送は大きな技術的ブレイクスルーとして期待されている。  The miniaturization of semiconductor devices does not stop, and the physical gate length has reached an area of several tens of nanometers. According to International Technology Roadmap for Semiconductors (ITRS), which is a future prediction by an international standardization body, the gate length will reach 10 nanometers in 2015, and will reach several nanometers thereafter. As for the transport of electrons or holes inside such an ultrafine semiconductor device, ballistic transport becomes more conspicuous than the conventional drift / diffusion type. If ballistic transport becomes significant, positive effects are expected on the current value and amplification factor. Ballistic transport is expected as a major technological breakthrough in future predictions by ITRS.

弾道輸送の解析と弾道輸送を活用した半導体デバイス設計手段として、これまでモンテカルロ法とボルツマン方程式に基づくエネルギー・運動量バランス方程式を解く方法の2つが知られている。モンテカルロ法はキャリアである電子または正孔を荷電粒子として扱い、荷電粒子1個1個の運動をニュートン力学により追跡するもので、散乱過程を乱数で記述する方法である。電子または正孔のミクロな運動に基づく理論体系をもつので、詳細な物理情報を得ることができる。しかし、膨大な計算時間を必要とする。このため、実際の構造設計には適用できず、基本的物理量を算出する利用に留っている。これに対し、エネルギー・運動量バランス方程式を解く方法は実用的計算時間に収めることができるが、基礎となる緩和係数をモンテカルロ法より求める必要がある。  As a semiconductor device design means utilizing ballistic transport analysis and ballistic transport, two methods have been known so far: a Monte Carlo method and a method for solving an energy-momentum balance equation based on the Boltzmann equation. The Monte Carlo method treats electrons or holes as carriers as charged particles, tracks the movement of each charged particle by Newtonian mechanics, and describes the scattering process with random numbers. Since it has a theoretical system based on microscopic movement of electrons or holes, detailed physical information can be obtained. However, enormous calculation time is required. For this reason, it cannot be applied to an actual structural design, and is only used for calculating a basic physical quantity. On the other hand, the method of solving the energy-momentum balance equation can be put into practical calculation time, but it is necessary to obtain the basic relaxation coefficient by the Monte Carlo method.

緩和係数を支配する因子としてバルクの半導体内部ではキャリアのエネルギーと基板不純物濃度の二つが重要因子である。更に、現代の最新型Large−Scale−Integration(LSI)の基本要素である金属−酸化物−半導体−電界効果型トランジスタ(Metal−Oxide−Semiconductor Field−Effect−Transistor:MOSFET)内部では、キャリアを輸送するソース・ドレイン方向電界と垂直な電界成分、所謂ゲート電界が電子または正孔の分布状態と散乱・緩和過程を決定する上で重要な因子となる。従って、MOSFETでは、バルク半導体における緩和係数決定因子であるエネルギーと基板不純物濃度にゲート電界を加え、合計3つの要素が重要因子として挙げられることになる。3要素の重要性は冨澤一隆著「半導体デバイスシミュレーション」(コロナ社)、並びに”Numerical solution of submicron semiconductor devices”(サブミクロン半導体デバイスの数値解析)Artech House,0−89006−620−5の中で記述されている。  Two factors that control the relaxation coefficient are the carrier energy and the substrate impurity concentration inside the bulk semiconductor. Further, carriers are transported inside metal-oxide-semiconductor field-effect-transistors (MOSFETs), which are the basic elements of the latest modern large-scale-integration (LSI). The electric field component perpendicular to the source / drain direction electric field, the so-called gate electric field, is an important factor in determining the electron or hole distribution and the scattering / relaxation process. Therefore, in the MOSFET, the gate electric field is added to the energy and substrate impurity concentration, which are determinants of the relaxation coefficient in the bulk semiconductor, and a total of three factors are listed as important factors. The importance of the three elements is “Semiconductor Device Simulation” written by Kazutaka Serizawa (Corona) and “Numerical solution of submicron semiconductor devices” (Numerical Analysis of Submicron Semiconductor Devices) Arttech House, 0-89006-620-5. It is described in.

エネルギー・運動量バランス方程式を具体的に解くには、緩和係数をエネルギー、不純物濃度、ゲート電界の3変数の関数として準備する必要がある。しかし、緩和係数を3変数に対し決定するには、しかも、計算精度向上のためには、各変数の変域を広くかつ、細かな分割で求める必要がある。これに対応するモンテカルロ法計算を実行することは、最新の計算機環境を用いても天文学的時間となり、事実上、実行困難である。
冨澤一隆著「半導体デバイスシミュレーション」(コロナ社 1996年) Kazutaka Tomizawa”Numerical solution of submicron semiconductor devices”(サブミクロン半導体デバイスの数値解析)Artech House,0−89006−620−5 K.Yamaguchi“A mobility model for carriers in the MOS inversion layer”,(MOS反転層内のキャリアに関する移動度モデル)IEEE Trans.Electron Devices,vol.ED−30,no.6,pp.658−663,June 1983
To specifically solve the energy / momentum balance equation, it is necessary to prepare a relaxation coefficient as a function of three variables of energy, impurity concentration, and gate electric field. However, in order to determine the relaxation coefficient for the three variables, and to improve the calculation accuracy, it is necessary to obtain the domain of each variable in a wide and fine division. Executing the Monte Carlo calculation corresponding to this becomes astronomical time even using the latest computer environment, and is practically difficult to execute.
"Semiconductor device simulation" by Kazutaka Serizawa (Corona Corp. 1996) Kazutaka Tomizawa “Numerical solution of submicron semiconductor devices” (Numerical analysis of submicron semiconductor devices) Arttech House, 0-89006-620-5 K. Yamaguchi “A mobility model for carriers in the MOS inversion layer”, (mobility model for carriers in MOS inversion layer) IEEE Trans. Electron Devices, vol. ED-30, no. 6, pp. 658-663, June 1983

本発明は上記従来の問題点を解決すべく各変数間の相関関係を物理的に精査・検討し、その物理的性質を明らかにすることにより、合理的で精度が高く、且つ、適用範囲の広い緩和係数を短時間で提供することを課題とする。  In order to solve the above-mentioned conventional problems, the present invention physically examines and examines the correlation between each variable, and clarifies the physical properties thereof. It is an object to provide a wide relaxation coefficient in a short time.

課題は、電子又は正孔が強いゲート電界を感じながらソースからドレインへ弾道輸送されるときに受ける散乱機構の成分の内、従属関係にある因子と比較的独立な関係にある因子を識別し、各因子の依存性に成分分離することで、精度向上と計算処理速度向上の両立を達成可能とする。  The problem is to identify factors that are relatively independent of the dependent factors among the components of the scattering mechanism that electrons or holes undergo when ballistically transported from the source to the drain while feeling a strong gate electric field, By separating the components into the dependency of each factor, it is possible to achieve both improvement in accuracy and calculation speed.

具体的には、静電ポテンシャル分布を求めるポアソン方程式とキャリア分布を求める電流連続の式を解く手段と、運動量・エネルギー状態を求める運動量・エネルギーバランス方程式を解く手段と、且つ、方程式群を自己無撞着に数値解法する手段と、運動量・エネルギーバランス方程式を解くための運動量・エネルギー緩和係数の数値データをエネルギー、不純物濃度、ゲート電界の3変数の関数として具備する計算処理方法において、3変数表現による運動量・エネルギー緩和係数の数値データ群を2変数ずつの組合せへ分解する計算処理方法が挙げられる。2変数データを更に1変数データの組合せへ分解することも挙げられる。  Specifically, there are means for solving the Poisson equation for obtaining the electrostatic potential distribution and current continuity for obtaining the carrier distribution, means for solving the momentum / energy balance equation for obtaining the momentum / energy state, In a calculation processing method comprising a numerical solution method for coherence and numerical data of momentum / energy relaxation coefficient for solving the momentum / energy balance equation as a function of three variables of energy, impurity concentration, and gate electric field. There is a calculation processing method that decomposes a numerical data group of momentum and energy relaxation coefficient into combinations of two variables. It is also possible to further decompose the two-variable data into a combination of one-variable data.

まず、MOSFET内部を運動するキャリアの散乱確率は上述したようにエネルギー、不純物濃度、ゲート電界強度の3変数により表現される。  First, as described above, the scattering probability of carriers moving inside the MOSFET is expressed by the three variables of energy, impurity concentration, and gate electric field strength.

f=f(w,N,E) (1)
ここに、wがエネルギー、Nが不純物濃度、Eがゲート電界強度を表し、fが散乱確率を表している。前記該3変数の成分分離は、理論的には、散乱プロセスの独立性により判断できる。独立な散乱過程が複数種類存在する時、全体の散乱確率は各散乱確率の和となることが一般的に知られている。該MOSFET内部のキャリアの散乱過程を考えてみると、反転層内におけるキャリアの散乱過程とバルク中でのキャリアの散乱過程は第1近似として独立事象と言える。そこで、エネルギー(w)、不純物濃度(N)、ゲート電界(E)の3変数により表現される全散乱確率は、エネルギーとバルク中の不純物濃度により決まる散乱確率[fBULK(w,N)]と、反転層内で起こるキャリア散乱の確率fInversion(w,E)の和として
f = f (w, N, E G ) (1)
Here, w denotes energy, N is the impurity concentration, E G is a gate electric field strength, f represents the scattering probability. The component separation of the three variables can theoretically be determined by the independence of the scattering process. It is generally known that when there are a plurality of independent scattering processes, the total scattering probability is the sum of the scattering probabilities. Considering the carrier scattering process in the MOSFET, it can be said that the carrier scattering process in the inversion layer and the carrier scattering process in the bulk are independent events as a first approximation. Therefore, the total scattering probability expressed by the three variables of energy (w), impurity concentration (N), and gate electric field (E G ) is a scattering probability [f BULK (w, N) determined by the energy and the impurity concentration in the bulk. ] And the probability of carrier scattering occurring in the inversion layer f Inversion (w, E G )

f(w,N,E)=fBULK(w,N)+fInversion(w,E) (2)
の2成分分離が可能となる。更に詳しく式(2)の意味するところを考えてみる。簡単のため、特定のエネルギー状態(w)を仮定して不純物濃度依存性とゲート電界依存性を考えてみると、不純物濃度Nとゲート電界Eの2変数で決まるような散乱確率は、該散乱過程が独立事象であれば、不純物濃度Nにより決定される散乱確率とゲート電界Eにより決定される散乱確率の和といえる。このことは、不純物濃度Nにより決定される散乱確率、実質的に1次元データと、また同様に1次元データであるゲート電界Eにより決定される散乱確率を知れば、2次元データを簡易に、精度良く推定できることを意味している。
f (w, N, E G ) = f BULK (w, N) + f Inversion (w, E G ) (2)
It is possible to separate the two components. Consider in more detail the meaning of equation (2). For simplicity, when consider the hypothetical impurity concentration dependency and gate field dependence of the specific energy state (w), the scattering probability as determined by two variables of the impurity concentration N and the gate electric field E G, said if scattering process is independent event, it can be said a sum of the scattering probability determined by the scattering probability and gate field E G determined by the impurity concentration N. This scattering probability determined by the impurity concentration N, substantially one-dimensional data, also knowing the scattering probability determined by the gate electric field E G are equally 1-dimensional data, two-dimensional data easily This means that it can be estimated with high accuracy.

キャリアの運動の散乱・緩和過程として、エネルギー緩和と運動量緩和の2種類が存在する。そこで、式(2)で記載される散乱確率をエネルギー緩和と運動量緩和のそれぞれの計算処理手段に分けて考えて行く。先ず第1にエネルギー緩和計算手段について考えてみる。不純物による散乱はキャリアの運動の方向を変化させるだけの弾性散乱であるため、エネルギーの授受が発生しない。従って、エネルギー緩和過程において、不純物濃度依存性は消失する。つまり、式(2)の右辺第1項から不純物濃度依存性が消失する。エネルギー緩和に関しては更に簡単化され、次のように表記可能となる。  There are two types of carrier motion scattering / relaxation processes: energy relaxation and momentum relaxation. Therefore, the scattering probability described by equation (2) will be considered separately for each of the energy relaxation and momentum relaxation calculation processing means. First, let us consider energy relaxation calculation means. Since scattering due to impurities is elastic scattering that only changes the direction of carrier motion, no energy is transferred. Therefore, the impurity concentration dependency disappears in the energy relaxation process. That is, the impurity concentration dependency disappears from the first term on the right side of the equation (2). Energy relaxation is further simplified and can be expressed as follows.

energy(w,N,E)=fenergy_BULK(w)+fenergy_Inversion(w,E) (3)
式(3)の右辺第1項はエネルギーのみに依存する完全1次元表記である。
f energy (w, N, E G ) = f energy_BULK (w) + f energy_Inversion (w, E G ) (3)
The first term on the right side of Equation (3) is a complete one-dimensional notation that depends only on energy.

次に運動量緩和の計算処理過程を考えてみる。この場合、バルク中の不純物散乱が低電界の移動度を決定する主要因であり、式(2)の形式は保存され  Next, consider the calculation process of momentum relaxation. In this case, impurity scattering in the bulk is the main factor determining the mobility of the low electric field, and the form of equation (2) is preserved.

moment(w,N,E)=fmoment_BULK(w,N)+fmoment_Inversion(w,E) (4)
となる。以上纏めると、エネルギー、不純物濃度、ゲート電界の3変数により標記されていた散乱確率はエネルギーと不純物濃度、エネルギーとゲート電界の2変数ずつの組合せにより表現可能なことがわかる。つまり、3次元データ構造が2次元データ構造2組に分解され、計算処理が簡単化される。
f moment (w, N, E G ) = f moment_BULK (w, N) + f moment_Inversion (w, E G ) (4)
It becomes. In summary, it can be seen that the scattering probability marked by the three variables of energy, impurity concentration, and gate electric field can be expressed by a combination of two variables of energy and impurity concentration and energy and gate electric field. That is, the three-dimensional data structure is decomposed into two sets of two-dimensional data structures, and the calculation process is simplified.

式(3)と式(4)の右辺第2項は共にエネルギーwとゲート電界Eの2変数構造となっている。ゲート電界依存項は基本的には状態密度依存性が主体であるので、あるエネルギーwにおける緩和係数がf(w,E(1))であるなら、f(w,E(2))はf(w,E(1))の定数倍(g倍)で、その定数は状態密度に比例、結果としてEに依存することになる。従って、gという項目はg(E)と書き下すことができる。つまり、任意のゲート電界に対する緩和係数は、或る出発点として選択したゲート電界E(0)における値f(w,E(0))を基に、g倍した値を加算し続けることに等価である。加算を継続することは乗算に他ならずThe second term of the right side of equation (3) and (4) are both a two variable structure of energy w and the gate electric field E G. Since the gate electric field dependence term is basically state density dependent, if the relaxation coefficient at a certain energy w is f (w, E G (1)), f (w, E G (2)). in f constant multiple (g times) of (w, E G (1)), the constant is proportional to the density of states, it will depend on the E G as a result. Therefore, the item g can be written down as g (E G ). That is, as the relaxation coefficient for an arbitrary gate electric field, the value multiplied by g is continuously added based on the value f (w, E G (0)) in the gate electric field E G (0) selected as a certain starting point. Is equivalent to Continuing addition is nothing but multiplication.

Inversion(w,E)=fInversion_energy(w)×g(E) (5)
となる。式(5)の意味するところは、エネルギーとゲート電界の2変数構造はエネルギーとゲート電界の各1次元データ構造の積に分解されることであり、計算処理手順が大幅に軽減される。式(5)がエネルギー緩和項、運動量緩和項のいずれにも適用可能なこと、即ち、式(3)と式(4)の右辺第2項のいずれにも適用可能なことは言うまでもない。
f Inversion (w, E G ) = f Inversion_energy (w) × g (E G ) (5)
It becomes. The meaning of equation (5) is that the two-variable structure of energy and gate electric field is decomposed into the product of each one-dimensional data structure of energy and gate electric field, and the calculation processing procedure is greatly reduced. Needless to say, the expression (5) can be applied to both the energy relaxation term and the momentum relaxation term, that is, the expression (3) and the second term on the right side of the expression (4).

運動量緩和項、つまり、式(4)の右辺第1項は、エネルギーwと不純物濃度Nの2変数構造となっている。半導体中のキャリア散乱は、エネルギー状態が低い場合の運動量緩和は主に不純物散乱が支配し、一方、エネルギー状態が高い場合は光学的フォノン散乱が支配的となり、不純物散乱からの寄与はほとんど無視することができる。つまり、低エネルギー側と高エネルギー側の緩和過程は独立事象ということができる。そこで、或る緩和係数fmoment_BULK(w,N)に対し、不純物濃度がNからNに変化したことによる緩和係数の変化分は、エネルギーwに依存した変化分The momentum relaxation term, that is, the first term on the right side of Equation (4) has a two-variable structure of energy w and impurity concentration N. In carrier scattering in semiconductors, momentum relaxation is mainly dominated by impurity scattering when the energy state is low, whereas optical phonon scattering is dominant when the energy state is high, and the contribution from impurity scattering is almost ignored. be able to. In other words, the relaxation processes on the low energy side and the high energy side can be said to be independent events. Therefore, for a certain relaxation coefficient f moment_BULK (w, N 1 ), the change in the relaxation coefficient due to the change in the impurity concentration from N 1 to N 2 is the change depending on the energy w.

Δfmoment_BULK(w) (6)
に定数倍した重みをつけて加算することに相当する。式(6)における変化分を全てのエネルギーについて加算操作を継続することは、結果として、乗算になり
Δf moment_BULK (w) (6)
It is equivalent to adding a weight multiplied by a constant to. Continuing the addition operation for all energies with the change in equation (6) results in multiplication.

Δfmoment_BULK(w,N)=h×Δfmoment_BULK_energy(w) (7)
となる。式(7)における重み係数hは不純物濃度Nに依存して決定されるので、h(N)と記述できる。従って、式(7)は
Δf moment_BULK (w, N) = h × Δf moment_BULK_energy (w) (7)
It becomes. Since the weighting factor h in Equation (7) is determined depending on the impurity concentration N, it can be described as h (N). Therefore, equation (7) is

Δfmoment_BULK(w,N)=h(N)×Δfmoment_BULK_energy(w) (8)
と拡張される。
Δf moment_BULK (w, N) = h (N) × Δf moment_BULK_energy (w) (8)
And expanded.

ここで、或る出発点となる不純物濃度Nにおける緩和係数をHere, the relaxation coefficient at the impurity concentration N 0 as a certain starting point is

moment_BULK(w,N) (9)
と記述しておくことにすると、任意の不純物濃度に対する運動量緩和係数は
f moment_BULK (w, N 0 ) (9)
The momentum relaxation coefficient for any impurity concentration is

moment_BULU(w,N)=fmoment_BULU(w,N)+h(N)×Δfmoment_BULK_energy(w) (10)
となる。式(10)の右辺第1項は特定の不純物濃度に対してのエネルギー依存性を記載したものであるから、式(10)全体をみて明らかな如く、エネルギーと不純物濃度依存性の2次元データ構造もまた、エネルギーと不純物濃度の各1次元データ構造に分解できる。結果、2次元構造は1次元構造数列の単純積と単純和で処理でき、数値処理手順が簡略化できることが分かる。
f moment_BULU (w, N) = f moment_BULU (w, N 0 ) + h (N) × Δf moment_BULK_energy (w) (10)
It becomes. Since the first term on the right side of the equation (10) describes the energy dependency with respect to a specific impurity concentration, as is apparent from the entire equation (10), two-dimensional data on the energy and impurity concentration dependency. The structure can also be broken down into one-dimensional data structures of energy and impurity concentration. As a result, it can be seen that the two-dimensional structure can be processed by a simple product and a simple sum of a one-dimensional structure sequence, and the numerical processing procedure can be simplified.

MOSFETの表面反転層におけるエネルギー、運動量緩和係数の展開式(式(5))において、エネルギーとゲート電界依存性に分解された右辺に式(10)の右辺第1項に相当する項目が無いのは、ゲート電界を小さくした極限では表面反転層が消失するためである。仮に、極端に低いゲート電界を選択した場合、反転層内におけるキャリアの散乱確率は、バルク中の該散乱確率よりも低い値となり、実質上、意味を持たない。従って、適当なゲート電界におけるfInversion(w)を決定した後、ゲート電界が大きくなればg(E)として1より大きい数を、逆に、ゲート電界が小さくなればg(E)として1より小さい数を設定すれば良い。In the expansion formula (formula (5)) of the energy and momentum relaxation coefficient in the surface inversion layer of the MOSFET, there is no item corresponding to the first term on the right side of formula (10) on the right side decomposed into energy and gate electric field dependence. This is because the surface inversion layer disappears in the limit where the gate electric field is reduced. If an extremely low gate electric field is selected, the carrier scattering probability in the inversion layer is lower than the scattering probability in the bulk, which is substantially meaningless. Therefore, after determining f Inversion (w) in an appropriate gate electric field, if the gate electric field is increased, g (E G ) is set to a number larger than 1, and conversely, if the gate electric field is decreased, g (E G ) is set. A number smaller than 1 may be set.

式(5)、式(10)に従って計算された値をデータベース化して保存、利用することは計算処理を高速に行う上で有利である。散乱確率である緩和係数をエネルギー、不純物濃度、ゲート電界の3変数の関数として事前に計算しデータベースとして外部記憶装置に格納、計算時に計算機本体へ移動し利用することによる高速処理も有効である。It is advantageous to store and use the values calculated according to the equations (5) and (10) in a database in order to perform the calculation process at high speed. It is also effective to perform high-speed processing by calculating a relaxation coefficient, which is a scattering probability, in advance as a function of three variables of energy, impurity concentration, and gate electric field, storing it in an external storage device as a database, moving to the computer main body at the time of calculation and using it.

ボルツマン方程式に基づくエネルギー・運動量バランス方程式を数値解析する方法において、エネルギーと運動量の緩和係数をモンテカルロ法数値実験により求める必要がある。MOSFET内の電子または正孔の運動に該方法を適用する場合、緩和係数はエネルギー、不純物濃度、ゲート電界の3つの因子について、それぞれの依存性をモンテカルロ法数値実験により求める必要がある。エネルギー、不純物濃度、ゲート電界の3次元構造のデータを得ることは、最新の計算機環境を用いても事実上不可能な計算時間である。本発明によれば、3つの因子に対し各依存性を記述する1次元データ構造を作成した後、それぞれの依存性同志の乗算処理と加算処理の組合せにより3次元データ構造作成を高速に、且つ高精度に実現可能となる。In the method of numerical analysis of the energy-momentum balance equation based on the Boltzmann equation, it is necessary to obtain the relaxation coefficient of energy and momentum by the Monte Carlo numerical experiment. When this method is applied to the movement of electrons or holes in the MOSFET, it is necessary to determine the dependence of the relaxation coefficient on the three factors of energy, impurity concentration, and gate electric field by Monte Carlo numerical experiments. Obtaining data of the three-dimensional structure of energy, impurity concentration, and gate electric field is a calculation time that is virtually impossible even with the latest computer environment. According to the present invention, after creating a one-dimensional data structure that describes each dependency on three factors, a three-dimensional data structure can be created at high speed by combining multiplication and addition processing of the respective dependencies. It can be realized with high accuracy.

図1は、本発明の原理を模式的に示すもので、縦軸にエネルギー緩和係数、横軸に電子エネルギーを例に取り計算処理手順を示すものである。バルク中の電子のエネルギー緩和係数はゲート電界がゼロの状態で定義され、かつ基板不純物濃度に殆ど依存しないので、エネルギー依存性は図1の中で1本の曲線にて示される。これは、式(3)の右辺第1項に相当する。一方、MOS反転層内における電子の緩和係数はゲート電界に依存するので、複数本描画してある。特定のゲート電界に着目すると、反転層内の擬2次元電子ガスの緩和係数はエネルギーに対し、急激な増加と減少からなり、高エネルギー側でバルク中の電子の緩和係数曲線と交差する。一般的にMOSFET内部の電子は陰極から放出され陽極に到達する間、ゲート電界が弱く、あるいはゼロまたはゼロに近いバルク的伝導形態とMOS反転層内の擬2次元電子ガス状態の2つの伝導形態を経緯する。緩和係数はキャリアの散乱確率そのものであるから、図1における2つの緩和係数、即ちバルク中の電子の緩和係数と反転層内の緩和係数のいずれか大きい方を採用すれば良いことが理解される。このように2つの曲線の外周線を採用する計算処理方法については非特許文献1で既に述べられている。  FIG. 1 schematically shows the principle of the present invention, and shows a calculation processing procedure taking an energy relaxation coefficient on the vertical axis and electron energy on the horizontal axis. Since the energy relaxation coefficient of electrons in the bulk is defined in a state where the gate electric field is zero and hardly depends on the substrate impurity concentration, the energy dependency is shown by a single curve in FIG. This corresponds to the first term on the right side of Equation (3). On the other hand, since the relaxation coefficient of electrons in the MOS inversion layer depends on the gate electric field, a plurality of lines are drawn. Focusing on a specific gate electric field, the relaxation coefficient of the quasi-two-dimensional electron gas in the inversion layer consists of a sudden increase and decrease with respect to energy, and intersects the relaxation coefficient curve of the electrons in the bulk on the high energy side. In general, electrons inside the MOSFET are emitted from the cathode and reach the anode. The gate electric field is weak, or two conduction forms, a bulk conduction form that is zero or close to zero and a quasi-two-dimensional electron gas state in the MOS inversion layer. The background. Since the relaxation coefficient is the carrier scattering probability itself, it is understood that the larger one of the two relaxation coefficients in FIG. 1, that is, the relaxation coefficient of electrons in the bulk or the relaxation coefficient in the inversion layer may be adopted. . Non-Patent Document 1 has already described a calculation processing method that employs the outer peripheral line of two curves in this way.

本特許の有用性はMOS反転層内の緩和係数がエネルギー依存とゲート電界依存の2変数関数、データ構造で言えば2次元データ構造となっていることを、ある特定のゲート電界に対するエネルギー依存性の定数倍とすることで、異なるゲート電界に対する該依存性を表現する計算処理手法、即ち、1次元データ構造で処理する方法である。図1に記載した如く、ある特定のゲート電界に対する緩和係数のエネルギー依存性が分かっていれば、異なるゲート電界における該依存性はゲート電界の増加に対し1より大きい係数を乗算することにより、小さければ1より小さい係数を乗算することで、任意のエネルギー、ゲート電界に対する該依存性を精度良く推定可能とし、2次元データを高速に演算処理することである。  The usefulness of this patent is that the relaxation coefficient in the MOS inversion layer is an energy-dependent and gate-field-dependent two-variable function, that is, a two-dimensional data structure in terms of data structure. This is a calculation processing method for expressing the dependence on different gate electric fields, that is, a method of processing with a one-dimensional data structure. As shown in FIG. 1, if the energy dependence of the relaxation coefficient for a particular gate field is known, the dependence on the different gate fields can be reduced by multiplying the increase in gate field by a factor greater than one. For example, by multiplying a coefficient smaller than 1, the dependence on arbitrary energy and gate electric field can be estimated with high accuracy, and two-dimensional data is processed at high speed.

本特許の実施例を、図2を用いて説明する。図2はモンテカルロ法計算により、異なるゲート電界に対するエネルギー緩和係数のエネルギー依存性を数値実験により求めたものである。ゲート電界として低い方から500kV/cm、1MV/cm、2MV/cm、3MV/cmの4種類について数値実験を行った。モンテカルロ法計算は電子1個1個を多数回追跡してエネルギー状態を求めるため、図2の表現では多数の点の集合となる。本数値実験では滑らかな曲線が得られるようにするため、計算点数を増加させ、最終的に得られた多数の点の集合を滑らかに内挿する曲線で結果を示している。上記4種類のゲート電界の内、基準値として1MV/cmの時のエネルギー依存性を選択してみると、ゲート電界の増加に対して、例えば2MV/cm、3MV/cmに対してそれぞれ1.58倍、2.09倍、逆にゲート電界の減少に対して、例えば500kV/cmに対しては0.62倍することにより、モンテカルロ法計算結果と一致する曲線を得ることができる。これは、ゲート電界をパラメータとして緩和係数のエネルギー依存性を詳細に数値実験する替わりに、特定のゲート電界について緩和係数のエネルギー依存性を1回求め、ゲート電界依存性については定数倍することで精度良く緩和係数をエネルギーとゲート電界の関数として求めることができることを示している。  An example of this patent will be described with reference to FIG. FIG. 2 shows the energy dependence of the energy relaxation coefficient with respect to different gate electric fields obtained by numerical experiments by Monte Carlo calculation. Numerical experiments were performed on four types of gate electric fields of 500 kV / cm, 1 MV / cm, 2 MV / cm, and 3 MV / cm from the lowest. Since the Monte Carlo calculation calculates the energy state by tracking each electron many times, the expression in FIG. 2 is a set of many points. In this numerical experiment, in order to obtain a smooth curve, the number of calculation points is increased, and the result is shown by a curve that smoothly interpolates a set of many finally obtained points. When the energy dependence at the time of 1 MV / cm is selected as the reference value among the above four types of gate electric fields, the increase in the gate electric field is, for example, 1 for each of 2 MV / cm and 3 MV / cm. By reducing the gate electric field by 58 times, 2.09 times, or by 0.62 for 500 kV / cm, for example, a curve that matches the Monte Carlo calculation result can be obtained. Instead of conducting detailed numerical experiments on the energy dependence of the relaxation coefficient using the gate electric field as a parameter, the energy dependence of the relaxation coefficient is obtained once for a specific gate electric field, and the gate electric field dependence is multiplied by a constant. It shows that the relaxation coefficient can be obtained with high accuracy as a function of energy and gate electric field.

エネルギー緩和係数のゲート電界依存性をモンテカルロ法計算により数値実験したのが図3である。式(5)の右辺で定義された係数g(E)を、図3に示した数値実験から求めるとFIG. 3 shows a numerical experiment of the dependence of the energy relaxation coefficient on the gate electric field by Monte Carlo calculation. When the coefficient g (E G ) defined on the right side of Equation (5) is obtained from the numerical experiment shown in FIG.

g(E)=A×(E (11)
と簡単な級数展開が可能となる。式(5)の左辺で記述される如くエネルギー緩和係数は本来、エネルギー、不純物濃度、ゲート電界の3変数関数表示、即ち、3次元データ構造を有しているが、本発明によれば、式(3)の右辺に展開した如く、第1項のエネルギー依存の1次元データ構造と第2項のエネルギーとゲート電界による2次元データ構造に分解、簡略化される。更に、右辺第2項の2次元データ構造も、式(5)に示される如く1次元データ構造2組の単純な乗算で表記可能となる。具体的にはエネルギー依存を表現する1次元データ構造と式(11)で表記されるようなゲート電界依存に対する1次元データ構造との乗算で2次元データを作成可能となる。これにより、モンテカルロ法計算回数の大幅な減少を図り、精度と計算速度の同時向上を図ることが可能となる。
g (E G ) = A × (E G ) B (11)
Simple series expansion is possible. As described on the left side of the equation (5), the energy relaxation coefficient originally has a three-variable function display of energy, impurity concentration, and gate electric field, that is, a three-dimensional data structure. As developed on the right side of (3), the energy-dependent one-dimensional data structure of the first term and the two-dimensional data structure by the energy of the second term and the gate electric field are simplified. Furthermore, the two-dimensional data structure of the second term on the right side can also be expressed by simple multiplication of two sets of one-dimensional data structures as shown in Equation (5). Specifically, two-dimensional data can be created by multiplying a one-dimensional data structure expressing energy dependence with a one-dimensional data structure for gate electric field dependence as expressed by Equation (11). As a result, the number of Monte Carlo calculations can be greatly reduced, and the accuracy and calculation speed can be improved simultaneously.

本特許の第2の実施例を、図4を用いて説明する。図4はモンテカルロ法数値実験により求めた運動量緩和係数のエネルギー依存性を表現したものである。図中、高エネルギー側で急激な増加傾向を示しているのがバルク中の電子の運動量緩和係数で、逆に、高エネルギー側で飽和傾向を示すのがMOS反転層内の擬2次元電子ガスの運動量緩和係数である。バルク中の電子の緩和係数はエネルギーレベルが低い状態で不純物濃度に依存するが、高エネルギー側では不純物濃度に殆ど依存しない。低エネルギー側では不純物散乱(弾性散乱)が優勢であるが、高エネルギー側では光学フォノン散乱が優勢となり不純物による弾性散乱効果が相対的に低減されるためである。従って、運動量緩和係数についても、エネルギー緩和の場合と同様、全緩和係数はバルク中の電子の緩和係数を表現する曲線と、反転層内の該緩和係数曲線の外周線をとればよいことが理解できる。  A second embodiment of this patent will be described with reference to FIG. FIG. 4 represents the energy dependency of the momentum relaxation coefficient obtained by the Monte Carlo numerical experiment. In the figure, the momentum relaxation coefficient of the electrons in the bulk shows a sharp increase trend on the high energy side, and conversely, the pseudo two-dimensional electron gas in the MOS inversion layer shows a saturation trend on the high energy side. Is the momentum relaxation coefficient. The relaxation coefficient of electrons in the bulk depends on the impurity concentration at a low energy level, but hardly depends on the impurity concentration on the high energy side. This is because impurity scattering (elastic scattering) is dominant on the low energy side, but optical phonon scattering is dominant on the high energy side, and the elastic scattering effect due to impurities is relatively reduced. Therefore, as for the momentum relaxation coefficient, as in the case of energy relaxation, it is understood that the total relaxation coefficient should be a curve expressing the relaxation coefficient of electrons in the bulk and the outer circumference of the relaxation coefficient curve in the inversion layer. it can.

反転層内の擬2次元電子ガスの運動量緩和係数はゲート電界に依存する。異なるゲート電界に対する運動量緩和係数のエネルギー依存性を数値実験により求めたものを、図2の場合と同様、ゲート電界として低い方から500kV/cm、1MV/cm、2MV/cm、3MV/cmの4種類について図4に示している。異なるゲート電界に対するデータ列は類似の形状を成し、相似形を形成していることが推測できる。  The momentum relaxation coefficient of the quasi-two-dimensional electron gas in the inversion layer depends on the gate electric field. The energy dependence of the momentum relaxation coefficient with respect to different gate electric fields was obtained by numerical experiments. As in the case of FIG. 2, the gate electric fields from the lowest one of 500 kV / cm, 1 MV / cm, 2 MV / cm, 3 MV / cm The types are shown in FIG. It can be inferred that the data strings for different gate electric fields have similar shapes and form similar shapes.

ゲート電界が1MV/cmにおける運動量緩和係数のエネルギー依存性を基準に設定し、異なるゲート電界に対してはゲート電界に依存した重みづけ係数を乗算することを試みた。結果を図5に示す。反転層内の擬2次元電子ガスのエネルギー依存性は低エネルギー側で発現するので、図5では高エネルギー側を削除している。特定のゲート電界におけるエネルギー依存性を示す曲線を、図中、縦方向へ一定比率で拡大・縮小することで、モンテカルロ法数値実験結果を良く近似できることが分かった。なお、モンテカルロ法数値実験データは多数の点の集合で表現されるため、点データを内挿する曲線を作成したところ、縦方向へ一定比率で拡大・縮小した該曲線と一致している。このため、図では各ゲート電界に対し1本の曲線となっている。従って、式(4)の右辺第2項もまた、式(5)と同一形式で分離表現できることが数値実験により確認される。ここでも、本発明により2次元データ構造が1次元データ構造同志の乗算で表現可能なことが実証されている。  The gate electric field was set based on the energy dependence of the momentum relaxation coefficient at 1 MV / cm, and different gate electric fields were tried to be multiplied by a weighting coefficient depending on the gate electric field. The results are shown in FIG. Since the energy dependence of the quasi-two-dimensional electron gas in the inversion layer appears on the low energy side, the high energy side is deleted in FIG. It was found that the Monte Carlo numerical experiment results can be well approximated by enlarging / reducing the curve indicating the energy dependence in a specific gate electric field at a constant ratio in the vertical direction in the figure. In addition, since the Monte Carlo method numerical experimental data is represented by a set of a large number of points, a curve for interpolating the point data is created, which coincides with the curve expanded and reduced at a constant ratio in the vertical direction. For this reason, in the figure, there is one curve for each gate electric field. Therefore, it is confirmed by numerical experiments that the second term on the right side of Equation (4) can also be expressed separately in the same format as Equation (5). Again, it has been demonstrated by the present invention that a two-dimensional data structure can be represented by multiplication of one-dimensional data structures.

反転層内の擬2次元電子ガスの運動量緩和係数のゲート電界依存性をモンテカルロ法数値実験により求めたものが図6である。図3の場合と同様、ゲート電界依存性は式(11)で示した級数で精度良く近似可能である。  FIG. 6 shows the gate electric field dependence of the momentum relaxation coefficient of the quasi-two-dimensional electron gas in the inversion layer obtained by a Monte Carlo numerical experiment. As in the case of FIG. 3, the gate electric field dependency can be approximated with high accuracy by the series shown in the equation (11).

本特許の第3の実施例として、バルク中の電子の運動量緩和係数について適用した場合を図7に示す。図はモンテカルロ法数値実験結果として、バルクシリコン中の電子の運動量緩和係数のエネルギー依存性について不純物濃度をパラメータとして示したものである。いずれも右上がりで、エネルギーレベルが低い領域で緩和係数に大きな差が出現し、高エネルギー側では不純物濃度に依存しなくなる。低エネルギー側では不純物散乱による弾性散乱が優勢で、高エネルギー側では光学フォノンによる非弾性散乱が優勢となった結果である。  FIG. 7 shows a case where the third embodiment of this patent is applied to the momentum relaxation coefficient of electrons in the bulk. The figure shows, as a result of Monte Carlo numerical experiments, the impurity concentration as a parameter for the energy dependence of the momentum relaxation coefficient of electrons in bulk silicon. In all cases, the difference increases to the right, a large difference appears in the relaxation coefficient in the region where the energy level is low, and the high energy side does not depend on the impurity concentration. The result is that elastic scattering due to impurity scattering is dominant on the low energy side, and inelastic scattering due to optical phonons is dominant on the high energy side.

図7の中で、基板不純物濃度が高い領域(1E19cm−3(ここに、1E19とは1x1019の略記形))と不純物濃度が実質的にゼロに近い低濃度領域(1E14cm−3)の2つの領域に着目してモンテカルロ法数値実験データから緩和係数対エネルギー曲線2本を引出し、且つ、その変化分(ΔN(1E19−1E14cm−3))を示したのが図8である。本発明によれば低濃度領域の値が式(9)に相当し、変化分が式(8)の右辺におけるエネルギー依存性を表現する部分である。In FIG. 7, a region having a high substrate impurity concentration (1E19 cm −3 (here, 1E19 is an abbreviation of 1 × 10 19 )) and a low concentration region (1E14 cm −3 ) whose impurity concentration is substantially close to zero. FIG. 8 shows two relaxation coefficient versus energy curves drawn from Monte Carlo numerical experimental data, and the change (ΔN (1E19-1E14 cm −3 )), focusing on one region. According to the present invention, the value of the low density region corresponds to Equation (9), and the amount of change is the portion expressing the energy dependence on the right side of Equation (8).

本発明によれば、式(8)の右辺に示した不純物濃度依存性を表現する重み関数h(N)を分離することで、多次元関数表現データ構成を1次元構造データの積で単純化可能となる。不純物濃度依存性が強く出現するのは低電界移動度であり、移動度と散乱確率は逆数関係にある。そこで、低電界移動度の不純物濃度依存性の逆数を重み関数h(N)として適用する方法が考えられる。実験的に得られている移動度の不純物濃度依存性を表現する近似式として  According to the present invention, by separating the weight function h (N) representing the impurity concentration dependency shown on the right side of the equation (8), the multidimensional function representation data structure is simplified by the product of the one-dimensional structure data. It becomes possible. Impurity concentration dependence appears strongly in low electric field mobility, and mobility and scattering probability are in an inverse relationship. Therefore, a method is conceivable in which the inverse of the dependence of the low electric field mobility on the impurity concentration is applied as the weighting function h (N). As an approximate expression that expresses the impurity concentration dependence of mobility obtained experimentally

Figure 0005235124
ける移動度の値、SとNは不純物濃度依存性を表現するパラメータで、電子と正孔のそれぞれについて、事前に与えることができる。上記文献に示されたSとNの値を式(12)へ代入し、計算結果を図示したのが図9である。右下がりの滑らかな曲線で、不純物濃度が1E16cm−3以下で移動度は殆ど変化しない。図7に示したモンテカルロ法数値実験データにおいても、不純物濃度が1E16cm−3以下で緩和係数のエネルギー依存性曲線はほぼ一致し、低濃度領域で移動度が一定値となることを示している。
Figure 0005235124
The mobility values S and Nr are parameters that express the impurity concentration dependence, and can be given in advance for each of electrons and holes. The values of S and N r indicated in the above document is substituted into Equation (12), was shown the calculation result is Fig. It is a smooth curve that descends to the right, and the mobility hardly changes when the impurity concentration is 1E16 cm −3 or less. Also in the Monte Carlo method numerical experimental data shown in FIG. 7, when the impurity concentration is 1E16 cm −3 or less, the energy dependence curves of the relaxation coefficients are almost the same, and the mobility becomes a constant value in the low concentration region.

重み関数h(N)は式(8)と式(10)を見ればわかるように、低濃度領域を基準値とした任意の不純物濃度までの変化分に対する重みであるので、h(N)は式(12)を基に  As can be seen from the equations (8) and (10), the weight function h (N) is a weight for the change up to an arbitrary impurity concentration with the low concentration region as a reference value. Based on formula (12)

Figure 0005235124
と設定可能である。式(13)は比例関係を表すもので、次に比例定数の決定方法について述べる。図8に示された変化分、或いは、式(8)の右辺で示した変化分は不純物濃度で言うと1E14(cm−3)から1E19(cm−3)への変化分であることを考えれば、h(N)は
Figure 0005235124
And can be set. Equation (13) expresses a proportional relationship. Next, a method for determining a proportional constant will be described. The change shown in FIG. 8 or the change shown on the right side of Equation (8) can be considered to be a change from 1E14 (cm −3 ) to 1E19 (cm −3 ) in terms of impurity concentration. H (N) is

Figure 0005235124
と設定できる。式(14)を図示したのが、図10である。図から分かるようにh(N)は不純物濃度1E14(cm−3)を基準点(h=0)とし、1E19(cm−3)において1となるように規格化された重み関数である。低濃度領域(1E14−1E16cm−3)でhの値はほぼゼロである。式(14)のNの値として1E18(cm−3)と1E17(cm−3)を代入してhの値を計算し、運動量緩和係数のエネルギー依存性を、本発明になる乗算処理を用いて処理したものが図11である。1次元関数の乗算処理による図11と、モンテカルロ法計算による図7の結果とは妥当な一致をみることができる。即ち、本発明による1次元処理による計算処理方法が精度を落とすことなく、高速にデータ処理可能なことが示されたわけである。
Figure 0005235124
Can be set. FIG. 10 illustrates Expression (14). As can be seen from the figure, h (N) is a weighting function that is standardized to be 1 at 1E19 (cm −3 ) with the impurity concentration 1E14 (cm −3 ) as a reference point (h = 0). The value of h is almost zero in the low concentration region (1E14-1E16 cm −3 ). The value of h is calculated by substituting 1E18 (cm −3 ) and 1E17 (cm −3 ) as the value of N in Equation (14), and the energy dependence of the momentum relaxation coefficient is used in the multiplication process according to the present invention. FIG. 11 shows the result of processing. A reasonable agreement can be seen between the result of FIG. 11 by the multiplication process of the one-dimensional function and the result of FIG. 7 by the Monte Carlo calculation. That is, it has been shown that the calculation processing method by the one-dimensional processing according to the present invention can process data at high speed without reducing accuracy.

ゲート電界が3MV/cm以上の高電界になると、MOS反転層内の擬2次元電子ガスのエネルギー準位はバルク中のそれと同等あるいは更に高くなり、擬2次元電子ガスの状態からのズレが発生する。従って、図3並びに図6にて示されたゲート電界依存性の妥当性は3MV/cm程度と考えるべきである。一方、このような高電界が発生するとゲート酸化膜の絶縁破壊が物理的に発生する。従って、ゲート高電界領域の取扱については3MV/cmにおける状態を最終値として用いることが実用上妥当、かつ高速処理可能な方法でもある。  When the gate electric field becomes a high electric field of 3 MV / cm or more, the energy level of the quasi-two-dimensional electron gas in the MOS inversion layer becomes equal to or higher than that in the bulk, and a deviation from the state of the quasi-two-dimensional electron gas occurs. To do. Therefore, the validity of the gate electric field dependency shown in FIGS. 3 and 6 should be considered to be about 3 MV / cm. On the other hand, when such a high electric field is generated, dielectric breakdown of the gate oxide film physically occurs. Therefore, regarding the handling of the gate high electric field region, using the state at 3 MV / cm as the final value is a practically appropriate method and capable of high-speed processing.

実施例ではシリコン表面反転層における擬2次元電子ガスについて本特許の有効性を述べてきたが、正孔についても有効な処理方法であることは言うまでもない。また、式(2)の右辺第2項では不純物濃度Nの依存性がほとんど無いのでNを除去する形式を採用したが、Nを追加しても式(8)に示す如き1次元データ構造分解が可能である。依存性の小さい項目を陽に導入することは、計算負荷の増加から、可能な限り排除することが望ましい。  In the embodiments, the effectiveness of this patent has been described for the quasi-two-dimensional electron gas in the silicon surface inversion layer, but it goes without saying that it is also an effective treatment method for holes. The second term on the right-hand side of equation (2) has almost no dependency on impurity concentration N, so the form of removing N is adopted. However, even if N is added, one-dimensional data structure decomposition as shown in equation (8) is adopted. Is possible. It is desirable to explicitly introduce items with low dependency from the increase in calculation load as much as possible.

本発明によるデータ構造の簡略化は数値計算手順を簡便にするだけでなく、使用するデータそのものが簡単な構造となるため、データベース化して保存、読み出し、内挿処理による処理手順もまた有効である。更に、適用する材料は実施例で述べたシリコン半導体のみならず、シリコンとゲルマの混晶系半導体、シリコンとゲルマの混晶とヘテロ接合を介したシリコン半導体、化合物半導体、複数の材料系から成るヘテロ接合を有する化合物半導体にも適用可能である。  The simplification of the data structure according to the present invention not only simplifies the numerical calculation procedure, but also the data itself to be used has a simple structure. Therefore, the processing procedure by storing, reading, and interpolating in a database is also effective. . Furthermore, the material to be applied is not only the silicon semiconductor described in the embodiment, but also a silicon-german mixed crystal semiconductor, a silicon-german mixed crystal and a silicon semiconductor via a heterojunction, a compound semiconductor, and a plurality of material systems. The present invention can also be applied to a compound semiconductor having a heterojunction.

本発明の原理を示す図。  The figure which shows the principle of this invention. シリコンMOS反転層擬2次元電子ガスのエネルギー緩和係数のエネルギー依存性に関するモンテカルロ数値実験と本発明の実施例を示す図。  The figure which shows the Monte Carlo numerical experiment regarding the energy dependence of the energy relaxation coefficient of a silicon MOS inversion layer quasi-two-dimensional electron gas, and the Example of this invention. シリコンMOS反転層擬2次元電子ガスのエネルギー緩和係数のゲート電界依存性に関するモンテカルロ数値実験と本発明の実施例を示す図。  The figure which shows the Monte Carlo numerical experiment regarding the gate electric field dependence of the energy relaxation coefficient of silicon MOS inversion layer quasi-two-dimensional electron gas, and the Example of this invention. シリコンMOS反転層擬2次元電子ガスの運動量緩和係数のエネルギー依存性に関するモンテカルロ数値実験結果を示す図。  The figure which shows the Monte Carlo numerical experiment result regarding the energy dependence of the momentum relaxation coefficient of a silicon MOS inversion layer quasi-two-dimensional electron gas. シリコンMOS反転層擬2次元電子ガスの運動量緩和係数のエネルギー依存性に関する本発明の実施例を示す図。  The figure which shows the Example of this invention regarding the energy dependence of the momentum relaxation coefficient of a silicon MOS inversion layer quasi-two-dimensional electron gas. シリコンMOS反転層擬2次元電子ガスの運動量緩和係数のゲート電界依存性に関するモンテカルロ数値実験と本発明の実施例を示す図。  The figure which shows the Monte Carlo numerical experiment regarding the gate electric field dependence of the momentum relaxation coefficient of a silicon MOS inversion layer quasi-two-dimensional electron gas, and the Example of this invention. バルクシリコン中の電子の運動量緩和係数のエネルギー依存性に関するモンテカルロ数値実験結果を示す図。  The figure which shows the Monte Carlo numerical experiment result regarding the energy dependence of the momentum relaxation coefficient of the electron in bulk silicon. バルクシリコン中の電子の運動量緩和係数のエネルギー依存性の近似曲線(モンテカルロ数値実験データを滑らかに内挿する近似曲線)と2本の曲線により示される2組のデータの変化分を示す図。  The figure which shows the change part of two sets of data shown by the approximate curve (approximate curve which interpolates a Monte Carlo numerical experimental data smoothly) and the two curves of the energy dependence of the momentum relaxation coefficient of the electron in bulk silicon. シリコン中の電子と正孔の低電界移動度の不純物濃度依存性を示す図。  The figure which shows the impurity concentration dependence of the low electric field mobility of the electron and hole in silicon. 運動量緩和係数の不純物濃度依存重み関数を示す図。  The figure which shows the impurity concentration dependence weight function of a momentum relaxation coefficient. 本発明になる運動量緩和係数の不純物濃度、エネルギー依存性を示す図。  The figure which shows the impurity concentration and energy dependence of the momentum relaxation coefficient which become this invention.

Claims (6)

半導体装置内部のキャリア輸送をコンピューターに求めさせるためのプログラムであって、前記コンピューターに実行されると、静電ポテンシャル分布を求めるポアソン方程式とキャリア分布を求める電流連続の式を解く手段と、キャリアのエネルギーを求めるエネルギーバランス方程式を解く手段と、運動量を求める運動量バランス方程式を解く手段と、キャリアの走行に対し垂直方向成分の電界成分を計算する手段と、エネルギー緩和係数が、前記キャリアエネルギーと、前記キャリアの走行に対し垂直方向成分の電界との関数で表現されることによる計算手段と、運動量緩和係数が、前記キャリアエネルギーと、前記キャリアの走行に対し垂直方向成分の電界と、基板不純物濃度との関数で表現されることによる計算手段と、前記エネルギー緩和係数が、前記キャリアの走行に対し垂直方向成分の電界成分が殆どないバルク状態におけるエネルギー緩和係数と、前記キャリアの走行に対し垂直方向成分の電界成分により強く影響を受ける擬2次元電子ガス状態におけるエネルギー緩和係数との和、またはいずれか大きい方で決定され、及び/または、前記運動量緩和係数が、前記キャリアの走行に対し垂直方向成分の電界成分が殆どないバルク状態における運動量緩和係数と、前記キャリアの走行に対し垂直方向成分の電界成分により強く影響を受ける擬2次元電子ガス状態における運動量緩和係数との和、またはいずれか大きい方で決定される計算手段として機能し前記バルク状態におけるエネルギー緩和係数と前記バルク状態における運動量緩和係数と前記キャリアエネルギーと前記基板不純物濃度の関数で決定され前記擬2次元電子ガス状態におけるエネルギー緩和係数と前記擬2次元電子ガス状態における運動量緩和係数と前記キャリアエネルギーと前記キャリアの走行に対し垂直方向成分の電界成分の関数で決定されることを特徴とするプログラム。 A program for causing a computer to calculate carrier transport inside a semiconductor device, and when executed by the computer, means for solving a Poisson equation for obtaining an electrostatic potential distribution and a current continuity equation for obtaining a carrier distribution , Means for solving an energy balance equation for obtaining energy; means for solving a momentum balance equation for obtaining momentum; means for calculating an electric field component of a vertical component with respect to carrier travel; and an energy relaxation coefficient comprising the carrier energy, The calculation means by being expressed as a function of the electric field of the vertical component with respect to the carrier traveling, the momentum relaxation coefficient, the carrier energy, the electric field of the vertical component with respect to the carrier traveling, and the substrate impurity concentration A calculation means represented by a function of The quasi-two-dimensional electron gas is strongly influenced by the energy relaxation coefficient in the bulk state in which the electric field component of the vertical component with respect to the traveling of the carrier is hardly present, and the electric field component of the vertical component with respect to the traveling of the carrier. And / or the momentum relaxation coefficient in the bulk state in which there is almost no electric field component in the vertical direction with respect to the traveling of the carrier. functions as the calculation means is determined the sum of the momentum relaxation coefficient in quasi-two-dimensional electron gas states strongly affected by the electric field component in the vertical direction component with respect to the travel of the carrier or any larger, the bulk wherein the momentum relaxation coefficient in the bulk state energy relaxation coefficient in the carrier Wherein the energy is determined as a function of the substrate impurity concentration, a direction perpendicular to the travel of the carrier and the carrier energy and momentum relaxation coefficient in the quasi-two-dimensional electron gas states with energy relaxation coefficient in the quasi-two-dimensional electron gas states A program characterized by being determined as a function of the electric field component of the component. 前記バルク状態におけるエネルギー緩和係数は前記キャリアエネルギーのみで表現されることを特徴とする請求項1記載のプログラム。 Claim 1, wherein the program energy relaxation coefficient in the bulk state is characterized by being represented only by the carrier energy. 前記バルク状態における運動量緩和係数が、前記キャリアエネルギーで記述される緩和係数と前記基板不純物濃度の関数で記述される係数との積、または前記キャリアエネルギーで記述される緩和係数と少なくとも一つ以上のエネルギー状態で記述される緩和係数との和、または該乗算と該和算の和で表現されることを特徴とする請求項1または2記載のプログラム。 Momentum relaxation factor in the bulk state, wherein the relaxation factor described by the carrier energy and coefficients described by a function of the substrate impurity concentration product or alleviating factor that is described by the carrier energy and at least one, 3. The program according to claim 1, wherein the program is expressed by a sum of relaxation coefficients described in energy states, or a sum of the multiplication and the summation. 前記擬2次元電子ガス状態におけるエネルギー緩和係数と運動量緩和係数とが、前記キャリアエネルギーで記述される緩和係数と前記キャリアの走行に対し垂直方向成分の電界成分の関数で記述される係数との積で表現されることを特徴とする請求項1または2または3記載のプログラム。 Product of the quasi-two-dimensional electron gas states in the energy relaxation factor and the momentum relaxation factor, factor described by a function of the electric field component in the vertical direction component with respect to the traveling relaxation factor and the carrier described by the carrier energy The program according to claim 1, 2, or 3. 少なくとも2つ以上の材料から成る半導体において装置内部のキャリア輸送をコンピューターに求めさせるためのプログラムであって、非混晶状態の前記エネルギー緩和係数と前記運動量緩和係数は、材料の混晶比率の関数で記述される値を乗算することにより決定されることを特徴とする請求項1〜4記載のプログラム。 A program for causing a computer to calculate carrier transport inside a device in a semiconductor composed of at least two or more materials, wherein the energy relaxation coefficient and the momentum relaxation coefficient in a non-mixed crystal state are functions of the mixed crystal ratio of the material. in claims 1 to 4, wherein the program characterized Rukoto determined by multiplying the stated values. 請求項1〜5記載のプログラムにより生成された前記エネルギー緩和係数と前記運動量緩和係数とをデータベース化し、前記データベース化された前記エネルギー緩和係数と前記運動量緩和係数とを元に輸送問題を解析することを特徴とするプログラム。 A database of the energy relaxation coefficient and the momentum relaxation coefficient generated by the program according to claim 1, and a transportation problem is analyzed based on the energy relaxation coefficient and the momentum relaxation coefficient that are stored in the database. A program characterized by
JP2008282458A 2008-10-07 2008-10-07 Semiconductor device simulation program Active JP5235124B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008282458A JP5235124B2 (en) 2008-10-07 2008-10-07 Semiconductor device simulation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008282458A JP5235124B2 (en) 2008-10-07 2008-10-07 Semiconductor device simulation program

Publications (2)

Publication Number Publication Date
JP2010092445A JP2010092445A (en) 2010-04-22
JP5235124B2 true JP5235124B2 (en) 2013-07-10

Family

ID=42255063

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008282458A Active JP5235124B2 (en) 2008-10-07 2008-10-07 Semiconductor device simulation program

Country Status (1)

Country Link
JP (1) JP5235124B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6118079B2 (en) * 2012-11-22 2017-04-19 アドバンスソフト株式会社 Semiconductor device simulation calculation processing method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06167926A (en) * 1992-12-01 1994-06-14 Sony Corp Energy transport type device simulating method
JP3569750B2 (en) * 1995-06-07 2004-09-29 沖電気工業株式会社 Semiconductor device simulation method and apparatus
JP2005228894A (en) * 2004-02-12 2005-08-25 Semiconductor Leading Edge Technologies Inc Semiconductor device simulation method
JP2006113749A (en) * 2004-10-13 2006-04-27 Fujitsu Ltd Semiconductor-designing model and semiconductor design device using it

Also Published As

Publication number Publication date
JP2010092445A (en) 2010-04-22

Similar Documents

Publication Publication Date Title
Nedjalkov et al. Wigner function approach
Tsuchiya et al. A quantum-corrected Monte Carlo study on quasi-ballistic transport in nanoscale MOSFETs
Martinez et al. NEGF simulations of a junctionless Si gate-all-around nanowire transistor with discrete dopants
Afzalian et al. ATOMOS: An ATomistic MOdelling Solver for dissipative DFT transport in ultra-scaled HfS 2 and Black phosphorus MOSFETs
Vyas Theoretical Study of Quantum Transport in Realistic Semiconductor Devices
Jin et al. Theoretical study of carrier transport in silicon nanowire transistors based on the multisubband Boltzmann transport equation
Tyaginov Physics-based modeling of hot-carrier degradation
Watling et al. The impact of interface roughness scattering and degeneracy in relaxed and strained Si n-channel MOSFETs
Sampedro et al. Multi-Subband Ensemble Monte Carlo simulation of bulk MOSFETs for the 32 nm-node and beyond
JP5235124B2 (en) Semiconductor device simulation program
Abdallah et al. A deterministic solver for a hybrid quantum-classical transport model in nanoMOSFETs
Macucci et al. Status and perspectives of nanoscale device modelling
Jeong Quantum-mechanical analysis of amorphous oxide-based thin-film transistors
Bescond Quantum transport in semiconductor nanowires
Fiori et al. Code for the 3D simulation of nanoscale semiconductor devices, including drift-diffusion and ballistic transport in 1D and 2D subbands, and 3D tunneling
Chen et al. Two-dimensional quantum mechanical simulation of electron transport in nano-scaled Si-based MOSFETs
Koba et al. Quantum transport simulation of nanoscale semiconductor devices based on Wigner Monte Carlo approach
Ezaki et al. Self-consistent quantum mechanical Monte Carlo MOSFET device simulation
Ferry et al. Physics of silicon nanodevices
Zhouyin et al. Automatic differentiable nonequilibrium Green's function formalism: An end-to-end differentiable quantum transport simulator
Heinz et al. Full quantum simulation of silicon-on-insulator single-electron devices
Soares Ensemble Monte Carlo simulation of hole transport in Si, Ge and SiGe Alloys
Marchetti et al. Spin decoherence in n-type GaAs: The effectiveness of the third-body rejection method for electron-electron scattering
Xiao-Yan et al. An analytical model for the electron effective mobility in a strained silicon inversion layer
Charles Efficient inelastic scattering in atomistic tight binding

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110915

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20110915

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20110915

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20111220

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20111220

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20121128

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20121204

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130201

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130226

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130325

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 5235124

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160405

Year of fee payment: 3

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250