| 查看: 1379 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者fyc801将赠送您 50 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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 |
回复此楼» 猜你喜欢
26药学专硕105500求调剂
已经有8人回复
-
297，工科调剂?
已经有11人回复
-
申博自荐
已经有4人回复
-
294求调剂
已经有9人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有5人回复
-
291求调剂
已经有11人回复
-
304求调剂
已经有8人回复
-
求调剂推荐
已经有9人回复
-
求调剂
已经有12人回复
-
291求调剂
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
ch0722
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 356.5
- 帖子: 42
- 在线: 94.1小时
- 虫号: 3196205
- 注册: 2014-05-10
- 性别: MM
- 专业: 金属材料的微观结构
3楼2021-04-27 19:02:04
好学又懒惰
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 596.4
- 散金: 10
- 帖子: 176
- 在线: 10.5小时
- 虫号: 23176721
- 注册: 2020-07-17
- 专业: 理论和计算化学
2楼2021-04-17 21:15:11
4楼2022-05-03 13:25:27
50