- Canonical ensemble
A

**canonical ensemble**instatistical mechanics is astatistical ensemble representing aprobability distribution of microscopic states of the system. The probability distribution is characterised by the proportion "p_{i}" of members of the ensemble which exhibit a measurable macroscopic state "i", where the proportion of microscopic states for each macroscopic state "i" is given by theBoltzmann distribution ,:$p\_i\; =\; frac\{1\}\{Z\}e^\{-\; E\_i/(kT)\}\; =\; e^\{-(E\_i\; -A)/(kT)\}$

where "E"

_{"i"}is the energy of state "i". It can be shown that this is the distribution which is most likely, if each system in the ensemble can exchange energy with aheat bath , or alternatively with a large number of similar systems. Equivalently, it is the distribution which hasmaximum entropy for a given average energy <"E"_{i}>.It is also referred to as an NVT ensemble: the number of particles ("N"), the volume ("V"), of each system in the ensemble are the same, and the ensemble has a well defined temperature ("T"), given by the temperature of the heat bath with which it would be in equilibrium.

The quantity "k" is

Boltzmann's constant , which relates the units of temperature to units of energy. It may be suppressed by expressing the absolute temperature usingthermodynamic beta , $eta\; =\; 1/(kT)$.The quantities "A" and "Z" are constants for a particular ensemble, which ensure that $Sigma\; p\_i$ is normalised to 1. "Z" is therefore given by

:$Z\; =\; sum\; e^\{-\; E\_i/(kT)\}\; =\; sum\; e^\{-eta\; E\_i\}$.

This is called the partition function of the canonical ensemble. Specifying this dependence of "Z" on the energies "E"

_{i}conveys the same mathematical information as specifying the form of "p"_{i}above.The canonical ensemble (and its partition function) is widely used as a tool to calculate thermodynamic quantites of a system under a fixed temperature. This article derives some basic elements of the canonical ensemble. Other related thermodynamic formulas are given in the partition function article. When viewed in a more general setting, the canonical ensemble is known as the

Gibbs measure , where, because it has theMarkov property of statistical independence, it occurs in many settings outside of the field of physics.**Deriving the Boltzmann factor from ensemble theory**Let $E\_i,$ be the energy of the microstate $i,$ and suppose there are $n\_i,$ members of the ensemble residing in this state. Further we assume the total number of systems in the ensemble, $mathcal\{N\},$, and the total energy of all systems of the ensemble, $mathcal\{E\},$, are fixed, i.e.,

:$mathcal\{N\}=\; sum\_i\; n\_i\; ,\; ,$

:$mathcal\{E\}=\; sum\_i\; n\_i\; E\_i\; ,.$

Since systems in the ensemble are distinguishable, for each set $\{n\_i\}\; ,$, the number of ways of shuffling systems is equal to

:$W\; (\{n\_i\})\; =\; mathcal\{N\}!/\; prod\_\{i\}\; n\_i!\; ,\; .$

So for a given $\{n\_i\},$, there are $W(\{n\_i\}),$ rearrangements that specify the same state of the ensemble.

The most probable distribution is the one that maximizes $W\; (\{n\_i\}),$. The probability for any other distribution to occur is extremely small in the limit $mathcal\{N\}\; ightarrow\; infty\; ,$.To determine this distribution, one should maximize $W\; (\{n\_i\}),$ with respect to the $n\_i,$'s, under two constraints specified above. This can be done by using two

