Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hierarchical control”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Self‐Limiting Polymerization‐Induced Crystallization‐Driven Self‐Assembly (SL‐PI‐CDSA) Enables Templated Synthesis of Chiral Plasmonic, Hybrid 2D Hexagonal Assemblies

Examples of self-regulating synthetic self-assembly are relatively few, with most known chemical systems relying on kinetic rather than thermodynamic control. Herein, we demonstrate the rapid generation (t ≤ 5 min) of size-tunable ultralow dispersity (Ð ≤ 1.01) 2D hexagonal nanosheets governed by self-limiting self-assembly (SLSA). Self-assembly in natural systems occurs with exquisite control of structure, function, and dimension. We demonstrate that key aspects of biological assembly can be rationally applied toward the development of bottom-up approaches for the construction of chiral nanomaterials. To this end, self-limiting polymerization-induced crystallization-driven self-assembly (SL-PI-CDSA) of modular and templating aryl isocyanide (AIC) monomers yields functional 2D assemblies permitting post-polymerization/assembly modifications. Detailed study of the internal and external structure of the hexagonal nanosheets reveals topological defects which offer mechanistic insights into both their assembly and subsequent utilization. Specifically, these features enable fabrication of chiral hybrid organic–inorganic nanomaterials incorporating chiral plasmonic metal nanoparticles (MNPs). Our results suggest that the synergistic interplay of template-driven confinement and hierarchical chirality induce symmetry breaking of in situ-generated gold MNPs. We anticipate that the platform presented will facilitate fabrication of new hybrid, chiral organic–inorganic nanostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Surfactant-Free Synthesis of Three-Dimensional Perovskite Titania-Based Micron-Scale Motifs Used as Catalytic Supports for the Methanol Oxidation Reaction

We synthesized and subsequently rationalized the formation of a series of 3D hierarchical metal oxide spherical motifs. Specifically, we varied the chemical composition within a family of ATiO 3 (wherein “A” = Ca, Sr, and Ba) perovskites, using a two-step, surfactant-free synthesis procedure to generate structures with average diameters of ~3 microns. In terms of demonstrating the practicality of these perovskite materials, we have explored their use as supports for the methanol oxidation reaction (MOR) as a function of their size, morphology, and chemical composition. The MOR activity of our target systems was found to increase with decreasing ionic radius of the “A” site cation, in order of Pt/CaTiO 3 (CTO) > Pt/SrTiO 3 (STO) > Pt/BaTiO 3 (BTO). With respect to morphology, we observed an MOR enhancement of our 3D spherical motifs, as compared with either ultra-small or cubic control samples. Moreover, the Pt/CTO sample yielded not only improved mass and specific activity values but also a greater stability and durability, as compared with both commercial TiO 2 nanoparticle standards and precursor TiO 2 templates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Cobalt-Free Cathodes for Next Generation Li-Ion Batteries

In this U.S. Department of Energy sponsored project Nexceris, in collaboration with project partners; The Ohio State University and Navitas Advanced Systems have advanced the technical maturity of a non-cobalt containing cathode for next-generation Li-ion batteries. The cathode is based on the lithium manganese nickel-titanium oxide, LiNi 0.5 Mn 1.5 TiO 4 (LNMTO) high voltage spinel. To address limitations with poor cycle and calendar life an microstucturally hierarchical LNMO/LNMTO core-shell cathode powder has been developed that enables the formation of a solid-electrolyte interface that effectively passivates the cathode surface. The microstructural enhancements of the cathode material focus on preferentially enriching the surface with titanium. In parallel, new, optimized binder and electrolyte chemistries have been incorporated to address degradation mechanisms associated with high-voltage systems. Single-layer pouch cell and large-format 2-Ah cell testing have shown that an optimized LNMO/LNMTO core-shell powder significantly improves initial cell capacity and cycle life compared to homogeneous LNMO powder. To support development and the fabrication of 2-Ah cells a novel Hybrid Alternative Wet-Chemical Synthesis (HAWCS) process have been developed. This low-cost, synthesis approach enables the excellent compositional and particle morphology control achieved with co-precipitation without the strict process controls and associated expensive process equipment.

25 ENERGY STORAGE↗
DOE OSTI2020

Layer size polydispersity in hydrated montmorillonite creates multiscale porosity networks

The aluminosilicate layers of the swelling clay mineral montmorillonite, and the saturated pores they delineate, control the mechanical properties and the transport of solutes in many natural and engineered environments. However, the structural basis of montmorillonite porosity remains poorly characterized due to the difficulty in visualizing hydrated samples in their native state. Here in this paper, we used cryogenic transmission electron microscopy (cryo-TEM) and cryo electron tomography (cryo-ET) to show that stacking defects within minimally altered, fully hydrated montmorillonite particles define multiscale porosity networks. Variations in layer lateral dimensions over tens to thousands of nanometers cause a range of topological and dynamic defects that generate pervasive curvature and introduce previously uncharacterized solute transport pathways. Observations of long-range rotational order between neighboring layers indicate that the layer-layer interactions that govern clay swelling involve three dimensional orienting forces that operate across nanoscale pores. These direct observations of the hierarchical structure of hydrated montmorillonite pore networks with nanoscale resolution reveal potentially general aspects of colloidal interactions in fluid-saturated clay minerals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rational Design of Novel Biomimetic Sequence-Defined Polymers for Mineralization Applications

Silica biomineralization is a naturally occurring process, wherein organisms use proteins and other biological structures to direct the formation of complex, hierarchical nanostructures. Discovery and characterization of such proteins and their underlying mechanisms spurred significant efforts to identify routes for biomimetic mineralization that reproduce the exquisite shapes and size selectivities found in nature. A common strategy has been the use of short peptide sequences with chemistry mimicking those found in natural systems, such as the use of the silaffin-derived R5 peptide. While progress has been made using this approach, there are many limitations that have prevented breakthroughs in biomimicry. To advance our ability to use charged macromolecules for silica formation, we propose to use sequence-defined synthetic polymers known as peptoids, or N-substituted polyglycines, which present significant capability for the precise tuning of sequence and structure beyond what can often be achieved with peptides alone. This study presents a computationally predicted design of these polymers that leads to the controlled formation of silica nanomaterials. We investigate surface adsorption and the mineralization process through analysis of binding mechanisms and energetics of the R5 system. Next, we synthesized two R5-inspired peptoids and validated our prediction in the design of mineralization polymers through characterization using surface plasmon resonance and electron microscopy. Here, this computationally guided study holds great promise for designing new sequences with unprecedented control of the placement of chemical functional groups, thus allowing for further unraveling of silicification mechanisms and the eventual design of sequence-defined synthetic polymers leading to the predictive synthesis of nanostructured functional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗
DOE OSTI2021

A Sparse-Grid Probabilistic Scheme for Approximation of the Runaway Probability of Electrons in Fusion Tokamak Simulation

Runaway electrons (RE) generated during magnetic disruptions present a major threat to the safe operation of plasma nuclear fusion reactors. A critical aspect of understanding RE dynamics is to calculate the runaway probability, i.e., the probability that an electron in the phase space will runaway on, or before, a prescribed time. Such probability can be obtained by solving the adjoint equation of the underlying Fokker-Planck equation that controls the electron dynamics. In this effort, we present a sparse-grid probabilistic scheme for computing the runaway probability. The key ingredient of our approach is to represent the solution of the adjoint equation as a conditional expectation, such that discretizing the differential operator reduces to the approximation of a set of integrals. Adaptive sparse grid interpolation is utilized to approximate the map from the phase space to the runaway probability. The main novelties of this effort are the integration of the sparse-grid method into the probabilistic numerical scheme for computing escape probability, and the application of the proposed method in computing RE probabilities. Two numerical examples are given to illustrate that the proposed method can achieve O(Δt) convergence, and that the local anisotropic adaptive refinement strategy (M. Stoyanov, Adaptive sparse grid construction in a context of local anisotropy and multiple hierarchical parents. In: Sparse Grids and Applications-Miami 2016, Springer, Berlin, 2018, pp. 175–199) can effectively handle the sharp transition layer between the runaway and non-runaway regions.

Yang, Minglei↗
DOE OSTI2024

Combining solution-, precipitation- and load-transfer strengthening in a cast Al-Ce-Mn- Sc -Zr alloy

