JP3959469B2 - 膜の屈折率及び厚さの測定方法、測定プログラム及びコンピュータ読取可能な記録媒体 - Google Patents

Description

本発明は、エリプソメトリーを用いて膜の屈折率及び厚さを測定する方法に関する。

光学フィルタのコーティングにおいて、設計が複雑になれば、積層工程における屈折率及び膜の成長速度の正確な制御が重要になる。エリプソメトリーは、形成される薄膜の屈折率、膜厚の測定精度が高いので、コーティング制御における有効な測定手段と考えられている。しかし、従来の光学モデルでは、一定の屈折率である均質な薄膜に関する測定方法しか開示されていない（下記非特許文献１参照）。

また多くの場合、リアルタイムに製膜の制御を行うことが望ましく、そのためには、ほぼ瞬時に積層物の特性を得ることが必要であり、最外層の屈折率及び厚さを、同時に決定するために、エリプソメトリーデータの「直接数値反転」（下記非特許文献２〜４参照）及びフィッティングプロセス（下記非特許文献５、６参照）が使用されるようになってきた。

この場合、「共通擬似基板近似」が有効であることが知られている。その理由は、多層又は基板の下層について完全に計算を行うことなく、最上層の特性評価が可能であるからである（下記非特許文献７参照）。しかし、この方法を拡張して、厚膜に適用することはできない。

一方、多層構造の下層に関する特性値を用いた繰り返し計算は、アルゴリズムが正確であるだけでなく、効率的である。その理由は、下層に関する特性値に基づき、次の外層に関する特性値が計算されるからである（下記非特許文献２〜４、６参照）。特に、２×２マトリックス法は、有効であり、反射及び透過の両方を扱うことができ、不均一な層の場合にも容易に拡張することができる。
M. Kildemo and B. Drevillon, "Real-time monitoring of the growth of transparent thin films by spectroscopic ellipsometry," Rev. Sci. Instrum., 67, 1956-1960(1996) I-F. Wu, J.B. Dottellis, and M. Dagenais, "Real-time in situ ellipsometric control of antireflection coatings for semiconductor laser amplifiers using SiOx," J. Vac. Sci. Technol. A, 11, 2398-2406(1993) G. M. W. Kroessen, G. S. Oehrlein, E. de Fresart, and M. Hverlag, "Depth profiling of the Ge concentration in SiGe alloys using in situ ellipsometry during reactive-ion etching," J. Appl. Phys., 73, 8017-8026(1993) D. Kouznetsov, A. Hofrichter, and B. Drevillon, "Direct numerical inversion method for kinetic ellipsometric data. I. Presentation of the method and numerical evaluation," Applied Optics, 41, 4510-4518(2002) M. Kildemo, "Real-time monitoring and growth control of Si-gradient-index structures by multiwavelength ellipsometry," Appl. Opt., 37, 113-124(1998) T. Heitz, A. Hofrichter, P. Bulkin, and B. Drevillon, "Real time control pf plasma deposited optical filters by multiwavelength ellipsometry," J. Vac. Sci. Technol. A, 18, 1303-1307(2000) D. E. Aspens, "Minimal-data approaches for determining outer layer dielectric responses of films from kinetic reflectometric and ellipsometric measurements," J. Opt. Soc. Am. A, 10, 974-983(1993)

しかし、従来報告されている２×２マトリックス法を使用したアルゴリズムでは、マトリックスの要素を２次項までのテイラー級数で近似するため、数ｎｍ程度以下の小さな膜厚変化にしか対応できず、急激な膜厚の変化が起こった場合計算不能に陥る欠点があった（上記非特許文献５参照）。また、屈折率と膜厚を求めるために高次（８次）の多項式を解かなければならなかった（上記非特許文献５参照）。また、フィッティング処理が必要である（上記非特許文献６参照）ため、計算に時間がかかりリアルタイム処理が困難であった。

一方、エリプソメトリーは、膜厚の増加とともに周期的に精度が大きく低下するという原理的な問題点を有しているため、膜厚の大きな積層膜の製膜過程を単一波長で測定することは困難であり、一般に多波長で測定が行われる。しかし、多層構造の下層に関する特性値を用いた繰り返し計算の多くは、単一の波長で測定したデータを用いる場合しか検討されていない（上記非特許文献２および３参照）。複数の波長に対応した手法についても、屈折率の波長依存性を適当な数式であらかじめ仮定するため、異なる波長における屈折率を精度よく求めることができなかった（上記非特許文献５および６参照）。

本発明は、以上の課題を解決すべく、エリプソメトリーを用い、近似や仮定を用いることなく、また高次多項式の数値解法やフィッティング処理を使用することなく、大きな膜厚変化にも対応して、単一または複数の波長において測定したエリプソメトリーの係数から、多層膜の最上位層の屈折率及び厚さを、高速に精度良く求めることができ、リアルタイム制御に適用可能な、膜の屈折率及び厚さの測定方法、測定プログラム及びコンピュータ読取可能な記録媒体を提供することを目的とする。

本発明の目的は、以下によって達成される。

