特開2023-180839成膜シミュレーション方法、成膜シミュレーションプログラム、成膜シミュレータおよび成膜装置
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023180839
(43)【公開日】2023-12-21
(54)【発明の名称】成膜シミュレーション方法、成膜シミュレーションプログラム、成膜シミュレータおよび成膜装置
(51)【国際特許分類】
C23C 16/52 20060101AFI20231214BHJP
C23C 14/34 20060101ALI20231214BHJP
H01L 21/31 20060101ALI20231214BHJP
【ＦＩ】
C23C16/52
C23C14/34 U
H01L21/31 B
H01L21/31 C
H01L21/31 D
【審査請求】未請求
【請求項の数】15
【出願形態】ＯＬ
(21)【出願番号】P 2022094465
(22)【出願日】2022-06-10
(71)【出願人】
【識別番号】316005926
【氏名又は名称】ソニーセミコンダクタソリューションズ株式会社
(74)【代理人】
【識別番号】110001357
【氏名又は名称】弁理士法人つばさ国際特許事務所
(72)【発明者】
【氏名】久保井 信行
【テーマコード（参考）】
4K029
4K030
5F045
【Ｆターム（参考）】
4K029CA01
4K029CA05
4K029DA03
4K029EA06
4K030HA12
4K030KA39
4K030KA41
5F045AA04
5F045AA05
5F045AA08
5F045AA15
5F045AA18
5F045AA19
5F045GB17
(57)【要約】
【課題】計算時間の短縮化および計算精度の向上を図ることの可能な成膜シミュレーション方法、成膜シミュレーションプログラムおよび成膜シミュレータと、このような成膜シミュレータを備えた成膜装置を提供する。
【解決手段】本開示の一実施の形態に係る成膜シミュレーション方法は、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積を計算するステップを含む。この成膜シミュレーション方法は、上記ステップにおいて、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積を計算することで成膜表面における膜の被覆性および膜質を表現することを含む。
【選択図】図１
【解決手段】本開示の一実施の形態に係る成膜シミュレーション方法は、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積を計算するステップを含む。この成膜シミュレーション方法は、上記ステップにおいて、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積を計算することで成膜表面における膜の被覆性および膜質を表現することを含む。
【選択図】図１
【特許請求の範囲】
【請求項1】
入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算する第１ステップを含み、
前記第１ステップにおいて、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現することを含む
成膜シミュレーション方法。
前記第１ステップにおいて、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現することを含む
成膜シミュレーション方法。
【請求項2】
前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して直接入射してくるガス成分と、前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して周囲構造に応じて入射してくるガス回り込み成分とにより、前記入射ラジカルフラックスを算出する第２ステップを更に含む
請求項１に記載の成膜シミュレーション方法。
請求項１に記載の成膜シミュレーション方法。
【請求項3】
前記第１ステップにおいて、前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して入射するイオンのフラックスおよびエネルギーに基づいて、前記代表粒子と前記成膜表面との結合および未結合の判断を行うことを更に含む
請求項１に記載の成膜シミュレーション方法。
請求項１に記載の成膜シミュレーション方法。
【請求項4】
前記第１ステップにおいて、前記状態情報を用いて、膜密度、透水性および前記成膜表面との密着性を計算することを更に含む
請求項１に記載の成膜シミュレーション方法。
請求項１に記載の成膜シミュレーション方法。
【請求項5】
前記第１ステップにおいて、複数種類のガスを用いて成膜を行う場合、同じ時間ステップ内でガスごとの前記代表粒子の生成個数を算出することを更に含む
請求項１に記載の成膜シミュレーション方法。
請求項１に記載の成膜シミュレーション方法。
【請求項6】
前記第１ステップにおいて、複数種類のガスを用いて成膜を行う場合、ガスごとに別々の時間ステップで前記代表粒子の生成個数を算出することを更に含む
請求項１に記載の成膜シミュレーション方法。
請求項１に記載の成膜シミュレーション方法。
【請求項7】
前記第１ステップにおいて、成膜後に、アニールを施した時の膜密度およびブリスターの少なくとも一方を計算することを更に含む
請求項１に記載の成膜シミュレーション方法。
請求項１に記載の成膜シミュレーション方法。
【請求項8】
成膜条件を取得する入力部と、
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレーションプログラム。
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレーションプログラム。
【請求項9】
前記入力部は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成され、
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
請求項８に記載の成膜シミュレーションプログラム。
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
請求項８に記載の成膜シミュレーションプログラム。
【請求項10】
成膜条件を取得する入力部と、
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレータ。
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレータ。
【請求項11】
前記入力部は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成され、
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
請求項１０に記載の成膜シミュレータ。
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
請求項１０に記載の成膜シミュレータ。
【請求項12】
成膜チャンバと、
前記成膜チャンバの動作を制御する制御部と、
入射ラジカルフラックスの値に応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算し、それにより得られた計算結果に基づいて、成膜プロセスの最適化条件を探索する最適化演算部と、
前記成膜チャンバの成膜条件が、前記最適化演算部で見つかった前記最適化条件となるのに必要なデータを生成し、前記制御部へ出力する出力部と
を備え、
前記最適化演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜装置。
前記成膜チャンバの動作を制御する制御部と、
入射ラジカルフラックスの値に応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算し、それにより得られた計算結果に基づいて、成膜プロセスの最適化条件を探索する最適化演算部と、
前記成膜チャンバの成膜条件が、前記最適化演算部で見つかった前記最適化条件となるのに必要なデータを生成し、前記制御部へ出力する出力部と
を備え、
前記最適化演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜装置。
【請求項13】
前記成膜チャンバ内の状態をモニタリングするモニタリング装置と、
前記モニタリング装置により得られたモニタリングデータを用いて前記入射ラジカルフラックスを計算するガスフラックス演算部と
を更に備えた
請求項１２に記載の成膜装置。
前記モニタリング装置により得られたモニタリングデータを用いて前記入射ラジカルフラックスを計算するガスフラックス演算部と
を更に備えた
請求項１２に記載の成膜装置。
【請求項14】
前記最適化演算部は、データベースを利用して予測した前記入射ラジカルフラックスの値に応じた代表粒子を生成する
請求項１２に記載の成膜装置。
請求項１２に記載の成膜装置。
【請求項15】
前記最適化演算部は、ウェハ毎、または、ロット毎に、前記最適化条件を探索する
請求項１２に記載の成膜装置。
請求項１２に記載の成膜装置。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、成膜処理の際の加工表面（成膜対象の成膜面）の形状をシミュレーションする成膜シミュレーション方法、この成膜シミュレーション方法を実行する成膜シミュレーションプログラムおよび成膜シミュレータに関する。また、本開示は、成膜シミュレータを備えた成膜装置に関する。
【背景技術】
【0002】
半導体デバイス製造のキープロセスの１つとして、パターン上に薄膜(ｎｍオーダー)を形成するないしは埋め込む成膜プロセスがある。近年の半導体デバイスの高機能化の要求の高まりにより、高低アスペクト比が混在する複雑な積層構造を有するようになってきており、成膜プロセスにおけるこれらのパターンの埋め込みの難易度も高くなってきている。そのため、プロセス条件に依存した膜の被覆性（カバレッジ）および膜質（例えば、密度、欠陥密度、透水性、および密着性など）の予測の重要性がますます増加している。このような状況において、予測技術の１つとして、成膜プロセスの数値シミュレーションが有用である。
【0003】
例えば、特許文献１では、モンテカルロ法を用いて、パターン表面へのガス粒子の付着および脱離、基板温度に依存したマイグレーション、ならびに表面凹凸の影響を考慮した粒子の堆積を計算し、膜のモフォロジーを予測する。さらに、あらかじめＭＤ(Molecular Dynamics)計算や第一原理計算で構築しておいたモフォロジーに紐づいた膜質のデータベースを用いて膜質分布を予測する方法を提案している。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特表２０１７－２０４７５７号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかし、特許文献１に記載の方法では、モンテカルロ法による照射ガス計算のため、モンテカルロ法に用いる粒子数に精度と計算時間が大きく依存してしまう。さらに、膜質計算には、モフォロジーとガスフラックスに依存したデータベースをＭＤ(Molecular Dynamics)計算や第一原理計算によって別途準備する必要がある。そのため、実際にシミュレーションツールとして用いるまでに、その準備にかなりの時間と手間がかかってしまう。従って、計算時間の短縮化および計算精度の向上を図ることの可能な成膜シミュレーション方法、成膜シミュレーションプログラムおよび成膜シミュレータと、このような成膜シミュレータを備えた成膜装置とを提供することが望ましい。
【課題を解決するための手段】
【0006】
本開示の一実施の形態に係る成膜シミュレーション方法は、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積を計算するステップを含む。この成膜シミュレーション方法は、上記ステップにおいて、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積を計算することで成膜表面における膜の被覆性および膜質を表現することを含む。
【0007】
本開示の一実施の形態に係る成膜シミュレーションプログラムは、入力部と、演算部と、出力部とを備えている。入力部は、成膜条件を取得する。演算部は、入力部で取得した成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積を計算する。出力部は、演算部での演算結果を出力する。演算部は、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積を計算することで成膜表面における膜の被覆性および膜質を表現する。
【0008】
本開示の一実施の形態に係る成膜シミュレータは、入力部と、演算部と、出力部とを備えている。入力部は、成膜条件を取得する。演算部は、入力部で取得した成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積を計算する。出力部は、演算部での演算結果を出力する。演算部は、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積を計算することで成膜表面における膜の被覆性および膜質を表現する。
【0009】
本開示の一実施の形態に係る成膜装置は、成膜チャンバと、成膜チャンバの動作を制御する制御部と、最適化演算部と、出力部とを備えている。最適化演算部は、入射ラジカルフラックスの値に応じた代表粒子を生成し、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積を計算する。最適化演算部は、それによって得られた計算結果に基づいて、成膜プロセスの最適化条件を探索する。出力部は、成膜チャンバの成膜条件が、最適化演算部で見つかった最適化条件となるのに必要なデータを生成し、制御部へ出力する。最適化演算部は、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積を計算することで成膜表面における膜の被覆性および膜質を表現する。
【図面の簡単な説明】
【0010】
【図1】本開示の一実施の形態に係る成膜シミュレーション方法における処理手順の一例を表す図である。
【図3】本開示の一実施の形態に係る成膜シミュレータの構成例を表す図である。
【図4】ガス輸送の概念を表す図である。
【図5】実際の成膜表面の概念を表す図である。
【図7】１Ｖｏｘｅｌでの代表粒子の生成について説明するための図である。
【図8】マイグレーションにおける実際の成膜表面の概念を表す図である。
【図9】マイグレーションにおけるＶｏｘｅｌモデルの概念を表す図である。
【図10】Ｖｏｘｅｌ空間内でのＣｒａｄｌｅの判定方法の一例を表す図である。
【図11】Ｖｏｘｅｌ空間内でのＣｒａｄｌｅの判定方法の一例を表す図である。
【図12】Ｖｏｘｅｌ空間内でのＣｒａｄｌｅの判定方法の一例を表す図である。
【図13】Ｖｏｘｅｌ空間内での膜密度および透水性の算出について説明するための図である。
【図14】Ｖｏｘｅｌ空間内での密着性の算出について説明するための図である。
【図15】Ｖｏｘｅｌ空間内でのアニール処理について説明するための図である。
【図16】Ｖｏｘｅｌ空間内でのアニール処理について説明するための図である。
【図17】Ｖｏｘｅｌ空間内でのブリスターについて説明するための図である。
【図18】Ｖｏｘｅｌ空間内でのブリスターについて説明するための図である。
【図19】Ｖｏｘｅｌ空間内でのブリスターについて説明するための図である。
【図20】Ｖｏｘｅｌ空間内でのブリスターについて説明するための図である。
【図21】ＣＶＤプロセスによるカバレッジと膜密度分布の計算結果を表す図である。
【図22】Ｖｏｘｅｌの結合状態/未結合状態の分布、ガスフラックス分布の結果を表す図である。
【図23】透水性分布の計算結果を表す図である。
【図24】密着性分布の計算結果を表す図である。
【図25】アニール後の膜密度分布の計算結果を表す図である。
【図26】アニール後の密着性分布の計算結果を表す図である。
【図27】アニール後の透水性分布の計算結果を表す図である。
【図28】ブリスター計算の概念を表す図である。
【図29】トレンチの交差構造例を表す図である。
【図31】本開示の一実施の形態に係る成膜シミュレーションプログラムの一構成例を表す図である。
【図32】本開示の一実施の形態に係る成膜装置の一構成例を表す図である。
【発明を実施するための形態】
【0011】
以下、本開示を実施するための形態について、図面を参照して詳細に説明する。
【0012】
＜実施の形態＞
＜概略＞
まず、本開示の一実施形態に係る成膜シミュレーション方法の概略について説明する。本実施形態に係る成膜シミュレーション方法では、加工表面（成膜表面）に対して、原料粒子を投射し、原料粒子からなる膜を形成する成膜方法を扱う。
＜概略＞
まず、本開示の一実施形態に係る成膜シミュレーション方法の概略について説明する。本実施形態に係る成膜シミュレーション方法では、加工表面（成膜表面）に対して、原料粒子を投射し、原料粒子からなる膜を形成する成膜方法を扱う。
【0013】
具体的には、本実施形態に係る成膜シミュレーション方法は、各種蒸着法を扱い、成膜される膜の被覆性および膜質を予測するものである。本実施形態に係る成膜シミュレーション方法で扱うことが可能な成膜方法としては、例えば、抵抗加熱蒸着法、電子ビーム蒸着法、分子線エピタキシー法、イオンプレーティング法、およびスパッタリング法などの物理蒸着法（Physical Vapor Deposition :PVD）ならびに熱またはプラズマ化学気相蒸着法、原子層堆積法（Atomic Layer Deposition :ADL）、有機金属気相成長法などの化学蒸着法(Chemical Vapor Deposition :CVD）などを例示することができる。
【0014】
原料粒子は、例えば、原子、分子、またはこれらを電離させたイオンである。原料粒子は、成膜チャンバに導入された原料ガスを熱またはプラズマ等を用いて分解または電離することで形成されてもよく、金属ターゲットに希ガス原子等を衝突させることで形成されてもよい。また、原料粒子は、１種であってもよく、２種以上であってもよい。すなわち、膜は、単一の原料から形成された膜であってもよく、複数の原料を反応させることで形成された膜であってもよい。
【0015】
成膜対象は、例えば、金属基板、半導体基板、ガラス基板、石英基板、または樹脂基板などである。成膜対象の表面の形状および材質は、特に限定されない。例えば、成膜対象の表面には、薄膜が形成されていてもよく、微細構造が形成されていてもよい。成膜対象に成膜される膜は、例えば、数マイクロメートル程度の膜厚の薄膜である。また、本実施形態に係る成膜シミュレーション方法にて扱うことが可能な領域の大きさは、例えば、一辺の長さが数マイクロメートル程度の領域である。
【0016】
本実施形態に係る成膜シミュレーション方法によれば、上述した成膜方法において、成膜される膜の被覆性（カバレッジ）および膜質を数十ｎｍの範囲にて予測することが可能である。
【0017】
＜成膜シミュレーション方法の流れ＞
続いて、図１、図２を参照して、本実施形態に係る成膜シミュレーション方法の流れの概要について説明する。図１は、本実施形態に係る成膜シミュレーション方法の流れの概要を示すフローチャート図である。図２は、図１に記載の成膜シミュレーション方法の流れをより詳細に示すフローチャート図である。図３は、この成膜シミュレーション方法を実現するための情報処理装置（成膜シミュレータ）の一構成例を示す図である。
続いて、図１、図２を参照して、本実施形態に係る成膜シミュレーション方法の流れの概要について説明する。図１は、本実施形態に係る成膜シミュレーション方法の流れの概要を示すフローチャート図である。図２は、図１に記載の成膜シミュレーション方法の流れをより詳細に示すフローチャート図である。図３は、この成膜シミュレーション方法を実現するための情報処理装置（成膜シミュレータ）の一構成例を示す図である。
【0018】
（成膜シミュレータの構成）
図３に示した成膜シミュレータ１は、入力部１１と、演算部１２と、出力部１３とを備えている。入力部１１は、成膜表面に対して所定の成膜処理を行う際の成膜条件を取得して演算部１２に入力するものである。入力部１１は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成されている。演算部１２は、入力部１１を介して入力された成膜条件に基づいて、後述する図１，図２に示したシミュレーション方法によって、成膜表面の形状進展や膜質の計算を行う。
図３に示した成膜シミュレータ１は、入力部１１と、演算部１２と、出力部１３とを備えている。入力部１１は、成膜表面に対して所定の成膜処理を行う際の成膜条件を取得して演算部１２に入力するものである。入力部１１は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成されている。演算部１２は、入力部１１を介して入力された成膜条件に基づいて、後述する図１，図２に示したシミュレーション方法によって、成膜表面の形状進展や膜質の計算を行う。
【0019】
なお、本実施の形態では、演算部１２をハードウェアで構成して後述する計算処理を実現してもよいが、所定のシミュレーションプログラム（ソフトウェア）を用いて計算処理を実行してもよい。この場合、演算部１２を、例えばＣＰＵ（Central Processing Unit）等の演算装置で構成し、シミュレーションプログラムを外部から読み込み、そのプログラムを実行することにより計算を実行する。
【0020】
シミュレーションプログラムは、例えば、図示しないデータベースや、別途設けられる例えばＲＯＭ（Read Only Memory）等の記憶部などに格納することができる。この際、シミュレーションプログラムを、例えばデータベースや別途設けられた記憶部等に予め実装した構成にしてもよいし、外部から例えばデータベースや別途設けられた記憶部等に実装する構成にしてもよい。なお、外部からシミュレーションプログラムを取得する場合には、シミュレーションプログラムを、光ディスクや半導体メモリなどの媒体から配布するようにしてもよいし、インターネットなどの伝送手段を介してダウンロードするようにしてもよい。
【0021】
出力部１３は、演算部１２での演算結果（演算部１２によって計算された所定の成膜処理のシミュレーション結果）を出力するものである。出力部１３は、演算部１２での演算結果（演算部１２によって計算された所定の成膜処理のシミュレーション結果）可視化するためのＧＵＩによって構成されている。なお、この際、出力部１３は、成膜処理のシミュレーション結果と共に、例えば、演算に用いた成膜条件およびパラメータ等の情報も出力してもよい。出力部１３は、例えば、シミュレーション結果を表示する表示装置、シミュレーション結果を印刷して出力する印刷装置、シミュレーション結果を記録する記録装置等の装置のいずれか、または、これらの装置を適宜組み合わせて構成される。なお、本実施の形態では、シミュレータが出力部１３を備える例を説明するが、本技術はこれに限定されず、出力部１３がシミュレータの外部に設けられていてもよい。
【0022】
このシミュレータはさらに、演算部１２における計算処理に必要となる各種パラメータを記憶するデータベース部を備えていてもよい。また、このようなデータベース部を、シミュレータの外部に設けてもよい。なお、計算処理に必要となる各種パラメータを外部から随時入力する場合には、データベース部を設けなくてもよい。
【0023】
（成膜シミュレーション方法の詳細）
本シミュレーションでは、成膜処理の予測技術として、フラックス法ベースのＶｏｘｅｌモデルが用いられる。一般にＶｏｘｅｌモデルでは、計算領域内に配置するＶｏｘｅｌは立方体となっている。本実施の形態では、Ｖｏｘｅｌには、膜がそこに存在するか否かの存在情報だけでなく、膜の被覆性（カバレッジ）および膜質（例えば、密度、欠陥密度、透水性、および密着性など）といった情報も含まれる。さらに、本実施の形態では、Ｖｏｘｅｌに流入するガスフラックス（入射フラックス）は、直接入射してくるガス成分（直接入射成分）だけでなく、周囲のパターンから廻りこんでくるガス成分（周囲構造に応じたガス廻りこみ成分）も考慮して計算する。さらに、入射フラックスに応じた代表粒子という概念を用いて、膜の存在や、被覆性、さらには膜質を予測する。これにより、モンテカルロ法による流体計算と比べて、少ない計算コストで、膜の存在や、被覆性、さらには膜質を数十ｎｍの範囲にて精度良く予測することが可能である。
本シミュレーションでは、成膜処理の予測技術として、フラックス法ベースのＶｏｘｅｌモデルが用いられる。一般にＶｏｘｅｌモデルでは、計算領域内に配置するＶｏｘｅｌは立方体となっている。本実施の形態では、Ｖｏｘｅｌには、膜がそこに存在するか否かの存在情報だけでなく、膜の被覆性（カバレッジ）および膜質（例えば、密度、欠陥密度、透水性、および密着性など）といった情報も含まれる。さらに、本実施の形態では、Ｖｏｘｅｌに流入するガスフラックス（入射フラックス）は、直接入射してくるガス成分（直接入射成分）だけでなく、周囲のパターンから廻りこんでくるガス成分（周囲構造に応じたガス廻りこみ成分）も考慮して計算する。さらに、入射フラックスに応じた代表粒子という概念を用いて、膜の存在や、被覆性、さらには膜質を予測する。これにより、モンテカルロ法による流体計算と比べて、少ない計算コストで、膜の存在や、被覆性、さらには膜質を数十ｎｍの範囲にて精度良く予測することが可能である。
【0024】
【0025】
まず、演算部１２は、成膜の初期条件を設定する（ステップＳ１０１、図１，図２）。成膜の初期条件とは、具体的には、成膜条件に関する情報、および下地層に関する情報などを含む。例えば、成膜方法がガスを原料とする蒸着方法である場合、ステップＳ１０１では、成膜条件に関する情報として、表面反応に関するモデルパラメータ、入射ガスフラックス、イオンの入射エネルギーおよび角度分布等が設定されてもよい。また、下地層に関する情報として、下地層の材質および形状が設定されてもよい。
【0026】
次に、演算部１２は、表面Ｖｏｘｅｌを選択する（ステップＳ２０１、図２）。１つ前の時間ステップで堆積された堆積膜のＶｏｘｅｌに隣接する、空気領域（Ａｉｒ）のＶｏｘｅｌ（Ａｉｒ Ｖｏｘｅｌ）が表面Ｖｏｘｅｌとして定義され、この表面Ｖｏｘｅｌが入射ラジカルフラックス計算および表面反応計算（付着、脱離、マイグレーション、堆積）の対象となる。
【0027】
次に、演算部１２は、入射ラジカルフラックスを計算する（ステップＳ１０２、図１、図２）。通常の流体計算を逐次行うとすると、計算コストが大きく、半導体プロセスで取り扱うμｍオーダーのパターンへの適用は現実的ではない。そこで本実施の形態では、演算部１２は、表面Ｖｏｘｅｌへ直接入射してくるガス成分（直接入射成分）と周囲のパターンから廻りこんでくるガス成分（廻りこみ成分）の２つに分けて計算を行う。
【0028】
ここで、ガス輸送の概念図を図４に示す。パターンＡはパターンＢよりも間口の大きいパターンである。ＣＢは直接入射成分（密度）である。ＣＢは、表面Ｖｏｘｅｌからパターン間口をながめたときの立体角に依存する。ただし、隣接パターンの立体角寄与は含まない。ＣＡは廻りこみ成分（密度）である。演算部１２は、以下の式（１）～（３）を用いてＣＡを算出する。連続の式（式（１））とベルヌーイの定理（式（２），（３））から、パターンＡ(間口の大きいパターン)からパターンＢ（間口の小さいパターン）へ流入する接続面でのガス量ＣＡを求めることができる。
【0029】
【数1】
【数2】
【数3】
【0030】
ここで、ＳＡはパターンＡの開口面積である。ＳＢはパターンＡとパターンＢとの接続面の面積である。ＶＡはパターンＡにおけるガス熱速度、ＶＢはパターンＢにおけるガス熱速度、ＶはＶＡおよびＶＢの総称、ＫＢはボルツマン定数、ＴＮはガス温度、Ｍはガス質量、ＰＡは直接入射成分の圧力、ＰＢは廻りこみ成分の圧力である。
【0031】
さらに、演算部１２は、ＣＡの該当Ｖｏｘｅｌへの寄与成分Ｃ（ｚ）を下記の式（４）から算出する。接続面ＳＢから、パターンＢ内の黒丸にかけての領域を近似的に部分的なトレンチ（ドットパターンで表現された領域）ととらえ、下記の式（４）から、近似的にＣ（ｚ）を求めることができる。ここで、Ｗはトレンチ幅、Ｌは接続面から反対側のパターンＢの境界までの距離である。該当ＶｏｘｅｌでのトータルフラックスＦは、下記の式（５）から算出することができる。このトータルフラックスＦが、該当Ｖｏｘｅｌでの入射ラジカルフラックスに相当する。パターンＢの周囲に複数のパターン（以下、「周囲パターン」と称する。）が存在する場合には、各周囲パターンについて求めたＣ（ｚ）に対して１／Ｌを掛けることにより、Ｃ（ｚ）に重み付けを行い、得られた重み付けの値（Ｃ（ｚ）×１／Ｌ）を足し合わせることにより、該当ＶｏｘｅｌでのトータルフラックスＦを求めることができる。
【0032】
【数4】
【数5】
【0033】
次に、演算部１２は、入射ラジカルフラックス（＝トータルフラックスＦ）に応じた代表粒子を生成する（ステップＳ１０３、図１、図２）。図４に、実際の成膜表面（例として、ALD-SiO2のBDEASステップ後のＯ２ステップ）の概念図を示す。図６に、図５の概念図をＶｏｘｅｌモデルに落とし込んだ概念図を示す。
【0034】
実際の表面では、入射粒子１個１個が表面粒子との結合・ポテンシャル変化を引き起こしながら、付着、マイグレーションおよび堆積を繰り返し、膜が形成されていく。本アルゴリズムでは、計算負荷を大きく低減させる工夫として、複数の入射粒子を１つの代表粒子として代表させることで、～１０１８[／ｃｍ２／Ｓ]もあるラジカルフラックスの取り扱いを簡便にし、その代表粒子について、確率を用いて、付着、脱離、マイグレーションおよび堆積を取り扱う。これは、いわゆる統計的アンサンブル手法の応用である。
【0035】
演算部１２は、Ｖｏｘｅｌごとに代表粒子の生成個数Ｎを算出する。代表粒子の生成個数Ｎについては、堆積膜の仮の密度ρ[個／ｃｍ３]を設定し、体積Ｌ３×１０－２１[ｃｍ３]の１Ｖｏｘｅｌ（例えば、図７の太枠で囲んだＶｏｘｅｌ）内に入射するラジカル粒子の個数（Ｆ×Ｌ２×１０－１４×ｄｔ[個]）を堆積膜の１Ｖｏｘｅｌ内の粒子の個数（ρ×Ｌ３×１０－２１[個]）で割ることで算出する（式（６）参照）。代表粒子の生成個数Ｎは、パターン内での入射ラジカルフラックスＦ[個／ｃｍ２／ｓ]に応じて変動する。
【数6】
【0036】
例えば、Ｆ＝１０１８[個／ｃｍ２／ｓ]、ｄｔ＝０．１[ｓ]、Ｌ＝５[ｎｍ]、ρ＝１０２２[個／ｃｍ３]としたとき、代表粒子の生成個数Ｎは２０個となる。複数種類のラジカル粒子が存在する場合には、演算部１２は、この計算をラジカル粒子の種類ごとに行う。ＣＶＤまたはＰＶＤによって成膜を行う場合には、演算部１２は、同じ時間ステップ内でガスごとの代表粒子の生成個数Ｎを算出する。ＡＬＤによって成膜を行う場合には、演算部１２は、ガスごとに別々の時間ステップで代表粒子の生成個数Ｎを算出する。
【0037】
次に、演算部１２は、各々の代表粒子に対して、付着および脱離を決定する（ステップＳ１０４、図１、図２）。具体的には、各々の代表粒子に対して、付着確率Ｙ（０≦Ｙ≦１）で付着するか、確率（１－Ｙ）で脱離するかを、乱数を用いて決定する。代表粒子が下地層と接する場合（成膜のごく初期）には、加工による下地ダメージＤａによって付着確率Ｙｓが変動する影響を加味している（式（７））。一方、堆積膜上では付着確率Ｙｄは一定とする（式（８））。ただし、ａ，ｂ，ｃはユーザ設定の定数である。Ｄａには、他の成膜シミュレータで計算された結果ないしは実験値が与えられる。脱離と判定された場合は、それ以上の表面反応計算を行わない。
【数7】
【数8】
【0038】
次に、演算部１２は、マイグレーションの範囲を決定する（ステップＳ１０５、図１、図２）。図８に、マイグレーションおよび堆積時における実際の成膜表面の概念図を示した。図９に、マイグレーションおよび堆積時におけるＶｏｘｅｌモデルの概念図を示した。
【0039】
基板温度のエネルギーを使って付着したラジカル粒子は表面をマイグレーションする。表面ポテンシャルの安定した領域でラジカル粒子が堆積膜（成膜表面）と結合し、膜が堆積される。この時、一部のラジカル粒子において未結合手も存在する（図８）。この現象のＶｏｘｅｌモデル化に際して以下を考慮する。
【0040】
まず、基板温度Ｔと活性化エネルギーＥｄに依存して算出されるマイグレーション長ＬＤ（式（９））の２ＬＤ立方範囲（マイグレーションの範囲）を決定する（ステップＳ１０５、図１、図２）。続いて、２ＬＤ立方範囲内でＣｒａｄｌｅ（曲率半径の小さい凹）を探索する（図９）。ここで、Ｄ０は拡散定数、τは時定数、ＫＢはボルツマン定数である。式（９）の右辺の√Ｄ０τおよびＥｄは、第一原理計算ないしは基板温度と膜密度の実測の相関関係（傾きＥｄ、ｙ切片√Ｄ０τ）等から値を決定するパラメータである。
【数9】
【0041】
表面ポテンシャルμは曲率半径１／Ｒを用いて式（１０）のように示される。表面ポテンシャルμを最小にするように、すなわち、曲率半径１／Ｒの小さい構造（Ｃｒａｄｌｅ）がなくなるようにラジカル粒子はマイグレーションするものとする。ここで、μ０は平坦膜での化学ポテンシャル、γは表面張力、Ωは原子体積、１／Ｒは曲率半径である。
【数10】
【0042】
次に、演算部１２は、Ｖｏｘｅｌ空間内でのＣｒａｄｌｅを識別する（ステップＳ２０２、図２）。具体的には、図１０に示したように、太枠で囲ったマイグレーション範囲内（２ＬＤ）で、堆積膜に接するＡｉｒ Ｖｏｘｅｌ（重心座標は、ｉ，ｊ，ｋ）に着目し、着目したＡｉｒ Ｖｏｘｅｌに隣接する上下左右のＶｏｘｅｌの材質を判断する。着目したＡｉｒ Ｖｏｘｅｌに隣接するＶｏｘｅｌが堆積膜に該当する場合には、着目したＡｉｒ Ｖｏｘｅｌの判定指標：Ｎｖ（ｉ，ｊ，ｋ）に１を加算する。着目したＡｉｒ Ｖｏｘｅｌを、判定指標Ｎｖの値に応じて次のように判定する。
【0043】
・Ｎｖ≧３ →着目したＡｉｒ ＶｏｘｅｌはＣｒａｄｌｅである判定する。
・Ｎｖ＝１または２ →着目したＡｉｒ ＶｏｘｅｌはＣｒａｄｌｅではないと判定する。
フラグ：Ｆ（ｉ，ｊ，ｋ）＝１とする。
・Ｎｖ＝０ →着目したＡｉｒ ＶｏｘｅｌはＡｉｒのままであると判定する。
・Ｎｖ＝１または２ →着目したＡｉｒ ＶｏｘｅｌはＣｒａｄｌｅではないと判定する。
フラグ：Ｆ（ｉ，ｊ，ｋ）＝１とする。
・Ｎｖ＝０ →着目したＡｉｒ ＶｏｘｅｌはＡｉｒのままであると判定する。
【0044】
例えば、図１０の太枠で囲った２つのＡｉｒ Ｖｏｘｅｌのうち、一方のＡｉｒ Ｖｏｘｅｌでは、Ｎｖ＝２となり、他方のＡｉｒ Ｖｏｘｅｌでは、Ｎｖ＝３となったとする。このとき、Ｎｖ＝２のＡｉｒ ＶｏｘｅｌについてはＣｒａｄｌｅではないと判定される。Ｎｖ＝３のＡｉｒ ＶｏｘｅｌについてはＣｒａｄｌｅである判定される。
【0045】
次に、演算部１２は、代表粒子を移動させる（ステップＳ２０３、図２）。具体的には、マイグレーション範囲内にＮｖ≧３のＡｉｒ Ｖｏｘｅｌがあれば、代表粒子(ガス１、例えば、ALD-SiO2プロセスにおけるBDEAS)はそのＶｏｘｅｌへマイグレーションするとして、そのＶｏｘｅｌにおいて、結合判定のためのフラグ：ＦＢ（ｉ，ｊ，ｋ）を１とする。マイグレーション範囲内に複数のＣｒａｄｌｅがあった場合には、代表粒子から最も距離の近いＣｒａｄｌｅにおいてＦＢ（ｉ，ｊ，ｋ）を１とする（図１１）。
【0046】
一方、マイグレーション範囲内にＮｖ≧３のＡｉｒ Ｖｏｘｅｌがない場合は、乱数を用いて（確率によって）、Ｆ（ｉ，ｊ，ｋ）＝１のＶｏｘｅｌの中から、代表粒子の移動先Ｖｏｘｅｌを決定する。その移動先Ｖｏｘｅｌにおいて、結合判定のためのフラグ：ＦＢ（ｉ，ｊ，ｋ）を１とする。（図１２）。演算部１２は、このようにして、代表粒子の堆積位置を選択（決定）する（ステップＳ１０６、図１、図２）。
【0047】
演算部１２は、生成された全ての代表粒子について、堆積位置が確定するまで上記のステップＳ１０４～Ｓ１０６を実行する（ステップＳ２０４；Ｎ）。
【0048】
なお、複数種類のガス（例えば、ガス１，ガス２）を用いて成膜を行う場合には、演算部１２は、例えば、ガス１（たとえば、BDEAS）の代表粒子の堆積位置を選択（決定）した後、ガス２(たとえば、Ｏ２)の代表粒子の付着・脱離～移動、堆積位置の選択（決定）について実行する（ステップＳ２０５，Ｓ２０６）。つまり、ＡＬＤによって成膜を行う。その際、演算部１２は、ガス１で既にＦＢ（Ｉ，Ｊ，Ｋ）＝１となっているＶｏｘｅｌにおいては（ステップＳ２０７；Ｙ）、乱数を用いて、堆積時の状態が結合状態（確率ＹＢ）なのか、未結合状態（ダングリングボンド系: 確率（１－ＹＢ）)なのかを決定する。一方、ＦＢ（Ｉ，Ｊ，Ｋ）＝０のＶｏｘｅｌについては、演算部１２は、堆積時の状態が未結合状態であると判定する（ステップＳ２０８）。なお、演算部１２は、Ｖｏｘｅｌの状態を決定する際に、そのＶｏｘｅｌにおいて、イオンが入射する場合（つまり、入射イオンフラックスΓｉがある場合）には、入射するイオンのフラックスおよびエネルギーに基づいてイオンによる結合増強の影響（結合および未結合の判断）を判断する。
【0049】
具体的な計算方法を以下に示す。すなわち、移動先Ｖｏｘｅｌ内（Ｖｏｘｅｌサイズ：Ｌ）に時間ステップｄｔ内に入射するイオン数は、式（１１）で示される。この式（１１）と、イオンエネルギー分布（IEDF：Ion Energy Distribution Function）とを用いて、該Ｖｏｘｅｌへの総入射エネルギーＥｔｏｔを算出できる。一方、１Ｖｏｘｅｌ内に存在する結合数は堆積膜の仮密度ρを用いると、ρ×Ｌ３－１となるので、結合１つにつき、式（１２）に示したエネルギーが加わる。
【0050】
【数11】
【数12】
【0051】
このようにして求めたエネルギーが、堆積膜の原子同士の結合エネルギーＥＢ（ユーザー設定値）と比較して大きい場合には、代表粒子は確率ＹＢで結合状態として堆積される。また、このようにして求めたエネルギーが、堆積膜の原子同士の結合エネルギーＥＢ（ユーザー設定値）と比較して小さい場合には、確率ＹＢ’（ユーザー設定値）の確率として結合状態として堆積される。演算部１２は、このようにして、Ｖｏｘｅｌの膜状態（結合、未結合）を決定する（ステップＳ２０９）。つまり、代表粒子と成膜表面との結合および未結合のいずれかの状態を付与したＶｏｘｅｌとして堆積を計算する。
【0052】
演算部１２は、マイグレーション範囲内の全てのＶｏｘｅｌにおいて膜状態（結合、未結合）を決定するまで上記のステップＳ２０５～Ｓ２０９を実行する（ステップＳ２１０；Ｎ）。演算部１２は、マイグレーション範囲内の全てのＶｏｘｅｌにおいて膜状態（結合、未結合）を決定した場合には（ステップＳ２１０；Ｙ）、形状進展を行う。具体的には、膜状態が結合となっているＶｏｘｅｌを堆積膜として認識し、膜状態が未結合となっているＶｏｘｅｌを未結合膜（欠陥膜）として認識する。このようにして、時間ステップごとに、膜が積まれていき、形状が変化していく。
【0053】
次に、演算部１２は、膜質計算を行う（ステップＳ１１０、図１、図２）。膜質として、膜密度、透水性、密着性の分布の計算を行う。図１２に示した膜モフォロジーに対して、各Ｖｏｘｅｌに付与した結合状態（Ｂ）および未結合状態（ＵＢ）ならびにＶｏｉｄ領域の情報を用いて、太枠で示した設定範囲内での平均化した物理量として、膜密度、透水性、成膜表面との密着性を計算する。つまり、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積を計算することで成膜表面における膜の被覆性および膜質を表現する。設定範囲は、ユーザによって設定されるパラメータであり、たとえば、１００ｎｍとなっている。
【0054】
（膜密度）
演算部１２は、図１３の設定領域内において、各Ｖｏｘｅｌの密度重みＷを与え、これらの総和平均を局所的な膜密度ＤＬ（ｌ，ｍ，ｎ）として定義する（式（１３））。密度重みＷは、結合ＶｏｘｅｌにおいてＷＢとし、未結合ＶｏｘｅｌにおいてＷＵＢとし、Ｖｏｉｄ領域のＶｏｘｅｌにおいてＷＶとする。なお、ｌ，ｍ，ｎは設定領域の重心座標、Ｑは設定領域内のＶｏｘｅｌ数である。ただし、ＤＬ０は、別途、典型的なプロセス条件によって実測から求められるか、または、ＭＤ計算等から求められる。ＤＬ０は、ユーザによって入力されるパラメータである。
演算部１２は、図１３の設定領域内において、各Ｖｏｘｅｌの密度重みＷを与え、これらの総和平均を局所的な膜密度ＤＬ（ｌ，ｍ，ｎ）として定義する（式（１３））。密度重みＷは、結合ＶｏｘｅｌにおいてＷＢとし、未結合ＶｏｘｅｌにおいてＷＵＢとし、Ｖｏｉｄ領域のＶｏｘｅｌにおいてＷＶとする。なお、ｌ，ｍ，ｎは設定領域の重心座標、Ｑは設定領域内のＶｏｘｅｌ数である。ただし、ＤＬ０は、別途、典型的なプロセス条件によって実測から求められるか、または、ＭＤ計算等から求められる。ＤＬ０は、ユーザによって入力されるパラメータである。
【数13】
【0055】
（透水性）
演算部１２は、図１３の設定領域内におけるＶｏｉｄ領域の体積Ｖｖと、図１３の設定領域の体積Ｖｔｏｔとの比を透水性ＷＰ(の指標)として算出する（式（１４））。なお、ｌ，ｍ，ｎは、設定領域の重心座標である。
演算部１２は、図１３の設定領域内におけるＶｏｉｄ領域の体積Ｖｖと、図１３の設定領域の体積Ｖｔｏｔとの比を透水性ＷＰ(の指標)として算出する（式（１４））。なお、ｌ，ｍ，ｎは、設定領域の重心座標である。
【数14】
【0056】
（密着性）
演算部１２は、図１４に示したように、下地層に接するＶｏｘｅｌを対象とし、図１４の設定領域内における未結合Ｖｏｘｅｌ（ＵＢ）の体積ＶＵＢと、図１４の設定領域の体積Ｖｔｏｔとの比を密着性Ａｄ(の指標)として算出する（式（１５））。ｌ，ｍ，ｎは、設定領域の重心座標である。
演算部１２は、図１４に示したように、下地層に接するＶｏｘｅｌを対象とし、図１４の設定領域内における未結合Ｖｏｘｅｌ（ＵＢ）の体積ＶＵＢと、図１４の設定領域の体積Ｖｔｏｔとの比を密着性Ａｄ(の指標)として算出する（式（１５））。ｌ，ｍ，ｎは、設定領域の重心座標である。
【数15】
【0057】
演算部１２は、このようにして膜質計算を行った後、計算中の時間を時間ステップｄｔだけ進める（ステップＳ２１１）。そして、計算時間が所定の処理時間に到達していない場合には（ステップＳ２１２；Ｎ）、演算部１２は、引き続き、上述のステップＳ２０１以降の各ステップを実行する。一方、計算時間が所定の処理時間に到達した場合には（ステップＳ２１２；Ｙ）、演算部１２は、成膜プロセスを終了し、続いて、アニールプロセスによる膜質計算を行う（ステップＳ１１１）。
【0058】
（アニール計算）
ここでは、演算部１２は、アニールによる膜の熱変性、つまり、アニールを施した時の膜密度を計算する。アニールによる膜モフォロジーと密度の変化の概念図を図１５に示す。まず、基板温度Ｔに応じた拡散長Ｌａを計算し、図１５において太枠で示した２Ｌａ範囲内にあるＮａ個の、Ａｉｒを含むＶｏｘｅｌへ密度ρ（Ｉ，Ｊ、Ｋ）を等分配する。これにより、図１５中の黒点のあるＶｏｘｅｌから、図１５において太枠で示した２Ｌａ範囲内にある他のＶｏｘｅｌへ代表粒子を拡散させる（図１６）。なお、拡散長Ｌａの計算には、式（８）と同じ式を用いる。ただし、式（８）におけるパラメータ値は、マイグレーション長ＬＤの計算のときのパラメータ値とは異なる。
ここでは、演算部１２は、アニールによる膜の熱変性、つまり、アニールを施した時の膜密度を計算する。アニールによる膜モフォロジーと密度の変化の概念図を図１５に示す。まず、基板温度Ｔに応じた拡散長Ｌａを計算し、図１５において太枠で示した２Ｌａ範囲内にあるＮａ個の、Ａｉｒを含むＶｏｘｅｌへ密度ρ（Ｉ，Ｊ、Ｋ）を等分配する。これにより、図１５中の黒点のあるＶｏｘｅｌから、図１５において太枠で示した２Ｌａ範囲内にある他のＶｏｘｅｌへ代表粒子を拡散させる（図１６）。なお、拡散長Ｌａの計算には、式（８）と同じ式を用いる。ただし、式（８）におけるパラメータ値は、マイグレーション長ＬＤの計算のときのパラメータ値とは異なる。
【0059】
さらに、周囲のＶｏｘｅｌ（Ｉ’，Ｊ’、Ｋ’）からも同時に代表粒子が流入してくることを考慮して、アニール後の密度ρａを以下のように表現する。
【数16】
【数17】
【0060】
さらに、熱による堆積膜の結晶性の改善を反映させるために、基板温度Ｔに依存した補正項βを式（１８）のρａに掛ける。ここで、Ａは係数、Ｅａは活性化エネルギー、ＫＢはボルツマン定数である。
【数18】
【0061】
上記の計算後、ＵＢ（未結合状態）のＶｏｘｅｌについては、アニール前のＢ（結合状態）のＶｏｘｅｌの密度よりも大きくなっている場合にはＢ（結合状態）のＶｏｘｅｌとして状態を変更する。また、ＡｉｒだったＶｏｘｅｌに関しては、密度設定閾値（ユーザー設定）よりも大きくなった場合には、堆積膜Ｖｏｘｅｌとし、その結合状態についても同様の処理を行う。以上の計算を、下地Ｖｏｘｅｌを除く全てのＶｏｘｅｌに対して行う。
【0062】
（ブリスター計算）
演算部１２は、成膜プロセスが終わった直後、あるいは、アニール直後（アニールを施した時）のブリスター（膜剥がれ）計算として、堆積膜中のＵＢ（未結合状態）のＶｏｘｅｌに対して以下に述べるＶｏｘｅｌ移動の処理を行う。
演算部１２は、成膜プロセスが終わった直後、あるいは、アニール直後（アニールを施した時）のブリスター（膜剥がれ）計算として、堆積膜中のＵＢ（未結合状態）のＶｏｘｅｌに対して以下に述べるＶｏｘｅｌ移動の処理を行う。
【0063】
まず、ＵＢ内のブリスター要因のガス密度ｎＢを算出する。成膜計算時のガス１とガス２のフラックス計算時に要因ガスのフラックスＦＢも導出し（SiO2成膜であれば、水素(Ｈ)）、Ｖｏｘｅｌ内に記憶させておく。このＦＢとＶｏｘｅｌサイズＬ、成膜時間ｔ、さらに、超粒子の個数ＮＳを用いると、ガス密度ｎＢは、以下の式（１９）で表される。
【数19】
【0064】
式１９を用いて得られたガス密度ｎＢのうちの一部（ＲＢ≦1、入射イオンエネルギーの依存）が脱ガスとしてＶｏｘｅｌから出ていくとすると、ガス圧のエネルギーは、以下の式（２０）で表される。
【数20】
【0065】
ここで、Ｔは成膜時の基板温度である。ＰＢ（Ｌ３）と、図１７中の太枠をまたぐ太い実線で示した、Ｖｏｘｅｌ隣接面の結合エネルギーの総和ＥＢｔｏｔとを比較し、ＰＢ（Ｌ３）が総和ＥＢｔｏｔよりも大きい場合には、図１７中の太枠で囲まれた各Ｖｏｘｅｌを１つ上に移動させ、元の場所はＡｉｒとする（図１８）。さらに、ＰＢ（Ｌ３）と、図１９中の太枠（図１７中の太枠の外側に隣接するＶｏｘｅｌを囲む太枠）をまたぐ太い実線で示した、Ｖｏｘｅｌ隣接面の結合エネルギーの総和ＥＢｔｏｔとを比較し、ＰＢ（Ｌ３）が総和ＥＢｔｏｔよりも大きい場合には、図１９中の太枠で囲まれた各Ｖｏｘｅｌを１つ上に移動させ、元の場所はＡｉｒとする（図２０）。総和ＥＢｔｏｔがＰＢ（Ｌ３）を超えるまで、このようなプロセルを繰り返し実行する。複数のＵＢ起因のＡｉｒ領域が互いに重なり合う場合には、互いに重なり合う複数のＡｉｒ領域を含む領域について、上記と同様の大小判定を行う。
【0066】
なお、総和ＥＢｔｏｔについては、隣接Ｖｏｘｅｌの面密度σと平均の結合エネルギーＥＢを用いて、近似的に次のように示す。また、ＲＢについては、膜堆積中の結合形成にイオンが影響していることを反映させるために、イオンエネルギーＥに依存した補正項：ｅｘｐ（－Ｅ０／Ｅ）×２．７１８を掛けた値を用いる。ここで、Ｅ０は、基準となるイオンエネルギーである。
【数21】
【0067】
＜効果＞
本実施の形態に係る成膜シミュレーション方法および成膜シミュレータ１では、入射ラジカルフラックスに応じた代表粒子が生成され、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積の計算が実行される。これにより、精度と計算時間が粒子数に依存しなくなる。また、本実施の形態では、上記計算において、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積が計算されることで、成膜表面における膜の被覆性および膜質が表現される。これにより、膜質計算のために、モフォロジーとガスフラックスに依存したデータベースを別途準備する必要がない。従って、計算時間の短縮化および計算精度の向上を図ることができる。また、本アルゴリズムの活用までのＴＡＴ短縮が期待される。
本実施の形態に係る成膜シミュレーション方法および成膜シミュレータ１では、入射ラジカルフラックスに応じた代表粒子が生成され、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積の計算が実行される。これにより、精度と計算時間が粒子数に依存しなくなる。また、本実施の形態では、上記計算において、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積が計算されることで、成膜表面における膜の被覆性および膜質が表現される。これにより、膜質計算のために、モフォロジーとガスフラックスに依存したデータベースを別途準備する必要がない。従って、計算時間の短縮化および計算精度の向上を図ることができる。また、本アルゴリズムの活用までのＴＡＴ短縮が期待される。
【0068】
本実施の形態では、成膜表面におけるＶｏｘｅｌに対して直接入射してくるガス成分と、成膜表面におけるＶｏｘｅｌに対して周囲構造に応じて入射してくるガス回り込み成分とにより、入射ラジカルフラックスが算出される。これにより、実構造を反映した計算精度の向上を図ることができる。
【0069】
本実施の形態では、成膜表面におけるＶｏｘｅｌに対して入射するイオンのフラックスおよびエネルギーに基づいて、代表粒子と成膜表面との結合および未結合の判断が行われる。これにより、イオンによる結合増強の影響を考慮した成膜プロセスを実現することができる。
【0070】
本実施の形態では、結合および未結合のいずれかの状態情報を用いて、膜密度、透水性および成膜表面との密着性が計算される。これにより、より精度の高い成膜プロセスを実現することができる。
【0071】
本実施の形態では、複数種類のガスを用いて成膜を行う場合、同じ時間ステップ内でガスごとの代表粒子の生成個数が算出される。これにより、ＣＶＤまたはＰＶＤによる成膜プロセスを表現することができる。
【0072】
本実施の形態では、複数種類のガスを用いて成膜を行う場合、ガスごとに別々の時間ステップで代表粒子の生成個数が算出される。これにより、ＡＬＤによる成膜プロセスを表現することができる。
【0073】
本実施の形態では、成膜後に、アニールを施した時の膜密度およびブリスターの少なくとも一方が計算される。このように、本実施の形態では、成膜プロセスだけでなく、アニールによる膜の変化を表現することができる。
【0074】
本実施の形態に係る成膜シミュレータ１では、入力部１１が成膜条件を設定するためのＧＵＩもしくはＣＵＩによって構成され、出力部１３が演算部１２での演算結果を可視化するためのＧＵＩによって構成されている。これにより、ユーザは、成膜シミュレータ１を比較的簡単に操作することができ、また、成膜シミュレータ１による演算結果を比較的簡単に把握することができる。
【0075】
＜実施例１＞
～基本のＣＶＤ計算（カバレッジ＋膜密度）～
本実施例に係るＣＶＤプロセスによるカバレッジと膜密度分布の計算を図２１に示す。基板温度が８００℃、ガス圧力は４００Ｐａ、BDEASとO2のガスによる下地Ｓｉトレンチに対する成膜である。下地Ｓｉトレンチは、線幅１５０ｎｍ、深さ２μｍ、奥行き４５０ｎｍの２つのトレンチが交差した構造となっている。図２１には、トレンチエッチング時の生成された加工ダメージを反映した付着確率（０．５＋０．５×Ｄａ）と入射イオン（最大５０ｅＶのエネルギー）による堆積時のＶｏｘｅｌ結合状態の確率（イオン入射がない場合には０．２、ある場合には１）を加味したシミュレーション結果が示されている。活性化エネルギー（マイグレーション、イオン入射による結合状態の判定）は一律０．４ｅＶと設定した。Ｖｏｘｅｌサイズは５ｎｍ、膜質計算で用いる平均化の設定範囲は２５ｎｍである。また、パターンへの入射イオンフラックスは４×１０１６／ｃｍ２／ｓ、ラジカル（BDEAS、O2）のフラックスは１×１０１８／ｃｍ２／ｓである。
～基本のＣＶＤ計算（カバレッジ＋膜密度）～
本実施例に係るＣＶＤプロセスによるカバレッジと膜密度分布の計算を図２１に示す。基板温度が８００℃、ガス圧力は４００Ｐａ、BDEASとO2のガスによる下地Ｓｉトレンチに対する成膜である。下地Ｓｉトレンチは、線幅１５０ｎｍ、深さ２μｍ、奥行き４５０ｎｍの２つのトレンチが交差した構造となっている。図２１には、トレンチエッチング時の生成された加工ダメージを反映した付着確率（０．５＋０．５×Ｄａ）と入射イオン（最大５０ｅＶのエネルギー）による堆積時のＶｏｘｅｌ結合状態の確率（イオン入射がない場合には０．２、ある場合には１）を加味したシミュレーション結果が示されている。活性化エネルギー（マイグレーション、イオン入射による結合状態の判定）は一律０．４ｅＶと設定した。Ｖｏｘｅｌサイズは５ｎｍ、膜質計算で用いる平均化の設定範囲は２５ｎｍである。また、パターンへの入射イオンフラックスは４×１０１６／ｃｍ２／ｓ、ラジカル（BDEAS、O2）のフラックスは１×１０１８／ｃｍ２／ｓである。
【0076】
成膜時間が進むにつれて、トレンチのスリット領域が閉塞され、トレンチ内部にはSiO2が薄く堆積してボイドが形成されている。また、イオン照射の影響を強く受けている（Ｖｏｘｅｌが結合状態と判定される確率が高い）平坦部での膜密度が高く、より影響の低いトレンチ内部では膜密度は低くなっている。Ｖｏｘｅｌの結合状態/未結合状態の分布、ガスフラックス分布の結果を合わせて、図２２に示す。
【0077】
＜実施例２＞
～基本のＣＶＤ計算（透水性）～
本実施例に係る透水性分布の計算結果を図２３に示す。実施例１の成膜条件と下地構造を用いた時の計算結果である。平坦部では透水性は低く、トレンチ内部で透水性がより高くなっている。
～基本のＣＶＤ計算（透水性）～
本実施例に係る透水性分布の計算結果を図２３に示す。実施例１の成膜条件と下地構造を用いた時の計算結果である。平坦部では透水性は低く、トレンチ内部で透水性がより高くなっている。
【0078】
＜実施例３＞
～基本のＣＶＤ計算（密着性）～
本実施例に係る密着性分布の計算結果を図２４に示す。実施例１の成膜条件と下地構造を用いた時の計算結果である。平坦部では密着性は高く、トレンチ内部で密着性がより低くなっている。
～基本のＣＶＤ計算（密着性）～
本実施例に係る密着性分布の計算結果を図２４に示す。実施例１の成膜条件と下地構造を用いた時の計算結果である。平坦部では密着性は高く、トレンチ内部で密着性がより低くなっている。
【0079】
＜実施例４＞
～基本のＣＶＤ計算（アニール）～
本実施例に係るアニール計算の結果を図２５、図２６、図２７に示す。実施例１の成膜後に、１１００℃アニールを施した時の膜密度分布（図２５）、密着性分布（図２６）、透水性分布（図２７）の変化であり、これら膜質がトレンチ内部にかけて改善する計算結果となっている。
～基本のＣＶＤ計算（アニール）～
本実施例に係るアニール計算の結果を図２５、図２６、図２７に示す。実施例１の成膜後に、１１００℃アニールを施した時の膜密度分布（図２５）、密着性分布（図２６）、透水性分布（図２７）の変化であり、これら膜質がトレンチ内部にかけて改善する計算結果となっている。
【0080】
＜実施例５＞
～基本のＣＶＤ計算（ブリスター）～
本実施例に係るブリスター計算の概念図を図２８に示す。実施例１の成膜後に、ブリスター計算を行ったものである。基板温度は８００℃、要因ガスフラックスは１×１０１８／ｃｍ２／ｓ、ＲＢ＝０．０１、結合エネルギーは５ｅＶと設定した場合、平坦部にＶｏｘｅｌ２層分のブリスター（膜剥がれ）が形成される結果となった。
～基本のＣＶＤ計算（ブリスター）～
本実施例に係るブリスター計算の概念図を図２８に示す。実施例１の成膜後に、ブリスター計算を行ったものである。基板温度は８００℃、要因ガスフラックスは１×１０１８／ｃｍ２／ｓ、ＲＢ＝０．０１、結合エネルギーは５ｅＶと設定した場合、平坦部にＶｏｘｅｌ２層分のブリスター（膜剥がれ）が形成される結果となった。
【0081】
＜実施例６＞
～応用ＣＶＤ計算（線幅が広いパターン＋線幅が狭いパターン）
実施例１の計算では２つの同じ線幅のトレンチの交差構造を例にとった計算であった。一方、本アルゴリズムによる計算はそれに限定されるものではなく、線幅が広いパターンと線幅が狭いパターンとが互いに混ざった任意の構造にも応用することができる。
～応用ＣＶＤ計算（線幅が広いパターン＋線幅が狭いパターン）
実施例１の計算では２つの同じ線幅のトレンチの交差構造を例にとった計算であった。一方、本アルゴリズムによる計算はそれに限定されるものではなく、線幅が広いパターンと線幅が狭いパターンとが互いに混ざった任意の構造にも応用することができる。
【0082】
例えば、図２９（Ａ）に示したように、円形状のパターンＢと、その周囲に連結された４つの方形状のパターンＡとによって構成されたパターンについて、本アルゴリズムによる計算を行うことができる。また、例えば、図２９（Ｂ）に示したように、方形状のパターンＢと、その周囲に連結された４つの方形状のパターンＡとによって構成されたパターンについて、本アルゴリズムによる計算を行うことができる。
【0083】
また、例えば、図２９（Ｃ）に示したように、円形状のパターンＢと、その周囲であって、かつ、その円形状のパターンＢを介して互いに対向する位置に連結された２つの方形状のパターンＡとによって構成されたパターンについて、本アルゴリズムによる計算を行うことができる。また、例えば、図２９（Ｄ）に示したように、方形状のパターンＢと、その周囲であって、かつ、その方形状のパターンＢを介して互いに対向する位置に連結された２つの方形状のパターンＡとによって構成されたパターンについて、本アルゴリズムによる計算を行うことができる。
【0084】
また、例えば、図２９（Ｅ）に示したように、円形状のパターンＢと、その周囲であって、かつ、互いの延在方向が９０°で交差する位置に連結された２つの方形状のパターンＡとによって構成されたパターンについて、本アルゴリズムによる計算を行うことができる。また、例えば、図２９（Ｆ）に示したように、方形状のパターンＢと、その周囲であって、かつ、互いの延在方向が９０°で交差する位置に連結された２つの方形状のパターンＡとによって構成されたパターンについて、本アルゴリズムによる計算を行うことができる。
【0085】
＜実施例７＞
～ＡＬＤ計算～
実施例１～６ではＣＶＤの計算を示した。一方で、本実施例では、図２でのガス１とガス２に関する計算を成膜終了時間まで設定時間ごと（たとえば、１０ｓ）に交互に繰り返すことで、ＡＬＤプロセスとしてカバレッジと膜質を計算することができる。
～ＡＬＤ計算～
実施例１～６ではＣＶＤの計算を示した。一方で、本実施例では、図２でのガス１とガス２に関する計算を成膜終了時間まで設定時間ごと（たとえば、１０ｓ）に交互に繰り返すことで、ＡＬＤプロセスとしてカバレッジと膜質を計算することができる。
【0086】
＜実施例８＞
～ＰＶＤ計算～
実施例１～６ではＣＶＤの計算を示した。一方で、本実施例では、図２でのガス１のみの計算を行うことで単一元素から構成される金属膜ターゲットのＰＶＤ計算とすることができる。また、ガス１とガス２、さらにはガス３以降の計算を追加することで、複数元素からなる化合物ターゲットのＰＶＤ計算とすることができる。
～ＰＶＤ計算～
実施例１～６ではＣＶＤの計算を示した。一方で、本実施例では、図２でのガス１のみの計算を行うことで単一元素から構成される金属膜ターゲットのＰＶＤ計算とすることができる。また、ガス１とガス２、さらにはガス３以降の計算を追加することで、複数元素からなる化合物ターゲットのＰＶＤ計算とすることができる。
【0087】
＜実施例９＞
～気相計算との連携～
本実施例を図３０に示す。図１の入射フラックス計算に関して、気相計算を行うことでチャンバ面内でのガス１とガス２の密度分布をシミュレーションする。これら密度分布を入力として、ガスの熱速度を掛けあわせることでフラックスを導出する。気相計算では、上部電極に印加するパワーと周波数、ガス種、ガス流量、ガス圧力、チャンバ壁状態（ガスの付着確率）、チャンバ構成、排気速度をインプットとする。あるいは、実測データベースに基づく補間計算でもよいし、機械学習（たとえば、ガウス過程回帰やディープラーニング等。ただし、手法は限定しない）を用いた計算でもよい。
～気相計算との連携～
本実施例を図３０に示す。図１の入射フラックス計算に関して、気相計算を行うことでチャンバ面内でのガス１とガス２の密度分布をシミュレーションする。これら密度分布を入力として、ガスの熱速度を掛けあわせることでフラックスを導出する。気相計算では、上部電極に印加するパワーと周波数、ガス種、ガス流量、ガス圧力、チャンバ壁状態（ガスの付着確率）、チャンバ構成、排気速度をインプットとする。あるいは、実測データベースに基づく補間計算でもよいし、機械学習（たとえば、ガウス過程回帰やディープラーニング等。ただし、手法は限定しない）を用いた計算でもよい。
【0088】
【0089】
（成膜シミュレーションプログラム２の構成）
図３１に示した成膜シミュレーションプログラム２は、入力部２１と、演算エンジン部２２と、シミュレーション結果を可視化するための表示部２３とを備えている。演算エンジン部２２は、ガスフラックス計算部２２ａ、カバレッジ計算部２２ｂと、膜質計算部２２ｃと、出力部２２ｄとを有している。
図３１に示した成膜シミュレーションプログラム２は、入力部２１と、演算エンジン部２２と、シミュレーション結果を可視化するための表示部２３とを備えている。演算エンジン部２２は、ガスフラックス計算部２２ａ、カバレッジ計算部２２ｂと、膜質計算部２２ｃと、出力部２２ｄとを有している。
【0090】
入力部２１は、初期条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成されている。入力部２１は、設定された初期条件の値を演算エンジン部２２に出力する。表示部２３は、出力部２２ｄから得られたデータを可視化するためのＧＵＩによって構成されている。
【0091】
ガスフラックス計算部２２ａは、入力部２１で設定された初期条件の値を用いて、上述のステップＳ１０２を実行することにより、表面Ｖｏｘｅｌへ入射する入射ラジカルフラックスを計算する。カバレッジ計算部２２ｂは、ガスフラックス計算部２２ａで計算された入射ラジカルフラックスの値を用いて、上述のステップＳ１０３～Ｓ１０９を実行することにより、膜の形状や状態を計算する。膜質計算部２２ｃは、カバレッジ計算部２２ｂで得られた膜の形状や状態の下、上述のステップＳ１１０、Ｓ１１１を実行することにより、膜質を計算する。出力部２２ｄは、カバレッジ計算部２２ｂおよび膜質計算部２２ｃによって計算された所定の成膜処理のシミュレーション結果を出力する。
【0092】
この成膜シミュレーションプログラム２の実行プラットフォームは、例えば、Ｗｉｎｄｏｗｓ（登録商標）、Ｌｉｎｕｘ（登録商標）、Ｕｎｉｘ（登録商標）、またはＭａｃ（登録商標）のいずれでもよい。また、入力部２１および表示部２３に用いられるＧＵＩは、ＯｐｅｎＧＬ、Ｍｏｔｉｆ、ｔｃｌ／ｔｋなど構成言語を問わない。演算エンジン部２２のプログラミング言語は、Ｃ、Ｃ＋＋、Ｆｏｒｔｒａｎ、ＪＡＶＡ（登録商標）などその種類を問わない。この成膜シミュレーションプログラム２は、他のシミュレータによる計算結果を取り込み演算エンジン部２２に渡す機能を有していてもよい。
【0093】
初期条件としては、例えば、レシピ情報、計算用パラメータ、パターン構造情報が入力される。カバレッジ計算部２２ｂおよび膜質計算部２２ｃでの計算終了後には、出力部２２ｄからは、各Ｖｏｘｅｌの材質と座標、各Ｖｏｘｅｌでの結合状態・膜質（膜密度、透水性、密着性）・表面フラックスをファイルないしは一時記憶領域に出力される。ＧＵＩによってこれらの結果の可視化が行われる。データ出力や可視化は、計算中にリアルタイムに行われても構わない。
【0094】
本実施例に係る成膜シミュレーションプログラム２では、入射ラジカルフラックスに応じた代表粒子が生成され、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積の計算が実行される。これにより、精度と計算時間が粒子数に依存しなくなる。また、本実施例では、上記計算において、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積が計算されることで、成膜表面における膜の被覆性および膜質が表現される。これにより、膜質計算のために、モフォロジーとガスフラックスに依存したデータベースを別途準備する必要がない。従って、計算時間の短縮化および計算精度の向上を図ることができる。また、本アルゴリズムの活用までのＴＡＴ短縮が期待される。
【0095】
本実施例では、入力部２１が成膜条件を設定するためのＧＵＩもしくはＣＵＩによって構成され、出力部２２ｄがカバレッジ計算部２２ｂおよび膜質計算部２２ｃでの演算結果を可視化するためのＧＵＩによって構成されている。これにより、ユーザは、成膜シミュレーションプログラム２を比較的簡単に実行することができ、また、成膜シミュレーションプログラム２による演算結果を比較的簡単に把握することができる。
【0096】
【0097】
（成膜装置３の構成）
成膜装置３の概念図を図３２に示す。成膜装置３は、成膜チャンバ３１、成膜シミュレーションシステム３２、制御システム３３、およびＦＤＣ／ＥＥＳ（Fault Detection and Classification/Equipment Engineering System）システム３４を備えている。
成膜装置３の概念図を図３２に示す。成膜装置３は、成膜チャンバ３１、成膜シミュレーションシステム３２、制御システム３３、およびＦＤＣ／ＥＥＳ（Fault Detection and Classification/Equipment Engineering System）システム３４を備えている。
【0098】
成膜チャンバ３１は、チャンバ内の状態をモニタリングするモニタリング装置３１Ａを有している。モニタリング装置３１Ａは、例えば、OES(Optical Emission Spectroscopy)を有している。ＯＥＳは、チャンバ内のプラズマからの発光をモニタリングする計測装置である。成膜時には、プラズマ中に存在するガス種ごとに特有の波長の光が発せられる。ＯＥＳは、その光を測定する。ＯＥＳは、波長ごとの発光強度を測定する。ＯＥＳは、測定された波長からガス種を特定し、特定したガス種についての情報をモニタリングデータとして出力する。モニタリング装置３１Ａは、例えば、チャンバ内の状態をモニタリングするシステムを有している。このシステムは、ガス圧力、流量、温度、パワー、バイアス、マッチャー容量、真空ポンプの開度等の時間変動を測定し、測定により得られたデータをモニタリングデータとして出力する。
【0099】
成膜チャンバ３１は、モニタリング装置３１Ａで得られたデータ（モニタリングデータ）を成膜シミュレーションシステム３２に送信する。成膜シミュレーションシステム３２は、成膜シミュレーションプログラム２を実行する成膜シミュレータ１を含んで構成されている。成膜シミュレーションシステム３２は、ガスフラックス演算部３２Ａ、最適化演算部３２Ｂおよび補正条件出力部３２Ｃを有している。ガスフラックス演算部３２Ａ、最適化演算部３２Ｂおよび補正条件出力部３２Ｃは、集積回路によって構成されていてもよいし、成膜シミュレーションプログラム２がロードされた演算装置によって構成されていてもよい。
【0100】
ガスフラックス演算部３２Ａは、モニタリング装置３１Ａから入力されたモニタリングデータを用いて、入射ラジカルフラックスを計算する。最適化演算部３２Ｂは、ガスフラックス演算部３２Ａで得られた入射ラジカルフラックスの値を用いて、上述のステップＳ１０３～Ｓ１０９を実行することにより、入射ラジカルフラックスの値に応じた代表粒子を生成し、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積を計算する。その結果、最適化演算部３２Ｂは、膜の形状や状態を予測することができる。最適化演算部３２Ｂは、さらに、得られた膜の形状や状態の下、上述のステップＳ１１０、Ｓ１１１を実行することにより、膜質を計算する。
【0101】
最適化演算部３２Ｂは、このようにして得られた計算値（予測値）と、所望のスペック（例えば、設定された膜厚情報Ｄ２）とを比較し、計算値（予測値）が許容範囲外となっている場合には、計算値（予測値）が許容範囲となるような成膜プロセスの最適化条件を探索する。具体的には、最適化演算部３２Ｂは、ガス流量、ガス圧力、基板温度、処理時間を変化させて、所定のアルゴリズムから解を見つけ出す。最適化演算部３２Ｂは、ウェハ毎、または、ロット毎に、最適化条件を探索する。
【0102】
その結果、最適化条件が見つからなかった場合には、最適化演算部３２Ｂは、ＦＤＣ／ＥＥＳシステム３４に異常信号を送信する。ＦＤＣ／ＥＥＳシステム３４は、最適化演算部３２Ｂから異常信号を受信したときは、制御システム３３に対して、成膜チャンバ３１の動作を停止させる信号を出力する。制御システム３３は、成膜チャンバ３１の動作を停止させる信号を受信すると、成膜チャンバ３１の動作を停止させる。
【0103】
一方で、最適化条件が見つかった場合には、最適化演算部３２Ｂは、見つかった最適化条件を補正条件出力部３２Ｃに出力する。補正条件出力部３２Ｃは、最適化条件が入力されると、成膜チャンバ３１の成膜条件が最適化条件となるのに必要なデータを生成し、制御システム３３に出力する。制御システム３３は、最適化条件とするのに必要なデータを受信すると、受信したデータに基づいて成膜チャンバ３１の動作を制御する。制御システム３３は、受信したデータに基づいてレシピ情報Ｄ１を修正する。
【0104】
最適化演算部３２Ｂに関しては、計算時間が実加工時間と同等のスケール以上の場合には、上記のようなオンラインで最適解を見出すのではなく、様々なプロセス条件に対してあらかじめ本アルゴリズムでシミュレーションしたデータベースを作成しておき、そのデータベースを利用して予測した入射ラジカルフラックスの値を用いて、上述のステップＳ１０３～Ｓ１０９を実行することにより、入射ラジカルフラックスの値に応じた代表粒子を生成する方法（オフラインによる方法）でも構わない。あるいは、データベースを使って機械学習モデルを構築し、そのモデルをオンラインの最適化演算部３２Ｂとしても構わない。
【0105】
本実施例に係る成膜装置３では、入射ラジカルフラックスに応じた代表粒子が生成され、確率によって成膜表面での各代表粒子の付着、脱離、マイグレーションおよび堆積の計算が実行される。これにより、精度と計算時間が粒子数に依存しなくなる。また、本実施例では、上記計算において、代表粒子と成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして堆積が計算されることで、成膜表面における膜の被覆性および膜質が表現される。これにより、膜質計算のために、モフォロジーとガスフラックスに依存したデータベースを別途準備する必要がない。従って、計算時間の短縮化および計算精度の向上を図ることができる。また、本アルゴリズムの活用までのＴＡＴ短縮が期待される。
【0106】
また、本実施例に係る成膜装置３では、上記計算によって得られた計算結果に基づいて、成膜プロセスの最適化条件が探索され、成膜チャンバの成膜条件が、見つかった最適化条件となるのに必要なデータが生成され、制御システム３３へ出力される。このようなフィードバックによって、より質の高い半導体デバイスを製造することが可能となる。
【0107】
本実施例では、成膜チャンバ内の状態がモニタリング装置３１Ａによってモニタリングされ、モニタリング装置３１Ａにより得られたモニタリングデータを用いて入射ラジカルフラックスが計算される。このような実測値を用いた計算を実行することにより、より質の高い半導体デバイスを製造することが可能となる。
【0108】
本実施例では、データベースを利用して予測した入射ラジカルフラックスの値に応じた代表粒子が生成される。これにより、計算時間が実加工時間と同等のスケール以上の場合であっても、より質の高い半導体デバイスを製造することが可能となる。
【0109】
本実施例では、ウェハ毎、または、ロット毎に、最適化条件が探索される。これにより、ウェハ毎、または、ロット毎に、より質の高い半導体デバイスを製造することが可能となる。
【0110】
以上、実施の形態およびその変形例を挙げて本開示を説明したが、本開示は上記実施の形態等に限定されるものではなく、種々変形が可能である。なお、本明細書中に記載された効果は、あくまで例示である。本開示の効果は、本明細書中に記載された効果に限定されるものではない。本開示が、本明細書中に記載された効果以外の効果を持っていてもよい。
【0111】
また、例えば、本開示は以下のような構成を取ることができる。
（１）
入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算する第１ステップを含み、
前記第１ステップにおいて、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現することを含む
成膜シミュレーション方法。
（２）
前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して直接入射してくるガス成分と、前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して周囲構造に応じて入射してくるガス回り込み成分とにより、前記入射ラジカルフラックスを算出する第２ステップを更に含む
（１）に記載の成膜シミュレーション方法。
（３）
前記第１ステップにおいて、前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して入射するイオンのフラックスおよびエネルギーに基づいて、前記代表粒子と前記成膜表面との結合および未結合の判断を行うことを更に含む
（１）または（２）に記載の成膜シミュレーション方法。
（４）
前記第１ステップにおいて、前記状態情報を用いて、膜密度、透水性および前記成膜表面との密着性を計算することを更に含む
（１）ないし（３）のいずれか１つに記載の成膜シミュレーション方法。
（５）
前記第１ステップにおいて、複数種類のガスを用いて成膜を行う場合、同じ時間ステップ内でガスごとの前記代表粒子の生成個数を算出することを更に含む
（１）ないし（３）のいずれか１つに記載の成膜シミュレーション方法。
（６）
前記第１ステップにおいて、複数種類のガスを用いて成膜を行う場合、ガスごとに別々の時間ステップで前記代表粒子の生成個数を算出することを更に含む
（１）ないし（３）のいずれか１つに記載の成膜シミュレーション方法。
（７）
前記第１ステップにおいて、成膜後に、アニールを施した時の膜密度およびブリスターの少なくとも一方を計算することを更に含む
（１）ないし（６）のいずれか１つに記載の成膜シミュレーション方法。
（８）
成膜条件を取得する入力部と、
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレーションプログラム。
（９）
前記入力部は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成され、
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
（８）に記載の成膜シミュレーションプログラム。
（１０）
成膜条件を取得する入力部と、
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレータ。
（１１）
前記入力部は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成され、
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
（１０）に記載の成膜シミュレータ。
（１２）
成膜チャンバと、
前記成膜チャンバの動作を制御する制御部と、
入射ラジカルフラックスの値に応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算し、それにより得られた計算結果に基づいて、成膜プロセスの最適化条件を探索する最適化演算部と、
前記成膜チャンバの成膜条件が、前記最適化演算部で見つかった前記最適化条件となるのに必要なデータを生成し、前記制御部へ出力する出力部と
を備え、
前記最適化演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜装置。
（１３）
前記成膜チャンバ内の状態をモニタリングするモニタリング装置と、
前記モニタリング装置により得られたモニタリングデータを用いて前記入射ラジカルフラックスを計算するガスフラックス演算部と
を更に備えた
（１２）に記載の成膜装置。
（１４）
前記最適化演算部は、データベースを利用して予測した前記入射ラジカルフラックスの値に応じた代表粒子を生成する
（１２）に記載の成膜装置。
（１５）
前記最適化演算部は、ウェハ毎、または、ロット毎に、前記最適化条件を探索する
（１２）ないし（１４）のいずれか１つに記載の成膜装置。
（１）
入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算する第１ステップを含み、
前記第１ステップにおいて、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現することを含む
成膜シミュレーション方法。
（２）
前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して直接入射してくるガス成分と、前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して周囲構造に応じて入射してくるガス回り込み成分とにより、前記入射ラジカルフラックスを算出する第２ステップを更に含む
（１）に記載の成膜シミュレーション方法。
（３）
前記第１ステップにおいて、前記成膜表面に隣接する、空気領域のＶｏｘｅｌに対して入射するイオンのフラックスおよびエネルギーに基づいて、前記代表粒子と前記成膜表面との結合および未結合の判断を行うことを更に含む
（１）または（２）に記載の成膜シミュレーション方法。
（４）
前記第１ステップにおいて、前記状態情報を用いて、膜密度、透水性および前記成膜表面との密着性を計算することを更に含む
（１）ないし（３）のいずれか１つに記載の成膜シミュレーション方法。
（５）
前記第１ステップにおいて、複数種類のガスを用いて成膜を行う場合、同じ時間ステップ内でガスごとの前記代表粒子の生成個数を算出することを更に含む
（１）ないし（３）のいずれか１つに記載の成膜シミュレーション方法。
（６）
前記第１ステップにおいて、複数種類のガスを用いて成膜を行う場合、ガスごとに別々の時間ステップで前記代表粒子の生成個数を算出することを更に含む
（１）ないし（３）のいずれか１つに記載の成膜シミュレーション方法。
（７）
前記第１ステップにおいて、成膜後に、アニールを施した時の膜密度およびブリスターの少なくとも一方を計算することを更に含む
（１）ないし（６）のいずれか１つに記載の成膜シミュレーション方法。
（８）
成膜条件を取得する入力部と、
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレーションプログラム。
（９）
前記入力部は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成され、
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
（８）に記載の成膜シミュレーションプログラム。
（１０）
成膜条件を取得する入力部と、
前記入力部で取得した前記成膜条件に基づいて、入射ラジカルフラックスに応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積計算する演算部と、
前記演算部での演算結果を出力する出力部と
を備え、
前記演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜シミュレータ。
（１１）
前記入力部は、成膜条件を設定するためのＧＵＩ（Graphical User Interface）もしくはＣＵＩ(Character-based User Interface)によって構成され、
前記出力部は、前記演算部での演算結果を可視化するためのＧＵＩによって構成されている
（１０）に記載の成膜シミュレータ。
（１２）
成膜チャンバと、
前記成膜チャンバの動作を制御する制御部と、
入射ラジカルフラックスの値に応じた代表粒子を生成し、確率によって成膜表面での各前記代表粒子の付着、脱離、マイグレーションおよび堆積を計算し、それにより得られた計算結果に基づいて、成膜プロセスの最適化条件を探索する最適化演算部と、
前記成膜チャンバの成膜条件が、前記最適化演算部で見つかった前記最適化条件となるのに必要なデータを生成し、前記制御部へ出力する出力部と
を備え、
前記最適化演算部は、前記代表粒子と前記成膜表面との結合および未結合のいずれかの状態情報を付与したＶｏｘｅｌとして前記堆積を計算することで前記成膜表面における膜の被覆性および膜質を表現する
成膜装置。
（１３）
前記成膜チャンバ内の状態をモニタリングするモニタリング装置と、
前記モニタリング装置により得られたモニタリングデータを用いて前記入射ラジカルフラックスを計算するガスフラックス演算部と
を更に備えた
（１２）に記載の成膜装置。
（１４）
前記最適化演算部は、データベースを利用して予測した前記入射ラジカルフラックスの値に応じた代表粒子を生成する
（１２）に記載の成膜装置。
（１５）
前記最適化演算部は、ウェハ毎、または、ロット毎に、前記最適化条件を探索する
（１２）ないし（１４）のいずれか１つに記載の成膜装置。
【符号の説明】
【0112】
１…成膜シミュレータ、２…成膜シミュレーションプログラム、３…成膜装置、１１…入力部、１２…演算部、１３…出力部、２１…入力部、２２…演算エンジン部、２２ａ…ガスフラックス計算部、２２ｂ…カバレッジ計算部、２２ｃ…膜質計算部、２２ｄ…出力部、２３…表示部、３１…成膜チャンバ、３１Ａ…モニタリング装置、３２…成膜シミュレーションシステム、３２Ａ…ガスフラックス演算部、３２Ｂ…最適化演算部、３２Ｃ…補正条件出力部、３３…制御システム、３４…ＦＤＣ／ＥＥＳシステム、Ｄ１…レシピ情報、Ｄ２…膜厚情報。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
【図31】
【図32】