In recent years there has been considerable study of the potential mechanisms underlying the electromechanical response of ionic-polymer-metal composites. The most recent models have been based on the response of the ion-containing clusters that are formed when the material is synthesized. Most of these efforts have employed assumptions of uniform ion distribution within spherical cluster shapes. This work investigates the impact of dispensing with these assumptions in order to better understand the parameters that impact cluster shape, size, and ion transport potential. A computational micromechanics model applying Monte Carlo methodology is employed to predict the equilibrium state of a single cluster of a solvated ionomeric polymer. For a constant solvated state, the model tracks the position of individual ions within a given cluster in response to ion-ion interaction, mechanical stiffness of the pendant chain, cluster surface energy, and external electric-field loading. Results suggest that cluster surface effects play a significant role in the equilibrium cluster state, including ion distribution; pendant chain stiffness also plays a role in ion distribution but to a lesser extent. Moreover, ion pairing is rarely complete even in cation-rich clusters; this in turn supports the supposition of the formation of anode and cathode boundary layers. (c) 2005 American Institute of Physics.