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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Algebraic Hastatic Order in One-Dimensional Two-Channel Kondo Lattice

The two-channel Kondo lattice likely hosts a rich array of phases, including hastatic order, a channel symmetry breaking heavy Fermi liquid. In this work, we revisit its one-dimensional phase diagram using density matrix renormalization group and, in contrast to previous work, find algebraic hastatic orders generically for stronger couplings. These are heavy Tomonaga-Luttinger liquids with nonanalyticities at Fermi vectors captured by hastatic density waves. We also find a predicted additional nonlocal order parameter due to interference between hastatic spinors, not present at large N, and residual repulsive interactions at strong coupling suggesting non-Fermi-liquid physics in higher dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Analytic first derivatives for a spin-adapted open-shell coupled cluster theory: Evaluation of first-order electrical properties

An analytic scheme is presented for the evaluation of first derivatives of the energy for a unitary group based spin-adapted coupled cluster (CC) theory, namely, the combinatoric open-shell CC (COSCC) approach within the singles and doubles approximation. The widely used Lagrange multiplier approach is employed for the derivation of an analytical expression for the first derivative of the energy, which in combination with the well-established density-matrix formulation, is used for the computation of first-order electrical properties. Derivations of the spin-adapted lambda equations for determining the Lagrange multipliers and the expressions for the spin-free effective density matrices for the COSCC approach are presented. Orbital-relaxation effects due to the electric-field perturbation are treated via the Z-vector technique. We present calculations of the dipole moments for a number of doublet radicals in their ground states using restricted open-shell Hartree-Fock (ROHF) and quasi-restricted HF (QRHF) orbitals in order to demonstrate the applicability of our analytic scheme for computing energy derivatives. We also report calculations of the chlorine electric-field gradients and nuclear quadrupole-coupling constants for the CCl, CH{sub 2}Cl, ClO{sub 2}, and SiCl radicals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Coherence Evolution and Transfer Supplemented by Sender’s Initial-State Restoring

The evolution of quantum coherences comes with a set of conservation laws provided that the Hamiltonian governing this evolution conserves the spin-excitation number. At that, coherences do not intertwist during the evolution. Using the transmission line and the receiver in the initial ground state we can transfer the coherences to the receiver without interaction between them, although the matrix elements contributing to each particular coherence intertwist in the receiver’s state. Therefore we propose a tool based on the unitary transformation at the receiver side to untwist these elements and thus restore (at least partially) the structure of the sender’s initial density matrix. A communication line with two-qubit sender and receiver is considered as an example of implementation of this technique.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Quantum chi-squared and goodness of fit testing

A quantum mechanical hypothesis test is presented for the hypothesis that a certain setup produces a given quantum state. Although the classical and the quantum problems are very much related to each other, the quantum problem is much richer due to the additional optimization over the measurement basis. A goodness of fit test for i.i.d quantum states is developed and a max-min characterization for the optimal measurement is introduced. We find the quantum measurement which leads both to the maximal Pitman and Bahadur efficiencies, and determine the associated divergence rates. We discuss the relationship of the quantum goodness of fit test to the problem of estimating multiple parameters from a density matrix. These problems are found to be closely related and we show that the largest error of an optimal strategy, determined by the smallest eigenvalue of the Fisher information matrix, is given by the divergence rate of the goodness of fit test.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Singlet, triplet, and pair density wave superconductivity in the doped triangular-lattice moiré system

