The Rigid Rotor


Another form of microscopic motion that shows interesting quantization effects is rotation. Rotation is a fundamental motion of molecules and atoms; for example, gaseous molecules are in a constant state of relatively free rotation, electrons in atoms orbit around the nucleus, and the nuclear particles themselves exhibit spin. Rotational motion is characterized by an angular momentum vector L. On a microscopic scale the angular momentum is quantitized in both length and direction by units of Panck's constant divided by 2¹

This quantization of angular momentum is a fundamental property of quantum mechanical systems. The central role of angular momentum in quantum mechanics is indicated by the fact that Planck's constant itself has units of angular momentum. In the early days of quantum mechanics Niels Bohr developed a successful model of the hydrogen atom spectrum by imposing a quantization condition on the electron's orbital angular momentum. Indeed many features of molecular and atomic spectra results from the fact that the total angular momentum can only take on certain values. Today, we see that quantization of angular momentum does not need to be imposed on a system but is a natural consequence of Schrödinger's model.

The simplest model of rotational motion that can be described quantum mechanically is a rigid rotor consisting of two masses held a fixed distance apart rotating about its center of mass. The rotation is assumed to be free of any outside potential energy. Angular momentum in quantum systems results from a direct correspondence of classical mechanical equations of rotational motion. In rotational motion the moment of inertia and angular velocity play the same role as the mass and velocity does in a straight moving body. In linear motion momentum and kinetic energy are given by p = mv and E =p2/2m , so a rigid rotor rotating with angular velocity w and having a moment of inertia I = µ R2 exhibits angular momentum |L| = I w; and an energy of E = L2/2I. To convert these classical equations to quantum equations we simply need to replace the momentum with its quantum operator and solve Schrödinger's eigenvalue equation.

In most respects the analysis of rotational systems is largely a generalization of the types of coordinates used to describe the system. Schrödinger's eigenvalue equation given above is very hard to solve in Cartesian coordinates because motions in the x,y, and z directions are not independent of each other. Polar coordinates most directly describe rotational motion and allow the Hamiltonian to be separated into independent coordinates. For example, the angular velocity is only dependent on the time derivative of the phi coordinate. To solve Schrödinger's equation we need to convert the Hamiltonian to polar coordinates. Chain rule differentiation provides the means for converting differential operators from Cartesian to polar coordinates.

Angular momentum is a vector quantity that results from the cross-product of the position vector r from the center of rotation with the linear momentum vector p of the particles in motion. Conversion of the angular momentum vector to polar coordinates is given in the following table.

For a rigid rotor r is a constant and the Hamiltonian becomes

The wavefunction for this Hamiltonian will be two-dimensional depending on both coordinates. To solve this problem we need to reduce the two-dimensional problem to two one-dimensional problems. Because the derivatives in the Hamiltonian are in respects to only one coordinate at a time the Hamiltonian is separable and we can look for a solution of the form of a product function:

Y(q,f) = Q(q)F(f)

separating like variables gives

For this equation to be true for all variations of both angles, it must be true that both sides are equal to a constant

It is easy to see by inspection that the first equation has a solution of

Differentiating twice shows the constant beta is equal to m2. Furthermore, because polar coordinates are cyclic, the wavefunction must be invariant on a rotation of 2¹.

F(f + p) = F(f)

substituting this condition into the wavefunction gives:

F(f + 2p) = Aeim(f + 2p) = A eimf eim2p = F(f) ei p2m = F(f) (-1)2m

Using this in the cyclic boundary condition gives:

(-1)2m = 1

Hence, 2m must be a positive or negative even integer restricting m to values: 0, ±1, ±2, ±3 .....

It is easy to see that this solution is also an eigenfunction of the z component of the angular momentum.

The remarkable implication of this equation is that the orientation of a rotating body is quantitized. A macroscopic object like a child's top starts out spinning rapidly in a vertical direction then begins to precess, its axis pointing in an ever increasing angle to the vertical until it wobbles to the floor. In quantum mechanical systems the angular momentum vector is not allowed to point in any direction. A quantum mechanical top would only precess at fixed angles to the vertical. Angular momentum is space quantitized such that z component of angular momentum has values of only m. The quantum number m is called the magnetic quantum number because it describes how spinning particles such as electrons or protons will align with an external magnetic field.

The second differential equation can now be used to solve for the energy.

In this program we will solve this differential equation numerically varing the energy until the cyclic boundary condition is satisfied. Note: this is not an easy differential equation to solve because the sine function cause singularities at 0 and ¹ . Small numerical errors that keep the solution from being zero at these points will be magnified to infinity. You might guess then that the analytical solutions would contain the sine function to keep these singularities in check.

This is a well known Legendre differential equation whose analytical solution are called associate Legendre functions. The cyclic boundary condition requires that the energy is equal to

The first 4 sets of associate Legendre functions are:

Q(q) = Pl|m|
|m| = 0|m| = 1 |m| = 2|m| = 3
l = 0P00 = 1
l = 1P10 = cosq P11 = sinq
l = 2P20 = (3cos2q -1)/2 P21 = 3 cosq sinq P22 = 3 sin2q
l = 3P30 = cosq(5cos2q -3)/2 P30 = 3sinq(5cos2q -1)/2 P31 = 15 cosq sin2q P33 = 15 sin3q

The associate Legendre functions are normalized in polar coordinate representation

The normalized product function solution of the full rotational Hamiltonian is called a spherical harmonic (vibrational modes of a sphere of jello) and is usually denoted as Ylm

The spherical harmonics are eigenfunctions of the L2 operator

The implication of this equation is again startling, the magnitude and hence the angular velocity is quantitized. In the microscopic world a child's yo-yo could only be set spinning a certain rates. One of these rates is zero. At first glance this seems to violate the uncertainty principle. Most systems do not have zero point energies that are zero. A zero velocity does not violate the uncertainty principle in this case because the angle and hence its position is completely undermined. Also note that there is a high degree of degeneracy. For a given l quantum number there are 2l+1 different m quantum numbers. For example when l=2 there will be five different m values corresponding to five different orientations of the angular momentum with respects to the z direction. Each will have the same rotation rate and energy because the length of the angular momentum vector has the same magnitude by different orientations.


An atom is also a spherical system where the total orbital angular momentum is quantitized. The spherical harmonics give the atomic orbitals their angular shape. For example, l=2 corresponds to a d orbital. Their are five d orbitals each corresponding to a different m quantum number. These five d orbitals are normally degenerate but in the presence of a magnetic field they will separate into five different energies corresponding to the orientation of the angular momentum vector with the field.