  Two methods are compared that are used in path integral simulations. Both
methods aim to achieve faster convergence to the quantum limit than the
so-called primitive algorithm (PA). One method, originally proposed by
Takahashi and Imada, is based on a higher-order approximation (HOA) of the
quantum mechanical density operator. The other method is based upon an
effective propagator (EPr). This propagator is constructed such that it
produces correctly one and two-particle imaginary time correlation functions in
the limit of small densities even for finite Trotter numbers P. We discuss the
conceptual differences between both methods and compare the convergence rate of
both approaches. While the HOA method converges faster than the EPr approach,
EPr gives surprisingly good estimates of thermal quantities already for P = 1.
Despite a significant improvement with respect to PA, neither HOA nor EPr
overcomes the need to increase P linearly with inverse temperature. We also
derive the proper estimator for radial distribution functions for HOA based
path integral simulations.
