Influence of Octahedral Cation Distribution in Montmorillonite on Interlayer Hydrogen Counter-Ion Retention Strength via First-Principles Calculations

Abstract

Although multiple types of adsorption sites have long been observed in montmorillonite, a consistent explanation about the chemical structure of these adsorption sites has not yet been established. Identifying the cation interlayer adsorption sites based on the octahedral cation distribution on montmorillonite was investigated in this study by using a Density Functional Theory (DFT) simulation. A clay structural model (H[Al6MgFe]Si16O40(OH)8) with a similar composition to Wyoming SWy-1 montmorillonite was built, where two octahedral Al were respectively substituted by Fe and Mg, and H+ was the charge compensating cation. This model had twenty-one different possible configurations as a function of the distribution of octahedral Al, Fe, and Mg cations. The DFT simulations of 15 of these different configurations showed no preference for the formation of any configuration with a specific octahedral Fe-Mg distance. However, the H+ adsorption energy was separated into three distinct groups based on the number of octahedral jumps from Fe to Mg atoms. The H+ adsorption energy significantly decreased with increasing number of octahedral jumps from Fe to Mg. Assuming an even probability of occurrence of 21 octahedral structures in montmorillonite, the percentages of these three groups are 43, 43, and 14%, respectively, which are very close to the three major sites on montmorillonite from published cation adsorption data. These DFT simulations offer an entirely new explanation for the location and chemical structure of the three major adsorption sites on montmorillonite, namely, all three sites are in the interlayer, and their adsorption strengths are a function of the number of octahedral jumps from Fe to Mg atoms.

New Approach To Understanding the Experimental 133Cs NMR Chemical Shift of Clay Minerals via Machine Learning and DFT-GIPAW Calculations

“…The Cs positions for the targeted chemical shifts are then output. Inverse analysis is applied to montmorillonite, and the resultant Cs positions are found to be consistent with previous results (Ohkubo, T.; et al J. Phys. Chem.…”

“…Computational approaches to the analysis of Cs adsorption on clay minerals have mainly been carried out using molecular dynamics (MD) simulations − and ab initio calculations. − The former directly simulates the dynamical behavior of Cs in water–clay systems, adsorption at edge sites, and exchange between inter- and intralayer spaces. Improved computational performance and highly accurate potential functions have made MD simulations involving many clay layers possible.…”

High-temperature transformations of Zr-pillared bentonite

“…The development of mullite, which is usually observed in thermally treated clays, was not observed on this sample. This can be explained by the low amount of Al in the clay mineral due to a high amount of isomorphic replacement of structural Al +3 with Mg +2 in the octahedral sheet [33,34], which favored instead the development of cordierite.…”

Effect of clay mineralogical composition on the structural and textural properties after hydrochloric acid treatment

“…The layer charge is one of the most important parameters of clay minerals. It determines the material properties, which are important for industrial applications (Yoshimoto et al, 2005;Sakuma et al, 2017;Li et al, 2019;Xiang et al, 2019).…”