[交流] 【讨论】对撞模拟－－速度设定 已有2人参与

echo both

dimension       3
boundary        p p p

units metal

atom_style    atomic
neighbor      2.0 bin
neigh_modify    every 10 delay 0 check no

lattice         fcc 3.52
region box block 0.0 100.0 0.0 5.0 0.0 5.0
create_box      1 box
mass            1 58.6934

region     p3 block 0.0 5.0 0.0 5.0 0.0 5.0
region     p4 block 95.0 100.0  0.0 5.0 0.0 5.0

create_atoms    1 region p4
create_atoms    1 region p3

group p3 region p3
group p4 region p4

timestep        0.001
thermo          1000

fix 1 all nve

pair_style      lj/cut 2.5
pair_coeff      * * 1.0 1.0 2.5

velocity p3 set -20.0 0.0 0.0
velocity p4 set 20.0 0.0 0.0
fix 3 all   deform 1000 x final 45.0 55.0 units lattice
dump            1 all xyz 50 dump.xyz

run 1000