Next Article in Journal
Synthesis and Physical Characteristics of Undoped and Potassium-Doped Cubic Tungsten Trioxide Nanowires through Thermal Evaporation
Next Article in Special Issue
First-Principles Study on Mechanical, Electronic, and Magnetic Properties of Room Temperature Ferromagnetic Half-Metal MnNCl Monolayer
Previous Article in Journal
Development of Carbon-Nanodot-Loaded PLA Nanofibers and Study of Their Barrier Performance for Medical Applications
Previous Article in Special Issue
Simulation of a Steep-Slope p- and n-Type HfS2/MoTe2 Field-Effect Transistor with the Hybrid Transport Mechanism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Silicate Dissolution Mechanism from Metakaolinite Using Density Functional Theory

by
Mohammadreza Izadifar
*,
Neven Ukrainczyk
* and
Eduardus Koenders
Institute of Construction and Building Materials, Technical University of Darmstadt, Franziska-Braun-Str. 3, 64287 Darmstadt, Germany
*
Authors to whom correspondence should be addressed.
Nanomaterials 2023, 13(7), 1196; https://doi.org/10.3390/nano13071196
Submission received: 22 February 2023 / Revised: 17 March 2023 / Accepted: 26 March 2023 / Published: 27 March 2023
(This article belongs to the Special Issue First-Principles Investigations of Low-Dimensional Nanomaterials)

Abstract

:
Metakaolin (MK) is a high-quality, reactive nanomaterial that holds promising potential for large-scale use in improving the sustainability of cement and concrete production. It can replace cement due to its pozzolanic reaction with calcium hydroxide and water to form cementitious compounds. Therefore, understanding the dissolution mechanism is crucial to fully comprehending its pozzolanic reactivity. In this study, we present an approach for computing the activation energies required for the dissolution of metakaolin (MK) silicate units at far-from-equilibrium conditions using the improved dimer method (IDM) and the transition-state theory (TST) within density functional theory (DFT). Four different models were prepared to calculate the activation energies required for breaking oxo-bridging bonds between silicate or aluminate units. Our results showed that the activation energy for breaking the oxo-bridging bond to a silicate neighbor is higher than that to an aluminate neighbor due to the ionic interaction. However, for complete silicate tetrahedra dissolution, a higher activation energy is required for breaking the oxo-bridging bond to the aluminate neighbor compared to the silicate neighbor. The findings provide methodology for missing input data to predict the mesoscopic dissolution rate, e.g., by the atomistic kinetic Monte Carlo (KMC) upscaling approach.

1. Introduction

