| 查看: 2384 | 回复: 3 | |||
木棉和我木虫 (职业作家)
|
[求助]
lammps运行时nvt下对应压力10391145bar，npt下对应体积扩大12倍，哪里出了问题？ 已有1人参与
|
| 诸位。采用lammps计算MgO-TiO2-B2O3体系的玻璃化转变或结晶转变，关注体积和能量变化。势函数采用Materials Explorer5.0软件势函数库中的Born-Mayer-Huggins potential parameters，文献来源为Hirao K, Kawamura K. Material design using personal computer[J]. Shokabo, Tokyo, 1994, 52。in文件的主要语句为：unit metal；timestep 0.001；velocity all create 3000 4928459 dist gaussian；fix 0 all nve；fix 1 all temp/rescale 1 3000 3000 5.0 1； run 8000；unfix 0；unfix 1；fix 2 all npt temp 3000 300 0.1 iso 10378167 10378167 1；run 27000；unfix 2。密度按照实验测定设置为4.00g/cm3。nvt下运行，显示压力为10391145bar，体积为39136.177埃3(对应于之前设置的密度)。npt下运行，压力设置为1bar，体积则为482515.33 埃3（与之前比扩大了12倍）。请问，哪里出了问题？ |
回复此楼» 猜你喜欢
越来越担心硕士毕业后没有竞争力了
已经有15人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
2楼2018-06-13 13:51:23
木棉和我
木虫 (职业作家)
- 应助: 5 (幼儿园)
- 金币: 9573.7
- 散金: 2362
- 红花: 3
- 帖子: 3424
- 在线: 686.1小时
- 虫号: 1852823
- 注册: 2012-06-08
- 性别: GG
- 专业: 冶金物理化学与冶金原理
3楼2018-06-13 14:37:23
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
4楼2018-06-16 19:37:12
50