10. Appendix
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
618 <strong>Appendix</strong> B<br />
approximation we can expand the term º(v) as a function of v about the equilibrium<br />
volume v0 and keep only terms up to the quadratic term:<br />
º(v) ∼ º(v0) (dº/dv)v0 (v v0) (1/2)(d 2 º/dv 2 2<br />
)V0 (v v0)<br />
(3.17.9)<br />
By definition, the linear term (dº/dv)V0 vanishes at the equilibrium volume so<br />
we can take the derivative of º(v) with respect to v and obtain:<br />
(dº/dv) ∼ (d 2 º/dv 2 )v0 (v v0) . (3.17.10)<br />
On the other hand, if we have started from (3.17.8) and take the derivative<br />
with respect to v at equilibrium volume then we would obtain:<br />
(dP/dv)V0 ∼ (d2 º/dv 2 )v0 (3.17.11)<br />
by noting that U is proportional to v. (3.17.11) can be simplified by introducing<br />
the isothermal compressibility Î0 (1/v0)(dv/dP)T. The bulk modulus B<br />
defined on Page 139 is related to Î by: B 1/Î. With this simplification we<br />
arrived at:<br />
(dP/dv)T 1/Î0v0 . (3.17.12)<br />
Combining (3.17.10) – (3.17.12) we obtain:<br />
(dº/dv) (dP/dv)V0 (v v0) (v v0)/Î0v0<br />
(3.17.13)<br />
On substituting (3.17.13) back into (3.17.8) we obtain the equation-of-state<br />
within the quasi-harmonic approximation:<br />
P (v v0)/Î0v0 〈Á〉U/v . (3.17.14)<br />
This equation was first derived by Mie and later by Grüneisen.<br />
To obtain the coefficient of thermal expansion ‚ (1/v)(v/T)P we take<br />
the partial derivative of (3.17.14) with respect to T while keeping P constant:<br />
0 ‚v<br />
<br />
〈Á〉 <br />
Î0v0<br />
‚U<br />
<br />
U<br />
‚ <br />
v v T<br />
1<br />
<br />
U<br />
. (3.17.15)<br />
v T v<br />
The quantities inside the braces can be expressed in terms of measurable<br />
quantities, such as the heat capacity. For example, the constant volume heat<br />
capacity Cv is equal to (U/T)V. By differentiating the equation (3.17.7) for<br />
U with respect to v, we can show that:<br />
(U/v)T (〈Á〉CvT/v) (〈Á〉U/v) . (3.17.16)<br />
Actually this calculation is very similar to the one used to derive (3.17.5) by<br />
taking the derivative of Fvib in (3.17.4) with respect to v. Substituting back<br />
(3.17.16) into (3.17.15) we obtain:<br />
‚V<br />
<br />
Î0v0<br />
〈Á〉‚U<br />
<br />
CVT U<br />
‚〈Á〉2 <br />
v v v<br />
〈Á〉Cv<br />
(3.17.17)<br />
v<br />
Since we have assumed a weak anharmonic term we expect that terms<br />
quadratic in 〈Á〉 can be neglected. This means that at P 0 (when v v0)