(57)【特許請求の範囲】
【請求項１】
処理ガスが導入された真空容器内に電界を形成して，未知の変動周期で変動するプラズマを生成し，前記真空容器内に配された被処理体に対してプラズマ処理を施すプラズマ処理方法において：
前記プラズマのプラズマ光を一定のサンプリング周期でサンプリングし，サンプリング信号のデータ列を求める工程と；
前記サンプリング信号のデータ列から複数の移動平均値算出期間における移動平均値を算出し，前記各移動平均値に対応する複数のデータ列を求める工程と；
前記各移動平均値のデータ列からそれぞれに対応する複数の近似式を求める工程と；
所定期間における前記各移動平均値のデータ列と前記各移動平均値のデータ列に対応する前記各近似式との変位量をそれぞれ求める工程と；
前記各変位量の中から最小の変位量を求める工程と；
前記最小の変位量に対応する前記移動平均値算出期間を前記プラズマの変動周期と判定する工程とを含み，
前記移動平均値算出期間は，前記サンプリング周期の整数倍に限定されず，前記サンプリング周期の整数倍ではない時刻における前記移動平均値の算出対象となるデータは，前記時刻の直前と直後の前記サンプリング信号データから算出される疑似サンプリングデータであることを特徴とする，プラズマ処理方法。
【請求項２】
前記移動平均値算出期間の最大値は，予想される前記変動周期の上限値であり，前記移動平均値算出期間の最小値は，予想される前記変動周期の下限値であることを特徴とする，請求項１に記載のプラズマ処理方法。
【請求項３】
前記所定の期間は，予想される前記変動周期の上限値以上であることを特徴とする，請求項１または２に記載のプラズマ処理方法。
【請求項４】
処理ガスが導入された真空容器内に電界を形成して，既知の変動周期で変動するプラズマを生成し，前記真空容器内に配された被処理体に対してプラズマ処理を施すプラズマ処理方法において：
前記プラズマのプラズマ光を一定のサンプリング周期でサンプリングし，サンプリング信号のデータ列を求める工程と；
前記サンプリング信号のデータ列から，移動平均値算出期間の移動平均値を算出し，移動平均値のデータ列を求める工程と；
前記移動平均値のデータ列からプラズマ処理の終点を判定する工程とを含み；
さらに，前記移動平均値算出期間は，前記サンプリング周期の整数倍に限定されないことを特徴とする，プラズマ処理方法。
【請求項５】
前記サンプリング周期の整数倍ではない時刻における前記移動平均値の算出対象となるデータは，前記時刻の直前と直後の前記サンプリングデータから算出される疑似サンプリングデータであることを特徴とする，請求項４に記載のプラズマ処理方法。
【請求項６】
前記移動平均値を算出するに際し，前記サンプリングデータは，前記疑似サンプリングデータに対して重み付けされていることを特徴とする，請求項５に記載のプラズマ処理方法。
【請求項７】
前記真空容器内にはさらに回転磁界が形成されることを特徴とする，請求項１〜６のいずれかに記載のプラズマ処理方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は，プラズマ処理方法に関する。
【０００２】
【従来の技術】
従来，半導体製造装置の分野においては，各種プラズマ源を用いたプラズマ処理装置の構成が提案されている。その内の一つに，処理室内に対向配置された上部電極と下部電極との間に電界を形成し，その電界により処理室内に導入された処理ガスをプラズマ化して，被処理体に対して所定のプラズマ処理を施すことが可能なプラズマ処理装置がある。さらに，上記プラズマ処理装置の中には，処理室内に回転磁界を形成可能な磁石を備えた，いわゆる磁場アシスト型のプラズマ処理装置がある。かかる構成を採用することにより，処理室内に形成された磁界によりプラズマ中の電子をトラップし，処理ガス粒子との衝突回数を増加させて，高密度プラズマを励起することが可能となり，さらに磁界を回転させることにより，プラズマ密度を均一化させ，高速かつ均一なプラズマ処理を実現することが可能である。
【０００３】
ところで，エッチングなどのプラズマ処理工程においては，プラズマ処理の終了時点を正確に判定して，プラズマ処理を遅滞なく終了することが重要である。そして，従来，処理室内のプラズマ中に含まれる特定物質のスペクトラム光の変化を検出し，その変化に基づいて終点検出を行う方法が提案されている。かかる方法は，被処理体に対するエッチングの進行とともにプラズマ中に含まれる成分も変化することに着目し，ある特定物質のスペクトラム光の強度変化を観測することにより，エッチング処理の終点を正確に実時間で検出しようとするものである。
【０００４】
しかし，上述のような磁場アシスト型プラズマ処理装置では，処理室内に形成される磁界方向に応じて，プラズマ中に密度の異なる領域が形成されるが，磁界の回転に応じて，プラズマの密度分布も変動してしまう。そのため，例えば，処理室の一壁に設けられた検出窓などから，プラズマ光を定点観測する場合には，回転磁界に従うプラズマの変動も考慮に入れる必要がある。
【０００５】
そのため，例えば，特開平４−３３８６６３号公報においては，磁石の回転と同期してパルスを発生するロータリエンコーダをエッチング装置に設けて，そのパルスに応じてプラズマ光をサンプリングして，磁石の回転周期に伴うノイズ成分を除去し，正確な終点検出を行おうとする技法が開示されている。
【０００６】
【発明が解決しようとする課題】
しかしながら，上記のようにハードウェア的に磁石の回転周期を求め，その回転周期に応じてサンプリングを行う構成では，ロータリエンコーダなどの装置を処理装置に付加せねば成らず，装置構成が複雑となるとともに，装置のイニシャルコストを押し上げるという問題があった。
【０００７】
また，上記のようにハードウェア的に磁石の回転周期を求め，その回転周期に応じてサンプリング行う構成では，サンプリング周期はその都度変動する変動周期である。しかし，終点検出用の信号処理ソフトウェアによって，サンプリング周期として固定周期を使用した方が効率的な処理を行える場合や，さらには固定周期を要求するようなものある。しかし，上記構成では，そのようなソフトウェア側の要求に柔軟に対応することができないという問題があった。
【０００８】
【課題を解決するための手段】
上記課題を解決するため，本発明の第１の観点によれば，処理ガスが導入された真空容器内に電界を形成して，未知の変動周期で変動するプラズマを生成し，真空容器内に配された被処理体に対してプラズマ処理を施すプラズマ処理方法において，プラズマのプラズマ光を一定のサンプリング周期でサンプリングし，サンプリング信号のデータ列を求める工程と，サンプリング信号のデータ列から複数の移動平均値算出期間における移動平均値を算出し，各移動平均値に対応する複数のデータ列を求める工程と，各移動平均値のデータ列からそれぞれに対応する複数の近似式を求める工程と，所定期間における各移動平均値のデータ列と各移動平均値のデータ列に対応する各近似式との変位量をそれぞれ求める工程と，各変位量の中から最小の変位量を求める工程と，最小の変位量に対応する移動平均値算出期間をプラズマの変動周期と判定する工程とを含むことを特徴とするプラズマ処理方法が提供される。
【０００９】
さらに本発明の別の目的は，相対的に間隔の広いサンプリング周期によるサンプリング信号から疑似サンプリング信号を作ることにより，実際にサンプリングする実サンプリング回数を増やすことなく，相対的に少ない演算量で，あたかも実サンプリング回数を増やしたかのようなより正確な移動平均値を求めることが可能な，新規かつ改良されたプラズマ処理方法を提供することである。
【００１０】
【課題を解決するための手段】
上記課題を解決するため，本発明の第１の観点によれば，請求項１に記載の発明のように，処理ガスが導入された真空容器内に電界を形成して，未知の変動周期で変動するプラズマを生成し，真空容器内に配された被処理体に対してプラズマ処理を施すプラズマ処理方法において，プラズマのプラズマ光を一定のサンプリング周期でサンプリングし，サンプリング信号のデータ列を求める工程と，サンプリング信号のデータ列から複数の移動平均値算出期間における移動平均値を算出し，各移動平均値に対応する複数のデータ列を求める工程と，各移動平均値のデータ列からそれぞれに対応する複数の近似式を求める工程と，所定期間における各移動平均値のデータ列と各移動平均値のデータ列に対応する各近似式との変位量をそれぞれ求める工程と，各変位量の中から最小の変位量を求める工程と，最小の変位量に対応する移動平均値算出期間をプラズマの変動周期と判定する工程とを含むことを特徴とするプラズマ処理方法が提供される。
【００１１】
かかる構成によれば，一定のサンプリング周期でプラズマ光のサンプリングをした場合であっても，各処理ごとに変化する可能性があるプラズマの変動周期を，ソフトウェア的な統計演算処理のみで算出することができる。従って，特段のハードウェア装置を付加せずとも，各処理ごとにプラズマの変動周期を正確に求めることができる。
【００１２】
また，例えば上記移動平均値算出期間の最大値を，予想される変動周期の上限値とし，移動平均値算出期間の最小値を，予想される変動周期の下限値とすることにより，予想される変動周期の全範囲からプラズマの変動周期を求めることができ，その変動周期を正確かつ確実に判定することができる。
【００１３】
さらに，上記移動平均値算出期間は，任意に設定することが可能であり，例えば一定のサンプリング周期の整数倍に限定されない。すなわち，本発明によれば，プラズマの変動周期がサンプリング周期の整数倍でない場合であっても，ソフトウェア的にプラズマの変動周期を判定することができる。すなわち，プラズマの変動周期がサンプリング周期の整数倍でない場合であっても，サンプリング周期の整数倍ではない時刻における移動平均値の算出対象となるデータは，時刻の直前と直後のサンプリング信号のデータから疑似サンプリングデータとして算出することができるので，その疑似サンプリングデータに基づいて，プラズマの変動周期をソフトウェア的に求めることができる。
【００１４】
さらに，所定期間にわたる各移動平均値のデータ列と各移動平均値のデータ列に対応する各近似式との変位量をそれぞれ求めるにあたっては，例えばその所定期間を予想される変動周期の上限値以上とすれば，プラズマの変動周期を判定するために必要な複数の変位量を全て求めることができる。
【００１５】
また，上記課題を解決するため，本発明の第２の観点によれば，処理ガスが導入された真空容器内に電界を形成して，既知の変動周期で変動するプラズマを生成し，真空容器内に配された被処理体に対してプラズマ処理を施すプラズマ処理方法において，プラズマのプラズマ光を一定のサンプリング周期でサンプリングし，サンプリング信号のデータ列を求める工程と，サンプリング信号のデータ列から移動平均値算出期間の移動平均値を算出し，移動平均値のデータ列を求める工程と，移動平均値のデータ列からプラズマ処理の終点を判定する工程とを含み，さらに移動平均値算出期間は，サンプリング周期の整数倍に限定されないことを特徴とするプラズマ処理方法が提供される。
【００１６】
かかる構成によれば，一定のサンプリング周期でプラズマ光のサンプリングをした場合であっても，各処理ごとに変化する可能性があるプラズマの変動周期を算出し，その変動分を補償した上で，プラズマ処理の終点をすることができる。従って，特段のハードウェア装置を付加せずとも，ソフトウェア的な処理のみで，各処理ごとに変動するプラズマ処理の終点を正確に求めることができる。
【００１７】
また，サンプリング周期の整数倍ではない時刻における移動平均値の算出対象となるデータは，例えば時刻の直前と直後のサンプリングデータから疑似サンプリングデータとして算出し，その疑似サンプリングデータから移動平均値を算出することができる。
【００１８】
さらに，移動平均値を算出するに際し，例えばサンプリングデータを疑似サンプリングデータに対して重み付けするように構成してもよい。移動平均値を算出する際に，プロセスに応じた変動値を適宜補償することが可能となり，終点検出の精度を高めることができる。
【００１９】
なお，真空容器内におけるプラズマの変動は，回転磁界によって生じることがあるが，このような回転磁界が存在する場合に，特に本発明を好適に適用することができる。
【００２０】
【発明の実施の形態】
以下に，添付図面を参照しながら，本発明にかかるプラズマ処理方法をエッチング処理の終点判定方法に適用した実施の一形態について詳細に説明する。
【００２１】
まず，かかる終点判定方法を適用可能なエッチング装置１００の装置構成について，図１を参照しながら説明する。
【００２２】
エッチング装置１００の処理室１０２は，導電性の真空容器１０４内に形成されている。処理室１０２内には，サセプタを構成する導電性の下部電極１０６と，接地された導電性の上部電極１０８が対向配置されている。下部電極１０６上には，処理時に，被処理体，例えば半導体ウェハ（以下，「ウェハ」と称する。）Ｗを載置固定することができる。また，下部電極１０６には，整合器１１６を介して高周波電源１１８が接続されており，処理時には，演算制御器１２０からの制御信号に応じて所定出力の高周波電力を印加可能である。
【００２３】
また，処理室１０２には，処理室１０２内に所定の処理ガスを供給するガス供給管１２４と，処理室１０２内の雰囲気を排気する排気管１２６が接続されている。さらに，真空容器１０４の上方には，処理室１０２内に所定の回転磁界を形成可能な永久磁石などの磁石１２８が配置されている。
【００２４】
また，処理室１０２の側壁には検出窓１３４が形成されており，この検出窓１３４を透過したプラズマ光が受光部１３６により検出される。この受光部１３６は，プラズマ光から特定物質に関する発光スペクトルのみを分光する不図示の分光器や，発光スペクトルの発光強度を検出する不図示の光電変換器や，検出されたアナログ信号をディジタル信号に変換するＡ／Ｄ変換器などを備えており，サンプリングされたデータは適宜演算制御器１２０に出力される。なお，図示の例では，受光部１３６は直接プラズマ光を検出構成となっているが，検出窓１３４と受光部１３６の間に，例えば光ファイバなどの光伝達手段を介装した構成としても良い。
【００２５】
次に，エッチング装置１００において，ウェハＷに対してエッチング処理を施す場合について説明する。まず，ウェハＷを下部電極１０６に載置して，ガス供給管１２４から処理室１０２内に処理ガスを導入するとともに，排気管１２６から処理室１０２内の雰囲気を排気して，処理室１０２内を所定の減圧雰囲気に維持する。次いで，真空容器１０４上方の磁石１２８を所定の回転数で回転させ，処理室１０２内に回転磁界を形成する。
【００２６】
しかる後，高周波電源１１８から下部電極１０６に高周波電力を印加すると，処理室１０２内に形成された電界により処理ガスが解離して高密度プラズマＰが生成され，所定のエッチング処理がウェハＷに施される。さらに，処理室１０２内に形成された回転磁界によって，プラズマＰは所定の回転周期で変動し，密度むらが平均化されて，ウェハＷに均一な処理を施すことができる。
【００２７】
次に，上述したエッチング装置１００に適用される本実施の形態にかかるエッチング処理の終点判定方法について，図２〜図５を参照しながら説明する。かかる終点判定方法は，プラズマの変動（回転）周期判定工程とエッチング処理の終点判定工程から構成されている。変動周期判定工程は，磁石１２８の回転に応じて変化する磁界の変動周期，すなわちプラズマの変動周期を判定する工程である。また，終点判定工程は，その判定された磁石１２８の回転周期に基づいて，プラズマの変動に伴って生じたサンプリング信号中のノイズ成分を除去した後，正確な終点判定を行う工程である。以下に，変動周期判定工程（Ａ）と終点判定工程（Ｂ）に分けて説明する。
【００２８】
（Ａ）変動周期判定工程
まず，図２〜図４を参照しながら回転周期判定工程について説明する。図２には，上記エッチング装置１００の受光部１３６においてサンプリングされたサンプリング信号Ｐ（ｔ）の時系列的変化が示されている。図示のように，プラズマ立ち上がり時には，サンプリング信号も不安定に推移するが，プラズマの安定とともにサンプリング信号も安定する。そこで，サンプリング信号が安定した後（図示の例では，時刻ａ以降）の期間Ａにおけるサンプリング信号に基づいて，プラズマの変動周期が求められる。
【００２９】
図３は，図２の期間Ａにわたるサンプリング信号の推移を拡大して示すものであるが，図示の例では０．１秒間隔でサンプリングが行われているものとする。ここで，処理室１０２内の磁界の変動周期は，ほぼ磁石１２８の回転周期に連動しているので，その大凡の周期についてはグローバルに推測することが可能である。例えば，本実施の形態にかかる方法では，プラズマは磁石１２８の回転周期に応じて３．０〜３．１秒周期で変動しているものとする。
【００３０】
したがって，３．０〜３．１秒周期の間において，プラズマの正確な変動周期を例えば０．０１秒のオーダーで求めようとする場合には，本来であれば，０．０１秒間でサンプリングを行う必要がある。しかし，単純にサンプリング周期を短くしたのでは，サンプリングデータの量がいたずらに増加し，データ処理に時間を要するという問題があった。この点，本実施の形態にかかる方法によれば，０．１秒間隔のサンプリングによっても，０．０１秒間隔，あるいはそれ以上に短い間隔でサンプリングを行った場合と同等の精度で，プラズマの変動周期を求めることが可能となる。以下，その方法について詳細に説明する。
【００３１】
まず，プラズマが３．０秒周期で変動しているものと仮定する。そして，従来の方法と同様に，サンプリングデータ中に含まれているノイズ成分を除去するために，サンプリングデータの移動平均値Ｈ３．００（ｔ）を下記の式（１）に基づいて求める。この場合，移動平均算出期間は，３．０秒である。
0003563949.tif 000002
【００３２】
次いで，プラズマが３．１秒周期で変動しているものと仮定する。そして，先の処理と同様に，サンプリングデータ中に含まれているノイズ成分を除去するために，サンプリングデータの移動平均値Ｈ３．１０（ｔ）を下記の式（２）に基づいて求める。この場合，移動平均算出期間は，３．１秒である。
0003563949.tif 000003
【００３３】
ここで，上記の場合には，０．１秒のサンプリング周期に対して，移動平均算出期間は，サンプリング周期の整数倍である３．０秒及び３．１秒にそれぞれ設定している。しかしながら，実際のプラズマの変動周期はサンプリング周期の整数倍とは限らない。例えば，実際のプラズマの変動周期が３．０３秒であるとすると，その移動平均算出期間は３．０３秒であり，０．１秒のサンプリング周期では測定不能であり，移動平均を求めるのであれば，サンプリング周期を０．０１秒で行う必要がある。しかしながら，このように，サンプリング周期を細かく設定したのでは，いたずらにデータ量が増えてしまい，処理の効率化を図ることができない。
【００３４】
この点，本実施の形態においては，以下に述べるような手順により，０．１秒のサンプリング周期でサンプリングしたデータを，あたかも０．０１秒のサンプリング周期でサンプリングしたデータであるかのように疑似し，実際のプラズマの周期変動に即した移動平均を可能にするものである。
【００３５】
例えば，プラズマが３．０３秒の周期で変動しているものと想定した場合の時刻ｔにおける移動平均値Ｈ３．０３（ｔ）は，
0003563949.tif 000004
で表される式（３）によって求めることができる。なお，式（３）中のＳ（ｔ−２．９−０．０１）や，Ｓ（ｔ−２．９−０．０２）や，Ｓ（ｔ−２．９−０．０３）で表されるデータは，実際にサンプリングされたデータではなく，実際にサンプリングされたデータからソフトウェア的に疑似されたデータであるので，本明細書中では疑似サンプリングデータと称することにする。これらの疑似サンプリングデータは，隣接するサンプリング信号のデータ間を直線又は近似曲線で結び，その直線又は近似曲線上における時刻ｔ−２．９−０．０１や，時刻ｔ−２．９−０．０２や，時刻ｔ−２．９−０．０３などの疑似サンプリングタイミングにおける対応値として求めることができる。このように，本明細書においては，疑似サンプリングデータを求めるために想定された時間間隔を疑似サンプリングタイミングと称している。
【００３６】
なお，隣接するサンプリング信号のデータ間を結ぶ近似曲線を求めるにあたっては，最小二乗法により１次以上の多項式で近似する方法や，ラグランジュ補間法や，スプライン補間法などの各種統計学手法を用いて，複数点のデータから求めることができる。
【００３７】
また，例えば隣接するサンプリング信号のデータ間を直線で結んだ場合には，上記Ｓ（ｔ−２．９−０．０１）は，
Ｓ（ｔ−２．９−０．０１）＝Ｓ（ｔ−２．９）＋（Ｓ（ｔ−３．０）−Ｓ（ｔ−２．９））／１０＊１    …（４）
で表される式（４）によって求めることができ，またＳ（ｔ−２．９−０．０２）は，
Ｓ（ｔ−２．９−０．０２）＝Ｓ（ｔ−２．９）＋（Ｓ（ｔ−３．０）−Ｓ（ｔ−２．９））／１０＊２    …（５）
で表される式（５）によって求めることができ，さらにＳ（ｔ−２．９−０．０３）は，
Ｓ（ｔ−２．９−０．０３）＝Ｓ（ｔ−２．９）＋（Ｓ（ｔ−３．０）−Ｓ（ｔ−２．９））／１０＊３    …（６）
で表される式（６）によって求めることができる。
【００３８】
このように，本実施の形態にかかる回転周期判定工程では，移動平均値Ｈ（ｔ）や，疑似サンプリングデータであるＳ（ｔ−２．９−０．０１）などを上述した式（１）〜（３）などや，式（４）〜式（６）など用いて算出することができる。かかる構成により，サンプリング周期の整数倍とはならない期間の移動平均値であっても，より少ないサンプリングデータから求めることができる。例えば上述の如く磁石１２８の回転周期を最大の回転周期である３．１秒と想定した場合でも，１１個のサンプリングデータであるＳ（ｔ−３．０），……，Ｓ（ｔ）から求めることができる。そして，本実施の形態によれば，結果的に，加重平均の手法により，疑似サンプリングデータを含めて３０３個のデータの移動平均値を少ない演算処理で算出することができる。
【００３９】
同様の方法により，磁石１２８の回転周期を３．０１秒〜３．０９秒と想定した場合の時刻ｔにおける移動平均Ｈ（ｔ）を求めることができる。
【００４０】
すなわち，図２中の期間Ａおいて，プラズマ光の信号を０．１秒間隔でサンプリングする度に，３．００秒，３．０１秒，３．０２秒，３．０３秒，３．０４秒，３．０５秒，３．０６秒，３．０７秒，３．０８秒，３．０９秒，３．１０秒の各移動平均値算出期間に対する移動平均値を並列的に，すなわち同時に算出する。このようにして，各移動平均値，Ｈ３．００（ｔ），Ｈ３．０１（ｔ），Ｈ３．０２（ｔ），Ｈ３．０３（ｔ），Ｈ３．０４（ｔ），Ｈ３．０５（ｔ），Ｈ３．０６（ｔ），Ｈ３．０７（ｔ），Ｈ３．０８（ｔ），Ｈ３．０９（ｔ），Ｈ３．１０（ｔ）を求めることができる。さらに，複数回のサンプリングを実行することにより，各サンプリング時刻に対応する各移動平均値Ｈ３．００（ｔ）〜Ｈ３．１０（ｔ）のデータ列を求めることができる。図３に，これら各移動平均値のうち，Ｈ３．０３（ｔ）をある期間にわたってプロットしたものをＨ（ｔ）として例示した。
【００４１】
以上のように，実サンプリング周期の１０倍の精度で移動平均値を求める方法について説明したが，本実施の形態はかかる構成に限定されるものではなく，例えば，実サンプリング周期の１００倍の精度で移動平均値を求めることができる。すなわち，例えばプラズマの回転周期を３．０３１秒を想定した場合の時刻ｔにおける移動平均値Ｈ３．０３１（ｔ）は，
0003563949.tif 000005
で表される式（７）によって求めることができる。
【００４２】
また，式（７）中のＳ（ｔ−２．９−０．０１）と，Ｓ（ｔ−２．９−０．０２）と，Ｓ（ｔ−２．９−０．０３）は，それぞれに対応する上記式（４）と，式（５）と，式（６）により求めることができる。さらに，同式中のＳ（ｔ−２．９−０．０３−０．００１）は，
Ｓ（ｔ−２．９−０．０３−０．００１）＝Ｓ（ｔ−２．９）＋（Ｓ（ｔ−３．０）−Ｓ（ｔ−２．９））／１００＊３１    …（８）
で表される式（８）により求めることができる。
【００４３】
次に，図４を参照しながら，上記工程により必要な数だけ求めた各移動平均値のデータ列から近似式，例えば一次近似式を求める工程について説明する。図４は，図３に示したＨ（ｔ）を拡大して示したものである。また，図４に示したＫ（ｔ）は，所定期間のＨ（ｔ）の値から，例えば最小二乗法で求めた一次近似式である。なお，回転周期を所定の精度で求めるためには，上記所定の期間は，少なくとも想定される最大の回転周期以上であることが好ましい。ここで，本実施の形態の場合には，プラズマの変動周期，すなわち，磁石１２８の回転周期が，予め，３．０秒〜３．１秒の範囲内であることがわかっているため，その所定の期間を３．１秒とする。そして，上述したＨ３．００（ｔ）〜Ｈ３．１０（ｔ）の各データ列に対応する一次近似式を求め，それら各一次近似式をそれぞれに対応してＫ３．００（ｔ）〜Ｋ３．１０（ｔ）とする。
【００４４】
次に，上記各工程により求められた各移動平均値のデータ列と，それら各移動平均値のデータ列に対応する各一次近似式の偏差値の絶対値の平均を求める工程について説明する。かかる偏差値の絶対値の平均Ａは，
0003563949.tif 000006
で表される式（９）によって求めることができる。なお，式（９）中の時刻ｔは，上述した移動平均値を求める際に用いた式（１）〜式（８）中の時刻ｔとは異なる時刻を表している。
【００４５】
上述した式（９）は，移動平均値を求める期間，すなわち想定される磁石１２８の回転周期が３．００秒から３．１０秒の間のいずれであっても，偏差値の絶対値を求めることができる。ただし，磁石１２８の回転周期を所定の精度で求めるためには，上記期間は，少なくとも予想される最大の磁石１２８の回転周期より長いことが好ましい。そして，上記式（９）により，移動平均値Ｈ３．００（ｔ）〜Ｈ３．１０（ｔ）のデータ列と，それら各データ列に対応する一次近似式Ｋ３．００（ｔ）〜Ｋ３．１０（ｔ）のそれぞれに対して，偏差値の絶対値の平均を求め，それぞれＡ３．００（ｔ）〜Ａ３．１０（ｔ）とする。
【００４６】
なお，上記偏差値の絶対値の平均Ａを求める式は，Ｈ３．００（ｔ）〜Ｈ３．１０（ｔ）のデータ列と，Ｋ３．００（ｔ）〜Ｋ３．１０（ｔ）の変位の度合いを示すものであれば良く，上記の式（９）に限定されるものではない。また，合計平均Ａの対象となる時刻は，上述の如くＡ３．００〜Ａ３．１０の全てで共通でなくても良く，例えばＡ３．０３を求める場合には，時刻ｔ−２．９，ｔ−２．８，ｔ−２．７，……，ｔ−０．１，ｔ，ｔ−２．９−０．０１，ｔ−２．９−０．０２，ｔ−２．９−０．０３におけるデータ値，すなわち移動平均値の算出を対象とした時刻におけるデータ値から求めてもよい。
【００４７】
次に，上述した工程で求めた偏差値の絶対値の平均から，プラズマの変動周期，すなわち磁石１２８の回転周期を算出する工程について説明する。上記式（９）により求めた偏差値の絶対値の平均Ａ３．００〜Ａ３．１０の中で，最小の値であるものを探し，その値に対する移動平均値算出期間を回転周期Ｔとする。すなわち，Ａ３．０３が最小値である場合には，回転周期Ｔは，３．０３秒に最も近いことになる。この理由としては，次のようなことが考えられる。
【００４８】
すなわち，図２中の期間Ａでは，処理室１０２内のプラズマＰの状態が相対的に安定しており，回転変動を含まないプラズマ光に基づく元波形は，相対的に滑らかな曲線となる。また，計算を行う期間は，微小であるため，その元波形は直線と近似することができる。さらに，移動平均値を算出する対象期間が磁石１２８の回転周期に完全に一致していれば，移動平均値は，回転変動を含まない元波形と一致する。従って，回転周期に近い期間で求めた移動平均値のデータ列は，そのデータ列から求めた一次近似式との変位（偏差）が小さくなることから，上記偏差値の絶対値の平均の最小の値が回転周期に最も近い値となる。
【００４９】
本実施の形態にかかる回転周期判定工程は，以上のように構成されており，変動するプラズマ光を一定周期でサンプリングした場合でも，プラズマ光の変動周期である磁石１２８の回転周期を正確に判定することができる。また，上述した演算で加重平均を用いているため，サンプリングデータ数を増やすことなく，かつ相対的に少ない演算処理で，磁石１２８の回転周期の判定に必要な所定の移動平均値を求めることができる。
【００５０】
（Ｂ）終点判定工程
次に，図５を参照しながら，プラズマ処理の終点判定工程について説明する。当該終点判定工程は，上述した変動周期判定工程で，プラズマの変動周期を求めた後，さらにプラズマ光の信号Ｐ（ｔ）をその磁石１２８の回転周期に基づいてサンプリングして，サンプリング信号から移動平均値のデータを求めてデータ列を記憶し，移動平均値のデータ列の変化を，例えばソフトウェア的に解析することによりエッチング処理の終点を判定する工程である。
【００５１】
図５は，移動平均値のデータ列を示したものであり，この移動平均値は，
0003563949.tif 000007
で表される式（１０）によって求めることができる。なお，この式（１０）は，上述した式（３）と同様の式である。また，式（１０）中のｍは（Ｔ／サンプリング周期）の商であり，ｎは（Ｔ／サンプリング周期）の余りの１０倍であり，０．１はサンプリング周期である。
【００５２】
式（１０）中の疑似サンプリングデータであるＳ（ｔ−０．１＊（ｍ−１）−０．０１）〜Ｓ（ｔ−０．１＊（ｍ−１）−０．０１＊ｎ）は，上述した式（４）〜式（６）と同様にして求めることができる。従って，式（１０）を用いて演算することにより，式（４）〜式（６）の如く，サンプリング周期の整数倍ではない期間の移動平均値を少ないサンプリングデータで求めることができると共に，疑似サンプリングデータを含む複数のデータの移動平均値を相対的に少ない演算処理で算出することができる。さらに，上記移動平均値を算出する際に，サンプリングデータを疑似サンプリングデータに対して重み付けすることにより，より一層正確な終点判定を行うことができる。
【００５３】
なお，このエッチング処理の終点判定工程は，上述した回転周期判定工程と組み合わせて実施される場合に限定されない。例えば，磁石１２８の回転周期が予めわかっている場合でも，その回転周期がサンプリング間隔の整数倍ではない場合に，上記式（１０）の如く終点判定工程と同様の方法で，移動平均値を求めることにより，磁石１２８の回転変動の影響を受けていないプラズマ光の信号を取り出すことができる。また，エッチング処理の終点の判定は，一般的に用いられている方法，例えば移動平均値のデータ列を微分し，微分値が所定値以上であればエッチング処理の終点と判定することにより行われる。
【００５４】
本実施の形態にかかる終点判定工程は，以上のように構成されており，上述した回転周期判定工程において判定された正確な磁石１２８の回転周期に基づき，その回転周期の期間の移動平均値のデータ列から得られる信号を用いて処理の終点を判定するため，正確な終点判定を行うことができ，エッチング処理を正確かつ確実に終了させることができる。
【００５５】
以上，本発明の好適な実施の一形態について，添付図面を参照しながら説明したが，本発明はかかる構成に限定されない。特許請求の範囲に記載された技術的思想の範疇において，当業者であれば，各種の変更例及び修正例に想到し得るものであり，それら変更例及び修正例についても本発明の技術的範囲に属するものと了解される。
【００５６】
例えば，上記実施の形態において，永久磁石から成る磁石１２８を回転させたときに生じるプラズマ光の信号の変動の変動周期を判定する方法について説明したが，本発明はかかる構成に限定されるものではなく，永久磁石に代えて電磁石を使用した場合や，磁石を回転ではなく往復運動させた場合にも本発明を適用することができる。さらに，本発明は，磁石を使用しなくてもプラズマに変動が生じる場合，例えば電極上で電界を印加する位置を順次周期的に変化させる場合にも適用することができる。
【００５７】
【発明の効果】
本発明によれば，プラズマが変動する際に，そのプラズマの変動周期やプラズマ処理の終点を，ソフトウェア的な演算処理のみで正確に判定することができるため，処理装置にその判定を行うためのハードウェアを形成する必要がない。その結果，処理装置の装置構成を容易にすることができると共に，既存の装置に対しても本発明を容易に実施することができる。また，本発明では，一定周期のサンプリング周期でサンプリングされたサンプリング信号に基づいて，上記変動周期や処理の終点を判定することができるため，装置のソフトウェアの都合上，固定周期をサンプリング周期として使用した方がよい場合や，固定周期を使用しなければならない場合などにも適用することができる。
【図面の簡単な説明】
【図１】本発明を適用可能なエッチング装置を示した概略的な断面図である。
【図２】図１に示したエッチング装置に適用される回転周期判定工程を説明するための概略的な説明図である。
【図３】図１に示したエッチング装置に適用される回転周期判定工程を説明するための概略的な説明図である。
【図４】図１に示したエッチング装置に適用される回転周期判定工程を説明するための概略的な説明図である。
【図５】図１に示したエッチング装置に適用される終点判定工程を説明するための概略的な説明図である。
【符号の説明】
１００            エッチング装置
１０２            処理室
１０６            下部電極
１０８            上部電極
１１８            高周波電源
１２０            演算制御器
１２８            磁石
１３４            検出窓
１３６            受光部
Ｗ                ウェハ
【図１】
0003563949.tif 000008
【図２】
0003563949.tif 000009
【図３】
0003563949.tif 000010
【図４】
0003563949.tif 000011
【図５】
0003563949.tif 000012