特許 6841980 装置診断装置、プラズマ処理装置及び装置診断方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6841980

(24)【登録日】2021年2月22日

(45)【発行日】2021年3月10日

(54)【発明の名称】装置診断装置、プラズマ処理装置及び装置診断方法

(51)【国際特許分類】

   H01L 21/3065 20060101AFI20210301BHJP

【ＦＩ】

   H01L21/302 103

   H01L21/302 101G

【請求項の数】8

【全頁数】12

(21)【出願番号】特願2020-532826(P2020-532826)

(86)(22)【出願日】2019年7月30日

(86)【国際出願番号】JP2019029762

(87)【国際公開番号】WO2020152889

(87)【国際公開日】20200730

【審査請求日】2020年6月12日

(73)【特許権者】

【識別番号】501387839

【氏名又は名称】株式会社日立ハイテク

(74)【代理人】

【識別番号】110001689

【氏名又は名称】青稜特許業務法人

(72)【発明者】

【氏名】梅田 祥太

(72)【発明者】

【氏名】玉置 研二

(72)【発明者】

【氏名】角屋 誠浩

(72)【発明者】

【氏名】石黒 正貴

【審査官】 鈴木 智之

(56)【参考文献】

【文献】 特開2014-022695（ＪＰ，Ａ）

【文献】 特開2002-270581（ＪＰ，Ａ）

【文献】 特開2013-041954（ＪＰ，Ａ）

【文献】 特開2003-347275（ＪＰ，Ａ）

【文献】 特開2010-165949（ＪＰ，Ａ）

【文献】 国際公開第2018/061842（ＷＯ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｈ０１Ｌ ２１／３０６５

(57)【特許請求の範囲】

【請求項1】

プラズマ処理装置の状態を診断する装置診断装置において、
第一のプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求め、前記予め求められた事前分布情報と、前記第一のプラズマ処理装置と異なる第二のプラズマ処理装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定し、前記推定された確率分布を用いて前記第二のプラズマ処理装置の状態を診断することを特徴とする装置診断装置。

【請求項2】

請求項１に記載の装置診断装置において、
前記第一のプラズマ処理装置は、複数であり、
前記第一のセンサの各々における予め求められた確率分布関数は、前記第一のプラズマ処理装置の各々に対して求められた確率分布関数の中で最も多くの前記第一のプラズマ処理装置に求められた確率分布関数であることを特徴とする装置診断装置。

【請求項3】

請求項２に記載の装置診断装置において、
前記第一のセンサの各々における予め求められた確率分布関数は、確率分布関数の候補の中から尤度を基に選択されることを特徴とする装置診断装置。

【請求項4】

プラズマ処理装置の状態を診断する装置診断装置において、
第一のプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求め、前記予め求められた事前分布情報と、第二のプラズマ処理装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定し、前記推定された確率分布に対する尤度である第一の尤度と、正規分布に対する尤度である第二の尤度とを比較し、
前記第一の尤度が前記第二の尤度より大きい場合、前記推定された確率分布を用いて前記第二のプラズマ処理装置の状態を診断し、
前記第二の尤度が前記第一の尤度より大きい場合、前記正規分布を用いて前記第二のプラズマ処理装置の状態を診断することを特徴とする装置診断装置。

【請求項5】

請求項１に記載の装置診断装置において、
前記確率分布は、マルコフ連鎖モンテカルロ法を用いて推定されることを特徴とする装置診断装置。

【請求項6】

請求項１に記載の装置診断装置において、
前記第二のプラズマ処理装置の状態の診断値として前記確率分布のプラズマ処理装置間差を出力するとともに前記確率分布の経時的な推移幅も出力することを特徴とする装置診断装置。

【請求項7】

試料がプラズマ処理される処理室と自装置の状態を診断する装置診断装置とを備えるプラズマ処理装置において、
前記装置診断装置は、自装置と異なるプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求め、前記予め求められた事前分布情報と、自装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定し、前記推定された確率分布を用いて前記自装置の状態を診断することを特徴とするプラズマ処理装置。

【請求項8】

プラズマ処理装置の状態を診断する装置診断方法において、
第一のプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求める工程と、
前記予め求められた事前分布情報と、第二のプラズマ処理装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定する工程と、
前記推定された確率分布を用いて前記第二のプラズマ処理装置の状態を診断する工程とを有することを特徴とする装置診断方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、装置診断装置、プラズマ処理装置及び装置診断方法に関する。

