Compare the classical and quantum rigid rotor in three dimensions
Demonstrate how to use the Separation of Variable technique to solve the 3D rigid rotor Schrödinger Equation
Identify and interpret the two quantum numbers for a 3D quantum rigid rotor including the range of allowed values
Describe the wavefunctions of the 3D quantum rigid rotor in terms of nodes, average displacements and most probable displacements
Describe the energies of the 3D quantum rigid rotor in terms of values and degeneracies
Rigid rotor means when the distance between particles do not change as they rotate. A rigid rotor only approximates a rotating diatomic molecular if vibration is ignored.
The Classical Rigid Rotor in 3D
The rigid rotor is a mechanical model that is used to explain rotating systems. The linear rigid rotor model consists of two point masses located at fixed distances from their center of mass. The fixed distance between the two masses and the values of the masses are the only characteristics of the rigid model. However, for many actual diatomics this model is too restrictive since distances are usually not completely fixed and corrections on the rigid model can be made to compensate for small variations in the distance. Even in such a case the rigid rotor model is a useful model system to master.
For a rigid rotor, the total energy is the sum of kinetic (\(T\)) and potential (\(V\)) energies
\[E_{tot} = T + V \label{5.8.2} \]
The potential energy, \(V\), is set to \(0\) because the distance between particles does not change within the rigid rotor approximation. However, In reality, \(V \neq 0\) because even though the average distance between particles does not change, the particles still vibrate. The rigid rotor approximation greatly simplifys our discussion.
Since \(V=0\) then \(E_{tot} = T\) and we can also say that:
Equation \ref{5.8.10} shows that the energy of the rigid rotor scales with increasing angular frequency (i.e., the faster is rotates) and with increasing moment of inertia (i.e, the inertial resistance to rotation). Also, as expected, the classical rotational energy is not quantized (i.e., all possible rotational frequencies are possible).
The Quantum Rigid Rotor in 3D
It is convenient to discuss rotation with in the spherical coordinate system rather than the Cartesian system (Figure 5.8.1
).
To solve the Schrödinger equation for the rigid rotor, we will separate the variables and form single-variable equations that can be solved independently. Only two variables \(\theta\) and \(\varphi\) are required in the rigid rotor model because the bond length, \(r\), is taken to be the constant \(r_0\). We first write the rigid rotor wavefunctions as the product of a theta-function depending only on \(\theta\) and a phi-function depending only on \(\varphi\)
Since \(r = r_0\) is constant for the rigid rotor and does not appear as a variable in the functions, the partial derivatives with respect to \(r\) are zero; i.e. the functions do not change with respect to \(r\). We also can substitute the symbol \(I\) for the moment of inertia, \(\mu r^2_0\) in the denominator of the left hand side of Equation \(\ref{5.8.13}\), to give
To begin the process of the Separating of Variables technique, multiply each side of Equation \(\ref{5.8.14}\) by \(\dfrac {2I}{\hbar ^2}\) and \(\dfrac {-\sin ^2 \theta}{\Theta (\theta) \Phi (\varphi)} \) to give
Note that this \(\lambda\) has no connection to a wavelength; it is merely being used as an algebraic symbol for the combination of constants shown in Equation \(\ref{5.8.16}\).
Inserting \(\lambda\), evaluating partial derivatives, and rearranging Equation \(\ref{5.8.15}\) produces
Carry out the steps leading from Equation \(\ref{5.8.15}\) to Equation \(\ref{5.8.17}\). Keep in mind that, if \(y\) is not a function of \(x\),
\[ \dfrac {dy}{dx} = y \dfrac {d}{dx} \nonumber \]
Equation \(\ref{5.8.17}\) says that the function on the left, depending only on the variable \(\theta\), always equals the function on the right, depending only on the variable \(\varphi\), for all values of \(\theta\) and \(\varphi\). The only way two different functions of independent variables can be equal for all values of the variables is if both functions are equal to a constant (review separation of variables). We call this constant \(m_J^2\) because soon we will need the square root of it. The two differential equations to solve are the \(\theta\)-equation
The partial derivatives have been replaced by total derivatives because only a single variable is involved in each equation.
Often \(m_J\) is referred to as just \(m\) for convenience.
Exercise 5.8.2
Show how Equations \(\ref{5.8.18}\) and \(\ref{5.8.21}\) are obtained from Equation \(\ref{5.8.17}\).
Solving the \(\varphi\) Equation
The \(\varphi\)-equation is similar to the Schrödinger Equation for the free particle. Since we already solved this previously, we immediately write the solutions:
where we introduce the number \(m\) to track how many wavelengths of the wavefunction occur around one rotation (similar to the wavelength description of the Bohr atom).
Exercise 5.8.3
Substitute Equation \(\ref{5.8.22}\) into Equation \(\ref{5.8.21}\) to show that it is a solution to that differential equation.
The range of the integral is only from \(0\) to \(2π\) because the angle \(\varphi\) specifies the position of the internuclear axis relative to the x-axis of the coordinate system and angles greater than \(2π\) do not specify additional new positions.
Exercise 5.8.4
Use the normalization condition in Equation \(\ref{5.8.23}\) to demonstrate that \(N = 1/\sqrt{2π}\).
Answer
We need to evaluate Equation \ref{5.8.23} with \(\psi(\varphi)=N e^{\pm i m J \varphi} \)
\[\begin{align*} \psi^{*}(\varphi) \psi(\varphi) &= N e^{+i m J \varphi} N e^{-i m J \varphi} \\[4pt] &=N^{2} \\[4pt] 1=\int_{0}^{2 \pi} N^* N d \varphi=1 & \\[4pt] N^{2} (2 \pi) =1 \\[4pt] N=\sqrt{1 / 2 \pi} \end{align*} \nonumber \]
Values for \(m\) are found by using a cyclic boundary condition. The cyclic boundary condition means that since \(\varphi\) and \(\varphi + 2\varphi \) refer to the same point in three-dimensional space, \(\Phi (\varphi)\) must equal \(\Phi (\varphi + 2 \pi )\), i.e.
Finding the \(\Theta (\theta)\) functions that are solutions to the \(\theta\)-equation (Equation \(\ref{5.8.18}\)) is a more complicated process. Solutions are found to be a set of power series called Associated Legendre Functions (Table M2), which are power series of trigonometric functions, i.e., products and powers of sine and cosine functions. The \(\Theta (\theta)\) functions, along with their normalization constants, are shown in the third column of Table 5.8.1
.
The solution to the \(\theta\)-equation requires that \(λ\) in Equation \(\ref{5.8.17}\) be given by
\[\lambda = J (J + 1) \label {5.8.28} \]
where
\[ J \ge |m_J| \label {5.8.29} \]
\(J\) can be 0 or any positive integer greater than or equal to \(m_J\). Each pair of values for the quantum numbers, \(J\) and \(m_J\), identifies a rotational state with a wavefunction (Equation \(\ref{5.8.11}\)) and energy (below). Equation \(\ref{5.8.29}\) means that \(J\) controls the allowed values of \(m_J\).
Each pair of values for the quantum numbers, \(J\) and \(m_J\), identifies a rotational state and hence a specific wavefunction with associated energy.
The combination of Equations \(\ref{5.8.16}\) and \(\ref{5.8.28}\) reveals that the energy of this system is quantized.
Using Equation \(\ref{5.8.30}\), you can construct a rotational energy level diagram (Figure 5.8.2
). For simplicity, use energy units of \(\dfrac {\hbar ^2}{2I}\).
\(J=0\): The lowest energy state has \(J = 0\) and \(m_J = 0\). This state has an energy \(E_0 = 0\). There is only one state with this energy, i.e. one set of quantum numbers, one wavefunction, and one set of properties for the molecule.
\(J=1\): The next energy level is \(J = 1\) with energy \(\dfrac {2\hbar ^2}{2I}\). There are three states with this energy because \(m_J\) can equal +1, 0, or ‑1. These different states correspond to different orientations of the rotating molecule in space. States with the same energy are said to be degenerate. The degeneracy of an energy level is the number of states with that energy. The degeneracy of the \(J = 1\) energy level is 3 because there are three states with the energy \(\dfrac {2\hbar ^2}{2I}\).
\(J=2\): The next energy level is for \(J = 2\). The energy is \(\dfrac {6\hbar ^2}{2I}\), and there are five states with this energy corresponding to \(m_J = +2, \,+1,\, 0,\, ‑1,\, ‑2\). The energy level degeneracy is five. Note that the spacing between energy levels increases as J increases. Also note that the degeneracy increases. The degeneracy is always \(2J+1\) because \(m_J\) ranges from \(+J\) to \(‑J\) in integer steps, including 0.
Each allowed energy of rigid rotor is \((2J+1)\)-fold degenerate. Hence, there exist \((2J+1)\) different wavefunctions with that energy.
Exercise 5.8.6
Compute the energy levels for a rotating molecule for \(J = 0\) to \(J = 5\) using units of \(\dfrac {\hbar ^2}{2I}\).
Answer
This rotating molecule can be assumed to be a rigid rotor molecule. From solving the Schrödinger equation for a rigid rotor we have the relationship for energies of each rotational eigenstate (Equation \ref{5.8.30}):
\[E = J(J+1)(ħ^2/2I) \nonumber \]
Using this equation, we can plug in the different values of the \(J\) quantum number so that
For J=0, \(E = (0)(1)(ħ^2/2I) = 0\) For J=1, \(E = (1)(2)(ħ^2/2I) = 2(ħ^2/2I)\)
This shows that as \(J\) increases, the energy levels get farther apart (Figure 5.8.2
).
For J=2, \(E = (2)(3)(ħ^2/2I) = 6(ħ^2/2I)\)
For J=3, \(E = (3)(4)(ħ^2/2I) = 12(ħ^2/2I)\)
For J=4, \(E = (4)(5)(ħ^2/2I) = 20(ħ^2/2I)\)
For J=5, \(E = (5)(6)(ħ^2/2I) = 30(ħ^2/2I)\)
Exercise 5.8.7
For \(J = 0\) to \(J = 5\), identify the degeneracy of each energy level and the values of the \(m_J\) quantum number that go with each value of the \(J\) quantum number. Construct a rotational energy level diagram including \(J = 0\) through \(J=5\). Label each level with the appropriate values for the quantum numbers \(J\) and \(m_J\). Describe how the spacing between levels varies with increasing \(J\).
Answer
This rotating molecule can be assumed to be a rigid rotor molecule. From solving the Schrödinger equation for a rigid rotor, we have: \[λ = 2IE/ħ^2 \nonumber \]
Where λ is an arbitrary parameter not related to wavelength. Additionally, λ is assigned a quantum relation so that \[λ = J(J+1) \nonumber \]
Thus, combining the two equations and solving for E yields \[E = J(J+1)(ħ^2/2I) \nonumber \]
Using this equation, we can plug in the different values of J so that
For J=0, \(E = (0)(1)(ħ^2/2I) = 0\)
For J=1, \(E = (1)(2)(ħ^2/2I) = 2(ħ^2/2I)\)
For J=2, \(E = (2)(3)(ħ^2/2I) = 6(ħ^2/2I)\)
For J=3, \(E = (3)(4)(ħ^2/2I) = 12(ħ^2/2I)\)
For J=4, \(E = (4)(5)(ħ^2/2I) = 20(ħ^2/2I)\)
For J=5, \(E = (5)(6)(ħ^2/2I) = 30(ħ^2/2I)\)
This shows that as \(J\) increases, the energy levels get closer and closer together.
Interpretation of Quantum Numbers for a Rigid Rotor
The \(m_J\) quantum number reflects the component of the angular momentum along the \(z\) direction (and hence is sometimes called the azimuthal quantum number). For a fixed value of \(J\), the different values of \(m_J\) reflect the different directions the angular momentum vector could be pointing – for large, positive \(m_J\) the angular momentum is mostly along +z; if \(m_J\) is zero the angular momentum is orthogonal to \(z\). Physically, the energy of the rotation does not depend on the direction, which is reflected in the fact that the energy depends only on \(J\) (Equation \(\ref{5.8.30}\)), which measures the length of the vector, not its direction given mb \(m_J\).
Example 5.8.7
: Molecular Oxygen
Calculate \(J = 0\) to \(J = 1\) rotational transition of the \(\ce{O2}\) molecule with a bond length of 121 pm.
convert from atomic units to kilogram using the conversion: 1 au = 1.66 x 10-27 kg. Plug and chug.
\[E = 5.71 \times 10^{-27} \;Joules \nonumber \]
Spherical Harmonics
A wavefunction that is a solution to the rigid rotor Schrödinger Equation (Equation \(\ref{5.8.11}\)) can be written as a single function \(Y(\theta, \varphi)\), which is called a spherical harmonic function.
The spherical harmonic wavefunction is labeled with \(m_J\) and \(J\) because its functional form depends on both of these quantum numbers. These functions are tabulated above for \(J = 0\) through \(J = 2\) and for \(J = 3\) in the Spherical Harmonics Table (M4) Polar plots of some of the \(\theta\)-functions are shown in Figure 5.8.3
.
The two-dimensional space for a rigid rotor is defined as the surface of a sphere of radius \(r_0\), as shown in Figure 5.8.2
.
The probability of finding the internuclear axis at specific coordinates \(\theta _0\) and \(\varphi _0\) within an infinitesimal area \(ds\) on this curved surface is given by
where the area element \(ds\) is centered at \(\theta _0\) and \(\varphi _0\).
Within the Copenhagen interpretation of wavefunctions, the absolute square (or modulus squared) of the rigid rotor wavefunction \(Y^{m_{J*}}_J (\theta, \varphi) Y^{m_J}_J (\theta, \varphi) \) gives the probability density for finding the internuclear axis oriented at \(\theta\) to the z-axis and \(\varphi\) to the x-axis. In spherical coordinates the area element used for integrating \(\theta\) and \(\varphi\) is
Use calculus to evaluate the probability of finding the internuclear axis of a molecule described by the \(J = 1\), \(m_J = 0\) wavefunction somewhere in the region defined by a range in \(\theta\) of 0° to 45°, and a range in of 0° to 90°. Note that a double integral will be needed. Sketch this region as a shaded area on Figure 5.8.1
.
Consider the significance of the probability density function by examining the \(J = 1\), \(m_J = 0\) wavefunction. The Spherical Harmonic for this case is
The polar plot of \(( Y^0_1)^2\) is shown in Figure 5.8.1
. For \(J = 1\) and \(m_J = 0\), the probability of finding the internuclear axis is independent of the angle \(\varphi\) from the x-axis, and greatest for finding the internuclear axis along the z‑axis, but there also is a probability for finding it at other values of \(\theta\) as well. So, although the internuclear axis is not always aligned with the z-axis, the probability is highest for this alignment. Also, since the probability is independent of the angle \(\varphi\), the internuclear axis can be found in any plane containing the z-axis with equal probability.
The \(J = 1\), \(m_J = 0\) function is 0 when \(\theta\) = 90°. Therefore, the entire xy-plane is a node. This fact means the probability of finding the internuclear axis in this particular horizontal plane is 0 in contradiction to our classical picture of a rotating molecule. In the classical picture, a molecule rotating in a plane perpendicular to the xy‑plane must have the internuclear axis lie in the xy‑plane twice every revolution, but the quantum mechanical description says that the probability of being in the xy-plane is zero. This conclusion means that molecules are not rotating in the classical sense, but they still have some, but not all, of the properties associated with classical rotation. The properties they retain are associated with angular momentum.
Exercise 5.8.9
For each state with \(J = 0\) and \(J = 1\), use the function form of the \(Y\) spherical harmonics and Figure 5.8.1
to determine the most probable orientation of the internuclear axis in a diatomic molecule, i.e., the most probable values for \(\theta\) and \(\theta\).
Exercise 5.8.10
Write a paragraph describing the information about a rotating molecule that is provided in the polar plot of \(Pr [\theta, \theta ] \) for the \(J = 1\), \(m_J = \pm 1\) state in Figure 5.8.1
. Compare this information to the classical picture of a rotating object.
Summary
There are two quantum numbers that describe the quantum behavior of a rigid rotor in three-deminesions: \(J\) is the total angular momentum quantum number and \(m_J\) is the z-component of the angular momentum. The spherical harmonics called \(Y_J^{m_J}\) are functions whose probability \(|Y_J^{m_J}|^2\) has the well known shapes of the s, p and d orbitals etc learned in general chemistry.