What comes to mind very quickly is the corresponding magnetostatics problem with a cylinder of uniform magnetization ## \vec{M} ##. In that case, ## \nabla \times \vec{M} =\vec{J}_m ## results in surface currents per unit length of ## \vec{K}_m=\vec{M} \times \hat{n} ## on the surface of the cylinder. ## \\ ## In the simplest case of uniform ## \vec{ P} ##, if you take ## \nabla \times \vec{P} ##, you will get places where ## \nabla \times \vec{P } ## diverges at parts of the surface/air interface. In general, if ## \vec{P} ## is uniform, the derivative ## \nabla \times \vec{P} ## vanishes, but this derivative can diverge when ## \vec{P} ## undergoes a discontinuity such as at the surface/air interface.