Variational and diffusion quantum Monte Carlo methods are employed to investigate the zero-temperature phase diagram of three-dimensional homogeneous electron gas at very low density. Fermi fluid body-centered cubic Wigner crystal ground-state energies determined using Slater-Jastrow-backflow Slater-Jastrow many-body wave functions different densities spin polarizations in finite simulation cel...