Dear lammps users,
I am curious that is the pair style lj/cut/dipole/long a way to simulate a polarizable force field described by point dipoles. The Doc of this pair style shows the expressions, which are exactly how electronic polarization can be accounted for due to the change of the environment of each atom. If it does so, then why it doesn’t ask for the polarizability of atoms, and how the command fix nve/sphere update dipole works, the Docs of this fix command says that the orientation of the dipole is updating, what about the magnitude. I appreciate for any ideas and suggestions. Thanks very much.