Recent experimental progress has established the twisted bilayer transition metal dichalcogenide (TMD) as a highly tunable platform for studying many-body physics. Particularly, the homobilayer TMDs under displacement field are believed to be described by a generalized triangular-lattice Hubbard model with a spin-dependent hopping phase θ. To explore the effects of θ on the system, we perform density matrix renormalization group calculations for the relevant triangular lattice t-J model. By changing θ at small hole doping, we obtain a region of quasi-long-range superconducting order coexisting with charge and spin density wave within 0 < θ < π/3. The superconductivity is composed of a dominant spin singlet d-wave and a subdominant triplet p-wave pairing. Intriguingly, the S z = ±1 triplet pairing components feature pair-density waves. In addition, we find a region of triplet superconductivity coexisting with charge-density wave and ferromagnetism within π/3 < θ < 2π/3, which is related to the former phase at smaller θ by a combined operation of spin-flip and gauge transformation. Furthermore, our findings provide insights and directions for experimental search for exotic superconductivity in twisted TMD systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Using Graph Partitioning for Scalable Distributed Quantum Molecular Dynamics

The simulation of the physical movement of multi-body systems at an atomistic level, with forces calculated from a quantum mechanical description of the electrons, motivates a graph partitioning problem studied in this article. Several advanced algorithms relying on evaluations of matrix polynomials have been published in the literature for such simulations. We aim to use a special type of graph partitioning to efficiently parallelize these computations. For this, we create a graph representing the zero–nonzero structure of a thresholded density matrix, and partition that graph into several components. Each separate submatrix (corresponding to each subgraph) is then substituted into the matrix polynomial, and the result for the full matrix polynomial is reassembled at the end from the individual polynomials. This paper starts by introducing a rigorous definition as well as a mathematical justification of this partitioning problem. We assess the performance of several methods to compute graph partitions with respect to both the quality of the partitioning and their runtime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Strongly modified four-wave mixing in a coupled semiconductor quantum dot-metal nanoparticle system

We study the four-wave mixing effect in a coupled semiconductor quantum dot-spherical metal nanoparticle structure. Depending on the values of the pump field intensity and frequency, we find that there is a critical distance that changes the form of the spectrum. Above this distance, the four-wave mixing spectrum shows an ordinary three-peaked form and the effect of controlling its magnitude by changing the interparticle distance can be obtained. Below this critical distance, the four-wave mixing spectrum becomes single-peaked; and as the interparticle distance decreases, the spectrum is strongly suppressed. The behavior of the system is explained using the effective Rabi frequency that creates plasmonic metaresonances in the hybrid structure. In addition, the behavior of the effective Rabi frequency is explained via an analytical solution of the density matrix equations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Noether’s theorem for dissipative quantum dynamical semi-groups

Noether’s theorem on constants of the motion of dynamical systems has recently been extended to classical dissipative systems (Markovian semi-groups) by Baez and Fong [J. Math. Phys. 54, 013301 (2013)]. We show how to extend these results to the fully quantum setting of quantum Markov dynamics. For finite-dimensional Hilbert spaces, we construct a mapping from observables to completely positive maps that leads to the natural analogue of their criterion of commutativity with the infinitesimal generator of the Markov dynamics. Using standard results on the relaxation of states to equilibrium under quantum dynamical semi-groups, we are able to characterise the constants of the motion under quantum Markov evolutions in the infinite-dimensional setting under the usual assumption of existence of a stationary strictly positive density matrix. In particular, the Noether constants are identified with the fixed point of the Heisenberg picture semi-group.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

π phase shift across stripes in a charge density wave system

Many strongly correlated materials are characterized by deeply intertwined charge and spin order. Besides their high superconducting transition temperatures, one of the central features of these complex patterns in cuprates is a phase shift which occurs across lines of decreased hole density. That is, when doped away from their antiferromagnetic phase, the additional charge is not distributed uniformly, but rather in “stripes.” The sublattices preferentially occupied by up and down spin are reversed across these stripes, a phenomenon referred to as a “π phase shift.” Many of the spin-charge patterns, including the π phase shift, are reproduced by density matrix renormalization group and quantum Monte Carlo calculations of simplified tight binding (repulsive Hubbard) models. In this paper we demonstrate that this sublattice reversal is generic by considering the corresponding phenomenon in the attractive Hubbard Hamiltonian, where a charge density wave phase forms at half filling. Here, we introduce charge stripes via an appropriate local chemical potential; measurements of charge correlation across the resulting lines of lowered density reveal a clear π phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Design and performance testing of a novel far field Gamma system to assay radioactive waste in 400 and 200 litre drums - 15339