Lagrange multipliers $alpha\; ,$ and $eta,$. (The assumption that $mathcal\{N\}\; ightarrow\; infty\; ,$ would be invoked in such calculation, which allows one to applyStirling's approximation .) The result is:$n\_i\; =\; e^\{-alpha\; -eta\; E\_i\}\; ,$.

This distribution is called the canonical distribution. To determine $alpha\; ,$ and $eta,$, it is useful to introduce the

partition function as a sum over microscopic states:$Z(eta)\; =\; sum\_j\; e^\{-eta\; E\_j\}\; .,$

Comparing with thermodynamic formulae, it can be shown that $eta,$, is related to the absolute temperature $T,$ as, $eta=1/k\_B\; T,$. Moreover the expression

:$-\; ln\; Z(eta)\; /eta,$

is identified as the

Helmholtz free energy $F$. A derivation is given here. Consequently, from the partition function we can obtain the average thermodynamic quantities for the ensemble. For example, the average energy among members of the ensemble is:$langle\; E\; angle\; =\; frac\{\; mathcal\{E\{\; mathcal\{N\}\; \}\; =\; -\; frac\{partial\}\{partial\; eta\; \}\; ln\; Z(eta)\; ,$.

This relation can be used to determine $eta,$. $alpha,$ is determined from

:$e^\{alpha\}\; =\; Z(eta)/\; mathcal\{N\},$.

**A derivation from heat-bath viewpoint**Define the following:

* S - the system of interest

* S′ - the heat reservoir in which S resides; S is small compared to S′

* S* - the system consisting of S and S′ combined together

* "m" - an indexing variable which labels all the available energy states of the system S

* "E"_{m}- the energy of the state corresponding to the index m for the system S

* "E"′ - the energy associated with the heat bath

* "E"* - the energy associated with S*

* "Ω"′("E") - denotes the number of microstates available at a particular energy E for the heat reservoir.It is assumed that the system S and the reservoir S′ are in thermal equilibrium. The objective is to calculate the set of probabilities "p"

_{"m"}that S is in a particular energy state "E"_{"m"}.Suppose S is in a microstate indexed by "m". From the above definitions, the total energy of the system S* is given by

:$E^ast\; =\; E\text{'}\; +\; E\_m\; ,$

Notice "E"* is constant, since the combined system S* is taken to be isolated.

Now, arguably the key step in the derivation is that "the probability of "S" being in the m-th state, $;\; p\_m$, is proportional to the corresponding number of microstates available to the reservoir when "S" is in the m-th state". Therefore,

:$p\_m\; =\; C\text{'}Omega\text{'}(E\text{'})\; ,$

for some constant $;\; C\text{'}$. Taking the logarithm gives

:$ln\; p\_m\; =\; ln\; C\text{'}\; +\; ln\; Omega\text{'}\; (E\text{'})\; =\; ln\; C\text{'}\; +\; ln\; Omega\text{'}\; (E^*\; -\; E\_m)\; ,$

Since "E"

_{"m"}is small compared to "E"*, aTaylor series expansion can be performed on the latter logarithm around the energy "E"*. A good approximation can be obtained by keeping the first two terms of theTaylor series expansion::$ln\; Omega\text{'}(E\text{'})\; =\; sum\_\{k=0\}^infty\; frac\{(E\text{'}\; -\; E^ast\; )^k\; \}\{k!\}\; frac\{d^k\; ln\; Omega\text{'}\; (E^ast)\}\{dE\text{'}^k\}\; approx\; ln\; Omega\text{'}(E^ast)\; -\; frac\{d\}\{dE\text{'}\}\; ln\; Omega\text{'}(E^ast)\; E\_m$

The following quantity is a constant which is traditionally denoted by "β", known as the

.thermodynamic beta :$eta\; =\; frac\{d\}\{dE\text{'}\}\; ln\; Omega\text{'}(E^ast)\; =\; left\; .\; frac\{d\}\{dE\text{'}\}\; ln\; Omega\text{'}(E\text{'})\; ight\; |\_\{E\text{'}=E^ast\}$

Finally,:$ln\; p\_m\; =\; ln\; C\text{'}\; +\; ln\; Omega\text{'}(E^ast)\; -\; eta\; E\_m\; ,$

Exponentiating this expression gives:$p\_m\; =\; C\text{'}\; Omega\text{'}(E^ast)\; e^\{-eta\; E\_m\}$

The factor in front of the exponential can be treated as a normalization constant "C", where:$C\; =\; C\text{'}\; Omega\text{'}(E^ast)\; ,$

From this:$p\_m\; =\; C\; e^\{-eta\; E\_m\}\; ,$

**Normalization to recover the partition function**Since probabilities must sum to 1, it must be the case that

:$sum\_m\; p\_m\; =\; 1\; =\; sum\_m\; C\; e^\{-eta\; E\_m\}\; =\; C\; sum\_m\; e^\{-eta\; E\_m\}\; iff\; C\; =\; frac\{1\}\{sum\_m\; e^\{-eta\; E\_m\; equiv\; frac\{1\}\{Z(eta)\}$

where $Z$ is known as the Partition function for the canonical ensemble.

**Note on derivation**As mentioned above, the derivation hinges on recognizing that the probability of the system being in a particular state is proportional to the corresponding multiplicities of the reservoir (the same can be said for the grand canonical ensemble). As long as one makes that observation, it is flexible as how one might proceed. In the derivation given, the logarithm is taken, then a linear approximation based on physical arguments is used. Alternatively, one can apply the thermodynamic identity for differential

entropy ::$d\; S\; =\; \{1\; over\; T\}\; (d\; U\; +\; P\; d\; V\; -\; mu\; d\; N)$

and obtain the same result. See the article on

Maxwell-Boltzmann statistics where this approach is employed.The canonical ensemble is also called the

**Gibbs ensemble**, in honor of J.W. Gibbs, widely regarded with Boltzmann as being one of the two fathers of statistical mechanics. In his definitive original book "Elementary Principles in Statistical Mechanics", Gibbs viewed an ensemble as a list of the allowed states of the system (each state appearing once and only once in the list) and the associated statistical weights. The states do not interact with each other, or with a reservoir, until Gibbs treats what happens when two complete ensembles at two different temperatures are allowed to interact weakly (Gibbs, pp 160). Gibbs writes that "...the distribution in phase..." (the phase space density in modern language) "... [is] called canonical... [if] the index of probability" (the logarithm of the statistical weight of the phase space density) "...is a linear function of the energy..." (Gibbs, Ch. 4). In Gibbs' formulation, this requirement (his equation 91, in modern notation:$P\; =\; e^\{frac\{E-A\}\{kT\}\; \}\; ,$

is taken to "define" the canonical ensemble and to "be" the fundamental postulate. Gibbs does show that a large collection of interacting microcanonical systems approaches the canonical ensemble, but this is part of his demonstration (Gibbs, pp 169-183) that the principle of equal a priori probabilities, therefore the microcanonical ensemble, are inferior to the canonical ensemble as an axiomatization of statistical mechanics, at every point where the two treatments differ.

Gibbs original formulation is still standard in modern mathematically rigorous treatments of statistical mechanics, where the canonical ensemble is defined as the

probability measure :$e^\{\; \{E\; -\; A\; over\; kT\}\; \}\; dp\; ,\; dq$with "p" and "q" being the canonical coordinates.**Characteristic state function**The

characteristic state function of the canonical ensemble is theHelmholtz free energy function, as the following relationship holds::$Z(T,V,N)\; =\; e^\{-\; eta\; A\}\; ,;$

**Quantum mechanical systems**By applying the canonical partition function, one can easily obtain the corresponding results for a canonical ensemble of quantum mechanical systems. A quantum mechanical ensemble in general is described by a

density matrix . Suppose the Hamiltonian "H" of interest is aself adjoint operator with onlydiscrete spectrum . The energy levels $\{\; E\_n\; \}$ are then the eigenvalues of "H", corresponding to eigenvector $|\; psi\; \_n\; angle$. From the same considerations as in the classical case, the probability that a system from the ensemble will be in state $|\; psi\; \_n\; angle$ is $p\_n\; =\; C\; e^\{-\; eta\; E\_n\}$, for some constant $C$. So the ensemble is described by the density matrix:$ho\; =\; sum\; p\_n\; |\; psi\; \_n\; angle\; langle\; psi\_n\; |\; =\; sum\; C\; e^\{-\; eta\; E\_n\}\; |\; psi\; \_n\; angle\; langle\; psi\_n$

(Technical note: a density matrix must be

trace-class , therefore we have also assumed that the sequence of energy eigenvalues diverges sufficiently fast.) A density operator is assumed to have trace 1, so:$operatorname\{Tr\}\; (\; ho)\; =\; Q\; =\; sum\; C\; e^\{-\; eta\; E\_n\}\; =\; 1$

, which means

:$C\; =\; frac\{1\}\{sum\; e^\{-\; eta\; E\_n\}\; \}\; =\; frac\{1\}\{Q\}.$

"Q" is the quantum-mechanical version of the canonical partition function. Putting "C" back into the equation for ρ gives

:$ho\; =\; frac\{1\}\{sum\; e^\{-\; eta\; E\_n\; sum\; e^\{-\; eta\; E\_n\}\; |\; psi\; \_n\; angle\; langle\; psi\_n|\; =frac\{1\}\{\; operatorname\{Tr\}(\; e^\{-\; eta\; H\}\; )\; \}\; e^\{-\; eta\; H\}\; .$

By the assumption that the energy eigenvalues diverge, the Hamiltonian "H" is an

unbounded operator , therefore we have invoked theBorel functional calculus to exponentiate the Hamiltonian "H". Alternatively, in non-rigorous fashion, one can consider that to be the exponential power series.Notice the quantity

:$operatorname\{Tr\}(\; e^\{-\; eta\; H\}\; )$

is the quantum mechanical counterpart of the canonical partition function, being the normalization factor for the mixed state of interest.

The density operator ρ obtained above therefore describes the (mixed) state of a canonical ensemble of quantum mechanical systems. As with any density operator, if "A" is a physical observable, then its expected value is

:$langle\; A\; angle\; =\; operatorname\{Tr\}(\; ho\; A\; ).$

**Relations with other ensembles**A generalization of this is the

grand canonical ensemble , in which the systems may share particles as well as energy. By contrast, in themicrocanonical ensemble , the energy of each individual system is fixed.**ee also****References****External links*** [

*http://www.iop.org/EJ/article/0295-5075/54/4/475/6545.html The Ornstein-Zernike equation in the canonical ensemble*]* [

*http://stp.clarku.edu/notes/ Draft Chapters of Thermal and Statistical Physics Textbook*]

*Wikimedia Foundation.
2010.*

### Look at other dictionaries:

**canonical ensemble**— ▪ physics in physics, a functional relationship for a system of particles that is useful for calculating the overall statistical and thermodynamic behaviour of the system without explicit reference to the detailed behaviour of particles.… … Universalium**canonical ensemble**— kanoninis ansamblis statusas T sritis fizika atitikmenys: angl. canonical assembly; canonical ensemble vok. kanonisches Ensemble, n rus. канонический ансамбль, m pranc. ensemble canonique, m … Fizikos terminų žodynas**Grand canonical ensemble**— In statistical mechanics, the grand canonical ensemble is a statistical ensemble (a large collection of identically prepared systems), where each system is in equilibrium with an external reservoir with respect to both particle and energy… … Wikipedia**grand canonical ensemble**— didysis kanoninis ansamblis statusas T sritis fizika atitikmenys: angl. grand canonical ensemble vok. große Gibbssche Gesamtheit, f; große kanonische Gesamtheit, f; großes kanonisches Ensemble, n rus. большой канонический ансамбль, m pranc.… … Fizikos terminų žodynas**Canonical**— is an adjective derived from . Canon comes from the Greek word kanon , rule (perhaps originally from kanna reed , cognate to cane ), and is used in various meanings. Basic, canonic, canonical : reduced to the simplest and most significant form… … Wikipedia**canonical assembly**— kanoninis ansamblis statusas T sritis fizika atitikmenys: angl. canonical assembly; canonical ensemble vok. kanonisches Ensemble, n rus. канонический ансамбль, m pranc. ensemble canonique, m … Fizikos terminų žodynas**ensemble canonique**— kanoninis ansamblis statusas T sritis fizika atitikmenys: angl. canonical assembly; canonical ensemble vok. kanonisches Ensemble, n rus. канонический ансамбль, m pranc. ensemble canonique, m … Fizikos terminų žodynas**Ensemble average**— In statistical mechanics, the ensemble average is defined as the mean of a quantity that is a function of the micro state of a system (the ensemble of possible states), according to the distribution of the system on its micro states in this… … Wikipedia**Canonical probability distribution**— In thermal physics, the canonical probability distribution is a statistical function which equates to the Boltzmann factor divided by the partition function. The function was introduced by Willard Gibbs in his 1901 Elementary Principles in… … Wikipedia**ensemble grand-canonique**— didysis kanoninis ansamblis statusas T sritis fizika atitikmenys: angl. grand canonical ensemble vok. große Gibbssche Gesamtheit, f; große kanonische Gesamtheit, f; großes kanonisches Ensemble, n rus. большой канонический ансамбль, m pranc.… … Fizikos terminų žodynas