Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential theory”

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 19 records

Stochastic Distribution Control Theory-Its Potential Application in Risk Management in Financial Systems

Stochastic Distribution Control (SDC) theory [1], originated by the author in 1996, aims at developing modeling and control strategies for dynamic and non-Gaussian stochastic systems by controlling the shape of the probability density functions of some concerned variables and parameters in stochastic systems. It generalizes the capability of standard stochastic differential equations and can therefore be applied to generic non-Gaussian systems. Since it was established in 1996, it has found a wide spectrum of applications in non-Gaussian stochastic system control, data mining, filtering and optimization for uncertain systems. In this short opinion article, discussions will be made on potential applications of SDC theory to financial systems in terms of risk analysis and management.

97 MATHEMATICS AND COMPUTING↗

Survey on stochastic distribution systems: A full probability density function control theory with potential applications

Complex systems seen either in general engineering practice or economics are subjected to ever increased uncertainties that are mostly represented as random variables or parameters, and the characteristics of random variables are represented by their probability density functions (PDFs). Controlling their PDFs means to shape their stochastic distributions and in general it would provide a full treatment for system analysis and operational control and optimization. This leads to the development of stochastic distribution control (SDC) systems theory in the past decades, where the original aim of the controller design is to realize a shape control of the distributions of certain random variables in their PDFs sense for some engineering processes. Indeed, once the PDFs of these random variables or parameters are used to describe their distribution characters, the control task is to obtain control signals so that the output PDFs of stochastic systems are made to follow their target PDFs. The subject of SDC was initially originated for non-Gaussian stochastic control systems design but has found a wide spectrum of applications in general systems in terms of data-driven modeling, analysis, signal processing (filtering), data mining via multivariable statistics, decision-making (optimization) for systems subjected to uncertainties and even in economics. In this context, SDC constitutes an effective primer tool for complex system analysis, control and operational optimizations. In this review paper, a detailed survey of the developments on the research of SDC systems will be made together with their wide spectrum applications and future perspectives.

42 ENGINEERING↗

Maximizing efficiency of dataset compression for machine learning potentials with information theory

Machine learning interatomic potentials (MLIPs) balance high accuracy and lower costs compared to density functional theory calculations, but their performance often depends on the size and diversity of training datasets. Large datasets improve model accuracy and generalization but are computationally expensive to produce and train on, while smaller datasets risk discarding rare but important atomic environments and compromising MLIP accuracy/reliability. Here, we develop an information-theoretical framework to quantify the efficiency of dataset compression methods and propose an algorithm that maximizes this efficiency. By framing atomistic dataset compression as an instance of the minimum set cover (MSC) problem over atom-centered environments, our method identifies the smallest subset of structures that contains as much information as possible from the original dataset while pruning redundant information. The approach is extensively demonstrated on the GAP-20 and TM23 datasets and validated on 64 varied datasets from the ColabFit repository. Across all cases, MSC consistently retains outliers, preserves dataset diversity, and reproduces the long-tail distributions of forces even at high compression rates, outperforming other subsampling methods. Furthermore, MLIPs trained on MSC-compressed datasets exhibit reduced error for out-of-distribution data even in low-data regimes. We explain these results using an outlier analysis and show that such quantitative conclusions could not be achieved with conventional dimensionality reduction methods. The algorithm is implemented in the open-source QUESTS package and can be used for several tasks in atomistic modeling, from data subsampling, outlier detection, and training improved MLIPs at a lower cost.

36 MATERIALS SCIENCE↗

Influence of strain and dislocations on GaSb/GaAs quantum dots: From nested to staggered band alignment

We investigate the influence of strain and dislocations on band alignment in GaSb/GaAs quantum dot systems. Composition profiles from cross-sectional scanning tunneling microscopy images are interpolated onto a finite element mesh in order to calculate the distribution of local elastic strain, which is converted to a spatially varying band alignment using deformation potential theory. Our calculations predict that dislocation-induced strain relaxation and charging lead to significant local variations in band alignment. Furthermore, misfit strain induces a transition from a nested (type I) to a staggered (type II) band alignment. Although dislocation-induced strain relaxation prevents the type I to type II transition, electrostatic charging at dislocations induces the staggered band alignment once again.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmark calculations of pure neutron matter with realistic nucleon-nucleon interactions

Here, we report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner–Bethe–Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne v 18 interaction and two of its simplified versions, and four of the new Norfolk Δ -full chiral effective field theory potentials. Primarily because of the advancements in the auxiliary-field diffusion Monte Carlo, we observe good agreement among the three many-body techniques up to nuclear saturation density—the maximum difference in the energy per particle is within 1.5 MeV for all the potentials we consider. At higher densities, the divergences become more important, and are mainly driven by the Fermi hypernetted chain/single-operator calculations. We also study the connection between nucleon-nucleon scattering data and the energy per particle of pure neutron matter. Our results suggest that fitting to higher-energy nucleon-nucleon scattering helps reduce the spread of energies among the models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Unified View of Topological Phase Transition in Band Theory

We develop a unified view of topological phase transitions (TPTs) in solids by revising the classical band theory with the inclusion of topology. Reevaluating the band evolution from an “atomic crystal” (a normal insulator (NI)) to a solid crystal, such as a semiconductor, we demonstrate that there exists ubiquitously an intermediate phase of topological insulator (TI), whose critical transition point displays a linear scaling between electron hopping potential and average bond length, underlined by deformation-potential theory. The validity of the scaling relation is verified in various two-dimensional (2D) lattices regardless of lattice symmetry, periodicity, and form of electron hoppings, based on a generic tight-binding model. Significantly, this linear scaling is shown to set an upper bound for the degree of structural disorder to destroy the topological order in a crystalline solid, as exemplified by formation of vacancies and thermal disorder. Our work formulates a simple framework for understanding the physical nature of TPTs with significant implications in practical applications of topological materials.

42 ENGINEERING↗

Transport and mechanical behavior in PEO-LLZO composite electrolytes

Abstract Composite solid electrolytes (CEs), wherein ion-conducting polymer and ceramic/glass is mixed, are promising candidates for all-solid-state batteries due to their promise of acceptable ionic conductivity and mechanical properties compared to their individual constituents. While numerous studies have focused on improving the performance of CEs, it is still unclear what the material targets are that can result in improved macroscopic performance especially in light of the coupled needs for high transport and high mechanical strength in these materials. In this study, a two-dimensional (2D) mathematical model is developed to investigate electrochemical and mechanical characteristics of CEs. The model is compared to CEs consisting of poly-ethylene-oxide (PEO) polymer and lithium lanthanum zirconium oxide (LLZO) ceramic material with examination of the impact of varying LLZO volume fractions. The potential drop at the PEO-LLZO interface is evaluated using the junction potential theory. Using experimental data from the literature, the model estimates the ionic conductivity, effective transference number, and mechanical stiffness of the CEs. While the mechanical stiffness improves with increasing volume fraction of LLZO, the impact on conductivity and transference number depends on interfacial resistance at the interface. Finally, the study reports CE’s potential to enhance Li-ion transport and mechanical properties to inhibit lithium (Li) dendrite growth.

25 ENERGY STORAGE↗

Membrane Degradation in PEM Fuel Cells: Part I. Modeling Gas Crossover and the Pt Band

Understanding chemical degradation of the proton-exchange membrane in fuel cells is crucial for extending their lifetimes. Herein, various degradation reactions reported in literature are organized and analyzed, including direct radical generation and an indirect (Fenton) pathway. To understand the transport of dissolved H 2 and O 2 crossover gases as they relate to membrane degradation, an agglomerate-scale model is introduced, treating gas, ionomer, and catalyst as discrete phases. The model reveals a key phenomenon: at working potentials, dissolved gases are mostly consumed at the interface between the catalyst layer and the membrane, leaving little gas to cross the membrane. Under open-circuit conditions, dissolved gases are not consumed and can then cross the membrane. This explains high H 2 O 2 concentrations and degradation rates seen in experiments but not captured in previous models. Following mixed-potential theory, crossover gases supply the hydrogen-oxidation and oxygen-reduction reactions, which occur simultaneously on individual Pt particles comprising the Pt band in the membrane, forming reactive species (H 2 O 2 , OH·). Results show crossover gas almost entirely reacts on the Pt band, allowing little to reach the opposite electrode. Furthermore, the micro-scale geometry of the catalyst-layer/membrane interface impacts the gas crossover at working potentials, indicating that cell construction affects membrane durability.

Johnson, Evan F. [Lawrence Berkeley National Labor↗

Comparison of integral equations used to study ${T}_{cc}^{+}$ for a stable D *

We perform a detailed comparison between three formalisms used in recent studies of DD* scattering at heavier-than-physical pion masses, which aim to understand the properties of the doubly-charmed tetraquark, ${T}_{cc}^{+}$ (3875). These methods are the three-particle relativistic field theory (RFT) formalism, the two-body Lippmann-Schwinger (LS) equation with chiral effective field theory potentials, and the two-particle relativistic framework proposed by Baião Raposo and Hansen (BRH approach). In a simplified single-channel setting, we derive the conditions under which the infinite-volume integral equations from the RFT and BRH approaches reduce to the LS form. We present numerical examples showing that differences between these methods can be largely removed by adjusting short-range couplings. We also address a number of technical issues in the RFT approach.

Hadronic Spectroscopy↗

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗

Performance estimation of resonance-enhanced dual-buoy wave energy converter using coupled time-domain simulation

This article presents the modeling methodology and performance evaluation of the resonance-enhanced dual-buoy WEC (Wave Energy Converter) by HEM (hydrodynamic & electro-magnetic) fully-coupled-dynamics time-domain-simulation program. The numerical results are systematically compared with the authors’ 1/6-scale experiment. With a direct-drive linear generator, the WEC consists of dual floating cylinders and a moon-pool between the cylinders, which can utilize three resonance phenomena from moon-pool dynamics as well as heave motions of inner and outer buoys. The contact and friction between the two buoys observed in the experiment are also properly modeled in the time-domain simulation by the Coulomb-friction model. Moon-pool resonance peaks significantly exaggerated in linear potential theory are empirically adjusted through comparisons with measured values. A systematic comparative study between the simulations and experiments with and without PTO (power-take-off) is conducted, and the relative heave displacements/velocities and power outputs are well matched. Then, parametric studies are carried out with the simulation program to determine optimum generator parameters. The performance with various wave conditions is also assessed.

16 TIDAL AND WAVE POWER↗

Annual Performance of the Second-Generation Variable-Geometry Oscillating Surge Wave Energy Converter

Recent studies in wave energy have highlighted the need for a structured innovation approach in wave energy converter (WEC) design because cost-of-energy estimates have remained high. One specific innovation being investigated by the National Renewable Energy Laboratory is WEC geometry control, which uses control surfaces in combination with an oscillating surge WEC (OSWEC) to increase device availability and power generation while limiting structural costs. This study performs the first analysis of annual performance for a novel OSWEC with geometry control to understand how the geometry control affects availability, annual power generation, and structural loadings like the power-take-off (PTO) torque or surge foundation force. Device hydrodynamic coefficients are calculated using linear potential theory for six geometry configurations. A nonideal PTO system is assumed and quadratic viscous damping losses are considered. Annual performance is evaluated and compared for three U.S. wave energy sites. The WEC geometry and PTO system are controlled on a sea-state basis to optimize for power capture while remaining under limits set on motion amplitude and structural loads. Results indicate that geometry control can increase availability up to 25 days in an average year depending on design limits, increase average power generation, and significantly reduce peak structural loads.

30 DIRECT ENERGY CONVERSION↗

Machine learning nonequilibrium electron forces for spin dynamics of itinerant magnets

Abstract We present a generalized potential theory for conservative as well as nonconservative forces for the Landau-Lifshitz magnetization dynamics. Importantly, this formulation makes possible an elegant generalization of the Behler-Parrinello machine learning (ML) approach, which is a cornerstone of ML-based quantum molecular dynamics methods, to the modeling of force fields in adiabatic spin dynamics of out-of-equilibrium itinerant magnetic systems. We demonstrate our approach by developing a deep-learning neural network that successfully learns the electron-mediated exchange fields in a driven s-d model computed from the nonequilibrium Green’s function method. We show that dynamical simulations with forces predicted from the neural network accurately reproduce the voltage-driven domain-wall propagation. Our work also lays the foundation for ML modeling of spin transfer torques and opens a avenue for ML-based multi-scale modeling of nonequilibrium dynamical phenomena in itinerant magnets and spintronics.

36 MATERIALS SCIENCE↗

Numerical simulation of parametric resonance in point absorbers using a simplified model

Abstract Parametric resonance is a non‐linear phenomenon in which a system can oscillate at a frequency different from its exciting frequency. Some wave energy converters are prone to this phenomenon, which is usually detrimental to their performance. Here, a computationally efficient way of simulating parametric resonance in point absorbers is presented. The model is based on linear potential theory, so the wave forces are evaluated at the mean position of the body. However, the first‐order variation of the body's centres of gravity and buoyancy is taken into account. This gives essentially the same result as a more rigorous approach of keeping terms in the equation of motion up to second order in the body motions. The only difference from a linear model is the presence of non‐zero off‐diagonal elements in the mass matrix. The model is benchmarked against state‐of‐the‐art non‐linear Froude–Krylov and computational fluid dynamics models for free decay, regular wave, and focused wave group cases. It is shown that the simplified model is able to simulate parametric resonance in pitch to a reasonable accuracy even though no non‐linear wave forces are included. The simulation speed on a standard computer is up to two orders of magnitude faster than real time.

16 TIDAL AND WAVE POWER↗

Reaching the equilibrium state of the frustrated triangular Ising magnet Ca 3 Co 2 O 6

Ca 3 Co 2 O 6 is a frustrated magnet consisting of a triangular arrangement of chains of Ising spins. It shows regular magnetization steps versus magnetic field every 1.2 T that are metastable with very slow dynamics. This has puzzled the community for many years and given rise to numerous potential theories. Here we approach the problem by seeking the elusive magnetic equilibrium state at T=2 K. To this end, we explore two approaches: (1) bypassing the slow dynamics produced by changing fields by instead field-cooling directly to the target temperature and (2) quantum annealing in transverse magnetic fields. While we observe no measurable effect of the quantum annealing in fields up to 7 T, which is likely due to the large Ising anisotropy of Co spins in this material, we find that for the field cooling in longitudinal fields we achieve the predicted equilibrium 1/3 magnetization. We perform Monte Carlo simulations of the ground state phase diagram and we also simulate the quantum annealing process and find good agreement between experiment and theory. Thus we present an investigation of the elusive ground state properties of the canonical frustrated triangular system Ca 3 Co 2 O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Strong increase in ultrasound attenuation below T c in Sr 2 RuO 4 : Possible evidence for domains

Recent experiments suggest that Sr 2 RuO 4 has a two-component superconducting order parameter (OP). A two-component OP has multiple degrees of freedom in the superconducting state that can result in low-energy collective modes or the formation of domain walls—a possibility that would explain a number of experimental observations including the smallness of the signature of time reversal symmetry breaking at T c and telegraph noise in critical current experiments. We use resonant ultrasound spectroscopy to perform ultrasound attenuation measurements across the superconducting T c of Sr 2 RuO 4 . We find that compressional sound attenuation increases by a factor of 7 immediately below T c , in sharp contrast with what is found in both conventional (s-wave) and high-T c (d-wave) superconductors. Our observations are most consistent with the presence of domain walls that separate different configurations of the superconducting OP. Here, the fact that we only observe an increase in sound attenuation for compressional strains, and not for shear strains, suggests an inhomogeneous superconducting state formed of two distinct, accidentally degenerate superconducting OPs that are not related to each other by symmetry. Whatever the mechanism, a factor of 7 increase in sound attenuation is a singular characteristic that must be reconciled with any potential theory of superconductivity in Sr 2 RuO 4 .

36 MATERIALS SCIENCE↗

Exploring NaCl-PuCl 3 molten salts with machine learning interatomic potentials and graph theory

Actinide molten salts are the basis of the liquid fuels used in molten salt reactors. Due to the inherent difficulties associated with high temperature and hazardous conditions, experimental investigations of fundamental properties of these materials are usually challenging. In this work, we describe the structure and transport of NaCl-PuCl 3 mixtures using computational techniques. Three compositions were considered (16, 25, and 36 mol% PuCl 3 ) over a temperature range (730 – 1257K) using ab initio molecular dynamics, which provided the necessary data sets for training machine learned interatomic potentials. Further, molecular dynamics simulations based on these potentials were then used to determine structure and transport properties. A substantial change was noted in the structure factor when increasing the PuCl 3 content from 25 to 36 mol%. This change is linked to the aggregation of larger Pu 3+ clusters. In addition, the similarity of the atomic environments of metal cations in molten salt systems to their solid states counterparts was investigated using an unsupervised learning technique. Finally, graph theory was employed to explore the structure and size of actinide networks. Consistent with the structure factor, a dense Pu 3+ intermolecular structure is observed within the 36 mol% PuCl 3 mixture. The structure of cation-cation inter-junctions is also discussed. In all cases, the diffusion of Pu 3+ is significantly lower than that of Na + and Cl - .

36 MATERIALS SCIENCE↗

Redox Chemistry of the Subphases of α-CsPbI2Br and β-CsPbI2Br: Theory Reveals New Potential for Photostability

The logic in the design of a halide-mixed APb(I 1-x Br x ) 3 perovskite is quite straightforward: to combine the superior photovoltaic qualities of iodine-based perovskites with the increased stability of bromine-based perovskites. However, even small amounts of Br doped into the iodine-based materials leads to some instability. In the present report, using first-principles computations, we analyzed a wide variety of α-CsPbI 2 Br and β-CsPbI 2 Br phases, compared their mixing enthalpies, explored their oxidative properties, and calculated their hole-coupled and hole-free charged Frenkel defect (CFD) formations by considering all possible channels of oxidation. Nanoinclusions of bromine-rich phases in α-CsPbI 2 Br were shown to destabilize the material by inducing lattice strain, making it more susceptible to oxidation. The uniformly mixed phase of α-CsPbI 2 Br was shown to be highly susceptible towards a phase transformation into β-CsPbI 2 Br when halide interstitial or halide vacancy defects were introduced into the lattice. The rotation of PbI 4 Br 2 octahedra in α-CsPbI2Br allows it either to transform into a highly unstable apical β-CsPbI 2 Br, which may phase-segregate and is susceptible to CFD, or to phase-transform into equatorial β-CsPbI 2 Br, which is resilient against the deleterious effects of hole oxidation (energies of oxidation >0 eV) and demixing (energy of mixing <0 eV). Thus, the selective preparation of equatorial β-CsPbI 2 Br offers an opportunity to obtain a mixed perovskite material with enhanced photostability and an intermediate bandgap between its constituent perovskites.

36 MATERIALS SCIENCE↗