Abstract:
We devise a hierarchy of computational algorithms to enumerate the microstates of a system
comprising N independent, distinguishable particles. An important challenge is to cope with
integers that increase exponentially with system size, and which very quickly become too large to
be addressed by the computer. A related problem is that the computational time for the most
obvious brute-force method scales exponentially with the system size which makes it difficult to
study the system in the large N limit. Our methods address these issues in a systematic and
hierarchical manner. Our methods are very general and applicable to a wide class of problems
such as harmonic oscillators, free particles, spin J particles, etc. and a range of other models for
which there are no analytical solutions, for example, a system with single particle energy spectrum
given by ε(p, q) = ε0(p2 +q4), where p and q are non-negative integers and so on. Working within
the microcanonical ensemble, our methods enable one to directly monitor the approach to the
thermodynamic limit (N ! 1), and in so doing, the equivalence with the canonical ensemble is
made more manifest. Various thermodynamic quantities as a function of N may be computed using
our methods; in this paper, we focus on the entropy, the chemical potential and the temperature.