arXiv:2302.14313

Finite Temperature Simulations of Strongly Correlated Systems

Chong Sun

Published Dec 18, 2020• 1 citation

Abstract

This thesis describes several topics related to finite temperature studies of strongly correlated systems: finite temperature density matrix embedding theory (FT-DMET), finite temperature metal-insulator transition, and quantum algorithms including quantum imaginary time evolution (QITE), quantum Lanczos (QLanczos), and quantum minimally entangled typical thermal states (QMETTS) algorithms. While the absolute zero temperature is not reachable, studies of physical and chemical problems at finite temperatures, especially at low temperature, is essential for understanding the quantum behaviors of materials in realistic conditions. Here we define low temperature as the temperature regime where the quantum effect is not largely dissipated due to thermal fluctuation. Treatment of systems at low temperature is specially difficult compared to both high temperature - where classical approximation can be applied - and zero temperature where only the ground state is required to describe the system of interest. FT-DMET is a wavefunction-based embedding scheme which can handle finite temperature simulations of a variety of strongly correlated problems. The "high-level in low-level" framework enables FT-DMET to tackle large bulk sizes and capture the majority of the entanglement at the same time. FT-DMET formulations and implementation details for both model systems and ab initio problems are provided in Chapter 2 and Chapter 3. Metal-insulator transition is a common but important phase transition in many strongly correlated materials. The widely accepted scheme to distinguish an insulator from a metal is band structure theory based on a single-particle picture. However, insulating phases caused by disorder or strong correlation cannot be explained merely with the band structure. In Chapter 4, we demonstrate that electron locality/mobility is a more general criteria to detect metal-insulator transition. We further introduce complex polarization as the order parameter to reflect the electron locality/mobility and provide a formalism based on thermofield theory to evaluate the complex polarization at finite temperature. Quantum algorithms are designed to perform simulations on a quantum device. The infrastructure of a quantum processing unit (QPU) utilizes the superposition property of quantum bits (qubits), and thus can potentially outplay the classical simulations in computational scaling for certain problems. In Chapter 5, we introduce the QITE algorithm, which can be applied to quantum simulations of both ground state and finite temperature problems. We further introduce a subspace method, QLanczos algorithm, and a a finite temperature quantum algorithm, QMETTS, where QITE is used as a building block for the two algorithms. We demonstrate above quantum algorithms with simulations on both classical computers and quantum computers.

DOI: 10.48550/arxiv.2302.14313

Related papers

Ab initio quantum embedding at finite temperature with density matrix embedding theory

Laurence W. Giordano, Y. Stanley Tan, Zhi‐Hao Cui +1 more

We present a finite-temperature extension of density matrix embedding theory (FT-DMET) for realistic crystalline systems. We describe a practical framework for constructing extended bath orbitals, solving the embedding problem, and performing DMET self-consistency at finite temperature. To reduce computational cost, we introduce strategies based on mutual-information-guided bath truncation, controlled treatment of the thermal electron number without explicit optimization, and the use of low-temperature impurity solvers and one-shot FT-DMET in the low-temperature regime. We apply this approach to periodic hydrogen chains and square lattices to characterize their finite-temperature phases. We observe the Pomeranchuk-like effect in one dimension and enhanced stability of long-range order in two dimensions.

The Journal of Chemical PhysicsApr 15, 20261 citation

Ab initio quantum embedding at finite temperature with density matrix embedding theory

Laurence W. Giordano, Y. Stanley Tan, Zhi‐Hao Cui +1 more

arXiv:2601.01641Jan 4, 20260 citations

Density Matrix Embedding Theory and Strongly Correlated Lattice Systems

Bo-Xiao Zheng

This thesis describes the development of the density matrix embedding theory (DMET) and its applications to lattice strongly correlated electron problems, including a review of DMET theory and algorithms (Ch 2), investigation of finite size scaling (Ch 3), Applications to high-temperature superconductivity (Ch 4-6), a framework for finite-temperature DMET (Ch 7).

