JP2012154868A - System for analyzing ion behavior in insulating film - Google Patents

System for analyzing ion behavior in insulating film Download PDF

Info

Publication number
JP2012154868A
JP2012154868A JP2011016016A JP2011016016A JP2012154868A JP 2012154868 A JP2012154868 A JP 2012154868A JP 2011016016 A JP2011016016 A JP 2011016016A JP 2011016016 A JP2011016016 A JP 2011016016A JP 2012154868 A JP2012154868 A JP 2012154868A
Authority
JP
Japan
Prior art keywords
insulating film
ion
ion density
density distribution
analysis system
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.)
Pending
Application number
JP2011016016A
Other languages
Japanese (ja)
Inventor
Kinya Kobayashi
金也 小林
Junpei Kusukawa
順平 楠川
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2011016016A priority Critical patent/JP2012154868A/en
Priority to PCT/JP2012/051602 priority patent/WO2012102321A1/en
Publication of JP2012154868A publication Critical patent/JP2012154868A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/22Details, e.g. general constructional or apparatus details
    • G01N29/24Probes
    • G01N29/2431Probes using other means for acoustic excitation, e.g. heat, microwaves, electron beams
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/02863Electric or magnetic parameters

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Testing Or Measuring Of Semiconductors Or The Like (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a system for analyzing ion behavior in an insulating film capable of predicting the lifetime of an insulating film by performing nondestructive measurement.SOLUTION: The system for analyzing ion behavior in an insulating film includes a pulse voltage generator, a rolling element, an arithmetic unit and a display unit. The system obtains ion density distribution in the insulating film by detecting ultrasonic waves which occur in the insulating film by pulse voltage applied by the pulse voltage generator, by the rolling element. The arithmetic unit obtains the ion density distribution in the insulating film at least at different times based on a detection result of ultrasonic waves, performs advective diffusion calculation based on the ion density distribution at least at different times to obtain the ion density distribution in the insulating film at arbitrary time. The display unit displays the ion density distribution in the insulating film at arbitrary time.

Description

本発明は、絶縁膜内のイオンの挙動を解析するシステムに関する。   The present invention relates to a system for analyzing behavior of ions in an insulating film.

絶縁膜は、装置の作動中に徐々に劣化していき、寿命に達すると絶縁破壊を起こす。モータ、インバータ、または変圧器等の電気機器に用いられる絶縁膜では、高電界化に伴い絶縁膜内のイオンマイグレーションによる絶縁破壊が起こりやすいので、絶縁破壊の抑制や絶縁膜の寿命の予測が必須である。   The insulating film gradually deteriorates during the operation of the apparatus, and causes dielectric breakdown when the life is reached. Insulating films used in electrical equipment such as motors, inverters, and transformers are prone to dielectric breakdown due to ion migration in the insulating film as the electric field increases, so it is essential to suppress dielectric breakdown and predict the life of the insulating film. It is.

絶縁膜の寿命を予測するには、絶縁膜内のイオン密度分布を求める必要があるが、絶縁膜内のイオンの挙動を評価するのは困難である。このため、従来は、絶縁マージンを大きくすることで絶縁破壊を抑制しており、モータ、インバータ、または変圧器等の電気機器の小型化と信頼性の両立が困難となっている。   In order to predict the lifetime of the insulating film, it is necessary to obtain the ion density distribution in the insulating film, but it is difficult to evaluate the behavior of ions in the insulating film. For this reason, conventionally, insulation breakdown is suppressed by increasing the insulation margin, and it is difficult to achieve both downsizing and reliability of electrical equipment such as a motor, an inverter, or a transformer.

特許文献1には、絶縁膜を破壊させて、絶縁膜の寿命を推定する方法が開示されている。特許文献2には、絶縁材料を破壊せずにイオンマイグレーションを検出することができるイオンマイグレーション発生評価方法が開示されている。非特許文献1には、非破壊の計測を行って絶縁膜内のイオン分布を計測する空間電荷測定装置が開示されている。   Patent Document 1 discloses a method for estimating the lifetime of an insulating film by destroying the insulating film. Patent Document 2 discloses an ion migration occurrence evaluation method capable of detecting ion migration without destroying an insulating material. Non-Patent Document 1 discloses a space charge measuring device that performs nondestructive measurement and measures the ion distribution in the insulating film.

特開平10−261681号公報Japanese Patent Laid-Open No. 10-261681 特開2002−156359号公報JP 2002-156359 A

ファイブラボ 株式会社、「PEA空間電荷測定装置」、[online]、[2010年12月6日検索]、インターネット〈URL:http://www.5lab.co.jp/products_pea.htm〉Fibravo, "PEA space charge measuring device", [online], [December 6, 2010 search], Internet <URL: http://www.5lab.co.jp/products_pea.htm>

特許文献1に記載の技術のように絶縁膜を破壊して絶縁膜の寿命を推定する方法では、破壊により露出した絶縁膜の断面が化学変化を起こして実際とは異なる状況になり、推定した寿命に誤差が生じる可能性がある。また、断面作成には研磨が必要であり、このために工数や時間がかかるというデメリットもある。このため、絶縁膜を破壊しないで絶縁膜の寿命を予測する方法が求められている。   In the method of estimating the lifetime of the insulating film by destroying the insulating film as in the technique described in Patent Document 1, the cross section of the insulating film exposed by the breakdown causes a chemical change, which is different from the actual situation. There may be an error in the service life. In addition, the cross-section needs to be polished, which has the demerit of man-hours and time. For this reason, there is a need for a method for predicting the life of an insulating film without destroying the insulating film.

特許文献2に記載されているようなイオンマイグレーション発生評価方法は、絶縁膜を破壊することなくイオンマイグレーションの発生の可能性を検出するものであるが、絶縁膜の寿命を予測するものではない。   The ion migration occurrence evaluation method described in Patent Document 2 detects the possibility of the occurrence of ion migration without destroying the insulating film, but does not predict the lifetime of the insulating film.

非特許文献1に開示されている空間電荷測定装置では、絶縁膜を破壊することなく、絶縁膜内のイオン分布を導出して表示する。さらに、求めたイオン分布から電界分布を導出して表示することが可能である。しかし、測定した時刻以外の時刻でのイオン分布を導出することはできない。したがって、現在より後の任意の時刻でのイオン分布を求めることができず、絶縁膜の寿命を予測することができない。   In the space charge measuring device disclosed in Non-Patent Document 1, the ion distribution in the insulating film is derived and displayed without destroying the insulating film. Furthermore, it is possible to derive and display the electric field distribution from the obtained ion distribution. However, an ion distribution at a time other than the measured time cannot be derived. Therefore, the ion distribution at an arbitrary time after the present time cannot be obtained, and the lifetime of the insulating film cannot be predicted.

本発明の目的は、非破壊の計測を行って絶縁膜内のイオン密度分布を求め、絶縁膜の寿命を予測することが可能な絶縁膜内イオン挙動解析システムを提供することである。   An object of the present invention is to provide an ion behavior analysis system in an insulating film that can perform nondestructive measurement to obtain an ion density distribution in the insulating film and predict the life of the insulating film.

本発明による絶縁膜内イオン挙動解析システムは、パルス電圧発生器と圧延素子と演算装置と表示装置を備え、前記パルス電圧発生器が印加したパルス電圧により絶縁膜内に発生した超音波を前記圧延素子で検出することにより、前記絶縁膜内のイオン密度分布を求める。前記演算装置は、少なくとも異なる2つの時刻における前記絶縁膜内のイオン密度分布を前記超音波の検出結果から求め、前記少なくとも異なる2つの時刻における前記イオン密度分布を基に移流拡散計算を行って、任意の時刻における前記絶縁膜内のイオン密度分布を求める。前記表示装置は、前記任意の時刻における前記絶縁膜内のイオン密度分布を表示する。   An ion behavior analysis system in an insulating film according to the present invention includes a pulse voltage generator, a rolling element, a calculation device, and a display device, and ultrasonic waves generated in the insulating film by a pulse voltage applied by the pulse voltage generator are rolled. By detecting with an element, the ion density distribution in the insulating film is obtained. The arithmetic device obtains ion density distribution in the insulating film at at least two different times from the detection result of the ultrasonic wave, performs advection diffusion calculation based on the ion density distribution at the at least two different times, An ion density distribution in the insulating film at an arbitrary time is obtained. The display device displays an ion density distribution in the insulating film at the arbitrary time.

本発明による絶縁膜内イオン挙動解析システムでは、非破壊の超音波計測を行って絶縁膜内のイオン密度分布を解析し、絶縁膜の寿命を予測することが可能である。   With the ion behavior analysis system in an insulating film according to the present invention, it is possible to predict the lifetime of the insulating film by analyzing the ion density distribution in the insulating film by performing nondestructive ultrasonic measurement.

実施例1による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。It is a figure which shows the processing flow for calculating | requiring the ion density distribution and electric field distribution inside the insulating film of a sample of the ion behavior analysis system in an insulating film by Example 1. FIG. 試料の平面図(上面図)である。It is a top view (top view) of a sample. 試料の断面図である。It is sectional drawing of a sample. 本実施例による絶縁膜内イオン挙動解析システムの構成を示す概略図である。It is the schematic which shows the structure of the ion behavior analysis system in an insulating film by a present Example. 実施例1により求めたイオン密度分布のグラフの例である。2 is an example of a graph of ion density distribution obtained by Example 1. FIG. 実施例1により求めた電界分布のグラフの例である。3 is an example of a graph of electric field distribution obtained by Example 1. 実施例2による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の寿命を予測するための処理フローを示す図である。It is a figure which shows the processing flow for predicting the lifetime of the insulating film of a sample of the ion behavioral analysis system in an insulating film by Example 2. FIG. 実施例3による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。It is a figure which shows the processing flow for calculating | requiring the ion density distribution and electric field distribution inside the insulating film of a sample of the ion behavior analysis system in an insulating film by Example 3. FIG. 実施例4による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の内部のイオン密度分布、電界分布、電流分布、及び電位分布を求めるための処理フローを示す図である。It is a figure which shows the processing flow for calculating | requiring the ion density distribution, electric field distribution, electric current distribution, and electric potential distribution inside the insulating film of a sample of the ion behavior analysis system in an insulating film by Example 4. FIG. 実施例5による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。It is a figure which shows the processing flow for calculating | requiring the ion density distribution and electric field distribution inside the insulating film of a sample of the ion behavior analysis system in an insulating film by Example 5. FIG. 実施例5と実施例6による絶縁膜内イオン挙動解析システムの、パルス電圧発生器と試料台を示す図である。It is a figure which shows the pulse voltage generator and sample stand of the ion behavior analysis system in the insulating film by Example 5 and Example 6. FIG. 実施例6による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の内部のイオン密度分布及び電界分布を求め、絶縁寿命が最も短い領域を再成型するための処理フローを示す図である。It is a figure which shows the processing flow for calculating | requiring the ion density distribution and electric field distribution inside the insulating film of a sample of the ion behavior analysis system in an insulating film by Example 6, and remolding the area | region with the shortest insulation lifetime. 実施例7による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。It is a figure which shows the processing flow for calculating | requiring the ion density distribution and electric field distribution inside the insulating film of a sample of the ion behavior analysis system in an insulating film by Example 7. FIG. 実施例8による絶縁膜内イオン挙動解析システムの、試料の絶縁膜の内部のイオン密度分布の温度依存性を求めるための処理フローを示す図である。It is a figure which shows the processing flow for calculating | requiring the temperature dependence of the ion density distribution inside the insulating film of a sample of the ion behavioral analysis system in an insulating film by Example 8. FIG.

本発明による絶縁膜内イオン挙動解析システムでは、非破壊の超音波計測を行って絶縁膜内のイオン密度分布と電界分布を解析し、絶縁膜の寿命を予測することが可能である。このため、寿命に達する前に絶縁膜を交換することができ、絶縁膜の絶縁破壊を防止し、モータ、インバータ、または変圧器等の絶縁膜を用いる電気機器の異常動作を防止することが可能である。   The ion behavior analysis system in an insulating film according to the present invention can predict the lifetime of an insulating film by performing nondestructive ultrasonic measurement to analyze the ion density distribution and electric field distribution in the insulating film. For this reason, it is possible to replace the insulation film before reaching the end of its life, prevent insulation breakdown of the insulation film, and prevent abnormal operation of electrical equipment using insulation films such as motors, inverters, or transformers It is.

対象の絶縁膜の例としては、エポキシ、酸化膜、マイカ、碍子、グラファイト、酸化アルミニューム、及びガラスが挙げられる。絶縁膜内のイオンの例としては、正イオンは銅、アルミニューム、及び銀などの金属イオンが主であり、負イオンは電子が主である。   Examples of the target insulating film include epoxy, oxide film, mica, insulator, graphite, aluminum oxide, and glass. As examples of ions in the insulating film, positive ions are mainly metal ions such as copper, aluminum, and silver, and negative ions are mainly electrons.

以下、図面を参照し、本発明による絶縁膜内イオン挙動解析システムの実施例について説明する。   Embodiments of an ion behavior analysis system in an insulating film according to the present invention will be described below with reference to the drawings.

図4は、本実施例による絶縁膜内イオン挙動解析システムの構成を示す概略図である。絶縁膜内イオン挙動解析システムは、試料201を載置する試料台401、試料台401の試料201の載置面とは反対の面に設けられた圧延素子403、試料201を覆うように試料台401の上方に設けられたパルス電圧発生器402、パルス電圧発生器402の電源408、試料201に電圧または電流を印加する直流バイアス電源410、アンプ(Amp)を介して圧延素子403と接続されたデジタルオシロスコープ404、GP−IBによりデジタルオシロスコープ404に接続されたパソコン(PC)405、PC405に接続された画像処理・計算用のパソコン(PC)406、及びPC405と画像処理・計算用のPC406に接続された表示装置であるモニタ407を備える。   FIG. 4 is a schematic diagram showing the configuration of the ion behavior analysis system in the insulating film according to the present embodiment. The ion behavior analysis system in the insulating film includes a sample table 401 on which the sample 201 is mounted, a rolling element 403 provided on the surface of the sample table 401 opposite to the mounting surface of the sample 201, and the sample table so as to cover the sample 201. 401 connected to the rolling element 403 via a pulse voltage generator 402 provided above 401, a power source 408 of the pulse voltage generator 402, a DC bias power source 410 for applying voltage or current to the sample 201, and an amplifier (Amp). Connected to the digital oscilloscope 404, the personal computer (PC) 405 connected to the digital oscilloscope 404 by the GP-IB, the personal computer (PC) 406 for image processing / calculation connected to the PC 405, and the PC 405 and the PC 406 for image processing / calculation. A monitor 407 which is a display device is provided.

パルス電圧発生器402は、キャパシタンス(C)や抵抗(R)を有し、ケースに囲まれており、例えばパルス電圧波形411に示されるようなパルス電圧を試料201に印加する。印加したパルス電圧により試料201内に発生した超音波を用いて、試料201を非破壊で計測することができる。   The pulse voltage generator 402 has a capacitance (C) and a resistance (R), and is surrounded by a case. For example, a pulse voltage as indicated by a pulse voltage waveform 411 is applied to the sample 201. The sample 201 can be measured nondestructively using ultrasonic waves generated in the sample 201 by the applied pulse voltage.

直流バイアス電源410は、高圧側電極(HV)とグランド端子(G)を有し、高圧側電極からパルス電圧発生器402のケースを介して、試料201に電圧または電流を印加する。   The DC bias power supply 410 has a high voltage side electrode (HV) and a ground terminal (G), and applies a voltage or current to the sample 201 from the high voltage side electrode through the case of the pulse voltage generator 402.

PC405と画像処理・計算用のPC406は、どちらも演算装置であり、PC405は主として計測用に用いられ、画像処理・計算用のPC406は主として移流拡散計算を実施する。画像処理・計算用のPC406は、PC405を兼ねていてもよい。PC405と画像処理・計算用のPC406は、メモリやハードディスク等の記憶装置を備える。   Both the PC 405 and the image processing / calculation PC 406 are arithmetic units, the PC 405 is mainly used for measurement, and the image processing / calculation PC 406 mainly performs advection diffusion calculation. The image processing / calculation PC 406 may also serve as the PC 405. The PC 405 and the image processing / calculation PC 406 include storage devices such as a memory and a hard disk.

絶縁膜内イオン挙動解析システムの計測データと計算データは、PC405または画像処理・計算用のPC406の記憶装置に、時系列的に格納する。計測データと計算データを時系列的にPC405または画像処理・計算用のPC406に格納することで、データがおかしい場合には、計測環境や計測過程に問題がないか確認できる
図2と図3を用いて、本実施例における試料201について説明する。図2は、試料201の平面図(上面図)であり、図3は、試料201の断面図である。試料201は、銅電極からなるアノード304、絶縁膜305、及び銅基板からなるカソード306で構成される。以下に述べる実施例では、絶縁膜内イオン挙動解析システムにより、絶縁膜305の内部のイオン密度分布及び電界分布を求め、絶縁膜305の寿命を予測する。
Measurement data and calculation data of the ion behavior analysis system in the insulating film are stored in time series in the storage device of the PC 405 or the PC 406 for image processing / calculation. By storing measurement data and calculation data in time series in the PC 405 or the image processing / calculation PC 406, it is possible to check whether there is any problem in the measurement environment or the measurement process when the data is incorrect. The sample 201 in this example will be described. FIG. 2 is a plan view (top view) of the sample 201, and FIG. 3 is a cross-sectional view of the sample 201. The sample 201 includes an anode 304 made of a copper electrode, an insulating film 305, and a cathode 306 made of a copper substrate. In the embodiment described below, the ion density distribution and the electric field distribution inside the insulating film 305 are obtained by the ion behavior analysis system in the insulating film, and the lifetime of the insulating film 305 is predicted.

本発明の第1の実施例について、図1、図5、及び図6を用いて説明する。本実施例では、超音波を用いた非破壊の計測により、予め指定した任意の時刻における、試料201の絶縁膜305の内部のイオン密度分布及び電界分布を求める。   A first embodiment of the present invention will be described with reference to FIG. 1, FIG. 5, and FIG. In this embodiment, the ion density distribution and the electric field distribution inside the insulating film 305 of the sample 201 at an arbitrary time designated in advance are obtained by nondestructive measurement using ultrasonic waves.

図1は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。   FIG. 1 is a diagram showing a processing flow for obtaining an ion density distribution and an electric field distribution inside the insulating film 305 of the sample 201 in the ion behavior analysis system in the insulating film according to the present embodiment.

第1ステップ101にて、試料201を試料台401に設置する。そして、直流バイアス電源410を用いて、試料201に電圧または電流を印加する。電圧と電流のどちらを印加するかは、ユーザが任意に決めることができる。   In the first step 101, the sample 201 is set on the sample stage 401. Then, a voltage or current is applied to the sample 201 using the DC bias power supply 410. The user can arbitrarily determine which voltage or current is applied.

第2ステップ103では、予め定めた指定の時刻t1まで、試料201に電圧または電流を印加する。時刻t1に到達したら、試料201への電圧または電流の印加を止める。   In the second step 103, a voltage or current is applied to the sample 201 until a predetermined designated time t1. When the time t1 is reached, the voltage or current application to the sample 201 is stopped.

第3ステップ105にて、パルス電圧発生器402を用いて、パルス電圧を試料201に印加する。すると、試料201内のイオンの静電応力により超音波が発生するので、この超音波を圧延素子403で検出する。圧延素子403で検出された信号の大きさから、デジタルオシロスコープ404とPC405により、イオンの密度とイオンの座標値を求め、絶縁膜305の内部のイオン密度分布を求める。さらに、イオン密度分布から、ポアソン方程式により、絶縁膜305の内部の電界分布を求める。イオン密度分布や電界分布は、従来の方法と同様の方法を用いて求めることができる。   In the third step 105, a pulse voltage is applied to the sample 201 using the pulse voltage generator 402. Then, since an ultrasonic wave is generated by the electrostatic stress of ions in the sample 201, this ultrasonic wave is detected by the rolling element 403. From the magnitude of the signal detected by the rolling element 403, the digital oscilloscope 404 and the PC 405 obtain the ion density and the ion coordinate value, and obtain the ion density distribution inside the insulating film 305. Further, the electric field distribution inside the insulating film 305 is obtained from the ion density distribution by the Poisson equation. The ion density distribution and the electric field distribution can be obtained using a method similar to the conventional method.

第4ステップ106では、PC405にて、時刻t1でのイオンの密度、イオンの座標値、イオン密度分布、及び電界分布のデータをテーブル化する。テーブル化した上記のデータは、画像処理・計算用のPC406の記憶装置に格納する。   In the fourth step 106, the PC 405 tabulates the ion density, ion coordinate value, ion density distribution, and electric field distribution data at time t1. The above table data is stored in the storage device of the PC 406 for image processing / calculation.

ここで、第2ステップ103に戻り、直流バイアス電源410を用いて、予め定めた指定の時刻t2まで、試料201に電圧または電流を印加する。時刻t2は、時刻t1と異なる時刻であり、時刻t1より後の時刻(t2>t1)とする。時刻t2に到達したら、試料201への電圧または電流の印加を止める。   Here, returning to the second step 103, a voltage or current is applied to the sample 201 using the DC bias power supply 410 until a predetermined designated time t 2. The time t2 is a time different from the time t1, and is a time after the time t1 (t2> t1). When time t2 is reached, application of voltage or current to the sample 201 is stopped.

第3ステップ105にて、パルス電圧発生器402を用いて、パルス電圧を試料201に印加する。そして、上述したように、圧延素子403、デジタルオシロスコープ404、及びPC405により、絶縁膜305の内部のイオンの密度、イオンの座標値、イオン密度分布、及び電界分布を求める。   In the third step 105, a pulse voltage is applied to the sample 201 using the pulse voltage generator 402. Then, as described above, the ion density, the ion coordinate value, the ion density distribution, and the electric field distribution inside the insulating film 305 are obtained by the rolling element 403, the digital oscilloscope 404, and the PC 405.

第4ステップ106では、PC405にて、時刻t2でのイオンの密度、イオンの座標値、イオン密度分布、及び電界分布のデータをテーブル化する。テーブル化した上記のデータは、画像処理・計算用のPC406の記憶装置に格納する。   In the fourth step 106, the PC 405 tabulates the ion density, ion coordinate value, ion density distribution, and electric field distribution data at time t2. The above table data is stored in the storage device of the PC 406 for image processing / calculation.

第5ステップ107にて、画像処理・計算用のPC406で、時刻t1と時刻t2でのイオン密度分布と電界分布を基に移流拡散計算を実施して、予め指定した任意の時刻t3での絶縁膜305の内部のイオン密度分布を求める。移流拡散計算は、イオンの拡散と電界による移流を考慮した過渡解析であり、移流拡散方程式を解くことでイオン密度分布を求めることができる。本実施例に用いる移流拡散計算には、従来の移流拡散計算に用いられている計算方法と同様の方法を用いることができる。   In the fifth step 107, the PC 406 for image processing / calculation performs advection diffusion calculation based on the ion density distribution and electric field distribution at time t1 and time t2, and insulation at an arbitrary time t3 designated in advance. An ion density distribution inside the film 305 is obtained. The advection diffusion calculation is a transient analysis in consideration of ion diffusion and advection due to an electric field, and the ion density distribution can be obtained by solving the advection diffusion equation. For the advection diffusion calculation used in the present embodiment, a method similar to the calculation method used for the conventional advection diffusion calculation can be used.

本実施例では、2つの異なる時刻(時刻t1と時刻t2)でのイオン密度分布を用いて、移流拡散方程式のパラメータを求める。時間をt、場所をz、イオン密度分布をρ(t,z)で表すと、イオン密度分布ρ(t,z)は、正のイオン密度分布ρh(t,z)と負のイオン密度分布ρe(t,z)での和で構成される。
ρ(t,z)=ρh(t,z)+ρe(t,z) (1)
ρh(t,z)はプラスで、ρe(t,z)はマイナスのイオン密度分布であるとする。
In this embodiment, the parameters of the advection diffusion equation are obtained using ion density distributions at two different times (time t1 and time t2). Expressing time as t, location as z, and ion density distribution as ρ (t, z), the ion density distribution ρ (t, z) is a positive ion density distribution ρh (t, z) and a negative ion density distribution. It is composed of the sum of ρe (t, z).
ρ (t, z) = ρh (t, z) + ρe (t, z) (1)
It is assumed that ρh (t, z) is positive and ρe (t, z) is negative ion density distribution.

図5は、イオン密度分布の例を示したものであり、横軸が負電極からの位置、縦軸はイオン密度である。この図5にて、イオン密度がプラスの領域の分布はρh(t,z)、負のイオン密度分布はρe(t,z)である。ここでρh(t、z)とρe(t、z)が満足する移流拡散方程式は、それぞれ式(2)、式(3)で表される。さらに、イオン密度分布ρ(t,z)はポアソン方程式(4)を満足する。
dρh(t,z)/dt−d(Dhdρp(t,z)/dz)/dz
+μpE(t,z)dρh(t,z)/dz+ρh(t,z)/τh=0 (2)
dρe(t,z)/dt−d(Dedρe(t,z)/dz)/dz
−μeE(t,z)dρe(t,z)/dz+ρe(t,z)/τe=0 (3)
dE(t,z)/dz = q ρ(t,z) (4)
ここで、Dh=μh×kB×T/q、De=μe×kB×T/q、とする。E(t,z)は時間tと場所zでの電界の値、kBはボルツマン定数、Tは温度、qは電荷である。さらに、μh、μe、は、それぞれ正イオン移流係数、負イオン移流係数である。そして、τh、τeは、それぞれ正イオン寿命定数、負イオン寿命定数である。移流係数μh、μeと寿命定数τh、τeは未知のパラメータであり、その他は既知のパラメータである。
FIG. 5 shows an example of ion density distribution, where the horizontal axis is the position from the negative electrode, and the vertical axis is the ion density. In FIG. 5, the distribution in the positive ion density region is ρh (t, z), and the negative ion density distribution is ρe (t, z). Here, the advection diffusion equations satisfying ρh (t, z) and ρe (t, z) are expressed by Equation (2) and Equation (3), respectively. Furthermore, the ion density distribution ρ (t, z) satisfies the Poisson equation (4).
dρh (t, z) / dt−d (Dhdρp (t, z) / dz) / dz
+ ΜpE (t, z) dρh (t, z) / dz + ρh (t, z) / τh = 0 (2)
dρe (t, z) / dt−d (Dedρe (t, z) / dz) / dz
−μeE (t, z) dρe (t, z) / dz + ρe (t, z) / τe = 0 (3)
dE (t, z) / dz = qρ (t, z) (4)
Here, Dh = μh × kB × T / q and De = μe × kB × T / q. E (t, z) is the value of the electric field at time t and place z, kB is the Boltzmann constant, T is the temperature, and q is the charge. Further, μh and μe are a positive ion advection coefficient and a negative ion advection coefficient, respectively. Τh and τe are a positive ion lifetime constant and a negative ion lifetime constant, respectively. The advection coefficients μh and μe and the life constants τh and τe are unknown parameters, and the others are known parameters.

ここで、電界E(t,z)は、ポアソン方程式(4)の右辺にρ(t,z)を代入し、試料の構造、誘電率、及び電極での印加電位を考慮して求める。未知のパラメータである式(2)の移流係数μhと寿命定数τhは、次のようにして求める。
1)移流係数μhと寿命定数τhに適当な値を設定する。
2)時刻t1での正イオン密度分布の実験値Rh(t1,z)を初期値として移流拡散方程式(2)を解いて、時刻t2でのρh(t2,z)を求める。
3)ρh(t2,z)と時刻t2での正イオン密度分布の実験値Rh(t2,z)との差の絶対値を求める、
4)求めた絶対値の最大値が予め定めた閾値以下になるまで、上記の手順1)から3)を繰り返す。
以上の方法により、未知のパラメータである移流係数μhと寿命定数τhが求められる。
Here, the electric field E (t, z) is obtained by substituting ρ (t, z) for the right side of the Poisson equation (4) and taking into account the structure of the sample, the dielectric constant, and the applied potential at the electrode. The advection coefficient μh and the life constant τh in Equation (2), which are unknown parameters, are obtained as follows.
1) Set appropriate values for the advection coefficient μh and life constant τh.
2) Solve the advection diffusion equation (2) using the experimental value Rh (t1, z) of the positive ion density distribution at time t1 as an initial value, and obtain ρh (t2, z) at time t2.
3) Find the absolute value of the difference between ρh (t2, z) and the experimental value Rh (t2, z) of the positive ion density distribution at time t2.
4) The above steps 1) to 3) are repeated until the maximum value of the obtained absolute value is equal to or less than a predetermined threshold value.
By the above method, the advection coefficient μh and the life constant τh, which are unknown parameters, are obtained.

