Search NASA⌕ Search

SEARCH · Search NASA

Results for “operator exponential”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Matrix Multiply Performance of GPUs on Exascale-class HPE/Cray Systems

The computation of dense matrix-matrix products (GEMMs) is central to many modeling and simulation workloads as well as AI/ML deep learning campaigns. In fact, millions of dollars are spent annually on computing GEMMs, and large model training demands are increasing exponentially. Specialized processors such as GPUs are designed to perform well for these operations. However, the performance of GEMMs on GPUs can exhibit complex behaviors depending on many factors, making it challenging to optimize the performance of GEMMs on these processors. In this study we undertake an examination of GEMM performance on several leading GPU models taken from product lines of GPUs to be deployed in forthcoming exascale computing systems. We show results to illustrate the many factors that can affect performance of GEMMs on GPUs. We then present data collected from a large number of test runs for an example GEMM operation to show the dependence behaviors of GEMM rate on matrix dimensions. Finally, we show results from machine learning-based performance models using novel feature engineering methods to fit the measured performance, providing a potential basis for GEMM performance tuning and autotuning methods for GPUs. Recommendations are also given for how to achieve high GEMM performance on modern GPUs.

Melesse Vergara, Veronica↗
DOE OSTI2024

Universal bounds on CFT Distance Conjecture

For any unitary conformal field theory in two dimensions with the central charge c, we prove that, if there is a nontrivial primary operator whose conformal dimension ∆ vanishes in some limit on the conformal manifold, the Zamolodchikov distance t to the limit is infinite, the approach to this limit is exponential ∆ = exp(−αt + O(1)), and the decay rate obeys the universal bounds c−1/2 ≤ α ≤ 1. In the limit, we also find that an infinite tower of primary operators emerges without a gap above the vacuum and that the conformal field theory becomes locally a tensor product of a sigma-model in the large radius limit and a compact theory. As a corollary, we establish a part of the Distance Conjecture about gravitational theories in three-dimensional anti-de Sitter space. In particular, our bounds on α indicate that the emergence of exponentially light states is inevitable as the moduli field corresponding to t rolls beyond the Planck scale along the steepest path and that this phenomenon can begin already at the curvature scale of the bulk geometry. We also comment on implications of our bounds for gravity in asymptotically flat spacetime by taking the flat space limit and compare with the Sharpened Distance Conjecture.

AdS-CFT Correspondence↗
DOE OSTI2020

Fixed-point quantum circuits for quantum field theories

Renormalization group ideas and effective operators are introduced to efficiently prepare ground states of massive lattice field theories on digital quantum devices. This is accomplished with a systematic approximation through localized unitaries that removes an exponentially costly barrier in the spatial volume of the quantum simulation. With these methods, classically computed ground states in a spatial volume L, containing a few Compton wavelengths, can be used to determine operators for preparing the ground state toward the thermodynamic limit with a precision improving as e –mL on beyond-classical quantum registers. Here, due to the exponential spatial decay of correlations in massive theories and the double exponential suppression of digitization artifacts in the number of qubits representing the scalar field, the derived fixed-point quantum circuits are expected to be relevant for simulations of quantum field theories throughout the evolution from small-scale near-term quantum devices to large-scale fault-tolerant quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Development of data‐driven exponential integrators with application to modeling of delay photocurrents

Abstract Radiation‐induced photocurrent effects represent a threat to microelectronic components operating in space, manmade and terrestrial radiation environments. Analysis of these threats by circuit simulations requires accurate and computationally efficient compact models. Most existing compact models are based on closed form analytic solutions of the governing equations and require empirical assumptions and idealizations that can limit their validity. In this paper we formulate an alternative numerical, data‐driven approach that learns a compact model from data representative of the type of measurements one can obtain in an experimental facility. To develop the model we start from a generic discrete‐time dynamical system and then use physics knowledge to refine its structure. Numerical studies demonstrate the potential of the model and establish some empirical guidelines for its training.

