求助：forcite退火运算总是优化的第一步就已经是能量最小

goosebay · 发表于 Post on 2021-9-14 20:23:37

体系：10条聚合度为10的葡萄糖模型，无定形晶胞初始密度：0.5g/cm3

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

AaronO_o · 发表于 Post on 2021-9-23 12:53:10

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

goosebay · 发表于 Post on 2021-9-26 22:01:22

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Material Studio] 求助：forcite退火运算总是优化的第一步就已经是能量最小

浏览过的版块