[go: up one dir, main page]

JP2010061190A - Model parameter extraction device and model parameter extraction method - Google Patents

Model parameter extraction device and model parameter extraction method Download PDF

Info

Publication number
JP2010061190A
JP2010061190A JP2008223247A JP2008223247A JP2010061190A JP 2010061190 A JP2010061190 A JP 2010061190A JP 2008223247 A JP2008223247 A JP 2008223247A JP 2008223247 A JP2008223247 A JP 2008223247A JP 2010061190 A JP2010061190 A JP 2010061190A
Authority
JP
Japan
Prior art keywords
model parameter
electrical characteristics
unit
approximate function
measurement point
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.)
Withdrawn
Application number
JP2008223247A
Other languages
Japanese (ja)
Inventor
Naoki Wakita
田 直 樹 脇
Kuniaki Utsumi
海 邦 朗 内
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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP2008223247A priority Critical patent/JP2010061190A/en
Publication of JP2010061190A publication Critical patent/JP2010061190A/en
Withdrawn legal-status Critical Current

Links

Images

Landscapes

  • Insulated Gate Type Field-Effect Transistor (AREA)

Abstract

【課題】回路シミュレーションの精度の向上。
【解決手段】測定点の電気特性を示す測定データを入力する入力部10と、入力された測定データが存在している測定点について、入力された測定データに基づいて、測定点の電気特性を近似する第1近似関数式を作成する第1近似関数式作成部201と、入力された測定データが不足している測定点について、第1近似関数式作成部201によって作成された第1近似関数式を修正して、測定点の電気特性を近似する第2近似関数式を作成する第2近似関数式作成部202と、第1および第2近似関数式作成部201,202によって作成された第1および第2近似関数式に基づいて、測定点の電気特性の近似値を算出する近似値算出部30と、算出された近似値を含むモデルパラメータセットを生成するモデルパラメータセット生成部80と、生成されたモデルパラメータセットを出力する出力部90と、を備える。
【選択図】図1
To improve the accuracy of circuit simulation.
An input unit for inputting measurement data indicating the electrical characteristics of a measurement point, and the electrical characteristics of the measurement point for the measurement point where the input measurement data exists are determined based on the input measurement data. A first approximate function formula creating unit 201 that creates a first approximate function formula to be approximated, and a first approximate function created by the first approximate function formula creating unit 201 for measurement points for which input measurement data is insufficient The second approximate function formula creating unit 202 that modifies the formula to create a second approximate function formula that approximates the electrical characteristics of the measurement point, and the first and second approximate function formula creating units 201 and 202 created by the first and second approximate function formula creating units 201 and 202. Based on the first and second approximate function formulas, an approximate value calculation unit 30 that calculates an approximate value of the electrical characteristic of the measurement point, and a model parameter set generation unit 8 that generates a model parameter set including the calculated approximate value Comprising the, an output unit 90 for outputting the generated model parameter set, a.
[Selection] Figure 1

Description

本発明は、モデルパラメータ抽出装置およびモデルパラメータ抽出方法に関し、特に、SPICE(Simulation Program with Integrated Circuit Emphasis)モデルのモデルパラメータを抽出するモデルパラメータ抽出装置およびモデルパラメータ抽出方法に関する。   The present invention relates to a model parameter extraction apparatus and a model parameter extraction method, and more particularly to a model parameter extraction apparatus and a model parameter extraction method for extracting model parameters of a SPICE (Simulation Program with Integrated Circuit Emphasis) model.

MOSFET(Metal Oxide Semiconducter Field Effect Transistor)のSPICEモデルパラメータセットには、グローバル(Global)モデルとビニング(binning)モデルとがある。グローバルモデルは、MOSFETのゲート長Lおよびゲート幅Wに対するL−W平面上の全領域の電気特性を1つのモデル式で表現するモデルである。一方、ビニングモデルは、L−W平面を複数の局所領域に分割して、MOSFETのゲート長Lおよびゲート幅Wに対するL−W平面上の各局所領域での電気特性を各局所領域のみで有効な局所モデル式で表現するモデルである。この各局所領域はビン(bin)と呼ばれ、複数のビンに分割されたモデルを作成することはビニングと呼ばれる。   The SPICE model parameter set of a MOSFET (Metal Oxide Semiconductor Field Effect Transistor) includes a global model and a binning model. The global model is a model that expresses the electrical characteristics of the entire region on the LW plane with respect to the gate length L and gate width W of the MOSFET by one model formula. On the other hand, the binning model divides the LW plane into a plurality of local regions, and the electrical characteristics in each local region on the LW plane with respect to the gate length L and gate width W of the MOSFET are effective only in each local region. This model is expressed by a local model expression. Each local region is called a bin, and creating a model divided into a plurality of bins is called binning.

しかしながら、一般的には、微細なMOSFETの電気特性は、ゲート長Lおよびゲート幅Wに対して非線形であり、且つ、互いに複雑な依存性を有する。したがって、グローバルモデルでは、MOSFETの電気特性を高精度に表現することは難しい。   However, in general, the electrical characteristics of a fine MOSFET are nonlinear with respect to the gate length L and the gate width W, and have complicated dependencies on each other. Therefore, in the global model, it is difficult to express the electrical characteristics of the MOSFET with high accuracy.

これに対して、ビニングモデルでは、各局所モデル式が各ビンのみを表現するので、MOSFETの電気特性を高精度に表現することができる。たとえば、回路シミュレーション実行時に、所定のゲート長Lおよびゲート幅Wの寸法に基づいて、適切な局所モデルを選択し、使用することによって、L−W平面上の全領域の電気特性を高精度に表現することができる。   On the other hand, in the binning model, each local model expression represents only each bin, so that the electrical characteristics of the MOSFET can be expressed with high accuracy. For example, when a circuit simulation is executed, an appropriate local model is selected and used on the basis of the dimensions of a predetermined gate length L and gate width W, so that the electrical characteristics of the entire region on the LW plane can be accurately determined. Can be expressed.

ここで、一般的なビニングモデルの作成手順について説明する。   Here, a general procedure for creating a binning model will be described.

はじめに、L/W領域全体を対象とするグローバルモデルを作成する。このグローバルモデルの精度は、高精度でなくても良い。次に、このグローバルモデルに基づいて、格子状に配置された各測定点の局所モデルを作成する。この局所モデルの精度は、できるだけ高いことが好ましい。次に、各ビンの局所モデルを作成する。この局所モデルは、各ビンを囲む4つの測定点のモデルパラメータを用いて所定の連立方程式を解くことによって作成される。   First, a global model for the entire L / W area is created. The accuracy of this global model need not be high. Next, based on this global model, a local model of each measurement point arranged in a grid is created. The accuracy of this local model is preferably as high as possible. Next, a local model for each bin is created. This local model is created by solving a predetermined simultaneous equation using model parameters of four measurement points surrounding each bin.

ここで、グローバルモデルから局所モデルを作成する方法について説明する。なお、モデルパラメータU00の例を説明する。各測定点のL/Wの寸法(L,W)が(L,W),(L,W),(L,W),(L,W)であり、各測定点のモデルパラメータU00の値が(U01,U02,U03,U04)であり、各測定点によって囲まれるビンが“bin_A”であるとする。これらの4つのモデルパラメータU00の値(U01,U02,U03,U04)を使用して、式1乃至式4の4元連立1次方程式を立てる。この連立方程式の解(U01,U02,U03,U04)が“bin_A”のモデルパラメータとなる。 Here, a method for creating a local model from a global model will be described. Incidentally, an example of a model parameter U 00. The L / W dimensions (L i , W i ) at each measurement point are (L 1 , W 1 ), (L 2 , W 1 ), (L 1 , W 2 ), (L 2 , W 2 ). Assume that the value of the model parameter U 00 at each measurement point is (U 01 , U 02 , U 03 , U 04 ), and the bin surrounded by each measurement point is “bin_A”. Using the values (U 01 , U 02 , U 03 , U 04 ) of these four model parameters U 00 , the quaternary simultaneous linear equations of Expressions 1 to 4 are established. The solution (U 01 , U 02 , U 03 , U 04 ) of the simultaneous equations is a model parameter of “bin_A”.

