Simulating the thermodynamics of charging in weak polyelectrolytes: the Debye-Huckel limit

The coil-globule transition in weak (annealed) polyelectrolytes involves a subtle balance of pH, charge strength, and solvation forces. In this work, we utilize a coarse-grained hybrid grand-canonical Monte Carlo and molecular dynamics approach to explore the swelling behavior of weak linear and star polyelectrolytes under different ionic screening conditions and pH. Importantly, we are able to quantify topology-dependent effects in charging which arise at the core of star polymers. Our results are suggestive of suppression of charging in star weak polyelectrolytes in comparison to linear weak polyelectrolytes. Furthermore, we characterize the coil-globule transition in linear and star weak polyelectrolyte through expanded ensemble density-of-states simulations which suggest a change from a first order to second order phase transition moving from linear to star polyelectrolytes. Lastly, we characterize the inhomogeneous charging across the weak star polyelectrolyte through observed shifts in Delta pK(o), and compare with experimental work. We discuss these results in relation to surfaces functionalized by weak polyelectrolyte brushes and weak polyelectrolyte-based drug delivery applications.