04.08.2013 Views

Statistical Mechanics - Physics at Oregon State University

Statistical Mechanics - Physics at Oregon State University

Statistical Mechanics - Physics at Oregon State University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

150 CHAPTER 7. CLASSICAL STATISTICAL MECHANICS.<br />

Similar to wh<strong>at</strong> we derived before, it is possible to show th<strong>at</strong> the last, exponential,<br />

factor in the logarithm does not give a contribution to the integral. We<br />

expand the logarithm to get:<br />

1<br />

S(U, V, N) = −kB<br />

N!h3N <br />

d XC 1<br />

ɛ √ (U−H)2<br />

e− ɛ<br />

π 2<br />

<br />

log(C) − log(ɛ √ (U − H)2<br />

π) −<br />

ɛ2 <br />

(7.34)<br />

The third term in the integral is very small when |U − H| ≫ ɛ due to the<br />

exponential. It is also very small when |U − H| ≪ ɛ because of the quadr<strong>at</strong>ic<br />

term. Hence it only contributes when |U − H| ≈ ɛ. In th<strong>at</strong> case the integrant<br />

is proportional to 1<br />

ɛ . But the volume in phase space th<strong>at</strong> corresponds to this<br />

value is proportional to ɛ6N−1 and hence the total contribution vanishes when<br />

ɛ goes to zero.<br />

The second factor in the logarithm gives a contribution kB log(ɛ √ π) to the<br />

entropy, and this factor vanishes in the thermodynamic limit. Hence we find<br />

with<br />

S(U, V, N) = −kB log(C) = kB log(D(U, V, N)) (7.35)<br />

D(U, V, N) = C −1 =<br />

1<br />

N!h3N ɛ √ <br />

π<br />

d (U−H)2<br />

−<br />

Xe ɛ2 (7.36)<br />

Here we can again take the limit ɛ → 0 to recover the expression we derived<br />

before. But note th<strong>at</strong> we have used the thermodynamic limit, and we did take<br />

it before the limit ɛ → 0. Th<strong>at</strong> is technically not correct. This was a topic of<br />

much discussion in the early days of the development of st<strong>at</strong>istical mechanics.<br />

For an ideal gas the density of st<strong>at</strong>es D follows from<br />

D(U, V, N) =<br />

1<br />

N!h3N <br />

dp 3 1 · · · dp 3 N δ(U −<br />

N<br />

i=1<br />

p 2 i<br />

2m )<br />

<br />

dr 3 1 · · · dr 3 N<br />

(7.37)<br />

The sp<strong>at</strong>ial integral is easy and gives a factor V N . The integral over the momenta<br />

is equal to the surface area of a 3N-dimensional hyper-sphere of radius<br />

√ 2mU. This leads to<br />

If N is even, Γ( 3N<br />

2<br />

D(U, V, N) =<br />

) = ( 3N<br />

2<br />

N V<br />

N!h3N 2π 3N<br />

√2mU3N−1 2<br />

Γ( 3N<br />

2 )<br />

(7.38)<br />

3N<br />

−1)! and if N is odd, Γ( 2 ) = √ 3N−1<br />

− π2 2 (3N −2)!!<br />

In both cases we can approxim<strong>at</strong>e log(Γ( 3N 3N−1<br />

2 )) by ( 2 ) log( 3N 3N<br />

2 ) − 2 .<br />

This assumes th<strong>at</strong> N is very large (thermodynamic limit) and th<strong>at</strong> we can<br />

replace N −1 by N. The entropy in this limit is (again ignoring terms not linear<br />

in N)

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

Saved successfully!

Ooh no, something went wrong!