| 查看: 2505 | 回复: 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 …… |
回复此楼» 猜你喜欢
298求调剂
已经有3人回复
-
336求调剂
已经有4人回复
-
296求调剂
已经有3人回复
-
266求调剂
已经有3人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有10人回复
-
初试 317
已经有4人回复
-
化学调剂
已经有3人回复
-
286求调剂
已经有7人回复
-
0703化学调剂
已经有3人回复
-
一志愿重庆大学085700资源与环境总分308求调剂
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- TIP3P 和 SPC/E 适用的体系
已经有7人回复
4楼2014-06-18 21:18:19
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
5楼2014-08-29 22:03:44