Metakaolinite is a type of nanomaterial due to its small particle size, typically in the range of 10 to 100 nanometers, which gives it unique reactivity properties. It is a derivative of kaolinite that is produced by heating kaolinite to temperatures between 500 and 800 degrees Celsius, which causes it to undergo a phase transformation. All this makes it useful in various applications, such as in the production of cementitious, geopolymer, and catalyst materials. Its small particle size also makes it a suitable candidate for use in nanocomposites and nanofillers. In this paper, the focus is on reactivity of metakaolin in cement- and geopolymer-based materials, which involves the dissolution mechanism in alkaline pore solutions. The annual production of more than 20 billion tons of concrete contributes 6–8% to the worldwide anthropogenic carbon dioxide production [1]; in other words, 7% of all greenhouse gases released worldwide [2]. Nevertheless, Portland clinkers, which produce high levels of CO2 emissions, can be partially substituted with supplementary cementitious materials (SCMs) to reduce those emissions by 30–40% without substantially affecting the strength, durability, performance, or cost of the material [3]. Supplementary cementitious materials (SCMs) that are used in the cement industry generally comprise industrial waste products, natural pozzolans, and activated clay minerals that exhibit either pozzolanic or hydraulic properties. Kaolin calcination produces millions of tons of metakaolinite annually, which serves as a SCM, geopolymer precursor, and concrete nano-additive. The major mineral phase, “ideal or ordered kaolinite” (Al2(OH)4Si2O5), dehydroxylates into metakaolinite (Al2Si2O7) [4] at 500–700 °C (Figure 11 reported in [5]). “Low-crystallinity” or “disordered kaolinites” with stacking faults dehydroxylate into metadiskaolinite at lower temperatures (400–700 °C, Figure 12 reported in [5]). Complete dehydroxylation at stationary points causes metakaolinite and metadiskaolinite supercells to expand by around 4% and 8%, respectively, compared to ideal and disordered kaolinite, as shown in Figure 18 from the study reported by Izadifar et al. [5].
Garg and Skibsted’s recent study [6] measured the MK initial dissolution rate for 0.10 mol/L NaOH and longer-term reaction behavior. The initial dissolution rate remained constant at 38 µmol/L/h for both Al and Si for the first 24 h, in far-from-equilibrium (highly diluted) conditions. Subsequently, the dissolution rate slowed as it approached equilibrium, highlighting the need to distinguish the fundamental contributions of forward (and reverse) dissolution reactions from other precipitation reactions. Briki et al. [7] found that the slag reaction slowed in both the early and late stages. They observed a non-steady state dissolution stage for aluminum, silicon, and sulfur ions in the first 30 min of a slag-NaOH solution, followed by a steady stage and then a decrease in release rates for aluminum and silicon ions due to precipitation. The release rates of aluminum and silicon were initially identical but started to deviate after two days. Slag reactivity was six times higher in NaOH solution than in cement paste after one day, and the slower reaction in cement paste was attributed to the presence of calcium ions from clinker phases. On two kaolinite types, N’Guessan et al. [6] studied aluminum and silicon leaching in different alkali hydroxide solutions. Dissolution of kaolinites was related to the hydroxide concentration, with Na+ inducing more dissociation than K+ due to its higher charge density. Werling et al. [8] reported the amount of different phases formed during geopolymerization by determining the solubility of metakaolinite and amorphous SiO2 in NaOH solution at different concentrations. The increment of solubility of metakaolinite calcined at 700 °C has been observed with the increasing concentration of NaOH. In metakaolinite, solubilities of Si < 5% were reached after 24 h and even 7 days in 0.01 and 0.1 mol/L NaOH. In 10.79 mol/L NaOH, the solubilities of Si for the same metakaolinite reached 65% and 80% after 24 h and 7 days, respectively.
Atomistic computational methods have emerged recently as a powerful approach for understanding the microstructure of cement [9,10,11,12] and glasses [13] and their relationship to reactivity. Molecular dynamic (MD) [14,15,16], classical MD, density functional theory (DFT), and ab initio MD computations are extensively used in chemistry and materials science to investigate properties such as reaction mechanisms, mechanical properties computations [17,18], and activation and bonding energies by utilizing force fields or electronic density, respectively, to derive the system’s total energy. Kurganskaya and Luttge [19] used the kinetic Monte Carlo approach to study the dissolution of phyllosilicates and predict the temporal surface evolution dynamics from the kinetics of elementary surface reactions. Their findings suggested that the dissolution of these minerals occurs through the formation of etch pits at opened screw dislocation cores via the movement and coalescence of stepwaves. More recently, Izadifar et al. [20] investigated the nucleation of alkaline aluminosilicate gels based on the binding energies of four different monomer species, including dissolution of metakaolinite [21], using a coarse-grained Monte Carlo (CGMC) approach. However, they did not consider activation energies. Morrow et al. [22] investigated dissolution mechanisms for Al- and Si-terminated sites in protonated, neutral, and deprotonated states in aluminosilicate minerals using density functional theory calculations. Schliemann and Churakov [23,24] employed an ab initio meta-dynamics simulation approach to comprehensively investigate the atomistic scale mechanism for pyrophyllite dissolution on the (010) surface and at the (110) edge surface. Tian et al. [25] conducted DFT calculations on the crystalline and amorphous forms of the anorhite calcium-aluminosilicate. Finally, Salha et al. [26] demonstrated that Al substitution promotes alumino-silicate chain growth in neutral conditions. Their findings suggest that substituting Si with Al-IV, V, and VI units can lower the energy barriers for generating trimers and pentamers.
Despite extensive research in the field of clays, no atomistic computational approach has been reported, to date, for studying the reactivity of protonated silicate tetrahedra in metakaolinite (calcined) clay while accounting for the influence of the existing neighboring aluminate and silicate units. Nevertheless, DFT computational silicate chemistry has been well researched. Xiao et al. studied gas-phase hydrolysis of silicic acid dimer in the presence of OH, finding that a 5-coordinated Si species forms with a 79 kJ/mol barrier, followed by dissociation into two silicic acid molecules with a 19 kJ/mol barrier. Pelmenschikov [27] found that Si–O–Si formation requires about 120 kJ/mol in the hydrolysis of silica in neutral and acidic solutions. Gong and White [13] used a force field molecular dynamics (MD) simulations approach to generate detailed structural representations for amorphous CaO-MgO-SiO2-Al2O3 and CaO-SiO2-Al2O3 glasses for their applications as SCMs in blended Portland cements. They have illustrated that the glass structural representations obtained through the MD simulation approach are in good agreement with the experimental results of X-ray and neutron pair distribution function data and literature data in terms of the nearest interatomic distance, coordination number, and degree of depolymerization.
After obtaining the atomistic activation energies as input data, Izadifar et al. examined dissolution of portlandite [11] and β-C2S clinker [12] through the atomistic kinetic Monte Karlo (KMC) upscaling approach [9,10]. However, the abovementioned aspects of the MK dissolution reaction mechanism have not been visited. Inspired by the pioneer works on silicates, we performed first-principles-based (DFT) calculations to highlight the mechanisms and kinetics for the MK dissolution, initially focusing on silicate-aluminate interactions in neutral aqueous solution. Hence, the main objective of this study is to compute the enthalpy activation energy (ΔH*) at far-from-equilibrium condition, which is based on the transition-state theory (TST) for the calculation of atomistic reaction rates for silicate tetrahedra dissolution in MK through the DFT computational approach. In general, the presented approach forms a basis for future DFT studies on aluminosilicate dissolution mechanisms. Thereby, our and future findings would provide missing data to feed mesoscopic forward dissolution rate computations, e.g., to run the atomistic kinetic Monte Karlo (KMC) upscaling approach.

2. Materials and Computational Approach

2.1. Structural Preparation

The initial periodic unit cell was expanded by a 1 × 1 × 3 ideal kaolinite supercell, as shown in Figure 1A. The unit cell parameters after expansion for B-vacant 1:1 layer were a = 5.15 Å, b = 8.94 Å, c = 22.18 Å, α = 91.92°, β = 105.04°, γ = 89.79°. Input structure was based on the transformation of ideal kaolinite Al12Si12O30(OH)24 into the partially dehydroxylated ideal kaolinite phase Al12Si12O39(OH)6 that took place after losing nine water molecules (Figure 1B), following Izadifar et al. [5], when the chemical potential of water was −1.93 eV. The final structure of the ideal MK implemented in this study is illustrated in Figure 1C. After losing three more water molecules, it had a water chemical potential of −2.67 eV and the chemical formula Al12Si12O42 with the space group P1.

2.2. Density Functional Theory (DFT) Calculation

The density functional theory (DFT) calculations [28] were carried out employing the Vienna ab initio simulation package (VASP) [29,30,31,32,33,34,35] with the projected-augmented wave (PAW) method [36,37] and pseudopotential were used to define electron-ion interaction. A well-converged plane-wave cutoff energy of 400 eV was employed for the structural relaxations. The electron exchange and correlation function were chosen in the generalized gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) parametrization [38]. The Brillouin zone was sampled using a well-converged k-sampling equivalent given by 1 × 1 × 1 Monkhorst-Pack k-points mech size for the total system [39]. The break condition of 10−5 eV was set for the convergence criterion for the electronic self-consistent cycles. In these calculations, the ionic positions, cell volume and cell shape were relaxed until all the forces were smaller than 10−3 eV/Å. To ensure that the stress had also converged, a separate convergence criterion was employed for the stress. Specifically, the “STRESS-TENSOR” was required to reach a value of 0.01 kbar. A three-dimensional visualization software of (VESTA) was also utilized for the structural analysis of our models [40].
The intention for this study was to find the reaction transition states for silicate dissolution in MK through the dimer method using the DFT computational approach. In this way, and through the dimer method [41], the saddle point associated with the reaction and, consequently, the activation energy computation for silicate tetrahedra dissolution based on the transition state theory (TST) can be obtained. In fact, the dimer method is a technique for the optimization of the transition states at the atomic scale for the computation of a minimum energy pathway and finding the transition state between a reactant and a product state. In VASP calculation, the method has been improved by Heyden et al. [42] and, accordingly, it is called the improved dimer method (IDM). The minimum energy pathway represents how the atoms would evolve (between the given initial and the final states), and where a maximum in the potential energy along that path represents the activation energy of the studied process.
Initially, the individual geometries of the initial (reactants) and the final (products) systems are optimized to minimize their energy. Secondly, the initial dimer direction was provided, and then the vibrational modes of the TS were computed to find the mode of decay direction to be given in the improved dimer calculation (IDM). Finally, and through IDM, an optimization step was taken, and the potential energy was maximized along the unstable direction (i.e., dimer axis) while it was minimized in all other directions. It is also worth mentioning that the ground state coordination structure of all four presented models at the Transition state and final products have been provided in the supplementary materials.

3. Results and Discussions

Dissolution of silicate tetrahedra in MK depends on the different neighbors, representing different bonding via the bridging oxygen between silicate and aluminate unit neighbors, as shown in Figure 2. In fact, dissolution of silicate tetrahedra can be described through four various scenarios (models), depending on which bridging oxygen is hydrolyzed (bond broken and saturated by a proton). As illustrated in Figure 2, the first scenario can take place when the bridging oxygen (i) bonded to the silicate neighbor (Si1) is broken down by protonation through water molecule contribution (Figure 3). The second scenario involves protonation of the bridging oxygen (j) bonded to the neighboring aluminum (Al1) (Figure 4). The third scenario concerns silicate tetrahedra dissolution because of protonation of the bridging oxygen, bonded to the aluminum neighbor (Al1) (Figure 5) for the structure, which was obtained from the 1st scenario. The fourth scenario also concerns silicate tetrahedra dissolution by cause of protonation of the bridging oxygen, bonded to the silicate neighbor (Si1) (Figure 6), from the structure obtained from the 2nd scenario.
Liu et al. [43] have recently investigated the oligomerization (the reverse process of hydrolysis) of poly-silicic acid molecules, proceeding through the lateral attacking and simultaneous proton transfer from the approaching molecule for the formation of a 5-coordinated Si species as the transition state, resulting in the ejection of a water molecule from the formed poly-silicic acid. Therefore, we benchmarked our IDM with the results published by Liu et al., confirming the excellent repeatability of their results. Thus, our IDM implementation has been validated as an accurate method to evaluate the configuration of the structure and the computation of the total energy of the new (MK) system at the saddle point. For this, four different models (scenarios) have been identified in this study, as shown in Figure 3, Figure 4, Figure 5 and Figure 6. Findings are summarized in Table 1, listing for each scenario (model) the results of the (activation) energy barrier (Ea) for the hydrolysis reaction, the energy change of the reaction enthalpy (Ee), the Owater-Ooxo distance between the attacking water and the oxo bridge, and the imaginary frequencies obtained from the first principles-based calculations for confirmation of the transition states. Moreover, all distances between the oxygen from water, Si, Al, and oxo-bridging in the reactants (A) and transition states (B) from Figure 3, Figure 4, Figure 5 and Figure 6 have been reported in Table 2 and Table 3, respectively.
The first scenario is representative of model 1, as shown in Figure 3. Figure 3A shows the ground-state structure of the reactant, including the water molecule, as absorbent with the distance of 3.99 Å (Table 2) from the bridging oxygen, bonded to the silicate neighbor (Si1) from the MK surface (00 1 ). By performing IDM, as shown in Figure 3B, the water molecule rotated 90 degrees out of the plane and approached closer to the bridging oxygen at the distance of 2.48 Å (1.51 Å closer compared to the reactant structure). The activation energy was computed to be 1.753 eV as the energy of the system reached from −178.155 eV for reactant structure to the −176.402 eV for the TS structure. Figure 3C shows that the ground state obtained as an optimized geometric structure of the product was a result of protonation of the bridging oxygen through proton transfer from water and making silicate tetrahedra configuration (4-coordinated Si) for the neighboring silicate by contribution of the hydroxyl group from water.
The second scenario is presented as model 2 (Figure 4), in which it is considered that water initially attacks the bridging oxygen bonded to the aluminum neighbor on the MK surface (001). For the ground-state structure of the reactant, the water molecule is initialized within 2.89 Å from the bridging oxygen, bonded to the aluminum neighbor (Al1) from the MK surface 001, as shown in Figure 4A. After implementation of IDM, water rotated slightly and approached closer to the oxygen bridging, bonded to the aluminum neighbor (Al1) at the distance of 2.33 Å.
The activation energy of 0.891 eV was computed when the energy of the system reached from −177.920 eV for the reactant structure to −177.029 eV for the TS structure. As can be observed from Figure 4C, the ground state was obtained as an optimized geometric structure of the product as a result of protonation of bridging oxygen through proton transfer from water and making silicate tetrahedra configuration (4-coordinated Si) for the neighboring silicate by contribution of a hydroxyl group from water. After protonation of the initial oxo-bridging, depending on the silicate (Si1) or aluminum (Al1) neighbor for two different scenarios represented as models 1 and 2, models 3 and 4 illustrate the silicate tetrahedra dissolution, depending on the protonation of oxo-bridging, also bonded to the silicate or aluminate neighbor units. Model 2(C), with one water molecule involvement, is used as the initial structure for the third scenario, presented as model 3, when the water molecule attacks the bridging oxygen, bonded to the silicate neighbor (Si1), for silicate tetrahedra dissolution on the surface (00 1 ), as shown in Figure 5.
Figure 5A shows the ground-state structure of the reactant presented as model 3 (obtained from model 2(C), as shown in Figure 4), including the water molecule as absorbent, with the distance of 3.74 Å from bridging oxygen, bonded to the silicate neighbor (Si1). When the IDM was implemented, the water molecule rotated almost 90 degrees and came closer to the bridging oxygen at the distance of 2.69 Å. The activation energy for model 3 was equal to 1.420 eV as the energy of the system reached from −192.843 eV for the reactant structure to −191.423 eV for the TS structure. Figure 5C shows the ground-state structure of the product, when the silicate tetrahedra dissolved after complete protonation of oxo-bridging. Concerning the fourth and the last scenario for dissolution of silicate tetrahedra, model 1(C) with the contribution of one water molecule was used as the initial structure (the reactant) for model 4, when the water molecule attacks the bridging oxygen bonded to the aluminum neighbor (Al1) on surface (001) as illustrated in Figure 6. Figure 6A displays that 4.66 Å was computed for the distance between water and bridging oxygen, bonded to the silicate neighbor (Si1) for the ground-state structure of the reactant.
After running the IDM optimization for model 4 (Figure 6), the water molecule reaches very close to the bridging oxygen at the distance of 2.39 Å. The computed activation energy of 1.641 eV was obtained as the energy difference between the reactant system structure (−192.664 eV) and the transition-state structure (ETS = −191.023 eV). The final configuration of Figure 6C is identical with that of Figure 5C for complete dissolution of silicate tetrahedra, representative of the ground-state structure of the product. The findings of this study can be discussed by the following points. The activation energy needed to initiate the reaction concerning model 1 for breaking the oxo-bridging bonded to the silicate neighbor is 0.862 eV (97%) higher than that of the aluminum neighbor through the ionic interaction (model 2). On the contrary, for complete silicate tetrahedra dissolution, 0.221 eV (13%) higher activation energy is needed for the breaking of the oxo-bridging, bonded through the ionic interaction to the aluminum neighbor (model 4), in comparison to the breaking of oxo-bridging, bonded through the covalent bond to the silicate neighbor (model 3). Moreover, model 3 for complete silicate tetrahedra dissolution needs 0.333 eV (23%) less activation energy in comparison to model 1 for protonation of oxo-bridging, bonded to the silicate neighbor. In contrast, model 4 for complete silicate tetrahedra dissolution needs 0.750 eV (84%) more activation energy in comparison to model 2 for protonation of the initial oxo-bridging, bonded to the aluminum neighbor.
Based on the analysis of activation energies, it can be concluded that breaking the initial oxo-bridge with an Al neighbor, as shown in model 2, requires almost half the energy (85.96 kJ/mol) compared to the other presented models, namely 1 (169.1 kJ/mol), 3 (137.009 kJ/mol), and 4 (158.33 kJ/mol). Additionally, the highest activation energy was observed in the breaking of the initial oxo-bridge with a Si neighbor, due to covalent bonding.
Comparing our results with previous studies, Chen and Brantley [44] reported almost the identical activation energy values of 88.9 ± 14.6, and 85.2 kJ/mol for the dissolution of albite at acidic and basic pH, respectively, when breaking the initial oxo-bridge with an Al neighbor, as presented in model 2. Morrow et al. [22] reported a lower ab initio barrier of 79 kJ/mol for breaking of Si-O-Al in the deprotonated state of aluminosilicate mineral.
Furthermore, Nangia and Garrison [45] computed values of 110, 159 (which is almost identical to our computation), and 69 kJ/mol for breaking of Si–O–Si in deprotonated, neutral, and protonated states of aluminosilicate mineral, respectively. Pelmenschikov et al. [27] reported that hydrolysis of the first Si–O–Si bond in a cluster with four Si–O–Si links requires a higher value of 36 kJ/mol (205.20 kJ/mol) in comparison to our computation result. It was also noted that a value of 73.34 kJ/mol is needed to break the last remaining bond.
In general, the presented approach forms a basis for future DFT studies on aluminosilicate dissolution mechanisms, identifying all other possible scenarios needed to provide a comprehensive database for upscaling the dissolution rate simulations, e.g., by developing atomistic kinetic Monte Carlo (KMC) models.

4. Conclusions

In this paper, an improved dimer method (IDM) for optimization of the transition states through a density functional theory (DFT) computational approach was employed to study silicate dissolution of metakaolin in water. For this, four different models (scenarios) have been identified to calculate the activation energies of breaking oxo-bridging, depending on the silicate or aluminate neighbor arrangements. Based on the results of this study, the following conclusions could be drawn.
The activation energy for breaking the oxo-bridging covalent bond between the two silicates is 97% higher than the ionic bond between the silicate and aluminate units. However, for final (complete) silicate tetrahedra dissociation, 13% higher activation energy is needed for the ionic bond with the (last) alumina neighbor compared to the covalent bond with the (last) silicate neighbor. In addition, 23% less activation energy is needed for model 3, for complete silicate tetrahedra dissociation, compared to model 1, for the protonation of the covalent bond between the two silicates. Nevertheless, for complete silicate tetrahedra dissociation, model 4 needs 84% more activation energy compared to model 2 for protonation of the initial oxo-bridging bond with the alumina.
The approach demonstrated in this paper is fundamental for computing activation energies, which could be used, e.g., for mesoscopic dissolution rates using the atomistic kinetic Monte Carlo (KMC) upscaling approach, to link meta-clay crystallographic structure with dissolution reactivity.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano13071196/s1.

Author Contributions

Writing—original draft preparation, visualization, methodology (DFT), software, investigations, formal analysis, validation, data curation, writing—review and editing, M.I.; conceptualization, writing, supervision, validation, methodology, investigation, project acquisition, project administration, N.U.; resources, project acquisition, E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Research Foundation (DFG) under project number 446266595 titled: ‘Multiscale modeling of advanced nano-reinforced geopolymer/CNTs materials.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors.

Acknowledgments

