| 查看: 634 | 回复: 0 | ||
wusongz金虫 (小有名气)
|
[求助]
lammps得到如图g(r)，为什么峰会劈裂？
|
本人初学者，实在搞不懂。 units metal atom_style charge boundary p p p variable x index 1250 print "--------Temperature=$x K---------" read_data data_amorphous_vasp.CU pair_style comb pair_coeff * * ffield.comb Cu Si neighbor 0.5 bin neigh_modify every 4 delay 0 check yes group CU id <> 1 3200 group SI id <> 3201 4000 compute A CU property/atom q compute Q CU reduce ave c_A compute B SI property/atom q compute W SI reduce ave c_B set group CU charge 0.0 set group SI charge 0.0 timestep 0.0001 #dump 1 all cfg 10 *.cfg id type xs ys zs id q vx vy vz fx fy fz #dump_modify 1 element Cu Si dump 2 all atom 100 CuSi_$x.lammpstrj thermo 1000 velocity all create $x 825577 thermo_style custom step temp epair etotal vol press pxx pyy pzz c_Q c_W min_modify dmax 0.05 line quadratic minimize 1.0e-14 1.0e-20 20 100 fix 1 all nvt temp $x $x 0.01 drag 0.2 fix 2 all qeq/comb 1 3e-3 run 10000 unfix 1 min_modify dmax 0.05 line quadratic minimize 1.0e-14 1.0e-20 100000 100000 fix 3 all box/relax x 0.0 y 0.0 z 0.0 vmax 0.001 fix 4 all npt temp $x $x 0.01 x 0.0 0.0 0.01 y 0.0 0.0 0.01 z 0.0 0.0 0.01 drag 0.2 fix 5 all qeq/comb 1 3e-3 run 2200000 clear next x jump in.comb.2 |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有7人回复
-
无聊看看filecode
已经有10人回复
-
无聊看看时间戳打发时间
已经有8人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
27年博士招生信息
已经有7人回复
-
求推荐期刊，重谢
已经有4人回复
-
有机能跨方向申博吗
已经有3人回复
-
求助！！
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
薄膜在lammps条件下的谐振
已经有5人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- lammps中固定键长
已经有4人回复
- 【求助】lammps如何测量模拟以后材料的晶格参数？
已经有7人回复
- 【讨论】lammps中关于fix ave/spatial命令的输出问题
已经有4人回复
- 【求助】lammps中多元化合物的晶格创建问题
已经有14人回复
- 【讨论】lammps运行eam/fs出错
已经有38人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
5