Dynamics of Confined Reactive Water in Smectite Clay–Zeolite Composites

First-principles investigation of positively charged and neutral oxygen vacancies in amorphous silica

The structural parameters, electron localization functions, electron paramagnetic resonance (EPR) parameters, formation energies, and thermodynamic transition levels of various oxygen vacancy defects in amorphous silica are comprehensively and integrally investigated by using density functional theory. The trends of changes in the oxygen vacancy defect structure and electron localization induced by the increase in distance between defective silicon atoms are clearly identified. It is shown that the dimer configuration may be the potential structure of the Eδ' center. For the back-projected unpuckered configuration and the puckered configuration, whose EPR parameters are more consistent with the experimental values of the Eγ' center, the unpaired electron localized on the sp3 hybridized silicon atom is a common feature. Due to the three-coordinated oxygen atom in the forward-oriented configuration, the EPR parameters are closest to those of the Eα' center. Transformations of oxygen vacancy defects under different charge states are studied by sequentially adding and removing electrons. The thermodynamic transition level analysis reveals that the dimer and forward configurations may behave as deep traps for electron accumulation. The back-projected puckered fourfold-coordinated and fivefold-coordinated configurations are comparatively stable and may be able to function as shallow traps for electron transport. The neutral double unpuckered, neutral back-projected puckered fourfold-coordinated, and neutral back-projected unpuckered configurations are more likely to lose electrons during hole trapping. As the bias voltage is repeatedly changed, the defect density of the puckered configuration may reduce, while that of the dimer and unpuckered configuration may take an opposite trend.

Managing Expectations and Imbalanced Training Data in Reactive Force Field Development: An Application to Water Adsorption on Alumina

ReaxFF is a computationally efficient model for reactive molecular dynamics simulations that has been applied to a wide variety of chemical systems. When ReaxFF parameters are not yet available for a chemistry of interest, they must be (re)optimized, for which one defines a set of training data that the new ReaxFF parameters should reproduce. ReaxFF training sets typically contain diverse properties with different units, some of which are more abundant (by orders of magnitude) than others. To find the best parameters, one conventionally minimizes a weighted sum of squared errors over all of the data in the training set. One of the challenges in such numerical optimizations is to assign weights so that the optimized parameters represent a good compromise among all the requirements defined in the training set. This work introduces a new loss function, called Balanced Loss, and a workflow that replaces weight assignment with a more manageable procedure. The training data are divided into categories with corresponding "tolerances", i.e., acceptable root-mean-square errors for the categories, which define the expectations for the optimized ReaxFF parameters. Through the Log-Sum-Exp form of Balanced Loss, the parameter optimization is also a validation of one's expectations, providing meaningful feedback that can be used to reconfigure the tolerances if needed. The new methodology is demonstrated with a nontrivial parametrization of ReaxFF for water adsorption on alumina. This results in a new force field that reproduces both the rare and frequent properties of a validation set not used for training. We also demonstrate the robustness of the new force field with a molecular dynamics simulation of water desorption from a γ-Al2O3 slab model.

Development and validation of a general-purpose ReaxFF reactive force field for earth material modeling

ReaxFF reactive force field bridges the gap between nonreactive molecular simulations and quantum mechanical calculations and has been widely applied during the past two decades. However, its application to earth materials, especially those under high T-P conditions relevant to Earth's interior, is still limited due to the lack of available parameters. Here, we present the development and validation of a ReaxFF force field containing several of the most common elements in Earth's crust, i.e., Si/Al/O/H/Na/K. The force field was trained against a large data set obtained from density functional theory (DFT) calculations, including charges, bond/angle distortion curves, equation of states, ion migration energy profiles, and condensation reaction energies. Different coordination environments were considered in the training set. The fitting results showed that the current force field can well reproduce the DFT data (the Pearson correlation coefficient, Rp, is 0.95). We validated the force field on mineral-water interfaces, hydrous melts/supercritical geofluids, and bulk crystals. It was found that the current force field performed excellently in predicting the structural, thermodynamic, and transport properties of various systems (Rp = 0.95). Moreover, possible applications and future development have been discussed. The results obtained in this study suggest that the current force field holds good promise to model a wide range of processes and thus open opportunities to advance the application of ReaxFF in earth material modeling.

Effects of -H and -OH Termination on Adhesion of Si-Si Contacts Examined Using Molecular Dynamics and Density Functional Theory

The contact between nanoscale single-crystal silicon asperities and substrates terminated with -H and -OH functional groups is simulated using reactive molecular dynamics (MD). Consistent with previous MD simulations for self-mated surfaces with -H terminations only, adhesion is found to be low at full adsorbate coverages, be it self-mated coverages of mixtures of -H and -OH groups, or just -OH groups. As the coverage reduces, adhesion increases markedly, by factors of ∼5 and ∼6 for -H-terminated surfaces and -OH-terminated surfaces, respectively, and is due to the formation of covalent Si-Si bonds; for -OH-terminated surfaces, some interfacial Si-O-Si bonds are also formed. Thus, covalent linkages need to be broken upon separation of the tip and substrate. In contrast, replacing -H groups with -OH groups while maintaining complete coverage leads to negligible increases in adhesion. This indicates that increases in adhesion require unsaturated sites. Furthermore, plane-wave density functional theory (DFT) calculations were performed to investigate the energetics of two Si(111) surfaces fully terminated by either -H or -OH groups. Importantly for the adhesion results, both DFT and MD calculations predict the correct trends for the relative bond strengths: Si-O > Si-H > Si-Si. This work supports the contention that prior experimental work observing strong increases in adhesion after sliding Si-Si nanoasperities over each other is due to sliding-induced removal of passivating species on the Si surfaces.

The inhibitory effect of excess calcium ions on the polymerization process of calcium aluminate silicate hydrate (CASH) gel

Calcium ion, as an essential component in CASH, affects the aggregation and formation process of CASH, thereby influencing its microstructure and mechanical properties. However, the mechanism by which calcium ions affect the polymerization process of CASH is not yet fully understood. In this study, the effects of calcium ions on the polymerization process, coagulation state, and microstructure of CASH are investigated via molecular dynamics simulation. The results indicate that the presence of a trace amount of Ca2+ attracts oligomers towards the calcium-rich region, thus speeding up the polymerization to some extent, but as the Ca2+ content increases, more Ca2+ binds to the oxygen atoms in silica-oxygen tetrahedra and aluminum-oxygen tetrahedra, forming tight ion pairs and occupying the hydroxyl binding sites required for the polycondensation reaction. This inhibits the continuous aggregation of CASH gel and slows down the rate of polymerization. Additionally, Ca2+ attracts oxygen atoms from free water molecules and free OH-, forming Ca(OH)2 dispersed in the spatial structure, which hinders the formation of larger clusters. As a result, the higher the Ca ion content in the system, the lower the overall polymerization degree of the CASH gel, resulting in a decrease in the conversion of the Q1 dimer to Q2 and Q3 chain structures, a shorter average chain length, poorer overall connectivity, and a transition from large clusters in a better-aggregated state to dispersed small clusters. This study sheds light on the polymerization reaction mechanism of CASH gels.

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.

Atomistic simulations of calcium aluminosilicate interfaced with liquid water

The dissolution behavior of calcium aluminosilicate based glass fibers, such as stone wool fibers, is an important consideration in mineral wool applications for both the longevity of the mineral wool products in humid environments and limiting the health impacts of released and inhaled fibers from the mineral wool product. Balancing these factors requires a molecular-level understanding of calcium aluminosilicate glass dissolution mechanisms, details that are challenging to resolve with experiment alone. Molecular dynamics simulations are a powerful tool capable of providing complementary atomistic insights regarding dissolution; however, they require force fields capable of describing not-only the calcium aluminosilicate surface structure but also the interactions relevant to dissolution phenomena. Here, a new force field capable of describing amorphous calcium aluminosilicate surfaces interfaced with liquid water is developed by fitting parameters to experimental and first principles simulation data of the relevant oxide-water interfaces, including ab initio molecular dynamics simulations performed for this work for the wüstite and periclase interfaces. Simulations of a calcium aluminosilicate surface interfaced with liquid water were used to test this new force field, suggesting moderate ingress of water into the porous glass interface. This design of the force field opens a new avenue for the further study of calcium and network-modifier dissolution phenomena in calcium aluminosilicate glasses and stone wool fibers at liquid water interfaces.

Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations

Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.

Mussel-Inspired Polydopamine Composite Mesoporous Bioactive Glass Nanoparticles: An Exploration of Potential Metal-Ion Loading Platform and In Vitro Bioactivity

Exploring new approaches to realize the possibility of incorporating biologically active elements into mesoporous silicate bioactive glass nanoparticles (MBG NPs) and guaranteeing their meso- structural integrity and dimensional stability has become an attractive and interesting challenge in biomaterials science. We present a postgrafting strategy for introducing different metal elements into MBG NPs. This strategy is mediated by polydopamine (PDA) coating, achieving uniform loading of copper or copper-cobalt on the particles efficiently and ensuring the stability of MBG NPs in terms of particle size, mesoporous structure, and chemical structure. However, the PDA coating reduced the ion-binding free energy of the MBG NPs for calcium and phosphate ions, resulting in the deposition of minimal CaP clusters on the PDA@MBG NP surface when immersed for 7 days in simulated body fluid, indicating the absence of hydroxyapatite mineralization.

Ultrahigh-rate and ultralong-life aqueous batteries enabled by special pair-dancing proton transfer

The design of Faradaic battery electrodes with high rate capability and long cycle life comparable to those of supercapacitors is a grand challenge. Here, we bridge this performance gap by taking advantage of a unique ultrafast proton conduction mechanism in vanadium oxide electrode, developing an aqueous battery with untrahigh rate capability up to 1000 C (400 A g^-1) and extremely long life of 0.2 million cycles. The mechanism is elucidated by comprehensive experimental and theoretical results. Instead of slow individual Zn²⁺ transfer or Grotthuss chain transfer of confined H⁺, the ultrafast kinetics and excellent cyclic stability are enabled by rapid 3D proton transfer in vanadium oxide via the special pair dance switching between Eigen and Zundel configurations with little constraint and low energy barriers. This work provides insight into developing high-power and long-life electrochemical energy storage devices with nonmetal ion transfer through special pair dance topochemistry dictated by hydrogen bond.

Microscopic insights into the variations of antibiotics sorption to clay minerals

Understanding the adsorption behavior of antibiotic molecules on minerals is crucial for determining the environmental fate and transport of antibiotics in soils and waters. However, the microscopic mechanisms that govern the adsorption of common antibiotics, such as the molecular orientation during the adsorption process and the conformation of sorbate species, are not well understood. To address this gap, we conducted a series of molecular dynamics (MD) simulations and thermodynamics analyses to investigate the adsorption of two typical antibiotics, tetracycline (TET) and sulfathiazole (ST), on the surface of montmorillonite. The simulation results indicated that the adsorption free energy ranged from - 23 to - 32 kJ·mol^-1, and - 9 to - 18 kJ·mol^-1 for TET and ST, respectively, which was consistent with the measured difference of sorption coefficient (K_d) for TET-montmorillonite of 11.7 L·g^-1 and ST-montmorillonite of 0.014 L·g^-1. The simulations also found that TET was adsorbed through dimethylamino groups (85% in probability) with a molecular conformation vertical to the montmorillonite's surface, while ST was adsorbed through sulfonyl amide group (95% in probability) with vertical, tilted and parallel conformations on the surface. The results confirmed that molecular spatial orientations could affect the adsorption capacity between antibiotics and minerals. Overall, the microscopic adsorption mechanisms revealed in this study provide critical insights into the complexities of antibiotics adsorption to soil and facilitate the prediction of adsorption capacity of antibiotics on minerals and their environmental transport and fate. This study contributes to our understanding of the environmental impacts of antibiotic usage and highlights the importance of considering molecular-level processes when assessing the fate and transport of antibiotics in the environment.

Nanoscale Insight into the Effect of Calcium on Early-Age Polymerization of CNASH Gels

Sodium-containing calcium-alumino-silicate-hydrate (CNASH) gels, the primary binder phase of alkali-activated materials (AAMs), significantly impact the performance of the AAM. Although the effect of the calcium content on the AAM has been extensively studied in the past, few studies focus on the effect of calcium on the structure and performance of gels at a molecular scale. As an important element in gels, the effect of calcium in gels on its atomic-scale properties remains unclear. This study establishes a molecular model of the CNASH gel via reactive molecular dynamics (MD) simulation and verifies the feasibility of the gel model. By employing the reactive MD, the effect of calcium on the physicochemical properties of gels in the AAM is investigated. The simulation highlights that the condensation process of the system containing Ca is accelerated dramatically. This phenomenon is explained from the perspective of thermodynamics and kinetics. The increased calcium content enhances the thermodynamic stability and reduces the energy barrier of the reaction. Then, the phenomenon is further analyzed through the nanosegregation in the structure. It is proved that this behavior is driven by the weaker affinity of calcium for aluminosilicate chains than the particles in the aqueous environment. The difference in affinity leads to nanosegregation in the structure, making Si(OH)₄ and Al(OH)₃ monomers and oligomers closer for better polymerization.

Molecular Insight into the Pozzolanic Reaction of Metakaolin and Calcium Hydroxide

The reaction mechanism of the pozzolanic reaction of metakaolin (MK) from the atomic point of view has not yet been explored. To explain the process and mechanism of the pozzolanic reaction from the atomic point of view, molecular insight into the pozzolanic reaction of MK and calcium hydroxide (CH) was analyzed through the reaction molecular dynamics (MD) simulation. The results show that the pozzolanic reaction of MK and CH can be essentially regarded as the CH decomposition and penetration into MK. Also, the structure evolution after the pozzolanic reaction shows that the water molecules cannot penetrate the MK structure till the participation of Ca²⁺ and OH^- ions of CH. The Ca²⁺ and OH^- ions have strong interaction with MK and drill into the MK structure, followed by the destruction of a part of the MK structure and water penetration. The final structure of CH removed by MK can be regarded as the precursor of the CASH gel structure.

A reactive molecular dynamics study of thermal pyrolysis behavior and mechanisms of lignin during the hydrothermal process: The function of the water molecules

The lignin hydrothermal processing is an important option but a full understanding of the role played by the water molecules in the depolymerization of lignin is still lacking. In order to clarify the role of the water molecules in the depolymerization of lignin, the evolution of chemical bonds, microstructural changes, and possible mechanisms of product generation were compared during the pyrolysis process under vacuum and water conditions using Reactive Molecular Dynamics Simulation. Compared with vacuum conditions, the role of water changes with temperature, identifying three stages: promotion (1200-1800 K)-inhibition (2100-2400 K)-promotion (2700-3000 K). Also compared with vacuum conditions, hydrothermal processing can promote the cleavage of the ether bonds while inhibiting the destruction of carbocycles. Water molecules promote the depolymerization of lignin into more C₄-molecules, thereby generating more combustible gas resources. Based on the research results, the pyrolysis conditions of lignin can be flexibly controlled to obtain solids, liquids or gases.

Critical role of water structure around interlayer ions for ion storage in layered double hydroxides

Water-containing layered materials have found various applications such as water purification and energy storage. The highly structured water molecules around ions under the confinement between the layers determine the ion storage ability. Yet, the relationship between the configuration of interlayer ions and water structure in high ion storage layered materials is elusive. Herein, using layered double hydroxides, we demonstrate that the water structure is sensitive to the filling density of ions in the interlayer space and governs the ion storage. For ion storage of dilute nitrate ions, a 24% decrease in the filling density increases the nitrate storage capacity by 300%. Quartz crystal microbalance with dissipation monitoring studies, combined with multimodal ex situ experiments and theoretical calculations, reveal that the decreasing filling density effectively facilitates the 2D hydrogen-bond networking structure in water around interlayer nitrate ions along with minimal change in the layered structure, leading to the high storage capacity.

Molecular Insights into the Reaction Process of Alkali-Activated Metakaolin by Sodium Hydroxide

When metakaolin (MK) is alkalized with an alkaline activator, it depolymerizes under the action of the alkali. However, the process of MK alkalinization is still unrevealed. Here, we supplied a molecular insight into the process of MK alkalinization through reaction molecular dynamics (MD) simulation. The structure, dynamics, and process of MK alkalinization are systematically investigated. The results showed that the layered structure of MK was destroyed and the silicates in MK were dissolved by sodium hydroxide solution during the alkalinization reaction of MK. The aluminates in MK are not dissolved, indicating that aluminates are more stable than silicates. Moreover, the equilibrium structures of MK with H₂O and MK with NaOH solution show that only when both sodium hydroxide and water are involved in the alkalinization reaction, the silicates in MK undergo depolymerization. Also, the observed final state of MK alkalinization can be recognized as the precursor of alkali-activated materials (AAMs).

Chemical and Reactive Transport Processes Associated with Hydraulic Fracturing of Unconventional Oil/Gas Shales

Hydraulic fracturing of unconventional oil/gas shales has changed the energy landscape of the U.S. Recovery of hydrocarbons from tight, hydraulically fractured shales is a highly inefficient process, with estimated recoveries of <25% for natural gas and <5% for oil. This review focuses on the complex chemical interactions of additives in hydraulic fracturing fluid (HFF) with minerals and organic matter in oil/gas shales. These interactions are intended to increase hydrocarbon recovery by increasing porosities and permeabilities of tight shales. However, fluid-shale interactions result in the dissolution of shale minerals and the release and transport of chemical components. They also result in mineral precipitation in the shale matrix, which can reduce permeability, porosity, and hydrocarbon recovery. Competition between mineral dissolution and mineral precipitation processes influences the amounts of oil and gas recovered. We review the temporal/spatial origins and distribution of unconventional oil/gas shales from mudstones and shales, followed by discussion of their global and U.S. distributions and compositional differences from different U.S. sedimentary basins. We discuss the major types of chemical additives in HFF with their intended purposes, including drilling muds. Fracture distribution, porosity, permeability, and the identity and molecular-level speciation of minerals and organic matter in oil/gas shales throughout the hydraulic fracturing process are discussed. Also discussed are analysis methods used in characterizing oil/gas shales before and after hydraulic fracturing, including permeametry and porosimetry measurements, X-ray diffraction/Rietveld refinement, X-ray computed tomography, scanning/transmission electron microscopy, and laboratory- and synchrotron-based imaging/spectroscopic methods. Reactive transport and spatial scaling are discussed in some detail in order to relate fundamental molecular-scale processes to fluid transport. Our review concludes with a discussion of potential environmental impacts of hydraulic fracturing and important knowledge gaps that must be bridged to achieve improved mechanistic understanding of fluid transport in oil/gas shales.

Ultrafast rectifying counter-directional transport of proton and metal ions in metal-organic framework-based nanochannels

Bioinspired control of ion transport at the subnanoscale has become a major focus in the fields of nanofluidics and membrane separation. It is fundamentally important to achieve rectifying ion-specific transport in artificial ion channels, but it remains a challenge. Here, we report a previously unidentified metal-organic framework nanochannel (MOF NC) nanofluidic system to achieve unidirectional ultrafast counter-directional transport of alkaline metal ions and proton. This highly effective ion-specific rectifying transport behavior is attributed to two distinct mechanisms for metal ions and proton, elucidated by theoretical simulations. Notably, the MOF NC exhibits ultrafast proton conduction stemming from ultrahigh proton mobility, i.e., 11.3 × 10^-7 m² /V·s, and low energy barrier of 0.075 eV in MIL-53-COOH subnanochannels. Furthermore, the MOF NC shows excellent osmotic power-harvesting performance in reverse electrodialysis. This work expects to inspire further research into multifunctional biomimetic ion channels for advanced nanofluidics, biomimetics, and separation applications.

Nanoarchitectonics of Illite-Based Materials: Effect of Metal Oxides Intercalation on the Mechanical Properties

Clay minerals inevitably interact with colloidal oxides (mainly iron and aluminum oxides) in the evolution of natural geomaterials. However, the interaction between the clay minerals and the colloidal oxides affecting the stability and the strength of geotechnical materials remains poorly understood. In the present work, the interaction between the clay minerals and the colloidal oxides was investigated by reaction molecular dynamics simulations to explore the mechanical properties of illite-based materials. It was found that the metal atoms of the intercalated amorphous iron and aluminum oxides interact with oxygen atoms of the silica tetrahedron at the interface generating chemical bonds to enhance the strength of the illite-based materials considerably. The deformation and failure processes of the hybrid illite-based structures illustrated that the Al–O bonds were more favorable to the mechanical properties’ improvement of the hybrid system compared with Fe–O bonds. Moreover, the anisotropy of illite was greatly improved with metal oxide intercalation. This study provides new insight into the mechanical properties’ improvement of clay-based materials through metal oxides intercalation.

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na⁺, Ca²⁺, and Mg²⁺ cations and the CO₃^2- anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me-O-C (Me = metal) three-body valence angle parameters and Me-C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na⁺ predominantly forms carbonates due to weaker angular strain, while Ca²⁺ and Mg²⁺ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. The nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca²⁺ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na⁺ and Mg²⁺ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

Condensation and Growth of Amorphous Aluminosilicate Nanoparticles by Aggregation Process

The precipitation of zeolite nanoparticles involves the initial formation of metastable precursors such as amorphous entities that crystalize through non-classical pathways. Here, using reactive force field based simulations, we reveal...

Transferable Gaussian Attractive Potentials for Organic/Oxide Interfaces

Organic/oxide interfaces play an important role in many areas of chemistry and in particular for lubrication and corrosion. Molecular dynamics simulations are the method of choice for providing complementary insight to experiments. However, the force fields used to simulate the interaction between molecules and oxide surfaces tend to capture only weak physisorption interactions, discarding the stabilizing Lewis acid/base interactions. We here propose a simple complement to the straightforward molecular mechanics description based on "out-of-the-box" Lennard-Jones potentials and electrostatic interactions: the addition of an attractive Gaussian potential between reactive sites of the surface and heteroatoms of adsorbed organic molecules, leading to the Gaussian Lennard-Jones (GLJ) potential. The interactions of four oxygenated and four amine molecules with the typical and widespread hematite and γ-alumina surfaces are investigated. The root mean square deviation (RMSD) for all probed molecules is only 5.7 kcal/mol, which corresponds to an error of 23% over hematite. On γ-alumina, the RMSD is 11.2 kcal/mol using a single parameter for all five chemically inequivalent surface aluminum atoms. Applying GLJ to the simulation of organic films on oxide surfaces demonstrates that the mobility of the surfactants is overestimated by the simplistic LJ potential, while GLJ and other qualitatively correct potentials show a strong structuration and slow dynamics of the surface films, as could be expected from the first-principles adsorption energies for model head groups.

Topological origin of phase separation in hydrated gels

HYPOTHESIS

Depending on their composition, hydrated gels can be homogeneous or phase-separated, which, in turn, affects their dynamical and mechanical properties. However, the nature of the structural features, if any, that govern the propensity for a given gel to phase-separate remains largely unknown. Here, we argue that the propensity for hydrated gels to phase-separate is topological in nature.

SIMULATIONS

We employ reactive molecular dynamics simulations to model the early-age precipitation of calcium-alumino-silicate-hydrate (CASH) gels with varying compositions, i.e., (CaO)_1.7(Al₂O₃)_x(SiO₂)_{1 -x}(H₂O)_{3.7 +x}. By adopting topological constraint theory, we investigate the structural origin of phase separation in hydrated gels.

FINDINGS

We report the existence of a homogeneous-to-phase-separated transition, wherein Si-rich (x ≤ 0.10) CASH gels are homogeneous, whereas Al-rich (x > 0.10) CASH gels tend to phase-separate. Furthermore, we demonstrate that this transition is correlated to a topological flexible-to-rigid transition within the atomic network. We reveal that the propensity for topologically-overconstrained gels to phase-separate arises from the existence of some internal stress within their atomic network, which acts as an energy penalty that drives phase separation.

Diffusion of confined fluids in microporous zeolites and clay materials

Fluids exhibit remarkable variation in their structural and dynamic properties when they are confined at the nanoscopic scale. Various factors, including geometric restriction, the size and shape of the guest molecules, the topology of the host, and guest-host interactions, are responsible for the alterations in these properties. Due to their porous structures, aluminosilicates provide a suitable host system for studying the diffusion of sorbates in confinement. Zeolites and clays are two classes of the aluminosilicate family, comprising very ordered porous or layered structures. Zeolitic materials are important due to their high catalytic activity and molecular sieving properties. Guest molecules adsorbed by zeolites display many interesting features including unidimensional diffusion, non-isotropic rotation, preferred orientation and levitation effects, depending on the guest and host characteristics. These are useful for the separation of hydrocarbons which commonly exist as mixtures in nature. Similarly, clay materials have found application in catalysis, desalination, enhanced oil recovery, and isolation barriers used in radioactive waste disposal. It has been shown that the bonding interactions, level of hydration, interlayer spacing, and number of charge-balancing cations are the important factors that determine the nature of diffusion of water molecules in clays. Here, we present a review of the current status of the diffusion mechanisms of various adsorbed species in different microporous zeolites and clays, as investigated using quasielastic neutron scattering and classical molecular dynamics simulation techniques. It is impossible to write an exhaustive review of the subject matter, as it has been explored over several decades and involves many research topics. However, an effort is made to cover the relevant issues specific to the dynamics of different molecules in microporous zeolites and clay materials and to highlight a variety of interesting features that are important for both practical applications and fundamental aspects.

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.

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.

Radiocesium interaction with clay minerals: Theory and simulation advances Post-Fukushima

Insights at the microscopic level of the process of radiocesium adsorption and interaction with clay mineral particles have improved substantially over the past several years, triggered by pressing social issues such as management of huge amounts of waste soil accumulated after the Fukushima Dai-ichi nuclear power plant accident. In particular, computer-based molecular modeling supported by advanced hardware and algorithms has proven to be a powerful approach. Its application can now generally encompass the full complexity of clay particle adsorption sites from basal surfaces to interlayers with inserted water molecules, to edges including fresh and weathered frayed ones. On the other hand, its methodological schemes are now varied from traditional force-field molecular dynamics on large-scale realizations composed of many thousands of atoms including water molecules to first-principles methods on smaller models in rather exacting fashion. In this article, we overview new understanding enabled by simulations across methodological variations, focusing on recent insights that connect with experimental observations, namely: 1) the energy scale for cesium adsorption on the basal surface, 2) progress in understanding the structure of clay edges, which is difficult to probe experimentally, 3) cesium adsorption properties at hydrated interlayer sites, 4) the importance of the size relationship between the ionic radius of cesium and the interlayer distance at frayed edge sites, 5) the migration of cesium into deep interlayer sites, and 6) the effects of nuclear decay of radiocesium. Key experimental observations that motivate these simulation advances are also summarized. Furthermore, some directions toward future solutions of waste soil management are discussed based on the obtained microscopic insights.

Enhancing Thermal Boundary Conductance of Graphite-Metal Interface by Triazine-Based Molecular Bonding

Thermal boundary conductance between graphite and metal plays an important role in developing thermally conductive composites and contacts for thermal management. On the basis of the premise that the thermal boundary conductance (TBC) correlates with interfacial bonding strength, we conducted triazine-based molecular-bonding process to improve interfacial adhesion forces between a-axis of highly oriented pyrolytic graphite and aluminum. The surface coverage of molecular bonding at the interface is estimated by the X-ray photoelectron spectroscopy and thermal boundary conductance is measured by the time-domain thermoreflectance method. It is found that the TBC is directly proportional to the surface coverage of covalently bonded triazine linkers, with the proportionality constant for their increment rates being about unity. The experimental finding is supported by the corresponding simulation using the atomic Green's function method, which exhibits the same linear dependence on the surface coverage.

Thermal Boundary Characteristics of Homo-/Heterogeneous Interfaces

The interface of two solids in contact introduces a thermal boundary resistance (TBR), which is challenging to measure from experiments. Besides, if the interface is reactive, it can form an intermediate recrystallized or amorphous region, and extra influencing phenomena are introduced. Reactive force field Molecular Dynamics (ReaxFF MD) is used to study these interfacial phenomena at the (non-)reactive interface. The non-reactive interfaces are compared using a phenomenological theory (PT), predicting the temperature discontinuity at the interface. By connecting ReaxFF MD and PT we confirm a continuous temperature profile for the homogeneous non-reactive interface and a temperature jump in case of the heterogeneous non-reactive interface. ReaxFF MD is further used to understand the effect of chemical activity of two solids in contact. The selected Si/SiO 2 materials showed that the TBR of the reacted interface is two times larger than the non-reactive, going from 1 . 65 × 10 - 9 to 3 . 38 × 10 - 9 m 2 K/W. This is linked to the formation of an intermediate amorphous layer induced by heating, which remains stable when the system is cooled again. This provides the possibility to design multi-layered structures with a desired TBR.

Mechanical size effects of amorphous polymer-derived ceramics at the nanoscale: experiments and ReaxFF simulations

Here we report an unprecedented mechanical size effect at the nanoscale in polymer-derived ceramic (PDC) nanofibers. Silicon oxycarbide (SiOC) PDCs were fabricated as micro- and nanofibers without the aid of fillers. By decreasing the size of SiOC ceramic fibers from 1.1 μm to 630 nm (reduction of 74%), the strength of nanofibers nearly tripled, going from ∼1 GPa to ∼3.3 GPa. This increase in strength exceeds the predictions of the Griffith theorem, which relies on the length-scale dependence of energy release rate during crack propagation, suggesting a reduction in flaw size more than proportional to sample size. Given the crosslinked and amorphous nature of SiOC PDCs, flaws are likely microcracks and voids, which form during polymer degassing as it is pyrolyzed to PDC nanofibers. A reduction in sample size may favor degassing via diffusion, preceding bubble and void formation. We developed a new reactive force field (ReaxFF) with parameters for Si/O/C/H/N to study the mechanics of PDCs in extreme cases where no void is present. The models and experiments compare favorably in terms of the elastic modulus. The simulations suggest a strength of ∼8.5 GPa for a "flawless" structure, which is in line with extrapolated experimental results, with C-C breakage as the root cause of failure. This work clearly shows the benefits of utilizing nanoscale components as building blocks of superstrong PDC structures.

Radiocesium interaction with clay minerals: Theory and simulation advances Post-Fukushima

Insights at the microscopic level of the process of radiocesium adsorption and interaction with clay mineral particles have improved substantially over the past several years, triggered by pressing social issues such as management of huge amounts of waste soil accumulated after the Fukushima Dai-ichi nuclear power plant accident. In particular, computer-based molecular modeling supported by advanced hardware and algorithms has proven to be a powerful approach. Its application can now generally encompass the full complexity of clay particle adsorption sites from basal surfaces to interlayers with inserted water molecules, to edges including fresh and weathered frayed ones. On the other hand, its methodological schemes are now varied from traditional force-field molecular dynamics on large-scale realizations composed of many thousands of atoms including water molecules to first-principles methods on smaller models in rather exacting fashion. In this article, we overview new understanding enabled by simulations across methodological variations, focusing on recent insights that connect with experimental observations, namely: 1) the energy scale for cesium adsorption on the basal surface, 2) progress in understanding the structure of clay edges, which is difficult to probe experimentally, 3) cesium adsorption properties at hydrated interlayer sites, 4) the importance of the size relationship between the ionic radius of cesium and the interlayer distance at frayed edge sites, 5) the migration of cesium into deep interlayer sites, and 6) the effects of nuclear decay of radiocesium. Key experimental observations that motivate these simulation advances are also summarized. Furthermore, some directions toward future solutions of waste soil management are discussed based on the obtained microscopic insights.

Role of Interfaces in Elasticity and Failure of Clay-Organic Nanocomposites: Toughening upon Interface Weakening?

Synthetic organic-inorganic composites constitute a new class of engineering materials finding applications in an increasing range of fields. The interface between the constituting phases plays a pivotal role in the enhancement of mechanical properties. In exfoliated clay-organic nanocomposites, individual, high aspect ratio clay sheets are dispersed in the organic matrix providing large interfaces and hence efficient stress transfer. In this study, we aim at elucidating molecular-scale reinforcing mechanisms in a series of model clay-organic composite systems by means of reactive molecular simulations. In our models, two possible locations of failure initiation are present: one is the interlayer space of the clay platelet, and the other one is the clay-organic interface. We systematically modify the cohesiveness of the interface and assess how the failure mechanism changes when the different model composites are subjected to a tensile test. Besides a change in the failure mechanism, an increase in the released energy at the interface (meaning an increased overall toughness) are observed upon weakening the interface by bond removal. We propose a theoretical analysis of these results by considering a cohesive law that captures the effect of the interface on the composite mechanics. We suggest an atomistic interpretation of this cohesive law, in particular, how it relates to the degree of bonding at the interface. In a broader perspective, this work sheds light on the importance of the orthogonal behavior of interfaces to nanocomposites.

Influence of Calcium on the Initial Stages of the Sol-Gel Synthesis of Bioactive Glasses

Understanding how calcium interacts with silica sources and influences their polycondensation in aqueous solutions is of central importance for the development of more effective biomaterials by sol-gel approaches. For this purpose, the atomic-scale evolutions of a calcium-containing precursor solution corresponding to a typical sol-gel bioactive glass and of a corresponding Ca-free solution were compared using reactive molecular dynamics simulations. The simulations highlight a significantly faster rate of condensation when calcium is present in the initial solution, resulting in the formation of large and ramified silica clusters within 5 ns, which are absent in the Ca-free system. This different behavior has been analyzed and interpreted in terms of the Ca-induced nanosegregation in calcium-rich and silica-rich regions, which promotes the condensation reactions within the latter. By identifying a possible mechanism behind the limited incorporation of calcium in the silica nanoclusters formed in the early stages of the sol-gel process, these results could guide further studies aimed at identifying favorable experimental conditions to enhance initial calcium incorporation and thus produce sol-gel biomaterials with improved properties.

Insight on Tricalcium Silicate Hydration and Dissolution Mechanism from Molecular Simulations

Hydration of mineral surfaces, a critical process for many technological applications, encompasses multiple coupled chemical reactions and topological changes, challenging both experimental characterization and computational modeling. In this work, we used reactive force field simulations to understand the surface properties, hydration, and dissolution of a model mineral, tricalcium silicate. We show that the computed static quantities, i.e., surface energies and water adsorption energies, do not provide useful insight into predict mineral hydration because they do not account for major structural changes at the interface when dynamic effects are included. Upon hydration, hydrogen atoms from dissociated water molecules penetrate into the crystal, forming a disordered calcium silicate hydrate layer that is similar for most of the surfaces despite wide-ranging static properties. Furthermore, the dynamic picture of hydration reveals the hidden role of surface topology, which can lead to unexpected water tessellation that stabilizes the surface against dissolution.

Tribochemical mechanism of amorphous silica asperities in aqueous environment: a reactive molecular dynamics study

Reactive molecular dynamics (ReaxFF) simulations are used to explore the atomic-level tribochemical mechanism of amorphous silica (a-SiO2) in a nanoscale, single-asperity contact in an aqueous environment. These sliding simulations are performed in both a phosphoric acid solution and in pure water under different normal pressures. The results show that tribochemical processes have profound consequences on tribological performance. Water molecules could help avoid direct adhesive interaction between a-SiO2 surfaces in pure water under low normal load. However, formation and rupture of interfacial siloxane bonds are obviously observed under higher normal load. In phosphoric acid solution, polymerization of phosphoric acid molecules occurs, yielding oligomers under lower load, and tribochemical reactions between the molecules and the sliding surfaces could enhance wear under higher load. The bridging oxygen atoms in silica play an important role in the formation of interfacial covalent bonds, and hydrogen is found to have a weakening effect on these bonds, resulting in the rupture during shear-related loading. This work sheds light on tribochemical reactions as a mechanism for lubrication and wear in water-based or other tribological systems.

ReaxFF molecular dynamics simulations of CO collisions on an O-preadsorbed silica surface

A quasiclassical trajectory dynamics study was performed for carbon monoxide collisions over an oxygen preadsorbed β-cristobalite (001) surface. A reactive molecular force field (ReaxFF) was used to model the potential energy surface. The collisions were performed fixing several initial conditions: CO rovibrational states (v = 0-5 and j = 0, 20, 35), collision energies (0.05 ≤ E(col) ≤ 2.5 eV), incident angles (θ(v) = 0°, 45°) and surface temperatures (T(surf) = 300 K, 900 K). The principal elementary processes were the molecular reflection and the non-dissociative molecular adsorption. CO₂ molecules were also formed in minor extension via an Eley-Rideal reaction although some of them were finally retained on the surface. The scattered CO molecules tend to be translationally colder and internally hotter (rotationally and vibrationally) than the initial ones. The present study supports that CO + O(ad) reaction should be less important than O + O(ad) reaction over silica for similar initial conditions of reactants, in agreement with experimental data.