Beyond the Average: Spatial and Temporal Fluctuations in Oxide Glass-Forming Systems.

Atomic structure dictates the performance of all materials systems; the characteristic of disordered materials is the significance of spatial and temporal fluctuations on composition-structure-property-performance relationships. Glass has a disordered atomic arrangement, which induces localized distributions in physical properties that are conventionally defined by average values. Quantifying these statistical distributions (including variances, fluctuations, and heterogeneities) is necessary to describe the complexity of glass-forming systems. Only recently have rigorous theories been developed to predict heterogeneities to manipulate and optimize glass properties. This article provides a comprehensive review of experimental, computational, and theoretical approaches to characterize and demonstrate the effects of short-, medium-, and long-range statistical fluctuations on physical properties (e.g., thermodynamic, kinetic, mechanical, and optical) and processes (e.g., relaxation, crystallization, and phase separation), focusing primarily on commercially relevant oxide glasses. Rigorous investigations of fluctuations enable researchers to improve the fundamental understanding of the chemistry and physics governing glass-forming systems and optimize structure-property-performance relationships for next-generation technological applications of glass, including damage-resistant electronic displays, safer pharmaceutical vials to store and transport vaccines, and lower-attenuation fiber optics. We invite the reader to join us in exploring what can be discovered by going beyond the average.

Effects of temperature and CO2 concentration on the early stage nucleation of calcium carbonate by reactive molecular dynamics simulations.

It is significant to investigate the calcium carbonate (CaCO3) precipitation mechanism during the carbon capture process; nevertheless, CaCO3 precipitation is not clearly understood yet. Understanding the carbonation mechanism at the atomic level can contribute to the mineralization capture and utilization of carbon dioxide, as well as the development of new cementitious materials with high-performance. There are many factors, such as temperature and CO2 concentration, that can influence the carbonation reaction. In order to achieve better carbonation efficiency, the reaction conditions of carbonation should be fully verified. Therefore, based on molecular dynamics simulations, this paper investigates the atomic-scale mechanism of carbonation. We investigate the effect of carbonation factors, including temperature and concentration, on the kinetics of carbonation (polymerization rate and activation energy), the early nucleation of calcium carbonate, etc. Then, we analyze the local stresses of atoms to reveal the driving force of early stage carbonate nucleation and the reasons for the evolution of polymerization rate and activation energy. Results show that the higher the calcium concentration or temperature, the higher the polymerization rate of calcium carbonate. In addition, the activation energies of the carbonation reaction increase with the decrease in calcium concentrations.

Water as a Modifier in a Hybrid Coordination Network Glass.

Chemical diversification of hybrid organic-inorganic glasses remains limited, especially compared to traditional oxide glasses, for which property tuning is possible through addition of weakly bonded modifier cations. In this work, it is shown that water can depolymerize polyhedra with labile metal-ligand bonds in a cobalt-based coordination network, yielding a series of nonstoichiometric glasses. Calorimetric, spectroscopic, and simulation studies demonstrate that the added water molecules promote the breakage of network bonds and coordination number changes, leading to lower melting and glass transition temperatures. These structural changes modify the physical and chemical properties of the melt-quenched glass, with strong parallels to the "modifier" concept in oxides. It is shown that this approach also applies to other transition metal-based coordination networks, and it will thus enable diversification of hybrid glass chemistry, including nonstoichiometric glass compositions, tuning of properties, and a significant rise in the number of glass-forming hybrid systems by allowing them to melt before thermal decomposition.

Electrosteric Control of the Aggregation and Yielding Behavior of Concentrated Portlandite Suspensions.

