07.04.2013 Views

Essentials of Computational Chemistry

Essentials of Computational Chemistry

Essentials of Computational Chemistry

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.

10.3 ENSEMBLE PROPERTIES AND BASIC STATISTICAL MECHANICS 359<br />

from the second to the third line, the sum has been changed so that it goes over discrete<br />

energy levels, rather than individual states, and gk is the degeneracy <strong>of</strong> level k. The quantity<br />

in brackets on the third line defines the molecular partition function q.<br />

A second consequence <strong>of</strong> the ideal gas assumption is that PV in Eq. (10.4) may be<br />

replaced by NkBT . In the special case where we are working with one mole <strong>of</strong> molecules,<br />

in which case N = NA (Avogadro’s number), we may replace PV with RT ,whereR is the<br />

universal gas constant (8.3145 J mol −1 K −1 ).<br />

10.3.2 Separability <strong>of</strong> Energy Components<br />

We have thus reduced the problem from finding the ensemble partition function Q to finding<br />

the molecular partition function q. In order to make further progress, we assume that the<br />

molecular energy ε can be expressed as a separable sum <strong>of</strong> electronic, translational, rotational,<br />

and vibrational terms, i.e.,<br />

q(V,T) =<br />

=<br />

=<br />

levels <br />

k<br />

gke −εk(V )/kBT<br />

levels <br />

gke<br />

k<br />

−[εelec+εtrans(V )+εrot+εvib]k/kBT<br />

<br />

elec<br />

gie<br />

i<br />

−εi/kBT<br />

⎡<br />

trans <br />

⎣ gje<br />

j<br />

−εj (V )/kBT<br />

⎤<br />

⎦<br />

rot<br />

<br />

gke −εk/kBT<br />

k<br />

vib<br />

<br />

gle −εl/kBT<br />

<br />

= qelec(T )qtrans(V, T )qrot(T )qvib(T ) (10.8)<br />

where again advantage is taken <strong>of</strong> the ability to express an exponential <strong>of</strong> sums as a product<br />

<strong>of</strong> sums <strong>of</strong> exponentials, and the separate lines make clear that the degeneracy <strong>of</strong> a total<br />

molecular energy level is simply the product <strong>of</strong> the degeneracies <strong>of</strong> each <strong>of</strong> its contributing<br />

components.<br />

Note in Eqs. (10.3) and (10.5), Q always appears as the argument <strong>of</strong> the natural logarithm<br />

function. Using Eqs. (10.7) and (10.8), our assumptions to this point allow us to write<br />

<br />

[qelec(T )qtrans(V, T )qrot(T )qvib(T )]<br />

ln[Q(N,V,T)] = ln<br />

N <br />

N!<br />

= N{ln[qelec(T )] + ln[qtrans(V, T )] + ln[qelec(T )]<br />

+ ln[qvib(T )]}−ln(N!)<br />

≈ N{ln[qelec(T )] + ln[qtrans(V, T )]<br />

+ ln[qelec(T )] + ln[qvib(T )]}−N ln N + N (10.9)<br />

where going from the second to the third equality makes use <strong>of</strong> Stirling’s approximation for<br />

ln(N!) when N is large. This separation <strong>of</strong> terms by the logarithm function makes evident<br />

l

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

Saved successfully!

Ooh no, something went wrong!