Spiers Memorial Lecture: Towards understanding of iontronic systems: electroosmotic flow of monovalent and divalent electrolyte through charged cylindrical nanopores.

In many practical applications, ions are the primary charge carrier and must move through either semipermeable membranes or through pores, which mimic ion channels in biological systems. In analogy to electronic devices, the "iontronic" ones use electric fields to induce the charge motion. However, unlike the electrons that move through a conductor, motion of ions is usually associated with simultaneous solvent flow. A study of electroosmotic flow through narrow pores is an outstanding challenge that lies at the interface of non-equilibrium statistical mechanics and fluid dynamics. In this paper, we will review recent works that use dissipative particle dynamics simulations to tackle this difficult problem. We will also present a classical density functional theory (DFT) based on the hypernetted-chain approximation (HNC), which allows us to calculate the velocity of electroosmotic flows inside nanopores containing 1 : 1 or 2 : 1 electrolyte solution. The theoretical results will be compared with simulations. In simulations, the electrostatic interactions are treated using the recently introduced pseudo-1D Ewald summation method. The zeta potentials calculated from the location of the shear plane of a pure solvent are found to agree reasonably well with the Smoluchowski equation. However, the quantitative structure of the fluid velocity profiles deviates significantly from the predictions of the Smoluchowski equation in the case of charged pores with 2 : 1 electrolyte. For low to moderate surface charge densities, the DFT allows us to accurately calculate the electrostatic potential profiles and the zeta potentials inside the nanopores. For pores with 1 : 1 electrolyte, the agreement between theory and simulation is particularly good for large ions, for which steric effects dominate over the ionic electrostatic correlations. The electroosmotic flow is found to depend very strongly on the ionic radii. In the case of pores containing 2 : 1 electrolyte, we observe a reentrant transition in which the electroosmotic flow first reverses and then returns to normal as the surface change density of the pore is increased.

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. Kosovan, Soft Matter, 2019, 15, 1155-1185.

In a recent review Landsgesell et al., Soft Matter 2019, 15, 1155 stated that pH - pKa is a "universal parameter" for titrating systems. We show that this is not the case. This broken symmetry has important implications for constant pH (cpH) simulations. In particular, we show that for concentrated suspensions the error resulting from the use of cpH algorithm described by Landsgesell et al. is very significant, even for suspension containing 1 : 1 electrolyte. We show how to modify the cpH algorithm to account for the grand-canonical nature of the cpH simulations and for the charge neutrality requirement.

Charge fluctuations in charge-regulated systems: dependence on statistical ensemble.

We investigate charge regulation of nanoparticles in concentrated suspensions, focusing on the effect of different statistical ensembles. We find that the choice of ensemble does not affect the mean charge of nanoparticles, but significantly alters the magnitude of its fluctuation. Specifically, we compared the behaviors of colloidal charge fluctuations in the semi-grand canonical and canonical ensembles and identified significant differences between the two. The choice of ensemble-whether the system is isolated or is in contact with a reservoir of acid and salt-will, therefore, affect the Kirkwood-Shumaker fluctuation-induced force inside concentrated suspensions. Our results emphasize the importance of selecting an appropriate ensemble that accurately reflects the experimental conditions when studying fluctuation-induced forces between polyelectrolytes, proteins, and colloidal particles in concentrated suspensions.

Canonical titration simulations.

We present a Monte Carlo approach for performing titration simulations in the canonical ensemble. The standard constant pH (cpH) simulation methods are intrinsically grand canonical, allowing us to study the protonation state of molecules only as a function of pH in the reservoir. Due to the Donnan potential between a system and an (implicit) reservoir of a semi-grand canonical simulation, the pH of the reservoir can be significantly different from that of an isolated system, for an identical protonation state. The new titration method avoids this difficulty by using the canonical reactive Monte Carlo algorithm to calculate the protonation state of macromolecules as a function of the total number of protons present inside the simulation cell. The pH of an equilibrated system is then calculated using a new surface insertion Widom algorithm, which bypasses the difficulties associated with the bulk Widom particle insertion for intermediate and high pH values. To properly treat the long range Coulomb force, we use the Ewald summation method, showing the importance of the Bethe potential for calculating the pH of canonical systems.

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.

