Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fundamental aspects of Spacetime and Quantum Fields

The proposal contained two goals: firstly, placing fundamental bound on thermalization in Quantum Field Theories (QFTs) and, secondly, developing our understanding of emergent spacetime from matrices through concrete models. Since the previous reporting period, in collaboration with Sean Hartnoll we have continued our study of entanglement edge modes in matrix quantum mechanics (MQM). This has resulted in two papers. The first applies our construction for the Matrix Quantum Hall system first to fuzzy sphere states known to correspond to stringy M2-branes in MQM. Entanglement in these states using machine learning methods have also been studied by Sean Hartnoll and Xizhi Han in previous work done under this grant. Our construction builds on this work, and further demonstrates how area laws on fuzzy geometries emerge from strongly coupled systems. The second paper generalizes this construction to all noncommutative geometries with curvature much larger than the noncommutativity parameter. We demonstrate that despite UV/IR mixing effects, the structure of entanglement edge mode irreducible representations is determined by the boundary area of subsystems. On manifolds without global symmetries, we have demonstrated that nonlocal effects inherent to noncommutative geometries resum into a change of frame of the metric structure, similar to the change from string frame to Einstein frame for entanglement entropies calculated in string theory. These advancements lay the groundwork for future progress in the understanding of emergent geometry from large-N theories. Using these techniques, we are currently working on applying our methods to noncom mutative geometries whose construction is not so well understood, such as the fuzzy 5-sphere. Despite their opacity these objects are quite important, as string physics in the bulk of holographic systems bears many features of noncommutative geometry. We have also laid the groundwork of applying our methods to tensor networks, one of the most powerful models for understanding how geometry emerges from entanglement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
DOE OSTI2025

professor

Professor is a tool to help you study complicated physical phenomena by providing tools to 1) fit machine learning models to 2D image arrays from simulations and 2) interactively explore these machine learning models in real time. Professor is most useful when studying ensembles of simulations. A typically workflow would look like: 1. A user is interested in how parameters A, B, C, & D influence some complicated physics 2. User setups up parameterized simulations to study ABCD and the results of these simulations to be image arrays (a 2d matrix of float32 values) 3. User runs an ensemble of simulations studying ABCD creating a dataset of image arrays 4. User runs `prof-trainer` to fit a machine learning model to learn the mapping from [A,B,C,D] to the image arrays 5. User then uses `prov-vis` to interactively explore the machine learning model in real time, gaining their insight into how those parameters influence the physics 6. Go profess your idea about ABCD!

Collis, HenryH [Lawrence Livermore National Labora↗
DOE OSTI2025

Strong CP and flavor in multi-Higgs theories

We introduce a class of multi-Higgs doublet extensions of the Standard Model that solve the strong problem with profound consequences for the flavor sector. The Yukawa matrices are constrained to have many zero entries by a “Higgs-flavor” symmetry, , that acts on Higgs and quark fields. The violation of both and occurs in the Higgs mass matrix so that, for certain choices of charges, the strong parameter is zero at tree level. Radiative corrections to are computed in this class of theories. They vanish in realistic two-Higgs doublet models with . We also construct realistic three-Higgs models with , where the one-loop results for are model-dependent. Requiring has important implications for the flavor problem by constraining the Yukawa coupling and Higgs mass matrices. Contributions to from higher-dimension operators are computed at one loop and can also be sufficiently small, although the hierarchy problem of this class of theories is worse than in the Standard Model.

Hall, Lawrence↗
DOE OSTI2025

Quantum Hamiltonian algorithms for maximum independent sets

ABSTRACT We compare two quantum Hamiltonian algorithms that address the maximum independent set problem: one based on the emergent non-Abelian gauge matrix in adiabatic evolution of an energetically isolated manifold of states; the other based on designed application of single-qubit operations. We demonstrate that they are mathematically equivalent in the sense that one is the other’s interaction picture. Despite their mathematical equivalence, our numerical simulations show significant differences between them in performance, which is explained analytically. Intriguingly, this equivalence unveils that the PXP model, recently prominent in quantum dynamics research, can be viewed as quantum diffusion over the median graph of all independent sets governed by the non-Abelian gauge matrix.

Science & Technology - Other Topics↗
DOE OSTI2021

Constructing tensor network influence functionals for general quantum dynamics

Here, we describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of its bond dimension and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation.

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

FLAG review 2024

We review lattice results related to pion, kaon, 𝐷-meson, 𝐵-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor 𝑓+⁡(0) arising in the semileptonic 𝐾 → 𝜋 transition at zero momentum transfer, as well as the decay-constant ratio 𝑓 𝐾 ⁡/𝑓 𝜋 and its consequences for the Cabibbo–Kobayashi–Maskawa (CKM) matrix elements 𝑉 𝑢⁢𝑠 and 𝑉 𝑢⁢𝑑 . We review the determination of the 𝐵 𝐾 parameter of neutral kaon mixing as well as the additional four 𝐵 parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for 𝑚 𝑐 and 𝑚 𝑏 as well as those for the decay constants, form factors, and mixing parameters of charmed and bottom mesons and baryons. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant 𝛼 𝑠 . We review the determinations of nucleon charges from the matrix elements of both isovector and flavor-diagonal axial, scalar and tensor local quark bilinears, and momentum fraction, helicity moment and the transversity moment from one-link quark bilinears. We also review determinations of scale-setting quantities. Finally, in this review we have added a new section on the general definition of the low-energy limit of the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Advanced modeling and simulation of research reactors using dynamic mode decomposition

Full text of publication follows. Due to the ever-increasing safety requirements, the current trend of nuclear reactor analysis is shifting towards high-fidelity multi-physics models, which have a very high computational cost and modelling complexity. As the cost of even a single model run makes it impossible to analyse the behaviour and performance of these models on large-scale commercial plants, it has become even more significant to provide suitable benchmarks to validate and test them extensively. In this sense, research reactors offer a promising solution for the initial validation of high-fidelity models, as they are significantly smaller than commercial reactors and their characteristics are well known. In particular, the reactors of the TRIGA family have been used to assess and validate models and methods for Generation-IV designs, as they have some similar features (such as the dominance of natural convection as cooling mechanism and the difficulties in performing sub-channel analysis using standard codes). Still, the computational requirements of high-fidelity models make them unsuitable for real-time analysis, even following their assessment on research reactors. In this sense, Model Order Reduction (MOR) techniques give an additional strategy to reduce the computational cost of high-fidelity models (whilst preserving sufficient accuracy). In particular, this work focuses on Dynamic Mode Decomposition (DMD), a non-intrusive MOR technique that aims at representing models with explicit temporal dynamics by extracting the time-varying characteristics and the governing structures based only on a set of available data, thus without needing any underlying knowledge of the governing equations. In addition, DMD also computes a low-dimensional surrogate of the dynamic matrix of the system, making it suited for stability analysis and real-time evaluations. This work focuses on the application and validation of the DMD method on the Computational Fluid-Dynamics (CFD) model TRIGA Mark II reactor, also discussing in detail the potentiality of this algorithm as an advanced modelling tool for nuclear reactor analysis. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Reinforcement Learning of Structured Stabilizing Control for Linear Systems With Unknown State Matrix

This paper delves into designing feedback control gains for a continuous-time linear quadratic regulator (LQR) problem that is constrained to certain predefined structure with unknown state matrix. We bring forth the ideas from reinforcement learning (RL) in conjunction with sufficient stability and performance guarantees in order to design these structured gains using the trajectory measurements of states and controls. Here we first formulate a model-based framework using dynamic programming (DP) to embed the structural constraint to the LQR gain computation in the continuous-time setting, and then subsequently, formulate a policy iteration RL algorithm that can alleviate the requirement of known state matrix in conjunction with maintaining the feedback gain structure. The design enables a distributed learning control design which is necessary for many large-scale cyber-physical systems. Theoretical guarantees are provided for stability and convergence of the structured reinforcement learning (SRL) algorithm. We validate our theoretical results with numerical simulations on a multi-agent networked linear time-invariant (LTI) dynamic system.

42 ENGINEERING↗
DOE OSTI2026

The deconstruction of flavor in the privately democratic Higgs sector

The Standard Model (SM) of particle physics fails to explain the observed hierarchy in fermion masses or the origin of fermion-flavor structure. We construct a model to explain these observations in the quark sector. We introduce a spectrum of new particles consisting of six of each — massive singlet vector-like quarks (VLQs), singlet scalars, and SU(2)-doublet scalars. SM quark masses are generated when the neutral components of the SU(2)-doublet scalars acquire non-zero vacuum expectation values (VEVs). We impose global symmetries to ensure that Yukawa couplings stay roughly flavor diagonal and democratic (of the same order), as well as to suppress tree-level flavor-changing neutral currents. Quark-mass hierarchy then follows from a hierarchy in scalar VEVs. The singlet scalars also acquire weak-scale VEVs. Together with the VLQs, they act as messengers between different generations of quarks in the SM. These messenger particles are responsible for generating the elements of the Cabibbo-Kobayashi-Masakawa (CKM) matrix which depend on the ratios of the singlet VEVs and VLQ masses. Constructed this way, the CKM matrix is found to be independent of the SM fermion masses. Using the measured values of the CKM matrix elements and assuming order-one couplings, we derive constraints on the masses of the VLQs and discuss prospects for probing our model in the near future.

CKM Parameters↗
DOE OSTI2020

High spatial resolution thermal conductivity mapping of SiC/SiC composites

Herein, silicon carbide (SiC) ceramic matrix composites are being investigated as a new generation of fuel cladding materials due to its higher accident tolerance compared to Zircaloy. Characterization of the thermal conductivity of SiC constituents and their interfaces and interphases in SiC/SiC composites is needed as inputs to models of cladding performance. We used time-domain thermoreflectance (TDTR) to map the thermal conductivity with a spatial resolution of 2 µm. The SiC/SiC composite is comprised of Hi-Nicalon Type S fibers in a matrix made by chemical vapor infiltration (CVI matrix). The interphase material is a pyrolytic carbon/SiC multilayer. We report thermal conductivity maps of the SiC fiber and SiC matrix at temperatures of 25 °C, 90 °C, 164 °C, and 250 °C. The fiber has a uniform and isotropic thermal conductivity of 22 W m -1 K -1 at 25 °C; the thermal conductivity of fiber is approximately independent of temperature. The matrix thermal conductivity is not uniform and varies from 50 to 120 W m -1 K -1 at 25 °C; the thermal conductivity of the matrix scales approximately with the inverse of the square root of temperature, i.e., ∝ 1/T 1/2 . The thermal conductivity of the interphase at room temperature is 6 W m -1 K -1 . The calculated spatially averaged thermal conductivity of this SiC/SiC composite at 1000°C is 29 W m -1 K -1 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Determination of Location-Specific Solidification Cracking Susceptibility for a Mixed Dissimilar Alloy Processed by Wire-Arc Additive Manufacturing

Solidification cracking is a major obstacle when joining dissimilar alloys using additive manufacturing. In this work, location-specific solidification cracking susceptibility has been investigated using an integrated computational materials engineering (ICME) approach for a graded alloy formed by mixing P91 steel and Inconel 740H superalloy. An alloy mixture of 26 wt.% P91 and 74 wt.% Inconel 740H, with high configurational and total entropy, was fabricated using wire arc additive manufacturing. Microstructure characterization revealed intergranular solidification cracks in the FCC matrix, which increased in length along with the enrichment of Nb (~27 to 56 wt.%) and Cu (~87 wt.%) in the middle and top regions. DICTRA simulations to model location-specific solidification cracking susceptibility showed that the top region with the highest cooling rate (270 K/s) has the highest solidification cracking susceptibility in comparison with the middle and bottom regions. This is in good agreement with the experimentally observed varying crack length. From Scheil simulations, it was deduced that enrichment of Nb and Cu affected the solidification range as high as ~77%, in comparison with the matrix composition. The overall solidification cracking susceptibility and freezing range was highest for the 26 wt.% P91 alloy amongst the mixed compositions between P91 steel and 740H superalloy, proving that solidification characteristics play a major role in alloy design for additive manufacturing.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Quantum entropy as a harbinger of factorizability

Deeply inelastic scattering (DIS) is a powerful probe for investigating the QCD structure of hadronic matter and testing the standard model (SM). DIS can be described through QCD factorization theorems which separate contributions to the scattering interaction arising from disparate scales — e.g ., with nonperturbative matrix elements associated with long distances and a perturbative hard scattering kernel applying to short-distance parton-level interactions. The fundamental underpinnings of factorization may be recast in the quantum-theoretic terms of entanglement, (de)coherence, and system localization in a fashion which sheds complementary light on the dynamics at work in DIS from QCD bound states. In this Letter, we propose and quantitatively test such a quantum-information theoretic approach for dissecting factorization in DIS and its domain of validity; we employ metrics associated with quantum entanglement such as a differential quantum entropy and associated Kullback-Leibler (KL) divergences in numerical tests. We deploy these methods on an archetypal quark-spectator model of the proton, for which we monitor quantum decoherence in DIS as underlying model parameters are varied. On this basis, we demonstrate quantitatively how factorization-breaking effects may be imprinted on quantum entropies in a kinematic regime where leading-twist factorization increasingly receives large corrections from finite- Q 2 effects; our findings suggest potential applications of quantum simulation to QCD systems and their interactions.

Deep inelastic scattering↗
DOE OSTI2021

Island for gravitationally prepared state and pseudo entanglement wedge

We consider spacetime initiated by a finite-sized initial boundary as a generalization of the Hartle-Hawking no-boundary state. We study entanglement entropy of matter state prepared by such spacetime. We find that the entanglement entropy for large subregion is given either by the initial state entanglement or the entanglement island, preventing the entropy to grow arbitrarily large. Consequently, the entanglement entropy is always bounded from above by the boundary area of the island, leading to an entropy bound in terms of the island. The island I is located in the analytically continued spacetime, either at the bra or the ket part of the spacetime in Schwinger-Keldysh formalism. The entanglement entropy is given by an average of complex pseudo generalized entropy for each entanglement island. We find a necessary condition of the initial state to be consistent with the strong sub-additivity, which requires that any probe degrees of freedom are thermally entangled with the rest of the system. We then find a large parameter region where the spacetime with finite-sized initial boundary, which does not have the factorization puzzle at leading order, dominates over the Hartle-Hawking no-boundary state or the bra ket wormhole. Due to the absence of a moment of time reflection symmetry, the island in our setup is a generalization of the entanglement wedge, called pseudo entanglement wedge. In pseudo entanglement wedge reconstruction, we consider reconstructing the bulk matter transition matrix on A ∪ I, from a fine-grained state on A. The bulk transition matrix is given by a thermofield double state with a projection by the initial state. We also provide an AdS/BCFT model by considering EOW branes with corners. We also find the exponential hardness of such reconstruction task using a generalization of Python’s lunch conjecture to pseudo generalized entropy.

2D gravity↗
DOE OSTI2020

Tuning crack-inclusion interaction with an applied T -stress

The interaction between cracks and inclusions plays an important role in the fracture behavior of particulate composites. It is commonly recognized that an inclusion stiffer than the matrix tends to deflect an approaching crack away while a softer inclusion attracts the crack. Here, we demonstrate by analytical modeling and numerical simulations that the crack-inclusion interaction can be tuned by an applied T-stress. Under a sufficiently large compressive applied T-stress, cracks can be attracted to stiffer inclusions while repelled by softer ones, therefore reversing the conventional trend. Potential applications of this work include composite electrodes in lithium-ion batteries and hydraulic fracturing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Assessment of CASL VERA for BWR analysis and application to accident tolerant SiC/SiC channel box

Application of the Virtual Environment for Reactor Applications (VERA) to BWR analysis is assessed in this paper by comparing results to those calculated using other widely-used modeling tools, namely the U.S. Nuclear Regulatory Commission’s PARCS/PATHS and the Serpent Monte Carlo particle transport code. Additionally, VERA is used to calculate 3-D temperature and fast neutron flux distributions in silicon carbide (SiC) fiber-reinforced, SiC matrix composite (SiC/SiC) BWR channel boxes, which are being studied as an Accident Tolerant Fuel core structural material concept. We found that the code-to-code comparisons were favorable, and the SiC/SiC channel box evaluation demonstrates the many advanced modeling features VERA provides while also highlighting the non-uniformity in fast neutron flux distributions that can play a role in potential SiC/SiC channel box deformation. Traditional BWR analysis tools do not have the calculation fidelity necessary for coupled assessment of flux and temperature gradients in a SiC/SiC channel box. VERA is a state-of-the-art modeling environment that was developed to increase the safety and economic competitiveness of nuclear power through improved modeling accuracy. While VERA has already been deployed in the nuclear industry for PWR applications, the current study is a vital initial step in the extensive development, validation, and verification that VERA must go through to be useful for BWR applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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↗