Some recent trends in computer simulations of aqueous double layers

Orientational dynamics of the water layer adjacent to Au surface accelerated by polarization effect

The orientation and rearrangement of water on a gold electrode significantly influences its physicochemical heterogeneous performance. Despite numerous experimental and theoretical studies aimed at uncovering the structural characteristics of interfacial water, the orientational behavior resulting from electrode-induced rearrangements remains a subject of ongoing debate. Here, we employed molecular dynamics simulations to investigate the adaptive structure and dynamics properties of interfacial water on Au(111) and Au(100) surfaces by considering a polarizable model for Au atoms in comparison with the non-polarizable model. Compared to the nonpolarizable systems, the polarization effect can enhance the interaction between water molecules and the gold surface. Unexpectedly, the rotational dynamics directly associated with the orientational behavior of water adjacent to the gold surface is accelerated, thereby reducing the hydrogen bond lifetime. The underlying mechanism for this anomalous phenomenon originates from the polarization effect, which induces the attraction of the positive hydrogen atoms to the surface by the negative image charge. This leads to a change in orientation that disrupts the hydrogen bonds in the first water layer and subsequently accelerates reorientation dynamics of water molecules adjacent to the gold surface. These results shed light on the intricate interplay between polarization effects and water molecule dynamics on metal surfaces, establishing the foundation for the rational regulation of the orientation of interfacial water.

Modeling interfacial electrochemistry: concepts and tools

This paper presents a grand canonical formalism and provides tools to investigate electrochemical effects at interfaces.

Polarizable Embedding with a Transferable H2O Potential Function II: Application to (H2O)n Clusters and Liquid Water

The incorporation of polarization in multiscale quantum-mechanics/molecular-mechanics (QM/MM) simulations is important for a variety of applications, for example, charge-transfer reactions. A recently developed formalism based on a density functional theory description of the QM region and a potential energy function for H₂O molecules that includes quadrupole as well as dipole polarizability of the MM region is used to simulate liquid water and water clusters. Analysis of the energy, atomic forces, MM polarization, and structure is presented. A quantitative assessment of the QM/MM-MM/MM interaction energy differences of all possible QM/MM configurations of (H₂O)_n clusters shows that the interquartile range of the distributions of the QM/MM binding energies is never more than 20 meV/molecule higher or lower than the binding energies produced with either of the single-model results. Comparing these interaction energy differences with the QM/MM induction differences show that they are not systematically caused by the induced MM moments of our polarizable embedding scheme. Optimized hexamer geometries as well as the liquid water structure are shown to be improved in comparison with results obtained using point-charge based embedding models neglecting polarization.

Off-center charge model revisited: Electrical double layer with multivalent cations

The off-center charge model of ions is a relatively simple model for introducing asymmetry in Coulomb interaction while retaining the simplicity and convenience of the spherical hard core geometry. A Monte Carlo simulation analysis of the planar electric double layer formed by this ionic model for 1+:1- valence systems [S. Lamperski et al., Langmuir 33, 11554-11560 (2017)] is extended to include solutions of multivalent (2+, 3+) hard spherical cations and single valence (1-) hard spherical anions near a uniformly charged, planar electrode. The solvent is modelled as a uniform dielectric continuum with a dielectric constant equal to that of the pure solvent, viz., the primitive model. Results are reported for the ion density, the cation charge profile, and the electrostatic potential profile at 1 mol/dm³ salt concentration. Additionally, the double layer potential drop, that is, the electrode potential, and the integral and the differential capacitances are computed as functions of the electrode surface charge density. The latter two quantities show an expected asymmetry as long as the cation valence is not too great and the charge of the off-center ion cannot approach too close to the electrode surface. It is unusual that the integral and differential capacitances are negative for high valence cations and a negatively charged electrode when the off-center charge is large and can be very near the surface of the electrode. The corresponding electrode potential versus surface charge density curve becomes non-monotonic and shows a change of slope, and thus the resultant integral and differential capacitances can become negative. This nonphysical result is the result of an incipient singularity when a large positive charge is too near a negatively charged electrode. Overall, the off-center charge model suggests a useful recipe to model electrical asymmetry within the broader context of the primitive model provided that the off-center charge is not too near the surface of the electrode.

Ion-Image Interactions and Phase Transition at Electrolyte-Metal Interfaces

