关于NVT模拟向NVE模拟的转变问题(已解决）

Deathaliveone · 发表于 Post on 2022-10-24 09:53:04

本帖最后由 Deathaliveone 于 2022-10-24 14:00 编辑

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

3163582842 · 发表于 Post on 2024-3-28 18:18:42

请问楼主，那可以用nvt平衡，然后使用nve使高低温平衡吗?nvt平衡了2ns.

Deathaliveone · 发表于 Post on 2022-10-24 13:58:45

问题已解决，可能原因是初始设置NVE+控温实现了NVT系综，后取消控温为NVE系综，系统无法立刻转变为NVE系综，执行单一NVE命令后一段时间，从近平衡态到达新的平衡态，能量先升高后趋于稳定。如图所示

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

[Lammps] 关于NVT模拟向NVE模拟的转变问题(已解决）