15 October 1996

Your browser must support super and subscripts to view this page!

Stathermo Calculations, Page 3


Previous Page

The integral approximation which was spectacularly successful for translation is only appropriate here for non-hydrogen rotations. With an atom as light as H, the moment of inertia, I, becomes quite small, making B large. Indeed, since BJ(J+1)/kT is the dimensionless exponent in the partition function, that means that B/k has units of temperature. In fact, it is called  Theta rot and measures the temperature at which significant rotational excitation can begin to occur. For the vast majority of molecules, such as carbon dioxide, B/k < 1 K, but this rotational onset temperature for hydrogen molecule is 88 K. That means that only a handful of rotational states are available at room temperature for hydrogen, making an integral approximation for z perilous.

But we'll make it anyway! In cases where B is too large, that only means we'll actually have to calculate the sum over states, but few states will contribute. In normal circumstances, the following will suffice instead. Same trick; different mode.

zrot =  Summation (2J+1)exp[-BJ(J+1)/kT]

~  Integral(0 to infinity) exp[-BJ(J+1)/kT] (2J+1)dJ

=  Integral(0 to infinity) exp[-BJ(J+1)/kT] d[J(J+1)]

=  Integral(0 to infinity) exp[-Bx/kT]dx where x=J(J+1)

= kT/B = T/ Theta rot

If we believe it for diatomic hydrogen, it says that since room temperature is only about 3 Theta rot, z=3 there. Certainly not enough states for a continuum over which to integrate! But for carbon dioxide under the same circumstances, it says z~600; good choice for integration.

Unfortunately, in both these cases, T/ Theta rot overestimates z by a factor of 2. This is because we've assumed that a 360° rotation gives us an indistinguishable configuration of the molecule. But both of these molecules become indistinguishable under only 180° of rotation! We've overcounted.

zrot = kT/ sigma B is the way this correction factor is normally inserted.  sigma is only appropriate for molecules whose symmetry is high enough to include rotational symmetries. So normal water has  sigma =2 but HOD doesn't. Standard molecular Symmetry Tables give  sigma directly as the sum of all real rotations (designated as C operations) plus the identity operation (designated E). Remember to include the number given prior to each C as this represents the number of different rotations of that type present in the molecule!

Not surprisingly, for example, ammonia has  sigma =3.

Since zrot has a single factor of T, as opposed to translation's T3/2, the derivation of Erot will yield RT instead of (3/2)RT. However, for non-linear molecule, like ammonia, there will be 3 rotation axes, not 2. Indeed, the general (integral approximation) rotational partition function for non-linear molecules will be given by

zrot = (1/ sigma ) (kT)3/2( pi /ABC)½ where A, B, and C are the  hbar 2/2I) for each separate rotation axis.


Vibration


The exact partition function for the Harmonic Oscillator (a not-unreasonable model for molecular vibrations near ground state) was given in the solution of the first quiz (fall '95), and was there found to be

zvib = 1/[1-exp(- Theta vib/T)]

where  Theta vib = h nu /k which is generally several times larger than room temperature. In other words, at room temperature, vibrations haven't turned on, and almost all of the molecules will be in their ground vibrational states for all vibration modes.

How many vibration modes can there be? Remember 3N-6 (or 3N-5 if a linear molecule)? That was a count of all the "degrees of freedom" of motion of the molecule (3N from x,y,z for each of the N atoms) less the Center of Mass translation (X,Y,Z of the molecule as a whole) less the 3 (2 if linear) rotational degrees of freedom. An analysis of how each symmetry operation in a Symmetry Table distorts each atomic displacement (dx, dy, dz) can tell you not only how many vibrations there are but also the degeneracy of each and which ones will be infrared or Raman active. See . . . Symmetry Tables are your friend.

If we substitute that partition function in the molar energy expression, we find

Evib = R Theta vib / [exp(+ Theta vib/T) - 1] for each vibration mode with its unique  Theta vib.

In the limit of T>> Theta vib, E reduces to RT for each vibration, but such a T is so high that the molecule will have disintegrated long since!


Return to the CHM 5414 Lecture Notes or Go To Next or Previous Lectures.

Chris Parr University of Texas at Dallas Programs in Chemistry, Room BE3.506 P.O. Box 830688 M/S BE2.6 (for snailmail) Richardson, TX 75083-0688
Voice: (972) 883-2485 Fax: (972) 883-2925 BBS: (972) 883-2168 (HST) or -2932 (V.32bis) Internet: parr@utdallas.edu (Click on that address to send Chris e-mail.)

Last modified 21 October 1996.