[解析事例] リバースマッピングによるアモルファス構造の作成
- 全原子MD / 粗視化MD
- 分子構造・親和性・溶解性
- マテリアルサイエンス
長鎖ポリマーモデルの緩和構造を効率的に作成
ポリマー物性をシミュレーションで評価する際、バルクモデルの作成はモデリングの第一ステップです。
J-OCTAでポリマーバルクモデルを作成するアモルファス作成機能には、通常版とリバースマップ版の2種類のオプションがあることをご存じでしょうか。通常版の機能では希薄状態から開始して、圧縮や複数の平衡化のためのMD計算によって平衡密度のバルク状態を作成します。リバースマップ版では密度と特性比を指定することで、全原子のモデルに対応するKremer-Grest (KG) モデルを作成し、平衡化MDをKGモデルで実施後にリバースマップを行います。(図1)
通常版は汎用性がありどんなモデルでも使える機能である一方で、ポリマー鎖を十分に平衡化するためには計算時間が長くなってしまう課題があります。リバースマップ版は粗視化モデルによって平衡化を行うために、この課題をクリアしてポリマーバルクが持つ構造特徴を再現することができます。
ここではポリエチレン(PE)、cisポリブタジエン(PB)、ポリスチレン(PS)の20、40、60量体をそれぞれの機能で各3水準ずつ作成し高分子鎖の構造特徴を比較しました。MD計算のための力場はGAFFを用いました。
図2には作成されたバルクモデルの密度、回転半径<Rg>を示しました。密度についてはどちらの方法でも同程度の値が得られていますが、 <Rg>については重合度の大きいモデルにおいて2種類の作成方法に差が確認されました。
次に3水準のモデルの<R2>の平均値の分子量依存性を図3に示しました。ポリマーの<R2>や <Rg2>は重合度に比例することが知られていますが、リバースマップを用いて作成したバルクモデルではこの傾向が明確に確認されました。また、参照値として<R2>/M [1]から計算される値とも比較し、リバースマップ版で作成される構造の妥当性を確認しました。
図1. アモルファス作成機能の処理フロー
(左)通常版(右)リバースマップ版
図2. アモルファス作成機能で作成されたポリマーモデル3水準の密度と<Rg>
(左)通常版(右)リバースマップ版
図3. 作成されたPB、PE、PSモデルの<R2>の重合度依存性の比較
○:通常版、☐:リバースマップ版、破線:参照値
またリバースマップ版においては、KGモデルを用いて緩和計算を行うことで絡み合った高分子鎖の構造も作成することができます。[2]
図4は本機能で作成されたポリカーボネート25量体の絡み合ったバルクモデルです。
図4: リバースマップ機能により作成されたポリカーボネート25量体の絡み合いモデル
(左)KGモデル、(右)全原子モデル
- *参考文献
-
- [1] L. J. Fetters, D. J. Lohse, and R. H. Colby, Physical properties of polymers handbook , 445?452 (2006).
- [2] H. Nitta, T. Ozawa and K. Yasuoka, J. Chem. Phys., 159, 194903 (2023).
事例一覧
-
- タイヤの耐摩耗性の向上
- COSMO法による物性推算
- VSOP-PSによる繊維配向材のシミュレーション
- 粗視化モデルによる脂質膜の解析
- 機械学習による沸点、屈折率、比誘電率の推算
- ナノトライボロジー(アブレシブ摩耗、ナノ加工)
- コバルト酸リチウムの基底状態と弾性率の解析(SIESTA事例ページへ)
- MD-GANによる固体電池内のLiイオン拡散解析
- リバースマッピングによるアモルファス構造の作成
- 電池電極の成形プロセス(カレンダリング)における圧力と空隙率の計算
- 粗視化分子動力学を用いた複屈折の解析
- 高分子膜の相分離プロセスシミュレーション
- 機械学習によるχパラメータの推定
- 水への溶解性評価
- フィラー樹脂複合材料の熱伝導率計算
- VSOP-PSによる繊維構造への樹脂含浸プロセス計算
- mol-inferを用いたQSPRの逆解析
- FMO-DPDを用いた高分子電解質の相分離構造
- FMO-DPDを用いた脂質膜とベシクルの形成
- カルサイト(方解石)の複屈折と光吸収(SIESTA事例ページへ)
- リチウムイオン電池のシミュレーション(SIESTA事例ページへ)
- 懸濁液の粘度の評価
- 量子補正を適用した固体の定積比熱の評価
- MD-GANを用いた長時間の分子運動の予測
- ポリマーの誘電緩和
- 金属錯体の吸収スペクトル(SIESTA事例ページへ)
- ゼオライトとガス分子の相互作用(SIESTA事例ページへ)
- 蒸着膜のシミュレーション
- 機械学習QSPRとマテリアルズ・インフォマティクス
- シミュレーション結果を用いた粘弾性マスターカーブの作成
- DPDを用いた粘弾性のシミュレーション
- 表面の再構成(SIESTA事例ページへ)
- フォノン分散を利用した剛性マトリクスの算出(SIESTA事例ページへ)
- Steered MDによるポリペプチドの自由エネルギー変化の解析
- 金属の電子比熱解析(SIESTA事例ページへ)
- 格子比熱の解析(SIESTA事例ページへ)
- J-OCTAによる蓄熱材の評価
- MDおよびMO/DFTを用いた比誘電率の評価
- 溶解度係数の算出
- 機械特性評価(SIESTA事例ページへ)
- 結晶の熱膨張(SIESTA事例ページへ)
- ゼオライトへのガス吸着
- GHz周波数領域における水の誘電分散(2)
- 反応のエネルギー変化(SIESTA事例ページへ)
- GHz周波数領域における水の誘電分散(1)
- 固体表面への分子の吸着エネルギー(SIESTA事例ページへ)
- 活性化エネルギーを用いたモンテカルロ判定によるエポキシ樹脂の架橋反応
- 樹脂 複屈折のためのマルチスケールシミュレーション
- 顔料の表面エネルギー(SIESTA事例ページへ)
- スラリー塗工プロセス
- 金属表面と分子の相互作用
- フィラー充填ゴムの繰り返し伸長
- 界面特性を考慮したCFRTPの破壊挙動の解析(Digimat事例ページへ)
- グラフェンシート添加によるCFRTPのマトリクス熱伝導特性改質(Digimat事例ページへ)
- 非平衡MDによる熱伝導率の計算
- ゴム材料のレオロジーシミュレーション
- 架橋フェノール樹脂の物性評価
- エポキシ樹脂架橋反応のシミュレーションからのガラス転移温度評価
- 電池用電解液の計算
- カーボンナノチューブの分散構造と非線形構造解析
- 粗視化ポテンシャルの評価
- DPDによる溶媒蒸発シミュレーション
- ガラス状態ポリマーの一軸伸長とクレーズ形成
- 高分子-固体界面の摩擦
- MDによる粘度の評価
- 相図からのχパラメータ推算
- 高分子-固体界面の剥離
- 気体の溶解係数と自由体積
- QSPRによる物性値の推算
- レオロジー特性
- 熱硬化性樹脂の架橋構造
- MDによる比誘電率の評価
- 複合材料の非線形力学特性(LS-DYNAとの連携)
- DPDによる界面張力の評価
- MDによる溶解度パラメータの評価
- DPDによる液滴のずり変形
- MDによる界面張力の評価
- 体積弾性率評価
- 界面活性剤の影響評価
- 燃料電池の高分子電解質膜解析
- 一軸伸張解析
- 光学特性評価
- ガラス転移温度評価
- ガス拡散解析
- ガス透過解析
- 架橋構造材料の特性評価
- コンポジット材料の特性解析
- ※記載されている製品およびサービスの名称は、それぞれの所有者の商標または登録商標です。