This article recapitulates the state of the art regarding simulations of ionization equilibria of weak polyelectrolyte solutions and gels.

In this work we study the titration behavior of weak polyelectrolytes by computer simulations. We analyze the local pH near the chains at various conditions and provide molecular-level insight which complements the recent experimental determination of this quantity. Next, we analyze the non-ideal titration behaviour of weak polyelectrolytes in solution, calculate the effective ionization constant and compare the simulation results with theoretical predictions. In contrast with the universal behaviour with respect to chain length, we find non-universality and deviations from theory with respect to polymer concentration and permittivity of the solvent. The latter we explain in terms of counterion condensation and ion correlation effects, which lead to reversal of the non-ideal titration behaviour at very low permittivities. We discuss the impact of these findings on the interpretation of experimental results.

We used computer simulations to explore the dissociative and conformational behaviour of branched weak polyelectrolytes with multivalent counterions.