Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic matrix completion”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phenomenology of the Prethermal Many-Body Localized Regime

The dynamical phase diagram of interacting disordered systems has seen substantial revision over the past few years. Theory must now account for a large prethermal many-body localized regime in which thermalization is extremely slow, but not completely arrested. We derive a quantitative description of these dynamics in short-ranged one-dimensional systems using a model of successive many-body resonances. Here, the model explains the decay timescale of mean autocorrelators, the functional form of the decay—a stretched exponential—and relates the value of the stretch exponent to the broad distribution of resonance timescales. The Jacobi method of matrix diagonalization provides numerical access to this distribution, as well as a conceptual framework for our analysis. The resonance model correctly predicts the stretch exponents for several models in the literature. Successive resonances may also underlie slow thermalization in strongly disordered systems in higher dimensions, or with long-range interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Matrix-Model Simulations Using Quantum Computing, Deep Learning, and Lattice Monte Carlo

Matrix quantum mechanics plays various important roles in theoretical physics, such as a holographic description of quantum black holes, and it underpins the only practical numerical approach to the study of complex high-dimensional supergravity theories. Understanding quantum black holes and the role of entanglement in a holographic setup is of paramount importance for the realization of a quantum theory of gravity. Moreover, a complete numerical understanding of the holographic duality and the emergence of geometric space-time features from microscopic degrees of freedom could pave the way for new discoveries in quantum information science. Euclidean lattice Monte Carlo simulations are the de facto numerical tool for understanding the spectrum of large matrix models and have been used to test the holographic duality. However, they are not tailored to extract dynamical properties or even the quantum wave function of the ground state of matrix models. Quantum computing and deep learning provide potentially useful approaches to study the dynamics of matrix quantum mechanics. If successful in the context of matrix models, these rapidly improving numerical techniques could become the new Swiss army knife of quantum gravity practitioners. In this paper, we perform the first systematic survey for quantum computing and deep-learning approaches to matrix quantum mechanics, comparing them to lattice Monte Carlo simulations. These provide baseline benchmarks before addressing more complicated problems. In particular, we test the performance of each method by calculating the low-energy spectrum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

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

Gauge-fixing quantum density operators at scale

We provide a theory, algorithms, and simulations of nonequilibrium quantum systems using a one-dimensional (1D) completely positive (CP), matrix-product (MP) density-operator (𝜌) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MP⁢𝜌 factorization naturally emerges. In this setting, we analytically and numerically examine the virtual gauge freedoms associated with the representation of quantum density operators. Based on this perspective, we simplify algorithms in certain limits to speed up the integration of the canonical-form master-equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance in the near term. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, because of loss from environmental interactions, nonlocal complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional nonequilibrium settings, beyond qubit engineering.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗
DOE OSTI2024

Flax fiber-reinforced fatty acid vitrimer biocomposite with enhanced chemical recyclability

Here, this study unveils a sustainable, easily recyclable biocomposite, leveraging the dynamic nature of covalently adaptive bonds in a vitrimer matrix. The fabrication involved a fatty acid-derived vitrimer as the polymer matrix and multi-layered, nonwoven flax mat as reinforcing scaffold. The incorporation of these fibers significantly improved the mechanical performance of the vitrimer matrix uniformly. The ester-based covalently adaptive network plays a crucial role in enabling exceptional fiber-matrix bonding, as well as recyclability. The vitrimer matrix dissolves in ethylene glycol through transesterification, facilitating complete material recovery and biocomposite recycling without compromising the original properties of the matrix and reinforcing fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Friction surface layer deposition of triple-phase Al 10 Cr 12 Fe 35 Mn 23 Ni 20 high entropy alloy: Process optimization and microstructural evolution

