(57)【要約】　　　（修正有）【課題】成膜等の反応機構解析を短期間、低コストで実現する、反応機構解析プログラム及び当該記録媒体並びに反応機構解析装置を提供する。【解決手段】反応器における反応機構を解析するためにコンピュータを、実験により得られた成膜形状の指針値を算出する実験値算出手段１、成膜形状を生成する成膜形状シミュレーション手段２、各成膜形状の指針値を算出するシミュレーション値算出手段３、各成膜形状の実験条件値および指針値に基づいて、実験条件と成膜形状指針との関係を表す実験条件−成膜形状指針モデルを生成するモデル化手段４、および　実験条件−成膜形状指針モデルを用いて算出される成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算する進化手段５、として機能させる。【選択図】　図１
2004119742.tif 000002
【特許請求の範囲】
【請求項１】
１種または複数種の反応器における化学反応の反応機構を解析するためにコンピュータを、
実験により得られた成膜形状の実験成膜形状指針値を算出する実験値算出手段、
数値シミュレーションによって、実験条件値を変化させながら複数のシミュレート成膜形状を生成する成膜形状シミュレーション手段、
各シミュレート成膜形状のシミュレート成膜形状指針値を算出するシミュレーション値算出手段、
各シミュレート成膜形状の実験条件値およびシミュレート成膜形状指針値に基づいて、実験条件と成膜形状指針との関係を表す実験条件−成膜形状指針モデルを生成するモデル化手段、および
進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算する進化手段、
として機能させることを特徴とする反応機構解析プログラム。
【請求項２】
１種または複数種の反応器における化学反応の反応機構を解析するためにコンピュータを、
実験により得られた成膜形状の実験成膜形状指針値を算出する実験値算出手段、
数値シミュレーションによって、実験条件値を変化させながら複数のシミュレート成膜形状を生成する成膜形状シミュレーション手段、
各シミュレート成膜形状のシミュレート成膜形状指針値を算出するシミュレーション値算出手段、
各シミュレート成膜形状の実験条件値およびシミュレート成膜形状指針値に基づいて、実験条件と成膜形状指針との関係を表す実験条件−成膜形状指針モデルを生成するモデル化手段、および
進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算する進化手段、
として機能させることを特徴とする反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項３】
１種または複数種の反応器における化学反応の反応機構を解析する装置であって、
実験により得られた成膜形状の実験成膜形状指針値を算出する実験値算出手段、
数値シミュレーションによって、実験条件値を変化させながら複数のシミュレート成膜形状を生成する成膜形状シミュレーション手段、
各シミュレート成膜形状のシミュレート成膜形状指針値を算出するシミュレーション値算出手段、
各シミュレート成膜形状の実験条件値およびシミュレート成膜形状指針値に基づいて、実験条件と成膜形状指針との関係を表す実験条件−成膜形状指針モデルを生成するモデル化手段、および
進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算する進化手段、
を備えることを特徴とする反応機構解析装置。
【請求項４】
進化手段は、進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値および前記実験時の実験条件から直接数値シミュレーションによって算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算することを特徴とする請求項１記載の反応機構解析プログラムまたは請求項２記載の反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体または請求項３記載の反応機構解析装置。
【請求項５】
反応機構が、
反応物が中間体を経て生成物にいたる気相および表面における定性的な反応経路、
反応物および中間体で構成される成膜種の付着確率、
成膜種の表面反応速度定数、
成膜種の気相反応速度定数、および
成膜種の濃度、
のうちの少なくともいずれかである請求項１記載の反応機構解析プログラムまたは請求項２記載の反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体または請求項３記載の反応機構解析装置。
【請求項６】
成膜形状指針が、基板上の任意の複数位置における膜厚の相互関係を表す関数である請求項１記載の成膜反応機構解析プログラムまたは請求項２記載の成膜反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体または請求項３記載の成膜反応機構解析装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この出願の発明は、反応機構解析プログラムおよび当該プログラムを記録したコンピュータ読み取り可能な記録媒体、ならびに反応機構解析装置に関するものである。さらに詳しくは、この出願の発明は、半導体デバイス製造等に有用な、成膜等の化学反応に関わる反応機構の解析を短期間で、且つ低コストで実現することのできる、新しい反応機構解析プログラムおよび当該プログラムを記録したコンピュータ読み取り可能な記録媒体、ならびに反応機構解析装置に関するものである。
【０００２】
【従来の技術】
従来、成膜等の化学反応に関わる反応機構の解析においては、たとえば、ＶＩＡホールやトレンチ等のテストストラクチャを持つ基板上にＣＶＤ（化学気相成長）などで成膜を行って成膜形状の実験値を得る一方で、モンテカルロシミュレーションなどで実験条件値を変化させながら様々な成膜形状のシミュレーション値を得て、それらを目視あるいは何らかの画像処理技術によって比較し、試行錯誤的に成膜実験を行いながら、妥当な反応機構を推測することがしばしば行われている（たとえば、特許文献１，２および非特許文献１，２，３参照）。
【０００３】
たとえば図１２は、トレンチへの成膜過程における１成膜種の付着確率ηを解析する場合の一従来例を示しており、トレンチの実際の成膜形状をＳＥＭ（走査型電子顕微鏡）で撮影する一方、トレンチ形状データおよび実験条件データ等を用いてモンテカルロシミュレーションにより成膜条件パラメータとしての付着確率ηを様々に変化させて複数のシミュレート成膜形状を算出し、そのシミュレート成膜形状と実際の成膜形状のＳＥＭ画像とを比較して正しい付着確率ηを試行錯誤的に決定している。
【０００４】
【特許文献１】
特開２００１−２７４１５２号公報
【特許文献２】
特開平９−２４６１８９号公報
【非特許文献１】
小宮山宏著，「ＣＲＥＡＴＩＶＥ　ＣＨＥＭＩＣＡＬ　ＥＮＧＩＮＥＥＲＩＮＧ　ＣＯＵＲＳＥ　３反応工学
反応装置から地球まで」，初版，培風館，１９９５年６月２６日，ｐ．１３０−１３１
【非特許文献２】
Ｙ．Ｓｈｉｍｏｇａｋｉ　ｅｔ　ａｌ．，　”Ａ　Ｋｉｎｅｔｉｃ　Ｓｔｕｄｙ　ｏｆ　ＷＳｉｘ−ＣＶＤ　Ｐｒｏｃｅｓｓｅｓ　ＢｙＭｏｎｔｅ　Ｃａｒｌｏ　Ｓｉｍｕｌａｔｉｏｎ”，　Ｃｏｎｆｅｒｅｎｃｅ　Ｐｒｏｃｅｅｄｉｎｇｓ　ＵＬＳＩ−ＶＩＩＩ，Ｕ．Ｓ．Ａ．，　１９９３，　ｐｐ．８８−９４
【非特許文献３】
Ｋ．Ｆｕｊｉｎｏ　ｅｔ　ａｌ．，　”Ｓｔｅｐ−Ｃｏｖｅｒａｇｅ　Ｓｉｍｕｌａｔｉｏｎ　ｆｏｒ　Ｔｅｔｒａｅｔｈｏｘｙｓｉｌａｎｅａｎｄ　Ｏｚｏｎｅ　Ａｔｍｏｓｐｈｅｒｉｃ　Ｐｒｅｓｓｕｒｅ　Ｃｈｅｍｉｃａｌ　Ｖａｐｏｒ　Ｄｅｐｏｓｉｔｉｏｎ”，Ｊｏｕｒｎａｌ　ｏｆ　ｔｈｅ　　Ｅｌｅｃｔｒｏｃｈｅｍｉｃａｌ　Ｓｏｃｉｅｔｙ，　Ｕ．Ｓ．Ａ，　Ａｕｇｕｓｔ　１９９３，Ｖｏｌ．１４０，　Ｎｏ．８，　ｐｐ．１２５−１２８
【０００５】
【発明が解決しようとする課題】
しかしながら、たとえば半導体デバイス製造においては、サブマイクロメータサイズのＶＩＡホールやトレンチ構造を薄膜によって”す”を残さずに埋め込む、つまり均一に成膜することが、断線や絶縁破壊防止、デバイスの長期的な信頼性確保などの観点から必要とされているところ、デバイスの高集積化によりＶＩＡホールやトレンチの高アスペクト比化、つまり幅が狭く穴の方向に深い状態が進むため、均一な成膜はますます難しくなっており、成膜プロセスの反応機構解析からより高アスペクト比に対応した新プロセスを、戦略的に、短期間で、且つ低コストで実現することが切望されている。
【０００６】
そこで、この出願の発明は、以上のとおりの事情に鑑み、従来技術の問題点を解消し、成膜等の化学反応の反応機構解析を短期間で、且つ低コストで実現することのできる、新しい反応機構解析プログラムおよび当該プログラムを記録したコンピュータ読み取り可能な記録媒体、ならびに反応機構解析装置を提供することを課題としている。
【０００７】
【課題を解決するための手段】
この出願の発明は、上記の課題を解決するものとして、第１には、１種または複数種の反応器における化学反応の反応機構を解析するためにコンピュータを、実験により得られた成膜形状の実験成膜形状指針値を算出する実験値算出手段、数値シミュレーションによって、実験条件値を変化させながら複数のシミュレート成膜形状を生成する成膜形状シミュレーション手段、各シミュレート成膜形状のシミュレート成膜形状指針値を算出するシミュレーション値算出手段、各シミュレート成膜形状の実験条件値およびシミュレート成膜形状指針値に基づいて、実験条件と成膜形状指針との関係を表す実験条件−成膜形状指針モデルを生成するモデル化手段、および進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算する進化手段、として機能させることを特徴とする反応機構解析プログラムを提供する。
【０００８】
また、この出願の発明は、１種または複数種の反応器における化学反応の反応機構を解析するためにコンピュータを、実験により得られた成膜形状の実験成膜形状指針値を算出する実験値算出手段、数値シミュレーションによって、実験条件値を変化させながら複数のシミュレート成膜形状を生成する成膜形状シミュレーション手段、各シミュレート成膜形状のシミュレート成膜形状指針値を算出するシミュレーション値算出手段、各シミュレート成膜形状の実験条件値およびシミュレート成膜形状指針値に基づいて、実験条件と成膜形状指針との関係を表す実験条件−成膜形状指針モデルを生成するモデル化手段、および進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算する進化手段、として機能させることを特徴とする反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体を提供する。
【０００９】
また、第３には、１種または複数種の反応器における化学反応の反応機構を解析する装置であって、実験により得られた成膜形状の実験成膜形状指針値を算出する実験値算出手段、数値シミュレーションによって、実験条件値を変化させながら複数のシミュレート成膜形状を生成する成膜形状シミュレーション手段、各シミュレート成膜形状のシミュレート成膜形状指針値を算出するシミュレーション値算出手段、各シミュレート成膜形状の実験条件値およびシミュレート成膜形状指針値に基づいて、実験条件と成膜形状指針との関係を表す実験条件−成膜形状指針モデルを生成するモデル化手段、および進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算する進化手段、を備えることを特徴とする反応機構解析装置を提供する。
【００１０】
そして、第４には、前記進化手段は、進化論的処理によって、実験条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値および前記実験時の実験条件から直接数値シミュレーションによって算出されるシミュレート成膜形状指針値が前記実験成膜形状指針値と一致するように、最適な反応機構を計算することを特徴とする請求項１記載の反応機構解析プログラムまたは請求項２記載の反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体または請求項３記載の反応機構解析装置、第５には、前記反応機構が、反応物が中間体を経て生成物にいたる気相および表面における定性的な反応経路、反応物および中間体で構成される成膜種の付着確率、成膜種の表面反応速度定数、成膜種の気相反応速度定数、および成膜種の濃度、のうちの少なくともいずれかであることを特徴とする成膜反応機構解析プログラムまたは成膜反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体または成膜反応機構解析装置、第６には、前記成膜形状指針が、基板上の任意の複数位置における膜厚の相互関係を表す関数であることを特徴とする成膜反応機構解析プログラムまたは成膜反応機構解析プログラムを記録したコンピュータ読み取り可能な記録媒体または成膜反応機構解析装置を提供する。
【００１１】
【発明の実施の形態】
上記のとおりの特徴を有するこの出願の発明は、コンピュータを用いて、１種または複数種の反応器における成膜等の化学反応の反応機構をコンピュータ演算によってシミュレーション解析するものであり、以下に、機能ブロックおよび処理フローを示した図１および図２に沿って、各処理について説明する。なお、図３はトレンチへの成膜過程における１成膜種の付着確率ηを解析する場合の一例を示したものであり、この図３をも適宜参酌して説明する。
【００１２】
＜実験値算出処理＞
まず、実験により実際にＶＩＡホールやトレンチ等を有する基板上へ成膜して得られた成膜形状に対して、その成膜形状指針の値（実験成膜形状指針値と呼ぶこととする）を算出する（実験値算出手段（１）、ステップＳ１）。
【００１３】
より具体的には、たとえば実験により得られた成膜形状をＳＥＭ等で撮影し、その画像からステップカバレッジ等の成膜形状の指針となる値を算出する。成膜形状画像からの成膜形状指針値の算出は公知手法で可能である。
【００１４】
たとえば図４に例示したような成膜形状におけるステップカバレッジ（段差被膜性とも呼ぶ）の値は、下記数１で求めることができる（たとえば、菅野卓雄　川西剛監修，「半導体大辞典」，初版，工業調査会，１９９９年，ｐ．１９３、Ｈｅｅ　Ｊｏｏｎ　Ｋｉｍ　他２名，「Ｍｉｃｒｏ−Ｔｒｅｎｃｈ法によるＡＰＣＶＤにおける表面反応速度定数の測定」，化学工学論文集，１９９１年，第１７巻，第６号，ｐ．１１７５−１１７７参照）。また、たとえば図５（トレンチ開口部に“ひさし”が出っ張ったもの）および図６（“す”が生じたもの）に例示したような成膜形状の場合においては、それぞれ下記数２および数３を成膜形状指針とすることができる。
【００１５】
【数１】
2004119742.tif 000003
【００１６】
【数２】
2004119742.tif 000004
【００１７】
【数３】
2004119742.tif 000005
【００１８】
これらの数式のように、ステップカバレッジ等の成膜形状指針は、たとえば基板上の任意の複数位置における膜厚の相互関係を表す関数とすることができる。なお、必要な各膜厚データについては、予め成膜形状画像からユーザが計測してデータ入力（入力手段（６））するようにしても、あるいは成膜形状画像データを入力（入力手段（６））して画像処理により自動算出するようにしてもよい。
【００１９】
＜成膜形状シミュレーション処理＞
数値シミュレーションによって、実験条件値を変化させながら、各実験条件値に対応するシミュレート成膜形状を生成する（成膜形状シミュレーション手段（２）、ステップＳ２）。
【００２０】
より具体的には、たとえば前記実験により得られた成膜形状のＳＥＭ画像等から基板上のＶＩＡホール形状やトレンチ形状等の基板形状データを決定し、成膜種の種類、濃度、付着確率（または表面反応速度）、平均自由行程の大きさなどといった実験条件の値を系統的に変化させてモンテカルロシミュレーション等の数値シミュレーションを行い、各実験条件値に対応する成膜形状を生成する。
【００２１】
実験条件値は、ユーザが任意に選択して入力（入力手段（６））するようにしても、コンピュータにより自動変化させるようにしてもよい。また、基板形状データは、予めユーザが決定した形状データを入力（入力手段（６））するようにしても、実際のＳＥＭ画像等の画像データを入力（入力手段（６））して画像処理により自動抽出するようにしてもよい。
【００２２】
図３の例では、予め成膜種の種類を一つと仮定しておき、ＳＥＭ像から決定したトレンチ形状データに対して１成膜種の付着確率ηを系統的に変化させてモンテカルロシミュレーションを行い、各付着確率ηに対応する成膜形状をトレンチ上に表した画像データを生成している。
【００２３】
＜シミュレーション値算出処理＞
続いて、上記のとおりに生成された各シミュレート成膜形状について、その成膜形状指針の値（シミュレート成膜形状指針値と呼ぶこととする）を算出する（シミュレーション値算出手段（３）、ステップＳ３）。
【００２４】
より具体的には、たとえば前記実験値算出処理（ステップＳ１）にて算出した実験成膜形状指針と同じ成膜形状指針の値を、各シミュレート成膜形状について算出する。ステップＳ１にて実験成膜形状指針値として数１によるステップカバレッジを算出していれば、本処理でも同じ数１によるステップカバレッジを算出する。
【００２５】
図３の例では、上述した各画像データについて、必要な膜厚データを自動算出して前記数１によりステップカバレッジＳｔを求めている。図中のグラフは、モンテカルロシミュレーションにて系統的に変化させた付着確率ηとそれに対応する各画像データごとのステップカバレッジＳｔとの関係を示している（グラフ中の○印）。
【００２６】
＜モデル化処理＞
ここで、上記各シミュレート成膜形状の実験条件値およびシミュレート成膜形状指針値に基づいて、実験条件と成膜形状指針との関係を表す成膜条件−成膜形状指針モデルを生成する（モデル化手段（４）、ステップＳ４）。
【００２７】
より具体的には、前記成膜形状シミュレーション処理（ステップＳ２）にて用いた実験条件値および上記シミュレーション値算出処理（ステップＳ３）にて求めたシミュレート成膜形状指針値とを基に、実験条件と成膜形状指針との関係を、多項式等の関数、解析解、ルックアップテーブル、ＮＮ（ニューラルネットワーク）などでモデル化する。このモデルは記憶手段（図示なし）に記憶しておき、後述の進化処理において随時読み出すようにする。
【００２８】
図３の例では、実験条件である１成膜種の付着確率ηと成膜形状指針であるステップカバレッジＳｔとの関係に対して下記数４のとおりの多項式モデルを得ており、グラフ中にてシミュレーション値（○印）とモデル曲線（実線）とが一致していることを示している。
【００２９】
【数４】
2004119742.tif 000006
【００３０】
＜進化処理＞
この後、進化論的処理によって、上記成膜条件−成膜形状指針モデルを用いて算出されるシミュレート成膜形状指針値が実験成膜形状指針値と一致するように、反応機構を進化させる（進化手段（５）、ステップＳ５・Ｓ６）。
【００３１】
より具体的には、遺伝的アルゴリズム等の「進化論的」（「遺伝的」とも呼べる）な探索技法によって、反応機構の情報を染色体（または遺伝子）に記述し、、成膜条件−成膜形状指針モデルにより算出されるシミュレート成膜形状指針値が上記実験成膜形状指針値と一致するように（たとえば一致度またはズレが許容範囲となるように）、染色体を進化させていく。そして、何回かの進化の末に両値が一致したときに、その世代の染色体が示す反応機構が正しい反応機構、つまり実際の成膜形状をもたらす反応機構として求められたことになるのである。
【００３２】
反応機構としては、反応物が中間体を経て生成物にいたる気相および表面における定性的な反応経路、反応物および中間体で構成される成膜種の付着確率、成膜種の表面反応速度定数、成膜種の気相反応速度定数、および成膜種の濃度などが挙げられる。
【００３３】
図３の例では、１成膜種の付着確率ηｓｉｍを染色体データとして設定し、まず複数の初代染色体を乱数を使って決め、各初代染色体の付着確率値ηｓｉｍ−１ごとにステップカバレッジＳｔｓｉｍ−１を算出して実測成膜形状のステップカバレッジ値Ｓｔｒｅａｌとの一致を自動判断（数値比較処理で可能）し、一致しない場合は初代染色体を次世代染色体へ進化させて、次世代染色体の付着確率値ηｓｉｍ−２で再びステップカバレッジＳｔｓｉｍ−２を算出してステップカバレッジ値Ｓｔｒｅａｌとの一致を自動判断し、一致すると判断されるまで進化を繰り返す。そして、一致したときの世代の染色体が示す付着確率値ηｓｉｍ−ｉが、実際のトレンチ上に成膜されたステップカバレッジＳｔｒｅａｌを得ることのできる反応機構ということになる。
【００３４】
なお、この出願の発明の発明者等は、既に反応機構解析に進化論的処理を用いた技術を提案しており（特願２００１−６９３９５号参照）、そこでは反応経路をトポロジカル染色体と呼ぶ染色体に表し、表面・気相反応速度定数を非トポロジカル染色体と呼ぶ染色体に表し、成膜種の濃度をテンポラリ染色体と呼ぶ染色体に表して、それら三種の染色体を入れ子状にして進化させていている。この出願の発明でもこの原理を利用することができ、たとえば上記図３の場合では成膜に関わる成膜種を予め１種類と仮定しているので反応経路を表すトポロジカル染色体は既に固定されているが、１成膜種の付着確率を非トポロジカル染色体とし、その濃度をテンポラリ染色体とすることで、入れ子状構造の各染色体を進化させることができる。
【００３５】
このような進化の具体的処理としては、選択、交叉、突然変異など、さらにはエリート戦略（エリート保存とも言う）、トーナメント方式選択、一様交叉などの各種公知の進化論的処理を用いることができ、交叉確率、突然変異確率、エリート戦略で残す染色体数などの処理条件は予め設定入力しておく（入力手段（６））。
【００３６】
以上のこの出願の発明によれば、ＶＩＡホールやトレンチ等を有する基板上への成膜過程に関わる様々な反応機構を、短期間且つ低コストで、しかも高精度でコンピュータ解析することができ、たとえば、１種の反応器としてのマクロキャビティ（Ｍａｃｒｏｃａｖｉｔｙ）での成膜反応機構はもちろんのこと、マクロキャビティとマイクロキャビティ（Ｍｉｃｒｏｃａｖｉｔｙ）という異なる２種の反応器を組み合わせた場合（マイクロ／マクロキャビティ：Ｍｉｃｒｏ／Ｍａｃｒｏｃａｖｉｔｙ）の成膜反応機構についても、コンピュータを用いた自動解析を実現することができるのである。
【００３７】
前述したように従来では、シミュレート成膜形状を実際の成膜形状に一致するまで変化させることをモンテカルロシミュレーションで行っていたため（図１２参照）、変化の度に計算コストの高いモンテカルロシミュレーションが必要となり、非常に非効率的であったのに対し、この出願の発明では、実際の成膜形状に一致するシミュレート成膜形状の探索は遺伝的アルゴリズム等の進化論的処理で行い、モンテカルロシミュレーション等の数値シミュレーションはその前処理としてのモデル化、つまりシミュレート成膜形状を計算するためのモデルを作成する処理に用いるだけであるため、モデル化に必要な少数のシミュレート成膜形状を生成する際にモンテカルロシミュレーション等の数値シミュレーションを行えばよく、全体の計算コストが低減され（たとえば１００分の１以下）、効率良く解析を実行できるのである。
【００３８】
またもちろん、この進化論的処理においては、実験条件−成膜形状指針モデルを用いるだけでなく、従来と同様な数値シミュレーションを行うことも併用でき、具体的には、実験条件−成膜形状指針モデルによって算出されるシミュレート成膜形状指針値（ステップカバレッジなど）と、実験条件つまり前記実験値算出処理において説明した実際の実験にて用いられた条件と同じ条件から直接数値シミュレーションによって算出されるシミュレート成膜形状指針値（膜厚分布など）との両値が、前記実験成膜形状指針値と一致するように、最適な反応機構を計算することになる。この処理は、進化手段（５）が全て行うようにしても、あるいは実験条件から直接数値シミュレーションによってシミュレート成膜形状指針値を算出する別途のシミュレート成膜形状指針値算出手段（図示ナシ）を設けてその算出値を進化論的処理時に進化手段（５）へ与えるように構築してもよい。このような併用処理の場合であっても従来と比較した計算コストの低減を図ることができ、より一層精度の高い解析を実現できる。
【００３９】
そして、以上のように優れたこの出願の発明は、成膜だけではなく様々な反応器における様々な化学反応についても適用できるものであり、同様に短時間、低コスト、高精度の反応機構解析を実現できる。
【００４０】
たとえば、上記付着確率η（または表面反応速度定数）に負の値をとることを許すという条件を設定すれば、成膜過程の一部にエッチング過程を含んだ場合の化学反応を扱え、また前記数１のＴｂ，Ｔｓおよびマクロキャビティの成膜速度に負の値をとることを許せば、全ての過程がエッチング過程である場合をも扱える。すなわち、成膜の形成方向に対して負の方向に寄与する処理を実行可能なものとすることで、エッチング現象を扱えるようになるのである。
【００４１】
【実施例】
［実施例１］
上述したこの出願の発明により、図７のマイクロ／マクロキャビティにおけるＣＶＤ法による成膜を仮想的に設定し、そのときの成膜に関わった反応機構をシミュレート解析した。マクロキャビティおよびマイクロキャビティの支配方程式はそれぞれ下記数５および数６で与えられる。
【００４２】
【数５】
2004119742.tif 000007
【００４３】
【数６】
2004119742.tif 000008
【００４４】
本実施例では、反応機構が気相および表面における反応物および中間体の一次反応のつながりで表現されるものと設定し、成膜形状指針として、マクロキャビティ中心部（ｘ＝０）におけるマイクロキャビティでのステップカバレッジＳｔ（ｘ）（数５および数６参照）を用い、またこの成膜形状指針と併用して、マクロキャビティ内部における膜厚分布を用いる。
【００４５】
図８内の表は、仮想的に設定した成膜形状（これが前記実験により得られた成膜形状になる）の成膜形状指針としての上記ステップカバレッジＳｔの仮想実験データを示しており、図８内のグラフは、これと併用されるマクロキャビティ内部における膜厚分布の仮想実験データを示している。
【００４６】
また、この仮想成膜形状の仮想実験データの作成に使用した仮想反応機構、ならびにマクロ／マイクロキャビティにおける各値（図中の表）は、図９に示したとおりである。この図９に示した仮想反応機構が求めるべき正解の仮想反応機構ということになり、上記仮想実験データとシミュレートされたステップカバレッジＳｔと膜厚分布値とが一致するように反応機構つまり染色体を遺伝的アルゴリズムにより進化させ、この進化により得られた反応機構が図９の仮想反応機構と十分に一致すれば、正しい反応機構を解析できたことになる。
【００４７】
なお、進化処理前のモデル化では、マクロ／マイクロキャビティにおけるトレンチの形状を矩形のものに設定しているので（図９の表参照）、下記数７で与えられるクヌーセン拡散域における拡散方程式を解くことで得られる解析解を実験条件−成膜形状指針モデルとして用いた。縦横比つまりアスペクト比ｌ／ｗは２０に設定している（図９の表参照）。
【００４８】
【数７】
2004119742.tif 000009
【００４９】
図１０および図１１はそれぞれ解析結果を示したものであり、図１０は染色体の各世代の評価値の最大値とメジアンを示し、図１１は２８世代目で最良（最大）の評価値を持つ染色体が示す反応機構情報を示している。図１１内の表は本実施例で探索対象とした各値の範囲である。この図１０に示したように２８世代目の染色体にて上記仮想実験データとシミュレートデータとのズレが許容範囲のものとなり、この２８世代目で最良の評価値を持つ染色体にて表された反応機構情報は図１１のとおりで、これと図９の仮想反応機構情報を比較すると両者が定性的にも定量的にも十分に一致していることがわかる。したがって、仮想成膜形状を実現できる正しい反応機構が正確に解析できたのである。
【００５０】
以上の結果から、この出願の発明により、図７のマイクロ／マクロキャビティにおけるＣＶＤ法による成膜の反応機構を高精度で解析できたことがわかる。もちろん計算量も従来に比べて格段に低減している。
【００５１】
この出願の発明は以上の実施形態および実施例に限定されるものではなく、細部については様々な態様が可能である。たとえばマイクロ／マクロキャビティにおけるＣＶＤ法による成膜についての実施例のみを示したが、その他の各種公知の反応器や成膜手法、さらには成膜以外の化学反応にも適用可能であることは言うまでもない。たとえば前記数５及び数６の支配方程式における付着確率ηまたは表面反応速度定数ｋｓに負の値を許容し、または成膜速度Ｇｒ（ｘ）および前記数１のＴｂ，Ｔｓに負の値を許容すれば、エッチング現象を扱うことができる。
【００５２】
【発明の効果】
以上詳しく説明したとおり、この出願の発明によって、様々な反応器における成膜等の化学反応の反応機構の解析を短期間で、且つ低コストで実現することのできる、新しい反応機構解析プログラムおよび当該プログラムを記録したコンピュータ読み取り可能な記録媒体、ならびに反応機構解析装置が提供される。
【図面の簡単な説明】
【図１】この出願の発明を説明するための機能ブロック図である。
【図２】この出願の発明を説明するための処理フロー図である。
【図３】この出願の発明によりトレンチへの成膜過程における成膜種の付着確率ηを解析する場合の一例を説明するための図である。
【図４】成膜形状指針を説明するための図である。
【図５】別の成膜形状指針を説明するための図である。
【図６】さらに別の成膜形状指針を説明するための図である。
【図７】実施例１で用いたマイクロ／マクロキャビティを説明するための図である。
【図８】実施例１で用いた仮想実験データを示した図である。
【図９】実施例１で用いた仮想反応機構を示した図である。
【図１０】実施例１における解析結果を示した図である。
【図１１】実施例１における解析結果を示した別の図である。
【図１２】従来手法によりトレンチへの成膜過程における成膜種の付着確率ηを解析する場合の一例を説明するための図である。
【符号の説明】
１　実験値算出手段
２　成膜形状シミュレーション手段
３　シミュレーション値算出手段
４　モデル化手段
５　進化手段
６　入力手段
７　出力手段
【図１】
2004119742.tif 000010
【図２】
2004119742.tif 000011
【図３】
2004119742.tif 000012
【図４】
2004119742.tif 000013
【図５】
2004119742.tif 000014
【図６】
2004119742.tif 000015
【図７】
2004119742.tif 000016
【図８】
2004119742.tif 000017
【図９】
2004119742.tif 000018
【図１０】
2004119742.tif 000019
【図１１】
2004119742.tif 000020
【図１２】
2004119742.tif 000021