Abstract The calculations of long‐range interaction energy are often based on multipole expansion. The truncated multipole expansion and interaction energy calculated with it are noninvariant with respect to an arbitrary choice of local coordinate systems. In this paper we show that truncated multipole expansion of form Σ C k R − k is “numerically” independent on a choice of local coordinate systems, if convergence conditions are satisfied.