Abstract
We summarize the essential ingredients, which enabled us to derive the path-integral for a system of harmonically interacting spin-polarized identical particles in a parabolic confining potential, including both the statistics (Bose–Einstein or Fermi–Dirac) and the harmonic interaction between the particles. This quadratic model, giving rise to repetitive Gaussian integrals, allows to derive an analytical expression for the generating function of the partition function. The calculation of this generating function circumvents the constraints on the summation over the cycles of the permutation group. Moreover, it allows one to calculate the canonical partition function recursively for the system with harmonic two-body interactions. Also, static one-point and two-point correlation functions can be obtained using the same technique, which make the model a powerful trial system for further variational treatments of realistic interactions