Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix theories”

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 91 records · Page 5

Square Root Statistics of Density Matrices and Their Applications

To estimate the degree of quantum entanglement of random pure states, it is crucial to understand the statistical behavior of entanglement indicators such as the von Neumann entropy, quantum purity, and entanglement capacity. These entanglement metrics are functions of the spectrum of density matrices, and their statistical behavior over different generic state ensembles have been intensively studied in the literature. As an alternative metric, in this work, we study the sum of the square root spectrum of density matrices, which is relevant to negativity and fidelity in quantum information processing. In particular, we derive the finite-size mean and variance formulas of the sum of the square root spectrum over the Bures–Hall ensemble, extending known results obtained recently over the Hilbert–Schmidt ensemble.

97 MATHEMATICS AND COMPUTING↗

Statistical uncertainty of fission matrix eigenvalues using perturbation theory

Eigenvalue search of high dominance ratio systems may be slow to converge. The fission matrix element is defined by its element (FM){sub ij}, which are the probability for a neutron born in cell i to create a fission in cell j for a spatial mesh of n{sub i}*n{sub j} elements. Fission matrices are used in Monte Carlo criticality simulations to enhance computing speed, but also to find higher order eigenvalues. However, few studies have been made on the link between statistical uncertainties of fission matrix elements and eigenvalues uncertainties. Thus, dominance ratio statistical uncertainties remain unknown. This paper uses a new generalized perturbation theory (GPT) method to estimate sensitivities of eigenvalues to fission matrix elements and then to calculate dominance ratio uncertainties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Holographic Quantum Simulation of Strongly Correlated Electron Systems

The project aimed to demonstrate a new holographic quantum simulation approach and co‐ designed quantum hardware to tackle three specific problems that fall within the broad umbrella of unraveling the physics of strongly correlated electron systems (SCES). These tasks were: (1) holographic preparation of ground‐ and thermal‐ states of correlated magnetic and electronic systems including quasi‐2d frustrated‐spin, Fermi‐Hubbard, and fractional quantum Hall (FQH) systems, (2) holographic‐simulation of long‐time out‐of‐equilibrium dynamics and (3) holographic analogs of embedding methods such as dynamical mean‐ field theory (DMFT) and density‐matrix embedding theory (DMET) to solve systems with complex structure or long‐range interactions. These tasks are prototypes for the kinds of material simulation problems of interest to BES, such as the simulation of multiferroic materials, perovskite photovoltaics and high‐temperature superconductors, that tax the capabilities of the most powerful classical supercomputers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the S-matrix of Ising field theory in two dimensions

We explore the analytic structure of the non-perturbative S-matrix in arguably the simplest family of massive non-integrable quantum field theories: the Ising field theory (IFT) in two dimensions, which may be viewed as the Ising CFT deformed by its two relevant operators, or equivalently, the scaling limit of the Ising model in a magnetic field. Our strategy is that of collider physics: we employ Hamiltonian truncation method (TFFSA) to extract the scattering phase of the lightest particles in the elastic regime, and combine it with S-matrix bootstrap methods based on unitarity and analyticity assumptions to determine the analytic continuation of the 2 → 2 S-matrix element to the complex s-plane. Focusing primarily on the “high temperature” regime in which the IFT interpolates between that of a weakly coupled massive fermion and the E 8 affine Toda theory, we will numerically determine 3-particle amplitudes, follow the evolution of poles and certain resonances of the S-matrix, and exclude the possibility of unknown wide resonances up to reasonably high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Infinite matrix product states for (1 + 1)-dimensional gauge theories

We present a matrix product operator construction that allows us to represent the lattice Hamiltonians of (abelian or non-abelian) gauge theories in a local and manifestly translation-invariant form. In particular, we use symmetric matrix product states and introduce link-enhanced matrix product operators (LEMPOs) that can act on both the physical and virtual spaces of the matrix product states. This construction allows us to study Hamiltonian lattice gauge theories on infinite lattices. As examples, we show how to implement this method to study the massless and massive one-flavor Schwinger model and adjoint QCD 2 .

confinement↗

Shear and bulk viscosity for a pure glue theory using an effective matrix model

At nonzero temperatures, the deconfining phase transition can be analyzed using an effective matrix model to characterize the change in holonomy. The model includes gluons and two-dimensional ghost fields in the adjoint representation, or “teens.” As ghosts, the teen fields are responsible for the decrease of the pressure as 𝑇 →𝑇 𝑑 , with 𝑇 𝑑 the transition temperature for deconfinement. Using the solution of this matrix model for a large number of colors, the parameters of the teen fields are adjusted so that the expectation value of the Polyakov loop is close to the values from the lattice. The shear, 𝜂, and bulk, 𝜁, viscosities are computed at nonzero holonomy to leading logarithmic order in weak coupling. In the pure glue theory, the value of the Polyakov loop is relatively large in the deconfined phase, ≈1/2 at 𝑇 𝑑 . Consequently, if 𝑠 is the entropy density, while 𝜂/𝑠 decreases as 𝑇 →𝑇 𝑑 , it is still well above the conformal bound. In contrast, 𝜁/𝑠 is largest at 𝑇 𝑑 , comparable to 𝜂/𝑠, then falls off rapidly with increasing temperature and is negligible by ∼2⁢𝑇 𝑑 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k · p two-band model

