# The Boltzmann Distribution Applied to a Harmonic Oscillator

In Chapter 3 of *Intermediate Physics for Medicine and Biology*, Russ Hobbie and I discuss the Boltzmann distribution. If you have a system with energy levels of energy populated by particles in thermal equilibrium, the Boltzmann distribution gives the probability of finding a particle in the *n*th level.

A classic example of the Boltzmann distribution is for the energy levels of a harmonic oscillator. These levels are equally spaced starting from a ground state and increasing without bound. To see the power of the Boltzmann distribution, solve this new homework problem.

Section 3.7

Problem 29½.Suppose the energy levels,En, of a system are given by

En=n ε, forn= 0, 1, 2, 3, …where

εis the energy difference between adjacent levels. Assume the probability of a particle occupying thenth energy level,Pn, obeys the Boltzmann distribution

whereAis a constant,Tis theabsolute temperature, andkis theBoltzmann constant.

(a) DetermineAin terms ofε,k, andT. (Hint:the sum of the probabilities over all levels is one.)

(b) Find theaverageenergyEof the particles. (Hint:E= ∑PnEn.)

(c) Calculate theheat capacityCof a system ofNsuch particles. (Hint:U=NEandC=dU/dT.)

(d) What is the limiting value ofCfor high temperatures (kT>>ε)? (Hint: use theTaylor seriesof theexponential.)

(e) What is the limiting value ofCfor low temperatures (kT<<ε)?

(f) Sketch a plot ofCversusT.

You may need theseinfinite series

1 +x+x^2 +x^3 + ⋯ = 1/(1−x)x

+ 2x^2 + 3x^3 + ⋯ =x/(1−x)^2

This is a somewhat advanced problem in statistical mechanics, so I gave several hints to guide the reader. The calculation contains much interesting physics. For instance, the answer to part (e) is known as the third law of thermodynamics. Albert Einstein was the first to calculate the heat capacity of a collection of harmonic oscillators (a good model for a crystalline solid). There s more physics than biology in this problem, because most of the interesting behavior occurs at cold temperatures but biology operates at hot temperatures.

If you’re having difficulty solving this problem, here’s one more hint:

e*^nx* = (e^*x*)*^n* .

Enjoy!

*Originally published at **http://hobbieroth.blogspot.com**.*