24小时热门版块排行榜    

查看: 991  |  回复: 0

河間

新虫 (初入文坛)

[求助] 【GULP错误报告】Coulomb subtraction not allowed

本人博一,岩土专业,在做Na-montmorillonite的分子模拟,研究hydration和力学性质。正常来讲,使用molecule关键词就可以去掉共价键原子间的库伦力。然而却遇到以下错误信息:

ERROR : Coulomb subtraction within periodic molecules is not allowed

体系主要状况是:cell里面有820个原子,包含粘土的晶体和水层以及水中的正离子。以下是.gin文件的主要部分:
# Keywords:
opti nosym full conp distance angle bond compare prop
# Structure(s):
cell
  20.692124  17.912804  24.272009  91.882561  98.519024  89.959525
fractional
。。。原子位置,种类,lj势函数的参数略。。。
# Molecular dynamic parameters
temperature    300.000
ensemble npt     0.1    0.1
timestep           0.001 ps
equilibration       2.0000 ps
production         8.0000 ps
sample             0.001 ps
write_MD           0.001 ps
lennard  epsilon geometric inter all       # 所有原子之间的lj势
0.0 9.0                                    # lj势的范围
harmonic bond intra                           # 水分子原子间的弹簧势
Oh H 48.0591  1
Hs H 48.0591  1
Ow Hw 48.0591  1
three bond intra regular
Ow    Hw  Hw  3.969515  109.47            # 水分子键角相关的势
cutp   9.0000   1.0000                     # 势函数cutoff radius
rtol 1.3                                   # bond length tolerance
switch_min rfo   gnorm     0.005000
。。。文件输出啥的略。。。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 河間 的主题更新
信息提示
请填处理意见