The dynamical vertex approximation and related methods represent a great leap forward in the description of nonlocal electronic correlations. It is based on a local two-particle irreducible vertex, from which the momentum-dependent self-energy is constructed via Bethe-Salpeter and Schwinger-Dyson equations. In this paper, the hitherto used one-shot calculations are improved by repeating these equations with feedback from the updated self-energy until self-consistency. The validity of the method is benchmarked by several numerical examples and, prominently, the formation of a pseudogap is demonstrated.
The dynamical vertex approximation and related methods represent a great leap forward in the description of nonlocal electronic correlations. It is based on a local two-particle irreducible vertex, from which the momentum-dependent self-energy is constructed via Bethe-Salpeter and Schwinger-Dyson equations. In this paper, the hitherto used one-shot calculations are improved by repeating these equations with feedback from the updated self-energy until self-consistency. The validity of the method is benchmarked by several numerical examples and, prominently, the formation of a pseudogap is demonstrated.
We suggest a new approach for treatment of local and non-local interactions in correlated electronic systems, which is based on the combination of the (extended) dynamical mean-field theory ((E)DMFT) and the two-particle irreducible functional renormaliztion-group (2PI-fRG) method. The considering approach uses self-energy and the two-particle irreducible vertices, obtained from (E)DMFT, as an input of 2PI-fRG approach. Using 2PI vertices may improve the applicability of various truncations and fulfilling conservation laws in comparison with one-particle irreducible approaches. In case of purely local interaction in a certain truncation of the DMFT+2PI-fRG equations, the obtained equation for the self-energy has a similar, although not identical, structure to that in the ladder dynamic vertex approximation (D$\Gamma $A). For the non-local interactions, in a simplest truncation we reproduce the results for the two-particle vertices/susceptibilities in the ladder approximation of the dual boson approach, but obtain different equation for the self-energy with a correct treatment of the one-particel reducible vertices of higher orders. The proposed scheme is rather general and can be applied to study various phenomena in strongly-correlated electronic systems, e.g. as a tool describing ab initio screening of the Coulomb interaction in strongly correlated systems.
A two-particle self-consistency is rarely part of mean-field theories. It is, however, essential for avoiding spurious critical transitions and unphysical behavior. We present a general scheme for constructing analytically controllable approximations with self-consistent equations for the two-particle vertices based on the parquet equations. We explain in detail how to reduce the full set of parquet equations so as not to miss quantum criticality in strong coupling. We further introduce a decoupling of convolutions of the dynamical variables in the Bethe-Salpeter equations to make them analytically solvable. We connect the self-energy with the two-particle vertices to satisfy the Ward identity and the Schwinger-Dyson equation and discuss the role of the one-particle self-consistency in making the approximations reliable in the whole spectrum of the input parameters. Finally, we demonstrate the general construction on the simplest static approximation that we apply to the Kondo behavior of the single-impurity Anderson model.
In a many-body perturbation treatment of electronic excitations, one has to solve the so-called Bethe–Salpeter equation (BSE) for the kernel (vertex) describing the interaction between electrons. In general, the BSE exhibits a non-local and frequency-dependent screened interaction and consequently it is extremely difficult to solve. We have developed a scheme that solves the BSE equation iteratively, including dynamically screening, by considering the change in the electron self-energy upon applying a time-dependent field. The BSE is obtained from the self-energy in the GW approximation (GWA) using the Schwinger functional technique. We apply our formalism to a simple model system and discuss briefly changes in the charge response, when dynamical screening is taken into account.
We propose an approach describing correlated electronic systems within the two-particle irreducible functional renormalization-group (2PI-fRG) approach, which uses initial self-energy and the two-particle irreducible vertices, obtained from (extended) dynamical mean-field theory [E(DMFT)]. Using the 2PI-fRG approach allows us to treat both local and nonlocal interactions. In the case of purely local interaction the corresponding equations have similar (although not identical) structure to the earlier developed DMFT+fRG (DMF2RG) approach. In a simplest truncation, neglecting scale dependence of the two-particle irreducible vertices, we reproduce the results for the two-particle vertices/susceptibilities in the ladder approximation of the dual-boson (DB) or the dynamic vertex approximation approach; in more sophisticated truncations the method allows us to consider nonlocal corrections to the self-energy as well as the interplay of charge and spin correlations. The proposed scheme is tested on the two-dimensional standard and extended $U\text{\ensuremath{-}}V$ half-filled Hubbard models. For the standard Hubbard model we obtain nonlocal self-energy, which is in agreement with numerical studies; for the extended Hubbard model we find the boundary of charge instability, which agrees well with the results of the DB approach. We also ascertain that the effect of spin correlations on electron interaction in the charge channel, not considered previously in the DB approach, only slightly reduces critical next-nearest-neighbor interaction of charge instability of the extended Hubbard model at the considered finite low temperature, yielding better agreement with dynamic cluster approximation. The considered method is rather general and can be applied to study various phenomena in strongly correlated electronic systems.
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.
We present a formalism for strongly correlated electron systems which consists in a local approximation of the dynamical three-leg interaction vertex. This vertex is self-consistently computed using a quantum impurity model with dynamical interactions in the charge and spin channels, similar to dynamical mean field theory approaches. The electronic self-energy and the polarization are both frequency and momentum dependent. The method interpolates between the spin-fluctuation or GW approximations at weak coupling and the atomic limit at strong coupling. We apply the formalism to the Hubbard model on a two-dimensional square lattice and show that as interactions are increased towards the Mott insulating state, the local vertex acquires a strong frequency dependence, driving the system to a Mott transition, while at low enough temperatures the momentum dependence of the self-energy is enhanced due to large spin fluctuations. Upon doping, we find a Fermi arc in the one-particle spectral function, which is one signature of the pseudogap state.