The arrangement of ions near a metallic electrode is crucial to energy storage in electrical double-layer capacitors. Classic Poisson-Boltzmann theory predicts that the charge stored in the double layer is a continuous function of applied voltage. However, recent experiments and simulations strongly suggest the presence of a voltage-induced first-order phase transition in the electrical double layer, leading to a hysteretic response: the capacitance-voltage relation is dependent on whether the voltage is increasing or decreasing. By developing a simple analytical model, we show that ion-image interaction could explain this phase transition. Moreover, our model shows that the presence of phase transition depends on the bulk energy of the ionic liquid. Our results justify mixing ionic liquids with solvents as a way to achieve large capacitance and avoid hysteresis.

Neural network molecular dynamics simulations of solid–liquid interfaces: water at low-index copper surfaces

Molecular dynamics simulation of the water–copper interface have been carried out using high-dimensional neural network potential based on density functional theory.

Atomistic Modeling of Corrosion Events at the Interface between a Metal and Its Environment

Atomistic simulation is a powerful tool for probing the structure and properties of materials and the nature of chemical reactions. Corrosion is a complex process that involves chemical reactions occurring at the interface between a material and its environment and is, therefore, highly suited to study by atomistic modeling techniques. In this paper, the complex nature of corrosion processes and mechanisms is briefly reviewed. Various atomistic methods for exploring corrosion mechanisms are then described, and recent applications in the literature surveyed. Several instances of the application of atomistic modeling to corrosion science are then reviewed in detail, including studies of the metal-water interface, the reaction of water on electrified metallic interfaces, the dissolution of metal atoms from metallic surfaces, and the role of competitive adsorption in controlling the chemical nature and structure of a metallic surface. Some perspectives are then given concerning the future of atomistic modeling in the field of corrosion science.

Water at an electrochemical interface--a simulation study

The results of molecular dynamics simulations of the properties of water in an aqueous ionic solution close to an interface with a model metallic electrode are described. In the simulations the electrode behaves as an ideally polarizable hydrophilic metal, supporting image-charge interactions with charged species, and it is maintained at a constant electrical potential with respect to the solution so that the model is a textbook representation of an electrochemical interface through which no current is passing. We show how water is strongly attracted to and ordered at the electrode surface. This ordering is different to the structure that might be imagined from continuum models of electrode interfaces. Further, this ordering significantly affects the probability of ions reaching the surface. We describe the concomitant motion and configurations of the water and ions as functions of the electrode potential, and we analyze the length scales over which ionic atmospheres fluctuate. The statistics of these fluctuations depend upon surface structure and ionic strength. The fluctuations are large--sufficiently so that the mean ionic atmosphere is a poor descriptor of the aqueous environment near a metal surface. The importance of this finding for a description of electrochemical reactions is examined by calculating, directly from the simulation, Marcus free-energy profiles for transfer of charge between the electrode and a redox species in the solution and comparing the results with the predictions of continuum theories. Significant departures from the electrochemical textbook descriptions of the phenomenon are found and their physical origins are characterized from the atomistic perspective of the simulations.

Insights from theory and simulation on the electrical double layer

Despite the fact that our conceptual understanding of the electrical double layer has advanced during the past few decades, the interpretation of experimental and applied work is still largely based on the venerable Poisson-Boltzmann theory of Gouy, Chapman and Stern. This is understandable since this theory is simple and analytic. However, it is not very accurate because the atomic/molecular nature of the ions/solvent and their correlations are ignored. Simulation and some theoretical studies by ourselves and others that have advanced our understanding are discussed. These studies show that the GCS theory predicts a narrow double layer with monotonic profiles. This is not correct. The double layer is wider, and there can be substantial layering that would be even more pronounced if explicit solvent molecules are considered. For many years, experimental studies of the double layer have been directed to the use of electrochemistry as an analytical tool. This is acceptable for analytic chemistry studies. However, the understanding of electrochemical reactions that typically occur at the electrode surface, where simulation and theory indicate that the GCS theory can have substantial errors, requires modern approaches. New, fundamental experimental studies that would lead to deeper insights using more novel systems would be desirable. Further, biophysics is an interesting field. Recent studies of the selectivity of ion channels and of the adsorption of ions in a binding sites of a protein have shown that the linearized GCS theory has substantial errors.

Liquid water confined in carbon nanochannels at high temperatures

