Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear operator 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Catalytic resonance theory for parametric uncertainty of programmable catalysis

Microkinetic models are useful tools for screening catalytic materials; however, errors in their input parameters can lead to significant uncertainty in model predictions of catalyst performance. Here, in this work, we investigate the impact of linear scaling and Brønsted-Evans-Polanyi relation parametric uncertainty on microkinetic predictions of programmable-catalyst performance. Two case studies are considered: a generic A-to-B prototype reaction and the oxygen evolution reaction (OER). The results show that error-unaware models can accurately predict trends and, for the prototype reaction, values of optimal waveform parameters. The specific model parameters driving output uncertainty are identified via variance-based global sensitivity analysis. However, predictions of dynamic rate enhancement can decrease when uncertainty is propagated into the models. In both cases, we identify operating conditions where the programmable catalyst achieves a rate enhancement of at least one order of magnitude despite parametric uncertainty in the model, supporting programmable catalysis as a viable strategy for exceeding the Sabatier limit.

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

Holographic thermal correlators: a tale of Fuchsian ODEs and integration contours

We analyze real-time thermal correlation functions of conserved currents in holographic field theories using the grSK geometry, which provides a contour prescription for their evaluation. We demonstrate its efficacy, arguing that there are situations involving components of conserved currents, or derivative interactions, where such a prescription is, in fact, essential. To this end, we first undertake a careful analysis of the linearized wave equations in AdS black hole backgrounds and identify the branch points of the solutions as a function of (complexified) frequency and momentum. All the equations we study are Fuchsian with only regular singular points that for the most part are associated with the geometric features of the background. Special features, e.g., the appearance of apparent singular points at the horizon, whence outgoing solutions end up being analytic, arise at higher codimension loci in parameter space. Using the grSK geometry, we demonstrate that these apparent singularities do not correspond to any interesting physical features in higher-point functions. We also argue that the Schwinger-Keldysh collapse and KMS conditions, implemented by the grSK geometry, continue to hold even in the presence of such singularities. For charged black holes above a critical charge, we furthermore demonstrate that the energy density operator does not possess an exponentially growing mode, associated with ‘pole-skipping’, from one such apparent singularity. Our analysis suggests that the connection between the scrambling physics of black holes and energy transport has, at best, a limited domain of validity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Predicting Critical Transitions in Multiscale Data

Predicting the dynamics of complex nonlinear systems remains a challenging problem both in dynamical systems theory as well as real world science and engineering applications. Data-driven methods utilizing the latest advances in machine learning (ML) provide a promising new paradigm for this task. Our work centered on Reservoir Computing (RC), which has shown itself to be capable of skillfully predicting chaotic dynamics in multiscale systems. In the first part of the work, the focus is on how to improve predictions of critical transitions in a class of slow-fast metastable systems in which the equations are known. An additional goal was to determine whether a relationship exists between RC and Koopman operator theory, to improve the efficiency and broaden the applicability of the approach. In the second part of this work, a variation on the RC model known as Reconstructive Reservoir Computing (RRC) is applied to real-world data to identify anomalies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Excitation spectra of quantum matter without quasiparticles. II. Random t - J models

Here, we present numerical solutions of the spectral functions of large dimension t -J models with random nearest neighbor exchange and global SU(M) spin rotation symmetry. The solutions are obtained from the saddle-point equations of the large dimension limit, followed by the large M limit. The same saddle point equations also apply to the model with all-to-all and random hopping and exchange. Such a t -J model realizes a deconfined critical state proposed as a model for the optimally doped cuprates. The large M theory involves Green’s functions for fractionalized spinons and holons carrying emergent U(1) gauge charges, obeying relations similar to those of the Sachdev-Ye-Kitaev (SYK) models. The low frequency spectral functions are compared with an analytic analysis of the operator scaling dimensions with good agreement. We also compute the low frequency and temperature behavior of experimentally observable gauge-invariant observables: the electron Green’s function, the local spin susceptibility and the optical conductivity, along with the temperature dependence of the d.c. resistivity. The time reparameterization soft mode (equivalent to the boundary graviton in holographically dual models of two-dimensional quantum gravity) makes important contributions to all observables and provides a linear-in-temperature contribution to the d.c. resistivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Derivative-based SINDy (DSINDy): Addressing the challenge of discovering governing equations from noisy data

