The effects of the Gutzwiller projection on a BCS wave function with varying particle number are considered. We show that a fugacity factor has to be introduced in these wave functions when they are Gutzwiller projected, and derive an expression for this factor within the Gutzwiller approximation. We examine the effects of the projection operator on BCS wave functions by calculating the average number of particles before and after projection. We also calculate particle number fluctuations in a projected BCS state. Finally, we point out the differences between projecting BCS wave functions in the micro- and grand-canonical schemes, and discuss the relevance of our results for variational Monte Carlo studies.
Gutzwiller-projected BCS wave function or the resonating-valence-bond (RVB) state in the two-dimensional (2D) extended $t\ensuremath{-}J$ model is investigated by using the variational Monte Carlo technique. We show that the results of ground-state energy and excitation spectra calculated in the grand-canonical scheme allowing particle number to fluctuate are essentially the same as previous results obtained by fixing the number of particle in the canonical scheme if the grand thermodynamic potential is used for minimization. To account for the effect of Gutzwiller projection, a fugacity factor proposed by Laughlin and Anderson a few years ago has to be inserted into the coherence factor of the BCS state. Chemical potential, particle number fluctuation, and phase fluctuation of the RVB state, difficult or even impossible to be calculated in the canonical ensemble, have been directly measured in the grand-canonical picture. We find that except for $La\ensuremath{-}214$ materials, the doping dependence of chemical potential is consistent with experimental findings on several cuprates. Similar to what has been reported by scanning tunneling spectroscopy experiments, the tunneling asymmetry becomes much stronger as doping decreases. We found a very large enhancement of phase fluctuation in the underdoped regime.
The Gutzwiller approximation (GA) for Gutzwiller-projected grand canonical wave functions with fugacity factors is investigated in detail. Our systems in general contain inhomogeneity and local magnetic moments. In deriving renormalization formulas, we also derive or estimate terms of higher powers of intersite contractions neglected in the conventional GA. We examine several different constraints, i.e., local/global spin-dependent/independent particle-number conservation. Out of the four, the local spin-dependent constraint seems the most promising at present. An improved GA derived from it agrees with the variational Monte Carlo method better than the conventional GA does. The corrections to the conventional GA can be interpreted as two-site correlation including the phase difference of configurations. Furthermore, projected quasiparticle excited states are orthogonal to each other within the GA. Using these states, spectral weights are calculated. We show that asymmetry between electron addition and removal spectra can appear by taking into account the higher powers of the intersite contractions in the case of the
d-wave superconductors and the Fermi sea; the addition is smaller than the removal. However, the asymmetry is quite weak especially near the Fermi level. In contrast, projected
s-wave superconductors can have the opposite asymmetry (addition larger than removal) especially near the Fermi level. In addition, formulas from the other three constraints are also derived, which may be useful depending on purposes.