Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix product state”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Scattering amplitudes for monopoles: pairwise little group and pairwise helicity

On-shell methods are particularly suited for exploring the scattering of electrically and magnetically charged objects, for which there is no local and Lorentz invariant Lagrangian description. In this paper we show how to construct a Lorentz-invariant S-matrix for the scattering of electrically and magnetically charged particles, without ever having to refer to a Dirac string. A key ingredient is a revision of our fundamental understanding of multi-particle representations of the Poincaré group. Surprisingly, the asymptotic states for electric-magnetic scattering transform with an additional little group phase, associated with pairs of electrically and magnetically charged particles. The corresponding “pairwise helicity” is identified with the quantized “cross product” of charges, e 1 g 2 - e 2 g 1 , for every charge-monopole pair, and represents the extra angular momentum stored in the asymptotic electromagnetic field. We define a new kind of pairwise spinor-helicity variable, which serves as an additional building block for electric-magnetic scattering amplitudes. We then construct the most general 3-point S-matrix elements, as well as the full partial wave decomposition for the 2 → 2 fermion-monopole S-matrix. In particular, we derive the famous helicity flip in the lowest partial wave as a simple consequence of a generalized spin-helicity selection rule, as well as the full angular dependence for the higher partial waves. Our construction provides a significant new achievement for the on-shell program, succeeding where the Lagrangian description has so far failed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

System Modeling and Robust Stability Region Analysis for Multi-inverters Based on VSG

As the number of inverter-interfaced distributed generators (IIDGs) increases, the decrease in the inertia of the microgrid affects the stability of microgrid voltage and current. The virtual synchronous generator (VSG) control method simulates the output characteristics of the synchronous machine to improve the inertia and damping of the microgrid. However, low-frequency oscillations can occur in systems with multi-VSGs connected in parallel, primarily when operating in island mode. Therefore, this article develops the small-signal state-space model for multi-inverters based on VSG and analyses the robust stability region (RSR) of the control parameters. Firstly, an extensible small-signal state-space model is established and simplified in the dq axis, which is a necessary preprocessing for the subsequent stability analysis. Secondly, the guardian map method based on bialternate product operation is constructed, and the RSR is established through translation mapping and rotation mapping. Therefore, the solution of RSR is transformed into the problem of identifying the Hurwitz matrix. Finally, the simulations and experiments verify that the calculated region can ensure stable operation during startup and meet the requirement to suppress low-frequency oscillations. In conclusion, the method of solving RSR proposed in this paper provides a basis for selecting controller parameters in practical engineering.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Efficient simulation of moiré materials using the density matrix renormalization group

Here we present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques - even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley, and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG and open up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moiré materials in future work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Evaluations for Nuclear Criticality Safety Program 12 C, 139 La, minor actinides, 235 U [Slides]

For light nuclei, preliminary work extends the evaluation from 6.5 MeV to ~ 10 MeV. For 139 La, the team delivered full evaluation in fast region to ORNL, including covariances. For sup>235 U, RPI data simulations, the team performed simulations and showed some improvement for neutrons below 5 MeV. Some of the changes needed for more improvement might not be supported by the current format. Some of the changes above 12 MeV to account for the angular distribution of preequilibrium neutrons require a change in the PFNS evaluation procedure.

235U re-evaluation↗
NASA NTRS2005

An Open Source Simulation System

An investigation into the current state of the art of open source real time programming practices. This document includes what technologies are available, how easy is it to obtain, configure, and use them, and some performance measures done on the different systems. A matrix of vendors and their products is included as part of this investigation, but this is not an exhaustive list, and represents only a snapshot of time in a field that is changing rapidly. Specifically, there are three approaches investigated: 1. Completely open source on generic hardware, downloaded from the net. 2. Open source packaged by a vender and provided as free evaluation copy. 3. Proprietary hardware with pre-loaded proprietary source available software provided by the vender as for our evaluation.

Slack, Thomas↗
NASA NTRS1976

Eigenvalue properties of structural mean-axis systems

