[lammps-users] polarizable core/shell model

I would like to use a recently published force field for Cadmium Sulfide:

Wright, K. and Gale, J.D., PHYSICAL REVIEW B, v.70, p. 035211 (2004)

For the sulphur atom it requires a core/shell model of polarization. Based on what I know so far, I think it should be possible to construct the core/shell model in LAMMPS using harmonically tethered dimers with fixed charges. My main concerns are numerical stability of the time integration and long-range electrostatic calculation. Has anyone tried to do a simulation like this in LAMMPS?

Thanks in advance,