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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamic Matrix Completion Based State Estimation in Distribution Grids

The power distribution network is undergoing tremendous transformation due to an increase in the penetration of renewable energy resources and electric vehicles. These changes have resulted in greater uncertainty and dynamics in the distribution grid states. Therefore, the ability to track and monitor system states has become a critical need for accurate and timely control actions. In this paper, we propose two dynamic sparsity-based state estimation approaches for distribution systems: (1) locally weighted matrix completion (LW-MC) and (2) Bayesian matrix completion with Kalman filter prediction (BMC-KF). The performance of the proposed dynamic state estimation strategies is compared with the classic/static matrix completion (static-MC) approach using the IEEE 37 and IEEE 123 bus test systems. Finally, results indicate that BMC-KF approach outperforms both LW-MC as well as static-MC even when 30% of the measurement data is available. Computational complexity associated with both approaches is quantified.

42 ENGINEERING↗
DOE OSTI2023

Toward more accurate adiabatic connection approach for multireference wavefunctions

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the “fixed-RDM” approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

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

Thermal Hydraulic Experimental Test Article - Fiscal Year 2023 (Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is a facility that is used to develop sodium components and instrumentation as well as acquire experimental data for validation of reactor thermal hydraulic and safety analysis codes. The facility simulates nominal conditions as well as protected/unprotected loss of flow accidents in a sodium-cooled fast reactor (SFR). High fidelity distributed temperature profiles of the developed flow field may be acquired with Rayleigh backscatter based optical fiber temperature sensors. The facility was designed in partnership with systems code experts to tailor the experiment to ensure the most relevant and highest quality data for code validation. THETA is comprised of a traditional primary coolant and secondary coolant system. The primary system is submerged in the pool of sodium and consists of a pump, electrically heated core, intermediate heat exchanger, and connected piping and thermal barriers (redan). The secondary system, located outside of the sodium pool, consists of a pump, sodium to air heat exchanger, and connected piping and valves. To date a test matrix has been completed utilizing the primary system of THETA. These tests, along with computational fluid dynamics and systems code models, determined the heat transfer across the core barrel and intermediate heat exchanger outlet was too great to effectively represent scaled thermal hydraulic phenomena of a liquid metal cooled reactor. Therefore, a significant effort was made to remove the primary system from the METL 28” test vessel #4 and clean the residual sodium from the primary system to facilitate upgrades. Thermal insulation was then incorporated in the core barrel and intermediate heat exchanger outlets. The primary system was then replaced, and a series of tests were performed to assess the performance of the thermal insulation. With primary system testing and upgrades complete, the secondary system could then be brought online. The tube side of the shell-and-tube intermediate heat exchanger was installed onto the primary system flange to begin installing the secondary system. The support structure for the secondary system was then erected on the METL mezzanine alongside the THETA primary system to facilitate installation of the secondary system components (sodium-to-air heat exchanger, flowmeter and pump). The piping and expansion tank were welded into the secondary system. Non-destructive examination of the secondary system welds was completed in order to satisfy ASME B31.3 pipe code for class M process fluids. The heating system and insulation were then added to prepare the system to be filled with sodium. The ancillary electrical equipment was installed which included the pump control box, blower VFD, pipe heater control system, etc. The secondary system will be filled, and a test matrix will be completed in early FY2024.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Ab initio ultrafast spin dynamics in solids

Spin relaxation and decoherence is at the heart of spintronics and spin-based quantum information science. Currently, theoretical approaches that can accurately predict spin relaxation of general solids including necessary scattering pathways and are capable of nanosecond to millisecond simulation time are urgently needed. We present a first-principles real-time density-matrix approach based on Lindblad dynamics to simulate ultrafast spin dynamics for general solid-state systems. Through the complete first-principles descriptions of pump, probe, and scattering processes including electron-phonon, electron-impurity, and electron-electron scatterings with self-consistent electronic spin-orbit couplings, our method can directly simulate the ultrafast pump-probe measurements for coupled spin and electron dynamics over nanoseconds at any temperatures and doping levels. We first apply this method to a prototypical system GaAs and obtain excellent agreement with experiments. We found that the relative contributions of different scattering mechanisms and phonon modes differ considerably between spin and carrier relaxation processes. In sharp contrast to previous work based on model Hamiltonians, we point out that the electron-electron scattering is negligible at room temperature but becomes dominant at low temperatures for spin relaxation in n-type GaAs. We further examine ultrafast dynamics in novel spin-valleytronic materials: monolayer and bilayer WSe 2 with realistic defects. We find that spin relaxation is highly sensitive to local symmetry and chemical bonds around defects. For the bilayer WSe 2 , we identify the scattering pathways in ultrafast dynamics and determine relevant dynamical properties, essential to its utilization of unique spin-valley-layer locking effects. In conclusion, our work provides a predictive computational platform for spin dynamics in solids, which has potential for designing new materials ideal for spintronics and quantum information technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

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

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ultrafast Molecular Frame Quantum Tomography

Here, we develop and experimentally demonstrate a methodology for a full molecular frame quantum tomography (MFQT) of dynamical polyatomic systems. We exemplify this approach through the complete characterization of an electronically nonadiabatic wave packet in ammonia (NH 3 ). The method exploits both energy and time-domain spectroscopic data, and yields the lab frame density matrix (LFDM) for the system, the elements of which are populations and coherences. The LFDM fully characterizes electronic and nuclear dynamics in the molecular frame, yielding the time- and orientation-angle dependent expectation values of any relevant operator. For example, the time-dependent molecular frame electronic probability density may be constructed, yielding information on electronic dynamics in the molecular frame. In NH 3 , we observe that electronic coherences are induced by nuclear dynamics which nonadiabatically drive electronic motions (charge migration) in the molecular frame. Here, the nuclear dynamics are rotational and it is nonadiabatic Coriolis coupling which drives the coherences. Interestingly, the nuclear-driven electronic coherence is preserved over longer timescales. In general, MFQT can help quantify entanglement between electronic and nuclear degrees of freedom, and provide new routes to the study of ultrafast molecular dynamics, charge migration, quantum information processing, and optimal control schemes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Bayesian High-Rank Hankel Matrix Completion for Nonlinear Synchrophasor Data Recovery

Phasor measurement units (PMUs) provide high temporal-resolution synchrophasor measurements for power system monitoring and control. The frequent data quality issues, such as missing and bad data, prevent the incorporation of synchrophasor data in real-time operations. Most existing data-driven data recovery methods assume the power system dynamics can be approximated by a linear dynamical system, and the recovery performance degrades significantly when the power system is experiencing nonlinear dynamics during significant events. Here, this paper proposes a data-driven Bayesian nonlinear synchrophasor data recovery method (Ba-NSDR) that can recover a consecutive time period of simultaneous data losses or errors across all channels, even when the underlying system is highly nonlinear. The idea is to lift the Hankel matrix of the spatial-temporal synchrophasor data to a higher dimension such that the lifted Hankel matrix is low-rank in that space and can be processed with the kernel trick. Our proposed Bayesian method then infers the probabilistic distributions of synchrophasor from the partial observations. Some distinctive features of Ba-NSDR include an uncertainty index to measure the accuracy of the recovery result and the robustness to parameter selections. Our method is verified on both synthetic and recorded event datasets.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗
DOE OSTI2022

Recent Advances toward Efficient Calculation of Higher Nuclear Derivatives in Quantum Chemistry

In this article, we provide an overview of state-of-the-art techniques that are being developed for efficient calculation of second and higher nuclear derivatives of quantum mechanical (QM) energy. Calculations of nuclear Hessians and anharmonic terms incur high costs and memory and scale poorly with system size. Three emerging classes of methods—machine learning (ML), automatic differentiation (AD), and matrix completion (MC)—have demonstrated promise in overcoming these challenges. We illustrate studies that employ unsupervised ML methods to reduce the need for multiple Hessian calculations in dynamics simulations and those that utilize supervised ML to construct approximate potential energy surfaces and estimate Hessians and anharmonic terms at reduced cost. By extension, if electronic structure operations could be written in a manner similar to functions underlying ML methods, rapid differentiation or AD routines can be employed to inexpensively calculate higher arbitrary-order derivatives. While ML approaches are typically black-box, we describe methods such as compressed sensing (CS) and MC, which explicitly leverage problem-specific mathematical properties of higher derivatives such as sparsity and low-rank, to complete higher derivative information using only a small, incomplete sample. The three classes of methods facilitate reliable predictions of observables ranging from infrared spectra to thermal conductivity and constitute a promising way forward in accurately capturing otherwise intractable higher-order responses of QM energy to nuclear perturbations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Kaon mixing beyond the standard model with physical masses

We present nonperturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit ( m u = m d ) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with N f = 2 + 1 flavors of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass m π ∼ 139 MeV for two different lattice spacings. We include data at three lattice spacings in the range a = 0.11 – 0.07 fm and with pion masses ranging from the physical value up to 450 MeV. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretization effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalize the lattice operators nonperturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract nonperturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the MS ¯ scheme at 3 GeV are B K ≡ B 1 = 0.5240 ( 17 ) ( 54 ) , B 2 = 0.4794 ( 25 ) ( 35 ) , B 3 = 0.746 ( 13 ) ( 17 ) , B 4 = 0.897 ( 02 ) ( 10 ) and B 5 = 0.6882 ( 78 ) ( 94 ) , where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to MS ¯ . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Open system approach to neutrinos propagating in an ultralight scalar background

We examine decoherence in neutrino oscillations induced by an ultralight scalar field coupled to neutrinos. The scalar induces time- and position-dependent shifts in the neutrino mass matrix. Neutrinos sample different field configurations throughout an experimental data-taking period, which leads to damping effects in the oscillation pattern in the form of decoherence. By recasting the neutrino-scalar dynamics within the open quantum systems framework, we establish a mapping between a complete model and phenomenological decoherence approaches. We find that the parameter driving decoherence scales as L 2 / E 2 , where L is the baseline and E is the neutrino energy, as opposed to L / E typically assumed in phenomenological studies of open system approaches to neutrino oscillations.

F.T. Airoldi, Lua [Sao Paulo U.; Fermilab] (ORCID:↗
DOE OSTI2025

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An effective matrix model for dynamical end of the world branes in Jackiw-Teitelboim gravity

We study Jackiw-Teitelboim gravity with dynamical end of the world branes in asymptotically nearly AdS2 spacetimes. We quantize this theory in Lorentz signature, and compute the Euclidean path integral summing over topologies including dynamical branes. The latter will be seen to exactly match with a modification of the SSS matrix model. The resolution of UV divergences in the gravitational instantons involving the branes will lead us to understand the matrix model interpretation of the Wilsonian effective theory perspective on the gravitational theory. We complete this modified SSS matrix model nonperturbatively by extending the integration contour of eigenvalues into the complex plane. Furthermore, we give a new interpretation of other phases in such matrix models. We derive an effective W(Φ) dilaton gravity, which exhibits similar physics semiclassically. In the limit of a large number of flavors of branes, the effective extremal entropy S0,eff has the form of counting the states of these branes.

2D Gravity↗
DOE OSTI2022

Virtual texture analysis to investigate the deformation mechanisms in metal microstructures at the atomic scale

Understanding the deformation behavior of metallic materials at high strain rates requires the characterization of plasticity contributors such as twins, phase transformed regions, and dislocations. However, predicting the contributions from phase transformation and twinning relies on a complete understanding of the selection of variants for various loading orientations and the evolution of their volume fractions. This manuscript presents a new virtual texture (VirTex) analysis approach to characterize phase transformation and twinning variants in deformed microstructures generated using molecular dynamics (MD) simulations. Furthermore, the VirTex method involves the construction of a rotation matrix to calculate the angle/axis pairs and misorientation angles for each atom in the microstructure. Any changes in the orientation angle from angle/axis pairs and/or structure types are analyzed to determine the nucleation and evolution of variants in the microstructure. The study uses shock deformed single-crystal Fe, Ta, and Cu to analyze the variant selections for phase transformation or twinning or both in BCC and FCC systems. In addition, the VirTex analysis is able to characterize the phase transformation and twinning variants in nanocrystalline Fe and Ta microstructures. Besides characterizing variants, orientation mapping also provides an accelerated and on-the-fly approach for quantifying twin fractions in MD microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

In-pixel automatic threshold calibration for the CMS Endcap Timing Layer readout chip

We present the implementation and verification of an in-pixel automatic threshold calibration circuit for the CMS Endcap Timing Layer (ETL) in the High-Luminosity LHC upgrade. The discriminator threshold of the ETL readout chip (ETROC) needs to be calibrated regularly to mitigate the circuit baseline change. Traditional methods need a lot of communication through a slow control system hence are time-consuming. This paper describes an in-pixel automatic scheme with improvements in operating time and usability. In this scheme, a sample-accumulation circuit is used to measure the average discriminator output. A binary successive approximation and linear combination scan are applied to find the equivalent baseline. The actual calibration procedure has been first implemented in FPGA firmware and tested with the ETROC front-end prototype chip (ETROC0). The calibration circuit has been implemented with Triple Modular Redundancy (TMR) and verified with Single Event Effects (SEEs) simulation. A complete calibration process lasts 35 ms with a 40 MHz clock. In the worst case, the dynamic and static power consumption are estimated to be 300 uW and 10.4 uW, respectively. The circuit design, implemented in a 65 CMOS technology, will be integrated into ETROC2, the next iteration of the ETROC with a 16 x16 pixel matrix.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Magnetic excitations, nonclassicality, and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations G(r,t), a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor S(k,ω). In this work, we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in S(k,ω) going from the noninteracting limit (U = 0) to strong coupling (U→∞), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group to find S(k,ω), from which QFI is then calculated. We find that QFI grows with U. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above U=2.5 (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth W(U) and the half-chain von Neumann entanglement entropy. We compute G(r,t) by Fourier transforming S(k,ω). The results indicate a crossover in the short-time short-distance dynamics at low U characterized by ferromagnetic light-cone wavefronts, to a Heisenberg-type behavior at large U featuring antiferromagnetic light cones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by U=3. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the G(r,t) dynamics found at low U.

1-dimensional systems↗