这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:dielec2020 [2021/08/06 16:33] – liu.jun | adf:dielec2020 [2024/05/31 11:19] (当前版本) – [小分子介电常数] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
======极化率、聚合物介电常数的计算====== | ======极化率、聚合物介电常数的计算====== | ||
- | 计算思路:截取聚合物的多个重复单元(原则上单元数越多,计算结果越可靠)通过恰当的加H饱和,使用ADF模块计算该分子其极化率,然后使用Clausius Mosotti方程,结合聚合物的密度,以及改分子的分子量,得到介电常数信息: | + | 直接计算整个聚合物的介电常数,当然是可以的,但是原子数太多,计算量过于巨大,没有实操性。 |
+ | |||
+ | 一种可操作的计算思路:截取聚合物的多个重复单元(原则上单元数越多,计算结果越可靠)通过恰当的加H饱和,使用ADF模块计算该分子其极化率,然后使用Clausius Mosotti方程,结合聚合物的密度,以及该分子的分子量,得到介电常数信息: | ||
{{ : | {{ : | ||
行 10: | 行 12: | ||
* M单位为 g/mol | * M单位为 g/mol | ||
* 真空介电常数8.854187817E-12 F/m | * 真空介电常数8.854187817E-12 F/m | ||
- | * 该公式使用上述单位后,等号右侧单位正好互相抵消,变为无量纲值,为了方便我们命名为A | + | * 该公式使用上述单位后,等号右侧单位正好互相抵消,变为无量纲值,为了方便我们将该值称为A |
行 24: | 行 26: | ||
* 不同泛函对结果影响很小,甚至GGA与杂化泛函之间,影响也很小,因此这里选择了计算量较小的GGA-PBE | * 不同泛函对结果影响很小,甚至GGA与杂化泛函之间,影响也很小,因此这里选择了计算量较小的GGA-PBE | ||
* 基组影响也较小,因此选择基本精度基组DZP。较重的元素可以选择TZP | * 基组影响也较小,因此选择基本精度基组DZP。较重的元素可以选择TZP | ||
- | * 计算了0.1eV~0.5eV之间3个频率的电场。其中频率与能量的变换关系是使用光子能量计算公式E=hγ。因为理论上无法计算静电场的激化,因此使用低频极限,也就是关心接近0.1eV甚至趋近于0.0eV的电场的极化率,通过计算可以发现0.1~0.5eV之间变化很小。所以,数据外推,认为0.0eV时极化率也是该数值。 | + | * 计算了0.1eV~0.5eV之间3个频率的电场。其中频率与能量的变换关系是使用光子能量计算公式E=hγ。因为理论上无法计算静电场的极化,因此使用低频极限,也就是关心接近0.1eV甚至趋近于0.0eV的电场的极化率,通过计算可以发现0.1~0.5eV之间变化很小。所以,数据外推,认为0.0eV时极化率也是该数值。 |
=====结果处理===== | =====结果处理===== | ||
行 36: | 行 38: | ||
以上,仅供参考。 | 以上,仅供参考。 | ||
+ | =====小分子介电常数===== | ||
+ | 这种算法,本身是以小分子的介电常数去近似聚合物的节点常数的,因此用来算小分子的介电常数当然也适用。 |