【背景技術】

【0002】

プラズマ処理装置では、半導体のウェハ上に微細形状を形成する。このため、物質をプラズマ化し、その物質の作用によりウェハ上の物質を除去するプラズマ処理を行う。

【0003】

プラズマ処理装置の製造時に複数の装置間で使用する共通部品に個体差（装置間差）があると、装置個体間のプラズマ状態の違いが生じる。このため、装置製造工場からの装置出荷前の最終検査時に原因部品を特定して交換あるいは調整するのに長時間を要するという問題があった。また、装置出荷後に半導体製造工場で想定寿命以上に使用した消耗部品の劣化が生じると同様にして装置間差が生じ得る。

【0004】

装置間差が生じると、所望の加工品質を得られなくなり、部品交換による予期せぬダウンタイムが生じる等の問題が起こる。そこで、調整の迅速化や予期せぬダウンタイム低減のために、装置間差を推定した上で、装置製造工程へのフィードバックや部品劣化診断を行う装置診断技術が求められる。

【0005】

このような装置診断技術として、例えば、特許文献１がある。特許文献１には、「異常検知装置は、観測値をまとめた要約値に対して統計モデリングを適用することにより、要約値からノイズを除去した状態を推測し、当該推測に基づき一期先の要約値を予測した予測値を生成する。異常検知装置は、予測値に基づき、監視対象装置の異常有無を検知する。」と記載されている。

【先行技術文献】

【特許文献】

【0006】

【特許文献1】ＷＯ２０１８／０６１８４２号公報

【発明の概要】

【発明が解決しようとする課題】

【0007】

特許文献１では、装置診断には一般に診断対象装置の状態センサの測定値（センサ値）の確率分布を用いるため、確率分布推定技術を要する。上述のように、プラズマ処理装置には装置間差が生じ得るため、特許文献１の方法では、高精度に異常検知するためには、装置毎にデータを大量に取得する必要がある。

【0008】

しかし、例えば、既存の工程に新規に装置を立上げる際や、新規工程に適用する際において、センサ値の大量取得は困難であるため、少量のセンサ値で確率分布を推定する方法が必要である。

【0009】

本発明の目的は、装置診断装置において、少量のセンサ値で確率分布を推定することにある。

【課題を解決するための手段】

【0010】

本発明の一態様の装置診断装置は、プラズマ処理装置の状態を診断する装置診断装置において、第一のプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求め、前記予め求められた事前分布情報と、前記第一のプラズマ処理装置と異なる第二のプラズマ処理装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定し、前記推定された確率分布を用いて前記第二のプラズマ処理装置の状態を診断することを特徴とする。

【0011】

本発明の一態様の装置診断装置は、プラズマ処理装置の状態を診断する装置診断装置において、第一のプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求め、前記予め求められた事前分布情報と、第二のプラズマ処理装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定し、前記推定された確率分布に対する尤度である第一の尤度と、正規分布に対する尤度である第二の尤度とを比較し、前記第一の尤度が前記第二の尤度より大きい場合、前記推定された確率分布を用いて前記第二のプラズマ処理装置の状態を診断し、前記第二の尤度が前記第一の尤度より大きい場合、前記正規分布を用いて前記第二のプラズマ処理装置の状態を診断することを特徴とする。

【0012】

本発明の一態様のプラズマ処理装置は、試料がプラズマ処理される処理室と自装置の状態を診断する装置診断装置とを備えるプラズマ処理装置において、前記装置診断装置は、自装置と異なるプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求め、前記予め求められた事前分布情報と、自装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定し、前記推定された確率分布を用いて前記自装置の状態を診断することを特徴とする。

【0013】

本発明の一態様の装置診断方法は、プラズマ処理装置の状態を診断する装置診断方法において、第一のプラズマ処理装置における第一のセンサにより取得された第一のセンサ値を用いて確率分布関数を含む事前分布情報を前記第一のセンサの各々に対して予め求める工程と、前記予め求められた事前分布情報と、第二のプラズマ処理装置における第二のセンサにより取得された第二のセンサ値とを基に前記第一のセンサの各々に対応する前記第二のセンサの各々における確率分布を推定する工程と、前記推定された確率分布を用いて前記第二のプラズマ処理装置の状態を診断する工程とを有することを特徴とする。