即ち、本発明に係る膜の屈折率及び厚さの測定方法（１）は、エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する方法であって、前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、反射光を測定することによって得られる反射型エリプソメトリーの計数Ψ_r、Δ_rを使用し、ln｜Ｘ｜＝０（ここで、Ｘは、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r−Ｂ₂ ^pＡ₂ ^s、ｂ＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^sＡ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、
ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−Ｂ₁ ^pＡ₁ ^sを係数とする２次方程式ａＸ²＋ｂＸ＋ｃ＝０の根であり、ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^p
は、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行
列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ
’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒体との境界面の特性行列
Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）、又は、(｜ｃ｜²−｜ａ｜²)²−｜ｃ^*ｂ−ｃｂ^*｜²＝０（ここで、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r−Ｂ₂ ^pＡ₂ ^s、ｂ
＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^sＡ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−Ｂ₁ ^pＡ₁ ^s、ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1
によって計算される散乱行列Ｓ_mの行列要素として求められ、前記薄膜が新しく形成され
る側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m
、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)m
Ｌ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）を評価関数として、非線形方程式の数値解法を用いて、前記波長に対する屈折率を計算する屈折率計算ステップを含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（２）は、エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する方法であって、前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、前記多層膜を透過した透過光を測定することによって得られる透過型エリプソメトリーの計数Ψ_t、Δ_tを使用し、ln｜Ｘ｜＝０（ここで、Ｘは、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^sＡ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^sを係数とする、１次方程式ａＸ＋ｂ＝０の根であり、ρ_t＝tanΨ_texp(−ｊΔ_t)であり、係数Ａ₁ ^s、
Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｋ^s、Ｋ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多層膜と前記
薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求め
られ、前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前
記媒体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素
として求められる）、又は、(ａ^*＋ｂ^*)―(ａ＋ｂ)＝０（ここで、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^s
Ａ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^s、ρ_t＝tanΨ_texp(−ｊΔ_t)であり、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｋ^s、Ｋ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多層膜と前記薄膜との
間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、
Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒体と
の境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から
、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）を評価関数として、非線形方程式の数値解法を用いて、前記波長に対する屈折率を計算する屈折率計算ステップを含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（３）は、上記の膜の屈折率及び厚さの測定方法（１）又は（２）において、前記屈折率計算ステップによって計算された前
記屈折率を用いて、ｄ＝Re｛ｊ(1/4π)λ(ｎ²−ｎ₀ ²sinφ₀)^-1/2lnＸ｝（ここで、ｎは前記薄膜の屈折率、ｎ₀は前記媒体の屈折率、λは前記波長、φ₀は前記光の入射角である）によって、前記薄膜の膜厚ｄを計算する膜厚計算ステップをさらに含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（４）は、上記の膜の屈折率及び厚さの測定方法（３）において、複数の波長の各々に対して、前記屈折率計算ステップによって前記屈折率を計算し、且つ前記膜厚計算ステップによって前記膜厚を計算し、計算された複数の前記屈折率及び前記膜厚の良否を判断し、精度良く前記薄膜の膜厚を決定する良否判断ステップをさらに含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（５）は、上記の膜の屈折率及び厚さの測定方法（４）において、前記良否判断ステップの後に、前記良否判断ステップにおいて不良と判断された屈折率を新たに計算する修正ステップをさらに含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（６）は、上記の膜の屈折率及び厚さの測定方法（５）において、前記良否判断ステップが、同じ波長に対する、前記屈折率計算ステップによって計算された前記屈折率、及び前記膜厚計算ステップによって計算された前記膜厚からなる複数の組の中から、良好な組を選別する選別ステップと、前記良好な組を、所定のグループ数のグループに分類する分類ステップと、分類された前記グループを用いて、信頼性の高い膜厚を含む最適グループを決定する最適グループ決定ステップと、前記最適グループに含まれる膜厚の平均値及び規格化標準偏差を用いて、前記良好な組の中から、良好な屈折率及び良好な膜厚を併せ持つ良好な組を再度選別する再選別ステップと、前記良好な膜厚の平均値を計算することによって精度良く前記薄膜の膜厚を決定する平均値計算ステップとを含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（７）は、上記の膜の屈折率及び厚さの測定方法（６）において、前記修正ステップが、前記再選別ステップにおいて再度選別された前記良好な屈性率の平均値を計算し、該平均値から屈折率の探索範囲を決定する範囲決定ステップと、非線形関数の極小探査法によって、前記探索範囲の中で、前記良否判断ステップにおいて不良と判断された屈折率を新たに計算する再計算ステップとを含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（８）は、上記の膜の屈折率及び厚さの測定方法（７）において、前記グループ数が、前記良好な組の数に応じて多くなるようにあらかじめ定めた規則に基づいて決定され、前記分類ステップが、各々のグループに属する膜厚の分散値の総和が最小になるように、前記良好な組を、前記グループ数のグループに分類するステップを含み、前記最適グループ決定ステップが、分類された前記グループの要素の数が２以上かつ最大である１つのグループ、および、そのグループを除く前記グループの要素の数が２以上かつ最大のグループを仮の最適グループとする仮最適グループ群決定ステップと、すべての仮最適グループＧ₁、Ｇ₂、Ｇ₃、…について、すべての
良好な組を対象として、仮最適グループＧ_iの要素ｄ_ijの対数値lnｄ_ijの平均値(lnｄ)_ai
（添字ａは平均値を表す）及び標準偏差σを用いて計算した規格化標準偏差（σ／(lnｄ)_ai）をσ_dNiとして、それぞれの仮最適グループＧ_iの(lnｄ)_aiに対する許容偏差δ_diがδ_di＝γσ_dNiによって表され、σ_dNiが大きくなるほど大きくなるようにあらかじめ定めたγの値に基づき、すべての良好な膜厚についての（lnｄ_k−(lnｄ)_ai）²の総和が最小になるように、前記良好な組の膜厚がabs（lnｄ_k−(lnｄ)_ai）＜δ_d（ここで、absは()内の絶対値を求める演算を表す）の条件を満たす良好な組（ｎ_k、ｄ_k）をそれぞれの仮最適グループに取り込む処理をＧ₁、Ｇ₂、Ｇ₃、…の順に行う取込ステップと、前記仮最適グルー
プ群へのデータの取込処理によって変化した要素の組み合わせに応じて、仮最適グループ群Ｇ_iのうち要素数最大のグループおよび２番目のグループを決定し、それぞれ最適グル
ープＧ₁および準最適グループＧ₂とする並べ替えステップと、前記最適および準最適グループＧ_iの要素数をＮ（Ｇ_i）、規格化標準偏差をσ_dN（Ｇ_i）（ｉ＝１又は２）とし、「
Ｎ（Ｇ₁）＝Ｎ（Ｇ₂）」および「σ_dN（Ｇ₂）＜σ_dN（Ｇ₁）」の条件を両方満たす場合、最適グループＧ₁と準最適グループＧ₂とを入れ替える処理を行い、最終的に最適グループＧ₁を決定する入替ステップとを含み、前記再選別ステップが、前記最適グループの屈折
率ｎ_iの平均値ｎ_a及び標準偏差σを用いて規格化標準偏差σ_nN（＝σ／ｎ_a）を計算し、
δ_n＝４σ_nNによって許容偏差δ_nを決定し、前記良好な組の屈折率の中で、abs（ｎ_k−ｎ_a）＜δ_nを満たす屈折率を前記良好な屈折率とし、前記最適グループの膜厚データｄ_iの
対数値の平均値(lnｄ)_a及び標準偏差σを用いて規格化標準偏差σ_dN（＝σ／(lnｄ)_a）を計算し、許容偏差δ_dがδ_d＝γσ_dNによって表され、σ_dNが大きくなるほど大きくなるようにあらかじめ定めたγの値に基づき、前記良好な組の膜厚がabs（lnｄ_k−(lnｄ)_ai）＜δ_dの条件を満たす膜厚を前記良好な膜厚とし、屈折率と膜厚が両方とも良好であること
を条件として、良好なデータの組（ｎ_k、ｄ_k）を決定するステップであることを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定方法（９）は、上記の膜の屈折率及び厚さの測定方法（８）において、前記範囲決定ステップが、前記良好な屈性率の平均値ｎ_gaと、あらかじめ定めた探索範囲の相対値ｒを用い、前記探索範囲の上限値を（１．０＋ｒ）ｎ_ga、下限値を（１．０−ｒ）ｎ_gaと決定するステップを含み、前記再計算ステップが、ｆ（ｎ）＝(Ｉ_S'−Ｉ_S)²−(Ｉ_C'−Ｉ_C)²（ここで、エリプソメメトリーの計数の測定値をΨ、Δとし、前記探索範囲内の屈折率ｎ'から計算されたエリプソメメトリーの計数を
Ψ'、Δ'とし、Ｉ_S＝sin２ΨsinΔ、Ｉ_C＝sin２ΨcosΔ、Ｉ_S'＝sin２Ψ'sinΔ'、Ｉ_C'＝sin２Ψ'cosΔ'である）を評価関数とする極小探査法を用いて、前記不良と判断された屈折率を新たに計算するステップを含むことを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定プログラム（１）は、コンピュータに、エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する機能を実現させるプログラムであって、前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、反射光を測定することによって得られた反射型エリプソメトリーの計数Ψ_r、Δ_rを、記録手段から読み出して、ln｜Ｘ｜＝０（ここで、Ｘは、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r
−Ｂ₂ ^pＡ₂ ^s、ｂ＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^sＡ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−
Ｂ₁ ^pＡ₁ ^sを係数とする２次方程式ａＸ²＋ｂＸ＋ｃ＝０の根であり、ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄
膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、前記薄膜が新しく形成される側と反対側の
前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前
記薄膜を通過する光がさらに透過する前記媒体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）、又は、(｜ｃ｜²−｜ａ｜²)²−｜ｃ^*ｂ−ｃｂ^*｜²＝０（ここで、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r−Ｂ₂ ^pＡ₂ ^s、ｂ＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^s
Ａ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−Ｂ₁ ^pＡ₁ ^s、ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、前記薄膜が新しく
形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁
、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、前記薄膜が新しく形成される側と反対側の前記他
方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜
を通過する光がさらに透過する前記媒体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）を評価関数として、非線形方程式の
数値解法を用いて、メモリ上で前記波長に対する屈折率を計算する機能を実現させることを特徴としている。

また、本発明に係る膜の屈折率及び厚さの測定プログラム（２）は、コンピュータに、エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する機能を実現させるプログラムであって、前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、前記多層膜を透過した透過光を測定することによって得られた透過型エリプソメトリーの計数Ψ_t、Δ_tを、記録手段から読み出して、ln｜Ｘ｜＝０（ここで、Ｘは、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^sＡ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^sを係数とする、１次方程式ａＸ＋ｂ＝０の根であり、ρ_t＝tanΨ_texp(−ｊΔ_t)であり、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、
Ｋ^s、Ｋ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多層膜と前記薄膜との間の境界面
の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、前記薄膜が新
しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒体との境界面の
特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）
、又は、(ａ^*＋ｂ^*)―(ａ＋ｂ)＝０（ここで、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^sＡ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^s、ρ_t＝tanΨ_texp(−ｊΔ_t)であり、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｋ^s、Ｋ^p
は、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行
列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ
’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒体との境界面の特性行列
Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）を評価関
数として、非線形方程式の数値解法を用いて、メモリ上で前記波長に対する屈折率を計算する機能を実現させることを特徴としている。

また、本発明に係るコンピュータ読取可能な記録媒体（１）は、上記の膜の屈折率及び厚さの測定プログラム（１）又は（２）を記録していることを特徴としている。

上記した課題解決手段によって、多層膜の最も外側に位置する薄膜の膜厚、及び測定に使用する光の波長に応じた屈折率を求めることができる。

また、従来の方法（上記非特許文献２および３参照）と対照的に、単一波長および多波長で測定したデータに対応しており、一時的な精度の低下のために膜の屈折率及び厚さの決定が単独ではできない波長についても、精度の高い他の波長のデータで補正することにより、その波長の屈折率を決定することができる。また、従来の方法（上記非特許文献５および６参照）のように屈折率の波長依存性をあらかじめ仮定する必要がないため、異なる波長における屈折率を精度よく求めることができる。すなわち、求められた膜厚、屈折率を良否判断した結果、特定の波長に対する屈折率が不良であると判断された場合でも、その不良な屈折率を修正することができ、測定に使用する光のほぼ全ての波長に対応する屈折率を精度良く求めることができる。即ち、測定に使用する波長に応じて周期的に発生する一時的な測定精度の低下を補償することができ、従って、測定精度が低下する領域を通過した後、不良データを生じた波長において、再び良データを得ることができる。このようにして、膜厚の大きな積層膜の場合でも、新たに積層される薄膜の膜厚及び屈折率を精度良く求めることができる。

また、本発明に係る測定方法では、近似を用いることなく評価関数を導出することができるため、従来の方法（例えば、マトリックスの行列要素を２次項までのテイラー級数で近似する方法（上記非特許文献５参照））で問題であった、急激な膜厚の変化が起こった場合に計算不能に陥る危険性が無く、高精度で膜の屈折率及び厚さを決定することができる。また、評価関数の導出に使用する方程式の次数が高々２次であることから、従来の方法（例えば、８次方程式を使用する方法（上記非特許文献５参照））と比較して、数値計算に要する時間が短時間になる。従って、本発明に係る膜の屈折率及び厚さの測定方法を薄膜の製造に適用した場合、製膜工程中において、膜厚をほぼリアルタイム且つ高精度に求めることができ、製膜のリアルタイム制御が可能となり、指定の膜厚の薄膜を高精度に形成することができる。