Recent advances in the field of data-driven dynamics allow for the discovery of ODE systems using state measurements. One approach, known as Sparse Identification of Nonlinear Dynamics (SINDy), assumes the dynamics are sparse within a predetermined basis in the states and finds the expansion coefficients through linear regression with sparsity constraints. This approach requires an accurate estimation of the state time derivatives, which is not necessarily possible in the high-noise regime without additional constraints. We present an approach called Derivative-based SINDy (DSINDy) that combines two novel methods to improve ODE recovery at high-noise levels. First, we denoise the state variables by applying a projection operator that leverages the assumed basis for the system dynamics. Second, we use a second order cone program (SOCP) to find the derivative and governing equations simultaneously. We derive theoretical results for the projection-based denoising step, which allow us to estimate the values of hyperparameters used in the SOCP formulation. This underlying theory helps limit the number of required user-specified parameters. Finally, we present results demonstrating that our approach leads to improved system recovery for the Van der Pol oscillator, the Duffing oscillator, the Rössler attractor, and the Lorenz 96 model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A freely falling graviton in the D1D5 CFT

We study a freely falling graviton propagating in AdS in the context of the D1D5 CFT, where we introduce an interaction by turning on a deformation operator. We start with one left and right moving boson in the CFT. After applying two deformation operators, the initial bosons split into three left moving and three right moving bosons. We compute the amplitude for various energies and extrapolate the result to the large energy region. At early times, the amplitude is linear in time. This corresponds to an infalling graviton becoming redshifted in AdS. At late times, the amplitude is periodic, which agrees with the fact that a freely falling graviton will not be thermalized.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Experimental evaluation of shear modulus scaling of dynamic strength at extreme pressures

Recent progress in the development of dynamic strength experimental platforms is allowing for unprecedented insight into the assumptions used to construct constitutive models operating in extreme conditions. In this work, we make a quantitative assessment of how tantalum strength scales with its shear modulus to pressures of hundreds of gigapascals through a cross-platform examination of three dynamic strength experiments. Specifically, we make use of Split–Hopkinson pressure bar and Richtmyer–Meshkov instability experiments to assess the low-pressure strain and strain rate dependence. Concurrent examination of magnetically driven ramp-release experiments up to pressures of 350GPa allows us to examine the pressure dependence. In this work, using a modern description of the shear modulus, validated against both ab initio theory and experimental measurements, we then assess how the experimentally measured pressure dependence scales with shear modulus. Furthermore, we find that the common assumption of scaling strength linearly with the shear modulus is too soft at high pressures and offer discussion as to how descriptions of slip mediated plasticity could result in an alternative scaling that is consistent with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Entropic order parameters in weakly coupled gauge theories

The entropic order parameters measure in a universal geometric way the statistics of non-local operators responsible for generalized symmetries. In this article, we compute entropic order parameters in weakly coupled gauge theories. To perform this computation, the natural route of evaluating expectation values of physical (smeared) non-local operators is prevented by known difficulties in constructing suitable smeared Wilson loops. We circumvent this problem by studying the smeared non-local class operators in the enlarged non-gauge invariant Hilbert space. This provides a generic approach for smeared operators in gauge theories and explicit formulas at weak coupling. In this approach, the Wilson and ’t Hooft loops are labeled by the full weight and co-weight lattices respectively. We study generic Lie groups and discuss couplings with matter fields. Smeared magnetic operators, as opposed to the usual infinitely thin ones, have expectation values that approach one at weak coupling. The corresponding entropic order parameter saturates to its maximum topological value, except for an exponentially small correction, which we compute. On the other hand, smeared ’t Hooft loops and their entropic disorder parameter are exponentially small. We verify that both behaviors match the certainty relation for the relative entropies. In particular, we find upper and lower bounds (that differ by a factor of 2) for the exact coefficient of the linear perimeter law for thin loops at weak coupling. This coefficient is unphysical/non-universal for line operators. We end with some comments regarding the RG flows of entropic parameters through perturbative beta functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Catalytic Resonance Theory: The Catalytic Mechanics of Programmable Ratchets

Catalytic reaction networks of multiple elementary steps operating under dynamic conditions via a programmed input oscillation are difficult to interpret and optimize due to reaction system complexity. To understand these dynamic systems, individual elementary catalytic reactions oscillating between catalyst states were evaluated to identify their three fundamental characteristics that define their ability to promote reactions away from equilibrium. First, elementary catalytic reactions exhibit directionality to promote reactions forward or backward from equilibrium as determined by a ratchet directionality metric comprised of the input oscillation duty cycle and the reaction rate constants. Second, catalytic ratchets are defined by the catalyst state of strong or weak binding that permit reactants to proceed through the transition state. Third, elementary catalytic ratchets exhibit a cutoff frequency which defines the transition in applied frequency for which the catalytic ratchet functions to promote chemistry away from equilibrium. All three ratchet characteristics are calculated from chemical reaction parameters including rate constants derived from linear scaling parameters, reaction conditions, and catalyst state. The characteristics of the reaction network’s constituent elementary catalytic reactions provided an interpretation of complex reaction networks and a method of predicting the behavior of dynamic surface chemistry on oscillating catalysts.

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