A Far Field Gamma ray Measurement system has been developed for measuring and assaying both 400 and 200 litre waste drums with a wide range of activity at nuclear power stations. These include radioactive waste drums with surface dose-rates of less than 2 mSv/h (200 mrem/h) and, at the higher range of activity, 400 litre drums with surface dose-rates of up to 50 mSv/h (5 rem/h). The assay system incorporates an electro-mechanically cooled high purity Germanium profile detector and advanced digital spectroscopic analysis electronics. A novel optimised conical-trapezoidal collimator has been incorporated to reduce the effect of the radioactive background in addition to the conventional detector lead shield with tin and copper graded lining to reduce the effect of lead X-rays. Two features have been incorporated to reduce detector dead time and maintain detector resolution when measuring high activity drums. The first is the inclusion of two automated tungsten filters of different thickness as an integral part of the detector collimator. The filters are used in conjunction with an automated rail system on which the detector platform is mounted. The detector can be positioned at different distances from the drum surface in order to reduce detector dead time. The deployment of the filters and the adjustment of the drum - detector position are automatically controlled based on a user defined dose-rate limit table. The data is provided by two Geiger-Muller dose-rate detectors, which are used to measure the drum surface dose-rate. This data is subsequently used to adjust the detector position and in cases of higher dose-rate deploy one or other of the tungsten filters. The drum rotation platform incorporates a load cell to determine drum weight and hence the drum density. This information is used to calculate an attenuation correction based on the assumption of uniform drum density and uniform distribution of activity within the drum. Where the waste drums have a more complicated and known regular internal structure, such as a small radioactive region surrounded by an annulus of shielding material, the analysis algorithms incorporated into the far field geometry and attenuation correction spectroscopic data analysis code are able to make a more accurate determination of drum activity than would be possible assuming a completely uniform drum. The detector energy and efficiency calibration is achieved using a point source with multiple gamma ray energy peaks. The response of the high purity Germanium detector is also modelled using the Monte Carlo Neutron Photon code and the model is benchmarked using the measured point source calibration data. The performance of the instrument has been determined for a range of uniform density matrix drums containing a set of Eu-152 line or rod sources located in re-entrant tubes positioned on an equal volume basis. When rotated, the test drums with volume distributed line sources simulate waste drums of different density each with a uniformly distributed radioactive source. The measurements have been validated using Monte Carlo simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Doped Mott insulators in the triangular-lattice Hubbard model

Here we investigate the evolution of the Mott insulators in the triangular lattice Hubbard Model, as a function of hole doping δ in both the strong and intermediate coupling limits. Using the advanced density matrix renormalization group (DMRG) method, at light hole doping δ ≲ 10%, we find a significant difference between strong and intermediate couplings. Notably, at intermediate coupling an unusual metallic state emerges, with short ranged spin correlations but long ranged spin-chirality order. Moreover, no clear Fermi surface or wave vector is observed, this chiral metal also exhibits staggered loop current, which breaks the translational symmetry. These features disappear on increasing interaction strength or on further doping. At strong coupling, the 120 degree magnetic order of the insulating magnet persists for light doping, and produces hole pockets with a well-defined Fermi surface. On further doping, δ ≈ 10%~20% SDW order and coherent hole Fermi pockets are found at both strong and intermediate couplings. At even higher doping δ ≳ 20%, the SDW order is suppressed and the spin-singlet Cooper pair correlations are simultaneously enhanced. We also briefly comment on the strong particle-hole asymmetry of the model.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Renormalization of quantum discord and Bell nonlocality in the XXZ model with Dzyaloshinskii–Moriya interaction