以下に、図面を参照して、本発明の本実施の形態について説明する。先ず、本発明の原理を説明する。

図１は、本実施の形態に係る膜の屈折率及び厚さの測定方法の原理を説明する図である。外入射型配置の場合、入射光学系ＩＮおよび射出光学系ＯＵＴは、光を入射する透明な媒体（雰囲気、層０）中に配置され、多層膜（層１〜層ｍ）は光が透過する媒体（基板、層ｍ＋１）の上に形成される。入射光学系ＩＮから所定波長の入射光Ｌ_inを入射角φ₀で
多層膜の表面に照射し、その反射光Ｌ_outを射出光学系ＯＵＴで観測する。内入射型配置
の場合、入射光学系ＩＮおよび射出光学系ＯＵＴは、基板の裏側に配置され、基板の裏側にプリズムを装着し、プリズムを通して基板の裏側から光を入射する。すなわち、多層膜（層１〜層ｍ）は光を入射する媒体（基板、層０）の下に形成される。入射光学系ＩＮから所定波長の入射光Ｌ_inを入射角φ₀で多層膜の裏面に照射し、その反射光Ｌ_outを射出光学系ＯＵＴで観測する。また、測定データはリアルタイムまたは測定後まとめて演算手段（例えば、コンピュータ）ＰＣに伝送される。

図１に示すように、層０と層ｍ＋１とに挟まれたｘｙ平面に平行な多層膜（層１、２、・・・、ｍ）からなる光学系を対象とする。層０と層１との間の境界面（０１界面）の層０側における電場Ｅ（ｚ₀−０）と、ｍ（ｍ＋１）界面の層ｍ＋１側における電場Ｅ（ｚ_m＋０）とは系の光学効果を表す散乱行列Ｓによって、次式（１）又は（２）のように関係づけられる。
Ｅ(ｚ₀−０)＝ＳＥ(ｚ_m＋０) (1)

ここで、Ｅ⁺は＋（プラス）方向に進む光の電場を表し、Ｅ^-は−（マイナス）方向に進む光の電場を表す。

新しい層は雰囲気と接する多層膜の最外層の外側に形成される。雰囲気が真空または空気の場合、雰囲気から光を入射し反射光を測定する外入射型配置によって測定が行われる。すなわち、図１において、雰囲気が層０、基板が層ｍ＋１となり、新しい層は０１界面に形成される。一方、雰囲気が液体または特殊気体の場合、基板から光を入射しその反射光を測定する内入射型が測定に適している。この場合、図１において、雰囲気が層ｍ＋１、基板が層０となり、新しい層はｍ（ｍ＋１）界面に形成される。いずれの場合も、ｍ層の積層膜からなる光学系全体の散乱行列Ｓ_mは、各層および各界面の光学効果を表す行列
の積であり、式（３）のように表される。

ここで、Ｌ_jは層ｊの特性行列であり、Ｉ_j(j+1)は層ｊと層（ｊ＋１）との間のｊ（ｊ
＋１）界面の特性行列であり、それぞれ、次式（４）で表される。

ここで、exp関数の指数部分の添字でないｊは、虚数（ｊ²＝−１）を表す（以下において同じ）。また、ｒ_j(j+1)、ｔ_j(j+1)は、それぞれフレネルの反射係数、透過係数であり、層ｊの屈折率ｎ_j、入射角φ_jを使用して、Ｓ偏光、Ｐ偏光に関して、次式（５）で表される。

また、β_jは、膜厚ｄ_j（複素数であり光学的な厚さを表す）、測定波長λを使用して、次式（６）のように表される。
β_j＝２π(ｄ_j/λ)(ｎ_j ²−ｎ₀ ²sinφ₀)^1/2 (6)
今、図１において、層１（またはｍ）が形成される前の、層２〜層ｍ（又は層１〜層（ｍ−１））のｍ−１個の積層膜に、新しい層１（またはｍ）が形成されるとする。この時、新しい層、即ち層１（またはｍ）を含む系全体の散乱行列Ｓ_mは次式（７）で表すこと
ができる。

ここで、Ｘ、βは、新しく形成される層１（またはｍ）の未知の膜厚ｄ、測定波長λ（既知）、層０の屈折率ｎ₀（既知）を使用して、
Ｘ＝exp(−ｊ２β) (8)
β＝２π(ｄ/λ)(ｎ²−ｎ₀ ²sinφ₀)^1/2 (9)
である。

また、式（７）の係数Ａ_i、Ｂ_i、Ｃ_i、Ｄ_i（ｉ＝１、２）、およびＫは、新しく形成される層１（またはｍ）の未知の屈折率ｎに加え、層０における光の入射角φ₀と、既存の
層全体の光学効果を現す擬散乱行列Ｓ’_m-1の行列要素と、外入射の場合、層０と既存の
層２の屈折率、内入射の場合、既存の層ｍ−１と層ｍ＋１の屈折率によって表すことができる。

このとき、系全体の反射率Ｒおよび透過率Ｔは、散乱行列Ｓ_mの行列要素を用いて、そ
れぞれ次式（１０）、（１１）のように決定される。

Ｓ偏光とＰ偏光に対する光学系全体の反射率Ｒの比ρ_rおよび透過率Ｔの比ρ_tは、エリプソメトリーの測定によって求められる計数Ψ、Δと次式（１２）、（１３）のような関係にある。

ここで、添字のｐ、ｓはそれぞれ、Ｐ偏光、Ｓ偏光に関する値であることを表す（以下において同じ）。

反射型エリプソメトリーによる場合は、式（１２）に式（１０）を代入して整理することにより、次式（１４）のＸについての２次方程式が導かれる。
ａＸ²＋ｂＸ＋ｃ＝０ (14)
ここで、ａ、ｂ、ｃは、

である。

一方、透過型エリプソメトリーによる場合は、式（１３）に式（１１）を代入し整理することにより、式（１６）のＸについての１次方程式が導かれる。
ａＸ＋ｂ＝０ (16)
ここで、ａ、ｂは、

である。

式（１５）および式（１７）に示す係数ａ、ｂ、ｃは未知数として新しい層１（又はｍ）の屈折率ｎのみを含む。そこで、屈折率ｎを決定するために、光を入射する層０と新しく形成される層１（またはｍ）は透明という仮定を用いる。その場合、層０の屈折率ｎ₀
および新しい層１（又はｍ）の屈折率ｎは実数となり、新しい層１（又はｍ）の膜厚ｄも幾何学的な長さを表す実数となる。従って、式（８）の指数部は虚数となり、Ｘの絶対値は１に等しくなり、次式（１８）が成り立つ。
ln｜Ｘ｜＝０ (18)
さらに、この場合は、反射型エリプソメトリーにおいて、式（１４）の複素共役は式（１９）のように表され、式（１４）と式（１９）が同一の解を有することから、式（２０）が成り立つ。
ｃ^*Ｘ²＋ｂ^*Ｘ＋ａ^*＝０ (19)
Ｇ＝(｜ｃ｜²−｜ａ｜²)²−｜ｃ^*ｂ−ｃｂ^*｜²＝０ (20)
また、透過型エリプソメトリーの場合は、式（１６）の複素共役は式（２１）のように表され、式（１６）と式（２１）が同一の解を有することから、式（２２）が成り立つ。ｂ^*Ｘ＋ａ^*＝０ (21)
Ｇ＝(ａ^*＋ｂ^*)−(ａ＋ｂ)＝０ (22)
従って、反射型エリプソメトリーの場合は式（１８）または（２０）を用い、透過型エリプソメトリーの場合は式（１８）または（２２）を用いることによって新たな層の屈折率ｎ（実数）を決定し、式（８）および（９）を用いて膜厚ｄ（実数）を求めることができる。

次に、図２、３に示した、フローチャートに基づいて、本実施の形態に係る膜の屈折率及び厚さの測定方法を説明する。なお、ここでは簡単のために、外入射型配置とする。

ここで、図１に示したように、入射光学系ＩＮから所定波長の光を入射角φ₀で多層膜
の表面に照射し、その反射光を射出光学系ＯＵＴで観測する。また、既に、基板層（層（ｍ＋１））の上に層ｍ〜層２が形成されており、これらの光学的特性が既知であり、層２の上に新たに形成される層１の屈折率及び膜厚が未知であると仮定する。即ち、層１の製膜を行う前に既に存在している既存の層を含む光学系の擬散乱行列Ｓ'_m-1は、式（３）を用いて、層ｍから開始して順にその時点の最外層である層２まで、各層の光学効果を現す伝達行列Ｌ_jＩ_j(j+1)を積算することによって、求められていると仮定する。

また、図２、３に示したフローチャートは、主に測定データを演算手段ＰＣで、即ち、演算手段ＰＣ内のＣＰＵで処理することによって実現される。従って、以下の説明においては、特に断らない限り、ＣＰＵが行うものとして説明する。また、各処理において、ＣＰＵは、記録手段に予め記録された測定データなどを読み出し、メモリに書き込み、メモリ上のワーク領域で所定の処理を行い、処理結果を、必要に応じてメモリ、記録手段などに記録し、モニターなどの表示手段に提示する。

