Grand-canonical Monte Carlo method for Donnan equilibria

Titration in Canonical and Grand-Canonical Ensembles

We discuss problems associated with the notion of pH in heterogeneous systems. For homogeneous systems, standardization protocols lead to a well-defined quantity, which, although different from Sørensen's original idea of pH, is well reproducible and has become accepted as the measure of the "hydrogen potential". On the other hand, for heterogeneous systems, pH defined in terms of the chemical part of the electrochemical activity is thermodynamically inconsistent and runs afoul of the Gibbs-Guggenheim principle that forbids splitting of the electrochemical potential into separate chemical and electrostatic parts, since only the sum of two has any thermodynamic meaning. The problem is particularly relevant for modern simulation methods which involve charge regulation of proteins, polyelectrolytes, nanoparticles, colloidal suspensions, and so forth. In this paper, we show that titration isotherms calculated using semigrand canonical simulations can be very different from the ones obtained using canonical reactive Monte Carlo simulations.

A new method for reactive constant pH simulations

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, and proteins. 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 the Donnan potential dramatically affects the pH of a suspension. Counterintuitively, we find that in suspensions with a large volume fraction of nanoparticles and 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-both systems being at exactly the same pH.

Widom insertion method in simulations with Ewald summation

We discuss the application of the Widom insertion method for calculation of the chemical potential of individual ions in computer simulations with Ewald summation. Two approaches are considered. In the first approach, an individual ion is inserted into a periodically replicated overall charge neutral system representing an electrolyte solution. In the second approach, an inserted ion is also periodically replicated, leading to the violation of the overall charge neutrality. This requires the introduction of an additional neutralizing background. We find that the second approach leads to a much better agreement with the results of grand canonical Monte Carlo simulation for the total chemical potential of a neutral ionic cluster.

Reactive Monte Carlo simulations for charge regulation of colloidal particles

We use a reactive Monte Carlo simulation method and the primitive model of electrolyte to study acid-base equilibrium that controls charge regulation in colloidal systems. The simulations are performed in a semi-grand canonical ensemble in which colloidal suspension is in contact with a reservoir of salt and strong acid. The interior of colloidal particles is modeled as a low dielectric medium, different from the surrounding water. The effective colloidal charge is calculated for different numbers of surface acidic groups, pH, salt concentrations, and types of electrolyte. In the case of potassium chloride, the titration curves are compared with the experimental measurements obtained using potentiometric titration. A good agreement is found between simulations and experiments. In the case of lithium chloride, the specific ionic adsorption is taken into account through the partial dehydration of lithium ion.

Improving the efficiency of Monte Carlo simulations of ions using expanded grand canonical ensembles

While ionic liquids have promising applications as industrial solvents, predicting their fluid phase properties and coexistence remains a challenge. Grand canonical Monte Carlo simulation is an effective method for such predictions, but equilibration is hampered by the apparent requirement to insert and delete neutral sets of ions simultaneously in order to maintain charge neutrality. For relatively high densities and low temperatures, previously developed methods have been shown to be essential in improving equilibration by gradual insertion and deletion of these neutral sets of ions. We introduce an expanded ensemble approach which may be used in conjunction with these existing methods to further improve efficiency. Individual ions are inserted or deleted in one Monte Carlo trial rather than simultaneous insertion/deletion of neutral sets. We show how charge neutrality is maintained and show rigorous quantitative agreement between the conventional and the proposed expanded ensemble approaches, but with up to an order of magnitude increase in efficiency at high densities. The expanded ensemble approach is also more straightforward to implement than simultaneous insertion/deletion of neutral sets, and its implementation is demonstrated within open source software.