24小时热门版块排行榜

返回列表

wusongz

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1077.3
红花: 1
帖子: 152
在线: 34.8小时
虫号: 839044
注册: 2009-09-02
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 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

回复此楼