| 查看: 3009 | 回复: 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 …… |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- TIP3P 和 SPC/E 适用的体系
已经有7人回复
chibang367
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 863.2
- 红花: 1
- 帖子: 197
- 在线: 78.3小时
- 虫号: 4191967
- 注册: 2015-11-02
- 性别: MM
- 专业: 工程热物理与能源利用
6楼2016-03-06 23:11:56
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
200