Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Excited states in hydrogenated single-layer MoS 2

Abstract Our calculations of the excitation spectrum of single-layer MoS 2 at several hydrogen coverages, using a density-matrix based time-dependent density-functional theory (TDDFT) show that the fully hydrogenated system is metallic, while at lower coverages the spectrum consists of spin-polarized partially filled localized mid-gap states. The calculated absorption spectrum of the system reveals standard excitonic peaks corresponding to the bound valence-band hole and conduction-band electron, as well as excitonic peaks that involve the mid-gap states. Binding energies of the excitons of the hydrogenated system are found to be relatively large (few tens of meV), making their experimental detection facile and suggesting hydrogenation as a knob for tuning the optical properties of single-layer MoS 2 . Importantly, we find hydrogenation to suppress visible light photoluminescence, in agreement with experimental observations. In contrast, both Li and Na atoms transform the system into an n-doped non-magnetic semiconductor that does not allow excitonic states.

Physics↗
DOE OSTI2021

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗
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 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 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 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 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 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 OSTI2022

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗
DOE OSTI2025

Superconductivity in the lightly doped Hubbard model on the cylindrical honeycomb lattice

Here, we have performed large-scale density-matrix renormalization group studies of the lightly doped Hubbard model on the honeycomb lattice on long three- and four-leg cylinders. We find that the ground state of the system upon lightly doping is consistent with that of a superconducting state with coexisting quasi-long-range superconducting and charge density wave orders. Both the superconducting and charge density wave correlations decay as a power law at long distances with corresponding exponents 𝐾 𝑠⁢𝑐 < 2 and 𝐾 𝑐 < 2. On the contrary, the spin-spin and single-particle correlations decay exponentially, although with relatively long correlation lengths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Majorana zero modes in Y -shape interacting Kitaev wires

Motivated by the recent experimental realization of minimal Kitaev chains using quantum dots, we investigate the Majorana zero modes (MZM) in Y-shape Kitaev wires. We solve the associated Kitaev models analytically at the sweet spot (t h = Δ) and derive the exact form of MZM wave-functions in this geometry. We observed exotic multi-site MZMs located near the junction center, on the nearby edge sites of each leg. This result is important for potential braiding of Majoranas and the performance of Y-junctions made from arrays of quantum dots. Furthermore, we study the stability of local (single-site) and multi-site MZMs modes in the presence of Coulomb repulsion, using density matrix renormalization group theory. Our local density-of-states calculation shows that these multi-site MZMs are as equally topologically protected as the single-site MZMs when in the presence of Coulomb repulsion or when away from the sweet-spot.

36 MATERIALS SCIENCE↗
DOE OSTI2020

TeraChem: Accelerating electronic structure and ab initio molecular dynamics with graphical processing units

Developed over the past decade, TeraChem is an electronic structure and ab initio molecular dynamics software package designed from the ground up to leverage graphics processing units (GPUs) to perform large-scale ground and excited state quantum chemistry calculations in the gas and the condensed phase. TeraChem’s speed stems from the reformulation of conventional electronic structure theories in terms of a set of individually optimized high-performance electronic structure operations (e.g., Coulomb and exchange matrix builds, one- and two-particle density matrix builds) and rank-reduction techniques (e.g., tensor hypercontraction). Recent efforts have encapsulated these core operations and provided language-agnostic interfaces. Finally, this greatly increases the accessibility and flexibility of TeraChem as a platform to develop new electronic structure methods on GPUs and provides clear optimization targets for emerging parallel computing architectures.

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

BCS-BEC crossover in a (t 2g ) 4 excitonic magnet

The condensation of spin-orbit-induced excitons in t 2 g 4 electronic systems is attracting considerable attention. At large Hubbard U , antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ( J eff = 1 ). Furthermore, we show that even at intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ( λ ) in one dimension using the density matrix renormalization group, while in two dimensions we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full λ versus U phase diagrams for both one- and two-dimensional lattices. Our main result is that at intermediate Hubbard U , increasing λ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O ( a ) and O ( 10 a ) in 1 d and 2 d , respectively, with a being the lattice spacing. Further increasing λ , the system eventually crosses over to the BEC limit (with r coh ≪ a ).

1-dimensional systems↗
DOE OSTI2021

Origin of the magnetic and orbital ordering in α-Sr 2 CrO 4

Motivated by recent experimental progress in transition metal oxides with the K 2 NiF 4 structure, in this study we investigate the magnetic and orbital ordering in α-Sr 2 CrO 4 . Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the ab initio band structure near the Fermi level. The unique reverse splitting of t 2g orbitals in α-Sr 2 CrO 4 , with the 3d 2 electronic configuration for the Cr 4+ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground state phase diagram for the two dimensional compound α-Sr 2 CrO 4 . We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of α-Sr 2 CrO 4 , we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound α-Sr 2 CrO 4 , in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with d xy driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the d xz and d yz orbitals driving the orbital order.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Triangular lattice Hubbard model physics at intermediate temperatures

Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Time-dependent approach to inelastic scattering spectroscopies in and away from equilibrium: Beyond perturbation theory

We propose a nonperturbative numerical approach to calculate the spectrum of a many-body Hamiltonian with time and momentum resolution by exactly recreating a scattering event using the time-dependent Schrödinger equation. Akin to an actual inelastic scattering experiment, we explicitly account for the incident and scattered particles (e.g., photons, neutrons, electrons, etc.) in the Hamiltonian and obtain the spectrum by measuring the energy and momentum lost by the particle after interacting with the sample. We illustrate the method by calculating the spin excitations of a Mott-insulating Hubbard chain after a sudden quench with the aid of the time-dependent density matrix renormalization group method. Our formalism can be applied to different forms of spectroscopies, such as neutron and Compton scattering and electron-energy-loss spectroscopy, as an example.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

High-pressure inelastic neutron scattering study of the anisotropic S = 1 spin chain [ Ni ( HF 2 ) ( 3 - Clpyradine ) 4 ] B F 4

[Ni(HF 2 )(3-Clpyradine) 4 ]BF 4 (NBCT) is a one-dimensional, S=1 spin-chain material that shows no long-range magnetic order down to thermometer temperatures of 0.1 K. Previous ambient pressure inelastic neutron scattering experiments identified NBCT to be in the large-D quantum paramagnetic phase of the D/J phase diagram, where D is the axial single-ion anisotropy and J is the intrachain superexchange. In this work, we extend the previous experiments to a hydrostatic pressure of 0.9 GPa. By comparing to density matrix renormalization group calculations, we find D/J increases from 1.5 to 3.2 as pressure increases from 0 GPa to 0.9 GPa, which pushes the system further into the large-D phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