24小时热门版块排行榜    

查看: 2545  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 硕博巨婴,也许才刚刚开始 +30 SNaiL1995 2024-06-12 65/3250 2024-06-16 02:10 by 青古
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] 面青地会评时间 +5 tanjydd 2024-06-15 5/250 2024-06-15 23:08 by andywei1028
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中! +5 恐龙爸爸 2024-06-14 5/250 2024-06-15 20:36 by wanghuawei
[教师之家] 请问事业编制和年薪制冲突吗? +6 ZHONGWU_U 2024-06-14 6/300 2024-06-15 20:16 by Ermito
[基金申请] 博后基金,博管会会提前知道消息吗? +3 yuyiang 2024-06-13 3/150 2024-06-15 19:52 by Lynn212
[教师之家] 我们学院常年位居 各学院 倒数第二。专业撤销的话,在编者有什么补偿? +13 河西夜郎 2024-06-09 14/700 2024-06-15 19:44 by LittleBush
[基金申请] 2024国社科通讯评审 +9 qsd10086 2024-06-13 14/700 2024-06-15 15:51 by thesuna
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +8 zju2000 2024-06-12 14/700 2024-06-15 13:59 by chemhua
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +9 味道很好啊 2024-06-13 15/750 2024-06-15 09:59 by yulan19992002
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[基金申请] 面上基金有一个(两个)C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[论文投稿] 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
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[基金申请] 连续两年医学口青年项目初审体会 +11 进击的荣耀 2024-06-09 18/900 2024-06-13 17:27 by 进击的荣耀
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[硕博家园] 求助 +6 LYWwrz 2024-06-09 9/450 2024-06-11 13:12 by powerhours
[教师之家] 公办双非,学生论文升学就业都不行。一本大学的论文升学就业没问题吧? +5 河西夜郎 2024-06-09 5/250 2024-06-10 17:59 by yyallen2003
信息提示
请填处理意见