arXiv:1803.10259Mar 27, 20184 citations

Algorithms for discrete and continuous quantum systems

Soumyodipto Mukherjee

This thesis is divided into three chapters. In the first chapter we outline a simple and numerically inexpensive approach to describe the spectral features of the single-impurity Anderson model. The method combines aspects of the density matrix embedding theory (DMET) approach with a spectral broadening approach inspired by those used in numerical renormalization group (NRG) methods. At zero temperature for a wide range of U , the spectral function produced by this approach is found to be in good agreement with general expectations as well as more advanced and complex numerical methods such as DMRG-based schemes. The theory developed here is simply transferable to more complex impurity problems The second chapter outlines the density matrix embedding methodology in the context of electronic structure applications. We formulate analytical gradients for energies obtained from DMET focusing on two scenarios: RHF-in-RHF embedding and FCI-in-RHF embedding. The former involves solving the small embedded system at Restricted Hartree-Fock (RHF) level of theory. This serves to check the validity of the formulas by reproducing the RHF results on the full system for energies and gradients. The latter scenario employs full configuration interaction (FCI) as a high level solver for the small embedded system. Our results show that only Hellmann-Feynman terms, which involve derivatives of one- and two-electron terms in the atomic orbital basis, are required to calculate energy gradients in both cases. We applied our methodology to the problem of H 10 ring dissociation where the analytical gradients matched those obtained numerically. The gradient formulation is applicable to geometry optimization of strongly correlated molecules and solids. It can also be used in ab initio molecular dynamics where forces on nuclei are obtained from DMET energy gradients. In the final chapter we focus on the study of finite-temperature equilibrium properties of quantum systems in continuous space. We formulate an expansion of the partition function in continuous-time and use Monte Carlo to sample terms in the resulting infinite series. Such a strategy has been highly successful in quantum lattice models but has found scant application in off-lattice systems. The Monte Carlo estimate of the average energy of quantum particles in continuous space subject to simple model potentials is found to converge with low statistical error to the exact solutions even when very high perturbation order is required. We outline two ways in which the algorithm can be applied to more complex problems. First, by drawing an analogy between the formulation in continuous-time with the discretized, Trotter factorized version of standard path integral Monte Carlo (PIMC). This allows one to use the suite of standard PIMC moves to carry out the position sampling required to obtain the weight of each time configuration. Finally, we propose an alternate route by fitting the many-particle potential with multidimensional Gaussians which provides an analytical form for the position integrals.

Jan 1, 20180 citations

Scalable Ab Initio Quantum Many-Body Methods for Crystalline Materials

Junjie Yang

Achieving materials-specific predictions for large, realistic molecules and solids with strong electron correlations remains a long-standing challenge in quantum chemistry. This dissertation develops scalable ab initio quantum many-body methods for crystalline materials, demonstrating that rigorous optimization of theoretical approximations and numerical algorithms enables predictive computation even in strongly correlated systems. The work tackles three interconnected challenges, namely describing material-specific properties of strongly correlated superconductors, achieving linear scaling with k-points for large-scale periodic calculations, and establishing reliable methods for electron-phonon interactions. The following chapters address each challenge in turn. Chapter 2 demonstrates that ab initio density matrix embedding theory, combined with symmetry-breaking strategies for superconducting order parameters, reproduces experimental trends in cuprate superconductivity, including pressure and layer effects on transition temperatures. By directly solving the electronic Schrödinger equation from material structures, these calculations capture multi-orbital covalency effects that simplified model Hamiltonians neglect, revealing their essential role in the spin fluctuations driving pairing. This predictive capability establishes a foundation for efficient computational screening of superconducting materials. Extending these calculations to larger k-point meshes and broader phase diagram explorations, however, demands more efficient computational algorithms. To overcome this computational bottleneck, Chapter 3 introduces interpolative separable density fitting framework. This development enables thermodynamic limit calculations with up to 1000 k-points, as validated on diamond, carbon dioxide, nickel monoxide, and cuprate systems. Integration with density matrix embedding theory and local natural orbital methods further yields converged ground-state energies from accurate correlated wavefunction methods such as CCSD(T) for crystalline materials in the thermodynamic limit. With computational efficiency established, the remaining challenge lies in extending the theoretical framework beyond purely electronic degrees of freedom. Chapter 4 addresses this need by establishing exponential ansatz wavefunctions as competitive approaches for electron-phonon systems. Unlike previous treatments based on perturbation theory, systematic benchmarking of coupled cluster theory and variational Lang-Firsov methods on the Holstein model demonstrates, for the first time, that these approaches accurately describe polaron formation across coupling regimes when combined with appropriate reference state optimization. These results provide a foundation for incorporating phonon effects into ab initio studies of superconducting materials. In addition to these methodological advances, GPU acceleration and MPI parallelization enhance computational throughput, enabling efficient large-scale calculations. Together, these developments constitute a computational framework for predictive ab initio quantum many-body calculations in crystalline materials, opening new avenues for understanding and designing strongly correlated quantum materials.

