Path integral molecular dynamics for anyons, bosons and fermions

In this article we develop a general method to numerically calculate physical properties for a system of anyons with path integral molecular dynamics. We provide a unified method to calculate the thermodynamics of identical bosons, fermions and anyons. Our method is tested and applied to systems of anyons, bosons and fermions in a two-dimensional harmonic trap. We also consider a method to calculate the energy for fermions as an application of the path integral molecular dynamics to simulate the anyon model.