Path integral hybrid Monte Carlo calculation of the bosonic oscillators

In this paper, various correlation functions for a bosonic many-body system have been calculated using the path integral hybrid Monte Carlo method proposed by the authors Chem. Phys. Lett. 308, 115 (1999)]. A system consisting of 32 independent particles in a harmonic well was chosen as a model problem. Calculations were carried out at two temperatures for the Boltzmann and Bose statistics. The large amplitude of the bosonic pair-correlation function around origin was found in comparison with the boltzmannoniac counterpart. This behavior can be interpreted by effective “attractive interaction” between particles arising from the Bose statistics.