Browsing by Author "Ferrero, Michel"

Now showing 1 - 3 of 3
  • No Thumbnail Available
    Item
    Analytical solution for time integrals in diagrammatic expansions: Application to real-frequency diagrammatic Monte Carlo
    Vučičević, Jakša; Stipsić, Pavle; Ferrero, Michel
    Recent years have seen a revived interest in the diagrammatic Monte Carlo (DiagMC) methods for interacting fermions on a lattice. A promising recent development allows one to now circumvent the analytical continuation of dynamic observables in DiagMC calculations within the Matsubara formalism. This is made possible by symbolic algebra algorithms, which can be used to analytically solve the internal Matsubara frequency summations of Feynman diagrams. In this paper, we take a different approach and show that it yields improved results. We present a closed-form analytical solution of imaginary-time integrals that appear in the time-domain formulation of Feynman diagrams. We implement and test a DiagMC algorithm based on this analytical solution and show that it has numerous significant advantages. Most importantly, the algorithm is general enough for any kind of single-time correlation function series, involving any single-particle vertex insertions. Therefore, it readily allows for the use of action-shifted schemes, aimed at improving the convergence properties of the series. By performing a frequency-resolved action-shift tuning, we are able to further improve the method and converge the self-energy in a nontrivial regime, with only 3-4 perturbation orders. Finally, we identify time integrals of the same general form in many commonly used Monte Carlo algorithms and therefore expect a broader usage of our analytical solution.
  • No Thumbnail Available
    Item
    Practical consequences of the Luttinger-Ward functional multivaluedness for cluster DMFT methods
    Vučičević, Jakša; Wentzell, Nils; Ferrero, Michel; Parcollet, Olivier
    The Luttinger-Ward functional (LWF) has been a starting point for conserving approximations in many-body physics for 50 years. The recent discoveries of its multivaluedness and the associated divergence of the two-particle irreducible vertex function Γ have revealed an inherent limitation of this approach. Here we demonstrate how these undesirable properties of the LWF can lead to a failure of computational methods based on an approximation of the LWF. We apply the nested cluster scheme (NCS) to the Hubbard model and observe the existence of an additional stationary point of the self-consistent equations, associated with an unphysical branch of the LWF. In the strongly correlated regime, starting with the first divergence of Γ, this unphysical stationary point becomes attractive in the standard iterative technique used to solve DMFT. This leads to an incorrect solution, even in the large cluster size limit, for which we discuss diagnostics.
  • No Thumbnail Available
    Item
    Real-frequency diagrammatic Monte Carlo at finite temperature
    Vučičević, Jakša; Ferrero, Michel
    Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculation ultimately requires the ill-defined analytical continuation from the imaginary-to the real-frequency domain. It was recently proposed [A. Taheridehkordi, Phys. Rev. B 99, 035120 (2019)2469-995010.1103/PhysRevB.99.035120] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy ς(ω) of the doped 32×32 cyclic square-lattice Hubbard model in a nontrivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.