Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Field redefinitions and infinite field anomalous dimensions

Field redefinitions are commonly used to reduce the number of operators in the Lagrangian by removing redundant operators and transforming to a minimal operator basis. We give a general argument that such field redefinitions, while leaving the S-matrix invariant and consequently finite, lead not only to infinite Green’s functions, but also to infinite field anomalous dimensions γ ϕ . These divergences cannot be removed by counterterms without reintroducing redundant operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Unveiling the sea: universality of the transverse momentum dependent quark distributions at small x

Within the Colour Glass Condensate effective theory, we demonstrate that back-to-back dijet correlations in dilute-dense collisions involving a small-quark from the nuclear target can be factorised in terms of universal transverse momentum dependent distributions (TMDs) for the sea quarks. Two building blocks are needed to construct all these TMDs at the operator level: the sea quark TMD operator which appears in semi-inclusive Deep-Inelastic Scattering (SIDIS) or in the Drell-Yan process and the elastic-matrix for a quark-antiquark dipole. Compared to SIDIS, the saturation effects are stronger for dijet production in forward proton-nucleus collisions, due to additional scattering in the initial and final state, effectively resulting in a larger value for the nuclear saturation momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Mitigation of DESI fiber assignment incompleteness effect on two-point clustering with small angular scale truncated estimators

We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We validate this approach on simulations reproducing the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) with and without fiber assignment. We show that we recover unbiased cosmological constraints using small angular scale truncated estimators from simulations with fiber assignment incompleteness, with respect to standard estimators from complete simulations. Additionally, we present an approach to remove the sensitivity of the fits to high k modes in the theoretical power spectrum, by applying a transformation to the data vector and window matrix. We find that our method efficiently mitigates the effect of fiber assignment incompleteness in two-point correlation function and power spectrum measurements, at low computational cost and with little statistical loss.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Utilizing integrated neutron diffraction and elastoplastic self-consistent crystal plasticity model to quantitatively assess the strengthening mechanism in Al–12.5Ce and Al–12.5Ce–0.4Mg alloys

An integrated in-situ neutron diffraction and elastic plastic self-consistent crystal plasticity (EPSC-CP) modeling scheme is performed on a binary Al–12Ce alloy and a ternary Al–12Ce–0.4Mg alloys. Using this scheme, the constitutive parameters, i.e. elastic constants and slip system parameters of individual phases can be calibrated which can be used in microstructure-based CP models to predict materials performance. From this study, it is shown that the elastic constants of Al 11 Ce 3 intermetallics calculated from density function theory calculation in the literature are rather accurate. When applied to the EPSC-CP model, the lattice strains of both the binary and ternary alloys are correctly predicted as compared with experiments, and large lattice strain differences between Al (100) plane and Al 11 Ce 3 (010) plane are demonstrated. The slip system parameters calibrated by the scheme shows that the addition of 0.4 wt% Mg in the alloy has little influence on the critical resolved shear stress of initial dislocation glide in the Al matrix which caused plastic yield in the material. This can be explained by the very dilute Mg solute content in the Al solid solution, causing large spacing of Al–Mg lattice misfit sites and little impact on resistance of dislocation glide at initial yield. The 0.4 wt% Mg addition, on the other hand, has a large influence on the hardening term in the slip system parameters, indicating those Al–Mg misfit sites do help dislocation accumulation during the deformation. The impact of dilute Mg addition on the Al slip system parameters is also reflected in the flow behavior of the ternary alloy: little impact on the yield stress, but a large impact on working hardening and tensile strength of the materials which is consistent with the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Analysis and Mitigation of Cascading Failures Using a Stochastic Interaction Graph with Eigen-analysis

In studies on complex network systems using graph theory, eigen-analysis is typically performed on an undirected graph model of the network. However, when analyzing cascading failures in a power system, the interactions among failures suggest the need for a directed graph beyond the topology of the power system to model directions of failure propagation. To accurately quantify failure interactions for effective mitigation strategies, this paper proposes a stochastic interaction graph model and associated eigen-analysis. Different types of modes on failure propagations are defined and characterized by the eigenvalues of a stochastic interaction matrix, whose absolute values are unity, zero, or in between. Finding and interpreting these modes helps identify the probable patterns of failure propagation, either local or widespread, and the participating components based on eigenvectors. Then, by lowering the failure probabilities of critical components highly participating in a mode of widespread failures, cascading can be mitigated. Here, the validity of the proposed stochastic interaction graph model, eigen-analysis and the resulting mitigation strategies is demonstrated using simulated cascading failure data on an NPCC 140-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A numerical study of bounds in the correlations of fractional quantum Hall states

