[lammps-users] Error

Dear all
I want to set the indentor on a grain boundary, and I have witten the following. However, when I run it, I face the error " Domain too large for neighbor bins". Could you say what is wrong with my code?

dimension 2

units metal

boundary s s p

atom_style atomic

neighbor 0.3 bin

neigh_modify delay 5

create geometry

lattice hex 1

region box block -20 20 0 20 -0.25 0.25

create_box 5 box

mass 1 1.0

mass 2 1.0

mass 3 1.0

mass 4 1.0

mass 5 1.0

LJ potentials

pair_style lj/cut 2.5

pair_coeff * * 1.0 1.0 2.5

define groups

region 1 block INF INF INF 1.25 INF INF

create_atoms 1 region 1

group lower region 1

region 2 block INF 0 1.25 INF INF INF

create_atoms 2 region 2

group left region 2

lattice hex 1.5

region 3 block 0.25 INF 1 INF INF INF

create_atoms 3 region 3

group rigth region 3

group mobile union rigth left

set group left type 2

set group rigth type 3

set group lower type 1

compute new mobile temp

velocity mobile create 0.2 482748 temp new

fix 1 all nve

fix 2 lower setforce 0.0 0.0 0.0

fix 3 all temp/rescale 100 0.1 0.1 0.01 1.0

#indentor

fix 4 all indent 1000.0 sphere 0 23 0 5.0 vel 0.0 -0.02 0.0

fix 5 all enforce2d

timestep 0.003

thermo 200

thermo_modify temp new

dump 1 all atom 500 dump.work1

run 5000

"Domain too large" typically occurs b/c you've blown
atoms out of the simulation box, meaning you have an
error. Are your thermodynamic energies reasonable? Have
you visualized what's happening in your simulation?

Steve