计算化学公社

标题: 关于NVT模拟向NVE模拟的转变问题(已解决） [打印本页]

作者
Author: Deathaliveone 时间: 2022-10-24 09:53
标题: 关于NVT模拟向NVE模拟的转变问题(已解决）
本帖最后由 Deathaliveone 于 2022-10-24 14:00 编辑

在lammps程序跑反应，目前确认力场没有问题，已经进行了多次恒温模拟。最近想进行绝热模拟，我的想法是取nve系综，先升温，再关闭控温命令。即关闭temp/beremdsen，继续nve系综。但是第一次这样做的时候发现体系的总能量Etotal持续增大，似乎违背了微正则系综中的能量不变的δ分布的基本假设。第二次尝试升温之后把berendsen控温和第一个nve系综全部关闭，重新设置nve系综，但是模拟输出的热力学量中，Eotal依然在变化，且表现为逐渐增大，温度仍然稳定控温命令附近。下面给出第二次系综设定和能量曲线。请问我的问题在哪？是否直接用NVT命令控温模拟再换成NVE命令可以i避免这种情况发生呢？
反应体系是做含能材料的热解。

(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 3)

作者
Author: Deathaliveone 时间: 2022-10-24 13:58
问题已解决，可能原因是初始设置NVE+控温实现了NVT系综，后取消控温为NVE系综，系统无法立刻转变为NVE系综，执行单一NVE命令后一段时间，从近平衡态到达新的平衡态，能量先升高后趋于稳定。如图所示 (, 下载次数 Times of downloads: 5)

作者
Author: 3163582842 时间: 2024-3-28 18:18
请问楼主，那可以用nvt平衡，然后使用nve使高低温平衡吗?nvt平衡了2ns.

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