Portlandite (calcium hydroxide: CH: Ca(OH)2) suspensions aggregate spontaneously and form percolated fractal aggregate networks when dispersed in water. Consequently, the viscosity and yield stress of portlandite suspensions diverge at low particle loadings, adversely affecting their processability. Even though polycarboxylate ether (PCE)-based comb polyelectrolytes are routinely used to alter the particle dispersion state, water demand, and rheology of similar suspensions (e.g., ordinary portland cement suspensions) that feature a high pH and high ionic strength, their use to control portlandite suspension rheology has not been elucidated. This study combines adsorption isotherms and rheological measurements to elucidate the role of PCE composition (i.e., charge density, side chain length, and grafting density) in controlling the extent of PCE adsorption, particle flocculation, suspension yield stress, and thermal response of portlandite suspensions. We show that longer side-chain PCEs are more effective in affecting suspension viscosity and yield stress, in spite of their lower adsorption saturation limit and fractional adsorption. The superior steric hindrance induced by the longer side chain PCEs results in better efficacy in mitigating particle aggregation even at low dosages. However, when dosed at optimal dosages (i.e., a dosage that induces a dynamically equilibrated dispersion state of particle aggregates), different PCE-dosed portlandite suspensions exhibit identical fractal structuring and rheological behavior regardless of the side chain length. Furthermore, it is shown that the unusual evolution of the rheological response of portlandite suspensions with temperature can be tailored by adjusting the PCE dosage. The ability of PCEs to modulate the rheology of aggregating charged particle suspensions can be generally extended to any colloidal suspension with a strong screening of repulsive electrostatic interactions.

Atomistic origin of kinetics in hydrated aluminosilicate gels upon precipitation.

Calcium-alumino-silicate-hydrate (CaO-Al2O3-SiO2-H2O, or C-A-S-H) gel, which is the binding phase of cement-based materials, greatly influences concrete mechanical properties and durability. However, the atomic-scale kinetics of the aluminosilicate network condensation remains puzzling. Here, based on reactive molecular dynamics simulations of C-A-S-H systems formation with varying Al/Ca molar ratios, we study the kinetic mechanism of the hydrated aluminosilicate gels upon precipitation. We show that the condensation activation energy decreases with the Al/Ca molar ratio, which suggests that the concentration of the Al polytopes has a great effect on controlling the kinetics of the gelation reaction. Significantly, we demonstrate that 5-fold Al atoms are mainly forming at high Al/Ca molar ratios since there are insufficient hydrogen cations or extra calcium cations to compensate the negatively charged Al polytopes at high Al/Ca molar ratios during accelerated aging.

New insights into the early stage nucleation of calcium carbonate gels by reactive molecular dynamics simulations.

The precipitation of calcium carbonate (CaCO3) is a key mechanism in carbon capture applications relying on mineralization. In that regard, Ca-rich cementitious binders offer a unique opportunity to act as a large-scale carbon sink by immobilizing CO2 as calcium carbonate by mineralization. However, the atomistic mechanism of calcium carbonate formation is still not fully understood. Here, we study the atomic scale nucleation mechanism of an early stage amorphous CaCO3 gel based on reactive molecular dynamics (MD) simulations. We observe that reactive MD offers a notably improved description of this reaction as compared to classical MD, which allows us to reveal new insights into the structure of amorphous calcium carbonate gels and formation kinetics thereof.

Interatomic potential parameterization using particle swarm optimization: Case study of glassy silica.

Classical molecular dynamics simulations of glassy materials rely on the availability of accurate yet computationally efficient interatomic force fields. The parameterization of new potentials remains challenging due to the non-convex nature of the accompanying optimization problem, which renders the traditional optimization methods inefficient or subject to bias. In this study, we present a new parameterization method based on particle swarm optimization (PSO), which is a stochastic population-based optimization method. Using glassy silica as a case study, we introduce two interatomic potentials using PSO, which are parameterized so as to match structural features obtained from ab initio simulations and experimental neutron diffraction data. We find that the PSO algorithm is highly efficient at searching for and identifying viable potential parameters that reproduce the structural features used as the target in the parameterization. The presented approach is very general and can be easily applied to other interatomic potential parameterization schemes.

Temperature-Induced Aggregation in Portlandite Suspensions.

