| 查看: 1699 | 回复: 3 | |||
[求助]
lammps模拟液态水，压力不合理，并且控压控不住，何解？ 已有2人参与
|
|
input script: units real dimension 3 boundary p p p atom_style full pair_style lj/charmm/coul/charmm 8.0 10.0 7.0 9.0 bond_style harmonic angle_style harmonic special_bonds charmm read_data data.droplet minimize 1.0e-4 1.0e-6 100 1000 reset_timestep 0 fix 1 all shake 0.0001 20 0 b 1 a 1 fix 2 all npt temp 298.0 298.0 100 iso 1.0 1.0 1000 velocity all create 298.0 123456 dist gaussian mom yes run 0 velocity all scale 298.0 timestep 1 thermo_style custom step temp press density vol thermo 1 run 20 log file: Step Temp Press Density Volume 0 298 2088.1099 1.0008098 343000 1 198.57626 901.42394 1.0007329 343026.35 2 195.61277 941.93247 1.0005897 343075.45 3 190.87646 989.45826 1.0003772 343148.32 4 184.67353 1052.2417 1.0000921 343246.15 5 177.3926 1122.761 0.99973009 343370.44 6 169.47103 1193.7707 0.99928675 343522.78 7 161.35931 1257.8854 0.99875821 343704.56 8 153.48719 1308.4637 0.99814209 343916.72 9 146.23467 1340.0094 0.99743833 344159.38 10 139.90957 1348.6115 0.99665001 344431.6 11 134.73268 1332.1385 0.99578397 344731.15 12 130.83089 1290.3488 0.99485107 345054.42 13 128.23776 1224.6234 0.99386601 345396.41 14 126.90087 1137.6957 0.99284661 345751.05 15 126.69476 1033.1263 0.99181251 346111.54 16 127.43816 914.95392 0.99078362 346470.96 17 128.91405 787.43923 0.98977849 346822.81 18 130.89062 654.85029 0.98881288 347161.49 19 133.14164 521.246 0.98789889 347482.68 20 135.46434 390.28164 0.98704467 347783.4 在step 0处，温度，密度都很合理，但是压力却与实际严重不符（实际应为1atm,这里却为2000多atm）,这里就不对，进一步控压感觉也没有意义了。而且之后的控压也是处于波动中，各位怎么看？ |
回复此楼» 猜你喜欢
J-J-W不为人知的一面
已经有19人回复
-
没消息就是被刷了呗
已经有14人回复
-
微信指数没变化，科研之友没阅读
已经有8人回复
-
时间戳他又来了
已经有12人回复
-
这种情况还有戏吗
已经有11人回复
-
同事接到电话了，我却没有
已经有4人回复
-
准备明年的基金了
已经有5人回复
-
今年的WR进展到哪一步了？
已经有6人回复
-
小木虫看见有人已经知道结果了
已经有15人回复
-
这个自发加氧反应的机理是什么？
已经有6人回复
jluchemlib
铁杆木虫 (正式写手)
- 应助: 119 (高中生)
- 金币: 11607.4
- 散金: 200
- 红花: 19
- 帖子: 521
- 在线: 356.1小时
- 虫号: 589656
- 注册: 2008-08-30
- 性别: GG
- 专业: 高分子物理与高分子物理化
2楼2017-03-07 08:16:07
3楼2017-04-03 21:58:05
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
4楼2017-04-06 01:26:44