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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Tensor renormalization group for fermions

Abstract We review the basic ideas of the tensor renormalization group method and show how they can be applied for lattice field theory models involving relativistic fermions and Grassmann variables in arbitrary dimensions. We discuss recent progress for entanglement filtering, loop optimization, bond-weighting techniques and matrix product decompositions for Grassmann tensor networks. The new methods are tested with two-dimensional Wilson–Majorana fermions and multi-flavor Gross–Neveu models. We show that the methods can also be applied to the fermionic Hubbard model in 1+1 and 2+1 dimensions.

Physics↗

DOE OSTI2025

Leveraging Hamiltonian simulation techniques to compile operations on bosonic devices

Circuit quantum electrodynamics enables the combined use of qubits and oscillator modes. Despite a variety of available gate sets, many hybrid qubit-boson (i.e. qubit-oscillator) operations are realizable only through optimal control theory, which is oftentimes intractable and uninterpretable. We introduce an analytic approach with rigorously proven error bounds for realizing specific classes of operations via two matrix product formulas commonly used in Hamiltonian simulation, the Lie–Trotter–Suzuki and Baker–Campbell–Hausdorff product formulas. We show how this technique can be used to realize a number of operations of interest, including polynomials of annihilation and creation operators, namely (a) p (a † ) q for integer p, q. We show examples of this paradigm including obtaining universal control within a subspace of the entire Fock space of an oscillator, state preparation of a fixed photon number in the cavity, simulation of the Jaynes–Cummings Hamiltonian, and simulation of the Hong-Ou-Mandel effect. This work demonstrates how techniques from Hamiltonian simulation can be applied to better control hybrid qubit-boson devices.

bosonic qubits↗

DOE OSTI2024

Test of CP-invariance of the Higgs boson in vector-boson fusion production and in its decay into four leptons

A search for CP violation in the decay kinematics and vector-boson fusion production of the Higgs boson is performed in the H → ZZ * → 4 ℓ ( ℓ = e, μ ) decay channel. The results are based on proton-proton collision data produced at the LHC at a centre-of-mass energy of 13 TeV and recorded by the ATLAS detector from 2015 to 2018, corresponding to an integrated luminosity of 139 fb –1 . Matrix element-based optimal observables are used to constrain CP-odd couplings beyond the Standard Model in the framework of Standard Model effective field theory expressed in the Warsaw and Higgs bases. Differential fiducial cross-section measurements of the optimal observables are also performed, and a new fiducial cross-section measurement for vector-boson-fusion production is provided. All measurements are in agreement with the Standard Model prediction of a CP-even Higgs boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Quantum embedding study of strain- and electric-field-induced Stark effects on the NV - center in diamond

The NV - color center in diamond has been demonstrated as a powerful nanosensor for quantum metrology due to the sensitivity of its optical and spin properties to external electric, magnetic, and strain fields. In view of these applications, we use quantum embedding to derive a many-body description of strain- and electric-field-induced Stark effects on the NV - center. Here, we quantify how strain longitudinal to the axis of NV - shifts the excited states in energy, while strain with a component transverse to the NV - axis splits the degeneracies of the 3 E and 1 E states. The largest effects are for the optically relevant 3 E manifold, which splits into E x and E y with transverse strain. From these responses we extract strain susceptibilities for the E x/y states within the quasilinear regime. Additionally, we study the many-body dipole matrix elements of the NV - and find a permanent dipole 1.6 D at zero strain, which is somewhat smaller than that obtained from recent density functional theory calculations. We also determine the transition dipole between the E x and E y and how it evolves with strain.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Mapping Rydberg States of H 2 with the Halfium R-Matrix Method

In this article, we use the Halfium R-matrix method to investigate the Rydberg states of the H 2 molecule up to n = 20, filling the gap above the low-lying bound states already calculated with configuration interaction packages. Moreover, we show that the use of Quantum Defect Theory scaling laws, allows for a comprehensive analysis of the regular patterns resulting from the coupling between Rydberg series and doubly excited states. The results should open the door for more efficient quasi-diabatization of the potential energy curves which is required for calculating cross sections and rate coefficients of the (e + H 2 + ) collisional processes, involved in the plasma modeling for fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