We numerically compute the guiding center static structure factor \bar S(k) S ‾ ( k ) of various fractional quantum Hall (FQH) states to \mathcal{O}(k\ell)^6 𝒪 ( k ℓ ) 6 where k k is the wavenumber and \ell ℓ is the magnetic length. Employing density matrix renormalization group on an infinite cylinder of circumference L_y L y , we study the two-dimensional limit using L_y/\xi \gg 1 L y / ξ ≫ 1 , where \xi ξ is the correlation length. The main findings of our work are: 1) the ground states that deviate away from the ideal conformal block wavefunctions, do not saturate the Haldane bound, and 2) the coefficient of O(k\ell)^6 O ( k ℓ ) 6 term appears to be bounded above by a value predicted by field theories proposed in the literature. The first finding implies that the graviton mode is not maximally chiral for experimentally relevant FQH states.

Kumar, Prashant (ORCID:0000000258002768)↗
DOE OSTI2025

Probing Nanorod Assembly and Dynamics in Polymer Nanocomposites in Equilibrium and Shear

Coarse-grained molecular dynamics simulations are used to examine the structure and dynamics of nanorod assemblies in polymer melts under equilibrium and simple shear. We show that as the concentration of nanorods increases, there is a transition from an isotropic phase to a two-phase region in which the nanorods phase separate into a dilute phase and dense bundles of hexagonally packed nanorods. The onset of the two-phase region is below that predicted by Onsager theory, which we attribute to an effective increase in the diameter of nanorods due to a layer of polymer bound to the rod surfaces. Equilibrium simulations show that increasing polymer chain length N enhances nanorod bundling at fixed nanorod concentration. Simulations of systems undergoing simple shear show that flow enhances nanorod alignment and bundling relative to those of equilibrium systems with similar properties. Finally, simulations reveal that increasing nanorod length enhances nanorod alignment under shear at equivalent shear rates but not at equivalent Péclet numbers. Overall, our simulations highlight that polymer matrix-nanorod attraction (i.e., chemistry), shear rate, and matrix chain length are desirable design variables to control the structure of nanorod-containing soft materials under simple shear.

diffusion↗
DOE OSTI2025

Electrolytic gold plating, stripping, and ion transport dynamics through a solid-state iodide perovskite

The pronounced electrochemical reactivity between halide perovskites and metal electrodes can introduce mobile extrinsic metal ions which can cause device instability or enable novel functionalities. Here we systematically investigate the kinetics of gold cation (Au + ) migration in indium tin oxide (ITO)/methylammonium lead triiodide (MAPbI 3 )/Au model devices under long-term potentiostatic biasing. Scanning electron microscopy (SEM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) analyses reveal that Au + ions, electrochemically generated at the Au anode, traverse the perovskite layer with diffusion coefficients on the order of 10 −11 to 10 −10 cm 2 s −1 and are subsequently reduced at the cathode as Au 0 clusters, resembling metal plating behavior in electrolytic cells and solid-state batteries during charging. Furthermore, reversing the applied bias strips the plated Au 0 , revealing reversibility suitable for bipolar resistive switching devices and providing direct evidence of the electrochemical and ionic nature of Au transport within the perovskite matrix. Quantitatively determining diffusion coefficients and ion concentrations provides foundational inputs for future drift-diffusion modelling opportunities and allows us to relate our findings to implications on long term operation of devices like photovoltaic modules. These results clearly demonstrate the solid-state electrochemical nature of perovskite devices, highlight methods to be more quantitative about ion transport properties, provide and emphasize the importance of disentangling electro-, photo-, photoelectrochemical processes for understanding device performance and unlocking new functionalities.

14 SOLAR ENERGY↗
DOE OSTI2024

Long-distance nuclear matrix elements for neutrinoless double-beta decay from lattice QCD

Neutrinoless double-beta ( 0 ν β β ) decay is a heretofore unobserved process which, if observed, would imply that neutrinos are Majorana particles. Interpretations of the stringent experimental constraints on 0 ν β β -decay half-lives require calculations of nuclear matrix elements. This work presents the first lattice quantum chromodynamics (LQCD) calculation of the matrix element for 0 ν β β decay in a multinucleon system, specifically the n n → p p e e transition, mediated by a light left-handed Majorana neutrino propagating over nuclear-scale distances. This calculation is performed with quark masses corresponding to a pion mass of m π = 806 MeV at a single lattice spacing and volume. The statistically cleaner Σ − → Σ + e e transition is also computed in order to investigate various systematic uncertainties. The prospects for matching the results of LQCD calculations onto a nuclear effective field theory to determine a leading-order low-energy constant relevant for 0 ν β β decay with a light Majorana neutrino are investigated. This work, therefore, sets the stage for future calculations at physical values of the quark masses that, combined with effective field theory and nuclear many-body studies, will provide controlled theoretical inputs to experimental searches of 0 ν β β decay. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Communication Lower Bounds and Optimal Algorithms for Symmetric Matrix Computations

