不同分子的相對(duì)介電常數(shù)計(jì)算

目的和方法

介電常數(shù)有三個(gè)分量：電子極化、離子極化和定向極化。在實(shí)驗(yàn)中，它們的總和被認(rèn)為是介電常數(shù)，但在模擬中進(jìn)行計(jì)算時(shí)，應(yīng)選擇合適的方法并對(duì)每種方法分別進(jìn)行計(jì)算。

分子動(dòng)力學(xué)計(jì)算(MD)

MD讓我們可以估測(cè)分子因振動(dòng)和取向產(chǎn)生的極化。相對(duì)介電常數(shù)可以由各個(gè)原子電荷偶極矩之和的時(shí)間波動(dòng)得到，公式如下：
μi : 分子i的偶極矩

qk : 原子電荷

rk : 原子空間

M : 分子偶極矩之和

ε r : 相對(duì)介電常數(shù)

T: 系統(tǒng)溫度

V: 系統(tǒng)體積
分子軌道法計(jì)算（MO）/密度泛函理論（DFT）計(jì)算

MO/DFT讓我們可以估測(cè)電子極化，由分子極化率計(jì)算相對(duì)介電常數(shù)。

εr : 相對(duì)介電常數(shù)

ε0 : 真空介電常數(shù)

N : 電偶極子的數(shù)密度

α : 極化率

模擬成果

圖2和表1給出了用MO和MD計(jì)算苯和丙酮的相對(duì)介電常數(shù)結(jié)果。其中MO估測(cè)值由高斯09測(cè)定的分子極化率和實(shí)驗(yàn)密度得出，MD估測(cè)值由300k和1百帕OPLS力場(chǎng)條件下計(jì)算液態(tài)體模型得出。在相同的OPLS力場(chǎng)中，丙酮的εMD也顯示為15 [2]。因此J-OCTA的計(jì)算是有效的。

由于苯具有極高的對(duì)稱性，幾乎沒(méi)有永久偶極子，使用MD來(lái)估測(cè)相對(duì)介電常數(shù)時(shí)因分子振動(dòng)和取向引起的極化是非常小的。這表明大部分實(shí)驗(yàn)結(jié)果是由電子極化得出。而丙酮?jiǎng)t相反，僅僅估測(cè)電子極化遠(yuǎn)遠(yuǎn)不夠，同時(shí)估測(cè)取向極化也非常重要。

圖1 仿真模型（左：苯環(huán) 右：丙酮）

圖2 相對(duì)介電常數(shù)估測(cè)值
表1 相對(duì)介電常數(shù)估測(cè)值和偶極矩

相對(duì)介電常數(shù) 偶極矩
同時(shí)MD和QSPR（定量構(gòu)效關(guān)系）也用來(lái)計(jì)算PVC聚合物的相對(duì)介電常數(shù)，其結(jié)果如表2和3所示。使用MD計(jì)算時(shí)，我們重復(fù)建10次建模過(guò)程并設(shè)置一個(gè)OPLS力場(chǎng)。我們把100個(gè)分子放入體系中，在300k和1百帕條件下釋放，并在2毫秒體系達(dá)到平衡時(shí)根據(jù)MD計(jì)算出相對(duì)介電常數(shù)。MD計(jì)算出平衡時(shí)密度為1.32 g/cm3，相對(duì)介電常數(shù)為2.92。

用QSPR計(jì)算相對(duì)介電常數(shù)時(shí)顯示密度為1.38g/cm3，相對(duì)介電常數(shù)為2.93，與MD結(jié)果相近。

圖3 PVC仿真模型
（白色：氫 灰色：碳 綠色：氯）
表2 MD計(jì)算的相對(duì)介電常數(shù)與QSPR計(jì)算的相對(duì)介電常數(shù)

· 參考文獻(xiàn)

· [1] J.Bicerano，《聚合物性能預(yù)測(cè)》，第三版，Marcel Dekker, 2002年

· [2] J. Chem.Eng. Data 2018, 63, 5, 1170

· [3]轉(zhuǎn)載自上海庭田信息科技有限公司

審核編輯 黃宇