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.

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 $\overline{1}00$, 010 or 0$\overline{1}$0, $\overline{1}10$ or $1\overline{1}0$, and 001 or 00$\overline{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$\overline{1}$0 of the crystal morphology. The steady-state dissolution rate for the most reactive facets (010 or 0$\overline{1}0$) was computed to be 1.0443 mol/(s cm²); however, 0.0032 mol/(s cm²) for $\overline{1}10$ or $1\overline{1}0$, 2.672 × 10⁻⁷ mol/(s cm²) for 001 or 00$\overline{1}$, and 0.31 × 10⁻¹⁶ mol/(s cm²) for 100 or $\overline{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.