I'm reading McQuarrie's Physical chemistry textbook [1]. On p. 281, the author explains:
If we use a more flexible trial function of the form in which $\psi(\mathbf{r_1},\mathbf{r_2})$ is a product of one-electron functions, or orbitals,
$$\psi(\mathbf{r_1},\mathbf{r_2})= \phi(\mathbf{r_1})\phi(\mathbf{r_2})\tag{8.15}$$
and allow $\phi(r)$ to be completely general, then we reach a limit that is both practical and theoretical. […]
This limiting value is the best value of the energy that can be obtained using trial function of the form of a product of one-electron wave equations (Equation 8.15). This limit is called the Hartree-Fock limit.
But why does Hartree–Fock limit exist? If we increase the number of parameters of the trial function we use, we may expect better calculation results. Nevertheless, why is the wave function multiplied by the two trial functions limited to the accuracy that can be achieved?
My guess is that the functional form does not properly consider the interelectronic correction term by treating the two trial functions independently. But can't this difference be overcome by using more parameters? I wonder if I can theoretically explain the existence of the Hartree–Fock limit.
In addition, if this limit exists, the function will have to be altered to produce better calculation results. I wonder what theory computational chemistry has developed about this.
Reference
McQuarrie, D. A.; Simon, J. D. Physical Chemistry: A Molecular Approach; University Science Books: Sausalito, California, 1997. ISBN 978-0-935702-99-6.