Noise-induced barren plateaus in variational quantum algorithms

Abstract Variational Quantum Algorithms (VQAs) may be a path to quantum advantage on Noisy Intermediate-Scale Quantum (NISQ) computers. A natural question is whether noise on NISQ devices places fundamental limitations on VQA performance. We rigorously prove a serious limitation for noisy VQAs, in that the noise causes the training landscape to have a barren plateau (i.e., vanishing gradient). Specifically, for the local Pauli noise considered, we prove that the gradient vanishes exponentially in the number of qubits n if the depth of the ansatz grows linearly with n . These noise-induced barren plateaus (NIBPs) are conceptually different from noise-free barren plateaus, which are linked to random parameter initialization. Our result is formulated for a generic ansatz that includes as special cases the Quantum Alternating Operator Ansatz and the Unitary Coupled Cluster Ansatz, among others. For the former, our numerical heuristics demonstrate the NIBP phenomenon for a realistic hardware noise model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A parallel strategy for density functional theory computations on accelerated nodes

Using the Löwdin orthonormalization of tall-skinny matrices as a proxy-app for wavefunction-based Density Functional Theory solvers, we investigate a distributed memory parallel strategy focusing on Graphics Processing Unit (GPU)-accelerated nodes as available on some of the top ranked supercomputers at the present time. Here we present numerical results in the strong limit regime, as it is particularly relevant for First-Principles Molecular Dynamics. We also examine how matrix product-based iterative solvers provide a competitive alternative to dense eigensolvers on GPUs, allowing to push the strong scaling limit of these computations to a larger number of distributed tasks. Our strategy, which relies on replicated Gram matrices and efficient collective communications using the NCCL library, leads to a time-to-solution under 0.5 s for the Löwdin orthonormalization of a tall-skinny matrix of 3000 columns on Summit at Oak Ridge Leadership Facility (OLCF). Given the similarity in computational operations between one iteration of a DFT solver and this proxy-app, this shows the possibility of solving accurately the DFT equations well under a minute for 3000 electronic wave functions, and thus perform First-Principles molecular dynamics of physical systems much larger than traditionally solved on CPU systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Simplified calculation of the area specific impedance for solid-state battery design

Simplified algebraic area specific impedance (ASI) correlations have been developed for solid-state composite battery electrodes made of a single ion conducting electrolyte, conductive additive, and intercalation active material. Two ASI expressions were developed, one for short times ( i.e. , pulsed power operation) and another for the pseudo steady state operation (i.e., sustained discharge for energy estimation). A full electrochemical model based on porous electrode theory was developed to examine the accuracy of the simplified ASI expressions. The simplified expressions agree favorably with full model results over a wide range of parameters (i.e. , electrode thicknesses, electrolyte conductivities, solid-state diffusion coefficients, specific surface areas, etc.) and conditions (i.e. , C-rates, states of charge, and pulse times). Under most conditions, the error between the full model and the correlations is well below 7 %. Higher errors were observed for the pseudo steady state expression at high/ low states of charge where the assumption of uniform reaction distributions loses validity. Here, the short time ASI has higher error at low states of charge due to the nonlinearity of the open circuit voltage equation, which is assumed linear in the formulation of the simplified algebraic expression.

25 ENERGY STORAGE↗
DOE OSTI2026

Combined effective field theory interpretation of Higgs boson, electroweak vector boson, top quark, and multijet measurements

Constraints on Wilson coefficients (WCs) corresponding to dimension-6 operators of the standard model effective field theory (SMEFT) are determined from a simultaneous fit to seven sets of CMS measurements probing Higgs boson, electroweak vector boson, top quark, and multijet production. Measurements of electroweak precision observables are also included and provide complementary constraints to those from the CMS experiment. The CMS measurements, using LHC proton-proton collision data at $\sqrt{s}=13\,\text {Te}\text {V} $, corresponding to integrated luminosities of 36.3 or 138$\,\text {fb}^{-1}$, are chosen to provide sensitivity to a broad set of operators, for which consistent SMEFT predictions can be derived. These are primarily measurements of differential cross sections which are parameterized as functions of the WCs. In measurements targeting ${\text {t}} (\bar{\textrm{t}})\text {X} $ production, SMEFT effects are modelled at the detector level. Individual constraints on 64 WCs, and constraints on 43 linear combinations of WCs, are obtained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quantum space, quantum time, and relativistic quantum mechanics

