Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Matrix”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Some properties of adiabatic blast waves in preexisting cavities

Cox and Anderson (1982) have conducted an investigation regarding an adiabatic blast wave in a region of uniform density and finite external pressure. In connection with an application of the results of the investigation to a study of interstellar blast waves in the very hot, low-density matrix, it was found that it would be desirable to examine situations with a positive radial density gradient in the ambient medium. Information concerning such situations is needed to learn about the behavior of blast waves occurring within preexisting, presumably supernova-induced cavities in the interstellar mass distribution. The present investigation is concerned with the first steps of a study conducted to obtain the required information. A review is conducted of Sedov's (1959) similarity solutions for the dynamical structure of any explosion in a medium with negligible pressure and power law density dependence on radius.

Cox, D. P.↗
DOE OSTI2022

Decoherence factor as a convolution: an interplay between a Gaussian and an exponential coherence loss

Abstract We identify and investigate the origin and nature of the transition between Gaussian and exponential forms of decoherence: the decoherence factor (that controls the time dependence of the off-diagonal terms of the density matrix expressed in the pointer basis representation) is the convolution of the Fourier transforms of the spectral density and of the overlap (between the eigenstates the environment with and without couplings to the system). Spectral density alone tends to lead to the (approximately) Gaussian decay of coherence while the overlap alone results in a (largely) exponential decay. We show that these two contributions combine as a convolution, their relative importance controlled by the strength of the system-environment coupling. The resulting decoherence factor in the strong and weak coupling limits leads to predominantly Gaussian or exponential decay, respectively, as is demonstrated with two paradigmatic examples of decoherence—a spin-bath model and the quantum Brownian motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Theory of x-ray scattering from laser-driven electronic systems

We describe, within the framework of quantum electrodynamics, an interaction between a nonresonant hard x-ray pulse and an electronic system in the presence of a temporally periodic laser field driving electron dynamics in this system. Here we apply Floquet theory to describe the laser-driven electronic system, and then obtain the scattering probability of an arbitrary nonresonant x-ray pulse from such a system employing the density-matrix formalism. We show that the scattering probability can be connected to the time-dependent electron density of the driven electronic system only under certain conditions, in particular, if the bandwidth of the probe x-ray pulse is sufficiently narrow to spectroscopically resolve transitions to different final states. A special focus is laid on application of the theory to laser-driven crystals in a strongly nonperturbative regime. We show how the time-dependent electron density of a crystal can be reconstructed from energy-resolved scattering patterns. This is illustrated by a calculation of a diffraction signal from a driven MgO crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Strongly Correlated States of Transition Metal Spin Defects: The Case of an Iron Impurity in Aluminum Nitride

We investigate the electronic properties of an exemplar transition metal impurity in an insulator, with the goal of accurately describing strongly correlated defect states. Here, we consider iron in aluminum nitride, a material of interest for hybrid quantum technologies, and we carry out calculations with quantum embedding methods, density matrix embedding theory (DMET) and quantum defect embedding theory (QDET), and with spin-flip time-dependent density functional theory (TDDFT). We show that both DMET and QDET accurately describe the ground state and low-lying excited states of the defect and that TDDFT yields photoluminescence spectra in agreement with experiments. In addition, we provide a detailed discussion of the convergence of our results as a function of the active space used in the embedding methods, thus defining a protocol to obtain converged data directly comparable with experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

The Frequency Detuning Correction and the Asymmetry of Line Shapes: The Far Wings of H2O-H2O

A far-wing line shape theory which satisfies the detailed balance principle is applied to the H2O-H2O system. Within this formalism, two line shapes are introduced, corresponding to band-averages over the positive and negative resonance lines, respectively. Using the coordinate representation, the two line shapes can be obtained by evaluating 11-dimensional integrations whose integrands are a product of two factors. One depends on the interaction between the two molecules and is easy to evaluate. The other contains the density matrix of the system and is expressed as a product of two 3-dimensional distributions associated with the density matrices of the absorber and the perturber molecule, respectively. If most of the populated states are included in the averaging process, to obtain these distributions requires extensive computer CPU time, but only have to be computed once for a given temperature. The 11-dimensional integrations are evaluated using the Monte Carlo method, and in order to reduce the variance, the integration variables are chosen such that the sensitivity of the integrands on them is clearly distinguished.

Ma, Q.↗
DOE OSTI2020

Complex Density Wave Orders and Quantum Phase Transitions in a Model of Square-Lattice Rydberg Atom Arrays

In this work, we describe the zero-temperature phase diagram of a model of a two-dimensional square-lattice array of neutral atoms, excited into Rydberg states and interacting via strong van der Waals interactions. Using the density-matrix renormalization group algorithm, we map out the phase diagram and obtain a rich variety of phases featuring complex density wave orderings, upon varying lattice spacing and laser detuning. While some of these phases result from the classical optimization of the van der Waals energy, we also find intrinsically quantum-ordered phases stabilized by quantum fluctuations. These phases are surrounded by novel quantum phase transitions, which we analyze by finite-size scaling numerics and Landau theories. Our work highlights Rydberg quantum simulators in higher dimensions as promising platforms to realize exotic many-body phenomena.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Scaling Up Data-Parallel Analytics Platforms: Linear Algebraic Operation Cases

