| 查看: 596 | 回复: 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 ] |
回复此楼» 猜你喜欢
青C资助名额大幅增加！
已经有13人回复
-
风把牡丹吹跑了
已经有7人回复
-
教学课件你会给同学吗
已经有8人回复
-
26应届毕业生考博求助
已经有3人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有7人回复
-
求助2，4-二氯-5-嘧啶甲醛的合成方法
已经有5人回复
-
本人最近太闲了，谁有问题可以提，每天会统一回复
已经有19人回复
-
申博自荐
已经有6人回复
-
这年头没有找到涵评专家，还有中面上的可能吗
已经有8人回复
-
材料类只有一篇综述能申博么
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
20
精华III: