Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theory - methods”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Entanglement entropy of a color flux tube in (2+1)D Yang-Mills theory

We construct a novel flux tube entanglement entropy (FTE 2 ), defined as the excess entanglement entropy relative to the vacuum of a region of color flux stretching between a heavy quark-anti-quark pair in pure gauge Yang-Mills theory. We show that FTE 2 can be expressed in terms of correlators of Polyakov loops, is manifestly gauge-invariant, and therefore free of the ambiguities in computations of the entanglement entropy in gauge theories related to the choice of the center algebra. Employing the replica trick, we compute FTE 2 for SU(2) Yang-Mills theory in (2+1)D and demonstrate that it is finite in the continuum limit. We explore the properties of FTE 2 for a half-slab geometry, which allows us to vary the width and location of the slab, and the extent to which the slab cross-cuts the color flux tube. Following the intuition provided by computations of FTE 2 in (1+1)D, and in a thin string model, we examine the extent to which our FTE 2 results can be interpreted as the sum of an internal color entropy and a vibrational entropy corresponding to the transverse excitations of the string.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Theory of strain phase equilibria and diagrams

Chemical phase equilibria and phase diagrams are well established and have served as indispensable guides for designing materials through chemical tuning. However, a general thermodynamic theory for strain phase equilibria remains elusive despite extensive studies which have revealed dramatic impacts of mechanical strain on the stability of phases and domain states in solid-state materials. Here, we establish the thermodynamic theory of strain equilibria and a general framework for efficiently constructing multidomain and multiphase diagrams of arbitrarily strained solids under incoherent conditions. As examples, we obtain temperature-strain phase diagrams of ferroelectric PbTiO 3 , strongly correlated VO 2 , and unconventional ferroelectric Hf 0.5 Zr 0.5 O 2 . We reveal the analogs of the Gibbs phase rules, multiple- and multi-critical points, common-tangent construction, and level rule in strain equilibria to those in the familiar chemical phase equilibria. Finally, our strain equilibria theory offers a powerful framework for predicting the thermodynamic stability of structural phases and domains in strained solids and for guiding the strain engineering of their functional properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fractonic gauge theory of smectics

Motivated by striped correlated quantum matter, and the recently developed duality between elasticity of a two-dimensional (2D) crystal and a gauge theory, we derive a dual coupled U(1) vector gauge theory for a two-dimensional (2D) quantum smectic, where the disclination is mapped onto the fractonic charge, that we demonstrate can only move transversely to smectic layers. This smectic gauge theory dual also emerges from a gauge dual of a quantum crystal through a Higgs transition corresponding to a single flavor of its dipole condensation, an anisotropic quantum melting via dislocation proliferation. A condensation of the second flavor of dislocations a corresponds to another Higgs transition describing the smectic-to-nematic melting. We also utilize the electrostatic limit of this duality to formulate a melting of a 2D classical smectic in terms of a higher derivative sine–Gordon model, demonstrating its instability to a nematic at any nonzero temperature. Generalizing this classical duality to a 3D smectic, gives formulation of a 3D nematic-to-smectic transition in terms of an anisotropic Abelian-Higgs model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Information theory optimization of signals from small-angle scattering measurements

Small-angle X-ray scattering (SAXS) of particles in solution informs on the conformational states and assemblies of biological macromolecules (bioSAXS) outside of cryo- and solid-state conditions. In bioSAXS, the SAXS measurement under dilute conditions is resolution limited, and through an inverse Fourier transform, the measured SAXS intensities directly relate to the physical space occupied by the particles via the P (r)-distribution. Yet, this inverse transform of SAXS data has been historically cast as an ill-posed, ill-conditioned problem requiring an indirect approach. Here, we show that through the applications of matrix and information theories, the inverse transform of SAXS intensity data is a well-conditioned problem. The so-called ill-conditioning of the inverse problem is directly related to the Shannon number. By exploiting the oversampling enabled by modern detectors, a direct inverse Fourier transform of the SAXS data is possible, provided the recovered information does not exceed the Shannon number. The Shannon limit corresponds to the maximum number of significant singular values that can be recovered in a SAXS experiment, suggesting this relationship is a fundamental property of band-limited inverse integral transform problems. This correspondence reduces the complexity of the inverse problem to the Shannon limit and maximum dimension. We propose a hybrid scoring function using an information theory framework that assesses both the quality of the model-data fit as well as the quality of the recovered P (r)-distribution. The hybrid score utilizes the Akaike information criteria and Durbin-Watson statistic that considers parameter-model complexity, i.e., degrees of freedom, and the randomness of the model-data residuals. The described tests and findings extend the boundaries for bioSAXS by completing the information theory formalism initiated by Peter B. Moore to enable a quantitative measure of resolution in SAXS, robustly determine maximum dimension, and more precisely define the best parameter model appropriately representing the observed scattering data.

