A NEW POTENTIAL FOR THE NUMERICAL SIMULATIONS OF ELECTROLYTE-SOLUTIONS ON A HYPERSPHERE

We propose a new way of performing numerical simulations of the restri cted primitive model of electrolytes-and related models-on a hypersphe re. In this new approach, the system is viewed as a single component f luid of charged bihard spheres constrained to move at the surface of a four dimensional sphere. A charged bihard sphere is defined as the ri gid association of two antipodal charged hard spheres of opposite sign s. These objects interact via a simple analytical potential obtained b y solving the Poisson-Laplace equation on the hypersphere. This new te chnique of simulation enables a precise determination of the chemical potential of the charged species in the canonical ensemble by a straig htforward application of Widom's insertion method. Comparisons with pr evious simulations demonstrate the efficiency and the reliability of t he method.