Polymorphism and phase transitions in Na2U2O7 from density functional perturbation theory.

Polymorphism and phase transitions in sodium diuranate, Na2U2O7, are investigated with density functional perturbation theory (DFPT). Thermal properties of crystalline α-, ß- and γ-Na2U2O7 polymorphs are predicted from DFPT phonon calculations, i.e., the first time for the high-temperature γ-Na2U2O7 phase (R3̄m symmetry). The standard molar isochoric heat capacities predicted within the quasi-harmonic approximation are for P21/a α-Na2U2O7 and C2/m ß-Na2U2O7, respectively. Gibbs free energy calculations reveal that α-Na2U2O7 (P21/a) and ß-Na2U2O7 (C2/m) are almost energetically degenerate at low temperature, with ß-Na2U2O7 becoming slightly more stable than α-Na2U2O7 as temperature increases. These findings are consistent with XRD data showing a mixture of α and ß phases after cooling of γ-Na2U2O7 to room temperature and the observation of a sluggish α → ß phase transition above ca. 600 K. A recently observed α-Na2U2O7 structure with P21 symmetry is also shown to be metastable at low temperature. Based on Gibbs free energy, no direct ß â†’ γ solid-solid phase transition is predicted at high temperature, although some experiments reported the existence of such phase transition around 1348 K. This, along with recent experiments, suggests the occurrence of a multi-step process consisting of initial ß-phase decomposition, followed by recrystallization into γ-phase as temperature increases.

Control of Structural Hydrophobicity and Cation Solvation on Interlayer Water Transport during Clay Dehydration.

Swelling clay hydration/dehydration is important to many environmental and industrial processes. Experimental studies usually probe equilibrium hydration states in an averaged manner and thus cannot capture the fast water transport and structural change in interlayers during hydration/dehydration. Using molecular simulations and thermogravimetric analyses, we observe a two-stage dehydration process. The first stage is controlled by evaporation at the edges: water molecules near hydrophobic sites and the first few water molecules of the hydration shell of cations move fast to particle edges for evaporation. The second stage is controlled by slow desorption of the last 1-2 water molecules from the cations and slow transport through the interlayers. The two-stage dehydration is strongly coupled with interlayer collapse and the coordination number changes of cations, all of which depend on layer charge distribution. This mechanistic interpretation of clay dehydration can be key to the coupled chemomechanical behavior in natural/engineered barriers.

Uranium carbonate complexes demonstrate drastic decrease in stability at elevated temperatures.

Quantitative understanding of uranium transport by high temperature fluids is crucial for confident assessment of its migration in a number of natural and artificially induced contexts, such as hydrothermal uranium ore deposits and nuclear waste stored in geological repositories. An additional recent and atypical context would be the seawater inundated fuel of the Fukushima Daiichi Nuclear Power Plant. Given its wide applicability, understanding uranium transport will be useful regardless of whether nuclear power finds increased or decreased adoption in the future. The amount of uranium that can be carried by geofluids is enhanced by the formation of complexes with inorganic ligands. Carbonate has long been touted as a critical transporting ligand for uranium in both ore deposit and waste repository contexts. However, this paradigm has only been supported by experiments conducted at ambient conditions. We have experimentally evaluated the ability of carbonate-bearing fluids to dissolve (and therefore transport) uranium at high temperature, and discovered that in fact, at temperatures above 100 °C, carbonate becomes almost completely irrelevant as a transporting ligand. This demands a re-evaluation of a number of hydrothermal uranium transport models, as carbonate can no longer be considered key to the formation of uranium ore deposits or as an enabler of uranium transport from nuclear waste repositories at elevated temperatures.

Structure-thermodynamics relationship of schoepite from first-principles.

The relationship between the structure and thermodynamic properties of schoepite, an important uranyl phase with formula [(UO2)8O2(OH)12]·12H2O formed upon corrosion of UO2, has been investigated within the framework of density functional perturbation theory (DFPT). Experimental crystallographic lattice parameters are well reproduced in this study using standard DFT. Phonon calculations within the quasi-harmonic approximation predict standard molar entropy and isobaric heat capacity of S0 = 179.60 J mol-1 K-1 and C0P = 157.4 J mol-1 K-1 at 298.15 K, i.e., ∼6% and ∼4% larger than existing DFPT-D2 calculations. The computed variation of the standard molar isobaric heat capacity with water content from schoepite (UO3·xH2O, x = 2.25) to dehydrated schoepite (x = 1) is predicted to be essentially linear along isotherms ranging from 100 to 500 K. These findings have important implications for the dehydration of layered uranyl corrosion phases and hygroscopic materials.

Gibbs Energy Minimization Model for Solvent Extraction with Application to Rare-Earths Recovery.

The emergence of technologies in which rare-earth elements provide critical functionality has increased the demand for these materials, with important implications for supply security. Recycling provides an option for mitigating supply risk and for creating economic value from the resale of recovered materials. While solvent extraction is a proven technology for rare-earth recovery and separation, its application often requires extensive trial-and-error experimentation to estimate parameter values and determine experimental design configurations. We describe a modeling strategy based on Gibbs energy minimization that incorporates parameter estimation for required thermodynamic properties as well as process design for solvent extraction and illustrate its applicability to rare earths separation. Visualization analysis during parameter estimation revealed a linear relationship between the standard enthalpies of the extractant and respective organo-metal complexes, analogous to the additivity principle for predicting molar volumes of organic compounds. Establishing this relationship reduced the size of the parameter estimation problem and yielded good agreement between model predictions and reported equilibrium extraction data, validating the property estimates for the organic phase species. Design exploration and optimization results map the space of feasible solvent extraction column configurations and identify the set of optimal design parameter values that meet recovery and purity targets.

Nonlinear dynamics and instability of aqueous dissolution of silicate glasses and minerals.

Aqueous dissolution of silicate glasses and minerals plays a critical role in global biogeochemical cycles and climate evolution. The reactivity of these materials is also important to numerous engineering applications including nuclear waste disposal. The dissolution process has long been considered to be controlled by a leached surface layer in which cations in the silicate framework are gradually leached out and replaced by protons from the solution. This view has recently been challenged by observations of extremely sharp corrosion fronts and oscillatory zonings in altered rims of the materials, suggesting that corrosion of these materials may proceed directly through congruent dissolution followed by secondary mineral precipitation. Here we show that complex silicate material dissolution behaviors can emerge from a simple positive feedback between dissolution-induced cation release and cation-enhanced dissolution kinetics. This self-accelerating mechanism enables a systematic prediction of the occurrence of sharp dissolution fronts (vs. leached surface layers), oscillatory dissolution behaviors and multiple stages of glass dissolution (in particular the alteration resumption at a late stage of a corrosion process). Our work provides a new perspective for predicting long-term silicate weathering rates in actual geochemical systems and developing durable silicate materials for various engineering applications.

Relationship between crystal structure and thermo-mechanical properties of kaolinite clay: beyond standard density functional theory.

The structural, mechanical and thermodynamic properties of 1 : 1 layered dioctahedral kaolinite clay, with ideal Al2Si2O5(OH)4 stoichiometry, were investigated using density functional theory corrected for dispersion interactions (DFT-D2). The bulk moduli of 56.2 and 56.0 GPa predicted at 298 K using the Vinet and Birch-Murnaghan equations of state, respectively, are in good agreement with the recent experimental value of 59.7 GPa reported for well-crystallized samples. The isobaric heat capacity computed for uniaxial deformation of kaolinite along the stacking direction reproduces calorimetric data within 0.7-3.0% from room temperature up to its thermal stability limit.

On the role of strong electron correlations in the surface properties and chemistry of uranium dioxide.

We report density functional calculations of the surface properties and chemistry of UO(2)(111) performed within the generalized gradient approximation corrected with an effective Hubbard parameter (GGA + U within Dudarev's formalism) to account for the strong on-site Coulomb repulsion between U 5f electrons. The variation of the properties of periodic slab models, with collinear ferromagnetic and antiferromagnetic arrangements of the uranium magnetic moments, was investigated while ramping up the effective Hubbard parameter from U(eff) = 0 eV, corresponding to standard density functional theory, up to U(eff) = 4 eV, the value that correctly reproduces the antiferromagnetic ground state of bulk UO(2). The chemical interactions of molecular water, dissociated water, dissociated oxygen and co-adsorbed molecular water and monatomic oxygen with the UO(2)(111) surface were also studied as functions of the U(eff) parameter. Calculations reveal that some of the key electronic and chemical properties controlling the surface reactivity are very sensitive to the value of this strong electron correlation parameter.

Ganglion dynamics and its implications to geologic carbon dioxide storage.

Capillary trapping of a nonwetting fluid phase in the subsurface has been considered as an important mechanism for geologic storage of carbon dioxide (CO(2)). This mechanism can potentially relax stringent requirements for the integrity of cap rocks for CO(2) storage and therefore can significantly enhance storage capacity and security. We here apply ganglion dynamics to understand the capillary trapping of supercritical CO(2) (scCO(2)) under relevant reservoir conditions. We show that, by breaking the injected scCO(2) into small disconnected ganglia, the efficiency of capillary trapping can be greatly enhanced, because the mobility of a ganglion is inversely dependent on its size. Supercritical CO(2) ganglia can be engineered by promoting CO(2)-water interface instability during immiscible displacement, and their size distribution can be controlled by injection mode (e.g., water-alternating-gas) and rate. We also show that a large mobile ganglion can potentially break into smaller ganglia due to CO(2)-brine interface instability during buoyant rise, thus becoming less mobile. The mobility of scCO(2) in the subsurface is therefore self-limited. Vertical structural heterogeneity within a reservoir can inhibit the buoyant rise of scCO(2) ganglia. The dynamics of scCO(2) ganglia described here provides a new perspective for the security and monitoring of subsurface CO(2) storage.

Structures of uranyl peroxide hydrates: a first-principles study of studtite and metastudtite.

The structures of the only known minerals containing peroxide, namely studtite [(UO(2))O(2)(H(2)O)(4)] and metastudtite [(UO(2))O(2)(H(2)O)(2)], have been investigated using density functional theory. The structure of metastudtite crystallizing in the orthorhombic space group Pnma (Z = 4) is reported for the first time at the atomic level and the computed lattice parameters, a = 8.45, b = 8.72, c = 6.75 Å, demonstrate that the unit cell of metastudtite is larger than previously reported dimensions (Z = 2) derived from experimental X-ray powder diffraction data.