skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Structure of Hydrated Gibbsite and Brucite Edge Surfaces: DFT Results and Further Development of the ClayFF Classical Force Field with Metal–O–H Angle Bending Terms

Abstract

Molecular scale understanding of the structure and properties of aqueous interfaces with clays, metal (oxy-) hydroxides, layered double hydroxides, and other inorganic phases is strongly affected by significant degrees of structural and compositional disorder of the interfaces. ClayFF was originally developed as a robust and flexible force field for classical molecular simulations of such systems. However, despite its success, multiple limitations have also become evident with its use. One of the most important limitations is the difficulty to accurately model the edges of finite size nanoparticles or pores rather than infinitely layered periodic structures. Here we propose a systematic approach to solve this problem by developing specific metal–O–H (M–O–H) bending terms for ClayFF, Ebend = k (θ – θ0)2 to better describe the structure and dynamics of singly protonated hydroxyl groups at mineral surfaces, particularly edge surfaces. On the basis of a series of DFT calculations, the optimal values of the Al–O–H and Mg–O–H parameters for Al and Mg in octahedral coordination are determined to be θ0,AlOH = θ0,MgOH = 110°, kAlOH = 15 kcal mol–1 rad–2 and kMgOH = 6 kcal mol–1 rad–2. Molecular dynamics simulations were performed for fully hydrated models of the basal and edge surfaces ofmore » gibbsite, Al(OH)3, and brucite, Mg(OH)2, at the DFT level of theory and at the classical level, using ClayFF with and without the M–O–H term. The addition of the new bending term leads to a much more accurate representation of the orientation of O–H groups at the basal and edge surfaces. Finally, the previously observed unrealistic desorption of OH2 groups from the particle edges within the original ClayFF model is also strongly constrained by the new modification.« less

Authors:
 [1];  [2];  [2]; ORCiD logo [1]
  1. Laboratoire SUBATECH (UMR 6457), Institut Mines-Télécom Atlantique, 44307 Nantes, France
  2. Geochemistry Department, Sandia National Laboratories, P.O Box 5800, MS 0754, Albuquerque, New Mexico 87185-0754, United States
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1369576
Alternate Identifier(s):
OSTI ID: 1367350
Report Number(s):
SAND-2017-6270J
Journal ID: ISSN 1932-7447
Grant/Contract Number:  
NA0003525; AC04-94AL85000
Resource Type:
Published Article
Journal Name:
Journal of Physical Chemistry. C
Additional Journal Information:
Journal Name: Journal of Physical Chemistry. C Journal Volume: 121 Journal Issue: 27; Journal ID: ISSN 1932-7447
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY

Citation Formats

