Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scaling laws of complex systems”

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

Scaling Laws of Exciton Recombination Kinetics in Low Dimensional Halide Perovskite Nanostructures

Carrier recombination is a crucial process governing the optical properties of a semiconductor. Although various theoretical approaches have been utilized to describe carrier behaviors, a quantitative understanding of the impact of defects and interfaces in low dimensional semiconductor systems is still elusive. Here, we develop a model system consisting of chemically tunable, highly luminescent halide perovskite nanocrystals to illustrate the role of carrier diffusion and material dimensionality on the carrier recombination kinetics and luminescence efficiency. Our advanced synthetic methods provide a well-controlled colloidal system consisting of nanocrystals with different aspect ratios, halide compositions, and surface conditions. Using this system, we reveal the scaling laws of photoluminescence quantum yield and radiative lifetime with respect to the aspect ratio of nanocrystals. The scaling laws derived herein are not only a phenomenological observation but proved a powerful tool disentangling the carrier dynamics of microscopic systems in a quantitative and interpretable manner. The investigation of our model system and theoretical formulation bring to light the dimensionality, as a hidden constraint on carrier dynamics, and identify the diffusion length as an important parameter that distinguishes nanoscale and macroscale carrier behaviors. Furthermore, the conceptual distinction in carrier dynamics in different dimensionality regimes informs new design rules for optical devices where complex microstructures are involved.

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

Roadmap to Iodine and Mercury Abatement Materials Selection in Nuclear Waste Processing Off-Gas Streams

This work provides a guide for candidate mercury (Hg) and iodine-129 (129I) abatement material identification, screening, evaluation, technical gap identification, and bench-scale testing prioritization while also conceptualizing a materials deployment roadmap for implementation of new materials in the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility at Hanford in Richland, Washington, and elsewhere in the U.S. Department of Energy (DOE) complex for similar applications. The study was prompted by the need to replace the Kombisorb BAT-37 due to uncertainties in Hg and 129I capture performance and future availability for use in the WTP LAW Facility secondary off-gas system, specifically the Carbon Adsorber units. However, replacement of this material could also mitigate two other issues with Kombisorb BAT-37 (and its successor BAT II 37): (1) fire safety risk due to exothermic heat generated by adsorption reactions between the carbon material and Hg and (2) the risk of low retention of 129I in glass and subsequent downstream impacts on secondary liquid waste treatment at the Effluent Treatment Facility.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Optimal Control of SOEC-Based Hydrogen Production Systems for Demand Response Using Deep Reinforcement Learning in Smart Grids

Solid oxide electrolysis cell (SOEC) hydrogen production technology can range in size from small, appliance-size equipment to large-scale, central production facilities that can be tied directly to renewable or non-greenhouse-gas-emitting forms of electricity production, making it an ideal resource for demand response (DR). The SOEC hydrogen production system is a complex integrated system that encompasses fluid dynamics, electrical dynamics, and electrochemical and thermal dynamics, all of which involve non-linearity and non-convexity. Proper control of the SOEC hydrogen production system is crucial to enable its participation in the DR program. Here, to overcome the difficulty of designing an explicit control law for such nonlinear systems with nonconvex optimization features in DR applications, deep reinforcement learning (DRL) is explored to achieve the optimal control of the SOEC system for DR participation. Specifically, a twin delayed deterministic policy gradient (TD3) control framework is applied to achieve optimal response performance during DR events by considering power tracking error and hydrogen production efficiency with a suitable reward function. Two case studies with grid connections for tracking different DR commands were investigated. The first case study involved operating conditions reaching the boundaries, while the second involved operating conditions within the boundaries. The results showed that the proposed DRL-based control for SOEC can track the DR signal in a timely manner while maintaining high energy efficiency.

08 HYDROGEN↗
DOE OSTI2025

Physics-informed latent neural operator for real-time predictions of time-dependent parametric PDEs

Deep operator network (DeepONet) has shown significant promise as surrogate models for systems governed by partial differential equations (PDEs), enabling accurate mappings between infinite-dimensional function spaces. However, when applied to systems with high-dimensional input-output mappings arising from large numbers of spatial and temporal collocation points, these models often require heavily overparameterized networks, leading to long training times. Latent DeepONet addresses some of these challenges by introducing a two-step approach: first learning a reduced latent space using a separate model, followed by operator learning within this latent space. While efficient, this method is inherently data-driven and lacks mechanisms for incorporating physical laws, limiting its robustness and generalizability in data-scarce settings. Here, in this work, we propose PI-Latent-NO, a physics-informed latent neural operator framework that integrates governing physics directly into the learning process. Our architecture features two coupled DeepONets trained end-to-end: a Latent-DeepONet that learns a low-dimensional representation of the solution, and a Reconstruction-DeepONet that maps this latent representation back to the physical space. By embedding PDE constraints into the training via automatic differentiation, our method eliminates the need for labeled training data and ensures physics-consistent predictions. The proposed framework is both memory and compute-efficient, exhibiting near-constant scaling with problem size and demonstrating significant speedups over traditional physics-informed operator models. We validate our approach on a range of parametric PDEs, showcasing its accuracy, scalability, and suitability for real-time prediction in complex physical systems.

Latent representations↗
DOE OSTI2022

Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review

This review considers the topological fermion condensation quantum phase transition (FCQPT) that explains the complex behavior of strongly correlated Fermi systems, such as frustrated insulators with quantum spin liquid and heavy fermion metals. The review contrasts theoretical consideration with recent experimental data collected on both heavy fermion metals (HF) and frustrated insulators. Such a method allows to understand experimental data. We also consider experimental data collected on quantum spin liquid in Lu 3 Cu 2 Sb 3 O 14 and quasi-one dimensional (1D) quantum spin liquid in both YbAlO 3 and Cu(C 4 H 4 N 2 )(NO 3 ) 2 with the aim to establish a sound theoretical explanation for the observed scaling laws, Landau Fermi liquid (LFL) and non-Fermi-liquid (NFL) behavior exhibited by these frustrated insulators. The recent experimental data on the heavy-fermion metal α-YbAl 1-x Fe x B 4 , with x=0.014, and on its sister compounds β-YbAlB 4 and YbCo 2 Ge 4 , carried out under the application of magnetic field as a control parameter are analyzed. We show that the thermodynamic and transport properties as well as the empirical scaling laws follow from the fermion condensation theory. We explain how both the similarity and the difference in the thermodynamic and transport properties of α-YbAl 1-x Fe x B 4 and in its sister compounds β-YbAlB 4 and YbCo 2 Ge 4 emerge, as well as establish connection of these (HF) metals with insulators Lu 3 Cu 2 Sb 3 O 14 , Cu(C 4 H 4 N 2 )(NO 3 ) 2 and YbAlO 3 . We demonstrate that the universal LFL and NFL behavior emerge because the HF compounds and the frustrated insulators are located near the topological FCQPT or are driven by the application of magnetic fields.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-driven prediction of the output composition of an atmospheric pressure plasma jet

Cold atmospheric plasma (CAP) in open air hosts numerous chemical species engaged in thousands of chemical reactions. Comprehensive diagnosis of its chemical composition is important across various fields from medicine, where reactive oxygen and nitrogen play key roles, to surface modification. In applications, a centimeter-scale helium–air jet operates for minutes, featuring micrometer-sized streamers and an atmospheric pressure-induced collision frequency in the hundreds of GHz range. To address this intricate multi-scale issue, we introduce a machine learning approach: using a physics-informed neural network (PINN) to tackle the multi-scale complexities inherent in predicting the complete list of species concentrations, gas temperature, and electron temperature of a CAP jet supplied with a mixture of helium and air. Experimental measurements of O 3 , N 2 O, and NO 2 concentrations downstream of the plasma jet, combined with fundamental physics laws, the conservation of mass and charge, constrain the PINN, enabling it to predict the concentrations of all species that are not available from the experiment, along with gas and electron temperatures. The results, therefore, obey all the physical laws we provided and can have a chemical balance with the measured concentrations. In conclusion, this methodology holds promise for describing and potentially regulating complex systems with limited experimental datasets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2022

Scalable algorithms for physics-informed neural and graph networks

