We calculate the off-diagonal density matrix of the homogeneous electron gas at zero temperature using unbiased Reptation Monte Carlo for various densities and extrapolate the momentum distribution, and the kinetic and potential energies to the thermodynamic limit. Our results on the renormalization factor allows us to validate approximate G0W0 calculations concerning quasiparticle properties over a broad density region (1 ≤ rs 10) and show that near the Fermi surface, vertex corrections and self-consistency aspects almost cancel each other out. The uniform electron gas (jellium) is one of the most fundamental models for understanding electronic properties in simple metals and semiconductors. Knowledge of its ground state properties, and, in particular, of modifications due to electron correlation are at the heart of all approximate approaches to the many-electron problem in realistic models. Quantum Monte Carlo methods (QMC) [1] have provided the most precise estimates of the correlation energy, electron pair density and structure factor of jellium; basic quantities for constructing and parameterizing the exchange-correlation energy used in density functional theory (DFT) [2].Correlations modify the momentum distribution, n k , of electrons, and introduce deviations from the ideal Fermi-Dirac step-function. The magnitude of the discontinuity at the Fermi surface (k F ), the renormalization factor Z, quantifies the strength of a quasi-particle excitation [3] and plays a fundamental role in Fermi liquid and many-body perturbation theory (GW) for spectral quantities. Whereas the momentum distribution (as well as other spectral information) is inaccessible in current Kohn-Sham DFT formulations, the reduced singleparticle density matrix -the Fourier transform of n k in homogeneous systems -is the basic object in the so-called density-matrix functional theory [4]; these theories rely on knowledge of n k of jellium. Inelastic x-ray scattering measurement of the Compton profile of solid sodium [5] have determined n k , but experiments for elements with different electronic densities are less conclusive.In this paper, we calculate n k for the electron gas (jellium) by QMC in the density region 1 ≤ r s ≤ 10. Here, r s = (4πna 3 B /3) −3 is the Wigner-Seitz density parameter, n is the density, and a B = 2 /me 2 is the Bohr radius. In contrast to previous calculations [6], our calculations are based on more precise backflow (BF) wave functions [7], and a careful extrapolation to the thermodynamic The momentum distribution (n k ) of the unpolarized electron gas for various densities extrapolated to the thermodynamic limit. The inset shows the extrapolation of n k for rs = 5 from a system with N = 54 electrons to the thermodynamic limit, N → ∞,leading to a significant reduction of the renormalization factor Z.limit [8,9]. Similar to the worm algorithm in finite temperature path-integral and lattice Monte Carlo [10, 11], we have extended Reptation Monte Carlo (RMC) [12] to include the off-diagonal density matrix in order to o...