| 查看: 2997 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】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 请大家帮忙解决下，谢谢！ |
回复此楼» 猜你喜欢
267一志愿南京工业大学0817化工求调剂
已经有7人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有3人回复
-
一志愿中国海洋大学，生物学，301分，求调剂
已经有4人回复
-
一志愿武理材料305分求调剂
已经有6人回复
-
材料专硕英一数二306
已经有3人回复
-
085600材料与化工调剂 324分
已经有7人回复
-
311求调剂
已经有3人回复
-
26调剂/材料/英一数二/总分289/已过A区线
已经有7人回复
-
295求调剂
已经有5人回复
-
354求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算界面结构时原子弛豫的问题
已经有12人回复
- 两个原子靠近成键的过程中，系统总能量会不会发生变化？
已经有3人回复
- NPT系综的平衡问题
已经有10人回复
- 为什么做表面弛豫时，真空层位置不同会对结果有显著影响？
已经有13人回复
- 【讨论】高温高压的体系需要注意的环节
已经有18人回复
- 【讨论】lammps中怎样检验模拟系统在NVT系综下弛豫是否达到平衡？
已经有6人回复
- 【求助】弱问lammps弛豫问题
已经有3人回复
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+10):thank you 2010-12-03 15:27:11
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+10):thank you 2010-12-03 15:27:11
|
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稳定了, 给大家上新图. |
4楼2010-12-03 11:26:48
2楼2010-08-30 15:17:27
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
3楼2010-12-03 11:16:09
5