Temperature is well known to affect the aggregation behavior of colloidal suspensions. This paper elucidates the temperature dependence of the rheology of portlandite (calcium hydroxide: Ca(OH)2) suspensions that feature a high ionic strength and a pH close to the particle's isoelectric point. In contrast to the viscosity of the suspending medium (saturated solution of Ca(OH)2 in water), the viscosity of Ca(OH)2 suspensions is found to increase with elevating temperature. This behavior is shown to arise from the temperature-induced aggregation of polydisperse Ca(OH)2 particulates because of the diminution of electrostatic repulsive forces with increasing temperature. The temperature dependence of the suspension viscosity is further shown to diminish with increasing particle volume fraction as a result of volumetric crowding and the formation of denser fractal structures in the suspension. Significantly, the temperature-dependent rheological response of suspensions is shown to be strongly affected by the suspending medium's properties, including ionic strength and ion valence, which affect aggregation kinetics. These outcomes provide new insights into aggregation processes that affect the temperature-dependent rheology of portlandite-based and similar suspensions that feature strong charge screening behavior.

How clay particulates affect flow cessation and the coiling stability of yield stress-matched cementing suspensions.

The remarkable increase in the flow resistance of dense suspensions can hinder 3D-printing processes on account of flow cessation in the extruder, and filament fragility/rupture following deposition. Understanding the nature of rheological changes that occur is critical to manipulate flow conditions or to dose flow modifiers for 3D-printing. Therefore, this paper elucidates the influences of clay particulates on controlling flow cessation and the shape stability of dense cementing suspensions that typically feature poor printability. A rope coiling method was implemented with varying stand-off distances to probe the buckling stability and tendency to fracture of dense suspensions that undergo stretching and bending during deposition. The contributions of flocculation and short-term percolation due to the kinetics of structure formation to deformation rate were deconvoluted using a stepped isostress method. It is shown that the shear stress indicates a divergence with a power-law scaling when the particle volume fraction approaches the jamming limit; φ → φj ≈ φmax. Such a power-law divergence of the shear stress decreases by a factor of 10 with increasing clay dosage. Such behavior in clay-containing suspensions arises from a decrease in the relative packing fraction (φ/φmax) and the formation of fractally-architected aggregates with stronger interparticle interactions, whose uniform arrangement controls flow cessation in the extruder and suspension homogeneity, thereby imparting greater buckling stability. The outcomes offer new insights for assessing/improving the extrudability and printability behavior during slurry-based 3D-printing process.

Dispersing nano- and micro-sized portlandite particulates via electrosteric exclusion at short screening lengths.

In spite of their high surface charge (zeta potential ζ = +34 mV), aqueous suspensions of portlandite (calcium hydroxide: Ca(OH)2) exhibit a strong tendency to aggregate, and thereby present unstable suspensions. While a variety of commercial dispersants seek to modify the suspension stability and rheology (e.g., yield stress, viscosity), it remains unclear how the performance of electrostatically and/or electrosterically based additives is affected in aqueous environments having either a high ionic strength and/or a pH close to the particle's isoelectric point (IEP). We show that the high native ionic strength (pH ≈ 12.6, IEP: pH ≈ 13) of saturated portlandite suspensions strongly screens electrostatic forces (Debye length: κ-1 = 1.2 nm). As a result, coulombic repulsion alone is insufficient to mitigate particle aggregation and affect rheology. However, a longer-range geometrical particle-particle exclusion that arises from electrosteric hindrance caused by the introduction of comb polyelectrolyte dispersants is very effective at altering the rheological properties and fractal structuring of suspensions. As a result, comb-like dispersants that stretch into the solvent reduce the suspension's yield stress by 5× at similar levels of adsorption as compared to linear dispersants, thus enhancing the critical solid loading (i.e., at which jamming occurs) by 1.4×. Significantly, the behavior of diverse dispersants is found to be inherently related to the thickness of the adsorbed polymer layer on particle surfaces. These outcomes inform the design of dispersants for concentrated suspensions that present strong charge screening behavior.

