Search NASASearch

SEARCH · Search NASA

Results for “Simulation and optimization”

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.

1,667 records · Page 16

Two-dimensional non-equilibrium melting of charged colloids

Thermodynamic two-dimensional melting has been extensively studied in experiments and simulations, and is well predicted by theory. For systems in equilibrium, this transition is well described by the Kosterlitz–Thouless–Halperin–Nelson–Young theory, where melting is directly linked to the unbinding of topological defects. For driven, non-equilibrium melting and other non-equilibrium phase transitions, the picture is less clear. Here, in this work, we study the two-dimensional melting of a crystal of charged colloids. By randomly replacing some charged colloids with magnetic colloids, we can melt our system by rotating a fraction of the particles to create non-equilibrium, hydrodynamic random flows and local stresses. We can also melt it thermally by changing the particle number density. We find that an effective temperature approach cannot explain the results of our driven system. Rather, in both experiments and simulations, we observe that plotting the hexatic order parameter and the hexatic correlation’s exponent versus the density of disclinations and dislocations, respectively, yields universal curves. This implies that in our systems, two-dimensional melting depends directly on the density of topological defects and is independent of whether thermal or non-equilibrium forces generate them.

critical phenomena

Mechanism of H 2 plasma-enabled reduction of hematite thin films

Hydrogen plasma is gaining significant interest as a promising pathway for direct iron ore reduction and for lowering process temperatures but the reduction mechanism remains poorly understood. In this work we analyzed the plasma and thermal reduction of thin-film hematite (Fe 2 O 3 ) at temperatures below 340 °C using X-ray diffraction and scanning electron microscopy with energy-dispersive X-ray spectroscopy. Plasma reduced the incubation period by an order of magnitude and increased the reduction rate by a factor of 2.6 compared to thermal reduction. Plasma-produced H-atoms facilitate the formation of numerous iron nucleation sites, bypassing the energetically unfavorable dissociative adsorption of H 2 on iron oxide. These iron nuclei can autocatalyze the reduction of the surrounding hematite via a hydrogen spillover mechanism. Our results demonstrate that plasma-derived H-atoms primarily impact the initial nucleation-limited stage. These new insights provide a mechanistic framework that can aid the implementation and optimization of hydrogen plasma-assisted iron oxide reduction at reduced temperatures.

08 HYDROGEN

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E

Cooperative and bifunctional Ga-Ca-Cr 2 O 3 @CaO structured monoliths as versatile platform for reactive capture of CO 2 and its subsequent conversion to ethylene

Cooperative and bifunctional materials (BFMs) that integrate adsorbents and catalysts offer a promising strategy for the reactive capture of CO 2 to produce valuable fuels and chemicals. In this study, we developed structured BFMs via 3D printing that combine CaO as an adsorbent with Ga–Ca–Cr 2 O 3 metal oxides as the catalyst for the reactive capture of CO 2 and its subsequent conversion to C 2 H 4 via the oxidative dehydrogenation of C 2 H 6 (CO 2 -ODHE). Three different Ga–Ca compositions were used to modify the catalyst surface characteristics and enhance C 2 H 4 selectivity. In these formulations, Ga ions stabilize the oxygen lattice of the BFM, while Ca ions interact strongly with Cr to form CaCrO 4 , thereby altering the oxygen species and enhancing the material’s basic properties. Under adsorption–reaction conditions at 600–650 °C, the optimal BFM achieved an excellent C 2 H 4 selectivity of 96.4 %, attributed to a balanced redox process and improved basicity that facilitate efficient C 2 H 6 conversion and rapid desorption of C 2 H 4 without excessive oxidation. Overall, this work provides new insights into the formulation of BFMs monoliths and highlights the critical role of catalytic surface modification in enhancing C 2 H 4 selectivity in the CO 2 -ODHE reactive capture process.

C2H4 production

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code

Humins-Derived Hard Carbon as a Low-Cost Material for Sodium-Ion Battery Anodes

The growing demand for sodium-ion batteries (SIBs) in grid storage underscores the need for electrode materials that balance performance and cost, including sustainable and robust carbon sources. The equitable and abundant distribution of materials for SIBs, along with their superior low-temperature performance, safety, and fast-charging capability, further distinguishes them from lithium-ion batteries (LIBs). Hard carbon (HC) is the state-of-the-art anode for SIBs, but current commercial HC production is localized mostly to one region of the world, raising concerns about supply chain vulnerability, critical material dependency, and environmental aspects. Here, the first demonstration of humins, an abundant biorefinery byproduct, as a precursor for HC anodes for sodium-ion storage is reported. Humins were carbonized at 1100 degrees C-1300 degrees C, and the resulting materials were subjected to comprehensive materials and electrochemical characterization. Among the temperatures studied, humins-derived hard carbon synthesized at 1200 degrees C delivers an initial reversible capacity 270mA h g-1 , with stable cycling performance up to 500 cycles and excellent rate capability, representing the optimal performance. This study establishes humins as a promising and low-cost carbon source that provides a route to mitigate supply chain risks and valorizes an underutilized biorefinery waste stream for high-performance SIB anodes.

25 ENERGY STORAGE

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE

Constraints on Nuclear Symmetry Energy Parameters

A review is made of constraints on the nuclear symmetry energy parameters arising from nuclear binding energy measurements, theoretical chiral effective field predictions of neutron matter properties, the unitary gas conjecture, and measurements of neutron skin thicknesses and dipole polarizabilities. While most studies have been confined to the parameters 𝑆 𝑉 and L, the important roles played by, and constraints on 𝐾 sym , or, equivalently, the neutron matter incompressibility 𝐾 𝑁 , are discussed. Strong correlations among 𝑆 𝑉 ,𝐿, and 𝐾 𝑁 are found from both nuclear binding energies and neutron matter theory. However, these correlations somewhat differ in the two cases, and those from neutron matter theory have smaller uncertainties. To 68% confidence, it is found from neutron matter theory that 𝑆 𝑉 = 32.0 ± 1.1 MeV, 𝐿 = 51.9 ± 7.9 MeV and 𝐾 𝑁 = 152.2 ± 38.1 MeV. Theoretical predictions for neutron skin thickness and dipole polarizability measurements of the neutron-rich nuclei 48 Ca, 120 Sn, and 208 Pb are compared to recent experimental measurements, most notably the CREX and PREX neutron skin experiments from Jefferson Laboratory. By themselves, PREX I+II measurements of 208 Pb and CREX measurement of 48 Ca suggest 𝐿 = 121 ± 47 MeV and 𝐿 = −5 ± 40 MeV, respectively, to 68% confidence. However, we show that nuclear interactions optimally satisfying both measurements imply 𝐿 = 53 ± 13 MeV, nearly the range suggested by either nuclear mass measurements or neutron matter theory, and is also consistent with nuclear dipole polarizability measurements. This small parameter range implies 𝑅 1.4 = 11.6 ± 1.0 km and Λ 1.4 = 228$^{+148}_{−90}$, which are consistent with NICER X-ray and LIGO/Virgo gravitational wave observations of neutron stars.

neutron skins

A Cross‐Linked Flexible Metaferroelectrolyte Regulated by 2D/2D Perovskite Heterostructures for High‐Performance Compact Solid‐State Sodium Batteries

Abstract To address the issues of limited ionic conductivity and poor interface stability at room and low temperatures in solid‐state electrolytes, a robust intrinsic ferroelectrolyte or nanoferroelectrolyte strategy for engineering solid‐state flexible ferroelectric composite electrolytes utilizing strongly coupled intrinsic ion conducting 2D/2D sodium‐rich anti‐perovskite (NaRAP)/ferroelectric perovskite heterostructures is introduced. Herein, highly scalable PVDF‐based metaferroelectrolytes with Na 2.99 Ba 0.005 OCl/Ca 2 Na 2 Nb 5 O 16 − (CNNO − ) nanosheets into a ferroelectric poly(vinylidene fluoride‐co‐hexafluoropropylene) (PVDF‐HFP) matrix, through an in situ cross‐linking and spontaneous bridging method, for compact solid‐state sodium batteries (SSBs), are reported. Benefiting from unique well‐dispersed 3D ferroelectric coupled network and the Na 2.99 Ba 0.005 OCl/CNNO − ‐induced PVDF‐HFP ferroelectric β phase, the Na + flux is regulated, thereby inhibiting Na dendrite growth at the interface. Notably, the optimized PH‐5% NC metaferroelectrolyte exhibits rapid ion transport (1.11 × 10 −4 S cm −1 at 25 °C), a wide electrochemical window (> 4.8V), superior conformal mechanical compatibility, improved flexibility, good elasticity and flame retardancy. The solid‐state Na 3 V 2 (PO 4 ) 3 /PH‐5% NC/Na batteries present a stable cycling performance (remaining 56.4 mAh g −1 after 500 cycles at 1 C) even at 0 °C, potential for cost‐effective, safe, stable and compact SSB energy storage over 600 Wh L −1 , vastly surpassing 365 Wh L −1 of the current commercial sodium‐ion liquid‐electrolyte batteries.

Chemistry

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality

Used Nuclear Fuel Management Using the Next Generation System Analysis Model

The U.S. Department of Energy (DOE) is leading the National effort to manage the back end of the nuclear fuel cycle, encompassing the safe transportation, storage/staging, and/or eventual disposal of used nuclear fuel (UNF) and high-level radioactive waste. The Next Generation System Analysis Model (NGSAM) is DOE’s discrete-event, agent-based simulation tool designed to model the full life cycle of UNF from reactor discharge to final disposal. NGSAM supports the DOE Office of Spent Fuel and High-Level Waste Disposition by enabling a detailed, scenario-based analysis of logistics, infrastructure, and shipping strategies. NGSAM replaces legacy models with a modern, flexible platform built on Repast Simphony and enhanced by the Process Analysis Tool. NGSAM simulates the movement and interaction of individual fuel assemblies with system components such as canisters, casks, railcars, and facilities. The model integrates with the Java Transportation Operations Model to plan and execute transportation scenarios, supporting both constrained and unconstrained resource allocation. Key features include customizable allocation and acceptance algorithms, detailed facility-level operations, and a Quick Edit tool for rapid scenario adjustments. NGSAM supports multimodal transportation modeling (e.g. rail, road, barge) and provides comprehensive cost, schedule, and infrastructure data. NGSAM utilizes data from sources such as DOE’s STANDARDS UNF database and DOE’s Stakeholder Tool for Assessing Radioactive Transportation, while also allowing user-defined inputs for scenario customization. NGSAM enables stakeholders to evaluate complex UNF management strategies, assess system performance under varying assumptions, and inform decision making for future infrastructure investments. Its modular architecture and integration with other Integrated Waste Management System tools make it a critical asset for planning the safe and efficient disposition of the Nation’s growing UNF inventory.

Craig, Brian [Argonne National Laboratory (ANL)]

Amphiphilic nanopores that condense undersaturated water vapor and exude water droplets

Condensation of water vapor in confined geometries, known as capillary condensation, is a fundamental phenomenon with far-reaching implications. While hydrophilic pores enable liquid formation from undersaturated vapor without energy input, the condensate typically remains confined, limiting practical utility. Here, we explore the use of amphiphilic nanoporous polymer-infiltrated nanoparticle films that condense and release liquid water under isothermal and undersaturated conditions. By tuning the polymer fraction and nanoparticle size, we optimize condensation and droplet formation. As vapor pressure increases, voids fill with condensate, which subsequently exudes onto the surface as microscopic droplets. This behavior, enabled by a balance of polymer hydrophobicity and capillarity, reveals how amphiphilic nanostructures can drive accessible water collection. Our findings provide design insights for materials supporting energy-efficient water harvesting and heat management without external input.

Science & Technology - Other Topics

Pinning ångström-size solid ionic channels for rare-earth element separation

