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.

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.

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.

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.

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.

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.

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.

Dynamics of self-reorganization explains passivation of silicate glasses.

Understanding the dissolution of silicate glasses and minerals from atomic to macroscopic levels is a challenge with major implications in geoscience and industry. One of the main uncertainties limiting the development of predictive models lies in the formation of an amorphous surface layer--called gel--that can in some circumstances control the reactivity of the buried interface. Here, we report experimental and simulation results deciphering the mechanisms by which the gel becomes passivating. The study conducted on a six-oxide borosilicate glass shows that gel reorganization involving high exchange rate of oxygen and low exchange rate of silicon is the key mechanism accounting for extremely low apparent water diffusivity (∼10-21 m2 s-1), which could be rate-limiting for the overall reaction. These findings could be used to improve kinetic models, and inspire the development of new molecular sieve materials with tailored properties as well as highly durable glass for application in extreme environments.

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.

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.

Free-Energy Landscape of the Dissolution of Gibbsite at High pH.

The individual elementary reactions involved in the dissolution of a solid into solution remain mostly speculative due to a lack of direct experimental probes. In this regard, we have applied atomistic simulations to map the free-energy landscape of the dissolution of gibbsite from a step edge as a model of metal hydroxide dissolution. The overall reaction combines kink formation and kink propagation. Two individual reactions were found to be rate-limiting for kink formation, that is, the displacement of Al from a step site to a ledge adatom site and its detachment from ledge/terrace adatom sites into the solution. As a result, a pool of mobile and labile adsorbed species, or adatoms, exists before the release of Al into solution. Because of the quasi-hexagonal symmetry of gibbsite, kink site propagation can occur in multiple directions. Overall, our results will enable the development of microscopic mechanistic models of metal oxide dissolution.

Boehmite and Gibbsite Nanoplates for the Synthesis of Advanced Alumina Products.

Boehmite (γ-AlOOH) and gibbsite (α-Al-(OH)3) are important archetype (oxy)hydroxides of aluminum in nature that also play diverse roles across a plethora of industrial applications. Developing the ability to understand and predict the properties and characteristics of these materials, on the basis of their natural growth or synthesis pathways, is an important fundamental science enterprise with wide-ranging impacts. The present study describes bulk and surface characteristics of these novel materials in comprehensive detail, using a collectively sophisticated set of experimental capabilities, including a range of conventional laboratory solids analyses and national user facility analyses such as synchrotron X-ray absorption and scattering spectroscopies as well as small-angle neutron scattering. Their thermal stability is investigated using in situ temperature-dependent Raman spectroscopy. These pure and effectively defect-free materials are ideal for synthesis of advanced alumina products.

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.

Effects of Oxygen-Containing Functional Groups on Supercapacitor Performance.

Molecular dynamics (MD) simulations of the interface between graphene and the ionic liquid 1-butyl-3-methylimidazolium trifluoromethanesulfonate (BMIM OTf) were carried out to gain molecular-level insights into the performance of graphene-based supercapacitors and, in particular, determine the effects of the presence of oxygen-containing defects at the graphene surface on their integral capacitance. The MD simulations predict that increasing the surface coverage of hydroxyl groups negatively affects the integral capacitance, whereas the effect of the presence of epoxy groups is much less significant. The calculated variations in capacitance are found to be directly correlated to the interfacial structure. Indeed, hydrogen bonding between hydroxyl groups and SO3 moieties prevents BMIM(+) and OTf(-) ions from interacting favorably in the interfacial layer and restrains the orientation and mobility of OTf(-) ions, thereby reducing the interfacial permittivity of the ionic liquid. The results of the simulations can facilitate the rational design of electrode materials for supercapacitors.