Exploring the landscape of Buckingham potentials for silica by machine learning: Soft vs hard interatomic forcefields.

Interatomic forcefields for silicate glasses often rely on partial (rather than formal) charges to describe the Coulombic interactions between ions. Such forcefields can be classified as "soft" or "hard" based on the value of the partial charge attributed to Si atoms, wherein softer forcefields rely on smaller partial charges. Here, we use machine learning to efficiently explore the "landscape" of Buckingham forcefields for silica, that is, the evolution of the overall forcefield accuracy as a function of the forcefield parameters. Interestingly, we find that soft and hard forcefields correspond to two distinct, yet competitive local minima in this landscape. By analyzing the structure of the silica configurations predicted by soft and hard forcefields, we show that although soft and hard potentials offer competitive accuracy in describing the short-range order structure, soft potentials feature a higher ability to describe the medium-range order.

On the equivalence of vapor-deposited and melt-quenched glasses.

Vapor deposition can yield glasses that are more stable than those obtained by the traditional melt-quenching route. However, it remains unclear whether vapor-deposited glasses are "allowable" or "forbidden," that is, if they are equivalent to glasses formed by cooling extremely slowly a liquid or if they differ in nature from melt-quenched glasses. Here, based on reactive molecular dynamics simulation of silica glasses, we demonstrate that the allowable or forbidden nature of vapor-deposited glasses depends on the temperature of the substrate and, in turn, is found to be encoded in their medium-range order structure.

Precipitation of calcium-alumino-silicate-hydrate gels: The role of the internal stress.

Concrete gains its strength from the precipitation of a calcium-alumino-silicate-hydrate (C-A-S-H) colloidal gel, which acts as its binding phase. However, despite concrete's ubiquity in the building environment, the atomic-scale mechanism of C-A-S-H precipitation is still unclear. Here, we use reactive molecular dynamics simulations to model the early-age precipitation of a C-A-S-H gel. We find that, upon gelation, silicate and aluminate precursors condensate and polymerize to form an aluminosilicate gel network. Notably, we demonstrate that the gelation reaction is driven by the existence of a mismatch of atomic-level internal stress between Si and Al polytopes, which are initially experiencing some local tension and compression, respectively. The polymerization of Si and Al polytopes enables the release of these competitive stresses.

Modifier clustering and avoidance principle in borosilicate glasses: A molecular dynamics study.

Oxide glasses are typically described as having a random, disordered skeleton of network-forming polyhedra that are depolymerized by network-modifying cations. However, the existence of local heterogeneity or clustering within the network-forming and network-modifying species remains unclear. Here, based on molecular dynamics simulations, we investigate the atomic structure of a series of borosilicate glasses. We show that the network-modifying cations exhibit some level of clustering that depends on composition-in agreement with Greaves' modified random network model. In addition, we demonstrate the existence of some mutual avoidance among network-forming atoms, which echoes the Loewenstein avoidance principle typically observed in aluminosilicate phases. Importantly, we demonstrate that the degree of heterogeneity in the spatial distribution of the network modifiers is controlled by the level of ordering in the interconnectivity of the network formers. Specifically, the mutual avoidance of network formers is found to decrease the propensity for modifier clustering.

Mesoscale texture of cement hydrates.

Strength and other mechanical properties of cement and concrete rely upon the formation of calcium-silicate-hydrates (C-S-H) during cement hydration. Controlling structure and properties of the C-S-H phase is a challenge, due to the complexity of this hydration product and of the mechanisms that drive its precipitation from the ionic solution upon dissolution of cement grains in water. Departing from traditional models mostly focused on length scales above the micrometer, recent research addressed the molecular structure of C-S-H. However, small-angle neutron scattering, electron-microscopy imaging, and nanoindentation experiments suggest that its mesoscale organization, extending over hundreds of nanometers, may be more important. Here we unveil the C-S-H mesoscale texture, a crucial step to connect the fundamental scales to the macroscale of engineering properties. We use simulations that combine information of the nanoscale building units of C-S-H and their effective interactions, obtained from atomistic simulations and experiments, into a statistical physics framework for aggregating nanoparticles. We compute small-angle scattering intensities, pore size distributions, specific surface area, local densities, indentation modulus, and hardness of the material, providing quantitative understanding of different experimental investigations. Our results provide insight into how the heterogeneities developed during the early stages of hydration persist in the structure of C-S-H and impact the mechanical performance of the hardened cement paste. Unraveling such links in cement hydrates can be groundbreaking and controlling them can be the key to smarter mix designs of cementitious materials.

