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.
The critical temperature for superconductivity for the two-dimensional attractive (negative- U ) Hubbard model on a square lattice is calculated as a function of the strength of the interaction and of the density of electrons. The critical temperature is obtained by calculating the helicity modulus associated with a wavelike distortion of the gap parameter in the BCS approximation as a function of temperature and subsequent comparison with the Kosterlitz-Thouless relation between critical temperature and helicity modulus. The adequacy of the approximations is discussed and a comparison with recent Monte Carlo calculations is made.
The slave-boson mean-field method is used to study the two-dimensional Hubbard model. A magnetic phase diagram allowing for paramagnetism, weak and strong ferromagnetism and antiferromagnetism is constructed and compared to the corresponding phase diagram using the Hartree-Fock approximation (HFA). Magnetically ordered regions are reduced by a factor of about three along both the t/U and density axes compared to the HFA. Using the spin-rotation-invariant formulation of the slave-boson method the helicity modulus is computed and for half filling is found practically to coincide with that found using variational Monte Carlo calculations with the Gutzwiller wave function. Off half filling the results can be used to compare with quantum Monte Carlo calculations of the effective hopping parameter. Contrary to the case of half filling, the slave-boson approach is seen to greatly improve the results of the HFA when off half filling.
We have considered the half-filled disordered attractive Hubbard model on a square lattice, in which the on-site attraction is switched off on a fraction
f of sites, while keeping a finite
U on the remaining ones. Through quantum Monte Carlo simulations for several values of
f and
U and for system sizes ranging from $8\ifmmode\times\else\texttimes\fi{}8$ to $16\ifmmode\times\else\texttimes\fi{}16$, we have calculated the configurational averages of the equal-time pair structure factor
Ps and, for a more restricted set of variables, the helicity modulus ${\ensuremath{\rho}}_{s}$, as functions of temperature. Two finite-size scaling Ans\"atze for
Ps have been used: one for zero temperature and the other for finite temperatures. We have found that the system sustains superconductivity in the ground state up to a critical impurity concentration
fc, which increases with
U, at least up to
U=4 (in units of the hopping energy). Also, the normalized zero-temperature gap as a function of
f shows a maximum near $f\ensuremath{\sim}0.07$ for $2\ensuremath{\lesssim}U\ensuremath{\lesssim}6$. Analyses of the helicity modulus and of the pair structure factor led to the determination of the critical temperature as a function of
f for
U=3, 4, and 6: they also show maxima near $f\ensuremath{\sim}0.07$, with the highest
Tc increasing with
U in this range. We argue that, overall, the observed behavior results from both the breakdown of charge-density-wave-superconductivity degeneracy and the fact that free sites tend to ``push'' electrons toward attractive sites; the latter effect being more drastic at weak couplings.
The field-theoretical renormalization-group techniques previously developed for the description of the mestastable phases of isotropic, homogeneous,
d-dimensional systems, which (close to their critical temperature) can be described in terms of the Landau-Ginzburg-Wilson functional integral with a two-component order parameter, and applied to obtain the scaling form of the helicity modulus, defined as a measure of the free-energy increment associated with a metastable, helical twisting of the order parameter. The general scaling form of this quantity is determined even in the presence of an external helical field, for arbitrary helix pitch and temperature. In the case of vanishing external field, this quantity represents the superfluid density and is a function of the temperature and the superfluid velocity (which corresponds to the inverse of the helix pitch). In the limit of small velocities this expression exhibits the well-known Josephson asymptotic form ${\ensuremath{\rho}}_{s}\ensuremath{\sim}{({T}_{c}\ensuremath{-}T)}^{\ensuremath{\nu}(d\ensuremath{-}2)}={({T}_{c}\ensuremath{-}T)}^{2\ensuremath{\beta}\ensuremath{-}\ensuremath{\nu}\ensuremath{\eta}}$. An explicit $\ensuremath{\epsilon}=4\ensuremath{-}d$ expansion for the helicity modulus is presented in the case of arbitrary external field to the one-loop approximation.
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.
Two-dimensional generalized XY spin model on a triangular lattice is studied by means of Monte-Carlo simulations. The critical temperatures of Berezinskii-Kosterlitz-Thouless (BKT) phase transition are obtained by the method of helicity modulus. It is found that the results are consistent with those obtained by other methods. The vortex density and the vortex-antivortex pair formation energy are also obtained. The result shows that the critical temperature decreases with the increase of the generalization parameter q . While the vortex-antivortex pair formation energy increases with the increase of q when q>1 .