Rigid rotor

Last updated

In rotordynamics, the rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water (asymmetric rotor), ammonia (symmetric rotor), or methane (spherical rotor).

Contents

Linear rotor

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. 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 point of departure (zeroth-order model).

Classical linear rigid rotor

The classical linear rotor consists of two point masses and (with reduced mass ) at a distance of each other. The rotor is rigid if 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 R3. In the physics convention the coordinates are the co-latitude (zenith) angle , the longitudinal (azimuth) angle and the distance . The angles specify the orientation of the rotor in space. The kinetic energy of the linear rigid rotor is given by

where and 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 )

The classical Hamiltonian function of the linear rigid rotor is

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, . In the center of mass reference frame, the moment of inertia is equal to:

where is the reduced mass of the molecule and 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:

where is the wave function and 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 reduced Planck constant 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 Schrödinger equation of the hydrogen atom after the radial part is separated off. The eigenvalue equation becomes The symbol represents a set of functions known as the spherical harmonics. Note that the energy does depend on through I. The energy is -fold degenerate: the functions with fixed and have the same energy.

Introducing the rotational constant, we write, In the units of reciprocal length the rotational constant is, with c the speed of light. If cgs units are used for , , and , is expressed in cm−1, or wave numbers, which is a unit that is often used for rotational-vibrational spectroscopy. The rotational constant depends on the distance . Often one writes where is the equilibrium value of (the value for which the interaction energy of the atoms in the rotor has a minimum).

A typical rotational absorption spectrum consists of a series of peaks that correspond to transitions between levels with different values of the angular momentum quantum number () such that , due to the selection rules (see below). Consequently, rotational peaks appear at energies with differences corresponding to an integer multiple of .

Selection rules

Rotational transitions of a molecule occur when the molecule absorbs a photon [a particle of a quantized electromagnetic (em) field]. Depending on the energy of the photon (i.e., the wavelength of the em field) this transition may be seen as a sideband of a vibrational and/or electronic transition. Pure rotational transitions, in which the vibronic (= vibrational plus electronic) wave function does not change, occur in the microwave region of the electromagnetic spectrum.

Typically, rotational transitions can only be observed when the angular momentum quantum number changes by . This selection rule arises from a first-order perturbation theory approximation of the time-dependent Schrödinger equation. According to this treatment, rotational transitions can only be observed when one or more components of the dipole operator have a non-vanishing transition moment. If is the direction of the electric field component of the incoming electromagnetic wave, the transition moment is,

A transition occurs if this integral is non-zero. By separating the rotational part of the molecular wavefunction from the vibronic part, one can show that this means that the molecule must have a permanent dipole moment. After integration over the vibronic coordinates the following rotational part of the transition moment remains,

Here is the z component of the permanent dipole moment. The moment is the vibronically averaged component of the dipole operator. Only the component of the permanent dipole along the axis of a heteronuclear molecule is non-vanishing. By the use of the orthogonality of the spherical harmonics it is possible to determine which values of , , , and will result in nonzero values for the dipole transition moment integral. This constraint results in the observed selection rules for the rigid rotor:

Non-rigid linear rotor

The rigid rotor is commonly used to describe the rotational energy of diatomic molecules but it is not a completely accurate description of such molecules. This is because molecular bonds (and therefore the interatomic distance ) are not completely fixed; the bond between the atoms stretches out as the molecule rotates faster (higher values of the rotational quantum number ). This effect can be accounted for by introducing a correction factor known as the centrifugal distortion constant (bars on top of various quantities indicate that these quantities are expressed in cm−1):

where

The non-rigid rotor is an acceptably accurate model for diatomic molecules but is still somewhat imperfect. This is because, although the model does account for bond stretching due to rotation, it ignores any bond stretching due to vibrational energy in the bond (anharmonicity in the potential).

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 R3, so that its energy consists only of rotational kinetic energy (and possibly 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 transitionsone usually classifies molecules (seen as rigid rotors) as follows:

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 the physical convention 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 framethe 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 axes), so that the body-fixed x, y, and z axes coincide with the space-fixed X, Y, and Z axis. Secondly, the body and its frame are rotated actively over a positive angle around the z-axis (by the right-hand 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 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 described here is known as the convention; it can be shown (in the same manner as in this article) 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 be the coordinate vector of an arbitrary point in the body with respect to the body-fixed frame. The elements of are the 'body-fixed coordinates' of . Initially is also the space-fixed coordinate vector of . Upon rotation of the body, 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 (in the physical convention).

Knowledge of the Euler angles as function of time t and the initial coordinates determine the kinematics of the rigid rotor.

Classical kinetic energy

The following text forms 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:

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 on the left hand side contains the components of the angular velocity of the rotor expressed with respect to the body-fixed frame. The angular velocity satisfies equations of motion known as Euler's equations (with zero applied torque, since by assumption the rotor is in field-free space). It can be shown that is not the time derivative of any vector, in contrast to the usual definition of velocity. [2]

The dots over the time-dependent Euler angles on the right hand side indicate time derivatives. Note that a different rotation matrix would result from a different choice of Euler angle convention used.

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 expressed in Euler anglesa 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. With respect to the body-fixed frame it has the components , and can be shown to be related to the angular velocity, This angular momentum is a conserved (time-independent) quantity if viewed from a stationary space-fixed frame. Since the body-fixed frame moves (depends on time) the 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 expressed in terms of the angular momentum 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 inverses of these functions, the position of these trigonometric functions in the classical Hamiltonian is arbitrary. After quantization the commutation does no longer hold 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 indicesin this case over the three Euler angles ):

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. (Note: The corresponding eigenvalue equation gives the Schrödinger equation for the rigid rotor in the form that it was solved for the first time by Kronig and Rabi [3] (for the special case of the symmetric rotor). This 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.)

Nowadays it is common to proceed as follows. 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 () is solved with the degenerate energy equal to .

The symmetric top (= symmetric rotor) is characterized by . It is a prolate (cigar shaped) top if . In the latter case we write the Hamiltonian as and use that Hence The eigenvalue is -fold degenerate, for all eigenfunctions with have the same eigenvalue. The energies with |k| > 0 are -fold degenerate. This exact solution of the Schrödinger equation of the symmetric top was first found in 1927. [3]

The asymmetric top problem () is not soluble analytically.

Direct experimental observation of molecular rotations

For a long time, molecular rotations could not be directly observed experimentally. Only measurement techniques with atomic resolution made it possible to detect the rotation of a single molecule. [4] [5] At low temperatures, the rotations of molecules (or part thereof) can be frozen. This could be directly visualized by scanning tunneling microscopy, i.e., the stabilization could be explained at higher temperatures by the rotational entropy. [5] The direct observation of rotational excitation at single molecule level was achieved recently using inelastic electron tunneling spectroscopy with the scanning tunneling microscope. The rotational excitation of molecular hydrogen and its isotopes were detected. [6] [7]

See also

Related Research Articles

<span class="mw-page-title-main">Pauli matrices</span> Matrices important in quantum mechanics and the study of spin

In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices that are traceless, Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

<span class="mw-page-title-main">Angular velocity</span> Direction and rate of rotation

In physics, angular velocity, also known as angular frequency vector, is a pseudovector representation of how the angular position or orientation of an object changes with time, i.e. how quickly an object rotates around an axis of rotation and how fast the axis itself changes direction.

<span class="mw-page-title-main">Law of sines</span> Property of all triangles on a Euclidean plane

In trigonometry, the law of sines, sine law, sine formula, or sine rule is an equation relating the lengths of the sides of any triangle to the sines of its angles. According to the law, where a, b, and c are the lengths of the sides of a triangle, and α, β, and γ are the opposite angles, while R is the radius of the triangle's circumcircle. When the last part of the equation is not used, the law is sometimes stated using the reciprocals; The law of sines can be used to compute the remaining sides of a triangle when two angles and a side are known—a technique known as triangulation. It can also be used when two sides and one of the non-enclosed angles are known. In some such cases, the triangle is not uniquely determined by this data and the technique gives two possible values for the enclosed angle.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition.

Unit quaternions, known as versors, provide a convenient mathematical notation for representing spatial orientations and rotations of elements in three dimensional space. Specifically, they encode information about an axis-angle rotation about an arbitrary axis. Rotation and orientation quaternions have applications in computer graphics, computer vision, robotics, navigation, molecular dynamics, flight dynamics, orbital mechanics of satellites, and crystallographic texture analysis.

An infinitesimal rotation matrix or differential rotation matrix is a matrix representing an infinitely small rotation.

A nonholonomic system in physics and mathematics is a physical system whose state depends on the path taken in order to achieve it. Such a system is described by a set of parameters subject to differential constraints and non-linear constraints, such that when the system evolves along a path in its parameter space but finally returns to the original set of parameter values at the start of the path, the system itself may not have returned to its original state. Nonholonomic mechanics is autonomous division of Newtonian mechanics.

<span class="mw-page-title-main">Bloch sphere</span> Geometrical representation of the pure state space of a two-level quantum mechanical system

In quantum mechanics and computing, the Bloch sphere is a geometrical representation of the pure state space of a two-level quantum mechanical system (qubit), named after the physicist Felix Bloch.

In linear algebra, a rotation matrix is a transformation matrix that is used to perform a rotation in Euclidean space. For example, using the convention below, the matrix

<span class="mw-page-title-main">Euler's rotation theorem</span> Movement with a fixed point is rotation

In geometry, Euler's rotation theorem states that, in three-dimensional space, any displacement of a rigid body such that a point on the rigid body remains fixed, is equivalent to a single rotation about some axis that runs through the fixed point. It also means that the composition of two rotations is also a rotation. Therefore the set of rotations has a group structure, known as a rotation group.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

<span class="mw-page-title-main">Velocity-addition formula</span> Equation used in relativistic physics

In relativistic physics, a velocity-addition formula is an equation that specifies how to combine the velocities of objects in a way that is consistent with the requirement that no object's speed can exceed the speed of light. Such formulas apply to successive Lorentz transformations, so they also relate different frames. Accompanying velocity addition is a kinematic effect known as Thomas precession, whereby successive non-collinear Lorentz boosts become equivalent to the composition of a rotation of the coordinate system and a boost.

A theoretical motivation for general relativity, including the motivation for the geodesic equation and the Einstein field equation, can be obtained from special relativity by examining the dynamics of particles in circular orbits about the Earth. A key advantage in examining circular orbits is that it is possible to know the solution of the Einstein Field Equation a priori. This provides a means to inform and verify the formalism.

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

The Wigner D-matrix is a unitary matrix in an irreducible representation of the groups SU(2) and SO(3). It was introduced in 1927 by Eugene Wigner, and plays a fundamental role in the quantum mechanical theory of angular momentum. The complex conjugate of the D-matrix is an eigenfunction of the Hamiltonian of spherical and symmetric rigid rotors. The letter D stands for Darstellung, which means "representation" in German.

In mathematics, the Schur orthogonality relations, which were proven by Issai Schur through Schur's lemma, express a central fact about representations of finite groups. They admit a generalization to the case of compact groups in general, and in particular compact Lie groups, such as the rotation group SO(3).

In classical mechanics, holonomic constraints are relations between the position variables that can be expressed in the following form:

In geometry, a plane of rotation is an abstract object used to describe or visualize rotations in space.

A proper reference frame in the theory of relativity is a particular form of accelerated reference frame, that is, a reference frame in which an accelerated observer can be considered as being at rest. It can describe phenomena in curved spacetime, as well as in "flat" Minkowski spacetime in which the spacetime curvature caused by the energy–momentum tensor can be disregarded. Since this article considers only flat spacetime—and uses the definition that special relativity is the theory of flat spacetime while general relativity is a theory of gravitation in terms of curved spacetime—it is consequently concerned with accelerated frames in special relativity.

In the field of mathematics known as complex analysis, the indicator function of an entire function indicates the rate of growth of the function in different directions.

References

  1. 1 2 Podolsky, B. (1928). "Quantum-Mechanically Correct Form of Hamiltonian Function for Conservative Systems". Phys. Rev. 32 (5): 812. Bibcode:1928PhRv...32..812P. doi:10.1103/PhysRev.32.812.
  2. Goldstein, Herbert; Poole, Charles P.; Safko, John L. (2002). Classical mechanics (3rd ed.). San Francisco: Addison Wesley. Chapter 4.9. ISBN   0-201-65702-3. OCLC   47056311.
  3. 1 2 R. de L. Kronig and I. I. Rabi (1927). "The Symmetrical Top in the Undulatory Mechanics". Phys. Rev. 29 (2): 262–269. Bibcode:1927PhRv...29..262K. doi:10.1103/PhysRev.29.262. S2CID   4000903.
  4. J. K. Gimzewski; C. Joachim; R. R. Schlittler; V. Langlais; H. Tang; I. Johannsen (1998), "Rotation of a Single Molecule Within a Supramolecular Bearing", Science (in German), vol. 281, no. 5376, pp. 531–533, Bibcode:1998Sci...281..531G, doi:10.1126/science.281.5376.531, PMID   9677189
  5. 1 2 Thomas Waldmann; Jens Klein; Harry E. Hoster; R. Jürgen Behm (2012), "Stabilization of Large Adsorbates by Rotational Entropy: A Time-Resolved Variable-Temperature STM Study", ChemPhysChem (in German), vol. 14, no. 1, pp. 162–169, doi:10.1002/cphc.201200531, PMID   23047526, S2CID   36848079
  6. Li, Shaowei; Yu, Arthur; Toledo, Freddy; Han, Zhumin; Wang, Hui; He, H. Y.; Wu, Ruqian; Ho, W. (2013-10-02). "Rotational and Vibrational Excitations of a Hydrogen Molecule Trapped within a Nanocavity of Tunable Dimension". Physical Review Letters. 111 (14): 146102. doi:10.1103/PhysRevLett.111.146102. ISSN   0031-9007.
  7. Natterer, Fabian Donat; Patthey, François; Brune, Harald (2013-10-24). "Distinction of Nuclear Spin States with the Scanning Tunneling Microscope". Physical Review Letters. 111 (17): 175303. arXiv: 1307.7046 . doi:10.1103/PhysRevLett.111.175303. ISSN   0031-9007.

General references