Structure, hydrogen bonding, electrostatics, dielectric, and dynamical properties of liquid water confined in flat graphene nanochannels are investigated by molecular dynamics simulations. A wide range of temperatures (between 20 and 360 degrees C) have been considered. Molecular structure suffers substantial changes when the system is heated, with a significant loss of structure and hydrogen bonding. In such case, the interface between adsorbed and bulk-like water has a marked tendency to disappear, and the two preferential orientations of water nearby the graphite layers at room temperature are essentially merging above the boiling point. The general trend for the static dielectric constant is its reduction at high temperature states, as compared to ambient conditions. Similarly, residence times of water molecules in adsorbed and bulk-like regions are significantly influenced by temperature, as well. Finally, we observed relevant changes in water diffusion and spectroscopy along the range of temperatures analyzed.

Surface charge density and electrokinetic potential of highly charged minerals: experiments and Monte Carlo simulations on calcium silicate hydrate

In this paper, we are concerned with the charging and electrokinetic behavior of colloidal particles exhibiting a high surface charge in the alkaline pH range. For such particles, a theoretical approach has been developed in the framework of the primitive model. The charging and electrokinetic behavior of the particles are determined by the use of a Monte Carlo simulation in a grand canonical ensemble and compared with those obtained through the mean field theory. One of the most common colloidal particles has been chosen to test our theoretical approach. That is calcium silicate hydrate (C-S-H) which is the main component of hydrated cement and is known for being responsible for cement cohesion partly due to its unusually high surface charge density. Various experimental techniques have been used to determine its surface charge and electrokinetic potential. The experimental and simulated results are in excellent agreement over a wide range of electrostatic coupling, from a weakly charged surface in contact with a reservoir containing monovalent ions to a highly charged one in contact with a reservoir with divalent ions. The electrophoretic measurements show a charge reversal of the C-S-H particles at high pH and/or high calcium concentration in excellent agreement with simulation predictions. Finally, both simulation and experimental results clearly demonstrate that the mean field theory fails not only quantitatively but also qualitatively to describe a C-S-H dispersion under realistic conditions.

Effect of the surface charge discretization on electric double layers: A Monte Carlo simulation study

The structure of the electric double layer in contact with discrete and continuously charged planar surfaces is studied within the framework of the primitive model through Monte Carlo simulations. Three different discretization models are considered together with the case of uniform distribution. The effect of discreteness is analyzed in terms of charge density profiles. For point surface groups, a complete equivalence with the situation of uniformly distributed charge is found if profiles are exclusively analyzed as a function of the distance to the charged surface. However, some differences are observed moving parallel to the surface. Significant discrepancies with approaches that do not account for discreteness are reported if charge sites of finite size placed on the surface are considered.

Testing one component plasma models on colloidal overcharging phenomena

In this paper, the mechanisms of overcharging of a colloidal macroion in the presence of multivalent counterions are investigated by means of Monte Carlo simulations. This computational technique appears as a powerful tool for probing the validity of semianalytical models developed for this issue. In particular, the simulations performed are compared with the predictions of two different models based on the one component plasma (OCP) theory. Therein, the multivalent ionic atmosphere confined at the macroion surface is approximated by a two-dimensional Wigner crystal. These kinds of models are largely used in the literature since (in some cases) they present quite simple equations to describe the electric double layer (EDL) of macroions with different geometries in the presence of much smaller (but still multivalent) ions. In this sense, charge inversion phenomena of membranes, polyelectrolytes, DNA molecules, etc., are straightforwardly predicted in terms of these expressions. Unfortunately, comparisons between these predictions and experimental results are scarce, mostly due to the difficulty to reproduce the experimental conditions in the laboratory. Accordingly, the goal of the present paper is to simulate EDLs under real conditions (in which overcharging phenomena are expected to happen) and use the results obtained in this way for comparing with those obtained from OCP models.

Molecular dynamics simulation of electro-osmotic flows in rough wall nanochannels

