Search NASASearch

SEARCH · Search NASA

Results for “characterization of rare earth elements (REE)”

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.

3,232 records · Page 18

Numerical Solution of Fluid Flow and Heat Tranfer Problems with Surface Radiation

This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

S Ahuja

Biocatalyst discovery and design for plastics deconstruction: A multi‐scale perspective

Plastic waste accumulation poses significant environmental challenges due to a lack of economical solutions for the molecular deconstruction of diverse synthetic polymers. Biological‐based degradation offers promise but is hindered by the crystallinity, hydrophobicity, and additive complexity of plastics, which restrict biocatalyst access and activity. To address these problems, we propose a multi‐scale framework that combines detailed materials characterization, optimization of plastic‐biomolecular interfacial interactions, and enhancement of biocatalytic kinetics to develop effective plastic‐deconstructing enzymes. This approach leverages principles from reaction kinetics, transport and interfacial phenomena, and enzyme engineering to systematically address barriers across diverse plastic types. Our framework aims to accelerate the discovery and optimization of biocatalysts capable of scalable, selective, and efficient deconstruction of plastic waste. These advances hold potential to enable sustainable biological recycling and upcycling pathways, contributing to global efforts in mitigating plastic pollution and promoting circular material economies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Optical Vibrational Spectroscopic Investigation of Natural and Synthetic Analogs of Uranyl Oxyhydroxyhydrate Minerals

Uranyl oxy-hydroxy-hydrate minerals are common alteration products of uraninite (UO2+x) which is chemically and structurally analogous to uranium dioxide nuclear fuel. Therefore, structural and spectroscopic investigations of these alteration minerals and their analogous anthropogenic counterparts can provide insight into the environmental behavior of nuclear fuel cycle materials. Previously, we compiled available vibrational spectroscopic data for uranyl minerals in the Compendium of Uranium Raman and Infrared Experimental Spectra (CURIES) and found that only 37% of known uranyl oxy-hydroxy-hydrate minerals had spectra readily available in the literature and existing databases for inclusion therein. Furthermore, no available infrared spectra for this mineral group were included in CURIES. To expand our understanding of the spectroscopic features of uranyl oxy-hydroxy-hydrates and the structural origins thereof, we collected, and now include in CURIES, Raman and infrared spectra for an additional 12 uranyl hydroxide phases. To better understand the impact of structural and compositional variations of these phases on their spectroscopic features, we compare Raman spectra of different anion sheet topological groups and of analogous phases hosting different counter cations. We identify spectroscopic variations related to differences in equatorial bonding and structural changes as a result of cation substitution. We also prepare a uranyl hydroxide phase via hydrolysis of uranyl fluoride (UO2F2) as an analog of hydrolysis reactions that occur in nuclear fuel cycle materials; and we find that the alteration product of UO2F2, despite chemical and structural similarities to uranyl oxy-hydroxy-hydrate minerals, is readily distinguishable from related mineral phases using Raman spectroscopy. In this work, we provide new insights into the structural origins of spectroscopic features in uranyl oxy-hydroxy-hydrate minerals, improve the average Raman spectrum for this group of minerals, and thereby improve capabilities for identifying these mineral species and related anthropogenic phases using Raman spectroscopy.

Barth, Brodie [ORNL] (ORCID:0000000256142601)

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT

Life cycle test results of a bipolar nickel hydrogen battery

A history is given of low Earth orbit (LEO) laboratory test data on a 6.5 ampere-hour bipolar nickel hydrogen battery designed and built at the NASA Lewis Research Center. The bipolar concept is a means of achieving the goal of producing an acceptable battery, of higher energy density, able to withstand the demands of low-Earth-orbit regimes. Over 4100 LEO cycles were established on a ten cell battery. It seems that any perturbation on normal cycling effects the cells performance. Explanations and theories of the battery's behavior are varied and widespread among those closely associated with it. Deep discharging does provide a reconditioning effect and further experimentation is planned in this area. The battery watt-hour efficiency is about 75 percent and the time averaged, discharge voltage is about 1.26 volts for all cells at both the C/4 and LEO rate. Since a significant portion of the electrode capacity has degraded, the LEO cycle discharges are approaching depths of 90 to 100 percent of the high rate capacity. Therefore, the low end-of-discharge voltages occur precipitously after the knee of the discharge curve and is more an indication of electrode capacity and is a lesser indicator of overall cell performance.

