计算化学公社

 找回密码 Forget password
 注册 Register
Views: 336|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] LAMMPS做GCMC吸水求助

[复制链接 Copy URL]

5

帖子

0

威望

105

eV
积分
110

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好,我目前在用lammps软件建造水化硅酸钙(CSH)模型,当前在GCMC吸水部分出现了问题,他一直在报错ERROR: All bond coeffs are not set (src/bond.cpp:84),我有点解决不了了,所以想请教一下各位老师,请大家指点一下,谢谢大家了
附页中是我的数据文件和日志文件,下面是我写的一个in文件:#C-S-H model
units                         real   
dimension           3                                                                  
boundary                p p p                                                                                    
atom_style                full                                                                                   
neighbor            2.0 nsq                                                                                                  
neigh_modify        every 1 delay 0 check yes  

pair_style          lj/cut/coul/long 10.0
pair_modify         mix arithmetic
bond_style          harmonic
angle_style         harmonic  
dihedral_style      harmonic              
kspace_style        ewald 0.0001

read_data           C-S-H-1.7.data  extra/atom/types 2  extra/bond/types 4 extra/angle/types 2   
molecule            155 H2O.txt

mass                1  40.080000 # cah
mass                2  28.085500 # st
mass                3  15.999400 # ob
mass                4   15.999  # OW      
mass                5   1.008   # HW

pair_coeff          4  4  0.1554000    3.16557      28.0   # OW      
pair_coeff          5  5  0.00000e+00  0.00000e+00  28.0   # HW  

group kln type 1:3
group water type 4:5
group molwater  molecule 155

fix shakew water shake 0.0001 20 5 b 1 a 1 mol 155
neigh_modify exclude group kln kln
compute_modify thermo_temp dynamic yes
thermo 100
fix 1 kln setforce 0.0 0.0 0.0

fix 2 molwater gcmc 10 1000 1000  0 2949499 300.0 -0.5 0.5 mol 155  maxangle 180 full_energy
thermo_style  custom step atoms temp press pe ke etotal density
run 300000



GCMC.in

1.57 KB, 下载次数 Times of downloads: 8

C-S-H-1.7.data

88.44 KB, 下载次数 Times of downloads: 8

H2O.txt

535 Bytes, 下载次数 Times of downloads: 4

log.lammps

3.22 KB, 下载次数 Times of downloads: 1

3

帖子

0

威望

101

eV
积分
104

Level 2 能力者

2#
发表于 Post on 2025-6-20 17:44:26 | 只看该作者 Only view this author
你水分子的bond coeffs没有设置

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-13 09:03 , Processed in 0.365647 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list