There is real statistical mechanics today. An ensemble is defined and an average calculated.
Microcanonical & Canonical Ensembles
It has been demonstrated that there is a huge number of microstates. It is possible to connect thermodynamics with the complexity of the microscopic world. Below are definitions.
Microstate
A microstate is a particular state of a system specified at the atomic level. This could be described by the many-body wavefunction. A system is something that over time fluctuates between different microstates. An example includes the gas from last time. There is immense complexity. Fix the variables <math>T, V,</math> and <math>N</math>. With only these variables known, there is no idea what microstate the system is in at the many body wavefunction level.
Consider a solid. It is possible to access different configurational states, and two are shown below. Diffusion is a process that results in changes in state over time. X-ray images may not be clear due to the diffusion of atoms and systems accessing different microstates. There may be vibrational excitations, and electronic excitations correspond to the excitation of electrons to different levels. These excitations specify what state the solid is in. Any combination of excitations specify the microstate of a system.
Unable to render embedded object: File (Two_configurational_states.PNG) not found.
Summary
A particular state of a system specified at the atomic level (many function wavebody level
[ \Psi_{\mbox{manybody}} ]
. The system over time fluctuates betwen different microstates
- There is immense complexity in a gas
- Solid
- Configurational states
- Vibrational excitations
- electronic excitations
- There is usually a combination of these excitations, which results in immense complexity
- Excitations specify the microstate of the system
Why Ensembles?
A goal is to find
[ P_v ]
. Thermodynamic variables are time averages. Sum over the state using the schrodinger equation to find energy and multiplying by probability. To facilitate averages, ensembles are introduced. Ensembles are collections of systems. Each is very large, and they are macroscopically identical. Look at the whole and see what states the system could be in. Below is a diagram of systems and ensembles. Each box represents a system, and the collection of systems is an ensemble. Each box could represent the class, and v could represent the sleep state. Look at the properties of the ensemble. There are
[ A ]
;macroscopically identical systems. Eventually
[ A ]
is taken to go to
[ \infty ]
. Each system evolves over time.
Unable to render embedded object: File (Ensemble_of_systems.PNG) not found.
Deriving P_v
The probability,
[ P_v ]
is defined for different kinds of boundary conditions. When looking at the probability that students in a class are asleep, it is possible to take an instantaneous snapshot. The term
[ a_v ]
is the occupation number and is defined as the number of systems that are in the state
[ v ]
at the time of the snapshot. The fraction of systems in state
[ v ]
is
[ P_v ]
. This is one approximation to get the probability, and it could be a bad approximation.
[ P_v \approx \frac
]
There are an additional definitions of
[ P_v ]
. It is equal to the probability of finding a system in state v at time t or identically it is the fraction of time spent by the system in state v. In the case of the sleep example, it is equal to the fraction of time that any one is in the sleep state.
The time average corresponds to looking at a class and seeing for what fraction of time does it find someone in the class asleep. The ensemble average corresponds to looking at the set of identical classes and seeing how many classes have at least one student asleep. There is a correlation between the state average and the time average. There is a need of boundary conditions. Take a picture of a large number of systems, look at everyone, and average.
Summary
Recap:
[ E = \sum_V E_V P_V ]
To facilitate averages, we introduce "ensembles" that we average over
- Averaging over many bodies rather than averaging over time
- Example: student = system, v =sleepstate
Ensemble of systems:
- 'A' (a very large number) macroscopically identical systems
- Each system evolves over time
Probability:
- Take an instantaneous snapshot
- Define
Unknown macro: {latex}= #of systems that are in state v at the time of snapshot
[ a_v ]
- Fraction of systems in state v is
Unknown macro: {latex}
[ \frac
Unknown macro: {a_v}Unknown macro: {A}\simeq P_v ]
- Probability to find a system in statevat timet
- Fraction of time spent in statev
Microcanonical Ensemble
The boundary conditions of all systems of the microcanonical ensemble are the same. The variablesN, V,andEcannot fluctuate. Each system can only fluctuate between states with fixed energy,E.
Unable to render embedded object: File (Microcanonical_ensemble.PNG) not found.
It is possible to get degeneracy from the Shrodinger equation.
[ \hat H \Psi = E \begin
\underbrace{ (\Psi_1, \Psi_2, ..., \Psi_
) }
\Omega(E) \end
]
Consider the example of the hydrogen atom. Below is an expression of the energy proportionality and the degeneracy when
[ n=1 ]
and
[ n=2 ]
<P> </P>
[ E \alpha \frac
]
<P> </P>
[ \Omega (n=1) = 1 ]
<P> </P>
[ \Omega (n=2) = 4 ]
What is Pv in microcanonical ensemble?
All states should be equally probable with variables
[ N, V, ]
and
[ E ]
[ fixed. The term
[ P_v ]
is the probability of being in any
[ E ]
state for a system, and it should be equal to a constant. An expression is below. Each state can be accessed, and one is not more favored. This is related to the principle of a priori probability. There is no information that states should be accessed with different probability.
[ P_v = \frac
]
Example
Consider an example of a box and gas. All the atoms are in one corner in the second box. Add to get complete degeneracy. The value of
[ \Omega_1 (E) ]
is large; there is enormous degeneracy.
Unable to render embedded object: File (Microcanonical_ensemble_II.PNG) not found.
[ \Omega_1(E) \gg \Omega_2(E) ]
Consider a poker hand. There is a lot of equivalence in bad hands. These are dealt most of the time and correspond to
[ \Omega_1(E) ]
. The royal flush corresponds to
[ \Omega_2(E) ]
. It is equally probable, but there are many fewer ways to get the royal flush. There are the same boundary conditions. In an isolated system, in which
[ N, V, ]
and
[ E ]
are fixed, it is equally probable to be in any of its
[ \Omega (E) ]
possible quantum states.
Summary
The variables
[ N, V, ]
and
[ E ]
are fixed
- Each system can only fluctuate between states with fixed energy
Unknown macro: {latex}(like from Schrodinger's equation)
[ E ]
- Unknown macro: {latex}
[ \hat H \Psi = E \Psi \rightarrow E[\Psi_2 .... \Psi_\omega \] ]
- Unknown macro: {latex}
[ \Omega(E) ]
- All states are equally probable, and are given equal weight.
Hydrogen atom
- Unknown macro: {latex}
[ E = \frac
Unknown macro: {1}Unknown macro: {n^2}]
- Unknown macro: {latex}
[ \Omega (n=1) =1 ]
- Unknown macro: {latex}
[ \Omega (n=2) =4 ]
Probability of being in any E state for a system
- employed the principle of equal a priory probabilities
- Unknown macro: {latex}
[ P_v = \mbox
Unknown macro: {constant}]
- Unknown macro: {latex}
[ P_v = \frac
Unknown macro: {1}Unknown macro: {Omega(E)}]
Systems
- Equally more probable
- Some accessed more times because of large degeneracy number
- There's just a lot more ways to get configuration left (like a craphand) than the right (like a straight flush)
- Unknown macro: {latex}
[ \Omega_1 (E) >> \Omega_2 (E) ]
An isolated system (
[ N, V, E= ]
fixed) is equally probable to be in any of its
[ \Omega (E) ]
possible quantum states.
Canonical Ensemble
There is a different set of boundary conditions in the canonical ensemble. There are heat conducting walls or boundaries of each system. Each of the
[ A ]
members of the ensemble find themselves in the heat bath formed by the
[ A-1 ]
members. Each system can fluctuate between different microstates. An energy far from average is unlikely. In the picture below, the energy of the ensemble on the right side of the ensemble is fixed, while the energy of a particular system is not fixed and can fluctuate.
Unable to render embedded object: File (Canonical_ensemble.PNG) not found.
Take another snapshot. There is interest in the distribution. The term
[ \overline
]
is equal to the number of systems in state
[ v ]
.
[
= \overline
]
Below is a table of microstates, energy, and occurence, and a graph. In the graph, equilibrium has occurred, but all states can be accessed. It is possible to access different states some distance from the average energy. The total energy,
[ \epsilon ]
, is fixed and is equal to the integral of the curve. As the number of systems increases, the curve becomes sharper.
| Unknown macro: {latex} [ \mbox Unknown macro: {microstate}
] |
1 |
2 |
3 |
Unknown macro: {latex} [ \nu ] |
|---|---|---|---|---|
| Unknown macro: {latex} [ \mbox Unknown macro: {energy}
] |
Unknown macro: {latex} [ E_1 ] |
Unknown macro: {latex} [ E_2 ] |
Unknown macro: {latex} [ E_3 ] |
Unknown macro: {latex} [ E_ Unknown macro: {nu}
] |
| Unknown macro: {latex} [ \mbox Unknown macro: {occurence}
] |
Unknown macro: {latex} [ a_1 ] |
Unknown macro: {latex} [ a_2 ] |
Unknown macro: {latex} [ a_3 ] |
Unknown macro: {latex} [ a_ Unknown macro: {nu}
] |
Unable to render embedded object: File (Occurence_versus_energy.PNG) not found.
Constraints
Below are constraints. The first is the sum of the occupation number. The second constraint is possible due to the system being isolated.
[ \sum_v a_v = A ]
<P> </P>
[ \sum_v a_v E_v = \epsilon ]
The term
[ P_v ]
is the probability of finding the system in state
[ v ]
. It is possible to use snapshot probability. There are many distributions that satisfy the boundary conditions. There is a better way to find
[ P_v ]
, and a relation is below. It corresponds to the average distribution. This is associated with a crucial insight.
[ P_v = \frac{\overline{a_v}}
]
Crucial Insight
An assumption is that the entire canonical ensemble is isolated. No energy can escape, and the energy
[ \epsilon ]
is constant. Every distribution of
[ \overline
]
that satisfies the boundary conditions is equally probable. it is possible to write many body wavefunction because the energy of the entire ensemble is fixed. The principle of equal a priori
probabilities is applied. Look at the whole distribution that satisfies the boundarty condition. Each distribution of occurance numbers must be given equal weights.
Summary
The variables
[ N, V, T ]
are fixed
- There are heat conducting bondaries of each system
- Each of the
Unknown macro: {latex}(= large number) members finds itself in a heat bath, formed by the
[ A ]
Unknown macro: {latex}other members[ (A - 1) ]
- Take snapshot; get distribution
Unknown macro: {latex}(= # of systems in state
[ a_v = \overline
Unknown macro: {a}]
Unknown macro: {latex})[ v ]
Constraints
The total energy,
[ \epsilon ]
, is fixed
[ \sum_v a_v = A \sum_v a_v E_v = \epsilon ]
(isolated!)
Probability
[ P_v \simeq \frac
]
is an approximation.
- Better to use
Unknown macro: {latex}, the averaged distribution
[ P_v = \frac{\overline{a_v}}
Unknown macro: {A}]
- There is an assumption that the whole canonical ensemble is isolated and that energy
Unknown macro: {latex}is constant. Every distribution of
[ \epsilon ]
Unknown macro: {latex}that satisfies the boundary conditions is equally probable.[ \overline
Unknown macro: {a}]
- We are applying the principle of equal a priori probabilities, and each distribution of occurance numbers must be given equal weights.
Some Math
Consider every possible distribution consistent with boundary conditions, and for each distribution consider every possible permutation. The term
[ w (\overline
) ]
is equal to the number of ways to obtain a distribution
[ \overline
]
, where
[ a_v ]
is the number of systems in state
[ v ]
. A bad hand in poker is defined by a large number of
[ a ]
. Use the multinomial distribution.
[ a_i = \overline
]
<P> </P>
[ w (\overline
) = \frac
]
<P> </P>
[ w (\overline
) = \frac
]
<P> </P>
[ a_1 ]
<P> </P>
[ \mbox
]
<P> </P>
[ a_
]
<P> </P>
[ \mbox
]
Below are expressions of the probability to be in a certain state. The term
[ a_v ]
is averaged over all possible distributions. Every distribution is given equal weight, and the one with the most permutations is the most favored.
[ P_v = \frac{\overline{a_v}}
]
<P> </P>
[ P_v = \frac
\frac{ \sum_{\overline {a}} \omega (\overline
) a_v (\overline
) }{ \sum_
\omega (\overline
) } ]
Example
Below is an example of four systems in an ensemble. The term
[ P_1 ]
is the probability of any system to be in state
[ 1 ]
.
| Unknown macro: {latex} [ \mbox Unknown macro: {state}
] |
Unknown macro: {latex} [ 1 ] |
Unknown macro: {latex} [ 2 ] |
|---|---|---|
| Unknown macro: {latex} [ \mbox Unknown macro: {energy}
] |
Unknown macro: {latex} [ E_1 ] |
Unknown macro: {latex} [ E_2 ] |
| Unknown macro: {latex} [ \mbox Unknown macro: {occupation}
] |
Unknown macro: {latex} [ a_1 ] |
Unknown macro: {latex} [ a_2 ] |
| Unknown macro: {latex} [ \mbox Unknown macro: {distribution}
] |
Unknown macro: {latex} [ a_0 ] |
Unknown macro: {latex} [a_1 ] |
[ \mbox
]
<P> </P>
[ w(A) = \frac
]
<P> </P>
[ w(A) = 1 ]
<P> </P>
[ \mbox
]
<P> </P>
[ w(B) = \frac
]
<P> </P>
[ w(B) = 4P_1 = \frac
\left ( \frac
\right ) ]
<P> </P>
[ P_1 = \frac
]
<P> </P>
Distribution of
[ w(a) ]
The term
[ w(a) ]
is the number of permutations for a particular distribution. As the number of systems increases, or asAincreases, the distribution becomes more peaked.
Consider the probability.
[ P_
= \frac
\frac{ \sum_{\overline {a}} \omega (\overline
) a_
(\overline
) }{ \sum_{\overline {a}} \omega (\overline
) } ]
[ P_
\approx \frac{ \frac
w ( \overline
^*) a_
^*}{w ( \overline
^*) } ]
[ P_
\approx \frac{a_
^*}
]
(Equation 1)
Look at the distribution that maximizes
[ w ( \overline
) ]
, the permutation number. To get
[ a_
]
, maximize
[ w ( \overline
) ]
subject to the constraints below.
[ \sum_
a_
- A = 0 ]
[ \sum_
a_
E_
- \epsilon = 0 ]
Use Lagrange multipliers, and maximize
[ \ln w ( \overline
) ]
in order to be able to use Stirling's approximation.
[ \frac
{\partial a_{\nu}} \left ( \ln w ( \overline
) - \alpha \sum_k a_k - \beta \sum_k a_k E_k \right ) = 0 ]
[ w ( \overline
) = \frac
]
[ \ln w ( \overline
) = \ln A! + \left ( - \ln \pi_k a_k! \right ) = \ln A! - \sum_k \ln a_k! ]
Use Stirling's approximation as
[ A ]
and the occupation number,
[ a_k ]
, go to infinity.
[ \sum_k \ln a_k! = \sum_k \left( a_k \ln a_k - a_k \right ) ]
[ = \sum_K a_K \ln a_K - \sum_K a_K ]
[ = \sum_K a_K \ln a_K - A\frac
{\partial a_{\nu}} \left ( \ln A! - \sum_k a_k \ln a_k + A - \alpha \sum_k a_k - \beta \sum_k a_k E_k \right ) = 0 ]
[ \left ( a_v \to x \mbox
\ln A! - \sum_k a_k \ln a_k + A \to -a_v \ln a_v \mbox
\alpha \sum_k a_k \to \alpha a_v \mbox
\beta \sum_k a_k E_k \to \beta a_v E_k \right )\ln a_v - 1 - \alpha - \beta E_
= 0 ]
[
The term
[ a_
]
is the occupation number that maximizes the expressione
[ \alpha ' \cdot e^{\beta E_{\nu}} ]
, where
[ \alpha ' = \alpha + 1 ]
. Use constraints to determine the Lagrange multipliers and determine the probability.
[ a_
=e^
\cdot e^{\beta E_
]
[ \sum_
a_
= A ]
[ \sum_
e^
\cdot e^{\beta E_
= A ]
[ e^
= \frac
\sum_
e^{\beta E_{\nu}} ]
The probability of being in a certain state
[ \nu ]
can be calculated, and it is still in terms of the second Lagrange multiplier. Plugging back into Equation 1:
[ P_v = \frac{a_{\nu}}
= \frac
{\sum_
e^{-\beta E_
}} \cdot \frac{e{\beta E_
}}
= \frac{e^{-\beta E_
}}{\sum_
e^{-\beta E_{\nu}} ]
Partition Function
The denominator is the partition function,
[ Q = \sum_
e^{-\beta E_
]
. It tells us how many states are accessible by the system. Determine
[ \beta ]
, a measure of thermally accessible states 
. Look at how the partition function connects to macroscopic thermodynamic variables. Find
[ \beta ]
and find
[ \overline
]
[ \overline
= \sum_
P_
E_
\overline
= \frac{\sum_
E_
e^{-\beta E_
}}
]
Consider the average pressure. The pressure for one microstate is
[ p_
]
.
[ p_
= \frac{-\partial E_{\nu}}
]
[ p_
= \sum_
P_
p_
]
[ p_
= \frac{ -\sum_
\left ( \frac{\partial E_{\nu}}
\right ) e^{-\beta E_
}}
]
Summary
In the case of a canonical ensemble, the energy of the entire ensemble is fixed. Each state is equally probable, and there is degeneracy. The probability is a function of how many ways to get the distribution. The distribution with the most permutation is the most probable. The graph can become very peaked. Once the distribution is known, do a maximization of
[ w(\overline
) ]
. Use Lagrange multipliers and two constraints. The term
[ a_
^* ]
is the distribution that maximizes an expression. This is what is most often found. Go back to the probability, get an expression, and give part of it a name. The term
[ \beta ]
is a measure of how many states are thermally accessible.
[ \overline
]
Consider every possible distribution
[ a_i = \overline
]
(consistent with boundary conditions)
- For each distribution, consider every possible permutation
- The number of ways to obtain
Unknown macro: {latex}is
[ \overline
Unknown macro: {a}]
Unknown macro: {latex}, where[ \omega (\overline
Unknown macro: {a}) = \frac
Unknown macro: { A! }Unknown macro: { a_1! a_2! a_3! ..... a_v! }= \frac
Unknown macro: { Pi_v a_v!}]
Unknown macro: {latex}is the number of systems in state[ a_v ]
Unknown macro: {latex}.[ v ]
Probability
[ P_v = \frac{\overline{a_v}}
= \frac
\frac{ \sum_
\omega (\overline
) a_v (\overline
) }{ \sum_{overline a}}\omega (\overline
) } ]
Averaging
[ a_v ]
over all possible distributions.
[ w(\overline a) ]
[ w(\overline a) ]
is very peaked around a specific distribution
Increase
[ A ]
and
[ \omega
]
becomes more peaked
[ a_v ]
To get
[a_v ]
, maximize
[ w (\overline a) ]
subject to constraints.
Find the partition function,
[ Q ]