Physics-informed machine learning (PIML) has emerged as a promising new approach for simulating complex physical and biological systems that are governed by complex multiscale processes for which some data are also available. In some instances, the objective is to discover part of the hidden physics from the available data, and PIML has been shown to be particularly effective for such problems for which conventional methods may fail. Unlike commercial machine learning where training of deep neural networks requires big data, in PIML big data are not available. Instead, we can train such networks from additional information obtained by employing the physical laws and evaluating them at random points in the space–time domain. Such PIML integrates multimodality and multifidelity data with mathematical models, and implements them using neural networks or graph networks. Here, we review some of the prevailing trends in embedding physics into machine learning, using physics-informed neural networks (PINNs) based primarily on feed-forward neural networks and automatic differentiation. For more complex systems or systems of systems and unstructured data, graph neural networks (GNNs) present some distinct advantages, and here we review how physics-informed learning can be accomplished with GNNs based on graph exterior calculus to construct differential operators; we refer to these architectures as physics-informed graph networks (PIGNs). We present representative examples for both forward and inverse problems and discuss what advances are needed to scale up PINNs, PIGNs and more broadly GNNs for large-scale engineering problems.

42 ENGINEERING↗
DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗
DOE OSTI2025

Numerical simulation of coupled THM behaviour of full-scale EBS in backfilled experimental gallery in the Horonobe URL

Bentonite-based engineered barrier system (EBS) is a key component of many repository designs for the geological disposal of high-level radioactive waste. Given the complexity and interaction of the phenomena affecting the barrier system, coupled thermo-hydro-mechanical (THM) numerical analyses are a potentially useful tool for a better understanding of their behaviour. In this context, a Task (the Horonobe EBS experiment) was undertaken to study, using numerical analyses, the thermo-hydro-mechanical (and thermo-hydro) interactions in bentonite based engineered barriers within the international cooperative project DECOVALEX 2023. One full-scale in-situ experiment and four laboratory experiments, largely complementary, were selected for modelling. The Horonobe EBS experiment is a temperature-controlled non-isothermal experiment combined with artificial groundwater injection. The Horonobe EBS experiment consists of the heating and cooling phases. Six research teams performed the THM or TH (depended on research team approach) numerical analyses using a variety of computer codes, formulations and constitutive laws. Finally, for each experiment, the basic features of the analyses are described and the comparison between calculations and laboratory experiments and field observations are presented and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Flue Gas CO 2 Capture via Electrochemically Mediated Amine Regeneration: Desorption Unit Design and Analysis

We analyzed the performance of an electrochemical desorption unit to be implemented in an electrochemically mediated amine regeneration (EMAR) process, which separates CO 2 and regenerates the amine upon modulation of a targeted species that complexes with the amine molecules. Using a simplified one-dimensional boundary layer model, we combined the kinetic expressions of the faradaic reactions with mass transfer effects to analyze the polarization behavior of the flow cell. Mass transfer limitations become dominant at high current densities, leading to increases in energy consumption. An understanding of scaling laws with key parameters, such as the Sherwood number, is required to achieve high power operation with minimal mass transfer losses. Because the EMAR chemistry is sensitive to temperature, heat distribution in the unit is also examined. Our results indicate that the desorption performance increases with addition of a high-temperature heat source. Additionally, incorporation of membrane crossover suggests energy losses because of diffusion and migration of active species. Furthermore, these findings inform on scaling laws for at-scale desorber design and provide a model for scaling up the EMAR system for larger scale operations.

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

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Blazar variability power spectra from radio up to TeV photon energies: Mrk 421 and PKS 2155−304

ABSTRACT We present the results of the power spectral density (PSD) analysis for the blazars Mrk 421 and PKS 2155−304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability time-scales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from Owens Valley Radio Observatory, optical and infrared (B, V, R, I, J, H, and Kbands), X-rays from the Swift and the Rossi X-ray Timing Explorer, high and very high energy (VHE) γ-rays from the Fermi and Very Energetic Radiation Imaging Telescope Array System as well as the High Energy Stereoscopic System. Our results are: (1) the power-law form of the variability power spectra at radio, infrared, and optical frequencies have slopes ∼1.8, indicative of random-walk-type noise processes; (2) the power-law form of the variability power spectra at higher frequencies, from X-rays to VHE γ-rays, however, have slopes ∼1.2, suggesting a flicker noise-type process; and (3) there is significantly more variability power at X-rays, high and VHE γ-rays on time-scales ≲ 100 d, as compared to lower energies. Our results do not easily fit into a simple model, in which a single compact emission zone is dominating the radiative output of the blazars across all the time-scales probed in our analysis. Instead, we argue that the frequency-dependent shape of the variability power spectra points out a more complex picture, with highly inhomogeneous outflow producing non-thermal emission over an extended, stratified volume.

Goyal, Arti↗
DOE OSTI2023

Quadruple Langmuir probe characterization of different fuel gases in a plasma deflagration accelerator

Astrophysical flows may be studied by reproducing similar conditions using a coaxial plasma accelerator operating in the deflagration regime (or plasma deflagration accelerator). This allows for the recreation and investigation of dynamics present in complex highly coupled plasma systems at the laboratory scale. Here we report on measurements of the plasma density, temperature, plasma potential and velocity found using a quadruple Langmuir probe (QLP) on such a deflagration accelerator in the form of the Stanford Coaxial High ENerGy (CHENG) device operating with multiple gases – specifically argon, nitrogen and hydrogen. Experiments show a general decrease in bulk plasma velocity with gas atomic mass from upwards of 120 km s –1 with hydrogen to less than 30 km s –1 with argon. There was an accompanying increase in peak plasma density with increasing atomic mass from ~3 × 10 20 m –3 with hydrogen to ~1.5 × 10 21 m –3 with argon. It was found that the momentum flux and internal energy density also generally increase with atomic mass while the particle flux is constant between shots. Further investigation is needed to understand these correlations and the underlying physics. Lastly, comparisons with scaling laws show that while the CHENG device may be operated in such a way as to simulate the effects of bulk solar wind movement, it may not properly capture the thermal effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hydrologic connectivity and dynamics of solute transport in a mountain stream: Insights from a long-term tracer test and multiscale transport modeling informed by machine learning

The movement of solutes in a watershed is a complex process with multiple interactions and feedbacks across spatial and temporal scales. Modeling the dynamics of solute transport along diverse hydrologic pathways within watersheds – from hillslopes to stream channels and in and out of the hyporheic zones – is challenging but critically important, as these processes integrate and contribute to the biogeochemical functioning of the river corridor up to the river network scale. Here we use results from a long-term network-scale tracer test at the H.J. Andrews experimental forest in western Cascade Mountains, Oregon, USA to inform a multiscale framework for transport in stream corridors. The framework uses a Lagrangian-based subgrid model to represent the effects of hyporheic exchange flow and advective transport at stream network scales. The spatially and temporally resolved stream discharge needed for the transport model is imputed across the river system by an entity-aware long short-term memory network. Modeled concentrations show good agreements with the observations and exhibit power scaling laws indicative of a very wide range of timescales over which hyporheic exchange flow occurs. Our results demonstrate a data-informed modeling framework that links dynamical processes occurring at small scales to a network context to help understand how changes at reach scale cascade into network-scale effects, providing a useful tool for sustainable river basin management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Discerning the thermodynamic feasibility of the spontaneous coexistence of multiple functional vegetation groups

Can the Second Law of Thermodynamics explain why ecosystems naturally organize into a complex structure composed of multiple vegetation species and functional groups? Ecosystem structure, which refers to the number and type of plant functional groups, is the result of self-organization, or the spontaneous emergence of order from random fluctuations. By considering ecosystems as open thermodynamic systems, we model and study these fluctuations of throughput signatures on short timescales to determine the drivers and characteristics of ecosystem structure. This diagnostic approach allows us to use fluxes of energy and entropy to calculate an ecosystem’s estimated work and understand the thermodynamic behavior of the system. We use a multi-layer canopy-root-soil model to calculate the energy and entropy fluxes of different scenarios for field sites across various climates. At each site, scenarios comprised of native individual plant functional groups and a coexisting multi-group composition scenario including all functional groups observed at the site are compared. Ecosystem-scale calculations demonstrate that entropy fluxes and work efficiency—the work performed for the amount of radiation entering the ecosystem—are greatest in the multi-group scenario when its leaf area is significantly larger than each of its individual functional groups. Thus, we conclude that ecosystems self-organize towards the vegetation structure with the greatest outgoing entropy flux and work efficiency, resulting in the coexistence of multiple functional groups and performing the maximum amount of work within the constraints of locally available energy, water, and nutrients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