The helicity modulus, which is the stiffness associated with a twisted order parameter, for the two-dimensional Hubbard model is calculated for the equivalent cases of (i) attractive on-site interaction (negative U) with arbitrary strength, arbitrary electron density, and zero magnetic field and (ii) repulsive on-site interaction (positive U) with arbitrary strength, at half-filling and in an arbitrary magnetic field. An explicit formula for the helicity modulus is derived using the Bogoliubov-Hartree-Fock approximation. An improved value for the helicity modulus is obtained by performing variational Monte Carlo calculations using a Gutzwiller projected trial wave function. To within a small correction term the helicity modulus is found to be given by -1/8 the average kinetic energy. The variational Monte Carlo calculation is found to increase the value of the helicity modulus by a small amount (about 5% for intermediate values of the interaction strength) compared to the results from the Bogoliubov-Hartree-Fock approximation. In the case of attractive interaction, from a comparison with the Kosterlitz-Thouless relation between critical temperature and helicity modulus, the critical temperature for a Kosterlitz-Thouless transition is calculated and a phase diagram is obtained. An optimal critical temperature is found for an intermediate value of U. We discuss connections of our results with results in the literature on the Hubbard model using the random-phase approximation and quantum Monte Carlo calculations.
We study the Hubbard model on the triangular lattice for general fillings by means of the slave-boson technique. In particular, we consider the stability regions of the paramagnetic, antiferromagnetic and linearly polarized spin-density-wave phases that are present at half-filling. The model presents a large anisotropy under particle or hole doping of the half-filled band, a consequence of the non-bipartite nature of the lattice. We also compute the helicity modulus in the antiferromagnetic phase, which characterizes the stiffness of the magnetic order.
We present a comprehensive study of the two-dimensional one-band Hubbard model applying the spin-rotation-invariant slave-boson method. We utilize a spiral magnetic mean field and fluctuations around a paramagnetic mean field to determine the magnetic phase diagram and find the two approaches to be in good agreement. Apart from the commensurate phases characterized by ordering wave vectors $\mathbf{Q}=(\ensuremath{\pi},\ensuremath{\pi})$, $(0,\ensuremath{\pi})$, and (0,0) we find incommensurate phases where the ordering wave vectors
Q=(Q,Q) and $(Q,\ensuremath{\pi})$ vary continuously with filling, interaction strength, or temperature. The mean-field quantities magnetization and effective mass are found to change discontinuously at the phase boundaries separating the
(Q,Q) and $(Q,\ensuremath{\pi})$ phases, indicating a first-order transition. The band structure and Fermi surface is shown in selected cases. The dynamic spin and charge susceptibilities as well as the structure factors are calculated and discussed, including the emergence of collective modes of the zero sound and Mott insulator type. The dynamical conductivity is calculated in dependence of doping, interaction strength, and temperature. Finally, a temperature-interaction strength phase diagram is established.
The velocity of long-wavelength spin waves for the two-dimensional Hubbard model on a square lattice at half-filling is calculated using the variational Monte Carlo method with a Gutzwiller-type wave function. We show that the result from the Random Phase Approximation (RPA) for the spin wave velocity can also be found as the square root of the ratio of the helicity modulus and the perpendicular susceptibility as calculated in the Hartree-Fock approximation (HFA). Since the Gutzwiller wave function improves upon the HFA, our result for the spin wave velocity improves upon the RPA result. For the on-site repulsion U between 3 and 12 times the hopping integral t the spin wave velocity increases by 34 to 15% compared to the RPA result.
The cluster slave-spin method is employed to investigate systematically the ground state properties of the Hubbard model on a square lattice with doping $\ensuremath{\delta}$ and coupling strength
U being its parameters. At half-filling, a relation between the staggered magnetization
M and the antiferromagnetic (AFM) gap ${\mathrm{\ensuremath{\Delta}}}_{\text{AFM}}$ is established in the small
U limit to compare with that from the Hartree-Fock theory, and a first-order metal-insulator Mott transition in the paramagnetic state is substantiated, which is characterized by discontinuities and hystereses at
UMott=10t. The interaction
Uc for the crossover in the AFM state, separating the weak- and strong-coupling regimes, is found to remain almost unchanged with large dopings, and smaller than
UMott at half-filling because of long range AFM correlations. Finally, an overall phase diagram in the $U\ensuremath{-}\ensuremath{\delta}$ plane is presented, which is composed of four regimes: the AFM insulator at half-filling, the AFM metal with the compressibility $\ensuremath{\kappa}>0$ or $\ensuremath{\kappa}<0$, and the paramagnetic metal, as well as three phase transitions: (i) From the AFM metal to the paramagnetic metal, (ii) between the AFM metal phases with positive and negative $\ensuremath{\kappa}$, and (iii) separating the AFM insulating phase at $\ensuremath{\delta}=0$ from the AFM metal phase for $\ensuremath{\delta}>0$.
Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.