No Arabic abstract
We discuss the phase diagram of the extended Hubbard model with both attractive and repulsive local and nonlocal interactions. The extended dynamical mean-field theory (EDMFT) and the dual boson method (DB) are compared. The latter contains additional nonlocal correlation effects that are not incorporated in EDMFT. We find that EDMFT and DB give almost identical results in the attractive $V$ regime, where phase separation occurs. This is quite a difference with the previously studied repulsive $V$ regime, where EDMFT and DB give very different phase boundaries for the checkerboard order phase, especially at small $U$.
The extended Hubbard model with an attractive density-density interaction, positive pair hopping, or both, is shown to host topological phases, with a doubly degenerate entanglement spectrum and interacting edge spins. This constitutes a novel instance of topological order which emerges from interactions. When the interaction terms combine in a charge-SU(2) symmetric fashion, a novel partially polarized pseudospin phase appears, in which the topological features of the spin degrees of freedom coexist with long-range $eta$-wave superconductivity. Thus, our system provides an example of an interplay between spontaneous symmetry breaking and symmetry-protected topological order that leads to novel and unexpected properties.
In this work, we adapt the formalism of the dynamical vertex approximation (D$Gamma$A), a diagrammatic approach including many-body correlations beyond the dynamical mean-field theory, to the case of attractive onsite interactions. We start by exploiting the ladder approximation of the D$Gamma$A scheme, in order to derive the corresponding equations for the non-local self-energy and vertex functions of the attractive Hubbard model. Second, we prove the validity of our derivation by showing that the results obtained in the particle-hole symmetric case fully preserve the exact mapping between the attractive and the repulsive models. It will be shown, how this property can be related to the structure of the ladders, which makes our derivation applicable for any approximation scheme based on ladder diagrams. Finally, we apply our D$Gamma$A algorithm to the attractive Hubbard model in three dimensions, for different fillings and interaction values. Specifically, we focus on the parameters region in the proximity of the second-order transition to the superconducting and charge-density wave phases, respectively, and calculate (i) their phase-diagrams, (ii) their critical behavior, as well as (iii) the effects of the strong non-local correlations on the single-particle properties.
A non-perturbative approach to the single-band attractive Hubbard model is presented in the general context of functional derivative approaches to many-body theories. As in previous work on the repulsive model, the first step is based on a local-field type ansatz, on enforcement of the Pauli principle and a number of crucial sum-rules. The Mermin-Wagner theorem in two dimensions is automatically satisfied. At this level, two-particle self-consistency has been achieved. In the second step of the approximation, an improved expression for the self-energy is obtained by using the results of the first step in an exact expression for the self-energy where the high- and low-frequency behaviors appear separately. The result is a cooperon-like formula. The required vertex corrections are included in this self-energy expression, as required by the absence of a Migdal theorem for this problem. Other approaches to the attractive Hubbard model are critically compared. Physical consequences of the present approach and agreement with Monte Carlo simulations are demonstrated in the accompanying paper (following this one).
The two-dimensional attractive Hubbard model is studied in the weak to intermediate coupling regime by employing a non-perturbative approach. It is first shown that this approach is in quantitative agreement with Monte Carlo calculations for both single-particle and two-particle quantities. Both the density of states and the single-particle spectral weight show a pseudogap at the Fermi energy below some characteristic temperature T*, also in good agreement with quantum Monte Carlo calculations. The pseudogap is caused by critical pairing fluctuations in the low-temperature renormalized classical regime $omega < T$ of the two-dimensional system. With increasing temperature the spectral weight fills in the pseudogap instead of closing it and the pseudogap appears earlier in the density of states than in the spectral function. Small temperature changes around T* can modify the spectral weight over frequency scales much larger than temperature. Several qualitative results for the s-wave case should remain true for d-wave superconductors.
We consider a modified extended Hubbard model (EHM) which, in addition to the on-site repulsion U and nearest-neighbor repulsion V, includes polarization effects in second-order perturbation theory. The model is equivalent to an EHM with renormalized U plus a next-nearest-neighbor repulsion term. Using a method based on topological quantum numbers (charge and spin Berry phases), we generalize to finite hopping t the quantum phase diagram in one dimension constructed by van den Brink et al. (Phys. Rev. Lett. 75, 4658 (1995)). At hopping t=0 there are two charge density-wave phases, one spin density-wave phase and one intermediate phase with charge and spin ordering, depending on the parameter values. At t eq 0 the nature of each phase is confirmed by studying correlation functions. However, in addition to the strong-coupling phases, a small region with bond ordering appears. The region occupied by the intermediate phase first increases and then decreases with increasing t, until it finally disappears for t of the order but larger than U. For small t, the topological transitions agree with the results of second order perturbation theory.