(57)【要約】【課題】  所望の物性値条件を満たす適合度の高い高分子材料を、短い時間であるいは労力を省いて効率的に見つけることができる高分子材料設計方法を提供する。【解決手段】  この高分子材料設計方法は、計算機を用いて所望の物性値条件を満たす高分子の構造を求めるものであって、高分子モデルの作成条件および物性値条件を入力するステップと、この作成条件の下でモノマーに基づく構成単位を高分子の設計単位とし、この設計単位を遺伝子として遺伝的アルゴリズムの個体となる高分子モデルを作成するステップと、作成された高分子モデルについて物性値を計算し、前記物性値条件に対する適合度を計算するステップとを有する。そして、遺伝的アルゴリズムを用いて、次世代以降の個体の作成および適合度計算のステップを繰り返し、適合度の大きい高分子の構造を探索することを特徴とする。【選択図】      図１
2007031541.tif 000002
【特許請求の範囲】
【請求項１】
計算機を用いて、所望の物性値条件を満たす高分子の構造を求める高分子材料設計方法であって、
高分子モデルの作成条件および物性値条件を入力するステップと、
前記作成条件の下でモノマーに基づく構成単位を高分子の設計単位とし、該設計単位を遺伝子として遺伝的アルゴリズムの個体となる高分子モデルを作成するステップと、
作成された高分子モデルについて物性値を計算し、前記物性値条件に対する適合度を計算するステップとを有し、
遺伝的アルゴリズムを用いて、次世代以降の個体の作成および適合度計算のステップを繰り返し、適合度の大きい高分子の構造を探索することを特徴とする高分子材料設計方法。
【請求項２】
前記高分子モデルの作成条件に、高分子モデルの原型の構造、高分子モデル当たりに配列される設計単位の数、高分子モデルを構成する設計単位の種類、の少なくとも１つを含むことを特徴とする請求項１に記載の高分子材料設計方法。
【請求項３】
立体配置上の異性体を有する構成単位を、該異性体毎に異なる設計単位として高分子モデルを作成することを特徴とする請求項１または請求項２に記載の高分子材料設計方法。
【請求項４】
前記遺伝的アルゴリズムにおける適合度は、所望の物性値条件を満たすほど大きな値を示す評価関数として定義されることを特徴とする請求項１ないし請求項３のいずれかに記載の高分子材料設計方法。
【請求項５】
前記物性値は、前記高分子モデルについて、分子動力学計算、あるいは分子動力学計算に分子力学計算、分子軌道計算、モンテカルロ計算およびバンド計算のうち少なくとも１つまたは複数を用いて求められることを特徴とする請求項１ないし請求項４のいずれかに記載の高分子材料設計方法。
【請求項６】
計算機を用いて、所望の物性値条件を満たす高分子の構造を求める高分子材料設計装置において、
高分子モデルの作成条件および物性値条件等を入力する装置と、
前記作成条件の下で設計単位を遺伝子として遺伝的アルゴリズムの個体となる高分子モデルを作成する装置と、
作成された高分子モデルについて物性値を計算する装置と、
物性値から適合度を計算し、遺伝的アルゴリズムを用いて適合度の大きい高分子の構造を探索する装置と、
探索結果を表示する装置と
を備えることを特徴とする高分子材料設計装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は高分子材料設計方法に関し、特に、遺伝的アルゴリズムを用いて高分子の構造を設計する高分子材料設計方法および装置に関する。
【背景技術】
【０００２】
高分子材料は、多種のものが、多種の用途に使用されている。従来、その物性値は、実際に製造したものに試験を行って求めていた。近年、計算機を用いた分子シミュレーションによって高分子材料の物性値を予測する技術が進んでいる。これにより、例えば特許文献１，２に示すように、高分子材料を実際に製造したり、材料試験を行ったりすることなく、高分子材料の評価を行うことができるようになり、設計・製造のコストや時間を低減させることができる。
【０００３】
そこで、このように精度が向上した分子シミュレーションを用いて、所定の要求性能を満たすような高分子材料を発見する、あるいは高分子材料を設計することが試みられている。これは、例えば、蓄積された過去のデータ等を基に、人間が試行錯誤的に高分子の構造を変えて選択した高分子モデルについて分子シミュレーションを用いて物性値を予測し、最も良いものを選ぶ、というような方法である。
【０００４】
しかしながら、このような方法では、高分子モデルを作成し、選択するプロセスに多くの手間や時間が掛かる。そして、この選択するプロセスにバイアスが有れば、算出される結果も適性が小さいものとなってしまう。つまり、最適の物性値をもつ高分子材料を見逃してしまう危険性もある。
【０００５】
【特許文献１】特開平９−１２７０３０号公報
【特許文献２】特開平１１−２８１６４２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
そこで、本発明の目的は、所望の物性値条件を満たす適合度の高い高分子材料を、短い時間であるいは労力を省いて効率的に見つけることができる高分子材料設計方法および装置を提供することである。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、請求項１に記載の高分子材料設計方法は、計算機を用いて所望の物性値条件を満たす高分子の構造を求める高分子材料設計方法であって、高分子モデルの作成条件および物性値条件を入力するステップと、前記作成条件の下でモノマーに基づく構成単位を高分子の設計単位とし、該設計単位を遺伝子として遺伝的アルゴリズムの個体となる高分子モデルを作成するステップと、作成された高分子モデルについて物性値を計算し、前記物性値条件に対する適合度を計算するステップとを有し、遺伝的アルゴリズムを用いて、次世代以降の個体の作成および適合度計算のステップを繰り返し、適合度の大きい高分子の構造を探索することを特徴とする。
【０００８】
この発明においては、入力された高分子モデルの作成条件の範囲内において、遺伝的アルゴリズムを用いて個体の作成および適合度計算のステップが繰り返し実行され、その結果、適合度の大きい高分子の構造が探索される。遺伝的アルゴリズムという、客観的で効率的な手法を用いることにより、少なくとも入力された高分子モデルの作成条件の範囲内において、最適な高分子構造が算出される。なお、ここで、「高分子モデル」とは、複数の設計単位が特定の配列をなしたもので、それに基づいて分子シミュレーションで物性値を算出することができるものを指す。
【０００９】
この発明において、「入力される高分子モデルの作成条件」は、解の最適度とそれを得るために要する時間やコストを左右する重要な因子である。「入力される高分子モデルの作成条件」が厳密であれば、解の最適度は高いが、それに要する計算量、つまりコンピュータを占有する時間やそのコストは大きくなる。このような問題は、分子シミュレーションという技術自体に既にあるものであるが、この発明においても引き継がれるものである。将来的には、コンピュータ技術のさらなる進歩によって時間やコストはより小さくなり、この問題の比重も小さくなる。
【００１０】
請求項２に記載の高分子材料設計方法は、請求項１に記載の発明において、前記高分子モデルの作成条件に、高分子モデルの原型の構造、高分子モデル当たりに配列される設計単位の数、高分子モデルを構成する設計単位の種類、の少なくとも１つを含むことを特徴とする。
【００１１】
高分子モデルを作成することは、上述したように、分子シミュレーションの手法自体の前提であるが、特に本発明においては、高分子モデルの作成条件を適当に設定して効率的に近似解を求めることが必要である。遺伝的アルゴリズムの機能は、基本的に同じ分岐構造・サイズの型の中で行われる。同じ分岐構造・サイズとは、高分子鎖の分岐構造が同じであり、高分子鎖１本当たりの全設計単位数および分岐している部分鎖それぞれの設計単位数が同じことを指す。分岐構造・サイズが異なれば、別途に遺伝的アルゴリズムを用いてその分岐構造・サイズにおける解を求め、それを比較することになる。
【００１２】
図３に示すように、原型として種々の骨組みが考えられ、（ａ）に示す１次元だけでなく、（ｂ）、（ｃ）のような２〜３次元的な構造も考えられる。しかしながら、多次元構造の高分子についても、分岐構造・サイズを規定することで、１次元構造の高分子と同じに扱うことができるので、高分子モデルの型は１次元つまり直線構造で考えれば充分である。
【００１３】
実際の高分子材料は、設計単位数が数千〜数万におよぶ長さであり、そのような高分子鎖が多数、相互に絡み合ってできている。一方、分子シミュレーションでは、従来より、計算時間の都合上、数を減らすことが必要である。高分子鎖１本あたりの設計単位数が少なすぎると、高分子鎖としての特徴が表現されないので、設計単位数は多い方が計算の精度上好ましいが、設計単位数が多すぎると、計算時間が長くなり、好ましくない。また、高分子鎖数が少ないと、高分子鎖同士の影響が表現されないので、高分子鎖数は多い方が計算の精度上好ましいが、高分子鎖数が多すぎると、計算時間が長くなり、好ましくない。そこで、高分子モデルとして数十〜数百の設計単位数から成る高分子鎖を数本ほど作成し、これを計算することが一般的に行われている。
【００１４】
本発明においても、これらの考え方はそのまま適用でき、原型のサイズ、すなわち、設計単位が配置されるサイトの数は５０〜５００程度、高分子鎖数ｍとしては２〜１０程度が好ましい。設計単位数と高分子鎖数の具体的な決定方法に、特に制限はない。決定方法としては、論文や学会等で発表されている計算事例を参照して決定する方法、設計単位数と高分子鎖数をそれぞれ変えて作成した高分子モデルについて実際に計算し、計算結果の各物性値を比較して決定する方法、および物性値の測定結果があれば、計算結果と測定結果を比較して決定する方法等が挙げられる。
【００１５】
高分子モデルを構成する設計単位は、予め採用可能なものを設計単位データベースに登録しておき、適宜に選択すればよい。設計単位の種類は、多くのものを選択すれば多種類の構造を作成することができるので、多種類の高分子モデルを探索することができる。その結果、他の条件と同様に、最適解に近づけるのに時間が掛かる、という相反関係にある。
【００１６】
請求項３に記載の高分子材料設計方法は、請求項１または請求項２に記載の発明において、立体配置上の異性体を有する構成単位を、該異性体毎に異なる設計単位として高分子モデルを作成することを特徴とする。
【００１７】
立体配置上の異性体とは、置換基の立体配置が異なる構成単位、配列の向きが異なる構成単位等、配列した際に立体配置が異なる高分子鎖が生成する構成単位である。モノマーに基づく構成単位には、分子式に加えて、置換基の立体配置、配列の向き等の立体配置の情報を持たせ、これらの情報が異なれば異なる設計単位として定義する。例えば、プロピレンに基づく構成単位ならば、構成単位を配置する際の向きを示す情報を、「Ｆｒｏｎｔ」および「Ｂａｃｋ」と記すとして、図４に示すように４個の異なる設計単位として定義し、プロピレン−（１）、プロピレン−（２）、プロピレン−（３）、プロピレン−（４）と名付ける。そして、これらの定義した設計単位を設計単位データベースに登録する。設計単位データベースは、図４のように設計単位名と設計単位の定義をデータとして保有する。設計単位データベースに登録される設計単位としては、エチレン−（１）、プロピレン−（１）、プロピレン−（２）、プロピレン−（３）、プロピレン−（４）、１−ブテン−（１）、１−ブテン−（２）、１−ブテン−（３）、１−ブテン−（４）、１−ペンテン−（１）、１−ペンテン−（２）、１−ペンテン−（３）、１−ペンテン−（４）、１−ヘキセン−（１）、１−ヘキセン−（２）、１−ヘキセン−（３）、１−ヘキセン−（４）、１−オクテン−（１）、１−オクテン−（２）、１−オクテン−（３）、１−オクテン−（４）、スチレン−（１）、スチレン−（２）、スチレン−（３）、およびスチレン−（４）等が挙げられる。なお、上記のプロピレン−（１）等の「−（１）」という表現は、高分子化学で使われている表現ではなく、異なる設計単位として定義するために、本明細書で独自に用いている表現である。
【００１８】
請求項４に記載の高分子材料設計方法は、請求項１ないし請求項３のいずれかに記載の発明において、前記遺伝的アルゴリズムにおける適合度は、所望の物性値条件を満たすほど大きな値を示す評価関数として定義されることを特徴とする。最適化の過程において、評価関数の最大となるような解を探索させるためである。
【００１９】
請求項５に記載の高分子材料設計装置は、請求項１ないし請求項４のいずれかに記載の発明において、前記物性値は、前記高分子モデルについて、分子動力学計算、あるいは分子動力学計算に分子力学計算、分子軌道計算、モンテカルロ計算およびバンド計算のうち少なくとも１つを組み合わせて求められることを特徴とする。
【００２０】
請求項６に記載の高分子材料設計方法は、計算機を用いて所望の物性値条件を満たす高分子の構造を求める高分子材料設計装置において、高分子モデルの作成条件や物性値条件等を入力する装置と、高分子モデルを作成する装置と、高分子モデルについて物性値を計算する装置と、物性値から適合度を計算し、遺伝的アルゴリズムを用いて適合度の大きい高分子の構造を探索する装置と、探索結果を表示する装置とを備えることを特徴とする。
【発明の効果】
【００２１】
請求項１ないし請求項６に記載の発明によれば、所望の物性値条件を満たす適合度の高い高分子材料を、短い時間であるいは労力を省いて効率的に見つけることができる。
【発明を実施するための最良の形態】
【００２２】
以下に、この発明の実施の形態について図面を参照して説明する。
まず、以下に、高分子材料設計装置について図１を用いて説明する。
高分子材料設計装置は、高分子モデルの作成条件や物性値条件等を入力する入力装置１０と、入力された高分子モデルの作成条件に基づいて、設計単位を構成単位とする高分子モデルを作成する高分子モデル作成装置２０と、高分子モデルについて物性値を計算する物性値計算装置３０と、物性値から適合度を計算し、遺伝的アルゴリズムを用いて適合度の大きい高分子の構造を探索する高分子構造探索装置４０と、および高分子構造探索装置で探索した結果を表示する表示装置５０から構成されている。
【００２３】
高分子モデル作成装置２０は、高分子モデル作成部２１、入力部２２、乱数発生部２３、および設計単位データベース２４を有する。また、物性値計算装置３０は、データ変換部３１、物性値計算部３２を有する。高分子構造探索装置４０は、適合度計算部４１、高分子モデル選択部４２、交叉発生部４３を有する。
なお、これらの各装置やそれを構成する各部は、コンピュータがこれらの機能を行なう部分を仮想的に表すものであって、ハードウエアとしての具体的な構成が限定されるものではない。
【００２４】
次に、１つの実施の形態の高分子材料設計方法について、図２のフロー図を用いて説明する。
まず、ステップ（１ａ）において、高分子モデルの作成条件を入力する。高分子モデルの作成条件としては、この実施の形態では、高分子モデルの原型、高分子モデル当たりに配列される設計単位の数、遺伝子として用いる設計単位の種類が挙げられる。
高分子モデルは、設計単位を構成単位として作成され、設計単位を複数個つなげることで高分子鎖を作成する。図３に示すように、原型として種々の骨組みが考えられるが、ここでは、２次元以上とする実益が無いので、１次元つまり直線構造を採用した。
【００２５】
次に、高分子鎖１本あたりの設計単位数ｎと、高分子鎖数ｍを入力する。設計単位数ｎに、特に制限はないが、高分子鎖１本あたりの設計単位数が少なすぎると、高分子鎖としての特徴が表現されないので、設計単位数は多い方が計算の精度上好ましい。従って、使用するコンピュータの性能や、許容される計算時間等を考慮して、可能な限り大きく設定する。高分子鎖数ｍに、特に制限はない。ただし、高分子鎖数が少ないと、高分子鎖同士の影響が表現されないので、高分子鎖数は多い方が計算の精度上好ましい。従って、使用するコンピュータの性能や、許容される計算時間等を考慮して、可能な限り大きく設定する。現状のコンピュータの性能では、設計単位数ｎとしては５０〜５００程度、高分子鎖数ｍとしては２〜１０程度が好ましい。
【００２６】
次に、ステップ（１ｂ）において、物性値条件を入力する。
高分子材料に求められる物性値を選択し、物性値条件を表現する適合度評価関数を定義し、入力する。物性値としては、密度、凝集エネルギー密度、溶解度パラメータ、ガラス転移温度、融点、体積弾性率、体積膨張率、定圧比熱、定積比熱、屈折率、固有複屈折率、配向複屈折率等が挙げられる。通常、これらの中から目的に応じた１つの物性値（Ｘ）を選択するが、場合によっては、複数の物性値（Ｘ，Ｙ，・・）を選択する。そして、選択した物性値条件を表現する適合度評価関数Ｆ＝ｆ（Ｘ）を定義し、入力する。
【００２７】
適合度評価関数Ｆは、物性値が、高分子材料に求められる条件を満たすほど大きくなるように定義する。例えば、探索する高分子材料の条件として、或る物性値Ｘが大きい程良いのであれば、
Ｆ＝Ｘ
と定義する。物性値Ｘが特定の値Ａに近い程良いのであれば、
Ｆ＝１／（Ｘ−Ａ）２
と定義する。
物性値が複数の場合には、適合度評価関数Ｆもより複雑になる。例えば、
Ｆ＝ＸａＹｂ
Ｆ＝１／｛ａ（Ｘ−Ａ）２＋ｂ（Ｙ−Ｂ）２｝
等となる。これらの式で、ａ，ｂは、複数の物性値Ｘ、Ｙに重み付けを与えるための定数である。
【００２８】
次に、ステップ（２ａ）において、高分子モデルを構成する設計単位を、設計単位データベースに登録されたものの中から選択する。登録された設計単位は、立体規則性高分子や、Ｈｅａｄ  ｔｏ  Ｔａｉｌ型構造、Ｈｅａｄ  ｔｏ  Ｈｅａｄ型構造も表現するため、置換基の立体配置、配列の向き等の立体配置の情報をそれぞれ持ち、これらの情報が異なれば異なる設計単位として定義する。例えば、プロピレンに基づく構成単位ならば、図４に示すように４個の異なる設計単位として定義し、プロピレン−（１）、プロピレン−（２）、プロピレン−（３）、プロピレン−（４）と名付ける。そして、これら定義した設計単位を設計単位データベースに登録しておく。
【００２９】
ここでは、エチレン、プロピレン、１−ブテン、１−ペンテン、および１−ヘキセンを共重合体した共重合体高分子を考える。このとき求められている或る物性値条件を満たす高分子材料として、エチレンに基づく構成単位、プロピレンに基づく構成単位、１−ブテンに基づく構成単位、１−ペンテンに基づく構成単位、および１−ヘキセンに基づく構成単位をそれぞれどれだけの比率で含有し、どのような順番に並んでいる構造がよいかを知りたいとする。その場合、選択する設計単位としては、エチレン−（１）、プロピレン−（１）、１−ブテン−（１）、１−ペンテン−（１）、および１−ヘキセン−（１）の５種類の設計単位が必要である。
【００３０】
あるいは、プロピレン−（１）の換わりにプロピレン−（２）、プロピレン−（３）、プロピレン−（４）のいずれかでもよく、１−ブテン−（１）の換わりに１−ブテン−（２）、１−ブテン−（３）、１−ブテン−（４）のいずれかでもよく、１−ペンテン−（１）の換わりに１−ペンテン−（２）、１−ペンテン−（３）、１−ペンテン−（４）のいずれかでもよく、１−ヘキセン−（１）の換わりに１−ヘキセン−（２）、１−ヘキセン−（３）、１−ヘキセン−（４）のいずれかでもよい。
【００３１】
さらに、立体規則性（タクティシティー）が違う構造も考えるならば、選択する設計単位としては、エチレン−（１）、プロピレン−（１）、１−ブテン−（１）、１−ペンテン（１）、１−ヘキセン−（１）、プロピレン−（２）、１−ブテン−（２）、１−ペンテン（２）、および１−ヘキセン−（２）の９種類の設計単位が必要である。
さらに、Ｈｅａｄ  ｔｏ  Ｔａｉｌ型構造やＨｅａｄ  ｔｏ  Ｈｅａｄ型構造も考えるならば、選択する設計単位としては、エチレン、プロピレン−（１）、１−ブテン−（１）、１−ペンテン（１）、１−ヘキセン−（１）、プロピレン−（２）、１−ブテン−（２）、１−ペンテン（２）、１−ヘキセン−（２）、プロピレン−（３）、１−ブテン−（３）、１−ペンテン（３）、１−ヘキセン−（３）、プロピレン−（４）、１−ブテン−（４）、１−ペンテン（４）、および１−ヘキセン−（４）の１７種類の設計単位が必要である。
【００３２】
次に、ステップ（２ｂ）において、遺伝的アルゴリズムにおける第１世代の個体集合となるＮ種類の高分子モデルを作成する。
高分子モデルの作成方法として、配列を使う。配列Ｓ（ｉ）が高分子鎖であり、配列要素ｉが設計単位である。
図３にあるような直鎖型の高分子鎖ならば、１次元配列Ｓ（ｉ）を使い、ステップ（２ａ）で選択した設計単位を１次元配列Ｓ（１），Ｓ（２），…，Ｓ（ｉ），…，Ｓ（ｎ）に入れる。これをＮ回繰り返し、図５のように、高分子モデルをＮ種類作成する。
【００３３】
配列への設計単位の入力方法としては、各番号の配列に入れる設計単位を、入力部２２を使ってユーザーが決めて入力する方法と、各番号の配列に入れる設計単位を、乱数発生部２３を使って乱数発生によって決めて入力する方法が挙げられる。ユーザーが決めて入力する方法では、既存の高分子材料や想定されている高分子材料の高分子モデルを作成することができ、乱数発生によって決めて入力する方法は、多様な高分子モデルを作成することに適している。
【００３４】
「Ｆｒｏｎｔ」および「Ｂａｃｋ」を接続点として、配列に入った設計単位を順次つなげる。ｎ（ｉ＝１，２，…，ｎ）個のうちのｉ番目の設計単位の「Ｂａｃｋ」と、ｉ＋１番目の設計単位の「Ｆｒｏｎｔ」を順次接続する。接続した時点で、これら「Ｂａｃｋ」と「Ｆｒｏｎｔ」は削除される。１番目の設計単位の「Ｆｒｏｎｔ」を水素原子「Ｈ」に置き換え、ｎ番目の設計単位の「Ｂａｃｋ」も同様に水素原子「Ｈ」に置き換える。以上の作業により、図６のように高分子鎖１本（ｍ＝１）の高分子モデルが作成される。
【００３５】
次に、ステップ（３ａ）において、高分子モデルの構造ファイルを作成する。
ステップ（２ｂ）で作成した高分子モデルを記述した構造ファイルを作成する。構造ファイルは、高分子モデルを構成する各原子の元素名、各原子の位置を表す３次元座標、および原子の結合のデータ等を有する。構造ファイルのフォーマットとして、図７に示したようなＭＤＬ−Ｍｏｌ形式をはじめとする汎用のフォーマットが挙げられる。構造ファイルのフォーマットの選び方としては、ステップ（３ｂ）で使用する各汎用ログラムで読み込めるフォーマットでなければならない。
【００３６】
次に、ステップ（３ｂ）において、Ｎ種類の高分子モデルについて、物性値Ｘをそれぞれ計算する。
物性値としては、密度、凝集エネルギー密度、溶解度パラメータ、ガラス転移温度、融点、体積弾性率、体積膨張率、定圧比熱、定積比熱、屈折率、固有複屈折率、配向複屈折率等が挙げられる。これら物性値は、分子動力学計算、あるいは分子動力学計算に分子力学計算、分子軌道計算、モンテカルロ計算およびバンド計算のうち少なくとも１つを組み合わせて求められる。これらの計算は、例えば、財団法人化学技術戦略推進機構のＯＣＴＡ、株式会社日本総合研究所のＪ−ＯＣＴＡ、Ａｃｃｅｌｒｙｓ  Ｉｎｃ．のＭＳ  Ｍｏｄｅｌｉｎｇ  ｐｒｏｄｕｃｔｓ等（いずれも商品名）の汎用プログラムを使って行うことができる。
【００３７】
このような汎用プログラムを起動させ、ステップ（３ａ）において作成した高分子モデルの構造ファイルを読み込む。読み込んだ高分子鎖１本の高分子モデルを図示させると、図８のようになる。高分子鎖ｍ本のアモルファス（非晶）状態を作成し、図示させると、図９のようになる。高分子モデルのアモルファス状態の安定構造を計算し、物性値Ｘを求める。以上の作業および計算は、汎用プログラムの機能を使って実施する。
【００３８】
次に、ステップ（４ａ）において、Ｎ種類の高分子モデルについて、適合度評価関数Ｆをそれぞれ計算する。適合度評価関数Ｆは、ステップ（１ｂ）で定義した関数式Ｆ＝ｆ（Ｘ）によって物性値Ｘに基づいて計算される。
【００３９】
次に、ステップ（４ｂ）おいて、Ｎ種類の高分子モデルから、適合度評価関数が大きいＭ種類の高分子モデルを選択し、Ｍ種類の高分子モデルを残す。このときＮ≧Ｍである。選択方式としては、ルーレット選択方式、トーナメント選択方式、エリート選択方式、およびこれらの組み合わせを使う。選択方式の組み合わせ方法に、特に制限はない。
ルーレット選択方式は、適合度評価関数の大きさに比例した割合で高分子モデルを選択する方式である。適合度評価関数の大きさに比例した面積でＮ個の領域に分割されているルーレットを作り、投げた玉の入った領域の高分子モデルを選択するという過程をＭ回繰り返して、Ｍ種類の高分子モデルを選択する方式である。
トーナメント選択方式は、高分子モデルをある数（トーナメントサイズ）だけランダムに選び出し、その中で適合度評価関数が最大の高分子モデルを選択するという過程をＭ回繰り返して、Ｍ種類の高分子モデルを選択する方式である。
エリート選択方式は、適合度評価関数が大きい高分子モデルを上位からＭ種類選択する方式である。
【００４０】
次に、ステップ（４ｃ）において、Ｍ種類の高分子モデルについて交叉を起こし、合計でＮ種類の高分子モデルを作成する。例えば、図１０に示したように高分子モデルを構成する設計単位（遺伝的アルゴリズムにおける遺伝子）の組み替えを行い、合計でＮ種類の高分子モデルを作成する。これが次世代の高分子となる。交叉のやり方として、特に制限はない。交叉数（１点交叉、２点交叉、さらに複数点の交叉をするのか）、交叉の位置に特に制限はない。
図９に例を示す。
【００４１】
次に、ステップ（４ｄ）において、ステップ（３ａ）からステップ（４ｃ）を第２世代から第Ｇ世代まで繰り返す。そして、最後に、ステップ（５）において、第１世代から第Ｇ世代までの（Ｎ×Ｇ）種類の高分子モデルのうち、適合度評価関数Ｆの値が大きい高分子モデルを複数表示する。これにより、所望の物性値条件を満たす適合度の高い高分子材料を、短い時間であるいは労力を省いて効率的に見つけることができる。
【図面の簡単な説明】
【００４２】
【図１】本発明の高分子材料設計装置の構成を示す図である。
【図２】本発明の高分子材料設計方法のフローチャートである。
【図３】本発明の高分子材料設計方法によって設計する高分子構造を概念的に示す図である。
【図４】本発明の高分子材料設計装置が備える設計単位データベースの記録内容を概念的に示す図である。
【図５】高分子モデルを配列で表現した図である。
【図６】高分子モデルを分子構造で表現した図である。
【図７】高分子モデルのＭＤＬ−Ｍｏｌ形式ファイルの一部を示す図である。
【図８】高分子鎖１本の高分子モデルの図である。
【図９】高分子鎖ｍ本の高分子モデルのアモルファス（非晶）状態の図である。
【図１０】交叉の図である。
【符号の説明】
【００４３】
１０  入力装置
２０  高分子モデル作成装置
２１  高分子モデル作成部
２２  入力部
２３  乱数発生部
２４  設計単位データベース
３０  物性値計算装置
３１  データ変換部
３２  物性値計算部
４０  高分子構造探索装置
４１  適合度計算部
４２  高分子モデル選択部
４３  交叉発生部
５０  表示装置
【図１】
2007031541.tif 000003
【図２】
2007031541.tif 000004
【図３】
2007031541.tif 000005
【図４】
2007031541.tif 000006
【図５】
2007031541.tif 000007
【図６】
2007031541.tif 000008
【図７】
2007031541.tif 000009
【図８】
2007031541.tif 000010
【図９】
2007031541.tif 000011
【図１０】
2007031541.tif 000012