Ab initio molecular dynamics study of the interlayer and micropore structure of aqueous montmorillonite clays

James L. Suter, Lara Kabalan, Mahmoud Khader, Peter V. Coveney

Research output: Contribution to journalArticle

15 Citations (Scopus)

Abstract

Ab initio molecular dynamics simulations have been performed to gain an understanding of the interfacial microscopic structure and reactivity of fully hydrated clay edges. The models studied include both micropore and interlayer water. We identify acidic sites through dissociation mechanisms; the resulting ions can be stabilized by both micropore and interlayer water. We find clay edges possess a complex amphoteric behavior, which depends on the face under consideration and the location of isomorphic substitution. For the neutral (110) surface, we do not observe any dissociation on the timescale accessible. The edge terminating hydroxyl groups participate in a hydrogen bonded network of water molecules that spans the interlayer between periodic images of the clay framework. With isomorphic substitutions in the tetrahedral layer of the (110) clay edge, we find the adjacent exposed apical oxygen behaves as a Brönsted base and abstracts a proton from a nearby water molecule, which in turn removes a proton from an AlOH<inf>2</inf> group. With isomorphic substitutions in the octahedral layer of the (110) clay edge the adjacent exposed apical oxygen atom does not abstract a proton from the water molecules, but increases the number of hydrogen bonded water molecules (from one to two). Acid treated clays are likely to have both sites protonated. The (010) surface does not have the same interfacial hydrogen bonding structure; it is much less stable and we observe dissociation of half the terminal SiOH groups (≡Si-O-H→≡SiO<sup>-</sup>+H<sup>+</sup>) in our models. The resulting anions are stabilized by solvation from both micropore and interlayer water molecules. This suggests that, when fully hydrated, the (010) surface can act as a Brönsted acid, even at neutral pH.

Original languageEnglish
Pages (from-to)17-29
Number of pages13
JournalGeochimica et Cosmochimica Acta
Volume169
DOIs
Publication statusPublished - 15 Nov 2015

Fingerprint

Bentonite
montmorillonite
Molecular dynamics
clay
Water
Molecules
Protons
substitution
Substitution reactions
water
hydrogen
Hydrogen
Oxygen
oxygen
Acids
acid
Solvation
Hydroxyl Radical
Anions
anion

ASJC Scopus subject areas

  • Geochemistry and Petrology

Cite this

Ab initio molecular dynamics study of the interlayer and micropore structure of aqueous montmorillonite clays. / Suter, James L.; Kabalan, Lara; Khader, Mahmoud; Coveney, Peter V.

In: Geochimica et Cosmochimica Acta, Vol. 169, 15.11.2015, p. 17-29.

Research output: Contribution to journalArticle

@article{21580f2cd5064ed5b0d1d0058d15b0a3,
title = "Ab initio molecular dynamics study of the interlayer and micropore structure of aqueous montmorillonite clays",
abstract = "Ab initio molecular dynamics simulations have been performed to gain an understanding of the interfacial microscopic structure and reactivity of fully hydrated clay edges. The models studied include both micropore and interlayer water. We identify acidic sites through dissociation mechanisms; the resulting ions can be stabilized by both micropore and interlayer water. We find clay edges possess a complex amphoteric behavior, which depends on the face under consideration and the location of isomorphic substitution. For the neutral (110) surface, we do not observe any dissociation on the timescale accessible. The edge terminating hydroxyl groups participate in a hydrogen bonded network of water molecules that spans the interlayer between periodic images of the clay framework. With isomorphic substitutions in the tetrahedral layer of the (110) clay edge, we find the adjacent exposed apical oxygen behaves as a Br{\"o}nsted base and abstracts a proton from a nearby water molecule, which in turn removes a proton from an AlOH2 group. With isomorphic substitutions in the octahedral layer of the (110) clay edge the adjacent exposed apical oxygen atom does not abstract a proton from the water molecules, but increases the number of hydrogen bonded water molecules (from one to two). Acid treated clays are likely to have both sites protonated. The (010) surface does not have the same interfacial hydrogen bonding structure; it is much less stable and we observe dissociation of half the terminal SiOH groups (≡Si-O-H→≡SiO-+H+) in our models. The resulting anions are stabilized by solvation from both micropore and interlayer water molecules. This suggests that, when fully hydrated, the (010) surface can act as a Br{\"o}nsted acid, even at neutral pH.",
author = "Suter, {James L.} and Lara Kabalan and Mahmoud Khader and Coveney, {Peter V.}",
year = "2015",
month = "11",
day = "15",
doi = "10.1016/j.gca.2015.07.013",
language = "English",
volume = "169",
pages = "17--29",
journal = "Geochmica et Cosmochimica Acta",
issn = "0016-7037",
publisher = "Elsevier Limited",

}

