24小时热门版块排行榜    

查看: 2638  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 大佬们,打的核磁氢谱与chemdraw预测的有偏差 +5 来了个 2024-06-26 5/250 2024-06-26 17:28 by Jasonlin7758
[基金申请] 基你太美 +13 lyfbangong 2024-06-24 15/750 2024-06-26 17:12 by lyfbangong
[考研] 刚当完兵回来想考研 +7 五官在线 2024-06-23 18/900 2024-06-26 17:05 by 1158057902
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[论文投稿] 我在写关于多维空间的论文,希望能与大家交流,准备在science上发表 +8 lgf519 2024-06-24 10/500 2024-06-26 16:25 by 梦渺岚烟
[基金申请] 博后面上今天下午会公布吗?有无消息? +26 hajkdfdf 2024-06-24 43/2150 2024-06-26 16:04 by kyukitu
[基金申请] 博后面上和特助今天出吗? +41 逗您玩 2024-06-21 78/3900 2024-06-26 16:00 by HAPPY_0225
[硕博家园] 回家两天,不想呆了 +6 368ghnf 2024-06-22 12/600 2024-06-26 13:48 by lyfbangong
[公派出国] CSC博士联培对将来就业有用吗 +4 也就这样 2024-06-22 4/200 2024-06-26 10:46 by lxs7405
[第一性原理] Vasp 版权问题 10+4 竹叶青9 2024-06-22 5/250 2024-06-25 14:58 by 无所谓109
[硕博家园] 数据不好 +5 Hetai 2024-06-23 7/350 2024-06-25 12:37 by 1591099
[考博] 申请25博士,可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[有机交流] 求助析晶问题 20+4 dengdawang 2024-06-24 5/250 2024-06-24 21:22 by cc116
[金属] EBSD的解析率只有10% +3 wallace6666 2024-06-20 7/350 2024-06-24 16:52 by wallace6666
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 说博后基金7月出的真打电话了吗? +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[公派出国] 博士csc联培会看重第一学历学校层次吗 +4 也就这样 2024-06-23 4/200 2024-06-24 08:18 by 晓目崇
[基金申请] 国自然青年基金,1A4B能上会吗?青年和面上的上会标准是一样的吗? +19 今晚推荐22 2024-06-20 32/1600 2024-06-23 23:17 by andywei1028
[基金申请] 工材口青年基金大概什么样能上会? +15 今晚推荐22 2024-06-20 21/1050 2024-06-22 23:04 by qbn0326
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
信息提示
请填处理意见