Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix renormalization group”

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

Design, Control and Application of Next Generation Qubits

Design, Control and Application of Next Generation Qubits Arun Bansil, Northeastern University (Principal Investigator) Claudio Chamon, Boston University (Co-Investigator) Adrian Feiguin, Northeastern University (Co-Investigator) Liang Fu, MIT (Co-Investigator) Eduardo Mucciolo, Univ. of Central Florida (Co-Investigator) Qimin Yan, Temple University (Co-Investigator) The quest for developing technologies for manipulating and storing information quantum mechanically is currently led by approaches that include Josephson-junctions, ion-traps, and qubits generated by defect spins in solids. Topological qubits, however, are inherently more robust to decoherence by environmental effects, and should be able to sprint ahead once practical barriers have been overcome. At the present stage of the development of the field, it is important to explore a variety of architectures and materials beyond the conventional paradigms in order to seed breakthroughs toward building a scalable quantum computer. Our comprehensive theoretical research program involved four interconnected thrusts as follows. • A materials discovery effort in two-dimensional compounds in search of materials to support Majorana zero modes and defect structures suitable as qubits. • Exploration of architectures for topological quantum computation by investigating both superconducting Majorana qubits, and robust platforms for braiding with new “meta-materials” built of arrays of Majorana qubits. • Investigation of properties of hybrid metal-organic qubits based on transition-metal centers in graphene, and molecular crystals of polyaromatic complexes with embedded transition-metal atoms. • Development of tensor-network and semiclassical approaches to study decoherence in the presence of random and dispersive spin baths, and NV centers in diamond. The full spectrum of theoretical and numerical approaches was used to address the goals of this project including first-principles, density-matrix-renormalization group, tensor networks, and data-driven high-throughput approaches using materials database and machine-learning.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A Dipolar Chiral Spin Liquid on the Breathed Kagome Lattice

Continuous control over lattice geometry, when combined with long-range interactions, offers a powerful yet underexplored tool to generate highly frustrated quantum spin systems. By considering long-range dipolar antiferromagnetic interactions on a breathed Kagome lattice, we demonstrate how these tools can be leveraged to stabilize a chiral spin liquid. We support this prediction with large-scale density-matrix renormalization group calculations and explore the surrounding phase diagram, identifying a route to adiabatic preparation via a locally varying magnetic field. At the same time, we identify the relevant low-energy degrees of freedom in each unit cell, providing a complementary language to study the chiral spin liquid. Finally, we carefully analyze its stability and signatures in finite-sized clusters, proposing direct, experimentally viable measurements of the chiral edge mode in both Rydberg atom and ultracold polar molecule arrays.

Machado↗
NASA NTRS1999

Phase Diagram of the Two-Chain Hubbard Model

We have calculated the charge gap and spin gap for the two-chain Hubbard model as a function of the on-site Coulomb interaction and the interchain hopping amplitude. We used the density matrix renormalization group method and developed a method to calculate separately the gaps numerically for the symmetric and antisymmetric modes with respect to the exchange of the chain indices. We have found very different behaviors for the weak and strong interaction cases. Our calculated phase diagram is compared to the one obtained by Balents and Fisher using the weak coupling renormalization group technique.

Park, Youngho↗
NASA NTRS1999

Charge and Spin Dynamics of the Hubbard Chains

We calculate the local correlation functions of charge and spin for the one-chain and two-chain Hubbard model using density matrix renormalization group method and the recursion technique. Keeping only finite number of states we get good accuracy for the low energy excitations. We study the charge and spin gaps, bandwidths and weights of the spectra for various values of the on-site Coulomb interaction U and the electron filling. In the low energy part, the local correlation functions are different for the charge and spin. The bandwidths are proportional to t for the charge and J for the spin respectively.

Park, Youngho↗
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 OSTI2019

Real-time density-matrix coupled-cluster approach for closed and open systems at finite temperature

