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

Predictive modelling of hydration and mechanical performance of low Ca composite cements: Possibilities and limitations from industrial perspective

We quantified hydration and performance development of composites cements containing metakaolin in the rage of 0–40% and investigated applicability of thermodynamic and micromechanical models to predict the evolution of hydrates assemblage and of performance, respectively. The reaction of metakaolin results in a formation of significant amounts of C-S-H phase for lower metakaolin content and of strätlingite for higher contents. This leads to very dense matrix with finer porosity allowing co-existence of thermodynamically incompatible phases, namely portlandite and strätlingite. Hence, thermodynamic modelling cannot be applied to correctly predict the phase assemblage for the systems with high metakaolin content. By introducing constraints into the thermodynamic model in a form of minimum phase content, the XRD-amorphous phases can be quantified. However, such constraints are only valid for the cases studied and limit the general applicability and predictive power of the model. Elastic properties and strength of mortars are influenced by the metakaolin additions in an opposite trend. It indicates that they are sensitive to different microstructural features/mechanisms. Continuum micromechanic-based model is capable to accurately predict the performance evolution based on the volume of hydrates predicted by the thermodynamic modelling. An additional step required for the accurate prediction - porosity is partitioned based on its saturation at 95% RH – is introduced and substantiated.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗
DOE OSTI2022

Mana and thermalization: Probing the feasibility of near-Clifford Hamiltonian simulation

Quantum hydrodynamics is the emergent classical dynamics governing transport of conserved quantities in generic strongly interacting quantum systems. Here, recent matrix product operator methods [1,2] have made simulations of quantum hydrodynamics in 1+1D tractable, but they do not naturally generalize to 2+1D or higher, and they offer limited guidance as to the difficulty of simulations on quantum computers. Near-Clifford simulation algorithms are not limited to one dimension, and future error-corrected quantum computers will likely be bottlenecked by non-Clifford operations. We therefore investigate the non-Clifford resource requirements for simulation of quantum hydrodynamics using mana, a resource theory of non-Clifford operations. For infinite-temperature starting states, we find that the mana of subsystems quickly approaches zero, while for starting states with energy above some threshold the mana approaches a nonzero value. Surprisingly, in each case the finite-time mana is governed by the subsystem entropy, not the thermal state mana; we argue that this is because mana is a sensitive diagnostic of finite-time deviations from canonical typicality.

1-dimensional spin chains↗
DOE OSTI2023

Spectral function of the $J_1$ – $J_2$ Heisenberg model on the triangular lattice

Spectral probes, such as neutron scattering, are crucial for characterizing excitations in quantum many-body systems and the properties of quantum materials. Among the most elusive phases of matter are quantum spin liquids, which have no long-range order even at zero temperature and host exotic fractionalized excitations with nontrivial statistics. These phases can occur in frustrated quantum magnets, such as the paradigmatic Heisenberg model with nearest- and next-nearest-neighbor exchange interactions on the triangular lattice, the so-called $J_1$ – $J_2$ model. In this work, we compute the spectral function using large-scale matrix product state simulations across the three different phases of this model's phase diagram, including a quantum spin-liquid phase at intermediate $J_2$/$J_1$. Despite a plethora of theoretical and experimental studies, the exact nature of this phase is still contested, with the dominant candidates being a gapped $\mathbb{Z}_2$, a gapless U(1) Dirac, and a spinon Fermi surface quantum spin-liquid state. We find a V-shaped spectrum near the center of the Brillouin zone ($Γ$ point), a key signature of a spinon Fermi surface, observed in prior neutron scattering experiments. However, we find a small gap near the $Γ$ point, ruling out such a phase. Furthermore, we find localized gapless excitations at the corner of the Brillouin zone boundary (K point) and the middle of the edge of the Brillouin zone boundary (M point), ruling out the gapped $\mathbb{Z}_2$ spin-liquid phase. Here, our results imply that the intermediate spin-liquid phase is a gapless $U(1)$ Dirac spin liquid, and provide clear signatures to detect this phase in future neutron scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Dissipative transport in superlattices within the Wigner function formalism

Here, we employ the Wigner function formalism to simulate partially coherent, dissipative electron transport in biased semiconductor superlattices. We introduce a model collision integral with terms that describe energy dissipation, momentum relaxation, and the decay of spatial coherences (localization). Based on a particle-based solution to the Wigner transport equation with the model collision integral, we simulate quantum electronic transport at 10 K in a GaAs/AlGaAs superlattice and accurately reproduce its current density vs field characteristics obtained in experiment.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Derivation of the Time Dependent Two Dimensional Focusing NLS Equation

