sobereva 发表于 2024-3-15 20:05
认真看此文并修改帖子
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sob ...

不好意思，老师。我使用lammps做得模拟，体系有224个正十六烷分子，我想要研究的是600K，30atm下的性质，我是fix npt all npt 1200 600 100*dt iso 500 500 1000*dt(就是在500atm下降温到600)，之后再 fix npt all npt 600 600 100*dt iso 500 30 1000*dt(在600K下降压到30atm)，这样跑完第二个动力学之后(也就是画的第二个npt的体积时间图)就出现了第一张图体积先减小后增大的趋势，而对于第一个npt后体积时间图像是降低的(我个人感觉没啥问题)。至于第二张图，我只是换了个体系(也是烷烃)，但是这次第一个npt跑完的体积时间图像就是第二张图片。关于模拟过程我不确定我是不是对的，但是对于出现这种趋势，我确实搞不明白，想要学习学习。

认真看此文并修改帖子
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）