请问LAMMPS软件使用 Buckingham Potential应该怎么设置atom_style

vaspMeng

我在使用buckingham势进行sio2的计算时候，出现报错ERROR: Pair style buck/coul/long requires atom attribute q (../pair_buck_coul_long.cpp:361)，查阅手册说The atom style defined does not have these attributes. 不知道是不是我的atom_style设置的不正确，求助各位老师，万分感谢。
我使用的in文件如下：
units      metal
atom_style    atomic
boundary     p p p
timestep     0.001
neighbor     2 bin
neigh_modify   every 1 delay 0 check yes
#读取data⽂件
read_data    sio2.data
#势函数设置
pair_style    tersoff
pair_coeff    * * SiO.tersoff O Si
#温度初始化
velocity all create 7000 98989
#⾼温弛豫
dump 1 all atom 1000 dump.xyz
fix 1 all npt temp 7000 7000 0.1 iso 100 100 1
run 50000
unfix 1
#降温
fix 1 all npt temp 7000 300 0.1 iso 1 1 1
run 50000
unfix 1
#300K弛豫
fix 1 all nvt temp 300 300 0.1
run 50000
#保存⾮晶⽂件
write_data sio2_a.data

likejade

ERROR: Pair style buck/coul/long requires atom attribute q  表明data文件中没有提供原子的电荷，问题在data 文件， atom style可以用charge 或者 full

另外，pair_style    tersoff 用的是tersoff势函数，并不是buckingham势

vaspMeng

likejade 发表于 2023-5-19 07:33
ERROR: Pair style buck/coul/long requires atom attribute q  表明data文件中没有提供原子的电荷，问题在 ...

非常感谢您的回复，另外还有一个问题想请教您，为什么我看到论文里使用buckingham势的时候，都没有给出同种原子的参数呢，如SiO2，只给了Si-O，O-O却没有Si-Si

thor

vaspMeng 发表于 2023-8-4 19:37
非常感谢您的回复，另外还有一个问题想请教您，为什么我看到论文里使用buckingham势的时候，都没有给出同 ...

Si-Si的是0 ，可以设置0 0.1 0

crychichi

likejade 发表于 2023-5-19 07:33
ERROR: Pair style buck/coul/long requires atom attribute q  表明data文件中没有提供原子的电荷，问题在 ...

您好请问我在使用buckingham势进行sio2的计算时候，文献里是qSi=2.4e, qO=−1.2e，但是我data文件里电荷都是0，这个应该怎么更改呢？报错是ERROR: Lost atoms: original 648 current 35 (src/thermo.cpp:481)