Cataldo, R. L.

Lunar and Planetary Science XXXIV

The 34th Lunar and Planetary Science Conference was held March 17-21, 2003. Topics included planetary exploration, crater research on Mars, Earth, Moon, and other planets or satellites, imaging techniques and image analysis, age determination, albedo studies, petrographic studies, isotope composition studies, instrument design, sampling methods, landform analysis, asteroids, impact analysis, impact melts, and related research.

Lunar and Planetary Institute

Influence of Powder Characteristics and Processing Methods on Creep Performance of Powder Metallurgy Hot Isostatic Pressed SS316

The U.S. nuclear energy expansion goals are driving the demand for manufacturing routes that can rapidly produce large, complex, near net shape components. Powder metallurgy hot isostatic pressing (PM HIP) is an advanced manufacturing technique that can be economically scaled-up, while alleviating the supply chain challenges that forging and casting face in terms of cost and lead time constraints. This makes PM-HIP a viable technology to aid and accelerate large-scale part manufacturing for nuclear applications. However, large-scale qualification and deployment of this technology require a thorough understanding of the influence of powder feedstock quality, powder handling history, hot isostatic pressing (HIP) parameters, and subsequent heat treatment on microstructural evolution and elevated temperature mechanical performance. The present work focusses on 316 austenitic stainless steel (SS316) which is most commonly used in high temperature environments for nuclear applications Results from this study show that PM HIPed SS316 meets ASME tensile requirements at room temperature and at elevated temperature. However, creep performance of PM-HIPed SS316 remains inferior to its wrought counterpart, demonstrating that tensile performance alone is not a reliable metric for long duration high temperature integrity. Further, this report delineates powder derived microstructural features that govern creep damage, with key evidences pointing to “microstructural inheritance” from gas atomized powder feedstocks. Commercial SS316 powders of varying chemical compositions and recycling histories were studies, and the results showed large differences in elemental segregation, oxide surface layers and secondary phase distributions. Multi-scale characterization revealed segregation of chromium, molybdenum, manganese and silicon at the boundaries and the precipitation of manganese-, silicon-, and molybdenum-oxides. During HIP consolidation, these surface oxides transform into decorated prior particle boundaries (PPBs) and grain boundary inclusions that persist through conventional post-HIP solution annealing treatment. The retained oxides in post-HIP microstructures were found to influence grain growth, precipitation behavior, and ultimately creep cavitation and fracture. Such post-HIP heat treatments are therefore limited by a complex trade-off between grain growth, and oxide coarsening which aggravate creep damage by acting as nucleation sites for cavities. The objective of this work is to establish an integrated processing–structure–property framework for PM-HIP 316 stainless steel by investigating the influence of powder feedstock characteristics in pre- and post-HIP processing as well as to understand the significance of post-HIP heat treatment on microstructural evolution and creep properties. The results from this report emphasize the significance of powder feedstock integrity in improving creep performance of PM-HIPed SS316, by highlight the effect of rapid solidification, elemental segregation, oxide formation and powder recycling on microstructural defect inheritance following HIP consolidation. Rather than considering HIP processing, solution annealing, and mechanical performance independently, this report treats powder production, HIP consolidation, post-HIP thermal processing, and creep deformation as interconnected stages within a continuous metallurgical process. The resulting framework will provide a scientific basis for developing feedstock engineering strategies capable of improving long-term reliability of PM-HIP stainless steels and accelerating their qualification for advanced nuclear applications.

Ajjarapu, Pavan [Oak Ridge National Laboratory (OR

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING

The Synthesis and Ring-Opening Metathesis Polymerization of Energetic Norbornene Materials

Energetic norbornenes are promising candidates toward the development of new energetic polymers due to the synthetic versatility of norbornene ring-opening metathesis polymerizations used in commercial applications. We report the synthesis of two energetic norbornene materials that can be made in two steps with modest yields, containing either trinitroethanol or fluorodinitroethanol substituents. The norbornene monomers were then polymerized, and the polymers were characterized by Fourier transform infrared spectroscopy (FT-IR), differential scanning calorimetry (DSC), contact angle measurements, and proton and carbon nuclear magnetic resonance spectroscopies ( 1 H and 13 C{ 1 H} NMR). Additionally, small-scale safety data consisting of electrostatic discharge (ESD), friction (FS), and impact (IS) sensitivities were measured for the norbornene monomers and their resulting polymers. These analyses revealed that energetic norbornene materials are relatively insensitive and have densities comparable to that of TNT (1.47–1.81 g·cm –3 ).

36 MATERIALS SCIENCE

Detection and characterization of detached tidal dwarf galaxies

Tidal interactions between galaxies often give rise to tidal tails, which can harbor concentrations of stars and interstellar gas resembling dwarf galaxies. Some of these tidal dwarf galaxies (TDGs) have the potential to detach from their parent galaxies and become independent entities, but their long-term survival is uncertain. In this study, we conducted a search for detached TDGs associated with a sample of 39 interacting galaxy pairs in the local Universe using infrared, ultraviolet, and optical images. We employed IR colors and UV/optical/IR spectral energy distributions to identify potential interlopers, such as foreground stars or background quasars. Through spectroscopic observations using the Boller and Chivens spectrograph at San Pedro Mártir Observatory, we confirmed that six candidate TDGs are at the same redshift as their putative parent galaxy pairs. We identified and measured emission lines in the optical spectra and calculated nebular oxygen abundances, which range from log(O/H) = 8.10 ± 0.01 to 8.51 ± 0.02. We have serendipitously discovered an additional detached TDG candidate in Arp72 using available spectra from SDSS. Utilizing the photometric data and the CIGALE code for stellar population and dust emission fitting, we derived the stellar masses, stellar population ages, and stellar metallicities for these detached TDGs. Compared to standard mass-metallicity relations for dwarf galaxies, five of the seven candidates have higher than expected metallicities, confirming their tidal origins. One of the seven candidates remains unclear due to large uncertainties in metallicity, and another has stellar and nebular metallicities compatible with those of a preexisting dwarf galaxy. The latter object is relatively compact in the optical relative to its stellar mass, in contrast to the other candidate TDGs, which have large diameters for their stellar masses compared to most dwarf galaxies. The derived stellar population ages range from 100 Myr to 900 Myr, while the inferred stellar masses are between 2 × 10 6 M ⊙ and 8 × 10 7 M ⊙ . Four of the six TDGs are associated with the gas-rich M51-like pair Arp 72, one TDG is associated with a second M51-like pair Arp 86, and another is associated with Arp 65, an approximately equal mass pair. In spite of the relatively low stellar masses of these TDGs, they have survived for at least 100–900 Myrs, suggesting that they are stable and in dynamical equilibrium. We conclude that encounters with a relatively low-mass companion (1/10th–1/4th of the mass of the primary) can also produce long-lasting TDGs.

Astronomy & Astrophysics

Multikilowatt Bipolar Nickel/Hydrogen Battery

High energy densities appear feasible. Nickel/hydrogen battery utilizing bipolar construction in common pressure vessel, addressing needs for multikilowatt storage for low-Earth-orbit applications, designed and 10-cell prototype model tested. Modular-concept-design 35-kW battery projected energy densities of 20 to 24 Wh/b (160 to 190 kj/kg) and 700 to 900 Wh/ft3 (90 to 110 MJ/m3) and incorporated significant improvements over state-of-the-art storage systems.

Source record

Characterization of ceramic fuel powder packing fractions to support INFLUX

Triply periodic minimal surface (TPMS)-based structures show marked potential in novel nuclear reactor fuel designs, as their high surface area-to-volume ratio increases the efficiency of heat transfer out of the fuel, enabling safer, more innovative reactor designs. This milestone report addresses the role of dUO 2 powder processing route on the fill behavior of TPMS-based cladding shells to understand and advance the feasibility of manufacturing TPMS-based nuclear fuel forms. dUO 2 powder was processed through either a dry granulation route, varying consolidation pressure, or through milling, varying milling time, milling method and milled size distribution. The lowest tapped bulk densities (TBD), but best powder flowabilities, were obtained when testing unprocessed dUO 2 powder which was prone to self-agglomeration and formed low-density spheroids. The highest TBD and lowest flowabilities were obtained when using powder produced by hammer-milling dUO 2 powder to pass through a 200-mesh sieve, which led to particles with angular morphologies. Powder produced by dry granulation exhibited TBD that varied according to the consolidation pressure used to form the initial pellets and exhibited improved flowabilities when compared to hammer-milled material. Because of the large span of granule sizes formed as well as the irregular shape associated with the granules, a packing fraction of 0.69 was achieved, exceeding the analytical solution for random close packing of mono-sized spheres. TPMS polymer shells were loaded with unprocessed, granulated, and hammer-milled dUO 2 powders, and their qualitative packing behaviors were analyzed using x-ray computed tomography (xCT). TBDs calculated after loading TPMS polymer shells were 10-20% lower when compared to tapped bulk density measurements taken in a glass graduated cylinder, indicating a non-trivial impact on the tapped bulk density of either the TPMS channel size, TPMS channel surface material, powder cohesiveness, or a combination of the two parameters. A metallic zircaloy-4 TPMS shell will be loaded with hammer-milled dUO 2 powder upon receipt of the shell from Oak Ridge National Laboratory (ORNL) and shipped to Idaho National Labs (INL) for subsequent hot isostatic pressing (HIP) densification experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Technology Readiness Assessment for Stirling Power Convertor

Stirling-based power conversion system could represent a critical enabling technology for future NASA missions, offering significantly enhanced power efficiency over traditional thermoelectric systems. To ensure the successful infusion of this technology into deep space and planetary surface missions, a rigorous evaluation of its maturity and associated risks is essential. The objective of this Technology Readiness Assessment (TRA) is to determine the current technology readiness level (TRL) of the Stirling power convertor and identify key technical risks and challenges associated with its maturation and potential infusion into a targeted flight system. The evaluation process will rigorously follow the NASA TRA procedure: (1) determining the requirements that underpin the design, function, and performance; (2) finding and listing all new technology elements (NTEs) and critical technology elements (CTEs); (3) identifying the level of integration assumed for the assessment; and (4) assessing the TRL for each NTE and CTE. The overall subassembly TRL will then be determined via a roll-up using the “weakest link” criterion, culminating in a final (5) assessment of risks to further progression of maturity. This paper will present the results of a comprehensive technology readiness assessment (TRA) conducted on the Stirling-based dynamic power conversion subassembly, excluding the radioisotope heat source. The paper will also describe a detailed breakdown of the TRL roll-up and a matrix of identified risks and recommended mitigation paths, along with identified key technical risks and challenges to achieve the required long-duration performance. Key Words: Radioisotope Power Conversion, Technology Readiness Assessment

Radioisotope Power Conversion

Molecular Inventory and Comparative Organic Profiling of A Homogenized Aggregate Sample From Asteroid (101955) Bennu

Evaluating the distribution and diversity of solvent-soluble organic molecules in extraterrestrial materials provide information about the physicochemical environments where prebiotic chemistry occurred and is essential for assessing the inventory of prebiotic compounds available to the early Earth. In this work, we conducted a comprehensive molecular characterization of solvent-soluble organic compounds in a homogenized aggregate sample (ID: OREX-800107-128) of unsorted regolith from asteroid (101955) Bennu returned by NASA’s OSIRIS-REx mission. Using two-dimensional gas chromatography coupled with high-resolution time-of-flight mass spectrometry (GC×GC-HRMS), we identified 87 organic compounds across several molecular classes, including alkanes, polycyclic aromatic hydrocarbons, and nitrogen-, oxygen-, and sulfur-bearing species. Comparison with the organic inventory of CI (Orgueil), C2-ungrouped (Tagish Lake and Tarda), and CM2 (Murchison) carbonaceous chondrites revealed that Bennu sample OREX-800107-128 is most similar to Tagish Lake and Tarda yet compositionally distinct from all of these meteorites. This Bennu sample exhibits a relatively high abundance of low-molecular-weight aromatics (e.g., toluene), pyridine derivatives, and sulfur species such as dimethyl sulfite, consistent with episodic low-temperature aqueous alteration. The presence of these organic compounds suggests that this material experienced chemically diverse processes in the parent body conducive to complex organic synthesis, including aldol condensation reactions responsible for N- and O-bearing compounds. We also compared the molecular distributions of this Bennu sample and previously analyzed samples from (162173) Ryugu to assess potential common origins between these two carbonaceous asteroids. Bennu shows more chemically diverse and nitrogen-rich organic materials than Ryugu. While both asteroids share broadly similar primitive compositions, the observed chemical differences indicate distinct evolutionary histories within a likely common formation environment though whether Bennu and Ryugu derive from a single compositionally heterogeneous parent body or from multiple closely related progenitors remains unresolved. Further detailed comparative analyses of material from Bennu, Ryugu, and specific carbonaceous meteorites (e.g., Tagish Lake and Tarda) will provide new insights into how the N-rich molecular profile observed in Bennu samples was formed and evolved in the outer Solar System.

sample return

Structural and spectroscopic characterization of a Tb(IV) polyoxometalate

There has been a renaissance in high valent lanthanide chemistry, which has resulted in the first examples of tetravalent praseodymium (Pr) and terbium (Tb) in molecular systems. These feats have been achieved with tailored ligands that facilitate tetravalent lanthanide stability in non-aqueous conditions. The next step in realizing the potential of high valent lanthanide chemistry is moving towards complexes that are stable in ambient and aqueous conditions. In this investigation, we advance this paradigm by obtaining definitive evidence of Tb(IV) in a molecular system under aqueous conditions utilizing the lacunary Wells-Dawson polyoxometalate, K 10 P 2 W 17 O 61 •20H 2 O. Herein we present the single crystal structure of K 16 Tb(IV)(P 2 W 17 O 61 ) 2 •39.85H 2 O (Tb(IV)W 34 ) as well as extensive spectroscopic evidence obtained via UV–Vis-NIR, X-ray absorption near edge spectroscopy, continuous wave X-band electron paramagnetic resonance spectroscopy and SQUID magnetometry measurements to confirm the oxidation state of Tb in W 34 complexes as +4.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Taking three-dimensional x-ray diffraction (3DXRD) from the synchrotron to the laboratory scale

Three-dimensional x-ray diffraction (3DXRD), a rotating x-ray diffraction technique, is a powerful tool for studying the micromechanical behavior of polycrystalline materials, capable of measuring the volume, position, orientation, and strain of thousands of grains simultaneously. However, its application has been historically limited to synchrotron facilities. Here, we present the first demonstration of laboratory-scale 3DXRD (Lab-3DXRD) using a liquid-metal-jet source. Lab-3DXRD achieves accuracy comparable to synchrotron-based 3DXRD, as validated against laboratory diffraction contrast tomography (LabDCT) and synchrotron-3DXRD. Over 96% of the grains detected with Lab-3DXRD are cross-validated, particularly for coarse grains (> ~60 μm), while the results suggest that finer grains should be accessible by taking advantage of high-efficiency detectors. We further demonstrate that its sensitivity to finer grains is enhanced by incorporating pre-characterization into the analysis. This study establishes Lab-3DXRD as a practical alternative to synchrotron techniques, making 3DXRD accessible to a wider range of academic and industrial researchers.

characterization and analytical techniques

MRV Challenge 3: velocity and passive scalar comparison in a 3D turbulent flow

The third iteration of a magnetic resonance velocimetry (MRV) challenge comparison activity consisting of 3D flow measurements for a turbulent flow in a water channel moving past centrally positioned cubic obstacles is reported. In this challenge iteration, MRV measurements are coupled with either temperature or concentration measurements to extend the diagnostic utility. Five research teams from around the world conducted the measurements on a single shared apparatus. The water channel included partial elements along the channel sidewalls that varied in height and precluded easy optical accessibility. In addition to a turbulent mainstream flow, a secondary flow entered the channel between the second and third cubic elements from a square injector hole at the channel bottom wall. The injector jet interacts with the mainstream flow and mixes turbulently as it advects downstream. For the selected flow regime and water solutions used as working fluids, the mixing of the higher temperature or concentration through the secondary flow with the mainstream flow satisfies the same dynamics so that both temperature and concentration fields behave as passive scalars and can be directly compared. The measurements are explained in detail by each participating team, and the results are interpolated onto a common grid and compared using line profiles, contour plots, and isosurfaces.

Benson, Mike [ORNL] (ORCID:000000023210116X)

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato