| 查看: 5071 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
对天空说lrj主管区长
|
[求助]
Nose-Hoover 热浴控温，lammps来讲就是指的fix nvt 吧？为什么温度会大增导致原子丢失 已有6人参与
|
||
|
我做热导率想先按照文献的方法模拟一遍验证自己的程序对不对，据文献里的表述“前50 万步，将整个系统置于NVE 系综中，对上下固壁和液体区域分别使用温度/速度调节法使系统在100 K 下达到平衡状态。后500 万步，利用Nose-Hoover 热浴来对上下固壁分别进行控温，给液体区域沿y 向施加一个温度梯度.” 对于上面描述的控温我使用如下的命令 fix NVE all nve fix 2 all temp/rescale 100 100 100 0.1 1.0 run 500000 unfix 2 fix NHlo Plow nvt temp 110 110 100 fix NHhi Pup nvt temp 90 90 100 run 5000000 第二个run的时候我没有对液相进行任何控温，Nose-Hoover 热浴控温，对于lammps来讲就是指的fix nvt 吧？但是这样运行之后前面还有nve,体系的温度就会不断增加，最后丢失原子，麻烦前辈看看我的问题出在哪里？ |
回复此楼» 猜你喜欢
271求调剂
已经有12人回复
-
322求调剂
已经有3人回复
-
环境工程 085701，267求调剂
已经有15人回复
-
材料0856 英一数二 323 求调剂
已经有3人回复
-
301求调剂
已经有3人回复
-
348求调剂
已经有3人回复
-
350求调剂
已经有3人回复
-
309分085801求调剂
已经有7人回复
-
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
dgwfly
实习版主
- 应助: 0 (幼儿园)
- 金币: 683.9
- 散金: 258
- 红花: 9
- 帖子: 543
- 在线: 202小时
- 虫号: 3454332
- 注册: 2014-10-04
- 性别: MM
- 专业: 数论
|
壁面加中间流体，并且最外层壁面原子固定不动，我也遇到类似的问题，跑NVE的时候，局部的壁面原子会跑出去，温度能量突然急剧上升。你把时间步长改为多少呐？ 发自小木虫Android客户端 |
12楼2016-12-16 23:06:40
3楼2015-07-29 02:37:49
对天空说lrj
管理员
- 应助: 1 (幼儿园)
- 金币: 1051.3
- 散金: 97
- 红花: 1
- 帖子: 383
- 在线: 182.3小时
- 虫号: 1466527
- 注册: 2011-10-29
- 性别: MM
- 专业: 化学反应工程
4楼2015-07-29 08:48:52
luzihen
兑换贵宾
- 应助: 47 (小学生)
- 金币: 833
- 散金: 265
- 红花: 29
- 帖子: 605
- 在线: 329.9小时
- 虫号: 2243418
- 注册: 2013-01-15
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
|
我做热导率想先按照文献的方法模拟一遍验证自己的程序对不对，据文献里的表述“前50 万步，将整个系统置于NVE 系综中，对上下固壁和液体区域分别使用温度/速度调节法使系统在100 K 下达到平衡状态。后500 万步，利用Nose-Hoover 热浴来对上下固壁分别进行控温，给液体区域沿y 向施加一个温度梯度.” ############################################################### 这里文献的描述不准确，NVE系宗能量是invariant。 ############################################################### 对于上面描述的控温我使用如下的命令 fix NVE all nve fix 2 all temp/rescale 100 100 100 0.1 1.0 run 500000 unfix 2 fix NHlo Plow nvt temp 110 110 100 fix NHhi Pup nvt temp 90 90 100 run 5000000 第二个run的时候我没有对液相进行任何控温，Nose-Hoover 热浴控温，对于lammps来讲就是指的fix nvt 吧？但是这样运行之后前面还有nve,体系的温度就会不断增加，最后丢失原子，麻烦前辈看看我的问题出在哪里？ ############################################################### 第二个run的时候你每步的对上下两个区域进行了两次积分，一个根据NVE，就是F=ma，另一次根据NVT，就是F=ma+一个量。我估计你体系出问题就是这个原因。 所以，在第二个run的时候，你应该对上下两部分NVT，中间NVE。不要有overlap。 还有一种方法是用郎之万thermalstat，lammps里的郎之万热浴不会直接进行积分，所以可以写成fix all nve+fix Nhi Langevin。 ############################################################### |
5楼2015-09-03 15:23:05
