Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Quantum Computing Algorithms and Applications for Coherent and Strongly Correlated Chemical Systems

This project advanced quantum algorithms, quantum information theory, strongly correlated electronic structure methods, molecular quantum materials, and exciton transport imaging in coherent condensed phase systems. Across the award period, the team developed new methods for open-quantum-system simulation, Hamiltonian learning, state tomography, reduced-density-matrix and contracted-quantum eigensolver approaches, and quantum diagnostics for device capability and openness.

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

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Polariton spectra under the collective coupling regime. I. Efficient simulation of linear spectra and quantum dynamics

We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein–Tavis–Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work.

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

Detecting superconductivity out of equilibrium

Recent pump-probe experiments on underdoped cuprates and similar systems suggest the existence of a transient superconducting state above Tc. This poses the question of how to reliably identify the emergence of long-range order, in particular superconductivity, out of equilibrium. We investigate this point by studying a quantum quench in an extended Hubbard model and by computing various observables, which are used to identify (quasi-)long-range order in equilibrium. Our findings imply that, in contrast to current experimental studies, it does not suffice to study the time evolution of the optical conductivity to identify superconductivity. In turn, we suggest to utilize time-resolved angle-resolved photoemission spectroscopy experiments to probe for the formation of a condensate in the two-particle channel.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Performance of wave function and Green's function methods for non-equilibrium many-body dynamics

Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime. Published by the American Physical Society 2025

Reeves, Cian C. (ORCID:0009000642581845)↗
DOE OSTI2023

Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory

The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density-matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies, relative to reference values from configuration-interaction-based calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2020

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Absolute electron density fluctuation reconstruction for two-dimensional hydrogen beam emission spectroscopy

Scrape-off layer (SOL) and edge plasma turbulence significantly contribute to the radial particle and heat transport, lowering the plasma confinement and increasing the heat load on the plasma facing components. SOL turbulence is predominantly intermittent, which manifests in the occurrence of isolated density filaments or blobs. Filaments propagate radially outward toward plasma facing components, limiting their lifetime by erosion and sputtering. To characterize this phenomenon in detail, few diagnostic techniques are available. Beam emission spectroscopy is a diagnostic capable of measuring plasma turbulence in both SOL and edge plasmas. Due to the finite lifetime of the excitation states during the beam–plasma interaction and the misalignment between the optics and the magnetic field, spatial smearing is introduced in the measurement. In this paper, a novel method is introduced to overcome this hindering effect by inverting the fluctuation response matrix on an optimally smoothed signal. We show that this method is fast and provides significantly more accurate absolute density fluctuation reconstruction than the direct inversion technique. Here, the presented method is usable for all types of beam emission diagnostics where the spatial resolution is higher than the combined smearing of the atomic physics and the observation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A Solubility Comparison of Cyclotetramethylenetetranitramine (HMX) and Potential Mock-HMX Candidates in Solvents Relevant to Slurry Formulating Polymer-Bonded Materials

Understanding solubility behavior of organic crystalline materials is critical in the development of slurry-formulated plastic bonded composites. In the field of high explosives (HE), a common method for obtaining high density, greater than 95 % theoretical maximum density, is to encapsulate HE in a polymer matrix resulting in the formation of plastic bonded explosive (PBX) prills. Prill formation is highly nuanced and requires understanding of particle size distributions, solvent, antisolvent, and binder system as well as interactions thereof. In this paper, in an effort to create mock-HE that mimics HE in density and mechanical strength, solubility of candidate molecules was compared to cyclotetramethylenetetranitramine (HMX) in binder solvents relevant to prill formation. The results indicate that the mock candidates chosen for this study have different solubilities compared to HMX in water and both organic solvents, 2-butanone and ethyl acetate, investigated. With these solubility differences known, prills to match PBXs based on HMX can be formulated, pressed and tested.

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

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

14 SOLAR ENERGY↗
DOE OSTI2021

Performance of the rigorous renormalization group for first-order phase transitions and topological phases

Expanding and improving the repertoire of numerical methods for studying quantum lattice models is an ongoing focus in many-body physics. While the density matrix renormalization group (DMRG) has been established as a practically useful algorithm for finding the ground state in one-dimensional systems, a provably efficient and accurate algorithm remained elusive until the introduction of the rigorous renormalization group (RRG) by Landau [Nat. Phys. 11, 566 (2015)1745-247310.1038/nphys3345]. In this paper, we study the accuracy and performance of a numerical implementation of RRG at first-order phase transitions and in symmetry-protected topological phases. Our study is motivated by the question of when RRG might provide a useful complement to the more established DMRG technique. In particular, despite its general utility, DMRG can give unreliable results near first-order phase transitions and in topological phases, since its local update procedure can fail to adequately explore (near-)degenerate manifolds. Furthermore, the rigorous theoretical underpinnings of RRG, meanwhile, suggest that it should not suffer from the same difficulties. We show this optimism is justified, and that RRG indeed determines well-ordered, accurate energies even when DMRG does not. Moreover, our performance analysis indicates that in certain circumstances seeding DMRG with states determined by coarse runs of RRG may provide an advantage over simply performing DMRG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Efficient simulation of low-temperature physics in one-dimensional gapless systems

Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O⁡(cN 0 ⁢log⁡β) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dimerization and spin decoupling in a two-leg Heisenberg ladder with frustrated trimer rungs

We study the antiferromagnetic spin-half Heisenberg ladder in the presence of an additional frustrating rung spin which is motivated and relevant also for the description of real two-dimensional materials such as the two-dimensional trimer magnet Ba 4 Ir 3 O 10 . We study the zero-temperature phase diagram, where we combine numerical and analytical methods into an overall consistent description. All numerical simulations are also accompanied by studies of the dynamical spin structure factor obtained via the density matrix renormalization group. Overall, we find in the regime of strong rung coupling a gapped dimerized phase related to competing symmetry sectors in Hilbert space that ultimately results in frustration-driven spin-Peierls transition. In the weak rung-coupling regime, the system is uniform, yet shows a gapped spinon continuum together with a sharp coherent low-energy branch which renders the system critical overall. In either case, the additional rung spin quickly get sidelined and nearly decouple in the regime when their bare coupling to the ladder is weakened

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Bootstrap embedding for interacting electrons in phonon coherent-state mean field

Here, we develop a Fermi–Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron–phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard–Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.

Islam, Shariful [North Carolina State University, ↗
DOE OSTI2026

Multireference Embedding and Fragmentation Methods for Classical and Quantum Computers: From Model Systems to Realistic Applications

One of the primary challenges in quantum chemistry is the accurate modeling of strong electron correlation. While multireference methods effectively capture such correlation, their steep scaling with system size prohibits their application to large molecules and extended materials. Quantum embedding offers a promising solution by partitioning complex systems into manageable subsystems. In this Review, we highlight recent advances in multireference density matrix embedding and localized active space self-consistent field approaches for complex molecules and extended materials. We discuss both classical implementations and the emerging potential of these methods on quantum computers. Here, by extending classical embedding concepts to the quantum landscape, these algorithms have the potential to expand the reach of multireference methods in quantum chemistry and materials.

Algorithms↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