Out-of-equilibrium Monte Carlo simulations of a classical gas with Bose-Einstein statistics
Algorithms to determine transition probabilities in Monte Carlo simulations are tested using a system of classical particles with effective interactions which reproduce Bose-Einstein statistics. The system is appropriate for testing different Monte Carlo simulation methods in out-of-equilibrium situations since non equivalent results are produced. We compare mobility numerical results obtained with transition probabilities derived from Glauber and Metropolis algorithms. Then, we compare these with a recent method, the interpolation algorithm, appropriate for non-equilibrium systems in homogeneous substrata and without phase transitions. The results of mobility obtained from the interpolation algorithm are qualitatively verified with molecular dynamics simulations for low concentrations.