2 νββ spectrum in chiral effective field theory

We investigate two-neutrino double beta decay (2νββ) in chiral effective field theory. We find contributions from weak magnetism and double-weak pion-exchange at next-to-leading-order in the chiral power counting. We discuss the impact of the chiral corrections on the electron spectra and find that they should be included in analyses of 2νββ decay that aim to uncover new physics signatures in the electron spectrum. We illustrate this point by revisiting the effect of sterile neutrinos and non-standard charged interactions. We also find that the pion-exchange contributions involve nuclear matrix elements that are related to those appearing in neutrinoless double beta decay (0νββ). We investigate whether the 0νββ nuclear matrix elements can be obtained from detailed measurements of the energy spectrum of the outgoing electrons in 2νββ transitions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗

DOE OSTI2024

Coherent spin states and emergent de Sitter quasinormal modes

As a toy model for the microscopic description of matter in de Sitter space, we consider a Hamiltonian acting on the spin-j representation of SU(2). This is a model with a finite-dimensional Hilbert space, from which quasinormal modes emerge in the large-spin limit. The path integral over coherent spin states can be evaluated at the semiclassical level and from it we find the single-particle de Sitter density of states, including 1/j corrections. Along the way, we discuss the use of quasinormal modes in quantum mechanics, starting from the paradigmatic upside-down harmonic oscillator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Variance reduction in lattice QCD observables via normalizing flows

Normalizing flows can be used to construct unbiased, reduced-variance estimators for lattice field theory observables that are defined by a derivative with respect to action parameters. This work implements the approach for observables involving gluonic operator insertions in the SU(3) Yang-Mills theory and two-flavor QCD in four space-time dimensions. Variance reduction by factors of 10–60 is achieved in glueball correlation functions and in gluonic matrix elements related to hadron structure, with demonstrated computational advantages. The observed variance reduction is found to be approximately independent of the lattice volume, so volume transfer can be utilized to minimize training costs.

Abbott, Ryan [Columbia U.; MIT, Cambridge, CTP; IA↗

DOE OSTI2025

Layer-dependent spin-resolved electronic structure of ferromagnetic triple-layered ruthenate Sr4Ru3O10

High-resolution angle- and spin-resolved photoemission spectroscopy (ARPES) of the triple-layered ruthenate Sr4Ru3O10 reveals features of the electronic structure that extend our understanding of the layered strontium ruthenates. The spectra near the Fermi energy are very different from the nonmagnetic analogues Sr2RuO4 and Sr3Ru2O7 with distinct Fermi surfaces for wide electronlike minority spin bands around the zone center and narrow holelike majority spin Fermi surface contours around the zone corners. The most dramatic results are two narrow spectral peaks ∼30 meV below the Fermi level, a spin-minority holelike band at the Brillouin zone center, and a spin-majority saddle-band van Hove singularity at the zone edge, which exhibits almost 100% spin polarization at low temperature, and a strong temperature dependent coherence-incoherence crossover attributed to Hund metal correlations. Quantitative comparison of the ARPES to spin-polarized density functional theory (DFT) calculations identify the specific antibonding and nonbonding orbital origins of the narrow bands, with a prediction of different spatial localization in the central and outer layers. This is shown to be consistent with experimental ARPES multizone matrix element intensity variations, and implicates outer-layer-specific control of the in-plane metamagnetism. The renormalization of the bands relative to the mean-field DFT, the demonstration of spin-polarized oxygen bands, and of spin-minority and spin-majority band-crossing hybridization provide a more complete picture of the magnetism which displays aspects of both delocalized and local moment behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗

DOE OSTI2026

Data from a four-day long microcosm experiment addressing the destabilization of artificial mineral-associated organic matter by model root exudates embedded in a soil matrix from the Rocky Mountain Biological Laboratory (Gothic, CO, USA), 2019

This dataset provides data collected during a four-day long laboratory soil microcosm experiment testing the efficacy of root exudate-driven mineral-associated organic matter destabilization. This dataset contains four data files in comma-separate values (*.csv). The files provide the metadata and the experimental results on microbial respiration, MAOM-derived respiration, and sequential mineral-extractions. This data was used to produce the figures in Bölscher et al., 2026. The results of the experiment can be found in the open access article Bölscher et al., 2026 (https://doi.org/10.1016/j.soilbio.2026.110276). Abstract: Mineral-associated organic matter (MAOM) is often considered stable, but root exudates can destabilize MAOM via various pathways. Theory and model system studies suggest that direct MAOM destabilization by strong ligands, like oxalic acid, or reducing agents, like catechol, is more effective than indirect, microbial-mediated MAOM destabilization, stimulated by less reactive compounds like glucose. Here, we demonstrate that the presence of a soil matrix alters the efficacy of exudate-driven MAOM destabilization pathways. Glucose and catechol destabilized significantly greater amounts of MAOM from ferrihydrite and aluminum hydroxide (Al (OH)3) embedded in a soil matrix than oxalic acid. Our findings indicate that indirect, microbial-mediated MAOM destabilization may play a larger role than direct MAOM destabilization in soil environments.

Destabilization↗

DOE OSTI2025

Three-dimensional imaging of pion using lattice QCD: generalized parton distributions

In this work, we report a lattice calculation of x-dependent valence pion generalized parton distributions (GPDs) at zero skewness with multiple values of the momentum transfer −t. The calculations are based on an N f = 2 + 1 gauge ensemble of highly improved staggered quarks with Wilson-Clover valence fermion. The lattice spacing is 0.04 fm, and the pion valence mass is tuned to be 300 MeV. We determine the Lorentz-invariant amplitudes of the quasi-GPD matrix elements for both symmetric and asymmetric momenta transfers with similar values and show the equivalence of both frames. Then, focusing on the asymmetric frame, we utilize a hybrid scheme to renormalize the quasi-GPD matrix elements obtained from the lattice calculations. After the Fourier transforms, the quasi-GPDs are then matched to the light-cone GPDs within the framework of large momentum effective theory with improved matching, including the next-to-next-to-leading order perturbative corrections, and leading renormalon and renormalization group resummations. We also present the 3-dimensional image of the pion in impact-parameter space through the Fourier transform of the momentum transfer −t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

nuclear-score-maximization v1.0

This software library presents efficient and multithreaded implementations of matrix low rank approximation via column selection in C++17 code. The algorithms are described in Fornace, Mark, and Michael Lindsey. "Column and row subset selection using nuclear scores: algorithms and theory for Nystro m approximation, CUR decomposition, and graph Laplacian reduction." arXiv preprint arXiv:2407.01698 (2024). The presented methods are by-and-large ver novel, have provable approximation guarantees, multiple use-cases, and exhibit higher quality approximations on a variety of studied examples.

Fornace, Mark↗

DOE OSTI2026

Scattering Processes from Quantum Simulation Algorithms for Scalar Field Theories

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Hardy, Andrew [Toronto U.] (ORCID:0000000235817382↗

DOE OSTI2026

Quantifying the phase diagram and Hamiltonian of S = 1/2 kagome antiferromagnets: bridging theory and experiment

Spin-1/2 kagome antiferromagnets are leading candidates for realizing quantum spin liquid (QSL) ground states. While QSL ground states are predicted for the pure Heisenberg model, understanding the robustness of the QSL to additional interactions that may be present in real materials is a forefront question in the field. Here we employ large-scale density-matrix renormalization group simulations to investigate the effects of next-nearest neighbor exchange couplings J 2 and Dzyaloshinskii-Moriya interactions D, which are relevant to understanding the prototypical kagome materials herbertsmithite and Zn-barlowite. By utilizing clusters as large as XC12 and extrapolating the results to the thermodynamic limit, we precisely delineate the scope of the QSL phase, which remains robust across an expanded parameter range of J 2 and D. Direct comparison of the simulated static and dynamic spin structure factors with inelastic neutron scattering reveals the parameter space of the Hamiltonians for herbertsmithite and Zn-barlowite, and, importantly, provides compelling evidence that both materials exist within the QSL phase. These results establish a powerful convergence of theory and experiment in this most elusive state of matter.

Jiang, Shengtao [SLAC National Accelerator Laborat↗

DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