We introduce regular series expansion for weakly and moderately correlated fermionic systems based on the fluctuating local field approach. The method relies on explicitly accounting for the leading fluctuating mode(s) and is therefore suitable for medium-sized lattices. It can also be used as a solver for the cluster approximations for infinite-size lattices. We introduce a classical fluctuating field coupled to a fermionic collective mode(s). In this way, fluctuations in selected modes are treated in a nonperturbative way. Other degrees of freedom are accounted for in the diagram expansion performed for each value of the fluctuating field. The method is benchmarked for the $U/t=1$ and $U/t=2$ Hubbard lattices at half filling. Results for susceptibility in the antiferromagnetic channel along with the spectral function are compared with the numerically exact quantum Monte Carlo data. Calculations up to the third order of the series expansion are performed and show a uniform convergence to the reference result for the susceptibility. This convergence is observed well below the dynamical mean-field theory (DMFT) Ne\'el temperature and makes our practically simple method applicable in a much wider temperature range than DMFT-based diagrammatic schemes.
We introduce regular series expansion for weakly and moderately correlated fermionic systems based on the fluctuating local field approach. The method relies on explicitly accounting for the leading fluctuating mode(s) and is therefore suitable for medium-sized lattices. It can also be used as a solver for the cluster approximations for infinite-size lattices. We introduce a classical fluctuating field coupled to a fermionic collective mode(s). In this way, fluctuations in selected modes are treated in a nonperturbative way. Other degrees of freedom are accounted for in the diagram expansion performed for each value of the fluctuating field. The method is benchmarked for the U/t=1 and U/t=2 Hubbard lattices at half filling. Results for susceptibility in the antiferromagnetic channel along with the spectral function are compared with the numerically exact quantum Monte Carlo data. Calculations up to the third order of the series expansion are performed and show a uniform convergence to the reference result for the susceptibility. This convergence is observed well below the dynamical mean-field theory (DMFT) Ne'el temperature and makes our practically simple method applicable in a much wider temperature range than DMFT-based diagrammatic schemes.
We establish a way to handle main collective fluctuations in correlated quantum systems based on a fluctuating local field concept. This technique goes beyond standard mean-field approaches, such as Hartree-Fock and dynamical mean-field theories (DMFT), as it includes a fluctuating classical field that acts on the leading order parameter of the system. Effective model parameters of this theory are determined from the variational principle, which allows one to resolve the Fierz ambiguity in decoupling of the local interaction term. In the saddle-point approximation for the fluctuating field our method reproduces the mean-field result. The exact numerical integration over this field allows one to consider nonlinear fluctuations of the global order parameter of the system while local correlations can be accounted for by solving the DMFT impurity problem. We apply our method to the magnetic susceptibility of finite Hubbard systems at half-filling and demonstrate that the introduced technique leads to a superior improvement of results with respect to parental mean-field approaches, without significant numerical complications. We show that the fluctuating local field method can be used in a very broad range of temperatures substantially below the N\'eel temperature of DMFT, which remains a major challenge for all existing theoretical approaches.
Dynamical mean-field theory (DMFT) provides an optimal local approximation for correlated lattice systems by mapping the lattice onto a self-consistent effective impurity model. To account for the missing long-range correlations, we propose a novel extended approach, which we term fluctuating dynamical mean-field theory (fDMFT). It incorporates collective fluctuations of auxiliary impurity models across different sites via functional integration. Technically, this method involves obtaining a family of DMFT solutions on a grid for a self-consistent auxiliary classical field applied to the lattice. While the result can, in principle, be improved diagrammatically, we find that the minimal version of the theory already yields accurate results, with lowest-order diagrammatic corrections offering only minor improvements. This consistent framework, based on our fluctuating local field concept, demonstrates superior performance for the nearly half-filled Hubbard model compared to other known diagrammatic extensions of DMFT.
The contributions from the non-local fluctuations to the dynamical mean field theory (DMFT) were studied using the recently proposed dual fermion approach. Straight forward cluster extensions of DMFT need the solution of a small cluster, where all the short-range correlations are fully taken into account. All the correlations beyond the cluster scope are treated in the mean-field level. In the dual fermion method, only a single impurity problem needs to be solved. Both the short and long-range correlations could be considered on equal footing in this method. The weak-coupling nature of the dual fermion ensures the validity of the finite order diagram expansion. The one and two particle Green's functions calculated from the dual fermion approach agree well with the Quantum Monte Carlo solutions, and the computation time is considerably less than with the latter method. The access of the long-range order allows us to investigate the collective behavior of the electron system, e.g. spin wave excitations.
The main focus of the present thesis lies on the development and application of numerical approaches based on quantum lattice models to describe collective excitations in strongly correlated materials. To this end, an efficient and balanced dual triply irreducible local expansion (D-TRILEX) approach that captures the major effects of spatial correlations in strongly correlated fermionic systems is presented, which is based on partial bosonization of the collective electronic fluctuations. Thereby, a specific choice of the bare interaction in the corresponding channel of instability allows a simultaneous description of collective excitations in different channels without facing the Fierz ambiguity problem which usually arises in methods based on partial bosonization. Compared to more elaborate methods, D-TRILEX does not require the computation of the numerically expensive four-point vertex. Furthermore, taking into account only the leading longitudinal contribution of the bosonic modes to the self-energy reduces the complexity of the diagrammatic structure drastically making the method applicable to more complex realistic multi-orbital systems. As the main result of this thesis, the D-TRILEX approach is introduced and its limit of applicability is studied by applying it to the two-dimensional single-band Hubbard model. Performing a comprehensive diagrammatic analysis by comparing D-TRILEX with methods, which consider different sets of diagrams and different approximations for the four-point vertex, we find that contributions that are not considered by the partially bosonized approximation of the four-point vertex have only a minor effect on electronic degrees of freedom and the longitudinal bosonic modes taking into account in D-TRILEX self-energy are the leading ones in a broad range of control parameters. After investigating the applicability of the D-TRILEX approach we apply the method to study collective electronic effects in the hole-doped InSe monolayers using a realistic electronic-structure model. Thereby, we find that due to the weakly screened long-range Coulomb interaction the system shows a charge density waves phase for the broad range of considered parameters, however, inside this regime we have found a coexisting ferromagnetic phase. D-TRILEX can be considered as an approximation of the more elaborate dual boson (DB) theory, which accounts for local correlations within the dynamical mean-field theory (DMFT) and nonlocal correlations perturbatively. As a further result of this thesis, we investigate the weak coupling regime of the two-dimensional single-band Hubbard model within the ladder DB approximation and compare our results to the numerically exact diagrammatic quantum Monte Carlo (DiagMC) method. We find that the DB method qualitatively captures all the different regimes. Even quantitatively DB shows a good agreement with the exact results for single- and two-particle quantities at high and moderate temperatures. However, the ladder approximation slightly underestimates the strong antiferromagnetic fluctuations that appear in the system at low temperatures. As the final result of this thesis, we perform a parametrization of the local four-point vertex in terms of the (three-point) Hedin vertex and screened interactions in the framework of single-boson exchange (SBE) decomposition. Similar to the D-TRILEX approach the computational costs are reduced by neglecting the interaction-irreducible part of the four-point vertex. The approximated vertex is used to calculate the DMFT susceptibility for different interaction strengths at different temperatures. We find that the approximation shows a good agreement in the limits of weak and strong couplings. Thereby, our study highlights the importance of the Hedin vertex for local vertex corrections.
The dual fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work, we compute the dual fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact diagrammatic determinant Monte Carlo simulations to systematically assess convergence of the dual fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT N\'eel transition, the dual fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4\ensuremath{\le}U/t\ensuremath{\le}12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find, however, that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime.
Collective excitations in fermionic systems play a crucial role in determining their physical properties. An important challenge is to develop efficient theoretical approaches for describing these excitations and their coupling to fermionic degrees of freedom. In this work, we revisit the problem of quantifying the contributions of individual bosonic modes of collective fluctuations to observable properties of correlated fermion systems within the framework of the Fluctuating Local Field (FLF) method. Whereas the auxiliary field in this method was previously considered only classically, we formulate its systematic extension termed Quantum FLF (Q-FLF) that incorporates selected bosonic Matsubara modes, thus tailoring it to description of quantum collective fluctuations. As a testbed, we apply the approach to a half-filled one-dimensional Hubbard chain and compute the Green's function, the total energy, and the antiferromagnetic susceptibility. Our results demonstrate that the proposed scheme enables an efficient and selective characterization of the contributions of individual bosonic modes. In particular, low Matsubara frequencies are found to have a quantitative impact on integrated observables such as total energy and antiferromagnetic susceptibility. At the same time, an accurate description of single-particle properties requires inclusion of higher-frequency bosonic modes.