Linear algebraic operations such as matrix manipulations form the kernel of many machine learning and other crucial algorithms. Scaling up as well as scaling out such algorithms are key to supporting large scale data analysis that require efficient processing over millions of data samples. To this end, we present, ARION, a hardware acceleration based approach for scaling-up individual tasks of Spark, a popular data-parallel analytics platform. We support both linear algebraic operations of between two dense matrices, and between sparse and dense matrices in distributed environments. ARION provides a flexible control of acceleration according to matrix density, along with efficient scheduling based on runtime resource utilization. We demonstrate the benefit of our approach for general matrix multiplication operations over large matrices with up to four billion elements by using Gramian matrix computation that is commonly used in machine learning. Experiments show that our approach achieves more than 2× and 1.5× end-to-end performance speedups for dense and sparse matrices, respectively, and up to 57.04× faster computation compared to MLlib, a state of the art Spark-based implementation. This work is sponsored in part by the NSF under the grants: CNS-1565314, CNS-1405697, and CNS-1615411. The manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan). This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725.

Xu, Luna↗
DOE OSTI2025

Instability of Nagaoka State and Quantum Phase Transition via Kinetic Frustration Control

We investigate the Nagaoka-Thouless (NT) ferromagnetic instability in the strongly interacting t-t' Hubbard model by continuously breaking particle-hole symmetry on a tunable square-triangular lattice geometry. We use an analytic approach to show that the fully spin-polarized state becomes unstable to a metastable spin-polaron when the kinetic frustration t'/t exceeds a critical, dimension-dependent value. Large-scale density matrix renormalization group (DMRG) simulations reveal a quantum phase transition from the NT ferromagnet to a spiral spin-density wave, which evolves continuously into the Haerter-Shastry antiferromagnet in the large-frustration limit. Remarkably, this transition remains robust at low but finite hole density, making it accessible in cold-atom and moiré Hubbard platforms under strong interactions. A variational analysis further captures the instability mechanism at finite density via frustration-induced magnon band deformation.

FOS: Physical sciences↗
DOE OSTI2017

Assessment of the Tao-Mo nonempirical semilocal density functional in applications to solids and surfaces

Recently, Tao and Mo developed a semilocal exchange-correlation density functional. The exchange part of this functional is derived from a density-matrix expansion corrected to reproduce the fourth-order gradient expansion of the exchange energy in the slowly-varying-density limit, while the correlation part is based on the Tao-Perdew-Staroverov-Scuseria (TPSS) correlation functional, with a modification for the low-density limit. In this work, the Tao-Mo (TM) functional is assessed by computing various properties of solids and jellium surfaces. This includes 22 lattice constants and bulk moduli, 30 band gaps, seven cohesive energies, and jellium surface exchange and correlation energies for the density parameter r s in the range from 2 to 3 bohr. Our calculations show that the TM approximation can yield consistently high accuracy for most properties considered here, with mean absolute errors (MAEs) of 0.025 Å for lattice constants, 7.0 GPa for bulk moduli, 0.08 eV/atom for cohesive energies, and 35 erg / cm 2 for surface exchange-correlation energies. The MAE in band gaps is larger than that of TPSS, but slightly smaller than the errors of the local spin-density approximation, Perdew-Burke-Ernzerhof generalized gradient approximation, and revised TPSS. However, band gaps are still underestimated, particularly for large-gap semiconductors, compared to the Heyd-Scuseria-Ernzerhof nonlocal screened hybrid functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Single-layer tensor network study of the Heisenberg model with chiral interactions on a kagome lattice

Here, we study the antiferromagnetic kagome Heisenberg model with an additional scalar-chiral interaction by using the infinite projected entangled-pair state (iPEPS) ansatz. We discuss in detail the implementation of the optimization algorithm in the framework of the single-layer tensor network based on the corner transfer matrix technique. Our benchmark based on the full-update algorithm shows that the single-layer algorithm is stable, which leads to the same level of accuracy as the double-layer ansatz but with much less computation time. We further apply this algorithm to study the nature of the kagome Heisenberg model with a scalar-chiral interaction by computing the bond dimension scaling of magnetization, bond energy difference, chiral order parameter, and correlation length. In particular, we find that for strong chiral coupling the correlation length, which is extracted from the transfer matrix, saturates to a finite value for large bond dimension, representing a gapped spin-liquid state. Further comparison with density matrix renormalization group results supports that our iPEPS faithfully represents the time-reversal symmetry breaking chiral state. Our iPEPS simulation results shed new light on constructing PEPS for describing gapped chiral topological states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

AGR-1, AGR-2 and AGR-3/4 Dimensional Change Data Analysis

