JP3937149B2 - 分光エリプソメータを用いた極薄膜２層構造の解析方法 - Google Patents
分光エリプソメータを用いた極薄膜２層構造の解析方法 Download PDFInfo
- Publication number
- JP3937149B2 JP3937149B2 JP2002110710A JP2002110710A JP3937149B2 JP 3937149 B2 JP3937149 B2 JP 3937149B2 JP 2002110710 A JP2002110710 A JP 2002110710A JP 2002110710 A JP2002110710 A JP 2002110710A JP 3937149 B2 JP3937149 B2 JP 3937149B2
- Authority
- JP
- Japan
- Prior art keywords
- analysis
- layer
- stage
- model
- film thickness
- 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.)
- Expired - Lifetime
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Spectrometry And Color Measurement (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、分光エリプソメータを用いて得たデータを、広範囲最極小値計算法（Extended Best Local Minimum Calculation 以下、Ｅｘｔｅｎｄｅｄ ＢＬＭＣ）を用いて処理する極薄膜２層構造の解析方法に関する。
【０００２】
【従来の技術】
分光エリプソメータを用いて入射光と反射光の偏光変化量を測定し、その結果から膜厚（ｄ）、複素屈折率Ｎ（Ｎ＝ｎ−ｉｋ）を算出することができる。
偏光変化量（ρ）はρ＝ｔａｎψｅｘｐ（ｉΔ）で表され、波長（λ）、入射角度（φ）、膜厚、複素屈折率等のパラメータに依存するので、その関係は次のようになる。
（ｄ，ｎ，ｋ）＝ｆ（Ψ，Δ，λ，φ）
【０００３】
入射角度を固定した場合、単一波長エリプソメータでは、（ｄ，ｎ，ｋ）の３つの未知数に対し、２つの独立変数しか測定できないので、ｄ，ｎ，ｋの内のいずれか一つを既知として固定する必要がある。
単一波長でも入射角度を変えると測定変数は増加する。しかしながら、入射角度（φ）の違いによる（Ψφ1 ，Δφ1 ）と（Ψφ2 ，Δφ2 ）に強い相関関係があるため、ｄ，ｎ，ｋを精度良く求めることは難しい。
【０００４】
【発明が解決しようとする課題】
偏光変化量は、光が通る体積、（位相角（β）×ビーム径の面積）で表すことができる。ビーム径が一定とすると偏光変化量は次のようになる。
偏光変化量∝膜厚（ｄ）×（複素屈折率Ｎ）×φ
ここにおいて、φは入射角である。膜厚（ｄ）が薄く、Ｎが低ければ、位相角（β）の変化が小さくなり、ｄとＮの間に相関関係が強くなる。
【０００５】
前述した問題は、極薄膜２層構造の場合、上記にさらに各層の間に相関関係が発生する可能性があり、偏光変化量を表している測定結果（ΨE ( λi ) ，ΔE ( λi ) ）から、各層のｎ，ｋ，ｄを求めることが困難となる。
【０００６】
前述したように、測定された多層薄膜の偏光変化を表している測定スペクトル（ΨE ( λi ) ，ΔE ( λi ) ）は、前記基板のｎ，ｋ情報、各層のｎ，ｋ，ｄの情報のすべてを含んでいるが、これから、前記基板のｎ，ｋ情報、各層のｎ，ｋ，ｄの情報の唯一の組み合わせを算出することはできない（基板のみの場合を除く）。この前記唯一の組み合わせを探す方法は、分光エリプソメータデータの解析と呼ばれる。解析を行うときは、前記基板のｎ，ｋ情報、各層のｎ，ｋ，ｄの情報を利用してモデルをたてる。この中の基板や各層のｎ，ｋ情報には、リファレンス（既知のテーブルデータ）や分散式、または似たような材質の単層薄膜の光学定数を使用する。
分散式とは、物質の誘電率の波長依存性を示す式であり、近赤外から紫外線領域では、この誘電率ε（λ）は材料の構成原子の結合様式から決定される。
分散式として、調和振動子をもとにした計算式、量子力学をもとにした計算式、経験式等が知られており、通常２つ以上のパラメータを含んでいる。
前述したモデルに含まれているすべての未知数（各層の厚さや、分散式パラメータ、混合比など）を変化させながら、測定データにあわせていく。
これをフィッティングといい、このフィッティングの結果、各層の膜厚や混合比などが求まり、分散式パラメータからは、材料の誘電率ε（λ）を計算することができる。材料の誘電率と屈折率は下記の関係にある。
ε＝Ｎ2
本発明では、様々な理由により１つの分散式では膜の誘電率の波長依存性を特定できないかまたは困難な場合についても問題にして、有効媒質論（Effective Medium Theory)を利用して、有効誘電率(Effective Dielectric Function) を計算する。一般的に、例えばホスト材料中に様々な誘電率を持つＮ個の物質（ゲスト）が混ざり合っている場合、有効誘電率（ε）は、下記の様に表される。
この時εh はホストの誘電率を、εj はｊ個目の誘電率を、ｋはスクリーニングファクターを示している。ここで、ホストの材料と、中に入っている材料がほぼ同じ量で混ざり合っているか、またはどちらがホストかゲストか分からない場合、ホストの材料自身が有効媒質材料と同じε＝εh になる。この条件の有効媒質論をブラッグマンの有効媒質近似（Bruggeman Effective Medium Approximation以下本願において単にＥＭＡ）と呼ばれている。３つの球状物質Ａ，Ｂ，Ｃが対照的に混在しているときの誘電率εは次の式で与えられる。
ｆa （εa −ε）／（εa ＋２ε）＋ｆb （εb −ε）／（εb ＋２ε）＋ｆc （εc −ε）／（εc ＋２ε）＝０
ここで、
ε : 求めようとする有効誘電率
εa , εb , εc : 球状物質Ａ，Ｂ，Ｃの誘電率
ｆa , ｆb , ｆc : 各物質の混合比（Volume Fraction 以下Vf）で、
ｆa ＋ｆb ＋ｆc ＝1
基板上の膜が不均一または不連続や、いくつかの材料が混ざり合っている場合、波長オーダーより十分小さく、物理的に混合している複数の物質から成る媒質については、有効媒質近似（ＥＭＡ）を利用してモデルをたてる。
物質Ａと物質Ｂと物質Ｃが混合している場合について説明する。このときの有効媒質近似（ＥＭＡ）は、物質Ａの混合比、物質Ｂの混合比、物質Ｃの混合比、ＡとＢとＣの混合層の膜厚、誘電率には分散式やリファレンスデータなどを推定してフィッティングを行い、評価する。
【０００７】
本発明の目的は、膜厚や複素屈折率などの組み合わせモデルを設定し、そのシミュレーションスペクトルを算出して、そのシミュレーションスペクトルと測定スペクトルとのフィッティングを広範囲最極小値計算法（Ｅｘｔｅｎｄｅｄ ＢＬＭＣ）を使用して行うことにより、極薄膜２層構造を決定する、分光エリプソメータを用いた極薄膜２層構造解析方法を提供することにある。
【０００８】
【課題を解決するための手段】
前記目的を達成するために、本発明による極薄膜２層構造を測定し解析する方法は基本的に、
まず、分光エリプソメータを用いて極薄膜２層構造を測定し、分光スペクトルを得る。
新しい解析方法は基本的には３つの段階から成立している。
解析第１段階の目的は、実際のサンプルに良く合うと思われるモデルを複数選択して初期値を決定する。
解析第２段階では、第１段階で得られた初期値をもとにＥｘｔｅｎｄｅｄ ＢＬＭＣを行う。
必要に応じて、解析第３段階では最後のフィッティング、確認、保存を行う。
【０００９】
本発明で利用する超薄膜および薄膜計測方法で用いる極小値計算法（ＢＬＭＣ）は、
極小値計算法（ＢＬＭＣ）による分光エリプソメータを用いた計測対象の基板表面の超薄膜および薄膜計測方法において、
計測対象の基板表面の薄膜を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE (λi )とΔE (λi )を得るΨE ，ΔE スペクトル測定ステップと、
前記基板の（Ｎ 0 （ｎ0 ,ｋ0 ）），基板上の薄膜の（ｄ,Ｎ（ｎ,ｋ））を分散式を用いて仮定し、さらに予想される範囲内にある複数の膜厚（ｄ±ｍΔｄ）および予想される範囲内にある複数の入射角（φ±ｍΔφ）を設定するステップと、
前記入射角と膜厚の組み合わせにもとづいて、分散式（ＤＳＰ）のパラメータ（εs ,ωt ）のフィッティングを行うステップと、
前記フィッティングにより得られた各ΨM (λi )とΔM (λi )の中から前記ΨE (λi )とΔE (λi )との差の最も少なくなる膜厚（ｄbest）と入射角（φbest）の組み合わせを設定したモデルのフィッティング結果（ＤＳＰbest）を選択するステップと、
前記選択するステップで選択された入射角（φbest）を確定値として、膜厚（ｄbest）と分散式（ＤＳＰbest）のフィッティングを行うステップで構成されている。
【００１０】
前記目的を達成するために本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
計測対象の基板上の極薄膜２層構造を、分光エリプソメータにより入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE(λi)とΔE (λi )を得るΨE ，ΔE スペクトル測定段階と、極薄膜２層構造の基板の複素屈折率（Ｎ0（ｎ0，ｋ0 ））および各薄膜の材料（Ｍａｔ1，Ｍａｔ2）の考えられる光学定数に係る複素屈折率（Ｎ1（ｎ1，ｋ1 ）），（Ｎ2 （ｎ2 ，ｋ2 ））、各層の膜厚（ｄ1，ｄ2）を利用し、複数のモデルをたてる第１ステップ、
前記各モデルごとに前記測定スペクトルとのフィッティングを行って平均二乗誤差（χ2 ）を得る第２ステップ、および
前記各モデルごとのフィッティングの結果、最低の平均二乗誤差（χ2 ）の 値をもつモデルまたは、あらかじめ設定した膜厚の最大、最小値の中に入っている平均二乗誤差（χ2 ）の値が最も低いモデルを選択し、膜厚（ｄ1(best) ，ｄ2(best) ）を含むそのモデルに係る各層の値を得る第３ステップ
から成る解析第１段階と、
解析第１段階で得られた値を、前記解析第１段階後の新しいモデルの初期値として設定する第１ステップ、
前記解析第１段階の第３ステップにて得られたいずれか一方の層の膜厚（ｄ1(best) またはｄ2(best) ）を中心値として、そのまわりの複数点ごとに、もう一方の層の膜厚（ｄ2(best)またはｄ1(best)）を変化させ最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いてフィッティングを行う第２ステップ、および
前記複数点ごとの最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いて行ったフィッティングの結果から最低の平均二乗誤差（χ2）の値をもつモデルまたは、あらかじめ設定した膜厚、分散式パラメータ、入射角のそれぞれ最大、最小値の中に入っている平均二乗誤差（χ2 ）の値が最も低いモデルを選択する第３ステップ
から成る解析第２段階と、
解析第２段階で選択されたモデルの複素屈折率及び膜厚を利用して最終的なフィッティングを行う第１ステップ、および
前記最終的なフィッティングで得られた結果の確認を行う第２ステップ
から成る解析第３段階と
から構成されている。
【００１１】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第２段階の第２ステップにおいて、前記いずれか一方の層は、２層構造中の材料において、光学定数がより分からない方の層とするものである。
【００１２】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第２段階の第２ステップにおいて、光学定数がより分かっている方の層の膜厚（ｄ1(best) またはｄ2(best)）を中心値として、その中心値の±１０％までの範囲の中で変化させた膜厚ごとに、もう一方の層について最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を行う。
また、本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第１段階の第１ステップは、１、２層どちらかの層が、複数の材料が混ざり合っている場合に、各材料の混合比（Vf1 ，Vf2）を利用し、有効媒質近似を用いてモデルをたてるものである。
【００１３】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第３段階の第１ステップにおいて、前記解析第２段階の第３ステップで選択したモデルの１、２層の光学定数を固定した１、２層目の膜厚のフィッティングまたは、モデルの１層の光学定数を固定した１、２層目の膜厚および２層目の光学定数のフィッティングを行うものである。
【００１４】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第３段階の第２ステップにおいて、フィッティングで得られた結果が、あらかじめ設定した膜厚と分散式パラメータでそれぞれの最大、最小値の間に入っているか否かを確認し、
入っていない場合は前記解析第１段階の第１ステップに戻り、極薄膜２層構造の基板の複素屈折率（Ｎ 0 （ｎ 0 ，ｋ 0 ））および各薄膜の材料（Ｍａｔ 1 ，Ｍａｔ 2 ）の考えられる光学定数に係る複素屈折率（Ｎ 1 （ｎ 1 ，ｋ 1 ）），（Ｎ 2 （ｎ 2 ，ｋ 2 ））、各層の膜厚（ｄ 1 ，ｄ 2 ）を利用し、複数のモデルをたてて前記解析第１段階の第２ステップ以降を行うものである。
【００１５】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
計測対象の基板上の極薄膜２層構造を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE(λi)とΔE (λi )を得るΨE ，ΔE スペクトル測定段階と、
１，２層どちらかの層が、複数の材料が混ざり合っている場合、有効媒質近似を用いてモデルをたて、極薄膜２層構造の基板の複素屈折率（Ｎ0（ｎ0 ，ｋ0 ））および各薄膜の材料（Ｍａｔ1 ，Ｍａｔ2 ）の考えられる複素屈折率（Ｎ1（ｎ1 ，ｋ1 ）），（Ｎ2（ｎ2 ，ｋ2 ））、材料が混ざり合っている場合の各材料の混合比（Vf1，Vf2）、各層の膜厚（ｄ1 ，ｄ2 ）を利用し、複数のモデルをたてる第１ステップ、
前記各モデルごとに前記測定スペクトルとのフィッティングを行う第２ステップ、および
前記各モデルごとのフィッティングの結果、最低の平均二乗誤差（χ2 ）の値をもつモデルまたはあらかじめ設定した膜厚と混合比のそれぞれの最大、最小値の中に入っている平均二乗誤差（χ2 ）の値の最も低いモデルに係る、膜厚（ｄ1(best) ，ｄ2(best) ），混合比Vf(best)を含む各層の値を得る第３ステップ
から成る解析第１段階と、
前記解析第１段階の第３ステップで得られた膜厚・混合比の値を、前記解析第１段階後の新しいモデルの初期値とし、未知の分散式が入っている他の層の前記解析第１段階の第３ステップで得られた膜厚の値を基に、予想される範囲内にある膜厚（（ｄ1±ｍΔｄ1 ）または（ｄ2 ±ｍΔｄ2 ））を設定し、もう一方の層の膜厚についても、前記解析第１段階で得られた膜厚の値を中心として、そのまわりの複数点（（ｄ2±ｍΔｄ2 ）または（ｄ1 ±ｍΔｄ1 ））を設定し、また、前記解析第１段階の第３ステップで得られた混合比の値を中心値として、そのまわりの複数点（Vf±ｍΔVf）を設定する第１ステップ、
混合比に係る複数点（Vf±ｍΔVf）と未知の分散式が入っている他の層のもう一方の層の膜厚に係る複数点（（ｄ2 ±ｍΔｄ2 ）または（ｄ1±ｍΔｄ1 ））との組み合わせのなかで、未知分散式が入っている他の層について最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いてフィッティングを行う第２ステップ、および
混合比と未知の分散式が入っている他の層のもう一方の層の膜厚との組み合わせごとに得られた結果の中から、最低の平均二乗誤差（χ2 ）の値またはあらかじめ設定した膜厚、分散式パラメータ、混合比、入射角のそれぞれ最大、最小値の間に入っている最低の平均二乗誤差（χ2 ）の値を持つ組み合わせのモデルを選択する第３ステップ
から成る解析第２段階と、
前記解析第２段階の第３ステップで選択されたモデルの複素屈折率及び膜厚を基に、両膜厚、混合比および、分散式のフィッティングまたは、両膜厚、混合比のフィッティングを行う第１ステップ、および
該第１ステップのフィッティングで得られた結果の確認を行う第２ステップ
から成る解析第３段階と
から構成されている。
【００１６】
【発明の実施の形態】
以下図面等を参照して本発明による分光エリプソメータを用いた極薄膜２層構造データの解析方法の実施の形態を説明する。
図１は、本発明方法で使用するエリプソメータの構成を示す略図である。
なお測定の対象の試料 (サンプル) ４の一部を拡大して示してある。
【００１７】
図１に示されている分光エリプソメータにより、後述する方法の分光測定データの獲得ステップが実行される。先ず、測定装置について簡単に説明する。Ｘｅランプ１は、多数の波長成分を含む、いわゆる白色光源である。このＸｅランプ１の発光は光ファイバ２を介して偏光子３に導かれる。偏光子３により偏光された光は、測定対象であるサンプル４の表面に特定の入射角（例えばφ＝７５．００°）で入射させられる。サンプル４は後述する基板表面に２層に薄膜を形成した測定試料である。
サンプル４からの反射は、光弾性変調器（ＰＥＭ）５を介して検光子６に導かれる。光弾性変調器（ＰＥＭ）５により５０ｋＨｚの周波数に位相変調されて、直線から楕円偏光までが作られる。そのため、数ｍ秒の分解能でΨ，Δを決定することができる。検光子６の出力は光ファイバ７を介して分光器８に接続される。分光器８の出力データがデータ取込部９に取り込まれ、分光測定データの獲得ステップを終了する。なお、ＰＥＭ５の位置は偏光子３の後か検光子６の前どちらでも可能とする。
【００１８】
図２は、本発明による薄膜計測方法の実施形態を示す流れ図である。
図３は、前記実施形態における解析第１段階を示す詳細説明図である。
図４は、前記実施形態における解析第２段階を示す詳細説明図である。
図５は、前記実施形態における解析第３段階を示す詳細説明図である。
次の段落に本発明で用いる記号を一括して示す。
【００１９】
（記号の定義）
Ｓｕｂ ： 基板（物理定数は既知，バルクとして取り扱い可能）
Ｍａｔ : 薄膜材料（物質の光学定数）
Ｍａｔij： j 番目のモデルのｉ 層目の材料（同上）
ｄ i : i 層目の膜厚
ｄ i (best): フィッティングによって得られたi 層目の膜厚
ｄij： j 番目のモデルのｉ 層目の膜厚
ｄij(best)： フィッティングで得たj 番目のモデルのｉ 層目の膜厚
X2 : 平均二乗誤差（χ2 ）の値
X2(j) : j 番目のモデルにおける平均二乗誤差（χ2 ）の値
Void : ｎ＝１，ｋ＝０の物質
Vf(ij) : j 番目のモデルのｉ 層目の混合比(Volume fraction)
Vf(ij)(best) : フィッティングで得られた j 番目のモデルのｉ 層目の混合比
【００２０】
（フィッティングの結果のχ2 の意味）
Ｎ個の測定データ対Ｅｘｐ（ｉ＝１，２．．．，Ｎ）と前記モデルの対応するＮ個のモデルの計算データ対Ｍｏｄ（ｉ＝１，２．．．，Ｎ）とし、測定誤差は正規分布をするとし、標準偏差をσi とすると、平均二乗誤差（χ2 ）は、次のようにして与えられる。
ここで、Ｐはパラメータの数である。χ2 が小さいということは測定結果とモデルの一致度が大きいということに他ならないから、複数のモデルについて比較するときに、χ2 の最も低いものがベストモデルということになる。
【００２１】
（測定段階データ測定ステップ）
測定は図１に示す装置で測定を行う。計測対象４の基板表面の２層構造の薄膜（図中の拡大図を参照されたい）を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE ( λi ) とΔE ( λi ) を得るΨE ，ΔE スペクトルを測定する。
（測定段階データ保存ステップ）
前のステップで測定したデータを保存して比較の対象データとする（図２参照）。
【００２２】
（解析第１段階ステップ１）
このステップで、基板の（Ｎ 0 （ｎ0 ,ｋ0 ））および各薄膜の材料（Ｍａｔ1 ，Ｍａｔ2 ）の考えられる複素屈折率（Ｎ 1 （ｎ1 ,ｋ1 ）），（Ｎ 2 （ｎ2 ,ｋ2 ）），各層の膜厚（ｄ1 ,ｄ2 ）を利用し、いくつかのモデルをたてる。
この実施例では、下記各モデル(1) 〜モデル(4) を設定したものとする。
図３の３１が示す部分に前記各モデルを略図示してある。
モデル(1) は、基板（Ｓｕｂ）の上に第１層(Ｍａｔ1 の光学定数，膜厚ｄ11)と第２層（Ｍａｔ2 の光学定数，膜厚ｄ21）を形成したものである。
モデル(2) は、基板（Ｓｕｂ）の上に第１層(Ｍａｔ1 の光学定数，膜厚ｄ12)と第２層(Ｍａｔ2 の光学定数＋Ｖｏｉｄ、膜厚ｄ22)を形成したものである。なおＶｏｉｄとは屈折率１の空間である。
モデル(3) は、基板（Ｓｕｂ）の上に第１層(Ｍａｔ1 の光学定数，膜厚ｄ13)と第２層(Ｍａｔ2 の光学定数＋Ｍａｔ3 の光学定数、膜厚ｄ23)を形成したものである。なお、（Ｍａｔ2 ＋Ｍａｔ3 ）は、材料２と材料３をある比率で混合したものを意味する。
モデル(4) は、基板（Ｓｕｂ）の上に第１層（Ｍａｔ1 の光学定数，膜厚ｄ14）と第２層（Ｍａｔ2 の光学定数の基準にする量(Special reference) ＋Ｖｏｉｄ、膜厚ｄ24)を形成したものである。なお、基準にする量とは似たような試料を用いて得られた光学定数である。
ここで、モデル(2) 〜(4) の２層目で設定されている混合比からは、有効媒質近似論を用いて、均質膜としての光学定数を求めることができる。なお前記各モデルは、基板Ｓｕｂの材料の光学定数は既知、第１層の光学定数は、ほぼ既知で、第２層の光学定数および、第１層の膜厚ｄ1 、第２層の膜厚ｄ2 は、未知（不確か）であることを前提にして前述の４個のモデルをたててある。
【００２３】
（解析第１段階ステップ２）
前記ステップ１で選定した４個のモデル(1) 〜(4) のそれぞれについて、前記測定スペクトルより得られた測定データΨE ，ΔE とのフィッティングを行う。
図３の３２の示す部分に各モデルのフィッティングの対象とフィッティングの結果得られたデータのχ2 値を示してある。
モデル(1) では、第１層目の膜厚ｄ11と第２層目の膜厚ｄ21をフィッティングして、その結果ｄ11(best)，ｄ21(best)とχ2 値 X2 (1) を得る。
モデル(2) では、第１層目の膜厚ｄ12，第２層目の膜厚ｄ22と第２層目の混合比のフィッティングをして、その結果ｄ12(best)，ｄ22(best)，Vf22(best)およびχ2 値 X2 (2) を得る。
モデル(3) では、第１層目の膜厚ｄ13，第２層目の膜厚ｄ23と第２層目の混合比のフィッティングをして、その結果ｄ13(best)，ｄ23(best)，Vf23(best)およびχ2 値 X2 (3) を得る。
モデル(4) では、第１層目の膜厚ｄ14，第２層目の膜厚ｄ24と第２層目の混合比のフィッティングをして、その結果ｄ14(best)，ｄ24(best)，Vf24(best)およびχ2 値 X2 (4) を得る。
（解析第１段階ステップ３）
前記複数組のフィッティングの結果から、最低χ2 値またはあらかじめ設定した膜厚の最大、最小値の間に入っている最低χ2 値のモデルの結果を選択するステップであり、これを図３の３３に示す。
【００２４】
（解析第２段階ステップ１）
図４の４１に解析第２段階で使用するモデル（初期値）が示されている。
この例では、Ｍａｔ1 の光学定数は、ほぼ既知となっているので、解析第１段階で使用した光学定数（Ｍａｔ1 ）をそのまま利用している。２層目の材料の光学定数（Ｍａｔ2 ）は未知のため、ここでは分散式を用いている。なお、膜厚（ｄ）の初期値は、解析第１段階ステップ３で得られた値を使用する。
（解析第２段階ステップ２）
このステップでは、図４の４２に示すように、前記モデル（初期値）の１層目膜厚（ｄ1(best) ）を中心値としてｄ1(best) ＋１０％，ｄ1(best) ＋ ５％，ｄ1(best) ，ｄ1(best) −５％，ｄ1(best) −１０％上下に変化させる。そして前記複数点（５点）に対して２層目の厚さ（ｄ2(best) ）を±１０％の範囲で変化させてＢＬＭＣを行う。
各モデルの２層目の厚さｄ2jと光学定数Ｍａｔ2jとχ2 値 X2 (j) を得る。
（解析第２段階ステップ３）
ステップ２で行われた結果から最低χ2 値またはあらかじめ設定した膜厚と分散式パラメータでそれぞれの最大、最小値の間に入っている最低χ2 値のモデルを選択する（図３、４３）。
【００２５】
（解析第３段階ステップ１）
解析第２段階ステップ３で選択したモデルの１，２層の光学定数を固定して１，２層目の膜厚のフィッティングを行う。または
１層の光学定数を固定して１，２層目の膜厚および２層目の光学定数をフィッティングする（図５、５１参照）。
（解析第３段階ステップ２）
前記ステップ１の結果があらかじめ設定した膜厚と分散式パラメータでそれぞれの最大、最小値の間に入っている最低χ2 値かどうかの確認をする（図５、５２参照）。なお、結果の確認が妥当でないときは解析第１段階に戻り、新しいモデルを設定して、フィッティング（ステップ２）を行う。
（解析第３段階ステップ３）
前記ステップ２で結果が妥当であるときは保存する（図５、５３参照）。
【００２６】
【実施例】
次に、前記同様に、１層目の光学定数は略既知、２層目の光学定数および１，２層目の膜厚が、未知の場合の実施例について説明する。図２に示した流れ図をそのまま使用することができる。この実施例は、基板ＳｕｂがSiで、１層目の材料をSiO2, ２層目の材料を SiNX としたものである。
前記試料を図１に示す装置で測定を行う。計測対象４の基板表面の２層構造の薄膜を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE ( λi ) とΔE ( λi ) を得るΨE ，ΔE スペクトルを測定して測定したデータを保存して比較の対象データとする。
【００２７】
図３の解析第１段階で、各薄膜の材料（SiO2 ， SiNx ）の考えられる複素屈折率（Ｎ 1 （ｎ1 ,ｋ1 ）），（Ｎ 2 （ｎ2 ,ｋ2 ））と膜厚（ｄ1 ,ｄ2 ）を利用し、モデルを準備する。この実施例では、下記各モデル(1) 〜モデル(4) を設定したものとする。
【００２８】
（解析第１段階ステップ１）
モデル(1) は、基板（Ｓｕｂ：Si）の上に第１層( SiO2の光学定数、膜厚ｄ11) と第２層（Si3N4 の光学定数，膜厚ｄ21）を形成したものである。
モデル(2) は、基板（Ｓｕｂ：Si）の上に第１層( SiO2の光学定数、膜厚ｄ12) と第２層( Si3N4 の光学定数＋Ｖｏｉｄ、膜厚ｄ22 )を形成したものである。なおＶｏｉｄとは屈折率１の空間である。
モデル(3) は、基板（Ｓｕｂ：Si）の上に第１層( SiO2の光学定数、膜厚ｄ13) と第２層( Si3N4 の光学定数＋ SiNX の光学定数、膜厚ｄ23) を形成したものである。
なお、（Si3N4 ＋SiNx）は、Si3N4 とSiNxをある比率で混合したものを意味する。
モデル(4) は、基板（Ｓｕｂ：Si）の上に第１層（SiO2の光学定数、膜厚ｄ14）と第２層（ SiNx の光学定数（既知の基準量）＋Ｖｏｉｄ、膜厚ｄ24 )を形成したものである。ここで、モデル(2) 〜(4) の２層目で設定されている混合比からは、有効媒質近似論を用いて、均質膜としての光学定数を求めることができる。なお前記各モデルは、基板Ｓｕｂの材料の光学定数は既知、第１層のSiO2の光学定数は、ほぼ既知で、第２層のSiNxの光学定数および、第１層の膜厚ｄ1,第２層の膜厚ｄ2 は、未知（不確か）であることを前提にしている。
【００２９】
（解析第１段階ステップ２）
前記ステップ１で選定した４個のモデル(1) 〜(4) のそれぞれについて、前記測定スペクトルより得られた測定データΨE ，ΔE とのフィッティングを行う。
図３の３１，３２の示す部分に各モデルのフィッティングの対象とフィッティングの結果得られたデータのχ2 値を示してある。
モデル(1) では、第１層目の膜厚ｄ11と第２層目の膜厚ｄ21をフィッティングしてその結果ｄ11(best)、ｄ21(best)とχ2 値 X2 (1) を得る。
モデル(2) では、第１層目の膜厚ｄ12、第２層目の膜厚ｄ22と第２層目の混合比のフィッティングしてその結果ｄ12(best)、ｄ22(best)、Vf22(best)およびχ2 値 X2 (2) を得る。
モデル(3) では、第１層目の膜厚ｄ13、第２層目の膜厚ｄ23と第２層目の混合比のフィッティングしてその結果ｄ13(best)、ｄ23(best)、Vf23(best)およびχ2 値 X2 (3) を得る。
モデル(4) では、第１層目の膜厚ｄ14、第２層目の膜厚ｄ24と第２層目の混合比のフィッティングしてその結果ｄ14(best)、ｄ24(best)、Vf24(best)およびχ2 値 X2 (4) を得る。
【００３０】
（解析第１段階ステップ３）
前記複数組のフィッティングの結果から、最低χ2 値またはあらかじめ設定した膜厚の最大、最小値の間に入っている最低χ2 値のモデルの結果を選択するステップであり、これを図３の３３に示す。
【００３１】
（解析第２段階ステップ１）〜（解析第３段階ステップ３）は前述したとおりである。
【００３２】
【発明の効果】
以上詳しく説明したように、本発明によれば、基板上に形成された極薄膜２層の厚さと光学定数を、分光エリプソメータを用いて得たデータを、広範囲最極小値計算法（Ｅｘｔｅｎｄｅｄ ＢＬＭＣ）を用いて処理することにより解析することができる。
前述したように基本的に３段階が含まれている手順は、２層以上の様々な構造の解析にも利用できる。
【００３３】
本発明方法によれば、
１．極薄膜多層構造でも、この方法を用いることで信頼できる膜厚（各層において）および光学定数（少なくとも１層）が得られる。
２．この手順の第１段階を行うことで、未知のパラメータの初期値範囲を狭めることができる。
３．Ｅｘｔｅｎｄｅｄ ＢＬＭＣを使用することで、好ましくないローカルミニマム（Ｌｏｃａｌ Ｍｉｎｉｍｕｍ）に落ちこむことは劇的に減少し、結果の信頼性があがる。
【００３４】
以上詳しく説明した実施例について、本発明の範囲内で種々の変形を施すことができる。理解を容易にするために、データの取得、モデルの設定に関連して、一貫してΨ，Δを用いて説明した。当業者には良く知られている以下のデータ対を用いても同様な、測定およびフィッティングが可能であり、本発明の技術的範囲に含まれるものである。
（ｎ，ｋ） （εr ，εi ） ( tan Ψ，cos Δ) 、 (Ｉs,Ｉc ）
【００３５】
基板上のＭａｔ1 ，Ｍａｔ2 には、極薄膜誘電体材料多層構造だけでなく、様々な厚さや材料など、幅広いアプリケーションにも使用することができる。また、実施例として光弾性変調器（ＰＥＭ）の例を示したが、ＰＥＭ以外のエリプソメータを使用することもできる。
【００３６】
上記手段の全部または１部を行う場合もあり、これも本発明の技術的範囲に含まれるものである。
【００３７】
基板もSiの他、ガラスや石英、化合物半導体なども同様に利用できる。また、基板の種類によらず、どんな平坦な基板でも、あれている基板でも使用することができる。
【００３８】
分散式には、クラシカル（古典力学）、アモルファス（量子力学）、経験式の他、様々な式・パラメータも使用可能であり、これも本発明の技術的範囲に含まれるものとする。
【００３９】
例ではあらゆる全てのパラメータを同時にフィッティングすると説明したが、別々にフィッティングする場合もあり、これも本発明の技術的範囲に含まれるものとする。
【００４０】
Ｅｘｔｅｎｄｅｄ ＢＬＭＣの一部であるＢＬＭＣでは、前述したとおり、入射角度をフィッティングする場合がある。また、手順では入射角度と様々なパラメータを同時にフィッティングするとしたが、別々にフィッティングする場合や、入射角度を固定する場合もあり、これらも全て本発明の技術的範囲に含まれるものとする。
【００４１】
入射角度はＢＬＭＣ以外でも、一般的なパラメータとしてフィッティングすることがあり、これも本発明の技術的範囲に含まれるものとする。
【００４２】
例では、解析第１段階第１ステップのモデル準備過程において、基板はSi、第１層目はSiO2、第２層目は(1) Si3N4 、(2) Si3N4 とＶｏｉｄ、(3) Si3N4 と SiNX 、および(4) SiNx とＶｏｉｄの(1) 〜(4) のモデルとしたが、製造プロセスになどによってモデルの種類、数が変わる場合もあり、これも本発明の技術的範囲に含まれるものとする。
【００４３】
上記のように、２ステップが含まれている手順は、２層以上の様々な構造にも対応可能である。
【００４４】
例ではＥＭＡを使用すると説明したが、他の有効媒質近似論の使用も可能であり、これも本発明の技術的範囲に含まれるものとする。
【図面の簡単な説明】
【図１】 本発明方法で使用するエリプソメータの構成を示す略図であり、測定の対象の試料 (サンプル) ４の一部を拡大して示してある。
【図２】 本発明による薄膜計測方法の実施形態を示す流れ図である。
【図３】 前記実施形態における解析第１段階を示す詳細説明図である。
【図４】 前記実施形態における解析第２段階を示す詳細説明図である。
【図５】 前記実施形態における解析第３段階を示す詳細説明図である。
【符号の説明】
１ Ｘｅランプ
２ 光ファイバ
３ 偏光子
４ サンプル
５ 光弾性変調器（ＰＥＭ）
６ 検光子
７ 光ファイバ
８ 分光器
９ データ取込部
３１ モデル複数組を準備する解析第１段階第１ステップ
３２ フィッティングを行う解析第１段階第２ステップ
３３ フィッティングの結果から、最低の平均二乗誤差（χ 2 ）の値のモデルの結果を選択する解析第１段階第３ステップ
４１ ３３で得られた結果を新しいモデルの初期値として設定する解析第２段階第１ステップ
４２ 一方の層の膜厚複数点ごとに、もう一方の層についてＢＬＭＣを行う解析第２段階第２ステップ
４３ 前記フィッティング結果から最低の平均二乗誤差（χ 2 ）の値のモデルを選択する解析第２段階第３ステップ
５１ 前記選択されたモデルの最終フィッティングを行う解析第３段階第１ステップ
５２ 結果の確認を行う解析第３段階第２ステップ
５３ 保存を行う解析第３段階第３ステップ
【発明の属する技術分野】
本発明は、分光エリプソメータを用いて得たデータを、広範囲最極小値計算法（Extended Best Local Minimum Calculation 以下、Ｅｘｔｅｎｄｅｄ ＢＬＭＣ）を用いて処理する極薄膜２層構造の解析方法に関する。
【０００２】
【従来の技術】
分光エリプソメータを用いて入射光と反射光の偏光変化量を測定し、その結果から膜厚（ｄ）、複素屈折率Ｎ（Ｎ＝ｎ−ｉｋ）を算出することができる。
偏光変化量（ρ）はρ＝ｔａｎψｅｘｐ（ｉΔ）で表され、波長（λ）、入射角度（φ）、膜厚、複素屈折率等のパラメータに依存するので、その関係は次のようになる。
（ｄ，ｎ，ｋ）＝ｆ（Ψ，Δ，λ，φ）
【０００３】
入射角度を固定した場合、単一波長エリプソメータでは、（ｄ，ｎ，ｋ）の３つの未知数に対し、２つの独立変数しか測定できないので、ｄ，ｎ，ｋの内のいずれか一つを既知として固定する必要がある。
単一波長でも入射角度を変えると測定変数は増加する。しかしながら、入射角度（φ）の違いによる（Ψφ1 ，Δφ1 ）と（Ψφ2 ，Δφ2 ）に強い相関関係があるため、ｄ，ｎ，ｋを精度良く求めることは難しい。
【０００４】
【発明が解決しようとする課題】
偏光変化量は、光が通る体積、（位相角（β）×ビーム径の面積）で表すことができる。ビーム径が一定とすると偏光変化量は次のようになる。
偏光変化量∝膜厚（ｄ）×（複素屈折率Ｎ）×φ
ここにおいて、φは入射角である。膜厚（ｄ）が薄く、Ｎが低ければ、位相角（β）の変化が小さくなり、ｄとＮの間に相関関係が強くなる。
【０００５】
前述した問題は、極薄膜２層構造の場合、上記にさらに各層の間に相関関係が発生する可能性があり、偏光変化量を表している測定結果（ΨE ( λi ) ，ΔE ( λi ) ）から、各層のｎ，ｋ，ｄを求めることが困難となる。
【０００６】
前述したように、測定された多層薄膜の偏光変化を表している測定スペクトル（ΨE ( λi ) ，ΔE ( λi ) ）は、前記基板のｎ，ｋ情報、各層のｎ，ｋ，ｄの情報のすべてを含んでいるが、これから、前記基板のｎ，ｋ情報、各層のｎ，ｋ，ｄの情報の唯一の組み合わせを算出することはできない（基板のみの場合を除く）。この前記唯一の組み合わせを探す方法は、分光エリプソメータデータの解析と呼ばれる。解析を行うときは、前記基板のｎ，ｋ情報、各層のｎ，ｋ，ｄの情報を利用してモデルをたてる。この中の基板や各層のｎ，ｋ情報には、リファレンス（既知のテーブルデータ）や分散式、または似たような材質の単層薄膜の光学定数を使用する。
分散式とは、物質の誘電率の波長依存性を示す式であり、近赤外から紫外線領域では、この誘電率ε（λ）は材料の構成原子の結合様式から決定される。
分散式として、調和振動子をもとにした計算式、量子力学をもとにした計算式、経験式等が知られており、通常２つ以上のパラメータを含んでいる。
前述したモデルに含まれているすべての未知数（各層の厚さや、分散式パラメータ、混合比など）を変化させながら、測定データにあわせていく。
これをフィッティングといい、このフィッティングの結果、各層の膜厚や混合比などが求まり、分散式パラメータからは、材料の誘電率ε（λ）を計算することができる。材料の誘電率と屈折率は下記の関係にある。
ε＝Ｎ2
本発明では、様々な理由により１つの分散式では膜の誘電率の波長依存性を特定できないかまたは困難な場合についても問題にして、有効媒質論（Effective Medium Theory)を利用して、有効誘電率(Effective Dielectric Function) を計算する。一般的に、例えばホスト材料中に様々な誘電率を持つＮ個の物質（ゲスト）が混ざり合っている場合、有効誘電率（ε）は、下記の様に表される。
ｆa （εa −ε）／（εa ＋２ε）＋ｆb （εb −ε）／（εb ＋２ε）＋ｆc （εc −ε）／（εc ＋２ε）＝０
ここで、
ε : 求めようとする有効誘電率
εa , εb , εc : 球状物質Ａ，Ｂ，Ｃの誘電率
ｆa , ｆb , ｆc : 各物質の混合比（Volume Fraction 以下Vf）で、
ｆa ＋ｆb ＋ｆc ＝1
基板上の膜が不均一または不連続や、いくつかの材料が混ざり合っている場合、波長オーダーより十分小さく、物理的に混合している複数の物質から成る媒質については、有効媒質近似（ＥＭＡ）を利用してモデルをたてる。
物質Ａと物質Ｂと物質Ｃが混合している場合について説明する。このときの有効媒質近似（ＥＭＡ）は、物質Ａの混合比、物質Ｂの混合比、物質Ｃの混合比、ＡとＢとＣの混合層の膜厚、誘電率には分散式やリファレンスデータなどを推定してフィッティングを行い、評価する。
【０００７】
本発明の目的は、膜厚や複素屈折率などの組み合わせモデルを設定し、そのシミュレーションスペクトルを算出して、そのシミュレーションスペクトルと測定スペクトルとのフィッティングを広範囲最極小値計算法（Ｅｘｔｅｎｄｅｄ ＢＬＭＣ）を使用して行うことにより、極薄膜２層構造を決定する、分光エリプソメータを用いた極薄膜２層構造解析方法を提供することにある。
【０００８】
【課題を解決するための手段】
前記目的を達成するために、本発明による極薄膜２層構造を測定し解析する方法は基本的に、
まず、分光エリプソメータを用いて極薄膜２層構造を測定し、分光スペクトルを得る。
新しい解析方法は基本的には３つの段階から成立している。
解析第１段階の目的は、実際のサンプルに良く合うと思われるモデルを複数選択して初期値を決定する。
解析第２段階では、第１段階で得られた初期値をもとにＥｘｔｅｎｄｅｄ ＢＬＭＣを行う。
必要に応じて、解析第３段階では最後のフィッティング、確認、保存を行う。
【０００９】
本発明で利用する超薄膜および薄膜計測方法で用いる極小値計算法（ＢＬＭＣ）は、
極小値計算法（ＢＬＭＣ）による分光エリプソメータを用いた計測対象の基板表面の超薄膜および薄膜計測方法において、
計測対象の基板表面の薄膜を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE (λi )とΔE (λi )を得るΨE ，ΔE スペクトル測定ステップと、
前記基板の（Ｎ 0 （ｎ0 ,ｋ0 ）），基板上の薄膜の（ｄ,Ｎ（ｎ,ｋ））を分散式を用いて仮定し、さらに予想される範囲内にある複数の膜厚（ｄ±ｍΔｄ）および予想される範囲内にある複数の入射角（φ±ｍΔφ）を設定するステップと、
前記入射角と膜厚の組み合わせにもとづいて、分散式（ＤＳＰ）のパラメータ（εs ,ωt ）のフィッティングを行うステップと、
前記フィッティングにより得られた各ΨM (λi )とΔM (λi )の中から前記ΨE (λi )とΔE (λi )との差の最も少なくなる膜厚（ｄbest）と入射角（φbest）の組み合わせを設定したモデルのフィッティング結果（ＤＳＰbest）を選択するステップと、
前記選択するステップで選択された入射角（φbest）を確定値として、膜厚（ｄbest）と分散式（ＤＳＰbest）のフィッティングを行うステップで構成されている。
【００１０】
前記目的を達成するために本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
計測対象の基板上の極薄膜２層構造を、分光エリプソメータにより入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE(λi)とΔE (λi )を得るΨE ，ΔE スペクトル測定段階と、極薄膜２層構造の基板の複素屈折率（Ｎ0（ｎ0，ｋ0 ））および各薄膜の材料（Ｍａｔ1，Ｍａｔ2）の考えられる光学定数に係る複素屈折率（Ｎ1（ｎ1，ｋ1 ）），（Ｎ2 （ｎ2 ，ｋ2 ））、各層の膜厚（ｄ1，ｄ2）を利用し、複数のモデルをたてる第１ステップ、
前記各モデルごとに前記測定スペクトルとのフィッティングを行って平均二乗誤差（χ2 ）を得る第２ステップ、および
前記各モデルごとのフィッティングの結果、最低の平均二乗誤差（χ2 ）の 値をもつモデルまたは、あらかじめ設定した膜厚の最大、最小値の中に入っている平均二乗誤差（χ2 ）の値が最も低いモデルを選択し、膜厚（ｄ1(best) ，ｄ2(best) ）を含むそのモデルに係る各層の値を得る第３ステップ
から成る解析第１段階と、
解析第１段階で得られた値を、前記解析第１段階後の新しいモデルの初期値として設定する第１ステップ、
前記解析第１段階の第３ステップにて得られたいずれか一方の層の膜厚（ｄ1(best) またはｄ2(best) ）を中心値として、そのまわりの複数点ごとに、もう一方の層の膜厚（ｄ2(best)またはｄ1(best)）を変化させ最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いてフィッティングを行う第２ステップ、および
前記複数点ごとの最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いて行ったフィッティングの結果から最低の平均二乗誤差（χ2）の値をもつモデルまたは、あらかじめ設定した膜厚、分散式パラメータ、入射角のそれぞれ最大、最小値の中に入っている平均二乗誤差（χ2 ）の値が最も低いモデルを選択する第３ステップ
から成る解析第２段階と、
解析第２段階で選択されたモデルの複素屈折率及び膜厚を利用して最終的なフィッティングを行う第１ステップ、および
前記最終的なフィッティングで得られた結果の確認を行う第２ステップ
から成る解析第３段階と
から構成されている。
【００１１】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第２段階の第２ステップにおいて、前記いずれか一方の層は、２層構造中の材料において、光学定数がより分からない方の層とするものである。
【００１２】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第２段階の第２ステップにおいて、光学定数がより分かっている方の層の膜厚（ｄ1(best) またはｄ2(best)）を中心値として、その中心値の±１０％までの範囲の中で変化させた膜厚ごとに、もう一方の層について最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を行う。
また、本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第１段階の第１ステップは、１、２層どちらかの層が、複数の材料が混ざり合っている場合に、各材料の混合比（Vf1 ，Vf2）を利用し、有効媒質近似を用いてモデルをたてるものである。
【００１３】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第３段階の第１ステップにおいて、前記解析第２段階の第３ステップで選択したモデルの１、２層の光学定数を固定した１、２層目の膜厚のフィッティングまたは、モデルの１層の光学定数を固定した１、２層目の膜厚および２層目の光学定数のフィッティングを行うものである。
【００１４】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
上述した分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第３段階の第２ステップにおいて、フィッティングで得られた結果が、あらかじめ設定した膜厚と分散式パラメータでそれぞれの最大、最小値の間に入っているか否かを確認し、
入っていない場合は前記解析第１段階の第１ステップに戻り、極薄膜２層構造の基板の複素屈折率（Ｎ 0 （ｎ 0 ，ｋ 0 ））および各薄膜の材料（Ｍａｔ 1 ，Ｍａｔ 2 ）の考えられる光学定数に係る複素屈折率（Ｎ 1 （ｎ 1 ，ｋ 1 ）），（Ｎ 2 （ｎ 2 ，ｋ 2 ））、各層の膜厚（ｄ 1 ，ｄ 2 ）を利用し、複数のモデルをたてて前記解析第１段階の第２ステップ以降を行うものである。
【００１５】
本発明による分光エリプソメータを用いた極薄膜２層構造の解析方法は、
計測対象の基板上の極薄膜２層構造を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE(λi)とΔE (λi )を得るΨE ，ΔE スペクトル測定段階と、
１，２層どちらかの層が、複数の材料が混ざり合っている場合、有効媒質近似を用いてモデルをたて、極薄膜２層構造の基板の複素屈折率（Ｎ0（ｎ0 ，ｋ0 ））および各薄膜の材料（Ｍａｔ1 ，Ｍａｔ2 ）の考えられる複素屈折率（Ｎ1（ｎ1 ，ｋ1 ）），（Ｎ2（ｎ2 ，ｋ2 ））、材料が混ざり合っている場合の各材料の混合比（Vf1，Vf2）、各層の膜厚（ｄ1 ，ｄ2 ）を利用し、複数のモデルをたてる第１ステップ、
前記各モデルごとに前記測定スペクトルとのフィッティングを行う第２ステップ、および
前記各モデルごとのフィッティングの結果、最低の平均二乗誤差（χ2 ）の値をもつモデルまたはあらかじめ設定した膜厚と混合比のそれぞれの最大、最小値の中に入っている平均二乗誤差（χ2 ）の値の最も低いモデルに係る、膜厚（ｄ1(best) ，ｄ2(best) ），混合比Vf(best)を含む各層の値を得る第３ステップ
から成る解析第１段階と、
前記解析第１段階の第３ステップで得られた膜厚・混合比の値を、前記解析第１段階後の新しいモデルの初期値とし、未知の分散式が入っている他の層の前記解析第１段階の第３ステップで得られた膜厚の値を基に、予想される範囲内にある膜厚（（ｄ1±ｍΔｄ1 ）または（ｄ2 ±ｍΔｄ2 ））を設定し、もう一方の層の膜厚についても、前記解析第１段階で得られた膜厚の値を中心として、そのまわりの複数点（（ｄ2±ｍΔｄ2 ）または（ｄ1 ±ｍΔｄ1 ））を設定し、また、前記解析第１段階の第３ステップで得られた混合比の値を中心値として、そのまわりの複数点（Vf±ｍΔVf）を設定する第１ステップ、
混合比に係る複数点（Vf±ｍΔVf）と未知の分散式が入っている他の層のもう一方の層の膜厚に係る複数点（（ｄ2 ±ｍΔｄ2 ）または（ｄ1±ｍΔｄ1 ））との組み合わせのなかで、未知分散式が入っている他の層について最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いてフィッティングを行う第２ステップ、および
混合比と未知の分散式が入っている他の層のもう一方の層の膜厚との組み合わせごとに得られた結果の中から、最低の平均二乗誤差（χ2 ）の値またはあらかじめ設定した膜厚、分散式パラメータ、混合比、入射角のそれぞれ最大、最小値の間に入っている最低の平均二乗誤差（χ2 ）の値を持つ組み合わせのモデルを選択する第３ステップ
から成る解析第２段階と、
前記解析第２段階の第３ステップで選択されたモデルの複素屈折率及び膜厚を基に、両膜厚、混合比および、分散式のフィッティングまたは、両膜厚、混合比のフィッティングを行う第１ステップ、および
該第１ステップのフィッティングで得られた結果の確認を行う第２ステップ
から成る解析第３段階と
から構成されている。
【００１６】
【発明の実施の形態】
以下図面等を参照して本発明による分光エリプソメータを用いた極薄膜２層構造データの解析方法の実施の形態を説明する。
図１は、本発明方法で使用するエリプソメータの構成を示す略図である。
なお測定の対象の試料 (サンプル) ４の一部を拡大して示してある。
【００１７】
図１に示されている分光エリプソメータにより、後述する方法の分光測定データの獲得ステップが実行される。先ず、測定装置について簡単に説明する。Ｘｅランプ１は、多数の波長成分を含む、いわゆる白色光源である。このＸｅランプ１の発光は光ファイバ２を介して偏光子３に導かれる。偏光子３により偏光された光は、測定対象であるサンプル４の表面に特定の入射角（例えばφ＝７５．００°）で入射させられる。サンプル４は後述する基板表面に２層に薄膜を形成した測定試料である。
サンプル４からの反射は、光弾性変調器（ＰＥＭ）５を介して検光子６に導かれる。光弾性変調器（ＰＥＭ）５により５０ｋＨｚの周波数に位相変調されて、直線から楕円偏光までが作られる。そのため、数ｍ秒の分解能でΨ，Δを決定することができる。検光子６の出力は光ファイバ７を介して分光器８に接続される。分光器８の出力データがデータ取込部９に取り込まれ、分光測定データの獲得ステップを終了する。なお、ＰＥＭ５の位置は偏光子３の後か検光子６の前どちらでも可能とする。
【００１８】
図２は、本発明による薄膜計測方法の実施形態を示す流れ図である。
図３は、前記実施形態における解析第１段階を示す詳細説明図である。
図４は、前記実施形態における解析第２段階を示す詳細説明図である。
図５は、前記実施形態における解析第３段階を示す詳細説明図である。
次の段落に本発明で用いる記号を一括して示す。
【００１９】
（記号の定義）
Ｓｕｂ ： 基板（物理定数は既知，バルクとして取り扱い可能）
Ｍａｔ : 薄膜材料（物質の光学定数）
Ｍａｔij： j 番目のモデルのｉ 層目の材料（同上）
ｄ i : i 層目の膜厚
ｄ i (best): フィッティングによって得られたi 層目の膜厚
ｄij： j 番目のモデルのｉ 層目の膜厚
ｄij(best)： フィッティングで得たj 番目のモデルのｉ 層目の膜厚
X2 : 平均二乗誤差（χ2 ）の値
X2(j) : j 番目のモデルにおける平均二乗誤差（χ2 ）の値
Void : ｎ＝１，ｋ＝０の物質
Vf(ij) : j 番目のモデルのｉ 層目の混合比(Volume fraction)
Vf(ij)(best) : フィッティングで得られた j 番目のモデルのｉ 層目の混合比
【００２０】
（フィッティングの結果のχ2 の意味）
Ｎ個の測定データ対Ｅｘｐ（ｉ＝１，２．．．，Ｎ）と前記モデルの対応するＮ個のモデルの計算データ対Ｍｏｄ（ｉ＝１，２．．．，Ｎ）とし、測定誤差は正規分布をするとし、標準偏差をσi とすると、平均二乗誤差（χ2 ）は、次のようにして与えられる。
【００２１】
（測定段階データ測定ステップ）
測定は図１に示す装置で測定を行う。計測対象４の基板表面の２層構造の薄膜（図中の拡大図を参照されたい）を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE ( λi ) とΔE ( λi ) を得るΨE ，ΔE スペクトルを測定する。
（測定段階データ保存ステップ）
前のステップで測定したデータを保存して比較の対象データとする（図２参照）。
【００２２】
（解析第１段階ステップ１）
このステップで、基板の（Ｎ 0 （ｎ0 ,ｋ0 ））および各薄膜の材料（Ｍａｔ1 ，Ｍａｔ2 ）の考えられる複素屈折率（Ｎ 1 （ｎ1 ,ｋ1 ）），（Ｎ 2 （ｎ2 ,ｋ2 ）），各層の膜厚（ｄ1 ,ｄ2 ）を利用し、いくつかのモデルをたてる。
この実施例では、下記各モデル(1) 〜モデル(4) を設定したものとする。
図３の３１が示す部分に前記各モデルを略図示してある。
モデル(1) は、基板（Ｓｕｂ）の上に第１層(Ｍａｔ1 の光学定数，膜厚ｄ11)と第２層（Ｍａｔ2 の光学定数，膜厚ｄ21）を形成したものである。
モデル(2) は、基板（Ｓｕｂ）の上に第１層(Ｍａｔ1 の光学定数，膜厚ｄ12)と第２層(Ｍａｔ2 の光学定数＋Ｖｏｉｄ、膜厚ｄ22)を形成したものである。なおＶｏｉｄとは屈折率１の空間である。
モデル(3) は、基板（Ｓｕｂ）の上に第１層(Ｍａｔ1 の光学定数，膜厚ｄ13)と第２層(Ｍａｔ2 の光学定数＋Ｍａｔ3 の光学定数、膜厚ｄ23)を形成したものである。なお、（Ｍａｔ2 ＋Ｍａｔ3 ）は、材料２と材料３をある比率で混合したものを意味する。
モデル(4) は、基板（Ｓｕｂ）の上に第１層（Ｍａｔ1 の光学定数，膜厚ｄ14）と第２層（Ｍａｔ2 の光学定数の基準にする量(Special reference) ＋Ｖｏｉｄ、膜厚ｄ24)を形成したものである。なお、基準にする量とは似たような試料を用いて得られた光学定数である。
ここで、モデル(2) 〜(4) の２層目で設定されている混合比からは、有効媒質近似論を用いて、均質膜としての光学定数を求めることができる。なお前記各モデルは、基板Ｓｕｂの材料の光学定数は既知、第１層の光学定数は、ほぼ既知で、第２層の光学定数および、第１層の膜厚ｄ1 、第２層の膜厚ｄ2 は、未知（不確か）であることを前提にして前述の４個のモデルをたててある。
【００２３】
（解析第１段階ステップ２）
前記ステップ１で選定した４個のモデル(1) 〜(4) のそれぞれについて、前記測定スペクトルより得られた測定データΨE ，ΔE とのフィッティングを行う。
図３の３２の示す部分に各モデルのフィッティングの対象とフィッティングの結果得られたデータのχ2 値を示してある。
モデル(1) では、第１層目の膜厚ｄ11と第２層目の膜厚ｄ21をフィッティングして、その結果ｄ11(best)，ｄ21(best)とχ2 値 X2 (1) を得る。
モデル(2) では、第１層目の膜厚ｄ12，第２層目の膜厚ｄ22と第２層目の混合比のフィッティングをして、その結果ｄ12(best)，ｄ22(best)，Vf22(best)およびχ2 値 X2 (2) を得る。
モデル(3) では、第１層目の膜厚ｄ13，第２層目の膜厚ｄ23と第２層目の混合比のフィッティングをして、その結果ｄ13(best)，ｄ23(best)，Vf23(best)およびχ2 値 X2 (3) を得る。
モデル(4) では、第１層目の膜厚ｄ14，第２層目の膜厚ｄ24と第２層目の混合比のフィッティングをして、その結果ｄ14(best)，ｄ24(best)，Vf24(best)およびχ2 値 X2 (4) を得る。
（解析第１段階ステップ３）
前記複数組のフィッティングの結果から、最低χ2 値またはあらかじめ設定した膜厚の最大、最小値の間に入っている最低χ2 値のモデルの結果を選択するステップであり、これを図３の３３に示す。
【００２４】
（解析第２段階ステップ１）
図４の４１に解析第２段階で使用するモデル（初期値）が示されている。
この例では、Ｍａｔ1 の光学定数は、ほぼ既知となっているので、解析第１段階で使用した光学定数（Ｍａｔ1 ）をそのまま利用している。２層目の材料の光学定数（Ｍａｔ2 ）は未知のため、ここでは分散式を用いている。なお、膜厚（ｄ）の初期値は、解析第１段階ステップ３で得られた値を使用する。
（解析第２段階ステップ２）
このステップでは、図４の４２に示すように、前記モデル（初期値）の１層目膜厚（ｄ1(best) ）を中心値としてｄ1(best) ＋１０％，ｄ1(best) ＋ ５％，ｄ1(best) ，ｄ1(best) −５％，ｄ1(best) −１０％上下に変化させる。そして前記複数点（５点）に対して２層目の厚さ（ｄ2(best) ）を±１０％の範囲で変化させてＢＬＭＣを行う。
各モデルの２層目の厚さｄ2jと光学定数Ｍａｔ2jとχ2 値 X2 (j) を得る。
（解析第２段階ステップ３）
ステップ２で行われた結果から最低χ2 値またはあらかじめ設定した膜厚と分散式パラメータでそれぞれの最大、最小値の間に入っている最低χ2 値のモデルを選択する（図３、４３）。
【００２５】
（解析第３段階ステップ１）
解析第２段階ステップ３で選択したモデルの１，２層の光学定数を固定して１，２層目の膜厚のフィッティングを行う。または
１層の光学定数を固定して１，２層目の膜厚および２層目の光学定数をフィッティングする（図５、５１参照）。
（解析第３段階ステップ２）
前記ステップ１の結果があらかじめ設定した膜厚と分散式パラメータでそれぞれの最大、最小値の間に入っている最低χ2 値かどうかの確認をする（図５、５２参照）。なお、結果の確認が妥当でないときは解析第１段階に戻り、新しいモデルを設定して、フィッティング（ステップ２）を行う。
（解析第３段階ステップ３）
前記ステップ２で結果が妥当であるときは保存する（図５、５３参照）。
【００２６】
【実施例】
次に、前記同様に、１層目の光学定数は略既知、２層目の光学定数および１，２層目の膜厚が、未知の場合の実施例について説明する。図２に示した流れ図をそのまま使用することができる。この実施例は、基板ＳｕｂがSiで、１層目の材料をSiO2, ２層目の材料を SiNX としたものである。
前記試料を図１に示す装置で測定を行う。計測対象４の基板表面の２層構造の薄膜を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE ( λi ) とΔE ( λi ) を得るΨE ，ΔE スペクトルを測定して測定したデータを保存して比較の対象データとする。
【００２７】
図３の解析第１段階で、各薄膜の材料（SiO2 ， SiNx ）の考えられる複素屈折率（Ｎ 1 （ｎ1 ,ｋ1 ）），（Ｎ 2 （ｎ2 ,ｋ2 ））と膜厚（ｄ1 ,ｄ2 ）を利用し、モデルを準備する。この実施例では、下記各モデル(1) 〜モデル(4) を設定したものとする。
【００２８】
（解析第１段階ステップ１）
モデル(1) は、基板（Ｓｕｂ：Si）の上に第１層( SiO2の光学定数、膜厚ｄ11) と第２層（Si3N4 の光学定数，膜厚ｄ21）を形成したものである。
モデル(2) は、基板（Ｓｕｂ：Si）の上に第１層( SiO2の光学定数、膜厚ｄ12) と第２層( Si3N4 の光学定数＋Ｖｏｉｄ、膜厚ｄ22 )を形成したものである。なおＶｏｉｄとは屈折率１の空間である。
モデル(3) は、基板（Ｓｕｂ：Si）の上に第１層( SiO2の光学定数、膜厚ｄ13) と第２層( Si3N4 の光学定数＋ SiNX の光学定数、膜厚ｄ23) を形成したものである。
なお、（Si3N4 ＋SiNx）は、Si3N4 とSiNxをある比率で混合したものを意味する。
モデル(4) は、基板（Ｓｕｂ：Si）の上に第１層（SiO2の光学定数、膜厚ｄ14）と第２層（ SiNx の光学定数（既知の基準量）＋Ｖｏｉｄ、膜厚ｄ24 )を形成したものである。ここで、モデル(2) 〜(4) の２層目で設定されている混合比からは、有効媒質近似論を用いて、均質膜としての光学定数を求めることができる。なお前記各モデルは、基板Ｓｕｂの材料の光学定数は既知、第１層のSiO2の光学定数は、ほぼ既知で、第２層のSiNxの光学定数および、第１層の膜厚ｄ1,第２層の膜厚ｄ2 は、未知（不確か）であることを前提にしている。
【００２９】
（解析第１段階ステップ２）
前記ステップ１で選定した４個のモデル(1) 〜(4) のそれぞれについて、前記測定スペクトルより得られた測定データΨE ，ΔE とのフィッティングを行う。
図３の３１，３２の示す部分に各モデルのフィッティングの対象とフィッティングの結果得られたデータのχ2 値を示してある。
モデル(1) では、第１層目の膜厚ｄ11と第２層目の膜厚ｄ21をフィッティングしてその結果ｄ11(best)、ｄ21(best)とχ2 値 X2 (1) を得る。
モデル(2) では、第１層目の膜厚ｄ12、第２層目の膜厚ｄ22と第２層目の混合比のフィッティングしてその結果ｄ12(best)、ｄ22(best)、Vf22(best)およびχ2 値 X2 (2) を得る。
モデル(3) では、第１層目の膜厚ｄ13、第２層目の膜厚ｄ23と第２層目の混合比のフィッティングしてその結果ｄ13(best)、ｄ23(best)、Vf23(best)およびχ2 値 X2 (3) を得る。
モデル(4) では、第１層目の膜厚ｄ14、第２層目の膜厚ｄ24と第２層目の混合比のフィッティングしてその結果ｄ14(best)、ｄ24(best)、Vf24(best)およびχ2 値 X2 (4) を得る。
【００３０】
（解析第１段階ステップ３）
前記複数組のフィッティングの結果から、最低χ2 値またはあらかじめ設定した膜厚の最大、最小値の間に入っている最低χ2 値のモデルの結果を選択するステップであり、これを図３の３３に示す。
【００３１】
（解析第２段階ステップ１）〜（解析第３段階ステップ３）は前述したとおりである。
【００３２】
【発明の効果】
以上詳しく説明したように、本発明によれば、基板上に形成された極薄膜２層の厚さと光学定数を、分光エリプソメータを用いて得たデータを、広範囲最極小値計算法（Ｅｘｔｅｎｄｅｄ ＢＬＭＣ）を用いて処理することにより解析することができる。
前述したように基本的に３段階が含まれている手順は、２層以上の様々な構造の解析にも利用できる。
【００３３】
本発明方法によれば、
１．極薄膜多層構造でも、この方法を用いることで信頼できる膜厚（各層において）および光学定数（少なくとも１層）が得られる。
２．この手順の第１段階を行うことで、未知のパラメータの初期値範囲を狭めることができる。
３．Ｅｘｔｅｎｄｅｄ ＢＬＭＣを使用することで、好ましくないローカルミニマム（Ｌｏｃａｌ Ｍｉｎｉｍｕｍ）に落ちこむことは劇的に減少し、結果の信頼性があがる。
【００３４】
以上詳しく説明した実施例について、本発明の範囲内で種々の変形を施すことができる。理解を容易にするために、データの取得、モデルの設定に関連して、一貫してΨ，Δを用いて説明した。当業者には良く知られている以下のデータ対を用いても同様な、測定およびフィッティングが可能であり、本発明の技術的範囲に含まれるものである。
（ｎ，ｋ） （εr ，εi ） ( tan Ψ，cos Δ) 、 (Ｉs,Ｉc ）
【００３５】
基板上のＭａｔ1 ，Ｍａｔ2 には、極薄膜誘電体材料多層構造だけでなく、様々な厚さや材料など、幅広いアプリケーションにも使用することができる。また、実施例として光弾性変調器（ＰＥＭ）の例を示したが、ＰＥＭ以外のエリプソメータを使用することもできる。
【００３６】
上記手段の全部または１部を行う場合もあり、これも本発明の技術的範囲に含まれるものである。
【００３７】
基板もSiの他、ガラスや石英、化合物半導体なども同様に利用できる。また、基板の種類によらず、どんな平坦な基板でも、あれている基板でも使用することができる。
【００３８】
分散式には、クラシカル（古典力学）、アモルファス（量子力学）、経験式の他、様々な式・パラメータも使用可能であり、これも本発明の技術的範囲に含まれるものとする。
【００３９】
例ではあらゆる全てのパラメータを同時にフィッティングすると説明したが、別々にフィッティングする場合もあり、これも本発明の技術的範囲に含まれるものとする。
【００４０】
Ｅｘｔｅｎｄｅｄ ＢＬＭＣの一部であるＢＬＭＣでは、前述したとおり、入射角度をフィッティングする場合がある。また、手順では入射角度と様々なパラメータを同時にフィッティングするとしたが、別々にフィッティングする場合や、入射角度を固定する場合もあり、これらも全て本発明の技術的範囲に含まれるものとする。
【００４１】
入射角度はＢＬＭＣ以外でも、一般的なパラメータとしてフィッティングすることがあり、これも本発明の技術的範囲に含まれるものとする。
【００４２】
例では、解析第１段階第１ステップのモデル準備過程において、基板はSi、第１層目はSiO2、第２層目は(1) Si3N4 、(2) Si3N4 とＶｏｉｄ、(3) Si3N4 と SiNX 、および(4) SiNx とＶｏｉｄの(1) 〜(4) のモデルとしたが、製造プロセスになどによってモデルの種類、数が変わる場合もあり、これも本発明の技術的範囲に含まれるものとする。
【００４３】
上記のように、２ステップが含まれている手順は、２層以上の様々な構造にも対応可能である。
【００４４】
例ではＥＭＡを使用すると説明したが、他の有効媒質近似論の使用も可能であり、これも本発明の技術的範囲に含まれるものとする。
【図面の簡単な説明】
【図１】 本発明方法で使用するエリプソメータの構成を示す略図であり、測定の対象の試料 (サンプル) ４の一部を拡大して示してある。
【図２】 本発明による薄膜計測方法の実施形態を示す流れ図である。
【図３】 前記実施形態における解析第１段階を示す詳細説明図である。
【図４】 前記実施形態における解析第２段階を示す詳細説明図である。
【図５】 前記実施形態における解析第３段階を示す詳細説明図である。
【符号の説明】
１ Ｘｅランプ
２ 光ファイバ
３ 偏光子
４ サンプル
５ 光弾性変調器（ＰＥＭ）
６ 検光子
７ 光ファイバ
８ 分光器
９ データ取込部
３１ モデル複数組を準備する解析第１段階第１ステップ
３２ フィッティングを行う解析第１段階第２ステップ
３３ フィッティングの結果から、最低の平均二乗誤差（χ 2 ）の値のモデルの結果を選択する解析第１段階第３ステップ
４１ ３３で得られた結果を新しいモデルの初期値として設定する解析第２段階第１ステップ
４２ 一方の層の膜厚複数点ごとに、もう一方の層についてＢＬＭＣを行う解析第２段階第２ステップ
４３ 前記フィッティング結果から最低の平均二乗誤差（χ 2 ）の値のモデルを選択する解析第２段階第３ステップ
５１ 前記選択されたモデルの最終フィッティングを行う解析第３段階第１ステップ
５２ 結果の確認を行う解析第３段階第２ステップ
５３ 保存を行う解析第３段階第３ステップ
Claims (7)
- 計測対象の基板上の極薄膜２層構造を、分光エリプソメータにより入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE(λi)とΔE (λi )を得るΨE ，ΔE スペクトル測定段階と、極薄膜２層構造の基板の複素屈折率（Ｎ0（ｎ0，ｋ0 ））および各薄膜の材料（Ｍａｔ1，Ｍａｔ2）の考えられる光学定数に係る複素屈折率（Ｎ1（ｎ1，ｋ1 ）），（Ｎ2 （ｎ2 ，ｋ2 ））、各層の膜厚（ｄ1，ｄ2）を利用し、複数のモデルをたてる第１ステップ、
前記各モデルごとに前記測定スペクトルとのフィッティングを行って平均二乗誤差（χ2 ）を得る第２ステップ、および
前記各モデルごとのフィッティングの結果、最低の平均二乗誤差（χ2 ）の 値をもつモデルまたは、あらかじめ設定した膜厚の最大、最小値の中に入っている平均二乗誤差（χ2 ）の値が最も低いモデルを選択し、膜厚（ｄ1(best) ，ｄ2(best) ）を含むそのモデルに係る各層の値を得る第３ステップ
から成る解析第１段階と、
解析第１段階で得られた値を、前記解析第１段階後の新しいモデルの初期値として設定する第１ステップ、
前記解析第１段階の第３ステップにて得られたいずれか一方の層の膜厚（ｄ1(best) またはｄ2(best) ）を中心値として、そのまわりの複数点ごとに、もう一方の層の膜厚（ｄ2(best)またはｄ1(best)）を変化させ最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いてフィッティングを行う第２ステップ、および
前記複数点ごとの最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いて行ったフィッティングの結果から最低の平均二乗誤差（χ2）の値をもつモデルまたは、あらかじめ設定した膜厚、分散式パラメータ、入射角のそれぞれ最大、最小値の中に入っている平均二乗誤差（χ2 ）の値が最も低いモデルを選択する第３ステップ
から成る解析第２段階と、
解析第２段階で選択されたモデルの複素屈折率及び膜厚を利用して最終的なフィッティングを行う第１ステップ、および
前記最終的なフィッティングで得られた結果の確認を行う第２ステップ
から成る解析第３段階と
から構成された分光エリプソメータを用いた極薄膜２層構造の解析方法。 - 請求項１記載の分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第２段階の第２ステップにおいて、前記いずれか一方の層は、２層構造中の材料において、光学定数がより分からない方の層とする分光エリプソメータを用いた極薄膜２層構造の解析方法。 - 請求項１または２記載の分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第２段階の第２ステップにおいて、光学定数がより分かっている方の層の膜厚（ｄ1(best) またはｄ2(best) ）を中心値として、その中心値の±１０％までの範囲の中で変化させた膜厚ごとに、もう一方の層について最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を行う分光エリプソメータを用いた極薄膜２層構造の解析方法。 - 請求項１乃至３のいずれか１つに記載の分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第１段階の第１ステップは、１、２層どちらかの層が、複数の材料が混ざり合っている場合に、各材料の混合比（Vf1 ，Vf2）を利用し、有効媒質近似を用いてモデルをたてる分光エリプソメータを用いた極薄膜２層構造の解析方法。 - 請求項１乃至４のいずれか１つに記載の分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第３段階の第１ステップにおいて、前記解析第２段階の第３ステップで選択したモデルの１、２層の光学定数を固定した１、２層目の膜厚のフィッティングまたは、モデルの１層の光学定数を固定した１、２層目の膜厚および２層目の光学定数のフィッティングを行う分光エリプソメータを用いた極薄膜２層構造の解析方法。 - 請求項１乃至５のいずれか１つに記載の分光エリプソメータを用いた極薄膜２層構造の解析方法において、
前記解析第３段階の第２ステップにおいて、フィッティングで得られた結果が、あらかじめ設定した膜厚と分散式パラメータでそれぞれの最大、最小値の間に入っているか否かを確認し、
入っていない場合は前記解析第１段階の第１ステップに戻り、極薄膜２層構造の基板の複素屈折率（Ｎ0 （ｎ0 ，ｋ0））および各薄膜の材料（Ｍａｔ1，Ｍａｔ2）の考えられる光学定数に係る複素屈折率（Ｎ1（ｎ1 ，ｋ1）），（Ｎ2 （ｎ2 ，ｋ2 ））、各層の膜厚（ｄ1，ｄ2 ）を利用し、複数のモデルをたてて前記解析第１段階の第２ステップ以降を行う分光エリプソメータを用いた極薄膜２層構造の解析方法。 - 計測対象の基板上の極薄膜２層構造を、入射光の波長を変えて各波長λi ごとの入射光と反射光の偏光の変化である測定スペクトルΨE(λi)とΔE (λi )を得るΨE ，ΔE スペクトル測定段階と、
１，２層どちらかの層が、複数の材料が混ざり合っている場合、有効媒質近似を用いてモデルをたて、極薄膜２層構造の基板の複素屈折率（Ｎ0（ｎ0 ，ｋ0 ））および各薄膜の材料（Ｍａｔ1 ，Ｍａｔ2 ）の考えられる複素屈折率（Ｎ1（ｎ1 ，ｋ1 ）），（Ｎ2（ｎ2 ，ｋ2 ））、材料が混ざり合っている場合の各材料の混合比（Vf1，Vf2）、各層の膜厚（ｄ1 ，ｄ2 ）を利用し、複数のモデルをたてる第１ステップ、
前記各モデルごとに前記測定スペクトルとのフィッティングを行う第２ステップ、および
前記各モデルごとのフィッティングの結果、最低の平均二乗誤差（χ2 ）の値をもつモデルまたはあらかじめ設定した膜厚と混合比のそれぞれの最大、最小値の中に入っている平均二乗誤差（χ2 ）の値の最も低いモデルに係る、膜厚（ｄ1(best) ，ｄ2(best) ），混合比Vf(best)を含む各層の値を得る第３ステップ
から成る解析第１段階と、
前記解析第１段階の第３ステップで得られた膜厚・混合比の値を、前記解析第１段階後の新しいモデルの初期値とし、未知の分散式が入っている他の層の前記解析第１段階の第３ステップで得られた膜厚の値を基に、予想される範囲内にある膜厚（（ｄ1±ｍΔｄ1 ）または（ｄ2 ±ｍΔｄ2 ））を設定し、もう一方の層の膜厚についても、前記解析第１段階で得られた膜厚の値を中心として、そのまわりの複数点（（ｄ2±ｍΔｄ2 ）または（ｄ1 ±ｍΔｄ1 ））を設定し、また、前記解析第１段階の第３ステップで得られた混合比の値を中心値として、そのまわりの複数点（Vf±ｍΔVf）を設定する第１ステップ、
混合比に係る複数点（Vf±ｍΔVf）と未知の分散式が入っている他の層のもう一方の層の膜厚に係る複数点（（ｄ2 ±ｍΔｄ2 ）または（ｄ1±ｍΔｄ1 ））との組み合わせのなかで、未知分散式が入っている他の層について最極小値計算法（ＢＬＭＣ:Best Local Minimum Calculation）を用いてフィッティングを行う第２ステップ、および
混合比と未知の分散式が入っている他の層のもう一方の層の膜厚との組み合わせごとに得られた結果の中から、最低の平均二乗誤差（χ2 ）の値またはあらかじめ設定した膜厚、分散式パラメータ、混合比、入射角のそれぞれ最大、最小値の間に入っている最低の平均二乗誤差（χ2 ）の値を持つ組み合わせのモデルを選択する第３ステップ
から成る解析第２段階と、
前記解析第２段階の第３ステップで選択されたモデルの複素屈折率及び膜厚を基に、両膜厚、混合比および、分散式のフィッティングまたは、両膜厚、混合比のフィッティングを行う第１ステップ、および
該第１ステップのフィッティングで得られた結果の確認を行う第２ステップ
から成る解析第３段階と
から構成された分光エリプソメータを用いた極薄膜２層構造の解析方法。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002110710A JP3937149B2 (ja) | 2002-04-12 | 2002-04-12 | 分光エリプソメータを用いた極薄膜２層構造の解析方法 |
| US10/488,596 US7196793B2 (en) | 2001-09-06 | 2002-09-06 | Method for analyzing thin-film layer structure using spectroscopic ellipsometer |
| KR1020047003286A KR100892743B1 (ko) | 2001-09-06 | 2002-09-06 | 분광타원계를 사용한 박막층 구조의 해석방법 |
| EP02798038A EP1435517B1 (en) | 2001-09-06 | 2002-09-06 | Method for analyzing thin-film layer structure using spectroscopic ellipsometer |
| PCT/JP2002/009080 WO2003023373A1 (fr) | 2001-09-06 | 2002-09-06 | Procede pour analyser une structure stratifiee en film mince au moyen d'un ellipsometre spectroscopique |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002110710A JP3937149B2 (ja) | 2002-04-12 | 2002-04-12 | 分光エリプソメータを用いた極薄膜２層構造の解析方法 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2003302334A JP2003302334A (ja) | 2003-10-24 |
| JP2003302334A5 JP2003302334A5 (ja) | 2005-12-15 |
| JP3937149B2 true JP3937149B2 (ja) | 2007-06-27 |
Family
ID=29393763
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002110710A Expired - Lifetime JP3937149B2 (ja) | 2001-09-06 | 2002-04-12 | 分光エリプソメータを用いた極薄膜２層構造の解析方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3937149B2 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5397693B2 (ja) | 2010-02-25 | 2014-01-22 | 大日本スクリーン製造株式会社 | 水素含有率取得装置および水素含有率取得方法 |
| US8666703B2 (en) | 2010-07-22 | 2014-03-04 | Tokyo Electron Limited | Method for automated determination of an optimally parameterized scatterometry model |
| US10663286B2 (en) * | 2017-08-22 | 2020-05-26 | Kla-Tencor Corporation | Measuring thin films on grating and bandgap on grating |
| CN114324184B (zh) * | 2021-12-30 | 2024-05-17 | 粤芯半导体技术股份有限公司 | 椭偏仪光谱浮动模型及建立方法 |
-
2002
- 2002-04-12 JP JP2002110710A patent/JP3937149B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2003302334A (ja) | 2003-10-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100892743B1 (ko) | 분광타원계를 사용한 박막층 구조의 해석방법 | |
| JP4633254B2 (ja) | 回折構造体、広帯域、偏光、エリプソメトリおよび下地構造の測定 | |
| US7239391B2 (en) | Method of analysis of multiple layer samples | |
| US5889592A (en) | Nondestructive optical techniques for simultaneously measuring optical constants and thicknesses of single and multilayer films | |
| KR100822419B1 (ko) | 시료 해석 방법 | |
| Elman et al. | Characterization of biaxially-stretched plastic films by generalized ellipsometry | |
| US5999267A (en) | Nondestructive optical techniques for simultaneously measuring optical constants and thicknesses of single and multilayer films | |
| CN112595673A (zh) | 一种单晶金刚石衬底光学常数测量方法 | |
| JP4136740B2 (ja) | 分光エリプソメータを用いた薄膜３層構造の解析方法 | |
| WO2002095372A1 (en) | Thin-film characteristic measuring method using spectroellipsometer | |
| JP3937149B2 (ja) | 分光エリプソメータを用いた極薄膜２層構造の解析方法 | |
| JP2003243467A (ja) | 偏光解析方法 | |
| US7463355B1 (en) | Nondestructive optical technique for simultaneously measuring optical constants and thickness of thin films | |
| JP3532165B2 (ja) | 分光エリプソメータを用いた薄膜計測方法 | |
| JP3928714B2 (ja) | 分光エリプソメータを用いた薄膜複数層構造の解析方法 | |
| JP3781245B2 (ja) | 半導体装置の製造方法 | |
| Stenzel et al. | A hybrid method for determination of optical thin film constants | |
| US7259850B2 (en) | Approach to improve ellipsometer modeling accuracy for solving material optical constants N & K | |
| JP3556183B2 (ja) | 基板上の化合物半導体層の組成決定方法 | |
| JP3613707B2 (ja) | 超薄膜および薄膜計測方法 | |
| JP3983093B2 (ja) | 分光エリプソメータを用いた多結晶化合物半導体の組成決定方法 | |
| Busch et al. | Application of ellipsometry for the accurate oxide layer measurement on silicon spheres | |
| JP4587690B2 (ja) | 超薄膜および薄膜計測方法 | |
| Edwards | Status and prospects for VUV ellipsometry (applied to high k and low k materials) | |
| Gibbons et al. | Dielectric functions of common YBCO substrate materials determined by spectroscopic ellipsometry |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20040601 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20040616 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050408 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051025 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060801 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060929 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20060929 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061212 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070131 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20070131 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070313 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070313 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 3937149 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130406 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130406 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130406 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140406 Year of fee payment: 7 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |