Monte Carlo simulations of oppositely charged macroions in solution

The structure and phase behavior of oppositely charged macroions in solution have been studied with Monte Carlo simulations using the primitive model where the macroions and small ions are described as charged hard spheres. Size and charge symmetric, size asymmetric, and charge asymmetric macroions at different electrostatic coupling strengths are considered, and the properties of the solutions have been examined using cluster size distribution functions, structure factors, and radial distribution functions. At increasing electrostatic coupling, the macroions form clusters and eventually the system displays a phase instability, in analogy to that of simple electrolyte solutions. The relation to the similar cluster formation and phase instability occurring in solutions containing oppositely charged polymers is also discussed.