24小时热门版块排行榜    

查看: 3081  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ustbmars

铜虫 (小有名气)

[交流] 【求助】NVT系综在MD模拟过程中的能量变化 已有5人参与

在对NVT系综的分子动力学模拟过程中,粒子数N、体积V均保持不变,温度T围绕设定值上下波动,那么能量E应该怎么变化呢?我在模拟过程中观察到能量随时间增加而不断下降,但下降趋势小于10meV/atom for every 1000 ionic steps. 请问在NVT系综中能量的这种变化正常吗?请大拿们帮助分析一下。谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weizhiyong

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
在NVT中,时间足够长时,能量应该也会上下波动,基本是一个定值,不会出现逐渐下降的情况。
因为温度和模拟盒子的体积恒定后,原子之间的相对距离不会变化很大,势能也基本不会改变。
3楼2010-07-06 15:10:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

老虎大王

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+1):谢谢 2010-07-06 10:12:33
你是在多少步以内观察到总能量或者总势能不断下降的?
2楼2010-07-06 09:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbmars

铜虫 (小有名气)

4000多步吧,是不是步数太少了?我用的是第一性原理分子动力学,几千步的计算量已经很大了。
引用回帖:
Originally posted by 老虎大王 at 2010-07-06 09:38:10:
你是在多少步以内观察到总能量或者总势能不断下降的?

4楼2010-07-06 22:32:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+3):谢谢 2010-07-07 10:23:16
AIMD,就不好说了。因为在第一原理计算中,MD方法经常作为一种优化方法,所以其中的能量值不断下降,大约也是正常的吧。但是严格来说,为了能够安心地进行平均,我想,还是要计算到能量大致不变才好吧。
5楼2010-07-07 07:59:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见