Rambo, Robert P. [Science and Technology Facilitie↗
DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Theory of subcycle time-resolved photoemission: Application to terahertz photodressing in graphene

Here, motivated by recent experimental progress we revisit the theory of pump–probe time- and angle-resolved photoemission spectroscopy (trARPES), which is one of the most powerful techniques to trace transient pump-driven modifications of the electronic properties. The pump-induced dynamics can be described in different gauges for the light–matter interaction. Standard minimal coupling leads to the velocity gauge, defined by linear coupling to the vector potential. In the context of tight-binding (TB) models, the Peierls substitution is the commonly employed scheme for single-band models. Multi-orbital extensions – including the coupling of the dipole moments to the electric field – have been introduced and tested recently. In this work, we derive the theory of time-resolved photoemission within both gauges from the perspective of nonequilibrium Green’s functions. This approach naturally incorporates the photoelectron continuum, which allows for a direct calculation of the observable photocurrent. Following this route we introduce gauge-invariant expressions for the time-resolved photoemission signal. The theory is applied to graphene pumped with short terahertz pulses, which we treat within a first-principles TB model. We investigate the gauge invariance and discuss typical effects observed in subcycle time-resolved photoemission. Our formalism is an ideal starting point for realistic trARPES simulations including scattering effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Core-to-core X-ray emission spectra from Wannier based multiplet ligand field theory

Recent advances using Density Functional Theory (DFT) to augment Multiplet Ligand Field Theory (MLFT) have led to ab-initio calculations of many formerly empirical parameters. This development makes MLFT more predictive instead of interpretive, thus improving its value for studies of highly correlated 3 d , 4 d , and ƒ-electron systems. Synchrotron time is always at a premium, and tools that provide predictive capabilities have clear value when it comes to planning studies. Here, in this work, we develop a DFT + MLFT based approach for core-to-core Kα x-ray emission spectra (XES) and evaluate its performance for a range of transition metal systems. We find good agreement between theory and experiment, as well as the ability to capture key spectral trends related to spin and oxidation state. We also discuss limitations of the model in the context of the remaining free parameters and suggest directions forward.

Ab initio theoretical treatment↗
DOE OSTI2021

Aqueous ion partitioning in Nafion: Applicability of Manning's counter-ion condensation theory

Nafion is a well-studied polymer of significant interest as a fuel cell membrane due to rapid transport of protons, which is often ascribed, at least in part, to its phase-separated morphology. The literature detailing ion sorption and transport in Nafion diverges considerably from that of water purification membranes, with contrasting analyses and models. To help bridge these gaps and facilitate comparisons of membranes with various structures, this study presents measurements of ion activity coefficients and sorption in Nafion 117 equilibrated with aqueous solutions of NaCl, MgCl 2 , and Na 2 SO 4 , which are relevant to water purification. The experimental data are compared to a theoretical framework based on Manning's counter-ion condensation theory and Donnan theory. Good agreement between the model and data is observed for NaCl and MgCl 2 , using no adjustable parameters. However, the model significantly underpredicts Na 2 SO 4 solubility, which is attributed to ion pairing effects. Furthermore, the results highlight the general applicability of counter-ion condensation theory to highly charged polymers, including those exhibiting nanoscale ordering, while revealing fundamental limitations related to ion pairing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the connection between perturbation theory and new semiclassical expansion in quantum mechanics

Here, it is shown that for the one-dimensional anharmonic oscillator with potential V(x) = ax 2 + bgx 3 + ... = $\frac{1}{g^2}$ $\hat{V}$ (gx), as well as for the radial oscillator V(r) = $\frac{1}{g^2}$ $\hat{V}$ (gr) and for the perturbed Coulomb problem V(r) = $\frac{a}{r}$ + βgr + ... = g $\tilde{V}$ (gr), the Perturbation Theory in powers of the coupling constant g (weak coupling regime) and the semiclassical expansion in powers of $\hbar$ 1/2 for the energies coincide. This is related to the fact that the dynamics developed in two spaces: x (r)-space and gx (gr)-space, lead to the same energy spectra. The equations which govern dynamics in these two spaces, the Riccati-Bloch equation and the Generalized Bloch equation, respectively, are presented. It is shown that the perturbation theory for the logarithmic derivative of the wavefunction in - space leads to (true) semiclassical expansion in powers of $\hbar$ 1/2 ; for the one-dimensional case this corresponds to the flucton calculus for the density matrix in the path integral formalism in Euclidean (imaginary) time proposed by one of the authors. Matching the perturbation theory in powers of g and the semiclassical expansion in powers of $\hbar$ 1/2 for the wavefunction leads to a highly accurate local approximation in the entire coordinate space, its expectation value for the Hamiltonian provides a prescription for the summation of the perturbative (trans)-series.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Time-domain Brillouin scattering theory for probe light and acoustic beams propagating at an angle and acousto-optic interaction at material interfaces

A theory has been developed to interpret time-domain Brillouin scattering (TDBS) experiments involving coherent acoustic pulse (CAP) and light pulse beams propagating at an angle to each other. It predicts the influence of the directivity pattern of their acousto-optic interaction on TDBS signals when heterodyne detection of acoustically scattered light is in backward direction to incident light. The theory reveals relationships between the carrier frequency, amplitude and duration of acoustically induced ”wave packets” in light transient reflectivity signals, and factors such as CAP duration, widths of light and sound beams, and their interaction angle. It describes the transient dynamics of these wave packets when the light and CAP encounter material interfaces, and how the light scattering by the incident CAP transforms into scattering by the reflected and transmitted CAPs. The theory suggests that single-point TDBS experiments can determine not only depth positions of buried interfaces but also their local inclinations/orientations.

42 ENGINEERING↗
DOE OSTI2023

Large-time correlation functions in bosonic lattice field theories

Large-time correlation functions have a pivotal role in extracting particle masses from Euclidean lattice field theory calculations, however little is known about the statistical properties of these quantities. In this work, the asymptotic form of the distributions of the correlation functions at vanishing momentum is determined for bosonic interacting lattice field theories with a unique gapped vacuum. It is demonstrated that the deviations from the asymptotic form at large Euclidean times can be utilized to determine the spectrum of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Cascading from $\mathscr{N}$ = 2 supersymmetric Yang–Mills theory to confinement and chiral symmetry breaking in adjoint QCD

We argue that adjoint QCD in 3 + 1 dimensions, with any SU(N) gauge group and two Weyl fermion flavors (i.e. one adjoint Dirac fermion), confines and spontaneously breaks its chiral symmetries via the condensation of a fermion bilinear. We flow to this theory from pure $\mathscr{N}$ = 2 SUSY Yang–Mills theory with the same gauge group, by giving a SUSY-breaking mass M to the scalars in the $\mathscr{N}$ = 2 vector multiplet. This flow can be analyzed rigorously at small M, where it leads to a deconfined vacuum at the origin of the $\mathscr{N}$ = 2 Coulomb branch. The analysis can be extended to all M using an Abelian dual description that arises from the N multi-monopole points of the $\mathscr{N}$ = 2 theory. At each such point, there are N −1 hypermultiplet Higgs fields h$^{i=1,2}_m$, which are SU(2) R doublets. We provide a detailed study of the phase diagram as a function of M, by analyzing the semi-classical phases of the dual using a combination of analytic and numerical techniques. The result is a cascade of first-order phase transitions, along which the Higgs fields h i m successively turn on, and which interpolates between the Coulomb branch at small M, where all h$^{i}_m$ = 0, and a maximal Higgs branch, where all h$^{i}_m$ ≠ 0, at sufficiently large M. We show that this maximal Higgs branch precisely matches the confining and chiral symmetry breaking phase of two-flavor adjoint QCD, including its broken and unbroken symmetries, its massless spectrum, and the expected large-N scaling of various observables. The spontaneous breaking pattern SU(2) R → U(1) R , consistent with the Vafa–Witten theorem, is ensured by an intricate alignment mechanism for the h$^{i}_m$ in the dual, and leads to a CP 1 sigma model of increasing radius along the cascade.

D’Hoker, Eric [Univ. of California, Los Angeles, C↗
DOE OSTI2025

Quantum complexity in gravity, quantum field theory, and quantum information science

Quantum complexity quantifies the difficulty of preparing a state or implementing a unitary transformation with limited resources. Applications range from quantum computation to condensed matter physics and quantum gravity. Here, we seek to bridge the approaches of these fields, which define and study complexity using different frameworks and tools. We describe several definitions of complexity, along with their key properties. In quantum information theory, we focus on complexity growth in random quantum circuits. In quantum many-body systems and quantum field theory (QFT), we discuss a geometric definition of complexity in terms of geodesics on the unitary group. In dynamical systems, we explore a definition of complexity in terms of state or operator spreading, as well as concepts from tensor-networks. We also outline applications to simple quantum systems, quantum many-body models, and QFTs including conformal field theories (CFTs). Finally, we explain the proposed relationship between complexity and gravitational observables within the holographic anti-de Sitter (AdS)/CFT correspondence.

Baiguera, Stefano [Istituto Nazionale di Fisica Nu↗
DOE OSTI2024

Resolving Discrepancies Between Theory and Experiment for the NCN + H Reaction

Flame modeling studies have highlighted the role of branching in the kinetics for the prompt NO switch reaction, NCN + H. In this reaction, there is strong competition between the CH + N2 (R1) and HCN + N (R2) product channels. Increased branching towards reaction R2 promotes the subsequent formation of Fenimore or Prompt NO. Recent direct shock tube studies on this complex reaction conclude that branching to reaction R2 is much more favorable than the predicted branching from theoretical studies. The experimentally predicted prompt NO switch temperature (TS), i.e., the temperature at which k2/(k1+k2) = 0.5, is at 1670 K, in contrast to theory that predicts Ts > 3200 K. The present theory/modeling work was initiated to highlight potential causes for this significant discrepancy between experiment and theory. Our analysis of the simulations reported in the experimental studies indicate that including a fall-off representation for C2H5 radical dissociation (with C2H5 produced from C2H5I dissociation and acting as an H-atom source) has a noticeable influence on the simulations for NCN decay and HCN formation in the shock tube studies. With the inclusion of a theory-based pressure fall-off representation for C2H5 dissociation, we show that this radical persists for much longer timescales even at the high-temperatures in the shock tube studies. This persistence then facilitates the rapid conversion of reactive H-atoms to CH3 radicals at the shock tube conditions via the H + C2H5 → CH3 + CH3 reaction. An updated theoretical analysis for this classic addition-elimination reaction is also provided in this work. The rapid formation of CH3 radicals necessitates the inclusion of additional reactions to adequately simulate the NCN and HCN data. Our simulations conclude that, with the inclusion of these additional reactions, the NCN and HCN profiles can be reasonably well simulated by the most recent theoretical predictions.

Ethyl Radicals↗
DOE OSTI2022

Quantum Perturbation Theory Using Tensor Cores and a Deep Neural Network

In this work, time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e., dense matrix–matrix multiplications, in mixed-precision arithmetics, which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree–Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low-precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

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

Computation of Auger Electron Spectra in Organic Molecules with Multiconfiguration Pair-Density Functional Theory

Efficient and accurate computation of molecular Auger electron spectra for larger systems is limited by the rapid increase in the number of doubly ionized final states as the system size grows. Here, in this work, we benchmark the application of multiconfiguration pair-density functional theory with a restricted active space (RAS) reference wave function for computing the carbon K-edge decay spectra of 20 organic molecules. Decay rates are computed within the one-center approximation. We evaluate the performance of different basis sets and on-top functionals and find that multiconfiguration pair-density functional theory achieves accuracy comparable to RAS followed by second-order perturbation theory, but at significantly lower computational cost.

Fouda, Adam E. A. [Argonne National Laboratory (AN↗
DOE OSTI2023

Local Molecular Field Theory for Coulomb Interactions in Aqueous Solutions

Coulomb interactions play a crucial role in a wide array of processes in aqueous solutions but present conceptual and computational challenges to both theory and simulations. We review recent developments in an approach addressing these challenges–local molecular field (LMF) theory. LMF theory exploits an exact and physically suggestive separation of intermolecular Coulomb interactions into strong short-range and uniformly slowly varying long-range components. This allows us to accurately determine the averaged effects of the long-range components on the short-range structure using effective single particle fields and analytical corrections, greatly reducing the need for complex lattice summation techniques used in most standard approaches. Furthermore, the simplest use of these ideas in aqueous solutions leads to the short solvent (SS) model, where both solvent–solvent and solute–solvent Coulomb interactions have only short-range components. Here we use the SS model to give a simple description of pairing of nucleobases and biologically relevant ions in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