Dear Dr. Kohlmeyer,
Sorry, my original input has some remarks and notes, so I delete all of them before I copy that onto the e-mail.
And I carelessly wiped out one important code " pair_coeff * * ".
I know what you mean. But why does the error still show the same message that tells me boost MAXTOUCH.
And according to your statement, there is no any way to solve the problem I bump into?
The model of my simulation is what I want to investigate!
I need to put many spheres in a narrow cylinder.
It is not so important for me that many contacts between spheres happen before they reach equilibrium.
Here is my correct input that can be executed:
units cgs
atom_style sphere
boundary f f f
newton off
communicate single vel yes
lattice fcc 1.0
region boxin1 cylinder z 0 0 5 0 250 units lattice
create_box 1 boxin1
create_atoms 1 random 10000 5 boxin1 units box
group bal region boxin1
set region boxin1 diameter 1
neighbor 0.2 bin
neigh_modify delay 0
pair_style gran/hooke/history 200000.0 NULL 100.0 NULL 0.5 1
pair_coeff * *
timestep 0.0001
fix wigz bal wall/gran 2000000.0 NULL 10000.0 NULL 0.5 1 zcylinder 5
fix wigx bal wall/gran 2000000.0 NULL 10000.0 NULL 0.5 1 zplane 0 250
fix 1 bal nve/sphere
fix g bal gravity 1000.0 spherical 0.0 -180.0
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 100
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
dump id all custom 500 dump.10k_1 x y z
run 100000
Thanks,
Hank