未知パラメータの移流係数μeと寿命定数τeに関しても、上記1)から4)と同様の手順で求めることができる。ただし、正イオン密度分布の実験値Rh(t1,z)を、負イオン密度分布の実験値Re(t1,z)に置き換え、移流拡散方程式(2)は移流拡散方程式(3)に置き換える。   The unknown parameter advection coefficient μe and life constant τe can also be obtained by the same procedure as in the above 1) to 4). However, the experimental value Rh (t1, z) of the positive ion density distribution is replaced with the experimental value Re (t1, z) of the negative ion density distribution, and the advection diffusion equation (2) is replaced with the advection diffusion equation (3).

これに加え、上記では場所zは一次元であったが、これを二次元、三次元に拡張できる。   In addition to this, the location z is one-dimensional in the above, but it can be expanded to two-dimensional and three-dimensional.

さらに、式(2)のρh(t,z)/τhと式(3)のρe(t,z)/τeの両者の項を、ρh(t,z)とρe(t,z)に依存する関数とすることが可能である。これらにより、方程式が現実のモデルに近くなり、ρhとρeを高精度に求めることができる。   Furthermore, both terms of ρh (t, z) / τh in equation (2) and ρe (t, z) / τe in equation (3) depend on ρh (t, z) and ρe (t, z). The function can be As a result, the equation becomes close to a real model, and ρh and ρe can be obtained with high accuracy.

