| 查看: 2666 | 回复: 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 |
» 猜你喜欢
无聊看看时间戳打发时间
已经有3人回复
基于自然哲学类比的风化壳型稀土矿
已经有14人回复
国自然申请五篇代表作大比拼,感觉这个是最重要的
已经有7人回复
评委有多少概率知道其他专家手中有哪些人的本子?
已经有6人回复
求推荐期刊,重谢
已经有3人回复
职称论文投稿
已经有7人回复
中!中!中!
已经有4人回复
E0414, 我的本子有没有希望?
已经有17人回复
青A35岁以下通知答辩了吗
已经有4人回复
小城的小雨
已经有3人回复
2楼2017-04-20 16:23:07
gtssongchi
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 4373.1
- 散金: 136
- 红花: 2
- 帖子: 157
- 在线: 99.5小时
- 虫号: 1137504
- 注册: 2010-11-02
- 性别: GG
- 专业: 凝聚态物性 II :电子结构
3楼2017-04-20 17:19:03
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
6楼2017-05-07 21:20:46











回复此楼