A numerically implementable multi-scale many-body approach to strongly correlated electron systems is introduced. An extension to quantum cluster methods, it approximates correlations on any given length-scale commensurate with the strength of the correlations on the respective scale. Short length-scales are treated explicitly, long ones are addressed at a dynamical mean-field level and intermediate length-regime correlations are assumed to be weak and are approximated diagrammatically. To illustrate and test this method, we apply it to the one-dimensional Hubbard model. The resulting multi-scale self-energy provides a very good quantitative agreement with substantially more numerically expensive, explicit quantum Monte Carlo calculations.
A numerically implementable multi-scale many-body approach to strongly correlated electron systems is introduced. An extension to quantum cluster methods, it approximates correlations on any given length-scale commensurate with the strength of the correlations on the respective scale. Short length-scales are treated explicitly, long ones are addressed at a dynamical mean-field level and intermediate length-regime correlations are assumed to be weak and are approximated diagrammatically. To illustrate and test this method, we apply it to the one-dimensional Hubbard model. The resulting multi-scale self-energy provides a very good quantitative agreement with substantially more numerically expensive, explicit quantum Monte Carlo calculations.
Many technologically important materials display complex behavior due to strong electronic correlations. These include magnets, magnetoresistive materials, lanthanides and actinides and high temperature superconductors. Some of these materials display complex phase diagrams with multiple competing ground states, spin and charge ordering and complex excitations like those associated with spin and charge separation. An unbiased treatment of these phenomena requires a theory able to treat the correlations over all relevant length scales. In part, due to the minus sign problem, numerically exact methods employed to study these systems scale exponentially. Thus, despite great recent progress such as the development of the dynamical mean field approximation and its cluster extensions, there is as yet no unbiased study of these complex phenomena. The Multi-Scale Many-Body (MSMB) formalism and algorithms address these problems by separating the correlations into different length scales and treating each with an appropriate approximation: strong correlations at short length scales are treated by numerically exact Quantum Monte Carlo (QMC) simulations, those at intermediate length scales are treated by Feynman diagrammatic methods, while the weak correlations at long length scales are treated by a (dynamical) mean field approximation with a self-consistently determined effective medium.
Abstract The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for their technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory and its cluster extension, the dynamical cluster approximation have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance aiming to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. We will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.
The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe a framework for a general multiscale method based on embedding a self-energy of a strongly correlated subsystem into a self-energy generated by a method able to treat large weakly correlated systems approximately. As an example, we present the embedding of an exact diagonalization self-energy into a self-energy generated from self-consistent second-order perturbation theory. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the two-dimensional Hubbard model, as a benchmark, we illustrate that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the method in multiple regimes of interaction strengths and dopings of the model. The general embedding framework we present avoids difficulties such as double counting corrections, frequency-dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.
Three different approaches to correlated electron systems are presented.(i)The single impurity Anderson model is solved by a semi-analytical approach based on the flow equation method.(ii)A numerically implementable Multi-scale Many-Body approach to strongly correlated electron systems is introduced.An extension to quantum cluster methods, it approximates correlations on any given length-scale commensurate with the strength of the correlations on the respective scale.(iii) The spectral properties of the one-dimensional Holstein and breathing polaron models using the self-consistent Born approximation are presented.4.10 a) Energy levels for zero coupling (g = 0).The solid line shows the non-interacting electron dispersion, the dashed ones are one phonon + one electron states and the dotted one two phonons + one electron, etc. k 1 is the momentum where the electron dispersion reaches ω 0 and has the most significant contribution to the scattering at the first phonon threshold.k 2 is the value of momentum with significant contribution to the scattering at the second phonon threshold.b) For the polaron with momentum k the scattering with momentum q has a significant contribution.Although the polaron momentum is close to k F , thus k a = k -k F is small, q = k a + k b is large because k b which connects two points on the Fermi surface can be large. . . . . . . . . . .
This thesis represents our effort to develop the next generation multi-scale quantum simulation methods suitable for strongly-correlated materials, where complicated phase-diagrams prevail, suggesting complicated underlying physics. We first give a detailed description of the parquet formalism. With its help, different approximate methods can be unified and a hierarchy of approximate methods with different accuracies and computational complexity can thus be designed. Next, we present a numerical solution of the parquet approximation. Results on the Hubbard model are compared to those obtained from Determinant Quantum Monte Carlo (DQMC), FLuctuation EXchange (FLEX), and self-consistent second-order approximation methods. The comparison shows a satisfactory agreement with DQMC and a significant improvement over the FLEX or the self-consistent second-order approximation. The parquet formalism can also be used to analyze the superconducting mechanism of the high-temperature superconductors. The dynamical cluster approximation (DCA) method is used to understand the proximity of the superconducting dome to the quantum critical point in the 2-D Hubbard model. At optimal doping, where Vd is revealed to be featureless, we find a power-law behavior of chi_0d(w=0), replacing the BCS logarithm behavior, and strongly enhanced T_c. After that we propose another multi-scale approach by combining the DCA and the recently introduced dual-fermion formalism. Within this approach, short and long length scale physics is addressed by the DCA cluster calculation, while intermediate length scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual fermionic Green functions scale as O(1/Lc), so perturbation theory on the dual lattice converges very quickly. Lastly, we study the responses to dynamical modulation of the optical lattice potential by analyzing properties of the repulsive fermionic Hubbard model in an optical lattice. We provide numerical evidence showing the modulations by on-site local interaction cannot be ignored, and instead can even strongly contribute to the dynamical behaviors of the system in highly-doped cases.
We study the Hubbard model using the cellular dynamical mean-field theory (CDMFT) with quantum Monte Carlo (QMC) simulations. We present the algorithmic details of CDMFT with the Hirsch-Fye QMC method for the solution of the self-consistently embedded quantum cluster problem. We use the one- and two-dimensional half filled Hubbard model to gauge the performance of CDMFT+QMC particularly for small clusters by comparing with the exact results and also with other quantum cluster methods. We calculate single-particle Green's functions and self-energies on small clusters to study their size dependence in one and two dimensions. It is shown that in one dimension, CDMFT with two sites in the cluster is already able to describe with high accuracy the evolution of the density as a function of the chemical potential and the compressibility divergence at the Mott transition, in good agreement with the exact Bethe ansatz result. With increasing U the result on small clusters rapidly approaches that of the infinite size cluster. Large scattering rate and a positive slope in the real part of the self-energy in one dimension suggest that the system is a non-Fermi liquid for all the parameters studied here. In two dimensions, at intermediate to strong coupling, even the smallest cluster $({N}_{c}=2\ifmmode\times\else\texttimes\fi{}2)$ accounts for more than 95% of the correlation effect of the infinite-size cluster in the single particle spectrum, suggesting that some of the important problems in strongly correlated electron systems may be studied highly accurately with a reasonable computational effort. Finally, as an application that is sensitive to details of correlations, we show that CDMFT+QMC can describe spin-charge separated Luttinger liquid physics in one dimension. The spinon and holon branches appear only for sufficiently large system sizes.