(57)【特許請求の範囲】
【請求項１】
パターンの形状，寸法，隣接パターンとの間隔をパラメタとしてＳＥＭシミュレーションを行い，該ＳＥＭシミュレーションにより得られたシミュレーション波形をライブラリに保存し、
計測対象パターンのＳＥＭ画像を取得し、
該取得したＳＥＭ画像を処理して計測対象パターンの特徴量を求め、
該求めた計測対象パターンの特徴量に基づいて前記ライブラリの中から前記計測対象パターンの特徴量に対応するシミュレーション波形のグループをマッチングに用いる波形として選択し、
該選択したマッチングに用いるシミュレーション波形のグループと計測対象パターンのＳＥＭ画像から得られる信号波形とのライブラリマッチングを行って最も一致したシミュレーション波形から前記計測対象パターンの形状情報を推定し、
該推定した計測対象パターンの形状情報に基づいて前記計測対象パターンの寸法を算出し、
前記推定した形状情報と算出したパターンの寸法の情報とを画面上にＳＥＭ画像と一緒に表示あるいは数値データとして出力する
ことを特徴とするＳＥＭ画像を用いたパターン寸法計測方法。
【請求項２】
前記ＳＥＭ画像を処理して求める計測対象パターンの特徴量は、前記計測対象パターンに隣接するパターンとの間隔又は前記計測対象パターンの寸法であることを特徴とする請求項１記載のＳＥＭ画像を用いたパターン寸法計測方法。
【請求項３】
パターンの形状，寸法，隣接パターンとの間隔をパラメタとしてＳＥＭシミュレーションを行い，該ＳＥＭシミュレーションにより得られたシミュレーション波形をライブラリに保存するステップと、
該シミュレーション結果から複数の画像特徴量を算出し，上記複数の形状入力パラメタと上記画像特徴量の関係を保存しておくステップと，
計測対象パターンのＳＥＭ画像を取得し、該取得したＳＥＭ画像を処理して計測対象パターンの特徴量を算出するステップと，
該算出された画像特徴量と上記形状入力パラメタと画像特徴量の関係を用いて計測対象パターンの形状パラメタを推定するステップと，
該推定された形状入力パラメタの情報を用いてライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めるステップと、
該シミュレーション条件により計測対象パターンの断面形状および寸法を推定するステップと
を有することを特徴とするＳＥＭ画像を用いたパターン寸法計測方法。
【請求項４】
前記ライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めるステップにおいて、前記推定した形状入力パラメタの情報を初期値として非線形最適化手法を用いて前記ライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めることを特徴とする請求項３記載のＳＥＭ画像を用いたパターン寸法計測方法。
【請求項５】
ＳＥＭ画像を用いたパターンの計測方法であって、
パターンの形状，寸法，隣接パターンとの間隔をパラメタとしてＳＥＭシミュレーションを行い、該ＳＥＭシミュレーションにより得られたシミュレーション波形をライブラリに保存するステップと、
前記計測対象パターンの断面概略形状を計測するステップと、
該計測された断面概略形状から形状入力パラメタを算出するステップと、
計測対象パターンのＳＥＭ画像を取得するステップと，
前記算出した形状入力パラメタの情報を用いてライブラリデータの中から前記取得したＳＥＭ画像の信号波形と最も波形形状が一致するシミュレーション条件を求めるステップと、
該シミュレーション条件から計測対象パターンの断面形状および寸法を推定するステップと
を有することを特徴とするＳＥＭ画像を用いたパターン寸法計測方法。
【請求項６】
前記計測対象パターンの断面概略形状を計測するステップにおいて、前記断面概略形状は、ＡＦＭあるいはスキャッタロメトリを用いて求めた計測対象パターンの断面概略形状であることを特徴とする請求項５記載のＳＥＭ画像を用いたパターン寸法計測方法。
【請求項７】
前記ライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めるステップにおいて、前記算出した形状入力パラメタの情報を初期値として非線形最適化手法を用いて前記ライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めることを特徴とする請求項５記載のＳＥＭ画像を用いたパターン寸法計測方法。
【請求項８】
前記ＳＥＭシミュレーションを、予め計測対象パターンの概略形状を数値データでモデル化し，該パターン形状の計測範囲に対応した複数の異なる側壁形状，寸法，隣接パターン間の距離の組み合せを形状入力パラメタとして実行することを特徴とする請求項１または３または５の何れかに記載のＳＥＭ画像を用いたパターン寸法計測方法。
【請求項９】
パターンの形状，寸法，隣接パターンとの間隔をパラメタとしてＳＥＭシミュレーションを行い，該ＳＥＭシミュレーションにより得られたシミュレーション波形をライブラリに保存するシミュレーション手段と、
ＳＥＭ装置で取得した計測対象パターンのＳＥＭ画像を処理して計測対象パターンの特徴量を求めるＳＥＭ画像処理手段と、
該ＳＥＭ画像処理手段で求めた計測対象パターンの特徴量に基づいて前記シミュレーション手段のライブラリの中から前記計測対象パターンの特徴量に対応するシミュレーション波形のグループをマッチングに用いる波形として選択するマッチング波形選択手段と、
該マッチング波形選択手段で選択したマッチングに用いるシミュレーション波形のグループと計測対象パターンのＳＥＭ画像から得られる信号波形とのライブラリマッチングを行って最も一致したシミュレーション波形から前記計測対象パターンの形状情報を推定するパターン形状情報推定手段と、
該パターン形状情報推定手段で推定した計測対象パターンの形状情報に基づいて前記計測対象パターンの寸法を算出するパターン寸法算出手段と、
前記パターン形状情報推定手段で推定した形状情報と前記パターン寸法算出手段で算出したパターンの寸法の情報とを画面上にＳＥＭ画像と一緒に表示あるいは数値データとして出力する出力手段と
を備えたことを特徴とするＳＥＭ画像を用いたパターン寸法計測装置。
【請求項１０】
前記ＳＥＭ画像処理手段でＳＥＭ画像を処理して求める計測対象パターンの特徴量は、前記計測対象パターンに隣接するパターンとの間隔又は前記計測対象パターンの寸法であることを特徴とする請求項９記載のＳＥＭ画像を用いたパターン寸法計測装置。
【請求項１１】
パターンの形状，寸法，隣接パターンとの間隔をパラメタとしてＳＥＭシミュレーションを行い，該ＳＥＭシミュレーションにより得られたシミュレーション波形をライブラリに保存するシミュレーション手段と、
該シミュレーション手段におけるシミュレーション結果から複数の画像特徴量を算出し，上記複数の形状入力パラメタと上記画像特徴量の関係を保存しておく画像特徴量算出手段と，
ＳＥＭ装置で取得した計測対象パターンのＳＥＭ画像における画像特徴量を算出し，該算出された画像特徴量と，上記形状入力パラメタと画像特徴量の関係を用いて，計測対象パターンの形状パラメタを推定するパターン形状パラメタ推定手段と，
該パターン形状パラメタ推定手段で推定された形状入力パラメタの情報を用いてライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めるシミュレーション条件決定手段と、
該シミュレーション条件決定手段で求めたシミュレーション条件により計測対象パターンの断面形状および寸法を推定するパターン断面形状・寸法推定手段と
を備えることを特徴とするＳＥＭ画像を用いたパターン寸法計測装置。
【請求項１２】
前記シミュレーション条件決定手段は、前記パターン形状パラメタ推定手段で推定した形状入力パラメタの情報を初期値として非線形最適化手法を用いて前記ライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めることを特徴とする請求項１１記載のＳＥＭ画像を用いたパターン寸法計測装置。
【請求項１３】
ＳＥＭ画像を用いたパターンの計測装置であって、
パターンの形状，寸法，隣接パターンとの間隔をパラメタとしてＳＥＭシミュレーションを行い、該ＳＥＭシミュレーションにより得られたシミュレーション波形をライブラリに保存するシミュレーション手段と、
前記計測対象パターンの断面概略形状情報を入力する断面概略形状情報入力手段と、
該断面概略形状情報入力手段に入力された断面概略形状情報から形状入力パラメタを算出する形状入力パラメタ算出手段と、
ＳＥＭ装置で取得した計測対象パターンのＳＥＭ画像を受けて前記形状入力パラメタ算出手段で算出した形状入力パラメタの情報を用いて前記シミュレーション手段のライブラリデータの中から前記計測対象パターンのＳＥＭ画像の信号波形と最も波形形状が一致するシミュレーション条件を求めるシミュレーション条件決定手段と、
該シミュレーション条件決定手段で求めたシミュレーション条件から計測対象パターンの断面形状および寸法を推定するパターン断面形状・寸法推定手段と
を備えたことを特徴とするＳＥＭ画像を用いたパターン寸法計測装置。
【請求項１４】
前記シミュレーション条件決定手段は、前記形状入力パラメタ算出手段で算出した形状入力パラメタの情報を初期値として非線形最適化手法を用いて前記ライブラリデータの中からＳＥＭ画像と最も波形形状が一致するシミュレーション条件を求めることを特徴とする請求項１３記載のＳＥＭ画像を用いたパターン寸法計測装置。
【請求項１５】
前記断面概略形状情報入力手段は、前記断面概略形状として、ＡＦＭあるいはスキャッタロメトリを用いて求めた計測対象パターンの断面概略形状を入力することを特徴とする請求項１３記載のＳＥＭ画像を用いたパターン寸法計測装置。
【請求項１６】
前記シミュレーション手段は、前記ＳＥＭシミュレーションを、予め計測対象パターンの概略形状を数値データでモデル化し，該パターン形状の計測範囲に対応した複数の異なる側壁形状，寸法，隣接パターン間の距離の組み合せを形状入力パラメタとして実行することを特徴とする請求項９または１１または１３の何れかに記載のＳＥＭ画像を用いたパターン寸法計測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は半導体の製造工程においてウェハ上に形成された回路パターンの加工形状の良否を，回路パターンの電子顕微鏡画像を用いて評価する方法及びそのシステムに関する。
【背景技術】
【０００２】
半導体ウェハの製造工程において，ウェハ上に多層に形成されるパターンの微細化が急速に進んでおり，それらのパターンが設計どおりにウェハ上に形成されているか否かを監視するプロセスモニタの重要性はますます増加している。中でもトランジスタゲート配線をはじめとする配線パターンは，その配線幅とデバイス動作特性に強い関連があり，その配線製造プロセスのモニタは特に重要である。
【０００３】
数十ナノメートルオーダの微細配線の線幅を測定する測長ツールとしては，それらの配線を10万から20万倍の拡大倍率で撮像可能な線幅測定用の走査型電子顕微鏡（測長ＳＥＭ(Scanning Electron Microscope)またはＣＤ(Critical dimension)ＳＥＭ）が従来から用いられている。このような走査型電子顕微鏡を用いた測長処理の一例が特開平１１−３１６１１５号公報（特許文献１）に記載されている。特許文献１の開示例においては，測定対象配線を撮像した画像内の局所領域から，配線の信号プロファイルを配線の長手方向に加算平均した投影プロファイルを作成し，このプロファイルにおいて検出した左右の配線エッジ間の距離として配線寸法を算出している。
【０００４】
しかし，非特許文献１のFig1に開示されているように，ＳＥＭの信号波形においては，その計測対象の形状が変化すると，それに応じて信号波形も変化するため計測誤差を生じるという問題がある。半導体パターンの微細化に伴い，これら計測誤差がプロセスモニタに与える影響は大きくなりつつある。このような計測誤差を低減する手法が非特許文献１および非特許文献２に開示されている。この手法では，パターン形状とＳＥＭ信号波形の関係を予めシミュレーションにより計算しておき，その結果を利用して対象形状に依存しない高精度計測を実現している。
【０００５】
また、パターン材料とパターンの下地材料とで二次電子放出効率が異なるようなサンプルを計測する場合に，事前に材料のパラメタを調整することで，ライブラリマッチングによる計測精度を向上させることが，非特許文献３に開示されている。
【０００６】
【特許文献１】特開平１１−３１６１１５号公報
【特許文献２】特開Ｐ２００７−２１８７１１号公報
【非特許文献１】J. S. Villarrubia， A. E. Vladar， J. R. Lowney， and M. T. Postek， “Scanning electron microscope analog of scatterometry，”Proc. SPIE 4689， pp. 304-312 (2002)
【非特許文献２】J. S. Villarrubia， A. E. Vladar， M. T. Postek， “A simulation study of repeatability and bias in the CD-SEM，” Proc. SPIE  5038， pp. 138-149， 2003.
【非特許文献３】M. Tanaka， J. S. Villarrubia and A. E. Vladar， “Influence of Focus Variation on Linewidth Measurements，” Proc. SPIE 5752， pp.144-155(2005)
【発明の開示】
【発明が解決しようとする課題】
【０００７】
背景技術に示されているように，測長ＳＥＭにより半導体パターンの寸法計測を行う際には，対象パターンの形状に依存した計測誤差を発生するという問題がある。これに対して，非特許文献１および非特許文献２に開示されている手法では，パターン形状とＳＥＭ信号波形の関係を予めＳＥＭシミュレーションにより計算しておき，その結果を利用して対象形状に依存しない高精度計測を実現している。パターン形状をパラメタにより数値化し，様々な形状のＳＥＭシミュレーション結果をライブラリとして保存し，実波形と比較することで，正確に形状と寸法を推定することが可能である。本明細書では，以後，本手法をモデルベース計測あるいはライブラリマッチング手法と呼ぶ。このようなモデルベース計測手法においては，いかに精度の良いシミュレーションを行うかが安定で高精度な計測を実現するために重要となる。
【０００８】
実画像とのマッチングに適したシミュレーション結果を得るためには，適切なパターン形状モデルの設定とシミュレーションパラメタの設定が重要である。しかしながら，最適なパターン形状モデルは，計測対象によって異なるが，これらを容易かつ適切に設定することは困難である。さらに，計測対象の材料特性などの物理パラメタを正確に求めるのは困難であるといった課題がある。
【０００９】
また，ライブラリとのマッチング処理においては，適切な形状モデルを用いていても，多くのパラメタの適切な組み合せを求めるのには多くの時間が必要になるという問題がある。多くの時間を要するパラメタ空間の全探索の代わりに，各種の非線形最適化手法を用いることができるが，この場合には，局所解におちいってしまい，正しい計測結果が得られなくなる。
【課題を解決するための手段】
【００１０】
本発明では，これらの課題に対して，計測対象であるサンプルの材料や形状のモデルおよびパラメタを，予め別の手段で推定した結果に基づいて推定することにより，非特許文献１および非特許文献２などに示されているようなシミュレーションを用いた計測手法の安定化，高速化，高精度化をはかる。
【００１１】
本発明では，寸法計測結果に与える影響が大きい隣接パターンとの距離を，実際のＳＥＭ画像を用いて従来法での計測結果に基づいて推定し，推定された隣接パターンとの距離の条件下でシミュレーションされた波形のみを用いてライブラリマッチングを行なうことで，高精度な計測を実現する。
【００１２】
また，パターン線幅についても，実際のＳＥＭ画像を用いて従来法での計測結果に基づいた推定を行い，推定された線幅の条件下でシミュレーションされた波形のみを用いてライブラリマッチングを行なうことで，高精度な計測を実現する。 
さらに，予めシミュレーションにより求めておいたパターン形状とＳＥＭ画像特徴量の関係に基づいて，計測時には，対象パターンのＳＥＭ画像から得られる画像特徴量を用いて，パターン形状を推定し，その推定結果をライブラリマッチングにおける初期値として用いることで，マッチング処理を高速化・安定化して，高精度かつ高速な計測を実現する。
【００１３】
さらに別の方法として，予めＡＦＭなどの別の計測装置で得られた情報に基づいて，ライブラリマッチングにおけるパターン形状パラメタの初期値を適切に設定することで，マッチング処理を高速化・安定化して，高精度かつ高速な計測を実現する。
【００１４】
また，形状モデルの最適化の課題に対しては，ＡＦＭや断面ＳＥＭなどにより得られる断面形状情報に基づいて，オペレータが容易に最適なモデルの選択とパラメタ設定が可能なグラフィカルユーザインターフェースを提供することにより，その結果，安定で高精度な計測を実現する。
【００１５】
複数種類の材料から構成されるパターンを計測する場合には，事前にＳＥＭ画像を取得し，シミュレーションによる各材料の信号量が，上記実ＳＥＭ画像のコントラストと一致するように，ライブラリ構築時の材料パラメタを設定することで，シミュレーションの精度を向上し，安定で高精度な計測を実現する。
【発明の効果】
【００１６】
本発明により，モデルベース計測手法に用いるシミュレーションの精度が向上でき，その結果モデルベース計測手法そのものの精度も向上する。また，材料パラメタや一部の形状パラメタを事前に設定することにより，推定パラメタの数を減らすことができるため，安定な推定が可能となり，また計測に要する計算時間を短縮することができる。さらに，推定が必要なパラメタについても，適切な初期値を用いることにより，安定で高速な推定が可能となり，ひいては計測の信頼性および速度を向上することが可能となる。
【発明を実施するための最良の形態】
【００１７】
本発明は，各種の荷電粒子線装置（ＳＥＭ，ＦＩＢ等）に適用可能であるが，以下の実施例では代表としてＳＥＭを用いる場合を例にとって説明する。
【実施例１】
【００１８】
第１の実施例では，計測対象パターン周辺のスペース，すなわち隣接パターンとの距離の違いによって生じるＳＥＭ計測の誤差を低減する方法について図１から図５を用いて説明する。本発明では予め従来手法により計測されたスペースの暫定計測結果に基づいて，ライブラリ波形マッチングにおいて用いるシミュレーション波形を，実際のパターンに合ったものに限定することで，マッチングの精度，ひいてはパターン計測精度を向上する。
【００１９】
図２は，従来手法の計測課題であるスペースに依存して発生する計測誤差を説明する図である。図２（ａ）のグラフに示すように，パターン２０１とパターン２０２とのスペース幅Ｓによってエッジ部のＳＥＭ信号波形Ｗは右側に拡大して表示したように、大きく変化する。この現象は，スペース部から発生する２次電子が側壁に衝突するため，検出信号量が減少することによるもので，パターンが高くスペースが狭いほど顕著である。図２（ｂ）はスペース幅Ｓを変化させたときに，従来法であるしきい値法（しきい値５０％）による計測値をシミュレーションにより計算したものである。
【００２０】
しきい値法とは，下地部分の信号量とエッジ部ピーク部分の信号量の間の任意の明るさを，下地部０，ピーク部１００％としてしきい値で指定し，この指定された閾値に相当する信号量をとる位置をパターンエッジ位置とする方法である。図２（ｂ）では，入力パターンの幅寸法Ｗ及び高さ寸法Ｈとは全て一定だが，計算結果はパターン２０１とパターン２０２との間のスペース幅Ｓによって変化している。誤差の範囲は計測対象パターンの寸法ＷとＨに依存するが，数ナノメートルから十数ナノメートルにおよぶ恐れがある。
【００２１】
一方，実際の計測対象パターンにおいては，図２（ｃ）の矢印示すように，そのレイアウトによって，隣接パターンとの距離は様々である。図２(ｃ)に示した例においては、パターン２１０，２１１，２１２のそれぞれの間の距離、Ｓ１〜Ｓ４に相当する。このことから，計測の高精度化にはスペースの影響を考慮することが重要であることがわかる。
【００２２】
この課題に対する解決手段として，図１に，本発明に係るパターン計測方法の第１の実施例のＳＥＭによる計測手順を示す。
【００２３】
計測開始前に，パターンの形状，寸法，隣接パターンとの距離を様々な値に設定したＳＥＭシミュレーションを行い，ＳＥＭシミュレーション波形をライブラリに保存しておく。シミュレーション条件の設定範囲は，計測対象パターンの製造プロセスに応じて適切な値に設定しておく。
【００２４】
計測時には，まず図１（ａ）に示したＳＥＭ００１にて，予め設定された撮像条件（倍率，照射ビームの加速電圧など）で試料のＳＥＭ像を取得する。具体的には、ＳＥＭ００１の電子銃１０１から発射された電子線１０２を集束レンズ１０３で収束させ、偏向器１０４でＸ方向及びＹ方向（図１において図面に垂直な平面内）に走査し、対物レンズ１０５で電子線の焦点を計測対象パターンが形成された試料１０６の表面に合わせて試料１０６の表面を走査して照射する。図１（ａ）では図示を省略してあるが、試料１０６はテーブルに載置されて平面内で移動可能になっており、試料１０６の表面の所望の領域が電子線１０２の照射領域に位置するように制御される。電子線１０２が照射された試料１０６の表面から発生した二次電子の一部は検出器１０７で検出され、電気信号に変換されて全体制御・画像処理部１０８に送られてＳＥＭ画像が作成され、演算部１０９でＳＥＭ画像を処理してパターンの寸法を算出し、結果が出力部１１０の画面上に表示される。全体制御・画像処理部１０８は、図示していない試料１０６を載置するテーブルを含めたＳＥＭ００１全体の制御も行う。
【００２５】
演算部１０９における処理手順を図１（ｂ）に示す。まず、上記に説明したように全体制御・画像処理部１０８でＳＥＭ００１を制御して計測対象パターンのＳＥＭ画像を取得する（Ｓ０００１）。次に，全体制御・画像処理部１０８で取得したＳＥＭ画像を受けて、演算部１０９でこのＳＥＭ画像を処理して，計測対象パターンの寸法および計測対象パターンと隣接パターンとの距離を暫定的に計測する（Ｓ０００２）。ステップ（Ｓ０００２）の計測はしきい値法など，従来のエッジ検出方法を用いて行えばよい。
【００２６】
次に，ステップ（Ｓ０００２）で計測された隣接パターンとの距離に基づいて，予め計算しておいたＳＥＭシミュレーション波形ライブラリの中から，ステップ（Ｓ０００２）で計測した実際のスペースに対応するシミュレーション波形をマッチングに用いる波形として選択する（Ｓ０００３）。
【００２７】
次に，ステップ（Ｓ０００３）で選択されたライブラリ波形のみを用いて，計測対象パターンの実波形とライブラリマッチングを行い，最も一致したシミュレーション波形の入力形状からパターン側壁エッジの位置と形状を推定する（Ｓ０００４）。マッチングにより得られた計測対象パターンのエッジ位置および形状推定結果に基づき，ユーザが予め指定した高さに相当するパターンの寸法を算出する（Ｓ０００５）。最後に得られた結果をＳＥＭ画像上への表示あるいは数値データとして出力部１１０の画面上に表示する（Ｓ０００６）。また、出力部１１０からの出力は、図示していない他のデータ処理装置又は記憶装置へ送信することも出来る。
以下，図１（ｂ）の各ステップの詳細について説明する。
（処理領域設定と従来法による初期計測）
ステップ（Ｓ０００２）での隣接パターンとの距離計測の詳細について図３を用いて説明する。隣接パターンとの距離計測は，演算部１０９の画像処理ユニット１０９１で行われ、予め設定された領域内の画像に対して処理を行う。図３（ａ）は計測対象のＳＥＭ画像００３に処理領域設定ウィンドウを設定する例である。画像中央の計測対象パターン００４に設定された寸法計測処理領域００５と，左右の隣接パターン００４１および００４２のエッジ間にまたがる領域に設定されたスペース計測用処理領域００６および００７の設定例を示している。これらの処理領域設定ウィンドウ００５〜００７の設定は，実際にＳＥＭ画像を見ながらオペレータが設定してもよいし，パターンの設計情報を用いれば自動的に設定することも可能である。これらのウィンドウは，計測レシピとして一度設定しておけば，他のウェハやショット上のサンプルに対しても，パターンの位置決めを行うことで自動的に設定することが可能である。
【００２８】
図３（ｂ）は計測ウィンドウに対応するｘ方向のＳＥＭ信号波形００８を示している。図３（ａ）のように設定された処理領域設定ウィンドウ内００５，００６，００７（図３（ｂ）では００５１，００６１，００７１が対応する領域）のＳＥＭ画像に対して処理を行うが，このとき，処理画像のＳ／Ｎを向上するための前処理として，複数のＳＥＭ信号波形を平均して用いてもよい。例えば，図３（ａ）のようなラインパターンであれば，局所的にはラインの長手方向の寸法やエッジ形状はあまり変わらないため，画像内の異なるｙ座標の複数の画素の画像データを平均して得られる波形を処理対象とすればよい。あるいは，非特許文献３に開示されているように，パターン寸法のばらつきに対応した平均処理を行ってももちろんよい。
【００２９】
このように，ＳＥＭ画像の処理領域が決まると，ステップＳ００２ではこれらの領域内において，従来法でのエッジ検出を行い，計測対象パターンの寸法および隣接パターンとの距離の暫定計測を行う。一例として，図３（ｂ）の波形に，しきい値法をしきい値５０％で適用したエッジ位置抽出結果を黒点で示している。このようにして抽出されたエッジ位置間の距離を計算することにより，図３（ｂ）のグラフに矢印で示すように，計測対象パターンの寸法および隣接パターンとの距離を暫定的に求めることができる。
【００３０】
ここでは，エッジ抽出手法としてしきい値法を用いたが，他のエッジ位置決定方法として，処理時間が十分に早いものであれば，信号量ピーク位置（しきい値１００％と等価），信号量の微分値の最大・最小値や，側壁部分のデータに直線を当てはめて下地明るさと交差する点など，他の一般的に用いられている計測手法を用いてももちろんよい。
【００３１】
また，図３（ａ）の例では左右に隣接パターンがある場合の例を示したが，ＳＥＭ画像の範囲内に隣接パターンが観測されない場合には，計測は行わず，代わりに設計値あるいは予め設定した値を用いて後の処理を行う。
【００３２】
図２（ｂ）に示したように，スペースに依存した計測誤差は，スペースが小さい程影響が大きい。このため，ある程度大きなスペースであれば，実際のスペース幅が正確にわからなくても，後の処理への影響は少なく問題ない。事前に，計測精度に影響を与えるスペース幅を調べておき，スペース幅に応じて，設計値を使うなど切り替えればよい。
【００３３】
例えば，図２（ｂ）において，スペース変化に対する計測誤差変動の感度が高い領域では必ずスペース計測ができる視野の画像を取得して正確なスペース値を用いればよい。感度が低い範囲であれば設計値をもちいればよいし，誤差が飽和する範囲であれば，事前に設定した固定値のシミュレーション結果を用いればよい。なお，スペースが狭い場合に設計値を用いないのは，スペースの変化に対する誤差感度が高いこととスペースが設計値どおりにできるとは限らないためである。
【００３４】
（ライブラリの構築，限定）
次に，図４を用いて本発明で用いるライブラリの構築方法と，ステップＳ０００３により暫定的に計測されたスペース幅を用いて使用するライブラリを制限する方法について説明する。このライブラリの構築は演算部１０９のシミュレーションユニット１０９２で行われる。
【００３５】
図４（ａ）は，パターンのスペース幅（x軸）と，同パターンをスペースの暫定計測時と同じ従来法で計測した結果（y軸）の関係を示しており，図２（ｂ）のグラフから導出できる。本グラフはシミュレーションライブラリの波形データに対して従来の計測処理を行えば容易に作成できる。図４（ａ）の関係がわかっていれば，ステップＳ０００２で得られた暫定計測結果(S_meas)から実際のスペース幅（S_est）をある程度の精度で推定できる。このようにして，実際のスペース幅を推定した後，シミュレーションライブラリから推定されたスペース幅を持つパターンに限定して計測用の波形マッチングを行う。
【００３６】
図４（ｂ）を用いて，本発明のライブラリマッチングおよびライブラリ限定方法を説明する。図４（ｂ）は，本発明で用いるシミュレーション波形ライブラリ００２の例を示している。シミュレーション波形ライブラリ００２にはシミュレーションの入力断面形状００９とそれに対応したＳＥＭシミュレーション波形０１０が記録されている。
【００３７】
側壁傾斜角θの変化を正しく計測するライブラリマッチング計測を行なう場合には，図４（ｂ）に示すように複数の異なる傾斜角θの形状についてシミュレーションを実施しておく。この例では簡単のため側壁傾斜角θが３種類の場合を示しているが，実際にはプロセス変動で生じうるパターン形状をカバーする範囲を，計測したい精度に見合う細かさでシミュレーションしておく。シミュレーションライブラリは図１の計測とは別に事前に作成しておく。本発明では，このライブラリに図４（ｂ）の横軸に示すように，スペースＳを変化させたシミュレーションも実施しておく。
【００３８】
なお，図４（ｂ）の例ではわかりやすいように２つの隣接するパターンの図で示したが，実際のライブラリ作成においては，スペース両側のシミュレーション結果は左右反転させただけの同じものになるため，隣接するパターンのどちらか片方のスペース側を計算しておけばよい。マッチング時には，マッチング対象のパターン方向によって必要に応じてシミュレーション信号を左右反転させればよい。
【００３９】
本発明では，このような形状パラメタの組み合わせに対して，ＳＥＭシミュレーションを実施し，形状情報００９とシミュレーション波形０１０を対応付けてライブラリ００２に保存する。なお，図４（ｂ）の例では，簡単のためスペースと傾斜角のみの例を表しているが，このほかにパターンのトップコーナＲｔやボトムコーナＲｂのラウンディングなどを変更したデータがライブラリに含まれていてももちろんよい（その場合には３次元以上の多次元空間となる）。
【００４０】
本発明では，このライブラリ００２に対して，図４（ａ）の処理により推定された実際のスペース幅を適用してライブラリ００２を限定する。図４（ｂ）では，推定されたスペース幅S_estがS3と一致した場合の例で，点線で囲ったように，実際のスペース幅S_estと一致するスペースで計算されたシミュレーション波形のみを選択した部分ライブラリ０１１をマッチング対象として用いることで，スペースに依存したＳＥＭ波形の変化を考慮した高精度な計測が実現可能となる。
【００４１】
図４（ｂ）の例では，簡単のため，単純にシミュレーション条件を選択する例を示したが，実際には事前に準備できるシミュレーション条件が限られていてライブラリ内の形状パラメタは離散値であるため，完全に一致するスペースのシミュレーションデータがあるとは限らない。このため，実際には，シミュレーションライブラリから推定値S_estに最も近いスペース条件を選択するか，あるいはS_estにおけるＳＥＭ波形を近い条件の波形を補間することで推定して用いればよい。
【００４２】
シミュレーション波形の補間は波形の補間で行うことができる。例えば，スペース50nmのシミュレーション波形f_s50(x)とスペース100nmのシミュレーション波形f_s100(x)がある場合に，スペース80nmの波形f_s80(x)は線形補間によりf_s80(x) = f_s50(x) + (f_s100(x)- f_s50(x))*(80-50)/(100-50)により求めることができる。
【００４３】
ライブラリのデータによっては，２つ以上のスペースパラメタの組み合わせを用いた非線形補間を行ってももちろんよい。これを利用すれば，ステップＳ０００３で推定スペース値に対応するライブラリが無い場合には，前後のスペース値の波形を用いて，当該スペース値（S_est）に対応する波形ライブラリを上記補間処理により新たに作成して，代わりに用いればよい。
【００４４】
（ライブラリマッチング）
次に，ステップＳ０００３により限定されたシミュレーションライブラリと実際の計測対象パターンの画像の比較により，ライブラリ内で最も一致度の高い波形を選択する。このライブラリマッチングは演算部１０９のマッチング波形選択ユニット１０９３で行われる。パターンの寸法計測では，パターン両側のエッジ間の距離が計測対象により変化する（このエッジ間距離が計測対象である寸法に相当する）。このため，図３（ｂ）の計測処理ウィンドウ００５内のＳＥＭ波形全体にマッチング処理を行うと，エッジ間の距離も変化させたシミュレーション波形が必要となり，マッチングに多くの計算量が必要となる。
【００４５】
そこで，本発明の計測方法では，図４（ｂ）に示すように，１つのパターンの評価領域００５１におけるＳＥＭ信号波形００８２のうち処理対象エッジ周辺の処理領域０１２のＳＥＭ信号波形００８３を局所的に切り出した後，エッジ毎にライブラリとの波形マッチング処理Ｓ０００４を行う。波形マッチングは，評価対象ウィンドウ内の実ＳＥＭ画像Ｓ００８とエッジ部の波形が最も一致するシミュレーション波形を前記の限定されたライブラリの中から探索することで行う。このとき，シミュレーション波形のｘ座標も同時に推定することで，形状だけでなく，それぞれのエッジ位置も正しく求めることができる。
【００４６】
マッチングの際の波形の一致度としては，例えば，波形間の差の自乗和を用いて，これを最小にするようなシミュレーション波形を選択すればよい。あるいは，離散的なライブラリに対して図４（ｂ）でスペースを限定した処理と同様に，補間処理を行うことで，ライブラリにないパラメタについてマッチングを行うことも可能である。補間処理により連続関数として扱えれば，Levenberg- Marquardt法などの非線形最適化手法を用いてマッチングするパラメタを推定することも可能である。
【００４７】
以上の手続きにより，実波形と最も一致度が高いシミュレーション条件を求めることで，各エッジの側壁断面形状およびその位置を高精度に推定することができる。図４（ｂ）の例であれば，マッチング結果から，例えば撮像されたＳＥＭ画像の断面形状がスペースＳ３と側壁傾斜角θ1に近い形状であると推定され，ライブラリマッチング領域０１２内の部分波形とマッチしたシミュレーション波形の位置ずれ量から，ＳＥＭ画像内におけるエッジ位置を求めることができる。推定されたエッジ部の側壁形状が得られるため，任意の高さにおけるパターン寸法を算出することが可能となる。この処理は演算部１０９のパターン寸法算出ユニット１０９４で行われる。
【００４８】
予め計測したいサンプル高さ（トップ，ボトム，中央、ボトムからパターン高さの１０％の位置など）を設定しておけば，得られた断面形状から，対応する箇所でのエッジ位置および対向するエッジ位置の差分から寸法を求めることができる（ステップＳ０００５）。上記のように，マッチングは左右エッジで別々に行なうため，それぞれのエッジについて，適切なスペース幅でのライブラリ制限を行なってから計測を行なうことで，左右非対称なパターン形状においても正確な計測が可能となる。
【００４９】
これにより、パターンの任意の位置の寸法データとパターン断面形状の情報とを出力部１１０の画面上に表示してユーザに提供することが出来る。また通信手段を介してデータサーバ１１１に出力部１１０からの情報を蓄積しておくことも出来る。
【００５０】
以上のように，本発明の第１の実施例を用いれば，パターンの形状や隣接パターン間の距離などによらず，パターンの正確な寸法および形状計測が可能となる。本実施例においては事前準備として多くのシミュレーションが必要であるが，ライブラリ構築は，各製品プロセスについて一度やっておけば後で計算しなおす必要はないので、特に量産ラインにおいてパターン寸法、形状の計測に本発明を適用してプロセスの管理に用いた場合、本発明の効果は顕著になる。
【実施例２】
【００５１】
本発明の第２の実施例について，図５を用いて説明する。図５（ａ）〜（ｃ）はパターン寸法が異なる場合のＳＥＭ信号波形の違いを模式的に示している。図５（ａ）はパターン５０１の寸法が大きい場合を示し、図５（ｂ）はパターン５０２の寸法が（ａ）の場合に比べて少し小さい場合、図５（ｃ）はパターン５０３の寸法が小さい場合を示す。図５（ａ）〜（ｃ）にかけてパターンの寸法が小さくなると，パターン上面部の信号量が５１１，５１２，５１３のように大きく変化する。これは，電子のサンプル材料内での拡散長に対して，パターン線幅が小さくなってしまうために，パターン中央に照射した電子が両側のエッジ部まで拡散し，多くの２次電子が発生してしまうためである。
【００５２】
このように，パターンの寸法，材料，電子線照射条件の組み合せによっては，実施例１のスペース幅だけではなく，パターン寸法によってもエッジ部のＳＥＭ信号波形が変化してしまう。このような状況は近年の半導体パターン微細化によって無視できない影響を計測精度に与える。そこで，このような条件における計測では，左右のエッジの干渉を考慮した計測が必要となる。
【００５３】
第２の実施例は，図５（ａ）〜（ｃ）に示したような、スペース幅だけではなく，パターン寸法によってもエッジ部のＳＥＭ信号波形が変化してしまうという課題に対応するものであり，図５（ｄ）を用いて説明する。なお，図５（ｄ）の手順は，図１（ｂ）で説明したステップＳ０００２，０００３がステップＳ００１０，Ｓ００１１と入れ替わる以外は第１の実施例の図１の手順と同じである。
【００５４】
第２の実施例では，第１の実施例でスペース幅に対して行なった処理と同様に，はじめに寸法幅を変更したデータを含むライブラリを用意しておき，パターン幅をあらかじめ従来法で暫定的に計測する（Ｓ００１０）。シミュレーション結果から事前に求めておいた実寸法と従来法での計測誤差の関係に基づき，暫定計測結果から実際の寸法値を予測，その結果に基づいて線幅のライブラリ条件を限定し（Ｓ００１１），限定されたライブラリによりＳＥＭ実波形とのマッチングを行い（Ｓ００１２），寸法計測を行ない（Ｓ００１３）、計測結果を出力する（Ｓ００１４）。この場合も，第１の実施例と同様に，寸法変動に対する波形変化および計測誤差の変化が小さい場合には設計値や予め指定した固定値を用いても問題ない。
【００５５】
以上のように，本発明の第２の実施例を用いれば，第１の実施例で説明した効果に加えて、ＳＥＭ画像内で左右のエッジが干渉するような微細パターンに対しても正確な寸法および形状計測が可能となる。また、第１の実施例で説明した場合と同様に、事前準備として多くのシミュレーションが必要であるが，ライブラリ構築は，各製品プロセスについて一度やっておけば後で計算しなおす必要はない。
【実施例３】
【００５６】
第３の実施例では，ライブラリとの波形マッチングに非線形最適化を用いた場合に，マッチングを安定かつ高速にする手段について図６から図８を用いて説明する。図６（ａ）は，ライブラリマッチングの問題点を説明する図であり，グラフはマッチングの際の誤差空間（形状パラメタとマッチング誤差の関係）を表している。簡単のために形状パラメタはp1， p2（例えば側壁傾斜角とトップコーナラウンディングなど）の２つのみで示しているが，実際にはさらに多くのパラメタからなる多次元空間であってよい。
【００５７】
等高線は計測対象の実ＳＥＭ画像と各パラメタセットにおけるシミュレーション波形のマッチング誤差（例えば，差の自乗和）を示しており，理想的なケースでは極小値をひとつだけ持ち，周辺に比べて誤差が十分に小さいことが望ましい。しかしながら，実際には画像ノイズやモデルの不完全さのため，図６（ａ）に示すように、複数の極小値を持つことがあり，このような場合には，非線形最適化によるマッチング処理によって，偽の解を選択してしまう恐れがある。場合によっては、偽の解の方が小さな誤差を持ってしまうケースもあり，この場合は数学的に正しい解を選択しても，計測結果は正しくなくなってしまう。
【００５８】
一般的に，非線形最適化においては，初期値の設定が重要であり，解に近い初期値が適切に設定されていれば，正しい解を求めることが容易となる。それに対し，不適切な初期値を与えてしまうと，誤った解を選択したり，収束に多くの時間を要したりという問題がある。そこで，第３の実施例では，この初期値を実際のパターンのＳＥＭ画像特徴量から推定することで適切な値に設定する。
【００５９】
第３の実施例を，図６（ｂ）を用いて説明する。なお，図６（ｂ）の手順は，ステップＳ０００２，０００３がステップＳ００２０，Ｓ００２１と入れ替わる以外の（Ｓ００２２〜Ｓ００２４）のステップは第1の実施例の図１の手順と同じである。本実施例では，ＳＥＭ画像取得後，計測対象パターンの形状をＳＥＭの画像特徴量を用いて推定する（Ｓ００２０）。次に，ステップＳ００２０で推定された断面形状を非線形最適化の初期値に設定し（Ｓ００２１），ＳＥＭの実画像とライブラリ波形のマッチング処理を行い（Ｓ００２２），計測結果を得（Ｓ００２３）、結果を出力する（Ｓ００２４）。
【００６０】
図７にステップＳ００２０で用いる画像特徴量の例を示す。図７（ａ）において，特徴量ｆ１はエッジピーク部分の幅である（以下ホワイトバンドと呼ぶ）。ホワイトバンド幅は，垂直上方からみた場合のエッジ部の見込み幅を反映する特徴量である。特徴量ｆ２は，ホワイトバンド部において，ピーク位置より外側部の平均幅であり，ボトム部の曲率の大きさを反映する特徴量である。特徴量ｆ３は，ホワイトバンド部において，ピーク位置より内側部の平均幅であり，トップ部の曲率の大きさを反映する特徴量である。特徴量ｆ４は信号強度の大きさであり，図７が示すようにテーパ角の大きさを反映する特徴量である。また，絶対信号量の評価が可能なシステムであれば，ピーク部分の絶対信号量ｆ６およびエッジ外側の最小絶対信号量ｆ７なども用いることができる。ｆ６は傾斜角効果により，テーパ角度に依存して変化する値であり，ｆ７はスペースによって変化する。
【００６１】
図７（ｂ）は別の特徴量の例である。エッジピーク部の1次微分を用いて，1次微分が極値を持つ点とゼロになる点の間の距離を特徴量Ｆ１，Ｆ２，Ｆ３としている。Ｆ１はトップコーナの曲率に応じて変化する値であり，Ｆ２は側壁傾斜角，Ｆ３は裾引きに相関がある。図７（ａ）（ｂ）に示すような種々の特徴量を，予めライブラリ内のＳＥＭシミュレーション波形に対して計算しておき，計測時には実際のＳＥＭ信号波形についても同様の特徴量計算を行い，これらの関係から計測対象パターンの形状推定を行い，結果を初期値としてライブラリマッチングを行なう。
【００６２】
次に，図８を用いて，これらの特徴量からパターン形状パラメタの推定を行なう方法について説明する。図８（ａ）は形状パラメタ空間における画像特徴量の分布を曲面８０１と８０２で示すグラフである。ｐ１，ｐ２は側壁傾斜角やコーナ曲率などのシミュレーション形状モデルの形状パラメタ，特徴量ｆＡ，ｆＢは図７に示したようなＳＥＭ画像特徴量を表している。説明を簡単にするため，形状パラメタ，特徴量ともに２つずつの例を示すが，実際にはさらに多くの形状パラメタおよび特徴量を用いてもちろんよい（その場合には，４次元以上の他パラメタ空間となる）。
【００６３】
図８（ａ）のグラフはシミュレーションライブラリの波形に対して，図７（ａ）（ｂ）に示すような種々の特徴量を算出することで得ることができるため，計測時に作成する必要はない。事前に，ライブラリと同時に作成しておく。次に，計測時の実ＳＥＭ画像から各画像特徴量ｆＡ＿ＳＥＭ，ｆＢ＿ＳＥＭを算出すると，図８（ｂ）に示す平面８０３と８０４が得られる。ライブラリの特徴量の曲面８０１又は８０２とこの平面８０３又は８０４の交点に対応する形状パラメタが計測対象パターン形状の候補となる。
【００６４】
そこで図８（ｃ）のように，この平面からの距離に応じた各パラメタの尤度関数を適当に与える。図８（ｃ）はそれぞれのパラメタが尤もらしい確立であるので，図（ｃ）のグラフを掛け合わせて，最大の確立をとる点を算出することで，実ＳＥＭ画像の形状パラメタを推定することができる。この推定された形状パラメタを初期値として用いて，ライブラリマッチングを行なえば，偽の解におちいることなく，安定で高速なマッチング処理が可能となる。ひいては，安定で高速に高精度なパターン形状および寸法の推定が可能となる。これらの初期値設定方法に加え，（特許文献２）に開示されているような手法により，装置の検出器ゲインやオフセット，ビーム径などの装置パラメタを事前に設定してももちろんよい。適切な値に装置パラメタを設定したうえで，パターン形状パラメタの初期値を適切に設定することで，より安定で高速なパターン計測が実現可能となる。
【００６５】
以上のように，本発明の第３の実施例を用いれば，安定で高速なパターン形状および寸法の計測が実現できる。
【実施例４】
【００６６】
第３の実施例では，パターン形状推定の初期値をＳＥＭ画像から得られる特徴量を用いて設定した。第４の実施例では，他の計測装置による計測結果を用いてこれらの初期値を設定する。本実施例ではＡＦＭ計測を組合せて用いる方法について説明する。ＡＦＭは比較的スループットが低いため，プロセスばらつきを考慮した大量計測には向かないが，パターンの概略形状を推定するだけであれば，数スキャン程度のデータ取得で済む。そこで，本実施例では，予め計測対象パターンをＡＦＭで計測し，得られた断面プロファイル形状に基づいて，形状パラメタの初期値あるいはマッチングを行なう形状パラメタの範囲を設定する。。
【００６７】
図９に手順を示す。なお，図９の手順は，図１で説明したステップＳ０００２，０００３がステップＳ００３０と入れ替わる以外のステップＳ００３１〜Ｓ００３３までは第1の実施例の図１に示したＳ０００４〜Ｓ０００６の手順と同じである。計測対象のＳＥＭ画像を取得した後，予め計測しておいたＡＦＭの結果から，シミュレーションモデルに用いた形状パラメタに相当する値を直接計測して，結果をライブラリマッチングの初期値として設定する（Ｓ００３０）。例えば，側壁傾斜角ならば，それぞれのエッジにおける高さ計測結果に直線フィッティングを行ない，その傾きから求めればよい。その後，これらの初期値を用いてＳＥＭ実波形とシミュレーションライブラリのマッチングを行い（Ｓ００３１），寸法計測を行ない（Ｓ００３２），計測結果を出力する（Ｓ００３３）。これにより，第３の実施例と同様に，安定で高速な計測が可能となる。
【００６８】
図９のステップＳ００３０では，ＡＦＭ計測結果に基づき，マッチングで推定される形状パラメタの初期値を決定したが，これとあわせて，マッチングを行なうパラメタの範囲を制限してもよい。パラメタ推定を行なう範囲が適切な範囲に制限されることで，マッチングの際に誤った解を選択してしまう可能性を減らすことができる。これにより，第３の実施例と同様に，安定で高速な計測が可能となる。
【００６９】
また，ＳＥＭ画像評価では，一般に，パターン高さの変化を知ることが困難である。そこで，ＡＦＭデータを用いれば，高さ情報を得ることができるため，第１の実施例のスペースと同様に，高さの異なるライブラリを準備し，得られた高さにライブラリを限定して計測を行なうことも可能である。この場合には，ステップＳ００３０では，初期値ではなく，高さパラメタを固定することになる。本発明により，パターン高さが変動するような場合においても，ＡＦＭとＳＥＭの組み合わせによって安定で，高精度な計測が実現可能となる。本実施例ではＡＦＭを用いる場合について述べたが，Scatterometryのような光学式のパターン形状装置を用いても同様の効果が得られる。
【実施例５】
【００７０】
次に，図１０を用いて，第５の実施例として，本発明のパターン計測方法において用いるシミュレーションライブラリの作成方法について説明する。
【００７１】
シミュレーションライブラリ作成時には，パターン形状のモデルを決めて，モデルに用いられるパラメタにより形状を数値化して表現する。例えば，図４の例では可変パラメタとして側壁の傾斜角θを使って形状を表現していた。図１０（ａ）はこの拡張であり，パターン形状モデル１００１のボトムコーナのラウンディングＲｂもパラメタに加えている。このパターン形状モデル１００１は，計測対象プロセスで起こりうる実際のパターン形状変動に対応したものであることが望ましい。
【００７２】
例えば，エッチングプロセスであれば，エッチング条件が変動しても，レジストで保護されるためトップコーナＲｔ形状はあまり変化しないが，側壁の形状は影響を受けやすく、パターン幅Ｗは変化し易い。特に，ストッパ層の選択比がよくない場合などは，エッチング条件をボトム付近で切り替えるため，切り替えに相当する高さにおいて，側壁の形状が変化することが多い。このような場合には，図１０（ｂ）に示すような台形１００２と１００３を積み重ねた形状モデルが実際のパターンを表現するのに適している。この場合，例えば全体の高さＨは成膜プロセスにより決まり，エッチングプロセスでは変化しないので，エッチングパターンの評価時にはHを設計値に固定するか，あるいは事前に膜厚計で計測した値に設定すればよい。また変化点の高さhはエッチングプロセス切り替え時の変化点であるから，実際のパターンの断面写真などから大体の値を一度決めておけばよい。
【００７３】
このように，形状モデルのうちいくつかのパラメタを評価対象プロセスに適した値に固定すれば，ライブラリマッチング時の推定パラメタを減らすことができ，安定かつ高速な計測が可能となる。図１０（ｃ）はレジストパターンに適した形状モデルの例であり，台形１００４とその頂点で側壁に接する楕円１００５の組み合わせを用いる。例えば，H，hを事前に指定しておけばＷは固定値としてθを設定するだけで楕円１００５が一意に決まるため，θのみを推定パラメタとして，トップ部分のラウンディングＲｔの推定が不要になる。
【００７４】
図１０（ｄ）はこのようなモデルの選択を，断面写真を見ながら設定できるグラフィカルユーザインターフェースの例である。例えば図１０（ｂ）のｈはパターンの断面ＳＥＭ写真を見れば，オペレータが容易に判断することができる。そこで，ライブラリ作成時に，事前に取得しておいた断面ＳＥＭ写真１０１１を画面１０１０に表示し，その上にモデル形状を重ねて表示すれば、画面１０１０を見ながらマウス操作により簡単に設定することができる。
【００７５】
図１０（ｄ）の例では，画面１０１０の左下に表示されたボタン１０１２で使用する断面画像を選択し，右上の画像１０１３の中から用いるモデルを選択する。その後，画面１０１０の右側中断の領域１０１４で、選択したモデルの形状パラメタの固定の有無や範囲を指定できる（図１０（ｄ）では、ボトム線幅１０１５、ピッチ１０１６、高さ１０１７を指定する場合を示す）。パラメタを選択して確定ボタン１０１８を押し、最後にマウス入力ボタン１０１９を押せば，画像上でマウスを用いて，対象パラメタを任意の値に設定できる。このようにして，実際のプロセスにあった適切なモデルの選定とパラメタの範囲設定に基づくシミュレーションライブラリが構築できる。
【００７６】
なお，図１０（ｄ）では断面ＳＥＭ写真の例を示したが，代わりにＳＴＥＭ写真やＡＦＭの計測結果を用いてももちろんよい。また，これらの断面情報は，実際のプロセスの変動範囲を考慮したものが望ましい。例えばレジストパターンであれば，フォーカスと露光エネルギー量を変化させて露光した複数のパターンを作成し，プロセスの正常範囲との境界条件でのパターンを用いれば，当該プロセスで発生しうるパターン形状をカバーしたシミュレーションライブラリが構築できる。
【００７７】
このように，プロセスに応じたモデルの選択と，推定すべきパラメタの選択を適切に行うことができれば，シミュレーション波形と実際のＳＥＭ波形の一致度が向上し，より正確で安定な計測が可能となる。また，全てのパラメタを推定するのではなく，プロセス変動によって変化しないパラメタを固定して，推定パラメタの数を減らすことにより，推定にかかる計算時間の短縮，推定結果の安定化などの効果が得られる。
【実施例６】
【００７８】
本発明のパターン計測方法において用いるシミュレーションライブラリの作成時の材料パラメタ設定方法について図１１を用いて説明する。
【００７９】
ＳＥＭで検出される２次電子信号量は，照射ビームの加速電圧やスキャン速度などのＳＥＭの撮像条件と計測対象の材料によって変化する。しかしながら，材料に依存した信号量の違いを正確にシミュレーションするのは難しい。例えば，同じ二酸化ケイ素のパターンであっても，その成膜方法が異なれば，得られるＳＥＭ信号量は異なるということがしばしば起こる。
【００８０】
一方で，同じ材料パターンでの形状変化に対する相対的なＳＥＭ信号波形の変化は，材料による違いが正確に再現できなくても，比較的正しくシミュレーションできるため，本発明のシミュレーション波形とのマッチングによる計測においても従来法に比べ高い計測精度を実現することができる。ところが，図５（ａ）〜（ｃ）に示したように，パターン部分と下地部分で材料が異なるサンプルのＳＥＭ画像を用いる場合には高精度な計測を行なうことが困難となる。
【００８１】
２次電子信号量は材料に応じて変化するため，例えば図１１（ａ）に示すようなパターン１１０１の材料Ａでできたパターントップ部分と材料Ｂでできた下地部分での信号量は異なる。例えば，図１１（ａ）のパターン１１０１のＳＥＭ画像の信号波形が図１１（ｂ）のような場合に，シミュレーションの結果は図１１（ｃ）の波形のようになってしまうことがある。この場合には，材料Ａ部分とＢ部分のコントラストが反転しているため，どんなにライブラリの形状パラメタを調整しても，実際の画像にマッチするシミュレーション波形を作ることはできないため，正しい計測結果を得ることが困難となる。
【００８２】
本発明では，このようなサンプルを計測する場合に，事前に材料のパラメタを調整することで，ライブラリマッチングによる計測精度を向上させることができる。使用されるモデルに依存するが，通常のＳＥＭシミュレータには材料を設定するためのパラメタがあるので，本発明ではそれを事前に調節する。例えば，非特許文献３に開示されているように，例えばＳＥＭシミュレータMONSELにはresidual energy loss rateという，電子の散乱によるエネルギーの減衰量を調節するパラメタがあり，これを調節すると発生する２次電子信号量が変化する。他のモデルを用いたシミュレータにも，通常上記に類似したパラメタある。ＳＥＭに絶対的な２次電子効率（１次電子の照射量に対する２次電子の割合）の計測が可能なハードウェアがある場合には，計測された２次電子効率にシミュレーション結果が合うようにライブラリ作成前に予め材料パラメタを調整しておく。
【００８３】
しかしながら，一般にはＳＥＭの絶対的な２次電子効率を計測するのは難しく，ＳＥＭシミュレーションの材料パラメタ調整も容易ではない。このような場合には基準材料を定め，それに対するＳＥＭ信号量の割合が実画像とシミュレーションで一致するようにパラメタを設定する。
【００８４】
例えば，図１１（ａ）の場合であれば，例えば下地材料Ｂを基準として，材料Ａの信号量を調整する。はじめに，スクライブライン部などの，比較的大きなパターン部分において，材料Ａ，Ｂそれぞれについてエッジから離れた平坦部のＳＥＭ画像を実際に取得する。次に，カラムと試料室の間のバルブを閉めるなどして，１次電子が試料に当たらないようにして，ＳＥＭ画像をとることで，センサに信号が当たっていない場合のＳＥＭ画像（暗電流）を取得する。このとき，材料Ａ，Ｂ，およびビーム照射の無い状態のＳＥＭ画像は全て必ず同じ検出条件（ビーム条件や検出器の感度など）で取得する。次に，それぞれの画像内の適当な領域内で平均明るさを計算し，材料Ａ，Ｂの明るさからそれぞれ暗電流の明るさを差し引く。その後，基準材料Ｂに対する材料Ａの相対信号量を算出し，シミュレーション結果がこれと一致するように材料Ａの材料パラメタを調整する。
【００８５】
図１１（ｄ）は調整の一例である。予め計測した実信号量を左のｙ軸に示している。この例では基準材料にくらべ材料Ａは６割程度の信号量となっている。これに対し，まずデフォルトのシミュレーション条件で材料Ａ，基準材料Ｂの２次電子効率を計算する（結果は白い丸および白い四角で表示）。図１１（ｄ）では，この基準材料Ｂの結果（白い四角）が実画像と一致するように右のｙ軸をあわせて表示している。このシミュレーション結果では，実画像と異なり，材料Ａの方が明るくなっている。これに対して，材料Ａの材料パラメタをいくつか変更して，２次電子効率を計算した結果が曲線あり，基準材料に対する材料Ａの明るさが実画像と一致する条件は黒丸で示したものとなる。
【００８６】
このようにして，材料間のコントラストが実際に一致するように，相対的な材料間の明るさを調整することで，完全ではないが，ライブラリマッチングが安定かつ高精度に行なえるようになる。基準材料としては，安定で特性のよく知られているものが望ましい。半導体ウェハであれば，基板のＳｉなどを用いることができる。
【００８７】
このように，材料のパラメタを事前に適切な値に設定しておくことで，シミュレーションとの整合性を向上することができ，その結果シミュレーションライブラリを利用したパターン計測手法においても，計測の安定性向上，精度向上が可能となる。
【００８８】
また，第４の実施例のように，他の計測装置によりパターンの概略形状を計測できる場合には，平らな部分の明るさに加えて，エッジ部の信号量も合わせこむ材料パラメタを設定することも可能となる。図１１（ａ）に示したパターンがある場合に，例えばまずＡＦＭで概略形状を計測する。その結果に基づいて，ＡＦＭ計測結果を入力形状として，いくつかの材料パラメタでシミュレーションを行い，材料Ａの平らな部分（上面部）と材料Ａの側壁部分の信号量の関係が実際のＳＥＭ画像に一致するパラメタを選択する。次に，これら材料Ａの平坦部分とエッジ部分の信号量に対する材料Ｂの信号量が実際のＳＥＭ画像と一致する材料Ｂのパラメタを決定すればよい。
【００８９】
本形態によれば，エッジ部分と平坦部における信号量の関係がより実際に近いものになり，さらなる計測の安定性向上，精度向上が可能となる。この方法を用いれば，基準サンプルがない場合においても材料パラメタを合わせこむことが可能となる。
【産業上の利用可能性】
【００９０】
本明細書に開示されているパターン計測技術は，電子顕微鏡あるいはそれに類似した荷電粒子線装置で画像取得が行える対象であれば適用が可能である。本明細書では半導体パターンの計測について説明を行ったが，MEMSや微細工業部品などに適用することができる。
【図面の簡単な説明】
【００９１】
【図１】本発明の第１の実施の形態の説明図である。
【図２】従来のライブラリマッチング手法の課題の説明図である。
【図３】本発明の第１の実施の形態における画像処理領域設定の説明図である。
【図４】本発明の第１の実施の形態におけるライブラリ制限手法の説明図である。
【図５】従来のライブラリマッチング手法の課題と本発明の第２の実施の形態の説明図である。
【図６】従来のライブラリマッチング手法の課題と本発明の第３の実施の形態の説明図である。
【図７】本発明の第３の実施の形態で用いるＳＥＭ画像特徴量の説明図である。
【図８】本発明の第３の実施の形態における形状パラメタ推定手法の説明図である。
【図９】本発明の第４の実施の形態の説明図である。
【図１０】本発明の第５の実施の形態の説明図である。
【図１１】本発明の第６の実施の形態の説明図である。
【符号の説明】
【００９２】
００１・・・ＳＥＭ    １０１・・・電子銃    １０３・・・集束レンズ    １０４・・・偏向器    １０５・・・対物レンズ    １０６・・・試料    １０７・・・検出器    １０８・・・全体制御・信号処理部    ００２・・・ＳＥＭ波形シミュレーションライブラリ    ００３・・・ＳＥＭ画像
００４・・・計測対象パターン    ００５・・・寸法計測処理領域ウィンドウ    ００６・・・左側スペース計測用処理領域ウィンドウ    ００７・・・右側スペース計測用処理領域ウィンドウ
００８・・・ＳＥＭ信号波形    ００９・・・シミュレーション入力断面形状    ０１０・・・ＳＥＭシミュレーション波形    ０１１・・・ＳＥＭシミュレーション部分ライブラリ    ０１２・・・ライブラリマッチング領域    ０１３・・・ＳＥＭ信号１次微分波形    ０１４・・・ＡＦＭ計測結果
【図１】
0005103219.tif 000002
【図２】
0005103219.tif 000003
【図３】
0005103219.tif 000004
【図４】
0005103219.tif 000005
【図５】
0005103219.tif 000006
【図６】
0005103219.tif 000007
【図７】
0005103219.tif 000008
【図８】
0005103219.tif 000009
【図９】
0005103219.tif 000010
【図１０】
0005103219.tif 000011
【図１１】
0005103219.tif 000012