The authors review some of the properties of the pseudoinverse and oblique pseudoinverse of a linear transformation T from one finite-dimensional inner-product space into another, and then use these properties and a theorem of Milne (1968), which states that the oblique pseudoinverse can be expressed in terms of a weak generalized inverse and two projection operators, in order to compute a mean-axis influence coefficient matrix for the dynamic analysis of an elastic body. Some eigenvalue invariance properties of the mean-axis structural dynamics equations are demonstrated, and on the simple example of a uniform beam, it is shown that the finite frequencies and mode shapes of the mean axis structural system are identical to the nonzero frequencies and mode shapes of the free structure.

Cavin, R. K., III↗
DOE OSTI2023

Production, isolation, and shipment of clinically relevant quantities of astatine-211: A simple and efficient approach to increasing supply

The alpha emitter astatine-211 ( 211 At) is a promising candidate for cancer treatment based on Targeted Alpha (α) Therapy (TAT). A small number of facilities, distributed across the United States, are capable of accelerating α-particle beams to produce 211 At. However, challenges remain regarding strategic methods for shipping 211 At in a form adaptable to advanced radiochemistry reactions and other uses of the radioisotope. Purpose: Our method allows shipment of 211 At in various quantities in a form convenient for further radiochemistry. Procedures: For this study, a 3-octanone impregnated Amberchrom® CG300M resin bed in a column cartridge was used to separate 211 At from the bismuth matrix on site at the production accelerator (Texas A&M) in preparation for shipping. Aliquots of 6 M HNO 3 containing up to ≈2.22 GBq of 211 At from the dissolved target were successfully loaded and retained on columns. Exempt packages (<370 MBq) were shipped to a destination radiochemistry facility, University of Texas MD Anderson Cancer Center, in the form of a convenient air-dried column. Type A packages have been shipped overnight to University of Alabama at Birmingham. Main findings: Air-dried column hold times of various lengths did not inhibit simple and efficient recovery of 211 At. Solution eluted from the column was sufficiently high in specific activity to successfully radiolabel a model compound, 4-(4,4,5,5-tetramethyl-1,3,2-dioxaborolan-2-yl)aniline (1), with 211 At. The method to prepare and ship 211 At described in this manuscript has also been used to ship larger quantities of 211 At a greater distance to University of Alabama at Birmingham. Principal conclusions: The successful proof of this method paves the way for the distribution of 211 At from Texas A&M University to research institutions and clinical oncology centers in Texas and elsewhere. Finally, use of this simple method at other facilities has the potential increase the overall availability of 211 At for preclinical and clinical studies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

On time-dependent backgrounds in 1 + 1 dimensional string theory

In perturbative string theory, one is generally interested in asymptotic observables, such as the S-matrix in flat spacetime, and boundary correlation functions in anti-de Sitter spacetime. However, there are backgrounds in which such observables do not exist. We study examples of such backgrounds in 1 + 1 dimensional string theory. In these examples, the Liouville wall accelerates and can become spacelike in the past and/or future. When that happens, the corresponding null infinity, at which the standard scattering states are defined, is shielded by the Liouville wall. We compute scattering and particle production amplitudes in these backgrounds in the region in parameter space where the wall remains timelike, and discuss the continuation of this picture to the spacelike regime. We also discuss the physics from the point of view of the dynamics of free fermions in backgrounds with a time-dependent Fermi surface.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Numerically exact configuration interaction at quadrillion-determinant scale

The combinatorial growth of configuration interaction (CI) has long limited this formally exact quantum chemistry method to only the smallest molecules. Here, we report a numerically exact CI calculation exceeding one quadrillion (10 15 ) determinants, made possible by a lossless categorical compression strategy within the small-tensor-product distributed active space (STP-DAS) framework. This approach overcomes the traditional memory bottlenecks of CI by a numerically exact compression of the wavefunction representation and reformulating the most computationally demanding matrix–vector operations. Using this method, we performed a fully relativistic CI calculation of the ground state of HBrTe with over 10 15 complex-valued determinants in just 34.5 h on 1000 computing nodes—the largest CI calculation ever reported. We further achieved fast computation for systems with hundreds of billions of determinants on only a few compute nodes. Extensive benchmarks confirm that the method retains full numerical exactness while cutting memory and computational cost by orders of magnitude. Compared to previous state-of-the-art CI calculations, this work achieves a 1000 times increase in CI space, a 10 6 -fold increase in floating-point operations performed, and a 10 6 -fold improvement in computational speed.