このように、複数の異なる時刻でのイオン密度分布を基にパラメータを求めると、パラメータの精度が向上し、移流拡散方程式の解(任意の時刻t3でのイオン密度分布)も精度良く求めることができる。また、移流拡散方程式のパラメータを求めるのに用いるイオン密度分布の数は、2つだけに限らない。パラメータを求めるのに用いるイオン密度分布の数が多ければ多いほど、パラメータの精度が高くなり、移流拡散方程式の解もより精度良く求めることができる。したがって、3つ以上の異なる時刻でのイオン密度分布を基にして、移流拡散方程式のパラメータを求めてもよい。   Thus, when parameters are obtained based on ion density distributions at a plurality of different times, the accuracy of the parameters is improved, and the solution of the advection diffusion equation (ion density distribution at an arbitrary time t3) can be obtained with high accuracy. it can. Further, the number of ion density distributions used for obtaining the parameters of the advection diffusion equation is not limited to two. The greater the number of ion density distributions used to determine the parameters, the higher the accuracy of the parameters and the more accurate the solution of the advection diffusion equation. Therefore, the parameters of the advection diffusion equation may be obtained based on the ion density distribution at three or more different times.

以下の実施例では、μは、μhとμeのいずれか一方、または両方を示すとする。さらに、τは、τhとτeのいずれか一方、または両方を示すとする。   In the following examples, μ represents one or both of μh and μe. Further, τ represents one or both of τh and τe.

第5ステップ107では、さらに、任意の時刻t3でのイオン密度分布を基に、ポアソン方程式の計算により、時刻t3での絶縁膜305の内部の電界分布を求める。   In the fifth step 107, the electric field distribution inside the insulating film 305 at time t3 is further obtained by calculation of the Poisson equation based on the ion density distribution at arbitrary time t3.

第6ステップ108にて、任意の時刻t3での絶縁膜305の内部のイオン密度分布と電界分布を、モニタ407に表示する。   In the sixth step 108, the ion density distribution and electric field distribution inside the insulating film 305 at an arbitrary time t3 are displayed on the monitor 407.

図5は、本実施例により求めたイオン密度分布のグラフの例である。モニタ407には、絶縁膜305の内部のイオン密度分布として、図5に示すようなグラフが表示される。図5の横軸は、試料201において、絶縁膜305のカソード306からの距離であり、縦軸は、絶縁膜305の内部のイオン密度である。時刻t1のイオン密度104と時刻t2のイオン密度112は、図1の第3ステップ105にて、パルス電圧を試料201に印加することにより求めたものである。時刻t3のイオン密度109は、図1の第5ステップ107にて、時刻t1のイオン密度分布と時刻t2のイオン密度分布を基に移流拡散計算を実施して、求めたものである。   FIG. 5 is an example of a graph of ion density distribution obtained by this example. The monitor 407 displays a graph as shown in FIG. 5 as the ion density distribution inside the insulating film 305. The horizontal axis in FIG. 5 is the distance from the cathode 306 of the insulating film 305 in the sample 201, and the vertical axis is the ion density inside the insulating film 305. The ion density 104 at time t1 and the ion density 112 at time t2 are obtained by applying a pulse voltage to the sample 201 in the third step 105 of FIG. The ion density 109 at time t3 is obtained by performing advection diffusion calculation based on the ion density distribution at time t1 and the ion density distribution at time t2 in the fifth step 107 of FIG.

