(57)【特許請求の範囲】
【請求項１】
互いにネットワークで接続された複数の電子顕微鏡装置の一つである第一の電子顕微鏡装
置を用いて寸法測定対象試料を走査し、前記寸法測定対象試料のパターンの画像を取得す
るステップと、
前記取得した前記寸法測定対象試料のパターンの画像の実信号波形と、予め格納された前
記複数の電子顕微鏡装置の装置特性パラメタのうち前記第一の電子顕微鏡装置の装置特性
パラメタを用いてシミュレーションにより作成したシミュレーション波形ライブラリと、
を比較して、前記寸法測定対象試料のパターンの推定断面形状を得るステップと、
前記得られた寸法測定対象試料のパターンの推定断面形状と、前記複数の電子顕微鏡装置
のうち前記第一の電子顕微鏡以外の複数の電子顕微鏡装置各々の装置特性パラメタと、を
用いて、前記第一の電子顕微鏡以外の複数の電子顕微鏡装置各々に対応する複数の模擬信
号波形を生成するステップと、
前記実信号波形に基づいて算出した前記寸法測定対象試料のパターンの寸法と前記複数の
模擬信号波形各々に基づいて算出した複数の前記寸法測定対象試料のパターンの寸法とを
用いて、前記第一の電子顕微鏡を含めた複数の電子顕微鏡装置各間の機差を算出し、算出
された前記複数の電子顕微鏡装置各間の機差を機差情報として機差データベースに格納す
るステップと、
前記複数の電子顕微鏡装置の少なくとも一つを用いて得られた前記寸法計測対象試料のパ
ターンの寸法計測結果を、前記機差データベースに格納された機差情報を用いて補正して
、前記寸法計測対象試料のパターンの寸法を計測するステップと、
を有することを特徴とするパターン寸法計測方法。
【請求項２】
請求項１記載のパターン寸法計測方法であって、
前記推定断面形状を得るステップで用いるシミュレーション波形ライブラリは、
前記予め格納された前記複数の電子顕微鏡装置の装置特性パラメタのうち前記第一の電子
顕微鏡装置の装置特性パラメタに加え、前記寸法測定対象試料のパターンの近隣パターン
情報を用いて作成されることを特徴とするパターン寸法計測方法。
【請求項３】
請求項１又は２記載のパターン寸法計測方法であって、
前記予め格納された前記複数の電子顕微鏡装置の装置特性パラメタは、前記複数の電子顕
微鏡装置各々の照射ビームのチルト角およびビーム径、並びに、試料からの２次電子の検
出効率を含むことを特徴とするパターン寸法計測方法。
【請求項４】
請求項３記載のパターン寸法計測方法であって、
前記予め格納された前記複数の電子顕微鏡装置の装置特性パラメタの一つである前記複数
の電子顕微鏡装置各々の照射ビームのチルト角は、装置特性評価用試料として異方性エッ
チングにより形成したピラミッドパターンを前記複数の電子顕微鏡装置各々により撮像し
、撮像して得た各撮像画像上で前記ピラミッドパターンの稜線を各々検出し、検出された
ピラミッドパターンの稜線の幾何学歪みからビームチルト角を計測して得たものであるこ
とを特徴とするパターン寸法計測方法。
【請求項５】
請求項３記載のパターン寸法計測方法であって、
前記予め格納された前記複数の電子顕微鏡装置の装置特性パラメタである前記複数の電子
顕微鏡装置各々の照射ビームのビーム径と検出効率として、装置特性評価用試料としてス
テップエッジサンプルを用い、ＳＥＭシミュレーションによって得たビーム径ゼロにおけ
るシミュレーション波形と前記ステップエッジサンプルのＳＥＭ像から得た実波形を合致
させるような，ビーム強度分布を表現する関数と，明るさ変換係数を用いることを特徴
とするパターン寸法計測方法。
【請求項６】
請求項１乃至５のいずれかに記載のパターン寸法計測方法であって、
前記寸法測定対象試料のパターンの推定断面形状を得るステップでは、
モデルベースライブラリ法を用いて前記推定断面形状を得ることを特徴とするパターン寸
法計測方法。
【請求項７】
互いにネットワークで接続された複数の電子顕微鏡装置と、
前記複数の電子顕微鏡装置各々の装置パラメタを格納する装置パラメタデータベースと、
前記複数の電子顕微鏡装置の一つである第一の電子顕微鏡装置を用いて寸法測定対象試料
を走査することにより取得した寸法測定対象試料のパターンの画像の実信号波形と、前記
装置パラメタデータベースに格納された前記第一の電子顕微鏡装置の装置特性パラメタを
用いてシミュレーションにより作成したシミュレーション波形ライブラリと、を比較して
、前記寸法測定対象試料のパターンの推定断面形状を得、前記得られた寸法測定対象試料
のパターンの推定断面形状と、前記複数の電子顕微鏡装置のうち前記第一の電子顕微鏡以
外の複数の電子顕微鏡装置各々の装置特性パラメタと、を用いて、前記第一の電子顕微鏡
以外の複数の電子顕微鏡装置各々に対応する複数の模擬信号波形を生成するコンピュータ
と、
前記コンピュータにより生成された前記複数の模擬信号波形各々に基づいて算出した複数
の前記寸法測定対象試料のパターンの寸法と前記実信号波形に基づいて算出した前記寸法
測定対象試料のパターンの寸法とを用いて算出された、前記第一の電子顕微鏡を含めた複
数の電子顕微鏡装置各間の機差を機差情報として格納する機差データベースと、
前記機差データベースに格納された機差情報を用いて寸法計測結果を補正する処理を含む
レシピが格納されたレシピデータベースと、
を有することを特徴とするパターン寸法計測システム。
【請求項８】
請求項７記載のパターン寸法計測システムであって、
前記機差情報を用いて寸法計測結果を補正する処理を含むレシピは、前記コンピュータの
レシピ作成のためのＧＵＩ上で、機差補正を実行する項目を選択することにより作成され
ることを特徴とするパターン寸法計測システム。
【請求項９】
装置特性評価用試料を電子顕微鏡装置により撮像して得たパターン画像を用いて装置特性
パラメタを定期又は不定期に算出するステップと、
前記計測された装置特性パラメタを算出した日時とリンクさせてデータベースに格納する
ステップと、
前記装置特性評価用試料のパターン画像の実信号波形と、前記装置特性パラメタを用いて
作成されたシミュレーション波形ライブラリとを比較して、前記装置特性評価用試料のパ
ターンの推定断面形状を得るステップと、
前記推定断面形状と前記装置特性パラメタを用いて模擬信号波形を生成するステップと、
前記生成された模擬信号波形を用いて前記装置特性評価用試料のパターンの寸法計測を行
い、当該寸法計測結果を格納するステップと、
前記格納された寸法計測結果を前記装置特性パラメタの経時変化と合わせて表示するステ
ップと、
を有することを特徴とする電子顕微鏡装置の経時変化のモニタ方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は半導体装置等の微少寸法を計測するのに用いられる電子顕微鏡装置のパターン寸法計測に関する。
【背景技術】
【０００２】
非特許文献１に述べられているように，半導体プロセスにおいてパターン寸法管理ツールとしても最も普及しているのは，走査電子顕微鏡（以下ＳＥＭと略す）を半導体計測専用に特化した測長ＳＥＭである。図２に測長ＳＥＭの原理を示す。電子銃９０１から放出された電子ビーム９１０は集束レンズ９０２で細く絞られ，対物レンズ９０４で焦点を試料９００の表面に合わせた状態で制御装置９０６で制御された偏向器９０３により試料９００上を２次元的に走査される。電子ビーム９１０の照射によって試料９００から発生した二次電子９２０を検出器９０５で捕らえて制御装置９０６で信号処理することで、ＣＲＴ９０７に表示されているような電子線像が得られる。二次電子はパターンエッジ部でより多く発生するため，電子線像は，ＣＲＴ９０７で表示されているように、パターンエッジに相当する部分が明るい画像となる。測長ＳＥＭにおいては，電子線像上でのエッジ間距離l（画素）に画素サイズp（nm/画素）を乗じる（l×p）ことにより寸法が求められる。
【０００３】
測長ＳＥＭにおける測長処理の一例が特許文献１に記載されている。特許文献１の開示例においては，測定対象配線を撮像した画像内の局所領域から，配線の信号波形を配線の長手方向に加算平均した投影波形を作成し，この波形において検出した左右の配線エッジ間の距離として配線寸法を算出している。
【０００４】
エッジ間距離lを自動計算するためのエッジ位置検出方法として種々の方法が提案されているが，ここでは，広く一般的に用いられている，しきい値法について述べる。
【０００５】
しきい値法は，例えば、特許文献２に開示されている。図３で示すように，左右のパターンエッジに相当する信号量の大きいピーク部分を，それぞれ左ホワイトバンド（左ＷＢ），右ホワイトバンド（右ＷＢ）と呼ぶことにする。しきい値法は，左右ＷＢそれぞれで，Ｍａｘ値，Ｍｉｎ値を求め，これらを所定の比率th（％）で内分するしきい値レベルを算出し，しきい値と信号波形の交差位置をエッジ位置と定義する。
【０００６】
半導体デバイスパターンの微細化により，計測精度に対する要求は年々厳しさを増している。測長ＳＥＭの場合，一つの半導体製造ラインに複数台の測長ＳＥＭが投入され，これらを併用することが多いため，単一の測長ＳＥＭの計測再現性のみならず，装置間の計測再現性，すなわち，複数台の測長ＳＥＭの計測値間の差（以下，機差と記す）が問題となる。機差に対する要求については、例えば、非特許文献２で示すように、DRAMハーフピッチ45nm世代での機差として＜0.18nm，同32nm世代での測長値機差として＜0.13nm，同22nm世代での機差として＜0.09nmが要求されている。
【０００７】
また，非特許文献２には機差の測定方法が開示されている。前述のように，機差は複数台のＳＥＭの計測値の差と定義されるが，測長SEMは，電子線照射によってコンタミネーション（試料の電子線を照射した部分にアモルファスカーボン膜が付着する現象），試料帯電，あるいは，シュリンク（電子線照射によるレジストパターンの縮み）といった試料変化を引き起こすため，同一計測対象を複数の装置で計測し、単純にそれらの計測値の差を求めても，機差を正しく測定したことにはならない。計測値の差には，試料変化の影響が含まれてしまうからである。非特許文献２にはこの問題に対処するための具体的な方法が示されている。ここでは，二つの代表的な手法，ABBA法とグレーティング法を説明する。
【０００８】
（１）ABBA法
装置Aに対する装置Bの機差をδ（nm），装置Aでの計測が引き起こすコンタミネーション，試料帯電，シュリンク起因の計測値変化量をc1（nm）と，装置Bでの計測が引き起こすによる計測値変化量をc2（nm）とする。次に示す，step1〜step3にて機差を求める。
Step1：試料上の複数箇所（通常は数十箇所）について，装置Aで計測した後に装置Bで計測を行い（ABシーケンス），装置Bの計測値の平均値から装置Aの計測値の平均値を減じた値δ１（nm）を求める。δ，δ１，C1の間には数式１の関係が成り立つ。
（数１）δ１＝δ＋C1
Step2：試料上の複数箇所（通常は数十箇所）について，装置Bで計測した後に装置Aで計測を行い（BAシーケンス），装置Aの計測値の平均値から装置Bの計測値の平均値を減じた値δ２（nm）を求める。δ，δ2，C2の間には数式２の関係が成り立つ。
（数２）δ２＝-δ＋C2
Step3：仮にC1＝C2が成り立つならば（装置Aと装置Bが引き起こす試料変化が同程度であるならば），数式１と数式２より，
（数３）δ＝（δ１-δ２）/2
にて，δが求められる。
ABBA法はc1＝c2の仮定が成り立ち（装置Aと装置Bが引き起こす試料変化が同程度とみなせるということ），かつ，個々の計測の計測誤差を考慮して，計測点数を十分に多くすれば正しい機差を求めることができる。
【０００９】
(2)グレーティング法
グレーティング法では，電子線照射による試料変化の影響を避けるため，同じ箇所を複数回計測せずに、図４のように、４台の装置（A，B，C，D）の機差を測定する際の計測箇所の配置をずらし、各装置は多数のフレッシュな箇所（一度も撮像されたことない箇所）を計測して，それらの平均値の差を機差とする。
この際、パターンの出来映えの面内分布の影響を受けづらいよう，各装置の計測箇所は偏らないように配置される。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開平１１−３１６１１５号公報
【特許文献２】特開昭５５−７２８０７号公報
【非特許文献】
【００１１】
【非特許文献１】社団法人  日本半導体製造装置協会（ＳＥＡＪ）平成１７年度半導体製造装置技術ロードマップ報告書  第５編計測
【非特許文献２】International SEMATECH Manufacturing Initiative， “Unified Advanced Critical Dimension Scanning Electron Microscope (CD-SEM) Specification for sub-90 nm Technology (2008 Version)，” ISMI Technology Transfer #04114595E-ENG
【発明の概要】
【発明が解決しようとする課題】
【００１２】
先に述べた機差に対する要求に応えるべく，測長SEMベンダ各社では機差低減のための努力が継続されているが，測長SEMの運用状況によっては，種々の型式の装置が混在していることもあるため，型式間に性能の違い（例えば分解能の違いなど）がある以上，機差の発生は避けがたい。こうした場合の対策として、本番の計測に先立ち機差を実測して（ABBA法やグレーティング法による）装置別機差として登録しておき，本番の計測時には各装置による寸法計測値から，上記の装置別機差を減ずることで機差を補正するという方法（以下，これを機差実測方式と呼ぶ）が採用されることもある。
【００１３】
しかしながら，ABBA法では、例えば、対象装置が３台（A，B，C）の場合は，ABBAシーケンスによって装置A，B間の機差を測定し、さらに，ACCAシーケンスによって装置A，C間の機差を測定するといった手順が必要であるため、機差測定には多大な手間を要する。また、グレーティング法では、計測箇所が偏らないようにしても，エッジラフネスの影響は残っているので，場所によってパターン寸法に差があり、この差に機差が埋もれないようにすべく、（エッジラフネスの程度によるが）計測点数を相当大きくすることが必要である。グレーティング法は，ABBA法での前提条件c1＝c2が成り立たない場合にも適用可能という利点はあるが，ABBA法と異なり，場所による寸法差の影響をダイレクトに受けるため，原理上，ABBA法よりずっと多くの計測点数を必要とする。図４では装置ごとの計測点数は２５点であるが，前述の機差に対する要求を考慮すると，１００点以上の計測点数が必要になることも少なくない。
このように、ABBA法にしてもグレーティング法にしても機差測定には多大な手間を要す上，ABBA法には試料変化という誤差要因が，グレーティング法にはパターンのエッジラフネスという誤差要因が存在するため，必ずしも正しい機差が測定されるとは限らない。
【００１４】
さらに，機差は，計測対象の断面形状（高さや側壁傾斜角など）やパターン密度によって異なるため，機差実測方式を厳密に適用しようとするならば，全ての計測対象パターンについて，機差を実測することが必要となるという非常にやっかいな問題がある。これは，機差が，装置間の入射ビーム特性や２次電子検出特性の違いによって生じるということに起因する。例として，図５（Ａ）〜（Ｄ）に示す各種のパターン１２０、１２１，１２２、１２３を，ビーム径が異なる２台の測長SEMを用いて計測した場合の機差のシミュレーション結果を図６〜図８に示す。
【００１５】
シミュレーションでは，ビーム径が５nm（σ=2.5nmのガウシアンビーム）の場合のSEM信号波形と，ビーム径が６nm（σ＝3nmのガウシアンビーム）の場合のSEM信号波形をモンテカルロシミュレーションにて求め，それぞれの波形に対してしきい値法（図３参照）を適用し寸法を算出し，ビーム径が6nmの場合の寸法計測値とビーム径が５nmの場合の寸法計測値の差を機差として求めた。
【００１６】
図６は，(a)パターン高さ80nm，側壁傾斜角０°の孤立ラインパターン（図５（Ａ）のパターン１２０のようなパターン）と，(b)パターン高さ80nm，側壁傾斜角1０°の孤立ラインパターン（図５（Ｂ）のパターン１２１のようなパターン）について，線幅が15nm〜45nmにおける機差を算出した結果である。この結果から，側壁傾斜角が異なると機差が異なるということが分かる。機差実測方式を適用するためには，側壁傾斜角別に機差を実測することが必要といえる。
【００１７】
図７は，上記(a)と，(c)パターン高さ80nm，側壁傾斜角０°のパターンピッチが45nmのライン・アンド・スペースパターン（図５（Ｃ）のパターン１２２のようなパターン）について，線幅が15nm〜45nmにおける機差を算出した結果である。この結果から，パターン密度が異なると機差が異なるということが分かる。なお，線幅が25nm以下では，(a)と（c）の機差の違い小さいのは，本シミュレーションではライン・アンド・スペースパターンでのパターンピッチを45nmに固定したため，線幅が小さくなると隣接ラインまでの距離が遠くなって，孤立ラインパターンに近い状況になるためである。線幅やパターンピッチに依存して，機差の変化傾向が異なることは明らかであり，機差実測方式を適用するためには，パターンピッチ，サイズ別に機差を実測することが必要といえる。
【００１８】
図８は，上記(a)と，(d)パターン高さ80nm，側壁傾斜角０°の孤立スペースパターン（図５（Ｄ）のパターン１２３のようなパターン）について，線幅が15nm〜45nmにおける機差を算出した結果である。ラインパターンの機差(a)と，スペースパターンの機差（d）の符号が逆転するのは、図９に示すように、エッジ部においては２次電子信号波形のピーク部が一定の広がりを有し，ビーム径大→ピーク部の広がり大という関係があるからである。スペース計測の他，ホールパターンの径を計測するような場合も，(d)と同様の状況となる。このように，計測対象のトポロジーによっても機差は異なる。機差実測方式を適用するためには，トポロジーに応じて機差を実測することが必要といえる。
【００１９】
図６〜図８は，装置間のビーム径の違いによって機差が生じている場合をシミュレートした結果であるが，他の機差要因，例えば，入射ビームのチルト角の装置間差，あるいは，２次電子の検出効率の装置間差が存在する場場合，やはり，計測対象の断面形状（高さや側壁傾斜角など）やパターン密度によって機差が異なるが，その異なり方は，上記のシミュレーション結果とは異なった傾向となる。
【００２０】
以上のように、機差実測方式は，単一のパターンですら大変な手間を要し，加えて，正確な測定が難しい機差を計測対象パターンごとに実施しなければならず，特に，計測対象パターンが多岐に及ぶような状況では，現実的な方法とは言い難い。
【００２１】
本発明は，型式が異なる等の理由で，測長SEMのハードウエハの調整では機差をなくすことが困難であり，また，機差実測方式では対応が困難な，計測対象が多岐に及ぶようなケースにも適用が可能な機差補正方法、パターン寸法計測方法、パターン寸法計測システムおよび電子顕微鏡装置の経時変化のモニタ方法を提供するものである。
【課題を解決するための手段】
【００２２】
本願において開示される発明のうち代表的なものの概要を簡単に説明すれば次のとおりである。
（１）本発明は，次の（Ａ）〜（Ｃ）のステップにて，機差を補正，すなわち，機差のない計測値を得るものである。
（Ａ）装置特性性をモデル化するステップ：各測長SEM（#1号機〜＃n号機）において，専用試料を撮像し，画像を解析することで装置特性パラメタを求める。例えば，次の（ａ）〜（ｃ）を行う。
（ａ）ビーム径：ナイフエッジパターンを撮像し、エッジ部におけるピーク波形のぼけの程度からのビーム径を算出する。（ｂ）ビームチルト角：ピラミッド形状パターン（四角錐形状パターン）を撮像し，画像上のピラミッド形状の歪みからビームチルト角を算出する。（ｃ）明るさ変換係数：断面形状が既知（原子間力顕微鏡による計測，あるいは，断面ＳＥＭ観察によって形状情報を得る）のサンプルを用い，その断面形状におけるSEM信号波形をSEMシミュレーションにて求め，これと，測長ＳＥＭによる実信号波形とを比較することで，シミュレーション波形と実波形の間の明るさ変換係数を算出する。
（Ｂ）機差を推定するステップ：基準装置（例えば，測長SEM#１号機）で計測対象パターンを撮像し，ＭＢＬ法（モデル・ベース・ライブラリ法，すなわち，計測対象パターンの電子線信号波形を予め作成しておいた試料の断面形状とＳＥＭ信号波形とを関連づけるライブラリに当てはめることによって，試料の断面形状を推定する方法）にて計測対象パターンの断面形状を推定する。この断面形状を入力として，装置パラメタ（上記（Ａ）で算出）を反映させたSEMシミュレーションを行って，#2号機〜＃n号機のSEM信号波形を求め，所定のエッジ検出手法（前記のしきい値法など）を適用して寸法計測を行い，これと，測長SEM#１号機の寸法計測結果を比較することで，各測長SEMの基準装置に対する機差を推定する。
（Ｃ）測長のステップ：各測長SEM（#1号機〜＃n号機）で通常通りに，計測対象パターンの計測を行い，非基準装置（#２号機〜＃n号機）については，得られた計測結果から，上記（Ｂ）で推定した機差を減ずる。
従来の手法である機差実測方式の場合，大変な手間と時間を要し，かつ，正確な測定が難しい機差測定を，全計測対象パターンについて実施しなければならないのに対し，本発明では，上記（Ａ）のステップで予め装置特性パラメタを計測しておけば，上記（Ｂ）のステップで機差が推定される。
また、本発明によれば，上記（Ｂ）で推定した機差を補正値として用いることによって，計測対象が多岐に及ぶようなケースにも適用が可能な機差補正方法を提供することが可能となる。これにより，型式が異なる測長SEMが混在している場合であっても，これらの混用が可能となる。 
（２）互いにネットワークで接続された複数の電子顕微鏡装置の一つである第一の電子顕微鏡装置を用いて寸法測定対象試料を走査し、前記寸法測定対象試料のパターンの画像を取得するステップと、前記取得した前記寸法測定対象試料のパターンの画像の実信号波形と、予め格納された前記複数の電子顕微鏡装置の装置特性パラメタのうち前記第一の電子顕微鏡装置の装置特性パラメタを用いてシミュレーションにより作成したシミュレーション波形ライブラリと、を比較して、前記寸法測定対象試料のパターンの推定断面形状を得るステップと、前記得られた寸法測定対象試料のパターンの推定断面形状と、前記複数の電子顕微鏡装置のうち前記第一の電子顕微鏡以外の複数の電子顕微鏡装置各々の装置特性パラメタと、を用いて、前記第一の電子顕微鏡以外の複数の電子顕微鏡装置各々に対応する複数の模擬信号波形を生成するステップと、前記実信号波形に基づいて算出した前記寸法測定対象試料のパターンの寸法と前記複数の模擬信号波形各々に基づいて算出した複数の前記寸法測定対象試料のパターンの寸法とを用いて、前記第一の電子顕微鏡を含めた複数の電子顕微鏡装置各間の機差を算出し、算出された前記複数の電子顕微鏡装置各間の機差を機差情報として機差データベースに格納するステップと、前記複数の電子顕微鏡装置の少なくとも一つを用いて得られた前記寸法計測対象試料のパターンの寸法計測結果を、前記機差データベースに格納された機差情報を用いて補正して、前記寸法計測対象試料のパターンの寸法を計測するステップと、を有することを特徴とするパターン寸法計測方法である。
（３）互いにネットワークで接続された複数の電子顕微鏡装置と、前記複数の電子顕微鏡装置各々の装置パラメタを格納する装置パラメタデータベースと、前記複数の電子顕微鏡装置の一つである第一の電子顕微鏡装置を用いて寸法測定対象試料を走査することにより取得した寸法測定対象試料のパターンの画像の実信号波形と、前記装置パラメタデータベースに格納された前記第一の電子顕微鏡装置の装置特性パラメタを用いてシミュレーションにより作成したシミュレーション波形ライブラリと、を比較して、前記寸法測定対象試料のパターンの推定断面形状を得、前記得られた寸法測定対象試料のパターンの推定断面形状と、前記複数の電子顕微鏡装置のうち前記第一の電子顕微鏡以外の複数の電子顕微鏡装置各々の装置特性パラメタと、を用いて、前記第一の電子顕微鏡以外の複数の電子顕微鏡装置各々に対応する複数の模擬信号波形を生成するコンピュータと、前記コンピュータにより生成された前記複数の模擬信号波形各々に基づいて算出した複数の前記寸法測定対象試料のパターンの寸法と前記実信号波形に基づいて算出した前記寸法測定対象試料のパターンの寸法とを用いて算出された、前記第一の電子顕微鏡を含めた複数の電子顕微鏡装置各間の機差を機差情報として格納する機差データベースと、前記機差データベースに格納された機差情報を用いて寸法計測結果を補正する処理を含むレシピが格納されたレシピデータベースと、を有することを特徴とするパターン寸法計測システムである。
（４）装置特性評価用試料を電子顕微鏡装置により撮像して得たパターン画像を用いて装置特性パラメタを定期又は不定期に算出するステップと、前記計測された装置特性パラメタを算出した日時とリンクさせてデータベースに格納するステップと、前記装置特性評価用試料のパターン画像の実信号波形と、前記装置特性パラメタを用いて作成されたシミュレーション波形ライブラリとを比較して、前記装置特性評価用試料のパターンの推定断面形状を得るステップと、前記装置特性パラメタの推定断面形状と前記装置特性パラメタを用いて模擬信号波形を生成するステップと、前記生成された模擬信号波形を用いて前記装置特性評価用試料のパターンの寸法計測を行い、当該寸法計測結果を格納するステップと、前記格納された寸法計測結果を前記装置特性パラメタの経時変化と合わせて表示するステップと、を有することを特徴とする電子顕微鏡装置の経時変化のモニタ方法である。
（５）（ａ１）複数の電子顕微鏡装置と，（ｂ１）各電子顕微鏡装置の装置特性パラメタを計測するシステムと，（ｃ１）各電子顕微鏡装置の装置特性パラメタを格納するデータベースと，（ｄ１）サンプルの断面形状を推定するシステムと，（ｅ１）該試料の推定断面形状と各電子顕微鏡装置の装置特性パラメタから各電子顕微鏡装置の模擬ＳＥＭ信号波形を生成するシステムと，（ｆ１）該模擬ＳＥＭ信号波形の寸法計測を行って，模擬ＳＥＭ信号波形間の寸法計測結果の差を格納するデータベースを有し，上記，複数の電子顕微鏡装置での寸法計測結果を，上記模擬ＳＥＭ信号波形間の寸法計測結果の差を格納するデータベースを参照して補正することを特徴とする機差補正方法である。
（６）（ａ２）電子顕微鏡装置と，（ｂ２）電子顕微鏡装置の装置特性パラメタを計測するシステムと，（ｃ２）電子顕微鏡装置の装置特性パラメタを上記装置特性パラメタを算出した日時とリンクさせて格納するデータベースと，（ｄ２）サンプルの断面形状を推定するシステムと，（ｅ２）該試料の推定断面形状と上記電子顕微鏡装置の装置特性パラメタから模擬ＳＥＭ信号波形を生成するシステムと，（ｆ２）該模擬ＳＥＭ信号波形の寸法計測を行い，寸法計測を格納するデータベースを有し，上記，（ｂ２），（ｃ２），（ｄ２），（ｅ２），（ｆ２）の各ステップの的体に行って，装置特性パラメタの経時変化と寸法計測結果を合わせて提示することを特徴とする電子顕微鏡装置の経時変化のモニタ方法である。
【発明の効果】
【００２３】
本発明によれば，型式が異なる等の理由で，測長SEMのハードウエハの調整では機差をなくすことが困難であり，また，機差実測方式では対応が困難な，計測対象が多岐に及ぶようなケースにも適用が可能な機差補正方法、パターン寸法計測方法、パターン寸法計測システムおよび電子顕微鏡装置の経時変化のモニタ方法を提供することが可能である。
【図面の簡単な説明】
【００２４】
【図１】全体の処理フロー
【図２】測長ＳＥＭの構成
【図３】測長ＳＥＭにおけるエッジ点検出方法の説明
【図４】機差の測定法の一つであるグレーティング法の説明図
【図５】種々の計測対象を示す図
【図６】側壁傾斜角の違いと機差の関係を説明する図
【図７】パターン密度の違いと機差の関係を説明する図
【図８】ラインパターンとスペースパターンの機差の違いを説明する図
【図９】２次電子の広がりを示す図
【図１０】ピラミッドサンプルを用いたビームチルト角計測の方法を示す図
【図１１】ビーム径と明るさ変換係数を計測する第１の方法を示す図
【図１２】ビーム径と明るさ変換係数を計測する第２の方法を示す図
【図１３】モデルベース計測法の原理図
【図１４】モデルベース計測法による断面形状推定の方法を示す図
【図１５】機差を予測する処理のフロー
【図１６】システム構成
【図１７】レシピ作成時のＧＵＩ
【図１８】ビームチルト角の違いと機差の関係を説明する図
【図１９】ホール計測への適用も可能であることを示す説明図
【図２０】装置の経時変化モニタへの適用を示す説明図
【発明を実施するための形態】
【００２５】
[第１の実施の形態：機差補正]
はじめに図１を参照して，本実施の形態に係る機差に対応したパターン寸法計測の全体フローを説明する。
本実施の形態は，装置特性モデル化のステップ，機差予測のステップ，測長のステップの３ステップを適宜有する。以下，各ステップについて説明する。
【００２６】
[装置特性モデル化のステップ]
本ステップでは装置特性評価用試料１０１にて，測長SEM２０１−１〜２０１−ｎの装置特性パラメタを，装置特性パラメタ算出部３０１を用いて算出し，この結果を装置特性パラメタデータベース４０１に保存する。
本実施の形態では機差補正を目的としているので，モデル化に望ましいのは，機差へのインパクトが大きい装置特性である。本実施の形態では、望ましい装置特性として、照射ビームのチルト角，照射ビームのビーム径，試料から放出された２次電子の検出効率の装置間差を対象とする。
【００２７】
このうち、照射ビームのビーム径の違いについては，これが機差を引き起こすこと，かつ，その機差が断面形状，パターン密度，さらに寸法にも依存して変化することは図６〜図８で説明した通りである。
【００２８】
次に、照射ビームのチルト角と機差の関係について、図１８（Ａ）〜（Ｃ）を用いて説明する。側壁傾斜角がθの台形の断面形状を有する試料１３０に対して，垂直上方からビームが照射された場合の信号波形を模式的に表したのが図１８（Ａ）の１５０である（図１０のように２次電子強度はエッジ部ほど大きい）。この信号波形に対して，図３に示すしきい値を適用すると，寸法計測値はＣＤaとなる。一方，図１８（Ｂ）の１５１はビームチルト角がφ１の場合の信号波形である。この信号波形に対して，同じくしきい値法を適用すると，寸法計測値はＣＤｂとなる。また、図１８（Ｃ）の１５２はビームチルト角がさらに大きいφ２の場合である。この時の寸法計測値はＣＤｃとなる。図から明らかなように，ＣＤａ＝ＣＤｂ，ＣＤａ＜ＣＤｃの関係にある。ビーム照射角がチルトしたとしても，そのチルト角が側壁傾斜角よりも小さい場合には，寸法計測値は変化しないが，チルト角が側壁傾斜角を上回ると寸法計測値が変化するという関係である。これは，装置間にビーム照射角の違いがある場合，その照射角の違いの程度と，計測対象の断面形状との関係によって，機差が生じる場合も，生じない場合もあるということを意味している。以上のように，ビームチルト角の装置間差は機差の要因となり得、照射ビーム径の装置間差と同様、その機差は計測対象の断面形状やパターン密度，さらに寸法にも依存して変化する（図１８は模式図につき，パターン密度や寸法に対する依存性は表れてこないが，実際はパターン密度や寸法によって信号波形は変化するので，チルト角変化のインパクトがこれらに応じて変化するのは明らかである）。
【００２９】
ここで，図１０を参照して，照射ビームのチルト角の計測について述べる。照射ビームのチルト角の計測は，例えば、特開２００５−１８３３６９に記述された方法が好適である。この方法では，結晶異方性エッチングの技術により作製した凹または凸のピラミッド形状（四角錐）が配列されたサンプルを用いる。図１０の１４０に凹型のピラミッド形状の例を示す。結晶異方性エッチング技術とは，単結晶に特定のエッチング溶液を作用させた場合に結晶面ごとにエッチング速度が異なることを利用して，結晶面を基準とした，三角山や段差構造を製作する技術である。結晶面に相当する側面（図１０においては４個の側面）の互いになす角が決まっている。通常の半導体プロセスによるパターニングに比べて形状精度が高いのが特徴である。ピラミッド形状のＳＥＭ像を１４１に示す。平坦部における２次電子強度よりも傾斜部における２次電子強度の方が大きいため，ＳＥＭ像上ではピラミッドの側壁部が平坦部よりも明るい画像となる。１４１はビーム照射角がチルトしていない場合であるが，チルトしている場合には，例えば、１４２のような状況となる。１４３のように，画像処理により稜線，及び，ピラミッド底面のエッジ線を検出し，ピラミッド稜線の交点のｘ方向のずれΔｘ，ｙ方向のずれΔｙを求める。Δｘはｘ方向のビームチルト角の変化φxに，Δｙはｙ方向のビームチルト角の変化φyを反映する。先に述べたようにピラミッド側面がなす角は既知なので，Δｘ，Δｙからビームチルト角φx，φyを算出することができる。
【００３０】
図１に戻ると，装置特性評価用試料として，例えば、上記のピラミッドサンプルを用い，ピラミッドサンプルの画像を各測長ＳＥＭ２０１−１〜２０１−ｎで取得し，装置特性パラメタ算出部３０１にて各測長ＳＥＭのビームチルト角を算出し，この結果を装置特性パラメタデータベース４０１に保存するということになる。
【００３１】
ピラミッドサンプルを用いる利点は，４個の側面のなす角が既知なので，所定の計算式を適用すればビーム傾斜角が計算できることにあるが，もう一つの利点として，この方法が，他の機差要因の影響を受けづらいことが挙げられる。測長ＳＥＭ間の照射ビームのビーム径の違いは画像上ではエッジ部のシャープさの違いとして現れるが，上記の方法はピラミッド形状の幾何学変形をもとにチルト角を計算しているので，エッジ部のシャープさの違いには殆ど影響されない。また，測長ＳＥＭ間の２次電子の検出効率の違いは画像上では明るさやＳ／Ｎの違いとして現れるが，この違いもチルト角の計測結果に影響を及ぼすことはない。
【００３２】
次に，図１１を参照して，照射ビームのビーム径，及び，２次電子の検出効率の計測方法について述べる。ここでは、図１１中に示した数式１７０のように，ビーム径がゼロの場合の信号波形Ｉｓｉｍ（ｘ）に対して，ビーム強度分布に相当するＧａｕｓｓ（ｘ）を畳み込み，ダイナミックレンジを調整する係数ｓｃａｌｅを乗じ，オフセットを調整する係数ｏｆｆｓｅｔを加算した結果が実際の信号波形Ｉｒｅａｌ（ｘ）波形になるとモデル化する。よって，装置特性パラメタとして決定すべきは，Ｇａｕｓｓ（ｘ），ｓｃａｌｅ，ｏｆｆｓｅｔである。Ｇａｕｓｓ（ｘ）については，ガウス関数のσ値が決定すべきパラメタである。
【００３３】
これらのパラメタ算出のために望ましい装置特性評価用試料としては，断面形状が既知のステップエッジ試料，例えば，米国ＶＬＳＩスタンダード社製のＮＣＤ線幅スタンダード等がある。図１１の１６１はステップエッジ、１６２はステップエッジの信号波形を示し、さらに１９０、１９１にステップエッジのＳＥＭ写真の例を示す。同試料は，シリコン酸化膜／シリコン／シリコン酸化膜からなる積層膜の小片をウェーハ上に垂直にマウントし，シリコン酸化膜部をエッチングすることで，シリコン部分をラインパターン状に残したもので，その作成原理上，垂直のエッジ（側壁傾斜角０°）を有するため、上記パラメタの算出に望ましい。
【００３４】
以下，図１１に示した各ステップ（Ｓ３１０〜Ｓ３１２）を順に説明する。
まず、ＳＥＭシミュレーションにより、ビーム径がゼロの場合のステップエッジにおける信号信号波形をＩｓｉｍ（ｘ）を計算する（ステップＳ３１０）。この際，図１１に示した方法で計測した照射ビームのチルト角の計測結果を反映する。すなわち，ＳＥＭシミュレーションにおいて，試料に対する電子の打ち込み角度をチルト角の計測結果に合わせる。
次に、ステップエッジの側長ＳＥＭ像１５１を撮像し，実信号波形Ｉｒｅａｌ（ｘ）を求める（ステップＳ３１１）。さらに、求めたＩｓｉｍ（ｘ）、Ｉｒｅａｌ（ｘ）を用いて、数式１７０を満たす、Ｇａｕｓｓ（ｘ），Ｓｃａｌｅ，Ｏｆｆｓｅｔを求める（ステップＳ３１２）。この際，これらを解析的に求めることはできないので，Ｇａｕｓｓ（ｘ），Ｓｃａｌｅ，Ｏｆｆｓｅｔを逐次変化させ，そのつど，右辺と左辺の乖離度（例えば差の二乗和）を求め，乖離度が所定の値以下になるまで，繰り返し演算を行うという方法をとる。例えば，ＬＭ法（レベンバーグ・マーカート法）などの非線形最適化手法が適用可能である。ステップＳ３１２の結果の一例を吹き出し内に示した。波形１８０は式１７０の左辺，波形１８１は決定されたＧａｕｓｓ（ｘ），Ｓｃａｌｅ，Ｏｆｆｓｅｔを代入した場合の右辺である。
【００３５】
図１にもどると，装置特性評価用試料として上記のＮＣＤ線幅スタンダードサンプルを用い，その画像を各測長ＳＥＭ２０１−１〜２０１−ｎで取得し，装置特性パラメタ算出部３０１にて各測長ＳＥＭごとに，Ｇａｕｓｓ（ｘ），Ｓｃａｌｅ，Ｏｆｆｓｅｔを算出し，この結果を装置特性パラメタデータベース４０１に保存するということになる。
【００３６】
ＮＣＤ線幅スタンダードサンプルを用いることの利点は，同試料が垂直のエッジ（側壁傾斜角０°）を有すことが保証されているので断面形状の実測が不要ということ，及び，垂直に近いエッジ形状の方が，テーパを有するエッジ形状よりも，ビーム径の変化（ここでは，ガウス関数のσ値の変化ということになる）に対する感度が大きいという２点が挙げられる。
【００３７】
ただし，パターンを実測するというステップを加えれば，ＮＣＤ線幅スタンダードサンプルという特殊な製法によるサンプルを用いずとも，Ｇａｕｓｓ（ｘ），Ｓｃａｌｅ，Ｏｆｆｓｅｔを求めることは可能である。この場合のパラメタ算出ステップを図１２に示す。図１２では，最初のステップＳ３２０において，試料の断面形状を計測する。Ｓ３２０の吹き出し内は，ｙ方向に伸びるラインパターンを原子間力顕微鏡で計測した結果を模式的に示したものである。原子間力顕微鏡でのスキャン数分の計測結果が得られるが，これらを平均して，平均的な断面形状を求める。次のステップＳ３２１では，先のS３２０で求めた断面形状を入力として，ＳＥＭシミュレーションを行い、求めた断面形状におけるビーム径がゼロの場合のシミュレーション波形Ｉｓｉｍ（ｘ）を計算する。以下、図１１で示したステップＳ３１１，Ｓ３１２と同様にして、ステップＳ３２２，Ｓ３２３を行い、Ｇａｕｓｓ（ｘ），Ｓｃａｌｅ，Ｏｆｆｓｅｔを求める。
【００３８】
[機差予測のステップ]
次に，図１を参照して，機差予測のステップについて述べる。本ステップでは，寸法計測対象である試料１０２を，ＳＥＭ＃１号機にて撮像し，ＭＢＬシステム３０２（詳細は後述）にて計測対象パターンの断面形状を推定し，ＳＥＭシミュレータ３０３にて，推定された断面形状と先に装置特性パラメタデータベース４０１に保存された各測長ＳＥＭ（＃１号機〜＃ｎ号機）の装置パラメタを用いて，各測長ＳＥＭの信号波形を生成する。そして，測長部３０４にて，各ＳＥＭの信号波形の寸法計測を行い、それらの差を機差算出部３０５で求め，この結果を，機差データベース４０２に保存する。
【００３９】
図１３と図１４にて，ＭＢＬシステムでの処理内容を説明する。ＭＢＬ法（モデル・ベース・ライブラリマッチング法）は，図１３に示すように，予め，ＳＥＭシミュレーションにより，種々の断面形状におけるＳＥＭ波形を計算しておき，これと実波形とのマッチングを行うことで，対象の断面形状を推定する方法である。なお，図１３は，文献「J. S. Villarrubia， A. E. Vladar， J. R. Lowney， and M. T. Postek， “Scanning electron microscope analog of scatterometry，” Proc. of the SPIE， Vol. 4689， pp. 304-312 (2002)」からの抜粋である。
【００４０】
図１３ではＳＥＭシミュレータに対して，形状バリエーションのみを入力しているが，本発明では，図１４に示すように，形状バリエーション３１０の他，装置特性パラメタ４０１，及び，近隣パターン情報３１６を入力する。前述のように，装置特性パラメタデータベースには，照射ビームのチルト角，および，Ｇａｕｓｓ（ｘ），ｓｃａｌｅ，ｏｆｆｓｅｔが保存されている。ここでは，ＳＥＭ＃１号機の装置パラメタをデータベースから取り出し，ＳＥＭシミュレーションにて信号波形を生成する際にこれらの値を代入する。なお，ＳＥＭシミュレータ３１１としては，図１１のＳ３１０あるいは，図１２のＳ３２１で使用するのと同じものを用いればよい。近隣パターン情報（例えばパターンピッチなど）を入力するのは，前述のように，パターン密度もＳＥＭ信号波形に影響を及ぼすからである。従って，ライブラリ３１２は，実際の装置状態と，近隣パターンの影響とが考慮されたものとなる。ライブラリ波形と，寸法計測対象のＳＥＭ像３１４から得られた実信号波形とをマッチングし（３１３），ライブラリ波形の選択、及び，エッジ位置を求めることで，寸法計測対象の断面形状３１５を推定する。
なお，図１４においては，説明を単純化するため，形状バリエーションとして側壁傾斜角θのみをパラメタとしたが，より正確にすべく、計測対象に応じて，その形状を表現するのに適したパラメタ（例えば，トップやボトムの丸まり，あるいは，側壁傾斜角が上方と下方とで異なるのであれば，上側側壁傾斜角と下側側壁傾斜角など）を適宜入力するようにしてもよい。また、近隣パターン情報としては寸法計測対象試料の設計データ等を適宜用いればよい。
【００４１】
従来のＭＢＬ法（図１３）と比べると，既に確定された装置パラメタを代入することで断面形状推定の安定性が増すと共に，近隣パターンの影響も加味されているので，より，正確な断面形状が推定可能という特長を有す。以上が，図１のＭＢＬシステム３０２における処理内容である。
【００４２】
次に，図１５を用いて，図１の測長部３０４，機差算出部３０５での処理内容を説明する。測長ＳＥＭ＃１号機（２０１−１）については、既に寸法計測対象パターンのＳＥＭ像３１４を撮像済みなので（図１の３０２のステップ，すなわち，ＭＢＬシステムにて断面形状を推定するステップにおいて用いるために取得済み），ＳＥＭ像３１４の実信号波形３３０に対して，しきい値法（図３参照）など所定のエッジ検出手法を適用して寸法３３１を算出する。測長ＳＥＭ＃２号機については，図示しない装置パラメタデータベースより，＃２号機の装置パラメタ４１０を取り出し，これと，先に推定した計測対象パターンの断面形状３１５をＳＥＭシミュレータ３０３への入力として，模擬信号波形３３２を生成し，これに対して同じくしきい値法などを算出して寸法３３３を算出する。そして、寸法３３１と寸法３３３の差を求め（３３４）この結果を機差データベース４０２に保存する。他の測長ＳＥＭ（３号機〜n号機）についても同様である。以上により，機差データベース４０２には、寸法計測対象パターンの予測機差が保存される。図１５では、側長ＳＥＭ＃１号機と他の側長ＳＥＭとの機差を求める例を示したが、これに限られず、各装置間の機差を全て求めて保存してもよい。
【００４３】
通常の測長ＳＥＭの運用では，レシピと呼ばれる撮像シーケンス（ウェーハアライメント，パターン位置認識，オートフォーカスなどの一連の動作）や，測長方法（ラインやホールと言った測長種，エッジ検出手法や，そのための諸条件）を記述したデータを予め作成しておき，本番の計測では，このデータを呼び出して自動で寸法計測を行うのが一般的である。機差を予測して，その結果を機差データベースに保存する処理は，レシピ作成の際に合わせて実施するのが望ましい。図１７は，コンピュータ３００の画面上に表示されるレシピ作成のためのＧＵＩの例である。側長種６０１、側長方式６０２、波形平滑化フィルタサイズ６０３、側長ボックスサイズ６０４の側長条件の項目の他に，機差補正の項目を設け、機差補正を適用する場合には，チェックマーク６１１をＯＮにし，６１２にて機差が記述されたデータベースのファイル名を入力する。さらに，６１０にて出力項目を指定する。例えば，機差補正前と機差補正後の寸法を両方出力する場合には，６１３と６１４の両方をＯＮにする。なお、ここで示したＧＵＩは一例に過ぎず、上記した機能を可能とする範囲において種々変更可能であることは言うまでもない。また、図示しないが、出力としては、前述の機差補正前後の寸法のほかに寸法計測対象試料のＳＥＭ画像等を同時に表示しても、画面を切り替えることにより表示するようにしてもよく、種々変更可能である。
【００４４】
チェック項目６１５〜６１７は，所望高さにおけるパターン寸法の計測を必要とする場合にＯＮとする項目である。本発明においては，ＭＢＬシステム（図１の３０２）にて計測対象パターンの断面形状を推定済みなので，推定断面形状から求めたボトム寸法，１／２高さ寸法，トップ寸法と，図１５の３３１の寸法計測結果との差を用いれば，所望の高さにおける寸法を出力することも可能である。すなわち，測長ＳＥＭ＃１号機の寸法計測結果とＭＢＬシステムにて推定した断面形状から求めたボトム寸法との差がΔbottom，測長ＳＥＭ＃ｉ号機での寸法計測結果がＣＤi，機差データベースに記述されたＳＥＭ＃１号機との機差がｄｉとすると，ＳＥＭ＃ｉ号機においてボトム寸法６１５をＯＮにした場合は，数式４にてボトム寸法が計算され，測長結果の一項目として出力される。
（数４）ボトム寸法＝ＣＤi−ｄｉ−Δbottom
１／２高さ寸法，トップ寸法についても同様である。なお，６１０は寸法計測対象パターンの断面形状を台形として扱う場合の項目であって，例えば，前述のように，上方の側壁傾斜角と下方の側壁傾斜角が異なる，二つの台形を積み重ねたような形状の場合は，６１０の項目をそれに即したものにすることは言うまでもない。
【００４５】
[測長のステップ]
最後に，図１を参照して，測長のステップについて述べる。各装置での寸法計測は，前述したようにレシピの記載内容に従って行われる。測長ＳＥＭ＃ｎ号機にて機差補正ありが選択されている場合（図１７の６１１がＯＮの場合)，測長ＳＥＭ＃ｎ号機２０１−ｎにて寸法計測対象試料の寸法計測が行われ（３０４−ｎ），その後，機差データベースよりＳＥＭ＃ｎ号機の機差が読み出しされ，寸法計測結果が補正されることで、最終的な寸法計測対象試料の寸法が計測される。３０４-ｎの段階では，寸法計測結果は機差を有するが，機差補正部３０６−ｎを経た後は，機差のない計測値となる。出力されるのは，図１７のＧＵＩにて指定された各項目等である。
【００４６】
[全体システム]
図１７を参照して、本発明に係る機差に対応したパターン寸法計測システムのシステム構成例を説明する。本システムは，ローカルエリアネットワーク５００に接続された複数の測長ＳＥＭ，２０１−１，２０１−２，・・・２０１−ｎ，測長ＳＥＭのレシピを格納するレシピデータベース５０１，測長ＳＥＭで撮像されたＳＥＭ像が格納されるデータベース４００，装置パラメタが格納されるデータベース４０１，機差が格納されるデータベース３０７，及び，コンピュータ３００を適宜用いて構成される。
【００４７】
コンピュータ３００は，図１に示した装置パラメタ算出（３０１），ＭＢＬシステムによる断面形状推定（３０２），ＳＥＭシミュレーションによるＳＥＭ信号波形生成（３０３）などの各処理を行う。これらは，全て共通のＳＥＭシミュレータを用いるものであり，また，予め測長ＳＥＭで取得したＳＥＭ像を画像データベース４００に格納しておけば，全て，オフラインで実行できるものであるため，本システム構成のように，専用の計算機としてコンピュータを設けることが望ましい。なお、入出力を行うコンピュータと、上記各処理を行う計算機とを個別に設けても構わない。
【００４８】
ここで，図５では計測対象としてライン状のパターンを示したが，本発明の対象は，ラインパターンに限定するものではない。図１９に示すようなホールパターンの計測や，それ以外の任意形状の計測についても同様の機差補正方法、機差に対応したパターン寸法計測方法が適用可能である。典型的なホール計測では，例えば、各方向（図１９では，０度，４５度，９０度，１３５度）について信号波形を求め，各信号波形に対してしきい値法などでエッジ検出を行って寸法を求め，全方向の寸法計測結果の平均値をホール径として出力する。このようなホール計測に対して本発明を適用する場合，図１４でライブラリを作成する際には，入力形状をホールとし、近隣パターン情報としてホール間隔などを入力してライブラリを作成し，波形マッチング３１３を行って，所定方向のホールの断面形状を求める。これを，各方向（図１９の計測方法の場合は０度，４５度，９０度，１３５度の４方向）について行う。そして，図１５においては，各方向の断面形状を入力としてＳＥＭシミュレーションを行って信号波形を求めて寸法計測を行い，全方向の平均値をホール径とすれば良い。その他、任意の形状でも、対象パターンを適宜分割するなどして同様の処理を行えば、寸法計測可能である。
【００４９】
[第２の実施の形態：装置の経時変化モニタ]
第１の実施の形態においては，複数台の測長ＳＥＭを対象とした機差補正方法、機差に対応したパターン寸法計測方法について述べたが，本発明は，各種装置の経時変化のモニタにも適用可能である。以下、第２の実施の形態で行う処理の詳細は、第１の実施形態で説明した処理と同様であるため、異なる点を主として説明する。
装置の経時変化のモニタのためには，装置特性パラメタの算出（図１の３０１のステップの実施）を定期的に行い，ＳＥＭシミュレータにその時々の装置パラメタを入力して，ＳＥＭ信号波形を生成して寸法計測を行う（図１の３０３及び３０４のステップを実施する）。この際、定期的に算出した装置特性パラメタは算出した日時とリンクさせてデータベースに格納しておく。図２０に示すように，横軸に日時を，縦軸に各装置パラメタと寸法計測結果を表示すれば，装置の経時変化の程度が確認できると共に，装置の経時変化の要因のヒントが得られる。例えば，ビームチルト角が変動していれば光軸のずれが，ビーム径が変化している場合は電子銃の劣化が，ｓｃａｌｅ値が変化している場合は，照射ビーム電流の変化が疑われる。なお、装置特性パラメタの算出は必ずしも定期的である必要はなく、不定期で実施するようにしても構わない。
【００５０】
装置特性パラメタの経時変化の監視だけでも意味はあるが，上記の方法によれば寸法計測結果へのインパクトを知ることができる。つまり，寸法計測結果に変化が現れていれば早急な対策が必要であるし，寸法計測結果に変化がなければ，早急な対策は不要といった判断を下すことが可能となる。
【００５１】
さらに，本実施の形態は，従来の経時変化モニタ手法と比較すると，精度の面でも有利である。従来，装置の経時変化は，背景の技術で述べたグレーティング法のような方法が用いられてきた（時間をさかのぼることはできないので，ＡＢＢＡ法は適用不可）。すなわち，試料上の多数箇所を計測して，その平均値を求め，その結果を過去の結果と比較することで，装置の経時変化をモニタしていた。しかし、グレーティング法は，前述のように、相当数の計測点数を用いない精度が得られないという問題を有する。加えて、表面に酸化膜が生じることで寸法が変化したり、その変化の程度が材質や試料の保存状態に依存することから、寸法計測結果が変化した場合，それが装置の経時変化によるのか，試料の経時変化によるのかの判別がつかないという致命的な問題がある。一方，本発明の場合，図１の装置特性モデル化のステップで用いる装置特性評価用試料１０１の経時変化の可能性は完全には否定できないものの，図１０に示したピラミッドサンプルにしても，図１１に示したステップエッジサンプルにしても，これらのサンプルを用いて抽出される装置特性パラメタに，表面酸化膜は殆ど影響しないと考えられる。ピラミッドサンプルの場合，既知情報として用いているのは，４つの側面がなす角のみなので，これが変化しなければビームチルト角の計測結果への影響はない。ステップエッジサンプルの場合も同様である。既知情報として用いているのは，側壁傾斜角が垂直ということだけなので，これが変化しなければ，装置特性パラメタ（Ｇａｕｓｓ（ｘ），ｓｃａｌｅ，ｏｆｆｓｅｔ）の算出結果に影響を及ぼすことはない。表面酸化膜の付着によって，ピラミッドサンプルの４つの側面のなす角が変化や，ステップエッジサンプルの側壁傾斜角の変化は考えづらい。そして，図２０で，出力する寸法計測結果は，装置状態が変化した場合の寸法計測結果の予測値なので，当然のことながら，装置の経時変化によるのか，試料の経時変化によるのかの判別がつかないという点は問題とならない。
【００５２】
以上、本発明者によってなされた発明を実施形態に基づき具体的に説明したが、本発明は上記実施形態に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能であることはいうまでもない。
例えば、上記の第１の実施の形態及び第２の実施の形態では，装置特性パラメタとして，ビームチルト角，ビーム径（Ｇａｕｓｓ（ｘ）），明るさ変換係数（ｓｃａｌｅ，ｏｆｆｓｅｔ）を計測対象としたが，これらに限られるものではなく、その他の寸法計測結果に影響を及ぼす装置特性パラメタが存在する場合には，それを計測対象に適宜用いても構わない。
【００５３】
また、本願において開示される発明のうち代表的なものによって得られる効果を簡単に説明すれば次の通りである。
本発明の第１の実施の形態によれば，装置特性を予めモデル化しておくことで，機差が予測され，その結果に基づいて機差補正が行われる。機差を実測して，その結果に基づき寸法計測値を補正するという従来の方法に比べると，大幅な手間と時間の削減可能である。この結果，型式が異なる測長SEMが混在している場合，例えば，旧型式の測長ＳＥＭとの混用も可能となるため，測長ＳＥＭの有効活用という観点でのメリットも大きい。
また，本発明の第２の実施の形態によれば，従来の経時変化モニタにおける最大の問題点である，試料の経時変化と装置の経時変化の判別ができないという問題が解決され，より正確な経時変化モニタが可能となる。
【符号の説明】
【００５４】
１０１・・・装置特性評価用資料
１０２・・・寸法計測対象試料
１２０〜１２３・・・パターン形状のバリエーション
１３０・・・パターン断面
１４０，１４１，１４２，１４３・・・ピラミッドサンプル
１５０，１５１，１５２・・・信号波形
１６１・・・ステップエッジサンプル
１６２，１８０，１８１・・・信号波形
１７０・・・モデル式
２０１−１〜２０１−ｎ・・・測長ＳＥＭ
３００・・・計算機
３０１，３０２，３０３，３０４，３０５，３０６・・・各処理部
３１０・・・形状バリエーション入力部
３１１・・・ＳＥＭシミュレータ
３１２・・・シミュレーション波形ライブラリ
３１３・・・波形マッチング処理部
３１５・・・推定断面形状
３１６・・・近隣パターン情報
３１４・・・ＳＥＭ像
３３０・・・実信号波形
３３２，３３５・・・模擬信号波形
３３１，３３３，３３６・・・寸法
３３４，３３７・・・機差算出部
４０１，４０２，３０７，５０１・・・データベース
４１０，４１１・・・装置パラメタ
５００・・・ＬＡＮ
６０１〜６１７・・・ＧＵＩのパーツ
Ｓ３１０〜Ｓ３１２・・・ビーム径と明るさ変換係数算出のステップ
Ｓ３２０〜Ｓ３２３・・・ビーム径と明るさ変換係数算出のステップ
【図１】
0005433522.tif 000002
【図２】
0005433522.tif 000003
【図３】
0005433522.tif 000004
【図４】
0005433522.tif 000005
【図５】
0005433522.tif 000006
【図６】
0005433522.tif 000007
【図７】
0005433522.tif 000008
【図８】
0005433522.tif 000009
【図９】
0005433522.tif 000010
【図１０】
0005433522.tif 000011
【図１１】
0005433522.tif 000012
【図１２】
0005433522.tif 000013
【図１３】
0005433522.tif 000014
【図１４】
0005433522.tif 000015
【図１５】
0005433522.tif 000016
【図１６】
0005433522.tif 000017
【図１７】
0005433522.tif 000018
【図１８】
0005433522.tif 000019
【図１９】
0005433522.tif 000020
【図２０】
0005433522.tif 000021