24小时热门版块排行榜    

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

wood1130

新虫 (初入文坛)

[求助] 初学者用LAMMPS做气体穿过多孔石墨烯遇到报错已有2人参与

各位大神我毕设做多孔石墨烯的气体分离模拟,由于包括CH4,准备将CH4和石墨烯碳原子这三种原子用AIREBO,其他原子用LJ势,在IN文件中pair style选用了hybrid,请问使用了AIREBO的C和H是不是不用再设置bond coeff和angle coeff等,只用设置pair coeff?还有我运算时总是报错,说我没有设置bond coeff,请各位大神帮助我一下谢谢大家!,我把我的in文件放这里,其中CA是石墨烯周边C,CS是将要挖成孔的C,CO是CO2的C,OT是氧原子。目前里面只有CO2和石墨烯,石墨烯CC间用AIREBO,其他用LJ,可能具体参数不对大家不要介意,只是在试是否能计算
units           metal
boundary        p p p
#communicate     single vel yes

atom_style      full
bond_style      harmonic
angle_style     harmonic
pair_style      hybrid airebo 3.0 1 0 lj/cut/coul/long 14.0 10.0
kspace_style    pppm 1.0e-4

read_data        co2merge.data

pair_coeff      * *  airebo ../potentials/CH.airebo C NULL C NULL
pair_coeff      1 2  lj/cut/coul/long 0.0557 3.079
pair_coeff      2 3  lj/cut/coul/long 0.0557 3.079
pair_coeff      2 4  lj/cut/coul/long 0.0945 2.895
pair_coeff      2 2  lj/cut/coul/long 0.0559 2.757
pair_coeff      1 4  lj/cut/coul/long 0.1102 3.292
pair_coeff      4 4  lj/cut/coul/long 0.1599 3.033
pair_coeff      3 4  lj/cut/coul/long 0.1102 3.292

bond_coeff      3 144.821 1.149
angle_coeff     5 0.045 180.0

group                 hole   type 3
group                 mem   type 1
group                 co2   type 2 4

neighbor        2.0 bin
neigh_modify    delay 5
neigh_modify    exclude type 1 1
neigh_modify    exclude type 1 3
fix             1 hole setforce 0.0 0.0 0.0
fix             2 mem setforce 0.0 0.0 0.0dump            dumpwall all dcd 10000000 dump.co2.dcd
dump_modify     dumpwall flush yes
thermo_style    custom step cpu temp press etotal ke
thermo          100
#thermo_modify   flush yes
timestep        2.0
velocity        co2 create 300.0 34239
run             60000
write_restart   restart.eq
日志文件如下:
LAMMPS (14 Mar 2015-ICMS)
WARNING: OMP_NUM_THREADS environment is not set. (../comm.cpp:89)
  using 1 OpenMP thread(s) per MPI task
Reading data file ...
  orthogonal box = (0 0 -100) to (50 50 100)
  1 by 1 by 1 MPI processor grid
  reading atoms ...
  1200 atoms
  scanning bonds ...
  2 = max bonds/atom
  scanning angles ...
  3 = max angles/atom
  scanning dihedrals ...
  8 = max dihedrals/atom
  scanning impropers ...
  3 = max impropers/atom
  reading bonds ...
  1595 bonds
  reading angles ...
  2910 angles
  reading dihedrals ...
  5541 dihedrals
  reading impropers ...
  920 impropers
Finding 1-2 1-3 1-4 neighbors ...
  3 = max # of 1-2 neighbors
  6 = max # of 1-3 neighbors
  18 = max # of 1-4 neighbors
  18 = max # of special neighbors
Reading potential file ../potentials/CH.airebo with DATE: 2011-10-25
10 atoms in group hole
998 atoms in group mem
192 atoms in group co2
AAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAA
PPPM initialization ...
  G vector (1/distance) = 0.168083
  grid = 10 10 32
  stencil order = 5
  estimated absolute RMS force accuracy = 0.00126951
  estimated relative force accuracy = 8.81624e-005
  using double precision FFTs
  3d grid and FFT values/proc = 8325 3200
WARNING: Using a manybody potential with bonds/angles/dihedrals and special_bond exclusions (../pair.cpp:207)
ERROR: All bond coeffs are not set (../bond.cpp:69)
DATA文件主要部分:
1200 atoms
1595 bonds
2910 angles
5541 dihedrals
920 impropers
4 atom types
4 bond types
5 angle types
6 dihedral types
4 improper types
0.00000 50.00000  xlo xhi
0.00000 50.00000  ylo yhi
-100.00000 100.00000  zlo zhi

# Pair Coeffs
#
# 1  CA
# 2  CO
# 3  CS
# 4  OT