We investigate the twisted bilayer graphene (TBG) model of Bistritzer and MacDonald (BM) [Bistritzer and MacDonald, Proc. Natl. Acad. Sci. 108, 12233 (2011)] to obtain an analytic understanding of its energetics and wave functions needed for many-body calculations. We provide an approximation scheme for the wave functions of the BM model, which first elucidates why the BM K M -point centered original calculation containing only four plane waves provides a good analytical value for the first magic angle (θ M ≈ 1°). The approximation scheme also elucidates why most of the many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a Γ M -point centered model of six plane waves, we analytically understand the reason for the small Γ M -point gap between the active and passive bands in the isotropic limit w 0 = w 1 . Furthermore, we analytically calculate the group velocities of the passive bands in the isotropic limit, and show that they are almost doubly degenerate, even away from the Γ M point, where no symmetry forces them to be. Furthermore, moving away from the Γ M and K M points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite the first magic angle is defined by only requiring a vanishing K M -point Dirac velocity. We derive analytically a connected “magic manifold” w 1 = $2\sqrt{1 + w^{2}_{0}}$ $-\sqrt{2 + 3w^2_0}$, on which the bands remain extremely flat as w 0 is tuned between the isotropic (w 0 = w 1 ) and chiral (w 0 = 0) limits. We analytically show why going away from the isotropic limit by making w 0 less (but not larger) than w 1 increases the Γ M -point gap between the active and the passive bands. Finally, by perturbation theory, we provide an analytic Γ M point k ∙ p two-band model that reproduces the TBG band structure and eigenstates within a certain w 0 , w 1 parameter range. Further refinement of this model are discussed, which suggest a possible faithful representation of the TBG bands by a two-band Γ M point k ∙ p model in the full w 0 , w 1 parameter range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The S-matrix of 2D type 0B string theory. Part I. Perturbation theory revisited

We study the perturbative S-matrix of closed strings in the two-dimensional type 0B string theory from the worldsheet perspective, by directly integrating correlation functions of N = 1 Liouville theory. The latter is computed numerically using recurrence relations for super-Virasoro conformal blocks. We show that the tree level 3- and 4-point amplitudes are in agreement with the proposed dual matrix quantum mechanics. The non-perturbative aspects of the duality will be analyzed in a companion paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Confinement/deconfinement transition in the D0-brane matrix model — A signature of M-theory?

We study the confinement/deconfinement transition in the D0-brane matrix model (often called the BFSS matrix model) and its one-parameter deformation (the BMN matrix model) numerically by lattice Monte Carlo simulations. Our results confirm general expectations from the dual string/M-theory picture for strong coupling. In particular, we observe the confined phase in the BFSS matrix model, which is a nontrivial consequence of the M-theory picture. We suggest that these models provide us with an ideal framework to study the Schwarzschild black hole, M-theory, and furthermore, the parameter region of the phase transition between type IIA superstring theory and M-theory. A detailed study of M-theory via lattice Monte Carlo simulations of the D0-brane matrix model might be doable with much smaller computational resources than previously expected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum Perturbation Theory Using Tensor Cores and a Deep Neural Network

In this work, time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e., dense matrix–matrix multiplications, in mixed-precision arithmetics, which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree–Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low-precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗

A linear surrogate for optimising functions of an orthogonal matrix with applications in wave function theory

The technique of surrogate optimisation is to use a simpler function to approximate a complex function that is time-consuming to evaluate. Here we show that the maximum of a special type of surrogate function f(U)=Tr(AU),UϵO(n) is at A T (AA T ) 1/2 , and that there is one and only one local maximum both in SO(n) and O(n)–SO(n). This function f(U) has been found to be useful in various aspects of electronic structure theory, including proving the Carlson-Keller theorem, and localising orbitals. As one other example, we apply it here to optimise the ground state of molecules using the Generalised Valence Bond wavefunction.

74 ATOMIC AND MOLECULAR PHYSICS↗

Structured ground states of randomly interacting bosons

Bosonic degrees of freedom and their emergence as part of complex quantum many-body dynamics, symmetries, collective behavior, clustering, and phase transitions play an important role in modern studies of quantum systems. In this paper, we present a systematic study of many-boson systems governed by random interactions. Our findings show that ground states of randomly interacting bosons are not random, being dominated by a few collective configurations containing condensates of clusters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Matrix ensembles with global symmetries and ’t Hooft anomalies from 2d gauge theory

The Hilbert space of a quantum system with internal global symmetry $\textit{G}$ decomposes into sectors labelled by irreducible representations of $\textit{G}$. If the system is chaotic, the energies in each sector should separately resemble ordinary random matrix theory. We show that such “sector-wise” random matrix ensembles arise as the boundary dual of two- dimensional gravity with a $\textit{G}$ gauge field in the bulk. Within each sector, the eigenvalue density is enhanced by a nontrivial factor of the dimension of the representation, and the ground state energy is determined by the quadratic Casimir. We study the consequences of ’t Hooft anomalies in the matrix ensembles, which are incorporated by adding specific topological terms to the gauge theory action. The effect is to introduce projective representations into the decomposition of the Hilbert space. Finally, we consider ensembles with $\textit{G}$ symmetry and time reversal symmetry, and analyze a simple case of a mixed anomaly between time reversal and an internal $\mathbb{Z}$ 2 symmetry.

2D Gravity↗