The problem of determination of globally optimal arrangements of N pairwise-interacting particles arises in a variety of biological, physical, and chemical applications. At the same time, the important related question of finding all, or many, local minima of the corresponding energy functions, and the study of structure of these minima, has received relatively little attention.A computational procedure is proposed to compute locally optimal and putative globally optimal arrangements of N particles constrained to a sphere. The procedure is able to handle a wide class of pairwise potentials, and can be generalized to other kinds of surfaces and interactions.As computational examples, locally and globally energy-minimizing arrangements of particles on the unit sphere, interacting via the Coulombic, logarithmic, and inverse square law potentials, are computed. We present new results for the logarithmic potential consisting of 45 new local minima for N≤65 and two new global minima (N=19,46), as well as results for the inverse square law potential which has not previously been studied. We provide comprehensive tables of all minima found, and exclude saddle points. The algorithm can perform computations exceeding N=100 with reasonable execution times. Program summaryProgram Title: EOPS 1.0 - Energy Optimizer for Particles on the SphereProgram Files doi:http://dx.doi.org/10.17632/cbn8jt2ffw.1Licensing provisions: GPLv3Programming language: MATLAB 2015b, C++98, MapleNature of problem: Computation of locally and globally optimal arrangements of N particles on the sphere for different pairwise potentials. This constitutes a constrained local optimization problem with 2N−3 degrees of freedom.Solution method: For N particles, the pairwise potential energy is minimized via steepest descent trajectory from a starting configuration generated from known putative (N−1)-particle optimal configurations.Restrictions: Spherical domain in R3 and pairwise potentials. The number of particles is limited by the computing power and memory of the machine.Unusual features: The programs are executed from MATLAB scripts which call C++ and Maple procedures which perform the bulk of the computations.