A unified treatment of the most probable method: from microcanonical to grand canonical ensembles

This paper systematically explores a unified framework for applying the most probable method across the microcanonical, canonical, and grand canonical ensembles in statistical mechanics. In many textbook presentations, the most probable method is primarily illustrated in the microcanonical ensemble, which may obscure its broader applicability across all major ensembles. By extending this method to the canonical and grand canonical ensembles through the minimization of the Helmholtz free energy and the grand thermodynamic potential, respectively, we derive the corresponding equilibrium distributions. A comparative analysis of the core functions, constraint conditions, and mathematical structures of the extremum problems in these three ensembles reveals their intrinsic consistency, as well as the formal differences arising from distinct constraints. The study demonstrates that the most probable method constitutes a unified principle permeating all major ensembles in statistical mechanics, thereby deepening the understanding of equilibrium statistical physics and offering a more coherent perspective for teaching.