Expansion coefficient and the Lennard-Jones potential

Consider a classical particle moving in a one-dimensional potential well u(x). The particle is in thermal equilibrium with a reservoir at temperature T, so the probabilities of its various states are determined by Boltzmann statistics.

(a) Show that the average position of the particle is given by

x =∫xe-βu(x) dx/∫e-βu(x) dx,

where each integral is over the entire x axis.

(b) If the temperature is reasonably low (but still high enough for classical mechanics to apply), the particle will spend most of its time near the bottom of the potential well. In that case we can expand u(x) in a Taylor series about the equilibrium point x0:

u(x) = u(x_0) + (x - x_0) du/dx |x_0 + ½ (x + x_0)2 d^2u/dx^2 |x_0 + 1/3! (x - x_0)^3 d^3u/dx^3 |x_0 + …

Show that the linear term must be zero, and that truncating the series after the quadratic term results in the trivial prediction x = x_0.

(c) If we keep the cubic term in the Taylor series as well, the integrals in the formula for x become difficult. To simplify them, assume that the cubic term is small, so its exponential can be expanded in a Taylor series (leaving the quadratic term in the exponent). Keeping only the smallest temperature-dependent term, show that in this limit x differs from x0 by a term proportional to kT. Express the coefficient of this term in terms of the coefficients of the Taylor series for u(x).

(d) The interaction of noble gas atoms can be modeled using the Lennard-Jones potential,

u(x) = u0[(x_0/x)^12 - 2(x_0/x)^6].

Sketch this function, and show that the minimum of the potential well is x = x0, with depth u0. For argon, x_0 = 3.9 Å and u_0 = 0.010 eV. Expand the Lennard-Jones potential in a Taylor series about the equilibrium point, and use the result of part (c) to predict the linear thermal expansion coefficient of a noble gas crystal in terms of u_0. Evaluate the result numerically for argon, and compare to the measured value α = 0.0007 K-1 (at 80 K).