【発明の効果】

【0014】

本発明に一態様によれば、装置診断装置において、少量のセンサ値で確率分布を推定することができる。

【図面の簡単な説明】

【0015】

【図1】実施例のプラズマ処理装置の装置診断装置の構成図である。

【図2】プラズマ処理装置のセンサ値記憶部に格納するデータの例を示す図である。

【図3】装置群Ａのセンサ値を用いた装置診断装置の共通部の処理を示すフローチャートである。

【図4】装置診断装置の事前分布記憶部に格納するデータの例を示す図である。

【図5】装置診断装置の個別制御部の処理を示すフローチャートである。

【図6A】確率分布決定前における装置診断装置の確率分布記憶部に格納するデータの例を示す図である。

【図6B】確率分布決定後における装置診断装置の確率分布記憶部に格納するデータの例を示す図である。

【図7】装置診断装置の装置診断値記憶部に格納するデータの例を示す図である。

【図8】装置診断結果の表示画面の例を示す図である。

【図9】確率分布推定結果の表示画面の例を示す図である。

【図10】実施例のプラズマ処理装置の構成図である。

【発明を実施するための形態】

【0016】

以下、実施例について、図面を参照しながら説明する。なお、実施例を説明するための全図において、同一部には原則として同一符号を付し、その繰り返しの説明は省略する。

【0017】

（１）プラズマ処理装置
図１を参照して、プラズマ処理装置１の構成について説明する。

【0018】

図１に示すように、プラズマ処理装置１は、処理部１０と記憶部１１を有する。処理部１０は設定した加工条件に従い、処理室１０４の内部にプラズマ１０１を発生させてウェハ（試料１０２）を加工する。

【0019】

記憶部１１は、センサ値記憶部１２と管理値記憶部１３を有する。センサ値記憶部１２は、ウェハ加工中に装置のセンサ１０３の測定値を時系列データとして設定した加工ステップ毎に記憶する。例えば、センサ１０３は、状態センサ群を構成し温度や圧力を測定する。

【0020】

さらに、センサ値記憶部１２は、装置診断において使用するために、時系列データから主要な統計値（例えば、平均値や標準偏差）を算出して記憶する。本実施例においては、本統計値をセンサ値として以後使用する。

【0021】

図２を参照して、センサ値記憶部１２に格納するデータの例について説明する。

【0022】

図２では、統計値として平均値を用いている。センサ値の各列はセンサの種類に対応し、各行はウェハ１枚１枚に対応している。センサ値と共に、装置ＩＤや加工ステップＩＤ、加工条件ＩＤ、ウェハＩＤ等の情報も格納されている。装置ＩＤは、加工を行ったプラズマ処理装置１を特定する情報である。加工ステップＩＤ、加工条件ＩＤはそれぞれ加工ステップ、加工条件を特定する情報であり、監視対象の指定に用いる。ウェハＩＤには、加工したウェハを特定する情報を格納する。管理値記憶部１３は、加工日時や加工条件等の管理値を記憶する。

【0023】

（２）装置診断装置
図１を参照して、装置診断装置２の構成について説明する。

【0024】

図１に示すように、装置診断装置２は、共通部２０と個別制御部２３と記憶部２６を有する計算機である。装置診断装置２とプラズマ処理装置１の装置群３はネットワークを介して接続されており、相互にデータ通信が可能である。

【0025】

本実施例においては、装置群３は、センサ値が多量に蓄積された既存の装置群Ａと、装置立上げ等の理由からセンサ値が多量に蓄積されておらず装置診断対象である装置群Ｂとに分けて記載している。

【0026】

基準となるプラズマ処理装置１があれば必ずしも装置群Ａは複数のプラズマ処理装置１を有する必要はない。また、診断対象である装置群Ｂも必ずしも複数のプラズマ処理装置１を有する必要はない。さらに、装置群Ａと装置群Ｂで重複するプラズマ処理装置がある構成としてもよい。

【0027】

