\[ \Psi_{\mbox{manybody}} \] 

 \[ P_v \] 

 \[ A \] 

 \[ A \] 

 \[ \infty \] 

 \[ P_v \] 

 \[ a_v \] 

 \[ v \] 

 \[ v \] 

 \[ P_v \] 

 \[ P_v \approx \frac {a_v}{A} \] 

 \[ P_v \] 

 \[ E

h1.  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.

h1. 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.

<p>
</p>

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.

<br>

<center>
!Two_configurational_states.PNG!
</center>

<br>

h2. Summary

<p>
</p>

A particular state of a system specified at the atomic level (many function wavebody level <math>\Psi_\mbox{manybody}</math> )

* 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

h1. Why Ensembles?

<p>
</p>

A goal is to find <math>P_v</math>.  Thermodynamic variables are time averages.  Sum over the state using the schrodinger equation to find energy and multiplying by probability.

<p>
</p>

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 <math>A</math> macroscopically identical systems.  Eventually <math>A</math> is taken to go to <math>\infty</math>.  Each system evolves over time.

<center>

<br>

!Ensemble_of_systems.PNG!

<br>

</center>

h3. Deriving <math>P_v</math>

The probability, <math>P_v</math> 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 <math>a_v</math> is the occupation number and is defined as the number of systems that are in the state <math>v</math> at the time of the snapshot.  The fraction of systems in state <math>v</math> is <math>P_v</math>.  This is one approximation to get the probability, and it could be a bad approximation.

<center>

<br>

<math>P_v \approx \frac{a_v}{A}</math>

<br>

</center>

There are an additional definitions of <math>P_v</math>.  It is equal to the probability of finding a system in state <math>v</math> at time <math>t</math> or identically it is the fraction of time spent by the system in state <math>v</math>.  In the case of the sleep example, it is equal to the fraction of time that any one is in the sleep state.  

<p>
</p>

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.

<p></p>

h2. Summary

<p>
</p>

Recap: 
<center>

<br>

<math>E = \sum_V E_V P_V</math>

<br>

</center>

To facilitate averages, we introduce V \] 

 \[ v</math>
** Probability to find a system in state <math>v</math> at time <math>t</math>
** Fraction of time spent in state <math>v</math>

h1. Microcanonical Ensemble

The boundary conditions of all systems of the microcanonical ensemble are the same.  The variables <math>N, V,</math> and <math>E</math> cannot fluctuate.  Each system can only fluctuate between states with fixed energy, <math>E</math>.

<center>

<br>

!Microcanonical_ensemble.PNG!

<br>

</center>

It is possible to get degeneracy from the Shrodinger equation.

<center>

<br>

<math>\hat H \Psi = E \begin {matrix} \underbrace{ (\Psi_1, \Psi_2, ..., \Psi_{\Omega}) } \\ \Omega(E) \end{matrix}</math>

<br>

</center>

Consider the example of the hydrogen atom.  Below is an expression of the energy proportionality and the degeneracy when <math>n=1</math> and <math>n=2</math>

<center>

<br>

<math>E \prop \frac \] 

 \[ n=1 \] 

 \[ n=2 \] 

 \[ N, V, \] 

 \[ E \] 

 \[ P_v \] 

 \[ E \] 

 \[ P(n=1) = 1</math>

<br>

<math>\Omega (n=2) = 4</math>

<br>

</center>

h2. What is <math>P_v</math> in microcanonical ensemble?

All states should be equally probable with variables <math>N, V,</math> and <math>E</math> fixed.  The term <math>P_v</math> is the probability of being in any <math>E</math> 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.

<center>

<br>

<math>P_v = \frac{1}{\Omega (E)} \] 

 \[ </math>

<br>

</center>

h2. 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 <math>\Omega_1 (E)</math> is  \] 

 \[  degeneracy.

<center>

<br>

!Microcanonical_ensemble_II.PNG!

<br>

<math>\Omega_1(E) \gg \Omega_2(E)</math>

<br>

</center>

Consider a poker hand. There is a lot of equivalence in bad hands.  These are  \] 

 \[ dealt most of the time and correspond to <math>\Omega_1(E)</math>. \] 

 \[  to <math>\Omega_2(E)</math>. \] 

 \[ N which <math>N, V,</math> and <math>E</math> are fixed, it is equally probable to be in any of \] 

 \[ E \] 

 \[  its <math>\Omega (E)</math> possible quantum states.

<p>
</p>

h2. Summary

<p>
</p>

The variables <math>N, V,</math> and <math>E</math> are fixed
* Each system can only fluctuate between states with fixed energy E (like from Schr��dinger's equation)
* <math> \] 

 \[ N, V, \] 

 \[ E \] 

 \[ N</math>

<p>
</p>

An isolated system (<math>N, V, E</math> = E= \] 

 \[  its <math>\Omega (E)</math> possible quantum states.

h1.  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 <math>A</math> members of the ensemble find themselves in the heat bath formed by the <math>A-1</math> 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.

<center>

<br>

!Canonical_ensemble.PNG!

<br>

</center>

Take another snapshot.  There is interest in the distribution.  The term <math>\overline{a}</math> is equal to the number of systems in state <math>v</math>.

<center>

<br>

<math> \] 

 \[ A \] 

 \[ A-1 \] 

 \[ \overline {a} \] 

 \[ v \] 

 \[ {a_v} = \overline{a} </math>

<br>

</center>

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, <math>\epsilon</math>, is fixed and is equal to the integral of the curve.  As the number of systems increases, the curve becomes sharper.

<center>

<br>

<table padding = 5>

<tr>

<td>
<math>\mbox{microstate}</math>
</td>

<td>
<math>1</math>
</td>

<td>
<math>2</math>
</td>

<td>
<math>3</math>
</td>

<td>
<math>\nu</math>
</td>

</tr>

<tr>

<td>
<math>\mbox{energy}</math>
</td>

<td>
<math>E_1</math>
</td>

<td>
<math>E_2</math>
</td>

<td>
<math>E_3</math>
</td>

<td>
<math>E_{\nu}</math>
</td>

</tr>

<tr>

<td>
<math>\mbox{occurence}</math>
</td>

<td>
<math>a_1</math>
</td>
 \] 

 \[ \epsilon \] 

 \[ 
<td>
<math>a_2</math>
</td>

<td>
<math>a_3</math>
</td>

<td>
<math>a_{\nu}</math>
</td>

</tr>

</table>

<br>

!Occurence_versus_energy.PNG!

<br>

</center>

h2. Constraints

Below are constraints.  The first is the sum of the occupation number.  The second constraint is possible due to the system being isolated.  

<center>

<br>

<math>\sum_v a_v = A</math>

<br>

<math>A \] 

 \[ P_v \] 

 \[ v \] 

 \[ P_v \] 

 \[ P \epsilon</math>

<br>

</center>

The term <math>P_v</math> is the probability of finding the system in state <math>v</math>.  It is possible to use snapshot probability.  There are many distributions that satisfy the boundary conditions.  There is a better way to find <math>P_v</math>, and a relation is below.  It corresponds to the average distribution.  This is associated with a crucial insight.

<center>

<br>

<math>P_v = \frac{\overline{a_v}}{A}</math>

<br>

</center>

h2. Crucial Insight

An assumption is that the entire canonical ensemble is isolated.  No energy can escape, and the energy <math>\epsilon</math> is constant. Every distribution of <math>\overline{a}</math> 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|http://en.wikipedia.org/wiki/A_priori_%28statistics%29] probabilities is applied.  Look at the whole distribution that satisfies the boundarty condition.  Each distribution of occurance numbers must be given equal weights.

<p>
</p>

h2. Summary

<p>
</p>

The variables <math>N, V, T</math> are fixed
* There are heat conducting bondaries of each system
* Each of the <math>A</math> (= large number) members finds itself in a heat bath, formed by the <math>(A - 1)</math> other members
* Take snapshot; get distribution <math>\{ a_v \} = \overline{a}</math> (= # of systems in state v)

<p>
</p>

Constraints
* The total  \] 

 \[ \epsilon \] 

 \[ \overline{a} \] 

 \[ N, V, T \] 

 \[ \epsilon \] 

 \[ energy, <math>\epsilon</math>, is fixed
* <math>\sum_v a_v = A</math>
*A <math>\sum_v a_v E_v = \epsilon \epsilon</math>] 

 \[ P

Probability
* <math>P_v \simeq \frac{a_v}{A}</math> is an approximation. 
* Better to use <math>P \] 

 \[ P_1 \] 

 \[ 1 \] 

 \[ \mbox{Distribution A} \] 

 \[ w(a) \] 

 \[ P</math>

<br>

<math>P_1 = \frac{1}{2}</math>

</center>

h2. Distribution of w(a)

<p>
</p>

The term <math>w(a)</math> is the number of permutations for a particular distribution.  As the number of systems increases, or as <math>A</math> increases, the distribution becomes more peaked.

<center>

<br>

!W%28a%29_versus_a.PNG!
!W%28a%29_versus_a_--_large_A.PNG!

<br>

</center>

Consider the probability.

<center>

<br>

<math>P_{\nu} = \frac{1}{A} \frac{  \sum_{\overline {a}} \omega (\overline{a}) a_{\nu} (\overline{a}  ) }{  \sum_{\overline {a}} \omega (\overline{a}) } </math>

<br>
 \] 

 \[ P
<math>P_{\nu} \approx \frac{ \frac{1}{A} w ( \overline{a}^*) a_{\nu}^*}{w ( \overline{a}^*) } </math>

<br>

<math>P \] 

 \[ P_{\nu} \approx \frac{a_{\nu}^*}{A} \] 

 \[ w</math>                          (Equation 1)

<br>

</center>

Look at the distribution that maximizes <math>w ( \overline{a} )</math>, the permutation number.  To get <math>a) \] 

 \[ a_{\nu} \] 

 \[ w</math>, maximize <math>w ( \overline{a} )</math> subject to the constraints below.

<center>

<br>

<math> \] 

 \[ \sum_{\nu} a_{\nu} - A = 0</math>

<br>

<math>0 \] 

 \[ \sum_{\nu} a_{\nu}E_{\nu} - \epsilon = 0 \] 

 \[  0</math>

<br>

</center>

Use Lagrange multipliers, and maximize <math>\ln w ( \overline{a} ) </math>\] 

 \[ 

<center>

<br>

<math>\frac{\partial}{\partial a_{\nu}} \left ( \ln w ( \overline{a} ) - \alpha \sum_k a_k - \beta \sum_k a_k E_k \right ) = 0</ 0 \] 

 \[ wmath>

<br>

<math>w ( \overline{a} ) = \frac{A!}{\pi_k a_k!} \] 

 \[ </math>

<br>

<math>\ln w ( \overline{a}) = \ln A! + \left ( - \ln \pi_k a_k! \right )  = \ln A! - \sum_k \ln a_k! </math>

<br>

</center>

Use  \] 

 \[ as A \] 

 \[ a_k \] 

 \[  <math>a_k</math>, go to infinity.

<center>

<br>

<math>\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_KK \] 

 \[ = \sum_K a_K \ln a_K - A </math>

<br>

<math>\frac{\partial}{\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</math>

<br>

<math>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 ) </math>

<br>

<math>- \ln a_v - 1 - \alpha - \beta E_{\nu} = 0 \] 

 \[ a 0</math>

<br>

</center>

The term <math>a_{\nu}^*</math> is the occupation number that maximizes the expression <math>e^{- \] 

 \[ \alpha '} \cdot e^{-\beta E_{\nu}} \] 

 \[ </math>, where <math>\alpha ' = \alpha + 1</math>.  use constraints to determine the Legrange multipliers and determine the probability.

<center>

<br>

<math>a1 \] 

 \[ a_{\nu}=e^{-\alpha '} \cdot e^{-\beta E_{\nu} \] 

 \[ </math>

<br>

<math>\sum_{\nu} a_{\nu} = A</math>

<br>

<math> A \] 

 \[ \sum_{\nu} e^{-\alpha '} \cdot e^{-\beta E_{\nu} = A</math>

<br>

<math>e^ A \] 

 \[ e^{\alpha '} = \frac{1}{A} \sum_{\nu} e^{-\beta E_{\nu}</math>

<br>

</center>

The probability of being in a certain state <math>\nu</math> can be calculated, and it is still in terms of the second Legrange multiplier. Plugging back into Equation 1:

<center>

<br>

<math>P} \] 

 \[ \nu \] 

 \[ P_v = \frac{a_{\nu}^*}{A}  = \frac{A}{\sum_{\nu} e^{-\beta E_{\nu}}} \cdot \frac{e^e{-\beta E_{\nu}}}{A}  =  \frac{e^{-\beta E_{\nu}}}{\sum_{\nu} e^{-\beta E_{\nu}} \]

 \[ Q</math>

<br>

</center>

h2. Partition Function

The denominator is the partition function, <math>Q = \sum_{\nu} e^{-\beta E_{\nu}</math>.  It tells us how many states are accessible by the system. Determine <math>\beta</math>, a measure of thermally accessible states (?).  Look at how the partition function connects to macroscopic thermodynamic variables.  Find <math>\beta</math> and find <math>\overline{E}</math>

<center>

<br>

<math> \] 

 \[ \beta \] 

 \[ \beta \] 

 \[ \overline{E} \] 

 \[ \overline{E} = \sum_{\nu} P_{\nu} E_{\nu}</math>

<br>

<math>\overline{E} = \frac{\sum_{\nu} E_{\nu} e^{-\beta E_{\nu}}}{Q} \] 

 \[ p</math>

<br>

</center>

Consider the average pressure.  The pressure for one microstate is <math>p_{\nu}</math>.

<center>

<br>

<math>p \] 

 \[ p_{\nu} = \frac{-\partial E_{\nu}}{\partial V}</math>

<br>

<math>p \] 

 \[ p_{\nu} = \sum_{\nu} P_{\nu} p_{\nu}</math>

<br>

<math>p \] 

 \[ p_{\nu} = \frac{ -\sum_{\nu} \left ( \frac{\partial E_{\nu}}{\partial V} \right ) e^{-\beta E_{\nu}}}{Q} \] 

 \[ w</math>

<br>

</center>

h2. 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 <math>w(\overline{a})</math>. \] 

 \[ a term <math>a_{\nu}^*</math> 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 <math>\beta</math> \] 

 \[ \beta \] 

 \[ \overline{a} \] 

 \[ a_i  is a measure of how many states are thermally accessible.

<br>

<math>\overline{a}</math>
* Consider every possible distribution <math>\{ a_i \} = \overline{a}</math>  \] 

 \[ P </math>, where<math> a_v</math> is the number of systems in state <math>v</math>.

<br>

Probability
* <math>P_v = \frac{\overline{a_v}}{A} = \frac{1}{A} \frac{  \sum_{\overline a} \omega (\overline{a}) a_v (\overline{a}  ) }{  \sum_{\overline a}} \omega (\overline{a}) } \] 

 \[ a_v \] 

 \[ </math>
* Averaging <math>a_v</math> over all possible distributions.

<br>

<math> \w(\overline a) \] 

</math>
* <math> \[ w(\overline a)</math> is very peaked around a specific distribution
* Increase <math>A</math> and <math> \] 

 \[ A \] 

 \[ \omega{\overline a} \] 

 \[ a_v \] 

 \[a_v \] 

 \[</math> becomes more peaked

<br>

<math>a_v^*</math>
* To get <math>a_v^*</math>, maximize<math> w (\overline a)</math> subject to constraints.
*Find the partition function, <math>Q</math>
 \] 

 \[ Q \]