Mar 5, 20260 citations

Quantum Embedding and Localization Methods for Strongly Correlated Molecules and Materials

Shreya Verma

Strongly correlated electronic systems underpin a wide range of emerging quantum technologies, including solid-state defects for quantum information processing, molecular magnets, spin qubits, and quantum sensing platforms. Predictive computational design of these materials requires electronic structure methods capable of simultaneously describing strong electron correlation, excited-state phenomena, spin-dependent interactions, and environmental effects. While multireference wave function methods provide the necessary accuracy, their computational cost severely limits their application to realistic molecular and condensed-phase systems. This dissertation addresses this challenge through the development and application of quantum embedding and fragmentation methodologies that extend multireference electronic structure theory to increasingly complex quantum materials. The first part of this dissertation develops embedding-based approaches for describing optical properties of strongly correlated defects in extended solids. Complete Active Space Density Matrix Embedding Theory (CAS-DMET) and NEVPT2-DMET are applied to oxygen vacancy defects in magnesium oxide, where a systematic double-extrapolation protocol is introduced to approach both the non-embedding and thermodynamic limits. This framework enables quantitatively reliable calculations of absorption and emission energies while significantly reducing the computational cost of conventional multireference treatments of periodic systems. The second part extends embedding methodologies to the description of spin-dependent electronic structure. Density Matrix Embedding Pair-Density Functional Theory (DME-PDFT) is developed to recover dynamic electron correlation beyond the embedded active space, providing an efficient route for studying transition metal complexes and other strongly correlated molecular systems. Building upon these developments, CAS-DMET is further extended to investigate spin-phonon relaxation in molecular magnets and molecular crystals, demonstrating that quantum embedding can accurately describe spin dynamics and decoherence processes while substantially reducing computational cost. Beyond classical electronic structure calculations, this dissertation explores fragmentation strategies for hybrid quantum-classical algorithms through the development of efficient methodologies within the Localized Active Space Unitary Selective Coupled Cluster (LAS-USCC) framework. These developments improve the scalability of quantum chemistry algorithms for large active spaces and establish connections between embedding concepts and emerging quantum computing approaches. Finally, electronic structure calculations based on time-dependent density functional theory are employed to investigate fluorescent protein chromophores as potential biocompatible quantum sensors. These studies demonstrate how first-principles calculations can provide insight into spin-dependent properties in biologically relevant systems and expand the scope of computational approaches for quantum materials beyond traditional inorganic and molecular platforms. Collectively, the work presented in this dissertation establishes embedding and fragmentation as versatile frameworks for accurately describing the optical and spin properties of strongly correlated quantum systems across multiple length scales. By combining multireference electronic structure methods with physically motivated subsystem decompositions,these approaches bridge the gap between chemical accuracy and computational scalability, providing practical computational tools for the study and design of quantum materials relevant to quantum information science, molecular spintronics, and quantum sensing.

Aug 1, 20260 citations