Hi All!
I am trying to simulate equal no. of positive and negative charged particles confined in a 3D box. I am using the SET type/ratio command to make half of the total no. of particles to either charge. But each time, I am getting the following error:
ERROR: Numeric index is out of bounds (…/set.cpp:634)
Last command: set type all type/ratio 2 0.5 58912
Below is my script. Kindly assist me here:
units LJ
dimension 3
atom_style charge
boundary p p p
lattice sc ${rhostar}
region box block -5 5 -5 5 -5 5
create_box 2 box
create_atoms 1 box
set type all type/ratio 2 0.5 58912
mass 1 1.0
mass 2 1.0
set type 1 charge 1.0
set type 2 charge -1.0
velocity all create 0.40 87287
I shall be highly obliged for this kind of act.
Thanking you