| 查看: 513 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者hzw82623将赠送您 15 个金币 | ||
[求助]
锂粒子碰撞c60，c60不稳定总是自己先炸了是什么原因？
|
||
|
我在做一个锂离子碰撞c60的模拟，应该是锂离子撞上c60后会与c结合，但我的c60自己因为不稳定就炸开了，锂离子还没接触。这是哪里出了问题呢？是我的势函数还是系综有问题？还是说我没有定义bond、二面角这些导致的？ in文件是这样 echo screen units metal dimension 3 boundary p p p atom_style full #bond_style harmonic read_data data.c60li pair_style hybrid airebo 3.0 lj/cut/coul/cut 10.0 8.0 yukawa 2.0 2.5 #pair_coeff 1 1 airebo CH.airebo C pair_coeff * * airebo CH.airebo C NULL pair_coeff 2 2 yukawa 0.05429457 0.03331 pair_coeff 1 2 lj/cut/coul/cut 2.0 3.5 9.0 #bond_coeff 5 80.0 1.2 dielectric 2.0 group 1 type 1 group 2 type 2 minimize 1.0e-4 1.0e-6 10000 100000 neighbor 3.0 bin neigh_modify every 1 delay 0 check yes dump 1 all atom 10 dump.lammpstrj velocity 1 set 0 0 0 velocity 2 set 0 -1 0 fix 1 1 nve run 30 fix 2 all nve timestep 0.001 run 15000 data文件 # LAMMPS data file written by OVITO 61 atoms 90 bonds 2 atom types 1 bond types -10.8150001 10.8150001 xlo xhi -10.989 20.7090001 ylo yhi -10.876 10.876 zlo zhi Masses 1 12.0110 2 6.941 Atoms # full 1 1 1 0.0 -0.732 -1.007 -3.332 2 1 1 0.0 0.732 -1.007 -3.332 3 1 1 0.0 1.184 0.385 -3.332 4 1 1 0.0 0.0 1.245 -3.332 5 1 1 0.0 -1.184 0.385 -3.332 6 1 1 0.0 -2.3039999 0.749 -2.6040001 7 1 1 0.0 -1.424 -1.96 -2.6040001 8 1 1 0.0 1.424 -1.96 -2.6040001 9 1 1 0.0 2.3039999 0.749 -2.6040001 10 1 1 0.0 0.0 2.4230001 -2.6040001 11 1 1 0.0 2.6080001 -1.576 -1.8350001 12 1 1 0.0 3.036 -0.259 -1.8350001 13 1 1 0.0 2.3039999 1.994 -1.8350001 14 1 1 0.0 1.184 2.8080001 -1.8350001 15 1 1 0.0 -1.184 2.8080001 -1.8350001 16 1 1 0.0 -2.3039999 1.994 -1.8350001 17 1 1 0.0 -3.036 -0.259 -1.8350001 18 1 1 0.0 -2.6080001 -1.576 -1.8350001 19 1 1 0.0 -0.693 -2.968 -1.8350001 20 1 1 0.0 0.692 -2.968 -1.8350001 21 1 1 0.0 -1.424 -3.2060001 -0.59 22 1 1 0.0 -2.6080001 -2.345 -0.59 23 1 1 0.0 1.424 -3.2060001 -0.59 24 1 1 0.0 2.6080001 -2.345 -0.59 25 1 1 0.0 3.4879999 0.364 -0.59 26 1 1 0.0 3.036 1.7560001 -0.59 27 1 1 0.0 0.732 3.431 -0.59 28 1 1 0.0 -0.732 3.431 -0.59 29 1 1 0.0 -3.036 1.7560001 -0.59 30 1 1 0.0 -3.4879999 0.364 -0.59 31 1 1 0.0 -0.732 -3.431 0.589 32 1 1 0.0 0.732 -3.431 0.589 33 1 1 0.0 3.036 -1.7560001 0.589 34 1 1 0.0 3.4879999 -0.364 0.589 35 1 1 0.0 2.6080001 2.345 0.589 36 1 1 0.0 1.424 3.2060001 0.589 37 1 1 0.0 -2.6080001 2.345 0.589 38 1 1 0.0 -1.424 3.2060001 0.589 39 1 1 0.0 -3.036 -1.7560001 0.589 40 1 1 0.0 -3.4879999 -0.364 0.589 41 1 1 0.0 1.184 -2.8080001 1.8350001 42 1 1 0.0 2.3039999 -1.994 1.8350001 43 1 1 0.0 -2.3039999 -1.994 1.8350001 44 1 1 0.0 -1.184 -2.8080001 1.8350001 45 1 1 0.0 -3.036 0.259 1.8350001 46 1 1 0.0 -2.6080001 1.576 1.8350001 47 1 1 0.0 -0.693 2.968 1.8350001 48 1 1 0.0 0.692 2.968 1.8350001 49 1 1 0.0 3.036 0.259 1.8350001 50 1 1 0.0 2.6080001 1.576 1.8350001 51 1 1 0.0 0.0 -2.4230001 2.6040001 52 1 1 0.0 -2.3039999 -0.749 2.6040001 53 1 1 0.0 -1.424 1.96 2.6040001 54 1 1 0.0 1.424 1.96 2.6040001 55 1 1 0.0 2.3039999 -0.749 2.6040001 56 1 1 0.0 -0.732 1.008 3.332 57 1 1 0.0 0.732 1.008 3.332 58 1 1 0.0 1.184 -0.385 3.332 59 1 1 0.0 0.0 -1.245 3.332 60 1 1 0.0 -1.184 -0.385 3.332 61 2 2 -1.0 -0.146 18.7090001 0.408 Bonds 1 1 1 2 2 1 1 5 3 1 1 7 4 1 2 3 5 1 2 8 6 1 3 4 7 1 3 9 8 1 4 5 9 1 4 10 10 1 5 6 11 1 6 16 12 1 6 17 13 1 7 18 14 1 7 19 15 1 8 11 16 1 8 20 17 1 9 12 18 1 9 13 19 1 10 14 20 1 10 15 21 1 11 12 22 1 11 24 23 1 12 25 24 1 13 14 25 1 13 26 26 1 14 27 27 1 15 16 28 1 15 28 29 1 16 29 30 1 17 18 31 1 17 30 32 1 18 22 33 1 19 20 34 1 19 21 35 1 20 23 36 1 21 22 37 1 21 31 38 1 22 39 39 1 23 24 40 1 23 32 41 1 24 33 42 1 25 26 43 1 25 34 44 1 26 35 45 1 27 28 46 1 27 36 47 1 28 38 48 1 29 30 49 1 29 37 50 1 30 40 51 1 31 32 52 1 31 44 53 1 32 41 54 1 33 34 55 1 33 42 56 1 34 49 57 1 35 36 58 1 35 50 59 1 36 48 60 1 37 38 61 1 37 46 62 1 38 47 63 1 39 40 64 1 39 43 65 1 40 45 66 1 41 42 67 1 41 51 68 1 42 55 69 1 43 44 70 1 43 52 71 1 44 51 72 1 45 46 73 1 45 52 74 1 46 53 75 1 47 48 76 1 47 53 77 1 48 54 78 1 49 50 79 1 49 55 80 1 50 54 81 1 51 59 82 1 52 60 83 1 53 56 84 1 54 57 85 1 55 58 86 1 56 57 87 1 56 60 88 1 57 58 89 1 58 59 90 1 59 60  [F4KJ(F]POPR)8GFFXKKELU.png  M_$DFG7)(77OBP`7XC(L6)F.png |
回复此楼» 猜你喜欢
投稿时忘记修改一作
已经有6人回复
-
竟然失眠了
已经有9人回复
-
硕博巨婴，也许才刚刚开始
已经有55人回复
-
34岁读博士晚吗
已经有20人回复
-
饶议：什么制度能保障大学普通教师不用为领导拎包，不用看领导脸色
已经有10人回复
-
山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元)
已经有3人回复
-
杭电、天津科技、青农和宁波工程学院如何选？
已经有6人回复
-
博士后创新人才支持计划公示
已经有11人回复
-
E12面上申请
已经有4人回复
-
车间生产，真空度很高，温度很高，但减压蒸馏速度很慢。
已经有31人回复
5