24小时热门版块排行榜    

查看: 3082  |  回复: 6

ustbmars

铜虫 (小有名气)

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

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

老虎大王

木虫 (著名写手)

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

weizhiyong

金虫 (正式写手)


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

10133202

铁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主,您好。我做AIMD模拟32个原子,SMASS=0;POTIM=3;PREC=low;K-point:1*1*1.结果和您的一样,F,T,Ek呈现振荡收敛趋势,但是总能E一直下降。您的问题解决了吗?
6楼2018-12-08 13:11:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖仅楼主可见
7楼2021-07-23 11:42:35
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 ustbmars 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见