Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Theory”

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 OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Target space entanglement in Matrix Models

We study target space entanglement in gauged multi-matrix models as models of entanglement between groups of D-branes separated by a planar entangling surface, paying close attention to the implementation of gauge invariance. We open with a review of target space entanglement between identical particles, which shares some important features (specifically a gauged permutation symmetry) with our main problem. For our matrix models, we implement a gauge fixing well-adapted to the entangling surface. In this gauge, we map the matrix model problem to that of entanglement of a U(1) gauge theory on a complete or all-to-all lattice. Matrix elements corresponding to open strings stretching across the entangling surface in the target space lead to interesting contributions to the entanglement entropy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Parametric matrix models

We present a general class of machine learning algorithms called parametric matrix models. In contrast with most existing machine learning models that imitate the biology of neurons, parametric matrix models use matrix equations that emulate physical systems. Similar to how physics problems are usually solved, parametric matrix models learn the governing equations that lead to the desired outputs. Parametric matrix models can be efficiently trained from empirical data, and the equations may use algebraic, differential, or integral relations. While originally designed for scientific computing, we prove that parametric matrix models are universal function approximators that can be applied to general machine learning problems. After introducing the underlying theory, we apply parametric matrix models to a series of different challenges that show their performance for a wide range of problems. For all the challenges tested here, parametric matrix models produce accurate results within an efficient and interpretable computational framework that allows for input feature extrapolation.

Computational science↗
DOE OSTI2024

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Two-photon above-threshold ionization of helium

Multiphoton ionization provides a clear window into the nature of electron correlations in the helium atom. In the present study, the final-state energy range extends up to the region near the N = 2 and N = 3 ionization thresholds, where two-photon ionization proceeds via continuum intermediate states above the lowest threshold. Our calculations are performed using multichannel quantum defect theory (MQDT) and the streamlined R-matrix method. The sum and integration over all intermediate states in the two-photon ionization amplitude is evaluated using the inhomogeneous R-matrix method developed by Robicheaux and Gao. The seamless connection of that method with MQDT allows us to present high-resolution spectra of the final-state Rydberg resonances. Our analysis classifies the resonances above the N = 2 threshold in terms of their group theory quantum numbers. Furthermore, their dominant decay channels are found to obey the previously conjectured propensity rule far more weakly for these even-parity states than was observed for the odd-parity states relevant to single-photon ionization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Equivalence of EFTs and timelike extra dimensions

A flat 3-brane probing a five-dimensional anti–de Sitter space has the same symmetries and symmetry breaking pattern as one probing a five-dimensional de Sitter space with two time directions. Despite the seemingly different physical setups, we show that the effective field theory of the brane bending mode is identical in both cases. In particular, despite “wrong” signs in the Dirac-Born-Infeld (DBI) action for the two-time de Sitter case, the theories have the same S -matrix. This is consistent with the expectation that effective field theories are determined solely by their degrees of freedom and pattern of symmetry breaking, even in the case of spacetime symmetries. We comment further on the equivalence between the Weyl/dilaton and DBI representations of the EFT of broken conformal symmetry. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Efficient simulation of low-temperature physics in one-dimensional gapless systems

Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O⁡(cN 0 ⁢log⁡β) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.

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

Matter couplings and equivalence principles for soft scalars

Scalar effective field theories with enhanced soft limits behave in many ways like gauge theories and gravity. In particular, symmetries fix the structure of interactions and the tree-level S-matrix in both types of theories. We explore how this analogy persists in the presence of matter by considering theories with additional fields coupled to the Dirac-Born-Infeld (DBI) scalar or the special galileon in a way that is consistent with their symmetries. Using purely on-shell arguments, we show that these theories obey analogues of the S-matrix equivalence principle whereby all matter fields must couple to the DBI scalar or the special galileon through a particular quartic vertex with a universal coupling. These equivalence principles imply the universality of the leading double soft theorems in these theories, which are scalar analogues of Weinberg’s gravitational soft theorem, and can be used to rule out interactions with massless higher-spin fields when combined with analogues of the generalized Weinberg-Witten theorem. We verify in several examples that amplitudes with external matter fields nontrivially exhibit enhanced single soft limits and we show that such amplitudes can be constructed using soft recursion relations when they have sufficiently many external DBI or special galileon legs, including amplitudes with massive higher-spin fields. As part of our analysis we construct a recently conjectured special galileon-vector effective field theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Renormalization of the Standard Model Effective Field Theory from geometry

S-matrix elements are invariant under field redefinitions of the Lagrangian. They are determined by geometric quantities such as the curvature of the field-space manifold of scalar and gauge fields. We present a formalism where scalar and gauge fields are treated together, with a metric on the combined space of both types of fields. Scalar and gauge scattering amplitudes are given by the Riemann curvature Rijkl of this combined space, with indices i, j, k, l chosen to be scalar or gauge indices depending on the type of external particle. One-loop divergences can also be computed in terms of geometric invariants of the combined space, which greatly simplifies the computation of renormalization group equations. We apply our formalism to the Standard Model Effective Field Theory (SMEFT), and compute the renormalization group equations for even-parity bosonic operators to mass dimension eight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Matrix-Based Process Modeling in Microsoft Excel

The purpose of this document is to describe the theory and methodology of matrix-based process modeling and demonstrate its use through a generic process model within Microsoft Excel. This type of model allows a modeler to quickly query various metrics, such as total production time, for a specified production run. These matrix models tend to be quick to build and use, but come with certain limitations and are more suited for smaller processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Misanthropic entropy and renormalization as a communication channel

