I eventually updated my LAMMPS from 12April version to 11Jul. And my scripts for hybrid granular system stopped working. It stops after:
atom_style hybrid sphere dipole
boundary f f f
communicate single vel yes
region reg block 0.00000E+00 0.10000E+02 0.00000E+00 0.10000E+02 0.00000E+00 0.20000E+02 units box
create_box 2 reg
neighbor 20.000 bin
neigh_modify delay 0
lattice sc 1.0000
pair_style hybrid/overlay gran/hertz/history 0.389503E+10 0.000000E+00 0.000000E+00 NULL 0.5 1 dipole/cut 2.5000
pair_coeff 1 1 gran/hertz/history
pair_coeff 1 1 dipole/cut 0.0 0.0 2.5000 2.5000
pair_coeff 1 2 gran/hertz/history
pair_coeff 1 2 dipole/cut 0.0 0.0 2.5000 2.5000
pair_coeff 2 2 gran/hertz/history
pair_coeff 2 2 dipole/cut 0.0 0.0 2.5000 2.5000
with an error:
ERROR: Illegal pair_style command (…/pair_gran_hertz_history.cpp:287)
Would be really grateful for the help.
please take a look at the note from June 6th here: http://lammps.sandia.gov/bug.html
Stan Moore (Sandia) has done some nice work to clean up and fix bugs in some code for both short- and long-range dipole/dipole and charge/dipole interactions. Two new pair_styles were also added: pair_style lj/cut/dipole/long and lj/long/coul/long. Long-range dipole interactions can now be computed by kspace_style ewald/disp. All of the options are described in the pertinent pair_style and kspace_style doc pages.
BACKWARD COMPATIBILITY: The names of two of the previous dipole pair styles were changed to be more compatible with the LAMMPS naming conventions. Specifically, dipole/cut was changed to lj/cut/dipole/cut, and dipole/sf was changed to lj/cut/dipole/sf.
You are using the pair style name that was changed between the two releases.
18.07.2013, 09:56, "Oksana Golovnya" <[email protected]>: