Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments.

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO3·nH2O systems, where M2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

Thin Water Films Enable Low-Temperature Magnesite Growth Under Conditions Relevant to Geologic Carbon Sequestration.

Injecting supercritical CO2 (scCO2) into basalt formations for long-term storage is a promising strategy for mitigating CO2 emissions. Mineral carbonation can result in permanent entrapment of CO2; however, carbonation kinetics in thin H2O films in humidified scCO2 is not well understood. We investigated forsterite (Mg2SiO4) carbonation to magnesite (MgCO3) via amorphous magnesium carbonate (AMC; MgCO3·xH2O, 0.5 < x < 1), with the goal to establish the fundamental controls on magnesite growth rates at low H2O activity and temperature. Experiments were conducted at 25, 40, and 50 °C in 90 bar CO2 with a H2O film thickness on forsterite that averaged 1.78 ± 0.05 monolayers. In situ infrared spectroscopy was used to monitor forsterite dissolution and the growth of AMC, magnesite, and amorphous SiO2 as a function of time. Geochemical kinetic modeling showed that magnesite was supersaturated by 2 to 3 orders of magnitude and grew according to a zero-order rate law. The results indicate that the main drivers for magnesite growth are sustained high supersaturation coupled with low H2O activity, a combination of thermodynamic conditions not attainable in bulk aqueous solution. This improved understanding of reaction kinetics can inform subsurface reactive transport models for better predictions of CO2 fate and transport.

Critical Water Coverage during Forsterite Carbonation in Thin Water Films: Activating Dissolution and Mass Transport.

In geologic carbon sequestration, CO2 is injected into geologic reservoirs as a supercritical fluid (scCO2). The carbonation of divalent silicates exposed to humidified scCO2 occurs in angstroms to nanometers thick adsorbed H2O films. A threshold H2O film thickness is required for carbonate precipitation, but a mechanistic understanding is lacking. In this study, we investigated carbonation of forsterite (Mg2SiO4) in humidified scCO2 (50 °C and 90 bar), which serves as a model system for understanding subsurface divalent silicate carbonation reactivity. Attenuated total reflection infrared spectroscopy pinpointed that magnesium carbonate precipitation begins at 1.5 monolayers of adsorbed H2O. At about this same H2O coverage, transmission infrared spectroscopy showed that forsterite dissolution begins and electrical impedance spectroscopy demonstrated that diffusive transport accelerates. Molecular dynamics simulations indicated that the onset of diffusion is due to an abrupt decrease in the free-energy barriers for lateral mobility of outer-spherically adsorbed Mg2+. The dissolution and mass transport controls on divalent silicate reactivity in wet scCO2 could be advantageous for maximizing permeability near the wellbore and minimize leakage through the caprock.

Immobilizing Pertechnetate in Ettringite via Sulfate Substitution.

Technetium-99 immobilization in low-temperature nuclear waste forms often relies on additives that reduce environmentally mobile pertechnetate (TcO4-) to insoluble Tc(IV) species. However, this is a short-lived solution unless reducing conditions are maintained over the hazardous life cycle of radioactive wastes (some ∼10,000 years). Considering recent experimental observations, this work explores how rapid formation of ettringite [Ca6Al2(SO4)3(OH)12·26(H2O)], a common mineral formed in cementitious waste forms, may be used to directly immobilize TcO4-. Results from ab initio molecular dynamics (AIMD) simulations and solid-phase characterization techniques, including synchrotron X-ray absorption, fluorescence, and diffraction methods, support successful incorporation of TcO4- into the ettringite crystal structure via sulfate substitution when synthesized by aqueous precipitation methods. One sulfate and one water are replaced with one TcO4- and one OH- during substitution, where Ca2+-coordinated water near the substitution site is deprotonated to form OH- for charge compensation upon TcO4- substitution. Furthermore, AIMD calculations support favorable TcO4- substitution at the SO42- site in ettringite rather than gypsum (CaSO4·2H2O, formed as a secondary mineral phase) by at least 0.76 eV at 298 K. These results are the first of their kind to suggest that ettringite may contribute to TcO4- immobilization and the overall lifetime performance of cementitious waste forms.

Iron Vacancies Accommodate Uranyl Incorporation into Hematite.

