Search NASA⌕ Search

SEARCH · Search NASA

Results for “Exponential Matrix”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A deterministic alternative to the full configuration interaction quantum Monte Carlo method

Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such approach, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, along with a stochastic projected wave function, to find the important parts of Hilbert space. Yet, the stochastic representation of the wave function is not required to search Hilbert space efficiently, and here we describe a highly efficient deterministic method that can achieve chemical accuracy for a wide range of systems, including the difficult Cr 2 molecule. We demonstrate for systems like Cr 2 that such calculations can be performed in just a few cpu hours which makes it one of the most efficient and accurate methods that can attain chemical accuracy for strongly correlated systems. Furthermore, our method also allows efficient calculation of excited state energies, which we illustrate with benchmark results for the excited states of C 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2012

Effect of disorder with long-range correlation on transport in graphene nanoribbon

Transport in disordered armchair graphene nanoribbons (AGR) with long-range correlation between quantum wire contacts is investigated by a transfer matrix combined with Landauer's formula. The metal-insulator transition is induced by disorder in neutral AGR. Therein, the conductance is one conductance quantum for the metallic phase and exponentially decays otherwise, when the length of AGR approaches infinity and far longer than its width. Similar to the case of long-range disorder, the conductance of neutral AGR first increases and then decreases while the conductance of doped AGR monotonically decreases, as the disorder strength increases. In the presence of strong disorder, the conductivity depends monotonically and non-monotonically on the aspect ratio for heavily doped and slightly doped AGR, respectively. For edge disordered graphene nanoribbon, the conductance increases with the disorder strength of long-range correlated disordered while no delocalization exists, since the edge disorder induces localization.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Dynamics of Nano-Confined Water under Pressure

We report a study of the effects of pressure on the diffusivity of water molecules confined in single- wall carbon nanotubes (SWNT) with average mean pore diameter of 16 A. The measurements were carried out using high-resolution neutron scattering, over the temperature range 220 T 260 K, and at two pressure conditions: ambient and elevated pressure. The high pressure data were collected at constant volume on cooling, with P varying from 1.92 kbar at temperature T = 260 K to 1.85 kbar at T = 220 K. Analysis of the observed dynamic structure factor S(Q, E) reveals the presence of two relaxation processes, a faster diffusion component (FC) associated with the motion of caged or restricted molecules, and a slower component arising from the free water molecules diffusing within the SWNT matrix. While the temperature dependence of the slow relaxation time exhibits a Vogel-Fulcher-Tammann law and is non-Arrhenius in nature, the faster component follows an Arrhenius exponential law at both pressure conditions. The application of pressure remarkably slows down the overall molecular dynamics, in agreement with previous observations, but most notably affects the slow relaxation. The faster relaxation shows marginal or no change with pressure within the experimental conditions.

Confined Water↗
DOE OSTI2019

Activation behavior for ion permeation in ion-exchange membranes: Role of ion dehydration in selective transport

Here we explored the mechanisms governing the selectivity of anion- and cation-exchange membranes for the transport of four monovalent anions (i.e., fluoride, chloride, bromide, and nitrate) and four monovalent cations (i.e., sodium, potassium, cesium, and ammonium), respectively. Our ion adsorption and transport tests with mixed ion solutions reveal that an ion with larger ionic radius and lower hydration energy is more favorably adsorbed onto the ion-exchange membrane but diffuses more slowly through the polymer matrix compared to an ion with smaller ionic radius and higher hydration energy. Individual anion (as sodium salt) or cation (as chloride salt) permeation tests at different temperatures were performed to evaluate the activation behavior of ion transport through the ion-exchange membranes by calculating the energy barrier and pre-exponential factor (i.e., the ion flux when the energy barrier is negligible) for ion transport from an Arrhenius-type equation. Our results show that an ion with smaller ionic radius and higher hydration energy experiences higher energy barrier (e.g., fluoride, 10.3 kcal mol -1 ) and possesses higher pre-exponential factor compared to an ion with larger ionic radius and lower hydration energy (e.g., bromide, 4.6 kcal mol -1 ). This correlation corroborates our main hypothesis that the activation behavior observed for ion transport is a result of ion dehydration at the water-membrane interface. Our proposed ion selectivity mechanism elucidates how ion dehydration governs the extent of ion permeation into the membrane and the subsequent transport through the charged polymer matrix. Future membrane design that promotes dehydration of target ions is challenging but can result in unprecedented ion selectivity.

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

On the thermal stability of physical vapor deposited oxide-hardened nanocrystalline gold thin films

