(57)【要約】【課題】経験的パラメータを用いずに分子レベルの性質を適切に取り込んだメソスケールのシミュレーション方法及びプログラムを提供する。【解決手段】本シミュレーション方法及びプログラムは、例えば材料のメソスケール情報を予測するために好適に利用される。具体的には、揺動散逸理論に基づいて、ミクロスケールの運動方程式からミクロ情報を削減した運動方程式を導出するステップを有する。これにより、計算時間及び計算資源を低減しつつ、メソスケール情報を得ることが可能となる。次に、ミクロ情報を削減した運動方程式中の係数を近似式により求めて、メソスケール運動方程式を導出する。本シミュレーション方法及びプログラムによれば、経験的パラメータを用いることなく、分子論的要因が注目する物性に対して及ぼす効果を精度良く予測することができる。【選択図】図１
2008052409.tif 000002
【特許請求の範囲】
【請求項１】
材料のメソスケールのシミュレーション方法であって、
前記材料の分子レベル情報を取得し、前記分子レベル情報に基づいてミクロスケールの運動方程式を作成するステップと、
揺動散逸理論に基づいて、前記ミクロスケールの運動方程式からミクロ情報を削減した運動方程式を導出するステップと、
前記ミクロ情報を削減した運動方程式中の係数を近似式により求めて、メソスケール運動方程式を導出するステップと、
前記メソスケール運動方程式を用いてシミュレーションを実行するステップと、を有することを特徴とするメソスケールのシミュレーション方法。
【請求項２】
前記ミクロ情報を削減した運動方程式を導出するステップは、前記ミクロスケールの運動方程式に対して射影演算するステップと、前記射影演算した後の運動方程式に対してミクロ時間の積分を行うステップと、を有することを特徴とする請求項１に記載のメソスケールのシミュレーション方法。
【請求項３】
材料のメソスケールのシミュレーションプログラムであって、コンピュータ上で実行されることにより、
前記材料の分子レベル情報を取得し、前記分子レベル情報に基づいてミクロスケールの運動方程式を作成するステップと、
揺動散逸理論に基づいて、前記ミクロスケールの運動方程式からミクロ情報を削減した運動方程式を導出するステップと、
前記ミクロ情報を削減した運動方程式中の係数を近似式により求めて、メソスケール運動方程式を導出するステップと、
前記メソスケール運動方程式を用いてシミュレーションを実行するステップと、を実行することを特徴とするメソスケールのシミュレーションプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、材料のメソスケールのシミュレーションを行う手法に関する。
【背景技術】
【０００２】
従来から、高分子材料などのメソスケールにおける物性を、シミュレーションにより予測することが行われている。例えば、特許文献１には、高分子材料の形態を、分子動力学を用いたシミュレーションにより予測する技術が記載されている。また、特許文献２には、高分子材料の相分離構造を予測する技術が記載されている。その他にも、本発明に関連のある技術が特許文献３に記載されている。
【０００３】
【特許文献１】特開２００３−１０５０９０号公報
【特許文献２】特開２００５−３２０５８号公報
【特許文献３】ＷＯ９７／０６４９６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
上記した特許文献１乃至３に記載された技術においては、基本的には、シミュレーションに適用する現象論的方程式（例えばLangevin方程式など）に含まれる粘性係数や相溶性などの物性パラメータ（以下、「経験的パラメータ」とも呼ぶ。）を、予め実験や文献などより取得しておく必要があった。しかしながら、既知の実験データが存在しなかったり、実験が困難であったりして、経験的パラメータを取得することができない場合があった。また、特許文献１乃至３に記載された技術では、注目する物性に対して、分子レベルの要因がどのように影響するかを明瞭にすることが困難であった。
【０００５】
本発明は、上記のような課題を解決するためになされたものであり、経験的パラメータを用いずに分子レベルの性質を適切に取り込んだメソスケールのシミュレーション方法及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の１つの観点では、材料のメソスケールのシミュレーション方法は、前記材料の分子レベル情報を取得し、前記分子レベル情報に基づいてミクロスケールの運動方程式を作成するステップと、揺動散逸理論に基づいて、前記ミクロスケールの運動方程式からミクロ情報を削減した運動方程式を導出するステップと、前記ミクロ情報を削減した運動方程式中の係数を近似式により求めて、メソスケール運動方程式を導出するステップと、前記メソスケール運動方程式を用いてシミュレーションを実行するステップと、を有することを特徴とする。
【０００７】
上記のメソスケールのシミュレーション方法において好適には、前記ミクロ情報を削減した運動方程式を導出するステップは、前記ミクロスケールの運動方程式に対して射影演算するステップと、前記射影演算した後の運動方程式に対してミクロ時間の積分を行うステップと、を有する。
【０００８】
本発明の他の観点では、材料のメソスケールのシミュレーションプログラムは、コンピュータ上で実行されることにより、前記材料の分子レベル情報を取得し、前記分子レベル情報に基づいてミクロスケールの運動方程式を作成するステップと、揺動散逸理論に基づいて、前記ミクロスケールの運動方程式からミクロ情報を削減した運動方程式を導出するステップと、前記ミクロ情報を削減した運動方程式中の係数を近似式により求めて、メソスケール運動方程式を導出するステップと、前記メソスケール運動方程式を用いてシミュレーションを実行するステップと、を実行することを特徴とする。
【０００９】
上記のシミュレーション方法及びプログラムは、例えば材料のメソスケール情報を予測などするために好適に利用される。具体的には、揺動散逸理論に基づいて、ミクロスケールの運動方程式からミクロ情報を削減した運動方程式を導出するステップを有する。これにより、計算時間及び計算資源を低減しつつ、メソスケール情報を得ることが可能となる。次に、ミクロ情報を削減した運動方程式中の係数を近似式により求めて、メソスケール運動方程式を導出する。上記のシミュレーション方法及びプログラムによれば、経験的パラメータを用いることなく、分子論的要因が注目する物性に対して及ぼす効果を精度良く予測することができる。即ち、注目する物性に対する分子レベルの各要因の効果やそのメカニズムなどを予測することが可能となる。
【発明を実施するための最良の形態】
【００１０】
以下、図面を参照して本発明の好適な実施の形態について説明する。
【００１１】
図１は、本発明によるメソスケールシミュレーションプログラムの概略フローチャートである。メソスケールシミュレーションプログラムは、材料のメソスケール情報などを求めるために好適に利用される。また、このプログラムは、例えばパーソナルコンピュータ（ＰＣ）などによって実行される。なお、メソスケールは、ナノメートル程度のスケールを示すものとする。
【００１２】
まず、ステップＳ１０１では、材料の分子レベルの情報を入力する。具体的には、構成分子などのミクロ要因、及び構成・境界条件や熱力学的条件などのマクロ要因を入力する。そして、ステップＳ１０２に進む。
【００１３】
ステップＳ１０２では、解析対象のハミルトニアンを記述し、注目する物性値に対するミクロスケール運動方程式を書き下す。このミクロスケール運動方程式は、分子論的運動方程式（Liouville方程式）に対応する。そして、ステップＳ１０３に進む。
【００１４】
ステップＳ１０３では、揺動散逸理論に基づいて、ステップＳ１０２で得られたミクロスケール運動方程式からミクロ情報を削減した運動方程式を導出する。ミクロスケール運動方程式は、分子論的条件を考慮して注目する物性値の挙動を記述しているが、これを解いて得られるのはミクロ情報であり、メソスケール情報は得られない。実際のスケールでこれを解くには、膨大な計算時間と計算資源（経験的パラメータなど）が必要となり、現実的ではない。したがって、本実施形態では、揺動散逸理論に基づいて、ミクロスケール運動方程式からミクロ情報を削減した運動方程式を導出する。具体的には、２段階にわたってミクロ情報の削減を行う。まず、注目する物性への射影演算を行うことによってミクロ情報を削減する。次に、こうして得られた運動方程式に対してミクロ時間の積分を行うことによって、ミクロ情報を更に削減する。このようなミクロ情報の削減を行うことによって、分子論的情報を取り込んでミクロスケールからメソスケールに情報が削減される。なお、揺動散逸理論とは、簡単に言うと、分子における揺らぎと散逸との間を結ぶ定理である。そして、ステップＳ１０４に進む。
【００１５】
ステップＳ１０４では、注目する物性値に対するメソスケール運動方程式を導出する。ここでは、ステップＳ１０３で得られた運動方程式の係数（いわゆる現象論的メソスケール運動方程式における経験的パラメータに対応する）を、計算可能な形の近似式により算出する。そして、算出された係数の近似式をモンテカルロ法などによって数値計算することによって、メソスケール運動方程式を導出する。このメソスケール運動方程式は、現象論的メソスケール運動方程式に対応する分子論的条件を考慮したものとなっている。ステップＳ１０４で行う処理は、数個の原子を１つの粒子として粗視化を行うＤＰＤ（Dissipative Particle Dynamics）法などに相当する。そして、ステップＳ１０５に進む。
【００１６】
ステップＳ１０５では、ステップＳ１０４で得られたメソスケール運動方程式を用いて、物性値挙動のシミュレーションを実行する。この後、ステップＳ１０５で得られたシミュレーション結果を解析することも可能である。
【００１７】
このように、本発明によるメソスケールシミュレーションプログラムによれば、経験的パラメータ（物性パラメータ）を用いることなく、分子論的要因が注目する物性に対して及ぼす効果を精度良く予測することができる。即ち、上記のプログラムによれば、注目する物性に対する分子レベルの各要因の効果やそのメカニズムなどを予測することが可能となる。
【実施例】
【００１８】
次に、本発明のメソスケールシミュレーションプログラムの好適な実施例について説明する。
【００１９】
図２は、本実施例においてメソスケールシミュレーションの解析対象とする高分子溶液を示す。図２では、符号１０が溶媒分子を示し、符号１１が高分子内原子を示し、符号１２が高分子鎖を示し、符号１３が高分子セグメントを示している。
【００２０】
図３は、上記した高分子溶液をシミュレーションするために用いるメソスケールシミュレーションプログラムのフローチャートである。ここでは、メソスケールシミュレーションを行うことによって、高分子の溶媒中での分散状態などを計算する。
【００２１】
まず、ステップＳ２０１では、ミクロ条件及びマクロ条件からなるシミュレーション条件を設定する。ミクロ条件としては、溶媒分子における溶媒相互の分子間ポテンシャル及び溶媒−高分子間の分子間ポテンシャルや、高分子構造や、セグメント内及びセグメント外の原子間ポテンシャルなどを設定する。また、マクロ条件としては、セル形状や分子数や境界条件などのシステム構成や、濃度や温度・圧力条件などを設定する。なお、図２に示す例では、溶媒分子は１粒子、高分子はセグメント（部分鎖）の集合として扱う。そして、ステップＳ２０２に進む。
【００２２】
ステップＳ２０２では、ミクロスケール運動方程式を記述する。具体的には、系のハミルトニアン（粒子間ポテンシャルに対応する）を記述し、高分子セグメントのミクロスケール運動方程式を書き下す。そして、ステップＳ２０３に進む。
【００２３】
ステップＳ２０３では、揺動散逸理論に基づいて、メソスケールシミュレーションのための情報の粗視化を行う。具体的には、２段階にわたってミクロ情報の削減を行う。まず、射影演算子を用いて、ミクロスケール運動方程式の注目する物性値以外の自由度（情報）を消去することによってミクロ情報を削減する。即ち、ミクロ自由度のセグメントへの射影演算を行う（ミクロ情報の削減（１））。次に、ミクロスケール運動方程式をミクロスケールな時間積分を行い、即ちミクロ時間の積分を行い、メソスケールな時間スケールの運動方程式を導出する（ミクロ情報の削減（２））。これにより、情報粗視化されたセグメントの運動方程式が得られる。そして、ステップＳ２０４に進む。
【００２４】
ステップＳ２０４では、セグメントに対するメソスケール運動方程式を導出する。具体的には、ステップＳ２０３で得られた運動方程式の係数（現象論的メソスケール運動方程式の経験的パラメータに対応する）を、計算可能な近似式に変形する。そして、この係数の近似式を、モンテカルロ法などによって数値計算することによって、メソスケール運動方程式を導出する。そして、ステップＳ２０５に進む。
【００２５】
ステップＳ２０５では、ステップＳ２０４で得られたセグメントのメソスケール運動方程式を用いて、注目する物性値挙動のシミュレーションを実行する。そして、ステップＳ２０６に進む。ステップＳ２０６では、シミュレーション結果を解析する。
【００２６】
上記のメソスケールシミュレーションプログラムによれば、経験的パラメータを用いることなく、分子論的要因が注目する物性に対して及ぼす効果を精度良く予測することが可能となる。
【００２７】
以下で、上記したプログラムのステップにおける主要部分について詳細に説明する。
【００２８】
（１）ミクロスケール運動方程式
ここでは、ステップＳ２０２で行うミクロスケール運動方程式について説明する。まず、粒子間ポテンシャルに対応するハミルトニアンは、式（１）で表される。
【００２９】
【数１】
2008052409.tif 000003
式（１）において、「ｉ」及び「ｊ」は溶媒分子の番号を示し、「Ｋ」及び「Ｌ」は高分子中のセグメントの番号を示し、「ｍ」及び「ｎ」はセグメント中の原子の番号を示し、「ｖｉ」はｉ−番溶媒分子の速度を示し、「ｖＫｍ」はＫセグメント中のｍ−番原子の速度を示し、「ｒｉ，ｊ」はｉ−番溶媒分子とｊ−番溶媒分子との相対位置ベクトルを示し、「ｒｉ，Ｋｍ」はｉ−番溶媒分子とＫセグメント中のｍ−番原子との相対位置ベクトルを示し、「ｒＫｍ，Ｌｎ」はＫセグメント中のｍ−番原子とＬセグメント中のｎ−番原子との相対位置ベクトルを示し、「ＶＡＡ(ｒｉ，ｊ)」は溶媒分子同士の相互作用エネルギーを示し、「ＶＡＢ(ｒｉ，Ｋｍ)」は溶媒分子とセグメント中の原子との相互作用エネルギーを示し、「ＶＢＢ(ｒＫｍ，Ｌｎ)」はセグメント中の原子同士の相互作用エネルギーを示している。
【００３０】
ここで、注目する物理量としてセグメントの重心座標ＲＧ、及びセグメント平方長Ｒ２は、式（２）で示すように表現される。
【００３１】
【数２】
2008052409.tif 000004
このとき、Ｋ−セグメント重心ＲＧＫのミクロスケール運動方程式は、式（３）で表され、Ｋ−セグメント平方長Ｒ２Ｋのミクロスケール運動方程式は、式（４）で表される。
【００３２】
【数３】
2008052409.tif 000005
【００３３】
【数４】
2008052409.tif 000006
（２）メソスケールシミュレーションのための情報の粗視化
次に、ステップＳ２０３で行う、メソスケールシミュレーションのための情報の粗視化について説明する。具体的には、２段階にわたって行うミクロ情報の削減を説明する。即ち、射影演算によるミクロ情報の削減と、ミクロ時間の積分によるミクロ情報の削減とについて説明する。
【００３４】
（ａ）注目する物理量への射影演算によるミクロ情報の削減
ＲＧＫへの射影、及びＲ２Ｋへの射影を、それぞれ以下の式（５）、式（６）で定義する。
【００３５】
【数５】
2008052409.tif 000007
【００３６】
【数６】
2008052409.tif 000008
ここで、式（５）及び式（６）中の＜・・・＞は、Γ−空間での統計平均を示しており、具体的には式（７）で表される。
【００３７】
【数７】
2008052409.tif 000009
次に、式（３）に示すＲＧに対するミクロスケール運動方程式に対してＲＧＫへの射影を行うと、式（８）が得られる。
【００３８】
【数８】
2008052409.tif 000010
式（８）で示す運動方程式は、３つの部分に分けて考えることができる。ここで、この３つの部分を、Liouville演算子ｉＬを用いて表す。具体的には、Liouville演算子ｉＬは式（９）で定義される。
【００３９】
【数９】
2008052409.tif 000011
式（８）の右辺の第１項は、ＲＧＫ’(ｔ)に働く他のセグメントとの相互作用による力に対応する運動部分である。その係数ｉΩＧＫＫ’は、式（１０）で表される。
【００４０】
【数１０】
2008052409.tif 000012
式（８）の右辺の第２項は、ＲＧＫ’(ｔ)に働く溶媒との相互作用による摩擦抵抗に対応する運動部分である。その摩擦係数ΦＧＫＫ’(ｔ)は、式（１１）で表される。
【００４１】
【数１１】
2008052409.tif 000013
式（８）の右辺の第３項は、ＲＧＫ’(ｔ)に働くランダムな揺動力部分である。その揺動力ＦＧＫ(ｔ)は、式（１２）で表される。
【００４２】
【数１２】
2008052409.tif 000014
同様にして、式（４）に示すＲ２に対するミクロスケール運動方程式に対してＲ２Ｋへの射影を行うと、式（１３）が得られる。
【００４３】
【数１３】
2008052409.tif 000015
この場合、式（１３）の右辺の第１項は、Ｒ２Ｋ’(ｔ)に働く他のセグメントとの相互作用による力に対応する運動部分である。その係数ｉΩ２ＫＫ’は、式（１４）で表される。
【００４４】
【数１４】
2008052409.tif 000016
式（１３）の右辺の第２項は、Ｒ２Ｋ’(ｔ)に働く溶媒との相互作用による摩擦抵抗に対応する運動部分である。その摩擦係数Φ２ＫＫ’(ｔ)は、式（１５）で表される。
【００４５】
【数１５】
2008052409.tif 000017
式（１３）の右辺の第３項は、Ｒ２Ｋ’(ｔ)に働くランダムな揺動力部分である。その揺動力Ｆ２Ｋ(ｔ)は、式（１６）で表される。
【００４６】
【数１６】
2008052409.tif 000018
上記した式（８）及び式（１３）が、ＲＧ及びＲ２に対するミクロ自由度（溶媒分子の情報など）を削減した運動方程式に対応する。
【００４７】
（ｂ）ミクロ時間の積分によるメソスケールな時間スケールの運動方程式
セグメントの溶媒分子との衝突による揺らぎ周期に比べて、セグメントのメソスケール的運動ＲＧはゆっくりである。また、ミクロ時間のランダム力の積分は、以下の式（１７）で表されるように、概ね「０」と考えることができる。
【００４８】
【数１７】
2008052409.tif 000019
よって、上記した射影演算によるミクロ情報の削減によって得られた式（８）をメソスケールな時間で積分すると、式（１８）によって表すことができる。
【００４９】
【数１８】
2008052409.tif 000020
同様に、Ｒ２Ｋ(ｔ)に対応する式（１３）をメソスケールな時間で積分すると、式（１９）によって表すことができる。
【００５０】
【数１９】
2008052409.tif 000021
これらの式（１８）及び式（１９）が、情報粗視化されたセグメントの運動方程式に対応する。この場合、式（１８）及び式（１９）中のセグメント間相互作用の係数を示す「ｉΩＧＫＫ’」及び「ｉΩ２ＫＫ’」（式（１０）、式（１４）参照）と、摩擦抵抗を示す「ΦＧＫＫ’」、「Φ２ＫＫ’」（式（１１）、式（１５）参照）とが計算できれば、メソスケール運動方程式が得られる。
【００５１】
（３）セグメントに対するメソスケール運動方程式の導出
次に、ステップＳ２０４で行う、メソスケール運動方程式の導出について説明する。ここでは、ステップＳ２０３で得られた運動方程式中の係数をモンテカルロ法によって数値計算する例を説明する。
【００５２】
まず、上記した式（１０）及び式（１４）で示したセグメント間相互作用の係数に対応する「ｉΩＧＫＫ’」及び「ｉΩ２ＫＫ’」は、式（７）で示す統計平均をモンテカルロ平均で置き換えることによって数値計算できると言える。
【００５３】
これに対して、以下の式（２０）及び式（２１）で表される摩擦抵抗を示す「ΦＧＫＫ’」及び「Φ２ＫＫ’」は、このままでは計算できない。
【００５４】
【数２０】
2008052409.tif 000022
【００５５】
【数２１】
2008052409.tif 000023
そのため、演算子の指数関数部をベキ展開して近似を行う。これにより、「ΦＧＫＫ’」は式（２２）で表され、「Φ２ＫＫ’」は式（２３）で表される。
【００５６】
【数２２】
2008052409.tif 000024
【００５７】
【数２３】
2008052409.tif 000025
更に、式（１８）中の「ΦＧＫＫ’」の積分部分を近似する。具体的には、０からｔまでの積分を、０から∞までの積分に変形する。これにより、式（２４）が得られる。
【００５８】
【数２４】
2008052409.tif 000026
同様に、式（１９）中の「Φ２ＫＫ’」の積分部分を近似することによって、式（２５）が得られる。
【００５９】
【数２５】
2008052409.tif 000027
このようにして得られた式（２４）及び式（２５）を用いて、モンテカルロ法により数値計算する。
【００６０】
次に、メソスケール運動方程式の係数に対する計算可能な表現が得られたので、シミュレーション対象の具体的モデルについて熱力学条件を設定したモンテカルロ法により数値計算を行う。ここで、対象をモデル化してモンテカルロ法により数値計算して得られた式（１０）、式（１４）、式（２４）、式（２５）の値を、それぞれ「ｉωＧＫＫ’」、「ｉω２ＫＫ’」、「φＧＫＫ’」、「φ２ＫＫ’」と表現する。こうすることにより、式（１８）及び式（１９）で示すミクロ自由度を消去したメソスケールな時間スケールの運動方程式は、それぞれ式（２６）及び式（２７）で示すメソスケール運動方程式によって表される。
【００６１】
【数２６】
2008052409.tif 000028
【００６２】
【数２７】
2008052409.tif 000029
このようにして得られた式（２６）及び式（２７）のメソスケール運動方程式を基礎方程式として、従来のＤＰＤ法などと同様に、シミュレーションを実行することができる。これにより、分子論的要因が注目する物性に対して及ぼす効果を精度良く予測することが可能となる。
【図面の簡単な説明】
【００６３】
【図１】本発明によるメソスケールシミュレーションプログラムの概略フローチャートを示す。
【図２】メソスケールシミュレーションの解析対象とする高分子溶液を示す。
【図３】高分子溶液をシミュレーションするために用いるプログラムのフローチャートである。
【符号の説明】
【００６４】
１０  溶媒分子
１１  高分子内原子
１２  高分子鎖
１３  高分子セグメント
【図１】
2008052409.tif 000030
【図２】
2008052409.tif 000031
【図３】
2008052409.tif 000032