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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
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↗
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↗
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.↗
DOE OSTI2023

Uranium Nitride as LWR TRISO Fuel: Thermodynamic Modeling of U-C-N and Thermomechanics

A new concept for light water reactor (LWR) nuclear fuel has recently been proposed that utilizes pellets of tristructural isotropic (TRISO) fuel particles embedded in a SiC or zirconium alloy matrix, replacing sintered urania pellets. The fuel is considered significantly more accident tolerant than urania as the SiC matrix or coated/protected zirconium alloy and TRISO particles would be highly resistant to oxidation and fission product release under beyond-design-basis accident conditions. As is well understood, failed fuel rods release gaseous and volatile fission products, and that is greatly exacerbated as urania oxidizes from UO 2+x to higher oxidations states in failed fuel pin with resulting expansion, loss of integrity, and significant further release of fission products.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Spectroscopy of states in 136 Ba using the 138 Ba(p,t) reaction

Background: The 136 Ba isotope is the daughter nucleus in 136 Xe $ββ$ decay. It also lies in a shape transitional region of the nuclear chart, making it a suitable candidate to test a variety of nuclear models. Purpose: To obtain spectroscopic information on states in 136 Ba, which will allow a better understanding of its low-lying structure. These data may prove useful to constrain future 136 Xe → 136 Ba neutrinoless $ββ$ decay matrix element calculations. Methods: Here, a 138 Ba(p, t) reaction was used to populate states in 136 Ba up to approximately 4.6 MeV in excitation energy. The tritons were detected using a high-resolution Q3D magnetic spectrograph. A distorted wave Born approximation analysis was performed for the measured triton angular distributions. Results: 102 excited states in 136 Ba were observed, out of which 52 are reported for the first time. Definite spin-parity assignments are made for 26 newly observed states, while previously ambiguous assignments for ten other states are resolved. Together with other available data, the results are used to determine level densities in 136 Ba. These were compared with theory predictions, obtained using shell model calculations with Hamiltonians previously used for 136 Xe neutrinoless $ββ$ decay matrix element evaluations. Conclusions: The shell model predicted level densities agree reasonably well for the two Hamiltonians. However the results for theory and experiment are found to agree only at lower energies, diverging from one another for the higher lying states, with the discrepancy increasing with energy. This is presumably because of lower production cross sections for a majority of the higher-lying predicted states and the experimental limitations in resolving a large number of nearly degenerate states predicted by the theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Constraints on anomalous Higgs boson couplings to vector bosons and fermions from the production of Higgs bosons using the 𝜏⁢𝜏 final state

A study of anomalous couplings of the Higgs boson to vector bosons and fermions is presented. The data were recorded by the CMS experiment at a center-of-mass energy of pp collisions at the LHC of 13 TeV and correspond to an integrated luminosity of 138 fb −1 . The study uses Higgs boson candidates produced mainly in gluon fusion or electroweak vector boson fusion at the LHC that subsequently decay to a pair of 𝜏 leptons. Matrix-element and machine-learning techniques were employed in a search for anomalous interactions. The results are combined with those from the four-lepton and two-photon decay channels to yield the most stringent constraints on anomalous Higgs boson couplings to date. The pure 𝐶⁢𝑃-odd scenario of the Higgs boson coupling to gluons is excluded at 2.4 standard deviations. The results are consistent with the standard model predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

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

In-situ ion irradiation of a spent UO2 fuel: evolution of fission products and nanograins with radiation dose

This study investigates the behavior of fission gas bubbles, five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) and nanograins in spent uranium dioxide (UO2) 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, 5MPs and nanograins in UO2 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. The study also reveals nanograin coarsening under irradiation. 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. This study provides critical insights into the dynamic behavior of fission products and nanograins under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