1
$\begingroup$

On p. 162 of Statistical Mechanics: R K Pathria, Paul D. Beale, the partition function for a rotating diatomic molecule is derived as

$$ Z_{rot}(T) = \frac{T}{\Theta_r} + \frac{1}{3} + \frac{1}{15}\frac{\Theta_r}{T} + \frac{4}{315}\bigg(\frac{\Theta_r}{T}\bigg)^2 + ...$$

In the text they say that the corresponding specific heat is

$$ C_V = Nk\bigg\{1 + \frac{1}{45}\bigg(\frac{\Theta_r}{T}\bigg)^2 + \frac{16}{945}\bigg(\frac{\Theta_r}{T}\bigg)^3 + ...\bigg\}$$

I want to know where this formula for $C_V$ came from? I tried using the relations $C_V = \big(\frac{\partial U}{\partial T}\big)_{N,V}$ and $U = \frac{kT^2}{N}\frac{\partial \ln Z}{\partial T}$ but just ended up with a long expression that didn't simplify.

$\endgroup$

1 Answer 1

2
$\begingroup$

Firstly, I think $N$ might need to be in the numerator of the formula for $U$ in terms of $Z$.

Secondly, let $x = T/\Theta_r$, then \begin{align} \frac{\partial}{\partial T}\left(T^2\frac{\partial\ln Z}{\partial T}\right) & = \frac{\partial}{\partial (T/\Theta_r)}\left((T/\Theta_r)^2\frac{\partial\ln Z}{\partial (T/\Theta_r)} \right) = \frac{\partial}{\partial x}\left(x^2\frac{\partial\ln Z}{\partial x} \right) = \frac{\partial}{\partial x}\left(\frac{x^2}{Z}\frac{\partial Z}{\partial x}\right)\\ &= \frac{\partial}{\partial x}\left[\frac{x\left(1 - \frac{1}{15}\frac{1}{x^2} - \frac{8}{315}\frac{1}{x^3}+\cdots\right)}{1+\frac{1}{3}\frac{1}{x}+\frac{1}{15}\frac{1}{x^2}+\frac{4}{315}\frac{1}{x^3}+\cdots}\right] \\ &= \frac{\partial}{\partial x}\left[x\left(1 - \frac{1}{15}\frac{1}{x^2} - \frac{8}{315}\frac{1}{x^3}+\cdots\right)\left(1-\frac{1}{3}\frac{1}{x}+\frac{2}{45}\frac{1}{x^2}-\frac{1}{189}\frac{1}{x^3}+\cdots\right)\right]. \end{align} In going from the second-to-last line to the last line, one needs to invert the power series to third order in the small quantity $1/x$ (see below for how to do this). Now simply multiply out the expression in brackets keeping only the terms up to order $1/x^2$, and take the remaining partial derivative to obtain the desired result.

Inversion of Power Series.

Consider the power series $1 + a_1 \varepsilon + a_2 \varepsilon^2 + a_3 \varepsilon^3 + \cdots$. We are looking for the power series $1 + b_1 \epsilon + b_2 \varepsilon^2 + b_3\varepsilon^3 + \cdots$ that is its reciprocal. This means that their product should be 1: $$ (1 + a_1 \varepsilon + a_2 \varepsilon^2 + a_3 \varepsilon^3 + \cdots)(1 + b_1 \varepsilon + b_2 \varepsilon^2 + b_3\varepsilon^3 + \cdots) = 1 $$ Now simply multiply out the terms to the desired order, set the coefficient of each order in $\epsilon$ to zero, and solve the resulting system of equations order-by-order to obtain the desired coefficients $b_1, b_2, b_3\, \dots$.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.