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 181 records · Page 10

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↗

Pseudospectral calculation of helium wave functions, expectation values, and oscillator strength

We show that the pseudospectral method is a powerful tool for finding precise solutions of Schroedinger's equation for two-electron atoms with general angular momentum. Realizing the method's full promise for atomic calculations requires special handling of singularities due to two-particle Coulomb interactions. We give a prescription for choosing coordinates and subdomains whose efficacy we illustrate by solving several challenging problems. One test centers on the determination of the nonrelativistic electric dipole oscillator strength for the helium 1 {sup 1} S{yields}2 {sup 1} P transition. The result achieved, 0.276 164 99(27), is comparable to the best in the literature. The formally equivalent length, velocity, and acceleration expressions for the oscillator strength all yield roughly the same accuracy. We also calculate a diverse set of helium ground-state expectation values, reaching near state-of-the-art accuracy without the necessity of implementing any special-purpose numerics. These successes imply that general matrix elements are directly and reliably calculable with pseudospectral methods. A striking result is that all the relevant quantities tested in this paper--energy eigenvalues, S-state expectation values and a bound-bound dipole transition between the lowest energy S and P states--converge exponentially with increasing resolution and at roughly the same rate. Each individual calculation samples and weights the configuration space wave function uniquely but all behave in a qualitatively similar manner. These results suggest that the method has great promise for similarly accurate treatment of few-particle systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

More on rotations as spin matrix polynomials

Any nonsingular function of spin j matrices always reduces to a matrix polynomial of order 2j. The challenge is to find a convenient form for the coefficients of the matrix polynomial. The theory of biorthogonal systems is a useful framework to meet this challenge. Central factorial numbers play a key role in the theoretical development. Explicit polynomial coefficients for rotations expressed either as exponentials or as rational Cayley transforms are considered here. Structural features of the results are discussed and compared, and large j limits of the coefficients are examined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Remarks on non-BPS string amplitudes and their all order α' contact interactions in IIB, IIA

Here, we explore the entire form of S-Matrix elements of a potential C n–1 Ramond-Ramond (RR) form field, a tachyon and two transverse scalar fields on both world volume and transverse directions of type IIB and IIA superstring theories. Apart from V C –2 V Φ 0 V Φ 0 V T 0 the other scattering amplitude, namely V C –1 V Φ –1 V Φ 0 V T 0 is also revealed. We then start to compare all singularity structures of symmetric and asymmetric analysis, generating all infinite singularity structures as well as all order α' contact interactions on the whole directions. This leads to deriving various new contact terms and several new restricted Bianchi identities in both type IIB and IIA. It is also shown that just some of the new couplings of type IIB (IIA) string theory can be re-verified in an Effective Field Theory (EFT) by pull-back of branes. To construct the rest of S-matrix elements one needs to first derive restricted world volume (or bulk) Bianchi identities and then discover new EFT couplings in both type IIB and IIA. Finally the presence of commutator of scalar fields inside the exponential of Wess-Zumino action for non-BPS branes has been confirmed as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Absence of Barren Plateaus and Scaling of Gradients in the Energy Optimization of Isometric Tensor Network States

Abstract Vanishing gradients can pose substantial obstacles for high-dimensional optimization problems. Here we consider energy minimization problems for quantum many-body systems with extensive Hamiltonians and finite-range interactions, which can be studied on classical computers or in the form of variational quantum eigensolvers on quantum computers. Barren plateaus correspond to scenarios where the average amplitude of the energy gradient decreases exponentially with increasing system size. This occurs, for example, for quantum neural networks and for brickwall quantum circuits when the depth increases polynomially in the system size. Here we prove that the variational optimization problems for matrix product states, tree tensor networks, and the multiscale entanglement renormalization ansatz are free of barren plateaus. The derived scaling properties for the gradient variance provide an analytical guarantee for the trainability of randomly initialized tensor network states (TNS) and motivate certain initialization schemes. In a suitable representation, unitary tensors that parametrize the TNS are sampled according to the uniform Haar measure. We employ a Riemannian formulation of the gradient based optimizations which simplifies the analytical evaluation.

Barthel, Thomas↗

Influence of interfacial oxide on the optical properties of single layer CdTe/CdS quantum dots in porous silicon scaffolds

Using a combination of continuous wave and time-resolved spectroscopy, we study the effects of interfacial conditions on the radiative lifetimes and photoluminescence intensities of sub-monolayer colloidal CdTe/CdS quantum dots (QDs) embedded in a three-dimensional porous silicon (PSi) scaffold. The PSi matrix was thermally oxidized under different conditions to change the interfacial oxide thickness. QDs embedded in a PSi matrix with ∼0.4 nm of interfacial oxide exhibited reduced photoluminescence intensity and nearly five times shorter radiative lifetimes (∼16 ns) compared to QDs immobilized within completely oxidized, porous silica (PSiO{sub 2}) frameworks (∼78 ns). The exponential dependence of QD lifetime on interfacial oxide thickness in the PSi scaffolds suggests charge transfer plays an important role in the exciton dynamics.

36 MATERIALS SCIENCE↗

Quantum algorithms for Gibbs sampling and hitting-time estimation

In this paper, we present quantum algorithms for solving two problems regarding stochastic processes. The first algorithm prepares the thermal Gibbs state of a quantum system and runs in time almost linear in √Nβ/Ζ and polynomial in log(1/ϵ), where N is the Hilbert space dimension, β is the inverse temperature, Ζ is the partition function, and ϵ is the desired precision of the output state. Our quantum algorithm exponentially improves the dependence on 1/ϵ and quadratically improves the dependence on β of known quantum algorithms for this problem. The second algorithm estimates the hitting time of a Markov chain. For a sparse stochastic matrix Ρ, it runs in time almost linear in 1/(ϵΔ 3/2 ), where ϵ is the absolute precision in the estimation and Δ is a parameter determined by Ρ, and whose inverse is an upper bound of the hitting time. Our quantum algorithm quadratically improves the dependence on 1/ϵ and 1/Δ of the analog classical algorithm for hitting-time estimation. Finally, both algorithms use tools recently developed in the context of Hamiltonian simulation, spectral gap amplification, and solving linear systems of equations.

