Hi, Steve:
Thanks for your suggestion. Now I’m testing the scheme. However, the removing atoms doesn’t seemt to work like I expected.
The test system is a slab of diamond atoms. All the atoms starts with a 10,000K initial temperature plus a -3A/ps initial speed in center of mass. Then the system is going NVE so that the atoms will move towards the left end of the box. I was expecting to see the atoms get removed when they hit the left end. However, the program stops when this happens.
I wish the program continues to run. What shall I do? Thanks,
Here’s the input of my lammps run:
variable x index 6
variable y index 6
variable z index 6
units metal
atom_style atomic
boundary p p f
lattice diamond 3.5657
region sbox block 0 $x 0 $y 0 $z
region mbox block 0 $x 0 $y 5 $z
create_box 1 sbox
create_atoms 1 region mbox basis 1 1
neighbor 0.5 bin
neigh_modify delay 1
pair_style rebo 0.0 0 0
pair_coeff * * potentials/ch.airebonew C
mass * 12.01
velocity all set NULL NULL -3 loop geom
velocity all create 10000.0 376847 sum yes loop geom
compute T_ALL all temp/com
fix LNVT all nve
timestep 0.0001
thermo 1
thermo_modify format float %.10f
dump Dump2 all custom 10000 xyzFxyz-*.bin tag x y z fx fy fz vx vy vz
run 50000