Summary: In this paper we study strictly positive definite functions on the unit sphere of the \(m\)-dimensional Euclidean space. Such functions can be used for solving a scattered data interpolation problem on spheres. Since positive definite functions on the sphere were already characterized by Schoenberg some fifty years ago, the issue here is to determine what kind of positive definite functions are actually strictly positive definite. The study of this problem was initiated recently by \textit{Y. Xu} and \textit{E. W. Cheney} [Proc. Am. Math. Soc. 116, No. 4, 977-981 (1992; Zbl 0787.43005)], where certain sufficient conditions were derived. A new approach, which is based on a critical connection between this problem and that of multivariate polynomial interpolation on spheres, is presented here. The relevant interpolation problem is subsequently analyzed by three different complementary methods. The first is based on the de Boor-Ron general ``least solution for the multivariate polynomial interpolation problem''. The second, which is suitable only for \(m=2\), is based on the connection between bivariate harmonic polynomials and univariate analytic polynomials, and reduces the problem to the structure of the integer zeros of bounded univariate exponentials. Finally, the last method invokes the realization of harmonic polynomials as the polynomial kernel of the Laplacian, thereby exploiting some basic relations between homogeneous ideals and their polynomial kernels.