Extending many-body numerical techniques which are powerful in the context of simple model calculations to the realm of realistic material simulations can be a challenging task. Realistic systems often involve multiple active orbitals, which increases the complexity and numerical cost because of the large local Hilbert space and the large number of interaction terms or sign-changing off-diagonal Green’s functions. The two-particle self-consistent approach (TPSC) is one such many-body numerical technique, for which multi-orbital extensions have proven to be involved due to the substantially more complex structure of the local interaction tensor. In this paper we extend earlier multi-orbital generalizations of TPSC by setting up two different variants of a fully self-consistent theory for TPSC in multi-orbital systems. We first investigate the strengths and limitations of the approach analytically and then benchmark both variants against dynamical mean-field theory (DMFT) and D-TRILEX results. We find that the exact behavior of the system can be faithfully reproduced in the weak-coupling regime, while at stronger couplings the performance of the two TPSC variants strongly depends on details of the system.
Extending many-body numerical techniques which are powerful in the context of simple model calculations to the realm of realistic material simulations can be a challenging task. Realistic systems often involve multiple active orbitals, which increases the complexity and numerical cost because of the large local Hilbert space and the large number of interaction terms or sign-changing off-diagonal Green’s functions. The two-particle self-consistent approach (TPSC) is one such many-body numerical technique, for which multi-orbital extensions have proven to be involved due to the substantially more complex structure of the local interaction tensor. In this paper we extend earlier multi-orbital generalizations of TPSC by setting up two different variants of a fully self-consistent theory for TPSC in multi-orbital systems. We first investigate the strengths and limitations of the approach analytically and then benchmark both variants against dynamical mean-field theory (DMFT) and D-TRILEX results. We find that the exact behavior of the system can be faithfully reproduced in the weak-coupling regime, while at stronger couplings the performance of the two TPSC variants strongly depends on details of the system.
Abstract One of the most challenging problems in solid state systems is the microscopic analysis of electronic correlations. A paramount minimal model that encodes correlation effects is the Hubbard Hamiltonian, which—regardless of its simplicity—is exactly solvable only in a few limiting cases and approximate many‐body methods are required for its solution. In this review, an overview on the non‐perturbative two‐particle self‐consistent method (TPSC), which was originally introduced to describe the electronic properties of the single‐band Hubbard model, is presented. A detailed derivation of the multi‐orbital generalization of TPSC is introduced here and particular features of the method on exemplary interacting models in comparison to dynamical mean‐field theory results are discussed.
In this paper, we present a multiorbital form of the two-particle self-consistent approach (TPSC) where the effective local and static irreducible interaction vertices are determined by means of the dynamical mean-field theory (DMFT). This approach replaces the approximate ansatz equations for the double-occupations $\ensuremath{\langle}{n}_{\ensuremath{\alpha},\ensuremath{\sigma}}^{}{n}_{\ensuremath{\beta},{\ensuremath{\sigma}}^{\ensuremath{'}}}^{}\ensuremath{\rangle}$ by sampling them directly for the same model using DMFT. Compared to the usual Hartree-Fock-like ansatz, this leads to more accurate local vertices in the weakly correlated regime and provides access to stronger correlated systems that were previously out of reach. This approach is extended by replacing the local component of the TPSC self-energy by the DMFT impurity self-energy, which results in an improved self-energy that incorporates strong local correlations but retains a nontrivial momentum dependence. We find that this combination of TPSC and DMFT provides a significant improvement over the multiorbital formulation of TPSC as it allows to determine the components of the spin vertex without artificial symmetry assumptions and opens the possibility to include the transversal particle-hole channel. The new approach is also able to remove unphysical divergences in the charge vertices in TPSC. We find a general trend that lower temperatures can be accessed in the calculation. Benchmarking single-particle quantities, such as the local spectral function with other many-body methods, we find significant improvement in the more strongly correlated regime.
The Emery model, or three-band Hubbard model, is a Hamiltonian that is thought to contain much of the physics of cuprate superconductors. This model includes two noninteracting p orbitals and one interacting d orbital per unit cell. Few methods that can solve multiorbital interacting Hamiltonians reliably and efficiently exist. Here, we introduce an application of the two-particle self-consistent (TPSC) approach to the Emery model. We construct this method within the framework of the TPSC + DMFT method, which can be seen as a way to introduce nonlocal corrections to dynamical mean-field theory (DMFT). We show that interacting orbital densities, rather than the noninteracting ones, must be used in the calculations. For the Emery model, we find that at constant bare interaction U, the vertex for spin fluctuations, Usp, decreases rapidly with filling. This may be one of the factors that contributes to electron-doped cuprates appearing less correlated than hole-doped ones. More generally, our work opens the road to the application of the TPSC approach to spin fluctuations in multiorbital models.
Nonlocal correlations play an essential role in correlated electron systems, especially in the vicinity of phase transitions and crossovers, where two-particle correlation functions display a distinct momentum dependence. In nonequilibrium settings, the effect of nonlocal correlations on dynamical phase transitions, prethermalization phenomena, and trapping in metastable states is not well understood. In this paper, we introduce a dynamical mean field theory (DMFT) extension to the nonequilibrium two-particle self-consistent (TPSC) approach, which allows to perform nonequilibrium simulations capturing short- and long-ranged nonlocal correlations in the weak-, intermediate-, and strong-correlation regimes. The method self-consistently computes local spin and charge vertices, from which a momentum-dependent self-energy is constructed. Replacing the local part of the self-energy by the DMFT result within this self-consistent scheme provides an improved description of local correlation effects. We explain the details of the formalism and the implementation, and demonstrate the versatility of DMFT+TPSC with interaction quenches and dimensional crossovers in the Hubbard model.
In systems of indistinguishable particles, correlations between the quantum degrees of freedom can give rise to collective modes where particles behave coherently together over some correlation length and time, thereby forming quasiparticles. One example of such quasiparticle would be the π-ton [66], describing the binding of two particle-hole pairs via the exchange of momentum corresponding to the antiferromagnetic (AFM) wave vector. These quasiparticles would emerge in the vicinity of an AFM ordering instability [66, 119, 120] and nonlocal correlations are crucial to describe them. Especially, in systems out of equilibrium, nonlocality can give rise to prethermal states wherein the distribution function can be momentum-dependent [92, 119], and nonthermal exotic physics can be observed and manipulated [40, 68, 79, 84, 93, 132, 147, 157–159]. To be able to model and simulate systems out of equilibrium is a challenge. A few numerical methods such as nonequilibrium Dynamical Mean Field Theory (DMFT) [44] have been used extensively [5]. Nonequilibrium DMFT is designed to treat nonperturbatively the local time-dependent single-particle correlations. Since phenomena such as quantum criticality, high-Tc superconductivity or the creation of π-tons via light-matter coupling stem from nonlocal two-particle correlations, one needs to resort to diagrammatic extensions to DMFT [103] and other methods to effectively describe them, such as the post-processing DMFT treatment of vertex corrections, which will be introduced. The Two-Particle Self-Consistent approach (TPSC) and its variants figure among those methods that include nonlocal two-particle and single-particle correlations in a self-consistent way [121, 122]. TPSC and its variants have the ability to capture the effect of the growing spin fluctuations (paramagnons) in the self-energy when close to amagnetic phase transition leading to the breakdown of the Fermi-liquid quasiparticles at the Fermi level: this would mark the onset of a pseudo-gap [145, 152]. It fulfils the Mermin-Wagner theorem in two dimensions, the Pauli principle, and various two-particle local sum-rules [152]. In the works embodied in this thesis, numerical methods have been designed to access nonequilibrium regimes, by changing model parameters across time, allowing for nonlocal correlations to be properly included. With these new algorithms at hand, the π-tons were studied in the optical conductivity and magnetic response [119, 120] as well as spin and charge fluctuations in the 2D and 3D single-band Hubbard model [121, 122]. It is demonstrated by means of those methods that nonlocal two-particle vertex corrections can underpin prethermal states and could generate transient nonthermal states such as low-energy charge excitations associated with nonthermal negative temperatures. It is also shown that states lying close to the Fermi level control the thermalization timescale and that lattice hopping quenches changing the unit-cell dimensionality could allow one to play around with the charge and spin scattering channels.
In this thesis I present results concerning realistic calculations of correlated fermionic many-body systems. One of the main objectives of this work was the implementation of a hybridization expansion continuous-time quantum Monte Carlo (CT-HYB) algorithm and of a flexible self-consistency loop based on the dynamical mean-field theory (DMFT). DMFT enables us to treat strongly correlated electron systems numerically. After the implementation and extensive testing of the program we investigated different problems to answer open questions concerning correlated systems and their numerical treatment.