A Monte Carlo Method for Calculating Lynden-Bell Equilibrium in Self-Gravitating Systems.

We present a Monte Carlo approach that allows us to easily implement Lynden-Bell (LB) entropy maximization for an arbitrary initial particle distribution. The direct maximization of LB entropy for an arbitrary initial distribution requires an infinite number of Lagrange multipliers to account for the Casimir invariants. This has restricted studies of Lynden-Bell's violent relaxation theory to only a very small class of initial conditions of a very simple waterbag form, for which the entropy maximization can be performed numerically. In the present approach, an arbitrary initial distribution is discretized into density levels which are then evolved using an efficient Monte Carlo algorithm towards the final equilibrium state. A comparison is also made between the LB equilibrium and explicit Molecular Dynamics simulations. We find that for most initial distributions, relaxation is incomplete and the system is not able to reach the state of maximum LB entropy. In particular, we see that the tail of the stationary particle distribution is very different from the one predicted by the theory of violent relaxation, with a hard cutoff instead of an algebraic decay predicted by LB's theory.

Theory of Charge Regulation of Colloidal Particles in Electrolyte Solutions.

We present a theory that enables us to (i) calculate the effective surface charge of colloidal particles and (ii) efficiently obtain titration curves for different salt concentrations. The theory accounts for the shift of pH of solution due to the presence of 1:1 electrolyte. It also accounts self-consistently for the electrostatic potential produced by the deprotonated surface groups. To examine the accuracy of the theory, we have performed extensive reactive Monte Carlo simulations, which show excellent agreement between theory and simulations without any adjustable parameters.

Reversal of Electroosmotic Flow in Charged Nanopores with Multivalent Electrolyte.

We study the reversal of electroosmotic flow in charged cylindrical nanopores containing multivalent electrolyte. Dissipative particle dynamics is used to simulate the hydrodynamics of the electroosmotic flow. The electrostatic interactions are treated using 3D Ewald summation, corrected for a pseudo-one-dimensional geometry of the pore. We observe that, for sufficiently large surface charge density, condensation of multivalent counterions leads to the reversal of the pore's surface charge. This results in the reversal of electroosmotic flow. Our simulations show that the Smoluchowski equation is able to quantitatively account for the electroosmotic flow through the nanopore, if the shear plane is shifted from the position of the Stern contact surface.

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.

Lattice-gas model of a charge regulated planar surface.

In this work, we consider a lattice-gas model of charge regulation with electrostatic interactions within the Debye-Hückel level of approximation. In addition to long-range electrostatic interactions, the model incorporates the nearest-neighbor interactions for representing non-electrostatic forces between adsorbed ions. The Frumkin-Fowler-Guggenheim isotherm obtained from the mean-field analysis accurately reproduces the simulation data points.

Osmotic stress and pore nucleation in charged biological nanoshells and capsids.

A model system is proposed to investigate the chemical equilibrium and mechanical stability of biological spherical-like nanoshells in contact with an aqueous solution with added dissociated electrolyte of a given concentration. The ionic chemical equilibrium across the permeable shell is investigated in the framework of an accurate Density Functional Theory (DFT) that incorporates electrostatic and hardcore correlations beyond the traditional mean-field (e.g., Poisson-Boltzmann) limit. The accuracy of the theory is tested by a direct comparison with Monte Carlo (MC) simulations. A simple analytical expression is then deduced which clearly highlights the entropic, electrostatic, and self-energy contributions to the osmotic stress over the shell in terms of the calculated ionic profiles. By invoking a continuum mean-field elastic approach to account for the shell surface stress upon osmotic stretching, the mechanical equilibrium properties of the shell under a wide variety of ionic strengths and surface charges are investigated. The model is further coupled to a continuum mechanical approach similar in structure to a Classical Nucleation Theory (CNT) to address the question of mechanical stability of the shells against a pore nucleation. This allows us to construct a phase diagram which delimits the mechanical stability of capsids for different ionic strengths and shell surface charges.

Charge regulation of colloidal particles in aqueous solutions.

