6.1: Separation of Contributions
- Page ID
- 285776
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Collective Degrees of Freedom
In Section [subsection:Einstein_Debye] we have seen that the treatment of condensed phases can be complicated by collective motion of particles. Such effects are absent in an ideal gas that consists of point particles, a model that is reasonable for noble gases far from condensation. For gases consisting of molecules, it does not suffice to consider only translational motion as in Maxwell’s kinetic theory of gases. We see this already when considering H\(_2\) gas, where each molecule can be approximated by a harmonic oscillator (Section [section:harmonic_oscillator]). Neglect of the vibrational degrees of freedom will lead to wrong results for internal energy, heat capacity, the partition function, and entropy, at least at high temperatures. In fact, an H\(_2\) molecule is not only an oscillator, it is also a rotor. As a linear molecule it has two rotational degrees of freedom, which also contribute to internal energy and to the partition function.
In principle, we could try to ignore all this and treat each atom as one particle. If the Hamiltonian includes the potentials that characterize interaction between the particles, our equations of motion would be correct. In practice, such a treatment is inconvenient and it is better to group the spatial degrees of freedom according to the type of motion. The H\(_2\) molecule has 3 translational degrees of freedom, 2 rotational degrees of freedom, and 1 vibrational degree of freedom in the collective motion picture. The sum is 6, as expected for two atoms with each of them having 3 translational degrees of freedom in an ’independent’ motion picture. In general, a molecule with \(n\) atoms has \(f_s=3n\) spatial degrees of freedom, 3 of which are translational, 3 are rotational, except for linear molecules, which have only 2 rotation degrees of freedom, and the rest are vibrational. We note that the number of degrees of freedom in phase space is \(f = 2 f_s\) because each spatial degree of freedom is also assigned a momentum degree of freedom.
These considerations take care of particle motion. Further contributions to internal energy and to the partition function can arise from spin. In both closed-shell and open-shell molecules, nuclear spin can play a role. This is indeed the case for H\(_2\), which can exist in ortho and para states that differ in correlation of the nuclear spins of the two hydrogen atoms. For open-shell molecules electron spin degrees of freedom must be considered. This is the case, for instance, for O\(_2\), which has a triplet ground state. In this case, rotational and spin degrees of freedom correspond to similar energies and couple. Finally, at sufficiently high temperatures electronic excitation becomes possible and then also makes a contribution to the partition function.
Factorization of Energy Modes
In many cases, the individual contributions are separable, i.e. the modes corresponding to different types of motions can be treated independently. Roughly speaking, this results from a separation of energy ranges (frequency bands) of the modes and a corresponding separation of time scales. Nuclear spin degrees of freedom have much lower energy than rotational degrees of freedom which usually have much lower energy than vibrational degrees of freedom which have much lower energies than electronic excitation. The independence of nuclear and electron motion is basis of the Born-Oppenheimer approximation and the independence of rotational and vibrational motion is invoked when treating a molecule as a rigid rotor. Separability of energy modes leads to a sum rule for the energy contributions for a single closed-shell molecule ,
\[\epsilon_j = \epsilon_{j,\mathrm{trs}} + \epsilon_{j,\mathrm{ns}} + \epsilon_{j,\mathrm{rot}} + \epsilon_{j,\mathrm{vib}} + \epsilon_{j,\mathrm{el}} \ , \label{eq:epsilon_mol}\]
where \(\epsilon_{j,\mathrm{trs}}\), \(\epsilon_{j,\mathrm{ns}}\), \(\epsilon_{j,\mathrm{rot}}\), \(\epsilon_{j,\mathrm{vib}}\), and \(\epsilon_{j,\mathrm{el}}\) are the translational, nuclear spin, rotational, vibrational, and electronic contributions, respectively. For a monoatomic molecule (atom) \(\epsilon_{j,\mathrm{rot}}\) and \(\epsilon_{j,\mathrm{vib}}\) vanish. If both the number of neutrons and of protons in the nucleus is even, the nucleus has spin \(I=0\). In that case the nuclear spin contribution vanishes for an atom, even in the presence of an external magnetic field. If all nuclei have spin zero, the nuclear spin contribution also vanishes for a diatomic or multi-atomic molecule.
If we assume the equipartition theorem to hold, or even more generally, the whole system to attain thermal equilibrium, there must be some coupling between the different modes. If we say that the energy modes are separable, we assume weak coupling, which means that for statistical purposes we can assume the modes to be independent of each other. The consequence for the computation of the partition function can be seen by considering a system of \(N\) particles with an \(\alpha\) mode associated with quantum number \(k\) and an \(\omega\) mode associated with quantum number \(r\). The total energy of a single molecule of this type is \(\epsilon_j = \epsilon_{j,\alpha k} + \epsilon_{j,\omega r}\) The molecular partition function is given by
\[Z = \sum_k \sum_r e^{-\beta \left( \epsilon_{\alpha k} + \epsilon_{\omega r}\right)} \ .\]
This sum can be rewritten as
\[\begin{align} Z & = \sum_k \sum_r e^{-\beta \epsilon_{\alpha k}} \cdot e^{-\beta \epsilon_{\omega r}} \\ & = \sum_k e^{-\beta \epsilon_{\alpha k}} \sum_r e^{-\beta \epsilon_{\omega r}} \\ & = \sum_k e^{-\beta \epsilon_{\alpha k}} Z_\omega = Z_\omega \sum_k e^{-\beta \epsilon_{\alpha k}} \\ & = Z_\alpha Z_\omega \ .\end{align}\]
We see that the total partition function is the product of the partition functions corresponding to the individual modes. This consideration can be extended to multiple modes. With Equation \ref{eq:epsilon_mol} it follows that
\[Z = Z_\mathrm{trs} \cdot Z_\mathrm{ns} \cdot Z_\mathrm{rot} \cdot Z_\mathrm{vib} \cdot Z_\mathrm{el} \ . \label{eq:z_factorization}\]
By considering Equation \ref{eq:molar_z_from_molecular_z} or Equation \ref{eq:z_indist} we see that we can also compute the partition function for a given mode for all \(N\) particles before multiplying the modes. We have already seen that we must set \(z_\mathrm{trs} = Z_\mathrm{trs}^N/N!\) to heal the Gibbs paradox. What about the other, internal degrees of freedom? If two particles with different internal states are exchanged, they must be considered distinguishable, exactly because their internal state ’tags’ them. Hence, for all the other modes we have \(z_\alpha = Z_\alpha^N\). Thus,
\[\begin{align} z & = \frac{1}{N!} \left( Z_\mathrm{trs} \cdot Z_\mathrm{ns} \cdot Z_\mathrm{rot} \cdot Z_\mathrm{vib} \cdot Z_\mathrm{el} \right)^N \\ & = \frac{Z_\mathrm{trs}^N}{N!} \cdot Z_\mathrm{ns}^N \cdot Z_\mathrm{rot}^N \cdot Z_\mathrm{vib}^N \cdot Z_\mathrm{el}^N \ . \label{eq:z_separation}\end{align}\]
Accordingly, we can consider each of the partition functions in turn. We also note that separability of the energies implies factorization of the molecular wavefunction,
\[\psi = \psi_\mathrm{trs} \cdot \psi_\mathrm{ns} \cdot \psi_\mathrm{rot} \cdot \psi_\mathrm{vib} \cdot \psi_\mathrm{el}\]