<p>I'm trying to run a gas phase cluster calculation in CP2K. I want to create an outside force to force water molecules that are in the simulation within a radius of a central ion. I want to create a spherical cell of a given radius λ, fix the ion in the center and then run QMMD using the waters.</p><p>Does anyone know how to create a spherical cell or a better way to impose this clustering constraint λ?</p>