(57)【要約】【課題】動的架橋熱可塑性エラストマーの分子モデルを生成する方法を提供する。【解決手段】第１ポリマー粒子２０が連なる第１ポリマー２と、第２ポリマー粒子３０が連なる第２ポリマー３と、架橋剤粒子４と、を有する混合モデルＭ１に相互作用を設定するステップＳＴ１と、混合モデルをせん断変形させる変形処理と、第１ポリマー粒子２０と架橋剤粒子４とを結合させる架橋処理と、を第２温度Ｔ２にて同時に実行するステップＳＴ４と、を含む。第２温度よりも低い第１温度がガラス転移点より大きくなる非結合相互作用を第１ポリマー粒子２０に設定し、第１温度Ｔ１がガラス転移点以下となり、第２温度Ｔ２がガラス転移点より大きくなる非結合相互作用を第２ポリマー粒子３０に設定し、第１ポリマー粒子２０と第２ポリマー粒子３０との非結合相互作用が、第１ポリマー粒子２０同士及び第２ポリマー粒子３０同士に作用する非結合相互作用よりも大きな斥力を有する。【選択図】図２
2020012075.tif 000002
【特許請求の範囲】
【請求項１】
１又は複数のプロセッサが実行する方法であって、
複数の第１ポリマー粒子が直鎖状又は分岐状に連なる複数の第１ポリマーと、複数の第２ポリマー粒子が直鎖状又は分岐状に連なる複数の第２ポリマーと、前記第１ポリマー同士を結合するための複数の架橋剤粒子と、を有する混合モデルに対して相互作用を設定するステップと、
前記混合モデルをせん断変形させる変形処理と、互いに所定距離内にある前記第１ポリマー粒子と前記架橋剤粒子とを所定確率で結合させる架橋処理と、を第２温度において同時に実行するステップと、を含み、
前記第１ポリマー粒子には、前記第２温度よりも低い第１温度がガラス転移点より大きくなる非結合相互作用を設定し、
前記第２ポリマー粒子には、前記第１温度がガラス転移点以下となり、前記第２温度がガラス転移点より大きくなる非結合相互作用を設定し、
前記第１ポリマー粒子と前記第２ポリマー粒子との間に作用する非結合相互作用が、前記第１ポリマー粒子同士及び前記第２ポリマー粒子同士に作用する非結合相互作用よりも大きな斥力を有するように設定する、熱可塑性エラストマーの分子モデルを生成する方法。
【請求項２】
前記第２ポリマー粒子同士に作用する非結合相互作用の強さを表すパラメータは、前記第１ポリマー粒子同士に作用する非結合相互作用の強さを表すパラメータよりも大きい、請求項１に記載の方法。
【請求項３】
前記第２ポリマー粒子と前記架橋剤粒子との間に作用する非結合相互作用の粒子の大きさを表すパラメータは、前記第１ポリマー粒子と前記架橋剤粒子との間に作用する非結合相互作用の粒子の大きさを表すパラメータよりも大きい、請求項１又は２に記載の方法。
【請求項４】
前記第１ポリマー粒子と前記第２ポリマー粒子との間に作用する非結合相互作用は、斥力のみを有する、請求項１〜３のいずれかに記載の方法。
【請求項５】
複数の第１ポリマー粒子が直鎖状又は分岐状に連なる複数の第１ポリマーと、複数の第２ポリマー粒子が直鎖状又は分岐状に連なる複数の第２ポリマーと、前記第１ポリマー同士を結合するための複数の架橋剤粒子と、を有する混合モデルに対して相互作用を設定する設定部と、
前記混合モデルをせん断変形させる変形処理を実行する変形処理実行部と、
互いに所定距離内にある前記第１ポリマー粒子と前記架橋剤粒子とを所定確率で結合させる架橋処理を実行する架橋処理実行部と、を備え、
前記変形処理実行部及び前記架橋処理実行部は、前記変形処理と前記架橋処理とを第２温度において同時に実行し、
前記設定部は、前記第１ポリマー粒子には、前記第２温度よりも低い第１温度がガラス転移点より大きくなる非結合相互作用を設定し、
前記第２ポリマー粒子には、前記第１温度がガラス転移点以下となり、前記第２温度がガラス転移点より大きくなる非結合相互作用を設定し、
前記第１ポリマー粒子と前記第２ポリマー粒子との間に作用する非結合相互作用は、前記第１ポリマー粒子同士及び前記第２ポリマー粒子同士に作用する非結合相互作用よりも大きな斥力を有するように設定する、熱可塑性エラストマーの分子モデルを生成するシステム。
【請求項６】
前記第２ポリマー粒子同士に作用する非結合相互作用の強さを表すパラメータは、前記第１ポリマー粒子同士に作用する非結合相互作用の強さを表すパラメータよりも大きい、請求項５に記載のシステム。
【請求項７】
前記第２ポリマー粒子と前記架橋剤粒子との間に作用する非結合相互作用の粒子の大きさを表すパラメータは、前記第１ポリマー粒子と前記架橋剤粒子との間に作用する非結合相互作用の粒子の大きさを表すパラメータよりも大きい、請求項５又は６に記載のシステム。
【請求項８】
前記第１ポリマー粒子と前記第２ポリマー粒子との間に作用する非結合相互作用は、斥力のみを有する、請求項５〜７のいずれかに記載のシステム。
【請求項９】
請求項１〜４のいずれかに記載の方法を１又は複数のプロセッサに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動的架橋熱可塑性エラストマーモデルを生成する方法、システム及びプログラムに関する。
【背景技術】
【０００２】
熱可塑性エラストマーは、高温で溶融して成形が容易となり、常温で固化してゴム状の性質を持つ材料であり、産業で広く用いられている。互いに異なる複数のポリマーの混合物にせん断変形を与えながら、一方のポリマーを架橋することで得られるものが、動的架橋熱可塑性エラストマーと呼ばれる
【０００３】
非特許文献１には、動的架橋熱可塑性エラストマーについての記載がある。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】動的架橋熱可塑性エラストマー、神品順二、日本ゴム協会誌、第７６巻第８号（２００３）、３１０〜３１５頁
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記文献は一般的な紹介程度にとどまり、動的架橋熱可塑性エラストマーの分子モデルを生成する方法については提案されていない。
【０００６】
動的架橋熱可塑性エラストマーの開発には、共連続相と呼ばれる相分離構造にせん断変形を加えたときの構造変化を知ることが重要と考えられるが、一般的な理論がなく、試行錯誤に多くの工数を要する。シミュレーションを用いて動的架橋熱可塑性エラストマーの力学特性を予測することが可能になれば、上記コストを軽減可能となる。
【０００７】
本発明の目的は、動的架橋熱可塑性エラストマーの分子モデルを生成する方法、システム及びプログラムを提供することである。
【課題を解決するための手段】
【０００８】
本発明の熱可塑性エラストマーの分子モデルを生成する方法は、
１又は複数のプロセッサが実行する方法であって、
複数の第１ポリマー粒子が直鎖状又は分岐状に連なる複数の第１ポリマーと、複数の第２ポリマー粒子が直鎖状又は分岐状に連なる複数の第２ポリマーと、前記第１ポリマー同士を結合するための複数の架橋剤粒子と、を有する混合モデルに対して相互作用を設定するステップと、
前記混合モデルをせん断変形させる変形処理と、互いに所定距離内にある前記第１ポリマー粒子と前記架橋剤粒子とを所定確率で結合させる架橋処理と、を第２温度において同時に実行するステップと、を含み、
前記第１ポリマー粒子には、前記第２温度よりも低い第１温度がガラス転移点より大きくなる非結合相互作用を設定し、
前記第２ポリマー粒子には、前記第１温度がガラス転移点以下となり、前記第２温度がガラス転移点より大きくなる非結合相互作用を設定し、
前記第１ポリマー粒子と前記第２ポリマー粒子との間に作用する非結合相互作用が、前記第１ポリマー粒子同士及び前記第２ポリマー粒子同士に作用する非結合相互作用よりも大きな斥力を有するように設定する。
【０００９】
このように、第１ポリマー及び第２ポリマーに非結合相互作用を設定すれば、動的架橋を行う第２温度よりも低い第１温度において第１ポリマーは溶融状態となるが第２ポリマーはガラス状態となり、第２温度では第１ポリマー及び第２ポリマーの両方とも溶融状態となり、熱可塑性エラストマーのマスターバッチが模擬可能となる。
また、第１ポリマー及び第２ポリマーが溶融状態となる第２温度においてせん断変形させつつ、第１ポリマーと架橋剤粒子とを架橋させるので、動的架橋可能となる。
さらに、第１ポリマー粒子と第２ポリマー粒子とに作用する非結合相互作用は、第１ポリマー粒子同士及び第２ポリマー粒子同士に作用する非結合相互作用よりも大きな斥力を有するように設定されているので、動的架橋熱可塑性エラストマーの特徴である第１ポリマーと第２ポリマーとの相分離を模擬可能となる。
したがって、動的架橋熱可塑性エラストマーの分子モデルを生成することができる。
【図面の簡単な説明】
【００１０】
【図１】熱可塑性エラストマーの分子モデルを生成するシステムを示すブロック図。
【図２】上記システムが実行する処理ルーチンを示すフローチャート。
【図３】複数の第１ポリマーと、複数の第２ポリマーと、複数の架橋剤粒子を有する混合モデルを示す説明図。
【図４】動的架橋熱可塑性エラストマーのモデルを示す図。
【図５】ＬＪ相互作用に関する説明図。
【図６】図４のモデルを単軸引張計算した結果を示す図。
【図７】第１温度における図４のモデルの粘弾性計算の結果を示す図。
【図８】第２温度における図４のモデルの粘弾性計算の結果を示す図。
【発明を実施するための形態】
【００１１】
以下、本発明の一実施形態を、図面を参照して説明する。
【００１２】
［熱可塑性エラストマーの分子モデルを生成するシステム］
本実施形態のシステム１は、分子動力学計算に用いられる、熱可塑性エラストマーの分子モデルを生成する。
【００１３】
図１に示すように、システム１は、設定部１０と、変形処理実行部１１と、架橋処理実行部１２と、を有する。これら各部１０〜１２は、ＣＰＵなどのプロセッサ、メモリ、各種インターフェイス等を備えたパソコン等の情報処理装置において予め記憶されている図２に示す処理ルーチンをプロセッサが実行することによりソフトウェア及びハードウェアが協働して実現される。本実施形態では、１つの装置におけるプロセッサが各部の処理を実行しているが、これに限定されない。例えば、ネットワークを用いて分散させ、複数のプロセッサが各部の処理を実行するように構成してもよい。
【００１４】
図１に示す設定部１０は、キーボードやマウス等の既知の操作部を介してユーザからの操作を受け付け、熱可塑性エラストマー（混合モデル）を構成する第１ポリマーモデル、第２ポリマーモデル、架橋剤粒子モデル、これらモデルに設定される相互作用に関するデータの設定、分子動力学計算に必要な解析条件（常温に相当する第１温度Ｔ１、動的架橋時（溶融時）の高温に相当する第２温度Ｔ２、圧力Ｐ）などの各種設定を実行し、これら設定値をメモリに記憶する。
【００１５】
熱可塑性エラストマーである混合モデルＭ１は、図３に示すモデル完成前の状態において、複数の第１ポリマー２と、複数の第２ポリマー３と、第１ポリマー２同士を結合するための複数の架橋剤粒子４と、を有する。熱可塑性エラストマーは、ゴムと熱可塑性樹脂とを溶融混錬すると同時に、せん断下でゴムを架橋することによって製造される。第１ポリマー２は、ゴムに相当し、第２ポリマー３は熱可塑性樹脂に相当する。第１ポリマー２は、図３に示すように、複数の第１ポリマー粒子２０が直鎖状又は分岐状に連なっている。第２ポリマー３は、複数の第２ポリマー粒子３０が直鎖状又は分岐状に連なっている。第１ポリマー粒子２０同士及び第２ポリマー粒子３０同士は、所定の結合相互作用によって結合されている。本実施形態では、所定の結合相互作用として、ＦＥＮＥ（Finitely extensible nonlinear elastic）−ＬＪ（レナード・ジョーンズ）型結合ポテンシャルを設定しているが、適宜変更可能である。ＦＥＮＥ−ＬＪ型結合ポテンシャルは、次の式（１）で表される。
【数１】
2020012075.tif 000003
Kは相互作用の強さを示すパラメータ、R０は結合距離の上限を示すパラメータ、εは斥力の強さを表すパラメータ、σは斥力が生じ始める距離を表すパラメータである。ｒは粒子間距離である。
【００１６】
設定部１０は、図３に示す混合モデルＭ１に対して相互作用（結合相互作用、非結合相互作用）を設定する。設定部１０は、第１ポリマー粒子２０には、第２温度Ｔ２よりも低い第１温度Ｔ１がガラス転移点より大きくなる非結合相互作用を設定する。設定部１０は、第２ポリマー粒子３０には、第１温度Ｔ１がガラス転移点以下となり、第２温度Ｔ２がガラス転移点より大きくなる非結合相互作用を設定する。このように設定することにより、常温である第１温度Ｔ１においては第１ポリマー２が溶融状態となり、第２ポリマー３がガラス状態となる。溶融混錬時の第２温度Ｔ２においては第１ポリマー２が溶融状態となり、第２ポリマー３も溶融状態となる。
【００１７】
設定部１０は、第１ポリマー粒子２０と第２ポリマー粒子３０との間に作用する非結合相互作用が、第１ポリマー粒子２０同士及び第２ポリマー粒子３０同士に作用する非結合相互作用よりも大きな斥力を有するように設定する。このように設定することにより、第１ポリマー２と第２ポリマー３とを相分離させることが可能となる。
【００１８】
非結合相互作用には、レナード・ジョーンズ相互作用（ポテンシャル）を使用している。この非結合相互作用は、次の式（２）で表される。
【数２】
2020012075.tif 000004
εは相互作用の強さを示すパラメータである。σは粒子の大きさを表すパラメータである。ｒは粒子間距離である。ｒｃは相互作用のカットオフである。
【００１９】
図１に示す変形処理実行部１１は、図３に示す混合モデルＭ１に対して第２温度Ｔ２にてせん断変形させる変形処理を実行する。具体的には、溶融状態となる第２温度Ｔ２において周期的なせん断変形を与え、分子動力学計算を実行する。
【００２０】
図１に示す架橋処理実行部１２は、第２温度Ｔ２において分子動力学計算を実行し、所定タイムステップが経過するたびに、互いに所定距離内にある第１ポリマー粒子２０と架橋剤粒子４とを所定確率で結合させる架橋処理を実行する。図３は、架橋剤粒子４が第１ポリマー２に結合していないが、架橋処理が完了すれば、架橋剤粒子４が第１ポリマー２に結合することになる。変形処理実行部１１による変形処理と架橋処理実行部１２による架橋処理とは第２温度Ｔ２で同時に実行する。全ての架橋剤粒子４が第１ポリマー粒子２０と結合すれば、変形処理及び架橋処理の実行を終了する。これにより、動的架橋熱可塑性エラストマーの分子モデルが完成する。その後、第２温度Ｔ２および第１温度Ｔ１の順序で平衡化処理を実行することが好ましい。平衡化処理では、図３に示す計算領域Ａｒ１すなわちモデルの体積がほぼ一定になる（体積変化が閾値以下になる）まで各分子の挙動を計算する。
【００２１】
図４に示す動的架橋熱可塑性エラストマーのモデルを示す図である。図４の例では、第１ポリマー粒子２０が５０個連なった第１ポリマー２を１００本、第２ポリマー粒子３０が５０個連なった第２ポリマー３を１００本、架橋剤粒子４を２００個、混合した混合モデルＭ１を生成し、第２温度Ｔ２を５．０（ＬＪ単位）にてせん断変形処理と架橋処理とを行い、架橋処理の完了後に、第１温度Ｔ１を１．０（ＬＪ単位）として平衡化処理した例である。密度は０．８５である。図中の黒色が架橋された第１ポリマー２を表し、図中の灰色が第２ポリマー３を表している。
【００２２】
図４の例の非結合相互作用について、第１ポリマー粒子２０同士の非結合相互作用は、ε＝１．０、σ＝１．０、カットオフｒｃ＝２．５としている。カットオフｒｃ＝２．５であるので、引力及び斥力の双方を有する。架橋剤粒子４同士の非結合相互作用も同じ値である。
【００２３】
第２ポリマー粒子３０同士の非結合相互作用は、ε＝３．０、σ＝１．０、カットオフｒｃ＝２．５である。第１ポリマー粒子２０同士のε＝１．０に対して、第２ポリマー粒子３０同士のε＝３．０とすることで、第２ポリマー３が第１温度Ｔ１にてガラス状態になるようにしている。
【００２４】
第１ポリマー粒子２０と架橋剤粒子４との間に作用する非結合相互作用は、ε＝３．０、σ＝１．０、カットオフｒｃ＝２．５である。ここでは引力と斥力の双方を有するが、少なくとも斥力が有ればよい。
【００２５】
第２ポリマー粒子３０と架橋剤粒子４との間に作用する非結合相互作用は、ε＝１．０、σ＝１．２、カットオフｒｃ＝２．５である。ここでは引力と斥力の双方を有するが、少なくとも斥力が有ればよい。本実施形態では、架橋反応の促進のために、σ＝１．２であるが、第１ポリマー粒子２０と架橋剤粒子４との間に作用する非結合相互作用のσ＝１．０と同じにしてもよい。
【００２６】
第１ポリマー粒子２０と第２ポリマー粒子３０との間に作用する非結合相互作用は、ε＝１．０、σ＝１．２、カットオフｒｃ＝１．１２である。ここでは斥力のみを有するが、斥力と引力の双方を有するようにしてもよい。カットオフｒｃ＝１．１２としているのは、２（１／６）≒１．１２２４６２０４であり、図５にＬＪ相互作用を示すように、斥力のみとするためである。第１ポリマー粒子２０と第２ポリマー粒子３０との間に作用する非結合相互作用のσ＝１．２が、第１ポリマー粒子２０同士及び第２ポリマー粒子３０同士の非結合相互作用のσ＝１．０により大きいため、大きな斥力を有し、相分離を実現可能となる。
【００２７】
図４に示す実施例が動的架橋熱可塑性エラストマーとしての特性を有するかを検証した。
【００２８】
図６は、生成した動的架橋熱可塑性エラストマーの分子モデルＭ１を単軸引張計算した結果である。伸長比が大きくなるに従い、応力が増していることから、ゴムの弾性を有することが理解できる。
【００２９】
図７は、常温（第１温度Ｔ１＝１．０）における粘弾性計算の結果である。四角が貯蔵弾性率を示し、三角が損失弾性率を示し、丸がｔａｎδを示す。図中の周波数が小さくなるほど、変形する速度が遅くなることを意味する。低周波数領域、すなわち変形速度が遅くなるに伴ってｔａｎδが一定になることから、ゴムの弾性を有することが理解できる。
【００３０】
図８は、高温（溶融時、第２温度Ｔ２＝５．０）における粘弾性計算の結果である。四角が貯蔵弾性率を示し、三角が損失弾性率を示し、丸がｔａｎδを示す。図中の周波数が小さくなるほど、変形する速度が遅くなることを意味する。低周波数領域、すなわち変形速度が遅くなるに伴ってｔａｎδが上昇している。ｔａｎδが上昇するということは、弾性体であったものが、引き延ばしに対して元に戻る力がなくなっていくことを示している。これにより、成形可能な熱可塑性であることが理解できる。
【００３１】
図６〜８の結果より、動的架橋熱可塑性エラストマーの分子モデルが生成できていることが理解できる。
【００３２】
［熱可塑性エラストマーの分子モデルを生成する方法］
図１に示すシステム１を用いて、熱可塑性エラストマーの分子モデルを生成する方法について、図２を用いて説明する。
【００３３】
まず、ステップＳＴ１において、設定部１０は、複数の第１ポリマー粒子が直鎖状又は分岐状に連なる複数の第１ポリマーと、複数の第２ポリマー粒子が直鎖状又は分岐状に連なる複数の第２ポリマーと、前記第１ポリマー同士を結合するための複数の架橋剤粒子と、を有する混合モデルＭ１に対して相互作用を設定する。
【００３４】
次のステップＳＴ２において、混合モデルＭ１の分子動力学計算を所定タイムステップ実行し、各粒子を分散させることが好ましい。
【００３５】
次に、ステップＳＴ４において、変形処理実行部１１による混合モデルＭ１をせん断変形させる変形処理と、架橋処理実行部１２による互いに所定距離内にある第１ポリマー粒子２０と架橋剤粒子４とを所定確率で結合させる架橋処理と、を第２温度において同時に実行する。ステップＳＴ４の処理は、全ての架橋剤粒子４が架橋完了するまで（ＳＴ３：ＹＥＳ）、繰り返し実行する。
【００３６】
次のステップＳＴ５において、分子動力学計算を実行し、平衡化することが好ましい。なお、ステップＳＴ２，ＳＴ５は省略可能である。
【００３７】
以上のように、本実施形態の熱可塑性エラストマーの分子モデルを生成する方法は、
１又は複数のプロセッサが実行する方法であって、
複数の第１ポリマー粒子２０が直鎖状又は分岐状に連なる複数の第１ポリマー２と、複数の第２ポリマー粒子３０が直鎖状又は分岐状に連なる複数の第２ポリマー３と、第１ポリマー２同士を結合するための複数の架橋剤粒子４と、を有する混合モデルＭ１に対して相互作用を設定するステップＳＴ１と、
混合モデルをせん断変形させる変形処理と、互いに所定距離内にある第１ポリマー粒子２０と架橋剤粒子４とを所定確率で結合させる架橋処理と、を第２温度Ｔ２において同時に実行するステップＳＴ４と、を含む。
第１ポリマー粒子２０には、第２温度Ｔ２よりも低い第１温度Ｔ１がガラス転移点より大きくなる非結合相互作用を設定し、
第２ポリマー粒子３０には、第１温度Ｔ１がガラス転移点以下となり、第２温度Ｔ２がガラス転移点より大きくなる非結合相互作用を設定し、
第１ポリマー粒子２０と第２ポリマー粒子３０との間に作用する非結合相互作用が、第１ポリマー粒子２０同士及び第２ポリマー粒子３０同士に作用する非結合相互作用よりも大きな斥力を有するように設定する。
【００３８】
本実施形態の熱可塑性エラストマーの分子モデルを生成するシステムは、
複数の第１ポリマー粒子２０が直鎖状又は分岐状に連なる複数の第１ポリマー２と、複数の第２ポリマー粒子３０が直鎖状又は分岐状に連なる複数の第２ポリマー３と、第１ポリマー２同士を結合するための複数の架橋剤粒子４と、を有する混合モデルＭ１に対して相互作用を設定する設定部１０と、
混合モデルＭ１をせん断変形させる変形処理を実行する変形処理実行部１１と、
互いに所定距離内にある第１ポリマー粒子２０と架橋剤粒子４とを所定確率で結合させる架橋処理を実行する架橋処理実行部１２と、を備える。
設定部１０は、第１ポリマー粒子２０には、第２温度Ｔ２よりも低い第１温度Ｔ１がガラス転移点より大きくなる非結合相互作用を設定し、
第２ポリマー粒子３０には、第１温度Ｔ１がガラス転移点以下となり、第２温度Ｔ２がガラス転移点より大きくなる非結合相互作用を設定し、
第１ポリマー粒子２０と第２ポリマー粒子３０との間に作用する非結合相互作用が、第１ポリマー粒子２０同士及び第２ポリマー粒子３０同士に作用する非結合相互作用よりも大きな斥力を有するように設定する。
【００３９】
このように、第１ポリマー２及び第２ポリマー３に非結合相互作用を設定すれば、動的架橋を行う第２温度Ｔ２よりも低い第１温度Ｔ１において第１ポリマー２は溶融状態となるが第２ポリマー３はガラス状態となり、第２温度Ｔ２では第１ポリマー２及び第２ポリマー３の両方とも溶融状態となり、熱可塑性エラストマーのマスターバッチが模擬可能となる。
まら、第１ポリマー２及び第２ポリマー３が溶融状態となる第２温度Ｔ２においてせん断変形させつつ、第１ポリマー２と架橋剤粒子４とを架橋させるので、動的架橋可能となる。
さらに、第１ポリマー粒子２０と第２ポリマー粒子３０とに作用する非結合相互作用は、第１ポリマー粒子２０同士及び第２ポリマー粒子３０同士に作用する非結合相互作用よりも大きな斥力を有するように設定されているので、動的架橋熱可塑性エラストマーの特徴である第１ポリマー２と第２ポリマー３との相分離を模擬可能となる。
したがって、動的架橋熱可塑性エラストマーの分子モデルを生成することができる。
【００４０】
本実施形態では、第２ポリマー粒子３０同士に作用する非結合相互作用の強さを表すパラメータ（ε＝３．０）は、前記第１ポリマー粒子２０同士に作用する非結合相互作用の強さを表すパラメータ（ε＝１．０）よりも大きい。
【００４１】
このようにすれば、、第１温度Ｔ１などの或る温度において、第１ポリマー２がガラス状でなく且つ第２ポリマー３がガラス状となるように、非結合相互作用を設定することが可能となる。
【００４２】
本実施形態では、第２ポリマー粒子３０と架橋剤粒子４との間に作用する非結合相互作用の粒子の大きさを表すパラメータ（σ＝１．２）は、第１ポリマー粒子２０と架橋剤粒子４との間に作用する非結合相互作用の粒子の大きさを表すパラメータ（σ＝１．０）よりも大きい。
【００４３】
このようにすれば、架橋剤粒子４と第２ポリマー粒子３０との間に斥力が発現する粒子間距離が、架橋剤粒子４と第１ポリマー粒子２０との間に斥力が発現する粒子間距離よりも長くなるので、架橋剤粒子４が第２ポリマー粒子３０から離れやすくなり、架橋剤粒子４が第１ポリマー粒子２０に集まりやすくなり、架橋反応を促進させることが可能となる。
【００４４】
本実施形態では、第１ポリマー粒子２０と第２ポリマー粒子３０との間に作用する非結合相互作用は、斥力のみを有する。
【００４５】
このように、第１ポリマー粒子２０と第２ポリマー粒子３０との間に作用する非結合相互作用に、引力がなく斥力のみが有れば、第１ポリマー粒子２０と第２ポリマー粒子３０とが反発するだけなので、第１ポリマー粒子２０と第２ポリマー粒子３０の相分離を促進および調節することが可能となる。
【００４６】
本実施形態に係るプログラムは、上記方法をコンピュータに実行させるプログラムである。このプログラムを実行することによっても、上記方法の奏する作用効果を得ることが可能となる。
【００４７】
以上、本発明の実施形態について図面に基づいて説明したが、具体的な構成は、これらの実施形態に限定されるものでないと考えられるべきである。本発明の範囲は、上記した実施形態の説明だけではなく特許請求の範囲によって示され、さらに特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれる。
【００４８】
例えば、図１に示す各部１０〜１２は、所定プログラムをコンピュータのプロセッサで実行することで実現しているが、各部を専用回路で構成してもよい。また、本実施形態では１つのコンピュータにおけるプロセッサが各部１０〜１２を実装しているが、少なくとも１又は複数のプロセッサに分散して実装してもよい。
【００４９】
上記の各実施形態で採用している構造を他の任意の実施形態に採用することは可能である。各部の具体的な構成は、上述した実施形態のみに限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形が可能である。
【符号の説明】
【００５０】
１…システム
１０…設定部
１１…変形処理実行部
１２…架橋処理実行部
【図１】
2020012075.tif 000005
【図２】
2020012075.tif 000006
【図３】
2020012075.tif 000007
【図４】
2020012075.tif 000008
【図５】
2020012075.tif 000009
【図６】
2020012075.tif 000010
【図７】
2020012075.tif 000011
【図８】
2020012075.tif 000012