Angular momentum (classical)

From Citizendium
Revision as of 23:37, 30 December 2007 by imported>David E. Volk (subpages)
Jump to navigation Jump to search
This article is developed but not approved.
Main Article
Discussion
Related Articles  [?]
Bibliography  [?]
External Links  [?]
Citable Version  [?]
 
This editable, developed Main Article is subject to a disclaimer.

In physics, angular momentum is a kinematic property of a system consisting of one or more point masses. Its importance derives from the fact that it is conserved in several physical circumstances, notably for objects in field-free space and for central-symmetric interactions. (Recall that a conserved property of a system is one that does not change over time.)

Before the advent of quantum mechanics, angular momentum did not play an important role in physics, as is witnessed by the 450 page book of Whittaker[1] in which three and a half page is devoted to it and the index contains only one reference to angular momentum. Cajori, in his well-known history of physics[2] gives exactly half a line to the conservation of angular momentum. Since angular momentum is of great importance in quantum mechanics, post quantum mechanics texts on classical mechanics pay much more attention to it. For instance Goldstein[3] has about 25 entries in its index and about a fifth of his book devoted to it. See this article for the definition and role of angular momentum in quantum mechanics.

Point mass m with velocity v has angular momentum Lm r × v with respect to O. L is a vector pointing towards the reader, perpendicular to the plane of drawing, which is spanned by r and v.

Definition

The angular momentum of a single point mass m is defined with respect to a point O. Denote the vector from O to m by r (see the figure). Let the mass have velocity v, then the angular momentum L of the point mass is defined as the cross product,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L} \equiv m \mathbf{r}\times \mathbf{v} = \mathbf{r}\times \mathbf{p}. }

It follows from the definition of cross product that the vector L is perpendicular to the plane of the figure and points towards the reader. Note that the definition pm v for the linear momentum is introduced here.

The definition of the angular momentum of a system of n point masses is the following simple generalization:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L}_\mathrm{sys} \equiv \mathbf{r}_1\times \mathbf{p}_1 + \mathbf{r}_2\times \mathbf{p}_2 + \cdots + \mathbf{r}_n\times \mathbf{p}_n = \sum_{i=1}^n \mathbf{r}_i\times \mathbf{p}_i. }

Conservation of angular momentum

For simplicity we consider first the case of one point mass, generalization to n point masses is straightforward. The following important result will be proved: a particle that moves in a centrally symmetric potential has a conserved angular momentum.

Using Newton's dot (fluxion) notation for time derivatives, we find

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{\dot{L}} = m\mathbf{\dot{r}} \times \mathbf{v} + m \mathbf{r} \times \mathbf{\dot{v}}. }

Use

where we invoked Newton's second law Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle m \mathbf{\dot{v}} \equiv m \mathbf{a} = \mathbf{F}. } Hence the time derivative of the angular momentum is equal to the torque N,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{\dot{L}} = \mathbf{r} \times \mathbf{F} \equiv \mathbf{N}. }

It follows that the angular momentum is conserved—i.e., its time derivative is zero—if the applied torque N is zero, which means that either F is zero, or r and F are parallel, because in that case the cross product vanishes.

The former case occurs when the mass moves uniformly in homogeneous space (no external forces). Then also the linear momentum of the point mass is conserved, and therefore one rarely considers the angular momentum in the case of uniform rectilinear motion.

The second case is more interesting. Let us assume that F is conservative, i.e., its curl vanishes, then F is the gradient of a potential V(x,y,z),

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{\nabla}\times \mathbf{F} = 0 \Longrightarrow \mathbf{F} = - \mathbf{\nabla} V(x,y,z), }

so that for a conservative force the time derivative of L is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{\dot{L}} = - \left( \mathbf{r} \times \mathbf{\nabla}\right) V . }

The quantity in brackets is known in group theory as a Lie derivative of the full rotation group SO(3). In quantum mechanics it is a well-known operator, namely the orbital angular momentum operator (except for the factor Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle i\hbar } ).

From the general properties of Lie derivatives—or equivalently orbital angular momenta—it follows that the time derivative of L vanishes whenever the potential V(x,y,z) is rotationally invariant (isotropic, central-symmetric). This means that the function V(x,y,z) depends on r only, when expressed in spherical polar coordinates (r, θ, φ). By expressing the gradient in spherical polar coordinates and writing the cross product as a determinant, we can prove the vanishing explicitly, using that the derivatives of V with respect to θ and φ are zero,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}\times \mathbf{\nabla}V = \begin{vmatrix} \mathbf{e}_r &\mathbf{e}_\theta & \mathbf{e}_\phi \\ r & 0 & 0 \\ \frac{\partial V}{\partial r} & \frac{1}{r}\frac{\partial V}{\partial \theta} & \frac{1}{r\sin \theta} \frac{\partial V}{\partial \phi} \\ \end{vmatrix} = -\mathbf{e}_\theta\, \left(\frac{1}{\sin \theta} \frac{\partial V}{\partial \phi}\right) + \mathbf{e}_\phi\, \left( \frac{\partial V}{\partial \theta}\right) = 0. }


