| 查看: 1177 | 回复: 1 | |||
[交流]
hbond/dreiding/lj相关计算
|
|
先运行oh.in，再运行ohre.in进行rerun。研究ohre.in和结果。 ohre.in文件如下： dimension 3 units metal boundary p p p atom_style full neighbor 3.0 bin neigh_modify delay 0 every 1 check yes read_data coh.data bond_style harmonic bond_coeff 1 13.878 1.410 bond_coeff 2 23.983 0.96 bond_coeff 3 13.878 1.410 bond_coeff 4 23.983 0.96 angle_style harmonic angle_coeff 1 2.1684 109.5 angle_coeff 2 2.3853 108.5 angle_coeff 3 2.1684 109.5 angle_coeff 4 2.3853 108.5 dihedral_style opls dihedral_coeff 1 0.1160 -0.1250 0.0445 0.0 dihedral_coeff 2 0.1160 -0.1250 0.0445 0.0 #dielectric 1.0 #kspace_style pppm 1.0e-6 pair_style hybrid/overlay airebo 2.5 1 1 lj/cut 10.0 hbond/dreiding/lj 4 3.3 3.5 150 pair_coeff * * airebo CH.airebo C C NULL NULL NULL NULL pair_coeff 1 2 none pair_coeff 1 2 lj/cut 0 3.50 pair_coeff 1 3 lj/cut 0.0046 3.31 pair_coeff 1 4 lj/cut 0 3.31 pair_coeff 2 3 lj/cut 0 3.31 pair_coeff 2 4 lj/cut 0.0046 3.31 pair_coeff 1 5 lj/cut 0.0016 2.1954 pair_coeff 1 6 lj/cut 0 2.1954 pair_coeff 2 5 lj/cut 0 2.1954 pair_coeff 2 6 lj/cut 0.0016 2.1954 pair_coeff 3 3 lj/cut 0.0074 3.12 pair_coeff 4 4 lj/cut 0.0074 3.12 pair_coeff 3 4 lj/cut 0 3.12 pair_coeff 3 5 lj/cut 0.0025 2.0054 pair_coeff 3 6 lj/cut 0 2.0054 pair_coeff 4 5 lj/cut 0 2.0054 pair_coeff 4 6 lj/cut 0.0025 2.0054 pair_coeff 5 5 lj/cut 8.6738e-4 0.8909 pair_coeff 6 6 lj/cut 8.6738e-4 0.8909 pair_coeff 5 6 lj/cut 0 0.8909 pair_coeff 3 4 hbond/dreiding/lj 5 i 0.412 2.75 4 3.3 3.5 150 pair_coeff 3 4 hbond/dreiding/lj 6 j 0.412 2.75 4 3.3 3.5 150 region r1 block INF INF INF INF INF 0.3 units box region rt block INF INF INF INF 0.3 INF units box group group2 region r1 group group1 region rt variable fg1x equal fcm(group1,x) compute hb all pair hbond/dreiding/lj variable nhb equal c_hb[1] variable ehb equal c_hb[2] compute gg2 group1 group/group group2 timestep 0.001 thermo_style custom v_fg1x c_gg2[1] c_hb[2] c_gg2 thermo 1 rerun coh.dump dump x y z vx vy vz ix iy iz box no 1、group1 包含原子类型 1,3,5。group2 包含原子类型 2,4,6。 2、bond/angle/dihedral 是 group1 或 group2 内部作用。 3、group1 和 group2 之间只有 hbond/dreiding/lj 作用。 4、无外力。 5、variable fg1x equal fcm(group1,x) 和 compute gg2 group1 group/group group2 的结果为什么不一样？ compute hb all pair hbond/dreiding/lj 和 compute gg2 group1 group/group group2 得到的能量有什么关系？是一部分吗？ fg1x gg2[1] hb[2] gg2 -0.26078056 0.068249515 -1.1930503 -0.78975531 -0.9591668 0.047818212 -1.4279895 -0.80196775 -0.98344781 0.029496808 -1.4659098 -0.80584598 -0.97878567 0.013287859 -1.4980585 -0.80407161 -0.93969166 0.019469678 -1.5237638 -0.80046352 -0.86897216 0.022923407 -1.5414305 -0.79845355 -0.77782396 0.021246764 -1.549718 -0.79996216 -0.68174755 0.015051052 -1.5484882 -0.80483706 -0.59494213 0.0049977361 -1.5388899 -0.81090889 -0.52521538 -0.0072240828 -1.5229889 -0.81482867 -0.4721573 -0.015303847 -1.5035062 -0.81353129 6、把 pair_coeff 3 4 hbond/dreiding/lj 5 i 0.412 2.75 4 3.3 3.5 150 注释掉， 结果是有氢键作用，但group/group结果为0。这是为什么？ fg1x gg2[1] hb[2] gg2 -0.19016095 0 -0.40329497 0 -0.88982777 0 -0.62602174 0 -0.91450677 0 -0.66006377 0 -0.90743726 0 -0.69398685 0 -0.86298851 0 -0.72330024 0 -0.78508871 0 -0.74297691 0 -0.6863355 0 -0.7497558 0 -0.58321373 0 -0.74365119 0 -0.49026671 0 -0.72798105 0 -0.41523376 0 -0.70816023 0 -0.35762304 0 -0.6899749 0 请大家指点。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : CH.airebo
- 附件 2 : coh.data
2014-08-25 19:22:52, 904.84 K
2014-08-25 19:22:52, 12.76 K
» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有24人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有18人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Ubuntu 编译Lammps 单机串行错误，帮帮忙
已经有5人回复
2楼2014-09-01 11:00:55
