27.04.2013 Views

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)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!