3.012 Fundamentals of Materials ScienceFall 2003

Lecture 22: 11.21.03 Connecting events at the atomic/molecular level to macroscopic thermodynamic behavior

Today:

Last time......

Temperature and the occupation of states......

The Boltzmann factor and Partition function......

All thermodynamic quantities can be calculated from the partition function......

The partition function of molecules vs. multi-molecular systems......

The Einstein Solid......

The complete partition function for an Einstein solid......

References......

Reading: Dill and Bromberg Ch. 10, ‘The Boltzmann Distribution Law,’ pp. 171-191

Supplementary Reading:-

Planning Notes:

HEAT CAPACITY REFLECTS THERMAL FLUCTUATIONS- DILL P. 228

  • distribution functions that describe the ensemble average under constraints of constant energy, constant temperature
  • a common laboratory situation: constant µ,V, T
  • Description of a magnetic material using a two-state energy model

Curie’s law of paramagnetism using the Schottky energy model (Dill p. 186)

HOMEWORK PROBLEMS:

Last time

Temperature and the occupation of states

  • As we described last time, the possible energies accessible by atoms and molecules in a material are quantized (we indexed them by the ‘quantum integer’ n), and thus the state of a material can be characterized by how the energy levels are occupied in the system. What we want to do now is determine how the microstates of a material are occupied when we consider the common experimental case of a system with fixed temperature instead of fixed energy.

The Boltzmann factor and Partition function

  • Last time we discussed the constant energy ensemble (fixed E, V, and N) for a Debye (harmonic oscillator) solid. Let’s now turn to the more interesting question of how to make calculations for systems where we fix T, V, and N.
  • The ensemble for fixed (T, V, and N) includes all possible microstates for the solid that have the same temperature; it is called the canonical ensemble (‘canonical’ because it is used so often to model real systems). With the temperature fixed, the internal energy and pressure may fluctuate in the real thermodynamic system- remember that the physical situation is that the system is surrounded by a large ‘heat bath’ which passes heat back and forth with the system to maintain a constant temperature. Thus, microstates with different total internal energies may be included in the allowed microstates.
  • We used the first postulate of statistical mechanics (equal a priori probabilities) to determine the state probabilities pj in the microcanonical ensemble. For the canonical ensemble, we will determine what values the pj’s take in order to satisfy the second law for a system at constant (T,V,N). Recall that for systems with constant (T,V,N), the second law is satisfied when the Helmholtz free energy (F = U - TS) is a minimum.
  • To determine the equilibrium probability pj for each individual state j, we simply calculate what values of pj(for each possible state j) minimize the Helmholtz free energy subject to the constraint that the pj’s act like a probability and sum to 1:

(Eqn 1)

  • At equilibrium at constant temperature and volume, F is minimized in order to satisfy the second law. We want to calculate the minimum in F with respect to pj for all possible states j:

(Eqn 2)for all states j

  • To satisfy the constraint that the pj sum to 1, we use the method of Lagrange multipliers. This simply requires that we minimize with the constraint included:

(Eqn 3)for all states j

  • We can fill in U and S with our formulas from last time:

(Eqn 4)

(Eqn 5)

  • Which gives us:

(Eqn 6)for each state j

for each state j

(Eqn 7)for each state j

  • We can determine the Lagrange multiplier a by referring back to our constraint equation:

(Eqn 8)

(Eqn 9)

  • We define the summation in the denominator as the partition function Q. The importance of this sum will soon be apparent (for now, at least it simplifies our notation!).

(Eqn 10)

  • Substituting (Eqn 9) into (Eqn 7), we arrive at the solution for the probability for the system to be in some arbitrary state j:

(Eqn 11)

  • The quantity exp(-Ei/kT) is known as the Boltzmann factor: it is the ‘thermal weighting factor’ that determines how many atoms access a given state of energy i. The Boltzmann factor indicates that states with high energies will not be accessible at low temperature, but may occur with a high frequency in the ensemble at high temperature:

(Dill)

All thermodynamic quantities can be calculated from the partition function

  • The Boltzmann factor and partition function are the two most important quantities for making statistical mechanical calculations. If we have a model for a material for which we can calculate the partition function, we know everything there is to know about the thermodynamics of that model of the material.
  • All thermodynamic quantities of interest can be derived using the partition function. Some examples are:

(Eqn 12)

(Eqn 13)

(Eqn 14)

  • …hence the utility of knowing Q for a given system.

The partition function of molecules vs. multi-molecular systems

  • Typically, it is straightforward to develop models at the molecular level for allowed energies/states, and to even write the partition function for individual molecules. But how do we handle the case when we have a mole of atoms in a system and we want to determine Q? It is not possible to enumerate all the possible states by hand.
  • One way to deal with this complexity is to assume the molecules of the solid behave independently. If the each molecule of the N molecules in a solid is independent, the energy of one molecule is not influenced by that of its neighbors, and the total energy is simply the sum of the energies of each individual molecule: FIX THIS DILL. P. 182

(Eqn 15)

  • If we then look at the total partition function for the multi-molecular system, we have:

(Eqn 16)

(Eqn 17)

  • This equation gives us a simple route to making calculations for macroscopic systems from molecular level-detailed models.

The Einstein Solid

  • Now that we have the formula for the probabilities in a system at constant temperature, we can start making some powerful predictions using our Einstein solid model.

The complete partition function for an Einstein solid

  • Recall that in the Einstein solid, the atoms are assumed to vibrate in a harmonic potential. The energy of this confined oscillation is quantized:

(Eqn 18)

  • We performed microcanonical calculations (fixed E,V,N) for a very imaginary 3-atom, 1D-oscillating solid. If we take the more realistic case of allowing each atom to oscillate in X, Y, and Z space, we have 3 quantized energies:

(Eqn 19)

(Eqn 20)

(Eqn 21)

  • The total energy in one microstate (characterized by one set of values nx, ny, nz) for one atom of the solid is:

(Eqn 22)

  • It follows that the molecular partition function for one atom of the solid is:

(Eqn 23)

  • Since the virbration in each of the 3 directions is equivalent (i.e. ), the three sums in Qatom are the same. We can therefore write:

(Eqn 24)

  • We showed above that the partition function for a system of N non-interacting distinguishable atoms or molecules is given by Q = (Qatom)N. Thus, for the partition function of the entire Einstein solid, we have:

(Eqn 25)with En = (n + 1/2)h, n = 0, 1, 2, …, 

  • The infinite sum looks messy, but we can simplify this partition function. To shorten the following notation, let’s use  = 1/kT:

(Eqn 26)

  • Using the approximation (for x2 < 1), we can obtain:

(Eqn 27)

  • Finally, substituting into (Eqn 25), we arrive at the simplified total partition function for the Einstein solid:

(Eqn 28)

References

Lecture 22 – Connectiing molecular events tu