  A variational method for many electron system is applied to momentum
distribution calculations. The method uses a generating two-electron geminal
and the amplitudes of the occupancies of one particle natural orbitals as
variational parameters. It introduces correlation effects beyond the free
fermion nodal structure.
