Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Integral Experiments”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Application of the Rossi-alpha method to simulations of HEU and organic scintillators

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook contains thousands of high-quality benchmark measurements that include detailed descriptions of experiment designs and exhaustive accounting of uncertainties. These benchmarks are used extensively by the criticality safety and nuclear data communities. A subcritical benchmark involving bare, highly enriched uranium (HEU) has yet to be accepted to the handbook. The Measurement of Uranium Subcritical and Critical (MUSiC) was proposed to address this gap. MUSiC will be performed at the National Criticality Experiments Research Center (NCERC) and will leverage bare HEU hemisphere shells that can be assembled in a variety of configurations. The measurement has the following key benefits: it will be the first set of fast, subcritical measurements of HEU to be submitted to ICSBEP; ten configurations will be measured spanning from nearly no multiplication to critical (e.g., it is a measurement of similar systems with different re activities); measurements of an all-HEU system will be of value for nuclear data validation, providing results that assist in determining biases in both measurements and simulations. A variety of detector systems will be used including a Neutron Multiplicity Array Detector (NoMAD) system (similar to the MC-15), four small volume 0.635 cm (#8960;) × 7.59 cm 3He detectors (ideal for measuring prompt neutron decay constants due to their fast recovery speed), and an array of eight 5.08 cm (#8960;) × 5.08 cm EJ-309 organic scintillator detectors. The focus of this work is the simulation of the subcritical MUSiC configurations including the EJ-309 detectors and the application of the Rossi-α method. The EJ-309 organic scintillator detector array is a new detection system that will be deployed at NCERC. Organic scintillators are ideal for measurements of fast, subcritical systems due to their short time response on the order of nanoseconds.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Neutron Absorber Plate Characterization Plan for Criticality Experiments Design

After being used in nuclear installations, depleted fuel can still be highly reactive and must be handled securely to prevent any radiological or criticality concerns. In particular, spent fuel from use in nuclear power reactors must be stored and transported in specifically designed containers using neutron absorber materials to prevent criticality. Various neutron absorber material types exist and are manufactured by various entities, as thoroughly described in the Handbook of Neutron Absorber Materials for Spent Nuclear Fuel Storage and Transportation Applications written by EPRI. Presently, one of the most modern and most widely used types of neutron absorber material contains particles of boron carbide, or B 4 C, embedded in aluminum matrix: Boralcan, manufactured by Rio Tinto. It is very important for the community to know as much as possible about such neutron absorber materials. Therefore, in the recent years, a US Department of Energy National Nuclear Security Administration–Nuclear Criticality Safety Program funded project initiated design of an experiment that places Boralcan neutron-absorbing plates in an established critical assembly using low-enriched uranium fuel at the Sandia Pulsed Reactor Facility/Critical Experiments (SPRF/CX) apparatus at Sandia National Laboratories. The goal of the experiment is to produce high-quality benchmark data to submit to the International Criticality Safety Benchmark Evaluation Project (ICSBEP), for use in validating calculational tools and nuclear data by criticality safety analysts. The project, named IER-554, is currently in its final design stage, following a successful preliminary design. In the work documented in the design study, ten critical configurations using Boralcan neutron absorber plates were designed, and the experiment was proven to be feasible, with a predicted low k eff uncertainty around 100 pcm. An overview of the modeled cutout of the critical assembly with a Boralcan plate is shown in Figure 1, representing one of the configurations planned for the critical experiments. Before the plates are inserted in the critical assembly, it is necessary to know more about their composition and uniformity. This summary focuses on the plate characterization plans. Each plate will undergo (1) neutron transmission measurements at different locations to determine the 10 B areal density and (2) an in-depth x-ray computed tomography (XCT) examination to obtain the exact Sizes and distribution of the B4C powder particles inside the plates. In parallel, plate modeling studies are performed with a goal to determine the validity of the currently used approximation of modeling the neutron absorber plates as a homogeneous mixture of Aluminum 1100 alloy and B4C— instead of explicitly modeling the B4C particles. By using the experimental 10 B areal density measurements, and the exact size and location of the B4C particles obtained by XCT, a plate model can theoretically be built that reproduces the plate with extremely high fidelity. The results of this modeling study could increase the confidence of the criticality safety community in its modeling methods when using this type of neutron absorber material, and the industry could use these validations to change the boron loading credit limits from the U.S. Nuclear Regulatory Commission standard review plan for dry cask storage of spent nuclear fuel. The modeling calculations are performed with SCALE 6.3.0 using the KENO V.a sequence for criticality calculations with the ENDF/B-VIII.0 continuous-energy cross section library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Preliminary RAM-RODD results for the MUSiC subcritical configurations

The Measurement of Uranium Subcritical and Critical (MUSiC) was performed at the DOE’s National Criticality Experiments Research Center (NCERC) located in the Nevada National Security Site (NNSS). The measurement utilized the Rocky Flats shells to perform benchmark measurements of similar highly enriched uranium (HEU) systems that span a wide range of reactivities. The Rocky Flats (RF) shells are 93.16% U-235 enriched metal hemishells that can be stacked concentrically. Ten configurations were measured with effective multiplication factors spanning between deeply subcritical (~ 0.64) through delayed critical. Details of the measured configurations are listed in Table 1. This unique set of measurements with its large span of reactivies is being used to determine the range over which neutron noise techniques such as Feynman variance-to-mean, Rossi-alpha, and pulsed neutron source techniques can be accurately employed for a bare HEU system. The results of the measurements will be published as a benchmark in The International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The results will support the growing amount of subcritical benchmark data, such as the SCRaP measurements, that is available to the community. The measurements were performed using several different detector systems for the purposes of cross-validation and obtaining detector independent results. Four detector systems were deployed, three by the NCERC team and one from the University of Michigan. The detector systems included a Neutron Multiplicity Array Detector (NoMAD) system (similar to the MC-15), four small volume 0.635 cm (Ø) × 7.59 cm 3 He detectors (ideal for measuring prompt neutron decay constants due to their fast recovery speed), the Rossi Alpha Measurements – Rapid Organic (n, γ) Discrimination Detector (RAM-RODD), and the Organic Scintillator Array (OSCAR). RAM-RODD is an array of eight 5.08 cm (Ø) × 5.08 cm EJ-309 organic scintillator detectors. OSCAR is a University of Michigan system and is an array of twelve 5.08 cm (Ø) × 5.08 cm stilbene detectors. Details of measurements performed with OSCAR will be discussed in a separate talk. The focus of this work is preliminary results obtained by RAM-RODD for the 8 subcritical configurations. Additional details on the measurements and on the setup and deployment of RAM-RODD will be discussed. Preliminary neutron noise analysis results including Rossi-alpha and Feynman-alpha will also be presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Unraveling the Highly Plastic Behavior of ALD-Aluminum Oxide Encapsulations by Small-Scale Tensile Testing

We present a study directly measuring the electron-beam-induced plasticity of amorphous Al 2 O 3 coatings. Core–shell nanostructures are employed as small-scale model systems for two-dimensional coatings made by atomic layer deposition (ALD). Copper nanowires (NWs) are used as substrates for ALD deposition, representing a model system for interconnects commonly found in integrated circuits. Experiments are performed in situ in a transmission electron microscope (TEM) and further analyzed with electron energy loss spectroscopy (EELS). Our in situ TEM tensile experiments reveal the highly plastic behavior of the ALD shell, which withstands a maximum strain of 188%. Comparable samples under beam-off conditions show a brittle fracture, which underlines the effect of electron irradiation. The electron-beam-activated bond switching within the amorphous network enables compensation of the applied tensile strain, leading to viscous flow. By incorporating an intermediate nanocrystalline layer within the Al 2 O 3 shell, the plasticity is suppressed and brittle fracture occurs. This work directly demonstrates the tuning of mechanical properties in amorphous ALD structures through electron irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Central Compact Finite‐Difference Scheme With High Spectral Resolution for KdV Equation

This work presents a combination of cell‐node and cell‐centered compact finite difference scheme for the approximation of third derivatives involved in Korteweg–de Vries (KdV) equations. This approach employs a half‐shifted derivative construction at cell centers, avoiding the need for compact interpolation, thereby removing transfer errors; hence, it improves spectral resolution and maintains high‐order accuracy. Fourier analysis is performed to show the spectral properties of the proposed formulation, which provides higher spectral resolutions as compared to node‐based compact schemes. A filtering strategy is incorporated to suppress high‐frequency oscillations without compromising the accuracy of the numerical scheme, and the total variation diminishing Runge Kutta (TVDRK3) method is applied for time integration. Numerical experiments on linear, nonlinear, and coupled KdV systems are conducted, and a comparative analysis with cell‐node compact schemes confirms that the proposed scheme consistently reduces errors by up to an order of magnitude and achieves high spectral resolution properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

On the Feynman-alpha method for reflected fissile assemblies

The Feynman-alpha method is a neutron noise technique that is used to estimate the prompt neutron period of fissile assemblies. The method and quantity are of widespread interest including in applications such as nuclear criticality safety, safeguards and nonproliferation, and stockpile stewardship; the prompt neutron period may also be used to infer the k eff multiplication factor. The Feynman-alpha method is predicated on time-correlated neutron detections that deviate from a Poisson random variable due to multiplication. Traditionally, such measurements are diagnosed with one-region point kinetics, but two-region models are required when the fissile assembly is reflected. This paper presents a derivation of the two-region point kinetics Feynman-alpha equations based on a double integration of the Rossi-alpha equations, develops novel propagation of measurement uncertainty, and validates the theory. Validation is achieved with organic scintillator measurements of weapons-grade plutonium reflected by various amounts of copper to achieve k eff values of 0.83–0.94 and prompt periods of 5–75 ns. The results demonstrate that Feynman-alpha measurements should use the two-region model instead of the one-region model. The simplified one-region model deviates from the validated two-region models by as much as 10% in the estimate of the prompt neutron period, and the two-region model reduces to the one-region model for small amounts of reflector. The Feynman-alpha estimates of the prompt neutron period are compared to those of the Rossi-alpha approach. The comparative results demonstrate that the Feynman-alpha method is more precise than the Rossi-alpha method and more accurate for k eff < 0.92, whereas the Rossi-alpha method is generally more accurate for higher multiplications. Here, the uncertainty propagation developed in this work should be used for all Feynman-alpha measurements and will therein improve fitting accuracy and appropriate precision estimates.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2026

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗
DOE OSTI2021

DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Enhancement of effective thermal conductivity of rGO/Mg nanocomposite packed beds

Engineering thermophysical properties of metal hydrides nanocomposites is crucial for effective thermal management during hydrogen storage reactions; however, the effect of microstructure on thermal transport mechanisms is still unclear. Here, we employed an integrated experiment-modeling approach to investigate microstructural factors that determine the effective thermal conductivity of individual reduced graphene oxide-magnesium (rGO/Mg) nanocomposites and their packed bed. Here the effective thermal conductivity of the rGO/Mg nanocomposite packed bed was measured by using guarded hot-plate method under various atmospheric conditions (i.e., vacuum, Ar and He). A microstructure-aware mesoscopic modeling revealed that anisotropy of the effective thermal conductivity of individual rGO/Mg nanocomposites plays an important role in determining the effective thermal conductivity of their packed bed. The validated mesoscopic model also disclosed a nontrivial interplay between the intrinsic rGO properties and the extrinsic composite structural features. Finally, quantitative sensitivity analysis based on the modeling framework is used to provide practical engineering guidance for controlling the thermal transport within nanocomposite packed beds.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interfacial characteristics of P91 steel - Inconel 740H bimetallic structure fabricated using wire-arc additive manufacturing

Although bimetallic structures can offer unique solutions to engineering problems with varying functionality, joining two dissimilar alloys using additive manufacturing involves challenges such as cracking and segregation. In this work, bimetallic structures of Inconel 740H superalloy and P91 steel are processed using wire-arc additive manufacturing (WAAM) for the first time. P91 steel is directly deposited over 740H with a cooling time of 20 minutes in between the two depositions. During microstructure characterization, a large gradient zone with coarse grains spanning for ~2200 μm is observed. Moreover, the microhardness of the gradient zone is the least in comparison with the P91 and 740H deposits. No deleterious phase formation is found; however, intergranular cracks are observed along the gradient zone. In order to deduce the reason for cracking, non-equilibrium solidification simulations are performed for the mixed compositions of P91 and 740H in the gradient zone. Here, the reasons for crack formation in the gradient zone are the sudden change in the volumetric coefficient of thermal expansion due to the formation of MC carbides as well as the development of local strains from the thermal residual stresses. Therefore, altering the deposition sequence and introducing gradient layers with mixed composition between the P91 and 740H deposits will be beneficial in achieving defect-free builds. Through the integration of experiments and simulation, this work gains valuable insights into the processability of dissimilar alloy manufacturing using the WAAM technique.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Rossi-alpha measurements of fast plutonium metal assemblies using organic scintillators

The reactivity and the k-effective multiplication factor (k eff ) of a fissionable assembly are quantities of widespread interest. These values can be inferred from Rossi-alpha measurements of the prompt neutron decay constant (or the inverse:α -1 ). It has been shown that 3 He-gas proportional counter-based detection systems are insensitive to α -1 of fast assemblies (much faster than tens of microseconds). Therefore, it is of interest to investigate fast detection systems such as those based on organic scintillation detectors. In this work, an array of 12 cylindrical, 5.08 cm x 5.08 cm diameter trans-stilbene organic scintillators was used to measure five subcritical assemblies. One assembly was a sphere of approximately 4.5 kg of alpha-phase, weapons-grade plutonium (k eff = 0.773, α -1 = 11.84 ns) known as the BeRP ball encased in a thin stainless-steel clad. The other assemblies used the same encased sphere and 7.62 cm of iron, nickel, copper, or tungsten reflectors (k eff = 0.884, 0.916, 0.924, 0.939, respectively and α -1 = 36.60, 41.56, 49.60, 70.32 ns, respectively). Finally, this work (1) validates Rossi-alpha measurements with organic scintillators by demonstrating good agreement between measurements and simulations; and (2), demonstrates that organic scintillator-based systems are sensitive to α -1 on the order of 10–100 ns.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Validation of the two-region Rossi-alpha model for reflected assemblies

Nuclear criticality safety, nonproliferation and safeguards, emergency response, and stockpile stewardship utilize estimates of the k eff multiplication factor. The value of k eff cannot be directly measured, but it can be inferred from the prompt neutron period. One modality of measuring the prompt neutron period is the Rossi-alpha method, which is predicated on fitting a histogram of detection time differences due to the nonrandom temporal distribution of same-fission-chain neutrons. Recent works have developed the motivation and theory to expand traditional one-region point kinetic models resulting in one-exponential histogram fits to two-region models that result in two-exponential fits. This paper validates the new two-region model using organic scintillator measurements of copper-reflected weapons-grade plutonium (0:83 ≤ k eff ≤ 0:94) and high-density- polyethylene-reflected highly enriched uranium (0:73 ≤ k eff ≤ 0:95). Furthermore, the results show that more thermal systems have shorter prompt neutron periods in the core region due to increased induced-fission probabilities for moderated neutrons. Here, a new parameter introduced by the two- region model is also shown to be correlated to the amount of reflection, and may be used to infer assembly properties such as type and amount of reflector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

SMR safety through HTTF modeling and benchmark efforts for code validation for gas-cooled reactor applications

Accurate modeling and simulation tools for thermal-hydraulics calculations are a key element needed to design and license new advanced reactors including Small Modular Reactors (SMR) and Microreactors. Uncertainties in modeling and simulation can have significant safety and economic implications. The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled prismatic-block nuclear reactors. High-quality measurement data is available from the HTTF that is suitable for a thermal-hydraulics code validation benchmark for gas-cooled reactor simulations. Here, this paper summarizes individual HTTF modeling efforts to date for tool validation at Idaho National Laboratory (INL), Argonne National Laboratory (ANL), Oregon State University (OSU) and Canadian Nuclear Laboratories (CNL) using system thermal-hydraulics codes, Computational Fluid Dynamics (CFD) codes and system-CFD code couplings. Also, the paper introduces the ongoing OECD Nuclear Energy Agency (NEA) High Temperature Gas Reactor Thermal-Hydraulics (HTGR T/H) benchmark that allows for better comparisons of results between different international modeling teams. The benchmark provides well defined computational problems that include code-to-code comparisons and comparisons to measured data. These problems provide an avenue for quantifying accuracy and identifying sources of uncertainty in thermal-hydraulics calculations, including in measured thermophysical properties, as part of validation for gas-cooled reactor simulation tools.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Ultralow-temperature cryogenic transmission electron microscopy using a new helium flow cryostat stage

Advances in cryogenic electron microscopy have opened new avenues for probing quantum phenomena in correlated materials. This study reports the installation and performance of a new side-entry condenZero cryogenic cooling system for JEOL (Scanning) Transmission Electron Microscopes (S/TEM), utilizing compressed liquid helium (LHe) and designed for imaging and spectroscopy at ultra-low temperatures. The system includes an external dewar mounted on a vibration-damping stage and a pressurized, low-noise helium transfer line with a remotely controllable needle valve, ensuring stable and efficient LHe flow with minimal thermal and mechanical noise. Performance evaluation demonstrates a stable base temperature of 4.37 K measured using a Cernox bare chip sensor on the holder with temperature fluctuations within ±0.004 K. Complementary in-situ electron energy-loss spectroscopy (EELS) via aluminum bulk plasmon analysis was used to measure the local specimen temperature and validate cryogenic operation during experiments. The integration of cryogenic cooling with other microscopy techniques, including electron diffraction and Lorentz TEM, was demonstrated by resolving charge density wave (CDW) transitions in NbSe2 using electron diffraction, and imaging nanometric magnetic skyrmions in MnSi via Lorentz TEM. In conclusion, this platform provides reliable cryogenic operation below 7 K, establishing a low-drift route for direct visualization of electronic and magnetic phase transformations in quantum materials.

Charge density wave↗
DOE OSTI2023

Trace Actinide Signatures of a Bulk Neptunium Sample

In nuclear forensic analyses, measurements of actinide elements in a sample can assist with identifying interdicted or unknown materials. While these radiochemical signatures have been extensively investigated in uranium materials, less is known about bulk neptunium samples. This paper describes the measurement of trace actinide concentrations and isotopic profiles in a 237 Np oxide sample. Uranium, plutonium, americium, and curium concentrations and isotopic profiles in the sample were determined and deemed potentially useful for distinguishing different sources of 237 Np. Furthermore, several different potential radiochronometry systems were also investigated; discordant results indicate that the Np sample was never completely purified of other actinide elements, or that subsequent contamination of the sample occurred. Few prior studies of neptunium materials have been reported, and these data suggest that trace actinide constituents could provide unique signatures to identify material out of regulatory control.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Synergetic Effects of Soil Organic Matter Components During Interactions with Minerals

Mineral-associated soil organic matter (SOM) is critical for stabilizing organic carbon and mitigating climate change. However, mineral-SOM interactions at the molecular scale, particularly synergetic adsorption through organic-organic interaction on the mineral surface known as organic multilayering, remain poorly understood. This study investigates the impact of organic multilayering on mineral-SOM interactions, by integrating macroscale experiments and molecular-scale simulations that assess the individual and sequential adsorption of major SOM compounds – lauric acid (lipid), pentaglycine (amino acid), trehalose (carbohydrate), and lignin onto soil minerals. Ferrihydrite, Al-hydroxide, and calcite are exposed to SOM compounds to determine adsorption affinities and binding energies. Results show that lauric acid has 20-40 times higher K d than pentaglycine, following the order K d (ferrihydrite) > K d (Al-hydroxide) >> K d (calcite). Molecular-scale simulations confirm that lauric acid has a higher binding energy (30.8 kcal/mol) on ferrihydrite than pentaglycine (6.0 kcal/mol), attributed to lipid hydrophobicity. The lower binding energy of pentaglycine results from its hydrophilic amide groups, facilitating partitioning into water. Sequential experiments examine how the first layer of lipid or amino acid affects the adsorption of carbohydrate/lignin, which show little or no individual adsorption affinities. Macroscale results reveal that lipid and amino acid adsorption induce ferrihydrite particle repulsion increasing reactive surface area and enhancing carbohydrate/lignin adsorption independently and synergistically through organic multilayering. Molecular-scale results reveal that amino acid adsorbed on ferrihydrite interacts more readily with lignin macroaggregates (preformed in solution) than with individual lignin units, indicating organic multilayering via H-bonding. Further, these findings reveal the molecular mechanisms of SOM-mineral interactions, crucial for enhancing soil carbon stabilization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