Search

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.

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.

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]