24小时热门版块排行榜    

查看: 2543  |  回复: 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

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 飞利白 at 2015-04-11 22:02:52
不好意思,我看错报错原因了,是bond不是pair,为什么你的bond_coeff      3 144.821 1.149键上来就是type 3了?bond类型不是1 2  3,我没用过这个力场,但是你要是键写出来就得设置bond_coeff吧,没有系数就直接不 ...

因为另外三个键我觉得都是包括在AIREBO中的,不知道我理解的有没有问题,因为官方帮助里面有这么一句话Only a single pair_coeff command is used with the airebo or rebo style which specifies an AIREBO potential
file with parameters for C and H.
因为AIREBO有三项相加,其中第一项就是短程作用力,我就理解为包含了键的作用。不知道对不对
6楼2015-04-11 23:35:01
已阅   回复此楼   关注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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中! +5 恐龙爸爸 2024-06-14 5/250 2024-06-15 20:36 by wanghuawei
[教师之家] 我们学院常年位居 各学院 倒数第二。专业撤销的话,在编者有什么补偿? +13 河西夜郎 2024-06-09 14/700 2024-06-15 19:44 by LittleBush
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[考博] 2025考博 +4 自强不息a?a 2024-06-15 6/300 2024-06-15 12:36 by 自强不息a?a
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +9 味道很好啊 2024-06-13 15/750 2024-06-15 09:59 by yulan19992002
[基金申请] 面上基金有一个(两个)C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[硕博家园] 关于硕博连读的一些疑问? +4 Lwenter 2024-06-14 4/200 2024-06-14 14:32 by ou0551
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[找工作] 应届大龄土博教职选择,石河子大学和广西大学 10+4 茅山老道士 2024-06-10 31/1550 2024-06-14 10:40 by ChmistryCHEN
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 最近写了一篇控制优化领域的文章,可以投哪里啊?有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[硕博家园] 机械研究生如何拿到年薪40+w +13 阿巴阿巴哦哦 2024-06-11 15/750 2024-06-13 15:40 by 113745685
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[找工作] 成都产品质量检测研究院 200+3 鲸鱼663 2024-06-11 8/400 2024-06-13 08:10 by 加纳居士
[基金申请] 博后特助这周出结果吗?往年都是啥时候啊? +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
[硕博家园] 申博 +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
信息提示
请填处理意见