Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix methods”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Where is the Quantum Spin Nematic?

We provide strong evidence of the spin-nematic state in a paradigmatic ferro-antiferromagnetic J 1 –J 2 model using analytical and density-matrix renormalization group methods. In zero field, the attraction of spin-flip pairs leads to a first-order transition and no nematic state, while pair repulsion at larger J 2 stabilizes the nematic phase in a narrow region near the pair-condensation field. A devil’s staircase of multipair condensates is conjectured for weak pair attraction. A suppression of the spin-flip gap by many-body effects leads to an order-of-magnitude contraction of the nematic phase compared to naïve expectations. Here, the proposed phase diagram should be broadly valid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl 2 and OsOCl 2

The van der Waals oxide dichlorides MOX 2 (M=V, Ta, Nb, Ru, and Os; X=halogen element), with different electronic densities, are attracting considerable attention. Ferroelectricity, spin-singlet formation, and orbital-selective Peierls phases were reported in this family with d 1 or d 2 electronic configurations, all believed to be caused by the strongly anisotropic electronic orbital degree of freedom. Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl 2 and OsOCl 2 with d 4 electronic configurations. Different from a previous study using VOI 2 with d 1 configuration, these systems with 4d4 or 5d4 do not exhibit a ferroelectric instability along the a axis. Due to the fully occupied dxy orbital in RuOCl 2 and OsOCl 2 , the Peierls instability distortion disappears along the b axis, leading to an undistorted Immm phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the a axis. For this reason, the materials of our focus can be regarded as “effective one-dimensional” systems even when they apparently have a dominant two-dimensional lattice geometry. The large crystal-field splitting energy (between dxz/yz and dxy orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the J = 0 singlet state in MOCl 2 (M=Ru or Os) with electronic density n = 4, resulting in a spin-1 system. Moreover, we find staggered antiferromagnetic order with π wave vector along the M-O chain direction (a axis) while the magnetic coupling along the b axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the t 2g orbitals for the Os atoms to construct a multiorbital Hubbard model for the M-O chains. Staggered AFM with ↑|-|↓-↑|-|↓ spin structure dominates in our density matrix renormalization group calculations, in agreement with density functional theory calculations. Our results for RuOCl 2 and OsOCl 2 provide guidance to experimentalists and theorists working on this interesting family of oxide dichlorides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spin-1 Kitaev-Heisenberg model on a honeycomb lattice

We study the Kitaev-Heisenberg model with spin-1 local degree of freedom on a honeycomb lattice numerically by the infinite density matrix renormalization group method on a cylinder geometry. By tuning the relative value of the Kitaev and Heisenberg exchange couplings, we obtain the phase diagram with two spin liquid phases and four symmetry-broken phases. Here, we identify that the spin liquid phases are gapless by calculating the central charge at the pure Kitaev points without the Heisenberg interactions. Comparing to its spin-1/2 counterpart, the position and number of gapless modes of the spin-1 case are quite different. Due to the approximate Z 2 local conservations, the expectation value of the Wilson loop operator measuring the flux of each plaquette stays near to 1, and the static spin-spin correlations remain short range in the entire spin liquid phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Block orbital-selective Mott insulators: A spin excitation analysis

We present a comprehensive study of the spin excitations—as measured by the dynamical spin structure factor S(q,ω)—of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system develops periodic ferromagnetic islands of various shapes and sizes, which are antiferromagnetically coupled. The 2×2 particular case was already found experimentally in the ladder material BaFe2Se3 that becomes superconducting under pressure. Here we discuss the electronic density as well as Hubbard and Hund coupling dependence of S(q,ω) using density matrix renormalization group method. Several interesting features were identified: (1) An acoustic (dispersive spin-wave) mode develops. (2) The spin-wave bandwidth establishes a new energy scale that is strongly dependent on the size of the magnetic island and becomes abnormally small for large clusters. (3) Optical (dispersionless spin excitation) modes are present for all block states studied here. In addition, a variety of phenomenological spin Hamiltonians have been investigated but none matches entirely our results that were obtained primarily at intermediate Hubbard U strengths. Our comprehensive analysis provides theoretical guidance and motivation to crystal growers to search for appropriate candidate materials to realize the block states, and to neutron scattering experimentalists to confirm the exotic dynamical magnetic properties unveiled here, with a rich mixture of acoustic and optical features.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

Superconductivity in a topological lattice model with strong repulsion

The highly tunable nature of synthetic quantum materials-both in the solid-state and cold atom contexts-invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moiré materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce the "double Hofstadter"model, a minimal 2D lattice model that incorporates exactly these features: Time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. Further, at the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the quantum spin Hall insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be p-wave-reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Intertwined charge, spin, and orbital degrees of freedom under electronic correlations in the one-dimensional Fe 3+ chalcogenide chain

Motivated by recent developments in the study of quasi-one-dimensional iron systems with Fe 2+ , we comprehensively study an Fe 3+ chalcogenide chain system. Based on first-principles calculations, the Fe 3+ chain has a similar electronic structure to that discussed before for the Fe 2+ chain, because of the similar Fe⁢𝑋 4 (𝑋 = S or Se) tetrahedron-chain geometry. Furthermore, a three-orbital electronic Hubbard model for this chain was constructed using the density matrix renormalization group method. A robust antiferromagnetic coupling was unveiled in the chain direction. In addition, in the intermediate electronic correlation 𝑈/𝑊 region, we found an interesting orbital-selective Mott phase with the coexistence of localized and itinerant electrons (𝑈 is the on-site Hubbard repulsion, while 𝑊 is the electronic bandwidth) based on the orbital-selective behavior observed in the charge fluctuations. Furthermore, we do not observe any obvious pairing tendency in the Fe 3+ chain in the electronic-correlation 𝑈/𝑊 region, where superconducting pairing tendencies were reported before in iron ladders. This suggests that superconductivity is unlikely to emerge in the Fe 3+ systems. Finally, our results clearly establish the similarities and differences between Fe 2+ and Fe 3+ iron chains, as well as iron ladders.

density matrix renormalization group↗
DOE OSTI2025

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗
DOE OSTI2023

Traces of electron-phonon coupling in one-dimensional cuprates

Abstract The appearance of certain spectral features in one-dimensional (1D) cuprate materials has been attributed to a strong, extended attractive coupling between electrons. Here, using time-dependent density matrix renormalization group methods on a Hubbard-extended Holstein model, we show that extended electron-phonon ( e–ph ) coupling presents an obvious choice to produce such an attractive interaction that reproduces the observed spectral features and doping dependence seen in angle-resolved photoemission experiments: diminished 3 k F spectral weight, prominent spectral intensity of a holon-folding branch, and the correct holon band width. While extended e–ph coupling does not qualitatively alter the ground state of the 1D system compared to the Hubbard model, it quantitatively enhances the long-range superconducting correlations and suppresses spin correlations. Such an extended e–ph interaction may be an important missing ingredient in describing the physics of the structurally similar two-dimensional high-temperature superconducting layered cuprates, which may tip the balance between intertwined orders in favor of uniform d -wave superconductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hilbert space multireference coupled cluster tailored by matrix product states

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants—the state universal one and the Brillouin–Wigner and Mukherjee’s state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO–LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Chemistry↗
DOE OSTI2025

Anyon superconductivity from topological criticality in a Hofstadter–Hubbard model

We argue that the combination of strong repulsive interactions and high magnetic fields can generate electron pairing and superconductivity. Inspired by the large lattice constants of moiré materials, which make large flux per unit cell accessible at laboratory fields, we study the triangular lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette, where previous literature has argued that a chiral spin liquid separates a weak-coupling integer quantum Hall phase and a strong-coupling topologically trivial antiferromagnetic insulator at a density of one electron per site. We argue that topological superconductivity emerges upon doping in the vicinity of the integer quantum Hall to chiral spin liquid transition. We employ exact diagonalization and density matrix renormalization group methods to examine this theoretical scenario and find that electronic pairing indeed occurs on both sides of criticality over a remarkably broad range of interaction strengths. On the chiral spin liquid side, our results provide a concrete model realization of the long-hypothesized mechanism of anyon superconductivity. Our study thus establishes a beyond-Bardeen-Cooper-Schrieffer route to electron pairing in a well-controlled limit, relying crucially on the interplay between electron correlations and band topology.

Hubbard model↗
DOE OSTI2026

Band-mixing effects in one-dimensional charge transfer insulators

The low-energy properties of transition metal oxides (TMOs) are governed by the electrons occupying strongly correlated $d$-orbitals that are hybridized with surrounding ligand oxygen $p$ orbitals to varying degrees. Their physics is thus established by a complex interplay between the transition-metal (TM)-ligand hopping $t$, charge transfer energy $\Delta_\mathrm{CT}$, and on-site TM Hubbard repulsion $U$. Here, we study the spectral properties of a one-dimensional (1D) analog of such a $pd$ system, with alternating TM $d$ and ligand anion $p$ orbitals situated along a chain. Using the density matrix renormalization group method, we study the model's single-particle spectral function, x-ray absorption spectrum, and dynamical spin structure factor as a function of $\Delta_\mathrm{CT}$ and $U$. In particular, we present results spanning from the Mott insulating ($\Delta_\mathrm{CT} > U$) to negative charge transfer regime $\Delta_\mathrm{CT} < 0$ to gain a better understanding of the ground and momentum-resolved excited state properties of these different regimes. Our results can guide new studies on TMOs that seek to situate them within the Mott-Hubbard/charge transfer insulator classification scheme.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Edge Reconstruction in a Quantum Spin Hall Insulator

We study interaction-driven edge reconstruction in a quantum spin Hall insulator described by the Bernevig-Hughes-Zhang model with Kanamori-Hubbard interactions using the real-space density matrix renormalization group method in both the grand-canonical and canonical ensembles. For a two-dimensional cylinder with a smooth edge, we identify discrete particle-number transitions that lead to a spin-polarized edge state stabilized by an emergent ferromagnetic exchange interaction. The reconstruction is orbital-selective, occurring predominantly in the 𝑠-orbital channel. Our results reveal a microscopic mechanism for emergent fluctuating moments at the edge that could compromise the topological protection of helical edge states by time-reversal symmetry.

Soni, Rahul [ORNL] (ORCID:0000000317714299)↗
DOE OSTI2020

Two-magnon bound states in the Kitaev model in a [111] field

It is now well established that the Kitaev honeycomb model in a magnetic field along the [111] direction harbors an intermediate gapless quantum spin liquid (QSL) phase sandwiched between a gapped non-Abelian QSL at low fields $H < H_{c1}$ and a partially polarized phase at high fields $H > H_{c2}$. Here, we analyze the low-field and high-field phases and phase transitions in terms of single- and two-magnon excitations using exact diagonalization and density matrix renormalization group methods. In this work, we find that the energy to create a bound state of two-magnon $Δ_p$ becomes lower than the energy to create a single spin flip $Δ_s$ near $H_{c2}$. In the entire Kitaev spin liquid, $Δ_p < Δ_s$ and both gaps vanish at $H_{c2}$. We make testable predictions for magnon pairing that could be observable in Raman scattering measurements on Kitaev QSL candidate materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spectroscopic signatures of next-nearest-neighbor hopping in the charge and spin dynamics of doped one-dimensional antiferromagnets

In this work, we study the impact of next-nearest-neighbor (nnn) hopping on the low-energy collective excitations of strongly correlated doped antiferromagnetic cuprate spin chains. Specifically, we use exact diagonalization and the density matrix renormalization group method to study the single-particle spectral function, the dynamical spin and charge structure factors, and the Cu L-edge resonant inelastic x-ray scattering (RIXS) intensity of the doped t-t'-J model for a set of t' values. We find evidence for the breakdown of spin-charge separation as |t'| increases and identify its fingerprints in the dynamical response functions. The inclusion of nnn hopping couples the spinon and holon excitations, resulting in the formation of a spin-polaron, where a ferromagnetic spin polarization cloud dresses the doped carrier. The spin-polaron manifests itself as additional spectral weight in the dynamical correlation functions, which appear simultaneously in the spin- and charge-sensitive channels. We also demonstrate that RIXS can provide a unique view of the spin-polaron, due to its sensitivity to both the spin and charge degrees of freedom.

1D spin chains↗
DOE OSTI2020

Coupled dimer and bond-order-wave order in the quarter-filled one-dimensional Kondo lattice model

Motivated by experiments on the organic compound (Per) 2 [Pt(mnt) 2 ], we study the ground state of the onedimensional Kondo lattice model at quarter filling with the density matrix renormalization group method. Here, we show a coupled dimer and bond-order-wave (BOW) state in the weak-coupling regime for the localized spins and itinerant electrons, respectively. The quantum phase transitions for the dimer and the BOW orders occur at the same critical coupling parameter J c , with the opening of a charge gap. The emergence of the combination of dimer and BOW order agrees with the experimental findings of the simultaneous Peierls and spin-Peierls transitions at low temperatures, which provides a theoretical understanding of such a phase transition. We also show that the localized spins in this insulating state have quasi-long-ranged spin correlations with collinear configurations, which resemble the classical dimer order in the absence of a magnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Peierls transition, ferroelectricity, and spin-singlet formation in monolayer VOI 2

Using ab initio density functional theory and single-orbital Hubbard model calculations via the density matrix renormalization group method, we systematically studied the monolayer VOI 2 with a 3d 1 electronic configuration. Our phonon calculations indicate that the orthorhombic Pmm2 FE-II phase is the most likely ground state, involving a ferroelectric (FE) distortion along the a axis and V-V dimerization along the b axis. Specifically, the “pseudo Jahn-Teller” effect caused by the coupling between empty V (d xz/yz and d 3z 2 -r 2 ) and O 2p states is proposed as the mechanism that stabilizes the FE distortion from the paraelectric phase. Moreover, the half-filled metallic d xy band displays a Peierls instability along the b axis, inducing a V-V dimerization. Additionally, we found very short-range antiferromagnetic coupling along the V-V chain due to the formation of nearly decoupled spin singlets in the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ferromagnetism and doublon localization in a Wannier-Hubbard chain

We derive a “Wannier-Hubbard” model consisting of an array of overlapping atomic orbitals interacting via a local Coulomb interaction. Transforming to an orthogonal Wannier basis set, the resulting Hamiltonian displays long-range hopping and interactions, with new terms such as correlated hopping and ferromagnetic direct exchange, among others. We numerically study the one-dimensional version of the model at half-filling using the density matrix renormalization group method, unveiling a rich phase diagram as a function of the interaction U and the overlap s with metallic and ferromagnetic phases separated by a ferrimagnetic region. Our results indicate a path toward understanding emergent phases under pressure and beyond standard model Hamiltonians.

36 MATERIALS SCIENCE↗