A well-known example is the motion of the earth around the sun. To a very good approximation the gravitation is central-symmetric (rotationally invariant) and hence the earth's orbit is characterized by a conserved angular momentum. Kepler's second law is a direct consequence of this conservation of angular momentum.

Extension to n point masses

When in a system of n point masses the particles do not interact mutually, the extension is simply:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{\dot{L}}_i = \mathbf{r}_i \times \mathbf{F}_i \Longrightarrow \mathbf{\dot{L}}_\mathrm{sys} = \sum_{i=1}^n \mathbf{r}_i \times \mathbf{F}_i =\mathbf{N}^\mathrm{ext}_\mathrm{sys}, }

where the forces F i are external. If the total external torque is zero, the total angular momentum of the system is conserved.

When the particles do interact, i.e., there are non-zero internal forces, then the result still holds, provided the internal forces are central-symmetric and additive. No internal torque occurs in that case. In order to show this, we recall that Newton's third law (action = −reaction) holds for central-symmetric forces. Indicating internal forces by a double suffix, we write Newton's law as follows,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{F}_{ij} = - \mathbf{F}_{ji}, \qquad i \ne j. }

where the left hand side is the force on particle i executed by particle j and the right hand side is i acting on j. Both forces act along the vector pointing form particle i to j,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_{ij} = \mathbf{r}_j - \mathbf{r}_i, }

so that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_{ij}\times \mathbf{F}_{ij} = 0, }

where

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_{ij}\times \mathbf{F}_{ij} = (\mathbf{r}_j - \mathbf{r}_i)\times \mathbf{F}_{ij} = -(\mathbf{r}_j \times \mathbf{F}_{ji} + \mathbf{r}_i\times \mathbf{F}_{ij}) . }

The total internal torque on particle i is,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{N}^\mathrm{int}_i = \mathbf{r}_{i}\times \Big(\sum^n_{{j=1 \atop j\ne i}} \mathbf{F}_{ij}\Big). }

Under assumption of additivity, the total internal torque can be written as

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{N}^\mathrm{int}_\mathrm{sys} = \sum_{i=1}^n \mathbf{N}^\mathrm{int}_i = \sum_{i=1}^n \sum^n_{j \ne i} \mathbf{r}_{i}\times \mathbf{F}_{ij} = \sum_{i=1}^n \sum_{j>i}^n \big( \mathbf{r}_{i}\times \mathbf{F}_{ij} + \mathbf{r}_{j}\times \mathbf{F}_{ji}\big) = - \sum_{i=1}^n \sum_{j>i}^n \mathbf{r}_{ij}\times \mathbf{F}_{ij} = 0, }

and

where F i (characterized by a single suffix) is external.

Separation of center of mass

The angular momentum of a system of point masses with respect to an arbitrary point O can be separated into two terms. The first term is the angular momentum with respect to O of a fictitious particle with total mass M positioned at the center of mass C of the system. The other term is the sum of the angular momenta of the particles with respect to C.

The center of mass C is determined by

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \overrightarrow{OC} = \mathbf{R} \equiv \frac{1}{M} \sum_{i=1}^n m_i \mathbf{r}_i \quad \hbox{with}\quad M \equiv \sum_{i=1}^n m_i. }

Introduce s i for particle i at point P i

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_i \equiv \overrightarrow{OP}_i = \overrightarrow{OC} + \overrightarrow{CP}_i \equiv \mathbf{R} + \mathbf{s}_i. }

Note that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sum_{i=1}^n m_i \mathbf{s}_i = \sum_{i=1}^n m_i \mathbf{r}_i - \mathbf{R} \sum_{i=1}^n m_i = \sum_{i=1}^n m_i \mathbf{r}_i -\left( \frac{1}{M} \sum_{i=1}^n m_i \mathbf{r}_i\right) M = 0, }

and hence also

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sum_{i=1}^n m_i \mathbf{\dot{s}}_i = 0. }