〔式1〕 U01=U+LU/L+WU/W+PU(L*W
〔式2〕 U02=U+LU/L+WU/W+PU(L*W
〔式3〕 U03=U+LU/L+WU/W+PU(L*W
〔式4〕 U04=U+LU/L+WU/W+PU(L*W
すなわち、回路シミュレーション実行時には、“bin_A”の任意のゲート長Lおよびゲート幅Wの寸法に対して、式5の解が算出される。
〔式5〕 U00=U+LU/L+WU/W+PU(L*W)
しかしながら、ビニングモデルでは、各ビンのモデルパラメータは、そのビンを囲む4つの測定点のモデルパラメータを用いて式1乃至式4に示されるような連立方程式を解くことによって作成される。したがって、ビンのモデルパラメータを作成するためには、そのビンを囲む4つの測定点の電気特性を示す測定データが必要である。すなわち、ビニングモデルには、全ての測定点の測定データが必要である。ところが、一般的には、TEG(Test Element Group)の欠損等によって、測定データが不足する(所定数以上の測定点にTEGが存在しない)場合がある。このような場合には、1つのビンの面積が大きくなり、各ビンのモデルパラメータとMOSFETの電気特性との間の誤差も大きくなる。その結果、回路シミュレーションの精度が低下するという問題が発生する。
[Expression 1] U 01 = U 0 + LU 0 / L 1 + WU 0 / W 1 + PU 0 (L 1 * W 1 )
[Formula 2] U 02 = U 0 + LU 0 / L 1 + WU 0 / W 2 + PU 0 (L 1 * W 2 )
[Equation 3] U 03 = U 0 + LU 0 / L 2 + WU 0 / W 1 + PU 0 (L 2 * W 1 )
[Equation 4] U 04 = U 0 + LU 0 / L 2 + WU 0 / W 2 + PU 0 (L 2 * W 2)
That is, at the time of executing the circuit simulation, the solution of Equation 5 is calculated for an arbitrary gate length L and gate width W of “bin_A”.
[Formula 5] U 00 = U 0 + LU 0 / L + WU 0 / W + PU 0 (L * W)
However, in the binning model, the model parameters of each bin are created by solving simultaneous equations as shown in Equations 1 to 4 using the model parameters of the four measurement points surrounding the bin. Therefore, in order to create a model parameter for a bin, measurement data indicating the electrical characteristics of four measurement points surrounding the bin is required. That is, the binning model requires measurement data at all measurement points. However, generally, measurement data may be insufficient (a TEG does not exist at a predetermined number or more measurement points) due to a lack of TEG (Test Element Group) or the like. In such a case, the area of one bin increases, and the error between the model parameters of each bin and the electrical characteristics of the MOSFET also increases. As a result, there arises a problem that the accuracy of circuit simulation is lowered.

また、測定データが充足する(所定数以上の測定点にTEGが存在する)場合であっても、測定点間で電気特性が急激に変化するような場合には、ビンのモデルパラメータとMOSFETの電気特性との間の誤差が大きくなる。その結果、回路シミュレーションの精度が低下するという問題が発生する。
特開2004−273903号公報
Even if the measurement data is satisfied (the TEG exists at a predetermined number or more of measurement points), if the electrical characteristics change suddenly between the measurement points, the model parameter of the bin and the MOSFET The error between the electrical characteristics increases. As a result, there arises a problem that the accuracy of circuit simulation is lowered.
JP 2004-273903 A

本発明の目的は、回路シミュレーションの精度を向上させるモデルパラメータ抽出装置およびモデルパラメータ抽出方法を提供することである。   An object of the present invention is to provide a model parameter extraction device and a model parameter extraction method that improve the accuracy of circuit simulation.

本発明の第1態様によれば、
トランジスタのゲート長方向およびゲート幅方向に対して格子状に配置された測定点によって囲まれるシミュレーションモデルのモデルパラメータを抽出するモデルパラメータ抽出装置であって、
前記測定点の電気特性を示す測定データを入力する入力部と、
前記入力部によって入力された測定データが充足している測定点について、前記入力部によって入力された測定データに基づいて、前記測定点の電気特性を近似する第1近似関数式を作成する第1近似関数式作成部と、
前記入力部によって入力された測定データが不足している測定点について、前記第1近似関数式作成部によって作成された第1近似関数式を修正して、前記測定点の電気特性を近似する第2近似関数式を作成する第2近似関数式作成部と、
前記第1および第2近似関数式作成部によって作成された第1および第2近似関数式に基づいて、前記測定点の電気特性の近似値を算出する近似値算出部と、
前記近似値算出部によって算出された近似値を含むモデルパラメータセットを生成するモデルパラメータセット生成部と、
前記モデルパラメータセット生成部によって生成されたモデルパラメータセットを出力する出力部と、を備えることを特徴とするモデルパラメータ抽出装置が提供される。
According to the first aspect of the present invention,
A model parameter extraction device for extracting model parameters of a simulation model surrounded by measurement points arranged in a grid pattern in a gate length direction and a gate width direction of a transistor,
An input unit for inputting measurement data indicating the electrical characteristics of the measurement point;
A first approximation function formula that approximates the electrical characteristics of the measurement point is generated based on the measurement data input by the input unit for the measurement point that is satisfied by the measurement data input by the input unit. An approximate function formula creation unit;
For a measurement point for which measurement data input by the input unit is insufficient, the first approximate function formula created by the first approximate function formula creation unit is modified to approximate the electrical characteristics of the measurement point. A second approximate function formula creating unit for creating a two approximate function formula;
An approximate value calculation unit that calculates an approximate value of the electrical characteristics of the measurement point based on the first and second approximate function formulas created by the first and second approximate function formula creation units;
A model parameter set generation unit for generating a model parameter set including the approximate value calculated by the approximate value calculation unit;
And a model parameter extracting apparatus comprising: an output unit that outputs the model parameter set generated by the model parameter set generating unit.

本発明の第2態様によれば、
トランジスタのゲート長方向およびゲート幅方向に対して格子状に配置された測定点によって囲まれるシミュレーションモデルのモデルパラメータを抽出するモデルパラメータ抽出装置であって、
前記ゲート幅または前記ゲート長方向の電気特性を所定間隔毎に読み取る読取部と、
前記読取部によって読み取られた電気特性が許容範囲外である場合に、その電気特性に対応する前記ゲート幅または前記ゲート長に前記測定点を設定する測定点設定部と、
前記測定点設定部によって設定された測定点の電気特性の推定値を算出する推定値算出部と、
前記推定値算出部によって算出された推定値を含むモデルパラメータセットを生成するモデルパラメータセット生成部と、
前記モデルパラメータセット生成部によって生成されたモデルパラメータセットを出力する出力部と、を備えることを特徴とするモデルパラメータ抽出装置が提供される。
According to a second aspect of the invention,
A model parameter extraction device for extracting model parameters of a simulation model surrounded by measurement points arranged in a grid in the gate length direction and gate width direction of a transistor,
A reading unit that reads electrical characteristics in the gate width or gate length direction at predetermined intervals;
A measurement point setting unit configured to set the measurement point to the gate width or the gate length corresponding to the electrical characteristic when the electrical characteristic read by the reading unit is outside an allowable range;
An estimated value calculation unit for calculating an estimated value of the electrical characteristics of the measurement point set by the measurement point setting unit;
A model parameter set generator for generating a model parameter set including the estimated value calculated by the estimated value calculator;
An output unit that outputs the model parameter set generated by the model parameter set generation unit.

本発明の第3態様によれば、
トランジスタのゲート長方向およびゲート幅方向に対して格子状に配置された測定点によって囲まれるシミュレーションモデルのモデルパラメータを抽出するモデルパラメータ抽出方法であって、
前記測定点の電気特性を示す測定データを入力し、
前記測定データが充足している測定点について、前記測定データに基づいて、前記測定点の電気特性を近似する第1近似関数式を作成し、
前記測定データが不足している測定点について、前記第1近似関数式を修正して、前記測定点の電気特性を近似する第2近似関数式を作成し、
前記第1および第2近似関数式に基づいて、前記測定点の電気特性の近似値を算出し、
前記近似値を含むモデルパラメータセットを生成し、
前記モデルパラメータセットを出力することを特徴とするモデルパラメータ抽出方法が提供される。
According to a third aspect of the invention,
A model parameter extraction method for extracting a model parameter of a simulation model surrounded by measurement points arranged in a grid pattern in a gate length direction and a gate width direction of a transistor,
Input measurement data indicating the electrical characteristics of the measurement point,
For the measurement points satisfied by the measurement data, based on the measurement data, create a first approximate function equation that approximates the electrical characteristics of the measurement points,
For the measurement point for which the measurement data is insufficient, the first approximate function formula is corrected to create a second approximate function formula that approximates the electrical characteristics of the measurement point,
Based on the first and second approximation function equations, an approximate value of the electrical characteristics of the measurement point is calculated,
Generating a model parameter set including the approximation,
A model parameter extraction method is provided that outputs the model parameter set.

本発明によれば、回路シミュレーションの精度を向上させることができる。   According to the present invention, the accuracy of circuit simulation can be improved.

以下、本発明の実施例について図面を参照して説明する。なお、以下の実施例は、本発明の実施の一形態であって、本発明の範囲を限定するものではない。   Embodiments of the present invention will be described below with reference to the drawings. The following examples are one embodiment of the present invention and do not limit the scope of the present invention.

はじめに、本発明の実施例1について説明する。本発明の実施例1は、TEGが存在していない測定点の電気特性の近似値や推定値を算出して、モデルパラメータを抽出するモデルパラメータ抽出装置の例である。   First, Example 1 of the present invention will be described. The first embodiment of the present invention is an example of a model parameter extracting apparatus that calculates an approximate value or an estimated value of an electrical characteristic at a measurement point where no TEG exists and extracts a model parameter.

まず、本発明の実施例1に係るモデルパラメータ抽出装置の構成について図1および図2を参照して説明する。図1は、本発明の実施例1に係るモデルパラメータ抽出装置1の構成を示すブロック図である。図2は、本発明の実施例1に係る測定データの概略を示す概略図である。   First, the configuration of the model parameter extracting apparatus according to the first embodiment of the present invention will be described with reference to FIG. 1 and FIG. FIG. 1 is a block diagram illustrating a configuration of a model parameter extracting apparatus 1 according to the first embodiment of the present invention. FIG. 2 is a schematic diagram showing an outline of measurement data according to Example 1 of the present invention.

図1に示すように、本発明の実施例1に係るモデルパラメータ抽出装置1は、入力部10と、近似関数式作成部20と、近似値算出部30と、読取部40と、測定点設定部50と、ビニングモデル作成部60と、推定値算出部70と、モデルパラメータセット生成部80と、出力部90と、を備えている。   As shown in FIG. 1, the model parameter extraction apparatus 1 according to the first embodiment of the present invention includes an input unit 10, an approximate function expression creation unit 20, an approximate value calculation unit 30, a reading unit 40, and a measurement point setting. Unit 50, binning model creation unit 60, estimated value calculation unit 70, model parameter set generation unit 80, and output unit 90.

図1に示すように、入力部10は、近似関数式作成部20および入力装置(図示せず)に接続されている。また、入力部10は、ユーザが入力装置(図示せず)を用いて入力した測定点の電気特性(たとえば、I−V特性カーブ、閾値電圧Vth、ドレイン電流Ion)を示す測定データを入力するようになっている。 As shown in FIG. 1, the input unit 10 is connected to an approximate function equation creation unit 20 and an input device (not shown). The input unit 10 also has measurement data indicating electrical characteristics (for example, an IV characteristic curve, a threshold voltage V th , and a drain current I on ) at a measurement point input by a user using an input device (not shown). It comes to input.

図2に示すように、測定データは、MOSFETのゲート幅W(=W1乃至W4)のそれぞれにおいて、MOSFETのゲート長L(=0.04乃至10)ごとにTEGが存在している測定点とTEGが存在していない測定点を含んでいる。なお、図2において、「o」は、TEGが存在している測定点を示し、「x」はTEGが存在していない測定点を示している。すなわち、図2は、ゲート幅W(=W)においてゲート長L(=0.1)のTEGが存在し、ゲート幅W(=W)においてゲート長L(=0.1)のTEGが存在していない。 As shown in FIG. 2, the measurement data includes the measurement points at which TEGs exist for each MOSFET gate length L (= 0.04 to 10) in each of the MOSFET gate widths W (= W1 to W4). The measurement point where TEG does not exist is included. In FIG. 2, “o” indicates a measurement point where the TEG exists, and “x” indicates a measurement point where the TEG does not exist. That is, in FIG. 2, there is a TEG having a gate length L (= 0.1) at a gate width W (= W 0 ), and a TEG having a gate length L (= 0.1) at a gate width W (= W 1 ). Does not exist.

ここで、図2の測定データは、ゲート幅Wが4段階であり、ゲート長Lが13段階である。したがって、図2の測定データに対応するビニングモデルを作成するためには、52(=4×13)個のTEGが必要である。しかしながら、図2の測定データでは、存在しているTEGの数は19個であるので、ビニングモデルを作成するには不足していることになる。   Here, in the measurement data of FIG. 2, the gate width W has four stages and the gate length L has thirteen stages. Therefore, in order to create a binning model corresponding to the measurement data of FIG. 2, 52 (= 4 × 13) TEGs are required. However, in the measurement data of FIG. 2, since there are 19 TEGs, it is insufficient to create a binning model.

図1に示すように、近似関数式作成部20は、第1近似関数式作成部201と、第2近似関数式作成部202と、を備えている。   As shown in FIG. 1, the approximate function formula creating unit 20 includes a first approximate function formula creating unit 201 and a second approximate function formula creating unit 202.

図1に示すように、第1近似関数式作成部201は、入力部10、第2近似関数式作成部202、近似値算出部30、および読取部40に接続されている。また、第1近似関数式作成部201は、入力部10によって入力された測定データのうち、同一のゲート幅Wで、異なるゲート長Lの多くのTEGが存在している測定点について、ゲート幅Wの電気特性をゲート長Lの関数として近似するn(nは2以上の整数)次多項式(以下、「第1近似関数式」という)を作成するようになっている。なお、nは、要求される近似の精度によって決定される値であって、大きな値ほど精度が高くなる。   As shown in FIG. 1, the first approximate function formula creating unit 201 is connected to the input unit 10, the second approximate function formula creating unit 202, the approximate value calculating unit 30, and the reading unit 40. In addition, the first approximate function formula creation unit 201 determines the gate width of the measurement data input by the input unit 10 for the measurement points where many TEGs with the same gate width W and different gate lengths L exist. An n-th order polynomial (hereinafter referred to as “first approximate function expression”) that approximates the electrical characteristics of W as a function of the gate length L is created. Note that n is a value determined by the required approximation accuracy, and the larger the value, the higher the accuracy.

図1に示すように、第2近似関数式作成部202は、第1近似関数式作成部201、近似値算出部30、および読取部40に接続されている。また、第2近似関数式作成部202は、入力部10によって入力された測定データのうちTEGが存在していない測定点について、第1近似関数式作成部201によって作成された第1近似関数式の係数を修正して、ゲート幅Wの電気特性をゲート長Lの関数として近似するn次多項式(以下、「第2近似関数式」という)を作成するようになっている。   As shown in FIG. 1, the second approximate function formula creating unit 202 is connected to the first approximate function formula creating unit 201, the approximate value calculating unit 30, and the reading unit 40. In addition, the second approximate function formula creating unit 202 creates the first approximate function formula created by the first approximate function formula creating unit 201 for the measurement points where the TEG does not exist in the measurement data input by the input unit 10. Is modified to create an n-order polynomial (hereinafter referred to as “second approximate function expression”) that approximates the electrical characteristics of the gate width W as a function of the gate length L.

図1に示すように、近似値算出部30は、近似関数式作成部20およびモデルパラメータセット生成部80に接続されている。また、近似値算出部30は、第1および第2近似関数式作成部201,202によって作成された第1および第2近似関数式に基づいて、測定点の電気特性の近似値を算出するようになっている。   As shown in FIG. 1, the approximate value calculation unit 30 is connected to the approximate function formula creation unit 20 and the model parameter set generation unit 80. In addition, the approximate value calculation unit 30 calculates the approximate value of the electrical characteristics of the measurement point based on the first and second approximate function formulas created by the first and second approximate function formula creation units 201 and 202. It has become.

図1に示すように、読取部40は、近似関数式作成部20および測定点設定部50に接続されている。また、読取部40は、第1および第2近似関数式作成部201,202によって作成された第1および第2近似関数式に基づいて、入力部10によって入力された測定データにTEGが存在していない点の電気特性を読み取り、後述のビニングモデルとの誤差が許容値ε以上である場合に、その点の座標および電気特性を保持するようになっている。   As shown in FIG. 1, the reading unit 40 is connected to the approximate function formula creation unit 20 and the measurement point setting unit 50. In addition, the reading unit 40 has a TEG in the measurement data input by the input unit 10 based on the first and second approximate function formulas created by the first and second approximate function formula creation units 201 and 202. The electrical characteristics of a point that has not been read are read, and the coordinates and electrical characteristics of the point are held when an error from a binning model described later is equal to or greater than an allowable value ε.

図1に示すように、測定点設定部50は、読取部40、ビニングモデル作成部60、および推定値算出部70に接続されている。また、測定点設定部50は、読取部40によって保持された座標に測定点を設定するようになっている。   As shown in FIG. 1, the measurement point setting unit 50 is connected to a reading unit 40, a binning model creation unit 60, and an estimated value calculation unit 70. The measurement point setting unit 50 sets the measurement point at the coordinates held by the reading unit 40.

図1に示すように、ビニングモデル作成部60は、測定点設定部50に接続されている。また、ビニングモデル作成部60は、入力部10によって入力された測定データと、近似値算出部30によって算出された近似値とを組み合わせたデータを用いて、ビニングモデルを作成する。   As shown in FIG. 1, the binning model creation unit 60 is connected to a measurement point setting unit 50. In addition, the binning model creation unit 60 creates a binning model using data obtained by combining the measurement data input by the input unit 10 and the approximate value calculated by the approximate value calculation unit 30.

図1に示すように、推定値算出部70は、測定点設定部50およびモデルパラメータセット生成部80に接続されている。また、推定値算出部70は、電気特性の推定値を算出するための式6乃至式9に示す近似関数式(後述する)に基づいて推定値を算出して、測定点設定部50によって設定された測定点の電気特性を推定値に置換するようになっている。   As shown in FIG. 1, the estimated value calculation unit 70 is connected to a measurement point setting unit 50 and a model parameter set generation unit 80. Further, the estimated value calculation unit 70 calculates an estimated value based on approximate function expressions (described later) shown in Expressions 6 to 9 for calculating the estimated value of the electrical characteristics, and is set by the measurement point setting unit 50. The electrical characteristics of the measured points are replaced with estimated values.

図1に示すように、モデルパラメータセット生成部80は、近似値算出部30、推定値算出部70、および出力部90に接続されている。また、モデルパラメータセット生成部80は、近似値算出部30によって算出された近似値または推定値算出部70によって算出された推定値を含むモデルパラメータセットを生成するようになっている。   As shown in FIG. 1, the model parameter set generation unit 80 is connected to the approximate value calculation unit 30, the estimated value calculation unit 70, and the output unit 90. The model parameter set generating unit 80 generates a model parameter set including the approximate value calculated by the approximate value calculating unit 30 or the estimated value calculated by the estimated value calculating unit 70.

図1に示すように、出力部90は、モデルパラメータセット生成部80および出力装置(図示せず)に接続されている。また、出力部90は、モデルパラメータセット生成部80によって生成されたモデルパラメータセットを出力装置(図示せず)に出力するようになっている。   As shown in FIG. 1, the output unit 90 is connected to a model parameter set generation unit 80 and an output device (not shown). The output unit 90 outputs the model parameter set generated by the model parameter set generation unit 80 to an output device (not shown).

次に、本発明の実施例1に係るモデルパラメータ抽出装置の処理について図3乃至図9を参照して説明する。図3は、本発明の実施例1に係るモデルパラメータ抽出処理の処理手順を示すフローチャートである。図4は、図2の測定データの概略を示す概略図である。図5は、近似関数式作成処理(図3のS302)の処理手順を示すフローチャートである。図6は、第1近似関数式の係数値を示す概略図である。図7は、第2近似関数式の係数値を示す概略図である。図8は、推定値算出処理(図3のS305)の処理手順を示すフローチャートである。図9は、推定値算出処理(図3のS305)の概略を示す概略図である。   Next, processing of the model parameter extracting apparatus according to the first embodiment of the present invention will be described with reference to FIGS. FIG. 3 is a flowchart illustrating a processing procedure of model parameter extraction processing according to the first embodiment of the present invention. FIG. 4 is a schematic diagram showing an outline of the measurement data of FIG. FIG. 5 is a flowchart showing a processing procedure of the approximate function formula creation processing (S302 in FIG. 3). FIG. 6 is a schematic diagram illustrating coefficient values of the first approximate function equation. FIG. 7 is a schematic diagram illustrating the coefficient values of the second approximate function equation. FIG. 8 is a flowchart showing the processing procedure of the estimated value calculation processing (S305 in FIG. 3). FIG. 9 is a schematic diagram showing an outline of the estimated value calculation process (S305 in FIG. 3).

はじめに、図3に示すように、測定データ入力工程(S301)が行われる。測定データ入力工程(S301)では、入力部10が、ユーザが入力装置(図示せず)を用いて入力した測定点の電気特性を示す測定データを入力する。たとえば、図4に示すように、測定データは、閾値電圧Vthの測定値および第1近似関数式の値である。 First, as shown in FIG. 3, a measurement data input step (S301) is performed. In the measurement data input step (S301), the input unit 10 inputs measurement data indicating the electrical characteristics of the measurement points input by the user using an input device (not shown). For example, as shown in FIG. 4, the measurement data is a measurement value of the threshold voltage Vth and a value of the first approximate function equation.

次に、図3に示すように、近似関数式作成処理(S302)が行われる。近似関数式作成処理(S302)では、はじめに、図5に示すように、第1近似関数式作成工程(S501)が行われる。第1近似関数式作成工程(S501)では、第1近似関数式作成部201が、測定データ入力工程(図3のS301)において入力された測定データが充足しているゲート幅W(たとえば、全てのゲート長LについてのTEGが存在しているゲート幅W)についての第1近似関数式(式6および式7を参照)を作成し、各係数A乃至AおよびB乃至Bの係数値を算出する。たとえば、図4の例では、第1近似関数式作成部201は、ゲート幅W(i=0)についての閾値電圧Vthをゲート長Lの関数として近似する第1近似関数式を作成し、図6に示すような係数A乃至Aの係数値を算出する。図6は、n=5のときの第1近似関数式の各係数A乃至Aの係数値を示している。

Figure 2010061190
Figure 2010061190
Next, as shown in FIG. 3, an approximate function formula creation process (S302) is performed. In the approximate function formula creation process (S302), first, as shown in FIG. 5, a first approximate function formula creation step (S501) is performed. In the first approximate function formula creation step (S501), the first approximate function formula creation unit 201 determines that the gate width W (for example, all of the measurement data input in the measurement data input step (S301 in FIG. 3) is satisfied. A first approximate function expression (see Expression 6 and Expression 7) for a gate width W i in which a TEG for the gate length L of the current exists is created, and each coefficient A 0 to A 5 and B 0 to B 5 is created. The coefficient value of is calculated. For example, in the example of FIG. 4, the first approximate function formula creating unit 201 creates a first approximate function formula that approximates the threshold voltage V th for the gate width W i (i = 0) as a function of the gate length L. Then, the coefficient values of the coefficients A 0 to A 5 as shown in FIG. 6 are calculated. FIG. 6 shows the coefficient values of the coefficients A 0 to A 5 of the first approximate function equation when n = 5.
Figure 2010061190
Figure 2010061190

次に、図5に示すように、第2近似関数式作成工程(S502)が行われる。第2近似関数式作成工程(S502)では、第2近似関数式作成部202が、測定データ入力工程(図3のS301)において入力された測定データが不足しているゲート幅W(たとえば、少なくとも1つのゲート長LについてのTEGが存在していないゲート幅W)についての第2近似関数式(式8および式9を参照)を作成し、各係数A乃至AおよびB乃至Bの係数値を算出する。このとき、第2近似関数式作成部202は、第1近似関数式作成工程(S501)において第1近似関数式が作成されたゲート幅Wのうち、ゲート幅Wに最も近いものを選択し、ゲート幅Wについての第2近似関数式のゲート長Lの数が、ゲート幅Wについての第1近似関数式のゲート長Lの数に合うように、ゲート幅WについてTEGが存在している数だけ第2近似関数式の係数を次数の低いものから修正して、各係数A乃至AおよびB乃至Bの係数値を算出する。たとえば、図4の例では、図7に示すように、第2近似関数式作成部202は、ゲート幅W(j=1乃至3)についての第2近似関数式を作成するときには、ゲート幅W(i=0)を選択し、ゲート幅W(j=1)については1個の係数(係数A)を修正し、ゲート幅W(j=2)については2個の係数(係数AおよびA)を修正し、ゲート幅W(j=3)については3個の係数(係数A乃至A)を修正する。

Figure 2010061190
Figure 2010061190
Next, as shown in FIG. 5, a second approximate function formula creation step (S502) is performed. In the second approximate function formula creating step (S502), the second approximate function formula creating unit 202 uses the gate width W (for example, at least the measurement data input in the measurement data input step (S301 in FIG. 3) is insufficient. create second approximate function expression for the gate width W j) the TEG is not present for one gate length L (see equation 8 and equation 9), the coefficients a 0 to a 5 and B 0 to B A coefficient value of 5 is calculated. At this time, the second approximate function expression creation unit 202 of the gate width W i of the first approximate function expression is created in the first approximate function expression generation step (S501), selects the one closest to the gate width W j and, the number of gate length L of the second approximate function expression for the gate width W j is to match the number of gate length L of the first approximate function expression for the gate width W i, TEG is the gate width W j The coefficients of the second approximate function formula are corrected from the lower order by the existing number, and the coefficient values of the coefficients A 0 to A 5 and B 0 to B 5 are calculated. For example, in the example of FIG. 4, as shown in FIG. 7, the second approximate function formula creating unit 202 creates the second approximate function formula for the gate width W j (j = 1 to 3). Select W i (i = 0), modify one coefficient (coefficient A 0 ) for gate width W j (j = 1), and two coefficients for gate width W j (j = 2) (Coefficients A 0 and A 1 ) are corrected, and three coefficients (Coefficients A 0 to A 2 ) are corrected for the gate width W j (j = 3).
Figure 2010061190
Figure 2010061190

図5に示すように、本発明の実施例1に係る近似関数式作成処理(図3のS302)は、第2近似関数式作成工程(S502)の後に終了する。   As shown in FIG. 5, the approximate function formula creation process (S302 in FIG. 3) according to Embodiment 1 of the present invention ends after the second approximate function formula creation step (S502).

次に、図3に示すように、近似値算出工程(S303)が行われる。近似値算出工程(S303)では、近似値算出部30が、第1および第2近似関数式作成工程(図5のS501,S502)において作成された第1および第2近似関数式に基づいて、測定点の電気特性の近似値を算出する。その結果、図4の例では、4個のゲート幅W(i=0乃至3)のそれぞれについて、13個のゲート長L(=0.04乃至10)の電気特性(すなわち、52個の電気特性)が得られる。 Next, as shown in FIG. 3, an approximate value calculation step (S303) is performed. In the approximate value calculating step (S303), the approximate value calculating unit 30 is based on the first and second approximate function formulas created in the first and second approximate function formula creating steps (S501 and S502 in FIG. 5). Approximate the electrical characteristics at the measurement point. As a result, in the example of FIG. 4, for each of the four gate widths W i (i = 0 to 3), the electrical characteristics of 13 gate lengths L (= 0.04 to 10) (that is, 52 gates) Electrical characteristics).

次に、図3に示すように、推定値を算出する場合には(S304−YES)、推定値算出処理(S305)が行われる。推定値算出処理(S305)では、はじめに、図8に示すように、ビニングモデル作成工程(S801)が行われる。ビニングモデル作成工程(S801)では、ビニングモデル作成部60が、測定データ入力工程(図3のS301)において入力された測定データと、近似値算出部30によって算出された近似値とを組み合わせたデータを用いてビニングモデルを作成する。このとき、ビニングモデル作成部60は、MOSFETのL−W平面を格子状に分割し、各格子点の中から任意の4点を選択し、それらの頂点によって囲まれる局所領域をビニングモデルとして作成する。また、ビニングモデル作成部60は、このようなビニングモデルを複数個作成する。   Next, as shown in FIG. 3, when an estimated value is calculated (S304—YES), an estimated value calculating process (S305) is performed. In the estimated value calculation process (S305), first, as shown in FIG. 8, a binning model creation step (S801) is performed. In the binning model creation step (S801), the binning model creation unit 60 combines the measurement data input in the measurement data input step (S301 in FIG. 3) and the approximate value calculated by the approximate value calculation unit 30. Create a binning model using. At this time, the binning model creation unit 60 divides the LW plane of the MOSFET into a lattice shape, selects any four points from each lattice point, and creates a local region surrounded by the vertices as a binning model. To do. Further, the binning model creation unit 60 creates a plurality of such binning models.

次に、図8に示すように、読取工程(S802)が行われる。読取工程(S802)では、読取部40が、第1および第2近似関数式作成工程(図5のS501,S502)において作成された第1および第2近似関数式のゲート長Lの電気特性をゲート長Lの最小値Lminから等間隔ΔLごとに読み取る(図9(a)を参照)。このとき、読取部40は、読み取られたゲート長Lmin+kΔL(k=0,1,2,・・・)の電気特性とビニングモデル作成工程(S801)において作成されたビニングモデルの電気特性とを比較する。 Next, as shown in FIG. 8, a reading step (S802) is performed. In the reading step (S802), the reading unit 40 obtains the electrical characteristics of the gate length L of the first and second approximate function formulas created in the first and second approximate function formula creation steps (S501 and S502 in FIG. 5). Read from the minimum value L min of the gate length L at regular intervals ΔL (see FIG. 9A). At this time, the reading unit 40 reads the electrical characteristics of the read gate length L min + kΔL (k = 0, 1, 2,...) And the electrical characteristics of the binning model created in the binning model creation step (S801). Compare

次に、図8に示すように、読取工程(S802)において読み取られたゲート長Lmin+kΔL(k=0,1,2,・・・)の電気特性とビニングモデル作成工程(S801)において作成途上のビニングモデルのゲート長Lmin+kΔL(k=0,1,2,・・・)の電気特性との比較結果(以下、「電気特性誤差」という)が所定の許容値ε以上である場合には(S803−YES)、保持工程(S804)が行われる。保持工程(S804)では、読取部40が、ゲート長Lmin+kΔL(k=0,1,2,・・・)の値およびその電気特性を保持する。 Next, as shown in FIG. 8, the electrical characteristics of the gate length L min + kΔL (k = 0, 1, 2,...) Read in the reading step (S802) and the binning model generation step (S801). When a comparison result (hereinafter referred to as “electric characteristic error”) of the gate length L min + kΔL (k = 0, 1, 2,...) Of the binning model on the way is equal to or greater than a predetermined allowable value ε (S803-YES), a holding process (S804) is performed. In the holding step (S804), the reading unit 40 holds the value of the gate length L min + kΔL (k = 0, 1, 2,...) And its electrical characteristics.

図8に示すように、読取工程(S802)乃至保持工程(S804)は、全てのゲート長Lについての電気特性と第1および第2近似関数式のゲート長Lについての電気特性との比較が完了するまで繰り返し行われる(S805−NO)。   As shown in FIG. 8, in the reading process (S802) to the holding process (S804), the electrical characteristics for all the gate lengths L and the electrical characteristics for the gate lengths L of the first and second approximate function formulas are compared. It is repeated until completion (S805-NO).

次に、図8に示すように、全てのゲート長Lについての電気特性と第1および第2近似関数式のゲート長Lについての電気特性との比較が完了した場合には(S805−YES)、図8に示すように、測定点設定工程(S806)が行われる。測定点設定工程(S806)では、測定点設定部50が、保持工程(S804)において保持されたゲート長Lmin+kΔL(k=0,1,2,・・・)に測定点を設定する。このとき、測定点設定部50は、そのゲート長Lmin+kΔL(k=0,1,2,・・・)より小さいものの中で最大となる測定点のゲート長Lとゲート長Lmin+kΔL(k=0,1,2,・・・)より大きいものの中で最小となる測定点のゲート長Lm+1に着目し、ゲート長Lとゲート長Lm+1との間に複数の候補が存在する場合には、任意の1つ(たとえば、電気特性誤差が最大となるもの)に測定点を設定する(図9(b)を参照)。 Next, as shown in FIG. 8, when the comparison between the electrical characteristics for all the gate lengths L and the electrical characteristics for the gate lengths L of the first and second approximate function formulas is completed (YES in S805). As shown in FIG. 8, a measurement point setting step (S806) is performed. In the measurement point setting step (S806), the measurement point setting unit 50 sets the measurement point to the gate length L min + kΔL (k = 0, 1, 2,...) Held in the holding step (S804). At this time, the measurement point setting unit 50 has the gate length L m and the gate length L min + kΔL at the maximum measurement point among those smaller than the gate length L min + kΔL (k = 0, 1, 2,...). Paying attention to the gate length L m + 1 of the smallest measurement point among those larger than (k = 0, 1, 2,...), There are a plurality of candidates between the gate length L m and the gate length L m + 1. In this case, the measurement point is set to an arbitrary one (for example, the electric characteristic error is maximized) (see FIG. 9B).

次に、図8に示すように、推定値算出工程(S807)が行われる。推定算出工程(S807)では、推定値算出部70が、測定点設定工程(S806)において設定された測定点の電気特性の推定値を所定の近似関数式に基づいて算出する。   Next, as shown in FIG. 8, an estimated value calculating step (S807) is performed. In the estimation calculation step (S807), the estimated value calculation unit 70 calculates the estimated value of the electrical characteristic of the measurement point set in the measurement point setting step (S806) based on a predetermined approximate function equation.

図8に示すように、読取工程(S802)乃至推定値算出工程(S807)は、全てのゲート幅Wの設定された測定点ついて推定値が算出されるまで繰り返し行われる(S808−NO)。また、本発明の実施例1に係る推定値算出処理(図3のS305)は、全てのゲート幅Wの設定された測定点について推定値が算出された後に終了する(S808−YES)。 As shown in FIG. 8, the reading step (S802) to the estimated value calculation step (S807) is repeatedly performed until the estimated value with the set measurement points of all the gate width W i is calculated (S808-NO) . Also, the estimated value calculation processing according to the first embodiment of the present invention (S305 in FIG. 3) is terminated after the estimated value is calculated for the set measurement points of all the gate width W i (S808-YES).

また、図3に示すように、推定値を算出しない場合(S304−NO)または推定値算出処理(S305)の後に、モデルパラメータセット作成工程(S306)が行われる。モデルパラメータセット作成工程(S306)では、モデルパラメータセット生成部80が、測定データ入力工程(S301)において入力された測定データのうちI−V特性カーブを用いて、L−W平面全体を1つのモデルで表現するグローバルモデルを作成し、式1乃至式4を解くことによって、各ゲート長Lまたは各ゲート幅Wについて図4に示す測定値または第1近似関数式の値に合うように、そのグローバルモデルの中の一部のパラメータ値を修正し、近似値算出工程(S303)において算出された電気特性の近似値を再現するような52個の格子点のモデルパラメータセットを作成して、その52個の格子点の中の任意の4個の格子点のモデルパラメータセットの値に基づいて、それらの4個の格子点によって形成されるビニングモデルを作成する。   Further, as shown in FIG. 3, when the estimated value is not calculated (S304-NO) or after the estimated value calculating process (S305), a model parameter set creating step (S306) is performed. In the model parameter set creation step (S306), the model parameter set generation unit 80 uses the IV characteristic curve of the measurement data input in the measurement data input step (S301) to make one entire LW plane. By creating a global model expressed by the model and solving Equations 1 to 4, each gate length L or each gate width W is adjusted so that it matches the measured value shown in FIG. 4 or the value of the first approximate function equation. A part of parameter values in the global model is corrected to create a model parameter set of 52 lattice points that reproduces the approximate value of the electrical characteristics calculated in the approximate value calculation step (S303). Based on the value of the model parameter set of any four lattice points among the 52 lattice points, the binaries formed by these four lattice points To create a Gumoderu.

次に、図3に示すように、出力工程(S307)が行われる。出力工程(S307)では、出力部90が、モデルパラメータセット作成工程(S306)において作成されたモデルパラメータセット群(以下、「ビニングモデルパラメータセット」という)を出力装置(図示せず)に出力する。   Next, as shown in FIG. 3, an output step (S307) is performed. In the output step (S307), the output unit 90 outputs the model parameter set group (hereinafter referred to as “binning model parameter set”) created in the model parameter set creation step (S306) to an output device (not shown). .

図3に示すように、本発明の実施例に係るモデルパラメータ抽出処理は、出力工程(S307)の後に終了する。   As shown in FIG. 3, the model parameter extraction process according to the embodiment of the present invention ends after the output step (S307).

なお、本発明の実施例1では、電気特性とゲート長Lの関係を与える近似関数式として、ゲート長Lを最大ゲート長LMAXで割って規格化し、対数を取り、それをn次多項式で表現した。ここで、n次多項式の各項の係数は、近似対象の電気特性を近似するための調整パラメータでありである。また、最大ゲート長LMAXは、使用が想定されるMOSFETのゲート長Lの最大値であり、図2の例では、最大ゲート長LMAX=10である。 In the first embodiment of the present invention, as an approximate function that gives the relationship between the electrical characteristics and the gate length L, the gate length L is divided by the maximum gate length L MAX and normalized, and the logarithm is taken. Expressed. Here, the coefficient of each term of the nth order polynomial is an adjustment parameter for approximating the electrical characteristics to be approximated. The maximum gate length L MAX is the maximum value of the gate length L of the MOSFET assumed to be used. In the example of FIG. 2, the maximum gate length L MAX = 10.

また、対数を取った理由は、ゲート長Lが小さいときの電気特性の急激な変化の再現性を高めるためである。   The reason why the logarithm is taken is to improve the reproducibility of a sudden change in electrical characteristics when the gate length L is small.

また、ゲート長Lを最大ゲート長LMAXで割って規格化した理由は、対数の中を1以下とすることによって、対数の評価値の符号を負に揃えるためである。 The reason why the gate length L is normalized by dividing it by the maximum gate length L MAX is to make the sign of the logarithmic evaluation value negative by setting the logarithm to 1 or less.

また、n次多項式を使用して各項の係数を調整する理由は、関数が表現する電気特性の自由度を高め、電気特性への合せ込み精度を高めるためである。   The reason for adjusting the coefficient of each term using an nth-order polynomial is to increase the degree of freedom of the electrical characteristics expressed by the function and increase the accuracy of fitting to the electrical characteristics.

使用する近似関数式は、電気特性を精度よく近似できるものであれば何でも良く、近似する電気特性によって適切なものが選択されることが好ましい。たとえば、他の近似関数式は、式6および式7に代えて、式10および式11、または式12および式13を用いても良い。このように、近似関数式として、ゲート長Lの関数で与えた中間変数(たとえば、log(L)を0乗、1乗、2乗、…、n乗して、それぞれに調整パラメータを掛けたものを加算した形式の関数か、または、そのような関数を内部に含む関数を用いることで、近似の精度を向上させることができる。

Figure 2010061190
The approximate function equation to be used is not particularly limited as long as the electrical characteristics can be approximated with high accuracy, and it is preferable to select an appropriate one depending on the electrical characteristics to be approximated. For example, instead of Equation 6 and Equation 7, Equation 10 and Equation 11 or Equation 12 and Equation 13 may be used as other approximate function equations. Thus, as an approximate function expression, an intermediate variable (for example, log (L)) given as a function of the gate length L is raised to the 0th power, the first power, the second power,. The accuracy of approximation can be improved by using a function in the form of adding the functions or a function including such a function inside.
Figure 2010061190

なお、本発明の実施例1では、ゲート長Lの全範囲を1つの関数で表したが、ゲート長Lの領域ごとに異なる関数で表しても良い。   In the first embodiment of the present invention, the entire range of the gate length L is represented by one function, but may be represented by a different function for each region of the gate length L.

また、本発明の実施例1では、第1および第2近似関数式作成部201は、入力部10によって入力された測定データのうちTEGが存在しているゲート長Lの電気特性をゲート幅Wの関数として禁じするn次の多項式を作成しても良い。この場合には、読取部40は、ゲート幅W方向の所定間隔毎に電気特性を読み取る。また、測定点設定部50は、読取部40によって読み取られた電気特性が許容範囲外を示すゲート幅Wに測定点を設定する。   Further, in the first embodiment of the present invention, the first and second approximate function equation creation unit 201 calculates the electrical characteristics of the gate length L where the TEG is present in the measurement data input by the input unit 10 as the gate width W. It is also possible to create an nth-order polynomial that is prohibited as a function of. In this case, the reading unit 40 reads the electrical characteristics at predetermined intervals in the gate width W direction. In addition, the measurement point setting unit 50 sets the measurement point to the gate width W where the electrical characteristic read by the reading unit 40 is outside the allowable range.

また、本発明の実施例1では、読取部40は、第1および第2近似関数式の一方について、ゲート長Lの電気特性をゲート長Lの最小値Lminから等間隔ΔLごとに読み取っても良い。 In Embodiment 1 of the present invention, the reading unit 40 reads the electrical characteristics of the gate length L from the minimum value L min of the gate length L at regular intervals ΔL for one of the first and second approximate function equations. Also good.

本発明の実施例1によれば、図3に示すように、測定データが不足している場合であっても、近似値や推定値を算出するので、高精度なモデルパラメータセットを抽出することができる。   According to the first embodiment of the present invention, as shown in FIG. 3, even if measurement data is insufficient, approximate values and estimated values are calculated, so that a highly accurate model parameter set can be extracted. Can do.

本発明の実施例1に係るモデルパラメータ抽出装置1の構成を示すブロック図である。It is a block diagram which shows the structure of the model parameter extraction apparatus 1 which concerns on Example 1 of this invention. 本発明の実施例1に係る測定データの概略を示す概略図である。It is the schematic which shows the outline of the measurement data which concern on Example 1 of this invention. 本発明の実施例1に係るモデルパラメータ抽出処理の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of the model parameter extraction process which concerns on Example 1 of this invention. 図2の測定データの概略を示す概略図である。It is the schematic which shows the outline of the measurement data of FIG. 近似関数式作成処理(図3のS302)の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of an approximate function formula creation process (S302 of FIG. 3). 第1近似関数式の係数値を示す概略図である。It is the schematic which shows the coefficient value of a 1st approximation function type | formula. 第2近似関数式の係数値を示す概略図である。It is the schematic which shows the coefficient value of a 2nd approximation function type | formula. 推定値算出処理(図3のS305)の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of an estimated value calculation process (S305 of FIG. 3). 推定値算出処理(図3のS305)の概略を示す概略図である。It is the schematic which shows the outline of an estimated value calculation process (S305 of FIG. 3).

符号の説明Explanation of symbols

1 モデルパラメータ抽出装置
10 入力部
20 近似関数式作成部
201 第1近似関数式作成部
201 第2近似関数式作成部
30 近似値算出部
40 読取部
50 測定点追加部
60 ビニングモデル作成部
70 推定値算出部
80 モデルパラメータセット生成部
90 出力部
DESCRIPTION OF SYMBOLS 1 Model parameter extraction apparatus 10 Input part 20 Approximate functional formula creation part 201 1st approximate function formula creation part 201 2nd approximate function formula creation part 30 Approximate value calculation part 40 Reading part 50 Measurement point addition part 60 Binning model creation part 70 Estimation Value calculation unit 80 Model parameter set generation unit 90 Output unit

Claims (5)

トランジスタのゲート長方向およびゲート幅方向に対して格子状に配置された測定点によって囲まれるシミュレーションモデルのモデルパラメータを抽出するモデルパラメータ抽出装置であって、
前記測定点の電気特性を示す測定データを入力する入力部と、
前記入力部によって入力された測定データが充足している測定点について、前記入力部によって入力された測定データに基づいて、前記測定点の電気特性を近似する第1近似関数式を作成する第1近似関数式作成部と、
前記入力部によって入力された測定データが不足している測定点について、前記第1近似関数式作成部によって作成された第1近似関数式を修正して、前記測定点の電気特性を近似する第2近似関数式を作成する第2近似関数式作成部と、
前記第1および第2近似関数式作成部によって作成された第1および第2近似関数式に基づいて、前記測定点の電気特性の近似値を算出する近似値算出部と、
前記入力部によって入力された測定データと、前記近似値算出部によって算出された近似値とを組み合わせたデータから、モデルパラメータセットを生成するモデルパラメータセット生成部と、
前記モデルパラメータセット生成部によって生成されたモデルパラメータセットを出力する出力部と、を備えることを特徴とするモデルパラメータ抽出装置。
A model parameter extraction device for extracting model parameters of a simulation model surrounded by measurement points arranged in a grid pattern in a gate length direction and a gate width direction of a transistor,
An input unit for inputting measurement data indicating the electrical characteristics of the measurement point;
A first approximation function formula that approximates the electrical characteristics of the measurement point is generated based on the measurement data input by the input unit for the measurement point that is satisfied by the measurement data input by the input unit. An approximate function formula creation unit;
For a measurement point for which measurement data input by the input unit is insufficient, the first approximate function formula created by the first approximate function formula creation unit is modified to approximate the electrical characteristics of the measurement point. A second approximate function formula creating unit for creating a two approximate function formula;
An approximate value calculation unit that calculates an approximate value of the electrical characteristics of the measurement point based on the first and second approximate function formulas created by the first and second approximate function formula creation units;
A model parameter set generation unit that generates a model parameter set from data obtained by combining the measurement data input by the input unit and the approximate value calculated by the approximate value calculation unit;
An output unit that outputs a model parameter set generated by the model parameter set generation unit.
前記第1および第2近似関数式作成部によって作成された第1および第2近似関数式の少なくとも一方に基づいて、前記ゲート幅または前記ゲート長方向の電気特性を所定間隔毎に読み取る読取部と、
前記読取部によって読み取られた電気特性と作成途上のビニングモデルの電気特性との誤差が許容範囲外である場合に、その電気特性に対応する前記ゲート幅または前記ゲート長に前記測定点を設定する測定点設定部と、
前記測定点設定部によって設定された測定点の電気特性の推定値を算出する推定値算出部と、をさらに備え、
前記モデルパラメータセット生成部は、前記推定値算出部によって算出された推定値を含むモデルパラメータセットを生成する請求項1に記載のモデルパラメータ抽出装置。
A reading unit that reads the electrical characteristics in the gate width or the gate length direction at predetermined intervals based on at least one of the first and second approximate function formulas created by the first and second approximate function formula creation units; ,
When the error between the electrical characteristic read by the reading unit and the electrical characteristic of the binning model being created is outside the allowable range, the measurement point is set to the gate width or the gate length corresponding to the electrical characteristic. A measurement point setting section;
An estimated value calculation unit that calculates an estimated value of the electrical characteristics of the measurement point set by the measurement point setting unit,
The model parameter extraction device according to claim 1, wherein the model parameter set generation unit generates a model parameter set including the estimated value calculated by the estimated value calculation unit.
前記読取部は、前記ゲート幅または前記ゲート長方向の電気特性を等間隔毎に読み取る請求項2に記載のモデルパラメータ抽出装置。   The model parameter extracting apparatus according to claim 2, wherein the reading unit reads the electrical characteristics in the gate width direction or the gate length direction at equal intervals. トランジスタのゲート長方向およびゲート幅方向に対して格子状に配置された測定点によって囲まれるシミュレーションモデルのモデルパラメータを抽出するモデルパラメータ抽出方法であって、
前記測定点の電気特性を示す測定データを入力し、
前記測定データが充足している測定点について、前記測定データに基づいて、前記測定点の電気特性を近似する第1近似関数式を作成し、
前記測定データが不足している測定点について、前記第1近似関数式を修正して、前記測定点の電気特性を近似する第2近似関数式を作成し、
前記第1および第2近似関数式に基づいて、前記測定点の電気特性の近似値を算出し、
前記近似値を含むモデルパラメータセットを生成し、
前記モデルパラメータセットを出力することを特徴とするモデルパラメータ抽出方法。
A model parameter extraction method for extracting a model parameter of a simulation model surrounded by measurement points arranged in a grid pattern in a gate length direction and a gate width direction of a transistor,
Input measurement data indicating the electrical characteristics of the measurement point,
For the measurement points satisfied by the measurement data, based on the measurement data, create a first approximate function equation that approximates the electrical characteristics of the measurement points,
For the measurement point for which the measurement data is insufficient, the first approximate function formula is corrected to create a second approximate function formula that approximates the electrical characteristics of the measurement point,
Based on the first and second approximation function equations, an approximate value of the electrical characteristics of the measurement point is calculated,
Generating a model parameter set including the approximation,
A model parameter extracting method, wherein the model parameter set is output.
前記第1および第2近似関数式に基づいて、前記ゲート幅または前記ゲート長方向の電気特性を所定間隔毎に読み取り、
前記読み取られた電気特性と作成途上のビニングモデルの電気特性との誤差が許容範囲外である場合に、その電気特性に対応する前記ゲート幅または前記ゲート長に前記測定点を設定し、
前記設定された測定点の電気特性の推定値を算出し、
前記算出された推定値を含むモデルパラメータセットを生成する請求項4に記載のモデルパラメータ抽出方法。
Based on the first and second approximate function equations, the electrical characteristics in the gate width or the gate length direction are read at predetermined intervals,
When the error between the read electrical characteristics and the electrical characteristics of the binning model being created is outside the allowable range, the measurement point is set to the gate width or the gate length corresponding to the electrical characteristics,
Calculate an estimated value of the electrical characteristics of the set measurement point,
The model parameter extraction method according to claim 4, wherein a model parameter set including the calculated estimated value is generated.
JP2008223247A 2008-09-01 2008-09-01 Model parameter extraction device and model parameter extraction method Withdrawn JP2010061190A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008223247A JP2010061190A (en) 2008-09-01 2008-09-01 Model parameter extraction device and model parameter extraction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008223247A JP2010061190A (en) 2008-09-01 2008-09-01 Model parameter extraction device and model parameter extraction method

Publications (1)

Publication Number Publication Date
JP2010061190A true JP2010061190A (en) 2010-03-18

Family

ID=42187951

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008223247A Withdrawn JP2010061190A (en) 2008-09-01 2008-09-01 Model parameter extraction device and model parameter extraction method

Country Status (1)

Country Link
JP (1) JP2010061190A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010067830A (en) * 2008-09-11 2010-03-25 Fujitsu Ltd Model parameter abstraction program
CN105373634A (en) * 2014-08-18 2016-03-02 中国科学院微电子研究所 A method and system for improving the extraction accuracy of specific device model parameters

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010067830A (en) * 2008-09-11 2010-03-25 Fujitsu Ltd Model parameter abstraction program
CN105373634A (en) * 2014-08-18 2016-03-02 中国科学院微电子研究所 A method and system for improving the extraction accuracy of specific device model parameters
CN105373634B (en) * 2014-08-18 2018-07-03 中国科学院微电子研究所 A method and system for improving the extraction accuracy of specific device model parameters

Similar Documents

Publication Publication Date Title
JP4882902B2 (en) Simulation method and program
JP5148354B2 (en) Integrated circuit modeling method and integrated circuit
JP5390154B2 (en) Circuit simulation apparatus, circuit simulation method, and program
US8954908B1 (en) Fast monte carlo statistical analysis using threshold voltage modeling
JP4874207B2 (en) Circuit simulation method, circuit simulation apparatus, and program
US8151240B2 (en) Effective gate length circuit modeling based on concurrent length and mobility analysis
KR19980064044A (en) Hot carrier degradation estimation method
US20060142987A1 (en) Circuit simulation method and circuit simulation apparatus
JP2020525872A (en) Influenza prediction model generation method, device, and computer-readable storage medium
JP2009026829A (en) Designing method for semiconductor integrated circuit, and mask data generating program
JP2012048706A (en) Constructing mapping between model parameters and electrical parameters
JP2008027940A (en) Semiconductor integrated circuit design method and circuit simulation method
US9646124B2 (en) Modeling transistor performance considering non-uniform local layout effects
CN113408231A (en) Method for establishing degradation simulation model
JPWO2012081158A1 (en) Circuit simulation method and semiconductor integrated circuit
US20190219625A1 (en) Predicting noise propagation in circuits
US8082532B1 (en) Placing complex function blocks on a programmable integrated circuit
JP2010061190A (en) Model parameter extraction device and model parameter extraction method
JP2008053692A (en) Simulation model of BT degradation of transistor and simulation modeling method
CN111579961B (en) Method, apparatus and computer-readable storage medium for determining electrical characteristics of a chip
TW200417887A (en) Circuit simulation for a circuit
US20160259871A1 (en) Model generation method and information processing apparatus
JP5560700B2 (en) Design support apparatus, design support method, and design support program
JP5227141B2 (en) Model parameter extracting apparatus and model parameter extracting program for semiconductor device model
JP2013232135A (en) Layout design device, layout design method, and program

Legal Events

Date Code Title Description
A300 Application deemed to be withdrawn because no request for examination was validly filed

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20111101