| 查看: 947 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者brickxj将赠送您 10 个金币 | ||
[求助]
LAMMPS 混合势能函数设置问题（二氧化碳水合物）
|
||
|
水分子选择TIP4P/ICE，二氧化碳分子选择TraPPE势能模型，请问在in文件的编写中，应该怎么写呢？ 如果均选择长程静电力的话，水分子应该用pppm/tip4p，而二氧化碳分子用pppm。会出现矛盾，程序无法运行。 如果用短程静电力来计算，结果明显不对。 应该怎么处理呢？ 以下是我的势能函数设置，结果是没有报错，但是程序算不动。 pair_style hybrid lj/cut/tip4p/cut 1 2 1 1 0.1546 10 lj/cut/coul/long 10 pair_modify tail yes pair_coeff 1 1 lj/cut/tip4p/cut 0.210838 3.1668 #水分子势能函数tip4p/ICE pair_coeff 2 2 lj/cut/tip4p/cut 0.0000 0.00 pair_coeff 1 2 lj/cut/tip4p/cut 0 0 pair_coeff 3 3 lj/cut/coul/long 0.053653 2.8 #二氧化碳分子势能函数EPM2 pair_coeff 4 4 lj/cut/coul/long 0.156986 3.05 pair_coeff 1 3 lj/cut/coul/long 0.106359 2.9834 pair_coeff 1 4 lj/cut/coul/long 0.196485 3.1084 pair_coeff 2 3 lj/cut/coul/long 0 0 pair_coeff 2 4 lj/cut/coul/long 0 0 pair_coeff 3 4 lj/cut/coul/long 0.091776 2.925 kspace_style pppm 1.0e-6 bond_style harmonic bond_coeff 1 0 0.9572 bond_coeff 2 0 1.16 angle_style harmonic angle_coeff 1 0 104.5 angle_coeff 2 0 180 哪位大神能够帮忙看看呀，不胜感激~ |
回复此楼» 猜你喜欢
面上提前没消息，有中的吗
已经有14人回复
-
这种情况还有戏吗
已经有12人回复
-
微信指数没变化，科研之友没阅读
已经有19人回复
-
准备明年的基金了
已经有7人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有24人回复
-
2026年国自然面上资助率
已经有21人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有10人回复
-
2027年申博
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有7人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
