Abstract We propose a systematic approach to the systems of correlated electrons, the so-called -DE-GWF method, based on reciprocal-space ( -resolved) diagrammatic expansion of the variational Gutzwiller-type wave function for parametrized models of correlated fermions. The present approach, in contrast to either variational Monte-Carlo (VMC), or the recently developed real-space diagrammatic expansion of the Gutzwiller-type wave function (direct-space DE-GWF technique), is applicable directly in the thermodynamic limit and thus is suitable for describing selected singular features of the wave-vector-dependent quantities. We employ the -DE-GWF method to extract the non-analytic part of the two leading moments of the fermion spectral-density function across the (two-dimensional) Brillouin zone for the Hubbard model and away from the half-filling. Those moments are used to evaluate the nodal quasiparticle velocities and their spectral weights in the correlated superconducting state. The two velocities determined in that manner exhibit scaling with the electron concentration qualitatively different from that obtained earlier for the excited states of the high- T c cuprates within the projected quasi-particle ansatz, and the results are in a very good quantitative agreement with experimental data if interpreted as those characterizing the spectrum below and above the observed kink. We provide a detailed discussion of the two gaps and two excitation branches (two velocities) appearing naturally within our DE-GWF approach. The two separate sets of characteristics distinguish the renormalized quasiparticle states very close to the Fermi surface from the deeper correlated-state properties. Also, an enhancement of the -dependent magnetic susceptibility is shown to contain a spin-fluctuation contribution within our language. Finally, the -DE-GWF approach is compared to both the VMC and real-space DE-GWF results for the cases of Hubbard and t - J - U models.
Abstract We propose a systematic approach to the systems of correlated electrons, the so-called -DE-GWF method, based on reciprocal-space ( -resolved) diagrammatic expansion of the variational Gutzwiller-type wave function for parametrized models of correlated fermions. The present approach, in contrast to either variational Monte-Carlo (VMC), or the recently developed real-space diagrammatic expansion of the Gutzwiller-type wave function (direct-space DE-GWF technique), is applicable directly in the thermodynamic limit and thus is suitable for describing selected singular features of the wave-vector-dependent quantities. We employ the -DE-GWF method to extract the non-analytic part of the two leading moments of the fermion spectral-density function across the (two-dimensional) Brillouin zone for the Hubbard model and away from the half-filling. Those moments are used to evaluate the nodal quasiparticle velocities and their spectral weights in the correlated superconducting state. The two velocities determined in that manner exhibit scaling with the electron concentration qualitatively different from that obtained earlier for the excited states of the high- T c cuprates within the projected quasi-particle ansatz, and the results are in a very good quantitative agreement with experimental data if interpreted as those characterizing the spectrum below and above the observed kink. We provide a detailed discussion of the two gaps and two excitation branches (two velocities) appearing naturally within our DE-GWF approach. The two separate sets of characteristics distinguish the renormalized quasiparticle states very close to the Fermi surface from the deeper correlated-state properties. Also, an enhancement of the -dependent magnetic susceptibility is shown to contain a spin-fluctuation contribution within our language. Finally, the -DE-GWF approach is compared to both the VMC and real-space DE-GWF results for the cases of Hubbard and t - J - U models.
A systematic diagrammatic expansion for Gutzwiller wavefunctions (DE-GWFs) proposed very recently is used for the description of the superconducting (SC) ground state in the two-dimensional square-lattice t – J model with the hopping electron amplitudes t (and ) between nearest (and next-nearest) neighbors. For the example of the SC state analysis we provide a detailed comparison of the methodʼs results with those of other approaches. Namely, (i) the truncated DE-GWF method reproduces the variational Monte Carlo (VMC) results and (ii) in the lowest (zeroth) order of the expansion the method can reproduce the analytical results of the standard Gutzwiller approximation (GA), as well as of the recently proposed 'grand-canonical Gutzwiller approximation' (called either GCGA or SGA). We obtain important features of the SC state. First, the SC gap at the Fermi surface resembles a wave only for optimally and overdoped systems, being diminished in the antinodal regions for the underdoped case in a qualitative agreement with experiment. Corrections to the gap structure are shown to arise from the longer range of the real-space pairing. Second, the nodal Fermi velocity is almost constant as a function of doping and agrees semi-quantitatively with experimental results. Third, we compare the doping dependence of the gap magnitude with experimental data. Fourth, we analyze the k -space properties of the model: Fermi surface topology and effective dispersion. The DE-GWF method opens up new perspectives for studying strongly correlated systems, as it (i) works in the thermodynamic limit, (ii) is comparable in accuracy to VMC, and (iii) has numerical complexity comparable to that of the GA (i.e., it provides the results much faster than the VMC approach).
The principal purpose of this topical review is to single out some of the universal features of high-temperature (high-Tc) and related strongly-correlated systems, which can be compared with experiment in a quantitative manner. The description starts with the concept of exchange-interaction-mediated (real-space) pairing, combined with strong correlations among narrow band electrons, for which the reference state is that of the Mott–Hubbard insulator. The physical discussion of concrete properties relies on variational approach, starting from generalized renormalized mean-field theory (RMFT) in the form of statistically-consistent Gutzwiller approximation (SGA), and its subsequent generalization, i.e., a systematic Diagrammatic Expansion of the Variational (Gutzwiller-type) Wave Function (DE-GWF). The solution leads to the two energy scales, one involving unconventional quasiparticles close to the Fermi energy, and the other reflecting the fully correlated state, involving electrons deeper below the Fermi surface. Those two regimes are separated by a kink in the dispersion relation, which is observed in photoemission. As a result, one obtains both the doping dependent properties in the correlated state, and effective renormalized quasiparticles. The reviewed ground-state characteristics for high-Tc systems encompass high-Tc superconductivity, nematicity, charge- (and pair-) density-wave effects, as well as non-BCS kinetic energy gain in the paired state. The calculated dynamic properties are: the universal Fermi velocity, the Fermi wave-vector, the effective mass enhancement, the pseudogap; and the d-wave gap magnitude. We discuss that, within the variational approach, the minimal realistic model is represented by the so-called t-J-U Hamiltonian. Inadequacy of the t-J and Hubbard models, particularly of their RMFT versions, is discussed explicitly. For heavy fermion systems, modeled by the Anderson lattice model and with the DE-GWF approach, we discuss the phase diagram, encompassing superconducting, Kondo insulating, ferro- and anti-ferromagnetic states. The superconducting state is then a two-d-wave gap system. If the orbital degeneracy of f-electrons is included, the coexistent ferromagnetic- (spin–triplet) superconducting phases appear and match those observed for UGe2 in a semiquantitative manner. Finally, in the second part, we generalize our approach to the collective spin and charge fluctuations in high-Tc systems, starting from variational approach, defining the saddle-point state, and combining it with 1/N expansion. The present scheme differs essentially from that starting from the saddle-point Hartree–Fock approximation, and incorporating the fluctuations in the random phase approximation (RPA). The spectrum of collective spin and charge excitations is determined for the Hubbard and t-J-U models, and subsequently compared quantitatively with recent experiments. The Appendices provide formal details to make this review self-contained.
A systematic diagrammatic expansion for Gutzwiller-wave functions (DE-GWF) is formulated and used for the description of superconducting (SC) ground state in the two-dimensional Hubbard model with electron-transfer amplitudes t (and ${t}^{\ensuremath{'}}$) between nearest (and next-nearest) neighbors. The method is numerically very efficient and allows for a detailed analysis of the phase diagram as a function of all relevant parameters (U, $\ensuremath{\delta}$, ${t}^{\ensuremath{'}}$) and a determination of the kinetic-energy driven pairing region. SC states appear only for substantial interactions, $U/t\ensuremath{\gtrsim}3$, and for not too large hole doping, $\ensuremath{\delta}\ensuremath{\lesssim}0.32$ for ${t}^{\ensuremath{'}}=0.25t$; this upper critical doping value agrees well with experiment for the cuprate high-temperature superconductors. We also obtain other important features of the SC state: (i) the SC gap at the Fermi surface resembles ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave only around the optimal doping and the corrections to this state are shown to arise from the longer range of the pairing; (ii) the nodal Fermi velocity is almost constant as a function of doping and agrees quantitatively with the experimental results; (iii) the SC transition is driven by the kinetic-energy lowering for low doping and strong interactions.
We study the superconducting phase of the Hubbard model using the Gutzwiller variational wave function (GWF) and the recently proposed diagrammatic expansion technique (DE-GWF). The DE-GWF method works on the level of the full GWF and in the thermodynamic limit. Here, we consider a finite-size system to study the accuracy of the results as a function of the system size (which is practically unrestricted). We show that the finite-size scaling used, e.g. in the variational Monte Carlo method can lead to significant, uncontrolled errors. The presented research is the first step towards applying the DE-GWF method in studies of inhomogeneous situations, including systems with impurities, defects, inhomogeneous phases, or disorder.
High-temperature superconductivity in copper oxides (cuprates, e.g. La2−xSrxCuO4, YBa2Cu3O6+δ or Bi2Sr2CaCu2O8+δ) remains among the most spectacular phenomena in condensed matter physics. Since its discovery in 1986, an enormous number (∼ 105) of papers on the subject have appeared. So far, there is no single, commonly accepted theory of high-temperature superconductivity. However, it is widely believed that a basic description of this phenomenon can be provided by a single-band Hubbard model or its derivative, the t-J model. The latter model is regarded as a minimal microscopic model, capable of describing the essential aspects of the complex physics of the cuprates. Unfortunately, in general case, neither of those two models can be solved exactly and therefore various approximate methods are used. Among them, the so-called mean-field methods provide a simple, yet fairly reasonable description of the cuprates. In particular, some of the main qualitative features of the phase diagram and the essential features of electronic spectrum are roughly reproduced. A standard mean-field approach to the t-J model, known under the name of renormalized mean-field theory (RMFT) goes beyond the Hartree-Fock approximation. Consequently, its fully consistent treatment requires a novel theoretical approach. This has been our original motivation to develop a general approach to the mean-field models, which is based on the maximum entropy (MaxEnt) principle. The method is presented in detail in Part II of this Thesis, and in Part III it is applied to study RMFT of the t-J model. First, we compare the results obtained within our formalism with those of the frequently used non-variational approach based entirely on the self-consistent equations. Also, various versions of RMFT are compared, and the most satisfactory of them is selected. This optimal version is subsequently used to study different versions of the original t-J Hamiltonian. As a result, upper critical concentration and doping dependence of the selected physical quantities (e.g. the superconducting gap and the Fermi velocity) is determined at low temperatures and in the absence of external magnetic field. We compare our findings both with theoretical results obtained from the Variational Monte Carlo (VMC) methods, as well as with the experimental data for selected cuprates. We show that the version of RMFT approach formulated in this Thesis provides a reasonable qualitative and in some cases semiquantitative rationalization of the principal characteristics of the hole-doped high-temperature superconductors at the optimal doping and in the overdoped regime. Possible extensions of the proposed analysis are mentioned at the end.
A nonzero-temperature technique for strongly correlated electron lattice systems, combining elements of both variational wave function (VWF) approach and expansion in the inverse number of fermionic flavors (1/Nf), is developed. The departure point, VWF method, goes beyond the renormalized mean-field theory and provides semiquantitative description of principal equilibrium properties of high-Tc superconducting cuprates. The developed here scheme of VWF+1/Nf, in the leading order provides dynamical spin and charge responses around the VWF solution, generalizing the weak-coupling spin-fluctuation theory to the regime of strong correlations. Thermodynamic corrections to the correlated saddle-point state arise systematically at consecutive orders. Explicitly, VWF+1/Nf is applied to evaluate dynamical response functions for the hole-doped Hubbard model and compared with available determinant quantum Monte Carlo data, yielding a good overall agreement in the regime of coherent collective-mode dynamics. The emergence of well-defined spin and charge excitations from the incoherent continua is explicitly demonstrated and a nonmonotonic dependence of the charge-excitation energy on the interaction magnitude is found. The charge mode energy saturates slowly when approaching the strong-coupling limit, which calls for a reevaluation of the $t\text{\ensuremath{-}}J$-model approach to the charge dynamics in favor of more general $t\text{\ensuremath{-}}J\text{\ensuremath{-}}U$ and $t\text{\ensuremath{-}}J\text{\ensuremath{-}}U\text{\ensuremath{-}}V$ models. The results are also related to recent inelastic resonant x-ray and neutron-scattering experiments for the high-Tc cuprates.