I was wondering whether a granular simulation with multisphere particles can be done on LAMMPS. Can fix rigid be used to serve the purpose? I am confused about the template to be used for the fix rigid. If I use a molecule template, it requires molecular atom_style for the fix rigid to work. But I need atom style to be spherical for granular simulation. Any suggestions?