Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systems 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quasi-Lindblad pseudomode theory for open quantum systems

Here, we introduce a new framework to study the dynamics of open quantum systems with linearly coupled Gaussian baths. Our approach replaces the continuous bath with an auxiliary discrete set of pseudomodes with dissipative dynamics, but we further relax the complete positivity requirement in the Lindblad master equation and formulate a quasi-Lindblad pseudomode theory. We show that this quasi-Lindblad pseudomode formulation directly leads to a representation of the bath correlation function in terms of a complex weighted sum of complex exponentials, an expansion that is known to be rapidly convergent in practice and thus leads to a compact set of pseudomodes. The pseudomode representation is not unique and can differ by a gauge choice. When the global dynamics can be simulated exactly, the system dynamics is unique and independent of the specific pseudomode representation. However, the gauge choice may affect the stability of the global dynamics, and we provide an analysis of why and when the global dynamics can retain stability despite losing positivity. We showcase the performance of this formulation across various spectral densities in both bosonic and fermionic problems, finding significant improvements over conventional pseudomode formulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quarkonium Polarization Kinetic Equation from Open Quantum Systems and Effective Field Theories

Recent measurements of polarization phenomena in relativistic heavy ion collisions have aroused a great interest in understanding dynamical spin evolution of the QCD matter. In particular, the spin alignment signature of J/ψ has been recently observed in Pb-Pb collisions at LHC, which may infer nontrivial spin transport of quarkonia in quark gluon plasmas. Motivated by this, we study the spin-dependent in-medium dynamics of quarkonia by using the potential nonrelativistic QCD (pNRQCD) and the open quantum system framework. By applying the Markovian approximation and Wigner transformation, we systematically derive the Boltzmann transport equation for vector quarkonia with polarization dependence in the quantum optical limit. As opposed to the previous study for the spin-independent case where the collision terms depend on chromoelectric correlators, the new kinetic equation incorporates gauge invariant correlators of chromomagnetic fields that determine the recombination and dissociation terms with polarization dependence at the order we are working in the multipole expansion. In the quantum Brownian motion limit, the Lindblad equation with new transport coefficients defined in terms of the chromomagnetic field correlators have also been derived. Our formalism is generic and valid for both weakly-coupled and strongly-coupled quark gluon plasmas. It may be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Theory of Two-Level Tunneling Systems in Superconductors

We develop a field theory formulation for the interaction of an ensemble of two-level tunneling systems (TLSs) with the electronic states of a superconductor. Predictions for the impact of two-level tunneling systems on superconductivity are presented, including T c and the spectrum of quasiparticle states for conventional BCS superconductors. We show that nonmagnetic TLS impurities in conventional s-wave superconductors can act as pair-breaking or pair-enhancing defects depending on the level population of the distribution of TLS impurities. We present calculations of the enhancement of superconductivity, both T c and the order parameter, for TLS defects in thermal equilibrium with the electrons and lattice. The scattering of quasiparticles by TLS impurities leads to subgap states below the bulk excitation gap, Δ, as well as resonances in the continuum above Δ. The energies and spectral weights of these states depend on the distribution of tunnel splittings, while the spectral weights are particularly sensitive to the level occupation of the TLS impurities. Under microwave excitation, or decoupling from the thermal bath, a nonequilibrium level population of the TLS distribution generates subgap quasiparticle states near the Fermi level that contribute to dissipation and thus degrade the performance of superconducting devices at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Two-dimensional non-equilibrium melting of charged colloids

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

critical phenomena↗
DOE OSTI2025

The mechanics of nucleation and growth and the surface tensions of active matter

Homogeneous nucleation, a textbook transition path for phase transitions, is typically understood on thermodynamic grounds through the prism of classical nucleation theory. However, recent studies have suggested the applicability of classical nucleation theory to systems far from equilibrium. Here, in this article, we formulate a purely mechanical perspective of homogeneous nucleation and growth, elucidating the criteria for the properties of a critical nucleus without appealing to equilibrium notions. Applying this theory to active fluids undergoing motility-induced phase separation, we find that nucleation proceeds in a qualitatively similar fashion to equilibrium systems, with concepts such as the Gibbs–Thomson effect and nucleation barriers remaining valid. We further demonstrate that the recovery of such concepts allows us to extend classical theories of nucleation rates and coarsening dynamics to active systems upon using the mechanically derived definitions of the nucleation barrier and surface tensions. Three distinct surface tensions—the mechanical, capillary, and Ostwald tensions—play a central role in our theory. While these three surface tensions are identical in equilibrium, our work highlights the distinctive role of each tension in the stability of active interfaces and the nucleation and growth of motility-induced phases.

Capillary waves↗
DOE OSTI2026

Global anomalies of Green's function zeros.

We study global anomalies of nonlocal effective theories proposed to describe symmetry-preserving Luttinger surfaces, i.e., the momentum-space manifolds of Green’s function zeros (GFZs) at zero energy, in strongly interacting fermionic systems. In particular, we focus on the simplest possible cases associated with a gapless Dirac zero, which is the counterpart of the gapless Dirac quasiparticle in weakly interacting systems. These theories may be derived by integrating out low-energy degrees of freedom that do not couple to the relevant gauge field. We discuss the global anomaly, the bulk-boundary correspondence, and the constraint on phases consistent with the anomaly, such as non-Fermi liquids and emergent gapless quasiparticles on Luttinger surfaces. Failing to avoid spontaneous symmetry breaking in the thermodynamical limit inevitably leads to unstable GFZs. We also provide some perspective on why the related nonlocal fermionic effective theory studied recently is not a suitable starting point for a symmetrically gapped phase

Su, Lei↗
DOE OSTI2025

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High-temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Theory of capillary tension and interfacial dynamics of motility-induced phases

The statistical mechanics of equilibrium interfaces has been well-established for over a half century. In the past decade, a wealth of observations have made increasingly clear that a new perspective is required to describe interfaces arbitrarily far from equilibrium. In this work, beginning from microscopic particle dynamics that break time-reversal symmetry, we derive the linear interfacial dynamics of coexisting motility-induced phases. Doing so allows us to identify the athermal energy scale that excites interfacial fluctuations and the nonequilibrium surface tension that resists these excitations. Our theory identifies that, in contrast to equilibrium fluids, this active surface tension contains contributions arising from nonconservative forces which act to suppress interfacial fluctuations and, crucially, is distinct from the mechanical surface tension of Kirkwood and Buff. Here we find that the interfacial stiffness scales linearly with the intrinsic persistence length of the constituent active particle trajectories, in agreement with simulation data. We demonstrate that at wavelengths much larger than the persistence length, the interface obeys surface-area minimizing Boltzmann statistics with our derived nonequilibrium interfacial stiffness playing a role identical to that of equilibrium systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Turbulent Heating in Collisionless Low-beta Plasmas: Imbalance, Landau Damping, and Electron–Ion Energy Partition

An understanding of how turbulent energy is partitioned between ions and electrons in weakly collisional plasmas is crucial for modeling many astrophysical systems. Using theory and simulations of a four-dimensional reduced model of low-beta gyrokinetics (the “Kinetic Reduced Electron Heating Model”), we investigate the dependence of collisionless heating processes on plasma beta and imbalance (normalized cross-helicity). These parameters are important because they control the helicity barrier, the formation of which divides the parameter space into two distinct regimes with remarkably different properties. In the first, at lower beta and/or imbalance, the absence of a helicity barrier allows the cascade of injected power to proceed to small (perpendicular) scales, but its slow cascade rate makes it susceptible to significant electron Landau damping, in some cases leading to a marked steepening of the magnetic spectra on scales above the ion Larmor radius. In the second, at higher beta and/or imbalance, the helicity barrier halts the cascade, confining electron Landau damping to scales above the steep “transition-range” spectral break, resulting in dominant ion heating. We formulate quantitative models of these processes that compare well to simulations in each regime, and combine them with results of previous studies to construct a simple formula for the electron–ion heating ratio as a function of beta and imbalance. This model predicts a “winner takes all” picture of low-beta plasma heating, where a small change in the fluctuations' properties at large scales (the imbalance) can cause a sudden switch between electron and ion heating.

Interplanetary turbulence↗
DOE OSTI2024

Response Functions of Correlated Systems within Green's Function Theory

Why use Green's functions as the fundamental variable? Wave-function (Psi) methods are king for high-fidelity and Density-functional (Rho) methods are very efficient (Kohn-Sham). Goldilocks principle: Green's function (G) methods straddle the Rho and Psi methods, intermediate in both accuracy and efficiency. Also, when interest lies in excitations & 2-particle properties: G-methods are natural - intrinsic to the theory.

DMFT↗
DOE OSTI2025

Assessing the Performance of Exchange‐Correlation Density Functionals in Describing the Iron‐Catalyzed Ammonia Synthesis System

Density functional theory (DFT) has been widely employed for elucidating mechanistic aspects of heterogeneous catalysis. However, the accuracy of DFT calculations relies heavily on selecting an exchange-correlation (XC) functional that correctly describes the electronic structure of materials involved in the reactions. This study assesses the accuracy of several XC density functionals for modeling the iron-catalyzed ammonia synthesis system. In the assessment of functional accuracy, experimental references are compared to DFT-calculated values for the formation energy of gas-phase ammonia and nitrogen, bulk Fe/Fe-nitride (γ′-Fe 4 N) lattice constants and cohesive/formation energies, and nitrogen and ammonia binding energies on Fe(100), Fe(111), Fe(110), and γ′-Fe 4 N(111). It is observed that the experimental value for each of these descriptors is accurately modeled by at least one functional. RPBE alone provides reliable estimates for both the lattice constant and cohesive energy of Fe and γ′-Fe 4 N. Temperature-programmed desorption experiments led to estimates for N and NH 3 adsorption across several Fe-based facets that are best captured by RPBE. These results highlight the importance of choosing an appropriate XC functional that accurately describes experimental systems and offer insights into effectively modeling the interactions between nitrogen and ammonia on Fe-based surfaces.

adsorption↗
DOE OSTI2025

Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems

Coupled-cluster theory (CC) is a highly accurate and versatile method for simulating complex interactions within quantum systems. The extension of CC theory to model mixed electron-photon processes with quantum electrodynamics (QED) has improved our capability to predict cavity-modified chemistry, a field where photons are used as cost-effective and eco-friendly alternatives to catalyze/inhibit chemical reactions. However, calculations with CC methods, even without incorporating QED effects, are often prohibitively expensive. Simulations of larger systems require scalable infrastructures that exist for traditional CC methods but not for QED-CC methods. As such, we present a GPU-enabled, high-performance, open-source implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) within the ExaChem quantum chemistry software package. ExaChem relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure: a parallel heterogeneous tensor library designed to achieve scalable performance on modern heterogeneous supercomputing platforms. Furthermore, we discuss theoretical foundations, algorithmic details, and numerical benchmarks to showcase the larger systems that ExaChem can simulate and how the integration of photonic degrees-of-freedom alters their ground-state properties.

Basis sets↗
DOE OSTI2025

Spacetime pq theory for AC and DC electric power systems

The 50/60 Hz alternating current (AC) electric power has been the standard and most flexible energy source powering our modern societies for one and a half centuries since the war of the currents: AC versus direct current (DC). A reactive power concept that was introduced at the beginning of the AC power was very useful for circuit/system analysis, design, control, optimization, and ultimately for more efficient and stable generation, transmission, distribution, and consumption. The initial reactive power theory was based on single-phase sinusoidal AC power to capture inductive and capacitive power that yields to net-zero average power over one fundamental cycle. Soon it was expanded to non-sinusoidal AC power and finally to instantaneous three-phase AC power. However, these reactive power theories remain separate and limited to special cases and have never been consolidated and made valid to all cases. Today, more widespread adoption of power electronics and renewable energy is bringing back DC power into the electric grids. The reactive power concept has never been applied to DC power systems. There is no reactive power in DC power systems according to the existing reactive power theories. Do DC power systems really have no reactive power? Capacitors and inductors are widely used in DC just like in AC power systems. Are they not reactive power components? Why are they different from their AC counterparts? Furthermore, are batteries active or reactive power components? What about active devices like power converters (or inverters) with AC (or DC) on one side and DC (or AC) on the other? Do they generate or consume reactive power? Finally, what about AC and DC hybrid power systems? How to define reactive power in such a complex power system that has a multitude of loads, buses, and sources? Is there reactive power between any two loads, any two buses, or any two sources in a power system and what is the total reactive power in such a complex power system as a whole? As the motivation and goal of this paper to answer the above basic questions, to unify the existing AC reactive power theories and to ultimately provide theoretical and insightful guidance for system analysis, design, control, efficiency, optimization, and operation of complex power systems, a concept of spacetime (both spatial and temporal) active and reactive power (pq) theory—the spatiotemporal aspect of active and reactive power—is developed for both AC and DC power systems. The theoretical definitions and physical meanings of the spacetime reactive power will be developed, and real applications and thought experiments/cases/exercises will be explored and discussed. The developed mathematics to define the active (or real) and reactive (or imaginary) power— p and q respectively by dot (scalar) and cross (vector) products of multi-dimension spacetime vectors and time-space mapping principle/law can have some fundamental implications as well.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2024

Analysis and Mitigation of Cascading Failures Using a Stochastic Interaction Graph with Eigen-analysis

In studies on complex network systems using graph theory, eigen-analysis is typically performed on an undirected graph model of the network. However, when analyzing cascading failures in a power system, the interactions among failures suggest the need for a directed graph beyond the topology of the power system to model directions of failure propagation. To accurately quantify failure interactions for effective mitigation strategies, this paper proposes a stochastic interaction graph model and associated eigen-analysis. Different types of modes on failure propagations are defined and characterized by the eigenvalues of a stochastic interaction matrix, whose absolute values are unity, zero, or in between. Finding and interpreting these modes helps identify the probable patterns of failure propagation, either local or widespread, and the participating components based on eigenvectors. Then, by lowering the failure probabilities of critical components highly participating in a mode of widespread failures, cascading can be mitigated. Here, the validity of the proposed stochastic interaction graph model, eigen-analysis and the resulting mitigation strategies is demonstrated using simulated cascading failure data on an NPCC 140-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