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 OSTI2021

Nonequilibrium optical response of a one-dimensional Mott insulator

We define, compute, and analyze the nonequilibrium differential optical conductivity of the one-dimensional extended Hubbard model at half-filling after applying a pump pulse, using the time-dependent density matrix renormalization group method. The melting of the Mott insulator is accompanied by a suppression of the local magnetic moment and ensuing photogeneration of doublon-holon pairs. The differential optical conductivity reveals (i) mid-gap states related to parity-forbidden optical states and (ii) strong renormalization and hybridization of the excitonic resonance and the absorption band, yielding a Fano resonance. We offer evidence and interpret such a resonance as a signature of nonequilibrium optical excitations resembling excitonic strings, (bi)excitons, and unbound doublon-holon pairs, depending on the magnitude of the intersite Coulomb repulsion. We discuss our results in the context of pump and probe spectroscopy experiments on organic Mott insulators.

1-dimensional systems↗
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy-particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). In this work, we apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

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

Intertwined charge, spin, and pairing orders in doped iron ladders

Motivated by recent experimental progress on iron-based ladder compounds, we study the doped two-orbital Hubbard model for the two-leg ladder BaFe2S3. The model is constructed by using ab initio hopping parameters and the ground state properties are investigated using the density matrix renormalization group method. We show that the (π,0) magnetic ordering at half filling, with ferromagnetic rungs and antiferromagnetic legs, becomes incommensurate upon hole doping. Moreover, depending on the strength of the Hubbard U coupling, other magnetic patterns, such as (0,π), are also stabilized. We found that the binding energy for two holes becomes negative for intermediate Hubbard interaction strength, indicating hole pairing. Due to the crystal-field split among orbitals, the holes primarily reside in one orbital, with the other one remaining half filled. This resembles orbital selective Mott states. The formation of tight hole pairs continues with increasing hole density, as long as the magnetic order remains antiferromagnetic in one direction. The study of pair-pair correlations indicates the dominance of the intraorbital spin-singlet channel, as opposed to other pairing channels. Although in a range of hole doping pairing correlations decay slowly, our results can also be interpreted as corresponding to a charge density wave made of pairs, a precursor of eventual superconductivity after interladder couplings are included. Such a scenario of intertwined orders has been extensively discussed before in the cuprates, and our results suggest a similar physics could exist in ladder iron-based superconductors. Finally, we also show that a robust Hund's coupling is needed for pairing to occur.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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)↗