(57)【特許請求の範囲】
【請求項１】
試料膜に対して測定光を照射し、膜表面で反射する反射光と、膜中を透過して反対側の境界面で反射する透過反射光とによる干渉光のスペクトルを取得し、該測定スペクトルに基づいて試料膜の膜厚を求める膜厚測定方法において、測定スペクトルの波形形状を近似する近似スペクトルの生成情報として、膜厚を変数として含み波数又は波長に応じて周期変化する項及びオフセットの項を含む関数又は演算式を予め用意しておき、測定スペクトルに基づいて膜厚を算出する際に、
未知である膜厚を想定し、その想定膜厚において前記生成情報を用いて生成される近似スペクトルと測定スペクトルとの誤差に基づいて最適近似を行い、該想定膜厚における最小二乗誤差を求める誤差取得処理と、
前記想定膜厚を変化させる毎に前記誤差取得処理を実行することにより、前記最小二乗誤差を順次求める誤差変動情報取得処理と、
該誤差変動情報から極小点を与える１乃至複数の膜厚を見い出す膜厚取得処理と、
を含むことを特徴とする膜厚測定方法。
【請求項２】
前記周期変化項の大きさが波数又は波長に応じて変化することを特徴とする請求項１に記載の膜厚測定方法。
【請求項３】
前記関数又は演算式は更に、波数又は波長に応じて線形に変化する項を含むことを特徴とする請求項１又は２に記載の膜厚測定方法。
【請求項４】
前記近似スペクトルの生成情報は、膜厚を変数として含む、基底スペクトルの線形和による近似演算式であることを特徴とする請求項１〜３のいずれかに記載の膜厚測定方法。
【請求項５】
前記近似スペクトルの生成情報には複数の膜厚をそれぞれ変数として含み、前記誤差取得処理では、複数の想定膜厚の組み合わせにおける最小二乗誤差を求め、誤差変動情報取得処理では、複数の想定膜厚の組み合わせを変化させる毎に前記誤差取得処理を実行することを特徴とする請求項１〜４のいずれかに記載の膜厚測定方法。
【請求項６】
前記試料膜は同一膜が複数の膜厚を有して成るものであることを特徴とする請求項１〜５のいずれかに記載の膜厚測定方法。
【請求項７】
試料膜に対して測定光を照射し、膜表面で反射する反射光と、膜中を透過して反対側の境界面で反射する透過反射光とによる干渉光のスペクトルを取得し、該測定スペクトルに基づいて試料膜の膜厚を求める膜厚測定装置において、
a)測定スペクトルの波形形状を近似する近似スペクトルの生成情報として、膜厚を変数として含み波数又は波長に応じて周期変化する項及びオフセットの項を含む関数、演算式又はそれに相当する情報を予め格納しておく記憶手段と、
b)未知である膜厚を想定し、その想定膜厚において前記生成情報を用いて生成される近似スペクトルと測定スペクトルとの誤差に基づいて最適近似を行い、該想定膜厚における最小二乗誤差を求める誤差取得手段と、
c)前記想定膜厚を変化させる毎に前記誤差取得手段による処理を実行することにより、前記最小二乗誤差を順次求める誤差変動情報取得手段と、
d)該誤差変動情報から極小点を与える１乃至複数の膜厚を見い出す膜厚取得手段と、
を備えることを特徴とする膜厚測定装置。
【請求項８】
前記周期変化項の大きさが波数又は波長に応じて変化することを特徴とする請求項７に記載の膜厚測定装置。
【請求項９】
前記関数、演算式又はそれに相当する情報は更に、波数又は波長に応じて線形に変化する項を含むことを特徴とする請求項７又は８に記載の膜厚測定装置。
【請求項１０】
前記近似スペクトルの生成情報は、膜厚を変数として含む、基底スペクトルの線形和による近似演算式であることを特徴とする請求項７〜９のいずれかに記載の膜厚測定装置。
【請求項１１】
前記近似スペクトルの生成情報には複数の膜厚をそれぞれ変数として含み、前記誤差取得手段は、複数の想定膜厚の組み合わせにおける最小二乗誤差を求め、前記誤差変動情報取得手段は、複数の想定膜厚の組み合わせを変化させる毎に前記誤差取得手段による処理を実行することを特徴とする請求項７〜１０のいずれかに記載の膜厚測定装置。
【請求項１２】
前記試料膜は同一膜が複数の膜厚を有して成るものであることを特徴とする請求項７〜１１のいずれかに記載の膜厚測定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、分光測定を利用して試料膜の膜厚を測定するための膜厚測定方法及び膜厚測定装置に関する。本発明に係る膜厚測定方法及び膜厚測定装置は、例えば、半導体製造工程などにおいてウエハ基板上に形成された各種薄膜の膜厚の検査等の各種分野で広く利用することができる。
【０００２】
【従来の技術】
紫外光、可視光又は赤外光を利用した分光光度計の応用分野の一つとして、薄膜状の試料に対する膜厚の測定がある。分光測定を用いた膜厚測定の基本原理は次の通りである。
【０００３】
図１１に示すように、膜状の試料Ｓに対して単一波長の測定光を入射すると、その一部は試料Ｓの表面Ｓ１で反射され、残りは試料Ｓ内部に入り込んで、その一部は光の入射面とは反対側の境界面Ｓ２で反射して試料Ｓ内部を再び戻り、試料Ｓの表面Ｓ１から外部へと出射する。前者の反射光と後者の透過反射光とでは光路差が生じるため、測定光の波長と膜厚とに応じた干渉が発生する。測定光の波長を所定範囲で走査したとき、波数（又は波長）を横軸に、干渉光の強度を縦軸にとってグラフを描くと、波状の干渉スペクトル波形が得られる。この干渉スペクトル波形は余弦関数で表すことができ、その余弦関数の周期は膜厚に対応したものとなる。
【０００４】
そこで、干渉スペクトル波形を利用して、そのピークの山又は谷に対応する波数を自動又は手動で読み取り、それらの波数間隔情報を最小二乗法などにより求め、試料Ｓに対する既知の屈折率ｎを利用して波数周期から膜厚を算出する。
【０００５】
【発明が解決しようとする課題】
しかしながら、分光測定によって得られる干渉スペクトル波形は、種々の要因のために、理想的な余弦波形となることは殆どない。こうした干渉スペクトル波形を乱す要因としては、例えば、干渉効率の波数依存性、光源のエネルギ分布の波数依存性、装置の各種ノイズなどが考え得る。従来の膜厚の計算方法ではこうした要因が考慮されておらず、理想的な余弦波形となることを前提として膜厚を算出しているため、膜厚の算出精度を高めることが困難であった。
【０００６】
また、従来は膜厚が或る１種類であることを前提として計算を行っている。しかしながら、例えば、ＣＭＰ（Chemical Mechanical Polishing）装置などへの応用を考えると、意図的に複数の膜厚を持つように形成された試料膜に対して、その複数の膜厚を精度よく求め、適宜の時点で研磨を停止することが必要となる。上記のような従来の方法では、こうした要求には対応することが不可能であった。
【０００７】
本発明はこのような点に鑑みて成されたものであり、その第１の目的とするところは、従来よりも膜厚の算出精度を大幅に向上させることができる膜厚測定方法及び膜厚測定装置を提供することにある。
【０００８】
また、本発明の第２の目的とするところは、試料上の測定領域内に複数種類の膜厚を有するような試料膜に対しても、それぞれの膜厚を精度よく算出することができる膜厚測定方法及び膜厚測定装置を提供することにある。
【０００９】
【課題を解決するための手段】
上記課題を解決するために成された第１発明は、試料膜に対して測定光を照射し、膜表面で反射する反射光と、膜中を透過して反対側の境界面で反射する透過反射光とによる干渉光のスペクトルを取得し、該測定スペクトルに基づいて試料膜の膜厚を求める膜厚測定方法において、測定スペクトルの波形形状を近似する近似スペクトルの生成情報として、膜厚を変数として含み波数又は波長に応じて周期変化する項及びオフセットの項を含む関数又は演算式を予め用意しておき、測定スペクトルに基づいて膜厚を算出する際に、
未知である膜厚を想定し、その想定膜厚において前記生成情報を用いて生成される近似スペクトルと測定スペクトルとの誤差に基づいて最適近似を行い、該想定膜厚における最小二乗誤差を求める誤差取得処理と、
前記想定膜厚を変化させる毎に前記誤差取得処理を実行することにより、前記最小二乗誤差を順次求める誤差変動情報取得処理と、
該誤差変動情報から極小点を与える１乃至複数の膜厚を見い出す膜厚取得処理と、
を含むことを特徴としている。
【００１０】
また、第２発明は上記第１発明に係る膜厚測定方法を具現化する装置であって、試料膜に対して測定光を照射し、膜表面で反射する反射光と、膜中を透過して反対側の境界面で反射する透過反射光とによる干渉光のスペクトルを取得し、該測定スペクトルに基づいて試料膜の膜厚を求める膜厚測定装置において、
a)測定スペクトルの波形形状を近似する近似スペクトルの生成情報として、膜厚を変数として含み波数又は波長に応じて周期変化する項及びオフセットの項を含む関数、演算式又はそれに相当する情報を予め格納しておく記憶手段と、
b)未知である膜厚を想定し、その想定膜厚において前記生成情報を用いて生成される近似スペクトルと測定スペクトルとの誤差に基づいて最適近似を行い、該想定膜厚における最小二乗誤差を求める誤差取得手段と、
c)前記想定膜厚を変化させる毎に前記誤差取得手段による処理を実行することにより、前記最小二乗誤差を順次求める誤差変動情報取得手段と、
d)該誤差変動情報から極小点を与える１乃至複数の膜厚を見い出す膜厚取得手段と、
を備えることを特徴としている。
【００１１】
上記膜厚測定方法及び膜厚測定装置において、前記周期変化項の大きさは波数又は波長に応じて変化してもよい。また、前記関数、演算式又はそれに相当する情報は更に、波数又は波長に応じて線形に変化する項を含むものであってもよい。
更に、上記膜厚測定方法及び膜厚測定装置において、前記近似スペクトルの生成情報は、膜厚を変数として含む、基底スペクトルの線形和による近似演算式であるものとすることができる。
【００１２】
また、特に、複数種類の膜厚を有している試料膜の各膜厚を精度よく算出するために、第１発明に係る膜厚測定方法において、前記近似スペクトルの生成情報には複数の膜厚をそれぞれ変数として含み、前記誤差取得処理では、複数の想定膜厚の組み合わせにおける最小二乗誤差を求め、誤差変動情報取得処理では、複数の想定膜厚の組み合わせを変化させる毎に前記誤差取得処理を実行するものとすると好ましい。
【００１３】
同様に、第２発明に係る膜厚測定装置において、前記近似スペクトルの生成情報には複数の膜厚をそれぞれ変数として含み、前記誤差取得手段は、複数の想定膜厚の組み合わせにおける最小二乗誤差を求め、前記誤差変動情報取得手段は、複数の想定膜厚の組み合わせを変化させる毎に前記誤差取得手段による処理を実行する構成とすることが好ましい。
【００１４】
【発明の実施の形態】
まず、第１発明に係る膜厚測定方法において、分光測定の結果として得られた干渉スペクトル波形に基づいて膜厚を算出する膜厚算出方法の原理を説明する。
【００１５】
〔第１の膜厚算出方法〕
いま、或る試料膜を測定した結果得られる干渉スペクトルの一例を図３に示す。この干渉スペクトルの波形から、次のようなことが把握できる。
(ｉ)一定周期の干渉波パターンが存在する
(ii)オフセットが存在する
(iii)右上がりにほぼ線形のドリフトが存在する
(iv)干渉効率のため、波数が大きいほど干渉波の振幅が小さくなる
こうした点を考慮すると、干渉波の周期ωが既知であれば、干渉スペクトル波形は次の(1)式で近似できると予測される。
ｆ(x)＝α0＋α1ｘ＋α2（１/ｘ）sin(ωｘ＋δ)    …(1)
上記(1)式において、右辺の第１項はオフセットを反映し、第２項はドリフトを反映し、第３項は波形状の周期性波形を反映している。更に詳しく言えば、第３項の中の（１/ｘ）の部分は波数増加に伴う振幅の減少を反映し、δは特に膜厚が大きい場合に顕著になる位相ずれを反映している。
【００１６】
ここで、
sin(ωｘ＋δ)＝sinωｘcosδ＋cosωｘsinδ
であるから、(1)式は(2)式に変形することができる。
ｆ(x)＝α0＋α1ｘ＋α2（１/ｘ）sinωｘ＋α3（１/ｘ）cosωｘ    …(2)
すなわち、４つの関数ｆ0(x)，ｆ1(x)，ｆ2(x)，ｆ3(x)を、
ｆ0(x)＝１
ｆ1(x)＝ｘ
ｆ2(x)＝（１/ｘ）sinωｘ
ｆ3(x)＝（１/ｘ）cosωｘ
と定義すれば、測定によって取得された測定スペクトルｇ(x)は、上記４つの関数の線形和である、式(3)のように表される関数ｆ(x) で近似表現できる筈である。
ｆ(x)＝α0ｆ0(x)＋α1ｆ1(x)＋α2ｆ2(x)＋α3ｆ3(x)    …(3)
【００１７】
上記説明では周期ωが既知であるとしたが、実際にはこのω自体が未知である（正確に把握できない）。上記のような近似による近似スペクトルｆ(x)が測定スペクトルｇ(x)に最適に近似された状態というのは、両者の二乗誤差が最小であることと同意である。そこで、周期ωを変数とし、周期ωを変化させる毎に二乗誤差が最小となるような近似スペクトルを見つけてそのときの最小二乗誤差を取得し、周期ωを変化させたときの最小二乗誤差の変化において、最小二乗誤差が最も小さくなったときの周期ωが求めようとする周期である可能性が極めて高い。
【００１８】
さて、ここで、試料膜の膜厚と干渉波形の周期との関係を考えてみる。従来から知られているように、屈折率ｎである試料膜に入射角θで測定光が入射したとき（図１１参照）に得られる反射光に基づいて膜厚ｄを算出する際の基本式は、次の(4)式である。
ｄ〔cm〕＝Δｍ /｛２・（ｎ2−sin2θ）1/2（１/λ2−１/λ1)｝    …(4)
ここで、Δｍは、波長λ1のピークを基点として短波長側へピークを順次数えていったときに、波長λ2のピークが何番目に現われるかを示す数であり、通常はΔｍ＝１とすればよい。
【００１９】
周期Ｔ〔cm-1〕＝１/λ2−１/λ1
を用いると(4)式は、
ｄ＝１/２・（ｎ2−sin2θ）1/2Ｔ
であり、これを変形すると、
Ｔ＝１/｛２・（ｎ2−sin2θ）1/2ｄ｝＝1×107/｛２・（ｎ2−sin2θ）1/2ｄ'｝
となる。但し、ｄ’は単位がcmであるｄを単位nmで表現したものである。
いまθ＝０とすれば、sinθの項は無視できるから、次の(5)式の関係が成り立つ。
Ｔ＝1×107/（２・ｎ・ｄ'）      …(5)
一例として、屈折率ｎ：1.466、膜厚300nmのときの周期Ｔを計算すると、
Ｔ＝1×107/２×1.466×300＝11368cm-1
である。
【００２０】
(5)式で表したように、周期Ｔは膜厚ｄ'の関数である。したがって、上述したように未知の周期ωを変数とすることは、未知の膜厚を変数とすることに置き換えられる。一方、上述したように、周期が既知であれば、測定スペクトルｇ(x)は近似スペクトルｆ(x)で良好に近似することができる。そこで、近似スペクトルの表現として基底スペクトルを導入し、測定スペクトルｇ(x)と近似スペクトルｆ(x)との誤差を表す誤差関数を、膜厚ｄ'を変数として求めることを考える。
【００２１】
いま、測定スペクトルと近似スペクトルとの測定点をそれぞれ、
【数１】
0003928478.tif 000002
とする。このＸ、Ｙを用い、誤差関数ε(d')を次の(6)式で定義する。
【数２】
0003928478.tif 000003
【００２２】
上記(6-2）式において、Ｔ(ｄ')は上記(5)式によるｄ'を変数としたときの周期Ｔである。すなわち、膜厚ｄ'が決まれば(6-2）式で表現される基底スペクトルが定まるから、(6)式による二乗誤差が最小になるように測定点を見い出すことができる。したがって、膜厚として或る値を仮定すれば、測定スペクトル（原波形）に対する二乗誤差が最小となるような近似スペクトルを見つけることができ、その膜厚に対する最小二乗誤差が求まる。
【００２３】
このようにして所定の膜厚範囲で所定膜厚間隔毎に、それぞれの膜厚に対して基底を定義し、その基底による最小二乗誤差を求め、膜厚に対して最小二乗誤差を順次プロットしてゆくと図４に示すようになる。図４に描出したカーブによれば、膜厚が約1750nmのところに最小二乗誤差の極小点（ピークトップ）が存在する。すなわち、膜厚が1750nmであるときに近似スペクトルｆ(x)は測定スペクトルｇ(x)に最も類似した形状になるから、求めるべき膜厚は1750nmであると結論付けることができる。
【００２４】
上記のような膜厚算出方法が確かであることを検証するために、1750nmの膜厚を与える周期を式(5)より求めると、Ｔ＝1950cm-1である。このときの近似スペクトルを測定スペクトルに重ねて描くと、図５に示すようになる。測定スペクトルの原波形はほぼ2000cm-1毎にピークを持っており、オフセットやドリフト、波形の振幅の減衰状態まで含めて、かなり高い精度で近似できていることがわかる。したがって、上記方法により算出される膜厚の精度は高いことがわかる。
【００２５】
上記説明は試料膜の厚さが均一（単一膜厚）であることを前提としているが、同一膜であって膜厚が２つ以上存在している場合には、干渉スペクトル波形のピークは周期的には出現せず、その出現は後述の図９に示すように不規則となる。例えば、図６に示すように、測定対象の試料膜の膜厚がｄ1、ｄ2の２つ存在する場合、上記のような解析処理を行って膜厚と最小二乗誤差との関係を求めると、図７に示すようになる。すなわち、膜厚と最小二乗誤差との関係を表すカーブには、２個の極小点が現れる。この２つのピークが目視等によって明らかに分離可能である場合には、それぞれのピークトップの位置から２つの膜厚を求めることができる。
【００２６】
一方、２つの膜厚ｄ1、ｄ2がかなり近い場合には、図８に示すように２つピークの一部が重なり合って、目視上ではピークトップの位置が正確に把握できなくなることがある。そのような場合には、所定の演算処理を実行して第１、第２なる２つのピークを分離し、それぞれのピークトップの位置から膜厚を求めればよい。
【００２７】
〔第２の膜厚算出方法〕
上述したように、上記第１の膜厚算出方法でも複数の膜厚を同時に求めることは可能である。しかしながら、(3)式及び(6-2)式は、元々、単一膜厚であるという前提の下での近似表現であるため、複数膜厚に対しては精度の高い近似は難しいという限界がある。図９は、測定スペクトルと上記第１の膜厚算出方法において求まった最適な近似スペクトルとを重ねて示した図である。このように近似スペクトルは、測定スペクトルを良好に近似しているとは言い難い。そこで、存在する膜厚の数が既知である場合には、この第２の膜厚算出方法を用いることにより、一層高い精度で複数の膜厚をそれぞれ求めることができる。
【００２８】
この第２の膜厚算出方法でも、基本的な解析処理の原理は上記第１の膜厚算出方法と同様である。異なる点は、測定スペクトルｇ(x)を予め複数の膜厚を想定した基底スペクトルの線形和で近似表現すること、そして、複数の膜厚をそれぞれ或る値に仮定した条件で、二乗誤差が最も小さくなるような近似スペクトルｆ(x)を見つけることである。
【００２９】
具体的に説明すると、いま膜厚が２つ存在することが既知である場合、それぞれの膜厚ｄ1、ｄ2に対応する周期ω1、ω2を用いれば、干渉スペクトル波形は次式で近似することができると予測できる。
ｆ(x)＝α0＋α1ｘ＋α2（１/ｘ）sin(ω1ｘ＋δ1)＋α3（１/ｘ）sin(ω2ｘ＋δ2)
したがって、上記第１の膜厚算出方法と同様に、sin項を分解すれば、
ｆ(x)＝α0＋α1ｘ＋α2（１/ｘ）sinω1ｘ＋α3（１/ｘ）cosω1ｘ＋α4（１/ｘ）sinω2ｘ＋α5（１/ｘ）cosω2
とすることができる。
【００３０】
すなわち、６つの関数ｆ0(x)，ｆ1(x)，ｆ2(x)，ｆ3(x)，ｆ4(x)，ｆ5(x)を、
ｆ0(x)＝１
ｆ1(x)＝ｘ
ｆ2(x)＝（１/ｘ）sinω1ｘ
ｆ3(x)＝（１/ｘ）cosω1ｘ
ｆ4(x)＝（１/ｘ）sinω2ｘ
ｆ5(x)＝（１/ｘ）cosω2ｘ
と定義すれば、測定スペクトルｇ(x)は、上記６つの関数の線形和で(7)式のように表される関数ｆ(x) で近似表現できる筈である。
ｆ(x)＝α0ｆ0(x)＋α1ｆ1(x)＋α2ｆ2(x)＋α3ｆ3(x)＋α4ｆ4(x)＋α5ｆ5(x)    …(7)
【００３１】
この場合にも、上記(6)、(6-1)、(6-2)式と同様の形式の基底スペクトルの線形和の形式での表現が可能である。但し、(6)、(6-1)、(6-2)式では変数は１つの膜厚であったのに対し、ここで、膜厚ｄ1、ｄ2がともに変数となる。したがって、測定スペクトルｇ(x)と近似スペクトルｆ(x)との二乗誤差を算出する際には、２つの膜厚ｄ1、ｄ2をそれぞれ個別に仮定する必要がある。すなわち、一方の膜厚、例えばｄ1を固定し、他方の膜厚ｄ2を順次所定間隔毎に変化させていったときに、それぞれ最小二乗誤差を求める作業が必要である。
【００３２】
こうした処理の結果、膜厚の組（ｄ1，ｄ2）に対して最小二乗誤差が求まるから、膜厚と最小二乗誤差との関係は第１の膜厚算出方法のように２次元グラフでは表すことができず、２つの膜厚を２軸とし、最小二乗誤差を他の１軸とする３次元空間内の表現となる。図１０はこうした膜厚と最小二乗誤差との関係を示す図である。図１０において膜厚ｄ1、膜厚ｄ2の２軸と直交する方向（紙面に垂直な方向）に二乗誤差を示す曲面が存在する。このようにして、ｄ1＝ｄ2以外の部分で極小となる（ｄ1，ｄ2）を探し、それを求める２つの膜厚とすればよい。
【００３３】
もちろん、この方法を更に拡張して、３以上の膜厚を同時に算出することも可能である。また、元々、存在する膜厚の数が不明である場合には、第１の膜厚算出方法により膜厚の数を調べ、膜厚の数を特定した後に第２の膜厚算出方法を適用すれば、正確な膜厚算出が行える。
【００３４】
次に、上記のような膜厚算出方法を採用した膜厚測定装置の一実施例について説明する。図１は本実施例による膜厚測定装置の概略構成図である。
【００３５】
この膜厚測定装置は、分光測定部として、光源１、分光部２、測定光学系３、光検出器４を含み、光検出器４による検出信号はスペクトル作成部５に与えられ、該スペクトル作成部５により作成された干渉スペクトルが演算処理部６に与えられて、後述するような所定の演算処理が実行されることにより膜厚が算出される。なお、スペクトル作成部５及び演算処理部６の実体は、ＣＰＵを中心に構成されるパーソナルコンピュータであって、該コンピュータ上で所定のプログラムを実行することにより演算処理が達成される。
【００３６】
上記構成の動作を概略的に説明すると、まず光源１から発した白色光の中から、分光部２により特定の波長を有する単色光が取り出され、測定光学系３を介して膜状の試料Ｓに測定光として照射される。試料Ｓの表面や裏面などで反射した光は測定光学系３を介して光検出器４に導入され、これら反射光の強度に応じた電気信号がスペクトル作成部５に送られる。後述するように試料Ｓからの反射光は干渉光となるから、スペクトル作成部５は、測定光の波長走査に対応して光検出器４で得られる信号に基づいて、横軸が波数（又は波長）、縦軸が相対強度である干渉スペクトルを作成する。これが上記測定スペクトルに相当する。
【００３７】
演算処理部６はこの干渉スペクトルを受け取り、上記のような第１又は第２の膜厚算出方法に基づく演算処理を実行する。ここでは第１の膜厚算出方法に基づく演算処理について説明するが、第２の膜厚算出方法に基づく演算処理も同様に具現化できることは明らかである。図２は、演算処理の具体的な処理手順を示すフローチャートである。
【００３８】
まず、演算条件として、外部より膜厚範囲ｄa〜ｄbと膜厚間隔Δｄを指示する（ステップＳ１）。予めおおよその膜厚が既知である場合には、その膜厚を含むように適宜の膜厚範囲ｄa〜ｄbを設定することにより、演算時間を短縮することが可能である。また、算出精度の点からは膜厚間隔Δｄを小さくすることが好ましいが、膜厚間隔Δｄを小さくするほど演算時間が長くなるから、必要な精度との兼ね合いで適宜に決めるとよい。
【００３９】
実際に演算処理が開始されると、まず計算上の想定膜厚ｄをｄaに設定し（ステップＳ２）、上述したような最小二乗法により、測定スペクトルｇ(x)に対する二乗誤差εが最小となるような近似スペクトルｆ(x)を探索する。そして、その近似スペクトルが与える最小二乗誤差を取得する（ステップＳ３）。
【００４０】
次に、そのときの想定膜厚ｄがｄb以上であるか否かを判定し（ステップＳ４）、ｄb未満であれば、ｄ＋Δｄを新たな想定膜厚ｄとし（ステップＳ５）、ステップＳ３へと戻る。そして、ステップＳ４で想定膜厚ｄがｄb以上であると判定されるまで、ステップＳ３、Ｓ４、Ｓ５なる処理を繰り返す。これにより、膜厚範囲ｄa〜ｄb内でΔｄ間隔毎に二乗誤差が最小になるような近似スペクトルが探索され、それに対応する最小二乗誤差が順次取得される。
【００４１】
その結果、図４に示したような、膜厚と最小二乗誤差との関係を示すカーブが得られる（ステップＳ６）。そして、そのカーブにおいて極小点を探し、極小点を与える膜厚を取得する（ステップＳ７）。このようにして、高い精度で膜厚を算出することができる。
【００４２】
膜厚が複数である場合でも、基本的に同様の手順で測定スペクトルから膜厚を算出できることは、既述の説明から明らかである。
【００４３】
なお、上記実施例は本発明の一例にすぎず、本発明の趣旨の範囲で適宜変更や修正を行えることは明らかである。
【００４４】
例えば、上記説明において近似スペクトルを表現するための基底のとり方は任意であって、上記記載のものに限らない。また、特に膜厚が薄い場合には、干渉スペクトル波形において明確な周期を持つピークが得られにくい傾向にある。このような場合には、基底としてsin項や１／ｘ項を除外するほうが良好な近似が行い易い。また、それ以外にも、(3)式又は(7)式に示した関数ｆ0(x)，ｆ1(x)，ｆ2(x)，ｆ3(x)，ｆ4(x)，ｆ5(x)の各項を適宜に変形又は加減することにより、近似精度を一層向上させることができる。
【００４５】
【発明の効果】
以上説明したように、第１発明に係る膜厚測定方法及び第２発明に係る膜厚測定装置によれば、様々な要因によって測定スペクトルに現れるオフセット、ドリフト、位相ずれ、周期性波形の振幅変化などの各種要素が考慮されるため、膜厚の算出精度が大幅に向上する。また、従来は、その測定領域内に複数種類の膜厚を有するような試料に対してその複数膜厚を同時に測定することは実質的に不可能であったが、このような測定も可能となり、しかも高い精度でもって複数の膜厚を算出することができる。
【図面の簡単な説明】
【図１】  本発明の一実施例である膜厚測定装置の概略構成図。
【図２】  本実施例の膜厚測定装置における膜厚算出の具体的な処理手順を示すフローチャート。
【図３】  或る試料膜を分光測定した結果得られる干渉スペクトル波形の一例を示す図。
【図４】  本発明の膜厚測定方法により得られる、膜厚とスペクトル近似の最小二乗誤差との関係の一例を示すグラフ。
【図５】  測定スペクトルと本発明の膜厚測定方法で求められる近似スペクトルとを重ねて示した図。
【図６】  ２つの膜厚を有する試料膜の断面構造の一例を示す図。
【図７】  試料膜が２つの膜厚を有する場合の膜厚とスペクトル近似の最小二乗誤差との関係の一例を示す図。
【図８】  ２つの膜厚が近い場合のピーク波形を示す図。
【図９】  試料膜が２つの膜厚を有する場合の測定スペクトルと近似スペクトルとの関係を示す図。
【図１０】  試料膜が２つの膜厚を有する場合の、膜厚と最小二乗誤差との関係を示す図。
【図１１】  分光測定を利用した膜厚測定の原理を説明するための図。
【符号の説明】
１…光源
２…分光部
３…測定光学系
４…光検出器
５…スペクトル作成部
６…演算処理部
Ｓ…試料
【図１】
0003928478.tif 000004
【図２】
0003928478.tif 000005
【図３】
0003928478.tif 000006
【図４】
0003928478.tif 000007
【図５】
0003928478.tif 000008
【図６】
0003928478.tif 000009
【図７】
0003928478.tif 000010
【図８】
0003928478.tif 000011
【図９】
0003928478.tif 000012
【図１０】
0003928478.tif 000013
【図１１】
0003928478.tif 000014