共通部２０は、共通分布関数選択部２１と事前分布設定部２２を有し、装置群Ｂの装置診断を実施する前に予め、装置群Ａに蓄積されたセンサ値からセンサ毎に事前分布情報を抽出して記憶部２６が有する事前分布記憶部２７に記憶する処理を行う。共通部２０の処理内容の例については、後述の（３）共通部の処理において説明する。

【0028】

個別制御部２３は、確率分布推定部２４と装置状態診断部２５を有する。確率分布推定部２４が、装置状態診断時に取得した装置群Ｂの個々のプラズマ処理装置のセンサ値と抽出済みの事前分布情報とから各センサ１０３の従う確率分布を事後分布として推定し、確率分布記憶部２８に記憶する。

【0029】

そして、推定したセンサ値の確率分布を基に、装置状態診断部２５が装置間差等の装置状態値を算出し、装置診断値記憶部２９に記憶する。個別処理部２３の処理内容の例については、後述の（４）個別処理部の処理において説明する。

【0030】

装置群３と装置診断装置２には出力部４０と入力部４１が接続されている。出力部４０は、例えばディスプレイやプリンタ等であり、記憶部２６の情報を基にユーザに対してグラフィカルに情報を出力する装置である。表示例については、後述の（５）出力部による表示例において説明する。入力部４１は、例えば、マウスやキーボード等の、ユーザの操作による情報入力を受け付ける入力装置である。

【0031】

（３）共通部の処理
図３を参照して、装置診断装置２の共通部２０で行われる装置群Ａのセンサ値を用いて事前分布を設定する処理の例について説明する。

【0032】

共通部２０の処理を実行する前には、予め装置群Ａの各プラズマ処理装置１において、監視対象の加工ステップにおけるプラズマ処理中の履歴であるセンサ値をセンサ値記憶部１２に格納しておく。

【0033】

監視対象の加工ステップとして、複数の工程で共通に行う処理部１０の状態を整えるためのプラズマ処理（例えば、エージング処理やクリーニング処理）を指定する。装置群Ａのセンサ値記憶部１２から、指定した加工ステップＩＤのセンサ値を取得する。（Ｓ１０１）
このように複数の工程で共通に行うプラズマ処理のセンサ値を用いることで、例えば装置群Ｂが新規工程への適用の際であっても、装置群Ａが蓄積したセンサ値を事前分布の設定に利用することができる。

【0034】

次に、事前分布設定部２２は、各センサ１０３、各確率分布関数候補に対してＳ１０３〜Ｓ１０４の処理を実行する。ここで、確率分布関数候補とは、プラズマ処理装置１の各センサ１０３が従い得る確率分布関数を予め候補として設定したものである。例えば、正規分布、歪正規分布、混合正規分布、コーシー分布等を確率分布関数候補として設定しておく（Ｓ１０２）。

【0035】

事前分布設定部２２は、当該センサ１０３のセンサ値に対して当該確率分布関数候補の確率分布パラメータを推定する。確率分布パラメータは、例えば正規分布であれば平均値と標準偏差に相当する値であり、確率分布関数毎に異なる種類を有する。確率分布パラメータの推定には、例えば、マルコフ連鎖モンテカルロ法（Ｍａｒｋｏｖ Ｃｈａｉｎ Ｍｏｎｔｅ Ｃａｒｌｏ ｍｅｔｈｏｄ、ＭＣＭＣ法）を利用する。ＭＣＭＣ法は、確率分布パラメータを確率変数とみなした上で、確率分布パラメータの事後分布に比例する確率分布パラメータの事前分布と尤度の積から乱数サンプルを大量に発生させて確率分布パラメータの事後分布を推定する手法である（Ｓ１０３）。

【0036】

推定した確率分布パラメータは確率分布として得られるが、例えば事後確率が最大となる値から一意に値を定める。このときの確率分布を用いて、事前分布設定部２２は得られている装置群Ａのセンサ値に対する当て嵌まりの度合いである対数尤度を算出する（Ｓ１０４）。

【0037】

最後に、共通分布関数選択部２１は、各センサ１０３に対して、装置群Ａで共通する確率分布関数として、対数尤度が最大となった確率分布関数を確率分布関数候補の中から選択する。また、選択した確率分布関数に関する確率分布パラメータの推定値を、装置群Ｂのセンサ値の確率分布を推定する際の事前分布として事前分布記憶部２７に格納する。また、正規分布に関する確率分布パラメータの推定値も対数尤度が最大となった確率分布関数と併せて事前分布記憶部２７に格納する（Ｓ１０５）。

