[lammps-users] About fix deposit

Dear all,

I have met a problem. I have used “fix deposit” command. However it seems impossible for me to add any atoms.

Although I have used the alternate command to create atoms, I still want to learn how to use the fix deposit command.

I will appreciate it if anyone can help me, the following is in.file:

ws 3d simulacion simple

units metal
dimension 3
boundary p p fs
newton on
atom_style molecular

read_data data.m-Bulk
region origbox block -0.056 44.164 -0.203 41.297 0.158 40.159 units box
group origbox region origbox

#add atom
lattice bcc 5
mass 1 69.72 #Ga
mass 2 14.01 #N
#pair parameter
pair_style tersoff
pair_coeff * * GaN.tersoff Ga N
neighbor 2 bin
neigh_modify delay 0
#region
region addbox block -0.056 44.164 -0.203 41.297 0.158 40.159 units box
group addbox region addbox
region fix block -0.056 44.164 -0.203 41.297 0.158 6.544 units box
group fix region fix
velocity fix set 0.0 0.0 0.0
fix 1 fix setforce 0.0 0.0 0.0
region mobile block -0.056 44.164 -0.203 41.297 7.00 40.159 units box
group mobile region mobile
velocity mobile create 1800.0 43454 dist gaussian
fix 2 mobile nvt 1800.0 1800.0 1 drag 0.0

compute mobile_temp all temp
compute_modify mobile_temp dynamic yes
#fix_modify 8 temp mobile_temp

region substrate block -0.056 44.164 -0.203 41.297 20.12 40.159 units box
group substrate region substrate
group substrate type 1 2

fix 3 substrate deposit 256 1 100 12345 region substrate vz -1.0 -1.0
fix 4 substrate deposit 512 2 100 12345 region substrate vz -1.0 -1.0

fix 5 substrate wall/reflect zhi

################################ injection1 ###########################################
thermo 100
dump 7 addbox xyz 100 dump.c2-T1
run 1000000