A high-strength Co-free triple-phase Al 10 Cr 12 Fe 35 Mn 23 Ni 20 high-entropy alloy (HEA) was successfully fabricated using Friction Surface Layer Deposition (FSLD), a bulk manufacturing method. Multiple single-layer deposits were produced by varying forging force (F) and traverse speeds (T r ) to optimize the process parameters. The optimized conditions (F = 40 kN & T r = 200 mm/min) were then applied to manufacture a scaled-up multi-layer specimen. The initial microstructure of the HEA consisted of coarse grains of the soft FCC-phase, long columnar dendrites of the hard BCC-phase, and small precipitates of the harder B2-phase within the BCC-dendrites. During FSLD, the FCC-matrix underwent continuous dynamic recrystallization due to high-temperature severe plastic deformation, forming finer equiaxed grains. Simultaneously, the BCC-dendrites fractured into smaller fragments, some of which experienced partial growth and coarsening under applied stress, resulting in an hourglass morphology. In contrast, the small B2-precipitates within the BCC-fragments dissolved during the elevated temperatures of FSLD and reprecipitated as substantially finer precipitates during continuous cooling post-FSLD. Additionally, the orientation relationships between the FCC and BCC/B2 phases were completely destroyed by the severe thermoplastic deformation inherent to FSLD. The microstructural refinements led to a substantial improvement in hardness from 177 HV to 283 HV, driven by Hall-Petch strengthening. The increased number of interfaces, including coherent BCC-B2 interfaces, potentially enhances the sink strength and radiation tolerance of the HEA, making it a promising candidate for nuclear applications. In conclusion, this study also highlights FSLD as a versatile technique for achieving tunable properties in HEAs, with detailed schematics illustrating the complex mechanisms of phase transformations during processing.

Additive Manufacturing↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spectral anomalies and broken symmetries in maximally chaotic quantum maps

Spectral statistics such as the level spacing statistics and spectral form factor (SFF) are widely expected to accurately identify “ergodicity,” including the presence of underlying macroscopic symmetries, in generic quantum systems ranging from quantized chaotic maps to interacting many-body systems. By studying various quantizations of maximally chaotic maps that break a discrete classical symmetry upon quantization, we demonstrate that this approach can be misleading and fail to detect macroscopic symmetries. Notably, the same classical map can exhibit signatures of different random matrix symmetry classes in short-range spectral statistics depending on the quantization. While the long-range spectral statistics encoded in the early time ramp of the SFF are more robust and correctly identify macroscopic symmetries in several common quantizations, we also demonstrate analytically and numerically that the presence of Berry-like phases in the quantization leads to spectral anomalies, which break this correspondence. Finally, we provide numerical evidence that long-range spectral rigidity remains directly correlated with ergodicity in the quantum dynamical sense of visiting a complete orthonormal basis.