We treat space and time as bona fide quantum degrees of freedom on an equal footing in Hilbert space. Motivated by considerations in quantum gravity, we focus on a paradigm dealing with linear, first-order Hamiltonian and momentum constraints that lead to emergent features of temporal and spatial translations. Unlike the conventional treatment, we show that Klein-Gordon and Dirac equations in relativistic quantum mechanics can be unified in our paradigm by applying relativistic dispersion relations to eigenvalues rather than treating them as operator-valued equations. With time and space being treated on an equal footing in Hilbert space, we show symmetry transformations to be implemented by unitary basis changes in Hilbert space, giving them a stronger quantum mechanical footing. Global symmetries, such as Lorentz transformations, modify the decomposition of Hilbert space; and local symmetries, such as U(1) gauge symmetry are diagonal in coordinate basis and do not alter the decomposition of Hilbert space. Here, we briefly discuss extensions of this paradigm to quantum field theory and quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Strategies for simulating the time evolution of Hamiltonian lattice field theories

Simulating the time evolution of quantum field theories given some Hamiltonian H requires developing algorithms for implementing the unitary operator e -iHt . A variety of techniques exist that accomplish this task, with the most common technique used so far being Trotterization, which is a special case of the application of a product formula. However, other techniques exist that promise better asymptotic scaling in certain parameters of the theory being simulated, the most efficient of which are based on the concept of block encoding. In this work we study the performance of such algorithms in simulating lattice field theories. We derive and compare the asymptotic gate complexities of several commonly used simulation techniques in application to Hamiltonian lattice field theories. Using the scalar $\hat{φ}$ 4 theory as a test, we also perform numerical studies and compare the gate costs required by product formulas and signal-processing-based techniques to simulate time evolution. For the latter, we use the linear combination of unitaries (LCU) construction augmented with the quantum Fourier transform circuit to switch between the field and momentum eigenbases, which leads to immediate order-of-magnitude improvement in the cost of preparing the block encoding. Further, this paper also includes a pedagogical review of the techniques used, in particular product formulas, LCU, qubitization, quantum signal processing, as well as the technique for simulating geometrically-local Hamiltonians developed by Haah, Hastings, Kothari, and Low.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A new and alternative look at nonlinear Alfvénic states

The formulation for studying nonlinear Alfvenic states, sustainable in Hall Magnetohydrodynamics (HMHD), becomes considerably simpler and more tractable when circularly polarized Beltrami vectors (the eigenstates of linear HMHD) are used as the basis functions. Nonlinear HMHD is, then, reduced to a rather simple looking set of scalar equations from which a model problem of three interacting Beltrami modes is formulated and analytically solved. The triplet interactions span a variety of familiar nonlinear processes leading to a redistribution as well as periodic exchange of energy. The energy exchange processes (whose strength is measured by an energy exchange/depletion time) will, perhaps, play a dominant role in determining the spectral content of an eventual Alfvenic state. All nonlinearities (sensitive functions of the interacting wave vectors) operate at par, and none is dominant over any substantial region of k-space; their intricate interplay prevents a “universal” picture from emerging; few generalizations on the processes that may, for instance, lead to a turbulent state, are possible. However, the theory can definitely claim: (1) the energy tends to flow from lower to higher k and (2) the higher k z (in the direction of the ambient magnetic field) components of a mode with a given k are depleted/oscillate faster—in some cases much faster. It is noteworthy that the mode coupling is the strongest (with the shortest depletion time) when the participating wave vectors are nearly perpendicular; perhaps, an expected consequence of the curl (cross product) nonlinearities. Numerical simulations will be necessary to help create a fully reliable picture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Cold ion beam in a storage ring as a platform for large-scale quantum computers and simulators: Challenges and directions for research and development

The purpose of this paper is to evaluate the possibility of constructing a large-scale storage-ring-type ion-trap system capable of storing, cooling, and controlling a large number of ions as a platform for scalable quantum computing (QC) and quantum simulations. In such a trap, the ions form a crystalline beam moving along a circular path with a constant velocity determined by the frequency and intensity of the cooling lasers. In this paper, we consider a large leap forward in terms of the number of ions that serve as qubits in QC, from fewer than 100 available in state of the art linear ion-trap devices today to an order of 10 5 crystallized ions in the storage-ring setup. This new trap design unifies two different concepts: the storage rings of charged particles and the linear ion traps used for QC and mass spectrometry. In this paper, we use the language of particle accelerators to discuss the ion state and dynamics. We outline the differences between the above concepts, analyze challenges of the large ring with a revolving chain of ions, and propose goals for the research and development required to enable future quantum computers with 1000 times more qubits than available today. The challenge of creating such a large-scale quantum system while maintaining the necessary coherence of the qubits and the high fidelity of quantum logic operations is significant. Performing analog quantum simulations may be an achievable initial goal for such a device. Quantum calculations and simulations of complex quantum systems will move forward both the fundamental science and the applied research. Nuclear and particle physics, many-body quantum systems, lattice gauge theories, and nuclear structure calculations are just a few examples in which a large-scale quantum simulation system will become a very powerful tool to move forward our understanding of nature.

36 MATERIALS SCIENCE↗