先ず、ステップＳ１において、上記したように、雰囲気中（層０）に配置された入射光学系ＩＮ及び射出光学系ＯＵＴを用いて、所定の波長λ_iでのエリプソメトリーの計数Ψ_i、Δ_iを、複数の波長λ_i（ｉ＝１〜Ｎ）に関して測定する。測定されたエリプソメトリーの計数Ψ_i、Δ_i（ｉ＝１〜Ｎ）、及びλ_i（ｉ＝１〜Ｎ）は、対応させて記録手段に記録
させることができる。また、測定後にまとめて演算を行う場合には、記録手段に記録した複数の波長λ_i（ｉ＝１〜Ｎ）におけるエリプソメトリーの計数Ψ_i、Δ_iを読み出す。

ステップＳ２において、複数の波長に関して計算するための繰り返し処理のカウンタｋ1に初期値として１を設定する。

ステップＳ３において、層１の屈折率の初期値ｎ（０）を設定し、繰り返し処理のカウンタｋ2に１を設定する。層１の屈折率の初期値ｎ（０）は任意であるが、約２．０を使
用することが望ましい。

ステップＳ４において、上記で説明した、反射型エリプソメトリーに関する式（１２）、（１４）、（１５）の条件の下で、（１８）を評価関数として、Newton-Raphson法を用いて、波長λ_k1に対する屈折率ｎ（ｋ2）（ここでは、繰り返しの初回でありｋ2＝１）を計算する。式（１４）、（１８）から、評価関数をより明示的に記載すれば、次式（２３
）となる。

ここで、ａ、ｂ、ｃは、式（１５）で定義されており、式（１５）の中で、ρ_rは、ステ
ップＳ１で測定され、または記録手段から読み出されたエリプソメトリーの計数Ψ_i、Δ_i（ｉ＝１〜Ｎ）を使用して、求めることができる。また、係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、新しく形成される層１を含む系の散乱行列Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1の行列要素として求めることができる。

上記したように、ａ、ｂ、ｃは、未知数として、層１の屈折率ｎ（ｋ2）（ここで、ｋ2＝１）のみを含む。従って、ｆ（ｘ）＝ln｜Ｘ（ｘ）｜、ｘ＝ｎ（ｋ2）（ｋ2＝１）として、次式（２４）によって、ｎ（０）、及びｎ（０）を微小変化させた値から新たな屈折率ｎ（１）を計算する。

ステップＳ５において、ステップＳ４において計算された屈折率ｎ（ｋ2）とｎ（ｋ2−１）との差が所定の値δｎよりも小さいか否かを判断することによって、計算の収束性を判断する。ここでは、ｋ2＝１であり、ｎ（１）−ｎ（０）＜δｎ か否かを判断する。
ｎ（１）−ｎ（０）＜δｎ でないと判断すれば、ステップＳ６に移行する。δｎは、例えば、０．００００１に設定することができる。

ステップＳ６において、繰り返し回数を制限するために、カウンタｋ2が所定の値より
も小さいか否かを判断する。ここでは、繰り返し回数の上限を２０とし、ｋ2が２０以上
であると判断すれば、繰り返し処理を止めて、計算時間短縮のためステップＳ１０に移行する。ｋ2が２０よりも小さいと判断すれば、ステップＳ７に移行する。

ステップＳ７において、カウンタｋ2を１だけ増大させて、ｎ（０）にステップＳ４で
計算した屈折率ｎ（ｋ2）（ｋ2＝１）を設定して、ステップＳ４に戻る。その後、繰り返し回数が上限値を超えない範囲で、ｎ（ｋ2）−ｎ（ｋ2−１）＜δｎ となるまで、ステップＳ４〜Ｓ７を繰り返す。

ステップＳ８において、以上の繰り返し処理によって求められた屈折率ｎ（ｋ2）を、
波長λ_k1に対する屈折率ｎ_k1として、記録手段に記録する。

ステップＳ９において、式（８）、（９）から得られる次式（２５）に、波長λ_k1に対する屈折率ｎ_k1を代入して、対応する層１の膜厚ｄ_k1を求める。
ｄ＝Re[ｊ(１/４π)λ(ｎ²−ｎ₀ ²sinφ₀)^-1/2lnＸ] (25)
ここで、Reは｛ ｝内の複素数の実数部分を表す。

ステップＳ１０において、繰り返しのカウンタｋ１が測定したデータ数Ｎ（波長λ_iの
数）よりも大きいか否かを判断し、ｋ１＜Ｎであると判断すれば、ステップＳ１１において、カウンタｋ１を１だけ増大させた後、ステップＳ３に戻る。これによって、測定データ全てに関して、ステップＳ３〜Ｓ１１の処理を繰り返し、ｉ＝１〜Ｎの各々の波長λ_i
に対する屈折率ｎ_i、膜厚ｄ_iを求める。

次に、ステップＳ１２において、以上で求められた屈折率ｎ_i、膜厚ｄ_i（ｉ＝１〜Ｎ）の良否を判断し、不良データの修正を行う。測定されるエリプソメトリーの計数Ψ、Δは、ある膜厚で周期的に精度が大きく低下することが知られている。そのため、複数の波長で測定を行うのが一般的であり、各々の波長λ_iに関する測定データについて良否を判断
することが重要である。測定波長λ_iが異なっていても膜厚は一定であることを基本原理
として、計算結果の良否を判断する方法を、以下に説明する。

図３は、ステップＳ１２における計算結果である屈折率、膜厚の組（ｎ_i、ｄ_i）（ｉ＝１〜Ｎ）の良否判断及び不良データの修正に関する処理を説明するフローチャートである。

ステップＳ２１において、後述する分類による判断を容易とするために、予備的に計算結果（ｎ_i、ｄ_i）（ｉ＝１〜Ｎ）を選別する。即ち、Ｓ４〜Ｓ７の繰り返し処理において、所定の回数以内に収束した（ステップＳ６でｋ2≧２０と判断されなかった）データの
組のうち、屈折率ｎ_i、膜厚ｄ_i（ｉ＝１〜Ｎ）について、ｎ_min＜ｎ_i＜ｎ_max、及びｄ_min＜ｄ_i＜ｄ_maxが成り立つ場合に、そのデータの組を「良」、即ち、物理的に妥当なデータとする。ここで、膜厚ｄの単位はｎｍ（ナノメートル：以下において同じ）であり、ｎ_min、ｎ_maxは、それぞれ屈折率の許容範囲の下限値、上限値である。真空の屈折率が１．０であり、透明材料の屈折率は高々２．５程度であることから、例えば、ｎ_min＝１．０、
ｎ_max＝３．０とすることができる。また、ｄ_min、ｄ_maxは、膜厚の許容範囲の下限値、
上限値である。これらの値は、物理的に意味のある限界値もしくは測定対象の薄膜の製膜条件で到達不可能である値を設定する。例えば、ｄ_min＝０（ｎｍ）、ｄ_max＝１０（ｎｍ）と設定することができる。

ステップＳ２１での予備選別で「良」と判断されたデータの組（ｎ_i、ｄ_i）について、さらにデータの分類と最適グループの決定を行うために、先ずステップＳ２２において、グループ数Ｎ_Gを決定する。ステップＳ２１で良データと選別された組の数をＮ_goodとし
て、グループ数Ｎ_Gを、例えば次のように決定する。
Ｎ_good＜５の場合、Ｎ_G＝floor（Ｎ_good／２）
５≦Ｎ_good＜８の場合、Ｎ_G＝floor（Ｎ_good／２）＋２
８≦Ｎ_good＜１０の場合、Ｎ_G＝floor（Ｎ_good／２）＋３
１０≦Ｎ_good＜１３の場合、Ｎ_G＝floor（Ｎ_good／２）＋４
Ｎ_good≧１３の場合、Ｎ_G＝floor（Ｎ_good／２）＋５
ここで、floor()は（ ）内の数値の小数部を切り捨てる演算を表す。

続くステップＳ２３〜２５において、Ｎ_good個の良データの組（ｎ_i、ｄ_i）を、膜厚ｄ_iの対数値（lnｄ_i）によってグループ化する。

ステップＳ２３において、ステップＳ２２で決定されたグループ数Ｎ_Gを用いて、各々
のグループに属する膜厚ｄ_iの分散値（σ²）の総和（Σσ²）が最小になるように、良デ
ータの組（ｎ_i、ｄ_i）をＮ_G個のグループに分類する。ここで、１組の良データのみを含
むグループについては、σ²＝０とする。このようにグループに分類する方法は周知であ
り、ここでは説明を省略する。本ステップでの処理の結果、良データの組（ｎ_i、ｄ_i）の
中で、他のデータとかけ離れた値が孤立、即ち単独でグループを形成することになる。

