24小时热门版块排行榜    

查看: 2994  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 299求调剂 +4 △小透明* 2026-03-17 4/200 2026-03-17 20:09 by peike
[考研] 0703化学336分求调剂 +4 zbzihdhd 2026-03-15 5/250 2026-03-17 17:33 by ruiyingmiao
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 5/250 2026-03-17 17:27 by Little-xue
[考研] 344求调剂 +4 knight344 2026-03-16 4/200 2026-03-17 17:27 by ruiyingmiao
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
请填处理意见