24小时热门版块排行榜    

查看: 600  |  回复: 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
回复此楼

» 猜你喜欢

----Whathappenedtoyou?----Life.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 wusongz 的主题更新
信息提示
请填处理意见