New insights into the sol-gel condensation of silica by reactive molecular dynamics simulations.

The sol-gel method is an attractive technique to synthesize homogeneous silicate glasses with high purity while relying on a lower synthesis temperature than in the melt-quench method. However, the mechanism and kinetics of the condensation of the silicate network in aqueous solution remain unclear. Here, based on reactive molecular dynamics simulations (ReaxFF), we investigate the sol-gel condensation kinetics of a silica glass. The influence of the potential parametrization and system size is assessed. Our simulation methodology is found to offer good agreement with experiments. We show that the aqueous concentration of the Si(OH)4 precursors and the local degree of polymerization of the Si atoms play a crucial role in controlling the condensation activation energy. Based on our simulations, we demonstrate that the gelation reaction is driven by the existence of some local atomic stress that gets released upon condensation.

The hydrophilic-to-hydrophobic transition in glassy silica is driven by the atomic topology of its surface.

The surface reactivity and hydrophilicity of silicate materials are key properties for various industrial applications. However, the structural origin of their affinity for water remains unclear. Here, based on reactive molecular dynamics simulations of a series of artificial glassy silica surfaces annealed at various temperatures and subsequently exposed to water, we show that silica exhibits a hydrophilic-to-hydrophobic transition driven by its silanol surface density. By applying topological constraint theory, we show that the surface reactivity and hydrophilic/hydrophobic character of silica are controlled by the atomic topology of its surface. This suggests that novel silicate materials with tailored reactivity and hydrophilicity could be developed through the topological nanoengineering of their surface.

New insights into the atomic structure of amorphous TiO2 using tight-binding molecular dynamics.

Amorphous TiO2 (a-TiO2) could offer an attractive alternative to conventional crystalline TiO2 phases for photocatalytic applications. However, the atomic structure of a-TiO2 remains poorly understood with respect to that of its crystalline counterparts. Here, we conduct some classical molecular dynamics simulations of a-TiO2 based on a selection of empirical potentials. We show that, on account of its ability to dynamically assign the charge of each atom based on its local environment, the second-moment tight-binding charge equilibration potential yields an unprecedented agreement with available experimental data. Based on these simulations, we investigate the degree of order and disorder in a-TiO2. Overall, the results suggest that a-TiO2 features a large flexibility in its local topology, which may explain the high sensitivity of its structure to the synthesis method being used.

Thermometer Effect: Origin of the Mixed Alkali Effect in Glass Relaxation.

Despite the dramatic increase of viscosity as temperature decreases, some glasses are known to feature room-temperature relaxation. However, the structural origin of this phenomenon-known as the "thermometer effect"-remains unclear. Here, based on accelerated molecular dynamics simulations of alkali silicate glasses, we show that both enthalpy and volume follow stretched exponential decay functions upon relaxation. However, we observe a bifurcation of their stretching exponents, with ß=3/5 and 3/7 for enthalpy and volume relaxation, respectively, in agreement with Phillips's topological diffusion-trap model. Based on these results, we demonstrate that the thermometer effect is a manifestation of the mixed alkali effect. We show that relaxation is driven by the existence of stressed local structural instabilities in mixed alkali glasses. This driving force is found to be at a maximum when the concentrations of each alkali atom equal each other, which arises from a balance between the concentration of each alkali atom and the magnitude of the local stress that they experience.

Topological Control on the Structural Relaxation of Atomic Networks under Stress.

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.