No Arabic abstract
We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalisation of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions U>=3t. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.
The nonequilibrium variational-cluster approach is applied to study the real-time dynamics of the double occupancy in the one-dimensional Fermi-Hubbard model after different fast changes of hopping parameters. A simple reference system, consisting of isolated Hubbard dimers, is used to discuss different aspects of the numerical implementation of the approach in the general framework of nonequilibrium self-energy functional theory. Opposed to a direct solution of the Euler equation, its time derivative is found to serve as numerically tractable and stable conditional equation to fix the time-dependent variational parameters.
Using the variational cluster approach (VCA), we study the transition from the antiferromagnetic to the superconducting phase of the two-dimensional Hubbard model at zero temperature. Our calculations are based on a new method to evaluate the VCA grand potential which employs a modified Lanczos algorithm and avoids integrations over the real or imaginary frequency axis. Thereby, very accurate results are possible for cluster sizes not accessible to full diagonalization. This is important for an improved treatment of short-range correlations, including correlations between Cooper pairs in particular. We investigate the cluster-size dependence of the phase-separation tendency that has been proposed recently on the basis of calculations for smaller clusters. It is shown that the energy barrier driving the phase separation decreases with increasing cluster size. This supports the conjecture that the ground state exhibits microscopic inhomogeneities rather than macroscopic phase separation. The evolution of the single-particle spectum as a function of doping is studied in addtion and the relevance of our results for experimental findings is pointed out.
We study the correlation-induced deformation of Fermi surfaces by means of a new diagrammatic method which allows for the analytical evaluation of Gutzwiller wave functions in finite dimensions. In agreement with renormalization-group results we find Pomeranchuk instabilities in two-dimensional Hubbard models for sufficiently large Coulomb interactions.
We applied the Recurrent Variational Approach to the two-leg Hubbard ladder. At half-filling, our variational Ansatz was a generalization of the resonating valence bond state. At finite doping, hole pairs were allowed to move in the resonating valence bond background. The results obtained by the Recurrent Variational Approach were compared with results from Density Matrix Renormalization Group.
The extended Hubbard model in the atomic limit (AL-EHM) on a square lattice with periodic boundary conditions is studied with use of the Monte Carlo (MC) method. Within the grand canonical ensemble the phase and order-order boundaries for charge orderings are obtained. The phase diagrams include three types of charge ordered phases and the nonordered phase. The system exhibits very rich structure and shows unusual multicritical behavior. In the limiting case of tij = 0, the EHM is equivalent to the pseudospin model with single-ion anisotropy 1/2U, exchange interaction W in an effective magnetic field (mu-1/2U-zW). This classical spin model is analyzed using the MC method for the canonical ensemble. The phase diagram is compared with the known results for the Blume-Capel model.