図6は、本実施例により求めた電界分布のグラフの例である。モニタ407には、絶縁膜305の内部の電界分布として、図6に示すようなグラフが表示される。図6の横軸は、試料201において、絶縁膜305のカソード306からの距離であり、縦軸は、絶縁膜305の内部の電界値の絶対値である電界強度に−1を掛けたものである。時刻t1の電界強度113と時刻t2の電界強度114は、図1の第3ステップ105にて、パルス電圧を試料201に印加して求めたイオン密度分布から、ポアソン方程式により求めたものである。時刻t3の電界強度115は、図1の第5ステップ107にて、時刻t3のイオン密度分布を基に、ポアソン方程式の計算により、求めたものである。   FIG. 6 is an example of a graph of the electric field distribution obtained by this example. The monitor 407 displays a graph as shown in FIG. 6 as the electric field distribution inside the insulating film 305. The horizontal axis in FIG. 6 is the distance from the cathode 306 of the insulating film 305 in the sample 201, and the vertical axis is the electric field strength that is the absolute value of the electric field value inside the insulating film 305 multiplied by −1. is there. The electric field intensity 113 at time t1 and the electric field intensity 114 at time t2 are obtained by Poisson's equation from the ion density distribution obtained by applying a pulse voltage to the sample 201 in the third step 105 in FIG. The electric field strength 115 at time t3 is obtained by calculation of the Poisson equation based on the ion density distribution at time t3 in the fifth step 107 of FIG.

以上のように、本実施例によれば、任意の時刻t3における絶縁膜305の内部のイオン密度分布及び電界分布を、時刻t3で測定せずに、導出できる。したがって、時刻t3を時刻t1や時刻t2より後の時刻(t3>t2>t1)とすると、現在より後の任意の時刻でのイオン密度分布を求めることができ、絶縁膜の寿命を予測することが可能である。   As described above, according to this embodiment, the ion density distribution and electric field distribution inside the insulating film 305 at an arbitrary time t3 can be derived without measuring at the time t3. Therefore, when the time t3 is a time after the time t1 or the time t2 (t3> t2> t1), the ion density distribution at an arbitrary time after the present can be obtained, and the lifetime of the insulating film is predicted. Is possible.

また、時刻t3が時刻t1や時刻t2より後の時刻(t3>t2>t1)の場合に、時刻t3での絶縁膜305の内部のイオン密度分布及び電界分布を知りたいときは、従来は時刻がt3になるまで待たなければならなかったが、本実施例では、時刻t1と時刻t2でのイオン密度分布及び電界分布から求めることができるので、時刻t3に到達するまでの時間(t3−t2)だけ時間を短縮できるという効果もある。   Further, when the time t3 is a time after the time t1 or the time t2 (t3> t2> t1), when it is desired to know the ion density distribution and the electric field distribution inside the insulating film 305 at the time t3, the conventional time is In this embodiment, since it can be obtained from the ion density distribution and the electric field distribution at time t1 and time t2, the time until reaching time t3 (t3-t2 ) Can also shorten the time.

なお、本実施例では時刻t1と時刻t2という異なる2つの時刻における測定から、任意の時刻t3でのイオン密度分布及び電界分布を求めた。上述したように、異なる3つ以上の時刻における測定から、時刻t3でのイオン密度分布及び電界分布を求めることもでき、この方がより精度良く時刻t3でのイオン密度分布及び電界分布を導出することができる。以下の実施例でも、本実施例と同様に、異なる時刻での測定回数が多ければ多いほど精度良く、時刻t3でのイオン密度分布及び電界分布を導出することができる。   In the present example, the ion density distribution and the electric field distribution at an arbitrary time t3 were obtained from the measurements at two different times, time t1 and time t2. As described above, the ion density distribution and electric field distribution at time t3 can also be obtained from measurements at three or more different times, which leads to more accurate derivation of the ion density distribution and electric field distribution at time t3. Can do. Also in the following examples, as in this example, the ion density distribution and electric field distribution at time t3 can be derived with higher precision as the number of measurements at different times increases.

本発明の第2の実施例について、図7を用いて説明する。本実施例では、予め指定した複数の任意の時刻における、試料201の絶縁膜305の内部のイオン密度分布及び電界分布を求め、電界分布の時間変化から試料201の絶縁膜305の寿命を予測する。   A second embodiment of the present invention will be described with reference to FIG. In this embodiment, the ion density distribution and the electric field distribution inside the insulating film 305 of the sample 201 are obtained at a plurality of predetermined times, and the lifetime of the insulating film 305 of the sample 201 is predicted from the time change of the electric field distribution. .

図7は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の寿命を予測する(予測寿命を求める)ための処理フローを示す図である。絶縁膜305の寿命とは、絶縁膜305が絶縁破壊を起こすまでの時間である。   FIG. 7 is a diagram showing a processing flow for predicting the life of the insulating film 305 of the sample 201 (determining the predicted life) of the ion behavior analysis system in the insulating film according to the present embodiment. The lifetime of the insulating film 305 is the time until the insulating film 305 causes dielectric breakdown.

第1ステップ701にて、試料201を試料台401に設置する。そして、直流バイアス電源410を用いて試料201に電圧または電流を印加する。電圧と電流のどちらを印加するかは、ユーザが任意に決めることができる。また、予測寿命tNを設定する。予測寿命tNは、初期値として任意の値を、画像処理・計算用のPC406に入力して設定する。さらに、画像処理・計算用のPC406には、絶縁膜305の絶縁破壊電界と、予測寿命tNの閾値を入力する。予測寿命tNの閾値は、予測寿命tNを求めるときの精度に関連する値であり、予め定めておく。   In the first step 701, the sample 201 is set on the sample stage 401. Then, a voltage or current is applied to the sample 201 using the DC bias power supply 410. The user can arbitrarily determine which voltage or current is applied. Also, a predicted life tN is set. The predicted life tN is set by inputting an arbitrary value to the PC 406 for image processing / calculation as an initial value. Furthermore, the dielectric breakdown electric field of the insulating film 305 and the threshold value of the predicted life tN are input to the PC 406 for image processing / calculation. The threshold value of the predicted life tN is a value related to the accuracy when obtaining the predicted life tN and is determined in advance.

第2ステップ703では、予め定めた指定の時間Δtが経過するまで、試料201に電圧または電流を印加する。時間Δtが経過したら、試料201への電圧または電流の印加を止める。   In the second step 703, voltage or current is applied to the sample 201 until a predetermined designated time Δt elapses. When the time Δt has elapsed, the application of voltage or current to the sample 201 is stopped.

第3ステップ705にて、パルス電圧発生器402を用いて、パルス電圧を試料201に印加する。そして、図1の第3ステップ105と同様にして、圧延素子403、デジタルオシロスコープ404、及びPC405により、絶縁膜305の内部のイオンの密度、イオンの座標値、イオン密度分布、及び電界分布を求める。   In the third step 705, a pulse voltage is applied to the sample 201 using the pulse voltage generator 402. In the same manner as in the third step 105 in FIG. 1, the ion density, ion coordinate value, ion density distribution, and electric field distribution inside the insulating film 305 are obtained by the rolling element 403, the digital oscilloscope 404, and the PC 405. .

第4ステップ706では、PC405にて、試料201へ電圧または電流を印加してから時間Δtが経過した時点でのイオンの密度、イオンの座標値、イオン密度分布、及び電界分布のデータをテーブル化する。テーブル化した上記のデータは、画像処理・計算用のPC406の記憶装置に格納する。   In the fourth step 706, the data of ion density, ion coordinate value, ion density distribution, and electric field distribution at the time when Δt has elapsed since the voltage or current was applied to the sample 201 by the PC 405 is tabulated. To do. The above table data is stored in the storage device of the PC 406 for image processing / calculation.

ここで、第2ステップ703に戻り、直流バイアス電源410を用いて、予め定めた指定の時間Δtが経過するまで、試料201に電圧または電流を印加する。時間Δtが経過したら、試料201への電圧または電流の印加を止める。   Here, returning to the second step 703, a voltage or current is applied to the sample 201 using the DC bias power supply 410 until a predetermined time Δt elapses. When the time Δt has elapsed, the application of voltage or current to the sample 201 is stopped.

第3ステップ705と第4ステップ706では、上述したのと同様の処理を行う。これにより、画像処理・計算用のPC406には、異なる2つの時刻における絶縁膜305の内部のイオンの密度、イオンの座標値、イオン密度分布、及び電界分布のデータが格納される。異なる2つの時刻における上記のデータが画像処理・計算用のPC406に格納されたら、第5ステップ707に進む。なお、第2ステップ703から第4ステップ706までの処理を3回以上繰り返し、異なる3つ以上の時刻における上記のデータを画像処理・計算用のPC406に格納してから第5ステップ707に進んでもよい。実施例1で述べたように、移流拡散方程式のパラメータを求めるのに用いるイオン密度分布の数が多ければ多いほど、パラメータの精度が高くなり、移流拡散方程式の解もより精度良く求めることができるからである。   In the third step 705 and the fourth step 706, the same processing as described above is performed. Thus, the image processing / calculation PC 406 stores the ion density, ion coordinate value, ion density distribution, and electric field distribution data inside the insulating film 305 at two different times. When the above data at two different times are stored in the image processing / calculation PC 406, the process proceeds to a fifth step 707. Note that the processing from the second step 703 to the fourth step 706 is repeated three times or more, and the above data at three or more different times are stored in the image processing / calculation PC 406 and then the processing proceeds to the fifth step 707. Good. As described in the first embodiment, as the number of ion density distributions used for obtaining the parameters of the advection diffusion equation increases, the accuracy of the parameters increases and the solution of the advection diffusion equation can be obtained more accurately. Because.

第5ステップ707では、画像処理・計算用のPC406で、今までに格納したイオン密度分布と電界分布を基に移流拡散計算を実施して、予め指定した複数の任意の時刻(t1,t2,・・・,2*tN)での絶縁膜305の内部のイオン密度分布を求める。移流拡散計算は、実施例1と同様に実施する。そして、時刻(t1,t2,・・・,2*tN)でのイオン密度分布を基に、ポアソン方程式の計算により、時刻(t1,t2,・・・,2*tN)での絶縁膜305の内部の電界分布をそれぞれ求める。   In the fifth step 707, the PC 406 for image processing / calculation performs the advection diffusion calculation based on the ion density distribution and the electric field distribution stored so far, and performs a plurality of arbitrary times (t1, t2, t2) designated in advance. .., 2 * tN), the ion density distribution inside the insulating film 305 is obtained. The advection diffusion calculation is performed in the same manner as in the first embodiment. Then, based on the ion density distribution at time (t1, t2,..., 2 * tN), the insulating film 305 at time (t1, t2,..., 2 * tN) is calculated based on the Poisson equation. The electric field distribution inside each is obtained.