ステップＳ２４において、ステップＳ２３で分類したグループから、仮の最適グループＧ_iの群を決定する。先ず、要素の数が２以上かつ最大である１つのグループを仮最適グ
ループＧ₁とし、Ｇ₁を除く要素数が２以上かつ最大のグループを仮最適グループＧ₂とす
る。このとき、要素数が等しくかつ最大のグループが複数ある場合には、そのうち１つのグループを選択し仮最適グループＧ₁とする。また、Ｇ₁を除くグループのうち、要素数が等しくかつ最大のグループが複数ある場合には、それら全てを仮最適グループとし、Ｇ₂
、Ｇ₃、…と順序付ける。また、Ｇ₁を除くグループのうち、要素数が２以上のグループがない場合は仮最適グループはＧ₁のみとする。なお、要素数２以上のグループが少なくと
も１つ存在できるように、グループ数Ｎ_GがステップＳ２２で決定される。

次に、すべての仮最適グループＧ₁、Ｇ₂、Ｇ₃、…について、すべての良データを対象
とし、仮最適グループＧ_iの要素ｄ_ijの対数値lnｄ_ijの平均値(lnｄ)_ai（添字ａは平均値
を表す）及び標準偏差σを用いて計算した規格化標準偏差（σ／(lnｄ)_ai）をσ_dNiとし
て、σ_dNiが大きくなるほど、良データの組に属する膜厚の対数値lnｄ_kと(lnｄ)_aiとの差の絶対値abs（lnｄ_k−(lnｄ)_ai）の許容値が大きくなるように、次の条件に合致する良データの組（ｎ_k、ｄ_k）を、仮最適グループＧ₁、Ｇ₂、Ｇ₃、…に取り込む。例えば、次の
ように条件を定める。
σ_dNi＝０．０ の場合、abs（lnｄ_k−(lnｄ)_ai）＜０．０２
σ_dNi＜０．０００００００２ の場合、abs（lnｄ_k−(lnｄ)_ai）＜２０００００σ_dni
σ_dNi＜０．００００００２ の場合、abs（lnｄ_k−(lnｄ)_ai）＜３００００σ_dni
σ_dNi＜０．０００００２ の場合、abs（lnｄ_k−(lnｄ)_ai）＜５０００σ_dni
σ_dNi＜０．０００００６ の場合、abs（lnｄ_k−(lnｄ)_ai）＜３０００σ_dni
σ_dNi＜０．００００２ の場合、abs（lnｄ_k−(lnｄ)_ai）＜１０００σ_dNi
σ_dNi＜０．００００５ の場合、abs（lnｄ_k−(lnｄ)_ai）＜３００σ_dNi
σ_dNi＜０．０００２ の場合、abs（lnｄ_k−(lnｄ)_ai）＜１００σ_dNi
σ_dNi＜０．０００５ の場合、abs（lnｄ_k−(lnｄ)_ai）＜４０σ_dNi
σ_dNi＜０．００２ の場合、abs（lnｄ_k−(lnｄ)_ai）＜２０σ_dNi
σ_dNi≧０．００２ の場合、abs（lnｄ_k−(lnｄ)_ai）＜１０σ_dNi
ここで、abs()は（ ）内の数値の絶対値を求める演算を表す。すべての良データについ
ての（lnｄ_k−(lnｄ)_ai）²の総和が最小になるように、この処理をＧ₁、Ｇ₂、Ｇ₃、…の
順に行う。なお、ここで示している条件の数および許容値は一例であり、エリプソメトリーの計数Ψ、Δを測定した条件に応じて異なる条件の数および許容値を使用することができる。また、許容値の決定法もここで示している方法に限らず、例えば、規格化標準偏差σ_dNiと許容値とを関係づける適当な数式によって算出しても良い。

次に、仮最適グループ群へのデータの取込処理によって変化した要素数に応じて、仮最適グループ群Ｇ_iのうち要素数最大のグループおよび２番目のグループを決定し、それぞ
れ最適グループＧ₁および準最適グループＧ₂とする。このとき、要素数が２以上で等しくかつ最大のグループが複数ある場合には、そのうち任意の１つのグループを選択し最適グループＧ₁とする。また、Ｇ₁を除くグループのうち、要素数が等しくかつ最大のグループが複数ある場合には、そのうち任意の１つのグループを選択し準最適グループＧ₂とする
。また、仮最適グループ群がＧ₁のみである場合は準最適グループＧ₂は決定しない。

さらに、最適および準最適グループＧ_iの要素数をＮ（Ｇ_i）、規格化標準偏差をσ_dN（Ｇ_i）（ｉ＝１又は２）とし、「Ｎ（Ｇ₁）＝Ｎ（Ｇ₂）」および「σ_dN（Ｇ₂）＜σ_dN（Ｇ₁）」の条件を両方満たす場合、最適グループＧ₁と準最適グループＧ₂とを入れ替える処
理を行い、最終的に最適グループＧ₁を決定する。

ステップＳ２５において、再度、データの良否を判断する。まず、ステップＳ２４で最終的に決定された最適グループＧ₁の要素の屈折率データｎ_jから、膜厚データｄ_jに関し
て上記で行ったのと同様に、屈折率データｎ_jの平均値ｎ_a及び標準偏差σを用いて計算した規格化標準偏差σ_nN（＝σ／ｎ_a）から、δ_n＝４σ_nNによって許容偏差δ_nを決定し、
ステップＳ２１で良と判断されたデータの組の屈折率データに関して、abs（ｎ_k−ｎ_a）
＜δ_nを満たすデータの組（ｎ_k、ｄ_k）を選別する。

次に、膜厚に関して、ステップＳ２４で決定された最適グループＧ₁の要素の膜厚デー
タｄ_jから、上記で行ったように、膜厚データｄ_jの対数値の平均値(lnｄ)_a及び標準偏差
σを用いて計算した規格化標準偏差σ_dN（＝σ／(lnｄ)_a）から、次のように許容偏差δ_dを決定し、屈折率データに関して選別したデータの組の膜厚データに関して、abs（lnｄ_k−(lnｄ)_a）＜δ_dを満たすデータの組（ｎ_k、ｄ_k）を良データとする。
σ_dN＜０．０００００００２ の場合、δ_d＝２０００００σ_dni
σ_dN＜０．００００００２ の場合、δ_d＝３００００σ_dni
σ_dN＜０．０００００２ の場合、δ_d＝５０００σ_dN
σ_dN＜０．０００００６ の場合、δ_d＝３０００σ_dni
σ_dN＜０．００００２ の場合、δ_d＝１０００σ_dN
σ_dN＜０．００００５ の場合、δ_d＝３００σ_dN
σ_dN＜０．０００２ の場合、δ_d＝１００σ_dN
σ_dN＜０．０００５ の場合、δ_d＝４０σ_dN
σ_dN＜０．００２ の場合、δ_d＝２０σ_dN
σ_dN≧０．００２ の場合、δ_d＝１０σ_dN
なお、ここで示している条件の数および許容偏差δ_dは一例であり、エリプソメトリーの
計数Ψ、Δを測定した条件に応じて異なる条件の数および許容偏差δ_dを使用することが
できる。また、許容偏差δ_dの決定法もここで示している方法に限らず、例えば、規格化
標準偏差σ_dNと許容偏差δ_dとを関係づける適当な数式によって算出しても良い。

ステップＳ２６〜Ｓ３７において、測定波長が異なっていても膜厚は一定であることを基本原理として、ステップＳ２１及びＳ２５で不良と判断されたデータを修正する。

ステップＳ２６において、ステップＳ２５で決定された良データ（ｄ_gi、ｎ_gi）から、各々の平均膜厚ｄ_ga、平均屈折率ｎ_gaを計算する。また、平均屈折率ｎ_gaから屈折率の探索範囲の上限を（１．０＋ｒ）ｎ_ga、下限を（１．０−ｒ）ｎ_gaに決定する。ここで、ｒは探索範囲の相対値であり、例えばその値を０．１とする。

ステップＳ２７において、不良データの総個数Ｎ_eを求め、後述する計算の収束性を判
断するための閾値ｆｓを、例えば０．００００１に設定し、繰り返し処理のカウンタｋ3
の初期値として１を設定する。

ステップＳ２８において、繰り返し処理のカウンタｋ4の初期値として１を設定する。

ステップＳ２９〜Ｓ３４において、ステップＳ２６で決定された平均膜厚ｄ_gaを用いて、非線形関数の極小探索法によって、前記探索範囲（０．９ｎ_ga〜１．１ｎ_ga）内で探索して、ステップＳ２１及びＳ２５で不良と判断された屈折率を新たに計算する。ここでは一例として、非線形関数の極小探索法としてBrent法を用いる場合を説明する。

次式（２６）の評価関数を使用する。

ｆ（ｎ）＝(Ｉ_S'−Ｉ_S)²−(Ｉ_C'−Ｉ_C)² (26)
ここで、Ｉ_sおよびＩ_cは、
Ｉ_S＝sin２ΨsinΔ (27)
Ｉ_C＝sin２ΨcosΔ (28)
であり、測定で得られる値である。Ｉ_sおよびＩ_cはΨおよびΔよりも周期的な精度の低下が著しくないと考えられる。Ｉ_s'及びＩ_c'を計算する方法は、下記のステップＳ３１で説明する。

ステップＳ２９において、探索範囲（０．９ｎ_ga〜１．１ｎ_ga）の範囲内で新たな屈折率ｎ（ｋ4）を、使用する極小探索法に特徴的な方法で計算する。

ステップＳ３０において、測定により得られたエリプソメトリーの計数Ψ、Δから式（２７）および（２８）を用いて、Ｉ_sおよびＩ_cを計算する。

ステップＳ３１において、ステップＳ２６で計算された膜厚平均値ｄ_ga、ステップＳ２９で計算された新しいｎ（ｋ4）、およびステップＳ４で計算されたＳ’_m-1を用いて、式（７）で表されるＳ_mの行列要素を算出し、式（１０）に代入して系全体の反射率ρ_rを求め、式（１２）からΨおよびΔを計算し、式（２７）および（２８）を用いてＩ_s'およびＩ_c'を計算する。計算されたＩ_s'、Ｉ_c'、及びステップＳ３０で計算されたＩ_s、Ｉ_cを使用して、式（２６）からｆ（ｎ）を計算する。

ステップＳ３２において、ステップＳ３１で計算されたｆ（ｎ）が、閾値ｆｓよりも小さいか否かを判断する。ｆ（ｎ）＜ｆｓと判断された場合、繰り返し処理を終了してステップＳ３５に移行する。ｆ（ｎ）≧ｆｓと判断された場合、ステップＳ３３に移行する。

ステップＳ３３において、繰り返し処理のカウンタｋ4が所定の値よりも小さいかを判
断する。ここでは、繰り返し回数の上限を２０とし、ｋ4≧２０の場合、繰り返し処理を
終了してステップＳ３６に移行するが、ｋ4＜２０の場合、ステップＳ３４に移行し、繰
り返しカウンタｋ4を１だけ増大させた後、ステップＳ２９に戻る。その後、繰り返しの
上限回数を超えない範囲内で、ｆ（ｎ）＜ｆｓと判断されるまで、ステップＳ２９〜Ｓ３４の処理を繰り返す。

ステップＳ３５において、以上の処理で求められたｎ（ｋ4）を、不良と判断されてい
た屈折率データｎ_ek3に設定する。

ステップＳ３６において、不良と判断されたデータの総数Ｎ_eに関して計算されたか否
かを判断するために、カウンタｋ3がＮ_eよりも小さいか否かを判断する。ｋ3＜Ｎ_eと判断すれば、ステップＳ３７に移行して、繰り返しカウンタｋ3を１だけ増大させた後、ステ
ップＳ２８に戻り、ｋ3≧Ｎ_eと判断されるまでステップＳ２８〜Ｓ３７を繰り返し、不良と判断されたデータ全てに関して上記した処理を行う。

以上で、測定に使用する波長によっては、測定精度が低下し、不良と判断されていた屈折率ｎ_eを精度良く、新たに決定することができる。通常の場合、不良データの修正はほ
とんどの場合成功し、光の波長に応じた屈折率を決定することができる。

以上のように、ステップＳ１２での計算結果の良否判断及び不良データの修正の後、ステップＳ１３において、終了の指示があったか否かを判断し、終了の指示が無ければ、ステップＳ１４に移行する。

ステップＳ１４では、新たに擬似散乱行列Ｓ'_mを計算した後、ステップＳ１に戻る。即ち、以上で決定された良データの膜厚平均値ｄ_ga、屈折率ｎ_i（良データ及び修正された
不良データ）、およびＳ'_m-1を用い、新しい層の光学効果を表す伝達行列Ｌ_jＩ_j(j+1)を
Ｓ’_m-1に積算して、Ｓ'_mを計算する。その後、ステップＳ１に戻る。

以上のように、ステップＳ２〜Ｓ１１の処理によって、最上位層の膜厚、及び測定に使用した光の波長に応じた屈折率を求めることができる。

また、ステップＳ２〜Ｓ１１の処理によって求められた複数組の膜厚、屈折率が、ステップＳ１２において良否判断された結果、特定の波長のデータ（膜厚、屈折率）が不良データであると判断された場合でも、同ステップＳ１２において、その不良な屈折率データを修正することができ、全ての測定波長に対応する屈折率を精度良く求めることができる。このようにして、従来の方法（上記非特許文献２および３参照）と対照的に、単一波長および多波長で測定したデータに対応しているため、一時的な精度の低下のために膜の屈折率及び厚さの決定が単独ではできない波長についても、精度の高い他の波長のデータで補正することにより、その波長の屈折率を決定することができる。また、従来の方法（上記非特許文献５および６参照）のように屈折率の波長依存性をあらかじめ仮定する必要がないため、異なる波長における屈折率を精度よく求めることができる。即ち、測定波長に応じて周期的に発生する一時的な測定精度の低下を補償することができ、従って、測定精度が低下する領域を通過した後、不良データを生じた波長において、ステップＳ２〜Ｓ１１の処理により再び良データを得ることができる。このようにして、膜厚の大きな積層膜の場合でも、新たに積層される薄膜の膜厚及び屈折率を精度良く求めることができる。

また、本発明の方法では、近似を用いることなく評価関数を導出することができるため、従来の方法（例えば、マトリックスの行列要素を２次項までのテイラー級数で近似する方法（上記非特許文献５参照））で問題であった、急激な膜厚の変化が起こった場合に計算不能に陥る危険性が無く、高精度で膜の屈折率及び厚さを決定することができる。また、方程式の次数が高々２次であることから、従来の方法（例えば、８次方程式を使用する方法（上記非特許文献５参照））と比較して計算に要する時間が短時間でよく、従って、本発明の方法を製膜工程に適用した場合、製膜中に膜厚をほぼリアルタイム且つ高精度に求めることができることから、指定膜厚の薄膜を高精度に形成することができる。

図４は、シミュレーションで求めたエリプソメトリーの計数Ψ、Δを測定データの代りに使用して、本実施の形態に係る膜の屈折率及び厚さの測定方法で、屈折率および膜厚を求めた結果を示す図である。

シミュレーション条件として、各々の薄膜１層の厚さが２ｎｍ、全体で１２５層（全体の厚さが２５０ｎｍ）からなる多層膜とし、積層された薄膜の屈折率が、所定の境界を挟んで階段状に変化するとした。また、５種類の波長（６００、６５０、７００、７５０、８００ｎｍ）を想定した。なお、シミュレーションで計算したΨ、Δには、ガウス分布を有する標準偏差０．００１°のランダムノイズを含ませている。このノイズの大きさは現在市販されている多くのエリプソメーターの性能に匹敵する値である。

図４のグラフは、縦軸を逆算した屈折率、横軸を逆算した膜厚として描いている。薄膜の屈折率が、所定の境界で階段状に変化している状態を良く再現できていることがわかる。

図５は、シミュレーションで求めたエリプソメトリーの計数Ψ、Δを測定データの代りに使用して、本実施の形態に係る膜の屈折率及び厚さの測定方法で、屈折率および膜厚を求めた結果を示す図である。

シミュレーション条件として、各々の薄膜１層の厚さが１ｎｍ、全体で２００層（全体の厚さが２００ｎｍ）からなる多層膜とし、各薄膜の屈折率が、積層される順序に応じて
、下層から上層に向かって一定勾配で増大するとした。また、３種類の波長（６００、７００、８００ｎｍ）を想定した。なお、シミュレーションで計算したΨ、Δには、ガウス分布を有する標準偏差０．００１°のランダムノイズを含ませている。

図５のグラフは、縦軸を逆算した屈折率、横軸を逆算した膜厚として描いている。薄膜の屈折率が、積層される順序に応じて、下層から上層に向かって一定勾配で増大している状態を良く再現できていることが分かる。

本実施の形態の方法での計算結果と、シミュレーションデータとを比較した計算精度は、図４の結果について、屈折率に関する平均絶対偏差が０．０３％以下、標準偏差が０．４％以下、膜厚に関する平均絶対偏差が−０．０１％以下、標準偏差が０．３％以下であった。また、クロック周波数が５００ＭＨｚの３２ビットＣＰＵ（PENTIUM（登録商標）
）を搭載したコンピュータを使用した図４についての計算では、１つの波長で逆算により屈折率および膜厚を求めるステップＳ４〜Ｓ９の繰り返し計算において、１回の数値計算に要した平均繰り返し回数は８．５回であり、平均時間は０．７２ｍｓであった。また、平均膜厚を用いて逆算結果が不良と判断された屈折率を新たに計算するステップＳ２９〜Ｓ３５の繰り返し計算において、１回の数値計算に要した平均繰り返し回数は１７．２回であり、平均時間は２．６６ｍｓであった。また、１つの最外層についての総所要時間の平均は１７．５ｍｓであった。従って、本実施の形態に係る方法を製膜工程に適用すれば、膜形成中の膜厚をほぼリアルタイムで測定することができ、膜厚の高精度の制御が可能となる。

以上の説明において、薄膜が新しく形成される側の透明媒体と多層膜との界面に光を入射する外入射型配置としたが、薄膜が新しく形成される側と反対側の透明媒体と多層膜との界面に光を入射する内入射型配置の場合も、系全体の散乱行列Ｓ_mが当該薄膜を除いた
多層膜と当該薄膜との間の境界面の特性行列Ｉ_(m-1)m、当該薄膜の特性行列Ｌ_m、当該薄
膜と当該薄膜を通過する光がさらに透過する媒体との境界面の特性行列Ｉ_m(m+1)、及び当該薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算されることを勘案し、同様に計算できる。

また、計算に用いるエリプソメトリーの計数はΨとΔとしたが、換算によってΨとΔが得られるのであれば、計算に用いるエリプソメトリーの計数の種類は特に限定されない。例えば、式（２７）および（２８）で表されるＩ_SおよびＩ_Cが測定値として得られる場合には、式（２７）および（２８）を用いてΨとΔを逆算し、これを用いて同様に計算を行うことができる。

また、ステップＳ４（図２）での屈折率の数値計算における評価関数として式（１８）を使用する場合を説明したが、評価関数として式（２０）を使用してもよい。

また、反射光を測定する反射型エリプソメトリーの場合を説明したが、多層膜を透過した透過光を測定する透過型エリプソメトリーの場合にも、評価関数として式（１８）又は式（２２）を使用して、同様に屈折率及び膜厚を計算することができる。この場合、方程式は、式（１４）の代りに、式（１７）を係数とする式（１６）を使用する。

また、繰り返し処理における収束性の判断は、上記した方法に限定されず、例えば、評価関数の値が設定した閾値よりも小さくなることを収束条件としてもよい。

また、ステップＳ４（図２）における数値計算は、Newton-Laphson法に限定されず、種々の数値解法を使用することができる。

また、ステップＳ２９（図３）における極小探索法は、Brent法に限定されず、種々の
極小探索法を使用することができる。

また、以上では、複数の波長で測定したエリプソメトリーの計数を処理対象とする場合を説明したが、ある測定波長で測定精度の低下が起こらないか、あるいは起こっても著しくないことが分かっている場合には、１つの波長でのみ測定したエリプソメトリーの計数を使用して、ステップＳ２〜Ｓ１０での処理を行う（Ｎ＝１の場合）ことによって、屈折率及び膜厚を、高速に且つ高精度に求めることができる。但しこの場合には、ステップＳ１２の処理は行わない。

また、図２、３に示したフローチャートの処理順序に限定されず、本発明の効果を生じる範囲内で、種々の変更、追加が可能である。

図６は、本実施の形態に係る膜の屈折率及び厚さの測定方法を、多層膜の製造に応用する場合の処理を示すフローチャートである。ステップＳ４１において、最初の第１層を形成する作業を行う。

製膜作業を中断することなく、図２、３に関して上記したのと同様に、ステップＳ４２において、複数の波長の光を使用してエリプソメトリーの計数を測定し、ステップＳ４３において、波長毎の屈折率及び膜厚を計算し、ステップＳ４４において、計算された屈折率、膜厚の良否判断及び不良データの修正を行い、良好な膜厚の平均値を求めると同時に、全ての波長に対応する屈折率を求める。ステップＳ４５において、ステップＳ４４で得られた良好な膜厚の平均値が、所望の値であるか否かを判断する。所望の膜厚でなければ、所望の膜厚になるまで、ステップＳ４６での製膜工程の制御、及びステップＳ４２〜Ｓ４５を繰り返す。

さらに、ステップＳ４７において、新たな薄膜を形成するか否かを判断し、新たな薄膜が形成される場合、ステップＳ４１に戻り、新たに形成される膜厚の形成を制御するために、上記した測定及び計算処理を繰り返す。これによって、所望の膜厚の層を複数備えた多層膜を製造することができる。

本発明の実施の形態に係る膜の屈折率及び厚さの測定方法の原理を説明する図である。 本発明の実施の形態に係る膜の屈折率及び厚さの測定方法を示すフローチャートである。 本発明の実施の形態に係る膜の屈折率及び厚さの測定方法における計算結果の良否判断及び不良データの修正処理を示すフローチャートである。 本発明の実施の形態に係る膜の屈折率及び厚さの測定方法によって計算した屈折率と厚さの関係を示す図である。 本発明の実施の形態に係る膜の屈折率及び厚さの測定方法によって計算した屈折率と厚さの関係を示す図である。 本発明の実施の形態に係る膜の屈折率及び厚さの測定方法を製膜工程に応用する場合を示すフローチャートである。

符号の説明

ＩＮ 入射光学系
ＯＵＴ 射出光学系
ＰＣ 演算手段
Ｌ_in 入射光
Ｌ_out 反射光
φ₀、φ_i 入射角
０ 雰囲気（外入射型配置）または基板（内入射型配置）
１〜ｍ 薄膜
ｍ＋１ 基板（外入射型配置）または雰囲気（内入射型配置）

Claims (12)

1. エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する方法であって、
   前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、反射光を測定することによって得られる反射型エリプソメトリーの計数Ψ_r、Δ_rを使用し、
   ln｜Ｘ｜＝０
   （ここで、Ｘは、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r−Ｂ₂ ^pＡ₂ ^s、ｂ＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^sＡ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−Ｂ₁ ^pＡ₁ ^sを係数とする２次方程式ａＸ²＋ｂＸ＋ｃ＝０
   の根であり、
   ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、
   係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記
   薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、
   前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
   体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
   から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）、又は、
   (｜ｃ｜²−｜ａ｜²)²−｜ｃ^*ｂ−ｃｂ^*｜²＝０
   （ここで、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r−Ｂ₂ ^pＡ₂ ^s、ｂ＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^sＡ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−Ｂ₁ ^pＡ₁ ^s、ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、
   係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記
   薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、
   前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
   体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
   から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）
   を評価関数として、非線形方程式の数値解法を用いて、前記波長に対する屈折率を計算する屈折率計算ステップを含むことを特徴とする膜の屈折率及び厚さの測定方法。
2. エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する方法であって、
   前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、前記多層膜を透過した透過光を測定することによって得られる透過型エリプソメトリーの計数Ψ_t、Δ_tを使用し、
   ln｜Ｘ｜＝０
   （ここで、Ｘは、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^sＡ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^sを係数とする、１次方程式ａＸ＋ｂ＝０の根であり、
   ρ_t＝tanΨ_texp(−ｊΔ_t)であり、
   係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｋ^s、Ｋ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多
   層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素
   として求められ、
   前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
   体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
   から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）、又は、
   (ａ^*＋ｂ^*)―(ａ＋ｂ)＝０
   （ここで、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^sＡ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^s、ρ_t＝tanΨ_texp(−ｊ
   Δ_t)であり、
   係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｋ^s、Ｋ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多
   層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素
   として求められ、
   前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
   体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
   から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）
   を評価関数として、非線形方程式の数値解法を用いて、前記波長に対する屈折率を計算する屈折率計算ステップを含むことを特徴とする膜の屈折率及び厚さの測定方法。
3. 前記屈折率計算ステップによって計算された前記屈折率を用いて、
   ｄ＝Re｛ｊ(1/4π)λ(ｎ²−ｎ₀ ²sinφ₀)^-1/2lnＸ｝
   （ここで、ｎは前記薄膜の屈折率、ｎ₀は前記媒体の屈折率、λは前記波長、φ₀は前記光の入射角である）
   によって、前記薄膜の膜厚ｄを計算する膜厚計算ステップをさらに含むことを特徴とする請求項１又は２に記載の膜の屈折率及び厚さの測定方法。
4. 複数の波長の各々に対して、前記屈折率計算ステップによって前記屈折率を計算し、且つ前記膜厚計算ステップによって前記膜厚を計算し、
   計算された複数の前記屈折率及び前記膜厚の良否を判断し、精度良く前記薄膜の膜厚を決定する良否判断ステップをさらに含むことを特徴とする請求項３に記載の膜の屈折率及び厚さの測定方法。
5. 前記良否判断ステップの後に、前記良否判断ステップにおいて不良と判断された屈折率を新たに計算する修正ステップをさらに含むことを特徴とする請求項４に記載の膜の屈折率及び厚さの測定方法。
6. 前記良否判断ステップが、
   同じ波長に対する、前記屈折率計算ステップによって計算された前記屈折率、及び前記膜厚計算ステップによって計算された前記膜厚からなる複数の組の中から、良好な組を選別する選別ステップと、
   前記良好な組を、所定のグループ数のグループに分類する分類ステップと、
   分類された前記グループを用いて、信頼性の高い膜厚を含む最適グループを決定する最適グループ決定ステップと、
   前記最適グループに含まれる膜厚の平均値及び規格化標準偏差を用いて、前記良好な組の中から、良好な屈折率及び良好な膜厚を併せ持つ良好な組を再度選別する再選別ステップと、
   前記良好な膜厚の平均値を計算することによって精度良く前記薄膜の膜厚を決定する平均値計算ステップとを含むことを特徴とする請求項５に記載の膜の屈折率及び厚さの測定方法。
7. 前記修正ステップが、
   前記再選別ステップにおいて再度選別された前記良好な屈性率の平均値を計算し、該平均値から屈折率の探索範囲を決定する範囲決定ステップと、
   非線形関数の極小探査法によって、前記探索範囲の中で、前記良否判断ステップにおいて不良と判断された屈折率を新たに計算する再計算ステップとを含むことを特徴とする請求項６に記載の膜の屈折率及び厚さの測定方法。