High-purity rare-earth elements are essential for modern technologies, yet current solvent extraction processes are energy-intensive and environmentally harmful because of inadequate selectivity and ligand toxicity. Although combining size exclusion and binding affinity can improve lanthanide separation, the role of long-range confinement remains underexplored. Here we report lanthanide separation in aqueous systems using extremely confined manganese oxide solid ionic channels with optimized layer spacing. Different lanthanides induce distinct solid-state phase transformations in manganese oxide, creating a strong driving force for separation. Two lanthanide groups, differing by ~1.4 Å in spacing, were identified and confirmed to be stable by density functional theory. The narrower confinement of heavier Group II lanthanides improves cross-group separation by increasing the dehydration barrier for lighter Group I lanthanides without inducing strong binding. Here, we further developed a strategy to pin the confinement dimensions and enhance same-group separation, increasing enrichment factors for La–Nd and La–Pr pairs from 1.6 ± 0.1 and 1.5 ± 0.1 to 5.4 ± 0.1 and 4.2 ± 0.1, respectively.

Chemical engineering

Compatibility of Type 304H Stainless Steel in Static Molten FLiBe for Inertial Confinement Fusion Reactors: Role of Impurities and Redox Control

Molten fluoride salts, such as FLiBe (LiF-BeF2), are promising candidates for tritium breeding and heat transfer in fusion reactors, but corrosion of structural materials remains a major challenge. This study investigates the corrosion behavior of austenitic stainless steel 304H in purified and NiF2-containing FLiBe at 500°C and 600°C, focusing on the effects of impurities and redox control. Exposure to purified FLiBe resulted in the concurrent depletion of Cr, Mn, and Fe, with corrosion at 500 °C dominated by the combined oxide formation and elemental dissolution, while at 600°C elemental depletion was predominant. The addition of a controlled NiF2 impurity significantly accelerated corrosion at both temperatures, demonstrating the sensitivity of 304H to the salt redox state. Beryllium additions were effective in mitigating corrosion for both baseline and NiF2-containing FLiBe; minimal depletion of Cr, Mn, and Fe occurred with Be additions as low as 2.5 mg (147 wppm), and no NiBe intermetallics formed at 5 mg (294 wppm), indicating that small Be inventories can provide substantial protection without deleterious phase formation. Thermodynamic equilibrium and coupled thermodynamic-kinetic analyses at the salt-alloy interface suggested low corrosion rates in systems with limited hydrogen fluoride (HF) generation, highlighting the importance of salt redox control. Estimates for a Be addition rate were calculated for the HYLIFE-II fusion reactor that can mitigate corrosion-induced degradation, assuming complete conversion of tritium to tritium fluoride (TF). Overall, 304H shows reasonable compatibility with FLiBe under optimized redox conditions. These results provide quantitative guidance for material selection and salt management in fusion blanket and heat exchanger systems and motivate validation under flowing, nonisothermal, and irradiated conditions.

Pillai, Rishi [ORNL] (ORCID:0000000243688197)

Selectivity mechanisms of ion intercalation in Prussian blue analogs

Prussian blue analogs (PBAs) are a family of materials with facile, reversible, and selective ion transport capability for various ions via electrochemical intercalation, owing to their vacancy structure. The large tunable compositional space of PBAs allows for manipulation of intercalation behavior and selectivity by controlling structural vacancy level through choice of transition metal centers and modifications to the synthesis process. However, a lack of understanding of the mechanisms of ion selectivity hinders the material’s design process. Here, for this work, we investigated the origins of ion selectivity using a model PBA, copper hexacyanoferrate, and focused on eight technologically and biologically prominent ions, for which we determined a sequence of selectivity: Rb + > K + > Na + > Ba 2+ > Sr 2+ ≈ Ca 2+ > Mg 2+ > Li + . We provide electrochemical, structural, and redox evidence of strong correlation between the ion identity, the dominant charge-compensating redox, and preferred occupancy site. Specifically, using synchrotron anomalous X-ray diffraction (AXRD), we reveal that monovalent ions exhibit significant association with the corner sites of the unit cell and iron redox, whereas divalent ions display affinity toward the center site with higher ratios of copper redox. Informed by selectivity results, we applied CuHCFe to Li purification and achieved 99.9% purity. Our findings demonstrate an approach to elucidating ion intercalation behavior in order to distinguish and manipulate material properties to optimize separation performance.

Prussian blue analog