第6ステップ708では、各時刻(t1,t2,・・・,2*tN)での絶縁膜305の内部の電界分布(電界強度)から、電界が絶縁膜305の絶縁破壊電界を超える最小の時刻を求める。そして、試料201に電圧または電流を印加した時刻からこの最小の時刻までの時間を、予測寿命tNとする。このようにして求めた予測寿命tNを、新たな予測寿命tNとする。例えば、自然数nで定められる時刻t(n)に対してt(n)<t(n+1)とし、時刻t(n)までは電界が絶縁破壊電界を下回っていたが、時刻t(n+1)以後では電界が絶縁破壊電界を超えた場合には、試料201に電圧または電流を印加した時刻から時刻t(n+1)までの時間を、予測寿命tNとする。このようにして、予測寿命tNの値は更新される。各時刻(t1,t2,・・・,2*tN)での絶縁膜305の内部のイオン密度分布と電界分布は、実施例1と同様に、モニタ407に表示する。このとき、予測寿命tNもモニタ407に表示する。予測寿命tNは、文字やグラフなど、テキスト表示やグラフィック表示により、表示することができる。   In the sixth step 708, from the electric field distribution (electric field strength) inside the insulating film 305 at each time (t1, t2,..., 2 * tN), the minimum electric field exceeds the breakdown electric field of the insulating film 305. Find the time. The time from the time when the voltage or current is applied to the sample 201 to the minimum time is defined as the predicted life tN. The predicted life tN obtained in this way is set as a new predicted life tN. For example, t (n) <t (n + 1) with respect to time t (n) determined by the natural number n, and the electric field was lower than the breakdown electric field until time t (n), but after time t (n + 1). Then, when the electric field exceeds the dielectric breakdown electric field, the time from the time when the voltage or current is applied to the sample 201 to the time t (n + 1) is defined as the predicted life tN. In this way, the value of the predicted life tN is updated. The ion density distribution and electric field distribution inside the insulating film 305 at each time (t1, t2,..., 2 * tN) are displayed on the monitor 407 as in the first embodiment. At this time, the predicted life tN is also displayed on the monitor 407. The predicted life tN can be displayed by text display or graphic display such as characters and graphs.

第7ステップ709では、前回の計算で求めた予測寿命tNと今回の計算で求めた予測寿命tNの差が、予測寿命tNの閾値以下であるかどうか判定する。差が閾値以下であれば、本処理を終了し、そうでなければ、第2ステップ703に戻り、第2ステップ703から第7ステップ709までの処理を繰り返す。このようにして、予測寿命tNの閾値で決められる精度で、予測寿命tNを求めることができる。   In a seventh step 709, it is determined whether or not the difference between the predicted life tN obtained in the previous calculation and the predicted life tN obtained in the current calculation is equal to or less than the threshold value of the predicted life tN. If the difference is less than or equal to the threshold value, the process ends. If not, the process returns to the second step 703 and the processes from the second step 703 to the seventh step 709 are repeated. In this way, the predicted life tN can be obtained with an accuracy determined by the threshold value of the predicted life tN.

本実施例によると、異なる複数の時刻における絶縁膜305の内部のイオン密度分布及び電界分布から、複数の任意の時刻における絶縁膜305の内部の電界分布を求めることができるので、絶縁破壊電界以上になるまで測定を続けなくても、予測寿命tNを推定することが可能である。   According to this embodiment, since the electric field distribution inside the insulating film 305 at a plurality of arbitrary times can be obtained from the ion density distribution and the electric field distribution inside the insulating film 305 at a plurality of different times, the dielectric breakdown electric field or more can be obtained. The estimated life tN can be estimated without continuing the measurement until.

なお、実施例1及び以下の実施例で述べた方法で求めた絶縁膜305の内部のイオン密度分布及び電界分布を用いても、本実施例で述べた方法と同様の方法により、絶縁膜305の予測寿命tNを推定することが可能である。   Note that even if the ion density distribution and the electric field distribution inside the insulating film 305 obtained by the method described in Embodiment 1 and the following embodiments are used, the insulating film 305 is obtained by the same method as that described in this embodiment. Can be estimated.

本発明の第3の実施例について、図8を用いて説明する。本実施例では、試料201に電圧または電流を印加したままで、絶縁膜305の内部のイオン密度分布及び電界分布を求める。なお、本明細書では、試料201に電圧または電流を印加したままの状態のことを、「オンライン」と称する。   A third embodiment of the present invention will be described with reference to FIG. In this embodiment, the ion density distribution and the electric field distribution inside the insulating film 305 are obtained while applying a voltage or current to the sample 201. Note that in this specification, a state where a voltage or a current is applied to the sample 201 is referred to as “online”.

図8は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。   FIG. 8 is a diagram showing a processing flow for obtaining the ion density distribution and electric field distribution inside the insulating film 305 of the sample 201 in the ion behavior analysis system in the insulating film according to the present embodiment.

第1ステップ801にて、試料201を試料台401に載置する。そして、直流バイアス電源410を用いて、試料201に電圧または電流を印加する。電圧と電流のどちらを印加するかは、ユーザが任意に決めることができる。   In the first step 801, the sample 201 is placed on the sample table 401. Then, a voltage or current is applied to the sample 201 using the DC bias power supply 410. The user can arbitrarily determine which voltage or current is applied.

第2ステップ803では、予め定めた指定の時刻t1まで、試料201に電圧または電流を印加する。   In the second step 803, voltage or current is applied to the sample 201 until a predetermined designated time t1.

第3ステップ805にて、時刻t1に到達したら、試料201に電圧または電流を印加したまま、すなわちオンラインで、パルス電圧発生器402を用いて、パルス電圧を試料201に印加する。そして、オンラインで、実施例1の第3ステップ105と同様にして、絶縁膜305の内部のイオン密度分布と電界分布を求める。   When the time t1 is reached in the third step 805, the pulse voltage is applied to the sample 201 using the pulse voltage generator 402 while applying the voltage or current to the sample 201, that is, online. Then, the ion density distribution and the electric field distribution inside the insulating film 305 are obtained online in the same manner as in the third step 105 of the first embodiment.

第4ステップ806から第6ステップ808までの処理は、実施例1の第4ステップ106から第6ステップ108までの処理と同様なので、説明を省略する。   Since the processing from the fourth step 806 to the sixth step 808 is the same as the processing from the fourth step 106 to the sixth step 108 of the first embodiment, the description thereof is omitted.

本実施例では、試料201を検査ごとに取り外したり、電圧または電流の印加を断続したりせず、試料201に電圧または電流を印加したままという実際の動作に近い条件の下でイオン密度分布及び電界分布を求めることができる。したがって、本実施例による絶縁膜内イオン挙動解析システムは、実施例2で述べた方法で絶縁膜305の予測寿命tNを推定することにより、実施例1で述べた効果に加えてより高精度に絶縁膜の寿命を予測することが可能であるという効果を有する。   In this embodiment, the sample 201 is not removed for each inspection, and the application of voltage or current is not interrupted, and the ion density distribution and the conditions are close to the actual operation of applying the voltage or current to the sample 201. The electric field distribution can be obtained. Therefore, the ion behavior analysis system in the insulating film according to the present embodiment estimates the predicted lifetime tN of the insulating film 305 by the method described in the second embodiment, so that the effect described in the first embodiment is increased with higher accuracy. It has the effect that it is possible to predict the lifetime of the insulating film.

本発明の第4の実施例について、図9を用いて説明する。本実施例では、試料201の絶縁膜305の内部のイオン密度分布及び電界分布だけでなく、試料201の絶縁膜305の内部の電流分布と電位分布も求める。   A fourth embodiment of the present invention will be described with reference to FIG. In this embodiment, not only the ion density distribution and electric field distribution inside the insulating film 305 of the sample 201 but also the current distribution and potential distribution inside the insulating film 305 of the sample 201 are obtained.

図9は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の内部のイオン密度分布、電界分布、電流分布、及び電位分布を求めるための処理フローを示す図である。   FIG. 9 is a diagram showing a processing flow for obtaining the ion density distribution, electric field distribution, current distribution, and potential distribution inside the insulating film 305 of the sample 201 in the ion behavior analysis system in the insulating film according to the present embodiment. .

第1ステップ901から第3ステップ905までの処理は、実施例3の第1ステップ801から第3ステップ805までの処理と同様なので、説明を省略する。   Since the processing from the first step 901 to the third step 905 is the same as the processing from the first step 801 to the third step 805 of the third embodiment, description thereof will be omitted.

第4ステップ906にて、試料201の電流の波形と電圧の波形を電流・電圧センサーやプローブを用いて測定する。そして、PC405にて、時刻t1と時刻t2での絶縁膜305内の電流の波形と電圧の波形を求める。   In the fourth step 906, the current waveform and voltage waveform of the sample 201 are measured using a current / voltage sensor or a probe. Then, the PC 405 obtains a current waveform and a voltage waveform in the insulating film 305 at time t1 and time t2.

第5ステップ907では、実施例3の第4ステップ806と同様の処理を行う。すなわち、PC405にて、時刻t1と時刻t2でのイオンの密度、イオンの座標値、イオン密度分布、及び電界分布のデータをテーブル化し、画像処理・計算用のPC406の記憶装置に格納する。さらに、本実施例では、第4ステップ906で求めた時刻t1と時刻t2での絶縁膜305内の電流の波形と電圧の波形もテーブル化し、画像処理・計算用のPC406の記憶装置に格納する。   In the fifth step 907, processing similar to that in the fourth step 806 of the third embodiment is performed. That is, the PC 405 tabulates the ion density, ion coordinate value, ion density distribution, and electric field distribution data at time t1 and time t2 and stores them in the storage device of the PC 406 for image processing / calculation. Further, in this embodiment, the current waveform and voltage waveform in the insulating film 305 at time t1 and time t2 obtained in the fourth step 906 are also tabulated and stored in the storage device of the PC 406 for image processing / calculation. .

第6ステップ908にて、画像処理・計算用のPC406で、時刻t1と時刻t2でのイオン密度分布と電界分布と電流の波形と電圧の波形を基に移流拡散計算を実施して、予め指定した任意の時刻t3での絶縁膜305の内部のイオン密度分布を求める。移流拡散計算は、実施例1と同様に実施する。移流拡散方程式のパラメータ(式1の移流係数μ、寿命定数τ、及び電界E)は、時刻t1と時刻t2でのイオン密度分布と電界分布と電流の波形と電圧の波形を用いて求める。そして、時刻t3でのイオン密度分布を基に、ポアソン方程式の計算により、時刻t3での絶縁膜305の内部の電界分布、電流分布、及び電位分布を求める。   In the sixth step 908, the PC 406 for image processing / calculation performs advection diffusion calculation based on the ion density distribution, electric field distribution, current waveform, and voltage waveform at time t1 and time t2, and designates in advance. The ion density distribution inside the insulating film 305 at the given time t3 is obtained. The advection diffusion calculation is performed in the same manner as in the first embodiment. The parameters of the advection diffusion equation (the advection coefficient μ, the life constant τ, and the electric field E in Equation 1) are obtained using the ion density distribution, electric field distribution, current waveform, and voltage waveform at time t1 and time t2. Then, based on the ion density distribution at time t3, the electric field distribution, current distribution, and potential distribution inside the insulating film 305 at time t3 are obtained by calculation of the Poisson equation.

第7ステップ909にて、任意の時刻t3での絶縁膜305の内部のイオン密度分布、電界分布、電流分布、及び電位分布をモニタ407に表示する。   In a seventh step 909, the ion density distribution, electric field distribution, current distribution, and potential distribution inside the insulating film 305 at an arbitrary time t3 are displayed on the monitor 407.

本実施例では、絶縁膜305の内部のイオン密度分布と電界分布だけではなく、電流分布と電位分布も求めるので、絶縁破壊電界だけではなく、絶縁破壊電流も求めることが可能である。したがって、本実施例による絶縁膜内イオン挙動解析システムは、実施例2で述べた方法で絶縁膜305の予測寿命tNを推定することにより、絶縁膜の寿命をより高精度に予測することが可能である。   In this embodiment, since not only the ion density distribution and electric field distribution inside the insulating film 305 but also the current distribution and potential distribution are obtained, it is possible to obtain not only the dielectric breakdown electric field but also the dielectric breakdown current. Therefore, the ion behavior analysis system in the insulating film according to the present embodiment can predict the life of the insulating film with higher accuracy by estimating the predicted life tN of the insulating film 305 by the method described in the second embodiment. It is.

なお、本実施例では時刻t1と時刻t2という異なる2つの時刻における測定から、任意の時刻t3での絶縁膜305の内部のイオン密度分布、電界分布、電流分布、及び電位分布を求めた。絶縁膜305の内部の電流分布と電位分布も、イオン密度分布や電界分布と同様に、異なる3つ以上の時刻における測定から、時刻t3での分布を求めることができ、この方がより精度良く電流分布と電位分布を導出することができる。   Note that in this example, the ion density distribution, electric field distribution, current distribution, and potential distribution inside the insulating film 305 at an arbitrary time t3 were obtained from measurements at two different times of time t1 and time t2. Similarly to the ion density distribution and the electric field distribution, the current distribution and the potential distribution inside the insulating film 305 can be obtained at the time t3 from the measurements at three or more different times, which is more accurate. Distribution and potential distribution can be derived.

本発明の第5の実施例について、図10と図11を用いて説明する。本実施例は、試料201が、パルス電圧発生器402より大きい場合の処理の例である。   A fifth embodiment of the present invention will be described with reference to FIGS. The present embodiment is an example of processing when the sample 201 is larger than the pulse voltage generator 402.

図11は、本実施例による絶縁膜内イオン挙動解析システムの、パルス電圧発生器402と試料台401を示す図である。試料台401には、パルス電圧発生器402より大きい試料201が載置されている。試料201は、絶縁膜305の表面の計測領域が、点線で表された複数の計測領域901に分割される。分割された計測領域901の大きさは、パルス電圧発生器402の大きさに基づいて定めることができる。図4に示したように、試料台401の試料201の載置面とは反対の面には、圧延素子403が設けられている。本実施例での絶縁膜内イオン挙動解析システムでは、試料台401に試料板稼動装置902が設けられている。試料板稼動装置902は、パルス電圧発生器402が試料201を走査し、絶縁膜305に設定された複数の計測領域901のそれぞれについて計測できるように、試料台401を前後左右に移動させる。   FIG. 11 is a diagram showing a pulse voltage generator 402 and a sample stage 401 of the ion behavior analysis system in an insulating film according to the present embodiment. A sample 201 larger than the pulse voltage generator 402 is placed on the sample table 401. In the sample 201, the measurement region on the surface of the insulating film 305 is divided into a plurality of measurement regions 901 represented by dotted lines. The size of the divided measurement region 901 can be determined based on the size of the pulse voltage generator 402. As shown in FIG. 4, a rolling element 403 is provided on the surface of the sample table 401 opposite to the surface on which the sample 201 is placed. In the insulating film ion behavior analysis system in the present embodiment, a sample plate operating device 902 is provided on a sample stage 401. The sample plate operating device 902 moves the sample stage 401 back and forth and right and left so that the pulse voltage generator 402 can scan the sample 201 and measure each of the plurality of measurement regions 901 set in the insulating film 305.

図10は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。   FIG. 10 is a diagram showing a processing flow for obtaining the ion density distribution and electric field distribution inside the insulating film 305 of the sample 201 in the ion behavior analysis system in the insulating film according to the present embodiment.

第1ステップ1001にて、試料201を試料台401に載置する。そして、直流バイアス電源410を用いて、試料201に電圧または電流を印加する。電圧と電流のどちらを印加するかは、ユーザが任意に決めることができる。   In the first step 1001, the sample 201 is placed on the sample table 401. Then, a voltage or current is applied to the sample 201 using the DC bias power supply 410. The user can arbitrarily determine which voltage or current is applied.

第2ステップ1003では、予め定めた指定の時刻t1まで、試料201に電圧または電流を印加する。時刻t1に到達したら、試料201への電圧または電流の印加を止める。   In the second step 1003, voltage or current is applied to the sample 201 until a predetermined designated time t1. When the time t1 is reached, the voltage or current application to the sample 201 is stopped.

第3ステップ1005にて、試料201の表面の計測領域を、複数の計測領域901に分割する。分割された全ての計測領域901にて、超音波計測が可能であるとする。試料台稼働装置902を用いて試料台401を動かし、全ての計測領域901にパルス電圧発生器402でパルス電圧を印加し、各計測領域901の内部のイオン密度分布と電界分布を求める。各計測領域901では、十分信号がとれるまで試料台401を固定して、信号を測定する。各計測領域901でのイオン密度分布と電界分布を求める方法は、実施例1の第3ステップ105と同様である。   In the third step 1005, the measurement area on the surface of the sample 201 is divided into a plurality of measurement areas 901. It is assumed that ultrasonic measurement is possible in all the divided measurement areas 901. The sample stage 401 is moved using the sample stage operating device 902, and a pulse voltage is applied to all the measurement areas 901 by the pulse voltage generator 402 to obtain the ion density distribution and electric field distribution inside each measurement area 901. In each measurement region 901, the sample stage 401 is fixed until a sufficient signal is obtained, and the signal is measured. The method for obtaining the ion density distribution and the electric field distribution in each measurement region 901 is the same as in the third step 105 of the first embodiment.

第4ステップ1006では、PC405にて、全ての計測領域901について、時刻t1でのイオンの密度、イオンの座標値、イオン密度分布、及び電界分布のデータをテーブル化する。テーブル化した上記のデータは、画像処理・計算用のPC406の記憶装置に格納する。   In the fourth step 1006, the PC 405 tabulates the data of the ion density, ion coordinate value, ion density distribution, and electric field distribution at time t1 for all measurement regions 901. The above table data is stored in the storage device of the PC 406 for image processing / calculation.

ここで、第2ステップ1003に戻り、直流バイアス電源410を用いて、予め定めた指定の時刻t2まで、試料201に電圧または電流を印加する。時刻t2は、時刻t1と異なる時刻であり、時刻t1より後の時刻(t2>t1)とする。時刻t2に到達したら、試料201への電圧または電流の印加を止める。   Here, returning to the second step 1003, a voltage or current is applied to the sample 201 using the DC bias power supply 410 until a predetermined designated time t 2. The time t2 is a time different from the time t1, and is a time after the time t1 (t2> t1). When time t2 is reached, application of voltage or current to the sample 201 is stopped.

そして、第3ステップ1005と第4ステップ1006を行い、全ての計測領域901について、時刻t2でのイオンの密度、イオンの座標値、イオン密度分布、及び電界分布を求め、画像処理・計算用のPC406の記憶装置に格納する。   Then, the third step 1005 and the fourth step 1006 are performed, and the ion density, the ion coordinate value, the ion density distribution, and the electric field distribution at the time t2 are obtained for all measurement regions 901, and image processing / calculation is performed. The data is stored in the storage device of the PC 406.

第5ステップ1007にて、画像処理・計算用のPC406で、時刻t1と時刻t2でのイオン密度分布と電界分布を基に、全ての計測領域901での移流拡散計算を実施する。これにより、予め指定した任意の時刻t3での全ての計測領域901での絶縁膜305の内部のイオン密度分布と電界分布を求める。移流拡散計算やイオン密度分布と電界分布を求める方法は、実施例1の第5ステップ107と同様である。   In a fifth step 1007, the advection diffusion calculation is performed in all measurement regions 901 on the image processing / calculation PC 406 based on the ion density distribution and electric field distribution at time t1 and time t2. Thereby, the ion density distribution and electric field distribution inside the insulating film 305 in all the measurement regions 901 at an arbitrary time t3 designated in advance are obtained. The advection diffusion calculation and the method for obtaining the ion density distribution and the electric field distribution are the same as in the fifth step 107 of the first embodiment.

第6ステップ1008にて、任意の時刻t3での全ての計測領域901での絶縁膜305の内部のイオン密度分布と電界分布を、モニタ407に表示する。   In a sixth step 1008, the ion density distribution and electric field distribution inside the insulating film 305 in all measurement regions 901 at an arbitrary time t3 are displayed on the monitor 407.

本実施例では、試料201が、パルス電圧発生器402より大きい場合でも、試料201の分割や破壊をせずに、イオン密度分布と電界分布を求めることができる。   In this embodiment, even when the sample 201 is larger than the pulse voltage generator 402, the ion density distribution and the electric field distribution can be obtained without dividing or destroying the sample 201.

さらに、超音波計測では深さ方向のイオン密度分布が計測できるので、試料201の絶縁膜305内のイオンの3次元分布を求めることが可能である。   Furthermore, since the ion density distribution in the depth direction can be measured by ultrasonic measurement, the three-dimensional distribution of ions in the insulating film 305 of the sample 201 can be obtained.

なお、本実施例では試料台401に試料板稼動装置902を設けて、試料台401を前後左右に移動させるようにしたが、試料台401を固定してパルス電圧発生器402と圧延素子403を前後左右に移動させるようにしてもよい。この場合には、パルス電圧発生器402と圧延素子403に稼動装置を設けて、パルス電圧発生器402が試料201を走査することができるようにする。   In this embodiment, the sample plate operating device 902 is provided on the sample table 401 and the sample table 401 is moved back and forth and left and right. However, the sample table 401 is fixed and the pulse voltage generator 402 and the rolling element 403 are connected. You may make it move back and forth and right and left. In this case, the pulse voltage generator 402 and the rolling element 403 are provided with operating devices so that the pulse voltage generator 402 can scan the sample 201.

本発明の第6の実施例について、図11と図12を用いて説明する。本実施例は、試料201が、パルス電圧発生器402より大きい場合の処理の例であり、絶縁膜305の寿命(絶縁寿命)を予測して、絶縁寿命が最も短い領域を再成型する例である。   A sixth embodiment of the present invention will be described with reference to FIGS. The present embodiment is an example of processing when the sample 201 is larger than the pulse voltage generator 402, and is an example in which the lifetime of the insulating film 305 (insulating lifetime) is predicted and the region having the shortest insulating lifetime is remolded. is there.

本実施例による絶縁膜内イオン挙動解析システムの、パルス電圧発生器402と試料台401は、図11に示しており、実施例5と同じなので説明を省略する。   The pulse voltage generator 402 and the sample stage 401 of the ion behavior analysis system in the insulating film according to the present embodiment are shown in FIG.

図12は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の内部のイオン密度分布及び電界分布を求め、絶縁寿命が最も短い領域を再成型するための処理フローを示す図である。   FIG. 12 shows a process flow for obtaining the ion density distribution and the electric field distribution inside the insulating film 305 of the sample 201 and re-molding the region with the shortest insulation lifetime in the ion behavior analysis system in the insulating film according to the present embodiment. FIG.

第1ステップ1201から第6ステップ1208までの処理は、実施例5の第1ステップ1001から第6ステップ1008(図10)までの処理と同一なので、説明を省略する。   Since the processing from the first step 1201 to the sixth step 1208 is the same as the processing from the first step 1001 to the sixth step 1008 (FIG. 10) of the fifth embodiment, description thereof is omitted.

第7ステップ1209にて、全ての計測領域901について、絶縁膜305の内部のイオン密度分布と電界分布から、絶縁膜305の寿命(絶縁寿命)を求める。絶縁膜305の寿命は、実施例2で説明した方法と同様にして求める。すなわち、予測寿命tNを予め設定し、複数の任意の時刻での絶縁膜305の内部のイオン密度分布と電界分布を求め、これらの値から絶縁膜305の寿命を求める。本実施例では、実施例5での時刻t3を、複数の異なる時刻(t3,t3’ ,t3’’,・・・)とする。全ての計測領域901について絶縁膜305の絶縁寿命を求めた後、絶縁寿命が最も短い計測領域409(図11参照)を求めて、モニタ407に表示する。絶縁寿命が最も短い計測領域409の表示方法としては、例えば、図11に示したように塗りつぶす色を変えて領域を強調表示するなど、他の領域と区別できるように表示する方法が考えられる。   In the seventh step 1209, the lifetime (insulation lifetime) of the insulating film 305 is obtained from the ion density distribution and the electric field distribution inside the insulating film 305 for all measurement regions 901. The lifetime of the insulating film 305 is obtained in the same manner as the method described in the second embodiment. That is, the predicted life tN is set in advance, the ion density distribution and the electric field distribution inside the insulating film 305 at a plurality of arbitrary times are obtained, and the life of the insulating film 305 is obtained from these values. In the present embodiment, the time t3 in the fifth embodiment is set to a plurality of different times (t3, t3 ′, t3 ″,...). After obtaining the insulation life of the insulating film 305 for all the measurement regions 901, the measurement region 409 (see FIG. 11) with the shortest insulation life is obtained and displayed on the monitor 407. As a display method of the measurement region 409 having the shortest insulation life, for example, a method of displaying the measurement region 409 so that it can be distinguished from other regions, such as highlighting the region by changing the color to be painted as shown in FIG.

第8ステップ1210にて、絶縁寿命が最も短い計測領域409を再成型する。再成型の方法の一例としては、絶縁寿命が最も短い計測領域409に電子線やイオンビーム、レーザを照射する等で表面処理をする方法が挙げられる。または、絶縁寿命が最も短い計測領域409を絶縁膜305から取り除き、取り除いた部分には新たな絶縁膜を設けてもよい。   In the eighth step 1210, the measurement region 409 having the shortest insulation life is remolded. As an example of the remolding method, there is a method in which a surface treatment is performed by irradiating an electron beam, an ion beam, or a laser to the measurement region 409 having the shortest insulation life. Alternatively, the measurement region 409 having the shortest insulation life may be removed from the insulating film 305, and a new insulating film may be provided in the removed portion.

本実施例によると、絶縁膜305の一部の寿命が短くなった場合でも、絶縁膜305の全てを廃棄したり置き換えたりすることなく、寿命が短くなった部分を再成型することができる。これにより、絶縁膜305の寿命の向上が、効率良く低コストで可能となる。   According to this embodiment, even when the lifetime of a part of the insulating film 305 is shortened, the portion having a shortened lifetime can be remolded without discarding or replacing the entire insulating film 305. Thereby, the lifetime of the insulating film 305 can be improved efficiently and at low cost.

本発明の第7の実施例について、図13を用いて説明する。本実施例は、絶縁膜305の内部のイオン密度分布及び電界分布を求めるのに、絶縁膜305の構造・材料データを用いる場合の処理の例である。   A seventh embodiment of the present invention will be described with reference to FIG. This embodiment is an example of processing when the structure / material data of the insulating film 305 is used to obtain the ion density distribution and electric field distribution inside the insulating film 305.

図13は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の内部のイオン密度分布及び電界分布を求めるための処理フローを示す図である。   FIG. 13 is a diagram showing a processing flow for obtaining the ion density distribution and electric field distribution inside the insulating film 305 of the sample 201 of the ion behavior analysis system in the insulating film according to the present embodiment.

