Angular momentum (classical): Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Paul Wormer
mNo edit summary
 
(10 intermediate revisions by 4 users not shown)
Line 1: Line 1:
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. (Recall that a conserved property of a system is one that does not change over time.)
{{subpages}}
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.)


[[Image:Angular momentum particle.png|right|200px|thumb|Point mass ''m'' with velocity '''v''' has angular momentum '''L''' ≡ ''m'' '''r''' × '''v''' with respect to ''O''. '''L''' is a vector pointing towards the reader]]
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<ref>E. T. Whittaker, ''A Treatise on the Analytical Dynamics of Particles and Rigid bodies'', Cambridge University Press, 1st edition (1904)</ref> 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<ref>F. Cajori, ''History of Physics'', Macmillan, New York (1929)</ref> 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<ref>H. Goldstein, ''Classical Mechanics'', Addison-Wesley, Reading, Massachusetts, 1st edition (1950)</ref> has about 25 entries in its index and about a fifth of his book devoted to it.  See [[angular momentum (quantum)|this article]] for the definition and role of angular momentum in quantum mechanics.
{{Image|Angular momentum particle.png|right|200px|Point mass ''m'' with velocity '''v''' has angular momentum '''L''' &equiv; ''m'' '''r''' &times; '''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==
==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 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
Line 16: Line 20:


==Conservation of angular momentum==
==Conservation of angular momentum==
For simplicity we consider 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.''
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 central-symmetric potential has a conserved angular momentum.''


Using Newton's dot (fluxion) notation for time derivatives, we find
Using Newton's dot (fluxion) notation for time derivatives, we find
Line 38: Line 42:
It follows that the angular momentum is conserved&mdash;i.e., its time derivative is zero&mdash;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.  
It follows that the angular momentum is conserved&mdash;i.e., its time derivative is zero&mdash;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 former case ('''F''' = '''0''') 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 (physics)|conservative]], i.e., its [[curl]] vanishes, then '''F''' is the [[gradient]] '''&nabla;''' of a potential ''V''(''x'',''y'',''z''),
The second case ('''r''' parallel to '''F''') is more interesting. In order to see what this means for the potential, we must assume that '''F''' is [[conservative (physics)|conservative]], i.e., its [[curl]] vanishes, otherwise the potential cannot be defined. If  '''F''' is conservative,
then '''F''' is minus the [[gradient]] '''&nabla;''' of a potential ''V''(''x'',''y'',''z''),
:<math>
:<math>
\mathbf{\nabla}\times \mathbf{F} = 0 \Longrightarrow \mathbf{F} = - \mathbf{\nabla} V(x,y,z),
\mathbf{\nabla}\times \mathbf{F} = 0 \Longrightarrow \mathbf{F} = - \mathbf{\nabla} V(x,y,z),
</math>
</math>
so that for a conservative force the time derivative of '''L''' is
so that for a conservative force the torque is
:<math>
:<math>
\mathbf{\dot{L}} = - \left( \mathbf{r} \times \mathbf{\nabla}\right) V .
\mathbf{N} = - \left( \mathbf{r} \times \mathbf{\nabla}\right) V .
</math>
</math>
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 <math>\scriptstyle i\hbar </math>).
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 <math>\scriptstyle i\hbar </math>).


From the general properties of Lie derivatives&mdash;or equivalently orbital angular momenta&mdash;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'', &theta;, &phi;). By expressing the gradient '''&nabla;''' in spherical polar coordinates and writing the cross product as a determinant, we can prove this explicitly, using that the derivatives of ''V'' with respect to &theta; and &phi; are zero,
From the general properties of Lie derivatives&mdash;or equivalently orbital angular momenta&mdash;it follows that the torque vanishes whenever the potential ''V''(''x'',''y'',''z'') is rotationally invariant (also known as an isotropic or central-symmetric potential). This means that the function ''V''(''x'',''y'',''z'') depends on ''r'' only, when expressed in [[spherical polar coordinates]] (''r'', &theta;, &phi;).  
 
We can prove the vanishing of the torque more explicitly by expressing the gradient '''&nabla;''' in spherical polar coordinates and writing the cross product as a determinant, using that the derivatives of ''V'' with respect to &theta; and &phi; are zero,
:<math>
:<math>
\mathbf{r}\times \mathbf{\nabla}V =  
\mathbf{r}\times \mathbf{\nabla}V =  
Line 58: Line 65:
\frac{\partial V}{\partial r} & \frac{1}{r}\frac{\partial V}{\partial \theta} & \frac{1}{r\sin \theta} \frac{\partial V}{\partial \phi} \\
\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}
\end{vmatrix}
= -\mathbf{e}_\theta\, \left(\frac{1}{\sin \theta} \frac{\partial V}{\partial \phi}\right)
=\begin{vmatrix}
+ \mathbf{e}_\phi\\left( \frac{\partial V}{\partial \theta}\right) = 0.
\mathbf{e}_r &\mathbf{e}_\theta & \mathbf{e}_\phi \\
r            & 0                &    0          \\
\frac{\partial V}{\partial r}  &    0 &0 \\
\end{vmatrix}
= 0.
</math>
</math>
Note that the third member in this equation shows that '''r''' (the second row) and &minus;'''F''' (the third row) are indeed parallel to each other in the case of an isotropic potential. Both vectors possess only a radial coordinate; '''r''' by definition and '''F''' because the gradient of ''V'' has only one non-vanishing component, the radial component.


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 laws|Kepler's]] second law is a direct consequence of this conservation of angular momentum.
A well-known example of a system with conserved angular momentum is the earth orbiting 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, i.e., the time derivative of the angular momentum (torque) is zero. [[Kepler's laws|Kepler's]] second law is a direct consequence of this conservation of angular momentum.
===Generalization to ''n'' point masses===
===Extension to ''n'' point masses===
When in a system of ''n'' point masses the particles do ''not'' interact mutually, the generalization is simply
When in a system of ''n'' point masses the particles do ''not'' interact mutually, the extension is simply:
:<math>
:<math>
\mathbf{\dot{L}}_i =  \mathbf{r}_i \times \mathbf{F}_i \Longrightarrow
\mathbf{\dot{L}}_i =  \mathbf{r}_i \times \mathbf{F}_i \Longrightarrow
Line 107: Line 119:


==Separation of center of mass==
==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 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 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
The center of mass ''C'' is determined by
Line 129: Line 141:
\sum_{i=1}^n m_i \mathbf{\dot{s}}_i = 0.
\sum_{i=1}^n m_i \mathbf{\dot{s}}_i = 0.
</math>
</math>
Now,
The two cross terms in the following equation vanish and the two  diagonal terms are non-vanishing,
:<math>
:<math>
\mathbf{L}_\mathrm{sys} = \sum_{i=1}^n m_i \mathbf{r}_i \times \mathbf{v}_i =   
\mathbf{L}_\mathrm{sys} = \sum_{i=1}^n m_i \mathbf{r}_i \times \mathbf{v}_i =   
Line 135: Line 147:
\Big(\mathbf{\dot{R}} + \mathbf{\dot{s}}_i \Big) =
\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{R} \times \mathbf{\dot{R}}  
\,+\, \sum_{i=1}^n m_i \mathbf{s}_i \times  \mathbf{\dot{s}}_i
\,+\, \sum_{i=1}^n m_i \mathbf{s}_i \times  \mathbf{\dot{s}}_i .
</math>  
</math>  
where the two cross terms vanished and the two  diagonal terms are non-vanishing. The two remaining terms are:
The two remaining terms are:
:<math>
:<math>
  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.
  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.
Line 150: Line 162:
==Circular motion around a fixed axis==
==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.  
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.  
[[Image:Circular motion.png|right|200px|thumb|Point mass ''m'' with instantaneous velocity '''v''' moves counterclock-wise on circle.]]
 
The coordinates of the particle and its time derivative are
The coordinates of the particle and its time derivative are
:<math>
:<math>
\mathbf{r} =  
\mathbf{r}(t) =  
\begin{pmatrix}
\begin{pmatrix}
r \cos\alpha \\
r \cos\alpha(t) \\
r \sin\alpha \\
r \sin\alpha(t) \\
0 \\
0 \\
\end{pmatrix}  
\end{pmatrix}  
\quad \hbox{and}\quad
\quad \hbox{and}\quad
\mathbf{v}\equiv \mathbf{\dot r} =\dot{\alpha}  
\mathbf{v}(t)\equiv \mathbf{\dot r}(t) =\dot{\alpha}(t)
\begin{pmatrix}
\begin{pmatrix}
- r\sin\alpha \\
- r\sin\alpha(t) \\
   r\cos\alpha \\
   r\cos\alpha(t) \\
0 \\
0 \\
\end{pmatrix},
\end{pmatrix},
Line 174: Line 186:
r\cos\alpha  & r\sin\alpha  &    0  \\
r\cos\alpha  & r\sin\alpha  &    0  \\
-r \dot{\alpha} \sin\alpha & r\dot{\alpha}\cos\alpha & 0 \\
-r \dot{\alpha} \sin\alpha & r\dot{\alpha}\cos\alpha & 0 \\
\end{vmatrix} = m r^2 \, \omega \mathbf{e}_z \quad\hbox{with}\quad \omega\equiv \dot\alpha .
\end{vmatrix} = m r^2 \, \omega(t) \mathbf{e}_z \quad\hbox{with}\quad \omega(t)\equiv \dot\alpha(t) .
</math>
</math>
 
{{Image|Circular motion.png|right|200px|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''<sup> 2</sup> is the [[inertia moment]] ''I'' of the present simple system (particle on a circle) and &omega; is the [[angular velocity]].  The angular velocity, as defined here,  is a scalar. Often the product &omega;'''e'''<sub>''z''</sub> is defined as angular velocity '''&omega;''', which is a vector.  
The quantity ''mr''<sup> 2</sup> is the [[inertia moment]] <math>I\,</math> of the present simple system (particle on a circle) and &omega; is the [[angular velocity]].  The angular velocity, as defined here,  is a scalar. Often the product &omega;'''e'''<sub>''z''</sub> is defined as angular velocity '''&omega;''', 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,
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,
Line 187: Line 199:


==Angular momentum of a rigid body; inertia tensor==
==Angular momentum of a rigid body; inertia tensor==
The following formula will be derived  for  the angular momentum of a [[rigid rotor]] (rigid body). The angular momentum is expressed with respect to a laboratory system of axes with origin at the center of mass of the rotor.  
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:
:<math>
:<math>
\mathbf{L} = \mathbb{I}(t) \boldsymbol{\omega},
\mathbf{L} = \mathbb{I}(t) \boldsymbol{\omega},
Line 200: Line 212:
\end{pmatrix},
\end{pmatrix},
</math>
</math>
and '''&omega;''' is the angular velocity of the rotor. The coordinates ''x''<sub>'' i''</sub>, ''y''<sub>'' i''</sub>, and ''z''<sub>'' i''</sub> of the particles depend on time.
and '''&omega;''' is the angular velocity of the rotor. The coordinates ''x''<sub>'' i''</sub>, ''y''<sub>'' i''</sub>, and ''z''<sub>'' i''</sub> 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 '''&omega;''', are in general not parallel, in contrast to the case with one fixed rotation axis&mdash;motion on a circle.
 
===Angular velocity of  a rigid body===
===Angular velocity of  a rigid body===
We consider a system of ''n'' point masses (particles) with fixed&mdash;time-independent&mdash;relative positions. The position vectors '''r'''<sub>'' i''</sub>(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,
We consider a system of ''n'' point masses (particles) with fixed&mdash;time-independent&mdash;relative positions. The position vectors '''r'''<sub>'' i''</sub>(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,
:<math>
:<math>
\dot{\mathbf{r}}_i(t) = \boldsymbol{\omega} \times  \mathbf{r}_i(t),
\dot{\mathbf{r}}_i(t) = \boldsymbol{\omega}(t) \times  \mathbf{r}_i(t),
\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad(1)
\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad(1)
</math>
</math>
where the vector '''&omega;''' is a 3-vector which is the same for all ''n'' particles, it is the ''angular velocity'' of the rigid body.
where the vector '''&omega;'''(''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
The time dependence of the position vectors of a rigid body has the following [[matrix]]-vector form
Line 214: Line 228:
\mathbb{R}(t)^T \mathbb{R}(t) = \mathbb{R}(t) \mathbb{R}(t)^T =\mathbb{I}.
\mathbb{R}(t)^T \mathbb{R}(t) = \mathbb{R}(t) \mathbb{R}(t)^T =\mathbb{I}.
</math>
</math>
The vectors  '''r'''<sub>'' i''</sub>(0) determine the shape of the rigid body at time zero. This shape is invariant since all inner products are time-independent. Indeed,
The vectors  '''r'''<sub>'' i''</sub>(0) determine the shape of the rigid body at time zero. This shape is invariant since all inner products are time-independent. Indeed,<ref>
:<math>
We use here two equivalent notations, valid for any two vectors of the same dimension,
\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).
</math>
We used here the equivalent notations, valid for any two vectors of the same dimension,
:<math>
:<math>
\mathbf{r} \cdot \mathbf{s} \equiv \mathbf{r}^T \mathbf{s} ,
\mathbf{r} \cdot \mathbf{s} \equiv \mathbf{r}^T \mathbf{s} ,
Line 225: Line 236:
:<math>
:<math>
\big( \mathbb{A}\mathbb{B} \big)^T = \mathbb{B}^T \mathbb{A}^T.
\big( \mathbb{A}\mathbb{B} \big)^T = \mathbb{B}^T \mathbb{A}^T.
</math>
</ref>
:<math>
\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).
</math>
</math>


Line 257: Line 272:
\boldsymbol{\omega} = \left(\omega_x, \, \omega_y,, \, \omega_z \right)
\boldsymbol{\omega} = \left(\omega_x, \, \omega_y,, \, \omega_z \right)
</math>
</math>
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.
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===
===Inertia moment===
Let us first consider one arbitrary particle, for the sake of clarity we drop the particle index.
Let us first consider one arbitrary particle, for the sake of clarity we drop the particle index.
Line 313: Line 328:
\mathbb{I}(t) \boldsymbol{\omega} .
\mathbb{I}(t) \boldsymbol{\omega} .
</math>
</math>
===Body-fixed angular momentum===
In quantum mechanical applications the angular momentum&mdash;a vector&mdash;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,
:<math>
\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],
</math>
where we used twice that for arbitrary vectors '''a''' and '''b''' and orthogonal matrix with unit determinant,
:<math>
\mathbb{R} \mathbf{a} \times \mathbb{R} \mathbf{b} = \mathbb{R}\left[ \mathbf{a} \times \mathbf{b} \right],
</math>
and where we introduced for temporary use the short-hand notation,
:<math>
\boldsymbol{\omega}' \equiv \mathbb{R}(t)^T \boldsymbol{\omega}.
</math>
In the same way as above, the expression in curly brackets can be worked out:
:<math>
\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}'.
</math>
From this and the earlier result follows parenthetically that the inertia tensor is indeed a tensor,
:<math>
\mathbb{I}(t)= \mathbb{R}(t)\, \mathbb{I}(0)\, \mathbb{R}(t)^T,
</math>
because this rotational relation is often taken as  defining a tensor.
Further,
:<math>
\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}
</math>
The vectors '''L''' and '''&omega;''' 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 <math>\scriptstyle \mathbb{I}(0)</math> is the inertia moment
in the body-fixed frame. Define further,
:<math>
\mathbf{L}^\mathrm{BF} \equiv  \mathbb{R}(t)^T \mathbf{L}
\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad(2)
</math>
and
:<math>
\boldsymbol{\omega}^\mathrm{BF} \equiv  \mathbb{R}(t)^T \boldsymbol{\omega} 
</math>
then we see that in the body-fixed frame angular momentum depends again linearly on angular velocity
through a tensor-vector expression,
:<math>
\mathbf{L}^\mathrm{BF}  = \mathbb{I}(0) \boldsymbol{\omega}^\mathrm{BF},
</math>
where the inertia tensor is constant.
We have shown above that '''L''', which contains the three  components of the physical vector ''angular momentum'' <math>\mathcal{L}</math> with respect to the stationary laboratory frame, is independent of time when there is no external torque.
Hence the physical vector <math>\mathcal{L}</math>, a quantity that exists irrespective of choice of frame, is itself time-independent in absence of torque. Now, '''L'''<sup>BF</sup> contains the components of <math>\mathcal{L}</math>, but  re-expressed with respect to a moving frame. Hence, the components of '''L'''<sup>BF</sup> are time-dependent.  Indeed,  its time-dependence is given by Eq. (2).


[[Category: CZ Live]] [[Category: Stub Articles]] [[Category: Physics Workgroup]]
==References==
<references />[[Category:Suggestion Bot Tag]]

Latest revision as of 16:00, 10 July 2024

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.

PD Image
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 central-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

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{r}} \times \mathbf{v} \equiv \mathbf{v} \times \mathbf{v}= 0 \quad \hbox{and}\quad m \mathbf{r} \times \mathbf{\dot{v}} = m \mathbf{r} \times \mathbf{a} = \mathbf{r} \times \mathbf{F} , }

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 (F = 0) 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 (r parallel to F) is more interesting. In order to see what this means for the potential, we must assume that F is conservative, i.e., its curl vanishes, otherwise the potential cannot be defined. If F is conservative, then F is minus 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 torque 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} = - \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 torque vanishes whenever the potential V(x,y,z) is rotationally invariant (also known as an isotropic or central-symmetric potential). This means that the function V(x,y,z) depends on r only, when expressed in spherical polar coordinates (r, θ, φ).

We can prove the vanishing of the torque more explicitly by expressing the gradient in spherical polar coordinates and writing the cross product as a determinant, 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} =\begin{vmatrix} \mathbf{e}_r &\mathbf{e}_\theta & \mathbf{e}_\phi \\ r & 0 & 0 \\ \frac{\partial V}{\partial r} & 0 &0 \\ \end{vmatrix} = 0. }

Note that the third member in this equation shows that r (the second row) and −F (the third row) are indeed parallel to each other in the case of an isotropic potential. Both vectors possess only a radial coordinate; r by definition and F because the gradient of V has only one non-vanishing component, the radial component.

A well-known example of a system with conserved angular momentum is the earth orbiting 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, i.e., the time derivative of the angular momentum (torque) is zero. 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

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}}_\mathrm{sys} = \mathbf{N}^\mathrm{ext}_\mathrm{sys} = \sum_{i=1}^n \mathbf{r}_i \times \mathbf{F}_i, }

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,

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{tot} \equiv \sum_{i=1}^n m_i \mathbf{s}_i \times \mathbf{\dot s}_i, }

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 is the inertia moment in the body-fixed frame. Define further,

and

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

where the inertia tensor is constant.

We have shown above that L, which contains the three components of the physical vector angular momentum with respect to the stationary laboratory frame, is independent of time when there is no external torque. Hence the physical vector , a quantity that exists irrespective of choice of frame, is itself time-independent in absence of torque. Now, LBF contains the components of , 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,
    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,