| 查看: 660 | 回复: 0 | ||
zzhwise金虫 (小有名气)
|
[求助]
有机分子与无机Fe基板模型无法运行！
|
|
各位大神： 我的模型是2块铁基板之间夹着长链有机分子. 为了测试模型，我分别对基板和有机分子进行了模拟 基板：利用in文件进行建模，计算，顺利通过 有机分子：自己写的data文件，计算，顺利通过 但是当自己讲Fe坐标写进有机分子的data再运行时，无法运行。具体文件如下： in文件 units metal dimension 3 boundary p p p #newton on atom_style full bond_style harmonic angle_style harmonic read_data data.fec24 pair_style hybrid eam/fs lj/cut 10.0 pair_coeff * * eam/fs Fe.eam.fs NULL NULL NULL Fe #只有铁原子间用eam势 pair_coeff 1 1 lj/cut 0.060 3.5 pair_coeff 1 2 lj/cut 0.060 3.5 pair_coeff 1 3 lj/cut 0.044 3.0 #使用LJ/COUL势 pair_coeff 1 4 lj/cut 0.249 3.0 pair_coeff 2 2 lj/cut 0.060 3.5 #不是使用杂化势时，不需要LJ/COUL pair_coeff 2 3 lj/cut 0.044 3.0 pair_coeff 2 4 lj/cut 0.249 3.0 pair_coeff 3 3 lj/cut 0.290 2.5 pair_coeff 3 4 lj/cut 0.165 2.5 group up molecule 122 group down molecule 121 neighbor 1.0 bin neigh_modify delay 5 min_style sd minimize 1.0e-4 1.0e-4 1000 1000 dump 1 all xyz 10 fec24.xyz # fix 1 all nvt temp 275.0 275.0 100.0 tchain 1 run 300 data文件 LAMMPS Description 20406 atoms 8760 bonds 17280 angles 0 dihedrals 0 impropers 4 atom types 2 bond types 3 angle types 0 dihedral types 0 improper types -5.0000 65.0000 xlo xhi -5.0000 50.0000 ylo yhi -35.0000 60.0000 zlo zhi Masses 1 12.0110 2 12.0110 3 1.00790 4 55.8450 Bond Coeffs 1 309.0 1.526 2 338.7 1.105 Angle Coeffs 1 44.1000 110.00000 2 46.300 110.50000 3 39.260 106.50000 Atoms 1 1 1 -0.18 24.109 30.702 1.682 2 1 2 -0.12 23.701 29.259 1.490 3 1 2 -0.12 24.912 28.438 1.784 4 1 2 -0.12 24.637 26.976 1.550 5 1 2 -0.12 26.005 26.295 1.739 6 1 2 -0.12 25.998 24.769 1.643 7 1 2 -0.12 27.032 23.634 2.129 。。。 20400 122 4 0.00 45.862 48.729 -14.8020 20401 122 4 0.00 48.729 48.729 -14.8020 20402 122 4 0.00 51.595 48.729 -14.8020 20403 122 4 0.00 54.462 48.729 -14.8020 20404 122 4 0.00 57.328 48.729 -14.8020 20405 122 4 0.00 60.194 48.729 -14.8020 20406 122 4 0.00 63.061 48.729 -14.8020 Bonds 1 2 8830 8880 2 2 8830 8879 data里没有Fe原子的成键对 只有有机分子的键长和键角信息 没有非键pair 请问有遇到这样问题的么？ [ Last edited by zzhwise on 2013-7-23 at 11:20 ] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
20