ERROR: Illegal delete_atoms command (…/delete_atoms.cpp:68)
I m working on LAMMPS. I m using it first time. my task is to create a crystal with different percentage of vacancy. I have written a input script and try to run the simulation. It is showing this error. Can you help me to write the delete command to create the random vancancy in the periodic crystal which contains 10976 atoms.
#input file for 0.1% vacancy
#initialization
processors 1 1 1
units metal
dimension 3
boundary p p p
atom_style atomic
Read data file
read_data CuStructure.data
#interatomic Potential and simulation settings
timestep 0.001
pair_style eam/alloy
pair_coeff * * Cu01.eam.alloy Cu
neighbor 2.0 bin
neigh_modify every 10
Create vacancies
region Traget block INF INF INF INF INF INF units box
delete_atoms random fraction 0.001 no all target 10946
#run a simualtion
minimize 1.0e-4 1.0e-6 100 1000
velocity all create 700.0 87827 mom yes rot yes dist gaussian
fix 1 all nvt temp 700.0 700.0 0.1
dump 1 all custom 1000 Cu_0.1pct_trajectory.lammpstrj id type x y z
dump_modify 1 sort id
compute msd all msd
fix msd_out all ave/time 1000 1 1000 c_msd[4] file msd_0.1pct.txt mode vector
thermo_style custom step temp pe etotal press vol
thermo 1000
run 100000
unfix 1
undump 1
reset_timestep 0
This is the input script I have tried to write by own. can you help me to wrie the delete command properly. I want to create a vacancy of 1% of total atoms. Also point out the mistakes and guide me for my future work.
Thank You