In this paper, we have investigated the dynamical behaviors of the two important quantum correlation witnesses, i.e. geometric quantum discord (GQD) and Bell–CHSH inequality in the XXZ model with DM interaction by employing the quantum renormalization group (QRG) method. The results have shown that the anisotropy suppresses the quantum correlations while the DM interaction can enhance them. Meanwhile, using the QRG method we have studied the quantum phase transition of GQD and obtained two saturated values, which are associated with two different phases: spin-fluid phase and the Néel phase. It is worth mentioning that the block–block correlation is not strong enough to violate the Bell–CHSH inequality in the whole iteration steps. Moreover, the nonanalytic phenomenon and scaling behavior of Bell inequality are discussed in detail. As a byproduct, the conjecture that the exact lower and upper bounds of Bell inequality versus GQD can always be established for this spin system although the given density matrix is a general X state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Nonlinear properties of gated graphene in a strong electromagnetic field

We develop a microscopic theory of a strong electromagnetic field interaction with gated bilayer graphene. Quantum kinetic equations for density matrix are obtained using a tight binding approach within second quantized Hamiltonian in an intense laser field. We show that adiabatically changing the gate potentials with time may produce (at resonant photon energy) a full inversion of the electron population with high density between valence and conduction bands. In the linear regime, excitonic absorption of an electromagnetic radiation in a graphene monolayer with opened energy gap is also studied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Enhanced superconductivity by near-neighbor attraction in the doped extended Hubbard model

A recent experiment has unveiled an anomalously strong electron-electron attraction in the one-dimensional copper-oxide chain Ba 2–x Sr x CuO 3+δ . While the effect of the near-neighbor electron attraction V in the one-dimensional extended Hubbard chain has been examined recently, its effect in the Hubbard model beyond the one-dimensional chain remains unclear. Here, in this work, we report a density matrix renormalization group study of the extended Hubbard model on long four-leg cylinders on the square lattice. We find that the near-neighbor electron attraction V can notably enhance the long-distance superconducting correlations while simultaneously suppressing the charge-density-wave (CDW) correlations. Specifically, for a modestly strong electron attraction, the superconducting correlations become dominant over the CDW correlations with a Luttinger exponent K sc ~1 and strong divergent superconducting susceptibility. Our results provide a promising way to realize long-range superconductivity in the doped Hubbard model. The relevance of our numerical results to cuprate materials is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Study of Thermal and Coherent A{sub 1g} Phonons in Bismuth Telluride

Totally symmetric A{sub 1g} phonons are studied for the equilibrium and coherent states of a Bi{sub 2}Te{sub 3} lattice. Equilibrium phonons were investigated in the frequency domain by the method of spontaneous Raman scattering, whereas coherent phonons were studied by the method of active femtosecond spectroscopy in the time domain. In the latter case, femtosecond laser pulses were used both for generating and detecting coherent A{sub 1g} phonons having a well-defined phase allowing the selective optical control of the lattice dynamics. A comparison of the results obtained in the frequency and time domains suggests that diagonal and nondiagonal elements of the density matrix of lattice excitations relax with the same characteristic time to the equilibrium and zero values, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Exact stochastic unraveling of an optical coherence dynamics by cumulant expansion

A numerically exact Monte Carlo scheme for calculation of open quantum system dynamics is proposed and implemented. The method consists of a Monte Carlo summation of a perturbation expansion in terms of trajectories in Liouville phase-space with respect to the coupling between the excited states of the molecule. The trajectories are weighted by a complex decoherence factor based on the second-order cumulant expansion of the environmental evolution. The method can be used with an arbitrary environment characterized by a general correlation function and arbitrary coupling strength. It is formally exact for harmonic environments, and it can be used with arbitrary temperature. Time evolution of an optically excited Frenkel exciton dimer representing a molecular exciton interacting with a charge transfer state is calculated by the proposed method. We calculate the evolution of the optical coherence elements of the density matrix and linear absorption spectrum, and compare them with the predictions of standard simulation methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