Structural Properties of Aqueous Solutions at the (100) and (101) Goethite Surfaces by Molecular Dynamics Simulation

Trace metal concentrations in soils and sediments are often controlled by adsorption to iron oxides such as goethite in both natural and contaminated systems. Because of goethite's importance as an adsorbent, its interaction with aqueous solutions has been studied extensively. Nonetheless, despite the use of numerous analytical and computational tools, the properties of goethite-aqueous solution interfaces are not fully understood. In this research, we investigate the interaction of water and aqueous NaCl, MgCl₂, and BaCl₂ solutions ranging in concentration from 0.1 to 4 M, with two goethite surfaces, (100) and (101), using classical molecular dynamics simulation. In the past, the (100) surface has been studied the most because of its simplicity; however, goethite crystals in the environment exhibit other prominent surfaces like the (101) surface which may exhibit very different adsorption properties than the (100) surface. The (100) surface has three surface sites; one is an under-coordinated Fe³⁺ which interacts strongly with water affecting the interfacial water structure, another site remains deprotonated and forms a hydrogen bond to the only hydroxylated surface site. The (101) surface is terminated with five hydroxyl groups that do not interact as strongly with water as the under-coordinated Fe³⁺ ion and that form a more corrugated surface structure. As a result, the (101) goethite-solution interface exhibits less water structure, weaker electric double layer oscillations, and more inner-sphere ion adsorption, especially for Cl^- and Ba²⁺ ions. The fundamental differences in interfacial properties for these surfaces suggest that the adsorption properties of one goethite surface cannot be averaged to represent goethite interfaces present in soils and sediments.

Force Field for Water over Pt(111): Development, Assessment, and Comparison

Metal/water interfaces are key in many natural and industrial processes, such as corrosion, atmospheric, or environmental chemistry. Even today, the only practical approach to simulate large interfaces between a metal and water is to perform force-field simulations. In this work, we propose a novel force field, GAL17, to describe the interaction of water and a Pt(111) surface. GAL17 builds on three terms: (i) a standard Lennard-Jones potential for the bonding interaction between the surface and water, (ii) a Gaussian term to improve the surface corrugation, and (iii) two terms describing the angular dependence of the interaction energy. The 12 parameters of this force field are fitted against a set of 210 adsorption geometries of water on Pt(111). The performance of GAL17 is compared to several other approaches that have not been validated against extensive first-principles computations yet. Their respective accuracy is evaluated on an extended set of 802 adsorption geometries of H₂O on Pt(111), 52 geometries derived from icelike layers, and an MD simulation of an interface between a c(4 × 6) Pt(111) surface and a water layer of 14 Å thickness. The newly developed GAL17 force field provides a significant improvement over previously existing force fields for Pt(111)/H₂O interactions. Its well-balanced performance suggests that it is an ideal candidate to generate relevant geometries for the metal/water interface, paving the way to a representative sampling of the equilibrium distribution at the interface and to predict solvation free energies at the solid/liquid interface.

Molecular dynamics study on ultrathin liquid water film sheared between platinum solid walls: Liquid structure and energy and momentum transfer

Molecular dynamics simulation has been performed on a liquid film that is sheared in between solid surfaces. As a shear is given to the liquid film, a Couette-like flow is generated in the liquid and energy conversion occurs from the macroscopic flow to the thermal energy, which is discharged back to the solid walls. In such a way, momentum and thermal energy fluxes are present simultaneously. And all these thermal and fluid phenomena take place in highly nonequilibrium state where thermal energy is not distributed equally to each degree of freedom of molecular motion in the vicinities of the solid-liquid interface. In the present paper, platinum and water are employed as solid and liquid, respectively. First, the structure and orientation of water molecules in the vicinities of the solid surfaces are analyzed and how these structure and orientation are influenced by the shear is considered. Based on this result, momentum and thermal energy transfer in the vicinities of and at the solid-liquid interfaces are investigated in detail. Results are compared with those of our previous study, in which monatomic and diatomic molecules are employed as liquid.