1. Introduction
Boltzmannian statistical mechanics (BSM) partitions the phase space of a system into cells consisting of macroscopically indistinguishable microstates. These cells correspond to the macrostates, and the largest cell is singled out as the equilibrium macrostate. The connection is not conceptual: there is nothing in the concept of equilibrium tying equilibrium to the largest cell. So what justifies the association of equilibrium with the largest cell?
After introducing BSM (sec. 2), we discuss three justificatory strategies based on Boltzmann’s combinatorial argument, the Maxwell-Boltzmann distribution, and maximum entropy considerations, respectively. We argue that all three fail because they either suffer from internal difficulties or are restricted to systems with negligible interparticle forces (sec. 3). This prompts the search for an alternative answer. In analogy with the standard thermodynamic definition of equilibrium, we characterize equilibrium as the macrostate in which the system spends most of its time. We then present a new mathematical theorem proving that such an equilibrium macrostate indeed corresponds to the largest cell (sec. 4). This result is completely general in that it is not based on any assumptions about the system’s dynamics or the nature of interactions within the system.
2. Boltzmannian Statistical Mechanics
Let us briefly introduce BSM.Footnote 1 Consider a system consisting of n particles that is isolated from the environment and in a bounded container. The system’s state is specified by a point x = (q, p) (the microstate) in its 6n-dimensional phase space Γ. The system’s dynamics is determined by its classical Hamiltonian H(x). Energy is preserved, and therefore the motion is confined to the 6n − 1 dimensional energy hypersurface ΓE defined by H(x) = E, where E is the energy value. The solutions of the equations of motion are given by the phase flow ϕt on ΓE, where ϕt(x) is the state into which x ∈ ΓE evolves after t time steps. The Lebesgue σ-algebra is ΣE and, intuitively speaking, consists of all relevant subsets of ΓE. Space Γ is endowed with the Lebesgue measure μ, which is preserved under ϕt. This measure can be restricted to a measure μE on ΓE, which is preserved as well and is normalized; that is, μE(ΓE) = 1. The quadruple (ΓE, ΣE, μE, ϕt) is a measure-preserving dynamical system.
Assume that the system can be characterized by a set {v 1, … , vk} of macrovariables (k ∈ ℕ). The vi assume values in
, and capital letters Vi denote the values of vi. A particular set of values {V 1, … , Vk} defines a macrostate
. We only write M rather than
if the specific Vi do not matter. A set of macrostates is complete if and only if it contains all states a system can be in.
A crucial posit of BSM is supervenience: a system’s microstate uniquely determines its macrostate. Every macrostate M is associated with a macroregion ΓM consisting of all x ∈ ΓE for which the system is in M. For a complete set of macrostates, the ΓM form a partition of ΓE (they do not overlap and jointly cover ΓE).Footnote 2
The Boltzmann entropy of a macrostate M is
, where kB is the Boltzmann constant. The Boltzmann entropy of a system at time t, SB(t), is the entropy of the system’s macrostate at t:
, where x(t) is the system’s microstate at t and Mx (t) is the macrostate supervening on x(t).
We denote the equilibrium macrostate by M eq and its macroregion by
. A crucial aspect of the standard presentation of BSM is that Γeq takes up most of
. To facilitate the discussion, we introduce the term β-dominance:
is β-dominant if and only if
for β ∈ (1/2, 1]. Often equilibrium is characterized as a state in which β is close to 1, but nothing in what follows depends on a particular choice of β.
The characterization of equilibrium as a β-dominant state goes back to Ehrenfest and Ehrenfest (1959, 30). While different versions of BSM explain the approach to equilibrium differently, β-dominance is a key factor in all of them. Those who favor an explanation based on ergodic theory have to assume that
takes up the majority of ΓE because otherwise the system would not spend most of the time in
(e.g., Frigg and Werndl Reference Frigg and Werndl2011, Reference Frigg and Werndl2012). Those who see the approach to equilibrium as the result of some sort of probabilistic dynamics assume that
takes up most of ΓE because they assign probabilities to macrostates that are proportional to μ(ΓM) and that equilibrium comes out as the most likely state only if the equilibrium macroregion is β-dominant (e.g., Boltzmann Reference Boltzmann1877). Proponents of the typicality approach see dominance as the key ingredient in explaining the approach to equilibrium and sometimes even seem to argue that systems approach equilibrium because the equilibrium region takes up nearly all of phase space (e.g., Goldstein and Lebowitz Reference Goldstein and Lebowitz2004). We do not aim to adjudicate between these different approaches. Our question is a more basic one: Why is the equilibrium state β-dominant?
3. Justificatory Strategies
A look at the literature reveals three justificatory strategies. In practice these are often pursued side by side and seen as providing mutual support to each other. We assess each of them and argue that none of them is conclusive.
3.1. The Largest Number of Microstates and the Combinatorial Argument
The leading idea of the first justificatory strategy is that equilibrium is the macrostate that is compatible with the largest number of microstates. This strategy is exemplified by Boltzmann’s (Reference Boltzmann1877) combinatorial argument.Footnote 3 The state of one particle is determined by a point in its six-dimensional state space Γμ, and the state of a system of n identical particles is determined by n points in this space. Since the system is confined to a finite container and has constant energy E, only a finite part of Γμ is accessible. Boltzmann partitions the accessible part of Γμ into cells of equal size δω whose dividing lines run parallel to the position and momentum axes. The result is a finite partition
, m ∈ ℕ. The cell in which a particle’s state lies is its coarse-grained microstate. The coarse-grained microstate of the entire gas, called an arrangement, is given by a specification of the coarse-grained microstate of each of particle.
The system’s macroproperties depend only on how many particles there are in each cell and not on which particles these are. A specification of the ‘occupation number’ of each cell is known as a distribution D = (n 1, n 2, … , nm), where ni is the number of particles whose state is in cell ωi. Since m and n are finite, there are only finitely many distributions D 1, … , Dk. Each distribution is compatible with several arrangements, and the number G(D) of arrangements compatible with a given distribution D = (n 1, n 2, … , nm) is
Every microstate x of ΓE is associated with exactly one distribution D(x). One then defines the set ΓD of all x that are associated with a distribution D:
Since macroproperties are fixed by the distribution, distributions are associated with macrostates. So we ask: Which of the distributions is the equilibrium distribution? Now Boltzmann’s main idea enters the scene: equilibrium is the macrostate that is compatible with the largest number of microstates. To determine the equilibrium distribution, Boltzmann assumed that the energy ei of particle i depends only on the cell in which it is located. Then the total energy is
He furthermore assumed that the number of cells in Ω is small compared to the number of particles (allowing him to use Stirling’s formula). With the further trivial assumption that
, Boltzmann shows that μE(ΓD) is maximal when
where γ and λ are parameters that depend on n and E. This is the discrete version of the Maxwell-Boltzmann distribution. Thus, the equilibrium macrostate corresponds to the Maxwell-Boltzmann distribution.Footnote 4
Its ingenuity notwithstanding, the combinatorial argument faces a number of important problems. The first is that it only applies to systems of noninteracting particles (Uffink Reference Uffink, Buttereld and Earman2007, 976–77). It provides a reasonable approximation for systems with negligible interparticle forces, but any other system is beyond its scope. Statistical mechanics (SM) ought to be a general theory of matter, and so this is a serious limitation.
The second problem is the absence of a conceptual connection between equilibrium in thermodynamics (TD) and the idea that the equilibrium macrostate is the one that is compatible with the largest number of microstates. In TD equilibrium is defined as the state to which isolated systems converge when left to themselves and which they never leave once they have reached it. This has very little, if anything, in common with the kind of considerations underlying the combinatorial argument. This is a problem for anyone who sees BSM as a reductionist enterprise. While the precise contours of the reduction of TD to SM remain controversial, we are not aware of any contributors who maintain radical antireductionism. Thus, the disconnect between the two notions of equilibrium is a serious problem.
Two replies come to mind. The first points out that since
is the largest subset of ΓE, systems approach equilibrium and spend most of their time in
. This shows that the BSM definition of equilibrium is a good approximation to the TD definition. This is not true in general. Whether a system spends most of its time in the β-dominant
depends on the dynamics. If, for instance, the dynamics is the identity function, it is not true that a system out of equilibrium approaches equilibrium and spends most of its time there. The second reply points out that we know for independent reasons that the equilibrium distribution is the Maxwell-Boltzmann distribution. This argument will be discussed in the next subsection, and our conclusion will be guarded.
Finally, the combinatorial argument (even if successful) shows that the equilibrium macrostate is larger than any other macrostate. However, as Lavis (Reference Lavis2005) points out, this need not imply that the equilibrium is β-dominant. There may be a large number of smaller macrostates that jointly take up a large part of ΓE. So the combinatorial argument does in fact not show that equilibrium is β-dominant.
3.2. The Maxwell-Boltzmann Distribution
According to the next justificatory strategy, a system is in equilibrium when its particles approximately satisfy the Maxwell-Boltzmann distribution (eq. [4]; e.g., Penrose Reference Penrose1989). This approach is misguided because the Maxwell-Boltzmann distribution is in fact the equilibrium distribution for a limited class of systems only, namely, for systems consisting of particles with negligible interparticle forces. For particles with nonnegligible interactions, different distributions correspond to equilibrium (Gupta Reference Gupta2003). Furthermore, for many simple models such as the Ising model (Baxter Reference Baxter1982) or the Kac ring (Lavis Reference Lavis2008) the equilibrium macrostate also does not correspond to the Maxwell-Boltzmann distribution.
This is no surprise given that the two common derivations of the distribution in effect assume that particles are noninteracting. Boltzmann’s (Reference Boltzmann1877) derivation is based on equation (3), the assumption that the total energy is the sum of the energy of the individual particles. This is true only if the particles are noninteracting (i.e., for ideal gases). While many expect that the argument also goes through for dilute gases (where this assumption holds approximately), the argument fails for nonnegligible interactions. In Maxwell’s 1860 derivation (see Uffink Reference Uffink, Buttereld and Earman2007) the noninteraction assumption enters via the postulate that the probability distributions in different spatial directions can be factorized, which is true only if there is no interaction between particles. For these reasons, the Maxwell-Boltzmann distribution is the equilibrium distribution for a limited class of systems only and cannot be taken as a general definition of equilibrium.
3.3. Maximum Entropy
A third strategy justifies dominance by maximum entropy considerations along the following lines:Footnote 5 we know from TD that, if left to itself, a system approaches equilibrium, and equilibrium is a maximum entropy state. Hence, the Boltzmann entropy of a macrostate SB is maximal in equilibrium. Since SB is a monotonic function, the macrostate with the largest Boltzmann entropy is also the largest macrostate, which is the desired conclusion.
There are serious problems with the understanding of TD in this argument as well as with its implicit reductive claims. First, that a system, when left to itself, reaches equilibrium where entropy is maximal is often taken to be a consequence of the second law of TD, but it is not. As Brown and Uffink (Reference Brown and Uffink2001) pointed out, that systems tend to approach equilibrium has to be added as an independent postulate, which they call the Minus First Law. But even if TD is amended with the Minus First Law, the conclusion does not follow. TD does not attribute an entropy to systems out of equilibrium. Thus, characterizing the approach to equilibrium as a process of entropy increase is meaningless from a TD point of view.
Even if all these issues could be resolved, there still would be a question why the fact that the TD entropy reaches a maximum in equilibrium would imply that the same holds for the Boltzmann entropy. To justify this inference, one would have to assume that the TD entropy reduces to the Boltzmann entropy. But it is far from clear that this is so. A connection between the TD entropy and the Boltzmann entropy has been established for ideal gases only, where the Sackur-Tatrode formula can be derived from BSM, which shows that both entropies have the same functional dependence on TD state variables. No such results are known for systems with interactions. Furthermore, there are well-known differences between the TD entropy and the Boltzmann entropy. Most importantly, the TD entropy is extensive while the Boltzmann entropy is not (Ainsworth Reference Ainsworth2012). But an extensive concept cannot reduce to a nonextensive concept (at least not without further qualifications). For these reasons, we conclude that maximum entropy considerations cannot be used to argue for the β-dominance of the equilibrium state.
4. Rethinking Equilibrium
The failure of standard justificatory strategies prompts the search for an alternative answer. In this section, we propose an alternative definition of equilibrium and introduce a new mathematical theorem proving that the equilibrium state thus defined is β-dominant.
The above strategies run into difficulties because there is no clear connection between the TD definition of equilibrium and β-dominance. Our aim is to provide the missing connection by taking as a point of departure the standard TD definition of equilibrium and then exploiting supervenience to ‘translate’ this macrodefinition into microlanguage.
The following is a typical TD textbook definition of equilibrium: “A thermodynamic system is in equilibrium when none of its thermodynamic properties are changing with time” (Reiss Reference Reiss1996, 3). In more detail: equilibrium is the state to which an isolated system converges when left on its own and which it never leaves once it has been reached (Callender Reference Callender2001; Uffink Reference Uffink2001). Equilibrium in TD is unique in the sense that the system always converges toward the same equilibrium state. This leads to the following definition (the qualification ‘strict’ will become clear later):
Definition 1: Strict BSM Equilibrium. Consider an isolated system S whose macrostates are specified in terms of the macrovariables {v 1, … , vk}, described as the measure-preserving dynamical system (ΓE, ΣE, μE, ϕt). Let M(x) be the macrostate that supervenes on microstate x ∈ Γ. Let
be the set of all microstates on which
supervenes. If there is a macrostate
satisfying the following condition, then it is the strict BSM equilibrium state of S: for all initial states x ∈ ΓE at t 0 there exists a time t* such that
for all t ≥ t*.Footnote 6 We then write
.
Note that this definition incorporates the Minus First Law of TD.
Before reflecting on this definition, we want to add a brief comment about reductionism. Reductive eliminativists may feel that a definition of equilibrium in SM that is based on ‘top down translation’ of its namesake in TD undermines the prospect of a reduction of TD to SM. They would argue that equilibrium has to be defined in purely mechanical terms and must then be shown to line up with the TD definition of equilibrium.
This point of view is not the only game in town, and reduction can be had even if equilibrium is defined ‘top down’ (as in the above definition).Footnote 7 First, whether the above definition undercuts a reduction depends on one’s concept of reduction. For someone with a broadly Nagelian perspective, there is no problem: the above definition provides a bridge law, which allows the derivation of the requisite macroregularities from the laws of the microtheory. And a similar argument can be made in the framework of New Wave Reductionism. Second, equilibrium is a macroconcept: when describing a system as being in equilibrium, we look at it in terms of macroproperties. From a micro point of view there are only molecules bouncing around. They always bounce—there is no such thing as a relaxation of particle motion to an immutable state. Hence, a definition of equilibrium in macroterms is no heresy.
Definition 1 is too rigid for two reasons. The first reason is Poincaré recurrence: as long as the ‘M’ in SM refers to a mechanical theory that conserves phase volume (and there is widespread consensus that this is the case),Footnote 8 any attempt to justify an approach to strict equilibrium in mechanical terms is doomed to failure. The system will at some point return arbitrarily close to its initial condition, violating strict equilibrium (Uffink Reference Uffink, Buttereld and Earman2007; Frigg Reference Frigg and Rickles2008). The second reason is that such a justification is not only unattainable but also undesirable. Experimental results show that equilibrium is not the immutable state that classical TD presents us with because systems exhibit fluctuations away from equilibrium (Wang et al. Reference Wang, Sevinck, Mittag, Searles and Evans2002). Thus, strict equilibrium is actually unphysical.
Consequently, strict definitions of equilibrium are undesirable both for theoretical and experimental reasons. So let us relax the condition that a system has to remain in equilibrium for all t ≥ t* by the weaker condition that it has to be in equilibrium most of the time:
Definition 2: BSM α-Equilibrium. Consider the same system as in definition 1. Let fM,x(t) be the fraction of time of the interval [t 0, t 0 + t] in which the system’s state is in M when starting in initial state x at t 0, and let α be a real number in [0.5, 1]. If there is a macrostate
satisfying the following condition, then it is the α-equilibrium state of S: for all initial states x ∈ ΓE,
in the limit t → ∞. We then write
.
An obvious question concerns the value of α. Often the assumption seems to be that α is close to 1. This is reasonable but not the only possible choice. For our purposes nothing hangs on a particular choice of α, and so we leave it open what the best choice would be.
One last step is needed to arrive at the definition of equilibrium suitable for BSM. It has been pointed out variously that in SM, unlike in TD, we should not expect every initial condition to approach equilibrium (see, e.g., Callender Reference Callender2001). Indeed, it is reasonable to allow for a set of very small measure ε for which the system does not approach equilibrium:
Definition 3: BSM α-ε-Equilibrium. Let S and fM,x(t) be as above. Let α be a real number in [0.5, 1], let 1 > ε ≥ 0 be a small real number, and let Y be a subset of ΓE such that
. If there is a macrostate
satisfying the following condition, then it is the α-ε-equilibrium state of S: for all initial states x ∈ Y,
in the limit t → ∞. We then write
.
Let us introduce the characteristic function of ΓM, 1M(x): 1M(x) = 1 for x ∈ ΓM and 0 otherwise. Definition 3 implies that, for all x ∈ Y,Footnote 9
An important assumption in this characterization of equilibrium is that μE (and not some other measure) is the relevant measure. It is often argued that μE can be interpreted as a probability or typicality measure (Frigg and Hoefer Reference Frigg, Hoefer, Dieks, Gonzales, Hartmann, Weber, Stadler and Übel2010; Werndl Reference Werndl2013). The condition then says that the system’s state spends more than a fraction α of its time in equilibrium with probability 1 − ε, or that typical initial conditions lie on trajectories that spend more than α of their time in equilibrium.
We contend that the relevant notion of equilibrium in BSM is α-ε-equilibrium. The central question then becomes: Why is the α-ε-equilibrium state β-dominant? Definition 3 in no way prejudges this question: it says nothing about the size of
, nor does it in an obvious sense imply anything about it.
That
is β-prevalent follows from the following theorem, which we prove in the appendix:
Equilibrium Theorem: If
is an α-ε-equilibrium of system S, then
.
We emphasize that the theorem is completely general in that no dynamical assumption is made (in particular, it is not assumed that the system is ergodic). So the theorem also applies to strongly interacting systems such as solids and liquids.
The equilibrium theorem is the centerpiece of our account. It shows in full generality that if the system S has an α-ε-equilibrium, then the equilibrium state is β-dominant for
.Footnote 10 This provides the sought-after justification of the β-dominance of the equilibrium state.
The equilibrium theorem makes the conditional claim that if there is an α-ε-equilibrium, then
. As with all conditionals, the crucial and often vexing question is whether, and under what conditions, the antecedent holds. Some systems do not have equilibria. For instance, if the dynamics is given by the identity function, then no approach to equilibrium takes place, and the antecedent of the conditional is wrong. By contrast, epsilon ergodicity allows for an equilibrium state to exist (Frigg and Werndl Reference Frigg and Werndl2011). This raises the question under which circumstances the antecedent is true, which is an important question for future research.
5. Conclusion
BSM partitions the phase space of a system into cells of macroscopically indistinguishable microstates. These cells are associated with the system’s macrostates, and the largest cell is identified with equilibrium. What justifies the association of equilibrium with the largest cell? We discussed three justificatory strategies that can be found in the literature: that equilibrium is the macrostate compatible with the largest number of microstates, that equilibrium corresponds to the Maxwell-Boltzmann distribution and that most states are characterized by that distribution, and that equilibrium is the maximum entropy state. We argued that none of them is successful. This prompted the search for an alternative answer. We characterized equilibrium as the state in which the system spends most of its time and presented a new mathematical theorem proving that such an equilibrium state indeed corresponds to the largest cell. This result is completely general in that it is not based on any assumptions about either the system’s dynamics or the nature of interactions within the system. It therefore provides the first fully general justification of the claim that the equilibrium state takes up most of the accessible part of the system’s phase space.
Appendix Proof of the Equilibrium Theorem
The proof appeals to the ergodic decomposition theorem (cf. Petersen Reference Petersen1983, 81), stating that for a dynamical system (ΓE, ΣE, μE, ϕt) the set ΓE is the disjoint union of sets X ω, each equipped with a σ-algebra
and a probability measure μω, and ϕt acts ergodically on each
. The indexing set is also a probability space (Ω, ΣΩ, P), and for any square integrable function f it holds that

Application of the ergodic decomposition theorem for
yields
For an ergodic system
, the long-run time average equals the phase average. Hence, for almost all x ∈ Xω
From requirement (5) and because ϕt acts ergodically on each
, for almost all x ∈ X ω, X ω ⊆ Y
Let us first consider the case ε = 0; that is, μE(Y) = 1. Here, from equation (A2)
Hence, if ε ≥ 0, it follows from equation (A2) that