Rigid rotor: Difference between revisions
imported>Paul Wormer (→Arbitrarily shaped rigid rotor: Worked through it again) |
mNo edit summary |
||
(12 intermediate revisions by 3 users not shown) | |||
Line 1: | Line 1: | ||
In [[classical mechanics]], | {{subpages}} | ||
In [[classical mechanics]] and [[quantum mechanics]], a '''rigid rotor''' is a 3-dimensional rigid body, such as a [[top]]—a children's toy. To describe the orientation of a 3-dimensional object with respect to a 3-dimensional coordinate system three angles are required. If the angles do not vary in time, the rigid body is standing still; when the angles vary in time the rigid body is rotating and is referred to as a '''rigid rotor''' (also known as '''rigid rotator'''). | |||
In [[molecular physics]], a molecule can often be regarded as a rigid rotor. An example is a [[diatomic molecule]] | A special rigid rotor is the linear rotor which is a 2-dimensional object (a pencil), requiring two angles to describe its orientation. Sometimes one considers rotation of a single particle around an axis fixed in space, then one angle is sufficient to describes the system fully. | ||
the water molecule (H<sub>2</sub>O) is an asymmetric rotor, the ammonia molecule (NH<sub>3</sub>) is a symmetric rotor, and methane (CH<sub>4</sub>) is spherical rotor. | |||
In [[molecular physics]], a molecule can often be regarded as a rigid rotor, which is why the (quantum mechanical) theory of rigid rotors plays an important role in molecular [[spectroscopy]]. An example of a linear rotor is a [[diatomic molecule]]; if one neglects its vibration, the diatomic molecule is a rigid linear rotor. More general molecules, too, can often be seen as rigid, i.e., often their vibration can be ignored. Rigid rotors can be classified by means of their [[inertia moment]]s, see [[classification of rigid rotors]]. For instance the water molecule (H<sub>2</sub>O) is an asymmetric rotor, the ammonia molecule (NH<sub>3</sub>) is a symmetric rotor, and methane (CH<sub>4</sub>) is spherical rotor. | |||
This article is restricted to the ''rotational kinematics'' of rigid bodies, that is, this article is about the kinetic energy of rotating bodies as function of time. It is shown that differentiation of the kinetic energy with respect to appropriate degrees of freedom gives rise to the velocity, acceleration, and angular momentum of the rigid rotor. ''Rotational dynamics'' (effect of forces, such as gravity, on the motion of the body) is not considered in this article. | |||
== The linear rotor == | == The linear rotor == | ||
Line 81: | Line 85: | ||
\tfrac{1}{2I}\; \hat{l}^2\;Y_\ell^m(\theta, \varphi ) = \tfrac{1}{2I}\; \hbar^2\ell(\ell+1) \;Y_\ell^m(\theta, \varphi ). | \tfrac{1}{2I}\; \hat{l}^2\;Y_\ell^m(\theta, \varphi ) = \tfrac{1}{2I}\; \hbar^2\ell(\ell+1) \;Y_\ell^m(\theta, \varphi ). | ||
</math> | </math> | ||
Here <math>Y_\ell^m (\theta, \varphi )</math> represents a set of functions known as the [[spherical harmonic]]s. Note that the energy does not depend on <math>m \,</math> | Here <math>Y_\ell^m (\theta, \varphi )</math> represents a set of functions known as the [[spherical harmonic]]s. Note that the energy does not depend on the magnetic quantum number <math>m \,</math>, in other words, the energy | ||
:<math> | :<math> | ||
E_\ell \equiv {\hbar^2 \over 2I} \ell \left (\ell+1\right ) | E_\ell \equiv {\hbar^2 \over 2I} \ell \left (\ell+1\right ) | ||
Line 94: | Line 98: | ||
:<math> \bar B \equiv \frac{B}{hc} = \frac{h}{8\pi^2cI}, </math> | :<math> \bar B \equiv \frac{B}{hc} = \frac{h}{8\pi^2cI}, </math> | ||
with ''c'' the speed of light. If [[cgs units]] are used for ''h'', ''c'', and ''I'', <math>\bar B</math> is expressed | with ''c'' the speed of light. If [[cgs units]] are used for ''h'', ''c'', and ''I'', <math>\bar B</math> is expressed | ||
in [[Wavenumber|wave numbers]], cm<sup> | in [[Wavenumber|wave numbers]], cm<sup>−1</sup>, a unit that is often used in rotational-vibrational spectroscopy. | ||
Obviously, the rotational constant <math>\bar B</math> depends on the distance ''R''. Often one writes <math> B_e = \bar B(R_e) </math> where ''R''<sub>''e''</sub> is the equilibrium value | Obviously, the rotational constant <math>\bar B</math> depends on the distance ''R''. Often one writes <math> B_e = \bar B(R_e) </math> where ''R''<sub>''e''</sub> is the equilibrium value | ||
Line 111: | Line 115: | ||
Different branches of physics and engineering use different coordinates for the description | Different branches of physics and engineering use different coordinates for the description | ||
of the kinematics of a rigid rotor. In molecular physics [[Euler angles]] are used almost exclusively. In quantum mechanical applications it is advantageous to use Euler | of the kinematics of a rigid rotor. In molecular physics [[Euler angles]] are used almost exclusively. In quantum mechanical applications it is advantageous to use Euler | ||
angles in a convention that is a simple extension of [[ | angles in a convention that is a simple extension of [[spherical polar coordinates]]. | ||
The first step is the attachment of a | The first step is the attachment of a ''right-handed'' orthonormal frame (3-dimensional system of orthogonal axes) to the rotor (a ''body-fixed frame'') . This frame can be attached arbitrarily to the body, but often one uses the principal axes frame—the normalized eigenvectors of the inertia tensor, which always can be chosen orthonormal, since the tensor is [[Symmetric matrix|symmetric]]. When the rotor possesses a symmetry-axis, it usually coincides with one of the principal axes. It is convenient to choose | ||
as body-fixed ''z''-axis the highest-order symmetry axis. | as body-fixed ''z''-axis the highest-order symmetry axis. | ||
One starts by aligning the body-fixed frame with a | One starts by aligning the body-fixed frame with a ''space-fixed frame'' | ||
(laboratory | (laboratory frame), so that the body-fixed ''x'', ''y'', and ''z'' axes coincide with the space-fixed ''X'', ''Y'', and ''Z'' axis. Then the body and its frame are rotated ''actively'' over a ''positive angle'' α around the ''z''-axis (by the [[right-hand screw rule]]), which moves the <math>y</math>- to the <math>y'</math>-axis. Thirdly, one rotates the body and its frame over a positive angle β around the <math>y'</math>-axis. The ''z''-axis of the body-fixed frame has after these two rotations the longitudinal angle α (commonly designated by φ) and the colatitude angle β (commonly designated by θ), both with respect to the space-fixed frame. If the rotor were cylindrical symmetric around its body-fixed ''z''-axis, like the linear rigid rotor, its orientation in space would be unambiguously specified at this point. | ||
If the body lacks cylinder (axial) symmetry, a last rotation around its ''z''-axis (which has polar coordinates β and α) is necessary to specify its orientation completely. Traditionally the last rotation angle is called γ. | If the body lacks cylinder (axial) symmetry, a last rotation around its ''z''-axis (which has polar coordinates β and α) is necessary to specify its orientation completely. Traditionally the last rotation angle is called γ. | ||
The [[Euler angles | The [[Euler angles|convention for Euler angles]] used here is known as the <math>z'-y'-z</math> convention; it can be shown that it is equivalent to the <math>z-y-z</math> convention in which the order of rotations is reversed. | ||
The total matrix of the three consecutive rotations is the product | The total matrix of the three consecutive rotations is the product | ||
Line 164: | Line 168: | ||
\end{pmatrix}, | \end{pmatrix}, | ||
</math> | </math> | ||
which shows the correspondence with the [[ | which shows the correspondence with the [[spherical polar coordinates]]. | ||
Knowledge of the Euler angles as function of time ''t'' and the initial coordinate vector '''r'''(0) determine the kinematics of the rigid rotor. | Knowledge of the Euler angles as function of time ''t'' and the initial coordinate vector '''r'''(0) determine the kinematics of the rigid rotor. | ||
Line 246: | Line 249: | ||
===== Angular momentum form ===== | ===== Angular momentum form ===== | ||
Often the kinetic energy is written as a function of the [[Angular momentum | Often the kinetic energy is written as a function of the [[Angular momentum (classical)|angular momentum]] <math>\vec{L}</math> of the rigid rotor. This vector is a conserved (time-independent) quantity. | ||
With respect to the body-fixed frame it has the components <math>\mathbf{L}</math>, which can be shown to | With respect to the body-fixed frame it has the components <math>\mathbf{L}</math>, which can be shown to | ||
be related to the angular velocity, | be related to the angular velocity, | ||
Line 346: | Line 349: | ||
The quantization rule is ''not'' sufficient to obtain the kinetic energy operator from | The quantization rule is ''not'' sufficient to obtain the kinetic energy operator from | ||
the classical Hamiltonian. Since classically <math>p_\beta</math> commutes with <math>\cos\beta</math> and <math>\sin\beta</math> and the reciprocals of these functions, the position of these trigonometric functions in the classical Hamiltonian is arbitrary. After | the classical Hamiltonian. Since classically <math>p_\beta</math> commutes with <math>\cos\beta</math> and <math>\sin\beta</math> and the reciprocals of these functions, the position of these trigonometric functions in the classical Hamiltonian is arbitrary. After | ||
quantization the commutation no longer holds and the order of operators and functions in the Hamiltonian (energy operator) becomes a point of concern. Podolsky<ref name="Podolsky" /> proposed in 1928 that the [[ | quantization the commutation no longer holds and the order of operators and functions in the Hamiltonian (energy operator) becomes a point of concern. Podolsky<ref name="Podolsky" /> proposed in 1928 that the [[Laplace-Beltrami operator]] (times <math>-\tfrac{1}{2}\hbar^2</math>) has the appropriate form for the quantum mechanical kinetic energy operator. This operator has the general form (summation convention: sum over repeated indices—in this case over the three Euler angles ''q''<sup>1</sup> ≡ α, | ||
(times <math>-\tfrac{1}{2}\hbar^2</math>) has the appropriate form for the quantum mechanical kinetic energy operator. This operator has the general form (summation convention: sum over repeated indices—in this case over the three Euler angles ''q''<sup>1</sup> ≡ α, | |||
''q''<sup>2</sup> ≡ β, ''q''<sup>3</sup> ≡ γ), | ''q''<sup>2</sup> ≡ β, ''q''<sup>3</sup> ≡ γ), | ||
Line 360: | Line 362: | ||
Given the inverse of the metric tensor above, the explicit form of the kinetic energy operator in terms of Euler angles follows by simple substitution. | Given the inverse of the metric tensor above, the explicit form of the kinetic energy operator in terms of Euler angles follows by simple substitution. | ||
It can be shown that <math>\hat{H}</math> can be expressed in body-fixed angular momentum operators (in this proof one must carefully commute differential operators with trigonometric functions). The result has the same appearance as the classical formula expressed in body-fixed | It can be shown that <math>\hat{H}</math> can be expressed in body-fixed angular momentum operators (in this proof one must carefully commute differential operators with trigonometric functions). The result has the same appearance as the classical formula expressed in body-fixed coordinates, | ||
:<math> | :<math> | ||
\hat{H} = \tfrac{1}{2}\left[ \frac{\mathcal{P}_x^2}{I_1}+ \frac{\mathcal{P}_y^2}{I_2}+ | \hat{H} = \tfrac{1}{2}\left[ \frac{\mathcal{P}_x^2}{I_1}+ \frac{\mathcal{P}_y^2}{I_2}+ | ||
Line 383: | Line 385: | ||
the [[Schrödinger equation]] as it was solved for the first time by Kronig and Rabi.<ref name="Kronig">R. de L. Kronig and I. I. Rabi, Phys. Rev., vol. 29, pp. 262-269 (1927).</ref> The symmetric top is one of the few cases where the Schrödinger equation can be solved analytically. All these cases were solved within a year of the formulation of the Schrödinger equation. | the [[Schrödinger equation]] as it was solved for the first time by Kronig and Rabi.<ref name="Kronig">R. de L. Kronig and I. I. Rabi, Phys. Rev., vol. 29, pp. 262-269 (1927).</ref> The symmetric top is one of the few cases where the Schrödinger equation can be solved analytically. All these cases were solved within a year of the formulation of the Schrödinger equation. | ||
Knowing | Knowing the form of the Hamiltonian in terms of angular momentum operators, and knowing the relation for Wigner D-matrices, | ||
:<math> | :<math> | ||
\mathcal{P}_z^2\, D^j_{m k}(\alpha,\beta,\gamma)^* = k^2\, D^j_{m k}(\alpha,\beta,\gamma)^* | \mathcal{P}_z^2\, D^j_{m k}(\alpha,\beta,\gamma)^* = \hbar^2 k^2\, D^j_{m k}(\alpha,\beta,\gamma)^*, | ||
</math> | </math> | ||
the [[Schrödinger equation]] does not any longer have to be solved from scratch, as Kronig and Rabi had to do. | |||
Indeed, | |||
:<math> | :<math> | ||
\hat{H}\,D^j_{m k}(\alpha,\beta,\gamma)^* = E_{jk} D^j_{m k}(\alpha,\beta,\gamma)^* | \hat{H}\,D^j_{m k}(\alpha,\beta,\gamma)^* = E_{jk} D^j_{m k}(\alpha,\beta,\gamma)^* | ||
\quad \hbox{with}\quad E_{jk} = \frac{j(j+1)}{2I_1} + k^2\left(\frac{1}{2I_3}-\frac{1}{2I_1}\right). | \quad \hbox{with}\quad E_{jk} = \hbar^2\frac{j(j+1)}{2I_1} + \hbar^2 k^2\left(\frac{1}{2I_3}-\frac{1}{2I_1}\right). | ||
</math> | </math> | ||
The eigenvalue ''E''<sub>''j''0</sub> is (2''j'' + 1)-fold degenerate, | The eigenvalue ''E''<sub>''j''0</sub> is (2''j'' + 1)-fold degenerate, because all eigenfunctions | ||
with ''m'' = − ''j'', −''j'' + 1, ..., ''j'' have the same eigenvalue. | with ''m'' = − ''j'', −''j'' + 1, ..., ''j'' have the same eigenvalue. Since ''k'' appears squared in the energy expression, the energies with k ≠ 0 are 2(2''j'' + 1)-fold degenerate. The complex conjugates of the Wigner D-matrix elements are the corresponding eigenfunctions, that is, they are the wave functions of the symmetric top. | ||
The asymmetric top problem (<math> I_1 \ne I_2 \ne I_3 </math>) is not exactly soluble. | The asymmetric top problem (<math> I_1 \ne I_2 \ne I_3 </math>) is not exactly soluble. | ||
==References== | ==References== | ||
{{reflist}}[[Category:Suggestion Bot Tag]] | |||
Latest revision as of 11:01, 12 October 2024
In classical mechanics and quantum mechanics, a rigid rotor is a 3-dimensional rigid body, such as a top—a children's toy. To describe the orientation of a 3-dimensional object with respect to a 3-dimensional coordinate system three angles are required. If the angles do not vary in time, the rigid body is standing still; when the angles vary in time the rigid body is rotating and is referred to as a rigid rotor (also known as rigid rotator).
A special rigid rotor is the linear rotor which is a 2-dimensional object (a pencil), requiring two angles to describe its orientation. Sometimes one considers rotation of a single particle around an axis fixed in space, then one angle is sufficient to describes the system fully.
In molecular physics, a molecule can often be regarded as a rigid rotor, which is why the (quantum mechanical) theory of rigid rotors plays an important role in molecular spectroscopy. An example of a linear rotor is a diatomic molecule; if one neglects its vibration, the diatomic molecule is a rigid linear rotor. More general molecules, too, can often be seen as rigid, i.e., often their vibration can be ignored. Rigid rotors can be classified by means of their inertia moments, see classification of rigid rotors. For instance the water molecule (H2O) is an asymmetric rotor, the ammonia molecule (NH3) is a symmetric rotor, and methane (CH4) is spherical rotor.
This article is restricted to the rotational kinematics of rigid bodies, that is, this article is about the kinetic energy of rotating bodies as function of time. It is shown that differentiation of the kinetic energy with respect to appropriate degrees of freedom gives rise to the velocity, acceleration, and angular momentum of the rigid rotor. Rotational dynamics (effect of forces, such as gravity, on the motion of the body) is not considered in this article.
The linear rotor
The classical as well as the quantum mechanical linear rigid rotor model consists of two point masses located at fixed distances from their center of mass. The fixed distance between the two point masses and the values of the masses are the only characteristics of the rigid model.
The classical linear rigid rotor
The two point masses are denoted by m1 and m2 [with reduced mass μ ≡ m1 m2/(m1 + m2)] and are kept at a fixed distance R. Classically, the rotor is rigid if R is independent of time. The kinematics of a linear rigid rotor is usually described by means of spherical polar coordinates, which form a coordinate system of 3-dimensional space. In the physics convention the coordinates are the colatitude (zenith) angle θ, the longitudinal (azimuth) angle φ and the distance R. The angles specify the orientation of the rotor in space. The kinetic energy T of the linear rigid rotor in Lagrangian form (in terms of coordinates and time derivatives of coordinates) is:
The classical Hamiltonian function of the linear rigid rotor is:
It is possible to write the Lagrangian in matrix form:
where hθ ≡ R and hφ ≡ Rsinθ are scale (or Lamé) factors.
Scale factors are of importance for quantum mechanical applications since they enter the Laplacian expressed in curvilinear coordinates. In the case at hand (constant R)
The quantum mechanical linear rigid rotor
The linear rigid rotor model can be used in quantum mechanics to predict the rotational energy of a diatomic molecule. The rotational energy depends on the moment of inertia for the system, I. In the center of mass reference frame, the moment of inertia is equal to:
where μ is the reduced mass of the molecule and R is the distance between the two atoms.
According to quantum mechanics, the energy levels of a system can be determined by solving the Schrödinger equation H Ψ = E Ψ, where Ψ is the wave function and H is the energy (Hamiltonian) operator. For the rigid rotor in a field-free space, the energy operator corresponds to the kinetic energy[1] of the system:
where is Planck's constant divided by and is the Laplacian. The Laplacian is given above in terms of spherical polar coordinates. The energy operator written in terms of these coordinates is:
This operator appears also in the (θ, φ)-part of the Schrödinger equation of the hydrogen-like atoms,
The eigenvalue equation becomes
Here represents a set of functions known as the spherical harmonics. Note that the energy does not depend on the magnetic quantum number , in other words, the energy
is -fold degenerate: the functions with fixed and have the same energy.
In microwave spectroscopy one prefers using the rotational constant B instead of the inertia moment, and one writes,
In the unit of reciprocal length the rotational constant is,
with c the speed of light. If cgs units are used for h, c, and I, is expressed in wave numbers, cm−1, a unit that is often used in rotational-vibrational spectroscopy.
Obviously, the rotational constant depends on the distance R. Often one writes where Re is the equilibrium value of R (the value for which the interaction energy of the atoms in the rotor has a minimum).
Arbitrarily shaped rigid rotor
An arbitrarily shaped rigid rotor is a rigid body of arbitrary shape with its center of mass fixed (or in uniform rectilinear motion) in field-free space , so that its energy consists only of rotational kinetic energy (plus possibly a constant translational energy that can be ignored). A rigid body can be (partially) characterized by the three eigenvalues of its moment of inertia tensor, which are real nonnegative values known as principal moments of inertia. In microwave spectroscopy—the spectroscopy based on rotational transitions—one usually classifies molecules (seen as rigid rotors) as follows:
- spherical rotors
- symmetric rotors
- oblate symmetric rotors
- prolate symmetric rotors
- asymmetric rotors
This classification depends on the relative magnitudes of the principal moments of inertia.
Coordinates of the rigid rotor
Different branches of physics and engineering use different coordinates for the description of the kinematics of a rigid rotor. In molecular physics Euler angles are used almost exclusively. In quantum mechanical applications it is advantageous to use Euler angles in a convention that is a simple extension of spherical polar coordinates.
The first step is the attachment of a right-handed orthonormal frame (3-dimensional system of orthogonal axes) to the rotor (a body-fixed frame) . This frame can be attached arbitrarily to the body, but often one uses the principal axes frame—the normalized eigenvectors of the inertia tensor, which always can be chosen orthonormal, since the tensor is symmetric. When the rotor possesses a symmetry-axis, it usually coincides with one of the principal axes. It is convenient to choose as body-fixed z-axis the highest-order symmetry axis.
One starts by aligning the body-fixed frame with a space-fixed frame (laboratory frame), so that the body-fixed x, y, and z axes coincide with the space-fixed X, Y, and Z axis. Then the body and its frame are rotated actively over a positive angle α around the z-axis (by the right-hand screw rule), which moves the - to the -axis. Thirdly, one rotates the body and its frame over a positive angle β around the -axis. The z-axis of the body-fixed frame has after these two rotations the longitudinal angle α (commonly designated by φ) and the colatitude angle β (commonly designated by θ), both with respect to the space-fixed frame. If the rotor were cylindrical symmetric around its body-fixed z-axis, like the linear rigid rotor, its orientation in space would be unambiguously specified at this point.
If the body lacks cylinder (axial) symmetry, a last rotation around its z-axis (which has polar coordinates β and α) is necessary to specify its orientation completely. Traditionally the last rotation angle is called γ.
The convention for Euler angles used here is known as the convention; it can be shown that it is equivalent to the convention in which the order of rotations is reversed.
The total matrix of the three consecutive rotations is the product
Let r(0) be the coordinate vector with respect to the body-fixed frame of an arbitrary point ℘ of the body. We call the three elements of r(0) the "body-fixed coordinates" of ℘. Initially r(0) is also the space-fixed coordinate vector of ℘. Upon rotation of the body (including its frame and its point ℘), the body-fixed coordinates of ℘ do not change, but the space-fixed coordinate vector of ℘ becomes,
In particular, if ℘ is initially on the space-fixed Z-axis, it has the space-fixed coordinates
which shows the correspondence with the spherical polar coordinates.
Knowledge of the Euler angles as function of time t and the initial coordinate vector r(0) determine the kinematics of the rigid rotor.
Classical kinetic energy
The following is a generalization of the well-known special case of the rotational energy of an object that rotates around one axis.
It will be assumed from here on that the body-fixed frame is a principal axes frame; it diagonalizes the instantaneous inertia tensor (expressed with respect to the space-fixed frame), i.e.,
where the Euler angles are time-dependent and in fact determine the time dependence of by the inverse of this equation. This notation implies that at the Euler angles are zero, so that at the body-fixed frame coincides with the space-fixed frame.
The classical kinetic energy T of the rigid rotor can be expressed in different ways:
- as a function of angular velocity
- in Lagrangian form
- as a function of angular momentum
- in Hamiltonian form.
Since each of these forms has its use and can be found in textbooks we will present all of them.
Angular velocity form
As a function of angular velocity T reads,
with
The vector ω = (ωx, ωy, ωz) contains the components of the angular velocity of the rotor expressed with respect to the body-fixed frame. It can be shown that ω is not the time derivative of any vector, in contrast to the usual definition of velocity. The dots over the time-dependent Euler angles indicate time derivatives.
Lagrange form
Backsubstitution of the expression of ω into T gives the kinetic energy in Lagrange form (as a function of the time derivatives of the Euler angles). In matrix-vector notation,
where is the metric tensor associated with Euler angles—a non-orthogonal system of curvilinear coordinates—
Angular momentum form
Often the kinetic energy is written as a function of the angular momentum of the rigid rotor. This vector is a conserved (time-independent) quantity. With respect to the body-fixed frame it has the components , which can be shown to be related to the angular velocity,
Since the body-fixed frame moves (depends on time) these components are not time independent. If we were to represent with respect to the stationary space-fixed frame, we would find time independent expressions for its components. The kinetic energy is given by
Hamilton form
The Hamilton form of the kinetic energy is written in terms of generalized momenta
where it is used that the is symmetric. In Hamilton form the kinetic energy is,
with the inverse metric tensor given by
This inverse tensor is needed to obtain the Laplace-Beltrami operator, which (multiplied by ) gives the quantum mechanical energy operator of the rigid rotor.
The classical Hamiltonian given above can be rewritten to the following expression, which is needed in the phase integral arising in the classical statistical mechanics of rigid rotors,
Quantum mechanical rigid rotor
As usual quantization is performed by the replacement of the generalized momenta by operators that give first derivatives with respect to its canonically conjugate variables (positions). Thus,
and similarly for and . It is remarkable that this rule replaces the fairly complicated function of all three Euler angles, time derivatives of Euler angles, and inertia moments (characterizing the rigid rotor) by a simple differential operator that does not depend on time or inertia moments and differentiates to one Euler angle only.
The quantization rule is sufficient to obtain the operators that correspond with the classical angular momenta. There are two kinds: space-fixed and body-fixed angular momentum operators. Both are vector operators, i.e., both have three components that transform as vector components among themselves upon rotation of the space-fixed and the body-fixed frame, respectively. The explicit form of the rigid rotor angular momentum operators is given here (but beware, they must be multiplied with ). The body-fixed angular momentum operators are written as . They satisfy anomalous commutation relations.
The quantization rule is not sufficient to obtain the kinetic energy operator from the classical Hamiltonian. Since classically commutes with and and the reciprocals of these functions, the position of these trigonometric functions in the classical Hamiltonian is arbitrary. After quantization the commutation no longer holds and the order of operators and functions in the Hamiltonian (energy operator) becomes a point of concern. Podolsky[1] proposed in 1928 that the Laplace-Beltrami operator (times ) has the appropriate form for the quantum mechanical kinetic energy operator. This operator has the general form (summation convention: sum over repeated indices—in this case over the three Euler angles q1 ≡ α, q2 ≡ β, q3 ≡ γ),
where is the determinant of the g-tensor:
Given the inverse of the metric tensor above, the explicit form of the kinetic energy operator in terms of Euler angles follows by simple substitution.
It can be shown that can be expressed in body-fixed angular momentum operators (in this proof one must carefully commute differential operators with trigonometric functions). The result has the same appearance as the classical formula expressed in body-fixed coordinates,
The action of the on the Wigner D-matrix is simple. In particular
so that the Schrödinger equation for the spherical rotor (I1 = I2 = I3) is solved with the (2j+1)2 degenerate energy equal to and Wigner D-matrix elements as eigenfunctions.
Schrödinger equation for symmetric top
The symmetric top (= symmetric rotor) is characterized by two equal principal moments of inertia. It is a prolate (cigar shaped) top if I3 < I1 = I2. In this case we write the Hamiltonian as
Introducing the angular momentum operators in their differential operator form, we obtain the Schrödinger equation as it was solved for the first time by Kronig and Rabi.[2] The symmetric top is one of the few cases where the Schrödinger equation can be solved analytically. All these cases were solved within a year of the formulation of the Schrödinger equation.
Knowing the form of the Hamiltonian in terms of angular momentum operators, and knowing the relation for Wigner D-matrices,
the Schrödinger equation does not any longer have to be solved from scratch, as Kronig and Rabi had to do.
Indeed,
The eigenvalue Ej0 is (2j + 1)-fold degenerate, because all eigenfunctions with m = − j, −j + 1, ..., j have the same eigenvalue. Since k appears squared in the energy expression, the energies with k ≠ 0 are 2(2j + 1)-fold degenerate. The complex conjugates of the Wigner D-matrix elements are the corresponding eigenfunctions, that is, they are the wave functions of the symmetric top.
The asymmetric top problem () is not exactly soluble.
References
- ↑ Jump up to: 1.0 1.1 B. Podolsky, Phys. Rev., vol. 32, p. 812 (1928)
- ↑ R. de L. Kronig and I. I. Rabi, Phys. Rev., vol. 29, pp. 262-269 (1927).