(57)【要約】【課題】硬化性組成物の膜を形成する処理における硬化性組成物の挙動をより短時間で計測するために有利な技術を提供する。【解決手段】第１部材の上に配置された硬化性組成物の複数の液滴と第２部材とを接触させ、前記第１部材と前記第２部材との間の空間に前記硬化性組成物の膜を形成する処理における前記硬化性組成物の挙動を予測するシミュレーション方法であって、前記硬化性組成物の複数の液滴のそれぞれについて、前記液滴の輪郭を、前記液滴の代表点、前記代表点から前記輪郭上の点に向かう方向、及び、前記代表点から前記輪郭上の点までの距離で定義し、前記硬化性組成物の複数の液滴のそれぞれについて、前記輪郭の内側の領域の面積を、前記第１部材と前記第２部材との間の距離の変化に応じて、前記液滴の体積及び前記第１部材と前記第２部材との間の距離から得られる前記液滴の面積に一致させるように前記代表点から前記輪郭上の点までの距離を求める、ことを特徴とするシミュレーション方法を提供する。【選択図】図３
2020205416.tif 000002
【特許請求の範囲】
【請求項１】
第１部材の上に配置された硬化性組成物の複数の液滴と第２部材とを接触させ、前記第１部材と前記第２部材との間の空間に前記硬化性組成物の膜を形成する処理における前記硬化性組成物の挙動を予測するシミュレーション方法であって、
前記硬化性組成物の複数の液滴のそれぞれについて、前記液滴の輪郭を、前記液滴の代表点、前記代表点から前記輪郭上の点に向かう方向、及び、前記代表点から前記輪郭上の点までの距離で定義し、
前記硬化性組成物の複数の液滴のそれぞれについて、前記輪郭の内側の領域の面積を、前記第１部材と前記第２部材との間の距離の変化に応じて、前記液滴の体積及び前記第１部材と前記第２部材との間の距離から得られる前記液滴の面積に一致させるように前記代表点から前記輪郭上の点までの距離を求める、
ことを特徴とするシミュレーション方法。
【請求項２】
前記第１部材と前記第２部材との間の距離の変化に応じて、前記硬化性組成物の複数の液滴のうち隣接する液滴が接合したかどうかを判定する、
ことを特徴とする請求項１に記載のシミュレーション方法。
【請求項３】
前記隣接する液滴のうちの第１液滴の輪郭上の点が前記隣接する液滴のうちの第２液滴の輪郭の内側に位置しているかを判定し、
前記第１液滴の輪郭上の点が前記第２液滴の輪郭の内側に位置している場合には、前記第１液滴と前記第２液滴とが接合したと判定し、
前記第１液滴の輪郭上の点が前記第２液滴の輪郭の内側に位置していない場合には、前記第１液滴と前記第２液滴とが接合していないと判定する、
ことを特徴とする請求項２に記載のシミュレーション方法。
【請求項４】
前記硬化性組成物の複数の液滴のうち隣接する液滴と接合した液滴について、当該液滴の接合した部分に対応する、前記代表点から前記輪郭上の点までの距離を固定する、
ことを特徴とする請求項２又は３に記載のシミュレーション方法。
【請求項５】
前記液滴の輪郭を定義するために設定される前記液滴の輪郭上の有限個の点を除く前記液滴の輪郭上の点を、前記代表点から前記輪郭上の前記有限個の点までの距離を用いて補間する、
ことを特徴とする請求項１乃至４のうちいずれか１項に記載のシミュレーション方法。
【請求項６】
前記液滴の輪郭を定義するために設定される前記液滴の輪郭上の有限個の点を除く前記液滴の輪郭上の点を、前記有限個の点を線でつなぐことで求める、
ことを特徴とする請求項１乃至４のうちいずれか１項に記載のシミュレーション方法。
【請求項７】
前記液滴の輪郭を定義するために設定される前記液滴の輪郭上の有限個の点は、前記第１部材と前記第２部材との間の距離の変化に応じて、隣接する液滴に接合するまで一定の方向に移動する、
ことを特徴とする請求項１乃至６のうちいずれか１項に記載のシミュレーション方法。
【請求項８】
前記代表点は、当該代表点に対応する液滴の内部に含まれる、
ことを特徴とする請求項１乃至７のうちいずれか１項に記載のシミュレーション方法。
【請求項９】
第１部材の上に配置された硬化性組成物の複数の液滴と第２部材とを接触させ、前記第１部材と前記第２部材との間の空間に前記硬化性組成物の膜を形成する処理における前記硬化性組成物の挙動を予測するシミュレーション装置であって、
前記硬化性組成物の複数の液滴のそれぞれについて、前記液滴の輪郭を、前記液滴の代表点、前記代表点から前記輪郭上の点に向かう方向、及び、前記代表点から前記輪郭上の点までの距離で定義し、
前記硬化性組成物の複数の液滴のそれぞれについて、前記輪郭の内側の領域の面積を、前記第１部材と前記第２部材との間の距離の変化に応じて、前記液滴の体積及び前記第１部材と前記第２部材との間の距離から得られる前記液滴の面積に一致させるように前記代表点から前記輪郭上の点までの距離を求める、
ことを特徴とするシミュレーション装置。
【請求項１０】
請求項１乃至８のうちいずれか１項に記載のシミュレーション方法をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、シミュレーション方法、シミュレーション装置及びプログラムに関する。
【背景技術】
【０００２】
基板上に硬化性組成物を配置し、硬化性組成物を型と接触させ、硬化性組成物を硬化させることによって基板上に硬化性組成物の硬化物からなる膜を形成する膜形成技術がある。このような膜形成技術は、インプリント技術や平坦化技術に適用される。インプリント技術では、パターンを有する型を用いて、基板上の硬化性組成物と型のパターンとを接触させて硬化性組成物を硬化させることによって基板上の硬化性組成物に型のパターンが転写される。平坦化技術では、平坦面を有する型を用いて、基板上の硬化性組成物と平坦面とを接触させて硬化性組成物を硬化させることによって平坦な上面を有する膜が形成される。
【０００３】
基板上には、硬化性組成物が液滴の状態で配置され、その後、硬化性組成物の液滴に型が押し付けられる。これにより、基板上の硬化性組成物の液滴が広がって硬化性組成物の膜が形成される。この際、厚さが均一な硬化性組成物の膜を形成することや膜に気泡が残存しないことなどが重要であり、これを実現するために、硬化性組成物の液滴の配置や硬化性組成物への型の押し付けの方法及び条件などが調整される。このような調整を、装置を用いた試行錯誤によって実現するためには、膨大な時間と費用とを必要とする。そこで、このような調整を支援するシミュレータの開発が望まれている。
【０００４】
特許文献１には、パターン形成面に配置された複数の液滴の濡れ広がり及び合一（液滴の接合）を予測するためのシミュレーション方法が開示されている。かかるシミュレーション方法では、パターン形成面をモデル化した解析面が複数の解析セルに分割され、また、液滴は、解析面上のドロップサイトごとに配置される。特許文献１では、ドロップサイトは、ｍ×ｎの格子状に分割した領域であると定義され、解析セルとは別個の概念であると説明されている。
【０００５】
通常、液滴の挙動を計算する場合、液滴の寸法（大きさ）よりも十分に小さい計算要素（解析セル）を定義する必要がある。しかしながら、このような小さい計算要素を定義しながら、例えば、１つのショット領域などの広い領域の全域にわたって液滴の挙動を計算することは、極めて現実性に乏しく、許容可能な時間内に計算結果を得ることができない可能性がある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第５５９９３５６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明は、このような従来技術の課題に鑑みてなされ、硬化性組成物の膜を形成する処理における硬化性組成物の挙動をより短時間で計測するために有利な技術を提供することを例示的目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、本発明の一側面としてのシミュレーション方法は、第１部材の上に配置された硬化性組成物の複数の液滴と第２部材とを接触させ、前記第１部材と前記第２部材との間の空間に前記硬化性組成物の膜を形成する処理における前記硬化性組成物の挙動を予測するシミュレーション方法であって、前記硬化性組成物の複数の液滴のそれぞれについて、前記液滴の輪郭を、前記液滴の代表点、前記代表点から前記輪郭上の点に向かう方向、及び、前記代表点から前記輪郭上の点までの距離で定義し、前記硬化性組成物の複数の液滴のそれぞれについて、前記輪郭の内側の領域の面積を、前記第１部材と前記第２部材との間の距離の変化に応じて、前記液滴の体積及び前記第１部材と前記第２部材との間の距離から得られる前記液滴の面積に一致させるように前記代表点から前記輪郭上の点までの距離を求める、ことを特徴とする。
【０００９】
本発明の更なる目的又はその他の側面は、以下、添付図面を参照して説明される実施形態によって明らかにされるであろう。
【発明の効果】
【００１０】
本発明によれば、例えば、硬化性組成物の膜を形成する処理における硬化性組成物の挙動をより短時間で計測するために有利な技術を提供することができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の実施形態における膜形成装置及びシミュレーション装置の構成を示す概略図である。
【図２】一般的な手法における計算格子を例示する図である。
【図３】本発明の一側面としてのシミュレーション方法を説明するためのフローチャートである。
【図４】硬化性組成物の液滴要素の概念を示す図である。
【図５】２０個の角度で定義された硬化性組成物の液滴要素を示す図である。
【図６】一般的な手法における計算格子を例示する図である。
【図７】型の運動による液滴要素の面積の変化を模式的に示す図である。
【図８】液滴要素の広がりの一例を示す図である。
【図９】液滴要素の広がりの一例を示す図である。
【図１０】隣接する液滴要素が接合しているかどうかを判定する処理を説明するための図である。
【図１１】複数の液滴要素の接合によって気泡の閉じ込めが発生した状態を示す図である。
【図１２】図１に示すシミュレーション装置によって硬化性組成物の液滴の挙動を計算した一例を示す図である。
【発明を実施するための形態】
【００１２】
以下、添付図面を参照して実施形態を詳しく説明する。なお、以下の実施形態は特許請求の範囲に係る発明を限定するものではない。実施形態には複数の特徴が記載されているが、これらの複数の特徴の全てが発明に必須のものとは限らず、また、複数の特徴は任意に組み合わせられてもよい。更に、添付図面においては、同一もしくは同様の構成に同一の参照番号を付し、重複した説明は省略する。
【００１３】
図１は、本発明の実施形態における膜形成装置ＩＭＰ及びシミュレーション装置１の構成を示す概略図である。膜形成装置ＩＭＰは、基板Ｓの上に配置された硬化性組成物ＩＭの複数の液滴と型Ｍとを接触させ、基板Ｓと型Ｍとの間の空間に硬化性組成物ＩＭの膜を形成する処理を実行する。膜形成装置ＩＭＰは、例えば、インプリント装置として構成されてもよいし、平坦化装置として構成されてもよい。ここで、基板Ｓと型Ｍとは相互に入れ替え可能であり、型Ｍの上に配置された硬化性組成物ＩＭの複数の液滴と基板Ｓとを接触させ、型Ｍと基板Ｓとの間の空間に硬化性組成物ＩＭの膜が形成されてもよい。従って、包括的には、膜形成装置ＩＭＰは、第１部材の上に配置された硬化性組成物ＩＭの複数の液滴と第２部材とを接触させ、第１部材と第２部材との間の空間に硬化性組成物ＩＭの膜を形成する処理を実行する装置である。本実施形態では、第１部材を基板Ｓとし、第２部材を型Ｍとして説明するが、第１部材を型Ｍとし、第２部材を基板Ｓとしてもよい。この場合、以下の説明における基板Ｓと型Ｍとを相互に入れ替えればよい。
【００１４】
インプリント装置では、パターンを有する型Ｍを用いて、基板Ｓの上の硬化性組成物ＩＭに型Ｍのパターンが転写される。インプリント装置では、パターンが設けられたパターン領域ＰＲを有する型Ｍが用いられる。インプリント装置では、基板Ｓの上の硬化性組成物ＩＭと型Ｍのパターン領域ＰＲとを接触させ、基板Ｓのパターンを形成すべき領域と型Ｍとの間の空間に硬化性組成物ＩＭを充填させ、その後、硬化性組成物ＩＭを硬化させる。これにより、基板Ｓの上の硬化性組成物ＩＭに型Ｍのパターン領域ＰＲのパターンが転写される。インプリント装置では、例えば、基板Ｓの複数のショット領域のそれぞれに硬化性組成物ＩＭの硬化物からなるパターンが形成される。
【００１５】
平坦化装置では、平坦面を有する型Ｍを用いて、基板Ｓの上の硬化性組成物ＩＭと型Ｍの平坦面とを接触させ、硬化性組成物ＩＭを硬化させることによって、平坦な上面を有する膜が形成される。平坦化装置では、基板Ｓの全域をカバーする寸法（大きさ）を有する型Ｍが用いられる場合、基板Ｓの全域に硬化性組成物ＩＭの硬化物からなる膜が形成される。
【００１６】
硬化性組成物としては、硬化用のエネルギーが与えられることにより硬化する材料が使用される。硬化用のエネルギーとしては、電磁波や熱などが用いられる。電磁波は、例えば、その波長が１０ｎｍ以上１ｍｍ以下の範囲から選択される光、具体的には、赤外線、可視光線、紫外線などを含む。このように、硬化性組成物は、光の照射、或いは、加熱により硬化する組成物である。光の照射により硬化する光硬化性組成物は、少なくとも重合性化合物と光重合開始剤とを含有し、必要に応じて、非重合性化合物又は溶剤を更に含有してもよい。非重合性化合物は、増感剤、水素供与体、内添型離型剤、界面活性剤、酸化防止剤、ポリマー成分などの群から選択される少なくとも一種である。硬化性組成物の粘度（２５℃における粘度）は、例えば、１ｍＰａ・ｓ以上１００ｍＰａ・ｓ以下である。
【００１７】
基板の材料としては、例えば、ガラス、セラミックス、金属、半導体、樹脂などが用いられる。必要に応じて、基板の表面に、基板とは別の材料からなる部材が設けられてもよい。基板は、例えば、シリコンウエハ、化合物半導体ウエハ、石英ガラスを含む。
【００１８】
本明細書及び添付図面では、基板Ｓの表面に平行な方向をＸＹ平面とするＸＹＺ座標系で方向を示す。ＸＹＺ座標系におけるＸ軸、Ｙ軸及びＺ軸のそれぞれに平行な方向をＸ方向、Ｙ方向及びＺ方向とし、Ｘ軸周りの回転、Ｙ軸周りの回転及びＺ軸周りの回転のそれぞれをθＸ、θＹ及びθＺとする。Ｘ軸、Ｙ軸、Ｚ軸に関する制御又は駆動は、それぞれ、Ｘ軸に平行な方向、Ｙ軸に平行な方向、Ｚ軸に平行な方向に関する制御又は駆動を意味する。また、θＸ軸、θＹ軸、θＺ軸に関する制御又は駆動は、それぞれ、Ｘ軸に平行な軸の周りの回転、Ｙ軸に平行な軸の周りの回転、Ｚ軸に平行な軸の周りの回転に関する制御又は駆動を意味する。また、位置は、Ｘ軸、Ｙ軸及びＺ軸の座標に基づいて特定される情報であり、姿勢は、θＸ軸、θＹ軸及びθＺ軸の値で特定される情報である。位置決めは、位置及び／又は姿勢を制御することを意味する。
【００１９】
膜形成装置ＩＭＰは、基板Ｓを保持する基板保持部ＳＨと、基板保持部ＳＨを駆動することで基板Ｓを移動させる基板駆動機構ＳＤと、基板駆動機構ＳＤを支持するベースＳＢとを有する。また、膜形成装置ＩＭＰは、型Ｍを保持する型保持部ＭＨと、型保持部ＭＨを駆動することで型Ｍを移動させる型駆動機構ＭＤとを有する。
【００２０】
基板駆動機構ＳＤ及び型駆動機構ＭＤは、基板Ｓと型Ｍとの相対位置が調整されるように、基板Ｓ及び型Ｍの少なくとも一方を移動させる相対移動機構を構成する。かかる相対移動機構による基板Ｓと型Ｍとの相対位置の調整は、基板Ｓの上の硬化性組成物ＩＭと型Ｍとの接触のための駆動、及び、基板Ｓの上の硬化した硬化性組成物ＩＭからの型Ｍの分離のための駆動を含む。また、相対移動機構による基板Ｓと型Ｍとの相対位置の調整は、基板Ｓと型Ｍとの位置合わせを含む。基板駆動機構ＳＤは、基板Ｓを複数の軸（例えば、Ｘ軸、Ｙ軸及びθＺ軸の３軸、好ましくは、Ｘ軸、Ｙ軸、Ｚ軸、θＸ軸、θＹ軸及びθＺ軸の６軸）に関して駆動するように構成されている。型駆動機構ＭＤは、型Ｍを複数の軸（例えば、Ｚ軸、θＸ軸及びθＹ軸の３軸、好ましくは、Ｘ軸、Ｙ軸、Ｚ軸、θＸ軸、θＹ軸及びθＺ軸の６軸）に関して駆動するように構成されている。
【００２１】
膜形成装置ＩＭＰは、基板Ｓと型Ｍとの間の空間に充填された硬化性組成物ＩＭを硬化させるための硬化部ＣＵを有する。硬化部ＣＵは、例えば、型Ｍを介して硬化性組成物ＩＭに硬化用のエネルギーを与えることによって、基板Ｓの上の硬化性組成物ＩＭを硬化させる。
【００２２】
膜形成装置ＩＭＰは、型Ｍの裏面側（基板Ｓに対面する面の反対側）に空間ＳＰを形成するための透過部材ＴＲを有する。透過部材ＴＲは、硬化部ＣＵからの硬化用のエネルギーを透過させる材料で構成され、基板Ｓの上の硬化性組成物ＩＭに対して硬化用のエネルギーを与えることを可能にする。
【００２３】
膜形成装置ＩＭＰは、空間ＳＰの圧力を制御することによって、型ＭのＺ軸方向への変形を制御する圧力制御部ＰＣを有する。例えば、圧力制御部ＰＣが空間ＳＰの圧力を大気圧よりも高くすることによって、型Ｍは、基板Ｓに向けて凸形状に変形する。
【００２４】
膜形成装置ＩＭＰは、基板Ｓの上に硬化性組成物ＩＭを配置、供給又は分配するためのディスペンサＤＳＰを有する。但し、膜形成装置ＩＭＰには、他の装置によって硬化性組成物ＩＭが配置された基板Ｓが供給（搬入）されてもよい。この場合、膜形成装置ＩＭＰは、ディスペンサＤＳＰを有していなくてもよい。
【００２５】
膜形成装置ＩＭＰは、基板Ｓ（又は基板Ｓのショット領域）と型Ｍとの位置ずれ（位置合わせ誤差）を計測するためのアライメントスコープＡＳを有していてもよい。
【００２６】
シミュレーション装置１は、膜形成装置ＩＭＰにおいて実行される処理における硬化性組成物ＩＭの挙動を予測する計算を実行する。具体的には、シミュレーション装置１は、基板Ｓの上に配置された硬化性組成物ＩＭの複数の液滴と型Ｍとを接触させ、基板Ｓと型Ｍとの間の空間に硬化性組成物ＩＭの膜を形成する処理における硬化性組成物ＩＭの挙動を予測する計算を実行する。
【００２７】
シミュレーション装置１は、例えば、汎用又は専用のコンピュータにシミュレーションプログラム２１を組み込むことによって構成される。なお、シミュレーション装置１は、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）などのＰＬＤ（Ｐｒｏｇｒａｍｍａｂｌｅ  Ｌｏｇｉｃ  Ｄｅｖｉｃｅ）によって構成されてもよい。また、シミュレーション装置１は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）によって構成されてもよい。
【００２８】
シミュレーション装置１は、本実施形態では、プロセッサ１０と、メモリ２０と、ディスプレイ３０と、入力デバイス４０とを有するコンピュータにおいて、メモリ２０にシミュレーションプログラム２１を格納することによって構成される。メモリ２０は、半導体メモリであってもよいし、ハードディスクなどのディスクであってもよいし、他の形態のメモリであってもよい。シミュレーションプログラム２１は、コンピュータによって読み取り可能なメモリ媒体に格納されて、又は、電気通信回線などの通信設備を介してシミュレーション装置１に提供されてもよい。
【００２９】
図２は、一般的な手法において、基板Ｓと型Ｍとの間の空間における硬化性組成物ＩＭの挙動をシミュレーションする際に定義される計算格子を例示する図である。一般的な手法において、計算格子は、計算のための最小単位である計算要素の集合体である。図２において、格子を構成するように配置された複数の微小な矩形の各々が計算要素である。基板Ｓの解析対象の領域、例えば、ショット領域に計算格子が定義される。一般的な手法では、硬化性組成物ＩＭの液滴の挙動を解析するために、硬化性組成物ＩＭの液滴の寸法よりも十分に小さい計算要素からなる計算格子が定義され、硬化性組成物ＩＭの液滴の体積が各々の計算要素の体積に対して占める割合として表現される。例えば、半導体製造における標準的な画角である２６ｍｍ×３３ｍｍの領域に、数十ｎｍ前後の膜厚となる硬化性組成物の液膜を、数ｐＬの液滴を配置することで形成する場合を考える。この場合、数万滴の液滴を同時に取り扱うことになるため、このような微細な計算要素からなる計算格子を定義すると、計算量が膨大なものとなり、許容可能な時間内に計算結果が得られることは期待できない。
【００３０】
基板Ｓと型Ｍとの間の空間に最終的に形成される硬化性組成物ＩＭの膜には、気泡（欠陥）が存在しないことが要求される。硬化性組成物ＩＭにおける気泡の有無をシミュレーションするためには、硬化性組成物ＩＭの液滴間に気泡が閉じ込められる際の液滴の挙動を正確に追跡する必要がある。ここで、一般的な手法では、図２に示すように、硬化性組成物ＩＭの液滴の形状を計算要素に対する液滴の体積分率の分布として表現する。従って、隣接する液滴同士の接合を精度よく計算するためには、液滴の輪郭を十分に表現（解像）できるような非常に小さい計算要素を定義しなければならず、計算コストの増大を招いてしまう。このような計算を１つのショット領域の全域にわたって行うことを考えると、許容可能な時間内に計算結果が得られることは期待できない。
【００３１】
図３を参照して、シミュレーション装置１によって実行されるシミュレーション方法を説明する。かかるシミュレーション方法は、工程として、Ｓ００１、Ｓ００２、Ｓ００３、Ｓ００４、Ｓ００５及びＳ００６を含む。Ｓ００１は、シミュレーションに必要な条件（シミュレーションの条件）を設定する工程である。Ｓ００２は、Ｓ００１で設定されたシミュレーション条件に基づいて、硬化性組成物ＩＭの初期状態を設定する工程である。Ｓ００１及びＳ００２は、それらを併せた１つの工程、例えば、準備工程として理解されてもよい。Ｓ００３は、型Ｍの運動を計算して、型Ｍの位置（基板Ｓと型Ｍとの間の距離）を更新（計算）する工程である。Ｓ００４は、Ｓ００３で更新された型Ｍの位置に基づいて、硬化性組成物ＩＭの複数の液滴のそれぞれについて、型Ｍで押し広げられる液滴の挙動（液滴の流動）を計算する。Ｓ００５は、Ｓ００４で計算された液滴の挙動に基づいて、硬化性組成物ＩＭの複数の液滴のうち隣接する液滴（液滴同士）が接合したかどうかを判定する工程である。Ｓ００６は、計算（シミュレーション）における時刻が終了時刻に達したかどうかを判定する。計算における時刻が終了時刻に達していなければ、時刻を次の時刻に進めて、Ｓ００３に移行する。一方、計算における時刻が終了時刻に達していれば、シミュレーション方法が終了する。シミュレーション装置１は、Ｓ００１、Ｓ００２、Ｓ００３、Ｓ００４、Ｓ００５及びＳ００６をそれぞれ実行するハードウェア要素の集合体として理解されてもよい。
【００３２】
以下、Ｓ００１、Ｓ００２、Ｓ００３、Ｓ００４、Ｓ００５及びＳ００６のそれぞれについて詳細に説明する。
【００３３】
Ｓ００１では、シミュレーションに必要な条件として各種のパラメータが設定される。パラメータは、基板Ｓの上における硬化性組成物ＩＭの液滴の配置、各液滴の体積、硬化性組成物ＩＭの物性値、型Ｍの表面の凹凸（例えば、パターン領域ＰＲのパターンの情報）に関する情報、基板Ｓの表面の凹凸に関する情報などを含む。また、パラメータは、型駆動部ＭＤが型Ｍに与える力の時間プロファイル、圧力制御部ＰＣが空間ＳＰ（型Ｍ）に与える圧力のプロファイルなどを含む。
【００３４】
Ｓ００２では、硬化性組成物ＩＭの複数の液滴のそれぞれの初期状態が設定される。かかる初期状態は、基板Ｓの上に配置された硬化性組成物ＩＭの各液滴が濡れ広がった際の各液滴の輪郭（の形状）及び高さを含む。かかる初期状態は、硬化性組成物ＩＭの物性値を用いて静的な釣り合い状態を仮定して算出することも可能である。また、硬化性組成物ＩＭの物性値に加えて、硬化性組成物ＩＭの液滴を基板Ｓの上に配置してからの経過時間などを入力とし、一般的な流体シミュレーションを実行することで、動的な濡れ広がり挙動から初期状態を算出することも可能である。
【００３５】
本実施形態におけるシミュレーション方法では、硬化性組成物ＩＭの各液滴は、図４に示すように、液滴要素ＤＲＰとしてモデル化される。図４は、硬化性組成物ＩＭの液滴要素ＤＲＰの概念を示す図である。図４を参照するに、計算領域内のｉ番目の液滴要素をＤＲＰｉと表記し、以後、下付き添え字ｉは、液滴要素ＤＲＰの番号を表すものとする。
【００３６】
まず、硬化性組成物ＩＭの液滴要素の内部に代表点が設定される。かかる代表点の座標をＣｉ（ｘ０，ｙ０）とする。硬化性組成物ＩＭの液滴要素の代表点は、液滴の重心としてもよいし、液滴の重心とは異なる点（位置）としてもよいが、液滴の輪郭の内側に設定する必要がある。また、硬化性組成物ＩＭの液滴要素の代表点を通る直線で規定される基準線が設定される。そして、硬化性組成物ＩＭの液滴要素の代表点から見て、角度θ（代表点と液滴の輪郭上の点とを結ぶ線と基準線とのなす角）の位置にある液滴要素の輪郭上（外周上）の点までの距離を半径ｒ（θ）と表現する。半径ｒ（θ）は、角度θごとに異なる値となる。ここで、液滴要素の輪郭上の各点が隣接する液滴要素に接合しているかどうかを示す情報を併せて保持する。隣接する液滴要素に接合した輪郭上の点は、その時点で位置が固定される。換言すれば、硬化性組成物ＩＭの複数の液滴のうち隣接する液滴と接合した液滴について、かかる液滴の接合した部分に対応する、代表点から輪郭上の点までの距離、即ち、半径ｒ（θ）を固定する。半径ｒ（θ）が固定された角度θの領域を、図４に太線で示すように、固定領域ＦＩＸｉとする。一方、半径ｒ（θ）が固定されていない角度θの領域を、図４に実線で示すように、自由領域ＦＲＥｉとする。硬化性組成物ＩＭの液滴の初期状態において、全ての角度θは自由領域に属する。
【００３７】
実際のプログラムとして本実施形態のシミュレーション方法を実装する場合には、角度θを有限個に分割して取り扱う（即ち、液滴の輪郭を定義するために、液滴の輪郭上に有限個の点が設定される）ことが考えられる。図５は、２０個の角度θｋ（ｋ＝１〜２０）で定義（分割）された硬化性組成物ＩＭの液滴要素を示す図である。この際、角度θｋは、３６０度を等分割して設定してもよいし、任意の角度を設定してもよい。有限個の角度で表される液滴要素の輪郭上の隣接する点同士の間の輪郭を求める際には、任意の補間を適用することができる。例えば、隣接する輪郭上の点を線で結んでもよいし、より高次の補間を適用することもできる。
【００３８】
このように、本実施形態では、硬化性組成物ＩＭの液滴要素の輪郭上の有限個の点の間の輪郭（形状）も精度よく表現することが可能である。一方、計算格子を定義する一般的な手法（従来技術）において、液滴要素の輪郭の表現精度を本実施形態と同等にするためには、図６に示すように、非常に多くの計算要素が必要となる。具体的には、図５に示す液滴要素を計算格子上で表現することを考えると、少なくとも２４個×２２個（＝５２８個）の計算要素が必要となり、計算要素と同数の情報を格納する必要があることがわかる。それに対して、本実施形態では、硬化性組成物ＩＭの液滴要素を２０個の情報で表現することができる。実際には、硬化性組成物ＩＭの複数の液滴を同時に扱う必要があるため、本実施形態と一般的な手法との差が更に増大する。本実施形態では、液滴の数に比例する情報量を扱うだけでよいが、一般的な手法では、液滴と液滴との間の領域にも計算要素を設定する必要があるため、更に多くの計算格子が必要となり、計算時間の増大を招いてしまう。このように、本実施形態は、計算格子を定義する一般的な手法と比較して、著しく少ない計算コストで、硬化性組成物ＩＭの液滴の輪郭を表現することが可能である。
【００３９】
Ｓ００３では、型Ｍの運動が計算され、型Ｍの位置が更新される。型Ｍの運動は、硬化性組成物ＩＭの液滴や液滴同士が接合した液膜が押しつぶされる際に発生する力、気体の流動に起因する力、型Ｍに印加する荷重、型Ｍの弾性変形の影響などを考慮した力学計算によって計算される。また、予め求めた計算における各時刻に対応する型Ｍの位置を用いて、現在の時刻に対応する型Ｍの位置を更新してもよい。
【００４０】
計算格子を定義する一般的な手法において、基板Ｓの上に配置された硬化性組成物ＩＭの液滴の広がりを計算するためには、液滴の流動に関する方程式を解く必要がある。このような硬化性組成物ＩＭの液滴の流動に関する方程式は、一般に、全計算要素に対する連立方程式となる。従って、硬化性組成物ＩＭの液滴の接合を精度よく判定可能な計算要素（サイズ）を用いた場合、計算コストが膨大になるため、現実的な時間の範囲内で計算結果を得ることが期待できない。
【００４１】
そこで、本実施形態では、一般的な手法とは異なり、硬化性組成物ＩＭの液滴の流動に関する方程式（連立方程式）を解かずに、液滴の挙動を計算する手法を提案し、計算コストの大幅な削減及び計算精度の両立を図っている。
【００４２】
Ｓ００４では、型Ｍで押し広げられる液滴要素ＤＲＰの挙動が計算される。Ｓ００４は、液滴要素ＤＲＰが型Ｍと接触しているかどうかを判定する工程を含む。Ｓ００２で得られた液滴要素ＤＲＰｉの高さｈｄｒｐ，ｉと、液滴要素ＤＲＰｉの代表点（ｘ０，ｙ０）における型Ｍと基板Ｓとの間の距離ｈｉとを比較し、以下の式（１）を満たす場合には、液滴要素ＤＲＰｉが型Ｍと接触したと判定する。
【００４３】
【数１】
2020205416.tif 000003
【００４４】
一方、式（１）を満たさない場合、計算における現在の時刻において、液滴要素ＤＲＰｉは、型Ｍと接触していないと判定する。この場合、液滴要素ＤＲＰｉの挙動の計算は行われない。
【００４５】
型Ｍと接触したと判定された液滴要素ＤＲＰｉについて、型Ｍの運動によって押し広げられる挙動を計算する。この過程において、硬化性組成物ＩＭの液滴の体積は保存（維持）される。従って、液滴要素ＤＲＰｉの体積Ｖｉと、現在の時刻における液滴要素位置における距離ｈｉとを用いて、現在の時刻における液滴要素ＤＲＰｉの面積Ｓｎｅｗは、以下の式（２）で表すことができる。
【００４６】
【数２】
2020205416.tif 000004
【００４７】
１つ前の時刻における液滴要素ＤＲＰｉの面積をＳｏｌｄとすると、１つ前の時刻から現在の時刻にいたる間の型Ｍの運動によって、液滴要素ＤＲＰｉの面積がΔＳだけ変化することになる。従って、液滴要素ＤＲＰｉの面積の変化ΔＳは、以下の式（３）で表される。図７は、型Ｍの運動による液滴要素ＤＲＰｉの面積の変化ΔＳを模式的に示す図である。
【００４８】
【数３】
2020205416.tif 000005
【００４９】
一方、液滴要素の面積Ｓは、以下の式（４）に示すように、半径ｒ（θ）と関連づけることができる。
【００５０】
【数４】
2020205416.tif 000006
【００５１】
本実施形態のシミュレーション方法では、液滴要素の広がりは、半径ｒ（θ）の増減によって表現される。隣接する液滴要素と接合した領域（部分）、即ち、固定領域ＦＩＸに含まれる角度方向に対しては、図８に示すように、液滴要素と隣接する液滴要素とが互いに広がろうとする（液滴要素がぶつかり合う）。この場合、液滴要素の広がりと隣接する液滴要素の広がりとが打ち消され、上述した角度方向には、それ以上広がらないとみなせる。従って、液滴要素の面積の変化ΔＳは、隣接する液滴と接合していない自由領域ＦＲＥに含まれる角度方向に広がることで現れる。
【００５２】
また、液滴要素の内部は型Ｍで一様に押しつぶされていることから、液滴要素は、自由領域ＦＲＥに属する角度方向に一様な大きさで広がるとみなすことができる。従って、１つ前の時刻からの液滴要素の面積の変化ΔＳは、以下の式（５）で表すことができる。
【００５３】
【数５】
2020205416.tif 000007
【００５４】
式（５）において、Δｒは、半径ｒ（θ）の変化（増分）であり、自由領域ＦＲＥに属する角度に対して共通の値である。
【００５５】
但し、基板Ｓに凹凸構造が存在する場合や型Ｍに微細パターンが形成されている場合には、基板Ｓの上に配置された硬化性組成物ＩＭの液滴（液滴要素）は、図９に示すように、異方的に広がることがわかっている。このような場合には、異方性を表す重みづけ関数ｆ（θ）を導入し、液滴要素の面積の変化ΔＳを、以下の式（６）で表せばよい。
【００５６】
【数６】
2020205416.tif 000008
【００５７】
重みづけ関数ｆ（θ）は、角度θと、型Ｍに形成された微細パターンなどによって生じる特徴方向（例えば、微細パターンの方向）とに依存して計算される。なお、基板Ｓの上に配置された硬化性組成物ＩＭの液滴（液滴要素）が等方的に広がる場合には、重みづけ関数ｆ（θ）は、以下の式（７）を満たす。
【００５８】
【数７】
2020205416.tif 000009
【００５９】
半径ｒ（θ）の変化Δｒは、以下の式（８）で求めることができる。
【００６０】
【数８】
2020205416.tif 000010
【００６１】
このようにして求めた半径ｒ（θ）の変化Δｒを用いれば、現在の時刻における半径ｒｎｅｗ（θ）は、１つ前の時刻における半径ｒｏｌｄ（θ）を用いて、以下の式（９）で表すことができる。
【００６２】
【数９】
2020205416.tif 000011
【００６３】
このように、硬化性組成物ＩＭの複数の液滴のそれぞれについて、型Ｍと基板Ｓとの間の距離の変化に応じて、液滴の体積及び型Ｍと基板Ｓとの間の距離から得られる液滴の面積に一致するように、液滴の代表点から輪郭上の点までの距離を求める。従って、本実施形態のシミュレーション方法によれば、硬化性組成物ＩＭの液滴の流動に関する方程式を解くことなく、液滴の挙動を表現（計算）することができる。
【００６４】
Ｓ００５では、隣接する液滴要素が接合したかどうかが判定される。Ｓ００４において、液滴要素の輪郭を計算した結果、自由領域ＦＲＥに属する角度θの輪郭上の点が、隣接する液滴要素の内部（輪郭の内側）に入ることが発生する。この場合、かかる角度θにおける半径ｒ（θ）が固定される（即ち、液滴の接合した部分に対応する、代表点から輪郭上の点までの距離を固定する）。換言すれば、かかる角度θは、固定領域ＦＩＸに含まれるようになり、それ以降の時刻において、かかる角度θの方向には、硬化性組成物ＩＭの液滴要素の広がり（流動）は発生しない。Ｓ００５では、隣接する液滴要素の全ての組について、上述したように、液滴同士が接合したかどうかを判定する。
【００６５】
図１０を参照して、隣接する液滴要素が接合しているかどうか、即ち、液滴の輪郭上の点が隣接する液滴の輪郭の内側に位置しているかどうかを判定する処理について説明する。まず、液滴要素ＤＲＰｉ（第１液滴）に着目し、液滴要素ＤＲＰｉの自由領域ＦＲＥｉに属する角度方向の輪郭上の点Ｐを考える。液滴要素ＤＲＰｉに隣接する液滴要素を液滴要素ＤＲＰｊ（第２液滴）とし、点Ｐと液滴要素ＤＲＰｊの代表点Ｃｊ（中心）とを結ぶ線分ＰＣｊの長さを求める。また、線分ＰＣｊと液滴要素ＤＲＰｊの基準線とのなす角（角度）θｊを求め、角度θｊにおける液滴要素ＤＲＰｊの半径ＱＣｊの長さを求める。そして、半径ＱＣｊの長さと線分ＰＣｊの長さとを比較して、半径ＱＣｊの長さが線分ＰＣｊの長さよりも長い場合に、液滴要素ＤＲＰｉの輪郭上の点Ｐが隣接する液滴要素ＤＲＰｊの輪郭の内側に位置している、即ち、接合していると判定する。一方、半径ＱＣｊの長さが線分ＰＣｊの長さよりも短い場合には、液滴要素ＤＲＰｉの輪郭上の点Ｐが隣接する液滴要素ＤＲＰｊの輪郭の内側に位置していない、即ち、接合していないと判定する。なお、図１０では、液滴要素ＤＲＰｉの輪郭上の点Ｐが隣接する液滴要素ＤＲＰｊの内部に大きく貫入しているように図示しているが、これは、本実施形態の特徴を強調して表現したものである。実際の計算においては、時刻の間隔を十分に細かくすることで、液滴要素ＤＲＰｉの輪郭上の点Ｐが隣接する液滴要素ＤＲＰｊの内部に貫入する貫入量を無視可能な大きさにすることができる。
【００６６】
隣接する液滴要素の間において、固定領域ＦＩＸの一部を共有する場合、両者は接合しているとみなすことができる。複数の液滴要素が接合していくと、図１１に示すように、液滴要素の間に囲まれた領域に気泡が閉じ込められ、気泡領域が形成される。図１１は、複数の液滴要素の接合によって気泡の閉じ込めが発生した状態を示す図である。このように、本実施形態のシミュレーション方法によれば、気泡が閉じ込められるタイミングを精度よく捉えることができる。また、液滴要素の間に囲まれた領域に閉じ込められた気泡の体積を簡便に算出することが可能である。具体的には、図１１に示すように、三角形Ｃ１Ｃ２Ｃ３の面積から、斜線で示す液滴要素の面積を差し引くことで、気泡領域の面積を算出することができる。また、気泡が閉じ込められた時刻における型Ｍと基板Ｓとの間の距離を、気泡領域の面積にかけることによって、液滴要素の間に囲まれた領域に閉じ込められた気泡の体積を算出することができる。
【００６７】
図１２は、本実施形態のシミュレーション方法を実装したシミュレーション装置１によって硬化性組成物ＩＭの液滴の挙動（広がり）を計算した一例を示す図である。型Ｍと基板Ｓとの間の距離は、図１２の中心ほど近く、図１２の中心から離れるほど遠くなっている。図１２を参照するに、基板Ｓの上の硬化性組成物ＩＭの液滴の配置に応じて、複雑な液滴同士の接合の様子などが表現できていることがわかる。
【００６８】
Ｓ００３、Ｓ００４及びＳ００５を含む計算工程は、予め設定された複数の時刻について実行される。複数の時刻は、例えば、型Ｍが初期位置から降下を開始する時刻から、複数の液滴と接触し、複数の液滴がつぶされながら広がり、複数の液滴が相互に結合し、最終的に１つの膜を形成し、硬化性組成物の硬化がなされるべき時刻までの期間内で任意に設定される。複数の時刻は、典型的には、一定の時間間隔で設定される。
【００６９】
Ｓ００６では、計算における時刻が終了時刻に達したかどうかが判定される。上述したように、計算における時刻が終了時刻に達していなければ、時刻を次の時刻に進めてＳ００３に移行し、計算における時刻が終了時刻に達していれば、シミュレーション方法が終了する。一例において、Ｓ００６では、現在の時刻が指定された時間刻み分だけ進められて、新たな時刻とされる。そして、新たな時刻が終了時刻に達した場合、計算が完了したと判定される。
【００７０】
Ｓ００３からＳ００５の順番は入れ替えてもよい。例えば、Ｓ００４、Ｓ００５、Ｓ００３の順番で実施する場合、硬化性組成物ＩＭの液滴の挙動の計算に用いる型Ｍの位置を、１つ前の時刻に対応する型Ｍの位置とすることによって、同様に処理することができる。
【００７１】
上述したように、本実施形態によれば、基板Ｓの上に配置された硬化性組成物ＩＭの各液滴が型Ｍによって押し広げられ、互いに接合していく挙動を、非常に少ない計算コストで実現可能である。また、硬化性組成物ＩＭの液滴同士の接合の過程において、液滴の間に閉じ込められる気泡の体積についても精度よく算出することができる。
【００７２】
なお、本実施形態においては、外周上の点を、代表点を通る基準線と代表点と外周上の点とを結ぶ線とのなす角、及び、代表点と外周上の点とを結ぶ線の長さで表現したが、これに限定されるものではない。例えば、上述した角度の代わりに計算空間のｘ軸やｙ軸、任意の基準線とのなす角を用いてもよい。また、上述した角度と長さを用いた表現の代わりに代表点から外周上に向かう、或いは、その逆向きのベクトル表現を用いたり、代表点と外周上の点の座標セットを用いたりしてもよい。このように、代表点から外周上の点までの向きと距離がわかる他の方法で表現してもよい。
【００７３】
本発明は、上述の実施形態の１つ以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサがプログラムを読出し実行する処理でも実現可能である。また、１つ以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【００７４】
発明は上記実施形態に制限されるものではなく、発明の精神及び範囲から離脱することなく、様々な変更及び変形が可能である。従って、発明の範囲を公にするために請求項を添付する。
【符号の説明】
【００７５】
１：シミュレーション装置    １０：プロセッサ    ２０：メモリ    ２１：シミュレーションプログラム    ３０：ディスプレイ    ４０：入力デバイス
【図１】
2020205416.tif 000012
【図２】
2020205416.tif 000013
【図３】
2020205416.tif 000014
【図４】
2020205416.tif 000015
【図５】
2020205416.tif 000016
【図６】
2020205416.tif 000017
【図７】
2020205416.tif 000018
【図８】
2020205416.tif 000019
【図９】
2020205416.tif 000020
【図１０】
2020205416.tif 000021
【図１１】
2020205416.tif 000022
【図１２】
2020205416.tif 000023