计算化学公社

标题: 不常见有机液相分子的静态介电常数怎么计算？ [打印本页]

作者
Author: yaoyuan0711 时间: 2021-11-11 15:38
标题: 不常见有机液相分子的静态介电常数怎么计算？
各位老师好！我在用GAFF力场计算一个不太常见的有机分子的静态介电常数（查不到，相似的结构也没有，如图所示），但是计算RESP2电荷（http://sobereva.com/531）的时候需要加该分子的溶剂模型，也需要静态介电常数。这对矛盾怎么解决？
以GAFF默认的AM1-BCC电荷计算出静态介电常数？或者计算RESP2电荷的时候预估一个静态介电常数？

(, 下载次数 Times of downloads: 9)

作者
Author: sobereva 时间: 2021-11-11 16:35
显然不可能用AM1-BCC，提这个没有任何意义，怎么想也不可能轮得到用AM1-BCC。

原理上，可以在真空下用B3LYP/def2-TZVP这种像样级别结合弱极性溶剂（比如eps=10）用Multiwfn算RESP电荷，然后跑动力学得到介电常数。或者真空下计算极化率，通过Clausius–Mossotti关系预测介电常数，论坛首页google框搜，以前回复过多次。

如果要求高，就需要迭代：在前述RESP电荷下跑MD得到介电常数，然后再用Multiwfn算RESP2电荷，再用RESP2电荷跑MD得到新的介电常数...反复交替更新RESP2电荷和介电常数直到介电常数收敛。

作者
Author: yaoyuan0711 时间: 2021-11-11 16:55

sobereva 发表于 2021-11-11 16:35
显然不可能用AM1-BCC，提这个没有任何意义，怎么想也不可能轮得到用AM1-BCC。

原理上，可以在真空下用B3 ...

通过Clausius–Mossotti关系预测介电常数也需要GMX模拟分子数密度，同样也得用RESP电荷吧？

迭代的思想很厉害，谢谢sob老师。

作者
Author: sobereva 时间: 2021-11-11 17:31

yaoyuan0711 发表于 2021-11-11 16:55
通过Clausius–Mossotti关系预测介电常数也需要GMX模拟分子数密度，同样也得用RESP电荷吧？

迭代的思 ...

是
密度相对来说是属于比较好算的属性，电荷不特别理想也不至于差太多。

作者
Author: yaoyuan0711 时间: 2021-11-11 17:36

sobereva 发表于 2021-11-11 17:31
是
密度相对来说是属于比较好算的属性，电荷不特别理想也不至于差太多。

好的，谢谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)