A new method for reactive constant pH simulations
arXiv:2305.10521 · doi:10.1063/5.0166840
Abstract
We present a simulation method that allows us to calculate the titration curves for systems undergoing protonation/deprotonation reactions -- such as charged colloidal suspensions with acidic/basic surface groups, polyelectrolytes, polyampholytes, proteins, etc. The new approach allows us to simultaneously obtain titration curves both for systems in contact with salt and acid reservoir (semi-grand canonical ensemble) and for isolated suspensions (canonical ensemble). To treat the electrostatic interactions, we present a new method based on Ewald summation -- which accounts for the existence of both Bethe and Donnan potentials within the simulation cell. We show that that the Donnan potential affects dramatically the pH of suspension. Counter intuitively, we find that for concentrated suspensions of low ionic strength, the number of deprotonated groups can be 100\% larger in an isolated system, compared to a system connected to a reservoir by a semi-permeable membrane -- with both systems being at exactly the same pH.
References in corpus (14)
- Electrostatic correlations: from Plasma to Biology
- Charge Regulation with Fixed and Mobile Charges
- Charge-regulation effects in nanoparticle self-assembly
- Simulations of Coulomb systems with slab geometry using an efficient 3d Ewald summation method
- Accelerated simulation method for charge regulation effects
- Theory of Charge Regulation of Colloidal Particles in Electrolyte Solutions
- Critical behavior of charge-regulated macro-ions
- Note: A pairwise form of the Ewald sum for non-neutral systems
- Charge regulation radically modifies electrostatics in membrane stacks
- Interaction between charge-regulated metal nanoparticles in an electrolyte solution
- Reactive Monte Carlo Simulations for Charge Regulation of Colloidal Particles
- Like-charge attraction at the nanoscale: ground-state correlations and water destructuring
- Sedimentation of colloidal plate-sphere mixtures and inference of particle characteristics from stacking sequences
- Comment on "Simulations of ionization equilibria in weak polyelectrolyte solutions and gels" by J. Landsgesell, L. Nová, O. Rud, F. Uhlík, D. Sean, P. Hebbeker, C. Holm and P. Ko\v sovan, Soft Matter, 2019,15, 1155-1185