A quantum Monte Carlo method of determining Jastrow-Slater and correlated m
ultideterminant wave functions for which the energy is stationary with resp
ect to variations in the single-particle orbitals is presented. A potential
is determined by a least-squares fitting of fluctuations in the energy wit
h a linear combination of one-body operators. This potential is used in a s
elf-consistent scheme for the orbitals whose solution ensures that the ener
gy of the correlated wave function is stationary with respect to variations
in the orbitals. The method is feasible for atoms, molecules, and solids a
nd is demonstrated for the beryllium, carbon, and neon atoms and for the so
lid diamond. (C) 2000 American Institute of Physics. [S0021-9606(00)30406-8
].