(57)【要約】【課題】    粒子・空孔の絶対量が少ないときであっても、高精度で粒子・空孔サイズ分布を解析し、かつ表面界面形状を評価する膜構造解析方法とその装置を提供すること。【解決手段】    単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析方法において、膜試料表面に対する入射角度と出射角度とが相関関係を有さないように入射角度と出射角度との組を複数設定し、それぞれの入射角度と出射角度との組についてＸ線強度を測定することにより得られる実測Ｘ線散乱曲線に対してフィッティングを行うことを特徴とする。【選択図】図１
2005121613.tif 000002
【特許請求の範囲】
【請求項１】
単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析方法において、膜試料表面に対する入射角度と出射角度とが相関関係を有さないように入射角度と出射角度との組を複数設定し、それぞれの入射角度と出射角度との組についてＸ線強度を測定することにより得られる実測Ｘ線散乱曲線に対してフィッティングを行うことを特徴とする膜構造解析方法。
【請求項２】
単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析方法において、オフセットスキャン測定、ロッキングスキャン測定、および、ディテクタースキャン測定の内の少なくとも２つの測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行うことを特徴とする膜構造解析方法。
【請求項３】
オフセットスキャン測定により得られる実測Ｘ線散乱曲線に対し、膜試料表面の界面粗さおよび膜試料内部の空孔サイズに関するパラメータを変化させて行うシミュレーション演算により得られるシミュレートＸ線散乱曲線をフィッティングして実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定することを特徴とする請求項２記載の膜構造解析方法。
【請求項４】
単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析装置において、膜試料表面に対する入射角度と出射角度とが相関関係を有さないように入射角度と出射角度との組を複数設定し、それぞれの入射角度と出射角度との組についてＸ線強度を測定することにより得られる実測Ｘ線散乱曲線に対してフィッティングを行うことを特徴とする膜構造解析装置。
【請求項５】
単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析装置において、オフセットスキャン測定、ロッキングスキャン測定、および、ディテクタースキャン測定の内の少なくとも２つの測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行うことを特徴とする膜構造解析装置。
【請求項６】
オフセットスキャン測定により得られる実測Ｘ線散乱曲線に対し、膜試料表面の界面粗さおよび膜試料内部の空孔サイズに関するパラメータを変化させて行うシミュレーション演算により得られるシミュレートＸ線散乱曲線をフィッティングして実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定することを特徴とする請求項５記載の膜構造解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この出願の発明は、膜構造解析方法およびその装置に関するものである。さらに詳しくは、この出願の発明は、単層膜あるいは多層膜における界面粗さおよび空孔分布について解析を行う膜構造解析方法およびその装置に関するものである。
【背景技術】
【０００２】
薄膜試料の表面・界面の構造解析方法として、電子顕微鏡などのマイクロスコープ技術が一般に用いられており、これらにより実空間観測という点で多くの情報を得ることができる。しかしながら、観測領域の制限や試料へのダメージ、また、埋もれた膜中における界面構造解析の困難などの問題点を抱えていることも事実である。
【０００３】
一方で、非破壊での膜構造解析方法としてＸ線反射率測定やＸ線散漫散乱測定が用いられている。
【０００４】
これまで、ポーラス膜等の密度不均一試料内の粒径分布についてＸ線反射率測定を応用して解析する方法が、この出願の発明の発明者らにより提案されている（特許文献１参照）。この解析方法は、Ｘ線の散漫散乱強度を測定し、その測定値に基づいて粒径分布を解析するものである。
【０００５】
Ｘ線を薄膜表面すれすれに入射すると、入射Ｘ線は反射および透過する。このとき、すれすれの入射角度が薄膜の全反射臨界角度近傍であれば、屈折の影響により（透過）Ｘ線は薄膜内部にほとんど進入できず、エバネッセント波として表面から数ｎｍの領域内において、薄膜表面へ移行方向に進行する。このように、Ｘ線入射角度が、全反射臨界角度に近いすれすれの入射条件であるとき、薄膜の表面構造解析が可能となる。
【０００６】
薄膜の表面界面は、試料表面法線方向におよそ数ｎｍのオーダーの電子密度揺らぎ（ラフネス）を持つ。そのため、上記の通り、試料表面にすれすれにＸ線を入射することにより、試料の表面界面の形状に対応したＸ線プロファイルを観測し、さらに、そのプロファイルを解析することで表面・界面形状を評価することを実現する密度不均一試料解析方法ならびにその装置およびシステムが提案されている（特許文献２参照）。
【特許文献１】特開２００１−３４９８４９号公報
【特許文献２】特開２００３−１４６６３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
これまでこの出願の発明者らによって提案されてきたＸ線を用いた密度不均一試料の解析方法は、優れた解析能力を実現するものであるが、試料の界面の形状と空孔サイズの同時解析が不可能であり、また、解析対象が単層膜である場合を想定された技術であることから、多層膜からなる試料に対しては適用することが困難であった。また、薄膜中の粒子・空孔の絶対量が少なくなると（単位体積あたりの粒子・空孔数や密度が小さくなると）、粒子・空孔等の散乱に対して、相対的に表面・界面粗さの散乱が顕著になる。表面・界面粗さによる散乱は、その電子密度が界面にのみ存在し、膜厚には依存しないからであるという問題がある。
【０００８】
この出願の発明は、以上のとおりの事情に鑑みてなされたものであり、試料の界面の形状と空孔サイズの同時解析を実現し、また、多層からなる試料に対しても適用可能であり、Ｘ線入射角度が全反射臨界角度より大きくなったときに粒子・空孔散乱が急激に大きくなる点、Ｘ線入射角度が０°から全反射臨界角度にかけて表面・界面散乱が顕著に観測される点という２つの特徴を用いることにより、粒子・空孔の絶対量が少ないときであっても、高精度で粒子・空孔サイズ分布を解析し、かつ表面界面形状を評価する膜構造解析方法とその装置を提供することを課題としている。
【課題を解決するための手段】
【０００９】
この出願の発明は、上記の課題を解決するものとして、第１に、〔請求項１〕
単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析方法において、膜試料表面に対する入射角度と出射角度とが相関関係を有さないように入射角度と出射角度との組を複数設定し、それぞれの入射角度と出射角度との組についてＸ線強度を測定することにより得られる実測Ｘ線散乱曲線に対してフィッティングを行うことを特徴とする膜構造解析方法を提供する。
【００１０】
また、この出願の発明は、第２に、〔請求項２〕
単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析方法において、オフセットスキャン測定、ロッキングスキャン測定、および、ディテクタースキャン測定の内の少なくとも２つの測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行うことを特徴とする膜構造解析方法を提供する。
【００１１】
この出願の発明における第３の態様として、上記の膜構造解析方法は、オフセットスキャン測定により得られる実測Ｘ線散乱曲線に対し、膜試料表面の界面粗さおよび膜試料内部の空孔サイズに関するパラメータを変化させて行うシミュレーション演算により得られるシミュレートＸ線散乱曲線をフィッティングして実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定することを特徴とする。
【００１２】
そして、この出願の発明は、第４に、単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析方法において、膜試料表面に対する入射角度と出射角度とが相関関係を有さないように入射角度と出射角度との組を複数設定し、それぞれの入射角度と出射角度との組についてＸ線強度を測定することにより得られる実測Ｘ線散乱曲線に対してフィッティングを行うことを特徴とする膜構造解析装置を、また、第５に、単層膜または多層膜からなる膜試料の表面に対する臨界角近傍の角度からのＸ線入射により得られる実測Ｘ線散乱強度曲線に対し、シミュレーション演算により得られるシミュレートＸ線散乱曲線を、少なくとも１個以上の試料の物性を示すパラメータを変化させることによりフィッティングを行い、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する膜構造解析方法において、オフセットスキャン測定、ロッキングスキャン測定、および、ディテクタースキャン測定の内の少なくとも２つの測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行うことを特徴とする膜構造解析装置を提供する。
【００１３】
上記の膜構造解析装置においては、この出願の発明の第６の態様として、オフセットスキャン測定により得られる実測Ｘ線散乱曲線に対し、膜試料表面の界面粗さおよび膜試料内部の空孔サイズに関するパラメータを変化させて行うシミュレーション演算により得られるシミュレートＸ線散乱曲線をフィッティングして実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定することを特徴とするものである。
【発明の効果】
【００１４】
この出願の発明の膜構造解析方法およびその装置により、試料の界面の形状と空孔サイズの同時解析が実現され、また、多層からなる試料に対しても適用可能となる。
【００１５】
以下に、この出願の発明である膜構造解析方法およびその装置による発明の効果を、具体的に示す。
１）  従来は、オフセットスキャン測定により得られる実測Ｘ線散乱曲線に対するフィッティングにより膜構造解析が行われていたが、薄膜試料表面に対する入射角度と出射角度とが相関関係を有さないように入射角度と出射角度との組を複数設定し、それぞれの入射角度と出射角度との組についてＸ線強度を測定することにより得られる実測Ｘ線散乱曲線に対してフィティングを行ったり、従来のスキャン測定を応用したオフセットスキャン測定、ロッキングスキャン測定、およびディテクタースキャン測定の内の少なくとも２つの測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィティングを行うことにより薄膜試料の物性を示すパラメータの最適値を求めることができる。この出願の発明である膜構造解析方法およびその装置においては、ロッキングスキャン測定により得られる実測Ｘ線散乱曲線を、オフセットスキャン測定により得られる実測Ｘ線散乱曲線と併せて解析を行うことにより、界面の形状と薄膜中の空孔サイズ分布解析を行うことが可能となる。
２）  この出願の発明である膜構造解析方法においては、多層膜中の空孔サイズ分布の評価を行うことが可能となる．
３）  ロッキングスキャン測定により得られる実測Ｘ線散乱曲線において、入射角度及び反射角度が全反射臨界角度より小さい領域においては、最表面界面からの散乱のみのシグナルが含まれており、複数種類のロッキングスキャン測定により得られる実測Ｘ線散乱曲線に対し、これらの角度領域を解析することにより最表面界面の形状を評価することが可能となる。
４）  界面散乱は、空孔散乱とは異なり、各界面間で生じる散乱が干渉し、膜厚に相当する振動が観測される。これらの振幅および干渉散乱の傾斜から表面以外の界面形状を評価することが可能となる。
５）  特に、界面散乱と空孔散乱とではロッキングスキャン測定により得られる実測Ｘ線散乱曲線に大きな違いが見られる。界面散乱においては、界面面内方向の電子密度揺らぎ相関距離および界面の形状に対応して鏡面反射の周りに強く散乱が観測される。一方、空孔散乱においては、鏡面反射近傍においてもほぼ一定の散乱強度が観測される。これらの違いを利用することにより空孔散乱と界面散乱の比率を推定することが可能となり、両散乱を同時に解析することが可能となる。これは、界面散乱が、面内と表面直行方向で形状に異方性があることに起因する。空孔散乱は、膜中にランダムに存在しているために、このような異方性はほとんど見られない。
６）  ロッキングスキャン測定において、空孔散乱は散乱体（空孔）が薄膜中に存在しているためＸ線が膜中に侵入できる条件、すなわち、全反射臨界角度を越える条件における観測が可能である。言い換えるならば、Ｘ線入射角度がポーラス膜の全反射臨界角度を越えると、急激に散乱が観測されることとなる。一方、表面界面ラフネスは、散乱体（ラフネス）が界面にあるため、界面散乱は界面における屈折Ｘ線の強度に依存し、Ｘ線入射角度およびＸ線反射角度が全反射臨界角度より小さい場合であっても測定される（Ｙｏｎｅｄａ−Ｗｉｎｇ）。また、屈折Ｘ線（Ｙｏｎｅｄａ−Ｗｉｎｇ）の入射角度依存性は連続的であり、急激に増加するものではない。つまり、ロッキングスキャン測定において、ポーラス膜の全反射臨界角度において散乱強度が急増しているか否か、または、散乱強度の増加量に着目することにより、界面散乱と空孔散乱とを区別することが可能となる。
【発明を実施するための最良の形態】
【００１６】
この出願の発明である膜構造解析方法においては、膜試料の表面に対し臨界角近傍の角度からのＸ線を入射することにより得られる実測Ｘ線散乱強度曲線に対し、少なくとも１個以上の試料の物性を示すパラメータを変化させて行うシミュレーション演算により得られるシミュレートＸ線散乱曲線をフィッティングして実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する。この出願の発明である膜構造解析方法においては、膜試料は単層膜または多層膜のどちらであってもよい。
【００１７】
この出願の発明である膜構造解析方法においては、膜試料表面に対するＸ線入射角度とＸ線出射角度とが相関関係を有さないように入射角度と出射角度との組を複数設定し、それぞれの入射角度と出射角度との組についてＸ線強度を測定することにより得られる実測Ｘ線散乱曲線に対してフィッティングを行う。すなわち、オフセットスキャン測定、ロッキングスキャン測定、および、ディテクタースキャン測定においては、膜試料表面に対するＸ線入射角度とＸ線出射角度とが特定の相関関係を有しているが、この出願の発明である膜構造解析方法においては、このようなＸ線入射角度とＸ線出射角度とが相関関係を有するスキャン測定を行わず、別個に任意に定められたＸ線入射角度とＸ線出射角度との組を複数設定し、それらの組について測定を行うことにより得られたＸ線強度に対してフィッティングを実施し、膜構造の解析を行うものである。
【００１８】
ディテクタースキャンは、薄膜試料表面に対するＸ線の入射角度を一定にし（試料回転軸θを固定）、検出器の回転軸２０をスキャンする方法で、Ｘ線入射角度を複数変化させたディテクタースキャンにより、すべての入射角度、出射角度のＸ線強度を測定することができる。
【００１９】
従来のスキャン測定を応用した本発明の膜構造解析方法としては、以下に示す４種類が知られている。
１）オフセットスキャン測定およびロッキングスキャン測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行い、併せて解析を行う方法
２）ロッキングスキャン測定およびディテクタースキャン測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行い、併せて解析を行う方法
３）オフセットスキャン測定およびディテクタースキャン測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行い、併せて解析を行う方法
４）オフセットスキャン測定とロッキングスキャン測定とディテクタースキャン測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングを行い併せて解析を行う方法
この出願の発明である膜構造解析方法においては、以上に示したＸ線入射角度とＸ線出射角度とが相関関係を持つスキャン方法を用いて測定した複数のＸ線強度のみによりフィッティングを行い、膜構造の解析をするのではなく、薄膜試料表面に対するＸ線入射角度とＸ線出射角度とを、全く別個に、かつ任意に決定し、複数のＸ線入射角度およびＸ線出射角度の組により測定された複数のＸ線強度を用いてフィッティングを行い、膜構造の解析を実施する。
【００２０】
もちろん、実測Ｘ線散乱強度曲線を取得する際には、オフセットスキャン測定、ロッキングスキャン測定、および、ディテクタースキャン測定の内の少なくとも２つの測定が行われ、それぞれの測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングが行ってもよい。
【００２１】
具体的な、フィッティングの手順としては、例えば、最小二乗法による解析が用いられ、Ｘ線反射率を求めている因子をパラメータ（変数）として、パラメータの値を少しずつ変化させながら反射率を計算して、実際の反射率データとの残差二乗和を最小にすることにより、測定データに最もフィッティングするような１組のパラメータを決定する。
【００２２】
最小二乗法による解析では、フィッティングの順番や方法、あるいは初期値の適正や不適正などによって、単純に残差二乗和を最小にすることが困難であり、発散してしまいフィッティングができないことが少なくない。通常の解析では、Ｘ線反射率の測定データへのシミュレーション演算データのフィッティングは以下の順番で実行することにより、解が収束することが知られている。
１）  基板を含む薄膜の層構造のモデルを作成する。
２）  測定データとシミュレーション演算データに対して、臨界角より低角度での最大強度を一致させる。
３）  測定データの強度とモデルのバックグランドの強度を一致させる。
４）  モデルのパラメータの中で固定するものを決定する。多くの場合は、基板ガラスは既知であるので、その密度を固定する。また、膜厚と密度などのその他のパラメータで既知のものがあれば定数化しておく。
５）  目視のレベルで測定データとシミュレーション演算データが一致するように、パラメータを手動で決定する。
６）  自動のフィッティングを実行する。
７）  測定データとシミュレーション演算データとの差異が最小となったときのパラメータの値を最適値として求められる。このとき、最小二乗和の残差が少数２桁以下（１０-2のオーダ）になれば、かなりの精度でフィッティングできていると判断される。
【００２３】
フィッティングにおいて実行されるシミュレーション演算の原理について以下に説明する。
量子論的散乱の記述
試料表面に対してすれすれにＸ線を入射した場合、入射したＸ線および界面で反射したＸ線は、互いに干渉するため、干渉縞が観測される。Ｘ線反射率測定においては、この干渉現象を利用し、薄膜の密度だけでなく、薄膜の膜厚および界面のラフネスを評価するのに利用されている。この出願の発明である膜構造解析方法においては、図１に示すように、試料の表面・界面における散乱は、すれすれのＸ線入射の条件でＸ線が試料の表面・界面の粗さによって散乱される場合を取り扱うため、散乱確率が反射・屈折条件を満足するＸ線の基底状態に、表面・界面粗さの散乱ポテンシャルを作用させた場合の散乱確率を計算することになる。量子力学の散乱理論によれば、入射Ｘ線の固有状態をΨi、散乱Ｘ線の固有状態をΨfと、試料の表面・界面の散乱ポテンシャルをＶとした場合、微分散乱面積は次式で与えられる。
【００２４】
【数１】
2005121613.tif 000003
【００２５】
式（Ｉ）に示したとおり、散乱を記述するためには、まず、入射Ｘ線の固有状態、散乱Ｘ線の固有状態、および、散乱ポテンシャルが必要となる。
薄膜中のＸ線の固有状態
多層膜中のＸ線の固有状態Ψの計算は、Ｐａｒｒａｔｔにより既に行われている。Ｆｒｅｎｅｌｌの式を用いると、第ｌ層の薄膜の深さＺiにおける入射Ｘ線の固有状態は次式にて与えられる。
【００２６】
【数２】
2005121613.tif 000004
【００２７】
また、散乱Ｘ線については、試料の表面・界面において散乱されたときの出射波について考える必要がある。そのような波動方程式の解は、入力Ｘ線の解（固有状態）を時間反転したものが考えられる。これは入射Ｘ線の解の複素共役を取り、さらにｋ→−ｋとすることで、次式のように得られるものである。
【００２８】
【数３】
2005121613.tif 000005
【００２９】
薄膜中の表面・界面散乱過程
薄膜中の多重反射内で生じる表面・界面散乱は、複雑な散乱過程を持つ。あるものは薄膜の全反射臨界角度近傍のみで強く観測されたり、また、あるものは広い角度領域で強く観測されたりと、様々な散乱過程が見られる。図２に、各層において観測される４成分の散乱過程を示す。
【００３０】
図２（ａ）は、界面ｌ−１において屈折したＸ線が散乱され、そのまま観測された場合を示しており、以下、散乱過程（ａ）と呼ぶ。図２（ｂ）は、界面ｌにおいて鏡面反射したＸ線が散乱され観測された場合を示しており、以下、散乱過程（ｂ）と呼ぶ。図２（ｃ）は、散乱されたＸ線が界面ｌで鏡面反射し観測された場合を示しており、以下、散乱過程（ｃ）と呼ぶ。図２（ｄ）は、界面ｌで鏡面反射したＸ線が散乱され、散乱されたＸ線が再び界面ｌで鏡面反射し観測された場合を示しており、以下、散乱過程（ｄ）と呼ぶ。図２（ａ）で示した場合を除き、図２（ｂ）〜（ｄ）で示した散乱過程は、界面ｌ層での反射の効果を含んでいるため、界面ｌが全反射を満足しているときのみ強く観測される。一方、図２（ａ）は、界面ｌ層の反射の寄与がないため、広い角度領域において観測される。
薄膜中の表面・界面散乱過程に依存する散乱ベクトル
図２に示したとおり、表面／界面散乱で生じる表面法線方向の散乱ベクトルｑsは２種類存在する。以下に、界面ｌ−１で生じる表面法線方向の散乱ベクトルの大きさを記述する。
散乱過程（ａ）および散乱過程（ｄ）の表面法線方向の散乱ベクトルの大きさ
【００３１】
【数４】
2005121613.tif 000006
【００３２】
散乱過程（ｂ）および散乱過程（ｃ）の表面法線方向の散乱ベクトルの大きさ
【００３３】
【数５】
2005121613.tif 000007
【００３４】
一方、面内方向の散乱ベクトルの大きさｑrは、各界面および散乱過程に依存せず、次式で与えられる。
【００３５】
【数６】
2005121613.tif 000008
【００３６】
微分散乱断面積の記述
第ｌ層における表面・界面の散乱ポテンシャルをＶlとした場合、図２に示した４つの散乱過程を考慮すれば、第ｌ層における散乱振幅は次式で与えられる。
【００３７】
【数７】
2005121613.tif 000009
【００３８】
全散乱振幅は、上記全ての層における散乱振幅の輪を取ればよく、次式により与えられる。
【００３９】
【数８】
2005121613.tif 000010
【００４０】
式（I）および式（VIII）より、散乱断面積は以下の式で求められる。
【００４１】
【数９】
2005121613.tif 000011
【００４２】
式（IX）の第２項目は、各界面間で生じる干渉性散乱に相当する。
表面・界面ラフネスによる散乱
理想的にフラットな界面状態に対して界面にラフネスが存在すると、界面での電子密度変化と界面ラフネスの形状に依存した散乱が観測される。散乱前、散乱後のＸ線の固有状態は、それぞれ式（II）または式（III）により厳密に与えられる。各界面でのラフネスによる散乱振幅は、式（VII）を用いて計算することができる。このとき、図２に示した散乱過程において、散乱ポテンシャルＶlは、次式で与えられる。
散乱過程（ａ）または散乱過程（ｃ）の場合
【００４３】
【数１０】
2005121613.tif 000012
【００４４】
散乱過程（ｂ）または散乱過程（ｄ）の場合
【００４５】
【数１１】
2005121613.tif 000013
【００４６】
ここで、ｎは各界面での屈折率であり、また、ｋ0は真空での波数ベクトルである。
【００４７】
この出願の発明で記述されるべき散乱式は、式（IX）の微分散乱断面積である。微分散乱断面積を記述するためには、各界面の各散乱過程の散乱振幅を求める必要がある。散乱振幅を記述できれば、微分散乱断面積は記述できたことに等しい。
【００４８】
まず、最初に界面ｌの散乱過程（ａ）において生じる散乱振幅の導出について、以下に示す。
【００４９】
【数１２】
2005121613.tif 000014
【００５０】
ここでは、表面散乱の計算を行うため、体積成分を表面積分に置き換える。
【００５１】
【数１３】
2005121613.tif 000015
【００５２】
ここで、ｄＳは表面に垂直な微分表面ベクトルである。また、Ａは任意の単位ベクトルである。ここでは、Ａを試料表面法線方向にとるものとする。そうすると、式（XIII）は次式のように変形可能である。
【００５３】
【数１４】
2005121613.tif 000016
【００５４】
さて、式（XIV）において表面が特定の形状を持っている場合については、計算を実行できるが、多くの場合、表面・界面の電子密度揺らぎが、ある相関関数に従うものと想定し、｜＜−ζl*｜Ｖl｜ηl＞｜2を計算することにする（Ｓｈｉｎｂａ、Ｐｈｙｓ。Ｒｅｖ。Ｂ３８、ｐ。２２９７  （１９９８）に詳しい）。
【００５５】
【数１５】
2005121613.tif 000017
【００５６】
ここで、σは表面・界面粗さ、また、ζは表面・界面粗さの面内方向の相関距離である。ｈは表面ラフネスの形状を与えるＨｕｒｓｔパラメータであり、ｈが１に近いほど滑らかな表面を与え、また、ｈが０に近いほどぎざぎざな表面形状を与える。
【００５７】
表面がある特定の形状を持たないため、表面・界面粗さの＜−ζl*｜Ｖl｜ηl＞を計算することはできないが、式（XV）より｜＜−ζl*｜Ｖl｜ηl＞｜2を求めることは可能である。式（XIV）および式（XV）により、界面ｌの散乱過程（ａ）における散乱振幅は以下のように求まる。
【００５８】
【数１６】
2005121613.tif 000018
【００５９】
同様に、界面ｌにおける散乱過程（ｂ）、散乱過程（ｃ）および散乱過程（ｄ）の散乱振幅も以下のとおりに求められる。
散乱振幅（ｂ）
【００６０】
【数１７】
2005121613.tif 000019
【００６１】
散乱振幅（ｃ）
【００６２】
【数１８】
2005121613.tif 000020
【００６３】
散乱振幅（ｄ）
【００６４】
【数１９】
2005121613.tif 000021
【００６５】
以上で導出した散乱振幅を用い、また、式（VII）〜（IX）を用いれば、微分散乱断面積が与えられる。
【００６６】
この出願の発明である膜構造解析方法においては、以上の原理を用いて、Ｘ線の界面散乱をシミュレーション演算が行われる。
【００６７】
この出願の発明である膜構造解析方法においては、前述の通り、実測Ｘ線散乱強度曲線を取得する際には、オフセットスキャン測定、ロッキングスキャン測定、および、ディテクタースキャン測定の３通りの測定が行われ、それぞれの測定により得られる実測Ｘ線散乱曲線のそれぞれに対してフィッティングが行われる。
【００６８】
オフセットスキャン測定により得られる実測Ｘ線散乱曲線に対しては、膜試料表面の界面粗さおよび膜試料内部の空孔サイズに関するパラメータを変化させてシミュレーション演算が行われ、このシミュレーション演算により得られるシミュレートＸ線散乱曲線を実測Ｘ線散乱曲線にフィッティングして実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する。
【００６９】
また、ロッキングスキャン測定は、複数の異なるスキャン角度を設定して行われ、それぞれのロッキングスキャン測定により得られた実測Ｘ線散乱曲線に対してシミュレートＸ線散乱曲線をフィッティングすることで、実測Ｘ線散乱曲線とシミュレートＸ線散乱曲線との差異が最小となったときのパラメータの値を最適値として求め、膜試料の構造を決定する。このとき、入射側の膜試料表面の全反射臨界角度から出射側の膜試料表面の全反射臨界角度までの間の実測Ｘ線散乱曲線に対してシミュレートＸ線散乱曲線をフィッティングすることにより、膜試料内部の空孔サイズに関するパラメータの最適値を求めることが可能である。また、入射側の膜試料面の全反射臨界角度以下あるいは出射側の膜試料表面の全反射臨界角度以上の実測Ｘ線散乱曲線に対してシミュレートＸ線散乱曲線をフィッティングすることにより、膜試料表面の界面粗さに関するパラメータの最適値を求めることが可能である。
【００７０】
以下に、この出願の発明である膜構造解析方法を実施するための装置構成および手順について説明する。
【００７１】
図３は、この出願の発明の膜構造解析方法を実施するための膜構造解析装置における光学系について示した概要図である。入射側においては、光源（３１）から照射される入射Ｘ線は、多層膜ミラー（３２）により反射され、さらに２体のナロースリット（３３）によりビーム径の絞り込みが行われる。試料（３４）に入射されたＸ線は反射・屈折・散乱し、受光側に設置された２体のナロースリット（３５）を通過して、検出手段（３６）により受光される。ナロースリット（３５）を用いることにより、検出手段（３６）側の分解能を高めることが可能となる。
【００７２】
この出願の発明である膜構造解析装置は以上の構成を最低限備えるものであって、実際には、上記の光学系の他に、光学系により得られた測定データを取り込み、取り込まれた測定データを用いて膜構造解析を実行するコンピュータなどの膜構造解析手段を備えている。
【００７３】
さらに、この出願の発明の膜構造解析方法を実施するための具体的手順について詳細に説明する。
【００７４】
まず、Ｘ線の入射角度と反射角度を一定角度だけずらしてオフセットスキャン測定を行うことにより、オフセット角度を見積もる。このとき、オフセット角度は、鏡面反射の影響の無い角度とする。オフセット角度は、検出器ステージ角度２θを固定し、試料ステージ角度θを一定角度分だけスキャンすることにより行われる鏡面反射のロッキングスキャンにより知ることができる。ロッキングスキャン測定によるＸ線散乱曲線の鏡面反射の半価全幅は、受光側のスリット条件により計算することができる。受光側の受光スリット幅をＲＳＷと、また、散乱抑制スリット幅をＳＳＷとし、さらに、ＲＳＷがＳＳＷに等しいものとした場合、試料中心から散乱抑制スリットまでの距離をＬとすれば、おおよその受光側の発散角度は次式で表される。
【００７５】
【数２０】
2005121613.tif 000022
【００７６】
例えば、ＲＳＷ＝０．２ｍｍ、ＳＳＷ＝０．２ｍｍ、Ｌ＝３４０ｍｍである場合、式（XX）を用いると受光側分解能はおよそ０．０３°となる。実際には、薄膜試料表面の影響等があるため、図４に示すような鏡面反射のロッキングスキャンによる測定結果より、鏡面反射の影響がないオフセット角度を設定する必要がある。図４に示した場合であれば、オフセット角度は０．１°と見積もられ、このように設定することにより、鏡面反射の影響をほとんど受けないことがわかる。
【００７７】
次いで、上記の通りに見積もられたオフセット角度を一定に保ちながら、検出器ステージ角度２θと試料ステージ角度θを２θ：θ＝２：１となるように移動しながら、オフセットスキャン測定を行う。オフセットスキャン測定においては、オフセット角度を可能な限り小さな値に設定していることから、面内の散乱ベクトルは試料表面法線方向の散乱ベクトルと比較して極めて小さくなる。そのため、界面散乱の情報としては、試料表面法線方向の電子密度揺らぎの相関距離、すなわち、ラフネスの情報を多く含むことになる。空孔散乱についても同様であり、試料表面方向の空孔サイズの情報を持つ。ただし、薄膜中の空孔の形状が、界面形状と異なり、面内と面法線方向で大きな異方性がなければ、薄膜中の平均の空孔サイズ情報を持っているものと考えられる。
【００７８】
このように、オフセットスキャン測定のＸ線散乱曲線から、界面のラフネス情報および薄膜中の平均空孔サイズに関する情報を知ることが可能となる。また、図５に示すように、界面のラフネス（試料表面法線方向の電子密度揺らぎの相関距離）が大きくなるにつれて、オフセットスキャンのＸ線散乱曲線の傾斜が大きくなる。同様に、図６に示すように、薄膜中の平均空孔サイズが大きくなるにつれて、オフセットスキャンのＸ線散乱曲線の傾斜が大きくなる。
【００７９】
続くロッキングスキャン測定は、前述の通り、検出器ステージ角度２θを一定に固定したまま、試料ステージ角度θを変化させることにより行う。オフセット角度の見積もりを行うステップにおいては、鏡面反射スペクトル近傍のみで測定を行っているが、このステップにおけるロッキングスキャン測定においては、図７に示したように、試料ステージ角度θを０°から検出器ステージの固定角度までスキャンする。
【００８０】
ロッキングスキャン測定においては、検出器の角度を固定することから、試料表面法線方向の散乱ベクトルはほぼ一定となる。一方、試料表面面内方向の散乱ベクトルは大きく変化している。そのため、界面散乱の情報としては、試料表面面内方向の電子密度揺らぎ相関距離および形状パラメータの情報が多く含まれる。
【００８１】
図８は、形状パラメータとしてＨｕｒｓｔパラメータを変化させた場合のロッキングスキャン測定により得られるシミュレーションＸ線散乱曲線であり、Ｈｕｒｓｔパラメータを大きく設定することによりＸ線散乱曲線はガウス型の曲線となる。また、Ｈｕｒｓｔパラメータを小さく設定することにより、鏡面反射近傍で急激に散乱強度が増加していることがわかる。
【００８２】
図９は、面内相関距離を変化させた場合のロッキングスキャン測定により得られるシミュレーションＸ線散乱曲線である。面内相関距離を大きく設定した場合、鏡面反射からわずかにずれるだけで、すなわち、面内の散乱ベクトルが０からわずかに変化するだけで、は大きく変化する。この傾斜の大きさは、面内の電子密度揺らぎの相関距離に対応し、傾きが大きいほど相関距離が大きくなる。一方、面内相関距離が小さい場合、もっと面内の散乱ベクトルが大きく変化する角度に検出器ステージ角度２θの固定角度を設定する必要がある。すなわち、面内の相関距離が小さい場合、面内相関距離に対応する傾斜が緩やかであるため、面内の散乱ベクトルが大きく変化するような条件で、ロッキングスキャン測定を行わなければならない。面内の相関距離が比較的小さい場合に検出器ステージ角度２θの固定角度を様々に変化させたロッキングスキャン測定のシミュレーションＸ線散乱曲線を図１０に示す。
【００８３】
実際に測定されるＸ線散乱曲線には、空孔から散乱と界面粗さからの散乱が同時に観測され、両者を実験的に分離できない。もし、界面ラフネスが大きく、また、薄膜中の平均空孔サイズが比較的小さい場合、高角度側においては界面粗さの散乱が、空孔からの散乱と比較して相対的に弱くなることが予想される。そのような場合には、検出器ステージ角度２θの固定角度を高角度側に設定しロッキングスキャン測定を行っても、界面粗さの情報はあまり含まれず、十分な解析を行うことは不可能である。逆に、空孔の平均サイズが大きく界面ラフネスが小さい場合に、検出器ステージ角度２θの固定角度を低角度側に設定してロッキングスキャン測定を行っても、界面粗さの情報を取得することをほとんど不可能といえる。そのため、本ステップのロッキングスキャン測定においては、検出器ステージ角度２θの固定角度を複数選択することが好ましい。
【００８４】
以上のロッキングスキャン測定により、試料表面内方向の散乱ベクトルを変化させることによって、面内方向の電子密度揺らぎ相関距離およびその形状パラメータであるＨｕｒｓｔパラメータの情報を知ることが可能であることを示した。また、上記の通り、無作為に条件を設定してロッキングスキャン測定を実施した場合には、界面粗さの散乱と空孔散乱とを分離して測定することは不可能である。しかし、以下で示すとおり、膜試料の最表面における界面散乱（真空と最表面薄膜との界面における散乱）をほぼ独立に求めることが可能である。
【００８５】
界面散乱は、膜試料の界面のみで散乱が生じる。一方、空孔散乱は膜試料内部において発生するため、Ｘ線が膜試料に侵入する条件または薄膜中から空孔散乱が反射される条件を満足するときのみ観測される。このような条件は、Ｘ線の入射角度及び反射角度と薄膜の全反射臨界角度との大小により決定される。すなわち、Ｘ線入射角度および反射角度が全反射臨界角度より小さい条件のもとでは、最表層の界面粗さの散乱のみが検出される（図１１）。また、界面散乱は界面において発生する散乱であることから、界面における屈折波の強度が大きいときに、非常に強く観測される（図１２）。
【００８６】
Ｘ線入射角度またはＸ線反射角度が膜試料の全反射臨界角度を越えたときに、Ｘ線はＰｏｒｏｕｓ膜内部に侵入され、また、Ｐｏｒｏｕｓ膜から膜外に放出されるため、空孔散乱が観測される（図１３）。
【００８７】
以上に、この出願の発明である膜構造解析方法における基本的な態様の一例を示したが、以上で示した形態に限定されることはなく、その細部について、様々な形態をとりうることが考慮されるべきであることは言うまでもない。
【００８８】
この出願の発明は、以上の特徴を持つものであるが、以下に実施例を示し、さらに具体的に説明する。
【実施例１】
【００８９】
上記の通り、オフセットスキャン測定により得られる実測Ｘ線散乱曲線およびロッキングスキャン測定により得られる実測Ｘ線散乱曲線には、それぞれ空孔サイズ分布の情報と界面形状の情報が含まれている。多層構造の解析を行う場合、界面形状パラメータ数は膜数に比例して増加するため各パラメータを高い精度で算出することは困難であるが、以下に示すような手順により多層構造膜試料の測定データを扱うことにより、高い精度で構造解析を実現することができる。
１）  界面粗さのパラメータであるラフネスは、Ｘ線反射率解析の結果からある程度推測することができる。ただし、Ｘ線反射率測定で得られた結果を初期値として用いる場合は、Ｘ線反射率解析が十分行われている必要がある。
２）  Ｘ線入射角度またはＸ線反射角度が、最表面の全反射臨界角度より小さい条件では、最表面の界面粗さの情報のみを取り出すことが可能である。複数種類測定されるロッキングスキャン測定データにおいて、入射角度もしくは反射角度が全反射臨界角度より小さい角度領域のデータを用いて、フィッティングを行うことにより表面形状のパラメータをある程度推定することができる。
３）  多層膜中に数十ｎｍの膜が存在する場合、実測Ｘ線散乱曲線にそれらの層の上下界面で発生した散乱が干渉し、振動構造が観測される。これらの干渉は、式（IX）の第２項によって記述されるものである。式（IX）より、干渉の振幅の強さは、上下界面の散乱振幅の大きさが等しいときに最も強くなる。例えば、基板上に４０ｎｍのキャップ層を仮定し、基板で生じる界面散乱の散乱振幅をＡ、キャップ層表面の散乱振幅をＢとしたとき、式（IX）によれば、全散乱強度は次式
【００９０】
【数２１】
2005121613.tif 000023
【００９１】
により与えられる。ベースとなる｜Ａ｜2＋｜Ｂ｜2の項に対して、２（Ａ*Ｂ＋Ｂ*Ａ）で示される振動成分は、散乱振幅ＡとＢの大きさが等しいときに最も顕著に振動が観測される。
【００９２】
図１４〜１７は、キャップ層のラフネスを変化させた場合の干渉散乱の変化を示したものである。
【００９３】
図１４は、以下の表に示すように基板とキャップ層のラフネスを同じに設定した場合における干渉散乱の変化である。傾斜が等しいため、高角度側まで振動振幅がほぼ一定で観測されることがわかる。
【００９４】
【表１】
2005121613.tif 000024
【００９５】
図１５は、以下の表に示すように、図１４における条件から基板とキャップ層のラフネスを少し変化させた場合における干渉散乱の変化である。各界面の散乱強度がほぼ同じ角度（θ＝１．５ｄｅｇ近傍）のとき、最も干渉の振幅が大きく、１．５ｄｅｇからずれるにしたがって振幅が小さくなる。
【００９６】
【表２】
2005121613.tif 000025
【００９７】
図１６は、以下の表に示すように、図１５と同様に各界面の散乱強度が等しい角度で最も干渉の振動が強く観測されている。両界面のラフネスの値に大きな差が出てきたため、それ以外の角度で干渉の振幅の減衰が大きくなっている。両界面のラフネスにこの程度の差異が生じることにより、各界面におけるＸ線散乱曲線の傾斜に大きな差が現れていることがわかる。
【００９８】
【表３】
2005121613.tif 000026
【００９９】
図１７は、以下の表に示すように、図１６の設定をさらに顕著したものであり、各界面の散乱強度が等しい角度での干渉の振幅の大きさは、図１４〜１６と変わらないが、それ以外の角度領域において、各界面の散乱強度に大きな差があるため干渉の振幅は小さい。図１６と同様に、各界面でラフネスにこの程度の差異を設けることにより、傾斜の大きさに顕著な差が生じている。
【０１００】
【表４】
2005121613.tif 000027
【０１０１】
図１４〜１７中に示した基板において発生する散乱に対応する曲線（Ｓｕｂｓｔｒａｔｅ）は上式の｜Ａ｜2に相当し、また、図中のキャップ層において発生する散乱に対応する曲線（Ｃａｐ  Ｌａｙｅｒ）は上式の｜Ｂ｜2に相当している。散乱振幅であるＡおよびＢは複素数であるため絶対値の２乗で示しているが、各界面の散乱強度が等しい角度のとき、すなわち、２つの曲線が重なったときに、干渉の振幅が最も大きくなっていることがわかる。当然、両者のラフネスが異なれば、前述したように、各界面によって生じるオフセットスキャン測定により得られる実測Ｘ線散乱強度曲線の傾斜が変わるため、図１４〜１７で示されるように、全散乱強度には傾斜の異なる２本の直線がベースとなって測定される。それらの傾きからラフネスを評価することが可能であるが、干渉の振幅の大きさがどこで変化しているかを注目することによってもラフネスを評価することが可能である。
４）  各界面の散乱振幅の比を変えた場合にも、実測Ｘ線散乱曲線にそれらの層の上下界面で発生した散乱が干渉し、振動構造が観測される。図１８〜２１に、基板とキャップ層のラフネスとを等しい設定にして、両界面の面内相間距離を変化させた場合における干渉の強さを示す。図１４〜１７で示した場合と同様に、キャップ層の面内相間距離を変化させ場合、すなわち、各界面間の散乱振幅比を変化させた場合、両界面の散乱振幅の差が小さくなるにつれて干渉の振幅が大きくなることがわかる。ラフネスの値がある程度決まれば、干渉の振幅の強さから面内の相関距離を見積もることも可能である。
５）  多層膜中の任意の層における膜厚が大きい場合、その上下界面で生じる界面散乱の干渉はその振動周期が小さいため、測定が困難となる。だからと言って、入射側の光学系または反射側の光学系の分解能を向上させる操作を行っても、界面散乱自体の散乱確率が小さいため、十分な強度を得られないため、測定が困難である。そのような場合は、反射率解析からある程度正確なラフネスの初期値を見積もることができても、面内の相関距離とラフネスの形状パラメータ（Ｈｕｒｓｔパラメータ）の初期値を予想することが困難になる。この出願の発明である膜構造解析方法により、ある程度、他の層のパラメータが決まれば、厚膜層の界面のパラメータを最小二乗法により求めることもできるが、情報量が少なく複数の局所解に陥る場合がある。そのような場合には、すべてのパラメータを完全に精密化することが不可能である。ただし、膜の成膜条件や設計値等から、ある層の上の界面と下の界面のパラメータを同じであると仮定するなど、パラメータを束縛して解析することによりこのような問題を解消できる場合がある。
６）  前述してきたように、界面散乱と空孔散乱とでは多層膜中の空孔サイズ分布解析が異なる。また、界面散乱は、面内方向と表面法線方向とでは、異なる形状を取るものと考えられる。それに対し、空孔散乱は材料にもよるが、ほとんどの場合において膜中で異方性がないものと考えられる。逆にいえば、そのような材料を対象とした解析になっている。特に、表面界面散乱は、Ｘ線入射角度またはＸ線反射角度が薄膜最表面の全反射臨界角度近傍になるにつれて連続的に増加するのに対し、空孔散乱は入射Ｘ線または反射Ｘ線が全反射臨界角度より大きくならない限りは観測されない。また、空孔散乱においては、全反射臨界角度を越えたとたんに急激に散乱強度が増加する。また、空孔散乱のロッキングスキャンは平均空孔径が極端に大きくならない限り、界面散乱とは異なり、フラットなＸ線散乱曲線を与える。このように、空孔散乱と界面散乱との異なる特性を利用することにより、空孔サイズ分布解析は難しくない。
【０１０２】
以上のとおり、多層膜試料に対してもこの出願の発明である膜構造解析方法は適用可能である。以下に、以上で詳しく説明したこの出願の発明である膜構造解析方法を多層膜試料の構造解析に適用した場合の特徴をまとめる。
（Ａ）  表面散乱は、表面屈折波の強度依存性をもつ特徴を有しているため、比較的容易に表面構造を知ることができる。
（Ｂ）  式（IX）の第２項に示した界面散乱の干渉成分を利用することにより多層膜試料の各界面の構造を評価することができる。
（Ｃ）  空孔散乱は空孔層の全反射臨界角度より大きくなったときに急激に散乱が強くなるという特性を利用することにより、多層膜中の空孔サイズを解析できる。
【実施例２】
【０１０３】
次いで、基板上に、下層から順に、ポーラス膜、キャップ層、および、酸化膜層が形成された多層膜を試料として、ポーラス膜内の空孔サイズおよび多層膜の界面形状解析結果を示す。
【０１０４】
まず、Ｘ線反射率測定を行い、図２２に示すようなＸ線反射率曲線が得られた。このＸ線反射率曲線を用いてＸ線反射率解析をした結果、以下の通りの密度および膜厚が得られた。
【０１０５】
【表５】
2005121613.tif 000028
【０１０６】
このＸ線反射率解析で得られた密度と膜厚をもとに、薄膜中のＸ線固有状態を記述し、Ｘ線散漫散乱測定データを解析することにより、多層膜の各界面の形状およびポーラス層中の空孔サイズ解析を行う。ここで、ラフネスの初期値はＸ線反射率解析の解析値を用いた。
入射角度と出射角度を０．１°だけずらしたオフセットスキャン測定により得られた実測Ｘ線散乱曲線と2θを０．６、０．８、１．０、１．２、１．５、１．８、２．０に固定して測定した７種類のロッキングスキャン測定により得られた実測Ｘ線散乱曲線を用い、シミュレーションＸ線散乱曲線をフィッティングすることにより構造解析を行った。それぞれの実測Ｘ線散乱曲線に対するフィッティングの結果を図２８〜３０に示す。また、解析されたポーラス層における空孔サイズ分布について、図３１に示す。フィッティングにより得られた解析結果、次表のとおりである。
【０１０７】
【表６】
2005121613.tif 000029
【０１０８】
従来、オフセットスキャン測定により得られる単一の実測Ｘ線散乱曲線を用いて空孔サイズ分布解析等の膜構造解析を行っていたが、界面粗さ散乱が無視できないほど強く観測される場合は、界面散乱を考慮する必要がある。界面散乱はロッキングスキャン測定より得られる実測Ｘ線散乱曲線に特徴を持つことから、この出願の発明である膜構造解析方法においては、オフセットスキャン測定による実測Ｘ線散乱曲線と複数種類のロッキングスキャン測定による実測Ｘ線散乱曲線を用いることにより空孔サイズ分布解析と界面形状評価を同時に行った。この実施例においては、Ｘ線反射率により得られたキャップ層のラフネスが、Ｘ線反射率および散漫散乱パターンに与える影響が小さいため、キャップ層の界面形状の解析結果に多少誤差が含まれていることが予想される。しかし、空孔散乱に関しては、上記の通り、空孔散乱と界面散乱とが異なる実測Ｘ線散乱曲線の形状を与えるという特性を活かしているため、空孔サイズ分布の解析結果に大きな誤差はないものと考えられる。
【産業上の利用可能性】
【０１０９】
この出願の発明により、試料の界面の形状と空孔サイズの同時解析を実現し、また、多層からなる試料に対しても適用可能な膜構造解析方法とその装置が提供される。
【０１１０】
この出願の膜構造解析方法とその装置により、従来のＸ線反射率測定装置の構成を変更することなく、多層膜試料を対象とした試料の界面の形状と空孔サイズの同時解析が可能となることから、各種の機能性材料の開発における特性評価等に貢献するものと考えられ、その実用化が強く期待される。
【図面の簡単な説明】
【０１１１】
【図１】入射波と散乱波のＸ線固有状態について示した概要図である。
【図２】薄膜中における界面散乱による散乱過程について示した概要図である。
【図３】この出願の発明である膜構造解析装置の光学系の構成について示した概要図である。
【図４】オフセット角度を見積もるために行われるロッキングスキャン測定により得られたＸ線散乱曲線について示したグラフである。
【図５】オフセットスキャンにより得られる界面散乱による成分のシミュレーションＸ線散乱曲線について示したグラフである。
【図６】オフセットスキャンにより得られる空孔散乱による成分のシミュレーションＸ線散乱曲線について示したグラフである。
【図７】ロッキングスキャン測定により得られる実測Ｘ線散乱曲線について示したグラフである。
【図８】形状パラメータとしてＨｕｒｓｔパラメータを変化させた場合のロッキングスキャン測定により得られるシミュレーションＸ線散乱曲線について示したグラフである。
【図９】面内相関距離を変化させた場合のロッキングスキャン測定により得られるシミュレーションＸ線散乱曲線について示したグラフである。
【図１０】検出器ステージ角度２θの固定角度を様々に変化させたロッキングスキャン測定のシミュレーションＸ線散乱曲線について示したグラフである。
【図１１】ロッキングスキャン測定により得られる実測Ｘ線散乱曲線について示したグラフである。
【図１２】最表層の界面散乱によるＸ線散乱曲線について示したグラフである。
【図１３】多層膜モデルに対するロッキングスキャン測定により得られるシミュレーションＸ線散乱曲線について示したグラフである。
【図１４】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図１５】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図１６】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図１７】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図１８】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図１９】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図２０】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図２１】多層膜の界面散乱について算出したシミュレーションＸ線散乱曲線について示したグラフである。
【図２２】この出願の発明の実施例において行われたＸ線反射率解析結果について示したグラフである。
【図２３】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図２４】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図２５】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図２６】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図２７】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図２８】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図２９】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図３０】この出願の発明の実施例において得られたロッキングスキャン測定による実測Ｘ線散乱曲線に対するフィッティング結果について示したグラフである。
【図３１】この出願の発明の実施例において解析されたポーラス層における空孔サイズ分布について示したグラフである。
【符号の説明】
【０１１２】
３１  光源
３２  多層膜ミラー
３３  ナロースリット
３４  試料
３５  ナロースリット
３６  検出手段
【図１】
2005121613.tif 000030
【図２】
2005121613.tif 000031
【図３】
2005121613.tif 000032
【図４】
2005121613.tif 000033
【図５】
2005121613.tif 000034
【図６】
2005121613.tif 000035
【図７】
2005121613.tif 000036
【図８】
2005121613.tif 000037
【図９】
2005121613.tif 000038
【図１０】
2005121613.tif 000039
【図１１】
2005121613.tif 000040
【図１２】
2005121613.tif 000041
【図１３】
2005121613.tif 000042
【図１４】
2005121613.tif 000043
【図１５】
2005121613.tif 000044
【図１６】
2005121613.tif 000045
【図１７】
2005121613.tif 000046
【図１８】
2005121613.tif 000047
【図１９】
2005121613.tif 000048
【図２０】
2005121613.tif 000049
【図２１】
2005121613.tif 000050
【図２２】
2005121613.tif 000051
【図２３】
2005121613.tif 000052
【図２４】
2005121613.tif 000053
【図２５】
2005121613.tif 000054
【図２６】
2005121613.tif 000055
【図２７】
2005121613.tif 000056
【図２８】
2005121613.tif 000057
【図２９】
2005121613.tif 000058
【図３０】
2005121613.tif 000059
【図３１】
2005121613.tif 000060