How colloidal particles interact with each other is one of the key issues that determines our ability to interpret experimental results for phase transitions in colloidal dispersions and our ability to apply colloid science to various industrial processes. The long-accepted theories for answering this question have been challenged by results from recent experiments. Herein we show from Monte-Carlo simulations that there is a short-range attractive force between identical macroions in electrolyte solutions containing divalent counterions. Complementing some recent and related results by others, we present strong evidence of attraction between a pair of spherical macroions in the presence of added salt ions for the conditions where the interacting macroion pair is not affected by any other macroions that may be in the solution. This attractive force follows from the internal-energy contribution of counterion mediation. Contrary to conventional expectations, for charged macroions in an electrolyte solution, the entropic force is repulsive at most solution conditions because of localization of small ions in the vicinity of macroions. Both Derjaguin-Landau-Verwey-Overbeek theory and Sogami-Ise theory fail to describe the attractive interactions found in our simulations; the former predicts only repulsive interaction and the latter predicts a long-range attraction that is too weak and occurs at macroion separations that are too large. Our simulations provide fundamental ''data'' toward an improved theory for the potential of mean force as required for optimum design of new materials including those containing nanoparticles.The potential of mean force between colloidal particles in electrolyte solutions plays a central role in describing the phase behavior and the kinetics of agglomeration in colloidal dispersions (1, 2). It is of fundamental importance for understanding the properties of inorganic materials (e.g., ceramics composed of nanoparticles), foods such as milk, and solutions of biomacromolecules including globular proteins (3-5). After decades of theoretical and experimental efforts, some aspects of colloidal interactions remain puzzling, in particular the issue of attractive electrostatic forces between like-charged colloidal particles in an electrolyte solution. Experimental evidence for such attraction has been indirect and largely complicated by boundary or polydispersity effects (6-8). Classical theories are not satisfactory because they are based on the mean-field approximation that neglects the effects of excluded volume and Coulombic correlations among small ions. Depending on simplifying assumptions, these theories have led to qualitatively different results. For example, theories based on the Derjaguin-Landau-Verwey-Overbeek (DLVO) approximation describe the electrostatic interaction between macroions of the same charge as screened repulsion (9-10, 12, 13), and others, represented by Sogami-Ise (SI) theory, predict a universal long-range attractive interaction (11,14). Likewise, the possibil...