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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Finding Matrix Product State Representations of Highly Excited Eigenstates of Many-Body Localized Hamiltonians

A key property of many-body localized Hamiltonians is the area law entanglement of even highly excited eigenstates. Matrix product states (MPS) can be used to efficiently represent low entanglement (area law) wave functions in one dimension. An important application of MPS is the widely used density matrix renormalization group (DMRG) algorithm for finding ground states of one-dimensional Hamiltonians. Here, we develop two algorithms, the shift-and-invert MPS (SIMPS) and excited state DMRG which find highly excited eigenstates of many-body localized Hamiltonians. Excited state DMRG uses a modified sweeping procedure to identify eigenstates, whereas SIMPS applies the inverse of the shifted Hamiltonian to a MPS multiple times to project out the targeted eigenstate. To demonstrate the power of these methods, we verify the breakdown of the eigenstate thermalization hypothesis in the many-body localized phase of the random field Heisenberg model, show the saturation of entanglement in the many-body localized phase, and generate local excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Spin-Orbital Density Wave and a Mott Insulator in a Two-Orbital Hubbard Model on a Honeycomb Lattice

Inspired by the recent discovery of correlated insulating states in twisted bilayer graphene, we study a two-orbital Hubbard model on the honeycomb lattice with two electrons per unit cell. Based on the real-space density matrix renormalization group simulation, we identify a metal-insulator transition around U c /t = 2.5–3. In the vicinity of U c , we find strong spin-orbital density wave fluctuations at commensurate wave vectors, accompanied by weaker incommensurate charge density wave fluctuations. Here, the spin-orbital density wave fluctuations are enhanced with increasing system sizes, suggesting the possible emergence of long-range order in the two-dimensional limit. At larger U, our calculations indicate a possible nonmagnetic Mott insulator phase without spin or orbital polarization. Our findings offer new insight into correlated electron phenomena in twisted bilayer graphene and other multiorbital honeycomb materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Superconductivity in the doped Hubbard model and its interplay with next-nearest hopping t ′

The Hubbard model is widely believed to contain the essential ingredients of high-temperature superconductivity. However, proving definitively that the model supports superconductivity is challenging. Here, we report a large-scale density matrix renormalization group study of the lightly doped Hubbard model on four-leg cylinders at hole doping concentration δ = 12.5%. We reveal a delicate interplay between superconductivity and charge density wave and spin density wave orders tunable via next-nearest neighbor hopping t ′. For finite t ′, the ground state is consistent with a Luther-Emery liquid with power-law superconducting and charge density wave correlations associated with half-filled charge stripes. In contrast, for t ′ = 0, superconducting correlations fall off exponentially, whereas charge density and spin density modulations are dominant. Our results indicate that a route to robust long-range superconductivity involves destabilizing insulating charge stripes in the doped Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Tunable plasmon-enhanced second-order optical nonlinearity in transition metal dichalcogenide nanotriangles

The development of nanomaterials with a large nonlinear susceptibility is essential for nonlinear nanophotonics. Here we show that transition metal dichalcogenide (TMD) nanotriangles have a large effective second-order susceptibility [χ (2) ] at midinfrared to near-infrared frequencies owing to their broken centrosymmetry. χ (2) is calculated within the density-matrix formalism that accounts for dissipation and screening. χ (2) peaks in the vicinity of both two-photon resonances (specified by the geometry) and plasmon resonances (tunable via the carrier density). Aligning the resonances yields the values of χ (2) as high as 10 –6 m/V. These findings underscore the potential of TMD nanotriangles for nonlinear nanophotonics, particularly second-harmonic generation.

0-dimensional systems↗
DOE OSTI2022

Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

