| 查看: 1588 | 回复: 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）,这里就不对，进一步控压感觉也没有意义了。而且之后的控压也是处于波动中，各位怎么看？ |
回复此楼» 猜你喜欢
最失望的一年
已经有16人回复
-
存款400万可以在学校里躺平吗
已经有33人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
求助一下有机合成大神
已经有3人回复
-
求推荐英文EI期刊
已经有5人回复
-
请教限项目规定
已经有4人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
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