ERROR on proc 0: Too many atom sorting bins

Vipul,

I suspect your problem is that your domain is enormous relative to the size of the particles and the neighbor size you've specified. Is your domain scaled properly?

Chris

With your region and bin settings, you will have around 200 million bins, which is just too much for LIGGGHTS to handle. In your case, increasing the binsize from 0.01 to 0.1 would surely resolve your issue.

Hey everyone ! I'm facing the similar error. In my case I have used only 5 particles and the binsize is also 0.1. The problem still persists

I'll attach the input file also

echo both
log ../DEM/log.liggghts
thermo_log ../DEM/post/thermo.txt

atom_style granular
atom_modify map array
communicate single vel yes

boundary p f f
newton off

units micro
processors 2 2 1

#read the restart file
read_restart ../DEM/post/restart/liggghts.restart

neighbor 0.1 bin
neigh_modify delay 0

# Smuts_Zirconia with Yukawa_Colloid
#------------------------------------------------------------------------------
#------------------------------------------------------------------------------
# Material properties required for granular pair styles_Zirconia
#------------------------------------------------------------------------------
#Variable definitions

variable YoungMod equal 88e5
variable PoissonsRat equal 0.33
variable CoeffResti equal 0.3
variable CoeffFrict equal 0.05

fix m1 all property/global youngsModulus peratomtype ${YoungMod}
fix m2 all property/global poissonsRatio peratomtype ${PoissonsRat}
fix m3 all property/global coefficientRestitution peratomtypepair 1 ${CoeffResti}
fix m4 all property/global coefficientFriction peratomtypepair 1 ${CoeffFrict}
fix m5 all property/global intermolecularDistance scalar 1e-4
fix m6 all property/global coefficientRollingFriction peratomtypepair 1 ${CoeffFrict}
fix m7 all property/global characteristicVelocity scalar 0.1
fix m8 all property/global cohesionEnergyDensity peratomtypepair 1 1.4e-2

#-------------------------------------------------------------------------------
# Pair Style definition_Yukawa_Colloid
#-------------------------------------------------------------------------------
# Variable definitions

variable ParticleDia equal 1.5
variable ParticleRad equal (${ParticleDia}/2)
variable ParticleDensity equal 6
variable allCutoffs equal 1.54 # 40nm + r1/2 +r2/2
variable NumParticle equal 5
variable KappaLength equal 1.0395e14 # Inverse of Debye Length
variable HamakerConst equal 9.24e-8
variable SigmaCC equal 2.78e-4
variable A_cc equal 3.6115e-4

#-------------------------------------------------------------------------------
#Physics definition
#-------------------------------------------------------------------------------
# Testing physics
#-------------------------------------------------------------------------------
#pair_style gran model hertz tangential history cohesion sjkr rolling_friction cdt

#pair_coeff * *
#-------------------------------------------------------------------------------
# Simulation Physics
#-------------------------------------------------------------------------------
#pair_style colloid 200.0
#colloid: pair_coeff i j hamaker sigma d1 d2 cutoff
#pair_coeff * * 144 1.0 0 0 100
#pair_coeff 1 1 39.5 278 150 150 200

# Yukawa/Colloid
pair_style yukawa/colloid ${KappaLength} ${allCutoffs}
pair_coeff * * ${A_cc} ${allCutoffs}

#-------------------------------------------------------------------------------
#timestep, gravity
#-------------------------------------------------------------------------------
variable timeCFD equal 0.025

timestep ${timeCFD}
fix gravi all gravity 0.00000981 vector 0.0 0.0 -1.0
#-------------------------------------------------------------------------------
#walls
#-------------------------------------------------------------------------------
#fix xwalls1 all wall/gran model hertz tangential history primitive type 1 xplane 0.0
#fix xwalls2 all wall/gran model hertz tangential history primitive type 1 xplane 100000
fix ywalls1 all wall/gran model hertz tangential history primitive type 1 yplane 0.0
fix ywalls2 all wall/gran model hertz tangential history primitive type 1 yplane 1000
fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0
fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 5000
#-------------------------------------------------------------------------------
#cfd coupling
#-------------------------------------------------------------------------------
fix cfd all couple/cfd couple_every 100 mpi
fix cfd2 all couple/cfd/force/implicit

#-------------------------------------------------------------------------------
# apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere
#fix ctg all check/timestep/gran 10 0.01 0.01
#-------------------------------------------------------------------------------

#center of mass
compute centerOfMass all com

#compute total dragforce
compute dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3]
#-------------------------------------------------------------------------------
#screen output
#-------------------------------------------------------------------------------
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol c_centerOfMass[3] c_dragtotal[1] c_dragtotal[2] c_dragtotal[3]
thermo 10
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

#insert the first particles so that dump is not empty
dump myDump all stl 1 post/dump_*.stl
#run 1
dump dmp all custom 10 ../DEM/post/dump*.liggghts_restart id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
undump myDump

run 0