We describe a correlation between electrical resistivity and grain size for PVD synthesized polycrystalline oxide-hardened metal-matrix thin films in oxide-dilute (<5 vol. % oxide phase) compositions. The correlation is based on the Mayadas-Shatzkes (M-S) electron scattering model, predictive of grain size evolution as a function of composition in the oxide-dilute regime for 2 μm thick Au-ZnO films. We describe a technique to investigate grain boundary (GB) mobility and the thermal stability of GBs based on in situ electrical resistivity measurements during annealing experiments, interpreted using a combination of the M-S model and the Michels et al. model describing solute drag stabilized grain growth kinetics. Using this technique, activation energy and pre-exponential Arrhenius parameter values of E{sub a} = 21.6 kJ/mol and A{sub o} = 2.3 × 10{sup −17} m{sup 2}/s for Au-1 vol. % ZnO and E{sub a} = 12.7 kJ/mol and A{sub o} = 3.1 × 10{sup −18} m{sup 2}/s for Au-2 vol. % ZnO were determined. In the oxide-dilute regime, the grain size reduction of the Au matrix yielded a maximum hardness of 2.6 GPa for 5 vol. % ZnO. A combined model including percolation behavior and grain refinement is presented that accurately describes the composition dependent change in electrical resistivity throughout the entire composition range for Au-ZnO thin films. The proposed correlations are supported by microstructural characterization using transmission electron microscopy and electron diffraction mapping for grain size determination.

36 MATERIALS SCIENCE↗
DOE OSTI2010

Simulation of classical thermal states on a quantum computer: A transfer-matrix approach

We present a hybrid quantum-classical algorithm to simulate thermal states of classical Hamiltonians on a quantum computer. Our scheme employs a sequence of locally controlled rotations, building up the desired state by adding qubits one at a time. We identified a class of classical models for which our method is efficient and avoids potential exponential overheads encountered by Grover-like or quantum Metropolis schemes. Our algorithm also gives an exponential advantage for two-dimensional Ising models with magnetic field on a square lattice, compared with the previously known Zalka's algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Communication: Active space decomposition with multiple sites: Density matrix renormalization group algorithm

We extend the active space decomposition method, recently developed by us, to more than two active sites using the density matrix renormalization group algorithm. The fragment wave functions are described by complete or restricted active-space wave functions. Numerical results are shown on a benzene pentamer and a perylene diimide trimer. It is found that the truncation errors in our method decrease almost exponentially with respect to the number of renormalization states M, allowing for numerically exact calculations (to a few μE{sub h} or less) with M = 128 in both cases. This rapid convergence is because the renormalization steps are used only for the interfragment electron correlation.

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

Linearity of holographic entanglement entropy

Here, we consider the question of whether the leading contribution to the entanglement entropy in holographic CFTs is truly given by the expectation value of a linear operator as is suggested by the Ryu-Takayanagi formula. We investigate this property by computing the entanglement entropy, via the replica trick, in states dual to superpositions of macroscopically distinct geometries and find it consistent with evaluating the expectation value of the area operator within such states. However, we find that this fails once the number of semi-classical states in the superposition grows exponentially in the central charge of the CFT. Moreover, in certain such scenarios we find that the choice of surface on which to evaluate the area operator depends on the density matrix of the entire CFT. This nonlinearity is enforced in the bulk via the homology prescription of Ryu-Takayanagi. We thus conclude that the homology constraint is not a linear property in the CFT. We also discuss the existence of entropy operators in general systems with a large number of degrees of freedom.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
DOE OSTI2014

Lignin-Derived Carbon Fiber as a Co-Product of Refining Cellulosic Biomass

Lignin by-products from biorefineries has the potential to provide a low-cost alternative to petroleum-based precursors to manufacture carbon fiber, which can be combined with a binding matrix to produce a structural material with much greater specific strength and specific stiffness than conventional materials such as steel and aluminum. The market for carbon fiber is universally projected to grow exponentially to fill the needs of clean energy technologies such as wind turbines and to improve the fuel economies in vehicles through lightweighting. In addition to cellulosic biofuel production, lignin-based carbon fiber production coupled with biorefineries may provide $2,400 to $3,600 added value dry Mg-1 of biomass for vehicle applications. Compared to producing ethanol alone, the addition of lignin-derived carbon fiber could increase biorefinery gross revenue by 30% to 300%. Using lignin-derived carbon fiber in 15 million vehicles per year in the US could reduce fossil fuel consumption by 2-5 billion liters year-1, reduce CO2 emissions by about 6.7 million Mg year-1, and realize fuel savings through vehicle lightweighting of $700 to $1,600 per Mg biomass processed. The value of fuel savings from vehicle lightweighting becomes economical at carbon fiber price of $6.60 kg-1 under current fuel prices, or $13.20 kg-1 under fuel prices of about $1.16 l-1.

09 BIOMASS FUELS↗
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 OSTI2019

Work Statistics, Loschmidt Echo and Information Scrambling in Chaotic Quantum Systems

Characterizing the work statistics of driven complex quantum systems is generally challenging because of the exponential growth with the system size of the number of transitions involved between different energy levels. We consider the quantum work distribution associated with the driving of chaotic quantum systems described by random matrix Hamiltonians and characterize exactly the work statistics associated with a sudden quench for arbitrary temperature and system size. Knowledge of the work statistics yields the Loschmidt echo dynamics of an entangled state between two copies of the system of interest, the thermofield double state. This echo dynamics is dictated by the spectral form factor. We discuss its relation to frame potentials and its use to assess information scrambling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