We extend the coupled-cluster method to correlated quantum dynamics of both closed and open systems at finite temperatures using the thermofield formalism. The approach expresses the time-dependent density matrix in an exponential ansatz and describes time-evolution along the Keldysh path contour. A distinct advantage of the approach is exact trace-preservation as a function of time, ensuring conservation of probability and particle number. Furthermore, the method avoids the computation of correlated bra-states, simplifying the computational implementation. We develop the method in a thermal quasiparticle representation, which allows seamless connection to the projection method and diagrammatic techniques of the traditional coupled-cluster formalism. For comparison, we also apply the thermofield framework to the density-matrix renormalization-group method to obtain reference results for closed and open systems at finite temperature. We test the singles and doubles approximation to the density-matrix coupled-cluster method on the correlated electronic dynamics of the single-impurity Anderson model, demonstrating that the new method successfully captures the correlated dynamics of both closed systems at finite temperature and driven-dissipative open systems. This encouraging performance motivates future applications to nonequilibrium quantum many-body dynamics in realistic systems.

Shushkov, Philip (ORCID:0000000266136117)↗
DOE OSTI2019

Time- and momentum-resolved tunneling spectroscopy of pump-driven nonthermal excitations in Mott insulators

We present a computational technique to calculate time- and momentum-resolved nonequilibrium spectral density of correlated systems using a tunneling approach akin to scanning tunneling spectroscopy. The important difference is that our probe is extended, basically a copy of the sample, allowing one to extract the momentum information of the excitations. We illustrate the method by measuring the spectrum of a Mott-insulating extended Hubbard chain after a sudden quench with the aid of time-dependent density-matrix renormalization-group calculations. In this work, we demonstrate that the system realizes a nonthermal state that is an admixture of spin and charge density wave states, with corresponding signatures that are recognizable as in-gap subbands. In particular, we identify a band of excitons and one of stable antibound states at high energies that gains enhanced visibility after the pump. We do not appreciate noticeable relaxation within the time scales considered, which is attributed to the lack of decay channels due to spin-charge separation. These ideas can be readily applied to study transient dynamics and spectral signatures of correlation-driven nonequilibrium processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quasi-Fermi liquid behavior in a one-dimensional system of interacting spinless fermions

We present numerical evidence for a paradigm in one-dimensional interacting fermion systems, whose phenomenology has traits of both Luttinger liquids and Fermi liquids. This state, dubbed a quasi-Fermi liquid, possesses a discontinuity in its fermion occupation number at the Fermi momentum. The excitation spectrum presents particlelike quasiparticles and absence of holelike quasiparticles, giving rise instead to edge singularities. Such a state is realized in a one-dimensional spinless fermion lattice Hamiltonian by fine-tuning the interactions to a regime where they become irrelevant in the renormalization group sense. We show, using uniform infinite matrix products states and finite-entanglement scaling analysis, that the system ground state is characterized by a Luttinger parameter K = 1 and a discontinuous jump in the fermion occupation number. We support the characterization with calculations of the spectral function that show a particle-hole asymmetry reflected in the existence of well-defined Landau quasiparticles above the Fermi level and edge singularities without the associated quasiparticles below. Furthermore, these results indicate that the quasi-Fermi liquid paradigm can be realized beyond the low-energy perturbative realm.

1-dimensional systems↗
DOE OSTI2022

Electronic structure, magnetic properties, and pairing tendencies of the copper-based honeycomb lattice Na 2 Cu 2 TeO 6

Spin-1/2 chains with alternating antiferromagnetic (AFM) and ferromagnetic (FM) couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density-matrix renormalization-group (DMRG) methods, we have systematically studied the dimerized chain system Na 2 Cu 2 TeO 6 with a d 9 electronic configuration. Near the Fermi level, in the nonmagnetic phase the dominant states are mainly contributed by the Cu 3d x 2 -y 2 orbitals highly hybridized with the O 2p orbitals, leading to an “effective” single-orbital low-energy model. By calculating the relevant hoping amplitudes, we explain the size and sign of the exchange interactions in Na 2 Cu 2 TeO 6 . In addition, a single-orbital Hubbard model is constructed for this dimerized chain system where the quantum fluctuations are taken into account. Both AFM and FM couplings (leading to an ↑|-|↓-|↓-↑| state) along the chain were found in our DMRG and Lanczos calculations, in agreement with density functional theory and neutron-scattering results. The hole pairing binding energy ΔE is predicted to be negative at Hubbard U~ 11eV, suggesting incipient pairing tendencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Higher Hall conductivity from a single wave function: Obstructions to symmetry-preserving gapped edge of (2+1)-dimensional topological order

