Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,387 records · Page 77

Hydrogen Diffusion in Slit Pores: Role of Temperature, Pressure, Confinement, and Roughness

Diffusion of hydrogen (H 2 ) is important to understand the leakage risk and transport behavior for H 2 geologic storage. We applied molecular dynamics simulations to investigate the influencing factors of H2 diffusion in the slit pores of calcite, hematite, and quartz, owing to their abundance. It is revealed that the H2 self-diffusion coefficient increases with the temperature, regardless of the type of pore minerals. The diffusion of H 2 in the 20 nm slit pores falls into the bulk diffusion regime when the pressure is 10 MPa. The self-diffusion of H 2 decreases with pressure in all three types of slit pores, following a power law model with the exponents ranging from -0.825 to -0.964. Furthermore, the impact of confinement on H 2 diffusion is more pronounced for the slit pores with stronger interactions with H 2 -like calcite. The role of surface roughness in H 2 diffusion depends on the slit aperture. The rough surface enhances H 2 diffusion in the larger slit pores due to the enlarged effective pore space, whereas it weakens H 2 diffusion in the small slit pores due to stronger adsorption. These findings will fill the knowledge gap on the coupling effect of different factors influencing H 2 diffusion.

08 HYDROGEN↗

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗

Real Time 3D Observations of Portland Cement Carbonation at CO2 Storage Conditions

Depleted oil reservoirs are considered a viable solution to the global challenge of CO 2 storage. A key concern is whether the wells can be suitably sealed with cement to hinder the escape of CO 2 . Under reservoir conditions, CO 2 is in its supercritical state, and the high pressures and temperatures involved make real-time microscopic observations of cement degradation experimentally challenging. Here, we present an in situ 3D dynamic X-ray micro computed tomography (μ-CT) study of well cement carbonation at realistic reservoir stress, pore-pressure, and temperature conditions. The high-resolution time-lapse 3D images allow monitoring the progress of reaction fronts in Portland cement, including density changes, sample deformation, and mineral precipitation and dissolution. By switching between flow and nonflow conditions of CO 2 -saturated water through cement, we were able to delineate regimes dominated by calcium carbonate precipitation and dissolution. For the first time, we demonstrate experimentally the impact of the flow history on CO 2 leakage risk for cement plugging. In-situ μ-CT experiments combined with geochemical modeling provide unique insight into the interactions between CO 2 and cement, potentially helping in assessing the risks of CO 2 storage in geological reservoirs.

58 GEOSCIENCES↗

Mobility of Radionuclides in Fractured Carbonate Rocks: Lessons from a Field-Scale Transport Experiment

Current research on radionuclide disposal is mostly conducted in granite, clay, saltstone, or volcanic tuff formations. These rock types are not always available to host a geological repository in every nuclear waste-generating country, but carbonate rocks may serve as a potential alternative. To assess their feasibility, a forced gradient cross-borehole tracer experiment was conducted in a saturated fractured chalk formation. The mobility of stable Sr and Cs (as analogs for their radioactive counterparts), Ce (an actinide analog), Re (a Tc analog), bentonite particles, and fluorescent dye tracers through the flow path was analyzed. The migration of each of these radionuclide analogs (RAs) was shown to be dependent upon their chemical speciation in solution, their interactions with bentonite, and their sorption potential to the chalk rock matrix. The brackish groundwater resulted in flocculation and immobilization of most particulate RAs. Nevertheless, the high permeability of the fracture system allowed for fast overall transport times of all aqueous RAs investigated. This study suggests that the geochemical properties of carbonate rocks may provide suitable conditions for certain types of radionuclide storage (in particular, brackish, high-porosity, and low-permeability chalks). Nevertheless, careful consideration should be given to high-permeability fracture networks that may result in high radionuclide mobility.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Quantification of CO 2 Mineralization at the Wallula Basalt Pilot Project

In 2013, the Pacific Northwest National Laboratory led a geologic carbon sequestration field demonstration where ~1000 tonnes of CO 2 was injected into several deep Columbia River Basalt zones near Wallula, Washington. Rock core samples extracted from the injection zone two years after CO 2 injection revealed nascent carbonate mineralization that was qualitatively consistent with expectations from laboratory experiments and reactive transport modeling. In this study, we report on a new detailed analysis of the 2012 pre-injection and 2015 post-injection hydrologic tests that capitalizes on the difference in fluid properties between scCO 2 and water to assess changes in near-field, wellbore, and reservoir conditions that are apparent approximately two years following the end of injection. This comparative hydrologic test analysis method provides a new way to quantify the amount of injected CO 2 that was mineralized in the field test. Modeling results indicate that approximately 60% of the injected CO 2 was sequestered via mineralization within two years, with the resulting carbonates occupying ~4% of the available reservoir pore space. The method presented here provides a new monitoring tool to assess the fate of CO 2 injected into chemically reactive basalt formations but could also be adapted for long-term monitoring and verification within more traditional subsurface carbon storage reservoirs.

54 ENVIRONMENTAL SCIENCES↗

Climate vs Energy Security: Quantifying the Trade-offs of BECCS Deployment and Overcoming Opportunity Costs on Set-Aside Land

Bioenergy with carbon capture and storage (BECCS) sits at the nexus of the climate and energy security. We evaluated trade-offs between scenarios that support climate stabilization (negative emissions and net climate benefit) or energy security (ethanol production). Our spatially explicit model indicates that the foregone climate benefit from abandoned cropland (opportunity cost) increased carbon emissions per unit of energy produced by 14–36%, making geologic carbon capture and storage necessary to achieve negative emissions from any given energy crop. The toll of opportunity costs on the climate benefit of BECCS from set-aside land was offset through the spatial allocation of crops based on their individual biophysical constraints. Dedicated energy crops consistently outperformed mixed grasslands. We estimate that BECCS allocation to land enrolled in the Conservation Reserve Program (CRP) could capture up to 9 Tg C year –1 from the atmosphere, deliver up to 16 Tg CE year –1 in emissions savings, and meet up to 10% of the US energy statutory targets, but contributions varied substantially as the priority shifted from climate stabilization to energy provision. Furthermore, our results indicate a significant potential to integrate energy security targets into sustainable pathways to climate stabilization but underpin the trade-offs of divergent policy-driven agendas.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Imprint of Anthropogenic Sources and Soil Removal on the Surface Concentration of H2 in the Contiguous US

Hydrogen (H2) is experiencing renewed interest throughout the world as a low carbon fuel alternative or complement to fossil fuels. Significant uncertainties remain regarding the environmental impact of increasing H2 usage, in part due to gaps in our understanding of the H2 atmospheric budget, including the H2 release from industrial activities and the H2 soil removal, the most important sink of H2. This study focuses on H2 dry air mole fractions measured by the NOAA Global Monitoring Laboratory in discrete ambient air samples collected every few days at sites located in the contiguous United States between 2010 and 2022. We take advantage of the long-term observations from this network to study the regional distribution of H2 sources using the potential source contribution function (PSCF). We find that H2 PSCF is consistent with a large anthropogenic source of atmospheric H2 from the incomplete combustion of fossil fuel. However, episodic air masses with high H2 dry air mole fractions (>700 ppb) recorded at some sampling locations in the Western and Southeastern US are not associated with elevated CO, suggesting significant noncombustion anthropogenic or geological sources of H2. Air masses depleted in H2 are recorded in rural or remote continental sites. Our analysis of observations from a tall tower in NE Colorado suggests that the degree of H2 depletion not only reflects differences in soil exposure but also the regional inhibition of the soil sink under arid conditions.

Paulot, Fabien↗

Locating Undocumented Wells Using Historical Oil and Gas Exploration Maps: A Case Study in Osage County, Oklahoma

Undocumented oil and gas wells lack reliable information about their locations and characteristics, making them difficult to identify. These wells can result in unanticipated delays and costs in the development of nearby surface and subsurface resources, and, if improperly plugged, can cause contamination. This study leverages historical petroleum exploration maps to locate such wells, focusing on Osage County, Oklahoma. Two sets of early 20th century oil and gas exploration maps by the United States Geological Survey were georeferenced and analyzed using a computer vision model to detect well symbols. The locations of detected wells were compared to the location of known wells in the database from the Bureau of Indian Affairs Osage Agency to identify potential undocumented wells. The analysis yielded over 500 potential undocumented wells, with dry holes constituting the largest fraction. Field verification confirmed the presence of some undocumented wells. Comparison with prior work revealed limited overlap, underscoring the complementary value of historical oil and gas maps for locating undocumented wells. This approach demonstrates the utility of integrating historical cartographic resources with modern geospatial and machine learning techniques to improve the identification and management of undocumented wells.

Energy - Petroleum↗

Sustainable Critical Minerals: Mining, Production, and Legacy

Critical minerals (CMs) are required in renewable energy systems, electric vehicles, electronics, and national defense. Global demand has increased rapidly due to energy and digital transitions, and the United States has become increasingly reliant on vulnerable supply chains. Geological scarcity, trade barriers, environmental hurdles, and processing challenges impede domestic CM production growth. Emerging strategies, including circular economy practices and advanced sustainable mining and processing techniques, can mitigate supply chain and economic risks. To advance these approaches, consistent research priorities, policy revisions, and technology investments are required to develop a resilient CM supply. Doing so with a focus on sustainability goals can simultaneously advance domestic economic and technological goals while limiting environmental impacts and supporting global environmental goals.

biomining↗

Nanoscale Mg-Depleted Layers Slow Carbonation of Forsterite (Mg 2 SiO 4 ) When Water Is Limited

Passivation of silicate surfaces by accumulated reaction products is an obstacle to efficient CO 2 mineralization. In this study, we investigate a unique passivation effect during the carbonation of the basalt mineral forsterite (Mg 2 SiO 4 ) in humid supercritical CO 2 (50 °C, 90 bar). Using in situ high-pressure infrared spectroscopy, we demonstrate that dissolution of forsterite into a thin water film slows significantly after reaction for ˜24 h, even under far-from-equilibrium conditions. Further, 29 Si magic angle spinning nuclear magnetic resonance spectroscopy detects a highly polymerized amorphous silica at this stage. On the basis of transmission electron microscopy and energy dispersive X-ray spectroscopy, we show that the silica is present as a Mg-depleted layer that is just 2–3 nm thick on the reacted forsterite particles. The decrease in the level of forsterite dissolution in the presence of an extraordinarily thin Mg-depleted layer can be strongly linked to properties of the thin fluid film at the surface, highlighting the importance of water during mineral carbonation. This study furthers our understanding of silicate mineral carbonation under select low-water, humidified fluid conditions relevant to basaltic geologic reservoirs, recovery of critical elements by carbonation of mafic ores, and sequestration of atmospheric CO 2 by enhanced rock weathering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of Mass Transfer and Sorption in CO 2 /Brine/Rock Systems via In Situ FT-IR

CO 2 geological storage in deep saline formations is considered a promising method to mitigate anthropogenic CO 2 emissions and, thereby, minimize changes to the Earth’s atmosphere. A fundamental understanding of CO 2 mass transfer and sorption phenomena in brine-saturated reservoir formations is necessary to understand the long-term fate of injected CO 2 as it is subjected to different (physical, dissolution, and mineral) trapping mechanisms. In this work, we investigate CO 2 sorption in brine-saturated and dry Mt. Simon sandstone samples via in-situ Fourier Transform infrared spectroscopy (FT-IR) at elevated pressures, ranging from 0.3 MPa to 8.3 MPa, at a temperature of 50 ⁰C. The FT-IR spectra of bulk-phase CO 2 were simultaneously recorded under the same conditions. For bulk-phase CO 2 , we observed, in agreement with past studies, a doublet peak at 2361 cm -1 and 2336 cm -1 and another peak (ν 2 bending mode) at 667 cm -1 . With increasing pressure, the position of the peak at 667 cm -1 remains invariant, however, when crossing into the supercritical region the doublet peak degenerates onto a single peak at 2336 cm -1 with a barely visible shoulder at 2361 cm -1 . The bulk CO 2 data provide a perfect fit for Beer’s law for the whole range of pressure conditions. For the dry sample, the IR spectrum is experimentally indistinguishable from the bulk CO 2 spectrum, signifying that if physical adsorption occurs to any significant extent, the adsorbed CO 2 molecules are not substantially more rotationally constrained than the dense bulk CO 2 molecules. For the brine-saturated sample, we observe a strong band centered at 2342 cm -1 and a small companion peak at 2360-2361 cm -1 that degenerates into a barely visible shoulder peak at the higher pressures. The 2342 cm -1 band has been previously observed by other investigators for CO 2 dissolved in bulk water/brine as well during its adsorption on a variety of other wet natural porous media. We observe no peaks corresponding to bicarbonate or carbonate bulk species, which correlates well with the prior literature on similar low-pH aqueous solutions. The integrated peak area for the CO 2 sorbed in the brine-saturated sample correlates linearly with its solubility in the same bulk brine, as measured separately via a PVT-cell approach. Here, this validates the accuracy of both techniques, and the potential of the FT-IR method to be used in the study of mass transfer and adsorption in such systems. To that effect, a simple mathematical model is presented to analyze the FT-IR data to determine the CO 2 effective diffusivity in the brine-saturated sandstone sample.

58 GEOSCIENCES↗

Hyperstoichiometric Uranium Dioxides: Rapid Synthesis and Irradiation-Induced Structural Changes

Uranium dioxide (UO 2 ), the primary fuel for commercial nuclear reactors, incorporates excess oxygen forming a series of hyperstoichiometric oxides. Thin layers of these oxides, such as UO 2.12 , form readily on the fuel surface and influence its properties, performance, and potentially geologic disposal. This work reports a rapid and straightforward combustion process in uranyl nitrate–glycine–water solutions to prepare UO 2.12 nanomaterials and thin films. We also report on the investigation of the structural changes induced in the material by irradiation. Despite the simple processing aspects, the combustion synthesis of UO 2.12 has a sophisticated chemical mechanism involving several exothermic steps. Raman spectroscopy and single-crystal X-ray diffraction (XRD) measurements reveal the formation of a complex compound containing the uranyl moiety, glycine, H 2 O, and NO 3 – groups in reactive solutions and dried combustion precursors. Combustion diagnostic methods, gas-phase mass spectroscopy, differential scanning calorimetry (DSC), and extracted activation energies from DSC measurements show that the rate-limiting step of the process is the reaction of ammonia with nitrogen oxides formed from the decomposition of glycine and uranyl nitrate, respectively. However, the exothermic decomposition of the complex compound determines the maximum temperature of the process. In situ transmission electron microscopy (TEM) imaging and electron diffraction measurements show that the decomposition of the complex compound directly produces UO 2 . The incorporation of oxygen at the cooling stage of the combustion process is responsible for the formation of UO 2.12 . Spin coating of the solutions and brief annealing at 670 K allow the deposition of uniform films of UO 2.12 with thicknesses up to 300 nm on an aluminum substrate. Irradiation of films with Ar 2+ ions (1.7 MeV energy, a fluence of up to 1 × 10 17 ions/cm 2 ) shows unusual defect-simulated grain growth and enhanced chemical mixing of UO 2.12 with the substrate due to the high uranium ion diffusion in films. As a result, the method described in this work allows the preparation of actinide oxide targets for fundamental nuclear science research and studies associated with stockpile stewardship.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Cation Short-Range Ordering of MgAl 2 O 4 and NiAl 2 O 4 Spinel Oxides at High Temperatures via In Situ Neutron Total Scattering

Complex oxides that adopt the isometric spinel structure (AB 2 O 4 ) are important for numerous technological applications and are relevant for certain geological processes, which involve exposure to extreme environments such as high pressures and temperatures. Recent studies have shown that the changes to the spinel structure caused by these environments are complex and depend on the material length scale under consideration. In this study we have expanded this approach to the behavior of spinels under high temperature. In situ neutron total scattering experiments, coupled with pair distribution function analysis, performed on two spinel compositions with various levels of pre-existing disorder (MgAl 2 O 4 and NiAl 2 O 4 ) revealed that both compositions trend to a state of maximum disorder where the A and B cations are randomly distributed amongst the two available sites. Temperature-induced cation inversion, conventionally understood as exchange of cations on the A and B sites, is locally expressed as an atomic rearrangement to a tetragonal symmetry, a correlation that is retained up to the maximum temperature studied (1000 °C). Furthermore, a complex thermal expansion behavior is revealed wherein the oxide materials expand heterogeneously at the level of coordination polyhedral with an apparent dependence on bond strength.

36 MATERIALS SCIENCE↗

Understanding the Hydronium Cation in the Solid State: A Study in Synthetic Hydronium Uranyl Phosphate and Arsenate Mineral Systems and Their Irradiation Stability

Here, we present neutron diffraction, electronic structure calculations, and optical vibrational spectroscopic characterization of hydronium uranyl phosphate (HUP), the synthetic analog of chernikovite ((H 3 O)(UO 2 )(PO 4 )·3H 2 O) to gain insight into structural features of the hydronium cation in the solid state. HUP crystallizes in P2 1 /c instead of previously reported P4/ncc 1 at room temperature and Pccn below 302 K. Lower symmetry was required due to hydronium and interstitial water positions within the structure. Crystallographic positions of hydronium cations determined from powder neutron diffraction data are bolstered by Raman and attenuated total reflectance infrared spectroscopic measurements and are further informed by density functional theory with phonon eigenvector analysis for spectral assignments. Finally, HUP and its arsenate analog (trogerite, (H 3 O)(UO 2 )(AsO 4 )·3H 2 O), (HUAs) were studied using He 2+ irradiation as an analog for α radiolysis to investigate the irradiation stability of these phases and the stability of hydronium cations in the solid state. Dose studies were employed wherein each sample was irradiated to 5, 10, 15, 25, and 50 MGy. Structural insight regarding irradiated materials is gained using optical vibrational spectroscopy and powder X-ray diffraction. Surprising irradiation stability of HUAs was found up to 50 MGy of dose, which may have important implications for understanding and modeling the geologic stability of legacy U.

actinides↗

NWPEsSe: an Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au_55, ligated Au_8^(2+), Au_8 supported on graphene oxide and defected rutile, and a large cluster assembly ?[Co?_6 Te_8 (PEt_3 )_6][C_60 ]_n, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems. The work described in this publication was performed at Pacific Northwest National Laboratory (PNNL), which is operated by Battelle for the United States Department of Energy (DOE) under Contract DE-AC05-76RL0180. J. Z. and V.-A. G. acknowledge support from DOE, Office of Science, Office of Basic Energy Sci-ences, Chemical, Geological and Biological Sciences Division and computing resources from PNNL’s Research Computing Facility and the National Energy Research Scientific Computing Center.

Zhang, Jun↗

Density Functional Tight-Binding Simulations Reveal the Presence of Surface Defects on the Quartz (101)–Water Interface

Understanding the structure and reactivity of quartz–water interfaces is critical for numerous applications in the geological, environmental, and biological sciences. However, disagreements on the atomic-level structure of the interfaces between experiments and simulations are hampering our ability to predict the surface reactivity. Here, we used density functional tight-binding (DFTB)-based molecular dynamics to simulate a series of quartz (101) surfaces having different types and densities of surface defects in water and compared them with the structures determined by X-ray reflectivity measurements. The DFTB simulations are able to reproduce previous classical and quantum mechanical predictions of the pristine quartz (101)–water interface that disagree with experimental observations. To remedy this situation, a set of defective quartz surfaces having various surface silicon (Si) vacancies were built as indicated by recent experimental studies. We found that the rotation of surface [SiO 4 ] tetrahedra near Si vacancies can lead to outward displacements of Si atoms similar to those observed in the experiments. The presence of additional surface Si vacancies caused inward relaxations of terminal oxygens through the formation of hydrogen bonds. Finally, the overall results indicate that the quartz (101)–water interface may include a mixture of geminal (≡Si–(OH) 2 )- and vicinal (≡Si–OH)-type silanol groups together with the presence of surface Si vacancies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Dynamics Simulations of Calcite Fracture in Water

Calcite (CaCO 3 ) is one of the most common minerals in geologic and engineered systems. It is often in contact with aqueous solutions, causing chemically assisted fracture that is critical to understanding the stability of subsurface systems and manmade structures. Calcite fracture was evaluated with reactive molecular dynamics simulations, including the impacts of crack tip geometry (notch), the presence of water, and surface hydroxyl groups. Chemo-mechanical weakening was assessed by comparing the loads where fracture began to propagate. Our analyses show that in the presence of a notch, the load at which crack growth begins is lower, compared to the effect of water or surface hydroxyls. Additionally, the breaking of two adjacent Ca–O bonds is the kinetic limitation for crack initiation, since transiently broken bonds can reform, not resulting in crack growth. In aqueous environments, fresh (not hydroxylated) calcite surfaces exhibited water strengthening. Manual addition of H + and/or OH – species on the (104) calcite surface resulted in chemo-mechanical weakening of calcite by 9%. Achieving full hydroxylation of the calcite surface was thermodynamically and kinetically limited, with only 0.17–0.01 OH/nm 2 surface hydroxylation observed on the (104) surface at the end of the simulations. In conclusion, the limited reactivity of pure water with the calcite surface restricts the chemo-mechanical effects and suggests that reactions between physiosorbed water and localized structural defects may be dominating the chemo-mechanical process in the studies where water weakening has been reported.

Calcite↗

Electronic and Vibrational Contributions to the Bulk Stabilities of Trivalent 3d Transition Metal Oxyhydroxides from Electronic Structure Calculations

Any understanding of the thermodynamic driving forces active in the global cycles of earth abundant metals requires attendant knowledge of the bulk mineral free energies. This can be particularly challenging for metal (oxyhydr)oxide phases which at low temperature aqueous conditions can be produced at a variety of hydration states and particle sizes. Here we present high quality density functional theory calculations of the relative bulk stabilities of a geologically important family of polymorphic trivalent metal oxyhydroxide crystal structures MOOH (M=Cr, Al, Fe). To improve understanding of their relative stabilities we permute the metal cations into the various possible canonical structure types including goethite (isostructural to diaspore and bracewillite), lepidocrocite (isosctructural to boehmite), guyaniite, and grimaldiite. Electronic, vibrational, and magnetic contributions are considered on a consistent theoretical footing using hybrid functionals and converged plane-wave bases. We find the considered bulk structures to be thermodynamically degenerate within the precision of the available theory, despite the fact that the vibrational and magnetic contributions vary substantially among the cations. A significant conclusion is that Cr3+ supports strong hydrogen bonds, low O-H stretching frequencies and facile H movement relative to the other cations.

Prange, Micah P.↗