Radiotoxic uranium contamination in natural systems and nuclear waste containment can be sequestered by incorporation into naturally abundant iron (oxyhydr)oxides such as hematite (α-Fe2O3) during mineral growth. The stability and properties of the resulting uranium-doped material are impacted by the local coordination environment of incorporated uranium. While measurements of uranium coordination in hematite have been attempted using extended X-ray absorption fine structure (EXAFS) analysis, traditional shell-by-shell EXAFS fitting yields ambiguous results. We used hybrid functional ab initio molecular dynamics (AIMD) simulations for various defect configurations to generate synthetic EXAFS spectra which were combined with adsorbed uranyl spectra to fit experimental U L3-edge EXAFS for U6+-doped hematite. We discovered that the hematite crystal structure accommodates a trans-dioxo uranyl-like configuration for U6+ that substitutes for structural Fe3+, which requires two partially protonated Fe vacancies situated at opposing corner-sharing sites. Surprisingly, the best match to experiment included significant proportions of vacancy configurations other than the minimum-energy configuration, pointing to the importance of incorporation mechanisms and kinetics in determining the state of an impurity incorporated into a host phase under low temperature hydrothermal conditions.

Incorporation Modes of Iodate in Calcite.

Iodate (IO3-) incorporation in calcite (CaCO3) is a potential sequestration pathway for environmental remediation of radioiodine-contaminated sites (e.g., Hanford Site, WA), but the incorporation mechanisms have not been fully elucidated. Ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure spectroscopy (EXAFS) were combined to determine the local coordination environment of iodate in calcite, the associated charge compensation schemes (CCS), and any tendency for surface segregation. IO3- substituted for CO32- and charge compensation was achieved by substitution of Ca2+ by Na+ or H+. CCS that minimized the I-Na/H distance or placed IO3- at the surface were predicted by density functional theory to be energetically favored, with the exception of HIO3, which was found to be metastable relative to the formation of HCO3-. Iodine K-edge EXAFS spectra were calculated from AIMD trajectories and used to fit the experimental spectrum. The best-fit combination consisted of a significant proportion of surface-segregated IO3- and charge compensation was predominantly by H+. Important implications are therefore that pH should strongly affect the extent of IO3- incorporation and that IO3- accumulated at the surface of CaCO3 particles may undergo mobilization under conditions that promote calcite dissolution. These impacts need to be considered in calcite-based iodate remediation strategies.

Electronic response of aluminum-bearing minerals.

Aluminum-bearing minerals show different hydrogen evolution and dissolution properties when subjected to radiation, but the complicated sequence of events following interaction with high-energy radiation is not understood. To gain insight into the possible mechanisms of hydrogen production in nanoparticulate minerals, we study the electronic response and determine the bandgap energies of three common aluminum-bearing minerals with varying hydrogen content: gibbsite (Al(OH)3), boehmite (AlOOH), and alumina (Al2O3) using electron energy loss spectroscopy, X-ray photoelectron spectroscopy, and first-principles electronic structure calculations employing hybrid density functionals. We find that the amount of hydrogen has only a small effect on the number and spectrum of photoexcitations in this class of materials. Electronic structure calculations demonstrate that low energy electrons are isotropically mobile, while holes in the valence band are likely constrained to move in layers. Furthermore, holes in the valence band of boehmite are found to be significantly more mobile than those in gibbsite, suggesting that the differences in radiolytic and dissolution behavior are related to hole transport.

Trace Uranium Partitioning in a Multiphase Nano-FeOOH System.

The characterization of trace elements in minerals using extended X-ray absorption fine structure (EXAFS) spectroscopy constitutes a first step toward understanding how impurities and contaminants interact with the host phase and the environment. However, limitations to EXAFS interpretation complicate the analysis of trace concentrations of impurities that are distributed across multiple phases in a heterogeneous system. Ab initio molecular dynamics (AIMD)-informed EXAFS analysis was employed to investigate the immobilization of trace uranium associated with nanophase iron (oxyhydr)oxides, a model system for the geochemical sequestration of radiotoxic actinides. The reductive transformation of ferrihydrite [Fe(OH)3] to nanoparticulate iron oxyhydroxide minerals in the presence of uranyl (UO2)2+(aq) resulted in the preferential incorporation of U into goethite (α-FeOOH) over lepidocrocite (γ-FeOOH), even though reaction conditions favored the formation of excess lepidocrocite. This unexpected result is supported by atomically resolved transmission electron microscopy. We demonstrate how AIMD-informed EXAFS analysis lifts the strict statistical limitations and uncertainty of traditional shell-by-shell EXAFS fitting, enabling the detailed characterization of the local bonding environment, charge compensation mechanisms, and oxidation states of polyvalent impurities in complex multiphase mineral systems.

Ab Initio Molecular Dynamics of Uranium Incorporated in Goethite (α-FeOOH): Interpretation of X-ray Absorption Spectroscopy of Trace Polyvalent Metals.

Incorporation of economically or environmentally consequential polyvalent metals into iron (oxyhydr)oxides has applications in environmental chemistry, remediation, and materials science. A primary tool for characterizing the local coordination environment of such metals, and therefore building models to predict their behavior, is extended X-ray absorption fine structure spectroscopy (EXAFS). Accurate structural information can be lacking yet is required to constrain and inform data interpretation. In this regard, ab initio molecular dynamics (AIMD) was used to calculate the local coordination environment of minor amounts of U incorporated in the structure of goethite (α-FeOOH). U oxidation states (VI, V, and IV) and charge compensation schemes were varied. Simulated trajectories were used to calculate the U LIII-edge EXAFS function and fit experimental EXAFS data for U incorporated into goethite under reducing conditions. Calculations that closely matched the U EXAFS of the well-characterized mineral uraninite (UO2), and constrained the S02 parameter to be 0.909, validated the approach. The results for the U-goethite system indicated that U(V) substituted for structural Fe(III) in octahedral uranate coordination. Charge balance was achieved by the loss of one structural proton coupled to addition of one electron into the solid (-1 H+, +1 e-). The ability of AIMD to model higher energy states thermally accessible at room temperature is particularly relevant for protonated systems such as goethite, where proton transfers between adjacent octahedra had a dramatic effect on the calculated EXAFS. Vibrational effects as a function of temperature were also estimated using AIMD, allowing separate quantification of thermal and configurational disorder. In summary, coupling AIMD structural modeling and EXAFS experiments enables modeling of the redox behavior of polyvalent metals that are incorporated in conductive materials such as iron (oxyhydr)oxides, with applications over a broad swath of chemistry and materials science.

Internal Domains of Natural Porous Media Revealed: Critical Locations for Transport, Storage, and Chemical Reaction.

Internal pore domains exist within rocks, lithic fragments, subsurface sediments, and soil aggregates. These domains, termed internal domains in porous media (IDPM), represent a subset of a material's porosity, contain a significant fraction of their porosity as nanopores, dominate the reactive surface area of diverse media types, and are important locations for chemical reactivity and fluid storage. IDPM are key features controlling hydrocarbon release from shales in hydraulic fracture systems, organic matter decomposition in soil, weathering and soil formation, and contaminant behavior in the vadose zone and groundwater. Traditionally difficult to interrogate, advances in instrumentation and imaging methods are providing new insights on the physical structures and chemical attributes of IDPM, and their contributions to system behaviors. Here we discuss analytical methods to characterize IDPM, evaluate information on their size distributions, connectivity, and extended structures; determine whether they exhibit unique chemical reactivity; and assess the potential for their inclusion in reactive transport models. Ongoing developments in measurement technologies and sensitivity, and computer-assisted interpretation will improve understanding of these critical features in the future. Impactful research opportunities exist to advance understanding of IDPM, and to incorporate their effects in reactive transport models for improved environmental simulation and prediction.

Solvation structure and transport properties of alkali cations in dimethyl sulfoxide under exogenous static electric fields.

A combination of molecular dynamics simulations and pulsed field gradient nuclear magnetic resonance spectroscopy is used to investigate the role of exogenous electric fields on the solvation structure and dynamics of alkali ions in dimethyl sulfoxide (DMSO) and as a function of temperature. Good agreement was obtained, for select alkali ions in the absence of an electric field, between calculated and experimentally determined diffusion coefficients normalized to that of pure DMSO. Our results indicate that temperatures of up to 400 K and external electric fields of up to 1 V nm(-1) have minimal effects on the solvation structure of the smaller alkali cations (Li(+) and Na(+)) due to their relatively strong ion-solvent interactions, whereas the solvation structures of the larger alkali cations (K(+), Rb(+), and Cs(+)) are significantly affected. In addition, although the DMSO exchange dynamics in the first solvation shell differ markedly for the two groups, the drift velocities and mobilities are not significantly affected by the nature of the alkali ion. Overall, although exogenous electric fields induce a drift displacement, their presence does not significantly affect the random diffusive displacement of the alkali ions in DMSO. System temperature is found to have generally a stronger influence on dynamical properties, such as the DMSO exchange dynamics and the ion mobilities, than the presence of electric fields.

Molecular dynamics simulations of uranyl and uranyl carbonate adsorption at aluminosilicate surfaces.

Adsorption at mineral surfaces is a critical factor controlling the mobility of uranium(VI) in aqueous environments. Therefore, molecular dynamics (MD) simulations were performed to investigate uranyl(VI) adsorption onto two neutral aluminosilicate surfaces, namely, the orthoclase (001) surface and the octahedral aluminum sheet of the kaolinite (001) surface. Although uranyl preferentially adsorbs as a bidentate inner-sphere complex on both surfaces, the free energy of adsorption on the orthoclase surface (-15 kcal mol(-1)) is significantly more favorable than that at the kaolinite surface (-3 kcal mol(-1)), which is attributed to differences in surface functional groups and the ability of the orthoclase surface to release a surface potassium ion upon uranyl adsorption. The structures of the adsorbed complexes compare favorably with X-ray absorption spectroscopy results. Simulations of the adsorption of uranyl complexes with up to three carbonate ligands revealed that uranyl complexes coordinated to up to two carbonate ions are stable on the orthoclase surface whereas uranyl carbonate surface complexes are unfavored at the kaolinite surface. Combining the MD-derived equilibrium adsorption constants for orthoclase with aqueous equilibrium constants for uranyl carbonate species indicates the presence of adsorbed uranium complexes with one or two carbonates under alkaline conditions, in support of current uranium(VI) surface complexation models.

Structure, kinetics, and thermodynamics of the aqueous uranyl(VI) cation.

In this work, molecular simulation techniques were employed to gain insight into the structural, kinetic, and thermodynamic properties of the uranyl(VI) cation (UO2(2+)) in aqueous solution. The simulations made use of an atomistic potential model (force field) derived in this work and based on the model of Guilbaud and Wipff [J. Mol. Struct. (THEOCHEM) 1996 , 366 , 55 - 63]. Reactive flux and thermodynamic integration calculations show that the derived potential model yields predictions for the water exchange rate and free energy of hydration, respectively, that are in agreement with experimental data. The water binding energies, hydration shell structure, and self-diffusion coefficient were also calculated and analyzed. Finally, a combination of metadynamics and transition path sampling simulations was employed to probe the mechanisms of water exchange reactions in the first hydration shell of the uranyl ion. These atomistic simulations indicate, based on two-dimensional free energy surfaces, that water exchanges follow an associative interchange mechanism. The nature and structure of the water exchange transition states were also determined. The improved potential model is expected to lead to more accurate predictions of uranyl adsorption energies at mineral surfaces using potential-based molecular dynamics simulations.

Local density changes and carbonate rotation enable Ba incorporation in amorphous calcium carbonate.

Incorporation of a Ba impurity in amorphous calcium carbonate (ACC) is shown with ab initio molecular dynamics simulations to have a long-range effect on its atomic-level structure and to be energetically favoured relative to incorporation in crystalline calcium carbonate polymorphs. The ability of carbonate ions to rotate and of ACC to undergo local density changes explain ACC's propensity for incorporating divalent metal impurities with a wide range of ionic radii. These findings provide an atomic-level basis for understanding the significant effects of low concentrations of impurities on the structure of ACC.

Diffusion and adsorption of uranyl carbonate species in nanosized mineral fractures.

Atomistic simulations were performed to study the diffusion and adsorption of Ca(2)UO(2)(CO3)3 and of some of its constituent species, i.e., UO(2)2+, CO(3)2­, and UO(2)CO3, in feldspar nanosized fractures. Feldspar is important to uranium remediation efforts at the U.S. Department of Energy Hanford site as it has been found in recent studies to host contaminants within its intragrain fractures. In addition, uranyl carbonate species are known to dominate U(VI) speciation in conditions relevant to the Hanford site. Molecular dynamics (MD) simulations showed that the presence of the feldspar surface diminishes the diffusion coefficients of all of the species considered in this work and that the diffusion coefficients do not reach their bulk aqueous solution values in the center of a 2.5 nm fracture. Moreover, the MD simulations showed that the rate of decrease in the diffusion coefficients with decreasing distance from the surface is greater for larger adsorbing species. Free energy profiles of the same species adsorbing on the feldspar surface revealed a large favorable free energy of adsorption for UO(2)2+ and UO(2)CO3, which are able to adsorb to the surface with their uranium atom directly bonded to a surface hydroxyl oxygen, whereas adsorption of CO(3)2­ and Ca(2)UO(2)(CO3)3, which attach to the surface via hydrogen bonding from a surface hydroxyl group to a carbonate oxygen, was calculated to be either only slightly favorable or unfavorable.

Atomistic simulations of uranium incorporation into iron (hydr)oxides.

Atomistic simulations were carried out to characterize the coordination environments of U incorporated in three Fe-(hydr)oxide minerals: goethite, magnetite, and hematite. The simulations provided information on U-O and U-Fe distances, coordination numbers, and lattice distortion for U incorporated in different sites (e.g., unoccupied versus occupied sites, octahedral versus tetrahedral) as a function of the oxidation state of U and charge compensation mechanisms (i.e., deprotonation, vacancy formation, or reduction of Fe(III) to Fe(II)). For goethite, deprotonation of first shell hydroxyls enables substitution of U for Fe(III) with a minimal amount of lattice distortion, whereas substitution in unoccupied octahedral sites induced appreciable distortion to 7-fold coordination regardless of U oxidation states and charge compensation mechanisms. Importantly, U-Fe distances of ∼3.6 Šwere associated with structural incorporation of U and cannot be considered diagnostic of simple adsorption to goethite surfaces. For magnetite, the octahedral site accommodates U(V) or U(VI) with little lattice distortion. U substituted for Fe(III) in hematite maintained octahedral coordination in most cases. In general, comparison of the simulations with available experimental data provides further evidence for the structural incorporation of U in iron (hydr)oxide minerals.

Low temperature and limited water activity reveal a pathway to magnesite via amorphous magnesium carbonate.

Forsterite carbonated in thin H2O films to magnesite via amorphous magnesium carbonate during reaction with H2O-bearing liquid CO2 at 25 °C. This novel reaction pathway contrasts with previous studies that were carried out at higher H2O activity and temperature, where more highly hydrated nesquehonite was the metastable intermediate.

Bond-valence constraints on liquid water structure.

The recent controversy about the structure of liquid water pits a new model involving water molecules in relatively stable "rings-and-chains" structures against the standard model that posits water molecules in distorted tetrahedral coordination. Molecular dynamics (MD) simulations, both classical and ab initio, almost uniformly support the standard model, but because none of them can yet reproduce all of the anomalous properties of water, they leave room for doubt. We argue that it is possible to evaluate these simulations by testing them against their adherence to the bond-valence model, a well-known and quantitatively accurate empirical summary of the behavior of atoms in the bonded networks of inorganic solids. Here we use the results of ab initio MD simulations of ice, water, and several solvated aqueous species to show that the valence sum rule (the first axiom of the bond-valence model) is followed in both solid and liquid bond networks. We then test MD simulations of water, employing several popular potential models against this criterion and the experimental O-O RDF. It appears that most of those tested cannot satisfy both criteria well, except TIP4P, TIP4P/2005, and TIP5P. If the valence sum rule really can be applied to simulated liquid structures, then it follows that the bonding behaviors of atoms in liquids are in some ways identical to those in solids. We support this interpretation by showing that the simulations produce O-H...O geometries that are completely consistent with the range of geometries available in solids, and the distributions of instantaneous valence sums reaching the atoms in both the ice and liquid water simulations are essentially identical. Furthermore, we show that none of the extant asymmetric water potentials that produce "rings-and-chains" structures can satisfy our geometric criteria. Taken together, this is powerful evidence in favor of the standard distorted tetrahedral model of liquid water structure.

Transition path sampling of water exchange rates and mechanisms around aqueous ions.

The rates and mechanisms of water exchange around two aqueous ions, namely, Na(+) and Fe(2+), have been determined using transition path sampling. In particular, the pressure dependence of the water exchange rates was computed to determine activation volumes. A common approach for calculating water exchange rates, the reactive flux method, was also employed and the two methods were compared. The water exchange rate around Na(+) is fast enough to be calculated by direct molecular dynamics simulations, thus providing a reference for comparison. Both approaches predicted exchange rates and activation volumes in agreement with the direct simulation results. Four additional sodium potential models were considered to compare the results of this work with the only activation volume for Na(+) previously determined from molecular simulation [D. Spangberg et al., Chem. Phys. Lett. 276, 114 (1997)] and provide the best possible estimate of the activation volume based on the ability of the models to reproduce known properties of the aqueous sodium ion. The Spangberg and Hermansson [D. Spangberg and K. Hermansson, J. Chem. Phys. 120, 4829 (2004)] and X-Plor/Charmm-22 [M. Patra and M. Karttunen, J. Comput. Chem. 25, 678 (2004)] models performed best and predicted activation volumes of -0.22 and -0.78 cm(3) mol(-1), respectively. For water exchange around Fe(2+), transition path sampling predicts an activation volume of +3.8 cm(3) mol(-1), in excellent agreement with the available experimental data. The potential of mean force calculation in the reactive flux approach, however, failed to sufficiently sample appropriate transition pathways and the opposite pressure dependence of the rate was predicted as a result. Analysis of the reactive trajectories obtained with the transition path sampling approach suggests that the Fe(2+) exchange reaction takes place via an associative interchange mechanism, which goes against the conventional mechanistic interpretation of a positive activation volume. Collectively, considerable insight was obtained not only for the exchange rates and mechanisms for Na(+) and Fe(2+) but also for identifying the most robust modeling strategy for these purposes.

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.