Calculations for this research were conducted on the Lichtenberg high-performance computer of the TU Darmstadt. The financial support by the Open Access Publishing Fund of Technische Universität Darmstadt is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Worrell, E.; Price, L.; Martin, N.; Hendriks, C.; Meida, L.O. Carbon dioxide emissions from the global cement industry. Annu. Rev. Energy Environ. 2001, 26, 303–329. [Google Scholar] [CrossRef]
  2. Davidovits, J. Geopolymers: Inorganic polymeric new materials. J. Therm. Anal. 1991, 37, 1633–1656. [Google Scholar] [CrossRef]
  3. Skibsted, J.; Snellings, R. Reactivity of supplementary cementitious materials (SCMs) in cement blends. Cem. Concr. Res. 2019, 124, 105799. [Google Scholar] [CrossRef]
  4. El-Diadamony, H.; Amer, A.A.; Sokkary, T.M.; El-Hoseny, S. Hydration and characteristics of metakaolin pozzolanic cement pastes. HBRC J. 2018, 14, 150–158. [Google Scholar] [CrossRef] [Green Version]
  5. Izadifar, M.; Thissen, P.; Steudel, A.; Kleeberg, R.; Kaufhold, S.; Kaltenbach, J.; Schuhmann, R.; Dehn, F.; Emmerich, K. Comprehensive examination of dehydroxylation of kaolinite, disordered kaolinite, and dickite: Experimental studies and density functional theory. Clays Clay Miner. 2020, 68, 319–333. [Google Scholar] [CrossRef]
  6. Garg, N.; Skibsted, J. Dissolution kinetics of calcined kaolinite and montmorillonite in alkaline conditions: Evidence for reactive Al(V) sites. J. Am. Ceram. Soc. 2019, 102, 7720–7734. [Google Scholar] [CrossRef]
  7. Briki, Y.; Zajac, M.; Ben Haha, M.; Scrivener, K. Factors affecting the reactivity of slag at early and late ages. Cem. Concr. Res. 2021, 150, 106604. [Google Scholar] [CrossRef]
  8. Werling, N.; Dehn, F.; Krause, F.; Steudel, A.; Schuhmann, R.; Emmerich, K. Solubility of precursors and carbonation of waterglass-free geopolymers. Clays Clay Miner. 2020, 68, 524–531. [Google Scholar] [CrossRef]
  9. Uddin, K.M.S.; Izadifar, M.; Ukrainczyk, N.; Koenders, E.; Middendorf, B. Dissolution of Portlandite in Pure Water: Part 1 Molecular Dynamics (MD) Approach. Materials 2022, 15, 1404. [Google Scholar] [CrossRef]
  10. Uddin, K.M.S.; Izadifar, M.; Ukrainczyk, N.; Koenders, E.; Middendorf, B. Dissolution of β-C2S Cement Clinker: Part 1 Molecular Dynamics (MD) Approach for Different Crystal Facets. Materials 2022, 15, 6388. [Google Scholar] [CrossRef]
  11. Izadifar, M.; Ukrainczyk, N.; Uddin, K.S.; Middendorf, B.; Koenders, E. Dissolution of Portlandite in Pure Water: Part 2 Atomistic Kinetic Monte Carlo (KMC) Approach. Materials 2022, 15, 1442. [Google Scholar] [CrossRef] [PubMed]
  12. Izadifar, M.; Ukrainczyk, N.; Uddin, K.M.S.; Middendorf, B.; Koenders, E. Dissolution of β-C2S Cement Clinker: Part 2 Atomistic Kinetic Monte Carlo (KMC) Upscaling Approach. Materials 2022, 15, 6716. [Google Scholar] [CrossRef]
  13. Gong, K.; White, C.E. Predicting CaO-(MgO)-Al2O3-SiO2 glass reactivity in alkaline environments from force field molecular dynamics simulations. Cem. Concr. Res. 2021, 150, 106588. [Google Scholar] [CrossRef]
  14. Izadifar, M.; Thissen, P.; Abadi, R.; Jam, A.N.; Gohari, S.; Burvill, C.; Rabczuk, T. Fracture toughness of various percentage of doping of boron atoms on the mechanical properties of polycrystalline graphene: A molecular dynamics study. Phys. E Low-Dimens. Syst. Nanostruct. 2019, 114, 113614. [Google Scholar] [CrossRef]
  15. Izadifar, M.; Abadi, R.; Shirazi, A.H.N.; Alajlan, N.; Rabczuk, T. Nanopores creation in boron and nitrogen doped polycrystalline graphene: A molecular dynamics study. Phys. E Low-Dimens. Syst. Nanostruct. 2018, 99, 24–36. [Google Scholar] [CrossRef]
  16. Izadifar, M.; Abadi, R.; Jam, A.N.; Rabczuk, T. Investigation into the effect of doping of boron and nitrogen atoms in the mechanical properties of single-layer polycrystalline graphene. Comput. Mater. Sci. 2017, 138, 435–447. [Google Scholar] [CrossRef]
  17. Mortazavi, B.; Cuniberti, G.; Rabczuk, T. Mechanical properties and thermal conductivity of graphitic carbon nitride: A molecular dynamics study. Comput. Mater. Sci. 2015, 99, 285–289. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, Q.; Mortazavi, B.; Aldakheel, F. Molecular Dynamics Modeling of Mechanical Properties of Polymer Nanocomposites Reinforced by C7N6 Nanosheet. Surfaces 2021, 4, 240–254. [Google Scholar] [CrossRef]
  19. Kurganskaya, I.; Luttge, A. A comprehensive stochastic model of phyllosilicate dissolution: Structure and kinematics of etch pits formed on muscovite basal face. Geochim. Cosmochim. Acta 2013, 120, 545–560. [Google Scholar] [CrossRef]
  20. Izadifar, M.; Valencia, N.C.; Xiao, P.; Ukrainczyk, N.; Koenders, E. 3D Off-Lattice Coarse-Grained Monte Carlo Simulations for Nucleation of Alkaline Aluminosilicate Gels. Materials 2023, 16, 1863. [Google Scholar] [CrossRef]
  21. White, C.E.; Provis, J.L.; Proffen, T.; van Deventer, J.S.J. Molecular mechanisms responsible for the structural changes occurring during geopolymerization: Multiscale simulation. AIChE J. 2012, 58, 2241–2253. [Google Scholar] [CrossRef]
  22. Morrow, C.P.; Nangia, S.; Garrison, B.J. Ab Initio Investigation of Dissolution Mechanisms in Aluminosilicate Minerals. J. Phys. Chem. A 2009, 113, 1343–1352. [Google Scholar] [CrossRef] [PubMed]
  23. Schliemann, R.; Churakov, S.V. Pyrophyllite dissolution at elevated pressure conditions: An ab initio study. Geochim. Cosmochim. Acta 2021, 307, 42–55. [Google Scholar] [CrossRef]
  24. Schliemann, R.; Churakov, S.V. Atomic scale mechanism of clay minerals dissolution revealed by ab initio simulations. Geochim. Cosmochim. Acta 2021, 293, 438–460. [Google Scholar] [CrossRef]
  25. Tian, K.V.; Mahmoud, M.Z.; Cozza, P.; Licoccia, S.; Fang, D.C.; Di Tommaso, D.; Chass, G.A.; Greaves, G.N. Periodic vs. molecular cluster approaches to resolving glass structure and properties: Anorthite a case study. J. Non-Cryst. Solids 2016, 451, 138–145. [Google Scholar] [CrossRef]
  26. Salha, M.S.; Yada, R.Y.; Farrar, D.H.; Chass, G.A.; Tian, K.V.; Bodo, E. Aluminium catalysed oligomerisation in cement-forming silicate systems. Phys. Chem. Chem. Phys. 2023, 25, 455–461. [Google Scholar] [CrossRef]
  27. Pelmenschikov, A.; Strandh, H.; Pettersson, L.G.M.; Leszczynski, J. Lattice Resistance to Hydrolysis of Si−O−Si Bonds of Silicate Minerals: Ab Initio Calculations of a Single Water Attack onto the (001) and (111) β-Cristobalite Surfaces. J. Phys. Chem. B 2000, 104, 5779–5783. [Google Scholar] [CrossRef]
  28. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef] [Green Version]
  29. Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558–561. [Google Scholar] [CrossRef]
  30. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  31. Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  32. Izadifar, M.; Königer, F.; Gerdes, A.; Wöll, C.; Thissen, P. Correlation between Composition and Mechanical Properties of Calcium Silicate Hydrates Identified by Infrared Spectroscopy and Density Functional Theory. J. Phys. Chem. C 2019, 123, 10868–10873. [Google Scholar] [CrossRef]
  33. Izadifar, M.; Dolado, J.S.; Thissen, P.; Ayuela, A. Interactions between Reduced Graphene Oxide with Monomers of (Calcium) Silicate Hydrates: A First-Principles Study. Nanomaterials 2021, 11, 2248. [Google Scholar] [CrossRef]
  34. Kresse, G.; Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 1994, 49, 14251–14269. [Google Scholar] [CrossRef] [PubMed]
  35. Izadifar, M.; Natzeck, C.; Emmerich, K.; Weidler, P.G.; Gohari, S.; Burvill, C.; Thissen, P. Unexpected Chemical Activity of a Mineral Surface: The Role of Crystal Water in Tobermorite. J. Phys. Chem. C 2022, 126, 12405–12412. [Google Scholar] [CrossRef]
  36. Blöchl, B.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [Green Version]
  37. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  38. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  40. Momma, K.; Izumi, F. VESTA: A three-dimensional visualization system for electronic and structural analysis. J. Appl. Crystallogr. 2008, 41, 653–658. [Google Scholar] [CrossRef]
  41. Henkelman, G.; Jónsson, H. A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives. J. Chem. Phys. 1999, 111, 7010–7022. [Google Scholar] [CrossRef]
  42. Heyden, A.; Bell, A.T.; Keil, F.J. Efficient methods for finding transition states in chemical reactions: Comparison of improved dimer method and partitioned rational function optimization method. J. Chem. Phys. 2005, 123, 224101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Liu, X.; Liu, C.; Meng, C. Oligomerization of Silicic Acids in Neutral Aqueous Solution: A First-Principles Investigation. IJMS 2019, 20, 3037. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Chen, Y.; Brantley, S.L. Temperature- and pH-dependence of albite dissolution rate at acid pH. Chem. Geol. 1997, 135, 275–290. [Google Scholar] [CrossRef]
  45. Nangia, S.; Garrison, B.J. Reaction Rates and Dissolution Mechanisms of Quartz as a Function of pH. J. Phys. Chem. A 2008, 112, 2027–2033. [Google Scholar] [CrossRef]
