3 years ago

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

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
Andrey G. Kalinichev, Randall T. Cygan, Maxime Pouvreau, Jeffery A. Greathouse
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 (Cygan, R. T.; Liang, J.-J.; Kalinichev, A. G. J. Phys. Chem. B 2004, 108, 1255–1266). 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. 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.

Publisher URL: http://dx.doi.org/10.1021/acs.jpcc.7b05362

DOI: 10.1021/acs.jpcc.7b05362

You might also like
Discover & Discuss Important Research

Keeping up-to-date with research can feel impossible, with papers being published faster than you'll ever be able to read them. That's where Researcher comes in: we're simplifying discovery and making important discussions happen. With over 19,000 sources, including peer-reviewed journals, preprints, blogs, universities, podcasts and Live events across 10 research areas, you'll never miss what's important to you. It's like social media, but better. Oh, and we should mention - it's free.

  • Download from Google Play
  • Download from App Store
  • Download from AppInChina

Researcher displays publicly available abstracts and doesn’t host any full article content. If the content is open access, we will direct clicks from the abstracts to the publisher website and display the PDF copy on our platform. Clicks to view the full text will be directed to the publisher website, where only users with subscriptions or access through their institution are able to view the full article.