[解析事例]MDおよびMO/DFTを用いた比誘電率の評価
- 量子化学・DFT
- 全原子MD
- 光学・電気・磁気
- マテリアルサイエンス
User Only
様々な分子の比誘電率の評価
目的と手法
比誘電率の種類は電子分極、イオン分極、配向分極があります。実験ではこれらの和が誘電率として求められますが、シミュレーションで評価する場合、それぞれ適切な手法を選択し計算を実施する必要があります。
MDによる評価
MDでは振動と配向に起因する分極を評価することができます。下記の式を用いて各原子の電荷から求められる双極子モーメント(µ)の和の時間揺らぎから求めます。
- µi: 分子iの双極子モーメント
- qk: 原子電荷
- rk: 原子位置
- M: 分子の双極子モーメントの総和
- εr: 比誘電率
- T: 原子電荷
- V: 原子位置
QMによる評価
QMでは電子分極を評価することができます。分子分極率αからClausius-Mossottiの式を用いて求めます。数密度Nは実験値やシミュレーション、QSPRなどを用いて別途求める必要があります。
- εr: 比誘電率
- ε0: 真空の誘電率
- N: 電気双極子の数密度
- α:分極率
解析結果
ベンゼンとアセトンについてMOとMDによる比誘電率の評価を行った結果を図2および表1に示します。QMはGaussianを用いて求めた分子分極率と実験値の密度を使用しました。MDはOPLS力場を用いた300K、1barの結果から算出しました。アセトンのεMDの値は同じOPLS力場で15との値が報告されており[2]、J-OCTAによる計算結果も妥当な結果であると言えます。
ベンゼンは対称性が高く永久双極子がほとんどないため、配向に起因する分極を評価するMDでは誘電率の値が非常に小さく、実験値のほとんどは電子分極に起因していることが分かります。一方、アセトンは電子分極を評価するだけでは十分ではなく、配向に起因する誘電率を評価することが重要であることが分かりました。
図1 解析モデル(左:ベンゼン、右:アセトン)
図2 比誘電率計算結果
εMO | εMD | εsim. | εexp. | |
---|---|---|---|---|
Benzene | 2.16 | 1.03 | 2.19 | 2.28 |
Acetone | 1.76 | 13.96 | 14.72 | 21.01 |
μMO[D] | μMD[D] | |
---|---|---|
Benzene | 0.013 | 0.190 |
Acetone | 3.196 | 3.110 |
表1 比誘電率計算結果
さらに、MDおよび定量的構造物性相関QSPRを用いてPVCポリマーの比誘電率を算出しました。結果を表2、3に示します。
MDでは繰り返し数10のPVCモデルを作成し、OPLS力場を設定しました。100分子を系に投入し、300K、1barで緩和を行い、平衡状態での2nsのMDの結果から比誘電率を算出しました。平衡状態での密度は1.32g/cm3、MDにより求められた比誘電率は2.92となりました。
QSPRによる比誘電率の推算では、密度1.3835 g/cm3、比誘電率2.9289となり、どちらもMDの結果と近い値を得られました。
図3 解析モデル(PVC)
(白:水素、灰色:炭素、緑:塩素)
ρMD[g/cm3] | εMD[-] | µMD[D] |
---|---|---|
1.32 | 2.92 | 2.01 |
ρMD[g/cm3] | ε298K[-] |
---|---|
1.3835 | 2.9289 |
- *参考文献
-
- [1] J.Bicerano, Prediction of Polymer Properties, 3rd Ed. Marcel Dekker, 2002
- [2] J. Chem. Eng. Data 2018, 63, 5, 1170
事例一覧
-
- タイヤの耐摩耗性の向上
- 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による界面張力の評価
- 体積弾性率評価
- 界面活性剤の影響評価
- 燃料電池の高分子電解質膜解析
- 一軸伸張解析
- 光学特性評価
- ガラス転移温度評価
- ガス拡散解析
- ガス透過解析
- 架橋構造材料の特性評価
- コンポジット材料の特性解析
- ※記載されている製品およびサービスの名称は、それぞれの所有者の商標または登録商標です。