Here, we develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems, including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input–output and the hierarchical equations of motion formalisms into a double hierarchy of density matrix equations. We show analytically that under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an N-photon Fock state input induces no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon with or without the influence of phonons. For short pulses (or, equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δ and an exciton–light coupling constant Γ. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyze the absorption probability in the opposite long pulse (narrow bandwidth) regime. We then derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective vs independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy equations.

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

One-dimensional Holstein model revisited

Here, we analyze the global ground-state (quantum) phase diagram of the one-dimensional spinful Holstein model at half filling as a function of the strength of the electron-phonon coupling (represented by the strength of the phonon-induced attraction, U) and the phonon frequency ω 0 . In addition to reanalyzing the various asymptotic regimes, we carry out density-matrix renormalization group simulations to correct previous inferences concerning the antiadiabatic (large ω 0 ) and strong-coupling (large U) regimes. There are two distinct phases—a fully gapped commensurate charge-density-wave phase and a spin-gapped Luther-Emery phase with a gapless charge mode—separated by a phase boundary, with a shape that reflects different microscopic physics in the weak- and strong-coupling limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Robust d-Wave Superconductivity in the Square-Lattice t–J Model

Unravelling competing orders emergent in doped Mott insulators and their interplay with unconventional superconductivity is one of the major challenges in condensed matter physics. Here, to explore the possible superconducting state in a doped Mott insulator, we study the square-lattice t-J model with both the nearest-neighbor and next-nearest-neighbor electron hoppings and spin interactions. By using the state-of-the-art density matrix renormalization group calculation with imposing charge U(1) and spin SU(2) symmetries on the six-leg cylinders, we establish a quantum phase diagram including three phases: a stripe charge density wave phase, a superconducting phase without static charge order, and a superconducting phase coexistent with a weak charge stripe order. Crucially, we demonstrate that the superconducting phase has a power-law pairing correlation that decays much slower than the charge density and spin correlations, which is a quasi-1D descendant of the uniform d-wave superconductor in two dimensions. These findings reveal that enhanced charge and spin fluctuations with optimal doping is able to produce robust d-wave superconductivity in doped Mott insulators, providing a foundation for connecting theories of superconductivity to models of strongly correlated systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Final Technical Report for "Cloud-based Low-Scaling Quantum Chemistry Simulations for Materials"

The objective of phase I was to develop the basic infrastructure for performing efficient hybrid DFT calculations on solid materials with a turnaround time amenable to the use in large-scale data-driven approaches. To fulfill the goal we have developed a novel algorithm that significantly reduces the cost of exchange matrix formation. The algorithm does so by making use of three ingredients (a) interpolative decomposition of the electron integrals (b) robust pseudospectral method and (c) occ-RI approach. Our published work demonstrates that the algorithm is orders of magnitude faster than any other hybrid-DFT method and is on the order of only 3-4 times slower than pure DFT. All steps of the algorithm developed during phase I can be accelerated to the point that, for large enough systems, the diagonalization of the Fock matrix will become the most expensive step. To treat such large systems we have created an interface to the ASCR-funded PEXSI library. The PEXSI method is used to compute the density matrix directly from the Fock matrix in a manner that preserves the sparsity of the local representation.

Shiozaki, Toru↗
DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Fractional charge and emergent mass hierarchy in diagonal two-leg t – J cylinders

Here, we define a class of “diagonal” t $-$ J ladders rotated by π / 4 relative to the canonical lattice directions of the square lattice, and study it using density matrix renormalization group. Here, we focus on the two-leg cylinder with a doped hole concentration near x = $\frac{1}{4}$ . At exactly x = $\frac{1}{4}$, the system forms a period 4 charge density wave and exhibits spin-charge separation. Slightly away from $\frac{1}{4}$ doping, we observe several topologically distinct types of solitons with well-defined fractionalized quantum numbers. Remarkably, given the absence of any obvious small parameter, the effective masses of the emergent solitons differ by several orders of magnitude.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1995

A far-wing line shape theory which satisfies the detailed balance principle

A far-wing theory in which the validity of the detailed balance principle is maintained in each step of the derivation is presented. The role of the total density matrix including the initial correlations is analyzed rigorously. By factoring out the rapidly varying terms in the complex-time development operator in the interaction representation, better approximate expressions can be obtained. As a result, the spectral density can be expressed in terms of the line-coupling functions in which two coupled lines are arranged symmetrically and whose frequency detunings are omega - 1/2(omega(sub ji) + omega (sub j'i'). Using the approximate values omega - omega(sub ji) results in expressions that do not satisfy the detailed balance principle. However, this principle remains satisfied for the symmetrized spectral density in which not only the coupled lines are arranged symmetrically, but also the initial and final states belonging to the same lines are arranged symmetrically as well.

Ma, Q.↗
DOE OSTI2024

Phase diagram of the one-dimensional t 1 - t 2 - J model: Ferromagnetism, triplet pairing, and charge and pair density waves

Here, we present a density matrix renormalization group (DMRG) study of an extended t-J model with hopping to the first and second neighbors -- the one dimensional t 1 -t 2 -J model. The full phase diagram as a function of the density n and exchange strength J, for both positive and negative values of t 2 , is obtained. For t 2 = -0.5 we observe that, in the strongly interacting region, Nagaoka Ferromagnetism (FM) is accompanied by a triplet pair density wave (PDW) upon doping. As the spin exchange J increases, a charge density wave (CDW) phase emerges and then gives way to singlet superconductivity (SC). This phase behaves as a singlet PDW with vanishing spacial average of the order parameter and a spin gap. When t 2 = 0.5, the physics is basically reminiscent of the conventional t-J model, undergoing a transition from a metallic to a SC phase as a function of J.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nonequilibrium spintronic transport through Kondo impurities

In this work we analyze the nonequilibrium transport through a quantum impurity (quantum dot or molecule) attached to ferromagnetic leads by using a hybrid numerical renormalization group–time-dependent density matrix renormalization group thermofield quench approach. For this, we study the bias dependence of the differential conductance through the system, which shows a finite zero-bias peak, characteristic of the Kondo resonance and reminiscent of the equilibrium local density of states. In the nonequilibrium settings, the resonance in the differential conductance is also found to decrease with increasing the lead spin polarization. The latter induces an effective exchange field that lifts the spin degeneracy of the dot level. Therefore, as we demonstrate, the Kondo resonance can be restored by counteracting the exchange field with a finite external magnetic field applied to the system. Lastly, we investigate the influence of temperature on the nonequilibrium conductance, focusing on the split Kondo resonance. Our work thus provides an accurate quantitative description of the spin-resolved transport properties relevant for quantum dots and molecules embedded in magnetic tunnel junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Interaction-driven quantum anomalous Hall insulator in a Dirac semimetal

Here, the interaction-driven quantum anomalous Hall (QAH) insulator has been sought for a long time in a Dirac semimetal with linear band touching points at the Fermi level. By combining exact diagonalization, density matrix renormalization group, and analytical methods, we study a spinless fermion system on the checkerboard lattice with twofold rotational symmetry, which realizes two Dirac band touching points in the absence of interaction. At weak coupling, the Dirac semimetal is stable. At a finite density-density repulsive interaction, we analyze possible symmetry broken states, and find that an QAH state is stabilized when the interaction strength exceeds the energy scale controlling the separation between the Dirac points. Through numerical simulations, we verify the existence of the QAH phase with spontaneous time-reversal symmetry breaking and quantized Chern number C = 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Exact matrix product state representation and convergence of a fully correlated electronic wavefunction in the infinite-basis limit

Here In this paper we present the exact representation of a fully correlated electronic wavefunction as the single-particle basis approaches completeness. It consists of a half-infinite chain of matrices of exponentially increasing size. The complete basis limit is illustrated numerically using the density-matrix renormalization-group method by computing the core-valence entanglement in the C 2 ground state in increasing subsets of cc-pVTZ and pVQZ bases until convergence is reached.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum Phases of Transition Metal Dichalcogenide Moiré Systems

Moiré systems provide a rich platform for studies of strong correlation physics. Recent experiments on heterobilayer transition metal dichalcogenide Moiré systems are exciting in that they manifest a relatively simple model system of an extended Hubbard model on a triangular lattice. Inspired by the prospect of the hetero-transition metal dichalcogenide Moiré system’s potential as a solid-state-based quantum simulator, we explore the extended Hubbard model on the triangular lattice using the density matrix renormalization group. Specifically, we explore the two-dimensional phase space spanned by the key tuning parameters in the extended Hubbard model, namely, the kinetic energy strength and the further-range Coulomb interaction strengths. We find competition between Fermi fluid, chiral spin liquid, spin density wave, and charge order. In particular, our finding of the optimal further-range interaction for the chiral correlation presents a tantalizing possibility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Interplay of Collective Excitations in Quantum Well Intersubband Resonances

Intersubband resonances in a semiconductor quantum well (QW) display some of the most fascinating features involving various collective excitations such as Fermi-edge singularity (FES) and intersubband plasmon (ISP). Using a density matrix approach, we treated many-body effects such as depolarization, vertex correction, and self-energy consistently for a two-subband system. We found a systematic change in resonance spectra from FES-dominated to ISP-dominated features, as QW- width or electron density is varied. Such an interplay between FES and ISP significantly changes both line shape and peak position of the absorption spectrum. In particular, we found that a cancellation of FES and ISP undresses the resonant responses and recovers the single-particle features of absorption for semiconductors with a strong nonparabolicity such as InAs, leading to a dramatic broadening of the absorption spectrum.

Li, Jian-Zhong↗
DOE OSTI2021

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