many body wavefunction and exchange correlation Everywhere I ready about HF or DFT the term exchange correlation functional comes up. I have a couple of fundamental questions about these:
1) Books say that the correlation energy is the difference between the exact energy (lets say we've found that somehow) and the hartree-fock energy and that the exchange term is accounted for exactly in the hartree-fock energy. While these claims depend on what exactly you include in the correlation energy, isn't the second part on getting the exact exchange totally wrong?? This is because the wave functions that are used to construct/evaluate the hartree fock energy are not the "exact" single electron wavefunctions to begin with so whatever mathematical construct we use to get the exchange energy term can only yield us a false answer. We don't know if there are things called single electron wavefunctions! Or am I wrong? I mean lets say we got the "exact" single electron wavefunctions from somehwhere. Then can we calculate the exchange energy the way its been calculated in hartree-fock method ? Is there a fundamental way to combine all the "exact" single electron wavefunctions to give us the total wavefunction?
2) Likewise, all books say that when we take into account coulombic repulsion, only using the composite wave function can only give us the same exact repulsion. If we use individual electron wavefunctions or densities, then this becomes a mean field approximation. I really do not understand this. In the end an electron is particle at least to the extent that it cannot feel repulsion from itself. So if it feels the repulsion from other charge densities (given we know what they are), why should this be a problem ?. I put forward the same proposal again - if i can get the "exact" single electronic wavefunctions from somewhere, then wouldn't it feel the repulsion of other electrons in the way this normal "averaged" coulomb term describes.
To be more specific, I provide this document http://www.physics.metu.edu.tr/~hande/teaching/741-lectures/lecture-05.pdf where on page 2, Eq. (11) the real accurate term has been described for calculating the exact repulsion from the many-body density. It says there "It can be proven that this (electron electron repulsion) term cannot be written in terms of the single-particle density but instead only in terms of the two-particle density". Can someone please direct me to the proof?
3) Are we in the end saying that there is no rule of physics which can tell how one exact single many body wavefunction for N electrons can be decomposed into N exact 1-body wavefunctions or vice versa. Is this really the challenge that we haven't been able to solve in order to get to the correct exchange-correlation functional? If yes please explain a bit further. Also, has this all got to do anything with the second quantization, whatever that means.
Thanks a lot for your attention. It would be great if you could give a long answer. with sources of the mathematical proofs wherever required.
 A: Let's take a usual non-relativistic many body bosonic (there is not much difference in taking bosons of fermions, at least for the purpose here) Hamiltonian for $N$ particles of mass $1/2$ on $L^2_{s}(\mathbb{R}^{Nd})$ (where $s$ stands for symmetric functions, wrt exchange of particles):
$$H=H_0+H_{int}=\sum_{j=1}^N -\Delta_{x_j}+\frac{1}{N}\sum_{j<i}V(x_i-x_j)\; ,$$
where the symmetric $V(\cdot)$ models the two-body interaction between particles (in the case of Coulomb, $V(x)=\pm \frac{1}{\lvert x\rvert}$).
It is pretty clear that a factorized structure for particles is not preserved by the interaction term: think of the simplest factorized symmetric function
$$\Phi_N(x_1,\dotsc,x_N)=\varphi(x_1)\dotsm\varphi(x_N)\; ,$$
where $\varphi\in L^2(\mathbb{R}^d)$ is a single particle wavefunction. 
While $H_0\Phi_N$ still gives a factorized structure (it suffices that $\varphi\in H^2(\mathbb{R}^d)\subset L^2(\mathbb{R}^d)$, and we get $H_0\Phi_N(x_1,\dotsc,x_N)=(-\Delta\varphi)(x_1)\dotsm(-\Delta\varphi)(x_N)$); $H_{int}\Phi_N$ is no more factorized, since in general (and specifically for Coulomb) it is not true that we can write $V(x_i-x_j)$ as $V'(x_i)V'(x_j)$.
So it is not a lack of knowledge, but a fact that systems of particles with Coulomb two-body interactions (and other potentials) are not describable by means of single-particle wavefunctions, because the interaction yields unavoidable correlations that affect the evolution of the system.
Now, said that, we would like to be able to describe the behavior of a single particle without studying the whole system, in suitable situations (because it is far more easy on a computational standpoint). This is done using the so-called mean field approximation, that is precise only in the limit $N\to\infty$ of many (infinite) particles. The huge number of particles makes possible to describe the motion of each particle by means of an effective self-interaction, that is an expression of the global action of all the other particles in the system (the mean field). The price to pay is that the mean field equation is no more linear. For the system above, it is the Hartree equation on $L^2(\mathbb{R}^d)$:
$$i\partial_t\varphi=-\Delta\varphi + (V*\lvert \varphi\rvert^2)\varphi\; .$$
However its validity holds only in the limit $N\to\infty$; if you have a very large $N$ but finite, you make an error using Hartree equation. However, the error can be quantified, and depends on how large $N$ is.
