Ｘ線回折測定解析方法及びプログラム
【課題】 既知情報及び実測データに基づく特定サテライトピークのフィッティングにより、十分な精度のプロファイル解析が行えるＸ線回折測定解析方法及びプログラムを実現し、また、一致状態を一般性を失うことなく大局的に見通しよく２次元表示する。
【解決手段】 実測プロファイルから周期Ｌと平均格子定数ａを求め、交互に積層された一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値を設定する設定段階と、Ｌ、ａ並びにＬ１及びａ１の初期値から、他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値を算出する変数算出段階と、Ｌ、ａ、Ｌ１及びａ１の初期値並びにＬ２、ａ２の計算値に基づいて解析プロファイルを得る解析プロファイル作成段階と、両プロファイルを比較して一致の度合いを調べる比較段階と、一致の度合いが一致基準を満たすか否かを判定し、満たすときはＬ１、ａ１、Ｌ２及びａ２の値を測定値とし、満たさないときは、Ｌ１及びａ１の初期値のいずれかを変更して再び変数算出段階に戻る判定段階とを含む。
【解決手段】 実測プロファイルから周期Ｌと平均格子定数ａを求め、交互に積層された一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値を設定する設定段階と、Ｌ、ａ並びにＬ１及びａ１の初期値から、他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値を算出する変数算出段階と、Ｌ、ａ、Ｌ１及びａ１の初期値並びにＬ２、ａ２の計算値に基づいて解析プロファイルを得る解析プロファイル作成段階と、両プロファイルを比較して一致の度合いを調べる比較段階と、一致の度合いが一致基準を満たすか否かを判定し、満たすときはＬ１、ａ１、Ｌ２及びａ２の値を測定値とし、満たさないときは、Ｌ１及びａ１の初期値のいずれかを変更して再び変数算出段階に戻る判定段階とを含む。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線回折測定解析方法及びプログラム、特に異種の薄膜を交互に積層した超格子構造の測定・解析に好適なＸ線回折測定解析方法及びプログラムに関する。
【０００２】
異種の薄膜を交互に積層した超格子構造を有する結晶構造は、通信分野で使用される高出力の半導体レーザの活性層を代表例として多くの用途に適用されている。
【背景技術】
【０００３】
近時、薄膜製造に対する要求精度が従前より更に厳しくなる中で、薄膜を測定・解析対象とする場合にＸ線回折を利用する高精度の測定解析法が多用されてきており、その中でも、異種の薄膜を交互に積層した超格子構造の測定・解析には、いわゆるロッキングカーブ測定・解析が利用されている。
【０００４】
ロッキングカーブ測定・解析を行なう従来のＸ線回折測定解析方法においては、基板に入射するＸ線の入射角を変化させていくと、その角度が基板結晶やその基板上に交互に積層された超格子構造の薄膜のブラッグ角にそれぞれ一致するところでＸ線が相互に強め合うように回折・出射することで、回折Ｘ線強度が回折角度の変化に応じてその基板及び超格子構造に由来する特定のパターンで変化する。したがって、回折角度を変えて結晶構造情報を含む回折Ｘ線強度変化のパターンを測定することで、そのプロファイル解析から超格子構造の異種薄膜の積層周期や膜厚や格子定数等の値を知ることができる。
【０００５】
前記測定に際しては、例えばＸ線源からのＸ線ビームを試料に照射し、回折Ｘ線の強度をＸ線カウンタでカウントすると、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされるＸ線の実測プロファイル中に、基板よりの回折Ｘ線に由来する基板ピークや、その基板上に積層された超格子構造に由来する複数のサテライトピークが現れる。
【０００６】
一方、Ｘ線回折現象による上記プロファイルの解析に際しては、基板上に超格子構造を持つ結晶構造体の既知の情報（基板の結晶構造情報や基板及び超格子構造の設計情報）に基づいて、Ｘ線のプロファイルを特徴付ける４つの変数をパラメータとして実測プロファイルに対応する解析プロファイルを理論計算で求め、前記４つの変数を変更しながらフィッティング処理を行なって、最終的に実測プロファイルに最も近い一致を見た時点の前記４つの変数を測定値としている。このように、Ｘ線回折で測定された実測プロファイルを解析することにより、結晶構造を数値的に把握することができる。
【０００７】
より具体的には、例えば、歪量子井戸半導体レーザの活性層を構成する量子井戸層の格子定数及び膜厚、並びにバリア層の格子定数及び膜厚の４つの値をパラメータとして、理論計算で解析プロファイルを求め、その解析プロファイルのパターンが実測プロファイルのパターンと、より良く一致するように理論計算のパラメータを変化させて、両プロファイルのフィッティングを行なう。そして、フィッティングにより最も良く一致する解析プロファイルの４つのパラメータ、すなわち量子井戸層の格子定数及び膜厚並びにバリア層の格子定数及び膜厚を測定値とする。なお、実測プロファイルと解析プロファイルを得る方法は、それぞれ既に確立されており学術界・産業界に広く普及している。
【０００８】
一方、上述のフィッティング処理は、対象によって異なるが、例えば最小二乗法近似によるフィッティングが良く知られている。この最小二乗法近似は、実測プロファイルと解析プロファイルの差の二乗誤差が最小になるようにフィッティングを行なうものであり、プロファイルの全角度範囲にわたって両プロファイルの差の二乗値を均等に加算し平方根をとりサンプル数で除算するという単純なアルゴリズムを採用でき、計算機を用いて容易に自動化することができる。また、パラメータの初期値を適切に設定すれば、コンピュータプログラムによって自動的にもっとも良く一致する点を見つけ出すことができるので、比較的広いパラメータ範囲を探索するような場合に適している。
【０００９】
しかし、実測プロファイルにはノイズ等が多く含まれており、パターン認識処理も容易ではないので、ソフトウェアにより最も一致した点を見つけ出すために、例えば、実測プロファイルデータに対してスムージング（平滑化）処理を施してノイズを含む実測プロファイルから超格子構造の特徴を示す主要なピーク（例えばサテライトピーク）についてのフィッティングを実行し、次いで、このスムージングを徐々に減らしていきながら、プロファイルの細かい特徴部分についてのフィッティングを行なうという手法が採られている（例えば、特許文献１、非特許文献１参照）。
【特許文献１】特表２００２−５３２７０３号公報
【非特許文献１】ＰＡＮａｌｙｔｉｃａｌ社の「Ｐｒｅｓｓ Ｒｅｌｅａｓｅｓ」、「New X'Pert Epitaxy software from Philips Analytical makes rocking-curve analysis easier than ever」［２００４年１２月１０日検索］、インターネットＵＲＬ＜http://www.panalytical.com/index.cfm?pid=21&itemid=153&contentitemid=12＞。
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、上記従来のＸ線回折測定解析方法及びプログラムにあっては、最小二乗法近似を用いているため、スムージングによる主要ピークについての第１段階のフィッティングとスムージングを減らしながらの第２段階のフィッティングとを行ったとしても、小さなピークやノイズまでフィッティングの対象としてしまうので、条件次第で本来１つしかない一致点が複数得られてしまうという結果を招きかねないという問題があった。
【００１１】
一方、基板の結晶に関する格子定数、並びに超格子構造の各種薄膜の配置、膜厚、格子定数等の設計値は既知であり、更に、実測プロファイルには、基板ピークと共に、超格子構造の平均格子定数に対応する０次のサテライトピークを中心に超格子構造の周期（Ｌ）、すなわち、交互に積層される薄膜の膜厚の和に相当する回折角度間隔を隔てて複数の他次数のサテライトピークが複数現れる。すなわち、実測プロファイルの形状的特徴とパラメータの間の特定の関係が明らかになっている。
【００１２】
したがって、このように実測プロファイルの形状的特徴と前記４つの変数（パラメータ）の明らかな関係を利用し、４つの変数を２つの独立変数に集約すれば、重要情報を有するプロファイル形状部分（特定のサテライトピーク）のみに着目してフィッティングを行なうことができる。このようにすれば、小ピーク、ノイズ等の微視的プロファイル部分の悪影響を抑えた簡便且つ確実な解析が可能になる。
【００１３】
そこで、本発明は、サテライトピーク間にある小さなピークやノイズはフィッティングに用いずに超格子構造を反映する主要サテライトピークについてフィッティングを行ない、既知の情報及び実測データに基づいて選択した特定のサテライトピークのフィッティングによって、必要にして十分な精度の解析が行なえるＸ線回折測定解析方法及びプログラムを実現することを目的とするものであり、併せて、４パラメータの一致状態を一般性を失うことなく大局的に見通しよく理解できるような２次元表示を実現することを目的とする。
【課題を解決するための手段】
【００１４】
本願の発明者は、解析プロファイルの４つの変数（パラメータ）を変化させて解析プロファイルを実測プロファイルに一致させることに代えて、まず、超格子物理の自然法則から、これら変数の従属関係を明確にすることで、４つの変数を２つの独立変数に集約できるとの知見を得、簡便、正確で局所最適値に陥る危険性のない優れた測定解析方法を創作したものである。また、これにより、２変数による鳥瞰図的な４変数の全体像を大局的に見通しよく表示することも可能となる。
【００１５】
すなわち、本発明は、上記目的達成のため、（１）基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、前記４つの変数を変化させながら比較し、最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析方法において、前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求め、且つ前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求めるとともに、前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１）の初期値を設定する設定段階と、前記超格子構造の周期（Ｌ）と前記超格子構造の平均格子定数（ａ）と前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を算出する変数算出段階と、前記超格子構造の周期（Ｌ）、前記超格子構造の平均格子定数（ａ）、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値、並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値に基づいて、前記理論計算により前記解析プロファイルを得る解析プロファイル作成段階と、前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較段階と、前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して再び変数算出段階に戻る判定段階と、を含むことを特徴とするものである。
【００１６】
本発明の方法では、まず、実測プロファイルから超格子構造の周期（Ｌ）と超格子構造の平均格子定数（ａ）が求められ、これらの値と既知の情報とに基づいて、第１又は第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値が４つの変数のうち２つの変数の初期値として設定され（設定段階）、これらの値を基に残りの２つの変数が算出されると（変数算出段階）、それまでに得られた情報に基づいて、理論計算により解析プロファイルが作成される（解析プロファイル作成段階）。そして、その解析プロファイル及び実測プロファイルが比較され（比較段階）、両プロファイルの一致の度合いが一致基準を満たすときには、そのときの４つの変数が測定値とされ、一致の度合いが前記基準を満たさないときには、２つの変数である設定値のいずれかを変更して再び変数算出段階に戻る（判定段階）。したがって、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、適切な２つの初期値の設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【００１７】
この発明の方法においては、（２）前記設定段階において、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイルから周期性のある孤立したピークを抽出して前記主要サテライトピークを特定するのが好ましい。これにより、小ピーク、ノイズ等の微視的プロファイル部分の悪影響をより有効に回避でき、しかも、処理を簡単化できることになる。
【００１８】
この発明の方法においては、さらに、（３）前記設定段階において、前記実測プロファイルから求めた前記超格子構造の周期（Ｌ）及び平均格子定数（ａ）と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値を設定するのがよい。実測プロファイルや設計情報等から得られる情報を基に、より得るべき測定値に近い初期値が設定できるからである。
【００１９】
また、この発明の方法においては、（４）前記比較段階において、前記複数のサテライトピークのうち前記超格子構造の平均格子定数（ａ）に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせ、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べるのが好ましい。これにより、４つの変数の変化に対する解析プロファイルの変化の特性を考慮して、一致度の比較処理を簡単化できることになる。
【００２０】
また、この発明の方法においては、（５）前記判定段階において、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのがよい。超格子構造の周期（Ｌ）は実測プロファイルから正確に得られるからであり、この周期（Ｌ）を固定することで、一致点への収束を早めることができる。
【００２１】
この場合、さらに、（６）前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのに先立って、前記一方の薄膜の格子定数（ａ１）を微調整するのが望ましい。膜厚（Ｌ１）と格子定数（ａ１）の変化における相関性を考慮して、一致点への収束を早めることができるからである。
【００２２】
また、この発明の方法においては、（７）前記判定段階で、前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶しておき、該記憶情報に基づいて前記解析プロファイル及び前記実測プロファイルの一致の度合いを示す分布図を作成するとともに、該分布図上に、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とを視覚的に区別される表示態様で表示し、前記分布図上の点が、該分布図上の前記一致領域及び環状分布領域のうちどの領域に属するかで前記一致の度合いを視覚的に識別可能に表示するのが好ましい。これにより、解析プロファイル及び実測プロファイルの一致の度合いがどの程度かを容易に見て取れることになり、ユーザーの迅速・的確な判断が可能になる。
【００２３】
本発明は、また、（８）基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、前記４つの変数を変化させながら比較し、最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、コンピュータにより、前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求める周期算出手段（１１）と、前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求める平均格子定数算出手段と、前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１）の初期値を設定する初期値設定手段と、前記超格子構造の周期（Ｌ）と前記超格子構造の平均格子定数（ａ）と前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値をそれぞれ算出する変数算出手段と、前記超格子構造の周期（Ｌ）、前記超格子構造の平均格子定数（ａ）、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値に基づいて、前記理論計算により前記解析プロファイルを作成する解析プロファイル作成手段と、前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較手段と、前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値、並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段と、の機能を実現させることを特徴とするものである。
【００２４】
このプログラムにより、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、適切な２つの初期値の設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【００２５】
本発明のプログラムにおいては、（９）前記初期値設定手段が、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイルから前記超格子構造の周期（Ｌ）に対応する周期性を有する孤立した複数のピークを抽出して前記主要サテライトピークを特定するのがよい。このようにすると、小ピーク、ノイズ等の微視的プロファイル部分の悪影響をより有効に回避でき、しかも、処理を簡単化できることになる。
【００２６】
また、この発明のプログラムにおいては、（１０）前記初期値設定手段が、前記実測プロファイルから求めた前記超格子構造の周期（Ｌ）及び平均格子定数（ａ）と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値を設定するのがよい。この構成により、実測プロファイルや設計情報等から得られる情報を基に、より得るべき測定値に近い初期値が設定可能となる。
【００２７】
また、この発明のプログラムにおいては、（１１）前記比較手段が、前記複数のサテライトピークのうち前記超格子構造の平均格子定数（ａ）に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせるシミュレーションを実行し、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べるのが好ましい。これにより、４つの変数の変化に対する解析プロファイルの変化の特性を考慮して、一致度の比較処理を簡単化できることになる。
【００２８】
この発明のプログラムにおいては、さらに、（１２）前記判定手段が前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのがよい。これにより、超格子構造の周期（Ｌ）は実測プロファイルから正確に得られるからであり、この周期（Ｌ）を固定することで、一致点への収束を早めることが可能となる。
【００２９】
この場合、さらに、（１３）前記判定手段が前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのに先立って、前記一方の薄膜の格子定数（ａ１）を微調整するのが望ましい。これにより、膜厚（Ｌ１）と格子定数（ａ１）の変化における相関性に基づき、一致点への収束を早めることが可能となる。
【００３０】
また、この発明のプログラムにおいては、（１４）前記コンピュータに、前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶する記憶手段と、該記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを示す分布図の画像情報を生成する分布図作成手段と、の機能を実現させるのがよい。このプログラムにより、複数の測定値の分布を、解析プロファイル及び実測プロファイルの一致の度合いの分布として視覚的に容易に把握可能となる。
【００３１】
この場合、さらに、（１５）前記分布図作成手段が、前記記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とを視覚的に区別される表示態様で表示するのが好ましい。これにより、解析プロファイル及び実測プロファイルの一致の度合いがどの程度かを容易に見て取れることになり、ユーザーの迅速・的確な大局的判断が可能になる。
【００３２】
本発明のプログラムは、あるいは、（１６）基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を異なる角度で入射したとき測定情報に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の実測プロファイルと、４つの変数として前記第１組成薄膜の格子定数（ａ１）及び膜厚（Ｌ１）並びに前記第２組成薄膜の格子定数（ａ２）及び膜厚（Ｌ２）を用いた理論計算に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の解析プロファイルとを、前記解析プロファイルの４つの変数を変化させながら比較し、最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、ｎを正の整数とするとき、ｎ個の前記角度の値（ｘ１，ｘ２，・・・ｘｎ）及び該角度に対応するｎ個の回折Ｘ線強度の測定値（ｇ１，ｇ２，・・・ｇｎ）を順次メモリに記憶させる実測データ記憶手段（４１）と、前記メモリに記憶されたデータに基づいて前記角度に応じて変化する回折Ｘ線強度の実測プロファイルを表示出力可能な所定のプロファイル表示形式で作成する実測プロファイル作成手段（４２）と、前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求める周期算出手段（４７）と、前記実測プロファイルの中に超格子構造の平均組成に対応して現れるサテライトピークと前記実測プロファイルの中に前記基板に由来して現れる基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求める平均格子定数算出手段（４７）と、前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１又はＬ２）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１又はａ２）の初期値を、前記４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）のうち２つの変数の初期値として設定する初期値設定手段（４４）と、前記超格子構造の周期（Ｌ）と、前記超格子構造の平均格子定数（ａ）と、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値とから、前記４つの変数（ａ１、Ｌ１、ａ２、Ｌ２）のうち残りの２つの変数である前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値をそれぞれ算出する変数算出手段（４８）と、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値と、前記変数算出手段（４８）で算出された前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値とを、前記４つの変数として、前記複数のサテライトピークのうち次数の異なる特定のサテライトピークに対応する角度範囲を特定する前記ｎより少ない正の整数であるｍ個の角度値（ｘｉ，ｘｊ，・・・ｘi+m-1；但し、ｉ、ｊはそれぞれｍより小さい正の整数）について、回折Ｘ線強度の値（ｆｉ，ｆｊ，・・・ｆi+m-1）を前記理論計算により算出して前記解析プロファイルを作成する解析プロファイル作成手段（４９）と、前記解析プロファイル及び前記実測プロファイルを前記複数のサテライトピークのうち前記次数の異なる特定のサテライトピークにつき比較して、両プロファイルの一致の度合いを調べる比較手段（５１）と、前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記解析プロファイルの４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）を前記一方の薄膜の膜厚及び格子定数並びに前記他方の薄膜の膜厚及び格子定数の測定値として出力し、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段（５２）と、の機能をコンピュータを使用して実現させることを特徴とするものである。
【００３３】
このプログラムにより、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、適切な２つの初期値の設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【発明の効果】
【００３４】
本発明によれば、サテライトピーク間にある小さなピークやノイズはフィッティングに用いることなく、既知の情報及び実測データに基づいて選択した特定のサテライトピークについて、適切な２つの変数の初期値設定とこれら初期値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行するので、超格子構造を反映する主要サテライトピークについて迅速で且つ有効なフィッティングを行ない、必要にして十分な精度の解析が行なえるＸ線回折測定解析方法及びプログラムを実現することができる。また、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、一致状態を一般性を失うことなく大局的に見通しよく理解できるような２次元表示を実現することができる。
【００３５】
すなわち、本発明は、解析プロファイルの４つの変数の従属関係を明確にすることで、４つの変数を２つの独立変数に集約させ、簡便、正確で局所最適値に陥る危険性のない優れた測定解析方法を創作したものであり、これにより、２変数による鳥瞰図的な４変数の全体像を大局的に見通しよく表示することも可能としたものである。
【発明を実施するための最良の形態】
【００３６】
以下、本発明の好ましい実施の形態について、図面を参照して説明する。
【００３７】
図１乃至図９は本発明の１つの実施の形態に係るＸ線回折測定解析方法を実施するための測定解析装置を示す図であり、この装置は本発明のＸ線回折測定解析プログラムに従って後述する処理を実行する。
【００３８】
まず、その構成について説明する。
【００３９】
図１及び図２に模式的に示すように、この測定装置では、第１回転台１に支持されたサンプル基板１０に対し、Ｘ線源２０からのＸ線ビームが照射されるとともに、サンプル基板１０で回折したＸ線を第２回転台２に支持されたＸ線カウンタ３０で検出し、そのＸ線強度を測定する光学系が構成されている。第１回転台１は、サンプル基板１０へのＸ線ビームの入射角度が変化するよう角速度ωで回転（往復回動）可能であり、第２回転台２は第１回転台１と同一中心軸回りに第１回転台１とは異なる角速度で回転可能である。これら回転台１、２は、それぞれに対応する公知の駆動手段３、４によって回転駆動されるとともに、回転中のサンプル基板１０の基板面と入射Ｘ線ビームのなす角度θ（視射角）に対して、Ｘ線ビームの照射軸線（図１に一点鎖線で示す）に対しＸ線カウンタ３０のＸ線受光軸線のなす角度が常に角度２θに保たれるように回転制御されるようになっている。すなわち、サンプル基板１０に対するＸ線ビームの入射角変化に対して常に回折出射角方向にＸ線カウンタ３０が配置された状態を保つようになっている。
【００４０】
サンプル基板１０は、基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造体で、例えばインジウム燐（以下、ＩｎＰという）基板上に異種のインジウムガリウム砒素燐（ＩｎＧａＡｓＰ）の薄膜を積層したもので、半導体レーザの活性層となるバリア層と量子井戸層の薄膜（例えば６０オングストローム程度）を積層した多重量子井戸（ＭＱＷ）構造を有している。
【００４１】
Ｘ線源２０は、例えばＸ線管２１、湾曲ミラー２２及びモノクロメータ２３によって構成されており、Ｘ線管２１で発生したＸ線束を湾曲ミラー２２によってサンプル基板１０の測定面上に収束させるとともに、モノクロメータ２３によって特定波長域に単色化するようになっている。また、Ｘ線カウンタ３０は、例えばシンチレーションカウンタで構成されており、回折Ｘ線の強度に応じたカウント値を出力するようになっている。
【００４２】
図２に示すように、Ｘ線源２０のＸ線管２１の駆動回路５並びに第１回転台１及び第２回転台２の駆動手段３、４はそれぞれ測定解析制御装置４０に接続されており、この測定解析制御装置４０によって上述のような測定系の動作制御がなされる。また、Ｘ線カウンタ３０の信号処理回路３１からのカウント情報並びに第１回転台１及び第２回転台２の回転位置検出器６、７の検出情報がそれぞれ測定解析制御装置４０に取り込まれるようになっており、測定解析制御装置４０はこれら入力情報及び設定・入力部４０６からの設定入力又は操作入力等に基づいて、所定のＸ線回折測定解析プログラムを実行し、公知のディスプレイデバイスからなる表示部４０７に横軸を回折角度（秒単位）とし縦軸を回折Ｘ線強度とする回折Ｘ線図形を表示出力したり、その表示画像情報を外部に出力したりするようになっている。
【００４３】
測定解析制御装置４０は、詳細のハードウェア構成を図示しないが、概略、例えばＣＰＵ（Central Processing Unit）４０１、ＲＯＭ（Read Only Memory）４０２、ＲＡＭ（Random Access Memory）４０３及び入出力インターフェース回路４０４に加え、ディスクドライブ等の補助記憶装置４０５を内蔵して構成され得るもので、前記設定・入力部４０６及び表示部４０７が付設されている。また、この測定解析制御装置４０は、ＣＰＵ４０１がＲＯＭ４０２及び補助記憶装置４０５に記憶された制御プログラムに従って、ＲＡＭ４０３との間でデータを授受しながら、Ｘ線カウンタ３０の信号処理回路３１からのカウント情報、回転位置検出器６、７の検出情報及び設定・入力部４０６からの入力情報等に基づいて、回折Ｘ線強度や回折角度を把握する演算処理を実行するとともに、図３に機能ブロックで示す各部の機能を発揮するようになっている。
【００４４】
図３に示すように、測定解析制御装置４０においては、回折Ｘ線強度とそのときの回折角度の情報が順次実測データ記憶部４１に記憶される。すなわち、上述のようにサンプル基板１０に対するＸ線ビームの入射角度を変化させながら、サンプル基板１０の結晶構造にＸ線を入射し回折されたＸ線をＸ線カウンタ３０で検出して得られるＸ線回折測定情報が、実測データ記憶部４１に順次記憶される。このとき、具体的には、例えば、ｎ（正の整数）個の角度の値ｘ１，ｘ２，・・・ｘｎ及びその角度に対応するｎ個の回折Ｘ線強度の測定値ｇ１，ｇ２，・・・ｇｎを順次メモリに記憶させることになる。
【００４５】
実測プロファイル作成部４２は、この実測データ記憶部４１の記憶データに基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる実測プロファイル、すなわち実測データから得られる回折Ｘ線図を公知の図形化処理によって作成する。この実測プロファイルは、Ｘ線回折測定結果として、横軸に角度、縦軸にＸ線強度をとったグラフで表示され、例えば図４に示すように、その横軸の回折角度が、サンプル基板１０のうち基板層の結晶の格子定数に対応した最大の鋭いピーク（基板ピーク）Ｐ10を中心に相対角度（秒）で示され、縦軸の回折Ｘ線強度が毎秒当たりのカウント数（counts／sec）で示されている。
【００４６】
実際の実測プロファイルは、図４に例示するように、サンプル基板１０の基板層よりの回折Ｘ線に由来する基板ピークＰ10の他に、そこに積層された超格子構造に由来する複数のサテライトピークＰ20（０次サテライトピーク）、Ｐ21（−１次サテライトピーク）、Ｐ22（＋１次サテライトピーク）、Ｐ23（＋２次サテライトピーク）等が超格子における異種薄膜の積層周期Ｌに対応する角度間隔で現れる。なお、図４中の他のピークは超格子構造に含まれる他の薄膜、例えば半導体レーザの光ガイド層を構成する複数層のインジウムガリウム砒素燐（ＩｎＧａＡｓＰ）の薄膜層（例えば１００オングストローム程度）からの回折Ｘ線等に基づくものである。
【００４７】
サンプル基板１０となる超格子構造の設計情報（交互に積層される異種薄膜の積層周期、各層の薄膜の膜厚及び格子情報、配置等）は、設計情報記憶部４３に予め設定記憶されている。
【００４８】
この設計情報記憶部４３の記憶情報は、実測プロファイル作成部からの実測プロファイル情報と共に特定ピーク選択部４５に取り込まれ、この特定ピーク選択部４５において、積層された超格子構造に関する設計値情報に基づき、実測プロファイル中で前記超格子構造の周期Ｌに対応する周期性を有する孤立した複数のピークの発生領域が特定される。
【００４９】
ここで特定のピークの発生領域を特定する条件は、ピーク特定条件記憶部４６に記憶されている。この条件は、例えば超格子を構成する異種の薄膜の格子定数とサンプル基板１０の基板層の結晶の格子定数との大小関係、並びに、それに応じて異なるサテライトピークの現れ方の特性に関する情報に基づいて、サテライトピークのうち孤立したサテライトピーク（例えば図４に示すサテライトピークＰ20、Ｐ21のように他の大きなピークとの干渉が少なく独立してピーク波形を形成しているもの）が発生する蓋然性が高い領域を回折角度範囲で特定する条件である。例えば、第１組成薄膜の格子定数が基板層（ＩｎＰ）の格子定数より大きく、第２組成薄膜の格子定数が基板層の格子定数と同等の設計仕様の場合には、通常、孤立したサテライトピークが基板ピークの左側（マイナス側角度域で超格子の積層周期に応じた所定角度範囲の特定領域内）に複数現れ、あるいは、第１組成薄膜の格子定数が基板層（ＩｎＰ）の格子定数と同等で、第２組成薄膜の格子定数が基板層の格子定数より小さい設計仕様の場合には、通常、孤立したサテライトピークが基板ピークの右側（プラス側角度域で超格子の積層周期に応じた所定角度範囲の特定領域内）に複数現れ、あるいはまた、第１組成薄膜の格子定数が基板層（ＩｎＰ）の格子定数より大きく、第２組成薄膜の格子定数が基板層の格子定数より小さい設計仕様の場合には、通常、孤立したサテライトピークが基板ピークの両側にそれぞれ複数現れる。そして、これらの孤立したサテライトピークがフィッティングに有効なサテライトピークとなる一般的傾向がある。このような条件に基づいて、フィッティングに有効な特定サテライトピークの発生領域を判定することができる。
【００５０】
また、特定ピーク選択部４５で孤立したサテライトピークを抽出する処理は、例えば特定のサテライトピークの発生領域内でのサテライトピークの間隔より十分に狭い所定角度範囲の単位で平均の回折Ｘ線強度を求めて各次数のサテライトピークを抽出することができる。この際、サテライトピークの間隔は設計情報に基づいて概ね予測でき、各次数のサテライトピークを特定（角度範囲で特定）するためのサンプリングは比較的少数で足りる。なお、０次や−１次、＋１次のサテライトピークが基板ピークや他の層に起因するピークに埋もれている場合は、他の孤立したサテライトピークを使用する。
【００５１】
特定ピーク選択部４５でこのようにピーク特定条件記憶部４６の記憶情報に基づいて特定された主要サテライトピークに関するプロファイル情報は、周期・平均格子定数算出部４７に送られ、ここで、第１組成薄膜及び第２組成薄膜の積層周期、すなわち多重量子井戸構造の井戸層及びバリア層からなる超格子の積層周期Ｌと多重積層された超格子構造の平均格子定数ａが算出される。
【００５２】
周期・平均格子定数算出部４７においては、第１組成薄膜及び第２組成薄膜の積層周期Ｌが、実測プロファイル中の主要なサテライトピークの角度間隔から、具体的には、例えば特定の次数の異なるサテライトピークＰ20、Ｐ21が現れる回折角度θｉ、θｊと、Ｘ線の波長とから、次式（１）により直接的に精度良く算出できる。ただし、ｉ、ｊはそれぞれサテライトピークの次数である。
【００５３】
Ｌ＝（ｉ＋ｊ）λ／（２（ｓｉｎθｉ−ｓｉｎθｊ）） ・・・（１）
ここで、第１組成薄膜及び第２組成薄膜の積層周期Ｌに対して、第１組成薄膜（たとえば井戸層）の膜厚をＬ１とし、第２組成薄膜（例えばバリア層）の膜厚をＬ２とすると、これらは次式（２）で表すことができる。
【００５４】
Ｌ＝Ｌ１＋Ｌ２ ・・・（２）
したがって、２つのパラメータＬ１、Ｌ２には従属関係があり、パラメータとしてのＬ２を定めると、次式（３）のようにＬ１が定まり、逆に、Ｌ１を定めると自動的にＬ２が定まる。
【００５５】
Ｌ１＝Ｌ−Ｌ２ ・・・（３）
周期・平均格子定数算出部４７においては、また、実測プロファイル中の複数のサテライトピークのうち、超格子構造の平均組成である平均格子定数ａに対応する０次サテライトピークＰ20が現れる回折角度と、サンプル基板のＩｎＰ基板層の結晶による基板ピークP10が現れる回折角度との差から、平均格子定数ａが精度良く直接的に算出できる。この平均格子定数ａは、また、次式（４）で表すことができる。
【００５６】
ａ＝（ａ１・Ｌ１＋ａ２・Ｌ２）／（Ｌ１＋Ｌ２） ・・・（４）
この（４）式中では、量子井戸層とバリア層の格子定数をそれぞれａ１、ａ２で示しており、パラメータとしてａ１を選択すると、次のように表せる。
【００５７】
ａ１＝（ａ・Ｌ−ａ２・Ｌ２）／（Ｌ−Ｌ２） ・・・（５）
このように、周期・平均格子定数算出部４７は、複数の主要なサテライトピークの間隔から超格子構造の周期Ｌを求める周期算出手段の機能と、実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークＰ20とＩｎＰ基板に由来して実測プロファイル中に現れる基板ピークＰ10との間隔から超格子構造の平均格子定数ａを求める平均格子定数算出手段の機能と、を併有している。
【００５８】
ところで、上記（４）式の右辺の値が全て既知であるとき、すなわち、格子定数ａ２が定められると、従属関係により、一義的に他の格子定数ａ１が定まる。勿論、逆も真なりで、格子定数ａ１を定めると、自動的に格子定数ａ２が定まる。すなわち、従属関係によって、４つのパラメータが２つのパラメータに集約され得るということができる。
【００５９】
そこで、初期値設定／変更部４４（初期値設定手段）では、４つのパラメータ（変数）Ｌ１、Ｌ２、ａ１、ａ２のうち、２つのパラメータ、すなわち、第１組成薄膜又は第２組成薄膜のうち一方の薄膜の膜厚Ｌ１の初期値と、第１組成薄膜又は第２組成薄膜のうち一方の薄膜の格子定数ａ１の初期値とを、それぞれ後述するように設定するとともに、その設定値を変更するようになっている。なお、これら膜厚Ｌ１の初期値と格子定数ａ１の初期値がそれぞれ設定されると、上述の従属関係から、第１組成薄膜又は第２組成薄膜のうち他方の薄膜の膜厚Ｌ２の初期と、第１組成薄膜又は第２組成薄膜のうち他方の薄膜の格子定数ａ２は一義的に定まる。
【００６０】
また、初期値設定／変更部４４は、実測プロファイルから求めた超格子構造の周期Ｌ及び平均格子定数ａの値と、これらの設計値並びに前記一方及び他方の薄膜の膜厚Ｌ１、Ｌ２及び格子定数ａ１、ａ２の設計値といった設計値情報とに基づいて、周期Ｌ及び平均格子定数ａの設計値からのずれをも考慮して、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値を例えば設計値と同一の相互比率を保つようにそれぞれ設定する。さらに、初期値設定／変更部４４は、後述する条件判定部５２（判定手段）からの更新指令を受けたとき、設定した初期値を所定の微調整量を単位として順次変更（増加又は減少）させる機能を有している。
【００６１】
変数算出部４８（変数算出手段）は、初期値設定／変更部４４で設定された前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の各初期値又は変更された膜厚Ｌ１及び格子定数ａ１の各値と、周期・平均格子定数算出部４７で算出された超格子構造の周期Ｌ及び超格子構造の平均格子定数ａの値とに基づいて、上述の（３）、（５）式に示すような条件で、超格子を構成する他方の薄膜の膜厚Ｌ２及び格子定数ａ２をそれぞれ算出する。
【００６２】
そして、上述のようにして得られた４つのパラメータ、すなわち、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値、並びに、前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値は、それぞれ解析プロファイル作成部４９に取り込まれ、この解析プロファイル作成部４９で、４つパラメータＬ１、Ｌ２、ａ１、ａ２に基づいて、これらの値を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルが作成される。そして、このとき、プロファイル作成領域が特定領域に特定されていることから、４つのパラメータＬ１、Ｌ２、ａ１、ａ２（変数）に基づき、複数の主要なサテライトピークのうち次数の異なる特定のサテライトピークに対応する角度範囲を特定する前記ｎより少ない正の整数であるｍ個の角度値ｘｉ，ｘｊ，・・・ｘi+m-1（但し、ｉ、ｊはそれぞれｍより小さい正の整数）について、回折Ｘ線強度の値ｆｉ，ｆｊ，・・・ｆi+m-1を前記理論計算により算出して解析プロファイルを作成することになる。
【００６３】
比較部５１では、解析プロファイル及び実測プロファイルのデータを実測プロファイル作成部４２及び解析プロファイル作成部４９からそれぞれ取り込み、両プロファイルを比較して、両プロファイルの一致の度合いを調べる。ここでは、次数の異なる特定のサテライトピークにつき比較して、両プロファイルの一致の度合いを調べるが、この比較は、例えば比較対象のサテライトピーク領域について平均二乗誤差Ｅを求めることで行われる。この平均二乗誤差Ｅは。比較対象として選択された特定のサテライトピークの個々のピーク形成角度範囲において、各サンプリング点の実測プロファイルの値、例えばｆｉ（ｉ＝１、２、・・・ｎ）と解析プロファイルの値ｇｉの差（ｆｉ−ｇｉ）の二乗値を均等に加算してその平方根をとり、サンプル数ｎで除算したものである。すなわち、Ｅ＝｛（ｆ１−ｇ１）２＋（ｆ２−ｇ２）２＋・・・＋（ｆｎ−ｇｎ）２｝１／２／ｎを比較結果とする。
【００６４】
この比較部５１での比較結果は、条件判定部５２（判定手段）に取り込まれ、この条件判定部５２において、各比較対象のサテライトピークについての平均二乗誤差Ｅが予め設定された許容値以下に達しているか否かが判定される。すなわち、前記一致の度合いが予め定められた一致基準（許容値以下）を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値、並びに前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値を、それぞれの膜厚測定値及び格子定数の測定値とする。一方、両プロファイルの一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれか１つの設定値を所定の微調整量単位で初期値設定／変更部４４に変更させて、変数算出部４８に再度変数の算出を実行させる。また、その結果が再度一致基準を満たさないときは、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうち残りの１つの設定値を所定の微調整量単位で初期値設定／変更部４４に変更させて、変数算出部４８に再度変数の算出を実行させる。それ以降に、一致基準を満たさないときは、まず、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれか１つの設定値を所定の微調整量単位で前回の設定値から変更させ、次いで、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうち残りの１つの前回の設定値を所定の微調整量単位で再度変更させて、変数の再算出を実行させることになる。また、このような変数の変更による解析プロファイルの再作成（シミュレーション）を所定回数実行してもなお一致基準を満たさない場合は、当該サンプル基板１０を不良と判定する。
【００６５】
また、本実施形態においては、条件判定部５２は、比較部５１と協働して、複数のサテライトピークのうち超格子構造の平均格子定数ａに対応する０次のサテライトピーク（ピークの点のみ）について、解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ実測プロファイルに合わせるシミュレーションを最初に実行するようになっており、これらによって、比較手段のシミュレーション処理部が構成されている。
【００６６】
そのため、条件判定部５２が前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して変数算出部４８に再度算出を実行させるとき、初期値設定／変更部４４及び変数算出部４８は、超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整し、更に、このように超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整するのに先立って、前記一方の薄膜の格子定数ａ１を微調整するようになっている。
【００６７】
また、比較部５１で比較対象として選択されるサテライトピークは、初回は０次のサテライトピークとこれとは次数の異なるサテライトピーク（例えばマイナス側の一次サテライトピーク）であり、０次のサテライトピークについての平均二乗誤差Ｅが所定値以下となるように設定値が変更された後は、次数の異なるサテライトピークのみが比較対象であってもよい。すなわち、比較部５１は、０次のサテライトピークについて両プロファイルを一致させるシミュレーション処理の後、主として他の次数のサテライトピークについて、前記一致の度合いを調べる誤差解析部となっている。
【００６８】
このように、測定解析制御装置４０は、補助記憶装置４０５に格納されたＸ線回折測定解析プログラムに従って、解析プロファイルを特定する４つのパラメータを２つに集約させつつ変化させながら実測プロファイルと解析プロファイルを比較し、両プロファイルが最も良く一致した状態における解析プロファイルの４つの変数を、第１、第２組成薄膜の格子定数ａ１、ａ２及び膜厚Ｌ１、Ｌ２のそれぞれの測定値として出力する。
【００６９】
また、測定解析制御装置４０は、複数のサンプル基板１０について上述のような測定解析が繰り返される場合に、条件判定部５２により前記一致の度合いが予め定められた一致基準を満たすと判定されたとき、あるいは不良判定がなされたときにも、そのときの前記一方又は他方の薄膜の膜厚Ｌ１又はＬ２と前記一方又は他方の薄膜の格子定数ａ１又はａ２の値を、格子定数及び膜厚を縦軸及び横軸とする一致度分布図、例えば図５に示すような等高線図様の分布図あるいは図６に示しすように一致度に応じて無段階に輝度、又は色及び色濃度が変化する分布図上の点（図５中に丸印で示す）として補助記憶装置４０５（分布データ記憶手段）に記憶蓄積させていき、複数のサンプル基板１０についての測定・解析がなされたとき、その分布データ記憶情報に基づいて、解析プロファイル及び実測プロファイルの一致の度合いの分布画像情報を生成する分布図作成手段の機能と、この分布図作成手段の機能により生成された分布画像情報を表示部４０７又は外部に出力する出力手段の機能とをそれぞれ有している。
【００７０】
図５に一例で示すように、分布図作成手段としての測定解析制御装置４０は、補助記憶装置４０５内の分布データ記憶情報に基づいて、解析プロファイル及び実測プロファイルの一致の度合いが最も高い一致領域１０１と、前記一致の度合いが一致領域１０１よりも低く一致領域１０１を取り囲むように分布する少なくとも１つの環状分布領域、例えば４つの環状分布領域１０２、１０３、１０４とが視覚的に区別可能な表示態様で生成するようになっている。ここで、視覚的に区別可能な表示態様とは、例えば異なる色、濃さ、１つ以上の一致度等高線による表示領域の分割、図示するような塗りつぶし条件の相違（例えば異なる網掛け）であり、図６に示すように、色を一致度が高い内方側から低い外方側へと徐々に変化させ、結果的に内外での一致どの相違が判定できる無段階的な変化の態様であってもよい。なお、図６では多数の分布データの点を個々に視認できるようには図示せず、分布の度合いを示す同色（例えば黄色）同輝度の環状分布領域をグラデーションにより無段階に多数表示しているが、データ自体はこの分布図中の一座標を持つ点として記憶されることになる。また、図５に示すように、一致領域１０１（図６では最も輝度の高い白い楕円形部分）の中心が解析・実測プロファイルの完全な一致点であり、この点から離れるほどずれが大きいことになる。したがって、一致領域１０１を取り囲むように分布する４つの環状分布領域１０２、１０３、１０４（図６では白い楕円形部分を囲む外側ほど徐々に黒くなる環状部分）は、両プロファイルが僅かにずれている状態を示す内側の領域１０２から、両プロファイルが部分的には近いものの他の部分に隙間があるような状態を示す中間の領域１０３、両プロファイル間に明らかに隙間があるような状態を示す外側の領域１０４へと順次外側になり、その外方では、両プロファイルが乖離して大きな隙間があるような状態（図６中では濃い黒の部分）となる。
【００７１】
なお、本実施の形態においては、一方の薄膜である井戸層の膜厚Ｌ１が決まると、他方の薄膜であるバリア層の膜厚Ｌ２が一義的に定まるため、図５に示すように、一致度分布図の一方の軸、例えば横軸に２つの薄膜の膜厚を併記している。
【００７２】
次に、上述の測定解析装置を用いて実施するＸ線回折測定解析方法の一実施形態について説明する。
【００７３】
図７は、本実施形態の測定解析装置で実行されるＸ線回折測定解析プログラムの概略の流れを示すフローチャートであり、このプログラムにより本発明の方法を実施することができる。
【００７４】
まず、設計情報記憶部４３への設計値情報の記憶やピーク特定条件記憶部４６への条件入力等の準備を行った後（ステップＳ１０）、測定を実行し（ステップＳ１１）、実測データ記憶部４１に実測データを記憶情報として順次格納する（ステップＳ１２）。
【００７５】
次いで、実測データ記憶部４１の記憶情報に基づいて、実測プロファイル作成部４２で実測プロファイルを作成する（ステップＳ１３）。
【００７６】
次いで、特定ピーク選択部４５でピーク特定条件記憶部４６の記憶条件に従って複数のうち孤立した次数の異なる特定のサテライトピークの発生領域を特定し（ステップＳ１４）、周期・平均格子定数算出部４７で超格子構造の平均組成に対応するこれら複数のサテライトピーク、例えばサテライトピークＰ20、Ｐ21の間隔から、超格子構造の積層周期Ｌを求めるともに（ステップＳ１５）、実測プロファイルの中の主要なサテライトピーク、例えば０次のサテライトピークＰ20と基板ピークＰ10との間隔から超格子構造の平均格子定数ａを求める（ステップＳ１６）。
【００７７】
次いで、初期値設定／変更部４４で前記４つの変数のうち第１組成薄膜又は第２組成薄膜のうち一方の薄膜の膜厚Ｌ１の初期値と、第１組成薄膜又は第２組成薄膜のうち一方の薄膜の格子定数ａ１の初期値とを、設計情報記憶部４３の記憶情報及び周期・平均格子定数算出部４７での算出結果Ｌ、ａに基づいて、それぞれ設定する（ステップＳ１７；設定段階）。この初期値は、実測プロファイルから求めた超格子構造の積層周期Ｌと平均格子定数ａが設計値に近い場合には、設計値と同一の値であるが、実測プロファイルから求めた超格子構造の積層周期Ｌと平均格子定数ａが設計値からずれている場合には、前記微調整量を単位としてそのずれ側に多少加減される。
【００７８】
次いで、変数算出部４８で、超格子構造の周期Ｌと、超格子構造の平均格子定数ａと、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値とから、第２組成薄膜又は第１組成薄膜の膜厚のうち他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値をそれぞれ算出する（ステップＳ１８；変数算出段階）。
【００７９】
次いで、解析プロファイル作成部４９で、超格子構造の周期Ｌ、超格子構造の平均格子定数ａ、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値、並びに前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値に基づいて、前記理論計算により解析プロファイルを作成する（ステップＳ１９；解析プロファイル作成段階）。
【００８０】
次いで、比較部５１で、作成された解析プロファイル及び実測プロファイルを比較して、両プロファイルの一致の度合いを調べる（ステップＳ２０；比較段階）。
【００８１】
次いで、条件判定部５２で、０次のサテライトピークについての一致の度合いが予め定められた一致基準を満たすか否かを判定し（ステップＳ２１）、一致の度合いが予め定められた一致基準を満たすときは、次数の異なるサテライトピーク、例えば−１次サテライトピークについての一致の度合いが予め定められた一致基準を満たすか否かを判定する（ステップＳ２２）。勿論、必要により、２次サテライトピーク（例えば−２次サテライトピーク）についての一致度合いを判定するステップ（ステップＳ２３）があってもよい。
【００８２】
ここで一致の度合いが予め定められた一致基準を満たす場合、そのときの前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値並びに前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値をそれぞれの膜厚測定値及び格子定数測定値として測定結果を表示部４０７に出力する（ステップＳ２４）。
【００８３】
一方、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれか１つの設定値を初期値設定／変更部４４で変更して（ステップＳ２６）、再び変数算出部４８での変数算出段階に戻る（ステップＳ１８）。
【００８４】
次いで、解析プロファイルを再度作成し（ステップＳ１９）、複数のサテライトピークのうち超格子構造の平均格子定数ａに対応する０次のサテライトピークＰ20について、比較部５１でのチェックを再度実行する（ステップＳ２１）。このように、最初に０次のサテライトピークＰ20のピーク点について、解析プロファイルを縦軸方向及び横軸方向でそれぞれ実測プロファイルに合わせるように、初期値設定／変更部４４での設定値の変更を実行させる。このとき、例えば、最初に０次ピーク強度を調整し、また、特定のサテライトピークＰ20、Ｐ21の角度間隔が解析プロファイルと実測プロファイルで一致していなければ、周期Ｌを加減してその間隔を近付ける。このような設定値変更は、図８中の左右両側に示すような解析プロファイルの変化をもたらす。すなわち、周期Ｌを増加させると、同図右側に示すようにサテライトピーク同士が近付き、周期Ｌを低減させると、同図左側に示すようにサテライトピーク同士が離れる。
【００８５】
そして、０次のサテライトピークＰ20のピーク点について一致基準を満たした場合に、主要サテライトピークのうち他の次数のサテライトピーク、例えば１次サテライトピークＰ21について、解析プロファイル及び実測プロファイルの一致の度合いを確認し（ステップＳ２２、Ｓ２３）、一致していなければ、初期値設定／変更部４４での設定値の変更を実行させる（ステップＳ２６）。
【００８６】
この間、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値変更が所定回数に達したか否かを判別し（ステップＳ２５）、所定回数に達した場合には不良と判定し、その判定結果を表示部４０７に出力する（ステップＳ２７）。
【００８７】
また、条件判定部５２では、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して再び変数算出部４８での算出を再度実行させる際、超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整する（ステップＳ２６）。さらに、このように超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整するのに先立って、０次のサテライトピークＰ20のピーク点について一致基準を満たすまで、前記一方の薄膜の格子定数ａ１を微調整する（ステップＳ２６）。このような設定値変更は、図８中の上下に示すような解析プロファイルの変化をもたらす。すなわち、格子定数を増加させると、同図上側に示すように複数のサテライトピークに接する包絡線ＥＮが同図中の反時計方向に傾き、格子定数を低減させると、同図下側に示すように複数のサテライトピークに接する包絡線ＥＮが同図中の時計方向に傾くように、解析プロファイルが変化する。したがって、格子定数を加減すると、その影響は次数の多いサテライトピークほど強く現れることになる。
【００８８】
このような膜厚Ｌ１及び格子定数ａ１一連の設定変更により、解析プロファイルは、図９に実践で模式的に示すように、同図に破線で概形を示す実測プロファイルに対して、同図（ａ）に示す最初の初期状態段階から、同図（ｂ）に示す０次サテライトピークの一致段階を経て、同図（ｃ）に示すような一致基準を満たす段階へと収束的に変化することになる。
【００８９】
一方、測定結果を表示部４０７に出力するか（ステップＳ２４）、あるいは、不良判定結果を表示部４０７に出力するか（ステップＳ２７）した場合、条件判定部５２で、前記一方又は他方の薄膜の膜厚Ｌ１又はＬ２の値と、前記一方又は他方の薄膜の格子定数ａ１又はａ２の値とを、図５に例示したような格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として、補助記憶装置４０５に記憶しておき（ステップＳ２９）、この記憶情報の蓄積に基づいて解析プロファイル及び実測プロファイルの一致の度合いを示す分布図を後から作成できるようにしておく。
【００９０】
そして、分布図の作成を指示する操作入力がなされたか否かをチェックし（ステップＳ３０）、指示された場合には、補助記憶装置４０５に記憶蓄積された一致度のデータを基に分布図を作成して、その分布図上に、解析プロファイル及び実測プロファイルの一致の度合いが最も高い一致領域１０１と、一致の度合いが一致領域１０１よりも低く一致領域１０１を取り囲むように分布する少なくとも１つの環状分布領域１０２〜１０４とを上述したような視覚的に区別可能な表示態様で表示させ、一致度分布図上の点が、その分布図上の一致領域１０１及び環状分布領域１０２〜１０４のうちどの領域に属するかで前記一致の度合いを視覚的に識別可能にする（ステップＳ３２）。
【００９１】
また、測定終了を指示する操作入力がなされたか否かをチェックし、指示された場合には、処理を終了する（ステップＳ３３）。また、終了指示がなければ、ステップＳ１１に戻って、次の測定開始のため待機する。
【００９２】
このように、本実施形態の回折Ｘ線測定解析方法及びプログラムにおいては、実測プロファイルから超格子構造の周期Ｌと超格子構造の平均格子定数ａが求められ、これらの値と既知の情報とに基づいて、第１又は第２組成薄膜のうち一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値が４つの変数のうち２つの変数の初期値として設定され、これらの値を基に残りの２つの変数が算出されると、それまでに得られた情報に基づいて、理論計算により解析プロファイルが作成され、その解析プロファイルと実測プロファイルが比較され、その比較結果に応じて、一致の度合いが一致基準を満たすときには、そのときの４つの変数が測定値とされ、基準を満たさないときには、２つの変数である前記初期値のいずれかを変更して再び変数算出段階に戻る。したがって、適切な２つの初期値設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【００９３】
また、本実施形態においては、設定段階において、積層された超格子構造に関する設計値情報に基づき、実測プロファイルから周期性のある孤立したピークを抽出して主要サテライトピークを特定するので、小ピーク、ノイズ等の微視的プロファイル部分の悪影響をより有効に回避でき、しかも、処理を簡単化することができる。
【００９４】
しかも、前記設定段階において、実測プロファイルから求めた超格子構造の周期Ｌ及び平均格子定数ａと、前記一方及び他方の薄膜の膜厚Ｌ１、Ｌ２及び格子定数ａ１、ａ２の設計値情報とに基づいて、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値を設定するので、実測プロファイルや設計情報等から得られる情報を基に、より得るべき測定値に近い初期値が設定できる。
【００９５】
また、本実施形態においては、比較段階において、複数のサテライトピークのうち超格子構造の平均格子定数ａに対応する０次のサテライトピークＰ20について、解析プロファイルを縦軸方向及び横軸方向でそれぞれ実測プロファイルに合わせ、主要サテライトピークのうち他の次数のサテライトピークＰ21等について、解析プロファイル及び実測プロファイルの一致の度合いを調べるので、４つの変数の変化に対する解析プロファイルの変化の特性を考慮して、一致度の比較処理を簡単化することができる。
【００９６】
また、判定段階において、一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して再び変数算出段階に戻る際、超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整するので、実測プロファイルから正確に得られる周期Ｌを固定することで、一致点への収束を早めることができる。
【００９７】
さらに、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して再び変数算出段階に戻る際、超格子構造の周期Ｌを一定にしつつ膜厚Ｌ１の設定値を微調整するのに先立って、前記一方の薄膜の格子定数ａ１を微調整するので、膜厚Ｌ１と格子定数ａ１の変化における相関性を考慮して、一致点への収束を早めることができる。
【００９８】
また、本実施形態においては、前記判定段階で、一方又は他方の薄膜の膜厚値Ｌ１又はＬ２と、一方又は他方の薄膜の格子定数値ａ１又はａ２とを、格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶しておき、その記憶情報に基づいて解析プロファイル及び実測プロファイルの一致の度合いを示す分布図を作成するとともに、その分布図上に、解析プロファイル及び実測プロファイルの一致の度合いが最も高い一致領域１０１と、一致の度合いが一致領域１０１よりも低く一致領域１０１を取り囲むように分布する環状分布領域１０２〜１０４とを視覚的に区別される表示態様で表示し、分布図上の点が、その分布図上の一致領域１０１及び環状分布領域１０２〜１０４のうちどの領域に属するかで一致の度合いを視覚的に識別可能に表示するので、解析プロファイル及び実測プロファイルの一致の度合いがどの程度かを容易に見て取れることになり、ユーザーの迅速・的確な大局的判断が可能になる。
【００９９】
なお、上述の実施形態においては、異種の薄膜を交互に積層した超格子構造を有する結晶構造として、通信分野で使用される高出力の半導体レーザの活性層を例としていたが、半導体レーザ以外の異種の薄膜を交互に積層した他の超格子構造にも適用できることはいうまでもない。
【産業上の利用可能性】
【０１００】
本発明は、以上説明したように、サテライトピーク間にある小さなピークやノイズはフィッティングに用いることなく、既知の情報及び実測データに基づいて選択した特定のサテライトピークについて、適切な２つの変数の設定値設定とこれら初期値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行するので、超格子構造を反映する主要サテライトピークについて迅速で且つ有効なフィッティングを行ない、必要にして十分な精度の解析が行なえるＸ線回折測定解析方法及びプログラムを実現することができ、さらに、一致度の分布表示による大局的判断を可能にできるＸ線回折測定解析方法及びプログラムを実現できるという効果を奏するものであり、異種の薄膜を交互に積層した超格子構造の測定・解析に用いるＸ線回折測定解析方法及びプログラム全般に有用である。
【図面の簡単な説明】
【０１０１】
【図１】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法を実施する装置の測定系の概略構成図である。
【図２】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法を実施する装置の全体の概略構成図である。
【図３】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法を実施する装置の制御装置の要部機能ブロック図である。
【図４】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で作成される解析・実測プロファイルの一例を例示する回折Ｘ線のプロファイル図で、縦軸は回折Ｘ線強度を、横軸は回折角度をそれぞれ示している。
【図５】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で作成される解析・実測プロファイルの一致度分布図で、縦軸は一方又は他方の薄膜の格子定数を、横軸はその薄膜の膜厚としている。
【図６】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で作成される解析・実測プロファイルの異なる態様の一致度分布図で、縦軸は一方又は他方の薄膜の格子定数を、横軸はその薄膜の膜厚としている。
【図７】本発明の第１の実施の形態に係るＸ線回折測定解析プログラムの概略の流れを示すフローチャートである。
【図８】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で行われるパラメータ変更とそれによる解析プロファイルの変化の関係を説明する説明図である。
【図９】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で行われるフィッティング処理とそれにより解析プロファイルが実測プロファイルに近付く状態を説明する説明図である。
【符号の説明】
【０１０２】
１ 第１回転台
２ 第２回転台
３、４ 駆動手段
６、７ 回転位置検出器
１０ サンプル基板（基板）
２０ Ｘ線源
２１ Ｘ線管
２２ 湾曲ミラー
２３ モノクロメータ
３０ Ｘ線カウンタ
４０ 測定解析制御装置（分布図作成手段、出力手段）
４１ 実測データ記憶部（実測データ記憶手段）
４２ 実測プロファイル作成部（実測プロファイル作成手段）
４３ 設計情報記憶部（設計情報記憶手段）
４５ 特定ピーク選択部
４６ ピーク特定条件記憶部（特定条件記憶手段）
４７ 周期・平均格子定数算出部（周期算出手段、平均格子定数算出手段）
４８ 変数算出部（変数算出手段）
４９ 解析プロファイル作成部（解析プロファイル作成手段）
５１ 比較部（シミュレーション処理部、誤差解析部）
５２ 条件判定部（判定手段、シミュレーション処理部）
１０１ 一致領域
１０２、１０３、１０４ 環状分布領域
４０１ ＣＰＵ
４０２ ＲＯＭ
４０３ ＲＡＭ
４０４ 入出力インターフェース回路
４０５ 補助記憶装置（分布データ記憶手段）
４０６ 設定・入力部
４０７ 表示部（出力手段）
【技術分野】
【０００１】
本発明は、Ｘ線回折測定解析方法及びプログラム、特に異種の薄膜を交互に積層した超格子構造の測定・解析に好適なＸ線回折測定解析方法及びプログラムに関する。
【０００２】
異種の薄膜を交互に積層した超格子構造を有する結晶構造は、通信分野で使用される高出力の半導体レーザの活性層を代表例として多くの用途に適用されている。
【背景技術】
【０００３】
近時、薄膜製造に対する要求精度が従前より更に厳しくなる中で、薄膜を測定・解析対象とする場合にＸ線回折を利用する高精度の測定解析法が多用されてきており、その中でも、異種の薄膜を交互に積層した超格子構造の測定・解析には、いわゆるロッキングカーブ測定・解析が利用されている。
【０００４】
ロッキングカーブ測定・解析を行なう従来のＸ線回折測定解析方法においては、基板に入射するＸ線の入射角を変化させていくと、その角度が基板結晶やその基板上に交互に積層された超格子構造の薄膜のブラッグ角にそれぞれ一致するところでＸ線が相互に強め合うように回折・出射することで、回折Ｘ線強度が回折角度の変化に応じてその基板及び超格子構造に由来する特定のパターンで変化する。したがって、回折角度を変えて結晶構造情報を含む回折Ｘ線強度変化のパターンを測定することで、そのプロファイル解析から超格子構造の異種薄膜の積層周期や膜厚や格子定数等の値を知ることができる。
【０００５】
前記測定に際しては、例えばＸ線源からのＸ線ビームを試料に照射し、回折Ｘ線の強度をＸ線カウンタでカウントすると、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされるＸ線の実測プロファイル中に、基板よりの回折Ｘ線に由来する基板ピークや、その基板上に積層された超格子構造に由来する複数のサテライトピークが現れる。
【０００６】
一方、Ｘ線回折現象による上記プロファイルの解析に際しては、基板上に超格子構造を持つ結晶構造体の既知の情報（基板の結晶構造情報や基板及び超格子構造の設計情報）に基づいて、Ｘ線のプロファイルを特徴付ける４つの変数をパラメータとして実測プロファイルに対応する解析プロファイルを理論計算で求め、前記４つの変数を変更しながらフィッティング処理を行なって、最終的に実測プロファイルに最も近い一致を見た時点の前記４つの変数を測定値としている。このように、Ｘ線回折で測定された実測プロファイルを解析することにより、結晶構造を数値的に把握することができる。
【０００７】
より具体的には、例えば、歪量子井戸半導体レーザの活性層を構成する量子井戸層の格子定数及び膜厚、並びにバリア層の格子定数及び膜厚の４つの値をパラメータとして、理論計算で解析プロファイルを求め、その解析プロファイルのパターンが実測プロファイルのパターンと、より良く一致するように理論計算のパラメータを変化させて、両プロファイルのフィッティングを行なう。そして、フィッティングにより最も良く一致する解析プロファイルの４つのパラメータ、すなわち量子井戸層の格子定数及び膜厚並びにバリア層の格子定数及び膜厚を測定値とする。なお、実測プロファイルと解析プロファイルを得る方法は、それぞれ既に確立されており学術界・産業界に広く普及している。
【０００８】
一方、上述のフィッティング処理は、対象によって異なるが、例えば最小二乗法近似によるフィッティングが良く知られている。この最小二乗法近似は、実測プロファイルと解析プロファイルの差の二乗誤差が最小になるようにフィッティングを行なうものであり、プロファイルの全角度範囲にわたって両プロファイルの差の二乗値を均等に加算し平方根をとりサンプル数で除算するという単純なアルゴリズムを採用でき、計算機を用いて容易に自動化することができる。また、パラメータの初期値を適切に設定すれば、コンピュータプログラムによって自動的にもっとも良く一致する点を見つけ出すことができるので、比較的広いパラメータ範囲を探索するような場合に適している。
【０００９】
しかし、実測プロファイルにはノイズ等が多く含まれており、パターン認識処理も容易ではないので、ソフトウェアにより最も一致した点を見つけ出すために、例えば、実測プロファイルデータに対してスムージング（平滑化）処理を施してノイズを含む実測プロファイルから超格子構造の特徴を示す主要なピーク（例えばサテライトピーク）についてのフィッティングを実行し、次いで、このスムージングを徐々に減らしていきながら、プロファイルの細かい特徴部分についてのフィッティングを行なうという手法が採られている（例えば、特許文献１、非特許文献１参照）。
【特許文献１】特表２００２−５３２７０３号公報
【非特許文献１】ＰＡＮａｌｙｔｉｃａｌ社の「Ｐｒｅｓｓ Ｒｅｌｅａｓｅｓ」、「New X'Pert Epitaxy software from Philips Analytical makes rocking-curve analysis easier than ever」［２００４年１２月１０日検索］、インターネットＵＲＬ＜http://www.panalytical.com/index.cfm?pid=21&itemid=153&contentitemid=12＞。
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、上記従来のＸ線回折測定解析方法及びプログラムにあっては、最小二乗法近似を用いているため、スムージングによる主要ピークについての第１段階のフィッティングとスムージングを減らしながらの第２段階のフィッティングとを行ったとしても、小さなピークやノイズまでフィッティングの対象としてしまうので、条件次第で本来１つしかない一致点が複数得られてしまうという結果を招きかねないという問題があった。
【００１１】
一方、基板の結晶に関する格子定数、並びに超格子構造の各種薄膜の配置、膜厚、格子定数等の設計値は既知であり、更に、実測プロファイルには、基板ピークと共に、超格子構造の平均格子定数に対応する０次のサテライトピークを中心に超格子構造の周期（Ｌ）、すなわち、交互に積層される薄膜の膜厚の和に相当する回折角度間隔を隔てて複数の他次数のサテライトピークが複数現れる。すなわち、実測プロファイルの形状的特徴とパラメータの間の特定の関係が明らかになっている。
【００１２】
したがって、このように実測プロファイルの形状的特徴と前記４つの変数（パラメータ）の明らかな関係を利用し、４つの変数を２つの独立変数に集約すれば、重要情報を有するプロファイル形状部分（特定のサテライトピーク）のみに着目してフィッティングを行なうことができる。このようにすれば、小ピーク、ノイズ等の微視的プロファイル部分の悪影響を抑えた簡便且つ確実な解析が可能になる。
【００１３】
そこで、本発明は、サテライトピーク間にある小さなピークやノイズはフィッティングに用いずに超格子構造を反映する主要サテライトピークについてフィッティングを行ない、既知の情報及び実測データに基づいて選択した特定のサテライトピークのフィッティングによって、必要にして十分な精度の解析が行なえるＸ線回折測定解析方法及びプログラムを実現することを目的とするものであり、併せて、４パラメータの一致状態を一般性を失うことなく大局的に見通しよく理解できるような２次元表示を実現することを目的とする。
【課題を解決するための手段】
【００１４】
本願の発明者は、解析プロファイルの４つの変数（パラメータ）を変化させて解析プロファイルを実測プロファイルに一致させることに代えて、まず、超格子物理の自然法則から、これら変数の従属関係を明確にすることで、４つの変数を２つの独立変数に集約できるとの知見を得、簡便、正確で局所最適値に陥る危険性のない優れた測定解析方法を創作したものである。また、これにより、２変数による鳥瞰図的な４変数の全体像を大局的に見通しよく表示することも可能となる。
【００１５】
すなわち、本発明は、上記目的達成のため、（１）基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、前記４つの変数を変化させながら比較し、最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析方法において、前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求め、且つ前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求めるとともに、前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１）の初期値を設定する設定段階と、前記超格子構造の周期（Ｌ）と前記超格子構造の平均格子定数（ａ）と前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を算出する変数算出段階と、前記超格子構造の周期（Ｌ）、前記超格子構造の平均格子定数（ａ）、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値、並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値に基づいて、前記理論計算により前記解析プロファイルを得る解析プロファイル作成段階と、前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較段階と、前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して再び変数算出段階に戻る判定段階と、を含むことを特徴とするものである。
【００１６】
本発明の方法では、まず、実測プロファイルから超格子構造の周期（Ｌ）と超格子構造の平均格子定数（ａ）が求められ、これらの値と既知の情報とに基づいて、第１又は第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値が４つの変数のうち２つの変数の初期値として設定され（設定段階）、これらの値を基に残りの２つの変数が算出されると（変数算出段階）、それまでに得られた情報に基づいて、理論計算により解析プロファイルが作成される（解析プロファイル作成段階）。そして、その解析プロファイル及び実測プロファイルが比較され（比較段階）、両プロファイルの一致の度合いが一致基準を満たすときには、そのときの４つの変数が測定値とされ、一致の度合いが前記基準を満たさないときには、２つの変数である設定値のいずれかを変更して再び変数算出段階に戻る（判定段階）。したがって、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、適切な２つの初期値の設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【００１７】
この発明の方法においては、（２）前記設定段階において、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイルから周期性のある孤立したピークを抽出して前記主要サテライトピークを特定するのが好ましい。これにより、小ピーク、ノイズ等の微視的プロファイル部分の悪影響をより有効に回避でき、しかも、処理を簡単化できることになる。
【００１８】
この発明の方法においては、さらに、（３）前記設定段階において、前記実測プロファイルから求めた前記超格子構造の周期（Ｌ）及び平均格子定数（ａ）と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値を設定するのがよい。実測プロファイルや設計情報等から得られる情報を基に、より得るべき測定値に近い初期値が設定できるからである。
【００１９】
また、この発明の方法においては、（４）前記比較段階において、前記複数のサテライトピークのうち前記超格子構造の平均格子定数（ａ）に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせ、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べるのが好ましい。これにより、４つの変数の変化に対する解析プロファイルの変化の特性を考慮して、一致度の比較処理を簡単化できることになる。
【００２０】
また、この発明の方法においては、（５）前記判定段階において、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのがよい。超格子構造の周期（Ｌ）は実測プロファイルから正確に得られるからであり、この周期（Ｌ）を固定することで、一致点への収束を早めることができる。
【００２１】
この場合、さらに、（６）前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのに先立って、前記一方の薄膜の格子定数（ａ１）を微調整するのが望ましい。膜厚（Ｌ１）と格子定数（ａ１）の変化における相関性を考慮して、一致点への収束を早めることができるからである。
【００２２】
また、この発明の方法においては、（７）前記判定段階で、前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶しておき、該記憶情報に基づいて前記解析プロファイル及び前記実測プロファイルの一致の度合いを示す分布図を作成するとともに、該分布図上に、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とを視覚的に区別される表示態様で表示し、前記分布図上の点が、該分布図上の前記一致領域及び環状分布領域のうちどの領域に属するかで前記一致の度合いを視覚的に識別可能に表示するのが好ましい。これにより、解析プロファイル及び実測プロファイルの一致の度合いがどの程度かを容易に見て取れることになり、ユーザーの迅速・的確な判断が可能になる。
【００２３】
本発明は、また、（８）基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、前記４つの変数を変化させながら比較し、最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、コンピュータにより、前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求める周期算出手段（１１）と、前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求める平均格子定数算出手段と、前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１）の初期値を設定する初期値設定手段と、前記超格子構造の周期（Ｌ）と前記超格子構造の平均格子定数（ａ）と前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値をそれぞれ算出する変数算出手段と、前記超格子構造の周期（Ｌ）、前記超格子構造の平均格子定数（ａ）、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値に基づいて、前記理論計算により前記解析プロファイルを作成する解析プロファイル作成手段と、前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較手段と、前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値、並びに前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段と、の機能を実現させることを特徴とするものである。
【００２４】
このプログラムにより、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、適切な２つの初期値の設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【００２５】
本発明のプログラムにおいては、（９）前記初期値設定手段が、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイルから前記超格子構造の周期（Ｌ）に対応する周期性を有する孤立した複数のピークを抽出して前記主要サテライトピークを特定するのがよい。このようにすると、小ピーク、ノイズ等の微視的プロファイル部分の悪影響をより有効に回避でき、しかも、処理を簡単化できることになる。
【００２６】
また、この発明のプログラムにおいては、（１０）前記初期値設定手段が、前記実測プロファイルから求めた前記超格子構造の周期（Ｌ）及び平均格子定数（ａ）と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の初期値を設定するのがよい。この構成により、実測プロファイルや設計情報等から得られる情報を基に、より得るべき測定値に近い初期値が設定可能となる。
【００２７】
また、この発明のプログラムにおいては、（１１）前記比較手段が、前記複数のサテライトピークのうち前記超格子構造の平均格子定数（ａ）に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせるシミュレーションを実行し、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べるのが好ましい。これにより、４つの変数の変化に対する解析プロファイルの変化の特性を考慮して、一致度の比較処理を簡単化できることになる。
【００２８】
この発明のプログラムにおいては、さらに、（１２）前記判定手段が前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのがよい。これにより、超格子構造の周期（Ｌ）は実測プロファイルから正確に得られるからであり、この周期（Ｌ）を固定することで、一致点への収束を早めることが可能となる。
【００２９】
この場合、さらに、（１３）前記判定手段が前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期（Ｌ）を一定にしつつ前記一方の薄膜の膜厚（Ｌ１）の設定値を微調整するのに先立って、前記一方の薄膜の格子定数（ａ１）を微調整するのが望ましい。これにより、膜厚（Ｌ１）と格子定数（ａ１）の変化における相関性に基づき、一致点への収束を早めることが可能となる。
【００３０】
また、この発明のプログラムにおいては、（１４）前記コンピュータに、前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶する記憶手段と、該記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを示す分布図の画像情報を生成する分布図作成手段と、の機能を実現させるのがよい。このプログラムにより、複数の測定値の分布を、解析プロファイル及び実測プロファイルの一致の度合いの分布として視覚的に容易に把握可能となる。
【００３１】
この場合、さらに、（１５）前記分布図作成手段が、前記記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とを視覚的に区別される表示態様で表示するのが好ましい。これにより、解析プロファイル及び実測プロファイルの一致の度合いがどの程度かを容易に見て取れることになり、ユーザーの迅速・的確な大局的判断が可能になる。
【００３２】
本発明のプログラムは、あるいは、（１６）基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を異なる角度で入射したとき測定情報に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の実測プロファイルと、４つの変数として前記第１組成薄膜の格子定数（ａ１）及び膜厚（Ｌ１）並びに前記第２組成薄膜の格子定数（ａ２）及び膜厚（Ｌ２）を用いた理論計算に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の解析プロファイルとを、前記解析プロファイルの４つの変数を変化させながら比較し、最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、ｎを正の整数とするとき、ｎ個の前記角度の値（ｘ１，ｘ２，・・・ｘｎ）及び該角度に対応するｎ個の回折Ｘ線強度の測定値（ｇ１，ｇ２，・・・ｇｎ）を順次メモリに記憶させる実測データ記憶手段（４１）と、前記メモリに記憶されたデータに基づいて前記角度に応じて変化する回折Ｘ線強度の実測プロファイルを表示出力可能な所定のプロファイル表示形式で作成する実測プロファイル作成手段（４２）と、前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求める周期算出手段（４７）と、前記実測プロファイルの中に超格子構造の平均組成に対応して現れるサテライトピークと前記実測プロファイルの中に前記基板に由来して現れる基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求める平均格子定数算出手段（４７）と、前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１又はＬ２）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１又はａ２）の初期値を、前記４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）のうち２つの変数の初期値として設定する初期値設定手段（４４）と、前記超格子構造の周期（Ｌ）と、前記超格子構造の平均格子定数（ａ）と、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値とから、前記４つの変数（ａ１、Ｌ１、ａ２、Ｌ２）のうち残りの２つの変数である前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値をそれぞれ算出する変数算出手段（４８）と、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値と、前記変数算出手段（４８）で算出された前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値とを、前記４つの変数として、前記複数のサテライトピークのうち次数の異なる特定のサテライトピークに対応する角度範囲を特定する前記ｎより少ない正の整数であるｍ個の角度値（ｘｉ，ｘｊ，・・・ｘi+m-1；但し、ｉ、ｊはそれぞれｍより小さい正の整数）について、回折Ｘ線強度の値（ｆｉ，ｆｊ，・・・ｆi+m-1）を前記理論計算により算出して前記解析プロファイルを作成する解析プロファイル作成手段（４９）と、前記解析プロファイル及び前記実測プロファイルを前記複数のサテライトピークのうち前記次数の異なる特定のサテライトピークにつき比較して、両プロファイルの一致の度合いを調べる比較手段（５１）と、前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記解析プロファイルの４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）を前記一方の薄膜の膜厚及び格子定数並びに前記他方の薄膜の膜厚及び格子定数の測定値として出力し、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段（５２）と、の機能をコンピュータを使用して実現させることを特徴とするものである。
【００３３】
このプログラムにより、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、適切な２つの初期値の設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【発明の効果】
【００３４】
本発明によれば、サテライトピーク間にある小さなピークやノイズはフィッティングに用いることなく、既知の情報及び実測データに基づいて選択した特定のサテライトピークについて、適切な２つの変数の初期値設定とこれら初期値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行するので、超格子構造を反映する主要サテライトピークについて迅速で且つ有効なフィッティングを行ない、必要にして十分な精度の解析が行なえるＸ線回折測定解析方法及びプログラムを実現することができる。また、４パラメータの一致状態を一般性を失うことなく２つのパラメータに集約することで、一致状態を一般性を失うことなく大局的に見通しよく理解できるような２次元表示を実現することができる。
【００３５】
すなわち、本発明は、解析プロファイルの４つの変数の従属関係を明確にすることで、４つの変数を２つの独立変数に集約させ、簡便、正確で局所最適値に陥る危険性のない優れた測定解析方法を創作したものであり、これにより、２変数による鳥瞰図的な４変数の全体像を大局的に見通しよく表示することも可能としたものである。
【発明を実施するための最良の形態】
【００３６】
以下、本発明の好ましい実施の形態について、図面を参照して説明する。
【００３７】
図１乃至図９は本発明の１つの実施の形態に係るＸ線回折測定解析方法を実施するための測定解析装置を示す図であり、この装置は本発明のＸ線回折測定解析プログラムに従って後述する処理を実行する。
【００３８】
まず、その構成について説明する。
【００３９】
図１及び図２に模式的に示すように、この測定装置では、第１回転台１に支持されたサンプル基板１０に対し、Ｘ線源２０からのＸ線ビームが照射されるとともに、サンプル基板１０で回折したＸ線を第２回転台２に支持されたＸ線カウンタ３０で検出し、そのＸ線強度を測定する光学系が構成されている。第１回転台１は、サンプル基板１０へのＸ線ビームの入射角度が変化するよう角速度ωで回転（往復回動）可能であり、第２回転台２は第１回転台１と同一中心軸回りに第１回転台１とは異なる角速度で回転可能である。これら回転台１、２は、それぞれに対応する公知の駆動手段３、４によって回転駆動されるとともに、回転中のサンプル基板１０の基板面と入射Ｘ線ビームのなす角度θ（視射角）に対して、Ｘ線ビームの照射軸線（図１に一点鎖線で示す）に対しＸ線カウンタ３０のＸ線受光軸線のなす角度が常に角度２θに保たれるように回転制御されるようになっている。すなわち、サンプル基板１０に対するＸ線ビームの入射角変化に対して常に回折出射角方向にＸ線カウンタ３０が配置された状態を保つようになっている。
【００４０】
サンプル基板１０は、基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造体で、例えばインジウム燐（以下、ＩｎＰという）基板上に異種のインジウムガリウム砒素燐（ＩｎＧａＡｓＰ）の薄膜を積層したもので、半導体レーザの活性層となるバリア層と量子井戸層の薄膜（例えば６０オングストローム程度）を積層した多重量子井戸（ＭＱＷ）構造を有している。
【００４１】
Ｘ線源２０は、例えばＸ線管２１、湾曲ミラー２２及びモノクロメータ２３によって構成されており、Ｘ線管２１で発生したＸ線束を湾曲ミラー２２によってサンプル基板１０の測定面上に収束させるとともに、モノクロメータ２３によって特定波長域に単色化するようになっている。また、Ｘ線カウンタ３０は、例えばシンチレーションカウンタで構成されており、回折Ｘ線の強度に応じたカウント値を出力するようになっている。
【００４２】
図２に示すように、Ｘ線源２０のＸ線管２１の駆動回路５並びに第１回転台１及び第２回転台２の駆動手段３、４はそれぞれ測定解析制御装置４０に接続されており、この測定解析制御装置４０によって上述のような測定系の動作制御がなされる。また、Ｘ線カウンタ３０の信号処理回路３１からのカウント情報並びに第１回転台１及び第２回転台２の回転位置検出器６、７の検出情報がそれぞれ測定解析制御装置４０に取り込まれるようになっており、測定解析制御装置４０はこれら入力情報及び設定・入力部４０６からの設定入力又は操作入力等に基づいて、所定のＸ線回折測定解析プログラムを実行し、公知のディスプレイデバイスからなる表示部４０７に横軸を回折角度（秒単位）とし縦軸を回折Ｘ線強度とする回折Ｘ線図形を表示出力したり、その表示画像情報を外部に出力したりするようになっている。
【００４３】
測定解析制御装置４０は、詳細のハードウェア構成を図示しないが、概略、例えばＣＰＵ（Central Processing Unit）４０１、ＲＯＭ（Read Only Memory）４０２、ＲＡＭ（Random Access Memory）４０３及び入出力インターフェース回路４０４に加え、ディスクドライブ等の補助記憶装置４０５を内蔵して構成され得るもので、前記設定・入力部４０６及び表示部４０７が付設されている。また、この測定解析制御装置４０は、ＣＰＵ４０１がＲＯＭ４０２及び補助記憶装置４０５に記憶された制御プログラムに従って、ＲＡＭ４０３との間でデータを授受しながら、Ｘ線カウンタ３０の信号処理回路３１からのカウント情報、回転位置検出器６、７の検出情報及び設定・入力部４０６からの入力情報等に基づいて、回折Ｘ線強度や回折角度を把握する演算処理を実行するとともに、図３に機能ブロックで示す各部の機能を発揮するようになっている。
【００４４】
図３に示すように、測定解析制御装置４０においては、回折Ｘ線強度とそのときの回折角度の情報が順次実測データ記憶部４１に記憶される。すなわち、上述のようにサンプル基板１０に対するＸ線ビームの入射角度を変化させながら、サンプル基板１０の結晶構造にＸ線を入射し回折されたＸ線をＸ線カウンタ３０で検出して得られるＸ線回折測定情報が、実測データ記憶部４１に順次記憶される。このとき、具体的には、例えば、ｎ（正の整数）個の角度の値ｘ１，ｘ２，・・・ｘｎ及びその角度に対応するｎ個の回折Ｘ線強度の測定値ｇ１，ｇ２，・・・ｇｎを順次メモリに記憶させることになる。
【００４５】
実測プロファイル作成部４２は、この実測データ記憶部４１の記憶データに基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる実測プロファイル、すなわち実測データから得られる回折Ｘ線図を公知の図形化処理によって作成する。この実測プロファイルは、Ｘ線回折測定結果として、横軸に角度、縦軸にＸ線強度をとったグラフで表示され、例えば図４に示すように、その横軸の回折角度が、サンプル基板１０のうち基板層の結晶の格子定数に対応した最大の鋭いピーク（基板ピーク）Ｐ10を中心に相対角度（秒）で示され、縦軸の回折Ｘ線強度が毎秒当たりのカウント数（counts／sec）で示されている。
【００４６】
実際の実測プロファイルは、図４に例示するように、サンプル基板１０の基板層よりの回折Ｘ線に由来する基板ピークＰ10の他に、そこに積層された超格子構造に由来する複数のサテライトピークＰ20（０次サテライトピーク）、Ｐ21（−１次サテライトピーク）、Ｐ22（＋１次サテライトピーク）、Ｐ23（＋２次サテライトピーク）等が超格子における異種薄膜の積層周期Ｌに対応する角度間隔で現れる。なお、図４中の他のピークは超格子構造に含まれる他の薄膜、例えば半導体レーザの光ガイド層を構成する複数層のインジウムガリウム砒素燐（ＩｎＧａＡｓＰ）の薄膜層（例えば１００オングストローム程度）からの回折Ｘ線等に基づくものである。
【００４７】
サンプル基板１０となる超格子構造の設計情報（交互に積層される異種薄膜の積層周期、各層の薄膜の膜厚及び格子情報、配置等）は、設計情報記憶部４３に予め設定記憶されている。
【００４８】
この設計情報記憶部４３の記憶情報は、実測プロファイル作成部からの実測プロファイル情報と共に特定ピーク選択部４５に取り込まれ、この特定ピーク選択部４５において、積層された超格子構造に関する設計値情報に基づき、実測プロファイル中で前記超格子構造の周期Ｌに対応する周期性を有する孤立した複数のピークの発生領域が特定される。
【００４９】
ここで特定のピークの発生領域を特定する条件は、ピーク特定条件記憶部４６に記憶されている。この条件は、例えば超格子を構成する異種の薄膜の格子定数とサンプル基板１０の基板層の結晶の格子定数との大小関係、並びに、それに応じて異なるサテライトピークの現れ方の特性に関する情報に基づいて、サテライトピークのうち孤立したサテライトピーク（例えば図４に示すサテライトピークＰ20、Ｐ21のように他の大きなピークとの干渉が少なく独立してピーク波形を形成しているもの）が発生する蓋然性が高い領域を回折角度範囲で特定する条件である。例えば、第１組成薄膜の格子定数が基板層（ＩｎＰ）の格子定数より大きく、第２組成薄膜の格子定数が基板層の格子定数と同等の設計仕様の場合には、通常、孤立したサテライトピークが基板ピークの左側（マイナス側角度域で超格子の積層周期に応じた所定角度範囲の特定領域内）に複数現れ、あるいは、第１組成薄膜の格子定数が基板層（ＩｎＰ）の格子定数と同等で、第２組成薄膜の格子定数が基板層の格子定数より小さい設計仕様の場合には、通常、孤立したサテライトピークが基板ピークの右側（プラス側角度域で超格子の積層周期に応じた所定角度範囲の特定領域内）に複数現れ、あるいはまた、第１組成薄膜の格子定数が基板層（ＩｎＰ）の格子定数より大きく、第２組成薄膜の格子定数が基板層の格子定数より小さい設計仕様の場合には、通常、孤立したサテライトピークが基板ピークの両側にそれぞれ複数現れる。そして、これらの孤立したサテライトピークがフィッティングに有効なサテライトピークとなる一般的傾向がある。このような条件に基づいて、フィッティングに有効な特定サテライトピークの発生領域を判定することができる。
【００５０】
また、特定ピーク選択部４５で孤立したサテライトピークを抽出する処理は、例えば特定のサテライトピークの発生領域内でのサテライトピークの間隔より十分に狭い所定角度範囲の単位で平均の回折Ｘ線強度を求めて各次数のサテライトピークを抽出することができる。この際、サテライトピークの間隔は設計情報に基づいて概ね予測でき、各次数のサテライトピークを特定（角度範囲で特定）するためのサンプリングは比較的少数で足りる。なお、０次や−１次、＋１次のサテライトピークが基板ピークや他の層に起因するピークに埋もれている場合は、他の孤立したサテライトピークを使用する。
【００５１】
特定ピーク選択部４５でこのようにピーク特定条件記憶部４６の記憶情報に基づいて特定された主要サテライトピークに関するプロファイル情報は、周期・平均格子定数算出部４７に送られ、ここで、第１組成薄膜及び第２組成薄膜の積層周期、すなわち多重量子井戸構造の井戸層及びバリア層からなる超格子の積層周期Ｌと多重積層された超格子構造の平均格子定数ａが算出される。
【００５２】
周期・平均格子定数算出部４７においては、第１組成薄膜及び第２組成薄膜の積層周期Ｌが、実測プロファイル中の主要なサテライトピークの角度間隔から、具体的には、例えば特定の次数の異なるサテライトピークＰ20、Ｐ21が現れる回折角度θｉ、θｊと、Ｘ線の波長とから、次式（１）により直接的に精度良く算出できる。ただし、ｉ、ｊはそれぞれサテライトピークの次数である。
【００５３】
Ｌ＝（ｉ＋ｊ）λ／（２（ｓｉｎθｉ−ｓｉｎθｊ）） ・・・（１）
ここで、第１組成薄膜及び第２組成薄膜の積層周期Ｌに対して、第１組成薄膜（たとえば井戸層）の膜厚をＬ１とし、第２組成薄膜（例えばバリア層）の膜厚をＬ２とすると、これらは次式（２）で表すことができる。
【００５４】
Ｌ＝Ｌ１＋Ｌ２ ・・・（２）
したがって、２つのパラメータＬ１、Ｌ２には従属関係があり、パラメータとしてのＬ２を定めると、次式（３）のようにＬ１が定まり、逆に、Ｌ１を定めると自動的にＬ２が定まる。
【００５５】
Ｌ１＝Ｌ−Ｌ２ ・・・（３）
周期・平均格子定数算出部４７においては、また、実測プロファイル中の複数のサテライトピークのうち、超格子構造の平均組成である平均格子定数ａに対応する０次サテライトピークＰ20が現れる回折角度と、サンプル基板のＩｎＰ基板層の結晶による基板ピークP10が現れる回折角度との差から、平均格子定数ａが精度良く直接的に算出できる。この平均格子定数ａは、また、次式（４）で表すことができる。
【００５６】
ａ＝（ａ１・Ｌ１＋ａ２・Ｌ２）／（Ｌ１＋Ｌ２） ・・・（４）
この（４）式中では、量子井戸層とバリア層の格子定数をそれぞれａ１、ａ２で示しており、パラメータとしてａ１を選択すると、次のように表せる。
【００５７】
ａ１＝（ａ・Ｌ−ａ２・Ｌ２）／（Ｌ−Ｌ２） ・・・（５）
このように、周期・平均格子定数算出部４７は、複数の主要なサテライトピークの間隔から超格子構造の周期Ｌを求める周期算出手段の機能と、実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークＰ20とＩｎＰ基板に由来して実測プロファイル中に現れる基板ピークＰ10との間隔から超格子構造の平均格子定数ａを求める平均格子定数算出手段の機能と、を併有している。
【００５８】
ところで、上記（４）式の右辺の値が全て既知であるとき、すなわち、格子定数ａ２が定められると、従属関係により、一義的に他の格子定数ａ１が定まる。勿論、逆も真なりで、格子定数ａ１を定めると、自動的に格子定数ａ２が定まる。すなわち、従属関係によって、４つのパラメータが２つのパラメータに集約され得るということができる。
【００５９】
そこで、初期値設定／変更部４４（初期値設定手段）では、４つのパラメータ（変数）Ｌ１、Ｌ２、ａ１、ａ２のうち、２つのパラメータ、すなわち、第１組成薄膜又は第２組成薄膜のうち一方の薄膜の膜厚Ｌ１の初期値と、第１組成薄膜又は第２組成薄膜のうち一方の薄膜の格子定数ａ１の初期値とを、それぞれ後述するように設定するとともに、その設定値を変更するようになっている。なお、これら膜厚Ｌ１の初期値と格子定数ａ１の初期値がそれぞれ設定されると、上述の従属関係から、第１組成薄膜又は第２組成薄膜のうち他方の薄膜の膜厚Ｌ２の初期と、第１組成薄膜又は第２組成薄膜のうち他方の薄膜の格子定数ａ２は一義的に定まる。
【００６０】
また、初期値設定／変更部４４は、実測プロファイルから求めた超格子構造の周期Ｌ及び平均格子定数ａの値と、これらの設計値並びに前記一方及び他方の薄膜の膜厚Ｌ１、Ｌ２及び格子定数ａ１、ａ２の設計値といった設計値情報とに基づいて、周期Ｌ及び平均格子定数ａの設計値からのずれをも考慮して、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値を例えば設計値と同一の相互比率を保つようにそれぞれ設定する。さらに、初期値設定／変更部４４は、後述する条件判定部５２（判定手段）からの更新指令を受けたとき、設定した初期値を所定の微調整量を単位として順次変更（増加又は減少）させる機能を有している。
【００６１】
変数算出部４８（変数算出手段）は、初期値設定／変更部４４で設定された前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の各初期値又は変更された膜厚Ｌ１及び格子定数ａ１の各値と、周期・平均格子定数算出部４７で算出された超格子構造の周期Ｌ及び超格子構造の平均格子定数ａの値とに基づいて、上述の（３）、（５）式に示すような条件で、超格子を構成する他方の薄膜の膜厚Ｌ２及び格子定数ａ２をそれぞれ算出する。
【００６２】
そして、上述のようにして得られた４つのパラメータ、すなわち、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値、並びに、前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値は、それぞれ解析プロファイル作成部４９に取り込まれ、この解析プロファイル作成部４９で、４つパラメータＬ１、Ｌ２、ａ１、ａ２に基づいて、これらの値を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルが作成される。そして、このとき、プロファイル作成領域が特定領域に特定されていることから、４つのパラメータＬ１、Ｌ２、ａ１、ａ２（変数）に基づき、複数の主要なサテライトピークのうち次数の異なる特定のサテライトピークに対応する角度範囲を特定する前記ｎより少ない正の整数であるｍ個の角度値ｘｉ，ｘｊ，・・・ｘi+m-1（但し、ｉ、ｊはそれぞれｍより小さい正の整数）について、回折Ｘ線強度の値ｆｉ，ｆｊ，・・・ｆi+m-1を前記理論計算により算出して解析プロファイルを作成することになる。
【００６３】
比較部５１では、解析プロファイル及び実測プロファイルのデータを実測プロファイル作成部４２及び解析プロファイル作成部４９からそれぞれ取り込み、両プロファイルを比較して、両プロファイルの一致の度合いを調べる。ここでは、次数の異なる特定のサテライトピークにつき比較して、両プロファイルの一致の度合いを調べるが、この比較は、例えば比較対象のサテライトピーク領域について平均二乗誤差Ｅを求めることで行われる。この平均二乗誤差Ｅは。比較対象として選択された特定のサテライトピークの個々のピーク形成角度範囲において、各サンプリング点の実測プロファイルの値、例えばｆｉ（ｉ＝１、２、・・・ｎ）と解析プロファイルの値ｇｉの差（ｆｉ−ｇｉ）の二乗値を均等に加算してその平方根をとり、サンプル数ｎで除算したものである。すなわち、Ｅ＝｛（ｆ１−ｇ１）２＋（ｆ２−ｇ２）２＋・・・＋（ｆｎ−ｇｎ）２｝１／２／ｎを比較結果とする。
【００６４】
この比較部５１での比較結果は、条件判定部５２（判定手段）に取り込まれ、この条件判定部５２において、各比較対象のサテライトピークについての平均二乗誤差Ｅが予め設定された許容値以下に達しているか否かが判定される。すなわち、前記一致の度合いが予め定められた一致基準（許容値以下）を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値、並びに前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値を、それぞれの膜厚測定値及び格子定数の測定値とする。一方、両プロファイルの一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれか１つの設定値を所定の微調整量単位で初期値設定／変更部４４に変更させて、変数算出部４８に再度変数の算出を実行させる。また、その結果が再度一致基準を満たさないときは、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうち残りの１つの設定値を所定の微調整量単位で初期値設定／変更部４４に変更させて、変数算出部４８に再度変数の算出を実行させる。それ以降に、一致基準を満たさないときは、まず、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれか１つの設定値を所定の微調整量単位で前回の設定値から変更させ、次いで、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうち残りの１つの前回の設定値を所定の微調整量単位で再度変更させて、変数の再算出を実行させることになる。また、このような変数の変更による解析プロファイルの再作成（シミュレーション）を所定回数実行してもなお一致基準を満たさない場合は、当該サンプル基板１０を不良と判定する。
【００６５】
また、本実施形態においては、条件判定部５２は、比較部５１と協働して、複数のサテライトピークのうち超格子構造の平均格子定数ａに対応する０次のサテライトピーク（ピークの点のみ）について、解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ実測プロファイルに合わせるシミュレーションを最初に実行するようになっており、これらによって、比較手段のシミュレーション処理部が構成されている。
【００６６】
そのため、条件判定部５２が前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して変数算出部４８に再度算出を実行させるとき、初期値設定／変更部４４及び変数算出部４８は、超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整し、更に、このように超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整するのに先立って、前記一方の薄膜の格子定数ａ１を微調整するようになっている。
【００６７】
また、比較部５１で比較対象として選択されるサテライトピークは、初回は０次のサテライトピークとこれとは次数の異なるサテライトピーク（例えばマイナス側の一次サテライトピーク）であり、０次のサテライトピークについての平均二乗誤差Ｅが所定値以下となるように設定値が変更された後は、次数の異なるサテライトピークのみが比較対象であってもよい。すなわち、比較部５１は、０次のサテライトピークについて両プロファイルを一致させるシミュレーション処理の後、主として他の次数のサテライトピークについて、前記一致の度合いを調べる誤差解析部となっている。
【００６８】
このように、測定解析制御装置４０は、補助記憶装置４０５に格納されたＸ線回折測定解析プログラムに従って、解析プロファイルを特定する４つのパラメータを２つに集約させつつ変化させながら実測プロファイルと解析プロファイルを比較し、両プロファイルが最も良く一致した状態における解析プロファイルの４つの変数を、第１、第２組成薄膜の格子定数ａ１、ａ２及び膜厚Ｌ１、Ｌ２のそれぞれの測定値として出力する。
【００６９】
また、測定解析制御装置４０は、複数のサンプル基板１０について上述のような測定解析が繰り返される場合に、条件判定部５２により前記一致の度合いが予め定められた一致基準を満たすと判定されたとき、あるいは不良判定がなされたときにも、そのときの前記一方又は他方の薄膜の膜厚Ｌ１又はＬ２と前記一方又は他方の薄膜の格子定数ａ１又はａ２の値を、格子定数及び膜厚を縦軸及び横軸とする一致度分布図、例えば図５に示すような等高線図様の分布図あるいは図６に示しすように一致度に応じて無段階に輝度、又は色及び色濃度が変化する分布図上の点（図５中に丸印で示す）として補助記憶装置４０５（分布データ記憶手段）に記憶蓄積させていき、複数のサンプル基板１０についての測定・解析がなされたとき、その分布データ記憶情報に基づいて、解析プロファイル及び実測プロファイルの一致の度合いの分布画像情報を生成する分布図作成手段の機能と、この分布図作成手段の機能により生成された分布画像情報を表示部４０７又は外部に出力する出力手段の機能とをそれぞれ有している。
【００７０】
図５に一例で示すように、分布図作成手段としての測定解析制御装置４０は、補助記憶装置４０５内の分布データ記憶情報に基づいて、解析プロファイル及び実測プロファイルの一致の度合いが最も高い一致領域１０１と、前記一致の度合いが一致領域１０１よりも低く一致領域１０１を取り囲むように分布する少なくとも１つの環状分布領域、例えば４つの環状分布領域１０２、１０３、１０４とが視覚的に区別可能な表示態様で生成するようになっている。ここで、視覚的に区別可能な表示態様とは、例えば異なる色、濃さ、１つ以上の一致度等高線による表示領域の分割、図示するような塗りつぶし条件の相違（例えば異なる網掛け）であり、図６に示すように、色を一致度が高い内方側から低い外方側へと徐々に変化させ、結果的に内外での一致どの相違が判定できる無段階的な変化の態様であってもよい。なお、図６では多数の分布データの点を個々に視認できるようには図示せず、分布の度合いを示す同色（例えば黄色）同輝度の環状分布領域をグラデーションにより無段階に多数表示しているが、データ自体はこの分布図中の一座標を持つ点として記憶されることになる。また、図５に示すように、一致領域１０１（図６では最も輝度の高い白い楕円形部分）の中心が解析・実測プロファイルの完全な一致点であり、この点から離れるほどずれが大きいことになる。したがって、一致領域１０１を取り囲むように分布する４つの環状分布領域１０２、１０３、１０４（図６では白い楕円形部分を囲む外側ほど徐々に黒くなる環状部分）は、両プロファイルが僅かにずれている状態を示す内側の領域１０２から、両プロファイルが部分的には近いものの他の部分に隙間があるような状態を示す中間の領域１０３、両プロファイル間に明らかに隙間があるような状態を示す外側の領域１０４へと順次外側になり、その外方では、両プロファイルが乖離して大きな隙間があるような状態（図６中では濃い黒の部分）となる。
【００７１】
なお、本実施の形態においては、一方の薄膜である井戸層の膜厚Ｌ１が決まると、他方の薄膜であるバリア層の膜厚Ｌ２が一義的に定まるため、図５に示すように、一致度分布図の一方の軸、例えば横軸に２つの薄膜の膜厚を併記している。
【００７２】
次に、上述の測定解析装置を用いて実施するＸ線回折測定解析方法の一実施形態について説明する。
【００７３】
図７は、本実施形態の測定解析装置で実行されるＸ線回折測定解析プログラムの概略の流れを示すフローチャートであり、このプログラムにより本発明の方法を実施することができる。
【００７４】
まず、設計情報記憶部４３への設計値情報の記憶やピーク特定条件記憶部４６への条件入力等の準備を行った後（ステップＳ１０）、測定を実行し（ステップＳ１１）、実測データ記憶部４１に実測データを記憶情報として順次格納する（ステップＳ１２）。
【００７５】
次いで、実測データ記憶部４１の記憶情報に基づいて、実測プロファイル作成部４２で実測プロファイルを作成する（ステップＳ１３）。
【００７６】
次いで、特定ピーク選択部４５でピーク特定条件記憶部４６の記憶条件に従って複数のうち孤立した次数の異なる特定のサテライトピークの発生領域を特定し（ステップＳ１４）、周期・平均格子定数算出部４７で超格子構造の平均組成に対応するこれら複数のサテライトピーク、例えばサテライトピークＰ20、Ｐ21の間隔から、超格子構造の積層周期Ｌを求めるともに（ステップＳ１５）、実測プロファイルの中の主要なサテライトピーク、例えば０次のサテライトピークＰ20と基板ピークＰ10との間隔から超格子構造の平均格子定数ａを求める（ステップＳ１６）。
【００７７】
次いで、初期値設定／変更部４４で前記４つの変数のうち第１組成薄膜又は第２組成薄膜のうち一方の薄膜の膜厚Ｌ１の初期値と、第１組成薄膜又は第２組成薄膜のうち一方の薄膜の格子定数ａ１の初期値とを、設計情報記憶部４３の記憶情報及び周期・平均格子定数算出部４７での算出結果Ｌ、ａに基づいて、それぞれ設定する（ステップＳ１７；設定段階）。この初期値は、実測プロファイルから求めた超格子構造の積層周期Ｌと平均格子定数ａが設計値に近い場合には、設計値と同一の値であるが、実測プロファイルから求めた超格子構造の積層周期Ｌと平均格子定数ａが設計値からずれている場合には、前記微調整量を単位としてそのずれ側に多少加減される。
【００７８】
次いで、変数算出部４８で、超格子構造の周期Ｌと、超格子構造の平均格子定数ａと、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値とから、第２組成薄膜又は第１組成薄膜の膜厚のうち他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値をそれぞれ算出する（ステップＳ１８；変数算出段階）。
【００７９】
次いで、解析プロファイル作成部４９で、超格子構造の周期Ｌ、超格子構造の平均格子定数ａ、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値、並びに前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値に基づいて、前記理論計算により解析プロファイルを作成する（ステップＳ１９；解析プロファイル作成段階）。
【００８０】
次いで、比較部５１で、作成された解析プロファイル及び実測プロファイルを比較して、両プロファイルの一致の度合いを調べる（ステップＳ２０；比較段階）。
【００８１】
次いで、条件判定部５２で、０次のサテライトピークについての一致の度合いが予め定められた一致基準を満たすか否かを判定し（ステップＳ２１）、一致の度合いが予め定められた一致基準を満たすときは、次数の異なるサテライトピーク、例えば−１次サテライトピークについての一致の度合いが予め定められた一致基準を満たすか否かを判定する（ステップＳ２２）。勿論、必要により、２次サテライトピーク（例えば−２次サテライトピーク）についての一致度合いを判定するステップ（ステップＳ２３）があってもよい。
【００８２】
ここで一致の度合いが予め定められた一致基準を満たす場合、そのときの前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値並びに前記他方の薄膜の膜厚Ｌ２及び格子定数ａ２の値をそれぞれの膜厚測定値及び格子定数測定値として測定結果を表示部４０７に出力する（ステップＳ２４）。
【００８３】
一方、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれか１つの設定値を初期値設定／変更部４４で変更して（ステップＳ２６）、再び変数算出部４８での変数算出段階に戻る（ステップＳ１８）。
【００８４】
次いで、解析プロファイルを再度作成し（ステップＳ１９）、複数のサテライトピークのうち超格子構造の平均格子定数ａに対応する０次のサテライトピークＰ20について、比較部５１でのチェックを再度実行する（ステップＳ２１）。このように、最初に０次のサテライトピークＰ20のピーク点について、解析プロファイルを縦軸方向及び横軸方向でそれぞれ実測プロファイルに合わせるように、初期値設定／変更部４４での設定値の変更を実行させる。このとき、例えば、最初に０次ピーク強度を調整し、また、特定のサテライトピークＰ20、Ｐ21の角度間隔が解析プロファイルと実測プロファイルで一致していなければ、周期Ｌを加減してその間隔を近付ける。このような設定値変更は、図８中の左右両側に示すような解析プロファイルの変化をもたらす。すなわち、周期Ｌを増加させると、同図右側に示すようにサテライトピーク同士が近付き、周期Ｌを低減させると、同図左側に示すようにサテライトピーク同士が離れる。
【００８５】
そして、０次のサテライトピークＰ20のピーク点について一致基準を満たした場合に、主要サテライトピークのうち他の次数のサテライトピーク、例えば１次サテライトピークＰ21について、解析プロファイル及び実測プロファイルの一致の度合いを確認し（ステップＳ２２、Ｓ２３）、一致していなければ、初期値設定／変更部４４での設定値の変更を実行させる（ステップＳ２６）。
【００８６】
この間、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の設定値変更が所定回数に達したか否かを判別し（ステップＳ２５）、所定回数に達した場合には不良と判定し、その判定結果を表示部４０７に出力する（ステップＳ２７）。
【００８７】
また、条件判定部５２では、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して再び変数算出部４８での算出を再度実行させる際、超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整する（ステップＳ２６）。さらに、このように超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整するのに先立って、０次のサテライトピークＰ20のピーク点について一致基準を満たすまで、前記一方の薄膜の格子定数ａ１を微調整する（ステップＳ２６）。このような設定値変更は、図８中の上下に示すような解析プロファイルの変化をもたらす。すなわち、格子定数を増加させると、同図上側に示すように複数のサテライトピークに接する包絡線ＥＮが同図中の反時計方向に傾き、格子定数を低減させると、同図下側に示すように複数のサテライトピークに接する包絡線ＥＮが同図中の時計方向に傾くように、解析プロファイルが変化する。したがって、格子定数を加減すると、その影響は次数の多いサテライトピークほど強く現れることになる。
【００８８】
このような膜厚Ｌ１及び格子定数ａ１一連の設定変更により、解析プロファイルは、図９に実践で模式的に示すように、同図に破線で概形を示す実測プロファイルに対して、同図（ａ）に示す最初の初期状態段階から、同図（ｂ）に示す０次サテライトピークの一致段階を経て、同図（ｃ）に示すような一致基準を満たす段階へと収束的に変化することになる。
【００８９】
一方、測定結果を表示部４０７に出力するか（ステップＳ２４）、あるいは、不良判定結果を表示部４０７に出力するか（ステップＳ２７）した場合、条件判定部５２で、前記一方又は他方の薄膜の膜厚Ｌ１又はＬ２の値と、前記一方又は他方の薄膜の格子定数ａ１又はａ２の値とを、図５に例示したような格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として、補助記憶装置４０５に記憶しておき（ステップＳ２９）、この記憶情報の蓄積に基づいて解析プロファイル及び実測プロファイルの一致の度合いを示す分布図を後から作成できるようにしておく。
【００９０】
そして、分布図の作成を指示する操作入力がなされたか否かをチェックし（ステップＳ３０）、指示された場合には、補助記憶装置４０５に記憶蓄積された一致度のデータを基に分布図を作成して、その分布図上に、解析プロファイル及び実測プロファイルの一致の度合いが最も高い一致領域１０１と、一致の度合いが一致領域１０１よりも低く一致領域１０１を取り囲むように分布する少なくとも１つの環状分布領域１０２〜１０４とを上述したような視覚的に区別可能な表示態様で表示させ、一致度分布図上の点が、その分布図上の一致領域１０１及び環状分布領域１０２〜１０４のうちどの領域に属するかで前記一致の度合いを視覚的に識別可能にする（ステップＳ３２）。
【００９１】
また、測定終了を指示する操作入力がなされたか否かをチェックし、指示された場合には、処理を終了する（ステップＳ３３）。また、終了指示がなければ、ステップＳ１１に戻って、次の測定開始のため待機する。
【００９２】
このように、本実施形態の回折Ｘ線測定解析方法及びプログラムにおいては、実測プロファイルから超格子構造の周期Ｌと超格子構造の平均格子定数ａが求められ、これらの値と既知の情報とに基づいて、第１又は第２組成薄膜のうち一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値が４つの変数のうち２つの変数の初期値として設定され、これらの値を基に残りの２つの変数が算出されると、それまでに得られた情報に基づいて、理論計算により解析プロファイルが作成され、その解析プロファイルと実測プロファイルが比較され、その比較結果に応じて、一致の度合いが一致基準を満たすときには、そのときの４つの変数が測定値とされ、基準を満たさないときには、２つの変数である前記初期値のいずれかを変更して再び変数算出段階に戻る。したがって、適切な２つの初期値設定とこれら設定値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行することができ、しかも、所要の測定精度が確保できることになる。
【００９３】
また、本実施形態においては、設定段階において、積層された超格子構造に関する設計値情報に基づき、実測プロファイルから周期性のある孤立したピークを抽出して主要サテライトピークを特定するので、小ピーク、ノイズ等の微視的プロファイル部分の悪影響をより有効に回避でき、しかも、処理を簡単化することができる。
【００９４】
しかも、前記設定段階において、実測プロファイルから求めた超格子構造の周期Ｌ及び平均格子定数ａと、前記一方及び他方の薄膜の膜厚Ｌ１、Ｌ２及び格子定数ａ１、ａ２の設計値情報とに基づいて、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１の初期値を設定するので、実測プロファイルや設計情報等から得られる情報を基に、より得るべき測定値に近い初期値が設定できる。
【００９５】
また、本実施形態においては、比較段階において、複数のサテライトピークのうち超格子構造の平均格子定数ａに対応する０次のサテライトピークＰ20について、解析プロファイルを縦軸方向及び横軸方向でそれぞれ実測プロファイルに合わせ、主要サテライトピークのうち他の次数のサテライトピークＰ21等について、解析プロファイル及び実測プロファイルの一致の度合いを調べるので、４つの変数の変化に対する解析プロファイルの変化の特性を考慮して、一致度の比較処理を簡単化することができる。
【００９６】
また、判定段階において、一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して再び変数算出段階に戻る際、超格子構造の周期Ｌを一定にしつつ前記一方の薄膜の膜厚Ｌ１の設定値を微調整するので、実測プロファイルから正確に得られる周期Ｌを固定することで、一致点への収束を早めることができる。
【００９７】
さらに、前記一方の薄膜の膜厚Ｌ１及び格子定数ａ１のうちいずれかの設定値を変更して再び変数算出段階に戻る際、超格子構造の周期Ｌを一定にしつつ膜厚Ｌ１の設定値を微調整するのに先立って、前記一方の薄膜の格子定数ａ１を微調整するので、膜厚Ｌ１と格子定数ａ１の変化における相関性を考慮して、一致点への収束を早めることができる。
【００９８】
また、本実施形態においては、前記判定段階で、一方又は他方の薄膜の膜厚値Ｌ１又はＬ２と、一方又は他方の薄膜の格子定数値ａ１又はａ２とを、格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶しておき、その記憶情報に基づいて解析プロファイル及び実測プロファイルの一致の度合いを示す分布図を作成するとともに、その分布図上に、解析プロファイル及び実測プロファイルの一致の度合いが最も高い一致領域１０１と、一致の度合いが一致領域１０１よりも低く一致領域１０１を取り囲むように分布する環状分布領域１０２〜１０４とを視覚的に区別される表示態様で表示し、分布図上の点が、その分布図上の一致領域１０１及び環状分布領域１０２〜１０４のうちどの領域に属するかで一致の度合いを視覚的に識別可能に表示するので、解析プロファイル及び実測プロファイルの一致の度合いがどの程度かを容易に見て取れることになり、ユーザーの迅速・的確な大局的判断が可能になる。
【００９９】
なお、上述の実施形態においては、異種の薄膜を交互に積層した超格子構造を有する結晶構造として、通信分野で使用される高出力の半導体レーザの活性層を例としていたが、半導体レーザ以外の異種の薄膜を交互に積層した他の超格子構造にも適用できることはいうまでもない。
【産業上の利用可能性】
【０１００】
本発明は、以上説明したように、サテライトピーク間にある小さなピークやノイズはフィッティングに用いることなく、既知の情報及び実測データに基づいて選択した特定のサテライトピークについて、適切な２つの変数の設定値設定とこれら初期値のいずれか一方又は双方を変更するだけのフィッティング処理を収束計算的に迅速に実行するので、超格子構造を反映する主要サテライトピークについて迅速で且つ有効なフィッティングを行ない、必要にして十分な精度の解析が行なえるＸ線回折測定解析方法及びプログラムを実現することができ、さらに、一致度の分布表示による大局的判断を可能にできるＸ線回折測定解析方法及びプログラムを実現できるという効果を奏するものであり、異種の薄膜を交互に積層した超格子構造の測定・解析に用いるＸ線回折測定解析方法及びプログラム全般に有用である。
【図面の簡単な説明】
【０１０１】
【図１】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法を実施する装置の測定系の概略構成図である。
【図２】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法を実施する装置の全体の概略構成図である。
【図３】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法を実施する装置の制御装置の要部機能ブロック図である。
【図４】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で作成される解析・実測プロファイルの一例を例示する回折Ｘ線のプロファイル図で、縦軸は回折Ｘ線強度を、横軸は回折角度をそれぞれ示している。
【図５】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で作成される解析・実測プロファイルの一致度分布図で、縦軸は一方又は他方の薄膜の格子定数を、横軸はその薄膜の膜厚としている。
【図６】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で作成される解析・実測プロファイルの異なる態様の一致度分布図で、縦軸は一方又は他方の薄膜の格子定数を、横軸はその薄膜の膜厚としている。
【図７】本発明の第１の実施の形態に係るＸ線回折測定解析プログラムの概略の流れを示すフローチャートである。
【図８】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で行われるパラメータ変更とそれによる解析プロファイルの変化の関係を説明する説明図である。
【図９】本発明の第１の実施の形態に係る回折Ｘ線測定解析方法で行われるフィッティング処理とそれにより解析プロファイルが実測プロファイルに近付く状態を説明する説明図である。
【符号の説明】
【０１０２】
１ 第１回転台
２ 第２回転台
３、４ 駆動手段
６、７ 回転位置検出器
１０ サンプル基板（基板）
２０ Ｘ線源
２１ Ｘ線管
２２ 湾曲ミラー
２３ モノクロメータ
３０ Ｘ線カウンタ
４０ 測定解析制御装置（分布図作成手段、出力手段）
４１ 実測データ記憶部（実測データ記憶手段）
４２ 実測プロファイル作成部（実測プロファイル作成手段）
４３ 設計情報記憶部（設計情報記憶手段）
４５ 特定ピーク選択部
４６ ピーク特定条件記憶部（特定条件記憶手段）
４７ 周期・平均格子定数算出部（周期算出手段、平均格子定数算出手段）
４８ 変数算出部（変数算出手段）
４９ 解析プロファイル作成部（解析プロファイル作成手段）
５１ 比較部（シミュレーション処理部、誤差解析部）
５２ 条件判定部（判定手段、シミュレーション処理部）
１０１ 一致領域
１０２、１０３、１０４ 環状分布領域
４０１ ＣＰＵ
４０２ ＲＯＭ
４０３ ＲＡＭ
４０４ 入出力インターフェース回路
４０５ 補助記憶装置（分布データ記憶手段）
４０６ 設定・入力部
４０７ 表示部（出力手段）
【特許請求の範囲】
【請求項１】
基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、
４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、
前記４つの変数を変化させながら比較し、
最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析方法において、
前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求め、且つ前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求めるとともに、前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１）の初期値を設定する設定段階と、
前記超格子構造の周期と前記超格子構造の平均格子定数と前記一方の薄膜の膜厚及び格子定数の設定値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を算出する変数算出段階と、
前記超格子構造の周期、前記超格子構造の平均格子定数、前記一方の薄膜の膜厚及び格子定数の設定値並びに前記他方の薄膜の膜厚及び格子定数の算出値に基づいて、前記理論計算により前記解析プロファイルを得る解析プロファイル作成段階と、
前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較段階と、
前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚及び格子定数の初期値並びに前記他方の薄膜の膜厚及び格子定数の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して再び変数算出段階に戻る判定段階と、を含むことを特徴とするＸ線回折測定解析方法。
【請求項２】
前記設定段階において、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイルから周期性のある孤立したピークを抽出して前記主要サテライトピークを特定することを特徴とする請求項１に記載のＸ線回折測定解析方法。
【請求項３】
前記設定段階において、前記実測プロファイルから求めた前記超格子構造の周期及び平均格子定数と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚及び格子定数の初期値を設定することを特徴とする請求項１に記載のＸ線回折測定解析方法。
【請求項４】
前記比較段階において、前記複数のサテライトピークのうち前記超格子構造の平均格子定数に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせ、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べることを特徴とする請求項１乃至３のいずれか１つに記載のＸ線回折測定解析方法。
【請求項５】
前記判定段階において、前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整することを特徴とする請求項１乃至４のいずれか１つに記載のＸ線回折測定解析方法。
【請求項６】
前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整するのに先立って、前記一方の薄膜の格子定数を微調整することを特徴とする請求項５に記載のＸ線回折測定解析方法。
【請求項７】
前記判定段階で、前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶しておき、該記憶情報に基づいて前記解析プロファイル及び前記実測プロファイルの一致の度合いを示す分布図を作成するとともに、
該分布図上に、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とを視覚的に区別される表示態様で表示し、
前記分布図上の点が、該分布図上の前記一致領域及び環状分布領域のうちどの領域に属するかで前記一致の度合いを視覚的に識別可能に表示することを特徴とする請求項１乃至６のいずれか１つに記載のＸ線回折測定解析方法。
【請求項８】
基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、
４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、
前記４つの変数を変化させながら比較し、
最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、
コンピュータにより、
前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期を求める周期算出手段（４７）と、
前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数を求める平均格子定数算出手段（４７）と、
前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数の初期値を設定する初期値設定手段（４４）と、
前記超格子構造の周期と前記超格子構造の平均格子定数と前記一方の薄膜の膜厚及び格子定数の設定値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚及び格子定数の値をそれぞれ算出する変数算出手段（４８）と、
前記超格子構造の周期、前記超格子構造の平均格子定数、前記一方の薄膜の膜厚及び格子定数の設定値並びに前記他方の薄膜の膜厚及び格子定数の算出値に基づいて、前記理論計算により前記解析プロファイルを作成する解析プロファイル作成手段（４９）と、
前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較手段（５１）と、
前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚（Ｌ１）及び格子定数の設定値並びに前記他方の薄膜の膜厚及び格子定数の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段（５２）と、の機能を実現させることを特徴とするＸ線回折測定解析プログラム。
【請求項９】
前記初期値設定手段が、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイル中で前記超格子構造の周期に対応する周期性を有する孤立した複数のピークの発生領域を特定する条件を記憶する特定条件記憶部を有し、該特定条件記憶部の記憶情報に基づいて前記主要サテライトピークを特定することを特徴とする請求項８に記載のＸ線回折測定解析プログラム。
【請求項１０】
前記初期値設定手段が、前記実測プロファイルから求めた前記超格子構造の周期及び平均格子定数と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚及び格子定数の初期値を設定することを特徴とする請求項８又は９に記載のＸ線回折測定解析プログラム。
【請求項１１】
前記比較手段が、前記複数のサテライトピークのうち前記超格子構造の平均格子定数に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせるシミュレーションを実行するシミュレーション処理部（５１、５２）と、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べる誤差解析部（５１）とを有することを特徴とする請求項８乃至１０のいずれか１つに記載のＸ線回折測定解析プログラム。
【請求項１２】
前記判定手段が前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整することを特徴とする請求項８乃至１１のいずれか１つに記載のＸ線回折測定解析プログラム。
【請求項１３】
前記判定手段が前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整するとともに、前記一方の薄膜の格子定数を微調整することを特徴とする請求項１２に記載のＸ線回折測定解析プログラム。
【請求項１４】
前記コンピュータに、
前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶する分布データ記憶手段（４０５）と、
該分布データ記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いの分布画像情報を生成する分布図作成手段（４０）と、
該分布図作成手段で生成された分布画像情報を出力する出力手段と、の機能を実現させることを特徴とする請求項８乃至１３のいずれか１つに記載のＸ線回折測定解析プログラム。
【請求項１５】
前記分布図作成手段が、前記分布データ記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とが視覚的に区別可能な表示画像の情報を生成することを特徴とする請求項１４に記載のＸ線回折測定解析プログラム。
【請求項１６】
基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を異なる角度で入射したとき測定情報に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の実測プロファイルと、
４つの変数として前記第１組成薄膜の格子定数（ａ１）及び膜厚（Ｌ１）並びに前記第２組成薄膜の格子定数（ａ２）及び膜厚（Ｌ２）を用いた理論計算に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の解析プロファイルとを、
前記解析プロファイルの４つの変数を変化させながら比較し、
最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、
ｎを正の整数とするとき、ｎ個の前記角度の値（ｘ１，ｘ２，・・・ｘｎ）および該角度に対応するｎ個の回折Ｘ線強度の測定値（ｇ１，ｇ２，・・・ｇｎ）を順次メモリに記憶させる実測データ記憶手段（４１）と、
前記メモリに記憶されたデータに基づいて前記角度に応じて変化する回折Ｘ線強度の実測プロファイルを表示出力可能な所定のプロファイル表示形式で作成する実測プロファイル作成手段（４２）と、
前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求める周期算出手段（４７）と、
前記実測プロファイルの中に超格子構造の平均組成に対応して現れるサテライトピークと前記実測プロファイルの中に前記基板に由来して現れる基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求める平均格子定数算出手段（４７）と、
前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１又はＬ２）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１又はａ２）の初期値を、前記４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）のうち２つの変数の初期値として設定する初期値設定手段（４４）と、
前記超格子構造の周期（Ｌ）と、前記超格子構造の平均格子定数（ａ）と、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値とから、前記４つの変数（ａ１、Ｌ１、ａ２、Ｌ２）のうち残りの２つの変数である前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値をそれぞれ算出する変数算出手段（４８）と、
前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値と、前記変数算出手段（４８）で算出された前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値とを、前記４つの変数として、前記複数のサテライトピークのうち次数の異なる特定のサテライトピークに対応する角度範囲を特定する前記ｎより少ない正の整数であるｍ個の角度値（ｘｉ，ｘｊ，・・・ｘi+m-1；但し、ｉ、ｊはそれぞれｍより小さい正の整数）について、回折Ｘ線強度の値（ｆｉ，ｆｊ，・・・ｆｉ＋ｍ−１）を前記理論計算により算出して前記解析プロファイルを作成する解析プロファイル作成手段（４９）と、
前記解析プロファイル及び前記実測プロファイルを前記複数のサテライトピークのうち前記次数の異なる特定のサテライトピークにつき比較して、両プロファイルの一致の度合いを調べる比較手段（５１）と、
前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記解析プロファイルの４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）を前記一方の薄膜の膜厚及び格子定数並びに前記他方の薄膜の膜厚及び格子定数の測定値として出力し、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段（５２）と、の機能をコンピュータを使用して実現させることを特徴とするＸ線回折測定解析プログラム。
【請求項１】
基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、
４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、
前記４つの変数を変化させながら比較し、
最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析方法において、
前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求め、且つ前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求めるとともに、前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１）の初期値を設定する設定段階と、
前記超格子構造の周期と前記超格子構造の平均格子定数と前記一方の薄膜の膜厚及び格子定数の設定値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値を算出する変数算出段階と、
前記超格子構造の周期、前記超格子構造の平均格子定数、前記一方の薄膜の膜厚及び格子定数の設定値並びに前記他方の薄膜の膜厚及び格子定数の算出値に基づいて、前記理論計算により前記解析プロファイルを得る解析プロファイル作成段階と、
前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較段階と、
前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚及び格子定数の初期値並びに前記他方の薄膜の膜厚及び格子定数の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して再び変数算出段階に戻る判定段階と、を含むことを特徴とするＸ線回折測定解析方法。
【請求項２】
前記設定段階において、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイルから周期性のある孤立したピークを抽出して前記主要サテライトピークを特定することを特徴とする請求項１に記載のＸ線回折測定解析方法。
【請求項３】
前記設定段階において、前記実測プロファイルから求めた前記超格子構造の周期及び平均格子定数と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚及び格子定数の初期値を設定することを特徴とする請求項１に記載のＸ線回折測定解析方法。
【請求項４】
前記比較段階において、前記複数のサテライトピークのうち前記超格子構造の平均格子定数に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせ、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べることを特徴とする請求項１乃至３のいずれか１つに記載のＸ線回折測定解析方法。
【請求項５】
前記判定段階において、前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整することを特徴とする請求項１乃至４のいずれか１つに記載のＸ線回折測定解析方法。
【請求項６】
前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して再び変数算出段階に戻る際、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整するのに先立って、前記一方の薄膜の格子定数を微調整することを特徴とする請求項５に記載のＸ線回折測定解析方法。
【請求項７】
前記判定段階で、前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶しておき、該記憶情報に基づいて前記解析プロファイル及び前記実測プロファイルの一致の度合いを示す分布図を作成するとともに、
該分布図上に、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とを視覚的に区別される表示態様で表示し、
前記分布図上の点が、該分布図上の前記一致領域及び環状分布領域のうちどの領域に属するかで前記一致の度合いを視覚的に識別可能に表示することを特徴とする請求項１乃至６のいずれか１つに記載のＸ線回折測定解析方法。
【請求項８】
基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を入射し回折されたＸ線を検出して得られるＸ線回折測定情報に基づき、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされ、前記基板よりの回折Ｘ線に由来する基板ピーク及び前記積層された超格子構造に由来する複数のサテライトピークを有する実測プロファイルと、
４つの変数として前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚を用いたＸ線回折現象の理論計算により、横軸をＸ線回折角度、縦軸をＸ線回折強度として表わされる解析プロファイルとを、
前記４つの変数を変化させながら比較し、
最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、
コンピュータにより、
前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期を求める周期算出手段（４７）と、
前記実測プロファイルの中の超格子構造の平均組成に対応するサテライトピークと前記基板ピークとの間隔から前記超格子構造の平均格子定数を求める平均格子定数算出手段（４７）と、
前記４つの変数のうち前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数の初期値を設定する初期値設定手段（４４）と、
前記超格子構造の周期と前記超格子構造の平均格子定数と前記一方の薄膜の膜厚及び格子定数の設定値とから、前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚及び格子定数の値をそれぞれ算出する変数算出手段（４８）と、
前記超格子構造の周期、前記超格子構造の平均格子定数、前記一方の薄膜の膜厚及び格子定数の設定値並びに前記他方の薄膜の膜厚及び格子定数の算出値に基づいて、前記理論計算により前記解析プロファイルを作成する解析プロファイル作成手段（４９）と、
前記解析プロファイル及び前記実測プロファイルを比較して、両プロファイルの一致の度合いを調べる比較手段（５１）と、
前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記一方の薄膜の膜厚（Ｌ１）及び格子定数の設定値並びに前記他方の薄膜の膜厚及び格子定数の値を膜厚測定値及び格子定数測定値とし、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段（５２）と、の機能を実現させることを特徴とするＸ線回折測定解析プログラム。
【請求項９】
前記初期値設定手段が、前記積層された超格子構造に関する設計値情報に基づき、前記実測プロファイル中で前記超格子構造の周期に対応する周期性を有する孤立した複数のピークの発生領域を特定する条件を記憶する特定条件記憶部を有し、該特定条件記憶部の記憶情報に基づいて前記主要サテライトピークを特定することを特徴とする請求項８に記載のＸ線回折測定解析プログラム。
【請求項１０】
前記初期値設定手段が、前記実測プロファイルから求めた前記超格子構造の周期及び平均格子定数と、前記一方及び他方の薄膜の膜厚（Ｌ１、Ｌ２）及び格子定数（ａ１、ａ２）の設計値情報とに基づいて、前記一方の薄膜の膜厚及び格子定数の初期値を設定することを特徴とする請求項８又は９に記載のＸ線回折測定解析プログラム。
【請求項１１】
前記比較手段が、前記複数のサテライトピークのうち前記超格子構造の平均格子定数に対応する０次のサテライトピークについて、前記解析プロファイルを前記縦軸方向及び横軸方向でそれぞれ前記実測プロファイルに合わせるシミュレーションを実行するシミュレーション処理部（５１、５２）と、前記主要サテライトピークのうち他の次数のサテライトピークについて、前記解析プロファイル及び前記実測プロファイルの一致の度合いを調べる誤差解析部（５１）とを有することを特徴とする請求項８乃至１０のいずれか１つに記載のＸ線回折測定解析プログラム。
【請求項１２】
前記判定手段が前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整することを特徴とする請求項８乃至１１のいずれか１つに記載のＸ線回折測定解析プログラム。
【請求項１３】
前記判定手段が前記一方の薄膜の膜厚及び格子定数のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させるとき、前記初期値設定手段が、前記超格子構造の周期を一定にしつつ前記一方の薄膜の膜厚の設定値を微調整するとともに、前記一方の薄膜の格子定数を微調整することを特徴とする請求項１２に記載のＸ線回折測定解析プログラム。
【請求項１４】
前記コンピュータに、
前記一方又は他方の薄膜の膜厚（Ｌ１又はＬ２）と前記一方又は他方の薄膜の格子定数（ａ１又はａ２）の値を、該格子定数及び膜厚を縦軸及び横軸とする等高線図様の分布図上の点として記憶する分布データ記憶手段（４０５）と、
該分布データ記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いの分布画像情報を生成する分布図作成手段（４０）と、
該分布図作成手段で生成された分布画像情報を出力する出力手段と、の機能を実現させることを特徴とする請求項８乃至１３のいずれか１つに記載のＸ線回折測定解析プログラム。
【請求項１５】
前記分布図作成手段が、前記分布データ記憶手段の記憶情報に基づいて、前記解析プロファイル及び前記実測プロファイルの一致の度合いが最も高い一致領域と、前記一致の度合いが前記一致領域よりも低く前記一致領域を取り囲むように分布する少なくとも１つの環状分布領域とが視覚的に区別可能な表示画像の情報を生成することを特徴とする請求項１４に記載のＸ線回折測定解析プログラム。
【請求項１６】
基板上に複数の第１組成薄膜と複数の第２組成薄膜とが交互に積層された超格子構造を有する結晶構造にＸ線を異なる角度で入射したとき測定情報に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の実測プロファイルと、
４つの変数として前記第１組成薄膜の格子定数（ａ１）及び膜厚（Ｌ１）並びに前記第２組成薄膜の格子定数（ａ２）及び膜厚（Ｌ２）を用いた理論計算に基づきプロファイル表示される、前記角度に応じて変化する回折Ｘ線強度の解析プロファイルとを、
前記解析プロファイルの４つの変数を変化させながら比較し、
最も良く一致した状態における前記解析プロファイルの４つの変数を、前記第１組成薄膜の格子定数及び膜厚並びに前記第２組成薄膜の格子定数及び膜厚のそれぞれの測定値とするＸ線回折測定解析プログラムであって、
ｎを正の整数とするとき、ｎ個の前記角度の値（ｘ１，ｘ２，・・・ｘｎ）および該角度に対応するｎ個の回折Ｘ線強度の測定値（ｇ１，ｇ２，・・・ｇｎ）を順次メモリに記憶させる実測データ記憶手段（４１）と、
前記メモリに記憶されたデータに基づいて前記角度に応じて変化する回折Ｘ線強度の実測プロファイルを表示出力可能な所定のプロファイル表示形式で作成する実測プロファイル作成手段（４２）と、
前記実測プロファイルの中の複数の主要なサテライトピークの間隔から前記超格子構造の周期（Ｌ）を求める周期算出手段（４７）と、
前記実測プロファイルの中に超格子構造の平均組成に対応して現れるサテライトピークと前記実測プロファイルの中に前記基板に由来して現れる基板ピークとの間隔から前記超格子構造の平均格子定数（ａ）を求める平均格子定数算出手段（４７）と、
前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の膜厚（Ｌ１又はＬ２）及び前記第１組成薄膜又は前記第２組成薄膜のうち一方の薄膜の格子定数（ａ１又はａ２）の初期値を、前記４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）のうち２つの変数の初期値として設定する初期値設定手段（４４）と、
前記超格子構造の周期（Ｌ）と、前記超格子構造の平均格子定数（ａ）と、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値とから、前記４つの変数（ａ１、Ｌ１、ａ２、Ｌ２）のうち残りの２つの変数である前記第１組成薄膜又は前記第２組成薄膜の膜厚のうち他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値をそれぞれ算出する変数算出手段（４８）と、
前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）の設定値と、前記変数算出手段（４８）で算出された前記他方の薄膜の膜厚（Ｌ２）及び格子定数（ａ２）の値とを、前記４つの変数として、前記複数のサテライトピークのうち次数の異なる特定のサテライトピークに対応する角度範囲を特定する前記ｎより少ない正の整数であるｍ個の角度値（ｘｉ，ｘｊ，・・・ｘi+m-1；但し、ｉ、ｊはそれぞれｍより小さい正の整数）について、回折Ｘ線強度の値（ｆｉ，ｆｊ，・・・ｆｉ＋ｍ−１）を前記理論計算により算出して前記解析プロファイルを作成する解析プロファイル作成手段（４９）と、
前記解析プロファイル及び前記実測プロファイルを前記複数のサテライトピークのうち前記次数の異なる特定のサテライトピークにつき比較して、両プロファイルの一致の度合いを調べる比較手段（５１）と、
前記一致の度合いが予め定められた一致基準を満たすか否かを判定し、前記一致の度合いが予め定められた一致基準を満たすときは、そのときの前記解析プロファイルの４つの変数（ａ１，Ｌ１，ａ２，Ｌ２）を前記一方の薄膜の膜厚及び格子定数並びに前記他方の薄膜の膜厚及び格子定数の測定値として出力し、前記一致の度合いが予め定められた一致基準を満たさないときは、前記一方の薄膜の膜厚（Ｌ１）及び格子定数（ａ１）のうちいずれかの設定値を変更して前記変数算出手段に再度算出を実行させる判定手段（５２）と、の機能をコンピュータを使用して実現させることを特徴とするＸ線回折測定解析プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図７】
【図８】
【図９】
【図６】
【図２】
【図３】
【図４】
【図５】
【図７】
【図８】
【図９】
【図６】
【公開番号】特開２００６−１７０７９１（Ｐ２００６−１７０７９１Ａ）
【公開日】平成１８年６月２９日（２００６．６．２９）
【国際特許分類】
【出願番号】特願２００４−３６３４４３（Ｐ２００４−３６３４４３）
【出願日】平成１６年１２月１５日（２００４．１２．１５）
【出願人】（００００００５７２）アンリツ株式会社 (838)
【Ｆターム（参考）】
【公開日】平成１８年６月２９日（２００６．６．２９）
【国際特許分類】
【出願日】平成１６年１２月１５日（２００４．１２．１５）
【出願人】（００００００５７２）アンリツ株式会社 (838)
【Ｆターム（参考）】
[ Back to top ]