Understanding, predicting and eventually improving the resistance to fracture of silicate materials is of primary importance to design new glasses that would be tougher, while retaining their transparency. However, the atomic mechanism of the fracture in amorphous silicate materials is still a topic of debate. In particular, there is some controversy about the existence of ductility at the nano-scale during the crack propagation. Here, we present simulations of the fracture of three archetypical silicate glasses using molecular dynamics. We show that the methodology that is used provide realistic values of fracture energy and toughness. In addition, the simulations clearly suggest that silicate glasses can show different degrees of ductility, depending on their composition.

The question of whether glass continues to relax at low temperature is of fundamental and practical interest. Here, we report a novel atomistic simulation method allowing us to directly access the long-term dynamics of glass relaxation at room temperature. We find that the potential energy relaxation follows a stretched exponential decay, with a stretching exponent β=3/5, as predicted by Phillips's diffusion-trap model. Interestingly, volume relaxation is also found. However, it is not correlated to the energy relaxation, but it is rather a manifestation of the mixed alkali effect.

The existence of nanoscale ductility during the fracture of silicate glasses remains controversial. Here, based on molecular dynamics simulations coupled with topological constraint theory, we show that nano-ductility arises from the spatial heterogeneity of the atomic network's rigidity. Specifically, we report that localized floppy modes of deformation in under-constrained regions of the glass enable plastic deformations of the network, resulting in permanent change in bond configurations. Ultimately, these heterogeneous plastic events percolate, thereby resulting in a non-brittle mode of fracture. This suggests that nano-ductility is intrinsic to multi-component silicate glasses having nanoscale heterogeneities.

Fly ash which encompasses a mixture of glassy and crystalline aluminosilicates is an abundant supplementary cementitious material (SCM), valuable for replacing ordinary portland cement (OPC) in the binder fraction in concrete. Because higher OPC replacement levels are desired, it is critically important to better understand and quantify fly ash reactivity. By combining molecular dynamics (MD) simulations and vertical scanning interferometry (VSI), this study establishes that the reactivity of the glassy fractions in a fly ash with water (i.e., their aqueous dissolution rate) is controlled by the number of constraints placed on atoms within the disordered aluminosilicate network. More precisely, an Arrhenius-like dependence of dissolution rates on the atomic network topology is observed. Such topological controls on fly ash reactivity are highlighted for a range of U.S. commercial fly ashes spanning CaO-enriched and SiO2-enriched compositions. The structure-property relationships reported herein establish an improved framework to control and estimate fly ash-cement interactions in concrete.

Concrete, used in the construction of nuclear power plants (NPPs), may be exposed to radiation emanating from the reactor core. Until recently, concrete has been assumed immune to radiation exposure. Direct evidence acquired on Ar(+)-ion irradiated calcite and quartz indicates, on the contrary, that, such minerals, which constitute aggregates in concrete, may be significantly altered by irradiation. More specifically, while quartz undergoes disordering of its atomic structure resulting in a near complete lack of periodicity, calcite only experiences random rotations, and distortions of its carbonate groups. As a result, irradiated quartz shows a reduction in density of around 15%, and an increase in chemical reactivity, described by its dissolution rate, similar to a glassy silica. Calcite however, shows little change in dissolution rate - although its density noted to reduce by ≈9%. These differences are correlated with the nature of bonds in these minerals, i.e., being dominantly ionic or covalent, and the rigidity of the mineral's atomic network that is characterized by the number of topological constraints (nc) that are imposed on the atoms in the network. The outcomes have major implications on the durability of concrete structural elements formed with calcite or quartz bearing aggregates in nuclear power plants.

Upon loading, atomic networks can feature delayed irreversible relaxation. However, the effect of composition and structure on relaxation remains poorly understood. Herein, relying on accelerated molecular dynamics simulations and topological constraint theory, we investigate the relationship between atomic topology and stress-induced structural relaxation, by taking the example of creep deformations in calcium silicate hydrates (C─S─H), the binding phase of concrete. Under constant shear stress, C─S─H is found to feature delayed logarithmic shear deformations. We demonstrate that the propensity for relaxation is minimum for isostatic atomic networks, which are characterized by the simultaneous absence of floppy internal modes of relaxation and eigenstress. This suggests that topological nanoengineering could lead to the discovery of nonaging materials.