Computational chemistry↗
DOE OSTI2021

A Method for Evaluating and Tracking Production-Sample Test (PST) Delay Risk

This report describes a method for evaluating risk from delayed Production Sample Testing (PST), and strategies for reducing test backlogs. Risk contributions to overall risk from delayed PST include a Steady State Risk (SSR) contribution that is inherent in the qualified design and manufacturing processes, and a Non-Homogenous Production Risk (NHPR) contribution from production variation. Both definitions and evaluation techniques for these risk contributions are provided along with correlation to the standard Quality Assurance Process, QA001.2, risk ranking matrix. Methods for tracking test backlogs, and reducing backlogs are also discussed. The delayed PST evaluation method described here includes the merger of both slightly quantitative methods for tracking defect levels and the propensity for production variation with general risk management tracking tools like the likelihood vs. consequence risk ranking matrix. This merger is important because it allows more effective evaluation and communication with regards to delayed PST.

42 ENGINEERING↗
DOE OSTI2020

Pyrazine-Fused Porous Graphitic Framework-Based Mixed Matrix Membranes for Enhanced Gas Separations

Membrane-based separations can mitigate the capital- and energy-intensive challenges associated with traditional thermally driven processes. To further push the boundary of gas separations, mixed matrix membranes (MMMs) have been extensively exploited; however, identifying an optimal nanofiller to boost the separation performance of MMMs beyond Robeson permeability-selectivity upper bounds remains an ongoing challenge. In this study, a new class of MMMs based on pyrazine-fused crystalline porous graphitic frameworks (PGFs) is reported. At a loading of 6 wt % PGFs, the MMMs surpass the current H 2 /CH 4 Robeson upper bound, ideally suited for applications such as H 2 regeneration. In addition, the fabricated MMMs exhibit appealing CO 2 separation performance, closely approaching the current Robeson upper bounds for CO 2 separation. Compared with the pristine polymeric membranes, the PGF-based MMMs display a record-high enhancement of gas permeability over 120% while maintaining intrinsic gas selectivities. Highlighting the crucial role of the crystallinity of nanofillers, this study demonstrates a facile and effective approach in formulating high-performance MMMs, complementing state-of-the-art membrane formation processes. The design principles open the door to energy-efficient separations of gas mixtures with enhanced productivity compatible with the current membrane manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A solid-state extrusion pathway extensible to upcycling aluminum scrap and end-of-life permanent magnets into ductile Al-6061/SmCo 5 structural magnetic composites

Permanent magnets offer strong potential for multifunctional structural components, but their inherent brittleness limits direct integration. Here, SolidStir® Extrusion (SSE) is used to embed SmCo5 powder into AA6061 in the solid state, and the resulting microstructure, interfacial chemistry, mechanical properties, and magnetic response are systematically characterized. SSE fragments the initial SmCo 5 powder into predominantly fine particles dispersed within an ultrafine Al matrix, with a volume fraction of ∼ 4–6 % in this proof-of-concept study. Interfacial reactions form a continuous Al–Co reaction layer, with local atomic-resolution indexing consistent with Al 13 Co 4 , together with Sm–Si-rich reaction products, while the matrix contains Sm–Co–Si-rich precipitates and Mg 2 Si, producing features from nanometers to micrometers. Compared with an Al extrudate, this architecture increases yield strength by ∼ 42 % and tensile strength by ∼ 31 %, while preserving uniform elongation, and it sustains work hardening in compression from 150 to 400 °C. The composite retained hard-magnetic behavior, with axial magnetization of 12.7 A·m 2 kg −1 at 6 T and intrinsic coercivity of 9.6 kOe. These findings demonstrate that SSE can produce Al-based load-bearing magnetic composites while preserving useful SmCo5 functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Macroscopic instructions vs microscopic operations in quantum circuits