Pouvreau, Maxime, Greathouse, Jeffery A., Cygan, Randall T., and Kalinichev, Andrey G. Structure of Hydrated Gibbsite and Brucite Edge Surfaces: DFT Results and Further Development of the ClayFF Classical Force Field with Metal–O–H Angle Bending Terms. United States: N. p., 2017. Web. https://doi.org/10.1021/acs.jpcc.7b05362.
Pouvreau, Maxime, Greathouse, Jeffery A., Cygan, Randall T., & Kalinichev, Andrey G. Structure of Hydrated Gibbsite and Brucite Edge Surfaces: DFT Results and Further Development of the ClayFF Classical Force Field with Metal–O–H Angle Bending Terms. United States. https://doi.org/10.1021/acs.jpcc.7b05362
Pouvreau, Maxime, Greathouse, Jeffery A., Cygan, Randall T., and Kalinichev, Andrey G. Wed . "Structure of Hydrated Gibbsite and Brucite Edge Surfaces: DFT Results and Further Development of the ClayFF Classical Force Field with Metal–O–H Angle Bending Terms". United States. https://doi.org/10.1021/acs.jpcc.7b05362.
@article{osti_1369576,
title = {Structure of Hydrated Gibbsite and Brucite Edge Surfaces: DFT Results and Further Development of the ClayFF Classical Force Field with Metal–O–H Angle Bending Terms},
author = {Pouvreau, Maxime and Greathouse, Jeffery A. and Cygan, Randall T. and Kalinichev, Andrey G.},
abstractNote = {Molecular scale understanding of the structure and properties of aqueous interfaces with clays, metal (oxy-) hydroxides, layered double hydroxides, and other inorganic phases is strongly affected by significant degrees of structural and compositional disorder of the interfaces. ClayFF was originally developed as a robust and flexible force field for classical molecular simulations of such systems. However, despite its success, multiple limitations have also become evident with its use. One of the most important limitations is the difficulty to accurately model the edges of finite size nanoparticles or pores rather than infinitely layered periodic structures. Here we propose a systematic approach to solve this problem by developing specific metal–O–H (M–O–H) bending terms for ClayFF, Ebend = k (θ – θ0)2 to better describe the structure and dynamics of singly protonated hydroxyl groups at mineral surfaces, particularly edge surfaces. On the basis of a series of DFT calculations, the optimal values of the Al–O–H and Mg–O–H parameters for Al and Mg in octahedral coordination are determined to be θ0,AlOH = θ0,MgOH = 110°, kAlOH = 15 kcal mol–1 rad–2 and kMgOH = 6 kcal mol–1 rad–2. Molecular dynamics simulations were performed for fully hydrated models of the basal and edge surfaces of gibbsite, Al(OH)3, and brucite, Mg(OH)2, at the DFT level of theory and at the classical level, using ClayFF with and without the M–O–H term. The addition of the new bending term leads to a much more accurate representation of the orientation of O–H groups at the basal and edge surfaces. Finally, the previously observed unrealistic desorption of OH2 groups from the particle edges within the original ClayFF model is also strongly constrained by the new modification.},
doi = {10.1021/acs.jpcc.7b05362},
journal = {Journal of Physical Chemistry. C},
number = 27,
volume = 121,
place = {United States},
year = {2017},
month = {6}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.1021/acs.jpcc.7b05362

Citation Metrics:
Cited by: 6 works
Citation information provided by
Web of Science

Save / Share:

Works referencing / citing this record:

Alumina: discriminative analysis using 3D correlation of solid-state NMR parameters
journal, January 2019

  • Chandran, C. Vinod; Kirschhock, Christine E. A.; Radhakrishnan, Sambhu
  • Chemical Society Reviews, Vol. 48, Issue 1
  • DOI: 10.1039/c8cs00321a

Distinguishing between bulk and edge hydroxyl vibrational properties of 2 : 1 phyllosilicates via deuteration
journal, January 2019

  • Harvey, Jacob A.; Johnston, Cliff T.; Criscenti, Louise J.
  • Chemical Communications, Vol. 55, Issue 24
  • DOI: 10.1039/c9cc00164f

Understanding methane/carbon dioxide partitioning in clay nano- and meso-pores with constant reservoir composition molecular dynamics modeling
journal, January 2019

  • Loganathan, Narasimhan; Bowers, Geoffrey M.; Ngouana Wakou, Brice F.
  • Physical Chemistry Chemical Physics, Vol. 21, Issue 13
  • DOI: 10.1039/c9cp00851a

Interfacial reactions of Cu( ii ) adsorption and hydrolysis driven by nano-scale confinement
journal, January 2020

  • Knight, Andrew W.; Ilani-Kashkouli, Poorandokht; Harvey, Jacob A.
  • Environmental Science: Nano, Vol. 7, Issue 1
  • DOI: 10.1039/c9en00855a

How do Nucleotides Adsorb Onto Clays?
journal, November 2018

  • Pedreira-Segade, Ulysse; Hao, Jihua; Razafitianamaharavo, Angelina
  • Life, Vol. 8, Issue 4
  • DOI: 10.3390/life8040059