高分子の物性値計算方法及び物性値計算装置
【課題】 高分子の物性値の算出を高精度かつ効率的に行う。
【解決手段】 本発明の物性値計算方法では、高分子に応じた分子のパラメータ算出用の全原子モデルを生成するための情報、及び高分子の物性値算出用の全原子モデルを生成するための情報を入力し（Ｓ０１）パラメータ算出用の全原子モデルを生成し（Ｓ０２）パラメータ算出用の全原子モデルを用いて分子動力学法によるシミュレーションを行い（Ｓ０３）当該シミュレーションの結果により示される原子の分布に基づいて物性値算出用の全原子モデル用のポテンシャルパラメータを算出し（Ｓ０４〜Ｓ０７）物性値算出用の全原子モデルを生成し（Ｓ０８）物性値算出用の全原子モデルを用いて、算出されたポテンシャルパラメータにて分子動力学法によるシミュレーションを行い（Ｓ０９）当該シミュレーションの結果に基づいて物性値を算出する（Ｓ１０）。
【解決手段】 本発明の物性値計算方法では、高分子に応じた分子のパラメータ算出用の全原子モデルを生成するための情報、及び高分子の物性値算出用の全原子モデルを生成するための情報を入力し（Ｓ０１）パラメータ算出用の全原子モデルを生成し（Ｓ０２）パラメータ算出用の全原子モデルを用いて分子動力学法によるシミュレーションを行い（Ｓ０３）当該シミュレーションの結果により示される原子の分布に基づいて物性値算出用の全原子モデル用のポテンシャルパラメータを算出し（Ｓ０４〜Ｓ０７）物性値算出用の全原子モデルを生成し（Ｓ０８）物性値算出用の全原子モデルを用いて、算出されたポテンシャルパラメータにて分子動力学法によるシミュレーションを行い（Ｓ０９）当該シミュレーションの結果に基づいて物性値を算出する（Ｓ１０）。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、高分子の物性値を算出する物性値計算方法及び物性値計算装置に関する。
【背景技術】
【０００２】
従来から、分子動力学法によるシミュレーションを用いて、高分子の密度、ガラス転移温度、ヤング率等の物性値を算出することが行われている。分子動力学法によるシミュレーションを行うためには、ポテンシャルパラメータを割り当てることが必要になる。例えば非特許文献１では、シミュレーションに用いられるポテンシャルパラメータを既存のものを用いるとしている。また、非特許文献２では、実験値と計算値とが合わない場合には上記のポテンシャルパラメータを、量子化学計算を用いて算出することとしている。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】D.Pavel et al., Journal of Polymer Science; Part B; Polymer Physics,Volume 37, 2334-2352, 1999
【非特許文献２】C.F. Fan et al., Macromolecules 1994, 27, 2383-2391
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、非特許文献１に記載された方法では、用いられるポテンシャルパラメータが正確な値でない場合、算出される物性値の精度が悪いものになってしまうという問題がある。また、非特許文献２に記載された方法では、算出される物性値の精度は良好であるものの、量子化学計算に要する時間が（分子動力学法によるシミュレーションに比べて）非常に長いという問題点がある。
【０００５】
本発明は、上記の問題点を解決するためになされたものであり、高分子の物性値の算出を、高精度かつ効率的に行うことができる物性値計算方法及び物性値計算装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記の目的を達成するために、本発明に係る高分子の物性値計算方法は、コンピュータにより構成され、高分子の物性値を算出する物性値計算装置によって実行される高分子の物性値計算方法であって、高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力ステップと、入力ステップにおいて入力された情報に基づいて、パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成ステップと、パラメータ算出用モデル生成ステップにおいて作成されたパラメータ算出用の全原子モデルを用いて、入力ステップにおいて入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーションステップと、パラメータ算出用シミュレーションステップにおいて行われたシミュレーションの結果により示される原子の分布に基づいて、物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出ステップと、入力ステップにおいて入力された情報に基づいて、物性値算出用の全原子モデルを生成する物性値算出用モデル生成ステップと、物性値算出用モデル生成ステップにおいて作成された物性値算出用の全原子モデルを用いて、パラメータ算出ステップにおいて算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーションステップと、物性値算出用シミュレーションステップにおいて行われたシミュレーションの結果に基づいて、物性値を算出して出力する物性値算出ステップと、を含むことを特徴とする。
【０００７】
本発明に係る高分子の物性値計算方法では、物性値を算出するために行われる物性値算出用の全原子モデルを用いたシミュレーションのためのポテンシャルパラメータが、物性値の算出対象の高分子に応じた分子の、パラメータ算出用の全原子モデルを用いたシミュレーションの結果に基づいて算出される。これにより、上記の物性値算出用の全原子モデルのポテンシャルパラメータを精度のよいものにすることができ、その結果、高分子の物性値の算出を高精度に行うことができる。その一方で、物性値算出用の全原子モデルのポテンシャルパラメータは、高分子に応じた分子の、パラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションによって算出されるため、量子化学計算等と比較して短時間で算出される。即ち、高分子の物性値の算出を効率的に行うことができる。
【０００８】
パラメータ算出用シミュレーションステップにおいて行われるパラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用シミュレーションステップにおいて行われる物性値算出用の全原子モデルを用いたシミュレーションの計算時間より少なくすることができる。パラメータ算出用の全原子モデルのシミュレーションは、物性値算出用の全原子モデル用のポテンシャルパラメータを算出することを目的としたものであり、物性値を算出することを目的とした物性値算出用の全原子モデルのシミュレーションとは異なり、短い計算時間で十分だからである。従って、この構成によれば、物性値算出用の全原子モデルのポテンシャルパラメータが更に効率的に算出されるため、これにより、高分子の物性値の算出を更に効率的に行うことができる。ここで、計算時間とは、シミュレーションにおける時間間隔とシミュレーションの繰り返しの回数との積によって決まる時間である。
【０００９】
入力ステップにおいて、高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報として、高分子の規模を小さくした分子の全原子モデルを生成するための情報を入力することができる。この構成によれば、高分子の規模を小さくした分子の、パラメータ算出用の全原子モデルを用いたシミュレーションによって、物性値算出用の全原子モデルのポテンシャルパラメータが算出されるため、更に短時間で算出される。即ち、高分子の物性値の算出を更に効率的に行うことができる。
【００１０】
ところで、本発明は、上記のように高分子の物性値計算方法の発明として記述できる他に、以下のように高分子の物性値計算装置の発明としても記述することができる。これはカテゴリが異なるだけで、実質的に同一の発明であり、同様の作用及び効果を奏する。
【００１１】
即ち、本発明に係る高分子の物性値計算装置は、コンピュータにより構成され、高分子の物性値を算出する物性値計算装置であって、高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力手段と、入力手段によって入力された情報に基づいて、パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成手段と、パラメータ算出用モデル生成手段によって作成されたパラメータ算出用の全原子モデルを用いて、入力手段によって入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーション手段と、パラメータ算出用シミュレーション手段によって行われたシミュレーションの結果により示される原子の分布に基づいて、物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出手段と、入力手段によって入力された情報に基づいて、物性値算出用の全原子モデルを生成する物性値算出用モデル生成手段と、物性値算出用モデル生成手段によって作成された物性値算出用の全原子モデルを用いて、パラメータ算出手段によって算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーション手段と、物性値算出用シミュレーション手段によって行われたシミュレーションの結果に基づいて、物性値を算出して出力する物性値算出手段と、を備えることを特徴とする。
【発明の効果】
【００１２】
本発明では、物性値を算出するために行われる物性値算出用の全原子モデルを用いたシミュレーションのためのポテンシャルパラメータが、物性値の算出対象の高分子に応じた分子の、パラメータ算出用の全原子モデルを用いたシミュレーションの結果に基づいて算出される。これにより、上記の物性値算出用の全原子モデルのポテンシャルパラメータを精度のよいものにすることができ、その結果、高分子の物性値の算出を高精度に行うことができる。その一方で、物性値算出用の全原子モデルのポテンシャルパラメータは、高分子に応じた分子の、パラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションによって算出されるため、量子化学計算等と比較して短時間で算出される。即ち、高分子の物性値の算出を効率的に行うことができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る高分子の物性値計算装置の構成図である。
【図２】本発明の実施形態に係る高分子の物性値計算方法を示すフローチャートである。
【発明を実施するための形態】
【００１４】
以下、図面と共に本発明に係る高分子の物性値計算方法及び物性値計算装置の好適な実施形態について詳細に説明する。なお、図面の説明においては同一要素には同一符号を付し、重複する説明を省略する。
【００１５】
本実施形態に係る高分子の物性値計算方法は、高分子の物性を示す値（物性値）を算出するものである。物性値の算出対象となる高分子は、通常、複数の高分子を含み、例えば、高分子を構成するモノマー数が１０００程度以上、又は分子量が１００００程度以上のものである。算出される物性値としては、例えば、高分子の密度、ガラス転移温度及びヤング率である。本実施形態に係る高分子の物性値計算方法では、コンピュータ上でモデル化した高分子に対して分子動力学法によるシミュレーションを行い、シミュレーション結果に基づいて演算を行うことによって物性値を算出する。
【００１６】
図１に本実施形態に係る高分子の物性値計算方法が実行される高分子の物性値計算装置１０を示す。物性値計算装置１０は、具体的には、ワークステーションやＰＣ（Personal Computer）等のコンピュータにより構成される。物性値計算装置１０は、例えばＣＰＵ（Central Processing Unit）やメモリ等のハードウェアにより構成されており、これらの構成要素が動作することにより後述する物性値計算装置１０としての機能が発揮される。なお、コンピュータに実行させるプログラムが物性値計算装置１０において実行されることにより、本実施形態に係る高分子の物性値計算方法が行われてもよい。
【００１７】
図１に示すように物性値計算装置１０は、記憶部１１と、入力部１２と、パラメータ算出用モデル生成部１３と、パラメータ算出用シミュレーション部１４と、パラメータ算出部１５と、物性値算出用モデル生成部１６と、物性値算出用シミュレーション部１７と、物性値算出部１８とを備えて構成される。また、物性値計算装置１０は、外部装置２０と接続されており、外部装置２０から情報が入力される。
【００１８】
記憶部１１は、高分子の物性値の計算に必要な情報を記憶する手段であり、物性値計算装置１０の各構成要素１２〜１８から必要に応じてアクセスされる。具体的には、記憶部１１は、入力部１２に入力された情報や物性値の計算の過程の演算結果等が記憶される。記憶部１１は、例えば、物性値計算装置１０が備えるメモリにより実現される。
【００１９】
入力部１２は、高分子の物性値の計算に必要な情報を入力する入力手段である。具体的には、入力部１２は、物性値計算装置１０におけるシミュレーションに用いられる高分子のモデルを生成するための情報が入力される。物性値計算装置１０におけるシミュレーションは、２つのモデルに対して行われる。２つのモデルは、高分子を構成する原子個々をモデルの構成要素として扱う全原子モデルであり、パラメータ算出用の全原子モデル、及び高分子の物性値算出用の全原子モデルである。
【００２０】
ここで、パラメータ算出用の全原子モデルは、物性値の算出対象となる高分子に応じた分子の全原子モデルであるとする。望ましくは、物性値の算出対象となる高分子よりも規模が小さい分子の全原子モデルであるとする。上記の全原子モデルは、単に算出対象となる高分子よりも規模が小さいだけで、原子の組成やその他の特徴は当該高分子と同じものであるとする。上記の物性値算出用の全原子モデルは、物性値の算出対象となる高分子をモデル化したものである。モデルを生成するための情報として、具体的には、高分子を構成する原子の組成、モノマー数及びポリマー数等がある。
【００２１】
高分子のパラメータ算出用の全原子モデルの生成は、例えば市販されている分子構造体モデルを用いて、モノマー構造を生成し、モノマー構造を基に行うこととしてもよい。セルの形状に関しては、立方体であることが、偏りのないモデル構造計算の観点から好ましい。パラメータ算出用の高分子の全原子モデルを構成するモノマー数については、特に制限はないが、少なすぎると実際の値との乖離が大きくなり好ましくない。モノマー数は多い方が精度上では好ましいが、多すぎるとシミュレーションを開始してから終了するまでに要する時間（シミュレーションに要する実時間、以下「実時間」と称することがある）が長くなりすぎる。上記を考慮して、具体的には、モノマー数としては５〜２００程度が好ましい。セルの境界に関しては、３次元方向に同じセルが繰り返して存在する３次元周期境界条件を用いることが、計算精度の観点から好ましい。高分子数についても、特に制限はないが、望ましくは高分子数は１〜１０である。高分子数を１〜１０とすることで、実時間を短縮化できる。上記の全原子モデルは、単に算出対象となる高分子よりも規模が小さいだけで、原子の組成やその他の特徴は当該高分子と同じものであるとする。
【００２２】
また、入力部１２は、上記のパラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションを行う際に設定される（割り当てられる）当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報が入力される。パラメータ算出用の全原子モデル用のポテンシャルパラメータは、例えば、文献等に記載されている既存のものが用いられる。
【００２３】
上記の情報の入力は、例えば、ユーザによる操作によって外部装置２０から行われる。また、予め高分子の物性値算出装置１０に記憶された情報を（ユーザによる操作等をトリガとして）入力することとしてもよい。入力部１２は、入力した各情報を記憶部１１に記憶させる。入力部１２は、各情報を入力するとその旨をパラメータ算出用モデル生成部１３に通知する。
【００２４】
パラメータ算出用モデル生成部１３は、情報を入力した旨が入力部１２から通知されると、入力部１２によって入力された情報に基づいて、上記のパラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成手段である。パラメータ算出用モデル生成部１３は、入力部１２によって入力されたパラメータ算出用の全原子モデルを生成するのに必要な情報を、記憶部１１から読み出してパラメータ算出用の全原子モデルの生成に利用する。パラメータ算出用の全原子モデルの生成自体は、従来の方法と同様に行われる。パラメータ算出用モデル生成部１３は、生成したパラメータ算出用の全原子モデルを示す情報を記憶部１１に記憶させる。パラメータ算出用モデル生成部１３は、パラメータ算出用の全原子モデルを生成するとその旨をパラメータ算出用シミュレーション部１４に通知する。
【００２５】
パラメータ算出用シミュレーション部１４は、パラメータ算出用の全原子モデルを生成した旨がパラメータ算出用モデル生成部１３から通知されると、パラメータ算出用モデル生成部１３によって作成されたパラメータ算出用の全原子モデルを用いて分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーション手段である。パラメータ算出用シミュレーション部１４は、上記のパラメータ算出用の全原子モデルに対して、入力部１２によって入力された情報により示されるポテンシャルパラメータを設定して上記のシミュレーションを行う。パラメータ算出用シミュレーション部１４は、パラメータ算出用の全原子モデルを示す情報、及び当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報を記憶部１１から読み出し、シミュレーションを行う。また、分子動力学法によるシミュレーションに用いられる原子や分子間の相互作用力を表す計算式等を、パラメータ算出用シミュレーション部１４に予め設定しておく。
【００２６】
パラメータ算出用の全原子モデルのシミュレーションにおける時間間隔（時間刻み）は、短いほど精度が高くなり好ましいが、実時間が長くなるため、０．１〜１０ｆｓ（フェムト秒）程度とするのがよい。例えば、時間間隔が１ｆｓの値に設定される（即ち、１時間単位が１ｆｓとされる）場合、計算時間（上記の時間間隔とシミュレーションの繰り返しの回数との積によって決まる時間）は特に制限はないが、１〜１００ｐｓ（ピコ秒）（１０３〜１０５時間単位）が好ましい。計算時間が短すぎると、物性値算出用の全原子モデルのポテンシャル精度が悪化し、計算時間が長すぎると、実時間も長くなるという観点から好ましくない。上記の時間間隔及び繰り返しの回数の情報については、パラメータ算出用シミュレーション部１４が予め記憶しており、当該記憶された情報に基づいてシミュレーションが行われる。
【００２７】
パラメータ算出用の全原子モデルのシミュレーションの計算時間は、後述する物性値算出用の全原子モデルのシミュレーションの計算時間よりも短いものである。即ち、パラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用の全原子モデルを用いたシミュレーションにおける計算時間よりも少ない。パラメータ算出用の全原子モデルのシミュレーションは、物性値算出用の全原子モデル用のポテンシャルパラメータを算出することを目的としたものであり、物性値を算出することを目的とした物性値算出用の全原子モデルのシミュレーションとは異なり、短い計算時間で十分だからである。但し、パラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用の全原子モデルを用いたシミュレーションにおける計算時間よりも必ずしも少なくなっている必要はない。なお、上記の分子動力学法によるシミュレーション自体は、従来の方法と同様に行われる。パラメータ算出用シミュレーション部１４は、実行したシミュレーションの結果を示す情報を記憶部１１に記憶させる。シミュレーションの結果を示す情報には、分子を構成する原子の（相対）位置を示す情報が含まれる。パラメータ算出用シミュレーション部１４は、シミュレーションの実行が完了するとその旨をパラメータ算出部１５に通知する。
【００２８】
パラメータ算出部１５は、シミュレーションの実行が完了した旨がパラメータ算出用シミュレーション部１４から通知されると、パラメータ算出用シミュレーション部１４によって行われたシミュレーションの結果により示される原子の分布に基づいて、物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出手段である。パラメータ算出部１５は、パラメータ算出用の全原子モデルを用いたシミュレーションの結果を示す情報を記憶部１１から読み出して、以下のように物性値算出用の全原子モデル用のポテンシャルパラメータを算出する。
【００２９】
まず、パラメータ算出部１５は、シミュレーションの結果により得られる原子の相対位置から原子間の間隔（長さ）や結合角等の分布を算出する。続いて、予め設定された分布（例えば、正規分布）に対する、算出した分布のフィッティングを行う。続いて、フィッティングした結果をポテンシャル関数に変換する。続いて、得られたポテンシャル関数から、物性値算出用の全原子モデル用のポテンシャルパラメータを算出する。なお、上記の一連のポテンシャルパラメータの算出の処理は、従来の方法と同様に行われる。パラメータ算出部１５は、算出したポテンシャルパラメータを示す情報を記憶部１１に記憶させる。パラメータ算出部１５は、ポテンシャルパラメータを算出するとその旨を物性値算出用モデル生成部１６に通知する。
【００３０】
物性値算出用モデル生成部１６は、ポテンシャルパラメータが算出された旨がパラメータ算出部１５から通知されると、入力部１２によって入力された情報に基づいて、上記の物性値算出用の全原子モデルを生成する物性値算出用モデル生成手段である。物性値算出用モデル生成部１６は、入力部１２によって入力された物性値算出用の全原子モデルを生成するのに必要な情報を、記憶部１１から読み出して物性値算出用の全原子モデルの生成に利用する。物性値算出用の全原子モデルの生成自体は、従来の方法と同様に行われる。物性値算出用モデル生成部１６は、生成した物性値算出用の全原子モデルを示す情報を記憶部１１に記憶させる。物性値算出用モデル生成部１６は、物性値算出用の全原子モデルを生成するとその旨を物性値算出用シミュレーション部１７に通知する。なお、物性値算出用モデル生成部１６による物性値算出用の全原子モデルの生成は、パラメータ算出部１５によるポテンシャルパラメータの算出処理を前提としていない。従って、物性値算出用の全原子モデルの生成は、必ずしもパラメータ算出部１５からの通知をトリガとして行なわれる必要はなく、記憶部１１に物性値算出用の全原子モデルを生成するのに必要な情報が記憶された後であればどのタイミングで行われてもよい。
【００３１】
上記の物性値算出用の全原子モデルは、物性値の算出対象となる高分子をモデル化したものである。物性値算出用の全原子モデルのモノマー数は、特に制限はないが少なすぎると実際の値との乖離が大きくなり好ましくない。モノマー数は多い方が精度上では好ましいが、多すぎると実時間が長くなりすぎる。上記を考慮して、具体的には、モノマー数としては５〜５００程度が好ましい。ポリマー数は、特に制限はなく１〜５０程度が好ましい。ポリマー数が少なすぎると精度が低下し、多すぎると、実時間が長くなるという観点から好ましくない。物性値算出用の全原子モデルは、上記のパラメータ算出用の全原子モデルと同数（同じサイズ）であってもよいが、パラメータ算出用の全原子モデルよりもモノマー数及びポリマー数の少なくともいずれかを多くすることで、物性値算出の精度を向上させることができる。
【００３２】
物性値算出用シミュレーション部１７は、物性値算出用の全原子モデルを生成した旨が物性値算出用モデル生成部１６から通知されると、物性値算出用モデル生成部１６によって作成されたパラメータ算出用の全原子モデルを用いて分子動力学法によるシミュレーションを行う物性値算出用シミュレーション手段である。物性値算出用シミュレーション部１７は、上記の物性値算出用の全原子モデルに対して、パラメータ算出部１５によって算出されたポテンシャルパラメータを設定して上記のシミュレーションを行う。物性値算出用シミュレーション部１７は、物性値算出用の全原子モデルを示す情報、及び当該物性値算出用の全原子モデル用のポテンシャルパラメータを示す情報を記憶部１１から読み出し、シミュレーションを行う。また、分子動力学法によるシミュレーションに用いられる原子や分子間の相互作用力を表す計算式等を、物性値算出用シミュレーション部１７に予め設定しておく。
【００３３】
物性値算出用の全原子モデルのシミュレーションにおける時間間隔は、短いほど精度が高くなり好ましいが、計算に要する時間が長くなるため、パラメータ算出用の全原子モデルのシミュレーションと同様に０．１〜１０ｆｓ程度とするのがよい。例えば、時間間隔が１ｆｓの値に設定される（即ち、１時間単位が１ｆｓとされる）場合、計算時間は特に制限はないが、１００ｐｓ〜１０ｎｓ（ナノ秒）（１０５〜１０７時間単位）が好ましい。計算時間が短すぎると、算出する高分子の物性の精度が悪化し、計算時間が長すぎると、実時間も長くなるという観点から好ましくない。上記の時間間隔及び繰り返しの回数の情報については、物性値算出用シミュレーション部１７が予め記憶しており、当該記憶された情報に基づいてシミュレーションが行われる。
【００３４】
この物性値算出用の全原子モデルのシミュレーションにおける計算時間は、上述したパラメータ算出用の全原子モデルのシミュレーションにおける計算時間よりも長いものである。物性値算出用の全原子モデルのシミュレーションは、物性値を算出することを目的としているためである。なお、上記の分子動力学法によるシミュレーション自体は、従来の方法と同様に行われる。物性値算出用シミュレーション部１７は、実行したシミュレーションの結果を示す情報を記憶部１１に記憶させる。物性値算出用シミュレーション部１７は、シミュレーションの実行が完了するとその旨を物性値算出部１８に通知する。
【００３５】
物性値算出部１８は、物性値算出用シミュレーション部１７によって行われたシミュレーションの結果に基づいて、物性値を算出して出力する物性値算出手段である。物性値算出部１８は、物性値算出用シミュレーション部１７によるシミュレーションの結果を示す情報を、記憶部１１から読み出して物性値の算出に利用する。物性値算出部１８は、算出する物性値毎に予め記憶している演算式を用いて、シミュレーションの結果に対して演算を行うことによって物性値を算出する。物性値算出部１８は、算出した物性値を示す情報を、例えば、物性値計算装置１０が備える表示装置に表示することによって出力する。なお、算出した物性値を示す情報の出力は、必ずしも表示が伴うものでなくてもよく、例えば、物性値を示す情報が別途用いられるように、自装置１０の記憶部１１や物性値計算装置１０に接続された別の装置に出力することとしてもよい。以上が、本実施形態に係る物性値計算装置１０の構成である。
【００３６】
引き続いて、図２のフローチャートを用いて、物性値計算装置１０により実行される処理である、本実施形態に係る物性値計算方法について説明する。
【００３７】
物性値計算装置１０では、まず、入力部１２によって、高分子の物性値の計算に必要な上述した情報が入力される（Ｓ０１、入力ステップ）。続いて、パラメータ算出用モデル生成部１３によって、入力部１２によって入力された情報を基にパラメータ算出用の全原子モデルが生成される。（Ｓ０２、パラメータ算出用モデル生成ステップ）。続いて、パラメータ算出用シミュレーション部１４によって、パラメータ算出用モデル生成部１３によって生成されたパラメータ算出用の全原子モデルが用いられて、入力部１２によって入力された情報に係るポテンシャルパラメータが設定されて、分子動力学法によるシミュレーションが行われる（Ｓ０３、パラメータ算出用シミュレーションステップ）。
【００３８】
続いて、パラメータ算出部１５によって、パラメータ算出用シミュレーション部１４によって行われたシミュレーションの結果により示される原子の分布に基づいて、以下のように物性値算出用の全原子モデル用のポテンシャルパラメータが算出される。シミュレーションに用いられる、原子や分子間の相互作用力を計算する計算式、及びその他のポテンシャルパラメータには、一般によく知られているもの（例えば川添良幸・三上益弘・大野かおる著『コンピュータ・シミュレーションによる物質化学』共立出版、ｐｐ．５５〜８２に記載されているもの）を用いることができる。分子動力学シミュレーションでは、結合力ポテンシャルと非結合ポテンシャルとを用いることができる。結合ポテンシャルは、結合距離を平衡値に保つ分子間の結合長ポテンシャル、結合角を平衡角に保つ結合角ポテンシャル、及び取り得る２面角を制御するトーション・ポテンシャルを含むことができる。非結合ポテンシャルは、分子内においてレナード・ジョーンズポテンシャルを用いることができる。シミュレーションの結果により得られる原子の相対位置から原子間の結合長、結合角、２面角、分子間距離を求める。分子間距離から非共有ポテンシャルを求めることができる。
【００３９】
まず、シミュレーションの結果により得られる原子の相対位置から原子間の間隔（長さ）や結合角等の分布が算出される（Ｓ０４、パラメータ算出ステップ）。続いて、予め設定された分布に対する、算出した分布のフィッティングが行われる（Ｓ０５、パラメータ算出ステップ）。続いて、フィッティングした結果がポテンシャル関数に変換される（Ｓ０６、パラメータ算出ステップ）。得られたポテンシャル関数から、物性値算出用の全原子モデル用のポテンシャルパラメータが算出される（Ｓ０７、パラメータ算出ステップ）。
【００４０】
続いて、物性値算出用モデル生成部１６によって、入力部１２によって入力された情報を基に物性値算出用の全原子モデルが生成される（Ｓ０８、物性値算出用モデル生成ステップ）。続いて、物性値算出用シミュレーション部１７によって、物性値算出用モデル生成部１６によって生成された物性値算出用の全原子モデルが用いられて、パラメータ算出部１５によって算出されたポテンシャルパラメータが設定されて、分子動力学法によるシミュレーションが行われる（Ｓ０９、物性値算出用シミュレーションステップ）。続いて、物性値算出部１８によって、物性値算出用シミュレーション部１７によって行われたシミュレーションの結果に基づいて、物性値が算出されて出力される（Ｓ１０、物性値算出ステップ）。以上が、本実施形態に係る物性値計算方法である。
【００４１】
上述したように、本実施形態では、物性値を算出するために行われる物性値算出用の全原子モデルを用いたシミュレーションのためのポテンシャルパラメータが、物性値の算出対象の高分子に応じた分子の、パラメータ算出用の全原子モデルを用いたシミュレーションの結果に基づいて算出される。このように、物性値算出用の全原子モデルのポテンシャルパラメータを算出すると、当該ポテンシャルパラメータを精度のよいものにすることができ、その結果、高分子の物性値の算出を高精度に行うことができる。その一方で、物性値算出用の全原子モデルのポテンシャルパラメータは、パラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションによって算出されるため、量子化学計算等と比較して短時間で算出される。即ち、高分子の物性値の算出を効率的に行うことができる。
【００４２】
また、本実施形態のように、物性値の算出対象の高分子の規模を小さくした分子の、パラメータ算出用の全原子モデルを用いて、物性値算出用の全原子モデルのポテンシャルパラメータを算出するためのシミュレーションを行うこととすれば、更に短時間で上記のポテンシャルパラメータが算出される。即ち、高分子の物性値の算出を更に効率的に行うことができる。但し、必ずしも、規模を小さくした分子の、パラメータ算出用の全原子モデルを用いる必要はなく、物性値の算出対象の高分子自体の全原子モデルを用いることとしてもよい。
【００４３】
また、本実施形態のようにパラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用の全原子モデルを用いたシミュレーションの計算時間より少なくすることとすれば、物性値算出用の全原子モデルのポテンシャルパラメータが更に効率的に算出されるため、これにより、高分子の物性値の算出を更に効率的に行うことができる。但し、必ずしも、上記の繰り返しの回数を少なくする必要はない。
【符号の説明】
【００４４】
１０…物性値計算装置、１１…記憶部、１２…入力部、１３…パラメータ算出用モデル生成部、１４…パラメータ算出用シミュレーション部、１５…パラメータ算出部、１６…物性値算出用モデル生成部、１７…物性値算出用シミュレーション部、１８…物性値算出部、２０…外部装置。
【技術分野】
【０００１】
本発明は、高分子の物性値を算出する物性値計算方法及び物性値計算装置に関する。
【背景技術】
【０００２】
従来から、分子動力学法によるシミュレーションを用いて、高分子の密度、ガラス転移温度、ヤング率等の物性値を算出することが行われている。分子動力学法によるシミュレーションを行うためには、ポテンシャルパラメータを割り当てることが必要になる。例えば非特許文献１では、シミュレーションに用いられるポテンシャルパラメータを既存のものを用いるとしている。また、非特許文献２では、実験値と計算値とが合わない場合には上記のポテンシャルパラメータを、量子化学計算を用いて算出することとしている。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】D.Pavel et al., Journal of Polymer Science; Part B; Polymer Physics,Volume 37, 2334-2352, 1999
【非特許文献２】C.F. Fan et al., Macromolecules 1994, 27, 2383-2391
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、非特許文献１に記載された方法では、用いられるポテンシャルパラメータが正確な値でない場合、算出される物性値の精度が悪いものになってしまうという問題がある。また、非特許文献２に記載された方法では、算出される物性値の精度は良好であるものの、量子化学計算に要する時間が（分子動力学法によるシミュレーションに比べて）非常に長いという問題点がある。
【０００５】
本発明は、上記の問題点を解決するためになされたものであり、高分子の物性値の算出を、高精度かつ効率的に行うことができる物性値計算方法及び物性値計算装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記の目的を達成するために、本発明に係る高分子の物性値計算方法は、コンピュータにより構成され、高分子の物性値を算出する物性値計算装置によって実行される高分子の物性値計算方法であって、高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力ステップと、入力ステップにおいて入力された情報に基づいて、パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成ステップと、パラメータ算出用モデル生成ステップにおいて作成されたパラメータ算出用の全原子モデルを用いて、入力ステップにおいて入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーションステップと、パラメータ算出用シミュレーションステップにおいて行われたシミュレーションの結果により示される原子の分布に基づいて、物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出ステップと、入力ステップにおいて入力された情報に基づいて、物性値算出用の全原子モデルを生成する物性値算出用モデル生成ステップと、物性値算出用モデル生成ステップにおいて作成された物性値算出用の全原子モデルを用いて、パラメータ算出ステップにおいて算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーションステップと、物性値算出用シミュレーションステップにおいて行われたシミュレーションの結果に基づいて、物性値を算出して出力する物性値算出ステップと、を含むことを特徴とする。
【０００７】
本発明に係る高分子の物性値計算方法では、物性値を算出するために行われる物性値算出用の全原子モデルを用いたシミュレーションのためのポテンシャルパラメータが、物性値の算出対象の高分子に応じた分子の、パラメータ算出用の全原子モデルを用いたシミュレーションの結果に基づいて算出される。これにより、上記の物性値算出用の全原子モデルのポテンシャルパラメータを精度のよいものにすることができ、その結果、高分子の物性値の算出を高精度に行うことができる。その一方で、物性値算出用の全原子モデルのポテンシャルパラメータは、高分子に応じた分子の、パラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションによって算出されるため、量子化学計算等と比較して短時間で算出される。即ち、高分子の物性値の算出を効率的に行うことができる。
【０００８】
パラメータ算出用シミュレーションステップにおいて行われるパラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用シミュレーションステップにおいて行われる物性値算出用の全原子モデルを用いたシミュレーションの計算時間より少なくすることができる。パラメータ算出用の全原子モデルのシミュレーションは、物性値算出用の全原子モデル用のポテンシャルパラメータを算出することを目的としたものであり、物性値を算出することを目的とした物性値算出用の全原子モデルのシミュレーションとは異なり、短い計算時間で十分だからである。従って、この構成によれば、物性値算出用の全原子モデルのポテンシャルパラメータが更に効率的に算出されるため、これにより、高分子の物性値の算出を更に効率的に行うことができる。ここで、計算時間とは、シミュレーションにおける時間間隔とシミュレーションの繰り返しの回数との積によって決まる時間である。
【０００９】
入力ステップにおいて、高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報として、高分子の規模を小さくした分子の全原子モデルを生成するための情報を入力することができる。この構成によれば、高分子の規模を小さくした分子の、パラメータ算出用の全原子モデルを用いたシミュレーションによって、物性値算出用の全原子モデルのポテンシャルパラメータが算出されるため、更に短時間で算出される。即ち、高分子の物性値の算出を更に効率的に行うことができる。
【００１０】
ところで、本発明は、上記のように高分子の物性値計算方法の発明として記述できる他に、以下のように高分子の物性値計算装置の発明としても記述することができる。これはカテゴリが異なるだけで、実質的に同一の発明であり、同様の作用及び効果を奏する。
【００１１】
即ち、本発明に係る高分子の物性値計算装置は、コンピュータにより構成され、高分子の物性値を算出する物性値計算装置であって、高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力手段と、入力手段によって入力された情報に基づいて、パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成手段と、パラメータ算出用モデル生成手段によって作成されたパラメータ算出用の全原子モデルを用いて、入力手段によって入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーション手段と、パラメータ算出用シミュレーション手段によって行われたシミュレーションの結果により示される原子の分布に基づいて、物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出手段と、入力手段によって入力された情報に基づいて、物性値算出用の全原子モデルを生成する物性値算出用モデル生成手段と、物性値算出用モデル生成手段によって作成された物性値算出用の全原子モデルを用いて、パラメータ算出手段によって算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーション手段と、物性値算出用シミュレーション手段によって行われたシミュレーションの結果に基づいて、物性値を算出して出力する物性値算出手段と、を備えることを特徴とする。
【発明の効果】
【００１２】
本発明では、物性値を算出するために行われる物性値算出用の全原子モデルを用いたシミュレーションのためのポテンシャルパラメータが、物性値の算出対象の高分子に応じた分子の、パラメータ算出用の全原子モデルを用いたシミュレーションの結果に基づいて算出される。これにより、上記の物性値算出用の全原子モデルのポテンシャルパラメータを精度のよいものにすることができ、その結果、高分子の物性値の算出を高精度に行うことができる。その一方で、物性値算出用の全原子モデルのポテンシャルパラメータは、高分子に応じた分子の、パラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションによって算出されるため、量子化学計算等と比較して短時間で算出される。即ち、高分子の物性値の算出を効率的に行うことができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る高分子の物性値計算装置の構成図である。
【図２】本発明の実施形態に係る高分子の物性値計算方法を示すフローチャートである。
【発明を実施するための形態】
【００１４】
以下、図面と共に本発明に係る高分子の物性値計算方法及び物性値計算装置の好適な実施形態について詳細に説明する。なお、図面の説明においては同一要素には同一符号を付し、重複する説明を省略する。
【００１５】
本実施形態に係る高分子の物性値計算方法は、高分子の物性を示す値（物性値）を算出するものである。物性値の算出対象となる高分子は、通常、複数の高分子を含み、例えば、高分子を構成するモノマー数が１０００程度以上、又は分子量が１００００程度以上のものである。算出される物性値としては、例えば、高分子の密度、ガラス転移温度及びヤング率である。本実施形態に係る高分子の物性値計算方法では、コンピュータ上でモデル化した高分子に対して分子動力学法によるシミュレーションを行い、シミュレーション結果に基づいて演算を行うことによって物性値を算出する。
【００１６】
図１に本実施形態に係る高分子の物性値計算方法が実行される高分子の物性値計算装置１０を示す。物性値計算装置１０は、具体的には、ワークステーションやＰＣ（Personal Computer）等のコンピュータにより構成される。物性値計算装置１０は、例えばＣＰＵ（Central Processing Unit）やメモリ等のハードウェアにより構成されており、これらの構成要素が動作することにより後述する物性値計算装置１０としての機能が発揮される。なお、コンピュータに実行させるプログラムが物性値計算装置１０において実行されることにより、本実施形態に係る高分子の物性値計算方法が行われてもよい。
【００１７】
図１に示すように物性値計算装置１０は、記憶部１１と、入力部１２と、パラメータ算出用モデル生成部１３と、パラメータ算出用シミュレーション部１４と、パラメータ算出部１５と、物性値算出用モデル生成部１６と、物性値算出用シミュレーション部１７と、物性値算出部１８とを備えて構成される。また、物性値計算装置１０は、外部装置２０と接続されており、外部装置２０から情報が入力される。
【００１８】
記憶部１１は、高分子の物性値の計算に必要な情報を記憶する手段であり、物性値計算装置１０の各構成要素１２〜１８から必要に応じてアクセスされる。具体的には、記憶部１１は、入力部１２に入力された情報や物性値の計算の過程の演算結果等が記憶される。記憶部１１は、例えば、物性値計算装置１０が備えるメモリにより実現される。
【００１９】
入力部１２は、高分子の物性値の計算に必要な情報を入力する入力手段である。具体的には、入力部１２は、物性値計算装置１０におけるシミュレーションに用いられる高分子のモデルを生成するための情報が入力される。物性値計算装置１０におけるシミュレーションは、２つのモデルに対して行われる。２つのモデルは、高分子を構成する原子個々をモデルの構成要素として扱う全原子モデルであり、パラメータ算出用の全原子モデル、及び高分子の物性値算出用の全原子モデルである。
【００２０】
ここで、パラメータ算出用の全原子モデルは、物性値の算出対象となる高分子に応じた分子の全原子モデルであるとする。望ましくは、物性値の算出対象となる高分子よりも規模が小さい分子の全原子モデルであるとする。上記の全原子モデルは、単に算出対象となる高分子よりも規模が小さいだけで、原子の組成やその他の特徴は当該高分子と同じものであるとする。上記の物性値算出用の全原子モデルは、物性値の算出対象となる高分子をモデル化したものである。モデルを生成するための情報として、具体的には、高分子を構成する原子の組成、モノマー数及びポリマー数等がある。
【００２１】
高分子のパラメータ算出用の全原子モデルの生成は、例えば市販されている分子構造体モデルを用いて、モノマー構造を生成し、モノマー構造を基に行うこととしてもよい。セルの形状に関しては、立方体であることが、偏りのないモデル構造計算の観点から好ましい。パラメータ算出用の高分子の全原子モデルを構成するモノマー数については、特に制限はないが、少なすぎると実際の値との乖離が大きくなり好ましくない。モノマー数は多い方が精度上では好ましいが、多すぎるとシミュレーションを開始してから終了するまでに要する時間（シミュレーションに要する実時間、以下「実時間」と称することがある）が長くなりすぎる。上記を考慮して、具体的には、モノマー数としては５〜２００程度が好ましい。セルの境界に関しては、３次元方向に同じセルが繰り返して存在する３次元周期境界条件を用いることが、計算精度の観点から好ましい。高分子数についても、特に制限はないが、望ましくは高分子数は１〜１０である。高分子数を１〜１０とすることで、実時間を短縮化できる。上記の全原子モデルは、単に算出対象となる高分子よりも規模が小さいだけで、原子の組成やその他の特徴は当該高分子と同じものであるとする。
【００２２】
また、入力部１２は、上記のパラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションを行う際に設定される（割り当てられる）当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報が入力される。パラメータ算出用の全原子モデル用のポテンシャルパラメータは、例えば、文献等に記載されている既存のものが用いられる。
【００２３】
上記の情報の入力は、例えば、ユーザによる操作によって外部装置２０から行われる。また、予め高分子の物性値算出装置１０に記憶された情報を（ユーザによる操作等をトリガとして）入力することとしてもよい。入力部１２は、入力した各情報を記憶部１１に記憶させる。入力部１２は、各情報を入力するとその旨をパラメータ算出用モデル生成部１３に通知する。
【００２４】
パラメータ算出用モデル生成部１３は、情報を入力した旨が入力部１２から通知されると、入力部１２によって入力された情報に基づいて、上記のパラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成手段である。パラメータ算出用モデル生成部１３は、入力部１２によって入力されたパラメータ算出用の全原子モデルを生成するのに必要な情報を、記憶部１１から読み出してパラメータ算出用の全原子モデルの生成に利用する。パラメータ算出用の全原子モデルの生成自体は、従来の方法と同様に行われる。パラメータ算出用モデル生成部１３は、生成したパラメータ算出用の全原子モデルを示す情報を記憶部１１に記憶させる。パラメータ算出用モデル生成部１３は、パラメータ算出用の全原子モデルを生成するとその旨をパラメータ算出用シミュレーション部１４に通知する。
【００２５】
パラメータ算出用シミュレーション部１４は、パラメータ算出用の全原子モデルを生成した旨がパラメータ算出用モデル生成部１３から通知されると、パラメータ算出用モデル生成部１３によって作成されたパラメータ算出用の全原子モデルを用いて分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーション手段である。パラメータ算出用シミュレーション部１４は、上記のパラメータ算出用の全原子モデルに対して、入力部１２によって入力された情報により示されるポテンシャルパラメータを設定して上記のシミュレーションを行う。パラメータ算出用シミュレーション部１４は、パラメータ算出用の全原子モデルを示す情報、及び当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報を記憶部１１から読み出し、シミュレーションを行う。また、分子動力学法によるシミュレーションに用いられる原子や分子間の相互作用力を表す計算式等を、パラメータ算出用シミュレーション部１４に予め設定しておく。
【００２６】
パラメータ算出用の全原子モデルのシミュレーションにおける時間間隔（時間刻み）は、短いほど精度が高くなり好ましいが、実時間が長くなるため、０．１〜１０ｆｓ（フェムト秒）程度とするのがよい。例えば、時間間隔が１ｆｓの値に設定される（即ち、１時間単位が１ｆｓとされる）場合、計算時間（上記の時間間隔とシミュレーションの繰り返しの回数との積によって決まる時間）は特に制限はないが、１〜１００ｐｓ（ピコ秒）（１０３〜１０５時間単位）が好ましい。計算時間が短すぎると、物性値算出用の全原子モデルのポテンシャル精度が悪化し、計算時間が長すぎると、実時間も長くなるという観点から好ましくない。上記の時間間隔及び繰り返しの回数の情報については、パラメータ算出用シミュレーション部１４が予め記憶しており、当該記憶された情報に基づいてシミュレーションが行われる。
【００２７】
パラメータ算出用の全原子モデルのシミュレーションの計算時間は、後述する物性値算出用の全原子モデルのシミュレーションの計算時間よりも短いものである。即ち、パラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用の全原子モデルを用いたシミュレーションにおける計算時間よりも少ない。パラメータ算出用の全原子モデルのシミュレーションは、物性値算出用の全原子モデル用のポテンシャルパラメータを算出することを目的としたものであり、物性値を算出することを目的とした物性値算出用の全原子モデルのシミュレーションとは異なり、短い計算時間で十分だからである。但し、パラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用の全原子モデルを用いたシミュレーションにおける計算時間よりも必ずしも少なくなっている必要はない。なお、上記の分子動力学法によるシミュレーション自体は、従来の方法と同様に行われる。パラメータ算出用シミュレーション部１４は、実行したシミュレーションの結果を示す情報を記憶部１１に記憶させる。シミュレーションの結果を示す情報には、分子を構成する原子の（相対）位置を示す情報が含まれる。パラメータ算出用シミュレーション部１４は、シミュレーションの実行が完了するとその旨をパラメータ算出部１５に通知する。
【００２８】
パラメータ算出部１５は、シミュレーションの実行が完了した旨がパラメータ算出用シミュレーション部１４から通知されると、パラメータ算出用シミュレーション部１４によって行われたシミュレーションの結果により示される原子の分布に基づいて、物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出手段である。パラメータ算出部１５は、パラメータ算出用の全原子モデルを用いたシミュレーションの結果を示す情報を記憶部１１から読み出して、以下のように物性値算出用の全原子モデル用のポテンシャルパラメータを算出する。
【００２９】
まず、パラメータ算出部１５は、シミュレーションの結果により得られる原子の相対位置から原子間の間隔（長さ）や結合角等の分布を算出する。続いて、予め設定された分布（例えば、正規分布）に対する、算出した分布のフィッティングを行う。続いて、フィッティングした結果をポテンシャル関数に変換する。続いて、得られたポテンシャル関数から、物性値算出用の全原子モデル用のポテンシャルパラメータを算出する。なお、上記の一連のポテンシャルパラメータの算出の処理は、従来の方法と同様に行われる。パラメータ算出部１５は、算出したポテンシャルパラメータを示す情報を記憶部１１に記憶させる。パラメータ算出部１５は、ポテンシャルパラメータを算出するとその旨を物性値算出用モデル生成部１６に通知する。
【００３０】
物性値算出用モデル生成部１６は、ポテンシャルパラメータが算出された旨がパラメータ算出部１５から通知されると、入力部１２によって入力された情報に基づいて、上記の物性値算出用の全原子モデルを生成する物性値算出用モデル生成手段である。物性値算出用モデル生成部１６は、入力部１２によって入力された物性値算出用の全原子モデルを生成するのに必要な情報を、記憶部１１から読み出して物性値算出用の全原子モデルの生成に利用する。物性値算出用の全原子モデルの生成自体は、従来の方法と同様に行われる。物性値算出用モデル生成部１６は、生成した物性値算出用の全原子モデルを示す情報を記憶部１１に記憶させる。物性値算出用モデル生成部１６は、物性値算出用の全原子モデルを生成するとその旨を物性値算出用シミュレーション部１７に通知する。なお、物性値算出用モデル生成部１６による物性値算出用の全原子モデルの生成は、パラメータ算出部１５によるポテンシャルパラメータの算出処理を前提としていない。従って、物性値算出用の全原子モデルの生成は、必ずしもパラメータ算出部１５からの通知をトリガとして行なわれる必要はなく、記憶部１１に物性値算出用の全原子モデルを生成するのに必要な情報が記憶された後であればどのタイミングで行われてもよい。
【００３１】
上記の物性値算出用の全原子モデルは、物性値の算出対象となる高分子をモデル化したものである。物性値算出用の全原子モデルのモノマー数は、特に制限はないが少なすぎると実際の値との乖離が大きくなり好ましくない。モノマー数は多い方が精度上では好ましいが、多すぎると実時間が長くなりすぎる。上記を考慮して、具体的には、モノマー数としては５〜５００程度が好ましい。ポリマー数は、特に制限はなく１〜５０程度が好ましい。ポリマー数が少なすぎると精度が低下し、多すぎると、実時間が長くなるという観点から好ましくない。物性値算出用の全原子モデルは、上記のパラメータ算出用の全原子モデルと同数（同じサイズ）であってもよいが、パラメータ算出用の全原子モデルよりもモノマー数及びポリマー数の少なくともいずれかを多くすることで、物性値算出の精度を向上させることができる。
【００３２】
物性値算出用シミュレーション部１７は、物性値算出用の全原子モデルを生成した旨が物性値算出用モデル生成部１６から通知されると、物性値算出用モデル生成部１６によって作成されたパラメータ算出用の全原子モデルを用いて分子動力学法によるシミュレーションを行う物性値算出用シミュレーション手段である。物性値算出用シミュレーション部１７は、上記の物性値算出用の全原子モデルに対して、パラメータ算出部１５によって算出されたポテンシャルパラメータを設定して上記のシミュレーションを行う。物性値算出用シミュレーション部１７は、物性値算出用の全原子モデルを示す情報、及び当該物性値算出用の全原子モデル用のポテンシャルパラメータを示す情報を記憶部１１から読み出し、シミュレーションを行う。また、分子動力学法によるシミュレーションに用いられる原子や分子間の相互作用力を表す計算式等を、物性値算出用シミュレーション部１７に予め設定しておく。
【００３３】
物性値算出用の全原子モデルのシミュレーションにおける時間間隔は、短いほど精度が高くなり好ましいが、計算に要する時間が長くなるため、パラメータ算出用の全原子モデルのシミュレーションと同様に０．１〜１０ｆｓ程度とするのがよい。例えば、時間間隔が１ｆｓの値に設定される（即ち、１時間単位が１ｆｓとされる）場合、計算時間は特に制限はないが、１００ｐｓ〜１０ｎｓ（ナノ秒）（１０５〜１０７時間単位）が好ましい。計算時間が短すぎると、算出する高分子の物性の精度が悪化し、計算時間が長すぎると、実時間も長くなるという観点から好ましくない。上記の時間間隔及び繰り返しの回数の情報については、物性値算出用シミュレーション部１７が予め記憶しており、当該記憶された情報に基づいてシミュレーションが行われる。
【００３４】
この物性値算出用の全原子モデルのシミュレーションにおける計算時間は、上述したパラメータ算出用の全原子モデルのシミュレーションにおける計算時間よりも長いものである。物性値算出用の全原子モデルのシミュレーションは、物性値を算出することを目的としているためである。なお、上記の分子動力学法によるシミュレーション自体は、従来の方法と同様に行われる。物性値算出用シミュレーション部１７は、実行したシミュレーションの結果を示す情報を記憶部１１に記憶させる。物性値算出用シミュレーション部１７は、シミュレーションの実行が完了するとその旨を物性値算出部１８に通知する。
【００３５】
物性値算出部１８は、物性値算出用シミュレーション部１７によって行われたシミュレーションの結果に基づいて、物性値を算出して出力する物性値算出手段である。物性値算出部１８は、物性値算出用シミュレーション部１７によるシミュレーションの結果を示す情報を、記憶部１１から読み出して物性値の算出に利用する。物性値算出部１８は、算出する物性値毎に予め記憶している演算式を用いて、シミュレーションの結果に対して演算を行うことによって物性値を算出する。物性値算出部１８は、算出した物性値を示す情報を、例えば、物性値計算装置１０が備える表示装置に表示することによって出力する。なお、算出した物性値を示す情報の出力は、必ずしも表示が伴うものでなくてもよく、例えば、物性値を示す情報が別途用いられるように、自装置１０の記憶部１１や物性値計算装置１０に接続された別の装置に出力することとしてもよい。以上が、本実施形態に係る物性値計算装置１０の構成である。
【００３６】
引き続いて、図２のフローチャートを用いて、物性値計算装置１０により実行される処理である、本実施形態に係る物性値計算方法について説明する。
【００３７】
物性値計算装置１０では、まず、入力部１２によって、高分子の物性値の計算に必要な上述した情報が入力される（Ｓ０１、入力ステップ）。続いて、パラメータ算出用モデル生成部１３によって、入力部１２によって入力された情報を基にパラメータ算出用の全原子モデルが生成される。（Ｓ０２、パラメータ算出用モデル生成ステップ）。続いて、パラメータ算出用シミュレーション部１４によって、パラメータ算出用モデル生成部１３によって生成されたパラメータ算出用の全原子モデルが用いられて、入力部１２によって入力された情報に係るポテンシャルパラメータが設定されて、分子動力学法によるシミュレーションが行われる（Ｓ０３、パラメータ算出用シミュレーションステップ）。
【００３８】
続いて、パラメータ算出部１５によって、パラメータ算出用シミュレーション部１４によって行われたシミュレーションの結果により示される原子の分布に基づいて、以下のように物性値算出用の全原子モデル用のポテンシャルパラメータが算出される。シミュレーションに用いられる、原子や分子間の相互作用力を計算する計算式、及びその他のポテンシャルパラメータには、一般によく知られているもの（例えば川添良幸・三上益弘・大野かおる著『コンピュータ・シミュレーションによる物質化学』共立出版、ｐｐ．５５〜８２に記載されているもの）を用いることができる。分子動力学シミュレーションでは、結合力ポテンシャルと非結合ポテンシャルとを用いることができる。結合ポテンシャルは、結合距離を平衡値に保つ分子間の結合長ポテンシャル、結合角を平衡角に保つ結合角ポテンシャル、及び取り得る２面角を制御するトーション・ポテンシャルを含むことができる。非結合ポテンシャルは、分子内においてレナード・ジョーンズポテンシャルを用いることができる。シミュレーションの結果により得られる原子の相対位置から原子間の結合長、結合角、２面角、分子間距離を求める。分子間距離から非共有ポテンシャルを求めることができる。
【００３９】
まず、シミュレーションの結果により得られる原子の相対位置から原子間の間隔（長さ）や結合角等の分布が算出される（Ｓ０４、パラメータ算出ステップ）。続いて、予め設定された分布に対する、算出した分布のフィッティングが行われる（Ｓ０５、パラメータ算出ステップ）。続いて、フィッティングした結果がポテンシャル関数に変換される（Ｓ０６、パラメータ算出ステップ）。得られたポテンシャル関数から、物性値算出用の全原子モデル用のポテンシャルパラメータが算出される（Ｓ０７、パラメータ算出ステップ）。
【００４０】
続いて、物性値算出用モデル生成部１６によって、入力部１２によって入力された情報を基に物性値算出用の全原子モデルが生成される（Ｓ０８、物性値算出用モデル生成ステップ）。続いて、物性値算出用シミュレーション部１７によって、物性値算出用モデル生成部１６によって生成された物性値算出用の全原子モデルが用いられて、パラメータ算出部１５によって算出されたポテンシャルパラメータが設定されて、分子動力学法によるシミュレーションが行われる（Ｓ０９、物性値算出用シミュレーションステップ）。続いて、物性値算出部１８によって、物性値算出用シミュレーション部１７によって行われたシミュレーションの結果に基づいて、物性値が算出されて出力される（Ｓ１０、物性値算出ステップ）。以上が、本実施形態に係る物性値計算方法である。
【００４１】
上述したように、本実施形態では、物性値を算出するために行われる物性値算出用の全原子モデルを用いたシミュレーションのためのポテンシャルパラメータが、物性値の算出対象の高分子に応じた分子の、パラメータ算出用の全原子モデルを用いたシミュレーションの結果に基づいて算出される。このように、物性値算出用の全原子モデルのポテンシャルパラメータを算出すると、当該ポテンシャルパラメータを精度のよいものにすることができ、その結果、高分子の物性値の算出を高精度に行うことができる。その一方で、物性値算出用の全原子モデルのポテンシャルパラメータは、パラメータ算出用の全原子モデルを用いた、分子動力学法によるシミュレーションによって算出されるため、量子化学計算等と比較して短時間で算出される。即ち、高分子の物性値の算出を効率的に行うことができる。
【００４２】
また、本実施形態のように、物性値の算出対象の高分子の規模を小さくした分子の、パラメータ算出用の全原子モデルを用いて、物性値算出用の全原子モデルのポテンシャルパラメータを算出するためのシミュレーションを行うこととすれば、更に短時間で上記のポテンシャルパラメータが算出される。即ち、高分子の物性値の算出を更に効率的に行うことができる。但し、必ずしも、規模を小さくした分子の、パラメータ算出用の全原子モデルを用いる必要はなく、物性値の算出対象の高分子自体の全原子モデルを用いることとしてもよい。
【００４３】
また、本実施形態のようにパラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、物性値算出用の全原子モデルを用いたシミュレーションの計算時間より少なくすることとすれば、物性値算出用の全原子モデルのポテンシャルパラメータが更に効率的に算出されるため、これにより、高分子の物性値の算出を更に効率的に行うことができる。但し、必ずしも、上記の繰り返しの回数を少なくする必要はない。
【符号の説明】
【００４４】
１０…物性値計算装置、１１…記憶部、１２…入力部、１３…パラメータ算出用モデル生成部、１４…パラメータ算出用シミュレーション部、１５…パラメータ算出部、１６…物性値算出用モデル生成部、１７…物性値算出用シミュレーション部、１８…物性値算出部、２０…外部装置。
【特許請求の範囲】
【請求項１】
コンピュータにより構成され、高分子の物性値を算出する物性値計算装置によって実行される高分子の物性値計算方法であって、
前記高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び前記高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力ステップと、
前記入力ステップにおいて入力された情報に基づいて、前記パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成ステップと、
前記パラメータ算出用モデル生成ステップにおいて作成された前記パラメータ算出用の全原子モデルを用いて、前記入力ステップにおいて入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーションステップと、
前記パラメータ算出用シミュレーションステップにおいて行われたシミュレーションの結果により示される原子の分布に基づいて、前記物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出ステップと、
前記入力ステップにおいて入力された情報に基づいて、前記物性値算出用の全原子モデルを生成する物性値算出用モデル生成ステップと、
前記物性値算出用モデル生成ステップにおいて作成された前記物性値算出用の全原子モデルを用いて、前記パラメータ算出ステップにおいて算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーションステップと、
前記物性値算出用シミュレーションステップにおいて行われたシミュレーションの結果に基づいて、前記物性値を算出して出力する物性値算出ステップと、
を含む高分子の物性値計算方法。
【請求項２】
前記パラメータ算出用シミュレーションステップにおいて行われる前記パラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、前記物性値算出用シミュレーションステップにおいて行われる前記物性値算出用の全原子モデルを用いたシミュレーションの計算時間より少ないことを特徴とする請求項１に記載の物性値計算方法。
【請求項３】
前記入力ステップにおいて、前記高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報として、前記高分子の規模を小さくした分子の全原子モデルを生成するための情報を入力することを特徴とする請求項１又は２に記載の物性値計算方法。
【請求項４】
コンピュータにより構成され、高分子の物性値を算出する物性値計算装置であって、
前記高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び前記高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力手段と、
前記入力手段によって入力された情報に基づいて、前記パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成手段と、
前記パラメータ算出用モデル生成手段によって作成された前記パラメータ算出用の全原子モデルを用いて、前記入力手段によって入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーション手段と、
前記パラメータ算出用シミュレーション手段によって行われたシミュレーションの結果により示される原子の分布に基づいて、前記物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出手段と、
前記入力手段によって入力された情報に基づいて、前記物性値算出用の全原子モデルを生成する物性値算出用モデル生成手段と、
前記物性値算出用モデル生成手段によって作成された前記物性値算出用の全原子モデルを用いて、前記パラメータ算出手段によって算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーション手段と、
前記物性値算出用シミュレーション手段によって行われたシミュレーションの結果に基づいて、前記物性値を算出して出力する物性値算出手段と、
を備える高分子の物性値計算装置。
【請求項１】
コンピュータにより構成され、高分子の物性値を算出する物性値計算装置によって実行される高分子の物性値計算方法であって、
前記高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び前記高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力ステップと、
前記入力ステップにおいて入力された情報に基づいて、前記パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成ステップと、
前記パラメータ算出用モデル生成ステップにおいて作成された前記パラメータ算出用の全原子モデルを用いて、前記入力ステップにおいて入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーションステップと、
前記パラメータ算出用シミュレーションステップにおいて行われたシミュレーションの結果により示される原子の分布に基づいて、前記物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出ステップと、
前記入力ステップにおいて入力された情報に基づいて、前記物性値算出用の全原子モデルを生成する物性値算出用モデル生成ステップと、
前記物性値算出用モデル生成ステップにおいて作成された前記物性値算出用の全原子モデルを用いて、前記パラメータ算出ステップにおいて算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーションステップと、
前記物性値算出用シミュレーションステップにおいて行われたシミュレーションの結果に基づいて、前記物性値を算出して出力する物性値算出ステップと、
を含む高分子の物性値計算方法。
【請求項２】
前記パラメータ算出用シミュレーションステップにおいて行われる前記パラメータ算出用の全原子モデルを用いたシミュレーションにおける計算時間は、前記物性値算出用シミュレーションステップにおいて行われる前記物性値算出用の全原子モデルを用いたシミュレーションの計算時間より少ないことを特徴とする請求項１に記載の物性値計算方法。
【請求項３】
前記入力ステップにおいて、前記高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報として、前記高分子の規模を小さくした分子の全原子モデルを生成するための情報を入力することを特徴とする請求項１又は２に記載の物性値計算方法。
【請求項４】
コンピュータにより構成され、高分子の物性値を算出する物性値計算装置であって、
前記高分子に応じた分子の、パラメータ算出用の全原子モデルを生成するための情報、当該パラメータ算出用の全原子モデル用のポテンシャルパラメータを示す情報、及び前記高分子の物性値算出用の全原子モデルを生成するための情報を入力する入力手段と、
前記入力手段によって入力された情報に基づいて、前記パラメータ算出用の全原子モデルを生成するパラメータ算出用モデル生成手段と、
前記パラメータ算出用モデル生成手段によって作成された前記パラメータ算出用の全原子モデルを用いて、前記入力手段によって入力された情報により示されるポテンシャルパラメータにて、分子動力学法によるシミュレーションを行うパラメータ算出用シミュレーション手段と、
前記パラメータ算出用シミュレーション手段によって行われたシミュレーションの結果により示される原子の分布に基づいて、前記物性値算出用の全原子モデル用のポテンシャルパラメータを算出するパラメータ算出手段と、
前記入力手段によって入力された情報に基づいて、前記物性値算出用の全原子モデルを生成する物性値算出用モデル生成手段と、
前記物性値算出用モデル生成手段によって作成された前記物性値算出用の全原子モデルを用いて、前記パラメータ算出手段によって算出されたポテンシャルパラメータにて、分子動力学法によるシミュレーションを行う物性値算出用シミュレーション手段と、
前記物性値算出用シミュレーション手段によって行われたシミュレーションの結果に基づいて、前記物性値を算出して出力する物性値算出手段と、
を備える高分子の物性値計算装置。
【図１】
【図２】
【図２】
【公開番号】特開２０１１−３８９７２（Ｐ２０１１−３８９７２Ａ）
【公開日】平成２３年２月２４日（２０１１．２．２４）
【国際特許分類】
【出願番号】特願２００９−１８８４８９（Ｐ２００９−１８８４８９）
【出願日】平成２１年８月１７日（２００９．８．１７）
【出願人】（０００００２０９３）住友化学株式会社 (8,981)
【公開日】平成２３年２月２４日（２０１１．２．２４）
【国際特許分類】
【出願日】平成２１年８月１７日（２００９．８．１７）
【出願人】（０００００２０９３）住友化学株式会社 (8,981)
[ Back to top ]