| 查看: 2739 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
gtssongchi木虫 (小有名气)
|
[求助]
fix nvt 温度设置为300K，可是结果却显示只有150K，是什么原因呢？ 已有2人参与
|
||
|
大家好，lammps初学，请问我在设置中出了什么问题，导致最后热力学输出的温度只有150K #lammps in file for melt of Cu echo screen units metal boundary p p p atom_style atomic lattice fcc 3.61 #建立一个晶格参数为3.61埃的fcc晶体 region box block 0 10 0 10 0 10 create_box 1 box create_atoms 1 box pair_style eam pair_coeff * * ./potentials/Cu_u3.eam #格式确保正确 timestep 0.003 #单位是ps neighbor 2.0 bin neigh_modify delay 0 every 5 check yes velocity all create 300 474747 dist gaussian #variable N equal abs{floor(x[100]+10000)} #velocity all create N dist gaussian #fix 1 all nvt temp 300 300 100.0 # 100是耦合系数，给出驰豫时间100×时间单位 fix 1 all nve fix 2 all temp/berendsen 300.0 300.0 100.0 thermo 100 #thermo_style 默认 #dump 1 all custom 100 dump.lammpstrj id type x y z run 10000 |
回复此楼» 猜你喜欢
fileCode有新解读？
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
gtssongchi
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 4373.1
- 散金: 136
- 红花: 2
- 帖子: 157
- 在线: 99.5小时
- 虫号: 1137504
- 注册: 2010-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2017-04-20 17:19:03
2楼2017-04-20 16:23:07
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
4楼2017-04-20 17:22:06
gtssongchi
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 4373.1
- 散金: 136
- 红花: 2
- 帖子: 157
- 在线: 99.5小时
- 虫号: 1137504
- 注册: 2010-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2017-04-20 21:12:20
100