# Bond Coeffs
#
# 1  CA-CA
# 2  CA-CS
# 3  CO-OT
# 4  CS-CS

# Angle Coeffs
#
# 1  CA-CA-CA
# 2  CA-CA-CS
# 3  CA-CS-CS
# 4  CS-CS-CS
# 5  OT-CO-OT

# Dihedral Coeffs
#
# 1  CA-CA-CA-CA
# 2  CA-CA-CA-CS
# 3  CA-CA-CS-CS
# 4  CA-CS-CS-CA
# 5  CA-CS-CS-CS
# 6  CS-CS-CS-CS

# Improper Coeffs
#
# 1  CA-CA-CA-CA
# 2  CA-CA-CA-CS
# 3  CA-CS-CS-CS
# 4  CS-CS-CS-CS
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wood1130

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 飞利白 at 2015-04-11 21:48:31
你的pair_coeff 1 3  没有设置,你的1和3分别使用airebo,那1类和3类碳原子之间也是吗,你先试试加上pair_coeff    1 3   lj/cut/coul/long 0.0557 3.079看看还报错不,如果不报错那就是1和3之间的没有设置

可是1 3之间难道不是也应该适用AIREBO吗,AIREBO不是包括了长程力难道不是任意两两原子之间的力吗,为什么还要设置成LJ呢
5楼2015-04-11 23:30:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

飞利白

新虫 (初入文坛)


wood1130(abinitio代发): 金币+1, 鼓励交流 2015-04-13 11:01:10
你的pair_coeff 1 3  没有设置,你的1和3分别使用airebo,那1类和3类碳原子之间也是吗,你先试试加上pair_coeff    1 3   lj/cut/coul/long 0.0557 3.079看看还报错不,如果不报错那就是1和3之间的没有设置
2楼2015-04-11 21:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞利白

新虫 (初入文坛)

不好意思,我看错报错原因了,是bond不是pair,为什么你的bond_coeff      3 144.821 1.149键上来就是type 3了?bond类型不是1 2  3,我没用过这个力场,但是你要是键写出来就得设置bond_coeff吧,没有系数就直接不写了。还是得看看你那个力场AIREBO,
3楼2015-04-11 22:02:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wood1130

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 飞利白 at 2015-04-11 21:48:31
你的pair_coeff 1 3  没有设置,你的1和3分别使用airebo,那1类和3类碳原子之间也是吗,你先试试加上pair_coeff    1 3   lj/cut/coul/long 0.0557 3.079看看还报错不,如果不报错那就是1和3之间的没有设置

可是1 3之间难道不是也应该适用AIREBO
4楼2015-04-11 23:28:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[教师之家] 关于2023的收入 +33 小龙虾2008 2024-06-10 34/1700 2024-06-15 23:01 by zeolitess
[教师之家] 每次骚扰女学生的都是院系领导,而不是普通教师,小编们要注意措辞正确 +7 zju2000 2024-06-15 9/450 2024-06-15 22:54 by zeolitess
[教师之家] 请问事业编制和年薪制冲突吗? +6 ZHONGWU_U 2024-06-14 6/300 2024-06-15 20:16 by Ermito
[教师之家] 我们学院常年位居 各学院 倒数第二。专业撤销的话,在编者有什么补偿? +13 河西夜郎 2024-06-09 14/700 2024-06-15 19:44 by LittleBush
[基金申请] 博士后创新人才支持计划公示 +9 aishida144 2024-06-14 13/650 2024-06-15 19:18 by 默默挥手
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +3 第二种态度 2024-06-11 5/250 2024-06-15 17:56 by 小懂事k
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了,能确定送审了吗? +5 潇洒怡惜 2024-06-12 10/500 2024-06-14 09:56 by 潇洒怡惜
[考博] 还能读上博士吗 +14 小伙走起 2024-06-09 24/1200 2024-06-14 06:32 by 18787065176
[有机交流] 旋蒸请教 +6 办公室摸鱼 2024-06-09 9/450 2024-06-13 21:31 by wuxiaominggo
[基金申请] 连续两年医学口青年项目初审体会 +11 进击的荣耀 2024-06-09 18/900 2024-06-13 17:27 by 进击的荣耀
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[基金申请] C口食品 +6 好好不过分 2024-06-09 19/950 2024-06-13 12:26 by zydemon525
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
[硕博家园] 求助 +6 LYWwrz 2024-06-09 9/450 2024-06-11 13:12 by powerhours
[论文投稿] water research状态咨询 5+3 Flyyawa 2024-06-10 6/300 2024-06-11 09:45 by bobvan
信息提示
请填处理意见