97 MATHEMATICS AND COMPUTING↗

Quantum Gauge Networks: A New Kind of Tensor Network

Although tensor networks are powerful tools for simulating low-dimensional quantum physics, tensor network algorithms are very computationally costly in higher spatial dimensions. We introduce quantum gauge networks: a different kind of tensor network ansatz for which the computation cost of simulations does not explicitly increase for larger spatial dimensions. We take inspiration from the gauge picture of quantum dynamics, which consists of a local wavefunction for each patch of space, with neighboring patches related by unitary connections. A quantum gauge network (QGN) has a similar structure, except the Hilbert space dimensions of the local wavefunctions and connections are truncated. We describe how a QGN can be obtained from a generic wavefunction or matrix product state (MPS). All 2k-point correlation functions of any wavefunction for M many operators can be encoded exactly by a QGN with bond dimension O(M k ). In comparison, for just k = 1, an exponentially larger bond dimension of 2 M/6 is generically required for an MPS of qubits. We provide a simple QGN algorithm for approximate simulations of quantum dynamics in any spatial dimension. The approximate dynamics can achieve exact energy conservation for time-independent Hamiltonians, and spatial symmetries can also be maintained exactly. We benchmark the algorithm by simulating the quantum quench of fermionic Hamiltonians in up to three spatial dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A nonperturbative light-front coupled-cluster method

The nonperturbative Hamiltonian eigenvalue problem for bound states of a quantum field theory is formulated in terms of Dirac's light-front coordinates and then approximated by the exponential-operator technique of the many-body coupled-cluster method. This approximation eliminates any need for the usual approximation of Fock-space truncation. Instead, the exponentiated operator is truncated, and the terms retained are determined by a set of nonlinear integral equations. These equations are solved simultaneously with an effective eigenvalue problem in the valence sector, where the number of constituents is small. Matrix elements can be calculated, with extensions of techniques from standard coupled-cluster theory, to obtain form factors and other observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pure Yang–Mills Solutions on dS{sub 4}

We consider pure SU(2) Yang–Mills theory on four-dimensional de Sitter space dS{sub 4} and construct smooth and spatially homogeneous classical Yang–Mills fields. Slicing dS{sub 4} as R×S{sup 3}, via an SU(2)-equivariant ansatz we reduce the Yang–Mills equations to ordinary matrix differential equations and further to Newtonian dynamics in a particular three-dimensional potential. Its classical trajectories yield spatially homogeneous Yang–Mills solutions in a very simple explicit form, depending only on de Sitter time with an exponential decay in the past and future. These configurations have not only finite energy, but their action is also finite and bounded from below. We present explicit coordinate representations of the simplest examples (for the fundamental SU(2) representation). Instantons (Yang–Mills solutions on the Wick-rotated S{sup 4} ) and solutions on AdS{sub 4} are also briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Subleading corrections in $ \mathcal{N} $ = 3 Gaiotto-Tomasiello theory

We study subleading corrections to the genus-zero free energy of the N = 3 Gaiotto-Tomasiello theory. In general, we obtain the endpoints and free energy as a set of parametric equations via contour integrals of the planar resolvent, up to exponentially suppressed corrections. In the particular case that the two gauge groups in the quiver are of equal rank, we find an explicit (perturbative) expansion for the free energy. If, additionally, both groups have equal levels, then we find the full expression for the genus-zero free energy, modulo exponentially suppressed corrections. We also verify our results numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Refinements for Bragg coherent X-ray diffraction imaging: electron backscatter diffraction alignment and strain field computation

Bragg coherent X-ray diffraction imaging (BCDI) allows the 3D measurement of lattice strain along the scattering vector for specific microcrystals. If at least three linearly independent reflections are measured, the 3D variation of the full lattice strain tensor within the microcrystal can be recovered. However, this requires knowledge of the crystal orientation, which is typically attained via estimates based on crystal geometry or synchrotron microbeam Laue diffraction measurements. Presented here is an alternative method to determine the crystal orientation for BCDI measurements using electron backscatter diffraction (EBSD) to align Fe–Ni and Co–Fe alloy microcrystals on three different substrates. The orientation matrix is calculated from EBSD Euler angles and compared with the orientation determined using microbeam Laue diffraction. The average angular mismatch between the orientation matrices is less than ∼6°, which is reasonable for the search for Bragg reflections. The use of an orientation matrix derived from EBSD is demonstrated to align and measure five reflections for a single Fe–Ni microcrystal via multi-reflection BCDI. Using this data set, a refined strain field computation based on the gradient of the complex exponential of the phase is developed. This approach is shown to increase accuracy, especially in the presence of dislocations. The results demonstrate the feasibility of using EBSD to pre-align BCDI samples and the application of more efficient approaches to determine the full lattice strain tensor with greater accuracy.

42 ENGINEERING↗

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. Finally, we observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Low-Rank QTT-based Finite Element Method for Elasticity Problems

We present an efficient and robust numerical algorithm for solving the linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. By integrating Z-ordering and subdomain concatenation, our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This efficiency is maintained while still guaranteeing exponential convergence with respect to the number of degrees of freedom. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented. This includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method’s performance in the numerical approximation of three representative test cases.

97 MATHEMATICS AND COMPUTING↗