计算化学公社

标题: 求助：forcite退火运算总是优化的第一步就已经是能量最小 [打印本页]

作者
Author: goosebay 时间: 2021-9-14 20:23
标题: 求助：forcite退火运算总是优化的第一步就已经是能量最小
体系：10条聚合度为10的葡萄糖模型，无定形晶胞初始密度：0.5g/cm3

按照论文里的步骤建立无定形晶包，然后NVT系综下从298到700K进行退火运算，然后NPT压缩和解压缩。在进行NVT退火运算时产生了20个过程结构文件，但是第一个就已经是能量最低的结构了，也就是说在进行退火之前的结构已经能量最小了，那后边的压缩和解压过程也是多余的，请问这究竟是怎么回事啊？

作者
Author: AaronO_o 时间: 2021-9-23 12:53
可以考虑检查一下退火前后的力场和电荷分配是否均合理。
库伦作用比vdW大了太多，且总势能为正，个人会怀疑力场、尤其是电荷是否合理。

作者
Author: goosebay 时间: 2021-9-26 22:01

AaronO_o 发表于 2021-9-23 12:53
可以考虑检查一下退火前后的力场和电荷分配是否均合理。
库伦作用比vdW大了太多，且总势能为正，个人会怀 ...

尝试了几个立场，compass、universial、dreiding等都试过，依然没有解决。请问电荷分配是在energy里的charges吗？我试过更改charges中的use current 为 QEq、Gasteiger等也都无济于事。是不是在建立聚合物模型的时候就要考虑电荷的问题呢？谢谢解答

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