In many experiments on microscopic quantum systems, it is implicitly assumed that when a macroscopic procedure or “instruction” is repeated many times – perhaps in different contexts – each application results in the same microscopic quantum operation. But in practice, the microscopic effect of a single macroscopic instruction can easily depend on its context. If undetected, this can lead to unexpected behavior and unreliable results. Here, we design and analyze several tests to detect context-dependence. They are based on invariants of matrix products, and while they can be as data intensive as quantum process tomography, they do not require tomographic reconstruction, and are insensitive to imperfect knowledge about the experiments. We also construct a measure of how unitary (reversible) an operation is, and show how to estimate the volume of physical states accessible by a quantum operation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Unifying Combinatorial and Graphical Methods in Artificial Intelligence

Recently, a new graph Laplacian, called the inner product Laplacian, was introduced which generalizes many existing Laplacians, including the normalized and combinatorial Laplacian and their weighted variants. The key observation behind the inner product Laplacian is that by defining appropriate inner product spaces on the vertices and edges, the standard Laplacians can be recovered as Hodge Laplacians over the simplicial complex formed by the edges and vertices. These inner product spaces form a natural way to incorporate non-combinatorial information into the definition of a domain-specific Laplacian. In particular, in contrast to current domain-specific weighting schemes which rely solely on edge weights, information regarding the similarity of non-adjacent vertices and arbitrary pairs of edges can be effectively incorporated into the Laplacian. In order to illustrate this approach we consider the problem of calculating the potential energy of an atomistic configuration using Graph Neural Networks. In comparison with start-of-the-art approaches, such as SchNet, our approach replaces a learned (via auto-encoder) representation of the atom types with an inner product space on atoms based on scientific knowledge (e.g., electronegativity). We will illustrate how this approach captures key chemical properties of the molecules and compare the energy calculations with state-of-the-art neural network approaches. However, to compute the resulting Laplacian involves a mixture of sparse and dense matrix computation and yields a dense matrix as the basis for the graph convolution. This dense convolutional kernel necessitates moving away from the standard message passing framework for graph neural networks and increases the computational cost of applying the kernel. In order to mitigate these costs we investigate means of leveraging the mixed sparse and dense computations to reduce the overall computational cost and how these approaches can be automatically transferred to energy efficient hardware (e.g., field programmable gate arrays (FPGAs)).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

An iterative technique to stabilize a linear time invariant multivariable system with output feedback

An iterative procedure for determining the constant gain matrix that will stabilize a linear constant multivariable system using output feedback is described. The use of this procedure avoids the transformation of variables which is required in other procedures. For the case in which the product of the output and input vector dimensions is greater than the number of states of the plant, general solution is given. In the case in which the states exceed the product of input and output vector dimensions, a least square solution which may not be stable in all cases is presented. The results are illustrated with examples.

Sankaran, V.↗
DOE OSTI2025

In-situ ion irradiation of fission products in a spent UO 2 fuel

This study investigates the behavior of fission gas bubbles and five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) in spent uranium dioxide (UO 2 ) fuel under various ion irradiation doses, temperatures, and flux conditions. Utilizing in-situ ion irradiation and advanced transmission electron microscopy, we analyzed the evolution of fission gas bubbles and 5MPs in UO 2 samples from the Belgium Reactor 3 (BR-3). Our findings reveal significant shrinkage of fission gas bubbles and 5MPs with increasing irradiation dose, accompanied by a decrease in pair density. We demonstrate that ion irradiation induces a homogeneous re-solution process where individual atoms are ejected from bubbles and precipitates, leading to their dissolution and subsequent re-precipitation in the matrix. In conclusion, this study provides critical insights into the dynamic behavior of fission products under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.

Fission products↗
DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