Revealing the Effect of Irradiation on Cement Hydrates: Evidence of a Topological Self-Organization

Role of steric repulsions on the precipitation kinetics and the structure of calcium-silicate-hydrate gels

The microstructure and properties of calcium-silicate-hydrate (C-S-H) gels are largely controlled by the physicochemical environment during their precipitation. However, the role of the steric repulsive environment induced by the pore solution chemistry on the kinetics, structure, and properties of C-S-H gels remains unclear. Here, we develop two potential formalisms, namely sinusoidal and polynomial, to simulate the role of steric repulsions in C-S-H. The results show excellent agreement with experimental observations of precipitation kinetics and elastic properties. We demonstrate that the repulsive interactions result in delayed precipitation and percolation, and an open and branched microstructure. Interestingly, the elastic properties (which are equilibrium properties) are also significantly affected by these second-neighbor interactions. Overall, the present study demonstrates that the kinetics, structure, and equilibrium properties of colloidal gels are controlled by the steric repulsions induced by the chemical environment.

The energy landscape governs ductility in disordered materials

Based on their structure, non-crystalline phases can fail in a brittle or ductile fashion. However, the nature of the link between structure and propensity for ductility in disordered materials has remained elusive. Here, based on molecular dynamics simulations of colloidal gels and silica glasses, we investigate how the degree of structural disorder affects the fracture of disordered materials. As expected, we observe that structural disorder results in an increase in ductility. By applying the activation-relaxation technique (an open-ended saddle point search algorithm), we demonstrate that the propensity for ductility is controlled by the topography of the energy landscape. Interestingly, we observe a power-law relationship between the particle non-affine displacement upon fracture and the average local energy barrier. This reveals that the dynamics of the particles upon fracture is encoded in the static energy landscape, i.e., before any load is applied. This relationship is shown to apply to several classes of non-crystalline materials (oxide and metallic glasses, amorphous solid, and colloidal gels), which suggests that it may be a generic feature of disordered materials.

Elucidating the constitutive relationship of calcium-silicate-hydrate gel using high throughput reactive molecular simulations and machine learning

Prediction of material behavior using machine learning (ML) requires consistent, accurate, and, representative large data for training. However, such consistent and reliable experimental datasets are not always available for materials. To address this challenge, we synergistically integrate ML with high-throughput reactive molecular dynamics (MD) simulations to elucidate the constitutive relationship of calcium–silicate–hydrate (C–S–H) gel—the primary binding phase in concrete formed via the hydration of ordinary portland cement. Specifically, a highly consistent dataset on the nine elastic constants of more than 300 compositions of C–S–H gel is developed using high-throughput reactive simulations. From a comparative analysis of various ML algorithms including neural networks (NN) and Gaussian process (GP), we observe that NN provides excellent predictions. To interpret the predicted results from NN, we employ SHapley Additive exPlanations (SHAP), which reveals that the influence of silicate network on all the elastic constants of C–S–H is significantly higher than that of water and CaO content. Additionally, the water content is found to have a more prominent influence on the shear components than the normal components along the direction of the interlayer spaces within C–S–H. This result suggests that the in-plane elastic response is controlled by water molecules whereas the transverse response is mainly governed by the silicate network. Overall, by seamlessly integrating MD simulations with ML, this paper can be used as a starting point toward accelerated optimization of C–S–H nanostructures to design efficient cementitious binders with targeted properties.

Dynamics of confined water and its interplay with alkali cations in sodium aluminosilicate hydrate gel: insights from reactive force field molecular dynamics

This paper presents the dynamics of confined water and its interplay with alkali cations in disordered sodium aluminosilicate hydrate (N-A-S-H) gel using reactive force field molecular dynamics. N-A-S-H gel is the primary binding phase in geopolymers formed via alkaline activation of fly ash. Despite attractive mechanical properties, geopolymers suffer from durability issues, particularly the alkali leaching problem which has motivated this study. Here, the dynamics of confined water and the mobility of alkali cations in N-A-S-H is evaluated by obtaining the evolution of mean squared displacements and Van Hove correlation function. To evaluate the influence of the composition of N-A-S-H on the water dynamics and diffusion of alkali cations, atomistic structures of N-A-S-H with Si/Al ratio ranging from 1 to 3 are constructed. It is observed that the diffusion of confined water and sodium is significantly influenced by the Si/Al ratio. The confined water molecules in N-A-S-H exhibit a multistage dynamic behavior where they can be classified as mobile and immobile water molecules. While the mobility of water molecules gets progressively restricted with an increase in Si/Al ratio, the diffusion coefficient of sodium also decreases as the Si/Al ratio increases. The diffusion coefficient of water molecules in the N-A-S-H structure exhibits a lower value than those of the calcium-silicate-hydrate (C-S-H) structure. This is mainly due to the random disordered structure of N-A-S-H as compared to the layered C-S-H structure. To further evaluate the influence of water content in N-A-S-H, atomistic structures of N-A-S-H with water contents ranging from 5-20% are constructed. Qⁿ distribution of the structures indicates significant depolymerization of N-A-S-H structure with increasing water content. Increased conversion of Si-O-Na network to Si-O-H and Na-OH components with an increase in water content helps explain the alkali-leaching issue in fly ash-based geopolymers observed macroscopically. Overall, the results in this study can be used as a starting point towards multiscale simulation-based design and development of durable geopolymers.