  The Neumann--Poincar\'e operator defined on a smooth surface has a sequence
of eigenvalues converging to zero, and the single layer potentials of the
corresponding eigenfunctions, called plasmons, decay to zero, i.e., are
localized on the surface, as the index of the sequence $j$ tends to infinity.
We investigate quantitatively the surface localization of the plasmons in three
dimensions. The results are threefold. We first prove that on smooth bounded
domains of general shape the sequence of plasmons converges to zero off the
boundary surface almost surely at the rate of $j^{-1/2}$. We then prove that if
the domain is strictly convex, then the convergence rate becomes $j^{-\infty}$,
namely, it is faster than $j^{-N}$ for any integer $N$. As a consequence, we
prove that cloaking by anomalous localized resonance does not occur on
three-dimensional strictly convex smooth domains. We then look into the surface
localization of the plasmons on the Clifford torus by numerical computations.
The Clifford torus is taken as an example of non-convex surfaces. The
computational results show that the torus exhibits the spectral property
completely different from strictly convex domains. In particular, they suggest
that there is a subsequence of plasmons on the torus which has much slower
decay than other entries of the sequence.
