Equilibrium Polyelectrolyte Bundles with Different Multivalent Counterion Concentrations

We present the results of molecular-dynamics simulations on the salt concentration dependence of the formation of polyelectrolyte bundles in thermodynamic equilibrium. Extending our results on salt-free systems we investigate here deficiency or excess of trivalent counterions in solution. Our results reveal that the trivalent counterion concentration significantly alters the bundle size and size distribution. The onset of bundle formation takes place at earlier Bjerrum length values with increasing trivalent counterion concentration. For the cases of 80%, 95%, and 100% charge compensation via trivalent counterions, the net charge of the bundles decreases with increasing size. We suggest that competition among two different mechanisms, counterion condensation and merger of bundles, leads to a nonmonotonic change in line-charge density with increasing Bjerrum length. The investigated case of having an abundance of trivalent counterions by 200% prohibits such a behavior. In this case, we find that the difference in effective line-charge density of different size bundles diminishes. In fact, the system displays an isoelectric point, where all bundles become charge neutral.