Ground state properties of fermionic systems with the Density Matrix Projection method