24小时热门版块排行榜    

查看: 1365  |  回复: 3
【悬赏金币】回答本帖问题,作者fyc801将赠送您 50 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fyc801

金虫 (正式写手)


[求助] Core/shell partner atom not found是什么原因?

大家好!我用lammps计算BaTiO3材料,用的core/shell模型。参考了手册的in文件写法,出现了错误“Core/shell partner atom not found”。手册的解释是Could not find one of the atoms in the bond pair. 我一直没想明白是哪里出了问题。寻求大家的帮助,十分感激。

data文件是

LAMMPS Description

        1250 atoms
         625 bonds

6 atom types
3 bond types

0   20.0500011445000      xlo xhi
0   20.0500011445000      ylo yhi
0   20.0500011445000      zlo zhi

Masses

1          47.867
2          0.0001
3          137.327
4          0.0001
5          15.9994
6          0.0001

Atoms

    1    1    1     4.7600     2.0050     2.0050     2.0050
    2    1    2    -1.5800     2.5050     2.5050     2.5050
    3    2    3     5.6200     0.0000     0.0000     0.0000
    4    2    4    -3.7600     0.5000     0.5000     0.5000
    5    3    5     0.9100     0.0000     2.0050     2.0050
    6    3    6    -2.5900     0.5000     2.5050     2.5050
(...)


Bonds

      1      1      1      2
      2      2      3      4
      3      3      5      6
      4      3      7      8
      5      3      9     10
(...)





in文件是

units                real
atom_style        full
boundary         p p p

#read_restart    restart.6523
read_data      BaTiO3-5x5x5.dat
run_style        verlet

group cores type 1 3 5
group shells type 2 4 6

compute CSequ all temp/cs cores shells
fix thermoberendsen all temp/berendsen 300 300 100   
fix thermostatequ all nve
fix_modify thermoberendsen temp CSequ


comm_modify vel yes
velocity all create 300.0 134 bias yes temp CSequ
velocity all scale 300.0 temp CSequ

pair_style     buck/coul/long/cs 10.0
pair_modify   shift yes
pair_coeff    * * 0.0 1.0 0.0
pair_coeff    2 6 104378.5 0.25239 0.0 #Ti-O
pair_coeff    4 6 19935.1 0.38729 0.0 #Ba-O  
pair_coeff    6 6 94604.04 0.29581 300.0 #O-O  

kspace_style pppm 1.0e-4   

bond_style  harmonic
bond_coeff      1 7401.852  0.0
bond_coeff      2 5787.741  0.0
bond_coeff      3 2335.158  0.0



thermo_style    custom  step  vol press pe  epair ebond ecoul  enthalpy temp
thermo_modify temp CSequ
thermo          100000
timestep        1




fix             1 all nvt temp 300.0 300.0 100 drag 2.0
dump            1 all xyz 100000 all.xyz
run             1000000
unfix  1
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

54疯癫

铁虫 (初入文坛)

壳核间距太远了?
4楼2022-05-03 13:25:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

好学又懒惰

新虫 (小有名气)

请问你的手册从哪下载的

发自小木虫Android客户端
2楼2021-04-17 21:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ch0722

银虫 (初入文坛)

请问楼主的问题解决了吗?我也想用lammps做BaTiO3材料,建模问题想请教一下楼主
3楼2021-04-27 19:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 0831生医工307分 求调剂 +3 小小怪wx 2026-03-08 3/150 2026-03-14 03:36 by JourneyLucky
[考研] 一志愿中国科学技术大学,086000,生物与医药,求调剂 +3 12341234789 2026-03-08 4/200 2026-03-14 03:14 by JourneyLucky
[考研] 296求调剂 +5 Xinyu Wu311 2026-03-09 5/250 2026-03-14 03:05 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-09 6/300 2026-03-14 02:14 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 泣血叩求调剂恩,愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 285求调剂 +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 304求调剂(085602一志愿985) +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[基金申请] 面上项目还需要AI说明吗? +3 liyundong 2026-03-08 3/150 2026-03-09 22:30 by kingkocxr
信息提示
请填处理意见