【0038】

図４を参照して、事前分布記憶部２７に格納するデータの例について説明する。

【0039】

図４に示すように、センサ毎に事前分布の情報を格納する形態となっている。確率分布関数の行は、図３のＳ１０５で選択した確率分布関数の名称を格納する。確率分布パラメータの行は、図３のＳ１０３で推定した確率分布パラメータの事後分布を用いて設定した各確率分布パラメータの事前分布を格納する。例として、事後分布の平均と標準偏差を有する正規分布を各確率分布パラメータの事前分布として設定している。対数尤度の行は、図３のＳ１０４で算出した対数尤度を格納している。

【0040】

このように、装置群Ａのセンサ値を用いて装置群で共通性が最大の確率分布関数を事前分布情報とすることで、装置群Ｂに対しても共通に使える事前分布を抽出でき、非正規分布であっても少量のセンサ値での確率分布推定を可能とする。

【0041】

（４）個別制御部の処理
図５を参照して、装置診断装置２の個別制御部２３で行われる確率分布の推定と装置診断の処理の例について説明する。

【0042】

装置群Ｂの診断対象とするプラズマ処理装置１のセンサ値記憶部１２から、図３のＳ１０１で指定した加工ステップＩＤのセンサ値を取得する（Ｓ２０１）。

【0043】

確率分布推定部２４は、取得したセンサ値の各センサに対して、Ｓ２０３〜Ｓ２０５の処理を実行する（Ｓ２０２）。

【0044】

まず、事前分布記憶部２７から当該センサに対応する事前分布情報を取得する（Ｓ２０３）。

【0045】

次に、取得した事前分布情報を事前分布として設定した上で、取得したセンサ値を用いてＭＣＭＣ法により、確率分布パラメータの事後分布を推定する。推定した確率分布パラメータを用いて、図３のＳ１０４と同様にしてセンサ値に対する対数尤度を算出する。本処理は、図３のＳ１０５で選択した確率分布関数の場合と、正規分布の場合両方に対して実施する（Ｓ２０４）。

【0046】

次に、図３のＳ１０５で選択した確率分布関数で推定した確率分布と正規分布で推定した確率分布とで対数尤度を比較し、対数尤度が大きい確率分布を当該センサに対する推定結果とし、確率分布記憶部２８に格納する（Ｓ２０５）。

【0047】

このように、装置群Ｂのセンサ値も併せて確率分布関数を決定することで推定結果のロバスト性を高めることができる。

【0048】

図６Ａ及び図６Ｂを参照して、各センサの確率分布決定前後における確率分布記憶部２８に格納するデータの例について説明する。

【0049】

図６Ａ及び図６Ｂに示すように、センサ名の各列に示すセンサ毎に確率分布情報を格納する。格納する確率分布情報は、図４と同様に、確率分布関数と確率分布パラメータ及び対数尤度である。

【0050】

確率分布パラメータは図４と同様に推定した事後分布を格納するものとしても良いし、図６Ａと図６Ｂのように例えば事後確率が最大となる値のように一意に定まる値を格納しても良い。図６Ａの各センサにおける正規分布を含む二つの確率分布関数のうち、対数尤度の大きい確率分布を決定した後が図６Ｂに相当する。

【0051】

最後に、確率分布記憶部２８に格納したデータを用いて装置診断値を算出する（Ｓ２０６）。例えば、装置間差を装置診断値として算出する場合は、まず、基準となるプラズマ処理装置１と診断対象のプラズマ処理装置１の確率分布情報を確率分布記憶部２８から取得する。取得した確率分布間の距離を装置間差診断値として装置状態診断値記憶部２９に格納する。確率分布間の距離指標としては、Ｋｕｌｌｂａｃｋ−ＬｅｉｂｌｅｒダイバージェンスやＪｅｎｓｅｎ−Ｓｈａｎｎｏｎダイバージェンス等を用いる。

【0052】

図７を参照して、装置間差を装置診断値として算出する場合における、装置診断値記憶部２９に格納するデータの例について説明する。

【0053】

図７に示すように、装置ＩＤ行には、比較する２台のプラズマ処理装置１を特定するＩＤを格納する。装置間差診断値行に、センサ毎に算出した装置間差診断値を格納する。