A series of Advanced Gas Reactor (AGR) experiments have been completed in the Advanced Test Reactor at Idaho National Laboratory in support of qualification and development of tristructural isotropic fuel. Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder. These capsules are instrumented with thermocouples embedded in the graphite, enabling temperature control. The fuel compacts are composed of fuel particles surrounded by a graphitic A3 matrix material. Dimensional change in AGR fuel compacts is vital because the swelling or shrinkage affects the size of the gas gaps that are used to control temperatures. Analysis of dimensional change in the AGR fuel compacts is needed to establish the variables directly relating to compact shrinkage. The variables initially identified for consideration were matrix density, compact density, fuel packing fraction, uranium loading, fuel particle diameter, cumulative fast neutron fluence, and volume average time average fuel temperature. In addition to the data available from the AGR experiments, the analysis included specimens formed from the same A3 matrix material used in Advanced Graphite Creep (AGC) experiments, which provide graphite creep data during irradiation for design and licensing purposes. The primary purpose of including the AGC specimens was to encompass dimensional behavior at zero packing fraction, zero uranium loading, and zero particle diameter. All possible combinations of first-order variable regressions were considered in the analysis. The study focused on identifying the best regression models for percent change in diameter, length, and volume. Bootstrap analysis was used to ensure the resulting regression models were robust and well-performing. The variables identified as very significant in predicting change in one or more dimensions (diameter, length, and volume) are volume average time average temperature, fast fluence, compact density, packing fraction and fuel particle diameter. Due to the presence of confounding effects between several variables, interpretation of these results is equivocal; the use of multiple statistical tests provides additional confidence in the conclusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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 OSTI2021

Precursor of pair-density wave in doping Kitaev spin liquid on the honeycomb lattice

Abstract We study the effects of doping the Kitaev model on the honeycomb lattice where the spins interact via the bond-directional interaction J K , which is known to have a quantum spin liquid as its exact ground state. The effect of hole doping is studied within the t - J K model on a three-leg cylinder using density-matrix renormalization group. Upon light doping, we find that the ground state of the system has a dominant quasi-long-range charge-density-wave correlations but short-range single-particle correlations. In the pairing channel, the even-parity superconducting correlation is dominant with d -wave-like symmetry, which oscillates in sign as a function of separation with a period equal to that of the spin-density wave and two times the charge-density wave. Although these correlations fall rapidly (possibly exponentially) at long distances, this is never-the-less the example where a pair-density wave is the leading instability in the pairing channel on the honeycomb lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Time- and momentum-resolved tunneling spectroscopy of pump-driven nonthermal excitations in Mott insulators

We present a computational technique to calculate time- and momentum-resolved nonequilibrium spectral density of correlated systems using a tunneling approach akin to scanning tunneling spectroscopy. The important difference is that our probe is extended, basically a copy of the sample, allowing one to extract the momentum information of the excitations. We illustrate the method by measuring the spectrum of a Mott-insulating extended Hubbard chain after a sudden quench with the aid of time-dependent density-matrix renormalization-group calculations. In this work, we demonstrate that the system realizes a nonthermal state that is an admixture of spin and charge density wave states, with corresponding signatures that are recognizable as in-gap subbands. In particular, we identify a band of excitons and one of stable antibound states at high energies that gains enhanced visibility after the pump. We do not appreciate noticeable relaxation within the time scales considered, which is attributed to the lack of decay channels due to spin-charge separation. These ideas can be readily applied to study transient dynamics and spectral signatures of correlation-driven nonequilibrium processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum Perturbation Theory Using Tensor Cores and a Deep Neural Network

In this work, time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e., dense matrix–matrix multiplications, in mixed-precision arithmetics, which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree–Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low-precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

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

Neutron elastic scattering on calcium isotopes from chiral nuclear optical potentials

We formulate microscopic neutron-nucleus optical potentials from many-body perturbation theory based on chiral two- and three-body forces. The neutron self-energy is first calculated in homogeneous matter to second order in perturbation theory, which gives the central real and imaginary terms of the optical potential. Here, the real spin-orbit term is calculated separately from the density matrix expansion using the same chiral interaction as in the self-energy. Finally, the full neutron-nucleus optical potential is derived within the improved local density approximation utilizing mean-field models consistent with the chiral nuclear force employed. We compare the results of the microscopic calculations to phenomenological models and experimental data up to projectile energies of E = 200 MeV. Experimental elastic differential scattering cross sections and vector analyzing powers are generally well reproduced by the chiral optical potential, but we find that total cross sections are overestimated at high energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Correlation matrix renormalization approximation for total-energy calculations of correlated electron systems

We generalized the commonly used Gutzwiller approximation for calculating the electronic structure and total energy of strongly correlated electron systems. In our method, the evaluation of one-body and two-body density matrix elements of the Hamiltonian is simplified using a renormalization approximation to achieve better scaling of the computational effort as a function of system size. To achieve a clear presentation of the concept and methodology, we describe the detailed formalism for a finite hydrogen system with minimal basis set. We applied the correlation matrix renormalization approximation approach to a H 2 dimer and H 8 cubic fragment with minimal basis sets, as well as a H 2 molecule with a large basis set. The results compare favorably with sophisticated quantum chemical calculations. We believe our approach can serve as an alternative way to build up the exchange-correlation energy functional for an improved density functional theory description of systems with strong electron correlations.

36 MATERIALS SCIENCE↗