SEARCH · Search NASA

Results for “Matrix 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Computer simulation of immobilized pH gradients at acidic and alkaline extremes - A quest for extended pH intervals

Computer simulations of the concentration profiles of simple biprotic ampholytes with Delta pKs 1, 2, and 3, on immobilized pH gradients (IPG) at extreme pH values (pH 3-4 and pH 10-11) show markedly skewed steady-state profiles with increasing kurtosis at higher Delta pK values. Across neutrality, all the peaks are symmetric irrespective of their Delta pK values, but they show very high contribution to the conductivity of the background gel and significant alteration of the local buffering capacity. The problems of skewness, due to the exponential conductivity profiles at low and high pHs, and of gel burning due to a strong electroosmotic flow generated by the net charges in the gel matrix, also at low and high pHs, are solved by incorporating in the IPG gel a strong viscosity gradient. This is generated by a gradient of linear polyacrylamide which is trapped in the gel by the polymerization process.

Mosher, Richard A.↗

DOE OSTI2026

Inverse problem in the large momentum effective theory framework

One proposal to compute parton distributions from first principles is the large momentum effective theory (LaMET), which requires the Fourier transform of matrix elements computed nonperturbatively. Lattice quantum chromodynamics (QCD) provides calculations of these matrix elements over a finite range of Fourier harmonics that are often noisy or unreliable in the largest computed harmonics. It has been suggested that enforcing an exponential decay of the missing harmonics helps alleviate this issue. Using nonperturbative data, we show that the uncertainty introduced by this inverse problem in a realistic setup remains significant without very restrictive assumptions, and that the importance of the exact asymptotic behavior is minimal for values of 𝑥 where the framework is currently applicable. We show that the crux of the inverse problem lies in harmonics of the order of 𝜆 = 𝑧⁢𝑃 𝑧 ∼ 5–15, where the signal in the lattice data is often barely existent in current studies, and the asymptotic behavior is not firmly established. We stress the need for more sophisticated techniques to account for this inverse problem, whether in the LaMET or related frameworks like the short-distance factorization. We also address a misconception that, with available lattice methods, the LaMET framework allows a “direct” computation of the 𝑥-dependence, whereas the alternative short-distance factorization only gives access to moments or fits of the 𝑥-dependence.

Dutrieux, Hervé [Aix-Marseille Université, Marseil↗

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

A knowledge gap analysis for transient CHF prediction within RELAP5-3D

Here we compare several experimental pool and flow boiling critical heat flux experiments with corresponding computational models developed with RELAP5-3D. Encompassed within this paper are several sensitivity analyses for the boiling heat transfer correlations, critical heat flux (CHF) prediction within RELAP5-3D and the thermal properties of the cladding for experiments in the Transient Reactor Test Loop (TRTL) and Transient Reactor Test (TREAT) facilities. The sensitivity analyses were consistent in determining that there is a large contribution from the critical heat flux multiplier which in this application, accounts for the effects of power transient critical heat flux. Comparison to experimental results of both the TREAT CHF-SERTTA experimental campaign and a University of New Mexico (UNM) pool boiling experimental campaign demonstrate both the under-prediction of transient CHF and nucleate boiling heat transfer. To address the gaps in knowledge regarding the transient heating effect for both boiling heat transfer and CHF prediction, we present an experimental test matrix for the TRTL facility that will isolate the power transient effect under PWR relevant operational conditions. The test matrix will incorporate 6 power pulses with a constant energy deposition and a timescale that spans over several orders of magnitude, including a steady state test with an exponential period of 10 s.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Gapless spin liquid and pair density wave of the Hubbard model on three-leg triangular cylinders

We study the ground state properties of the Hubbard model on three-leg triangular cylinders using large-scale density-matrix renormalization group simulations. At half-filling, we identify an intermediate gapless spin liquid phase, which has one gapless spin mode and algebraic spin–spin correlations but exponential decay scalar chiral–chiral correlations, between a metallic phase at weak coupling and Mott insulating dimer phase at strong interaction. Upon light doping the gapless spin liquid, the system exhibits power-law charge-density-wave (CDW) correlations but short-range single-particle, spin–spin, and chiral–chiral correlations. Similar to CDW correlations, the superconducting correlations also decay in power-law but oscillate in sign as a function of distance, which is consistent with the striped pair-density wave. When further doping the gapless spin liquid phase or doping the dimer order phase, another phase takes over, which has similar CDW correlations but all other correlations decay exponentially.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Quantum chaos on edge

Recently, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum [the Sachdev-Ye-Kitaev (SYK) model being a renowned representative] as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of “sparse” and the near edge of “dense” chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and invariant random matrix ensembles or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss. Published by the American Physical Society 2024

Altland, Alexander (ORCID:0000000229914805)↗

DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗

NASA NTRS1994

Control theory and splines, applied to signature storage

In this report the problem we are going to study is the interpolation of a set of points in the plane with the use of control theory. We will discover how different systems generate different kinds of splines, cubic and exponential, and investigate the effect that the different systems have on the tracking problems. Actually we will see that the important parameters will be the two eigenvalues of the control matrix.

Enqvist, Per↗

NASA NTRS1998

Attitude/attitude-rate estimation from GPS differential phase measurements using integrated-rate parameters

A sequential filtering algorithm is presented for attitude and attitude-rate estimation from Global Positioning System (GPS) differential carrier phase measurements. A third-order, minimal-parameter method for solving the attitude matrix kinematic equation is used to parameterize the filter's state, which renders the resulting estimator computationally efficient. Borrowing from tracking theory concepts, the angular acceleration is modeled as an exponentially autocorrelated stochastic process, thus avoiding the use of the uncertain spacecraft dynamic model. The new formulation facilitates the use of aiding vector observations in a unified filtering algorithm, which can enhance the method's robustness and accuracy. Numerical examples are used to demonstrate the performance of the method.

Oshman, Yaakov↗

NASA NTRS1997

Vector Observation-Aided/Attitude-Rate Estimation Using Global Positioning System Signals

Oshman, Yaakov↗

DOE OSTI2021

Island for gravitationally prepared state and pseudo entanglement wedge

We consider spacetime initiated by a finite-sized initial boundary as a generalization of the Hartle-Hawking no-boundary state. We study entanglement entropy of matter state prepared by such spacetime. We find that the entanglement entropy for large subregion is given either by the initial state entanglement or the entanglement island, preventing the entropy to grow arbitrarily large. Consequently, the entanglement entropy is always bounded from above by the boundary area of the island, leading to an entropy bound in terms of the island. The island I is located in the analytically continued spacetime, either at the bra or the ket part of the spacetime in Schwinger-Keldysh formalism. The entanglement entropy is given by an average of complex pseudo generalized entropy for each entanglement island. We find a necessary condition of the initial state to be consistent with the strong sub-additivity, which requires that any probe degrees of freedom are thermally entangled with the rest of the system. We then find a large parameter region where the spacetime with finite-sized initial boundary, which does not have the factorization puzzle at leading order, dominates over the Hartle-Hawking no-boundary state or the bra ket wormhole. Due to the absence of a moment of time reflection symmetry, the island in our setup is a generalization of the entanglement wedge, called pseudo entanglement wedge. In pseudo entanglement wedge reconstruction, we consider reconstructing the bulk matter transition matrix on A ∪ I, from a fine-grained state on A. The bulk transition matrix is given by a thermofield double state with a projection by the initial state. We also provide an AdS/BCFT model by considering EOW branes with corners. We also find the exponential hardness of such reconstruction task using a generalization of Python’s lunch conjecture to pseudo generalized entropy.

2D gravity↗

NASA NTRS1991

Mutual impedance of parallel and perpendicular coplanar surface monopoles

One dimensional integral formulas are derived for mutual impedance of arbitrary size, coplanar, parallel, and perpendicular surface monopoles. The integrals in formulas are expressed as exponential integrals where possible. The mutual impedance expression for parallel monopoles is a summation of exponential integrals and one-dimensional integrals. For perpendicular monopoles, the mutual impedance is in closed form, containing exponential integrals only. The final expressions are in a form suitable for numerical computation. Since the expressions contain at most one-dimensional integrals, they can be utilized to reduce the matrix filling time in the moment method formulations, especially when inhomogeneous sectioning is preferred. Additionally, they can be used in rectangular surface patch modeling of conducting surfaces with edges which are at an angle to the surface patches, providing the angle is small. To this end, the expressions were utilized in the moment method analysis of linearly tapered slot antennas. Very good accuracy was obtained with a reduction in computer time.

Koksal, Adnan↗

DOE OSTI2020

Faster Johnson–Lindenstrauss transforms via Kronecker products

The Kronecker product is an important matrix operation with a wide range of applications in signal processing, graph theory, quantum computing and deep learning. In this work, we introduce a generalization of the fast Johnson–Lindenstrauss projection for embedding vectors with Kronecker product structure, the Kronecker fast Johnson–Lindenstrauss transform (KFJLT). The KFJLT reduces the embedding cost by an exponential factor of the standard fast Johnson–Lindenstrauss transform’s cost when applied to vectors with Kronecker structure, by avoiding explicitly forming the full Kronecker products. Here, we prove that this computational gain comes with only a small price in embedding power: consider a finite set of $p$ points in a tensor product of $d$ constituent Euclidean spaces $\bigotimes _{k=d}^{1}{\mathbb{R}}^{n_k}$, and let $N = \prod _{k=1}^{d}n_k$. With high probability, a random KFJLT matrix of dimension $m \times N$ embeds the set of points up to multiplicative distortion $(1\pm \varepsilon )$ provided $m \gtrsim \varepsilon ^{-2} \, \log ^{2d - 1} (p) \, \log N$. We conclude by describing a direct application of the KFJLT to the efficient solution of large-scale Kronecker-structured least squares problems for fitting the CP tensor decomposition.

Kronecker structure↗

DOE OSTI2022

Radiation induced amorphization of carbides in additively manufactured and conventional ferritic-martensitic steels: In-situ experiments on extraction replicas

In-situ irradiations using 1 MeV Kr 2+ ions in a transmission electron microscope were performed on extraction replica samples containing M 23 C 6 carbides and MX carbonitrides from two ferritic-martensitic (FM) steels: (i) 9%Cr-1%W-TaV based Eurofer97 and (ii) 9%Cr-1%Mo-VNb (all in wt.%) based additively manufactured (AM) Grade91 steel. The irradiations were performed between 100 and 773 K up to a maximum dose of ~2.4 displacements per atom (dpa). The M 23 C 6 carbides are highly susceptible to radiation induced amorphization (RIA), while the MX type nanoprecipitates are highly amorphization resistant across the entire irradiation temperature range. Between 100 and 423 K, RIA of M 23 C 6 carbides occurs very rapidly with critical amorphization doses ranging between ~0.35 and 0.9 dpa, increasing to higher doses at 573 K. Complete amorphization of the M 23 C 6 carbides up to doses of ~2.4 dpa is not possible at 773 K under the present irradiation conditions. The critical amorphization dose of M 23 C 6 carbides increases nearly exponentially with irradiation temperature. The critical temperature for the crystalline-to-amorphous phase transformation (T c→a ) of M 23 C 6 carbides irradiated as replica samples, i.e. without the surrounding metal matrix, was estimated to be ~812 K. Comparing the present results with neutron irradiation data on bulk samples reveals a decrease in the critical amorphization dose for M 23 C 6 and an increase of T c→a , highlighting the effect of dose rate on the amorphization behavior that is qualitatively consistent with literature on other non-metals. Changes in the minor chemistry of M 23 C 6 , such as presence or absence of W, V, Nb, Mo, seems to have little effect on the amorphization behavior. Here, in-situ irradiations on extraction replica samples provide a novel pathway to explore radiation tolerance of nanoprecipitates in nuclear structural materials

36 MATERIALS SCIENCE↗

DOE OSTI2024

Explicit Quantum Circuits for Block Encodings of Certain Sparse Matrices

Many standard linear algebra problems can be solved on a quantum computer by using recently developed quantum linear algebra algorithms that make use of block encodings and quantum eigenvalue/singular value transformations. A block encoding embeds a properly scaled matrix of interest A in a larger unitary transformation U that can be decomposed into a product of simpler unitaries and implemented efficiently on a quantum computer. Although quantum algorithms can potentially achieve exponential speedup in solving linear algebra problems compared to the best classical algorithm, such a gain in efficiency ultimately hinges on our ability to construct an efficient quantum circuit for the block encoding of A, which is difficult in general, and not trivial even for well structured sparse matrices. Here, in this paper, we give a few examples on how efficient quantum circuits can be explicitly constructed for some well structured sparse matrices and discuss a few strategies used in these constructions. We also provide implementations of these quantum circuits in MATLAB.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Soft factorisation and exponentiation from Schwinger-space geometry

Infrared divergences in Quantum Field Theory govern the low-energy dynamics of many physical theories, and their understanding is a crucial ingredient in predicting the outcomes of collider experiments. We present a novel approach to deriving the structure of these divergences by employing the Schwinger parametrization of Feynman integrals. After using tropical geometry to identify divergent limits, we study the all-orders asymptotic properties of Feynman diagrams via matrix manipulations of graph Laplacians, which allows us to analyse their IR behaviour systematically. We explicitly demonstrate the soft-hard factorization of the integrand for a broad class of diagrams, and reveal that when written in terms of worldline distances, topologically distinct diagrams asymptote to the same integrand at leading order in the soft limit. In particular, for the case of Quantum Electrodynamics (with massive fermions), we use this fact to show how ladder-type diagrams combine in Schwinger-parameter space to yield the correct exponentiated soft anomalous dimension. This framework provides a foundation for extending these methods to more complex theories like Quantum Chromodynamics and offers a pathway towards a systematic understanding of infrared divergences in perturbative amplitudes.

Factorization↗

DOE OSTI2020

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Many-Body Level Statistics of Single-Particle Quantum Chaos

We consider a noninteracting many-fermion system populating levels of a unitary random matrix ensemble (equivalent to the q = 2 complex Sachdev-Ye-Kitaev model)—a generic model of single-particle quantum chaos. We study the corresponding many-particle level statistics by calculating the spectral form factor analytically using algebraic methods of random matrix theory, and match it with an exact numerical simulation. Despite the integrability of the theory, the many-body spectral rigidity is found to have a surprisingly rich landscape. In particular, we find a residual repulsion of distant many-body levels stemming from single-particle chaos, together with islands of level attraction. These results are encoded in an exponential ramp in the spectral form factor, which we show to be a universal feature of nonergodic many-fermion systems embedded in a chaotic medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