Estimating ground state entanglement entropy using path integral molecular dynamics