Shou, Laura [Univ. of Maryland, College Park, MD (↗
DOE OSTI2022

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
DOE OSTI2022

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence O - Raw Data

Sequence O: Sin AOA, Parked (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the absence of rotational influences by varying blade pitch angle. This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speed ranged from 18.9 m/s to 39.3 m/s, and the yaw angle was 0° throughout the sequence. The rotor was parked with the instrumented blade fixed at 0° azimuth, and the rotor lock was installed (see Appendix A). Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the OSU wind tunnel (Reuss Ramsey et al. 1995) as well as 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation O, followed by two digits for the span location (30, 47, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗
DOE OSTI2021

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

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

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence N - Raw Data

Sequence N: Sin AOA, Rotating (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the presence of rotational influences by varying blade pitch angle. Test sequence N used an upwind, rigid turbine with a 0° cone angle. The wind speed was 15 m/s, and the yaw angle was 0° throughout the sequence. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the Ohio State University (OSU) wind tunnel (Reuss Ramsey et al. 1995) and the 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation N, followed by two digits for the span location (30, 47, 63, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, most of the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗
DOE OSTI2023

Simulating noise on a quantum processor: interactions between a qubit and resonant two-level system bath

Material defects fundamentally limit the coherence times of superconducting qubits, and manufacturing completely defect-free devices is not yet possible. Therefore, understanding the interactions between defects and a qubit in a real quantum processor design is essential. We build a model that incorporates the standard tunneling model, the electric field distributions in the qubit, and open quantum system dynamics, and draws from the current understanding of two-level system (TLS) theory. Specifically, we start with one million TLSs distributed on the surface of a qubit and pick the 200 systems that are most strongly coupled to the qubit. We then perform a full Lindbladian simulation that explicitly includes the coherent coupling between the qubit and the TLS bath to model the time dependent density matrix of resonant TLS defects and the qubit. We find that the 200 most strongly coupled TLSs can accurately describe the qubit energy relaxation time. This work confirms that resonant TLSs located in areas where the electric field is strong can significantly affect the qubit relaxation time, even if they are located far from the Josephson junction (JJ). Similarly, a strongly-coupled resonant TLS located in the JJ does not guarantee a reduced qubit relaxation time if a more strongly coupled TLS is far from the JJ. In addition to the coupling strengths between TLSs and the qubit, the model predicts that the geometry of the device and the TLS relaxation time play a significant role in qubit dynamics. Our work can provide guidance for future quantum processor designs with improved qubit coherence times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Transplatformer: translating toxicogenomic profiles between generations of platforms

Background Transcriptomic profiling technologies have advanced the analysis of biological and toxicological responses. However, substantial differences in probe design, dynamic range, gene coverage, and preprocessing pipelines across platforms introduce artifacts that limit cross-study integration and hinder the reuse of historical datasets. We aim to develop computational methods for accurate cross-platform translation to maximize the value of legacy resources. Results We present TransPlatformer a deep learning framework for translating gene expression profiles across heterogeneous toxicogenomics platforms. TransPlatformer employs a novel attention-based architecture to map high-dimensional fold-change vectors from legacy microarray technologies to current platforms. Models are trained and evaluated using DrugMatrix, spanning three technological generations. We investigate mixed-tissue, single-tissue, and cross-tissue training paradigms and benchmark performance against multilayer perceptron and matrix-completion baselines. In mixed-tissue training, TransPlatformer achieves a greater than 50% reduction in mean absolute error (0.043 vs. 0.09) and nearly doubles Pearson correlation ( ≈ 0.71 vs. 0.37) relative to baseline methods. Importantly, TransPlatformer preserves rare but biologically meaningful over- and under-expressed signals, with mean absolute error below 0.22. Single-tissue models yield further improvements for well-represented organs, such as a 10% reduction in liver mean absolute error, while underscoring the need for data augmentation strategies in low-sample tissues.ra Conclusions TransPlatformer provides an effective and scalable computational solution for cross-platform transcriptomic translation. By enabling biologically faithful harmonization of gene expression data, the proposed approach facilitates the reuse of legacy toxicogenomics datasets, enhances downstream biomarker discovery, and supports more reproducible predictive modeling in toxicology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Minimal dark 𝑆⁢𝑈⁡(2) origin of a massless Dirac neutrino

We propose a gauge-symmetry origin of a rank-two Dirac neutrino mass matrix that enforces one exactly massless neutrino, while being consistent with the oscillation data, as well as cosmological constraints. The mechanism relies on a minimal dark 𝑆⁢𝑈⁢(2) 𝐷 gauge symmetry under which one right-handed neutrinolike Weyl fermion is charged, thereby forbidding its Standard Model Yukawa coupling. Quantum consistency then fixes the minimal dark-sector completion: Cancellation of the Witten anomaly requires a second fermionic 𝑆⁢𝑈⁢(2) 𝐷 doublet, while a discrete 𝑍 4 symmetry that forbids Majorana masses allows the two dark doublets to form a vectorlike pair. This anomaly-free completion gives rise to a secluded, confining dark sector that can contain a potential dark matter candidate, linking the protected neutrino texture to dark infrared dynamics.

confinement↗
DOE OSTI2024

Transverse momentum-dependent heavy-quark fragmentation at next-to-leading order

The transverse momentum-dependent fragmentation functions (TMD FFs) of heavy (bottom and charm) quarks, which we recently introduced, are universal building blocks that enter predictions for a large number of observables involving final-state heavy quarks or hadrons. They enable the extension of fixed-order subtraction schemes to quasi-collinear limits, and are of particular interest in their own right as probes of the nonperturbative dynamics of hadronization. In this paper we calculate all TMD FFs involving heavy quarks and the associated TMD matrix element in heavy-quark effective theory (HQET) to next-to-leading order in the strong interaction. Our results confirm the renormalization properties, large-mass, and small-mass consistency relations predicted in our earlier work. We also derive and confirm a prediction for the large-z behavior of the heavy-quark TMD FF by extending, for the first time, the formalism of joint resummation to capture quark mass effects in heavy-quark fragmentation. Our final results in position space agree with those of a recent calculation by another group that used a highly orthogonal organization of singularities in the intermediate momentum-space steps, providing a strong independent cross check. As an immediate application, we present the complete quark mass dependence of the energy-energy correlator (EEC) in the back-to-back limit at $\mathcal{O}\left({\alpha}_s\right)$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