We performed equilibrium and nonequilibrium molecular dynamics simulation to study electro-osmotic flows inside charged nanochannels with different types of surface roughness. We modeled surface roughness as a sequence of two-dimensional subnanoscale grooves and ridges (step function-type roughness) along the flow direction. The amplitude, spatial period, and symmetry of surface roughness were varied. The amplitude of surface roughness was on the order of the Debye length. The walls have uniform negative charges at the interface with fluids. We included only positive ions (counterions) for simplicity of computation. For the smooth wall, we compared our molecular dynamics simulation results to the well-known Poisson-Boltzmann theory. The density profiles of water molecules showed "layering" near the wall. For the rough walls, the density profiles measured from the wall are similar to those for the smooth wall except near where the steps are located. Because of the layering of water molecules and the finite size effect of ions and the walls, the ionic distribution departs from the Boltzmann distribution. To further understand the structure of water molecules and ions, we computed the polarization density. Near the wall, its z component dominates the other components, indicating the preferred orientation ("ordering") of water molecules. Especially, inside the groove for the rough walls, its maximum is 10% higher (stronger ordering) than for the smooth wall. The dielectric constant, computed with a Clausius-Mosotti-type equation, confirmed the ordering near the wall and the enhanced ordering inside the groove. The residence time and the diffusion coefficient, computed using the velocity autocorrelation function, showed that the diffusion of water and ions along the direction normal to the wall is significantly reduced near the wall and further decreases inside the groove. Along the flow direction, the diffusion of water and ions inside the groove is significantly lowered while it is similar to the bulk value elsewhere. We performed nonequilibrium molecular dynamics simulation to compute electro-osmotic velocities and flow rates. The velocity profiles correspond to those for overlapped electric double layers. For the rough walls, velocity inside the groove is close to zero, meaning that the channel height is effectively reduced. The flow rate was found to decrease as the period of surface roughness decreases or the amplitude of surface roughness increases. We defined the zeta potential as the electrostatic potential at the location of a slip plane. We computed the electrostatic potential with the ionic distribution and the dielectric constant both from our molecular dynamics simulation. We estimated the slip plane from the velocity profile. The zeta potential showed the same trend as the flow rate: it decreases with an increasing amplitude and a decreasing period of surface roughness.

Simulation of electric double layers undergoing charge inversion: mixtures of mono- and multivalent ions

In this paper, the electric double layer (EDL) of a charged plane in the presence of mixtures of 1:1 and 3:1 electrolytes has been investigated through Monte Carlo (MC) simulations using a nonrestrictive primitive model of EDL. In particular, the charge inversion in colloids (attributable to an accumulation of counterions on the surface) can be better understood by means of the simulations performed here. Moreover, two mechanisms proposed for charge inversion are probed: The formation of a strongly correlated layer (SCL) of multivalent counterions and excluded volume effects (to which we will also refer as ion size correlations). Our results are in agreement with the behavior found experimentally for some model colloids with increasing the concentration of monovalent salt in the presence of trivalent ions, which clearly supports the relevance of ion size correlations. In contrast, certain disagreement with predictions of SCL theories is reported.

Relaxation of the electrical double layer after an electron transfer approached by Brownian dynamics simulation

In this paper, the dynamical properties of the electrochemical double layer following an electron transfer are investigated by using Brownian dynamics simulations. This work is motivated by recent developments in ultrafast cyclic voltammetry which allow nanosecond time scales to be reached. A simple model of an electrochemical cell is developed by considering a 1:1 supporting electrolyte between two parallel walls carrying opposite surface charges, representing the electrodes; the solution also contains two neutral solutes representing the electroactive species. Equilibrium Brownian dynamics simulations of this system are performed. To mimic electron transfer processes at the electrode, the charge of the electroactive species are suddenly changed, and the subsequent relaxation of the surrounding ionic atmosphere are followed, using nonequilibrium Brownian dynamics. The electrostatic potential created in the center of the electroactive species by other ions is found to have an exponential decay which allows the evaluation of a characteristic relaxation time. The influence of the surface charge and of the electrolyte concentration on this time is discussed, for several conditions that mirror the ones of classical electrochemical experiments. The computed relaxation time of the double layer in aqueous solutions is found in the range 0.1 to 0.4 ns for electrolyte concentrations between 0.1 and 1 mol L(-1) and surface charges between 0.032 and 0.128 C m(-2).

Simulation of electric double layers with multivalent counterions: Ion size effect

In this paper, the structure of the electric double layer in the presence of (mostly) multivalent counterions is investigated through Monte Carlo simulations. Unlike previous similar studies addressing this matter, the difference of this study lies in the use of realistic hydrated ion sizes. Additionally, two different methods for calculating energies in the Metropolis algorithm are applied. The obtained results show that the conclusions of preceding papers must be revised. In particular, our simulations suggest the existence of certain ion layering effects at high surface charge densities, which are not accounted for by integral equation theories in the case of divalent counterions. These layering effects could justify why the overcharging phenomena due to ion size correlations are hardly observable in real colloids with divalent counterions. The existence of charge inversion due to ion size correlations (and without requiring specific counterion adsorption) is probed for trivalent counterions. Moreover, the hypernetted-chain/mean-spherical-approximation is tested under conditions not studied yet.