We present a microscopic derivation of the two-dimensional focusing cubic nonlinear Schrödinger equation starting from an interacting N-particle system of Bosons. The interaction potential we consider is given by W{sub β}(x)=N{sup −1+2β}W(N{sup β}x) for some spherically symmetric and compactly supported potential W∈L{sup ∞}(ℝ{sup 2},ℝ). The class of initial wave functions is chosen such that the variance in energy is small. Furthermore, we assume that the Hamiltonian H{sub W{sub β,t}}=−∑{sub j=1}{sup N}Δ{sub j}+∑{sub 1≤j 2, or in trace norm, for more general external potentials. For trapping potentials of the form A(x)=C|x|{sup s} ,C>0, the condition H{sub W{sub β,t}}≥− CN can be fulfilled for a certain class of interactions W{sub β}, for all 0<β<((s+1)/(s+2)), see Lewin et al. (Proc Am Math Soc 145:2441–2454, 2017).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Some Mathematical Problems of Control of Quantum Systems

Currently, quantum technologies are actively developing; these technologies are based on quantum effects in individual quantum systems—atoms or molecules. The mathematical study of problems of control of quantum systems is of particular importance. In this paper, we consider certain problems associated with control of quantum systems: extrema of target functionals for population transfer problems and generation of unitary processes and incoherent control and generation of arbitrary density matrices for open quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Distributed out-of-memory NMF on CPU/GPU architectures

We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Entropy evolution law in a laser process

For the first time, we obtain the entropy variation law in a laser process after finding the Kraus operator of the master equation describing the laser process with the use of the entangled state representation. The behavior of entropy is determined by the competition of the gain and damping in the laser process. The evolution formula for the number of photons is also obtained. -- Highlights: •We apply the thermo-field dynamics method to the master equation of a laser. •We find the operator-sum (Kraus) representation for the density operator. •We find both the normally ordered and compact forms of ρ(t) for ρ{sub 0}=|z〉〈z|. •We find the exact expression of the laser’s entropy. •Our results reveal quantitatively how a laser beam can be generated in a laser.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Nonequilibrium steady state in open quantum systems: Influence action, stochastic equation and power balance

The existence and uniqueness of a steady state for nonequilibrium systems (NESS) is a fundamental subject and a main theme of research in statistical mechanics for decades. For Gaussian systems, such as a chain of classical harmonic oscillators connected at each end to a heat bath, and for classical anharmonic oscillators under specified conditions, definitive answers exist in the form of proven theorems. Answering this question for quantum many-body systems poses a challenge for the present. In this work we address this issue by deriving the stochastic equations for the reduced system with self-consistent backaction from the two baths, calculating the energy flow from one bath to the chain to the other bath, and exhibiting a power balance relation in the total (chain + baths) system which testifies to the existence of a NESS in this system at late times. Its insensitivity to the initial conditions of the chain corroborates to its uniqueness. The functional method we adopt here entails the use of the influence functional, the coarse-grained and stochastic effective actions, from which one can derive the stochastic equations and calculate the average values of physical variables in open quantum systems. This involves both taking the expectation values of quantum operators of the system and the distributional averages of stochastic variables stemming from the coarse-grained environment. This method though formal in appearance is compact and complete. It can also easily accommodate perturbative techniques and diagrammatic methods from field theory. Taken all together it provides a solid platform for carrying out systematic investigations into the nonequilibrium dynamics of open quantum systems and quantum thermodynamics. -- Highlights: •Nonequilibrium steady state (NESS) for interacting quantum many-body systems. •Derivation of stochastic equations for quantum oscillator chain with two heat baths. •Explicit calculation of the energy flow from one bath to the chain to the other bath. •Power balance relation shows the existence of NESS insensitive to initial conditions. •Functional method as a viable platform for issues in quantum thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Quantum dynamics at finite temperature: Time-dependent quantum Monte Carlo study

In this work we investigate the ground state and the dissipative quantum dynamics of interacting charged particles in an external potential at finite temperature. The recently devised time-dependent quantum Monte Carlo (TDQMC) method allows a self-consistent treatment of the system of particles together with bath oscillators first for imaginary-time propagation of Schrödinger type of equations where both the system and the bath converge to their finite temperature ground state, and next for real time calculation where the dissipative dynamics is demonstrated. In that context the application of TDQMC appears as promising alternative to the path-integral related techniques where the real time propagation can be a challenge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Field’s entropy in the atom–field interaction: Statistical mixture of coherent states

We study the atom–field interaction when the field is in a mixture of coherent states. We show that in this case it is possible to calculate analytically the field entropy for times of the order of twice the collapse time. Such analytical results are done with the help of numerical analysis. We also give an expression in terms of Chebyshev polynomials for power of density matrices. - Highlights: • We calculate the field entropy for times of the order of twice the collapse time. • We give a relation between powers of the density matrices of the subsystems. • Entropy operators for both subsystems are obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

On the reduced dynamics of a subset of interacting bosonic particles

Highlights: Discusses open quantum systems where system and environment are indistinguishable. Perturbative, microscopic treatment of interacting bosonic many-body systems. First-order approximation: BBGKY-hierarchy and, as a special case, Gross–Pitaevskii equation. Second-order approximation: Nonlinear master equation with interaction-induced decoherence. - Abstract: The quantum dynamics of a subset of interacting bosons in a subspace of fixed particle number is described in terms of symmetrized many-particle states. A suitable partial trace operation over the von Neumann equation of an N-particle system produces a hierarchical expansion for the subdynamics of M≤N particles. Truncating this hierarchy with a pure product state ansatz yields the general, nonlinear coherent mean-field equation of motion. In the special case of a contact interaction potential, this reproduces the Gross–Pitaevskii equation. To account for incoherent effects on top of the mean-field evolution, we discuss possible extensions towards a second-order perturbation theory that accounts for interaction-induced decoherence in form of a nonlinear Lindblad-type master equation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Benchmarking the multiconfigurational Hartree method by the exact wavefunction of two harmonically trapped bosons with contact interaction

Highlights: • Comparison of exact and numerical single-particle density matrices. • Benchmark for multiconfigurational Hartree. • Convergence for large number of orbitals. We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, as well as overestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Transformations of Na, Al, Si and Fe species in red mud during synthesis of one-part geopolymers

Reaction mechanism of one-part geopolymer formation using red mud (RM) was explored in this study. Leaching test, X-ray photoelectron spectroscopy (XPS), Mössbauer spectroscopy and nitrogen isothermal adsorption (NIA) were conducted to study transformations of Na, Al, Si, Fe species, binding energies of Al−O and Si−O, and pore structure, respectively. The alkali-thermal-activated RM could dissolve in water to release large amounts of soluble Na, Al and Si species for geopolymerization. But more dissolved Si is needed from added silica fume (SF) to form stable geopolymer gels. The binding energies of Al−O and Si−O bonds in the geopolymer increased with the polymerization of reactive SiO{sub 2} and Al{sub 2}O{sub 3}. The coordinated Fe{sup 3+} in RM played a role like Al{sup 3+} by replacing it in the aluminosilicates structure, and exerted no obvious impact on geopolymerization. The dissolved SF took part in the geopolymerization to form dense geopolymer matrices in long-term curing. - Highlights: • Large amounts of Na, Al and Si dissolved from alkali-thermal-activated RM in water. • Low Si/Al geopolymers will dissolve and release soluble Al in long-term curing. • Dissolved Si from SF participates in the polymerization of soluble Na, Al and Si. • The binding energies of Al−O and Si−O will increase after geopolymerization. • The coordinated Fe{sup 3+} replaces Al{sup 3+} in aluminosilicate structure of geopolymer.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Bayesian estimation of Karhunen–Loève expansions; A random subspace approach

One of the most widely-used statistical procedures for dimensionality reduction of high dimensional random fields is Principal Component Analysis (PCA), which is based on the Karhunen-Lo eve expansion (KLE) of a stochastic process with finite variance. The KLE is analogous to a Fourier series expansion for a random process, where the goal is to find an orthogonal transformation for the data such that the projection of the data onto this orthogonal subspace is optimal in the L 2 sense, i.e, which minimizes the mean square error. In practice, this orthogonal transformation is determined by performing an SVD (Singular Value Decomposition) on the sample covariance matrix or on the data matrix itself. Sampling error is typically ignored when quantifying the principal components, or, equivalently, basis functions of the KLE. Furthermore, it is exacerbated when the sample size is much smaller than the dimension of the random field. In this paper, we introduce a Bayesian KLE procedure, allowing one to obtain a probabilistic model on the principal components, which can account for inaccuracies due to limited sample size. The probabilistic model is built via Bayesian inference, from which the posterior becomes the matrix Bingham density over the space of orthonormal matrices. We use a modified Gibbs sampling procedure to sample on this space and then build a probabilistic Karhunen-Lo eve expansions over random subspaces to obtain a set of low-dimensional surrogates of the stochastic process. We illustrate this probabilistic procedure with a finite dimensional stochastic process inspired by Brownian motion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Polyethylene Grafted Silica Nanoparticles Prepared via Surface-Initiated ROMP

Polyethylene and nanosilica represent the most ubiquitous commodity plastic and nanocomposite filler, respectively. Despite their possible utility, few examples exist in the literature of successfully combining these two materials to form polyethylene nanocomposites. Synthesizing well-defined polyethylene grafted to a surface is a significant challenge in the nanocomposites community. Presented in this paper is a synthetic approach toward polyethylene grafted nanoparticles with controllable graft density and molecular weight of the grafted polymer. The variably grafted nanoparticles were then incorporated into a commercial high density polyethylene matrix. The synthesis, characterization, and challenges in making these materials are further discussed.

36 MATERIALS SCIENCE↗