Molecular simulation of the magnetite-water interface

This paper reports molecular dynamics simulations of the magnetite (001)-water interface, both in pure water and in the presence of a 2.3 molal solution of NaClO4. The simulations are carried out using a potential model designed to allow the protonation states of the surface functional groups to evolve dynamically through the molecular dynamics trajectory. The primary structural quantities investigated are the populations of the surface functional groups, the distribution of electrolyte in the solution, and the surface hydrogen bonding relationships. The surface protonation states are dominated by extensive hydrolysis of interfacial water molecules, giving rise to a dipolar surface dominated by FeOH2+-OH2-OH− arrangements. Triply coordinated, more deeply buried, surface sites are inert, probably due to the relative lack of solvent in their vicinity. The electrolyte distribution is oscillatory, arranging preferentially in layers defined by the solvating water molecules. The presence of electrolyte has a negligible effect on the protonation states of the surface functional groups. Steady-state behavior is obtained for the protonation states of the surface functional groups and hydrogen-bonding network. Although the overall structure of the electrolyte distribution is fairly well established, the electrolyte distribution has not fully equilibrated, as evidenced by the asymmetry in the distribution from the top to the bottom of the slab.

[1]  L. Curtiss,et al.  Nonadditivity of ab initio pair potentials for molecular dynamics of multivalent transition metal ions in water , 1987 .

[2]  S. I. Yi,et al.  Surface structure of MBE-grown α-Fe 2 O 3 (0001) by intermediate-energy X-ray photoelectron diffraction , 1999 .

[3]  H. Jacobs,et al.  Ungewöhnliche H‐Brückenbindungen in Natriumhydroxidmonohydrat: Röntgen‐ und Neutronenbeugung an NaOH · H2O bzw. NaOD · D2O , 1991 .

[4]  E. Bylaska,et al.  Ab initio investigation of the structures of NaOH hydrates and their Na+ and OH– coordination polyhedra , 2003 .

[5]  F. Gygi,et al.  The solvation of Na+ in water: First-principles simulations , 2000 .

[6]  Michiel Sprik,et al.  A polarizable model for water using distributed charge sites , 1988 .

[7]  Venema,et al.  Intrinsic Proton Affinity of Reactive Surface Groups of Metal (Hydr)oxides: The Bond Valence Principle , 1996, Journal of colloid and interface science.

[8]  J. Rustad.  Molecular Models of Surface Relaxation, Hydroxylation, and Surface Charging at Oxide-Water Interfaces , 2001 .

[9]  The hematite (Alpha-Fe_2O_3)(0001) surface: Evidence for domains of distinct chemistry , 1998, cond-mat/9807202.

[10]  D. Dixon,et al.  Intrinsic acidity of aluminum, chromium (III) and iron (III) μ3-hydroxo functional groups from ab initio electronic structure calculations , 2000 .

[11]  W. Rudziński,et al.  A Combined Temperature−Calorimetric Study of Ion Adsorption at the Hematite−Electrolyte Interface: II. Models of a Heterogeneous Oxide Surface , 2000 .

[12]  J. W. Halley,et al.  Molecular dynamics simulation of iron(III) and its hydrolysis products in aqueous solution , 1995 .

[13]  D. A. Palmer,et al.  Magnetite surface charge studies to 290°C from in situ pH titrations , 2000 .

[14]  S. Sjöberg,et al.  Modeling proton binding at the goethite (α-FeOOH)–water interface , 2001 .

[15]  Bedzyk,et al.  Electrical Double-Layer Structure at the Rutile-Water Interface as Observed in Situ with Small-Period X-Ray Standing Waves. , 2000, Journal of colloid and interface science.

[16]  Shuichi Nosé,et al.  MOMENTUM CONSERVATION LAW IN THE CAR-PARRINELLO METHOD , 1999 .

[17]  V. Lobo.  Electrolyte solutions : literature data on thermodynamic and transport properties , 1981 .

[18]  W. Bleam.  On Modeling Proton Affinity at the Oxide/Water Interface , 1993 .

[19]  J. W. Halley,et al.  Ewald methods for polarizable surfaces with application to hydroxylation and hydrogen bonding on the (012) and (001) surfaces of α-Fe2O3 , 1997, cond-mat/9704070.

[20]  U. Diebold,et al.  Surface reconstruction of Fe3O4(001) , 2000 .

[21]  K. Rosso,et al.  Structures and energies of AlOOH and FeOOH polymorphs from plane wave pseudopotential calculations , 2001 .

[22]  Dennis R. Salahub,et al.  Optimization of Gaussian-type basis sets for local spin density functional calculations. Part I. Boron through neon, optimization technique and validation , 1992 .

[23]  W. Rudziński,et al.  A combined temperature-calorimetric study of ion adsorption at the hematite-electrolyte interface: I. Model of a homogeneous oxide surface , 2000 .

[24]  A. Navrotsky,et al.  ENERGETICS OF STABLE AND METASTABLE LOW-TEMPERATURE IRON OXIDES AND OXYHYDROXIDES , 1998 .

[25]  T. Healy,et al.  Adsorption of hydrolyzable metal ions at the oxide—water interface. III. A thermodynamic model of adsorption , 1972 .

[26]  J. Rustad,et al.  Molecular statics calculations of proton binding to goethite surfaces: A new approach to estimation of stability constants for multisite surface complexation models , 1996 .

[27]  Y. Kim,et al.  Selective growth and characterization of pure, epitaxial α-Fe2O3(0001) and Fe3O4(001) films by plasma-assisted molecular beam epitaxy , 1997 .

[28]  J. Rustad,et al.  MOLECULAR STATICS CALCULATIONS OF PROTON BINDING TO GOETHITE SURFACES : THERMODYNAMIC MODELING OF THE SURFACE CHARGING AND PROTONATION OF GOETHITE IN AQUEOUS SOLUTION , 1998 .

[29]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[30]  D. Dixon,et al.  Trivalent Ion Hydrolysis Reactions: A Linear Free-Energy Relationship Based on Density Functional Electronic Structure Calculations , 1999 .

[31]  S. Ziemniak.  Magnetite Surface Charge Studies to 290 degrees C from In-Situ pH Titrations , 1999 .

[32]  E. Wasserman,et al.  The Magnetite (001) Surface: Insights from Molecular Dynamics Calculations , 2000 .

[33]  J. Rustad,et al.  Interaction of water with the (1×1) and (2×1) surfaces of α-Fe2O3(012) , 1998 .

[34]  Frank H. Stillinger,et al.  Polarization model for water and its ionic dissociation products , 1978 .

[35]  C. Dellago,et al.  Autoionization in Liquid Water , 2001, Science.

[36]  James R. Rustad,et al.  A polarizable, dissociating molecular dynamics model for liquid water , 1993 .

[37]  W. Stumm.  Chemistry of the solid-water interface , 1992 .

[38]  J. W. Halley,et al.  Molecular dynamics studies of complexing in binary molten salts with polarizable anions: MAX4 , 1988 .