Interactions between Reduced Graphene Oxide with Monomers of (Calcium) Silicate Hydrates: A First-Principles Study

Graphene is a two-dimensional material, with exceptional mechanical, electrical, and thermal properties. Graphene-based materials are, therefore, excellent candidates for use in nanocomposites. We investigated reduced graphene oxide (rGO), which is produced easily by oxidizing and exfoliating graphite in calcium silicate hydrate (CSHs) composites, for use in cementitious materials. The density functional theory was used to study the binding of moieties, on the rGO surface (e.g., hydroxyl-OH/rGO and epoxide/rGO groups), to CSH units, such as silicate tetrahedra, calcium ions, and OH groups. The simulations indicate complex interactions between OH/rGO and silicate tetrahedra, involving condensation reactions and selective repairing of the rGO lattice to reform pristine graphene. The condensation reactions even occurred in the presence of calcium ions and hydroxyl groups. In contrast, rGO/CSH interactions remained close to the initial structural models of the epoxy rGO surface. The simulations indicate that specific CSHs, containing rGO with different interfacial topologies, can be manufactured using coatings of either epoxide or hydroxyl groups. The results fill a knowledge gap, by establishing a connection between the chemical compositions of CSH units and rGO, and confirm that a wet chemical method can be used to produce pristine graphene by removing hydroxyl defects from rGO

Comprehensive examination of dehydroxylation of kaolinite, disordered kaolinite, and dickite: Experimental studies and density functional theory

Kaolins and clays are important rawmaterials for production of supplementary cementitious materials and geopolymer precursors through thermal activation by calcination beyond dehydroxylation (DHX). Both types of clay contain different polytypes and disordered structures of kaolinite but little is known about the impact of the layer stacking of dioctahedral 1:1 layer silicates on optimum thermal activation conditions and following reactivity with alkaline solutions. The objective of the present study was to improve understanding of the impact of layer stacking in dioctahedral 1:1 layer silicates on the thermal activation by investigating the atomic structure after dehydroxylation. Heating experiments by simultaneous thermal analysis (STA) followed by characterization of the dehydroxylated materials by nuclear magnetic resonance spectroscopy (NMR) and scanning electron microscopy (SEM) together with first-principles calculations were performed. Density functional theory (DFT) was utilized for correlation of geometry-optimized structures to thermodynamic stability. The resulting volumes of unit cells were compared with data from dilatometry studies. The local structure changes were correlated with experimental results of increasing DHX temperature in the following order: disordered kaolinite, kaolinite, and dickite, whereupon dickite showed two dehydroxylation steps. Intermediate structures were found that were thermodynamically stable and partially dehydroxylated to a degree of DHX of 75% for kaolinite, 25% for disordered kaolinite, and 50% for dickite. These thermodynamically stable, partially dehydroxylated intermediates contained AlV while metakaolinite and metadickite contained only AlIV with a strongly distorted coordination shell. These results indicate strongly the necessity for characterization of the structure of dioctahedral 1:1 layer silicates in kaolins and clays as a key parameter to predict optimized calcination conditions and resulting reactivity

Dissolution of β-C₂S Cement Clinker: Part 2 Atomistic Kinetic Monte Carlo (KMC) Upscaling Approach

Cement clinkers containing mainly belite (β-C₂S as a model crystal), replacing alite, offer a promising solution for the development of environmentally friendly solutions to reduce the high level of CO₂ 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 β-C₂S, 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 β-C₂S 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

Dissolution of β-C₂S Cement Clinker: Part 1 Molecular Dynamics (MD) Approach for Different Crystal Facets

A major concern in the modern cement industry is considering how to minimize the CO₂ footprint. Thus, cements based on belite, an impure clinker mineral (CaO)₂SiO₂ (C₂S in cement chemistry notation), which forms at lower temperatures, is a promising solution to develop eco-efficient and sustainable cement-based materials, used in enormous quantities. The slow reactivity of belite plays a critical role, but the dissolution mechanisms and kinetic rates at the atomistic scale are not known completely yet. This work aims to understand the dissolution behavior of different facets of β-C₂S providing missing input data and an upscaling modeling approach to connect the atomistic scale to the sub-micro scale. First, a combined ReaxFF and metadynamics-based molecular dynamic approach are applied to compute the atomistic forward reaction rates (RD) of calcium (Ca) and silicate species of (100) facet of β-C₂S considering the influence of crystal facets and crystal defects. To minimize the huge number of atomistic events possibilities, a generalized approach is proposed, based on the systematic removal of nearest neighbors’ crystal sites. This enables us to tabulate data on the forward reaction rates of most important atomistic scenarios, which are needed as input parameters to implement the Kinetic Monte Carlo (KMC) computational upscaling approach. The reason for the higher reactivity of the (100) facet compared to the (010) is explained

3D Off-Lattice Coarse-Grained Monte Carlo Simulations for Nucleation of Alkaline Aluminosilicate Gels

This work presents a 3D off-lattice coarse-grained Monte Carlo (CGMC) approach to simulate the nucleation of alkaline aluminosilicate gels, their nanostructure particle size, and their pore size distribution. In this model, four monomer species are coarse-grained with different particle sizes. The novelty is extending the previous on-lattice approach from White et al. (2012 and 2020) by implementing a full off-lattice numerical implementation to consider tetrahedral geometrical constraints when aggregating the particles into clusters. Aggregation of the dissolved silicate and aluminate monomers was simulated until reaching the equilibrium condition of 16.46% and 17.04% in particle number, respectively. The cluster size formation was analyzed as a function of iteration step evolution. The obtained equilibrated nano-structure was digitized to obtain the pore size distribution and this was compared with the on-lattice CGMC and measurement results from White et al. The observed difference highlighted the importance of the developed off-lattice CGMC approach to better describe the nanostructure of aluminosilicate gels

Silicate Dissolution Mechanism from Metakaolinite Using Density Functional Theory

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

Interactions between Reduced Graphene Oxide with Monomers of (Calcium) Silicate Hydrates: A First-Principles Study

Graphene is a two-dimensional material, with exceptional mechanical, electrical, and thermal properties. Graphene-based materials are, therefore, excellent candidates for use in nanocomposites. We investigated reduced graphene oxide (rGO), which is produced easily by oxidizing and exfoliating graphite in calcium silicate hydrate (CSHs) composites, for use in cementitious materials. The density functional theory was used to study the binding of moieties, on the rGO surface (e.g., hydroxyl-OH/rGO and epoxide/rGO groups), to CSH units, such as silicate tetrahedra, calcium ions, and OH groups. The simulations indicate complex interactions between OH/rGO and silicate tetrahedra, involving condensation reactions and selective repairing of the rGO lattice to reform pristine graphene. The condensation reactions even occurred in the presence of calcium ions and hydroxyl groups. In contrast, rGO/CSH interactions remained close to the initial structural models of the epoxy rGO surface. The simulations indicate that specific CSHs, containing rGO with different interfacial topologies, can be manufactured using coatings of either epoxide or hydroxyl groups. The results fill a knowledge gap, by establishing a connection between the chemical compositions of CSH units and rGO, and confirm that a wet chemical method can be used to produce pristine graphene by removing hydroxyl defects from rGO

Interactions between reduced graphene oxide with monomers of (calcium) silicate hydrates: a first-principles study

Graphene is a two-dimensional material, with exceptional mechanical, electrical, and thermal properties. Graphene-based materials are, therefore, excellent candidates for use in nanocomposites. We investigated reduced graphene oxide (rGO), which is produced easily by oxidizing and exfoliating graphite in calcium silicate hydrate (CSHs) composites, for use in cementitious materials. The density functional theory was used to study the binding of moieties, on the rGO surface (e.g., hydroxyl-OH/rGO and epoxide/rGO groups), to CSH units, such as silicate tetrahedra, calcium ions, and OH groups. The simulations indicate complex interactions between OH/rGO and silicate tetrahedra, involving condensation reactions and selective repairing of the rGO lattice to reform pristine graphene. The condensation reactions even occurred in the presence of calcium ions and hydroxyl groups. In contrast, rGO/CSH interactions remained close to the initial structural models of the epoxy rGO surface. The simulations indicate that specific CSHs, containing rGO with different interfacial topologies, can be manufactured using coatings of either epoxide or hydroxyl groups. The results fill a knowledge gap, by establishing a connection between the chemical compositions of CSH units and rGO, and confirm that a wet chemical method can be used to produce pristine graphene by removing hydroxyl defects from rGO.This research was funded by Karlsruhe House of Young Science (KHYS), and Deutsche Forschungsgemeinschaft (DFG). A.A. and J.S.D. also acknowledge funding from the Spanish Ministry of Science and Innovation (grants nos. PID2019-105488GB-I00 and PCI2019-103657 grant), the Gobierno Vasco UPV/EHU (project no. IT-1246-19), and the European Commission NRG-STORAGE project (project no. GA 870114). Research conducted in the scope of the Transnational Common Laboratory (LTC) Aquitaine-Euskadi Network in Green Concrete and Cement-based Materials.Peer reviewe

Dissolution of Portlandite in Pure Water: Part 2 Atomistic Kinetic Monte Carlo (KMC) Approach

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 01¯0, 1¯10 or 11¯0, and 001 or 001¯) 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 01¯0 of the crystal morphology. The steady-state dissolution rate for the most reactive facets (010 or 01¯0) was computed to be 1.0443 mol/(s cm2); however, 0.0032 mol/(s cm2) for 1¯10 or 11¯0, 2.672 × 10−7 mol/(s cm2) for 001 or 001¯, 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