Rotational partition function explained

In chemistry, the rotational partition function relates the rotational degrees of freedom to the rotational part of the energy.

Definition

Z of a system of N identical, indistinguishable, noninteracting atoms or molecules can be divided into the atomic or molecular partition functions \zeta :^[1] $$Z\; =\; \backslash frac$$with:$$\backslash zeta\; =\; \backslash sum\_j\; g\_j\; e^,$$where g _{j} is the degeneracy of the j-th quantum level of an individual particle, k _{B} is the Boltzmann constant, and T is the absolute temperature of system.For molecules, under the assumption that total energy levels E _{j} can be partitioned into its contributions from different degrees of freedom (weakly coupled degrees of freedom)^[2] $$E\_j\; =\; \backslash sum\_i\; E\_j^i\; =\; E\_j^\backslash text\; +\; E\_j^\backslash text\; +\; E\_j^\backslash text\; +\; E\_j^\backslash text\; +\; E\_j^\backslash text$$and the number of degenerate states are given as products of the single contributions$$g\_j\; =\; \backslash prod\_i\; g\_j^i\; =\; g\_j^\backslash text\; g\_j^\backslash text\; g\_j^\backslash text\; g\_j^\backslash text\; g\_j^\backslash text,$$where "trans", "ns", "rot", "vib" and "e" denotes translational, nuclear spin, rotational and vibrational contributions as well as electron excitation, the molecular partition functions$$\backslash zeta\; =\; \backslash sum\_j\; g\_j\; e^$$can be written as a product itself$$\backslash zeta\; =\; \backslash prod\_i\; \backslash zeta^i\; =\; \backslash zeta^\backslash text\; \backslash zeta^\backslash text\; \backslash zeta^\backslash text\; \backslash zeta^\backslash text\backslash zeta^\backslash text.$$

Linear molecules

Rotational energies are quantized. For a diatomic molecule like CO or HCl, or a linear polyatomic molecule like OCS in its ground vibrational state, the allowed rotational energies in the rigid rotor approximation are$$E\_J^\backslash text\; =\; \backslash frac\; =\; \backslash frac\; =\; J(J+1)B.$$J is the quantum number for total rotational angular momentum and takes all integer values starting at zero, i.e.,

J = 0 , 1 , 2 , \ldots , B = (\hbar ^{2})/(2 I) is the rotational constant, and I is the moment of inertia. Here we are using B in energy units. If it is expressed in frequency units, replace B by hB in all the expression that follow, where h is Planck's constant. If B is given in units of cm ^{-1} , then replace B by hcB where c is the speed of light in vacuum.

For each value of J, we have rotational degeneracy,

g _{j} = (2J+1), so the rotational partition function is therefore$$\backslash zeta^\backslash text\; =\; \backslash sum\_^\backslash infty\; g\_j\; e^\; =\; \backslash sum\_^\backslash infty\; (2J+1)\; e^.$$

For all but the lightest molecules or the very lowest temperatures we have

B \ll k _{B} T . This suggests we can approximate the sum by replacing the sum over J by an integral of J treated as a continuous variable.$$\backslash zeta^\backslash text\; \backslash approx\; \backslash int\_0^\; (2J+1)e^\; dJ\; =\; \backslash frac\; .$$

This approximation is known as the high temperature limit. It is also called the classical approximation as this is the result for the canonical partition function for a classical rigid rod.

Using the Euler–Maclaurin formula an improved estimate can be found^[3] $$\backslash zeta^\backslash text\; =\; \backslash frac\; +\; \backslash frac\; +\; \backslash frac\; \backslash left(\backslash frac\; \backslash right)\; +\; \backslash frac\; \backslash left(\backslash frac\; \backslash right)^2\; +\; \backslash frac\; \backslash left(\backslash frac\; \backslash right)^3\; +\; \backslash cdots\; .$$

For the CO molecule at

T = 300K , the (unit less) contribution \zeta ^{rot} to \zeta turns out to be in the range of 10 ^{2} .

The mean thermal rotational energy per molecule can now be computed by taking the derivative of

\zeta ^{rot} with respect to temperature T . In the high temperature limit approximation, the mean thermal rotational energy of a linear rigid rotor is k _{B} T .

Quantum symmetry effects

For a diatomic molecule with a center of symmetry, such as

\rm H _{2,} N _{2,} CO _{2,} or H _{2} C _{2} (i.e. D _{infty} point group), rotation of a molecule by \pi radian about an axis perpendicular to the molecule axis and going through the center of mass will interchange pairs of equivalent atoms. The spin–statistics theorem of quantum mechanics requires that the total molecular wavefunction be either symmetric or antisymmetric with respect to this rotation depending upon whether an even or odd number of pairs of fermion nuclear pairs are exchanged. A given electronic & vibrational wavefunction will either be symmetric or antisymmetric with respect to this rotation. The rotational wavefunction with quantum number J will have a sign change of ( -1 ) ^{J} . The nuclear spins states can be separated into those that are symmetric or antisymmetric with respect to the nuclear permutations produced by the rotation. For the case of a symmetric diatomic with nuclear spin quantum number I for each nucleus, there are (I +1 )( 2 I +1 ) symmetric spin functions and I ( 2 I +1 ) are antisymmetric functions for a total number of nuclear functions g ^{ns} = ( 2 I + 1 ) ^{2} . Nuclei with an even nuclear mass number are bosons and have integer nuclear spin quantum number, I. Nuclei with odd mass number are fermions and had half integer I. For the case of H₂, rotation exchanges a single pair of fermions and so the overall wavefunction must be antisymmetric under the half rotation. The vibration-electronic function is symmetric and so the rotation-vibration-electronic will be even or odd depending upon whether J is an even or odd integer. Since the total wavefunction must be odd, the even J levels can only use the antisymmetric functions (only one for I = 1/2) while the odd J levels can use the symmetric functions (three for I = 1/2). For D2, I = 1 and thus there are six symmetric functions, which go with the even J levels to produce an overall symmetric wavefunction, and three antisymmetric functions that must go with odd J rotational levels to produce an overall even function. The number of nuclear spin functions that are compatible with a given rotation-vibration-electronic state is called the nuclear spin statistical weight of the level, often represented as g _{J} . Averaging over both even and odd J levels, the mean statistical weight is ( 1 / 2 ) ( 2 I +1 ) ^{2} , which is one half the value of g ^{ns} expected ignoring the quantum statistical restrictions. In the high temperature limit, it is traditional to correct for the missing nuclear spin states by dividing the rotational partition function by a factor \sigma = 2 with \sigma known as the rotational symmetry number which is 2 for linear molecules with a center of symmetry and 1 for linear molecules without.

Nonlinear molecules

A rigid, nonlinear molecule has rotational energy levels determined by three rotational constants, conventionally written

A, B, and C , which can often be determined by rotational spectroscopy. In terms of these constants, the rotational partition function can be written in the high temperature limit as ^[4] $$\backslash zeta^\backslash text\; \backslash approx\; \backslash frac\; \backslash sqrt$$with \sigma again known as the rotational symmetry number ^[5] which in general equals the number ways a molecule can be rotated to overlap itself in an indistinguishable way, i.e. that at most interchanges identical atoms. Like in the case of the diatomic treated explicitly above, this factor corrects for the fact that only a fraction of the nuclear spin functions can be used for any given molecular level to construct wavefunctions that overall obey the required exchange symmetries. Another convenient expression for the rotational partition function for symmetric and asymmetric tops is provided by Gordy and Cook:$$\backslash zeta^\backslash text\; \backslash approx\; \backslash frac\; \backslash sqrt$$where the prefactor comes from $$\backslash sqrt\; =\; 5.34\; \backslash times\; 10^6$$when A, B, and C are expressed in units of MHz. ^[6]

The expressions for

\zeta ^{rot} works for asymmetric, symmetric and spherical top rotors.

References

1. Donald A. McQuarrie, Statistical Mechanics, Harper & Row, 1973
2. Donald A. McQuarrie, ibid
3. G. Herzberg, Infrared and Raman Spectra, Van Nostrand Reinhold, 1945, Equation (V,21)
4. G. Herzberg, ibid, Equation (V,29)
5. G. Herzberg, ibid; see Table 140 for values for common molecular point groups
6. Book: Cook . Robert L. . Gordy . Walter . Microwave Molecular Spectra . 1970 . Interscience Pub. . 0-471-08681-9 . 56–57 . 2.

See also

• Translational partition function
• Vibrational partition function
• Partition function (mathematics)