ERROR on proc 0: Neighbor list overflow, boost neigh_modify one or page

Hi everybody,
Could you please tell me why I am getting following error to the code attached below.
ERROR on proc 0: Neighbor list overflow, boost neigh_modify one or page
thanks in advance.
Samantaray

#code starts here

units si
boundary s s s
atom_style granular
atom_modify sort 0 0.0
newton off
communicate single vel yes

lattice fcc 0.0352 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 origin 0 0 0
region reg block 0 16.0 0 10.0 0 28.28427
create_box 1 reg

create_atoms 1 box

pair_style gran/hooke/history 1 0
pair_coeff * *

# Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 2.e+5
fix m2 all property/global poissonsRatio peratomtype 0.20
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.1
fix m4 all property/global coefficientFriction peratomtypepair 1 0.4
fix m6 all property/global characteristicVelocity scalar 2.
fix m7 all property/global cohesionEnergyDensity peratomtypepair 1 1.9e+4

timestep 0.0001
run 1

Forums:

LIGGGHTS® - User Forum

ckloss_ | Fri, 08/12/2011 - 22:51

try changing the lattice

try changing the lattice spacing and/or the cohesion parameter. I guess you system is either exploding or imploding due to a wrong setting in one of these params

Christoph

prasantud | Wed, 08/17/2011 - 08:06

how to calculate the no. of neighbors

Hi Cristoph,
I am getting too much neighbors for my following file.
I am using the fcc lattice, so it will have have co-ordination no 12.
But the output is giving me co-ordination no as 655.
It creates 655 atoms

###### code starts here ###############
units si
boundary s s s
atom_style granular
atom_modify map array sort 0 0.0
newton off
communicate single vel yes

lattice fcc 0.002
region reg block 0 1.6e-2 0 1.0e-2 0 0.006 units box
create_box 1 reg

create_atoms 1 box units box

pair_style gran/hooke/history 1 0
pair_coeff * *

fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.001

neighbor 0.0003 bin
neigh_modify delay 0

compute rad all coord/gran
variable b equal c_rad[1]
fix extra all print 1 "$b" #screen yes

timestep 0.000000000001

run 3

############### output

Created orthogonal box = (0 0 0) to (0.016 0.01 0.006)
Created 655 atoms

Total # of neighbors = 356712
Ave neighs/atom = 544.598
Neighbor list builds = 0

Step Temp E_pair E_mol TotEng Press Volume
0 0 0 0 0 3.4842622e+12 1.024488e-06
654

ckloss_ | Wed, 08/17/2011 - 18:25

you do not set the diameter

you do not set the diameter for the particles you insert, so each particle will touch each other as they will be inserted with default diameter of 1 m, you must use the set command

Christoph

Sounik | Wed, 12/05/2018 - 19:43

Try changing the lattice

This comment was really helpful. Thanks a lot Christoph

anandmds | Tue, 06/23/2015 - 08:07

Time step and neigh_modify command

What is the combination for timestep to be entered and the bin size and neighbor command ? because I'm getting "ERROR on proc 0: Neighbor list overflow, boost neigh_modify one or page " inccrease the time step from 0.0000001 to 0.0001 and so on, and I'm getting 0 atoms in the simulation if I reduce the bin size form 0.3 to 0.03. What combination should I try ?