【0054】

（５）出力部による表示例
出力部４０は、記憶部１１や記憶部２６に格納された情報を用いて、装置状態の診断結果や確率分布の推定結果を表示する。

【0055】

図８を参照して、装置間差診断結果の表示画面Ｄ１００の例について説明する。

【0056】

図８に示すように、Ｄ１０２に注目する装置ＩＤを入力する。これにより、装置状態値記憶部２９に格納された該当する装置ＩＤの装置間差診断値を取得し、Ｄ１０４のように、センサ毎にグラフで表示する。この結果、基準となるプラズマ処理装置１との装置間差が大きいセンサを判断することができ、調整等に利用することができる。

【0057】

センサ１０３は通常、測定対象部品や測定対象項目等によって複数の群に分類することができる。センサ群として予め登録しておき、Ｄ１０３のように、センサ群に属する各センサの装置間差診断値を積算し、センサ群毎に装置間差診断値を表示する。こうすることで、調整対象の部品が明確となる。

【0058】

全センサの装置間差診断値を積算することで、Ｄ１０１のように、プラズマ処理装置１毎の比較も表示する。

【0059】

図９を参照して、センサ値の分布と確率分布記憶部２８に格納された確率分布の推定結果を確認する画面Ｄ１０５の例について説明する。

【0060】

図９に示すように、装置ＩＤやセンサ名を入力し、確認する対象を決定した後、グラフを表示する。複数の装置ＩＤを入力することで、比較表示することができる。ヒストグラムがセンサ値記憶部１２に格納されたセンサ値の測定値のヒストグラムであり、実線が確率分布記憶部２８に格納された確率分布推定結果である。

【0061】

例えば、図９では装置ＩＤがＣ１、Ｃ４の装置のセンサ名Ｘ２のセンサ値のヒストグラムと確率分布の推定結果を表示している。これにより、使用者は推定された確率分布が妥当なものであるかどうかを確認することができる。また、例えば、比較する二つの確率分布が平行移動している場合であればセンサ１０３の初期化を検討する等、装置診断結果に対する対策を考案することに活用することができる。

【0062】

このように、図８及び図９に示すように、出力部４０は、プラズマ処理装置１の状態の診断値として確率分布のプラズマ処理装置間差を出力するとともに確率分布の経時的な推移幅も出力する。

【0063】

上記実施例によれば、プラズマ処理装置において、例えば装置立上げ時や新規工程への適用時や出荷前の調整時に、診断対象装置のセンサ値が多量に取得できていない場合においても、既存の装置群で共通する事前分布情報を抽出した上で、診断対象の装置で新規に取得したセンサ値と併せて非正規分布を含めた確率分布を推定することで装置診断を可能とする。

【0064】

以上、実施例について説明したが、本発明は上記実施例に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能である。

【0065】

例えば、図１では、装置診断装置２とプラズマ処理装置１がネットワークを介して接続された構成を示したが、本発明は上記構成に限定されず、図１０に示すように、プラズマ処理装置１が装置診断装置２を備える構成であっても良い。この場合、プラズマ処理装置１は、試料１０２がプラズマ処理される処理室１０４と自装置の状態を診断する装置診断装置２とを備える。

【0066】

図１０に示すように、プラズマ処理装置１は、図１に示す処理部１０と記憶部１１を有する。処理部１０は設定した加工条件に従い、処理室１０４の内部にプラズマ１０１を発生させてウェハ（試料１０２）を加工する。

【0067】

また、装置診断装置２は、図１に示す構成と同様に、共通部２０と個別制御部２３と記憶部２６を有する。

【符号の説明】

【0068】

１ プラズマ処理装置
２ 装置診断装置
３ 装置群
１０ 処理部
１１ 記憶部
１２ センサ値記憶部
１３ 管理値記憶部
２０ 共通部
２１ 共通分布関数選択部
２２ 事前分布設定部
２３ 個別制御部
２４ 確率分布推定部
２５ 装置状態診断部
２６ 記憶部
２７ 事前分布記憶部
２８ 確率分布記憶部
２９ 装置診断値記憶部
４０ 出力部
４１ 入力部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6A】

【図6B】

【図7】

【図8】

【図9】

【図10】