| 查看: 2577 | 回复: 9 | |||
[求助]
SPC水分子模型，求密度有些许问题 已有1人参与
|
|
RT，我用MS建立了一个6000分子的SPC水分子模型。但是得到的在293K的密度总是0.96这样，距离实际值或者SPC水分子的实验值都有0.2～0.3的差距（虽然看上去不算大，不过这没关系吗？） 我想知道这是什么导致的，density=mass(all)/(参数*vol)，但是mass是不变的，密度的体积就只跟体积有关了。为什么在温度到达293附近的时候波动还比较大（270～300上下波动），密度的波动（0.95～0.97上下波动）。请大牛们帮我分析一下到底哪里有问题。 in文件和data文件如下： in文件： units real dimension 3 boundary p p p atom_style full pair_style lj/cut/coul/long 9.0 9.0 bond_style harmonic angle_style harmonic kspace_style pppm 1.0e-4 read_data water.lammps05 neighbor 2.0 bin neigh_modify delay 5 every 1 velocity all create 298 4928459 dist uniform variable density equal mass(all)*10/(6.0221417930*vol) timestep 0.01 thermo 2000 thermo_style custom step temp press vol etotal v_density fix fixShake all shake 0.0001 100 0 b 1 a 1 fix 1 all npt temp 293.0 293.0 100 iso 1.0 1.0 1000 drag 0.2 run 100000 unfix 1 fix 2 all nve dump 1 all custom 1000 dumpnve.lammpstrj id type xu yu zu run 10000 data文件： LAMMPS 2005 data file for water 6000 atoms 4000 bonds 2000 angles 0 dihedrals 0 impropers 2 atom types 1 bond types 1 angle types 0.000000000 39.107300000 xlo xhi 0.000000000 39.107300000 ylo yhi 0.000000000 39.107300000 zlo zhi Masses 1 15.999400 2 1.007970 Pair Coeffs 1 0.1553000000 3.1660000000 2 0.0000000000 0.0000000000 Bond Coeffs 1 450 1 Angle Coeffs 1 47.0000 109.4700 Atoms 1 1 1 -0.82 7.597767833 9.408946512 5.391758138 2 1 2 0.41 6.641116667 9.263121029 5.557195421 3 1 2 0.41 7.756749235 8.739511125 4.691137645 …… |
回复此楼» 猜你喜欢
求调剂
已经有11人回复
-
材料工程085601，270求调剂
已经有44人回复
-
085600材料与化工329分求调剂
已经有25人回复
-
药学305求调剂
已经有3人回复
-
085600材料与化工349分求调剂
已经有17人回复
-
农学0904 312求调剂
已经有4人回复
-
08工学 309分求调剂
已经有9人回复
-
考研求调剂
已经有4人回复
-
22408 312求调剂
已经有8人回复
-
279学硕食品专业求调剂院校
已经有26人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- TIP3P 和 SPC/E 适用的体系
已经有7人回复
liuasia
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 102.6
- 散金: 200
- 帖子: 69
- 在线: 122.6小时
- 虫号: 3001880
- 注册: 2014-02-27
- 性别: GG
- 专业: 艺术美学
2楼2014-04-09 16:10:49
laughingjudy
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 185.8
- 帖子: 33
- 在线: 31.5小时
- 虫号: 2738781
- 注册: 2013-10-20
- 专业: 工程热物理相关交叉领域
3楼2014-06-13 14:59:07
4楼2014-06-18 21:18:19
5楼2014-08-29 22:03:44
chibang367
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 863.2
- 红花: 1
- 帖子: 197
- 在线: 78.3小时
- 虫号: 4191967
- 注册: 2015-11-02
- 性别: MM
- 专业: 工程热物理与能源利用
6楼2016-03-06 23:11:56
cuso4852cuoh
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 223.4
- 散金: 20
- 帖子: 170
- 在线: 49.6小时
- 虫号: 1481571
- 注册: 2011-11-08
- 专业: 工程热物理相关交叉领域
7楼2018-11-15 08:37:03
8楼2018-12-12 14:08:21
cuso4852cuoh
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 223.4
- 散金: 20
- 帖子: 170
- 在线: 49.6小时
- 虫号: 1481571
- 注册: 2011-11-08
- 专业: 工程热物理相关交叉领域
9楼2018-12-27 21:18:42
cuso4852cuoh
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 223.4
- 散金: 20
- 帖子: 170
- 在线: 49.6小时
- 虫号: 1481571
- 注册: 2011-11-08
- 专业: 工程热物理相关交叉领域
10楼2019-03-09 13:57:29
5