Figure 1. (A) The initial periodic 1 × 1 × 3 ideal kaolinite supercell, Al12Si12O30(OH)24. (B) Partially dehydroxylated ideal kaolinite, Al12Si12O39(OH)6. (C) Metakaolinite, Al12Si12O42. Silicate tetrahedrons are depicted in dark blue, aluminate octahedrons in light blue, oxygen atoms in red, and protons in white.
Figure 1. (A) The initial periodic 1 × 1 × 3 ideal kaolinite supercell, Al12Si12O30(OH)24. (B) Partially dehydroxylated ideal kaolinite, Al12Si12O39(OH)6. (C) Metakaolinite, Al12Si12O42. Silicate tetrahedrons are depicted in dark blue, aluminate octahedrons in light blue, oxygen atoms in red, and protons in white.
Nanomaterials 13 01196 g001
Figure 2. Illustration of the initial scenario for dissolution of silicate tetrahedra (Si(OH)4) depending on the bridging of oxygen to the silicate and aluminate units. Silicate atoms are shown in blue; aluminum in light blue; oxygen in red; and hydrogen protons in white.
Figure 2. Illustration of the initial scenario for dissolution of silicate tetrahedra (Si(OH)4) depending on the bridging of oxygen to the silicate and aluminate units. Silicate atoms are shown in blue; aluminum in light blue; oxygen in red; and hydrogen protons in white.
Nanomaterials 13 01196 g002
Figure 3. Model 1: (A) Ground state obtained as an optimized geometric structure of reactant, including water molecule as absorbent on MK surface (00 1 ) for attacking and breaking the bridging oxygen, bonded to the silicate neighbor (Si1). (B) Transition state at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen to the silicate neighbor. The structural building units are described in caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Figure 3. Model 1: (A) Ground state obtained as an optimized geometric structure of reactant, including water molecule as absorbent on MK surface (00 1 ) for attacking and breaking the bridging oxygen, bonded to the silicate neighbor (Si1). (B) Transition state at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen to the silicate neighbor. The structural building units are described in caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Nanomaterials 13 01196 g003
Figure 4. Model 2: (A) Ground state obtained as an optimized geometric structure of reactant, including water molecule as absorbent on MK surface (001) for breaking the bridging oxygen, bonded to the aluminum neighbor (Al1). (B) Transition state (TS) at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen, bonded to the aluminum neighbor. The structural building units are described in the caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Figure 4. Model 2: (A) Ground state obtained as an optimized geometric structure of reactant, including water molecule as absorbent on MK surface (001) for breaking the bridging oxygen, bonded to the aluminum neighbor (Al1). (B) Transition state (TS) at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen, bonded to the aluminum neighbor. The structural building units are described in the caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Nanomaterials 13 01196 g004
Figure 5. Model 3: (A) Ground state obtained as an optimized geometric structure of reactant, including water molecule as absorbent on MK surface (00 1 ) obtained from model 2(C) for breaking the bridging oxygen, bonded to the aluminum neighbor (Al1). (B) TS at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen, bonded to the aluminum neighbor. The structural building units are described in caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Figure 5. Model 3: (A) Ground state obtained as an optimized geometric structure of reactant, including water molecule as absorbent on MK surface (00 1 ) obtained from model 2(C) for breaking the bridging oxygen, bonded to the aluminum neighbor (Al1). (B) TS at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen, bonded to the aluminum neighbor. The structural building units are described in caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Nanomaterials 13 01196 g005
Figure 6. Model 4: (A) Ground state obtained as an optimized geometric structure of reactant including water molecule as absorbent on MK surface (00 1 ) obtained from model 1(C) for breaking the bridging oxygen, bonded to the silicate neighbor (Si1). (B) TS at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen, bonded to the silicate neighbor. The structure building units are described in caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Figure 6. Model 4: (A) Ground state obtained as an optimized geometric structure of reactant including water molecule as absorbent on MK surface (00 1 ) obtained from model 1(C) for breaking the bridging oxygen, bonded to the silicate neighbor (Si1). (B) TS at the saddle point using improved dimer method. (C) Ground state obtained as an optimized geometric structure of product after protonation of the bridging oxygen, bonded to the silicate neighbor. The structure building units are described in caption of Figure 2. E_reactant, E_transition state, and E_product have been collected in Table 1.
Nanomaterials 13 01196 g006
Table 1. Results for each scenario (model): energy barrier for the hydrolysis reaction (ΔEa = E_transition state − E_reactant), the energy change of reaction enthalpy (ΔEe = E_product − E_reactant), the Owater–Ooxo distance between the attacking water and the oxo bridge, and the imaginary frequencies obtained from first principles-based calculations for confirmation of the transition states.
Table 1. Results for each scenario (model): energy barrier for the hydrolysis reaction (ΔEa = E_transition state − E_reactant), the energy change of reaction enthalpy (ΔEe = E_product − E_reactant), the Owater–Ooxo distance between the attacking water and the oxo bridge, and the imaginary frequencies obtained from first principles-based calculations for confirmation of the transition states.
ModelsΔEa (eV)ΔEa (kJ/mol)E_Reactant (eV)E_Transition State (eV)E_Product (eV)ΔEe (eV)Owater–Ooxo Distance (Transition State) (Å)Imaginary Frequencies (THz)
(1)1.753169.138−178.155−176.402−177.9830.1722.480−74.480
(2)0.89185.968−177.920−177.029−178.044−0.1242.330−74.106
(3)1.420137.009−192.843−191.423−191.9000.9432.690−76.330
(4)1.641158.330−192.664−191.023−191.9000.7642.390−81.593
Table 2. The distances between the oxygen from water, Si, Al, and oxo-bridging in the reactants (A), as shown in Figure 3, Figure 4, Figure 5 and Figure 6.
Table 2. The distances between the oxygen from water, Si, Al, and oxo-bridging in the reactants (A), as shown in Figure 3, Figure 4, Figure 5 and Figure 6.
ModelsO1-i (Å)O1-Si1 (Å)O1-Si2 (Å)Si1-i (Å)Si2-i (Å)O2-j (Å)O2-Si1 (Å)O2-Si2 (Å)O2-Al1 (Å)Al1-j (Å)Si2-j (Å)
(1)3.994.233.791.611.67----1.761.62
(2)---1.621.662.895.583.394.021.781.64
(3)3.744.093.861.601.66------
(4)-----4.666.353.375.971.751.61
Table 3. The distances between the oxygen from water, Si, Al, and oxo-bridging in the transition states (B) as shown in Figure 3, Figure 4, Figure 5 and Figure 6.
Table 3. The distances between the oxygen from water, Si, Al, and oxo-bridging in the transition states (B) as shown in Figure 3, Figure 4, Figure 5 and Figure 6.
ModelsO1-i (Å)O1-Si1 (Å)O1-Si2 (Å)Si1-i (Å)Si2-i (Å)O2-j (Å)O2-Si1 (Å)O2-Si2 (Å)O2-Al1 (Å)Al1-j (Å)Si2-j (Å)
(1)2.482.543.151.651.67----1.781.62
(2)---1.611.662.335.092.793.441.741.63
(3)2.692.893.491.641.67------
(4)-----2.395.112.163.331.791.69
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Izadifar, M.; Ukrainczyk, N.; Koenders, E. Silicate Dissolution Mechanism from Metakaolinite Using Density Functional Theory. Nanomaterials 2023, 13, 1196. https://doi.org/10.3390/nano13071196

AMA Style

Izadifar M, Ukrainczyk N, Koenders E. Silicate Dissolution Mechanism from Metakaolinite Using Density Functional Theory. Nanomaterials. 2023; 13(7):1196. https://doi.org/10.3390/nano13071196

Chicago/Turabian Style

Izadifar, Mohammadreza, Neven Ukrainczyk, and Eduardus Koenders. 2023. "Silicate Dissolution Mechanism from Metakaolinite Using Density Functional Theory" Nanomaterials 13, no. 7: 1196. https://doi.org/10.3390/nano13071196

APA Style

Izadifar, M., Ukrainczyk, N., & Koenders, E. (2023). Silicate Dissolution Mechanism from Metakaolinite Using Density Functional Theory. Nanomaterials, 13(7), 1196. https://doi.org/10.3390/nano13071196

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop