Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A matrix completion algorithm for efficient calculation of quantum and variational effects in chemical reactions

This work examines the viability of matrix completion methods as cost-effective alternatives to full nuclear Hessians for calculating quantum and variational effects in chemical reactions. The harmonic variety-based matrix completion (HVMC) algorithm, developed in a previous study (https://doi.org/10.1063/5.0018326), exploits the low-rank character of the polynomial expansion of potential energy to recover, using a small sample, vibrational frequencies (square roots of nuclear Hessian eigenvalues) constituting the reaction path. Furthermore, these frequencies are essential for calculating rate coefficients using variational transition state theory with multidimensional tunneling (VTST-MT). HVMC performance is examined for four SN2 reactions and five hydrogen transfer reactions, with each H-transfer reaction consisting of at least one vibrational mode strongly coupled to the reaction coordinate. HVMC is robust and captures zero-point energies, vibrational free energies, zero-curvature tunneling, and adiabatic ground state and free energy barriers as well as their positions on the reaction coordinate. For medium to large reactions involving H-transfer, with the exception of the most complex Ir catalysis system, less than 35% of total eigenvalue information is necessary for accurate recovery of key VTST-MT observables.

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

Carbon nanotube substrates enhance SARS-CoV-2 spike protein ion yields in matrix-assisted laser desorption–ionization mass spectrometry

Nanostructured surfaces enhance ion yields in matrix-assisted laser desorption–ionization mass spectrometry (MALDI-MS). The spike protein complex, S1, is one fingerprint signature of Sars-CoV-2 with a mass of 75 kDa. Here, we show that MALDI-MS yields of Sars-CoV-2 spike protein ions in the 100 kDa range are enhanced 50-fold when the matrix–analyte solution is placed on substrates that are coated with a dense forest of multi-walled carbon nanotubes, compared to yields from uncoated substrates. Nanostructured substrates can support the development of mass spectrometry techniques for sensitive pathogen detection and environmental monitoring.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Multi-site reaction dynamics through multi-fragment density matrix embedding

The practical description of disordered chemical reactions, where the reactions involve multiple species at multiple sites, is presently a challenge using correlated electronic structure methods due to their high computational cost and steep scaling. Here, we describe the gradient theory of multi-fragment density matrix embedding theory, which potentially provides a minimal computational framework to model such processes at the correlated electron level. We present the derivation and implementation of the gradient theory, its validation on model systems and chemical reactions using density matrix embedding, and its application to a molecular dynamics simulation of proton transport in a small water cluster, a simple example of multi-site reaction dynamics.

Chemistry↗
DOE OSTI2024

Matrix-formation dynamics dictate methyl nitrite conformer abundance

Methyl nitrite has two stable conformational isomers resulting from rotation about the primary C–O–N–O dihedral angle: cis-CH 3 ONO and trans-CH 3 ONO, with cis being more stable by ∼5 kJ/mol. Here, the barrier to rotational interconversion (∼45 kJ/mol) is too large for isomerization to occur under ambient conditions. This paper presents evidence of a change in conformer abundance when dilute CH 3 ONO is deposited onto a cold substrate; the relative population of the freshly deposited cis conformer is seen to increase compared to its gas-phase abundance, measured by in situ infrared spectroscopy. We observe abundance changes depending on the identity of the bath gas (N 2 , Ar, and Xe) and deposition angle. The observations indicate that the surface properties of the growing matrix influence conformer abundance—contrary to the widely held assumption that conformer abundance in matrices reflects gas-phase abundance. We posit that differences in the angle-dependent host-gas deposition dynamics affect the growing surfaces, causing changes in conformer abundances. Quantum chemistry calculations of the binding energies between CH 3 ONO and a single bath-gas component reveal that significant energetic stabilization is not observed in 1:1 complexes of N 2 :CH 3 ONO, Ar:CH 3 ONO, or Xe:CH 3 ONO. From our results, we conclude that the growing surface plays a significant role in trapping cis-CH 3 ONO more effectively than trans-CH 3 ONO, likely because cis-CH 3 ONO is more compact. Taken together, the observations highlight the necessity for careful characterization of conformers in matrix-isolated systems, emphasizing a need for further study into the deposition dynamics and surface structure of chemically inert matrices.

Hockey, Emily K.↗
DOE OSTI2025

Interface instabilities in hafnium hydride entrained iron metal matrix composites

The chemical interactions in Fe–HfH 2 metal matrix composites (MMCs) are studied across multiple length scales to elucidate the decomposition of the parent phases and corresponding reaction zone physics during direct current sintering. Fe–HfH 2 composites were synthesized with increasing as-mixed hydride contents of Fe–25% HfH 2 , Fe–40% HfH 2 , Fe–55% HfH 2 , and Fe–70% HfH 2 (all in vol. %) to demonstrate the ability to achieve sintered MMCs with target hydride contents. Samples were probed across multiple length scales through a multi-modal workflow employing x-ray diffraction, scanning electron microscopy and segmentation analysis, and synchrotron techniques including hard x-ray fluorescence mapping and nanoprobe x-ray absorption near-edge structure measurements. Under the selected sintering temperature and pressure conditions, hydrogen evolution is seen to evolve through parallel paths: thermal decomposition from during the transformation of HfH 2 to HfH x<2 and through subsequent reaction with the Fe matrix leading to intermetallic phase formation. Specifically, HfFe and HfFe 2 intermetallic formation accelerates the release of hydrogen with a subsequent HfO 2 phase forming at grain boundaries. For this MMC, the consumption or loss of hydrogen can be considerable in compacts with initial hydride loading of 25%–40% HfH 2 approaching 83% hydrogen loss for the lower volume fraction composites. Increasing the volume fraction of HfH 2 to 70% enhanced the retained hydrogen content to 53% and attributed to the reduced interfacial area intrinsic to the increased HfH 2 loading in this MMC.

36 MATERIALS SCIENCE↗
DOE OSTI2026

4D beam matrix reconstruction in particle accelerators

Transverse beam parameters in particle accelerators are commonly described using the Twiss parameters, which are experimentally accessible yet inherently limited because they neglect correlations between different transverse coordinates. Such correlations frequently arise from uncompensated cathode magnetic fields or misaligned focusing quadrupoles, affecting beam quality and accelerator performance. To address this limitation, we propose and validate a novel diagnostic method for the complete four-dimensional (4D) transverse beam matrix. Our method involves passing the beam through a beamline comprising both conventional and skew quadrupole magnets, followed by downstream measurements of the resulting two-dimensional (2D) beam profiles. These measurements represent distinct 2D projections of the underlying 4D transverse phase–space distribution. By systematically varying quadrupole strengths, multiple independent projections of the beam phase space are obtained. We reconstruct the original 4D beam matrix from these measured projections using an optimization-based least-square fit, providing fast and robust reconstruction regardless of the specific beamline configuration. Through extensive numerical simulations and realistic particle-tracking studies, we demonstrate the diagnostic’s accuracy, robustness, and capability to achieve reconstruction uncertainties smaller than measurement errors, particularly when employing sufficient numbers of quadrupole scans. This method presents a powerful and flexible approach for comprehensive beam characterization and accelerator tuning.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

A New Galerkin Quadrature Method Not Requiring a Matrix Inverse

We derive a new Galerkin quadrature (GQ) method for S 𝑛 calculations that differs from the two methods preceding it in that a matrix inverse for an 𝑁 𝑑 × 𝑁 𝑑 matrix, where 𝑁𝑑 is the number of directions in the quadrature set, is no longer required. Galerkin quadrature methods are designed for calculations with highly anisotropic scattering. Such methods are not simply special angular quadratures but also are methods for representing the S 𝑛 scattering source that offers several advantages relative to the standard scattering source representation when highly truncated Legendre cross-section expansions must be used. Galerkin quadrature methods are also useful when the scattering is moderately anisotropic, but the quadrature being used is not sufficiently accurate for the order of the scattering source expansion that is required. Furthermore, we derive the new method and present computational results showing that its performance for two challenging problems is comparable to those of the two GQ methods that preceded it.

Galerkin quadrature↗
DOE OSTI2026

Hanford Tank Waste Matrix Impact on Ion Exchange Performance Using Crystalline Silicotitanate

The removal of radiocesium from Hanford tank waste supernate is a critical step in preparing feed for low-activity waste immobilization. This study evaluated cesium ion exchange performance using crystalline silicotitanate (CST) media in a series of tests designed to evaluate the influence of waste matrix variability on capacity and kinetics. Tank waste supernate subsampled from five Hanford double-shell tanks encompassed a range of sodium, hydroxide, nitrate, and nitrite concentrations in order to assess the impact of feed variability on the performance of the ion exchange system. Both equilibrium and dynamic ion exchange tests were conducted to quantify cesium distribution coefficients and breakthrough behavior under prototypic operating conditions. Results indicated that effective cesium capacity varied by up to a factor of five across the matrices tested, with higher sodium concentrations significantly reducing uptake. Kinetic behavior was similarly matrix-dependent, with solution viscosity contributing to a twofold variation in mass-transfer rates. These results demonstrate the strong dependence of CST ion exchange performance on waste composition and must be incorporated into predictive models for future treatment system design and optimization.

Westesen, Amy M.↗
DOE OSTI2024

Three point amplitudes in matrix theory

We compute the three graviton amplitude in the Banks-Fischler-Shenker-Susskind matrix model for M-theory. Even though the three point amplitude is determined by super Poincare invariance in eleven dimensional M-theory, it requires a non-trivial computation in the matrix model. We consider a configuration where all three gravitons carry non-zero longitudinal momentum. To simplify the problem, we compactify one additional dimension and relate the amplitude to a supersymmetric index computation. Here, we find agreement with the expected answer even at finite values of N.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Classification of fragile topology enabled by matrix homotopy

Flat bands in twisted materials have attracted considerable attention due to the emergence of correlated phases that can be associated with the non-Wannier-representable nature of its single-particle states. Specifically, these bands can exhibit a class of topology that can be nullified by the addition of trivial bands, termed fragile topology, which has required an expansion of prior classification schemes. However, existing approaches for predicting fragile topology rely on momentum-space methods, e.g., Wilson loops, presenting a fundamental challenge for using fragile topology as a predictor of correlated phases in aperiodic systems, such as incommensurate twist angles in moiré materials. Here, we develop a ℤ 2 energy-resolved topological marker for classifying fragile phases using a system’s position-space description, enabling the direct classification of finite, disordered, and aperiodic materials. By translating the physical symmetries protecting the system’s fragile topological phase into matrix symmetries of the system’s Hamiltonian and position operators, we use matrix homotopy to construct our topological marker while simultaneously yielding a quantitative measure of topological robustness. We demonstrate our framework’s effectiveness in both a low-energy tight-binding model and a continuum photonic crystal model of 𝐶 2 ⁢𝒯-symmetric systems, and find that fragile topology can both persist under strong disorder and even exhibit disorder-induced reentrant phase transitions. Our photonic crystal results also demonstrate the robustness of fragile topology, and the applicability of our approach, to heterostructures lacking a bulk spectral gap. Overall, our framework serves as an efficient tool for elucidating fragile topology, offering guidance for the prediction and discovery of correlated phases in both crystalline and aperiodic materials.

Lee, Ki Young [Sandia National Laboratories (SNL-N↗
DOE OSTI2026

Three-nucleon lepton-number-violating potentials in chiral effective field theory and their matrix elements in light nuclei

Here, we derive the three-nucleon neutrinoless double-𝛽 decay potential in a Δ-full chiral effective field theory through next-to-next-to-next-to leading order in Weinberg's power counting. The matrix elements of the resulting operators are computed in light nuclei using variational Monte Carlo with wave functions constructed from the Norfolk family of nuclear interactions. We find that three-nucleon corrections induce a modest quenching of the total nuclear matrix elements. We discuss model dependencies and the potential impact of these corrections on the sensitivity of experimental programs to probe lepton number violating parameters. These results provide a benchmark for many-body methods capable of reaching heavier nuclei of experimental interest.

Chambers-Wall, Graham [Washington University, St. ↗
DOE OSTI2026

Comment on “Nonlocal nucleon matrix elements in the rest frame”

In a recent paper, “Nonlocal nucleon matrix elements in the rest frame” [1], it was observed that the next-to-leading order calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice quantum chromodynamics (QCD) rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. We argue that perturbative QCD breaks down at such a distance scale after resumming the associated large logarithms, while the Ansatz used in the analysis there is not justified in perturbation theory. Besides, we explain the observation in Ref. [1] and demonstrate that the Ansatz fails to describe the data for 𝑧 > 0.3 fm, showing an opposite trend. Finally, although Ref. [1] proposes multiplying the Ansatz by a Gaussian correction model, which is shown to reduce the discrepancy with the data, this does not legitimize the use of perturbative QCD at such distance scales.

Lattice QCD↗
DOE OSTI2021

Random circuit block-encoded matrix and a proposal of quantum LINPACK benchmark

The LINPACK benchmark reports the performance of a computer for solving a system of linear equations with dense random matrices. Although this task was not designed with a real application directly in mind, the LINPACK benchmark has been used to define the list of TOP500 supercomputers since the debut of the list in 1993. We propose that a similar benchmark, called the quantum LINPACK benchmark, could be used to measure the whole machine performance of quantum computers. The success of the quantum LINPACK benchmark should be viewed as the minimal requirement for a quantum computer to perform a useful task of solving linear algebra problems, such as linear systems of equations. We propose an input model called the Random Circuit Block-Encoded Matrix (RACBEM), which is a proper generalization of a dense random matrix in the quantum setting. The RACBEM model is efficient to be implemented on a quantum computer and can be designed to optimally adapt to any given quantum architecture, with relying on a black-box quantum compiler. Besides solving linear systems, the RACBEM model can be used to perform a variety of linear algebra tasks relevant to many physical applications, such as computing spectral measures, time series generated by a Hamiltonian simulation, and thermal averages of the energy. We implement these linear algebra operations on IBM Q quantum devices as well as quantum virtual machines, and demonstrate their performance in solving scientific computing problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Capturing electronic decoherence in quantum-classical dynamics using the ring-polymer-surface-hopping–density-matrix approach

Simulations of coupled electronic and nuclear dynamics in molecules can be quite challenging due to the involved interplay of the many degrees of freedom. Because a full quantum treatment of both electrons and nuclei is computationally very demanding, it is generally restricted to model systems or rather small molecules and short timescales. Mixed quantum-classical dynamics methods such as Tully's fewest switches surface hopping (FSSH) can be used to overcome this limitation. However, FSSH is known to poorly describe electronic coherences and decoherence phenomena. Here, we present an approach that combines FSSH with ring-polymer molecular dynamics (RPMD) in a specific way that aims to alleviate the coherence problem. Termed the ring-polymer-surface-hopping–density-matrix approach, this method uses an electronic density-matrix formulation to calculate surface hopping rates. Additionally, this incorporates decoherence effects into FSSH in a natural way by taking into account the spatial spreading of the ring polymer that mimics the width of a nuclear wave packet in each RPMD trajectory. By applying our method to Tully's one-dimensional model system, we demonstrate that this method captures a crucial decoherence mechanism that is missing in FSSH. Furthermore, our method turns out to be superior at describing electronic coherences compared with earlier attempts at combining RPMD and FSSH.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗
DOE OSTI2020

Benchmark neutrinoless double- β decay matrix elements in a light nucleus

We compute nuclear matrix elements of neutrinoless double-β decay mediated by light Majorana-neutrino exchange in the $\textit{A}$ = 6 system. The goal for this work is to benchmark two many-body approaches, the no-core shell model and the multireference in-medium similarity renormalization group (SRG). We use a SRG-evolved next-to-next-to-next-to-leading chiral order two-body potential for the nuclear interaction, and make the approximation that the isospin is conserved. We compare the results of the two approaches as a function of the cutoff on the many-body basis space. Although differences are seen in the predicted nuclear radii, the ground-state energies and neutrinoless double-β decay matrix elements produced by the two approaches show significant agreement. We discuss the implications for calculations in heavier nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Diabatic Hamiltonian matrix elements made simple

With a view to applying the generator coordinate method to large configuration spaces, we propose a simple approximate formula to compute diabatic many-body matrix elements without having to evaluate two-body interaction matrix elements. Here, the method is illustrated with two analytically solvable Hamiltonians based on the harmonic oscillator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

𝑅-matrix analysis of 22 Ne states populated in 18 O ⁡(𝛼,𝛼) resonant elastic scattering

The properties of resonances in the 18 O + 𝛼 interaction are important for analyses of astrophysical processes and for understanding of the development of the 𝛼-cluster structure with addition of neutrons to 4N nuclei. Previous high-energy resolution measurements have not been treated using the contemporary 𝑅-matrix approach. We have performed an 𝑅-matrix analysis of the resonant data [D. Powers, J. K. Bair, J. L. C. Ford, Jr., and H. B. Willard, Phys. Rev. 134, B1237 (1964)] and obtained new data on 15 levels in the excitation region of 11.7–12.5 MeV 22 Ne. It appears that the model calculations did not predict the observed tendency in the 20–22 Ne 𝛼 clustering. This paper is the first step of our program of investigation of the 22 Ne spectrum in resonant reactions.

6 ≤ A ≤ 19↗