Search

In the pursue of environmentally-friendly binders for the construction industry, metakaolin (MK) has emerged as promising material, with its hardening performance primarily driven by pozzolanic reactions. However, in systems containing MK and calcium hydroxide (CH), transformation reactions from metastable to stable phases, particularly in excessive CH conditions, can adversely affect material properties. To anticipate these processes, this study introduces a stoichiometry-based reaction modeling approach for pozzolanic MK reactions, encompassing both short-term kinetics and long-term transformation processes. The primary pozzolanic reactions of MK are briefly outlined, highlighting two sequential, partially overlapping reactions forming C4AH13 and C2ASH8. Short-term reaction kinetics are modeled using isothermal calorimetry measurements and deconvoluting the two reaction peaks. The model is in good agreement with experimental quantitative X–ray diffraction (qXRD), thermogravimetric analysis (TGA) and helium pycnometer (density, i.e. solid volume) results. Model limitations are discussed based on qualitative scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM/EDX) analysis. Leveraging the proposed model, the temperature dependency of pozzolanic MK reactions is analyzed, revealing an activation energy for the primary reaction of 84 kJ/mol. The model is intended to lay a foundation for designing innovative binder systems based on metakaolin, while paving the way for sustainable construction in the future.

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.

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. [ABSTRACT FROM AUTHOR]

Cement clinkers containing mainly belite (β-C2S as a model crystal), replacing alite, offer a promising solution for the development of environmentally friendly solutions to reduce the high level of CO2 emissions in the production of Portland cement. However, the much lower reactivity of belite compared to alite limits the widespread use of belite cements. Therefore, this work presents a fundamental atomistic computational approach for comprehending and quantifying the mesoscopic forward dissolution rate of β-C2S, applied to two reactive crystal facets of (100) and ( 1 ¯ 00 ). For this, an atomistic kinetic Monte Carlo (KMC) upscaling approach for cement clinker was developed. It was based on the calculated activation energies (ΔG*) under far-from-equilibrium conditions obtained by a molecular dynamic simulation using the combined approach of ReaxFF and metadynamics, as described in the Part 1 paper in this Special Issue. Thus, the individual atomistic dissolution rates were used as input parameters for implementing the KMC upscaling approach coded in MATLAB to study the dissolution time and morphology changes at the mesoscopic scale. Four different cases and 21 event scenarios were considered for the dissolution of calcium atoms (Ca) and silicate monomers. For this purpose, the (100) and ( 1 ¯ 00 ) facets of a β-C2S crystal were considered using periodic boundary conditions (PBCs). In order to demonstrate the statistical nature of the KMC approach, 40 numerical realizations were presented. The major findings showed a striking layer-by-layer dissolution mechanism in the case of an ideal crystal, where the total dissolution rate was limited by the much slower dissolution of the silicate monomer compared to Ca. The introduction of crystal defects, namely cutting the edges at two crystal boundaries, increased the overall average dissolution rate by a factor of 519. [ABSTRACT FROM AUTHOR]

Portlandite, as a most soluble cement hydration reaction product, affects mechanical and durability properties of cementitious materials. In the present work, an atomistic kinetic Monte Carlo (KMC) upscaling approach is implemented in MATLAB code in order to investigate the dissolution time and morphology changes of a hexagonal platelet portlandite crystal. First, the atomistic rate constants of individual Ca dissolution events are computed by a transition state theory equation based on inputs of the computed activation energies (ΔG*) obtained through the metadynamics computational method (Part 1 of paper). Four different facets (100 or 1 ¯ 00 , 010 or 0 1 ¯ 0, 1 ¯ 10 or 1 1 ¯ 0 , and 001 or 00 1 ¯ ) are considered, resulting in a total of 16 different atomistic event scenarios. Results of the upscaled KMC simulations demonstrate that dissolution process initially takes place from edges, sides, and facets of 010 or 0 1 ¯ 0 of the crystal morphology. The steady-state dissolution rate for the most reactive facets (010 or 0 1 ¯ 0 ) was computed to be 1.0443 mol/(s cm2); however, 0.0032 mol/(s cm2) for 1 ¯ 10 or 1 1 ¯ 0 , 2.672 × 10−7 mol/(s cm2) for 001 or 00 1 ¯ , and 0.31 × 10−16 mol/(s cm2) for 100 or 1 ¯ 00 were represented in a decreasing order for less reactive facets. Obtained upscaled dissolution rates between each facet resulted in a huge (16 orders of magnitude) difference, reflecting the importance of crystallographic orientation of the exposed facets. [ABSTRACT FROM AUTHOR]

[Display omitted] • Stoichiometric model combines kinetics of pozzolanic metakaolin reactions with long-term transformation processes. • Primary pozzolanic metakaolin reactions are divided into C 4 AH 13 and C 2 ASH 8 contributions. • Activation energy for the main (C 2 ASH 8) pozzolanic metakaolin reaction is determined as 84 kJ/mol. • KOH accelerates reaction processes and reduces the activation energy of the pozzolanic metakaolin reaction to 69 kJ/mol. In the pursue of environmentally-friendly binders for the construction industry, metakaolin (MK) has emerged as promising material, with its hardening performance primarily driven by pozzolanic reactions. However, in systems containing MK and calcium hydroxide (CH), transformation reactions from metastable to stable phases, particularly in excessive CH conditions, can adversely affect material properties. To anticipate these processes, this study introduces a stoichiometry-based reaction modeling approach for pozzolanic MK reactions, encompassing both short-term kinetics and long-term transformation processes. The primary pozzolanic reactions of MK are briefly outlined, highlighting two sequential, partially overlapping reactions forming C 4 AH 13 and C 2 ASH 8. Short-term reaction kinetics are modeled using isothermal calorimetry measurements and deconvoluting the two reaction peaks. The model is in good agreement with experimental quantitative X–ray diffraction (qXRD), thermogravimetric analysis (TGA) and helium pycnometer (density, i.e. solid volume) results. Model limitations are discussed based on qualitative scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM/EDX) analysis. Leveraging the proposed model, the temperature dependency of pozzolanic MK reactions is analyzed, revealing an activation energy for the primary reaction of 84 kJ/mol. The model is intended to lay a foundation for designing innovative binder systems based on metakaolin, while paving the way for sustainable construction in the future. [ABSTRACT FROM AUTHOR]