The two cross terms in the following equation vanish and the two diagonal terms are non-vanishing,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L}_\mathrm{sys} = \sum_{i=1}^n m_i \mathbf{r}_i \times \mathbf{v}_i = \sum_{i=1}^n m_i \Big(\mathbf{R}+\mathbf{s}_i \Big) \times \Big(\mathbf{\dot{R}} + \mathbf{\dot{s}}_i \Big) = \sum_{i=1}^n m_i \,\; \mathbf{R} \times \mathbf{\dot{R}} \,+\, \sum_{i=1}^n m_i \mathbf{s}_i \times \mathbf{\dot{s}}_i . }

The two remaining terms are:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle M \mathbf{R} \times \mathbf{\dot{R}} \quad\hbox{and}\quad \sum_{i=1}^n m_i \mathbf{s}_i \times \mathbf{\dot{s}}_i. }

If the total system moves in homogeneous (field free) space, the first term is conserved and is usually not of much interest. The second term, the sum of the angular momenta of the n particles with respect to the center of mass, is usually designated the total angular momentum of the system,

where the position vectors s i are defined with respect to a system of axes with origin in the center of mass C of the system.

Circular motion around a fixed axis

A system often encountered in physics is a particle (point mass) that moves along a circle in a fixed plane. Let us choose a system of axes such that the plane of motion is the x-y plane, with the z-axis perpendicular to the plane.

The coordinates of the particle and its time derivative are

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}(t) = \begin{pmatrix} r \cos\alpha(t) \\ r \sin\alpha(t) \\ 0 \\ \end{pmatrix} \quad \hbox{and}\quad \mathbf{v}(t)\equiv \mathbf{\dot r}(t) =\dot{\alpha}(t) \begin{pmatrix} - r\sin\alpha(t) \\ r\cos\alpha(t) \\ 0 \\ \end{pmatrix}, }

where we used that the radius r of the circle is time-independent. The angular momentum is given by

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L} = m \begin{vmatrix} \mathbf{e}_x & \mathbf{e}_y & \mathbf{e}_z \\ r\cos\alpha & r\sin\alpha & 0 \\ -r \dot{\alpha} \sin\alpha & r\dot{\alpha}\cos\alpha & 0 \\ \end{vmatrix} = m r^2 \, \omega(t) \mathbf{e}_z \quad\hbox{with}\quad \omega(t)\equiv \dot\alpha(t) . }
Point mass m moves counterclock-wise on a circle. Instantaneous velocity v is tangent to the circle. Angular velocity and angular momentum point toward the reader.

The quantity mr 2 is the inertia moment Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle I\,} of the present simple system (particle on a circle) and ω is the angular velocity. The angular velocity, as defined here, is a scalar. Often the product ωez is defined as angular velocity ω, which is a vector.

In summary, for a particle moving on a circle around a fixed axis, the angular momentum is proportional to the angular velocity with the inertia moment as proportionality factor,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L} = I \boldsymbol{\omega}. }

It is of some interest to point out that this simple relation becomes a vector-tensor relation for more general motions of more general systems (rigid rotors) for which the rotational axis is not fixed in space: L and ω remain vectors, but the inertia moment generalizes to a second rank, symmetric, tensor.

Angular momentum of a rigid body; inertia tensor

The angular momentum of a finite size rigid body (rigid rotor) will be introduced. The angular momentum is expressed with respect to a laboratory system of axes with origin at the center of mass of the rotor. The following formula will be derived:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L} = \mathbb{I}(t) \boldsymbol{\omega}, }

where the inertia tensor is given by

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbb{I}(t) \equiv \sum_{i=1}^n m_i \begin{pmatrix} y_i^2+z_i^2 & -x_iy_i & -x_iz_i \\ -y_ix_i & x_i^2+z_i^2 & -y_iz_i \\ -z_ix_i & -z_iy_i & x_i^2+y_i^2 \\ \end{pmatrix}, }

and ω is the angular velocity of the rotor. The coordinates x i, y i, and z i of the particles depend on time. In order not cloud the notation this dependence is not shown explicitly. Note that, since the inertia tensor is in general not (a constant times) the identity matrix, the two vectors, L and ω, are in general not parallel, in contrast to the case with one fixed rotation axis—motion on a circle.

Angular velocity of a rigid body

We consider a system of n point masses (particles) with fixed—time-independent—relative positions. The position vectors r i(t) are all with respect to a laboratory-fixed frame with origin in the center of mass of the system. The lengths of the n vectors and the angles between them are all time-independent, because the system is rigid. We will prove the following expression for the time derivatives of the position vectors,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\mathbf{r}}_i(t) = \boldsymbol{\omega}(t) \times \mathbf{r}_i(t), \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad(1) }

where the vector ω(t) is a 3-vector which is the same for all n particles, it is the angular velocity of the rigid body.

The time dependence of the position vectors of a rigid body has the following matrix-vector form

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_i(t) = \mathbb{R}(t)\, \mathbf{r}_i(0),\quad (i=1,2,\ldots,n) \qquad\hbox{with}\quad \mathbb{R}(t)^T \mathbb{R}(t) = \mathbb{R}(t) \mathbb{R}(t)^T =\mathbb{I}. }

The vectors r i(0) determine the shape of the rigid body at time zero. This shape is invariant since all inner products are time-independent. Indeed,[4]

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_i(t) \cdot \mathbf{r}_j(t) = \mathbf{r}_i(0)^T \mathbb{R}(t)^T \mathbb{R}(t) \mathbf{r}_j(0) = \mathbf{r}_i(0) \cdot \mathbf{r}_j(0). }

Differentiation with respect to time gives

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{0} = \dot\mathbb{I} = \dot{\mathbb{R}}(t) \mathbb{R}(t)^T + \mathbb{R}(t) \dot{\mathbb{R}}(t)^T, }

and hence the matrix product Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \dot{\mathbb{R}}(t) \mathbb{R}(t)^T } is equal to an antisymmetric matrix. A 3×3 antisymmetric matrix is parametrized by three real numbers and has the following form,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{\mathbb{R}}(t)\mathbb{R}(t)^T = -\Big( \dot{\mathbb{R}}(t) \mathbb{R}(t)^T\Big)^T = \begin{pmatrix} 0 & - \omega_z & \omega_y \\ \omega_z & 0 & - \omega_x \\ -\omega_y & \omega_x & 0 \\ \end{pmatrix}. }

Differentiation of a position vector with respect to time, followed by insertion of the identity matrix gives,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \dot{\mathbf{r}}_i(t) &= \dot{\mathbb{R}}(t)\, \mathbf{r}_i(0) = \dot{\mathbb{R}}(t)\, \mathbb{R}(t)^T \mathbb{R}(t)\, \mathbf{r}_i(0) = \dot{\mathbb{R}}(t) \mathbb{R}(t)^T \, \mathbf{r}_i(t)\\ &= \begin{pmatrix} 0 & - \omega_z & \omega_y \\ \omega_z & 0 & - \omega_x \\ -\omega_y & \omega_x & 0 \\ \end{pmatrix} \mathbf{r}_i(t) = \boldsymbol{\omega} \times \mathbf{r}_i(t) , \end{align} }

where the vector ω is composed of the matrix elements of the antisymmetric matrix,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \boldsymbol{\omega} = \left(\omega_x, \, \omega_y,, \, \omega_z \right) }

The fact that the matrix-vector multiplication may be written as a cross product follows most easily by direct verification. This proves the desired result, equation (1).

Inertia moment

Let us first consider one arbitrary particle, for the sake of clarity we drop the particle index. The angular momentum consists of terms of the type

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}(t)\times \mathbf{\dot{r}}(t) = \mathbf{r}(t)\times \big(\boldsymbol{\omega}\times \mathbf{r}(t)\big) }

where we used Eq. (1). From vector algebra we know that (we drop temporarily the time-dependence)

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}\times \big(\boldsymbol{\omega}\times \mathbf{r}\big) = r^2 \boldsymbol{\omega} - \mathbf{r} (\boldsymbol{\omega}\cdot \mathbf{r}), }

which written out in components is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} & \begin{pmatrix} r^2\omega_x - x(x\omega_x + y\omega_y + z\omega_z) \\ r^2\omega_y - y(x\omega_x + y\omega_y + z\omega_z) \\ r^2\omega_z -z(x\omega_x + y\omega_y + z\omega_z) \\ \end{pmatrix} = \begin{pmatrix} r^2 - x^2 & -xy & -xz \\ -yx & r^2 -y^2 & -yz \\ -zx & -zy & r^2-z^2 \\ \end{pmatrix} \begin{pmatrix} \omega_x \\ \omega_y \\ \omega_z \\ \end{pmatrix} \\ &= \begin{pmatrix} y^2+z^2 & -xy & -xz \\ -yx & x^2+z^2 & -yz \\ -zx & -zy & x^2+y^2 \\ \end{pmatrix} \begin{pmatrix} \omega_x \\ \omega_y \\ \omega_z \\ \end{pmatrix} \end{align} }

We define the inertia moment of the rigid body by

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbb{I}(t) \equiv \sum_{i=1}^n m_i \begin{pmatrix} y_i^2+z_i^2 & -x_iy_i & -x_iz_i \\ -y_ix_i & x_i^2+z_i^2 & -y_iz_i \\ -z_ix_i & -z_iy_i & x_i^2+y_i^2 \\ \end{pmatrix}. }

The coordinates of the particle are time-dependent and expressed with respect to a laboratory frame. Insertion into the expression of the angular momentum gives the final result

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L} = \sum_{i=1}^n m_i \mathbf{r}_i(t)\times \mathbf{\dot{r}}_i(t) = \mathbb{I}(t) \boldsymbol{\omega} . }

Body-fixed angular momentum

In quantum mechanical applications the angular momentum—a vector—is often expressed with respect to a system of axes that is attached to the body and rotates with it. We will find an expression for it. Consider to that end first the equation,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}(t) \times \big(\boldsymbol{\omega} \times \mathbf{r}(t)\big) \equiv \mathbb{R}(t)\mathbf{r}(0) \times \big(\boldsymbol{\omega} \times \mathbb{R}(t)\mathbf{r}(0) \big) = \mathbb{R}(t)\left[ \mathbf{r}(0) \times\big(\boldsymbol{\omega}' \times \mathbf{r}(0) \big) \right], }

where we used twice that for arbitrary vectors a and b and orthogonal matrix with unit determinant,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbb{R} \mathbf{a} \times \mathbb{R} \mathbf{b} = \mathbb{R}\left[ \mathbf{a} \times \mathbf{b} \right], }

and where we introduced for temporary use the short-hand notation,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \boldsymbol{\omega}' \equiv \mathbb{R}(t)^T \boldsymbol{\omega}. }

In the same way as above, the expression in curly brackets can be worked out:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L} =\sum_{i=1}^n m_i \mathbf{r}_i(t)\times \big(\boldsymbol{\omega}\times \mathbf{r}_i(t)\big) = \mathbb{R}(t)\Big\{\sum_{i=1}^n m_i \mathbf{r}_i(0) \times\big( \boldsymbol{\omega}' \times \mathbf{r}_i(0) \big) \Big\} = \mathbb{R}(t) \, \mathbb{I}(0)\, \boldsymbol{\omega}'. }

From this and the earlier result follows parenthetically that the inertia tensor is indeed a tensor,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbb{I}(t)= \mathbb{R}(t)\, \mathbb{I}(0)\, \mathbb{R}(t)^T, }

because this rotational relation is often taken as defining a tensor.

Further,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L} = \mathbb{R}(t)\, \mathbb{I}(0)\, \mathbb{R}(t)^T \,\boldsymbol{\omega} \Longrightarrow \mathbb{R}(t)^T\, \mathbf{L} = \mathbb{I}(0)\, \mathbb{R}(t)^T \,\boldsymbol{\omega} }

The vectors L and ω in this relation are expressed with respect to a laboratory-fixed frame. If we now make the assumption that this frame coincides at time zero with the body-fixed frame, it follows that Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \mathbb{I}(0)} is the inertia moment in the body-fixed frame. Define further,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L}^\mathrm{BF} \equiv \mathbb{R}(t)^T \mathbf{L} \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad(2) }

and

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \boldsymbol{\omega}^\mathrm{BF} \equiv \mathbb{R}(t)^T \boldsymbol{\omega} }

then we see that in the body-fixed frame angular momentum depends again linearly on angular velocity through a tensor-vector expression,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{L}^\mathrm{BF} = \mathbb{I}(0) \boldsymbol{\omega}^\mathrm{BF}, }

where the inertia tensor is constant.

We have shown above that L, which contains the three components of the physical vector angular momentum Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{L}} with respect to the stationary laboratory frame, is independent of time when there is no external torque. Hence the physical vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{L}} , a quantity that exists irrespective of choice of frame, is itself time-independent in absence of torque. Now, LBF contains the components of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{L}} , but re-expressed with respect to a moving frame. Hence, the components of LBF are time-dependent. Indeed, its time-dependence is given by Eq. (2).

References

  1. E. T. Whittaker, A Treatise on the Analytical Dynamics of Particles and Rigid bodies, Cambridge University Press, 1st edition (1904)
  2. F. Cajori, History of Physics, Macmillan, New York (1929)
  3. H. Goldstein, Classical Mechanics, Addison-Wesley, Reading, Massachusetts, 1st edition (1950)
  4. We use here two equivalent notations, valid for any two vectors of the same dimension,
    Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r} \cdot \mathbf{s} \equiv \mathbf{r}^T \mathbf{s} , }
    where the notation on the right hand side is based on the fact that column and row vectors are rectangular matrices. Further we used the matrix transposition rule,
    Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \big( \mathbb{A}\mathbb{B} \big)^T = \mathbb{B}^T \mathbb{A}^T. }