The Fermi hypernetted-chain (FHNC) theory is applied to spin-density-wave states combined with a Jastrow-type correlation operator. The commonly used approximation scheme FHNC/n is found to be unsuitable for dealing with this class of correlated wave functions. A new scheme, named FHNC-n, is proposed, and numerically tested for the one-dimensional (1D) Hubbard model and the Gutzwiller correlation operator. The results obtained are in very good agreement with the available variational Monte Carlo calculation. This scheme can be used to study 2D and 3D Hubbard models and can easily be extended to treat more general correlation operators, which include triplet and quartet correlations.
The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid
He3, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$.
F is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the
f(rij) embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements
Hmn,
Nmn of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the
Hmn are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (
Hmn) and (
Nmn), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (
Hmn) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (
Hmn). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (
Hmn). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any
F in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{},
Q-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid
He3 has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.
The Fermi-hypernetted-chain (FHNC) theory and the effective hypernetted-chain method are applied to study the composite-fermion (CF) states of the fractional quantum Hall effect. Using this theory we compute, in the thermodynamic limit, the correlation energy, radial distribution function, and static structure factor for all unprojected CF wave functions. The unprojected excitation gaps for $\ensuremath{\nu}=1/3,1/5$ were obtained by adopting in the FHNC a scheme previously used to compute nuclear matter excitation spectra. The results obtained so far are consistent with Monte Carlo simulations and small-number exact diagonalizations.
A simple alternative derivation of conventional cluster expansion formulas of Iwamoto and Yamada and others for the expectation value ≪H > of the energy of a many-Fermion system is offered. The assumed trial function is of the general form FΦp/(FΦp, FΦp)1/2, where Φp is a Slater determinant of appropriate single-particle functions and the factor F is attached to account for the possibility of strong interparticle correlations. The proposed scheme is applicable to finite as well as infinite systems and is easily generalized to the calculation of non-diagonal as well as diagonal matrix elements, with respect to the members of the set {FΦp/(FΦp, FΦp)1/2}. Explicit results for ≪H > are displayed to second cluster order, contributions being classified according to their assigned order of magnitude in a correlation parameter ξ=|Nω/Ω| of the sort introduced by Iwamoto and Yamada and Clark and Westhaus.