Bochev, Pavel↗
DOE OSTI2020

Formulations and Valid Inequalities for Optimal Black Start Allocation in Power Systems

The restoration of a power system after a blackout starts around units with enhanced technical capabilities, referred to as black start units (BSUs). We examine the planning problem of optimally allocating these units on the grid subject to a budget constraint. We present a mixed integer programming model based on current literature in power systems. Binary variables are associated with the allocation of BSUs and with the energization state of buses, branches, and generators of the power system over a time horizon. We extract a substructure of the feasible region which imposes the requirement that each island that appears during the restoration process must have at least one operational generator. We discuss three equivalent reformulations for this requirement. We introduce a family of exponentially many, polynomially separable, valid inequalities to strengthen the formulation. Under simplifying assumptions, we show that the convex hull of the feasible region is a full-dimensional polyhedron, and prove that some of the constraints we introduced are facet-defining. We perform experiments to examine the difference in strength between the formulations as well as the computational times to solve the problem to near optimality for synthetic instances of the IEEE-39, IEEE-118, Illinois-200, WECC-225, IEEE-300, South Carolina-500, and Texas-2000 power systems. We illustrate a use case of the model. We conclude by suggesting extensions of the current work for future research.

Black Start Allocation↗
DOE OSTI2022

Learning Distributed Geometric Koopman Operator for Sparse Networked Dynamical Systems

Koopman operator theory provides an alternative to study nonlinear networked dynamical systems by mapping the state space to an abstract higher dimensional space where the system evolution is linear. Recent works show the application of graph neural networks (GNNs) to learn state to object-centric embeddings and achieve centralized block-wise computation of Koopman operator (KO) under additional assumptions on the underlying agents properties and constraints on the KO structure. However, the computational complexity of learning the Koopman increases exponentially for networked systems where the number of possible system states grows in a combinatorial fashion with the number of nodes. The learning challenge is further amplified for sparse networks by two factors: 1) sample sparsity for learning the Koopman operator in the non-linear space, and 2) the divergence in the dynamics of individual nodes or from one subgraph to another. Our work aims to address these challenge by formulating the representation learning of networked dynamical systems into a multi-agent paradigm and learning the Koopman operator in a distributive manner. The computational as well as performance advantages of distributed Koopman is predominant for sparse networks whereas for fully connected networks, it is shown to coincide with the centralized one. The empirical study on rope system, network of oscillators and a synthetic power system show comparable and superior performance along with computational benefits with the state-of-the-art methods.

Mukherjee, Sayak↗
DOE OSTI2020

Calculating the divided differences of the exponential function by addition and removal of inputs

We introduce a method for calculating the divided differences of the exponential function by means of addition and removal of items from the input list to the function. Our technique exploits a new identity related to divided differences recently derived by F. Zivcovich. We show that upon adding an item to or removing an item from the input list of an already evaluated exponential, the re-evaluation of the divided differences can be done with only O(sn) floating point operations and O(sn) bytes of memory, where [z 0 ,...,z n ] are the inputs and s ∝ max i,j |z i – z j |. We demonstrate our algorithm’s ability to deal with input lists that are orders-of-magnitude longer than the maximal capacities of the current state-of-the-art. Here, we discuss in detail one practical application of our method: the efficient calculation of weights in the off-diagonal series expansion quantum Monte Carlo algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Covariance operator estimation: Sparsity, lengthscale, and ensemble Kalman filters

This paper investigates covariance operator estimation via thresholding. For Gaussian random fields with approximately sparse covariance operators, we establish non-asymptotic bounds on the estimation error in terms of the sparsity level of the covariance and the expected supremum of the field. We prove that thresholded estimators enjoy an exponential improvement in sample complexity compared with the standard sample covariance estimator if the field has a small correlation lengthscale. As an application of the theory, we study thresholded estimation of covariance operators within ensemble Kalman filters.

Covariance operator estimation↗
DOE OSTI2024

Spectral scheme for atomic structure calculations in density functional theory