In this article, we focus on the communication costs of three symmetric matrix computations: (i) multiplying a matrix with its transpose, known as a symmetric rank-k update (SYRK) (ii) adding the result of the multiplication of a matrix with the transpose of another matrix and the transpose of that result, known as a symmetric rank-2k update (SYR2K) (iii) performing matrix multiplication with a symmetric input matrix (SYMM). All three computations appear in the Level 3 Basic Linear Algebra Subroutines (BLAS) and have wide use in applications involving symmetric matrices. We establish communication lower bounds for these kernels using sequential and distributed-memory parallel computational models, and we show that our bounds are tight by presenting communication-optimal algorithms for each setting. Our lower bound proofs rely on applying a geometric inequality for symmetric computations and analytically solving constrained nonlinear optimization problems. As a result, the symmetric matrix and its corresponding computations are accessed and performed according to a triangular block partitioning scheme in the optimal algorithms.

Al Daas, Hussam [Rutherford Appleton Laboratory, D↗
DOE OSTI2024

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

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

Vibrational Spin-Orbit Coupling Contributions to Excited State Decay of Ligand-to-Ligand Charge Transfer States

Controlling excited state relaxation processes is important in a variety of photochemical and photophysical processes, including the generation of ground and excited state spin polarization for quantum information science applications. Here, we address how specific static distortions – based on vibrational spin-orbit active modes at C2v symmetry determined by group theory – in a series of low-symmetry ligand-to-ligand charge transfer complexes enable direct spin-orbit coupling contributions to T1 → S0 excited state decay. These results are used to address spin-vibronic coupling contributions to T1 → S0 decay in a high-symmetry (tBu2bpy)Pt(S,S) (tBu2bpy = 4,4’-di-tert-butyl-2,2’-bipyridine and S,S = benzene-1,2-dithiolate) ligand-to-ligand charge transfer complex with effective C2v symmetry, where T1 relaxation is both spin- and orbitally forbidden due to the direct spin-orbit coupling matrix element being zero by symmetry. Low-frequency vibrations that involve a pyridine-pyridine twisting motion within the bpy ligand generate large ∂/∂Qi values that will contribute significantly to T1 → S0 relaxation. The work advances ligand design strategies for the generation of tailored T1 → S0 relaxation rates, which can be utilized to optimize the generation of electron spin polarization in radical-elaborated ligand-to-ligand charge transfer complexes.

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

Measurement of t -channel production of single top quarks and antiquarks in pp collisions at 13 TeV using the full ATLAS Run 2 data sample

The production of single top quarks and top antiquarks via the t-channel exchange of a virtual W boson is measured in proton-proton collisions at a centre-of-mass energy of 13 TeV at the LHC using 140 fb -1 of ATLAS data. The total cross-sections are determined to be σ(tq) = 137$^{+8}_{-8}$ pb and σ($\bar{t}q$) = 84$^{+6}_{-5}$ pb for top-quark and top-antiquark production, respectively. The combined cross-section is found to be σ(tq + $\bar{t}q$) = 221$^{+13}_{-13}$ pb and the cross-section ratio is R t = σ(tq) / σ($\bar{t}q$) = 1.636$^{+0.036}_{-0.034}$. The predictions at next-to-next-to-leading-order in quantum chromodynamics are in good agreement with these measurements. The predicted value of R t using different sets of parton distribution functions is compared with the measured value, demonstrating the potential to further constrain the functions when using this result in global fits. The measured cross-sections are interpreted in an effective field theory approach, setting limits at the 95% confidence level on the strength of a four-quark operator and an operator coupling the third quark generation to the Higgs boson doublet: -0.37 < C$^{3,1}_{Qq}$ /Λ 2 < 0.06 and -0.87 < C$^{3}_{ΦQ}$ / Λ 2 < 1.42. The constraint |V tb | > 0.95 at the 95% confdence level is derived from the measured value of σ(tq + $\bar{t}q$), assuming that the Wtb interaction is a left-handed weak coupling and that |V tb | $\gg$ |V td |, |V ts |. In a more general approach, pairs of CKM matrix elements involving top quarks are simultaneously constrained, leading to confdence contours in the corresponding two-dimensional parameter spaces.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Influence of gadolinium doping on structural, optical, and electronic properties of polymeric graphitic carbon nitride

Polymeric graphitic carbon nitride (gCN) materials have received great attention in the fields of photo and electrocatalysis due to their distinct properties in metal-free systems with high physicochemical stability. Nevertheless, the activity of undoped gCN is limited due to its relatively low specific surface area, low conductivity, and poor dispersibility. Doping Gd atoms in a gCN matrix is an efficient strategy to fine-tune its catalytic activity and its electronic structure. Herein, the influence of various wt% of gadolinium (Gd) doped in melon-type carbon nitride was systematically investigated. Gadolinium-doped graphitic carbon nitride (GdgCN) was synthesized by adding gadolinium nitrate to dicyandiamide during polymerization. The X-ray diffraction (XRD) and transmission electron microscopy (TEM) results revealed that the crystallinity and the morphological properties are influenced by the % of Gd doping. Furthermore, X-ray photoelectron spectroscopy (XPS) studies revealed that the gadolinium ions bonded with nitrogen atoms. Complementary density functional theory (DFT) calculations illustrate possible bonding configurations of Gd ions both in bulk material and on ultrathin melon layers and provide evidence for the corresponding bandgap modifications induced by gadolinium doping.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

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

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