Dear LAMMPS users
Can anyone explain why in the following input script the values of all angular velocities are zero?
It was working well before upgrading to a newer version of LAMMPS thus I belive there is a bug in the newer version.
Unfortunately, the older versions have another problem so I should work with this newer version.
units real
variable T equal 86
dimension 3
boundary p p p
lattice fcc 5.6
region box block -3 3 -3 3 -3 3
create_box 2 box
create_atoms 1 box
region sphere1 sphere 0 0 0 1
delete_atoms region sphere1
lattice fcc 3.62
create_atoms 2 region sphere1
group Ar type 1
group Cu type 2
mass 1 39.948
mass 2 63.546
pair_style lj/cut 10
pair_coeff 1 1 0.2381 3.405
pair_coeff 2 2 9.3525 2.3377
#_______________________________________________________________ minimizeation
minimize 1e-6 1e-9 1000 100000