Here, a cast Al-9Ce-0.75Mn-0.18Sc-0.12Zr (wt%) alloy is designed to combine three strengthening phases: (i) micron-scale Al 11 Ce 3 platelets formed during eutectic solidification, (ii) nano-scale L1 2 -Al 3 (Sc,Zr) precipitates formed during aging, and (iii) Mn in solid solution in the α-Al matrix. Microstructural analyses by SEM, TEM, and atom-probe tomography reveal that Mn remains in solid solution in the as-cast alloy, providing solution strengthening with no influence on the eutectic Al-Al 11 Ce 3 microstructure, which provides precipitation- and load-transfer strengthening. During long-term over-aging at 400 °C, Mn-rich precipitates grow at the Al-Al 11 Ce 3 interface, with no effect on the microhardness. However, after short aging at 350 °C, a high number density of fine L1 2 -Al 3 (Sc,Zr) nanoprecipitates form in the Al matrix (with a coarser size at the Al-Al 11 Ce 3 interface), providing precipitation strengthening. The synergistic combination of the three strengthening mechanisms (solution, precipitation, and load transfer) in our Al-Ce-Mn-Sc-Zr alloy results in higher microhardness after aging at 350 and 400 °C, and higher creep resistance at 300 °C, as compared to alloys with two strengthening mechanisms: an Al-10Ce-0.93Mn control alloy (without precipitation strengthening from Sc and Zr), Al-Ce-Sc-Zr (without solution strengthening from Mn), and Al-Mn-Zr-Er (without load-transfer strengthening from Ce). Furthermore, these dual-strengthened alloys are more creep resistant than alloys with a single strengthening mechanism (Al-Ce, Al-Mn, and Al-Sc-Zr), confirming that the three mechanisms can be combined in pairs or all together.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From Cell to System: Accelerated hpc Simulations of BESS Aging under Frequency Regulation and Arbitrage use cases

Lithium-ion battery energy storage systems (BESS) packs have emerged as a leading solution for grid-scale energy storage, enhancing resiliency and balancing load fluctuations. Yet, experimental characterization of large-format LIB packs-particularly to assess performance and degradation over hundreds of cycles - demands substantial hardware investment and multi-year testing campaigns. In this work, we couple a hierarchical, physics-based modeling framework agnostic to electrode chemistries with high-performance computing to accelerate systems level evaluation by upto two orders of magnitude. Building on the open-source liionpack platform, we implement cell, module, and pack-scale electrochemical models enriched with mechanistic aging mechanisms and deploy them on an HPC cluster to simulate 150−200kWh systems over 500 - 1,000 cycles with in days. We subject these virtual B ESS to both constant-current cycling and realistic grid service profiles spanning frequency regulation, ramp-rate support, and energy arbitrage-and quantify the resulting degradation patterns. Our results reveal that localized cell aging can induce substantial nonuniformity at module and pack levels, with service-specific cycling protocols driving distinct aging modes. This rapid, multiscale modeling approach provides a powerful design-space exploration tool for optimizing electrical architecture, control strategies, and operational schedules to prolong pack lifetime and lower total cost of ownership.

Ayalasomayajula, Surya [ORNL] (ORCID:0009000860788↗
DOE OSTI2021

The pervasive impact of critical fluctuations in liquid–liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Altogether, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

An adaptive, data-driven multiscale approach for dense granular flows

The accuracy of coarse-grained continuum models of dense granular flows is limited by the lack of high-fidelity closure models for granular rheology. One approach to addressing this issue, referred to as the hierarchical multiscale method, is to use a high-fidelity fine-grained model to compute the closure terms needed by the coarse-grained model. The difficulty with this approach is that the overall model can become computationally intractable due to the high computational cost of the high-fidelity model. In this work, we describe a multiscale modeling approach for dense granular flows that utilizes neural networks trained using high-fidelity discrete element method (DEM) simulations to approximate the constitutive granular rheology for a continuum incompressible flow model. Our approach leverages an ensemble of neural networks to estimate predictive uncertainty that allows us to determine whether the rheology at a given point is accurately represented by the neural network model. Additional DEM simulations are only performed when needed, minimizing the number of additional DEM simulations required when updating the rheology. This adaptive coupling significantly reduces the overall computational cost of the approach while controlling the error. In addition, the neural networks are customized to learn regularized rheological behavior to ensure well-posedness of the continuum solution. We first validate the approach using two-dimensional steady-state and decelerating inclined flows. We then demonstrate the efficiency of our approach by modeling three-dimensional sub-aerial granular column collapse for varying initial column aspect ratios, where our multiscale method compares well with the computationally expensive computational fluid dynamics (CFD)-DEM simulation.

Dense granular flows↗
DOE OSTI2020

Non-Equilibrium Synthesis of Highly Active Nanostructured, Oxygen-Incorporated Amorphous Molybdenum Sulfide HER Electrocatalyst

Molybdenum sulfide emerged as promising hydrogen evolution reaction (HER) electrocatalyst thanks to its high intrinsic activity, however its limited active sites exposure and low conductivity hamper its performance. To address these drawbacks, the non-equilibrium nature of pulsed laser deposition (PLD) is exploited to synthesize self-supported hierarchical nanoarchitectures by gas phase nucleation and sequential attachment of defective molybdenum sulfide clusters. Furthermore, the physics of the process are studied by in situ diagnostics and correlated to the properties of the resulting electrocatalyst. The as-synthesized architectures have a disordered nanocrystalline structure, with nanodomains of bent, defective S-Mo-S layers embedded in an amorphous matrix, with excess sulfur and segregated molybdenum particles. Oxygen incorporation in this structure fosters the creation of amorphous oxide/oxysulfide nanophases with high electrical conductivity, enabling fast electron transfer to the active sites. The combined effect of the nanocrystalline pristine structure and the surface oxidation enhances the performance leading to small overpotentials, very fast kinetics (35.1 mV dec -1 Tafel slope) and remarkable long-term stability for continuous operation up to -1 A cm -2 . This work shows possible new avenues in catalytic design arising from a non-equilibrium technique as PLD and the importance of structural and chemical control to improve the HER performance of MoS-based catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Caldera_Grid

Caldera Grid is part of Caldera software platform, a suite of collective, open-source tools that was developed to improve the state of the art in modeling the impacts of Electric Vehicle (EV) charging on the grid. Caldera Grid is a co-simulation framework implemented using Hierarchical Engine for Large scale Infrastructure Co-Simulation (HELICS). The framework facilitates the co-simulation of EV charging and Smart Charge Management (SCM) strategies in Caldera Infrastructure Charge Model (ICM) with distribution level grid models in OpenDSS. Vehicle energy needs and charge session requirements generated using Caldera Charge Decision Module (CDM) are fed in as input to Caldera Grid. The EV charging models in Caldera ICM simulates both the uncontrolled charging loads as well as loads modified by the SCM strategies to develop distributed load profiles for each grid node hosting an Electric Vehicle Supply Equipment (EVSE) – also known as chargers. These loads were simulated in OpenDSS alongside existing distribution feeder loads. Each of these models are co-simulated in a HELICS federate. The HELICS co-simulation framework facilitates communication and synchronization between the federates. By co-simulating EV loads and distribution feeder loads, Caldera Grid can assess the potential grid impacts of EV charging on distribution feeders under various grid conditions. A control strategy federate is implemented with an interface where control strategies based on feedback from EV charging status and grid conditions can be developed and implemented. The platform can also support multiple control strategies in a single co-simulation.

Sundarrajan, ManojKumar Cebol↗
DOE OSTI2021

Enhancing the inherent catalytic activity and stability of TiO 2 supported Pt single-atoms at CeO x –TiO 2 interfaces

Single-atoms (SAs) with atomically coordinated reaction centers are considered the next generation of catalysts that can exhibit exceptional catalytic efficiency. However, the general concern about thermodynamic vulnerabilities of SAs questions their practical value. Moreover, whether the inherent catalytic nature of SAs is superior compared with that of larger nanoparticles is still under debate. Here, we address two controversies by a comparative study using two catalysts: Pt/TiO 2 and Pt/CeO x –TiO 2 . Based on a hierarchical study of density functional theory, time-resolved catalysis performance test, in situ infrared spectroscopy, and operando X-ray absorption spectroscopy, we could unveil the catalytic nature of Pt-SAs and their stability. By utilizing the heterogeneous interface formed between TiO 2 supporting particles and CeO x clusters formed on the surface of TiO 2 , we preferentially synthesized Pt-SAs pinned at the CeO x –TiO 2 interfaces on CeO x –TiO 2 hybrid-oxide supports. The strong electronic coupling between the Pt-SAs and the Ce ions at the CeO x –TiO 2 interfaces enhanced the catalytic activity toward CO oxidation of Pt-SAs and improved the long-term stability under CO oxidation conditions. The CO oxidation activity of Pt-SAs stabilized at CeO x –TiO 2 was improved by 13.5 times at 200 °C compared with the Pt-SAs on TiO 2 . The results present how to easily improve the activity and stability of Pt-SAs using a simple interface control method. Moreover, we demonstrate that the catalytic activity and the stability of Pt-SAs can be monitored through the chemical state of the interfaces. Finally, our study provides comprehensive understanding about the catalytic nature as well as a novel strategy toward applications of Pt-SAs, enabling sustainable use of Pt in heterogeneous catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

KGML-ag: a modeling framework of knowledge-guided machine learning to simulate agroecosystems: a case study of estimating N<sub>2</sub>O emission using data from mesocosm experiments

Abstract. Agricultural nitrous oxide (N2O) emission accounts for a non-trivial fraction of global greenhouse gas (GHG) budget. To date, estimating N2O fluxes from cropland remains a challenging task because the related microbial processes (e.g., nitrification and denitrification) are controlled by complex interactions among climate, soil, plant and human activities. Existing approaches such as process-based (PB) models have well-known limitations due to insufficient representations of the processes or uncertainties of model parameters, and due to leverage recent advances in machine learning (ML) a new method is needed to unlock the “black box” to overcome its limitations such as low interpretability, out-of-sample failure and massive data demand. In this study, we developed a first-of-its-kind knowledge-guided machine learning model for agroecosystems (KGML-ag) by incorporating biogeophysical and chemical domain knowledge from an advanced PB model, ecosys, and tested it by comparing simulating daily N2O fluxes with real observed data from mesocosm experiments. The gated recurrent unit (GRU) was used as the basis to build the model structure. To optimize the model performance, we have investigated a range of ideas, including (1) using initial values of intermediate variables (IMVs) instead of time series as model input to reduce data demand; (2) building hierarchical structures to explicitly estimate IMVs for further N2O prediction; (3) using multi-task learning to balance the simultaneous training on multiple variables; and (4) pre-training with millions of synthetic data generated from ecosys and fine-tuning with mesocosm observations. Six other pure ML models were developed using the same mesocosm data to serve as the benchmark for the KGML-ag model. Results show that KGML-ag did an excellent job in reproducing the mesocosm N2O fluxes (overall r2=0.81, and RMSE=3.6 mgNm-2d-1 from cross validation). Importantly, KGML-ag always outperforms the PB model and ML models in predicting N2O fluxes, especially for complex temporal dynamics and emission peaks. Besides, KGML-ag goes beyond the pure ML models by providing more interpretable predictions as well as pinpointing desired new knowledge and data to further empower the current KGML-ag. We believe the KGML-ag development in this study will stimulate a new body of research on interpretable ML for biogeochemistry and other related geoscience processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗
DOE OSTI2023

Roles of Hydrogen Bonds and Alignment in Oriented Attachment of Gibbsite Nanoparticles: Insights from Molecular Dynamics

Oriented attachment (OA) of nanoparticles is an important crystal growth pathway in the synthesis of hierarchical structures. Although a significant understanding of OA has been made, the effect of atomistic misalignment and the roles of solvent/particle and particle/particle interactions on the structure–energy relationship during an OA remain elusive. In this study, we perform molecular dynamics simulations to calculate the potential of mean force (PMF) profile for gibbsite particle translation on a gibbsite slab with 1 or 2 intervening water layers (1W or 2W). The structures of the gibbsite surfaces and the confined water are analyzed to determine how the number and type of hydrogen bonds (H-bonds) influence the free energy profile during the translation. The PMF profile exhibits a periodicity of length 5.078 Å, consistent with the gibbsite unit cell size along the translation direction. The changes in the surface–water and water–water hydrogen bond network and water and surface OH groups’ orientations during the translation are strongly coupled with the changes in the PMF profile in the 1W case. However, when increasing the number of intervening water layers from 1W to 2W, the particle/slab misalignment becomes a dominant factor controlling the behavior of the PMF profile. We also establish a method to quantify misalignment between the particle and the slab, which exhibits a strong correlation with the free energy for the 2W case. These results shed more light into the roles of particle/slab misalignment and hydrogen bond network in the OA of mineral particles in aqueous solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