In this study, we present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number Ζ = 1 to 83.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Chaotic roots of the modular multiplication dynamical system in Shor's algorithm

Shor's factoring algorithm, believed to provide an exponential speedup over classical computation, relies on finding the period of an exactly periodic quantum modular multiplication operator. This exact periodicity is the hallmark of an integrable system, which is paradoxical from the viewpoint of quantum chaos, given that the classical limit of the modular multiplication operator is a highly chaotic system that occupies the “maximally random” Bernoulli level of the classical ergodic hierarchy. In this work, we approach this apparent paradox from a quantum dynamical systems viewpoint, and consider whether signatures of ergodicity and chaos may indeed be encoded in such an “integrable” quantization of a chaotic system. We show that Shor's modular multiplication operator, in specific cases, can be written as a superposition of quantized 𝐴-baker's maps exhibiting more typical signatures of quantum chaos and ergodicity. This work suggests that the integrability of Shor's modular multiplication operator may stem from the interference of other “chaotic” quantizations of the same family of maps, and paves the way for deeper studies on the interplay of integrability, ergodicity, and chaos in and via quantum algorithms.

Dynamical systems↗
DOE OSTI2022

Performance Analysis and Optimization for Scientific Data Workloads

Scientific data generated at experimental and observational facilities are increasingly being processed on large-scale compute systems. Most of the experimental data analysis workflows are not designed or implemented to run on large scale environments and take full advantage of HPC compute and storage resources. These applications are unlike the traditional tightly-coupled scientific applications and hence face significant performance and scalability challenges as the volume of data increases exponentially. In this paper, we conduct a performance and scalability analysis for experimental analysis applications and workflows operating on data from light sources. Our analysis detects and quantifies I/O performance, scalability and runtime bottlenecks for three data analysis applications that run on NERSC resources. Based on our analysis we propose and implement a set of optimizations that lead to reducing the amount of time spent on I/O operations by almost 90%.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nonperturbative quantum gravity in a closed Lorentzian universe

We study how meaningful physical predictions can arise in nonperturbative quantum gravity in a closed Lorentzian universe. In such settings, recent developments suggest that the quantum gravitational Hilbert space is one-dimensional and real for each α-sector, as induced by spacetime wormholes. This appears to obstruct the conventional quantum-mechanical prescription of assigning probabilities via projection onto a basis of states. While previous approaches have introduced external observers or augmented the theory to resolve this issue, we argue that quantum gravity itself contains all the necessary ingredients to make physical predictions. We demonstrate that the emergence of classical observables and probabilistic outcomes can be understood as a consequence of partial observability: physical observers access only a subsystem of the universe. Tracing out the inaccessible degrees of freedom yields reduced density matrices that encode classical information, with uncertainties exponentially suppressed by the environment’s entropy. We develop this perspective using both the Lorentzian path integral and operator formalisms and support it with a simple microscopic model. Our results show that quantum gravity in a closed universe naturally gives rise to meaningful, robust predictions without recourse to external constructs.

AdS-CFT Correspondence↗
DOE OSTI2024

Understanding poromechanical response of a biogenic coalbed methane reservoir

Abstract Biogenic coalbed methane (BCBM) reservoirs aim to produce methane from in situ coal deposits following microbial conversion of coal. Success of BCBM reservoirs requires economic methane production within an acceptable timeframe. The work reported here quantifies the findings of previously published qualitative work, where it was found that bioconversion induces strains in the pore, matrix and bulk scales. Using imaging and dynamic strain monitoring techniques, the bioconversion induced strain is quantified here. To understand the effect of these strains from a reservoir geomechanics perspective, a corresponding poromechanical model is developed. Furthermore, findings of imaging experiments are validated using core-flooding flow experiments. Finally, expected field-scale behavior of the permeability response of a BCBM operation is modeled and analyzed. The results of the study indicated that, for Illinois coals, bioconversion induced strains result in a decrease in fracture porosity, resulting in a detrimental permeability drop in excess of 60% during bioconversion, which festers itself exponentially throughout its producing life. Results indicate that reservoirs with high initial permeability that will support higher Darcian flowrates, would be better suited for coal bioconversion, thereby providing a site-selection criteria for BCBM operations.

Pandey, Rohit (ORCID:0009000415483436)↗
DOE OSTI2020

The soft quark Sudakov

There has been recent interest in understanding the all loop structure of the subleading power soft and collinear limits, with the goal of achieving a systematic resummation of subleading power infrared logarithms. Most of this work has focused on subleading power corrections to soft gluon emission, whose form is strongly constrained by symmetries. In this paper we initiate a study of the all loop structure of soft fermion emission. In N = 1 QCD we perform an operator based factorization and resummation of the associated infrared logarithms using the formalism introduced in, and prove that they exponentiate into a Sudakov due to their relation to soft gluon emission. We verify this result through explicit calculation to O($α$ $^{3}_{s}$). We show that in QCD, this simple Sudakov exponentiation is violated by endpoint contributions proportional to (C A –C F ) n which contribute at leading logarithmic order. Combining our N = 1 result and our calculation of the endpoint contributions to O($α$ $^{3}_{s}$), we conjecture a result for the soft quark Sudakov in QCD, a new all orders function first appearing at subleading power, and give evidence for its universality. Our result, which is expressed in terms of combinations of cusp anomalous dimensions in different color representations, takes an intriguingly simple form and also exhibits interesting similarities to results for large-x logarithms in the off diagonal splitting functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Celestial amplitudes from UV to IR

Celestial amplitudes represent 4D scattering of particles in boost, rather than the usual energy-momentum, eigenstates and hence are sensitive to both UV and IR physics. We show that known UV and IR properties of quantum gravity translate into powerful constraints on the analytic structure of celestial amplitudes. For example the soft UV behavior of quantum gravity is shown to imply that the exact four-particle scattering amplitude is meromorphic in the complex boost weight plane with poles confined to even integers on the negative real axis. Would-be poles on the positive real axis from UV asymptotics are shown to be erased by a flat space analog of the AdS resolution of the bulk point singularity. The residues of the poles on the negative axis are identified with operator coefficients in the IR effective action. Far along the real positive axis, the scattering is argued to grow exponentially according to the black hole area law. Exclusive amplitudes are shown to simply factorize into conformally hard and conformally soft factors. The soft factor contains all IR divergences and is given by a celestial current algebra correlator of Goldstone bosons from spontaneously broken asymptotic symmetries. The hard factor describes the scattering of hard particles together with the boost-eigenstate clouds of soft photons or gravitons required by asymptotic symmetries. These provide an IR safe S-matrix for the scattering of hard particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2021

Renewable Energy for Industrial Environmental Management

Costs for renewable energy technologies have declined rapidly in the past decade and their use for residential, commercial, and utility scale electricity has grown exponentially as they become cost competitive. Simultaneously, industrial and manufacturing processes have been increasingly seeking ways to reduce emissions and operational costs in highly competitive sectors. With these combined drivers of lower cost and reduced environmental impact, renewable energy may become a viable energy provider for industrial processes such as oil and gas, mining, chemical refining, food production, and manufacturing. Renewable energy technologies may also partner with other reduced emission energy sources, such as small modular nuclear reactors and carbon capture and utilization, to create cleaner and circular industrial systems for reduced resource use. The Joint Institute for Strategic Energy Analysis (JISEA), which is a partnership of the National Renewable Energy Laboratory (NREL) and five universities and others, has been studying to potential for application of clean energy technologies to the heterogenous energy demands in industry. Dr. Jill Engel-Cox will present an overview of NREL and JISEA, the status and potential future of renewable energy technologies, and collaborations with the oil and gas industry and other industrial sectors to improve their environmental performance and reduce operational costs.

ENERGY PLANNING, POLICY, AND ECONOMY↗