TY - JOUR

T1 - Ab initio molecular dynamics study of the interlayer and micropore structure of aqueous montmorillonite clays

AU - Suter, James L.

AU - Kabalan, Lara

AU - Khader, Mahmoud

AU - Coveney, Peter V.

PY - 2015/11/15

Y1 - 2015/11/15

N2 - Ab initio molecular dynamics simulations have been performed to gain an understanding of the interfacial microscopic structure and reactivity of fully hydrated clay edges. The models studied include both micropore and interlayer water. We identify acidic sites through dissociation mechanisms; the resulting ions can be stabilized by both micropore and interlayer water. We find clay edges possess a complex amphoteric behavior, which depends on the face under consideration and the location of isomorphic substitution. For the neutral (110) surface, we do not observe any dissociation on the timescale accessible. The edge terminating hydroxyl groups participate in a hydrogen bonded network of water molecules that spans the interlayer between periodic images of the clay framework. With isomorphic substitutions in the tetrahedral layer of the (110) clay edge, we find the adjacent exposed apical oxygen behaves as a Brönsted base and abstracts a proton from a nearby water molecule, which in turn removes a proton from an AlOH2 group. With isomorphic substitutions in the octahedral layer of the (110) clay edge the adjacent exposed apical oxygen atom does not abstract a proton from the water molecules, but increases the number of hydrogen bonded water molecules (from one to two). Acid treated clays are likely to have both sites protonated. The (010) surface does not have the same interfacial hydrogen bonding structure; it is much less stable and we observe dissociation of half the terminal SiOH groups (≡Si-O-H→≡SiO-+H+) in our models. The resulting anions are stabilized by solvation from both micropore and interlayer water molecules. This suggests that, when fully hydrated, the (010) surface can act as a Brönsted acid, even at neutral pH.

AB - Ab initio molecular dynamics simulations have been performed to gain an understanding of the interfacial microscopic structure and reactivity of fully hydrated clay edges. The models studied include both micropore and interlayer water. We identify acidic sites through dissociation mechanisms; the resulting ions can be stabilized by both micropore and interlayer water. We find clay edges possess a complex amphoteric behavior, which depends on the face under consideration and the location of isomorphic substitution. For the neutral (110) surface, we do not observe any dissociation on the timescale accessible. The edge terminating hydroxyl groups participate in a hydrogen bonded network of water molecules that spans the interlayer between periodic images of the clay framework. With isomorphic substitutions in the tetrahedral layer of the (110) clay edge, we find the adjacent exposed apical oxygen behaves as a Brönsted base and abstracts a proton from a nearby water molecule, which in turn removes a proton from an AlOH2 group. With isomorphic substitutions in the octahedral layer of the (110) clay edge the adjacent exposed apical oxygen atom does not abstract a proton from the water molecules, but increases the number of hydrogen bonded water molecules (from one to two). Acid treated clays are likely to have both sites protonated. The (010) surface does not have the same interfacial hydrogen bonding structure; it is much less stable and we observe dissociation of half the terminal SiOH groups (≡Si-O-H→≡SiO-+H+) in our models. The resulting anions are stabilized by solvation from both micropore and interlayer water molecules. This suggests that, when fully hydrated, the (010) surface can act as a Brönsted acid, even at neutral pH.

UR - http://www.scopus.com/inward/record.url?scp=84940197400&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84940197400&partnerID=8YFLogxK

U2 - 10.1016/j.gca.2015.07.013

DO - 10.1016/j.gca.2015.07.013

M3 - Article

VL - 169

SP - 17

EP - 29

JO - Geochmica et Cosmochimica Acta

JF - Geochmica et Cosmochimica Acta

SN - 0016-7037

ER -