A central physical question is the extent to which infrared (IR) observations are sufficient to reconstruct a candidate ultraviolet (UV) completion. We recast this question as a problem of communication, with messages encoded in field configurations of the UV being transmitted to the IR degrees of freedom via a noisy channel specified by renormalization group (RG) flow, with noise generated by coarse graining/decimation. We present an explicit formulation of these considerations in terms of lattice field theory, where we show that the “misanthropic entropy” — the mutual information obtained from decimating neighbors — encodes the extent to which information is lost in marginalizing over/tracing out UV degrees of freedom. Our considerations apply both to statistical field theories as well as density matrix renormalization of quantum systems, where in the quantum case, the statistical field theory analysis amounts to a leading-order approximation. As a concrete example, we focus on the case of the 2D Ising model, where we show that the misanthropic entropy detects the onset of the phase transition at the Ising model critical point.

Physics↗
DOE OSTI2023

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗
DOE OSTI2023

In situ high-temperature 3D imaging of the damage evolution in a SiC nuclear fuel cladding material

Silicon carbide (SiC)-based nuclear fission fuel rod cladding has been considered as one of the possible designs for accident tolerant fuels. It is in the form of a SiC fibre reinforced SiC matrix composite tube (SiC f -SiC m ) with monolithic SiC outer and/or inner coating layers. This study focuses on the deformation and fracture processes in this material using in situ X-ray micro-computed tomography (XCT) at room temperature (RT) and 1200 °C in an inert gas environment in a C-ring compression loading configuration. Prior to testing, local properties and residual stresses were characterised using nanoindentation and Raman spectroscopy since they can impact the mechanical behaviour of the material. The 3D strain distribution, crack formation and propagation processes including the toughening mechanisms (e.g., crack deflection, micro-cracking, crack bridging and bifurcation) are investigated in the coating and underlying composites at RT and 1200 °C. There is no particular sequence which toughening mechanism occurs first – this is very different from the conventional toughening theory in ceramic-matrix composites under uniaxial tension loading. Indeed, no evidence of fibre pull-out or fibre fracture was observed in this SiC f -SiC m nuclear cladding material in the current C-ring compression configuration. The correlation between local measurements and bulk mechanical behaviour are discussed.

3D imaging↗
DOE OSTI2022

Electronic structure, magnetic properties, and pairing tendencies of the copper-based honeycomb lattice Na 2 Cu 2 TeO 6

Spin-1/2 chains with alternating antiferromagnetic (AFM) and ferromagnetic (FM) couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density-matrix renormalization-group (DMRG) methods, we have systematically studied the dimerized chain system Na 2 Cu 2 TeO 6 with a d 9 electronic configuration. Near the Fermi level, in the nonmagnetic phase the dominant states are mainly contributed by the Cu 3d x 2 -y 2 orbitals highly hybridized with the O 2p orbitals, leading to an “effective” single-orbital low-energy model. By calculating the relevant hoping amplitudes, we explain the size and sign of the exchange interactions in Na 2 Cu 2 TeO 6 . In addition, a single-orbital Hubbard model is constructed for this dimerized chain system where the quantum fluctuations are taken into account. Both AFM and FM couplings (leading to an ↑|-|↓-|↓-↑| state) along the chain were found in our DMRG and Lanczos calculations, in agreement with density functional theory and neutron-scattering results. The hole pairing binding energy ΔE is predicted to be negative at Hubbard U~ 11eV, suggesting incipient pairing tendencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Quantum Monte Carlo calculation of 𝛿 NS in 10 C using an effective field theory approach

Here, we compute radiative corrections to the superallowed 𝛽 decay of 10 C in an effective field theory approach using nuclear matrix elements obtained from quantum Monte Carlo calculations. These corrections are an important ingredient in the extraction of the Cabibbo-Kobayashi-Masakawa quark mixing matrix element 𝑉 𝑢⁢𝑑 , and the role of this work is to illuminate the uncertainties arising from nuclear structure. Our results provide good agreement with both the traditional extraction of 𝑉 𝑢⁢𝑑 , as well as with a more recent evaluation performed using the no-core shell model and a dispersion formalism. The dominant uncertainty in this approach is the presence of two unknown low-energy constants that enter into the relevant nuclear matrix elements. Future determinations of these low-energy constants—either from QCD or modeling them with two nucleon amplitudes—would improve the precision of the extraction in this formalism.

beta decay↗
DOE OSTI2024

Tritium β decay and proton-proton fusion in pionless effective field theory

The Gamow-Teller and Fermi matrix elements, 〈 GT 〉 and 〈 F 〉 , respectively, for tritium β decay are calculated to next-to-leading order (NLO) in pionless effective field theory in the absence of Coulomb interactions and isospin violation giving the leading order predictions 〈 GT 〉 0 = 0.9807 and 〈 F 〉 0 = 1 . Using an experimentally determined value for the tritium β decay GT matrix element, the two-body axial current low energy constant is fixed at NLO yielding L 1 , A = 6.01 ± 2.08 fm 3 at the renormalization scale of the physical pion mass, which agrees with predictions based on naive dimensional analysis. The impact of L 1 , A on proton-proton fusion is also discussed. Finally, the consequences of Wigner-SU(4) spin-isospin symmetry are considered for the Gamow-Teller matrix element. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