(57)【要約】【課題】ウエハの上方から観察した回路パターンの平面像のみから、回路パターンの断面形状又は３次元形状を推定し評価する。【解決手段】本発明は、立体的な構造体が上面に形成された基板の主面に略垂直な方向から収束エネルギー線を照射するとともに上記基板の上面を走査させ、上記基板及び上記構造体から発生した２次エネルギー線又は上記基板及び上記構造体により反射若しくは散乱されたエネルギー線の強度を検出及び／又は測定し、上記構造体の上面観察像を取得する処理と、上記上面観察像における収束エネルギー線の照射位置と測定された上記強度から上記構造体の表面の凹凸形状による散乱強度の不確かさ情報を求める処理と、求めた上記不確かさ情報に基づいて上記構造体の表面の傾斜角度θを求める処理と、求めた傾斜角度θに基づいて上記構造体の立体形状を推定する処理とを有する。【選択図】図１４
2014130077.tif 000002
【特許請求の範囲】
【請求項１】
計算機に、
立体的な構造体が上面に形成された基板の主面に略垂直な方向から収束エネルギー線を照射するとともに上記基板の上面を走査させ、上記基板及び上記構造体から発生した２次エネルギー線又は上記基板及び上記構造体により反射若しくは散乱されたエネルギー線の強度を検出及び／又は測定し、上記構造体の上面観察像を取得する処理と、
上記上面観察像における収束エネルギー線の照射位置と測定された上記強度から上記構造体の表面の凹凸形状による散乱強度の不確かさ情報を求める処理と、
求めた上記不確かさ情報に基づいて上記構造体の表面の傾斜角度θを求める処理と、
求めた傾斜角度θに基づいて上記構造体の立体形状を推定する処理と
を実行させるパターン形状評価方法。
【請求項２】
請求項１に記載のパターン形状評価方法において、
上記収束エネルギー線は収束電子線であり、上記２次エネルギー線は２次電子である
ことを特徴とするパターン形状評価方法。
【請求項３】
請求項１に記載のパターン形状評価方法において、
上記不確かさ情報を求める処理は、
上記２次エネルギー線の強度の異なる複数のレベルに対応して複数のエッジ点列を抽出し、上記複数のエッジ点列の各座標の設計座標からのずれを算出してエッジ点の変動値を算出するサブ処理と、
上記変動値に基づいて、上記異なる複数のレベルに対応する上記上面観察像上の各位置における前記構造体の表面の傾斜角度を求めるサブ処理と
を有することを特徴とするパターン形状評価方法。
【請求項４】
請求項３に記載のパターン形状評価方法において、
上記エッジ点の変動値を算出するサブ処理は、
上記複数のエッジ点列の空間周波数に依存しない成分を求める工程、
上記エッジ点列からそのパワースペクトルが空間周波数の２乗に反比例する成分を除去する工程、又は、
検出されたノイズ成分を除去する工程
を有することを特徴とするパターン形状評価方法。
【請求項５】
請求項１に記載のパターン形状評価方法において、
上記傾斜角度θを求める処理は、
上記上面観察像のある所定の領域内において、上記不確かさ情報σｘの関数として、上記領域の前記基板の主面に対する傾斜角度θを求める
ことを特徴とするパターン形状評価方法。
【請求項６】
請求項１に記載のパターン形状評価方法において、
上記傾斜角度θを求める処理は、
上記上面観察像の上記強度の分布に対して異なる複数のしきい値を設定し、各しきい値に対して複数のエッジ点列を抽出することにより求めた傾斜角度θを、上記複数のエッジ点列の各エッジ位置における傾斜角度θとし、
上記立体形状を推定する処理は、
上記上面観察像の各点における傾き角度分布を求め、上記傾き角度分布を積分することにより立体形状を推定する
ことを特徴とするパターン形状評価方法。
【請求項７】
請求項１に記載のパターン形状評価方法において、
上記不確かさ情報を求める処理は、
前記上面観察像内の局所領域における２次エネルギー線の強度変動のエッジに垂直及び平行な方向に対する空間周波数特性の差に基づいて、上記局所領域における前記構造体の表面の傾斜角度を求めるサブ処理を含む
ことを特徴とするパターン形状評価方法。
【請求項８】
請求項３に記載のパターン形状評価方法において、
上記エッジは曲線である
ことを特徴とするパターン形状評価方法。
【請求項９】
計算機に、
半導体装置を含む立体的な構造体が上面に形成された基板の主面に略垂直な方向から収束エネルギー線を照射するとともに上記基板の上面上を走査させ、上記基板及び上記構造体から発生した２次エネルギー線又は上記基板及び上記構造体により反射若しくは散乱されたエネルギー線の強度を検出及び／又は測定し、上記構造体の上面観察像を取得する処理と、
上記上面観察像における収束エネルギー線の照射位置と測定された上記強度から上記構造体の表面の凹凸形状による散乱強度の不確かさ情報を求める処理と、
求めた上記不確かさ情報に基づいて上記構造体の表面の傾斜角度θを求める処理と、
求めた傾斜角度θに基づいて上記構造体の立体形状を推定する処理と
推定された上記立体形状の特徴に基づいて上記半導体装置の製造工程における製造条件を推定する処理と
を実行させる半導体装置の製造方法。
【請求項１０】
立体的な構造体が上面に形成された基板の主面に略垂直な方向から収束エネルギー線を照射するとともに上記基板の上面上を走査させ、上記基板及び上記構造体から発生した２次エネルギー線又は上記基板及び上記構造体により反射若しくは散乱されたエネルギー線の強度を検出及び／又は測定し、上記構造体の上面観察像を取得するデータ処理部と、
上記上面観察像における収束エネルギー線の照射位置と測定された上記強度から上記構造体の表面の凹凸形状による散乱強度の不確かさ情報を求めるデータ処理部と、
求めた上記不確かさ情報に基づいて上記構造体の表面の傾斜角度θを求めるデータ処理部と、
求めた傾斜角度θに基づいて上記構造体の立体形状を推定するデータ処理部と
を有するパターン形状評価装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば走査型顕微鏡を用いた非破壊観測及び画像処理による詳細な形状評価方法、そのための装置、及び、当該技術を採用する半導体装置の製造方法に関する。
【背景技術】
【０００２】
半導体集積回路（LSI）は、回路パターンの微細化による高性能化・高集積化が進んでいる。現在、最先端LSIの最小回路パターンの線幅は30nm（ナノメートル）以下であり、LSIの性能を確保するためには、これらの回路寸法を厳密に（例えば寸法のバラツキ変動の許容値が設計値の１割以下となるような精度で）管理する必要がある。上記回路寸法の測定には現在広く走査型電子顕微鏡（SEM）が用いられている。非特許文献１には、ウエハを上面から観察する回路寸法計測専用の電界放出型電子顕微鏡（CD-SEM: Critical Dimension SEM）について記載されている。
【０００３】
上記CD-SEMは、半導体回路の線幅測定の他、様々な特徴量を計測するために用いられている。例えば、上記回路パターンのエッジにはラインエッジラフネス（LER）と呼ばれる凹凸が存在し、回路性能に悪影響を及ぼすことが知られている。CD-SEMは上記LERの計測に広く用いられており、例えば特許文献１にはその測定方法が記載されている。
【０００４】
その一方で、ウエハ上に積層・パターニングにより形成された素子構造の立体的な形状を知りたいという要求がある。特に、LSIの量産工程では、断面観察によらず非破壊でこれを行うことが望ましいが、このための手法として、一般にAFM又は光学的方法（Scatterometry）が用いられている。
【０００５】
AFMは、微細な先端をもつ探針により、探針先端と試料表面の間の原子間力が一定になるように走査することにより、試料表面の凹凸形状を計測する方法である。その詳細は、例えば特許文献２に記載されている。
【０００６】
Scatterometryは、周期的立体構造をもつパターンに光を入射してその反射回折光の波長又は回折角依存性を測定し、これをあらかじめ様々な断面形状に対して計算により求めた回折角依存性と比較することにより立体構造の断面形状を推定する。Scatterometryについては、例えば非特許文献２に記載されている。
【０００７】
Scatterometryと類似の手法であって、SEMを用いて断面形状を推定する方法に、MBL（model-based library）法がある。MBL法は、試料に収束電子ビームを走査して得られる２次電子検出信号強度分布と、予め様々な断面形状に対して計算により求めた２次電子信号強度分布との比較により、立体構造の断面形状を推定する。MBLについては、例えば特許文献３又は非特許文献３に記載されている。
【０００８】
また、SEMを用いて立体構造を測定する方法として、tilt-SEMがある。この方法はウエハに対して異なる角度から電子線を入射して得られた複数の像からステレオ画像の原理により３次元形状を推測する。tilt-SEMについては、例えば特許文献４に記載されている。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２００６−２１５０２０号公報
【特許文献２】特開２００９−２５７９３７号公報
【特許文献３】特開２００７−２２７６１８号公報
【特許文献４】特開２００５−１８３３６９号公報
【非特許文献】
【００１０】
【非特許文献１】Hitachi Review vol.60 (2011), No.5 pp.203-209
【非特許文献２】Solid State Technology, Vol.54, Issue 8, (2011)
【非特許文献３】Proceedings of SPIE Vol. 5375 (SPIE, Bellingham, WA, 2004) Dimensional Metrology of Resist Lines using a SEM Model-Based Library Approach
【発明の概要】
【発明が解決しようとする課題】
【００１１】
半導体集積回路（LSI）の製造工程では、回路パターン又はこれを形成するためのレジストパターン等の断面形状、もしくは３次元形状を所定の設計範囲内に管理することが重要である。これらの形状が設計値通りとなるように製造条件を設定したり、形状が劣化した場合には、いち早くこれを検知し、製造プロセスへフィードバックすることで品質劣化を抑制する必要がある。
【００１２】
しかし、現在まで、任意形状の微細な回路パターンの断面形状又は３次元形状を、非破壊、高精度かつ簡便に計測できる手法が存在しないという技術課題があった。例えばCD-SEMは、任意形状の微細な回路パターンの寸法を非破壊、高精度かつ簡便に計測できるが、ウエハの上方から回路の平面形状を観察するために、断面形状の推定が困難であるという技術課題を有している。
【００１３】
そこで一般には、ウエハを破壊して回路パターンの断面を露出し、その断面形状を電子顕微鏡で観察する手法が用いられる。ところが、製品ウエハに対し、この手法を適用することは難しい。このように、従来、CD-SEMを用いた断面形状の推定は困難であると考えられており、前述したような各種の方法が試みられてきた。
【００１４】
しかし、そのいずれの方法においても、以下のような技術上の課題がある。
例えばAFMによる方法では、回路パターンの寸法が縮小するのに伴い、探針がパターン間に入り込むことができず、形状測定が困難になるという課題がある。
【００１５】
また、scatterometry又はMBLによる方法では、様々な断面形状に対して予想される測定結果を予め計算により求めてライブラリとして用意する必要がある。しかし、これには膨大な計算を要し、適用可能な形状が台形等の比較的単純形状に限定されるという課題がある。特に、scatterometryによる方法では、被計測パターンが、数十ミクロン角の広い領域内に一様に存在する周期パターンに限定されるため、論理LSIの回路のような不規則的なパターンの断面形状を推定することが困難である。また、同方法は、計測パターンのために専用の広い領域を確保する必要があるという課題がある。
【００１６】
また、tilt-SEMを用いる方法では、収束電子線の入射角を変化させるための特殊な電子光学系が必要とされ、一般に、分解能等の性能が劣化する等の課題がある。また、ステージを傾けるために装置が大型化するという課題、計測に時間を要するという課題等も存在する。
【００１７】
そこで、本発明は、荷電粒子線装置を用いて取得した基板の上面からの観察像だけを用いつつ、基板の上面に形成された任意の構造体の断面形状を非破壊かつ高精度で推定できる技術を提供する。
【課題を解決するための手段】
【００１８】
上記課題を解決するために、例えば特許請求の範囲に記載の構成を採用する。一例としての発明は、(a) 立体的な構造体が上面に形成された基板の主面に略垂直な方向から収束エネルギー線を照射するとともに上記基板の上面を走査させ、上記基板及び上記構造体から発生した２次エネルギー線又は上記基板及び上記構造体により反射若しくは散乱されたエネルギー線の強度を検出及び／又は測定し、上記構造体の上面観察像を取得する処理と、(b) 上記上面観察像における収束エネルギー線照射位置と測定された上記強度から上記構造体の表面の凹凸形状による散乱強度の不確かさ情報を求める処理と、(c) 求めた上記不確かさ情報に基づいて上記構造体の表面の傾斜角度θを求める処理と、(d) 求めた傾斜角度θに基づいて上記構造体の立体形状を推定する処理とを有する。
【発明の効果】
【００１９】
本発明によれば、基板の上面観察像のみから、基板の上面に形成された構造体の断面形状又は３次元形状を推定し評価することができる。前述した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。
【図面の簡単な説明】
【００２０】
【図１】電子顕微鏡の観察像形成過程を説明する概念図である。
【図２】被観察構造体の特徴を示す概念図である。
【図３】本発明の原理を説明するための概念図である。
【図４】第１の実施形態を説明するための試料の断面形状の特徴を示す模式図である。
【図５】図４に示す各試料に対する計測結果を模式的に示す特性図である。
【図６】図４に示す各試料に対する計測結果を模式的に示す別の特性図である。
【図７】画像ノイズの影響を示す特性図である。
【図８】図６に示す計測結果に対する解析結果を示す特性図である。
【図９】構造体表面の傾斜角度とエッジ検出点変動幅の関係を示す原理図である。
【図１０】構造体表面の傾斜角度とエッジ検出点変動幅の関係を求める別のモデルを説明するための原理図である。
【図１１】構造体側面の変動が側面から比較的離れた点に入射した電子線に対する検出信号強度に及ぼす影響を示す模式図である。
【図１２】図６に示す計測結果に対する解析結果を示す別の特性図である。
【図１３】図４に示す各試料に対する断面形状の推定結果を示すための特性図である。
【図１４】第１の実施形態における処理手順を説明するフローチャートである。
【図１５】第１実施例で用いる装置の構成例を示す模式図である。
【図１６】第２実施例による３次元形状推定方法の模式図である。
【図１７】第３実施例によるプロセスモニターの模式図である。
【図１８】被観察構造体の特徴例を示す別の概念図である。
【図１９】第２の実施形態における処理手順を説明するフローチャートである。
【図２０】第２の実施形態による解析結果を説明するための特性図である。
【図２１】第５実施例を説明するための概念図である。
【図２２】第５実施例を説明するための概念図である。
【発明を実施するための形態】
【００２１】
以下、図面に基づいて、本発明の実施の形態を説明する。なお、本発明の実施態様は、後述する形態例に限定されるものではなく、その技術思想の範囲において、種々の変形が可能である。例えば以下の説明では電子顕微鏡を用いて試料を観察する場合について説明するが、集束イオンビーム装置その他の荷電粒子線装置を用いて試料を観察する場合や加工する場合についても同様に適用することができる。
【００２２】
（３次元形状観察像におけるエッジ変動要因とその成分分解）
まず、本発明で用いる電子顕微鏡による観察像の形成過程を、図１を用いて説明する。単純化のため、基板面上に形成された略立方体状の構造体（例えば半導体やレジストパターン等）を考える。基板面をｘ−ｙ平面、構造体のエッジの延びる方向（エッジ方向）をy方向とする。上記構造体の特徴寸法より十分に細く収束された電子線を、上記基板面と略垂直な方向（ｚ方向）から基板上に照射すると共に、基板上を上記エッジ方向と略垂直な方向（ｘ方向）に走査する。上記基板又は構造体に入射した電子は、基板又は構造体内部で散乱を受けて２次電子を放出し、又は直接反射（又は後方散乱）され、その一部を基板又は構造体の外部へ放出する。上記２次電子又は反射電子（以下、「２次電子等」という）の放出量は、電子線が凸形状の凸部（又は角の上部）に入射した場合に増大する。
【００２３】
そこで、電子線の走査中に放出された上記２次電子等を検出器で検出し、その検出強度を走査方向における電子線の照射位置ｘに対してプロットすると、図１の下段に示すような２次電子等の検出信号強度分布が得られる。一般に、検出信号強度分布を検出信号の最大値で規格化したときの所定のしきい値レベルで切った時の位置、又は、検出信号強度分布の傾斜が最大となる位置等に基づいてパターン及び非パターンの境界を決定する。２つのパターン境界位置間の距離からパターン寸法が計測される。また、走査をエッジ方向（y方向）の異なる位置で行い、求めたパターン／非パターン境界を結ぶことによりパターンのエッジ形状が得られる。通常、このようにして得られたエッジ形状は、エッジ方向に沿って凹凸状を呈する。この凹凸の大きさをラインエッジラフネスと呼ぶ。本発明者は、上記ラフネスの原因を考察することにより、パターンエッジに平行な方向及び垂直な方向を含む平面（以下、「基板面」という）内の検出信号強度分布から、上記基板面に垂直な方向のパターンの立体形状を推定する方法を考案した。
【００２４】
上述したラインエッジラフネスは、一般的に、図２の上段図（ａ）に示すような構造体のエッジ位置のばらつきに起因するものと考えられている。実際の構造体は立方体ではなく、その側壁は図２の上段図（ａ）に示すように傾斜しており、又は、図２の中段図（ｂ）に示すように表面に凹凸（表面ラフネス）を有する。検出信号強度分布は、傾斜角の変化や表面の凹凸の影響を受けると考えられる。
【００２５】
y方向について一定のサンプリング間隔で規定された開始点からx方向に電子線を走査して２次電子検出信号強度を計測すると、図２の下段図（ｃ）に示すような２次電子検出信号強度の２次元分布が得られる。ここで、２次電子検出信号強度の２次元分布は、（１）構造体の立体形状及び材料特性に起因する強度変化、（２）エッジ方向に沿った構造体のｘ方向位置変動による強度分布の変動、（３）構造体表面の凹凸に起因する強度分布の変動の３つの要因の影響を受けると考えられる。
【００２６】
そこで、これらの３つの要因を分解して考える。まず、ｙ座標を固定して、ｘｚ平面でパターンを切断した時の断面形状と、当該形状に対応する１次元信号強度分布について考える。この場合、図２の中段図（ｂ）の構造体は、図３の上段図（ａ）の左側に模式的に示すようなマクロなパターン構造の表面に、図３の中断図（ｂ）の左側に模式的に示すようなミクロな表面の凹凸構造が重なった形状を有するものと考えられる。図３の下段図（ｃ）の左側にマクロパターンとミクロパターンが重畳したイメージを示す。
【００２７】
前述したように、２次電子放出量は、電子線が凸部に入射すると増大し、逆に電子線が凹部に入射すると減少する。従って、パターン全体のマクロな構造を見ると、パターンの上面角部は凸形状なので、上面角部へ向かって信号強度が増大する（図３の上段図（ａ）の右側）。一方、ミクロな構造を見ると、パターン表面に存在する凹凸の凸部に入射したとき、信号強度は増大する（図３の中段図（ｂ）の右側）。従って、これに対応して２次電子信号強度分布は、マクロなパターン構造に対応した大きなピーク状分布（図３の上段図（ａ）の右側）に、表面凹凸構造に対応した細かい信号変動（図３の中段図（ｂ）の右側）が重畳した分布（図３の下段図（ｃ）の右側）となる。但し、構造体の側面では凹凸に対する電子線の入射角度が変化するため、凹凸による信号変動も変化することに留意する。
【００２８】
ここで、ｙ方向のある位置からｘ方向についての２次電子信号強度の分布を考える。上記表面凹凸はランダムに生じると考えられるので、対応する信号強度分布も図３の下段図（ｃ）の右側における実線から一点鎖線のように不確実性をもってばらつく。測長SEMでは、通常、得られた２次電子信号強度分布をあるしきい値で切ったときのｘ座標をエッジ座標として検出する。
【００２９】
ところで、しきい値を変えると、構造体の異なる高さ位置ｚに対応する複数のエッジが検出される。上記エッジ座標は不確実性をもって変動し、その変動量σの特徴はしきい値の関数となる。すなわち、構造体の異なる高さ位置ｚ（又はエッジと垂直方向の位置x）の関数となる。
【００３０】
後述する第１の実施形態では、上記関数の特徴から構造体の立体（断面）形状z(x)を推定する手法について説明する。また、第２の実施形態では、２次元信号強度分布を用い、２次元立体形状を推定する手法について説明する。なお、第２の実施形態の詳細は、後述する実施例４で説明する。
【００３１】
（第１の実施形態）
以下、本発明の第１の実施形態による構造体の立体（断面）形状推定方法の詳細について説明する。以下では、図４の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）に示す断面プロファイルを有する３種類の試料Ａ、Ｂ及びＣを解析対象とする。
【００３２】
まず、CD-SEMにより、前記３種類の断面プロファイルを有する構造体の上面観察画像を取得し、上記構造体のエッジ部を解析領域に指定する。対象エッジには、設計上の直線部分を選択し、上記直線に沿った方向をｙ方向とする。
【００３３】
次に、エッジ検出時のしきい値Ｔを様々に変えてエッジを検出し、Ｔの関数としてエッジ点列 {xi(T)} (ただし、i=1,2,  ,n)を求め、上記点列の平均エッジ位置<x(T)>からの偏差σ2 (T) ＝ Σ[xi(T)-<x(T)>]2 を計算した。次に、σ(T)をしきい値Ｔに対してプロットすると、図５の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）が得られる。前述したように、各段に示す点列は、図４の同じ段の断面プロファイルを有する構造体に対応する。
【００３４】
一般に、検出信号強度のｘ方向分布I(x)は、構造体のエッジ近傍にピークをもつ。例えばレジストラインパターン等、ある幅Ｗで凸構造をもついわゆるラインパターンを観察すると、図４に示す各パターンの左右エッジの各々に対応して２つの強度ピークが現れる。つまり、片方のエッジを観察すると、同じしきい値Ｔに対して検出信号強度の分布ピークの両側に２つのエッジ検出点が得られる。よって、σ(T)は、上記ピークの構造体の内側と外側の各々に対して別々に求める。
【００３５】
一般に、構造体の外側（スペース側）では、外側から強度が徐々に増大して所定のしきい値に達した点のｘ座標をエッジ位置とする。一方、構造体の内側（ライン側）では、内側から強度が徐々に増大して所定のしきい値に達した点のｘ座標をエッジ位置とし、これらの値からσ(T)を求める。なお、パターン幅Ｗが極端に小さくなると、左右のエッジに対する上記ピークが重なり、１本のラインパターンに対して１個のピークが観察されるようになる。この場合、しきい値とｘ座標は１対１に対応する。
【００３６】
次に、解析領域内で、検出信号強度の分布I(x)のy方向の平均分布<I>(x)を求め、各しきい値Ｔを、検出信号強度が<I>(x)＝Tとなる特定のｘ座標に対応させる。これにより、しきい値Ｔの関数としての偏差σ(T)を、ｘの関数としての偏差σ(x)に変換する。このとき、上記ピークの構造体の外側/内側のxに対する偏差σ(T)は、各々の構造体の外側/内側のxに対応させる。偏差σ(x)を検出信号強度分布<I>(x)とともにプロットすると、図６に示す特性図となる。図６の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）に示す特性図は、いずれも図４の同じ段の断面プロファイルを有する構造体に対応する。図６の各図は、位置ｘに入射した電子線に対して得られる平均的なエッジ位置における検出エッジ変動を表わす。
【００３７】
膜減りが大きく、かつ、順テーパ形状を有する試料Ａの場合（図４の上段図（ａ））、偏差σ(x)は構造体の外側下部の平坦面で大きく、入射電子線がエッジに近づくと急激に減少して最小値となり、その後、構造体の中心に向けて徐々に増大し、構造体の上部の平坦面近傍で極大をとる。
【００３８】
一方、ほぼ垂直又はやや逆テーパ形状の側壁を有し、かつ、膜減りが少ない試料Ｂ（図４の中段図（ｂ））、試料Ｃ（図４の下段図（ｃ））では、偏差σ(x)は、最小値を取った後、急激に増大して台地状のピークを取り、いったん減少した後、構造体の中心に向けて徐々に増大し、構造体の上部の平坦面近傍で極大となる。上記ピークの高さは、逆テーパ状側壁をもつ試料Ｃの方がやや大きい。
このような分布形状の差異は、断面プロファイルの違いを反映したものと考えられる。
【００３９】
次に、偏差σ(x)の形状のもつ意味について説明する。
まず、解析領域内部において、平均検出信号強度分布<I>(x)に対応して、便宜的に平均断面プロファイル<Z>(x)を考え、実際の構造体の3次元形状Ｚ(x,y)を次式で表わす。
【００４０】
【数１】
2014130077.tif 000003
【００４１】
すなわち、実際の３次元形状Ｚ(x,y)は、平均断面プロファイル<Z>(x)をエッジ点位置yに応じてｘ方向にΔx(y)シフトし、さらに断面プロファイルに変動ΔZ(x,y)を与えたものと考える。ここで、Δx(y)は、各y毎に、例えば|ΔZ(x,y)|のx方向積分値が最小になるように決める。
【００４２】
次に、上記表現に対応して、実際に観察される検出信号強度の３次元分布も、次のように表わすことができる。
【数２】
2014130077.tif 000004
【００４３】
ここで、ΔI(x, y)は、ΔZ(x,y)に対応して断面プロファイルの平均からのずれによって生じる項である。Δnoise(x,y)は、実際の検出画像に重畳するランダムな検出ノイズであり、その振幅は位置・画像に依らない。Δx(y)は、各y毎に、例えば|ΔI(x,y)|2のx方向のエッジ近傍における積分値が最小になるように決める。
【００４４】
さらに、上記表現に対応して、検出されるエッジ点列の平均からの偏差σ(x)も、上記各項に起因する成分に分割することができる。
【数３】
2014130077.tif 000005
【００４５】
ここで、σ_measuredは計測結果である。σ_yは、各y毎に、その断面を平均断面プロファイルで最適フィットしたときの構造体のエッジ位置のバラツキに起因する成分である。σ_xzは、構造体の断面形状の変化に起因する成分である。σ_noiseは、検出（画像）ノイズに起因するエッジ検出誤差に起因する成分である。
【００４６】
（エッジ変動成分の分解方法の詳細）
次に、計測されたエッジ変動値σ_measuredを、上述した４つの成分に分解する方法について説明する。以下、各成分について説明する。
【００４７】
構造体のエッジ位置のｙ方向の変動は、一般に、その空間周波数特性が1/f特性に従うのに対し、表面凹凸に起因する計測ばらつき及び検出ノイズは、その空間周波数特性がランダムであると考えられる。そこで、σ_yは、ｙ方向の空間周波数fに対して1/f特性を有する成分（便宜上「真のＬＥＲ」と呼ぶ）と考える。
【００４８】
なお、エッジに沿って測定されたｘ方向のエッジ点検出位置のばらつきを、1/f特性を有する成分と、それ以外の成分に分解する方法については、特許文献１に説明されている。ここで、ＬＥＲ成分の分解は、各しきい値毎に行ってもよいし、又は、ある代表的なしきい値Tに対して行ってもよい。
【００４９】
真のＬＥＲは、構造体のエッジ位置自体のｙ方向の変動である。よって、しきい値Tを変えて求めた計測結果に対してはほぼ共通と考えることができる。このため、後者の考えも合理的である。
【００５０】
σ_noiseは、信号強度分布の傾斜に依存する。平滑平坦面に対して得られた画像に、1次元(x)方向に強度勾配を加え、平均レベルのしきい値に対して得られたエッジ点の変動を求めた結果を図７に示す。図７に示すように、検出信号の強度勾配が増えると、σ_noiseは小さくなる。σ_noise(x)は、平均検出強度分布の位置xにおける勾配に対するσ_noiseを図７から求めることにより計算することができる。
【００５１】
計測結果σ_measuredから、σ_y及びσ_noiseの寄与を除くと、σ_xzを求めることができる。
ここで、図４の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）に示すような断面プロファイルを有する試料Ａ、Ｂ及びＣに対してσ_y及びσ_noiseを計算し、想定結果σ_measuredを、σ_y、σ_noise、σ_xzに分解した結果を図８の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）に示す。図８から分かるように、断面プロファイルが逆テーパ形状である場合、σ_xzに顕著なピークが現れる。
【００５２】
以下では、顕著なピークが発生する原因を考察するために、表面の凹凸が検出信号強度分布に及ぼす影響を次の２つに分類する。
【００５３】
［分類１］
この分類では、電子線は、最初に順テーパ面に入射すると考える。このとき、構造体内部で散乱された後、再び構造体外に出て検出される電子の数は、順テーパ面上の入射位置のごく近傍の凹凸パターンの影響を受けると考えられる。
【００５４】
［分類２］
この分類では、電子線は、エッジ点から離れた構造体上部における比較的平坦な面に入射すると考える。このとき、構造体内部で散乱された電子の一部が構造体の側面又は逆テーパ面から構造体の外部へ抜け出し、検出される。このとき、外部に抜け出す電子は、側面の凹凸パターンの影響を受ける。このように入射面とは異なる面から外部に抜け出す電子は、その絶対数が少ない。しかし、表面散乱により検出される電子数も比較的少ないので、側面の凹凸パターンの影響を受けた電子が順テーパ面に対する検出結果に及ぼす影響は無視できないと考えられる。
【００５５】
この影響は、構造体内部の電子線進入長の範囲で、垂直又は逆テーパ面の面積に依存し、垂直又は逆テーパ面高さが大きいほど大きい。また、その影響は、構造体のエッジから内部寄りの比較的広い領域に入射した電子線に及ぶ。
【００５６】
図８の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）に見られるσ_xzのピークは、垂直又は逆テーパの側面を有する構造体において、電子線が構造体の上部平坦面に入射した場合に見られる。このことから、構造体の上部平坦面に出現するσ_xzのピークは、分類２に起因すると考えられる。
【００５７】
そこで、σ_xzを、以下の２つの成分に分解した（式４）で表す
(1)電子線入射点近傍の順テーパ面上の凹凸パターンの影響
：σ_xz_near
(2)電子線入射点から比較的離れた略垂直又は逆テーパ面上の凹凸パターンの影響
：σ_xz_far
【００５８】
【数４】
2014130077.tif 000006
【００５９】
以下、上記(1)及び(2)の成分の詳細についてさらに説明する。
まず、凹凸パターンを有する順テーパ面に入射した電子線の散乱に及ぼす凹凸の影響（上記(1)の成分）について考える。ある一定角度の傾斜面に正弦波状の凹凸が存在し、凹部又は凸部の中心（図９の上段図（ａ）の点P）にビームが入射したとする。ビームの入射位置に対して凹凸の位相を180度変化させると、検出位置が、位相の±方向に１周期分変動する。例えば、ビームが凹部の中心に入射した場合、検出位置は位相の＋方向に最大のずれが生じ、ビームが凸部の中心に入射した場合、検出位置は位相の−方向に最大のずれが生じる。
【００６０】
以下では、説明を簡単にするため、図９の上段図（ａ）に示すように、平均的な傾斜角θが一定であり、かつ、その表面に沿って凹凸パターンをもつ側壁形状を有する構造体の平均側壁面上の点Ｐを考える。点Ｐ上に代表周期Ｌの表面凹凸の凸部中心が存在する場合、電子が点Ｐに向かって入射したとき、検出される信号強度の増加分が最大となり、エッジ検出位置のシフト量が−方向に最大となるとする。逆に、点Ｐ上に表面凹凸の凹部中心が存在する場合、検出される信号強度の減少分が最大となり、エッジ検出位置のシフト量が＋方向に最大となるとする。
【００６１】
このとき、凹凸によるエッジの検出位置のシフト量は、凸部中心のｘ座標位置により決まり、その変動振幅Δｘは、上記２つのケースにおける凸部中心間のｘ方向距離となる。従って、側壁面上の凹凸周期をLとし（すなわち、上記ｘ方向の距離はL/2であり）、検出位置が図９の中段図（ｂ）で決まるとすると、エッジの検出位置の変動幅ΔXは、凹凸周期L、側壁傾斜面の角度（傾斜角）θにより、次式のように推定することができる。
【００６２】
【数５】
2014130077.tif 000007
【００６３】
これに対し、検出位置が図９の下段図（ｃ）で決まるとすると、エッジの検出位置の変動幅ΔXは、高さH、側壁傾斜面の角度（傾斜角）θにより、次式のように推定することができる。
【００６４】
【数６】
2014130077.tif 000008
【００６５】
一般に、傾斜角θが比較的小さい領域では式５が支配的となり、傾斜角θが比較的大きい領域では式６が支配的となる。実際には、傾斜角θは、側面の凹凸周期と高さの両方に依存すると考えられる。従って、両式を加算した次式により、エッジ検出位置の変動幅ΔXを推定することも可能である。
【００６６】
【数７】
2014130077.tif 000009
【００６７】
さらに、エッジ検出位置の変動幅ΔXと傾斜面の傾斜角θとの関係は、表面凹凸による電子線信号の強度変動量が表面に対する電子線の入射角に依存するという特性を用いて求めることもできる。
【００６８】
例えば表面に適当な周期L又は高さHの凹凸を有する平面（図１０の上段図（ａ１））に、様々な入射角度の電子線を照射・走査したときの検出信号強度分布をシミュレーションにより計算する（図１０の上段図（ａ２））。このとき、検出強度の変動幅ΔIは、入射角θの関数となる（図１０の中段図（ｂ））。
【００６９】
一方、強度変動幅ΔI(θ)と、エッジ検出位置の変動幅ΔXと、平均信号強度分布I(x)の関係は、図１０の下段図（ｃ）に示すように、次式で表すことができる。
【００７０】
【数８】
2014130077.tif 000010
【００７１】
ただし、式８は、検出信号の強度分布のピークにおいて勾配が０（zero）になると、値が発散する。このため、式８は、検出信号の強度分布のピークにおいて勾配が０（zero）の場合には用いることが出来ないことに注意が必要である。実際には、エッジシフト量の振幅ΔXは、側壁表面凹凸によるエッジ検出位置の計測ばらつきから推定される統計量であり、σ_xz_nearがこれに相当する。また、側壁表面の凹凸パターンはランダム量であるので、その周期L及び高さHに関しても統計的な代表値を用いる。例えば便宜上、平坦面の傾斜角θが０（zero）となるように選んでもよい。
【００７２】
次に、上記(2)の成分（ビームの入射位置から離れた表面位置の凹凸パターンの影響）について説明する。
平均的な側壁面に対し、図１１の上段図（ａ）のＱ点に電子線が入射したときの検出強度をIQとする場合に、同図内の曲線Ａ及びＢに示すように側壁面のパターンが変動した場合を考える。このとき、図１１の下段図（ｂ）に示すように、検出強度は検出強度IQに対してΔIA及びΔIBだけ変化し、検出エッジはそれぞれΔxA、ΔxBだけずれる。
【００７３】
このとき、側壁面の大きさに対して電子線の広がりが同程度であるとすると、検出エッジに出現するずれ量ΔxA及びΔxBは、入射位置（Ｑ点）に依らず、ほぼ一定と考えることができる。
【００７４】
そこで、電子線の入射位置が、構造体のエッジ位置から構造体中の散乱電子が側壁に及ぶ範囲の間では、前述した(2)の影響がほぼ一定であると仮定し、σ_xz_far に対して台形状の分布形状として近似する。
【００７５】
台形の高さは、側壁の高さと角度に依存するが、例えば式４に従って、σ_xzの最小値から構造体側において、σ_xzからσ_xz_farを引いた後のσ_xz_nearの分布が、０(zero)からほぼ線形に増大するような高さに設定することが考えられる。図４の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）に対応する各試料について、σ_xzを上記２つの成分に分解した結果を図１２の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）に示す。
【００７６】
（断面形状の推定方法）
次に、不確かさを表す偏差σを分解した結果から、構造体の断面形状を推定する方法について説明する。
【００７７】
実際の構造体の断面形状は、丸みを帯びていたり、裾を引いていたり、様々な形状を有するものと考えられる。つまり、一般的に、側壁角θは、構造体の高さ方向又はｘ方向の位置に対して一定ではない。本実施形態に係る手法を用いれば、このような構造体であっても、その断面形状を適切に推定することができる。
【００７８】
以下、解析的形状推定法とモデルベースによる形状推定法の２つのアプローチによる断面推定法について説明する。
（解析的推定法）
解析的推定法では、しきい値又はｘ方向の位置に対し、式５、式６、式７、式８等を用いて、構造体の対応する高さにおける構造体表面の局所的な角度を求め、これをｘ方向に積分することにより構造体の断面形状を求める。すなわち、次式により、断面形状を求める。ただし、次式において積分範囲は０(zero)からxまでである。
【００７９】
【数９】
2014130077.tif 000011
【００８０】
ここでは、式５と式８を用い、図１２のσ_xz_nearより推定したθ(x)の分布と、推定される断面形状Z(x)を図１３の上段図（ａ）、中段図（ｂ）及び下段図（ｃ）にそれぞれ示す。推定された断面形状（実線で示す形状）は、別途、断面観察より得られた形状（図４の上段図（ａ）、中段図（ｂ）、下段図（ｃ）の形状）とよく一致していることが分かる。
【００８１】
一方、モデルベースによる形状推定法では、様々な断面形状に対して、その表面の凹凸（例えばその位相）を変化させたときの検出信号強度の分布を求め、各々に対して検出されるエッジ位置間のバラツキ（不確かさ）を、しきい値又はｘ方向の位置の関数として予め求めておく。次に、このバラツキ（不確かさ）と実際の測定結果とをマッチングさせ、最も近い断面形状を求める又は内挿・外挿により断面形状を推定する。
【００８２】
検出信号の強度分布自体の計算結果ライブラリと計測結果とのマッチングを求めることにより断面形状を推定する従来のMBL法は、測定結果が検出系や検出強度変動に影響を受けるという課題があるのに対し、本方法は、これらの影響を受け難い。もっとも、本方法と従来のMBL法とを併用することもできる。また、マッチングには最尤法等を用いてもよい。
【００８３】
（構造体立体形状推定処理手順）
図１４に、第１の実施形態に対応する、構造体の立体形状を推定するための方法のフローチャートを示す。なお、後述する一連の処理は、計算機で実行されるプログラムに基づいて実現される。
【００８４】
まず、電子顕微鏡により２次元画像を取得し、解析領域を指定する（ステップ１４０１）。この後、指定された解析領域における平均的な信号強度分布Ｉ（ｘ）を算出する（ステップ１４０２）。信号強度は、指定範囲における最大強度で規格化する。この際、上記範囲における平均的なエッジ方向がｙ方向となるように画像を調整することが望ましい。
【００８５】
次に、しきい値を指定の最小値から最大値まで所定の間隔毎（例えば５％から１００％まで５％毎）に設定し、各しきい値毎にパターンエッジを検出し、しきい値Ｔの関数としてＬＥＲを求める（ステップＳ１４０３〜１４０８）。ここでは、Ｔ＝Ｉ（ｘ）の関係を用いてＬＥＲをｘの関数に変換した後、これを各成分に分解する。さらに、式５、式６、式７、式８等に従って側壁角度θを求め（ステップS１４０９、１４１０）、側壁角度θをｘ方向に積分して断面形状を得る（ステップＳ１４１１）。ここまでの説明では、構造体の表面の凹凸パターンは、表面において等方的であると仮定した。局所的なレジストの溶解性ばらつき等に由来する凹凸に関しては、上記の仮定は、ほぼ正しい。一方、この仮定が成り立たないケースについては、後述する実施例５に述べる。
【００８６】
［実施例１］
本実施例では、前述した推定手法をCD-SEMに適用した実施例について説明する。
（装置構成）
図１５に、本実施例で用いたCD-SEMのハードウェア構成の模式図を示す。本実施例のCD-SEMは、主として電子光学カラム（SEMカラム）と試料室からなる走査型電子線顕微鏡の筐体1801と、走査型電子線顕微鏡の制御系1811と、情報処理装置1812により構成される。
【００８７】
情報処理装置1812には、得られた走査電子画像や解析に必要なCADデータなどを格納するデータ記憶装置1813が接続されている。もっとも、データ記憶装置1813は、情報処理装置1812内に格納されていても良い。図示されてはいないが、情報処理装置1812は、CD-SEMの操作者がデータ処理のために必要な情報を装置に対して入力するための情報入力端末及び取得される走査電子画像を表示するための画像表示手段を備えている。情報入力端末の具体例には、キーボード、マウス、画像表示手段上に表示されるGUI画面等がある。
【００８８】
電子光学カラムは、電子銃1802、収束レンズ1804、偏向器1805、対物レンズ1806、検出器1810などにより構成される。試料室は、被検査対象である観察ウエハ1807を載置するステージ1808を備える。電子銃1802から観察ウエハ1807に対して照射される電子線1803により発生した二次電子1809は検出器1810によって検出され、制御系1811によってデジタルデータに変換された後、情報処理装置1812に転送され、解析に使用する画像データが生成される。
【００８９】
本実施例では、CD-SEMに備えられた走査電子顕微鏡を用いてパターン観察を行なって被検査対象の画像データを取得した。得られた画像データをデータ記憶装置1813に保存し、観察終了後、情報入力端末を操作して画像データ解析を実行し、ラフネス指標及び断面構造推定の解析を行った。解析処理は、情報処理装置1812が実行する。
【００９０】
（画像取得工程）
まず、制御系1811（情報処理装置1812）は、ArＦレジストのラインパターンに対し、視野の左上から右下までの走査を32回行って得られた二次電子信号強度を平均化し、２次電子信号強度の２次元分布画像を取得する。必要に応じ、情報処理装置1812は、取得した画像をCD-SEMのモニタ画面上に表示する。ここで、観察画像の画素数は、縦・横方向に1500画素、１画素の１辺を1nm（ナノメートル）、視野内の観察画像の長さを縦横1.5μｍ（マイクロメートル）とする。
【００９１】
情報処理装置1812は、２次電子信号強度の２次元分布画像のうち解析したいエッジを含む領域に、それぞれ縦1024画素、横50画素の矩形形状の検査領域を設定する。また、エッジ抽出のためのしきい値Tの最小値Tmin、最大値Tmax、増分値ΔT、データ抽出のy方向のサンプリング間隔Δy、x方向のノイズ低減パラメータ、y方向の平均化パラメータS等、データ系列抽出に必要な情報を設定する。なお、データ抽出におけるy方向のサンプリング間隔ではなく、検出点数を設定することも可能である。これらの検査領域及びデータ系列抽出の条件は、例えばCD-SEMのGUI画面等を通じ、予め設定することが望ましい。
【００９２】
次いで、情報処理装置1812は、領域内のエッジラフネスのデータ系列を抽出するタスクを実行する。すなわち、情報処理装置1812は、設定された抽出開始点やサンプリング間隔に従って、検査領域内の画素データからサンプリング位置でのy座標に対応する信号強度分布を算出し、更に、最小値Tminから最大値Tmaxまで増分値ΔT毎に設定したしきい値Ｔに応じ、上記信号強度分布からエッジ点のx座標データを算出する。これらの処理を、設定されたサンプリング間隔に従って次々にy方向へ移動しながら実行し、最終的に、しきい値Tの関数として検査領域内でのエッジラフネスのデータ系列X(T)＝{Δxi(T)}＝{Δx1(T)、Δx2(T)、…}を取得する。
【００９３】
具体的には、情報処理装置1812は、データ抽出開始点のy座標として、上記検査領域の下辺に対応するy座標を設定し、ｙ方向のサンプリング間隔として1nm（ナノメートル）を設定し、検査領域内のエッジ点として1nm（ナノメートル）置きに1024個の点の位置（x1(T), y1(T)）,…（xi(T), yi(T)）,…(x1024(T), y1024(T))を抽出した。
【００９４】
次に、情報処理装置1812は、抽出された点を以下の直線で近似し、フィッティングパラメータであるαとβの値を求めた。
【数１０】
2014130077.tif 000012
【００９５】
次に、全てのしきい値Tに対するエッジ点の座標に対し、情報処理装置1812は、以下の式１１に従ってエッジ点の上記直線からのずれ量Δxiを求め、エッジラフネス系列X(T)＝{Δxi}Tを生成する。
【数１１】
2014130077.tif 000013
【００９６】
（画像処理工程）
次に、情報処理装置1812は、エッジラフネス系列X(T)＝{Δxi(T)}からしきい値毎にＬＥＲを求め、求めたＬＥＲを真のＬＥＲによる測定ばらつきσ_yと、表面凹凸による計測ばらつきσ_xzとに分解する。エッジラフネス系列から、真のＬＥＲによる測定ばらつきσ_yと表面凹凸による計測ばらつきσ_xzを求める方法については、例えば特許文献１に記載されている方法を用いることができるが、その代表的な手段について述べるならば以下の通りである。
【００９７】
前述したように、測定されたＬＥＲには、パワースペクトル密度が空間周波数fの２乗に反比例する成分（真のＬＥＲ）と、それ以外の高周波成分（ノイズ）とが重畳している。測定されたＬＥＲに平均化処理を施すと、後者の成分が減少する。このため、平均化処理の程度を表すパラメータ値を増やすのに伴い、高周波領域でのパワースペクトル密度分布がｆの２乗に反比例するようになる。具体的には、異なるｙ座標におけるｘ方向の１次元信号強度分布を、ｙ方向にS本平均化することを考えると、ランダムノイズの強度は、平均化によって1/Sに低減される。すなわち、平均化パラメータSの増加に伴って、パワースペクトル密度の周波数依存性が1/f2に近づいていく。このとき得られるパワースペクトルの密度が真のＬＥＲとなる。
【００９８】
ここで、平均化パラメータSで平均化処理されたデータから得られるラインエッジラフネス指標の値をσm(S)として、そのS依存性を次式にフィッティングする。
【数１２】
2014130077.tif 000014
【００９９】
ただし、σm(1)は、平均化処理前のデータから求めたラインエッジラフネス計測値であり、Δyはエッジ点のｙ方向抽出間隔であり、Ａはフィッティングパラメータである。S及びエッジ点の抽出間隔Δyは、2SΔy＜1/f0 [nm（ナノメートル）]を満たすことが望ましい。ここで、f0は、通常のレジストパターンの場合、スペクトルの屈曲点において、0.008nm-1以下であることが多い。
【０１００】
ここでは、実験結果をよく説明するSの中で最も小さいSに対して求められたσ0＝σを真のＬＥＲとした。この場合、空間周波数によらない計測ばらつき成分σｂは、次式により求められる。
【０１０１】
【数１３】
2014130077.tif 000015
【０１０２】
さらに、情報処理装置1812は、測定されたＬＥＲから計測誤差σ_noiseを求め、これをσｂ（Ｔ）から取り除くことにより、側壁表面の凹凸パターンがｘｙ平面に投影されることによる計測ばらつき（不確定性）を示すσ_xz（Ｔ）を求めた。
【０１０３】
【数１４】
2014130077.tif 000016
【０１０４】
ここで、構造体のしきい値Ｔに対応する高さにおける側壁角度θは、次式により求められる。
【０１０５】
【数１５】
2014130077.tif 000017
ここで、σｓは、表面凹凸の空間周期の代表値である。
【０１０６】
一方、ｙ方向の計測範囲に対して、ｘ方向の信号強度分布の平均<Ｉ>（ｘ）を算出した。しきい値Tに対応するエッジ座標ｘは、ｘ＝<Ｉ>-1（Ｔ）より求めることができる。ただし、<Ｉ>-1は、<Ｉ>の逆関数である。
【０１０７】
ところで、構造体の断面形状ｚ（＝Ｚ（ｘ））は、前述したように、次式により求めることができる。
【０１０８】
【数１６】
2014130077.tif 000018
ただし、積分範囲は、Ｔ＝Ｉa（０）からＴ＝Ｉa（ｘ）である。
【０１０９】
このようにして求めた断面形状を、ウエハのパターン部を切断してSEMで観察した断面形状と比較を行ったところ、両者はよく一致することが確認された。また、ＡＦＭによる計測結果との比較も行ったが、こちらもよく一致することも確認された。
【０１１０】
なお、上の説明ではσｓとσ_noiseは物理的な意味を持つ値として、各々独立に計測して求めるとしたが、これらの量はフィッティングパラメータとして考えてもよい。すなわち、他の方法により観察した断面形状と、本発明により推定した断面形状がよく一致するような、σｓとσ_noiseを用いればよい。
【０１１１】
［実施例２］
本実施例においては、１次元マスクパターンの断面構造だけでなく、２次元マスクパターンに対しても立体的構造を推定できる方法の一例を説明する。本実施例で用いる走査形電子顕微鏡の構成は、実施例１と同様なので説明を省略する。
【０１１２】
実施例１の手法を２次元マスクパターンに拡張する方法には次の２つの方法が考えられる。
［第１の方法］
この方法では、実施例１の方法により１次元パターンに対してしきい値とパターン高さの関係を求め、この関係をそのまま２次元パターンに対して適用する。
【０１１３】
［第２の方法］
この方法では、実施例１の方法を２次元エッジ情報に対して拡張する。
まず、第１の方法について説明する。情報処理装置1812は、実施例１の方法と同様に取得した２次元画像（図１６の上段左図）に対し、その１次元パターン部分に実施例１の方法を適用して断面形状ｚ＝Ｚ（x）を求める（推定する）。
【０１１４】
次に、情報処理装置1812は、上記ｚ＝Ｚ（x）と、上記形状の導出に用いた信号強度分布Ｔ＝Ｉ（ｘ）の２つの関係からｘを消去することにより、構造体の高さｚとしきい値Ｔの関係ｚ＝Ｚ’（T）を得る。図１６の上段左図の２次元マスクパターンに対するSEM観察画像に対し、しきい値Tを変えてエッジを抽出すると、図１６の上段中図のような等高線が得られる。
【０１１５】
次に、情報処理装置1812は、上記構造体の高さｚとしきい値Tの関係に従い、図１６の上段中図のTをｚに変換し、図１６の上段右図のような２次元パターンの高さ情報を得る。なお、Tとｚの間には、図１６の下段に示す関係が成立する。この関係を与えるデータは、例えばデータ記憶装置1813に予め格納しておく。
【０１１６】
本方法の問題は、実施例１の方法により求めた構造体の高さｚとしきい値Tの関係ｚ＝Ｚ（T）が、２次元マスクパターンの構造体において成立する保証が必ずしもないことである。一般に、高さと信号強度の関係は一意的ではない。しかしながら、ライン幅等のパターン特性が近いパターンに対しては、本方法は、比較的よい近似を与えることができる。
【０１１７】
次に、第２の方法について説明する。本方法は、同一のマスク上の異なる複数の位置に同一の設計パターンが存在する場合に適用することができる。まず、情報処理装置1812は、マスク上の複数の同一２次元パターンについて形成された、ウエハ上の複数の２次元パターンを同一条件でSEM観察し、複数のSEM画像（信号強度分布）を得る。
【０１１８】
これら複数のSEM画像に対し、情報処理装置1812は、あるしきい値Tで切った時のエッジ座標列を求め、各座標点を結んだ等高線を考える。この際、計測精度を確保するため、エッジの抽出は、等高線に略垂直な方向に沿った信号強度分布を求め、これに対してしきい値を適用することにより求めることが望ましい。
【０１１９】
次に、情報処理装置1812は、上記複数の画像に対して求めた等高線間の距離偏差の総和が最小となるように、各画像を平行方向にシフトする。シフト位置の基準は、適宜設定してよい。さらに、情報処理装置1812は、平行シフト後の各画像の等高線に対して平均形状を求め、各画像の等高線に対して等高線の接線と垂直な方向に対する平均形状との差分（距離）を算出する。
【０１２０】
情報処理装置1812は、各エッジ点に対し、差分の分布の標準偏差を求め、これを実施例１におけるσｂ（Ｔ）に相当する値と仮定する。これにより、実施例１と同様、式５、式６、式７、式８等を用いて求めた傾斜角θを、当該位置におけるしきい値Tに対応する高さの側壁角度θとする。しきい値Tを変えてエッジを抽出した方向に沿って、式９と同様の積分を行うことにより、上記方向に沿った縦構造を求めることができる。
【０１２１】
第２の方法で求めたσｂ（Ｔ）は、真のＬＥＲに相当する成分を含むことが懸念されるが、平行シフトを行うことにより、比較的長周期のＬＥＲに起因するエッジシフトは相殺することができると期待される。さらに、真のＬＥＲの短周期成分と計測誤算成分は、上記σｂ（Ｔ）から適宜除去することが望ましい。
【０１２２】
本実施例の場合にも、実施例１の場合と同様に、σｓとσｅは、物理的な意味を持つ値として、各々独立に計測して求めても良いし、フィッティングパラメータとして考えても良い。発明者は、このような手法により推定した構造体の断面形状が、同じ構造体を切断してSEM画像として観察した断面形状とよく一致することを確認した。
【０１２３】
［実施例３］
本実施例では、本発明を、半導体集積回路等の製造に用いられる光リソグラフィによるレジストパターン形成プロセス、又はそれを用いて形成した半導体集積回路の立体的構造の品質評価、並びに上記プロセスモニターに適用する例について説明する。
【０１２４】
本実施例では、形成されたパターンの品質指標として、推定したパターンの高さと設計上のパターンの高さの偏差の絶対値を評価領域の全面について積分した値、又は、その値の２乗平均値を用いる（図１７の上段図（ａ）の右図にハッチングで示す部分の面積に対応する量)。
【０１２５】
なお、より単純に、前述の実施例の手法を適用して推定した断面形状（例えば図１７の上段図（ａ）の左図における曲線下部）の断面積Spと理想的設計形状（例えば図１７の上段図（ａ）の左図における長方形部分）の断面積Siの比Sp／Si（１次元マスクパターンの場合）、又は、前述の実施例の手法により推定した形状の体積Vpと理想的設計形状の体積Viの比Vp／Vi（２次元マスクパターンの場合）を、パターン品質指標として用いても良い。ただし、この場合には、設計上のパターンが本来存在しないはずの領域に存在するパターンの面積又は体積を負の値で扱うことが好ましい。これにより、設計上のパターンが本来存在しないはずの領域にパターンが存在する場合にも、品質劣化として検知することができる。
【０１２６】
２次元マスクパターンに対する理想的な設計形状には、回路パターンの設計形状、上記回路設計形状に対して光学的近接効果補正等を施したマスクパターンに対する光学シミュレーションの結果、様々な実際の露光結果に対して前述の実施例で説明した手法を適用して推定した体積のうちの最大値等を用いることができる。
【０１２７】
パターン形成用露光装置のフォーカス及び露光量条件をウエハ上の露光ショット毎にマトリクス状に様々に変化させて露光したウエハに対し、上記露光ショット毎にマスク上の同一パターンに対するレジストパターンをCD-SEMで観察し、その観察像からショットごとにパターン寸法を測定すると共に、実施例１に記載した方法によりパターン断面形状を推定し、上記パターン品質指標を求めた。
【０１２８】
図１７の中段図（ｂ）に示すように、パターンの寸法は、露光量に対して単調減少したが、フォーカス設定値に対するパターン寸法変化は小さく、パターン寸法変化からフォーカス設定値を推定するのは難しかった。これに対して、図１７の下段図（ｃ）に示すように、パターン品質指標は、露光量及びフォーカスの両方の設定値に対してほぼ単調に変化し、フォーカス設定値のモニタとして適用できた。ただし、パターン品質指標は、露光量に対しても変化する。そこで、まずパターン寸法変化から露光量を推定し、上記露光量に対して、パターン品質指標を用いてフォーカス設定値を推定することが望ましい。
【０１２９】
本実施例のように、推定された断面形状を、形成されたパターンの品質指標に用いれば、簡便かつ高速に最適な製造条件を判定したり、品質低下を検知することができる。また、当該判定結果等を製造プロセスにフィードバックすれば、形成されるパターンの品質劣化を抑制し、半導体装置をはじめとする各種素子の性能を向上させることができる。
【０１３０】
［実施例４］
本実施例では、σｘ、σｓを求める別の方法として、ｘｙ方向の２次元信号強度分布を周波数分析する方法について説明する。すなわち、第２の実施形態による２次源立体形状の推定方法について説明する。
【０１３１】
ここでは、図１８の上段図（ａ）のような立体構造体の側壁の全面に、表面に沿って等方的な表面凹凸が存在すると仮定し、その凹凸情報のみを、立体構造体の上方から観察する場合を想定する。このとき、図１８の下段図（ｂ）に示すパターンが凹凸情報として取得されると考えられる。
【０１３２】
ここで、凹凸のｘ方向の周期をＬｘとし、ｙ方向の周期をＬｙとすると、傾斜角θは、次式により求められる。
【数１７】
2014130077.tif 000019
【０１３３】
このとき、観察される２次電子信号の２次元強度分布は、第１にｘ方向のマクロな強度分布、第２にｙ方向の真のＬＥＲによる強度変動、第３に表面凹凸による局所的な強度変動という、３つの変動要因を含む。従って、観察像から第１及び第２の変動要因を除去することにより、第３の表面凹凸情報のみを抽出することができる。具体的な手順の１例を図１９を用いて説明する。
【０１３４】
まず、情報処理装置1812は、SEM像を取得した後、解析領域を指定し、上記領域内で２次元強度分布を取得する（ステップＳ１９０１）。次に、情報処理装置1812は、取得された２次元強度分布に対して２次元空間周波数フィルタリングを適用し、第１の変動成分による部分及び第２の変動成分の高周波成分を除いた部分の和と、第２の変動成分の高周波成分と第３の変動要因による部分の和に分解する（ステップＳ１９０２）。後者は、真のＬＥＲによる変動の高周波成分と表面凹凸による変動を含む。
【０１３５】
一般に、構造体の側壁角θはｘ方向に変化するので、強度分布変動のｘ方向空間周期はｘ方向に変化する。従って、空間周期は、ｘ方向に対して局所的に評価することが望ましい。このような解析を行う一般的な手法として、例えばウエーブレット解析がある。情報処理装置1812は、ｘ方向にウエーブレット解析を行うことにより、ｘ方向に沿った空間周波数特性の変化を検出する（ステップＳ１９０３）。
【０１３６】
この方法により得られた空間周波数特性（パワースペクトル密度PSDの空間周波数依存性）のｘ依存性の一例を図２０に示す。空間周波数特性は、ｘに依存して変化し、特性にピークが得られる場合にはピーク位置の周波数を代表空間周波数とし、ピークが明確でない場合には空間周波数の分布広がり（例えば半値幅）を代表空間周波数とし、その逆数によりｘ方向の空間周期σｘを求める。同様に、情報処理装置1812は、ｙ方向に空間周波数解析を行うことにより求めた代表空間周波数からσｙを求める。更に、情報処理装置1812は、ｘｙ両方向に２次元のウエーブレット解析を行ってもよい。なお、この過程において、必要に応じ、情報処理装置1812は、第２の変動成分の高周波成分を推定し、後者から除去する。このようにして、情報処理装置1812は、σｘとσｙの２次元分布を求める。
【０１３７】
情報処理装置1812は、側壁角の２次元分布θ（ｘ、ｙ）を次式より求める（ステップＳ１９０４）。
【数１８】
2014130077.tif 000020
【０１３８】
さらに、情報処理装置1812は、上式を位置（ｘ、ｙ）で積分することにより２次元立体形状を推定する（ステップＳ１９０５）。
【０１３９】
［実施例５］
本実施例では、構造体の表面に形成された凹凸パターンが必ずしも等方的でない場合について説明する。構造体表面の凹凸パターンが等方的であるという前提が成り立たないケースには、例えばレジスト構造体の側壁形状が、レジスト膜内における露光光の干渉による定在波の影響を受ける場合がある。この場合、２次電子信号強度分布は、図２１に示すように、構造体の１本のエッジに対して、平行する複数のエッジとして観察される。この場合、複数エッジ間の側壁傾斜角度は、エッジ間隔をΔＬxとして、次式で表される。
【０１４０】
【数１９】
2014130077.tif 000021
【０１４１】
ただし、λは当該レジストパターンの露光に用いられた光の波長であり、ｎは上記波長の光に対する上記レジスト材料の屈折率である。定在波は、レジスト露光時にレジスト下地膜からの反射が大きいと出現するが、この場合、レジスト膜厚の変動によりレジスト寸法が敏感に変動する。これは、実際の生産工程では望ましいことではないので、通常これを抑制するために十分な反射防止策がなされる。このため、実際には上記のような定在波が影響することは少ない。
【０１４２】
表面凹凸が等方的である前提が必ずしも成立しない別のケースとして、等方的な表面凹凸を有するレジストパターンをエッチングして被エッチング構造を転写する場合がある。レジストパターン表面の凹凸は、エッチングにより被エッチング膜の構造体表面に転写されるが、エッチング転写後の構造体のラフネス解析によれば、被エッチング構造体のラフネスは、レジスト構造体のラフネスをほぼ反映したものであることが多い。
【０１４３】
この場合、図２２に示すように、凹凸のエッジ方向（ｙ方向）の周期は変わらないが、縦方向の寸法は変化する。ここで、レジスト材料のエッチング速度Vrと被エッチング膜のエッチング速度Veの比をＲetch＝Ve/Vrとして、エッチングが完全に異方的であるとすると、凹凸の縦方向寸法はおよそＲetch倍となると考えられる。この場合、側壁角度が比較的急峻であると仮定すると、式５におけるLもＲetch倍すればよい。
【０１４４】
実際には、エッチングはエッジ方向と垂直な方向（ｘ方向）にも進行する（エッチング）。この場合、凹凸の縦横寸法比は必ずしもＲetchに等しくない。従って、例えば他の方法により観察した断面形状と、前述した手法により推定した断面形状がよく一致するように、式５におけるＬをフィッティングパラメータとして最適化する等することが好ましい。
【０１４５】
なお、本発明の最も一般的な概念は、構造体を上面から観察した２次元画像の信号強度分布の局所的なばらつきの情報から、構造体の立体形状を推定することである。従って、推定アルゴリズムは上記に示した方法に限らない。
【０１４６】
（他の実施例）
なお、本発明は上述した実施例に限定されるものでなく、様々な変形例を含んでいる。例えば上述した実施例は、本発明を分かりやすく説明するために、一部の実施例について詳細に説明したものであり、必ずしも説明した全ての構成を備える必要は無い。また、ある実施例の一部を他の実施例の構成に置き換えることが可能であり、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成を追加、削除又は置換することも可能である。
【０１４７】
また、上述した各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路その他のハードウェアとして実現しても良い。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することにより実現しても良い。すなわち、ソフトウェアとして実現しても良い。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリやハードディスク、SSD（Solid State Drive）等の記憶装置、ICカード、SDカード、DVD等の記憶媒体に格納することができる。
【０１４８】
また、制御線や情報線は、説明上必要と考えられるものを示すものであり、製品上必要な全ての制御線や情報線を表すものでない。実際にはほとんど全ての構成が相互に接続されていると考えて良い。
【符号の説明】
【０１４９】
1801…走査型電子線顕微鏡の筐体
1802…電子銃
1803…電子線
1804…収束レンズ
1805…偏向器
1806…対物レンズ
1807…観察ウエハ
1808…ステージ
1809…二次電子
1810…検出器
1811…制御系
1812…情報処理装置
1813…データ記憶装置
【図１】
2014130077.tif 000022
【図２】
2014130077.tif 000023
【図３】
2014130077.tif 000024
【図４】
2014130077.tif 000025
【図５】
2014130077.tif 000026
【図６】
2014130077.tif 000027
【図７】
2014130077.tif 000028
【図８】
2014130077.tif 000029
【図９】
2014130077.tif 000030
【図１０】
2014130077.tif 000031
【図１１】
2014130077.tif 000032
【図１２】
2014130077.tif 000033
【図１３】
2014130077.tif 000034
【図１４】
2014130077.tif 000035
【図１５】
2014130077.tif 000036
【図１６】
2014130077.tif 000037
【図１７】
2014130077.tif 000038
【図１８】
2014130077.tif 000039
【図１９】
2014130077.tif 000040
【図２０】
2014130077.tif 000041
【図２１】
2014130077.tif 000042
【図２２】
2014130077.tif 000043