We study the charge regulation of colloidal particles inside aqueous electrolyte solutions. To stabilize a colloidal suspension against precipitation, colloidal particles are synthesized with either acidic or basic groups on their surface. On contact with water, these surface groups undergo proton transfer reactions, resulting in colloidal surface charge. The charge is determined by the condition of local chemical equilibrium between hydronium ions inside the solution and at the colloidal surface. We use a model of Baxter sticky spheres to explicitly calculate the equilibrium dissociation constants and to construct a theory which is able to quantitatively predict the effective charge of colloidal particles with either acidic or basic surface groups. The predictions of the theory for the model are found to be in excellent agreement with the results of Monte Carlo simulations. This theory is further extended to treat colloidal particles with a mixture of both acidic and basic surface groups.

Simulations of electrolyte between charged metal surfaces.

We present a new method for simulating ungrounded charged metal slabs inside an electrolyte solution. The ions are free to move between the interior and exterior regions of the slab-electrolyte system. This leads to polarization of both sides of each slab, with a distinct surface charge induced on each surface. Our simulation method is based on the exact solution of the Poisson equation using periodic Green functions. To efficiently perform the calculations, we decouple the electrostatic energy due to surface polarization from that of purely Coulomb interaction between the ions. This allows us to combine a fast 3D Ewald summation technique with an equally fast calculation of polarization. As a demonstration of the method, we calculate ionic density profiles inside an electrolyte solution and explore charge neutrality violation in between charged metal slabs.

Like-Charge Attraction between Metal Nanoparticles in a 1∶1 Electrolyte Solution.

We calculate the force between two spherical metal nanoparticles of charge Q_{1} and Q_{2} in a dilute 1∶1 electrolyte solution. Numerically solving the nonlinear Poisson-Boltzmann equation, we find that metal nanoparticles with the same sign of charge can attract one another. This is fundamentally different from what is found for like-charged, nonpolarizable, colloidal particles, the two-body interaction potential for which is always repulsive inside a dilute 1∶1 electrolyte. Furthermore, the existence of like-charge attraction between spherical metal nanoparticles is even more surprising in view of the result that such attraction is impossible between parallel metal slabs, showing the fundamental importance of curvature. To overcome a slow convergence of the numerical solution of the full nonlinear Poisson-Boltzmann equation, we developed a modified Derjaguin approximation which allows us to accurately and rapidly calculate the interaction potential between two metal nanoparticles or between a metal nanoparticle and a phospholipid membrane.

Charge Regulation of Colloidal Particles: Theory and Simulations.

To explore charge regulation (CR) in physicochemical and biophysical systems, we present a model of colloidal particles with sticky adsorption sites which account for the formation of covalent bonds between the hydronium ions and the surface functional groups. Using this model and Monte Carlo simulations, we find that the standard Ninham and Parsegian (NP) theory of CR leads to results which deviate significantly from computer simulations. The problem with the NP approach is traced back to the use of a bulk equilibrium constant to account for surface chemical reactions. To resolve this difficulty we present a new theory of CR. The fundamental ingredient of the new approach is the sticky length, which is nontrivially related to the bulk equilibrium constant. The theory is found to be in excellent agreement with computer simulations, without any adjustable parameters. As an application of the theory we calculate the effective charge of colloidal particles containing carboxyl groups, as a function of pH and salt concentration.

Isothermal adsorption of polyampholytes on charged nanopatterned surfaces.

We investigate the adsorption of neutral polyampholytes on charged nanopatterned surfaces. The surfaces have charged domains but are overall neutral. To perform efficient simulations, we use an approach which combines the explicit form of the interaction potential between the polyampholyte monomers and the surface with a 3d Ewald summation method. We observe that the amount of adsorption and the structure of the adsorbed polyampholytes depend strongly on the surface pattern, the relative size of the surface domains, and the charge distribution along the polyampholyte backbone.

Nonequilibrium Statistical Mechanics of Two-Dimensional Vortices.

It has been observed empirically that two-dimensional vortices tend to cluster, forming a giant vortex. To account for this observation, Onsager introduced the concept of negative absolute temperature in equilibrium statistical mechanics. In this Letter, we show that in the thermodynamic limit a system of interacting vortices does not relax to the thermodynamic equilibrium but becomes trapped in a nonequilibrium stationary state. We show that the vortex distribution in this nonequilibrium stationary state has a characteristic core-halo structure, which can be predicted a priori. All the theoretical results are compared with explicit molecular dynamics simulations.