计算化学公社

标题: 请问LAMMPS软件使用 Buckingham Potential应该怎么设置atom_style [打印本页]

作者
Author: vaspMeng 时间: 2023-5-18 21:16
标题: 请问LAMMPS软件使用 Buckingham Potential应该怎么设置atom_style
我在使用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

作者
Author: likejade 时间: 2023-5-19 07:33
ERROR: Pair style buck/coul/long requires atom attribute q 表明data文件中没有提供原子的电荷，问题在data 文件， atom style可以用charge 或者 full

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

作者
Author: vaspMeng 时间: 2023-8-4 19:37

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

作者
Author: thor 时间: 2024-7-22 16:17

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

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

作者
Author: crychichi 时间: 3 day ago

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)

欢迎光临计算化学公社 (http://bbs.keinsci.com/)