Energy estimators and calculation of energy expectation values in the path integral formalism

A recently developed method systematically improved the convergence of generic path integrals for transition amplitudes [A. Bogojević, A. Balaž, A. Belić, Phys. Rev. Lett. 94 (2005) 180403, A. Bogojević, A. Balaž, A. Belić, Phys. Rev. B 72 (2005) 064302, A. Bogojević, A. Balaž, A. Belić, Phys. Lett. A 344 (2005) 84]. This was achieved by analytically constructing a hierarchy of N  -fold discretized effective actions SN(p) labeled by a whole number p   and starting at p=1p=1 from the naively discretized action in the mid-point prescription. The derivation guaranteed that the level p   effective actions lead to discretized transition amplitudes differing from the continuum limit by a term of order 1/Np1/Np. Here we extend the applicability of the above method to the calculation of energy expectation values. This is done by constructing analytical expressions for energy estimators of a general theory for each level p  . As a result of this energy expectation values converge to the continuum as 1/Np1/Np. Finally, we perform a series of Monte Carlo simulations of several models, show explicitly the derived increase in convergence, and the ensuing speedup in numerical calculation of energy expectation values of many orders of magnitude.