A (2+1)D topologically ordered phase with U(1) symmetry may or may not have a symmetric gapped edge state, even if both thermal and electric Hall conductivity are vanishing. It has recently been discovered that there are “higher” versions of Hall conductivity valid for fermionic fractional quantum Hall (FQH) states that obstruct symmetry-preserving gapped edge states beyond thermal and electric Hall conductivity. In this paper, we show that one can extract higher Hall conductivity from a single wave function of an FQH state, by evaluating the expectation value of the “partial rotation” unitary, which is a combination of partial spatial rotation and a U(1) phase rotation. This result is verified numerically with the fermionic Laughlin state with 𝜈=1/3 and 1/5, as well as the non-Abelian Moore-Read state. Together with topological entanglement entropy, we prove that the expectation values of the partial rotation completely determine if a bosonic/fermionic Abelian topological order with U(1) symmetry has a symmetry-preserving gappable edge state or not. We also show that thermal and electric Hall conductivity of Abelian topological order can be extracted by partial rotations. Even in non-Abelian FQH states, partial rotation provides the Lieb-Schultz-Mattis type theorem constraining the low-energy spectrum of the bulk-boundary system. The generalization of higher Hall conductivity to the case with Lie group symmetry is also presented.

2-dimensional systems↗
DOE OSTI2019

Broken mirror symmetry, incommensurate spin correlations, and B 2 g nematic order in iron pnictides

Motivated by recent experiments in the extremely hole doped iron pnictide compounds AFe 2 As 2 (A = K, Rb, Cs), we consider spin-driven nematic order for incommensurate magnetic fluctuations. We classify the nematic order parameters by broken mirror symmetries of the tetragonal D 4h point group, and use this scheme to construct a general Ginzburg-Landau theory that links the nematic order to spatial pattern of magnetic fluctuations. Our analysis points to incommensurate (q,q) magnetic fluctuations as underlying a B 2 g nematic order in AFe 2 As 2 . We substantiate this idea by microscopic calculations based on three-sublattice (2π/3,2π/3) spin correlations in an extended bilinear-biquadratic Heisenberg model. Our classification scheme provides symmetry-based understanding for quasidegeneracy of several nematic channels. The proposed mechanism resolves recently emerged experimental puzzles. Furthermore, we suggest ways to further test it in future experiments and discuss the implications of our results for iron-based high temperature superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantum phases in the honeycomb-lattice J 1 – J 3 ferro-antiferromagnetic model

The significant role of quantum effects in magnets with competing interactions remains at the forefront of condensed matter physics for over 50 years, inspiring a multitude of quests for exotic states, models that can realize them, and real materials that can host them. The elusive spin-liquid states with strong entanglement are but one example; others include valence-bond phases with spatial symmetry breaking, quantum spin nematics that are quantum analogues of liquid crystals, and an especially extensive class of unconventional magnetically ordered phases that do not appear in the classical solutions of the underlying spin models. It is the latter group of phenomena that creates a broader context for the research effort in the present work concerning the search for and an identification of the novel ground states of quantum magnets. In this work, we have studied one of the paradigmatic models in quantum magnetism that is also attracting a significant recent interest because it appears to be providing a tantalizingly close description for many of the newly synthesized materials. We demonstrated that the phase diagram of the quantum ferro-antiferromagnetic model on the honeycomb lattice differs dramatically from the classical one. It hosts the double-zigzag and Ising-z phases as unexpected intermediaries between ferromagnetic and zigzag states that are also extended beyond their classical regions of stability. In broad agreement with quantum order-by-disorder arguments, these collinear phases completely supersede the classical spiral state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Recent developments in the PySCF program package

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. Here, we explain the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

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

A multilayer multi-configurational approach to efficiently simulate large-scale circuit-based quantum computers on classical machines

Here, a multilayer multi-configurational theory framework is adapted to simulate circuit-based quantum computers. Quantum addition of superpositions of an exponential number of summands is performed in polynomial time with high accuracy. We demonstrate numerically accurate calculations including up to one million qubits for entangling benchmarks. Simulation cost can be assessed by entropy-based entanglement measures. For the considered systems, we show that the entanglement only grows weakly with the system size. The present simulations demonstrate how quantum algorithms in low-entropy regimes can be used efficiently on classically simulated quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

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

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

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

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