8. 前記グループ数が、前記良好な組の数に応じて多くなるようにあらかじめ定めた規則に基づいて決定され、
   前記分類ステップが、各々のグループに属する膜厚の分散値の総和が最小になるように、前記良好な組を、前記グループ数のグループに分類するステップを含み、
   前記最適グループ決定ステップが、
   分類された前記グループの要素の数が２以上かつ最大である１つのグループ、および、そのグループを除く前記グループの要素の数が２以上かつ最大のグループを仮の最適グループとする仮最適グループ群決定ステップと、
   すべての仮最適グループＧ₁、Ｇ₂、Ｇ₃、…について、すべての良好な組を対象として
   、仮最適グループＧ_iの要素ｄ_ijの対数値lnｄ_ijの平均値(lnｄ)_ai（添字ａは平均値を表
   す）及び標準偏差σを用いて計算した規格化標準偏差（σ／(lnｄ)_ai）をσ_dNiとして、
   それぞれの仮最適グループＧ_iの(lnｄ)_aiに対する許容偏差δ_diがδ_di＝γσ_dNiによって表され、σ_dNiが大きくなるほど大きくなるようにあらかじめ定めたγの値に基づき、す
   べての良好な膜厚についての（lnｄ_k−(lnｄ)_ai）²の総和が最小になるように、前記良好な組の膜厚がabs（lnｄ_k−(lnｄ)_ai）＜δ_d（ここで、absは()内の絶対値を求める演算を表す）の条件を満たす良好な組（ｎ_k、ｄ_k）をそれぞれの仮最適グループに取り込む処理をＧ₁、Ｇ₂、Ｇ₃、…の順に行う取込ステップと、
   前記仮最適グループ群へのデータの取込処理によって変化した要素の組み合わせに応じて、仮最適グループ群Ｇ_iのうち要素数最大のグループおよび２番目のグループを決定し
   、それぞれ最適グループＧ₁および準最適グループＧ₂とする並べ替えステップと、
   前記最適および準最適グループＧ_iの要素数をＮ（Ｇ_i）、規格化標準偏差をσ_dN（Ｇ_i
   ）（ｉ＝１又は２）とし、「Ｎ（Ｇ₁）＝Ｎ（Ｇ₂）」および「σ_dN（Ｇ₂）＜σ_dN（Ｇ₁）」の条件を両方満たす場合、最適グループＧ₁と準最適グループＧ₂とを入れ替える処理を行い、最終的に最適グループＧ₁を決定する入替ステップとを含み、
   前記再選別ステップが、
   前記最適グループの屈折率ｎ_iの平均値ｎ_a及び標準偏差σを用いて規格化標準偏差σ_nN（＝σ／ｎ_a）を計算し、δ_n＝４σ_nNによって許容偏差δ_nを決定し、前記良好な組の屈
   折率の中で、abs（ｎ_k−ｎ_a）＜δ_nを満たす屈折率を前記良好な屈折率とし、
   前記最適グループの膜厚データｄ_iの対数値の平均値(lnｄ)_a及び標準偏差σを用いて規
   格化標準偏差σ_dN（＝σ／(lnｄ)_a）を計算し、許容偏差δ_dがδ_d＝γσ_dNによって表さ
   れ、σ_dNが大きくなるほど大きくなるようにあらかじめ定めたγの値に基づき、前記良好な組の膜厚がabs（lnｄ_k−(lnｄ)_ai）＜δ_dの条件を満たす膜厚を前記良好な膜厚とし、
   屈折率と膜厚が両方とも良好であることを条件として、良好なデータの組（ｎ_k、ｄ_k）を決定するステップであることを特徴とする請求項７に記載の膜の屈折率及び厚さの測定方法。
9. 前記範囲決定ステップが、前記良好な屈性率の平均値ｎ_gaと、あらかじめ定めた探索範囲の相対値ｒを用い、前記探索範囲の上限値を（１．０＋ｒ）ｎ_ga、下限値を（１．０−ｒ）ｎ_gaと決定するステップを含み、
   前記再計算ステップが、ｆ（ｎ）＝(Ｉ_S'−Ｉ_S)²−(Ｉ_C'−Ｉ_C)²（ここで、エリプソメメトリーの計数の測定値をΨ、Δとし、前記探索範囲内の屈折率ｎ'から計算されたエリ
   プソメメトリーの計数をΨ'、Δ'とし、Ｉ_S＝sin２ΨsinΔ、Ｉ_C＝sin２ΨcosΔ、Ｉ_S'＝sin２Ψ'sinΔ'、Ｉ_C'＝sin２Ψ'cosΔ'である）を評価関数とする極小探査法を用いて、前記不良と判断された屈折率を新たに計算するステップを含むことを特徴とする請求項８に記載の膜の屈折率及び厚さの測定方法。
10. コンピュータに、エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する機能を実現させるプログラムであって、
    前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、反射光を測定することによって得られた反射型エリプソメトリーの計数Ψ_r、Δ_rを、記録手段から読み出して、
    ln｜Ｘ｜＝０
    （ここで、Ｘは、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r−Ｂ₂ ^pＡ₂ ^s、ｂ＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^sＡ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−Ｂ₁ ^pＡ₁ ^sを係数とする２次方程式ａＸ²＋ｂＸ＋ｃ＝０
    の根であり、
    ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、
    係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記
    薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素として求められ、
    前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
    体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
    から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）、又は、
    (｜ｃ｜²−｜ａ｜²)²−｜ｃ^*ｂ−ｃｂ^*｜²＝０
    （ここで、ａ＝Ｂ₂ ^sＡ₂ ^pρ_r−Ｂ₂ ^pＡ₂ ^s、ｂ＝(Ｂ₁ ^sＡ₂ ^p＋Ｂ₂ ^sＡ₁ ^p）ρ_r−Ｂ₁ ^pＡ₂ ^s＋Ｂ₂ ^pＡ₁ ^s、ｃ＝Ｂ₁ ^sＡ₁ ^pρ_r−Ｂ₁ ^pＡ₁ ^s、ρ_r＝tanΨ_rexp(−ｊΔ_r)であり、
    係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｂ₁ ^s、Ｂ₂ ^s、Ｂ₁ ^p、Ｂ₂ ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜
    を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_m
    の行列要素として求められ、
    前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
    体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
    から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）
    を評価関数として、非線形方程式の数値解法を用いて、メモリ上で前記波長に対する屈折率を計算する機能を実現させる膜の屈折率及び厚さの測定プログラム。
11. コンピュータに、エリプソメトリーの計数を用いて膜の屈折率及び厚さを測定する機能を実現させるプログラムであって、
    前記膜が、２つの透明な媒体とそれらに挟まれた多層膜からなる光学系において、２つの前記媒体と前記多層膜との２つの界面のどちらか一方の界面に新しく形成される薄膜であって、前記一方の界面の側または他方の界面の側から少なくとも１つの波長の光を入射し、前記多層膜を透過した透過光を測定することによって得られた透過型エリプソメトリーの計数Ψ_t、Δ_tを、記録手段から読み出して、
    ln｜Ｘ｜＝０
    （ここで、Ｘは、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^sＡ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^sを係数とする、１次方程式ａＸ＋ｂ＝０の根であり、
    ρ_t＝tanΨ_texp(−ｊΔ_t)であり、
    係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｋ^s、Ｋ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多
    層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素
    として求められ、
    前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
    体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
    から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）、又は、
    (ａ^*＋ｂ^*)―(ａ＋ｂ)＝０
    （ここで、ａ＝Ｋ^pＡ₂ ^pρ_t−Ｋ^sＡ₂ ^s、ｂ＝Ｋ^pＡ₁ ^pρ_t−Ｋ^sＡ₁ ^s、ρ_t＝tanΨ_texp(−ｊ
    Δ_t)であり、
    係数Ａ₁ ^s、Ａ₂ ^s、Ａ₁ ^p、Ａ₂ ^p、Ｋ^s、Ｋ^pは、前記薄膜が新しく形成される前記一方の界面の側から光を入射する外入射型配置については、前記薄膜に入射する光が通過する前記媒体と前記薄膜との境界面の特性行列Ｉ₀₁、前記薄膜の特性行列Ｌ₁、前記薄膜を除いた多
    層膜と前記薄膜との間の境界面の特性行列Ｉ₁₂、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1から、Ｓ_m＝Ｉ₀₁Ｌ₁Ｉ₁₂Ｓ’_m-1によって計算される散乱行列Ｓ_mの行列要素
    として求められ、
    前記薄膜が新しく形成される側と反対側の前記他方の界面の側から光を入射する内入射型配置については、前記薄膜を除いた多層膜と前記薄膜との間の境界面の特性行列Ｉ_(m-1)m、前記薄膜の特性行列Ｌ_m、前記薄膜と前記薄膜を通過する光がさらに透過する前記媒
    体との境界面の特性行列Ｉ_m(m+1)、及び前記薄膜を除いた多層膜の擬似散乱行列Ｓ’_m-1
    から、Ｓ_m＝Ｓ’_m-1Ｉ_(m-1)mＬ_mＩ_m(m+1)によって計算される散乱行列Ｓ_mの行列要素として求められる）
    を評価関数として、非線形方程式の数値解法を用いて、メモリ上で前記波長に対する屈折率を計算する機能を実現させる膜の屈折率及び厚さの測定プログラム。
12. 請求項１０又は１１に記載の膜の屈折率及び厚さの測定プログラムを記録したコンピュータ読取可能な記録媒体。

Images