第1ステップ1301にて、試料201の絶縁膜305の構造・材料データを、PC405に格納する。本実施例では、構造・材料データの一例として、絶縁膜305の材料の移流係数μ、拡散係数D、及び絶縁寿命定数τ(または消滅定数(1/τ))のうち少なくとも1つを格納する。絶縁膜305の材料としては、例えば、エポキシ、炭素系ポリマー、SiO、Si、SiN、またはAlを用いることができる。絶縁膜305の構造・材料データ(移流係数μ、拡散係数D、及び絶縁寿命定数τ(または消滅定数(1/τ))のうち少なくとも1つ)は、移流拡散計算に使う移流拡散方程式のパラメータとして使用する。 In the first step 1301, the structure / material data of the insulating film 305 of the sample 201 is stored in the PC 405. In this embodiment, at least one of the advection coefficient μ, the diffusion coefficient D, and the insulation lifetime constant τ (or the extinction constant (1 / τ)) of the material of the insulating film 305 is stored as an example of the structure / material data. . As a material of the insulating film 305, for example, epoxy, carbon polymer, SiO 2 , Si, SiN, or Al 2 O 3 can be used. The structure / material data of the insulating film 305 (at least one of the advection coefficient μ, the diffusion coefficient D, and the insulation lifetime constant τ (or the extinction constant (1 / τ)) is a parameter of the advection diffusion equation used for the advection diffusion calculation. Use as

第2ステップ1303から第5ステップ1307までの処理は、実施例1の第1ステップ101から第4ステップ106(図1)までの処理と同一なので、説明を省略する。   Since the processing from the second step 1303 to the fifth step 1307 is the same as the processing from the first step 101 to the fourth step 106 (FIG. 1) of the first embodiment, description thereof is omitted.

第6ステップ1308にて、絶縁膜305の構造・材料データを、PC405から読み出す。   In the sixth step 1308, the structure / material data of the insulating film 305 is read from the PC 405.

第7ステップ1309と第8ステップ1310の処理は、実施例1の第5ステップ107と第6ステップ108(図1)の処理と同一なので、説明を省略する。ただし、第7ステップ1309で移流拡散計算を実施する際には、絶縁膜305の構造・材料データを移流拡散方程式のパラメータとして使用する。   Since the processes of the seventh step 1309 and the eighth step 1310 are the same as the processes of the fifth step 107 and the sixth step 108 (FIG. 1) of the first embodiment, description thereof is omitted. However, when the advection diffusion calculation is performed in the seventh step 1309, the structure / material data of the insulating film 305 is used as a parameter of the advection diffusion equation.

本実施例によると、絶縁膜305の構造・材料データも考慮して移流拡散計算を実施するので、絶縁膜305の内部のイオン密度分布及び電界分布をより高精度に求めることができる。したがって、本実施例による絶縁膜内イオン挙動解析システムは、実施例2で述べた方法で絶縁膜305の予測寿命tNを推定することにより、絶縁膜305の寿命もより高精度に予測することが可能となる。   According to this embodiment, since the advection diffusion calculation is performed in consideration of the structure / material data of the insulating film 305, the ion density distribution and electric field distribution inside the insulating film 305 can be obtained with higher accuracy. Therefore, the ion behavior analysis system in the insulating film according to the present embodiment can predict the lifetime of the insulating film 305 with higher accuracy by estimating the predicted lifetime tN of the insulating film 305 by the method described in the second embodiment. It becomes possible.

本発明の第8の実施例について、図14を用いて説明する。本実施例では、試料201の絶縁膜305の内部のイオン密度分布の温度依存性を求める。   An eighth embodiment of the present invention will be described with reference to FIG. In this embodiment, the temperature dependence of the ion density distribution inside the insulating film 305 of the sample 201 is obtained.

図14は、本実施例による絶縁膜内イオン挙動解析システムの、試料201の絶縁膜305の内部のイオン密度分布の温度依存性を求めるための処理フローを示す図である。   FIG. 14 is a diagram showing a processing flow for determining the temperature dependence of the ion density distribution inside the insulating film 305 of the sample 201 in the ion behavior analysis system in the insulating film according to the present embodiment.

第1ステップ1401にて、試料201の温度を変化させる。試料201の温度を変化させる方法としては、例えば、熱ヒータにより試料201へ熱を印加する方法や、赤外線レーザ装置により試料201に赤外光を照射する方法が挙げられる。このとき、必要に応じて、パルス電圧発生器402を試料201の上方から移動させて、熱を印加したり赤外光を照射したりする。   In the first step 1401, the temperature of the sample 201 is changed. Examples of a method for changing the temperature of the sample 201 include a method of applying heat to the sample 201 with a thermal heater and a method of irradiating the sample 201 with infrared light using an infrared laser device. At this time, if necessary, the pulse voltage generator 402 is moved from above the sample 201 to apply heat or irradiate infrared light.

第2ステップ1403から第5ステップ1407までの処理は、実施例1の第1ステップ101から第4ステップ106(図1)までの処理と同一なので、説明を省略する。   Since the processing from the second step 1403 to the fifth step 1407 is the same as the processing from the first step 101 to the fourth step 106 (FIG. 1) of the first embodiment, description thereof is omitted.

第6ステップ1408では、画像処理・計算用のPC406で、時刻t1と時刻t2でのイオン密度分布と電界分布を基に移流拡散計算を実施して、絶縁膜305の内部のイオン密度分布の温度依存性を求める。温度依存性は、第1ステップ1401で変化させた試料201の各温度でのイオン密度分布と電界分布から求める。そして、イオンの移流係数と拡散係数の温度依存性を導出し、表示する。イオンの移流係数と拡散係数の温度依存性は、試料201の温度を変化させたときのイオン密度分布と、これに対応するようにフィッテイングした移流拡散方程式における移流係数μ(μhまたはμe)と寿命係数τ(τhまたはτe)から、移流係数μ(μhまたはμe)と寿命係数τ(τhまたはτe)の温度依存性を導出して求める。イオンの移流係数と寿命係数の温度依存性は、例えば、横軸を温度、縦軸をμ、τとしたグラフで表示することができる。   In a sixth step 1408, the advancing diffusion calculation is performed based on the ion density distribution and the electric field distribution at the time t1 and the time t2 on the image processing / calculation PC 406, and the temperature of the ion density distribution inside the insulating film 305 is calculated. Dependency is sought. The temperature dependence is obtained from the ion density distribution and the electric field distribution at each temperature of the sample 201 changed in the first step 1401. Then, the temperature dependence of the ion advection coefficient and diffusion coefficient is derived and displayed. The temperature dependence of the advection coefficient and diffusion coefficient of ions is determined by the ion density distribution when the temperature of the sample 201 is changed, and the advection coefficient μ (μh or μe) in the advection diffusion equation fitted so as to correspond to this. The temperature dependence of the advection coefficient μ (μh or μe) and the life coefficient τ (τh or τe) is derived from the life coefficient τ (τh or τe). The temperature dependence of the ion advection coefficient and lifetime coefficient can be displayed, for example, as a graph with the horizontal axis representing temperature and the vertical axis representing μ and τ.

本実施例によると、試料201の絶縁膜305の内部のイオン密度分布の温度依存性を求めることができる。したがって、本実施例による絶縁膜内イオン挙動解析システムは、実施例2で述べた方法で絶縁膜305の予測寿命tNを推定することにより、イオン密度分布の温度依存性を考慮して絶縁膜の寿命を予測することが可能である。   According to this embodiment, the temperature dependence of the ion density distribution inside the insulating film 305 of the sample 201 can be obtained. Therefore, the ion behavior analysis system in the insulating film according to the present embodiment estimates the predicted lifetime tN of the insulating film 305 by the method described in the second embodiment, so that the temperature dependence of the ion density distribution is taken into consideration. It is possible to predict the lifetime.

104…時刻t1のイオン密度、109…時刻t3のイオン密度、112…時刻t2のイオン密度、113…時刻t1の電界強度、114…時刻t2の電界強度、115…時刻t3の電界強度、201…試料、304…アノード、305…絶縁膜、306…カソード、401…試料台、402…パルス電圧発生器、403…圧延素子、404…デジタルオシロスコープ、405…パソコン(PC)、406…画像処理・計算用のPC、407…モニタ、408…パルス電圧発生器の電源、409…絶縁寿命が最も短い計測領域、410…直流バイアス電源、411…パルス電圧波形、901…分割された計測領域、902…試料板稼動装置。   104 ... ion density at time t1, 109 ... ion density at time t3, 112 ... ion density at time t2, 113 ... electric field intensity at time t1, 114 ... electric field intensity at time t2, 115 ... electric field intensity at time t3, 201 ... Sample 304 ... Anode, 305 ... Insulating film, 306 ... Cathode, 401 ... Sample stand, 402 ... Pulse voltage generator, 403 ... Rolling element, 404 ... Digital oscilloscope, 405 ... PC (PC), 406 ... Image processing / calculation PC, 407 ... monitor, 408 ... pulse voltage generator power supply, 409 ... measurement region with the shortest insulation life, 410 ... DC bias power supply, 411 ... pulse voltage waveform, 901 ... divided measurement region, 902 ... sample Plate operating device.

Claims (7)

パルス電圧発生器と圧延素子と演算装置と表示装置を備え、前記パルス電圧発生器が印加したパルス電圧により絶縁膜内に発生した超音波を前記圧延素子で検出することにより、前記絶縁膜内のイオン密度分布を求める絶縁膜内イオン挙動解析システムにおいて、
前記演算装置は、少なくとも異なる2つの時刻における前記絶縁膜内のイオン密度分布を前記超音波の検出結果から求め、前記少なくとも異なる2つの時刻における前記イオン密度分布を基に移流拡散計算を行って、任意の時刻における前記絶縁膜内のイオン密度分布を求め、
前記表示装置は、前記任意の時刻における前記絶縁膜内のイオン密度分布を表示する、
ことを特徴とする絶縁膜内イオン挙動解析システム。
A pulse voltage generator, a rolling element, an arithmetic unit, and a display device are provided, and the ultrasonic wave generated in the insulating film by the pulse voltage applied by the pulse voltage generator is detected by the rolling element. In the ion behavior analysis system in the insulating film to obtain the ion density distribution,
The arithmetic device obtains ion density distribution in the insulating film at at least two different times from the detection result of the ultrasonic wave, performs advection diffusion calculation based on the ion density distribution at the at least two different times, Obtain the ion density distribution in the insulating film at an arbitrary time,
The display device displays an ion density distribution in the insulating film at the arbitrary time;
Insulation film ion behavior analysis system characterized by this.
請求項1記載の絶縁膜内イオン挙動解析システムにおいて、
前記演算装置は、少なくとも異なる2つの時刻における前記絶縁膜内のイオン密度分布から、複数の前記任意の時刻における前記絶縁膜内のイオン密度分布を求め、複数の前記任意の時刻における前記絶縁膜内のイオン密度分布から、複数の前記任意の時刻における前記絶縁膜内の電界分布を求め、前記絶縁膜内の電界分布の時間変化から前記絶縁膜の寿命を予測し、
前記表示装置は、前記演算装置が予測した前記絶縁膜の寿命を表示する絶縁膜内イオン挙動解析システム。
In the insulating film ion behavior analysis system according to claim 1,
The arithmetic unit obtains ion density distributions in the insulating film at the plurality of arbitrary times from ion density distributions in the insulating film at at least two different times, and calculates the inside of the insulating film at the plurality of arbitrary times. Obtaining the electric field distribution in the insulating film at a plurality of arbitrary times from the ion density distribution of, and predicting the lifetime of the insulating film from the time change of the electric field distribution in the insulating film,
The said display apparatus is an ion behavior analysis system in an insulating film which displays the lifetime of the said insulating film which the said arithmetic unit estimated.
請求項1または2記載の絶縁膜内イオン挙動解析システムにおいて、
電源装置を用いて前記絶縁膜に電圧または電流を印加しながら前記絶縁膜にパルス電圧を印加して発生した超音波を前記圧延素子で検出して、前記任意の時刻における前記絶縁膜内のイオン密度分布を求める絶縁膜内イオン挙動解析システム。
In the insulating film ion behavior analysis system according to claim 1 or 2,
An ultrasonic wave generated by applying a pulse voltage to the insulating film while applying voltage or current to the insulating film using a power supply device is detected by the rolling element, and ions in the insulating film at the arbitrary time are detected. Ion behavior analysis system for insulating film to obtain density distribution.
請求項3記載の絶縁膜内イオン挙動解析システムにおいて、
前記演算装置は、前記少なくとも異なる2つの時刻における前記絶縁膜内の電流の波形と電圧の波形と前記イオン密度分布とを基に移流拡散計算を行って、前記任意の時刻における前記絶縁膜内のイオン密度分布と電流分布を求め、
前記表示装置は、前記任意の時刻における前記絶縁膜内のイオン密度分布と電流分布を表示する絶縁膜内イオン挙動解析システム。
In the insulating film ion behavior analysis system according to claim 3,
The arithmetic unit performs advection diffusion calculation based on the current waveform, the voltage waveform, and the ion density distribution in the insulating film at the two different times, and Find the ion density distribution and current distribution,
The display device is an ion behavior analysis system in an insulating film that displays an ion density distribution and a current distribution in the insulating film at the arbitrary time.
請求項1または2記載の絶縁膜内イオン挙動解析システムにおいて、
前記絶縁膜内に設定された複数の領域のそれぞれでのイオン密度分布と電界分布を求めるために、前記絶縁膜または前記パルス電圧発生器と前記圧延素子とを移動させる稼動装置を備え、
前記演算装置は、前記任意の時刻における前記複数の領域のそれぞれでのイオン密度分布と電界分布を求め、
前記表示装置は、前記任意の時刻における前記複数の領域のそれぞれでのイオン密度分布と電界分布を表示する絶縁膜内イオン挙動解析システム。
In the insulating film ion behavior analysis system according to claim 1 or 2,
In order to obtain an ion density distribution and an electric field distribution in each of a plurality of regions set in the insulating film, an operating device for moving the insulating film or the pulse voltage generator and the rolling element is provided.
The arithmetic unit obtains an ion density distribution and an electric field distribution in each of the plurality of regions at the arbitrary time,
The display device is an ion behavior analysis system in an insulating film that displays an ion density distribution and an electric field distribution in each of the plurality of regions at the arbitrary time.
請求項1または2記載の絶縁膜内イオン挙動解析システムにおいて、
前記絶縁膜の移流係数、拡散係数、及び絶縁寿命定数のうち少なくとも1つを前記移流拡散計算に使う移流拡散方程式のパラメータとして用いる絶縁膜内イオン挙動解析システム。
In the insulating film ion behavior analysis system according to claim 1 or 2,
An ion behavior analysis system in an insulating film that uses at least one of the advection coefficient, diffusion coefficient, and insulation lifetime constant of the insulating film as a parameter of an advection diffusion equation used for the advection diffusion calculation.
請求項1または2記載の絶縁膜内イオン挙動解析システムにおいて、
前記演算装置は、記憶装置を備え、計測データと計算データを時系列的に前記記憶装置に格納する絶縁膜内イオン挙動解析システム。
In the insulating film ion behavior analysis system according to claim 1 or 2,
The said arithmetic unit is provided with a memory | storage device, The ion behavior analysis system in an insulating film which stores measurement data and calculation data in the said memory | storage device in time series.
JP2011016016A 2011-01-28 2011-01-28 System for analyzing ion behavior in insulating film Pending JP2012154868A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2011016016A JP2012154868A (en) 2011-01-28 2011-01-28 System for analyzing ion behavior in insulating film
PCT/JP2012/051602 WO2012102321A1 (en) 2011-01-28 2012-01-26 System for analyzing behavior of ions in insulation film

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011016016A JP2012154868A (en) 2011-01-28 2011-01-28 System for analyzing ion behavior in insulating film

Publications (1)

Publication Number Publication Date
JP2012154868A true JP2012154868A (en) 2012-08-16

Family

ID=46580885

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011016016A Pending JP2012154868A (en) 2011-01-28 2011-01-28 System for analyzing ion behavior in insulating film

Country Status (2)

Country Link
JP (1) JP2012154868A (en)
WO (1) WO2012102321A1 (en)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05216860A (en) * 1992-02-04 1993-08-27 Hitachi Ltd Simulation system
JP2008232906A (en) * 2007-03-22 2008-10-02 Ngk Insulators Ltd Method for measuring characteristic of piezo-electric element
JP4599439B2 (en) * 2008-08-07 2010-12-15 株式会社日立製作所 Control device deterioration diagnosis system
JP5416938B2 (en) * 2008-09-09 2014-02-12 日東電工株式会社 Space charge distribution measuring apparatus, space charge distribution measuring method using the apparatus, and high temperature insulating material measured using the method

Also Published As

Publication number Publication date
WO2012102321A1 (en) 2012-08-02

Similar Documents

Publication Publication Date Title
US7005851B2 (en) Methods and apparatus for inspection utilizing pulsed eddy current
JP6616302B2 (en) Lock-in thermography method and system for hot spot location determination
US8618810B2 (en) Identifying fuel cell defects
JP5701961B2 (en) Electrical conductivity measurement method and electrical conductivity measurement system using the same
Saguy et al. Flaw detection in metals by the ACPD technique: Theory and experiments
JP2015094717A (en) Thermal fatigue life prediction device, thermal fatigue life prediction method, and program
Liu et al. Investigation of thermal imaging sampling frequency for eddy current pulsed thermography
EP2807496A1 (en) Insulation defect detection of high voltage generator stator core
CN112180177B (en) Power frequency electromagnetic field evaluation method and system fusing measured data
JP2015230171A (en) Temperature measurement method of lightning protection element in arrester constituted of lightning protection element and porcelain tube covering the same
JP2013019841A (en) Defect evaluation method for structure
WO2012102321A1 (en) System for analyzing behavior of ions in insulation film
JP2018115874A (en) Inspection device, inspection method, inspection program, storage medium and inspection system
Yusa et al. Four-terminal measurement of the distribution of electrical resistance across stress corrosion cracking
Cheng et al. Development of a flexible capacitive sensor for concrete structure health monitoring
JP7237516B2 (en) Deterioration estimation device, deterioration estimation system, deterioration estimation method, and computer program
JP2008205230A (en) Mos semiconductor device of trench structure, lifetime evaluation device, and lifetime evaluation method
Yeh et al. Numerical simulations of electric potential field for alternating current potential drop associated with surface cracks in low-alloy steel nuclear material
JP2009180625A (en) Coil inspection apparatus and coil inspection method
JP6899644B2 (en) Damage inspection system, damage inspection method, and program
JP2010127698A (en) Apparatus and method for estimation of radiated electromagnetic wave frequency
Bekbaev et al. On the possibilities of dynamic evaluation of contact surface temperature under impulse-current loads
Majidnia et al. A method of determining the eddy current instrument response from models
Abdul Wahab et al. Initial Study of Multiple Excitation Source for Electrical Resistance Tomography in Steel Pipe Application
Bussolino et al. Virtual Experiments by Pulse Heating Techniques: Cylindrical Tungsten Specimens