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

好的，谢谢！

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

迭代的思 ...

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

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

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

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

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

显然不可能用AM1-BCC，提这个没有任何意义，怎么想也不可能轮得到用AM1-BCC。

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

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