A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of p H and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups

Competition effects in cation binding to humic acid: Conditional affinity spectra for fixed total metal concentration conditions

Information on the Pb and Cd binding to a purified Aldrich Humic Acid (HA) is obtained from the influence of different fixed total metal concentrations on the acid- base titrations of this ligand. NICA (Non-Ideal Competitive Adsorption) isotherm has been used for a global quantitative description of the binding, which has then been interpreted by plotting the Conditional Affinity Spectra of the H+ binding at fixed total metal concentrations (CAScTM). This new physicochemical tool, here introduced, allows the interpretation of binding results in terms of distributions of proton binding energies. A large increase in the acidity of the phenolic sites as the total metal concentration increases, especially in presence of Pb, is revealed from the shift of the CAScTM towards lower affinities. The variance of the CAScTM distribution, which can be used as a direct measure of the heterogeneity, also shows a significant dependence on the total metal concentration. A discussion of the factors that influence the heterogeneity of the HA under the conditions of each experiment is provided, so that the smoothed pattern exhibited by the titration curves can be justified