  We construct a quantum Monte Carlo algorithm for interacting fermions using
the two-body density as the fundamental quantity. The central idea is mapping
the interacting fermionic system onto an auxiliary system of interacting
bosons. The correction term is approximated using correlated wave functions for
the interacting system, resulting in an effective potential that represents the
nodal surface. We calculate the properties of 3He and find good agreement with
experiment and with other theoretical work. In particular, our results for the
total energy agree well with other calculations where the same approximations
were implemented but the standard quantum Monte Carlo algorithm was used
