We discuss the calculation of the double occupancy using dynamical mean-field theory in finite dimensions. The double occupancy can be determined from the susceptibility of the auxiliary impurity model or from the lattice susceptibility. The former method typically overestimates, whereas the latter underestimates the double occupancy. We illustrate this for the square-lattice Hubbard model. We propose an approach for which both methods lead to identical results by construction and which resolves this ambiguity. This self-consistent dual boson scheme results in a double occupancy that is numerically close to benchmarks available in the literature.
We discuss the calculation of the double occupancy using dynamical mean-field theory in finite dimensions. The double occupancy can be determined from the susceptibility of the auxiliary impurity model or from the lattice susceptibility. The former method typically overestimates, whereas the latter underestimates the double occupancy. We illustrate this for the square-lattice Hubbard model. We propose an approach for which both methods lead to identical results by construction and which resolves this ambiguity. This self-consistent dual boson scheme results in a double occupancy that is numerically close to benchmarks available in the literature.
The dual fermion (DF) method allows for calculating corrections due to nonlocal correlations relative to an effective impurity model. Choosing the impurity as that of a dynamical mean-field theory (DMFT) solution at self-consistency is popular, and the corrections from dual fermion theory are physically meaningful. We investigate the effect of choosing the impurity instead in a self-consistent manner and find for the two-dimensional Hubbard model an exponential increase of the correlation length and susceptibility at low temperatures. There are pronounced differences for the two self-consistency schemes that are discussed in the literature; the self-consistent DF solution can even be more metallic than the DMFT solution.
In this thesis, we use the Dual Fermion approximation (DF) to study the Hubbard model on the 2D square lattice by using two-particle quantities such as susceptibilities and full vertex functions. In the first chapter we give an introduction to the strongly correlated systems including their definitions and features. Then, in the next chapter, to solve our model we present some numerical methods. The numerical methods that we used for our study are the dynamical mean-field theory (DMFT), the continuous-time auxiliary field impurity solver (CTAUX), Fluctuation diagnostics and Dual Fermion approximation (DF). The DF method provides high-resolution results for single-particle quantities such as the density of states (DOS) and Self Energy as well two particles spin susceptibility and vertex functions. With these results, we can present high-resolution results of the fluctuation diagnostics method with minimal computational expense. We examine the full Q-vector dependence of the extended self-energy at the Nodal and Antinodal k-points and provide interpretation for the claim that spin fluctuations are responsible for the observed metal-to-insulator behaviour.
Since the first investigation of the Hubbard model in the limit of infinite dimensions by Metzner and Vollhardt, dynamical mean-field theory (DMFT) has become a very powerful tool for the investigation of lattice models of correlated electrons. In DMFT the lattice model is mapped on an effective quantum impurity model in a bath which has to be determined self-consistently. This approach led to significant progress in our understanding of typical correlation problems such as the Mott transition; furthermore, the combination of DMFT with ab-initio methods now allows for a realistic treatment of correlated materials. The focus of these lecture notes is on the relation between quantum impurity physics and the physics of lattice models within DMFT. Issues such as the observability of impurity quantum phase transitions in the corresponding lattice models are discussed in detail.
Extensions of dynamical mean-field theory (DMFT) make use of quantum impurity models as nonperturbative and exactly solvable reference systems which are essential to treat the strong electronic correlations. Through the introduction of retarded interactions on the impurity, these approximations can be made two-particle self-consistent. This is of interest for the Hubbard model because it allows to suppress the antiferromagnetic phase transition in two dimensions in accordance with the Mermin-Wagner theorem, and to include the effects of bosonic fluctuations. For a physically sound description of the latter, the approximation should be conserving. In this paper, we show that the mutual requirements of two-particle self-consistency and conservation lead to fundamental problems. For an approximation that is two-particle self-consistent in the charge and longitudinal spin channels, the double occupancy of the lattice and the impurity is no longer consistent when computed from single-particle properties. For the case of self-consistency in the charge and longitudinal as well as transversal spin channels, these requirements are even mutually exclusive so that no conserving approximation can exist. We illustrate these findings for a two-particle self-consistent and conserving DMFT approximation.
Disorder is ubiquitous in real materials and influences the physical properties like the conductivity to varying degrees. \nIf electron-electron interactions are strong, theoretical and numerical treatment of these systems becomes challenging. In this thesis a numerical approach is developed to address these systems, \ntreating both interactions and disorder on equal footing. \nThe approach is based on the dual fermion approach for interacting systems developed by Rubtsov et al. \nTerletska et al. applied the ideas of the dual fermion approach to disordered non-interacting systems. In this approach, the replica trick is used to integrate out the disorder in favor of an effective electron-electron interaction. \nWe extended the approach from Terletska et al. to treat disordered interacting systems. \n \nDual Fermions allow to take into account non-local fluctuations by means of a perturbative expansion around an impurity problem. The impurity reference system is determined self-consistently, analogously to the dynamical mean-field \ntheory. \nThe perturbative expansion is expected to yield good results for small and large values of interaction strength and disorder. A priori, it is not clear what to expect for intermediate values, \nbut experience shows that oftentimes good results are obtained for this region. \n \nAn advantage of the dual fermion approach is that there is no sign-problem for a single orbital model if quantum Monte Carlo is used to solve the interacting reference system. \nAdditionally, perturbation theory is usually numerically much cheaper than fully solving an interacting lattice or cluster problem. \nThus, the dual fermion approach allows to address regions of parameter space that are not accessible to lattice quantum Monte Carlo calculations or cluster extension of dynamical mean-field theory. \nCluster extensions of the dynamical mean-field theory are for example the dynamical cluster approximation or the cellular dynamical mean-field theory. \n \nThe new approach is benchmarked with respect to the Anderson-Falicov-Kimball model. The Anderson-Hubbard model is the main application. \nFor both models, the algorithm is tested in one dimension for the single particle Green function and compared to the dynamical cluster approximation to estimate the quality of the results. \nThe Anderson-Falicov-Kimball model is studied further in two dimensions and the phase diagram on the disorder-interaction plane is calculated from both the single particle Green function and the dc conductivity. \n \nFor the Anderson-Hubbard model, the antiferromagnetic phase transition, the Mott transition and the Anderson transition are investigated in three dimensions. The dynamical mean-field theory results are \nqualitatively confirmed for the most part. Quantitative corrections are in agreement with dynamical cluster approximation results, thus confirming the validity of the results. \nWe found a qualitative difference between the dual fermion approach and dynamical mean-field theory for the temperature dependence of the hysteresis of the double occupancy. The dual fermion result \nshows a decreasing lower critical interaction strength for decreasing temperature, which is in contrast with the dynamical mean-field result that shows an increase of the lower critical interaction strength. \nThis behavior happens for small temperatures that are not accessible using for example the dynamical cluster approximation, thus highlighting the usefulness of the dual fermion approach. \nWe found this behavior for all values of the disorder that we considered.
Numerical results for ground-state and excited-state properties (energies, double occupancies, and Matsubara-axis self-energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary-field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed-node approximation, unrestricted coupled cluster theory, and multireference projected Hartree-Fock methods. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic-limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.