| 查看: 6486 | 回复: 7 | ||
[求助]
lammps新手用tip4p模型模拟水
|
|
现在想用tip4p模型模拟液态水，计算时总是出现错误ERROR on proc:TIP4P hydrogen has incorrect atom type 一开始用pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5就会出错，试着换了个势函数 pair_style lj/cut 8.5 就没出错了，所以我觉得可能是势函数的问题。不过我比较想用tip4p。有没有什么办法呢？本人是非常新的新手，希望能讲的详细一点点。 下面是我的in文件： units metal atom_style full dimension 3 boundary p p p newton on read_data ye.data pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5 pair_coeff 1 1 0.1521 3.154000 8.5 pair_coeff * 2 0.000000 0.000000 8.5 bond_style harmonic bond_coeff 1 450 0.9572 angle_style harmonic angle_coeff 1 55 104.52 kspace_style pppm/tip4p 0.0001 neigh_modify delay 0 every 1 neighbor 2.0 bin timestep 0.0002 velocity all create 400 72489 mom yes rot yes fix 5 all box/relax iso 10 vmax 0.001 min_style cg minimize 1.0e-12 1.0e-12 10000 10000 unfix 5 reset_timestep 0 fix 2 all npt temp 300 300 0.02 iso 10 10 0.2 pchain 8 drag 0.5 thermo 500 #(每200步输出热力学结果) thermo_style custom step temp pe etotal press vol density enthalpy xlo xhi ylo yhi zlo zhi dump 1 all cfg 1000 dump2.temp_*.cfg mass type xs ys zs x y z run 10000 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有19人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有15人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有19人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有19人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有24人回复
-
面上提前没消息，有中的吗
已经有11人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MS粗粒化模型如何导入lammps
已经有10人回复
- lammps的in文件错误求助
已经有6人回复
- MS模型怎么导出lammps能读取的文件？win7系统
已经有8人回复
- 修改lammps程序的neighbor命令
已经有4人回复
- 求助LAMMPS两次运算结果不一致问题！
已经有11人回复
- 急求！！！MS中怎么画水分子模型啊？
已经有6人回复
- MS生成Lammps的data文件，文件中原子的类型和顺序并不是自己想要的，请问如何修改？
已经有8人回复
- pair_style lj/cut/coul/cut command的使用
已经有5人回复
- 【求助】关于LAMMPS建模
已经有9人回复
- 经过MS模型导入到LAMMPS，在运算时出现的错误求助
已经有22人回复
- 求助LAMMPS关于MSD计算的问题
已经有17人回复
- lammps data file中水的H，O排列
已经有4人回复
- lammps能用TIP4P/ICE来算冰吗
已经有5人回复
- 在LAMMPS中怎么建立图示的模型
已经有6人回复
2楼2014-12-11 15:03:13
3楼2015-01-12 15:37:17
chinaboy798
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 299.3
- 散金: 20
- 红花: 8
- 帖子: 177
- 在线: 69.6小时
- 虫号: 2007979
- 注册: 2012-09-17
- 性别: GG
- 专业: 金属材料的力学行为
4楼2015-02-07 14:30:01
mileyfan
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 221.8
- 帖子: 53
- 在线: 42.4小时
- 虫号: 3333003
- 注册: 2014-07-21
- 专业: 气体放电与放电等离子体技
5楼2015-07-17 11:36:06
mileyfan
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 221.8
- 帖子: 53
- 在线: 42.4小时
- 虫号: 3333003
- 注册: 2014-07-21
- 专业: 气体放电与放电等离子体技
★ ★ ★ ★ ★
abinitio: 金币+5, 感谢分享经验 2015-07-23 16:57:47
abinitio: 金币+5, 感谢分享经验 2015-07-23 16:57:47
|
帖子是很久以前的了，不知道楼主后来怎么办。分享一点我的经验。 根据mannual里面的语句介绍 pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5 中第一个“1”表示O原子atom type =1, 后面的“2”表示H原子atom type=2，这个要与输入的data文件对应，可能楼主的data文件里是反的，H的atom type=1了。 在 pair_style lj/cut 8.5语句中，没有关于这类参数的定义，只设定截断距离，因此不存在与data文件对应问题，所以不报错吧。 希望能给以后遇到同样问题的木虫一点小小的帮助~ |
6楼2015-07-19 09:35:50
chibang367
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 863.2
- 红花: 1
- 帖子: 197
- 在线: 78.3小时
- 虫号: 4191967
- 注册: 2015-11-02
- 性别: MM
- 专业: 工程热物理与能源利用
7楼2016-03-06 22:33:39
chibang367
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 863.2
- 红花: 1
- 帖子: 197
- 在线: 78.3小时
- 虫号: 4191967
- 注册: 2015-11-02
- 性别: MM
- 专业: 工程热物理与能源利用
8楼2016-03-06 23:38:00
50