Hybridization expansion Monte Carlo simulation of multi-orbital quantum impurity problems: matrix product formalism and improved Monte Carlo sampling

Abstract

We explore two complementary modifications of the hybridization-expansion continuous-time Monte Carlo method, aiming at large multi-orbital quantum impurity problems. One idea is to compute the imaginary-time propagation using a matrix product states representation. We show that bond dimensions considerably smaller than the dimension of the Hilbert space are sufficient to obtain accurate results, and that this approach scales polynomially, rather than exponentially with the number of orbitals. Based on scaling analyses, we conclude that a matrix product state implementation will outperform the exact-diagonalization based method for quantum impurity problems with more than 12 orbitals. The second idea is an improved Monte Carlo sampling scheme which is applicable to all variants of the hybridization expansion method. We show that this so-called sliding window sampling scheme speeds up the simulation by at least an order of magnitude for a broad range of model parameters, with the largest improvements at low temperature.

arXiv comment24 pages, 8 figures

References in corpus (19)

• The density-matrix renormalization group in the age of matrix product states
• Electronic Structure Calculations with Dynamical Mean-Field Theory: A Spectral Density Functional Approach
• Continuous-time Monte Carlo methods for quantum impurity models
• Quantum Cluster Theories
• A continuous-time solver for quantum impurity models
• Thermodynamic limit of the density matrix renormalization for the spin-1 Heisenberg chain
• Continuous Time Quantum Monte Carlo method for fermions
• Quantum Monte Carlo Impurity Solver for Cluster DMFT and Electronic Structure Calculations in Adjustable Base
• Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models
• Screening and Non-local Correlations in the Extended Hubbard Model from Self-Consistent Combined GW and Dynamical Mean Field Theory
• Markov Chain Monte Carlo Method without Detailed Balance
• Dynamical Screening in Correlated Electron Materials
• Conserved quantities of SU(2)-invariant interactions for correlated fermions and the advantages for quantum Monte Carlo simulations
• Density-matrix renormalisation group approach to quantum impurity problems
• Krylov implementation of the hybridization expansion impurity solver and application to 5-orbital models
• High Energy Dynamics of the Single Impurity Anderson Model
• Spectral Densities from Dynamic Density-Matrix Renormalization
• Spectral functions for single- and multi-Impurity models using DMRG
• Spectral Density of the Two-Impurity Anderson Model

Cited by in corpus (16)

• w2dynamics: Local one- and two-particle quantities from dynamical mean field theory
• Matrix Product State applications for the ALPS project
• $i$QIST: An open source continuous-time quantum Monte Carlo impurity solver toolkit
• Imaginary-time matrix product state impurity solver for dynamical mean-field theory
• Continuous-time hybridization expansion quantum impurity solver for multi-orbital systems with complex hybridizations
• Lazy skip lists, a new algorithm for fast hybridization-expansion quantum Monte Carlo
• A spin-freezing perspective on cuprates
• Symmetric Improved Estimators for Continuous-time Quantum Monte Carlo
• Efficient Continuous-time Quantum Monte Carlo Method for the Ground State of Correlated Fermions
• A Tensor Train Continuous Time Solver for Quantum Impurity Models
• Recommender Engine for Continuous Time Quantum Monte Carlo Methods
• State- and superstate-sampling in hybridization-expansion continuous-time quantum Monte Carlo
• Double-expansion impurity solver for multiorbital models with dynamically screened U and J
• Accelerated impurity solver for DMFT and its diagrammatic extensions
• Paramagnon dispersion and damping in doped Na${}_x$Ca${}_{2-x}$CuO${}_2$Cl${}_2$
• CT-X: an efficient continuous-time quantum Monte Carlo impurity solver in Kondo Regime