24小时热门版块排行榜    

查看: 3305  |  回复: 3

微风breeze

铁虫 (初入文坛)

[交流] 【求助】NVT弛豫问题 已有2人参与

现在在NVT弛豫这块遇到了问题:
  1、NVT弛豫之前需不需要minimize
  2、用NVT弛豫是不是就是fix NVT然后run一段时间,直到势能和总能在某一个值上下波动,即平衡?
  3、给定初始速度应该在哪里给定,minimize之前,还是minimize与fix NVT之间,还是别的?还有关于velocity命令的迷惑,velocity create XX seed(自己随便写还是用别的方法产生)temp temp_compute(这个计算的温度和前面的XX问题有什么区别,它对于velocity这个命令到底有什么作用)
   4、为什么我fix NVT之后,跑了很长时间还是无法达到我指定的温度,比如
    fix NVT 293.0 293.0 5 drag 0.2
      run 120000之后,得到的图像是这样的,
   
      好像没达到平衡,不知道是哪出错了,温度也不能达到293K
    请大家帮忙解决下,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微风breeze

铁虫 (初入文坛)

小木虫的人呢?。。。。没人回复自己先顶下
2楼2010-08-30 15:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

引用回帖:
Originally posted by 微风breeze at 2010-08-30 15:17:27:
小木虫的人呢?。。。。没人回复自己先顶下

小虫在吃菜.
Life, Love, Laugh.
3楼2010-12-03 11:16:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+10):thank you 2010-12-03 15:27:11
引用回帖:
Originally posted by 微风breeze at 2010-08-25 18:34:52:
现在在NVT弛豫这块遇到了问题:
  1、NVT弛豫之前需不需要minimize
  2、用NVT弛豫是不是就是fix NVT然后run一段时间,直到势能和总能在某一个值上下波动,即平衡?
  3、给定初始速度应该在哪里给定,minimiz ...

1、NVT弛豫之前需不需要minimize

当然需要!! minimization的主要作用是消除原子间的 overlap, 不做其他的事情.


  2、用NVT弛豫是不是就是fix NVT然后run一段时间,直到势能和总能在某一个值上下波动,即平衡?

对的.

只需要看总能量就可以了, 因为N和T是定值, 那么kinetic energy是定值, 所以势能和总能是一个概念了.




  3、给定初始速度应该在哪里给定,minimize之前,还是minimize与fix NVT之间,还是别的?还有关于velocity命令的迷惑,velocity create XX seed(自己随便写还是用别的方法产生)temp temp_compute(这个计算的温度和前面的XX问题有什么区别,它对于velocity这个命令到底有什么作用)

参见 1. 回答.  minimization 不需要你初速度的信息, 所以放哪里都无所谓, 你可以尝试的. 便于理解的话, 你就这样放在minize后面好了.

velocity的问题, 手册讲得很清楚了, 再去读一下吧.



   4、为什么我fix NVT之后,跑了很长时间还是无法达到我指定的温度,比如
    fix NVT 293.0 293.0 5 drag 0.2

显然你需要找到一个优化值, 而不是就用0.2, 参见这个帖子:

http://muchong.com/bbs/viewthread.php?tid=2435901&fpage=1

虽然fix NVT 用的是Nose-Hoover 和 berendsen thermostat 实现方法上有些不同, 但是基本的道理还是一样的.


好了, 注意了这些以后, 等你的T稳定了, 给大家上新图.
Life, Love, Laugh.
4楼2010-12-03 11:26:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 微风breeze 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +3 process2012 2026-08-23 4/200 2026-08-23 12:54 by huixian257
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见