In Taylor, read sections 10.5 to 10.6 for today, and 10.7 to 10.8 for Friday.

Since the inertia tensor is symmetric, we can find an orthogonal set of axes where it is diagonal.
These axes are known as principal axes and the diagonal values of the inertia tensor are the principal moments, λ_{i}.

We can easily calculate rotational kinetic energy by exploiting the principle axes and principal moments.

T = ½∑_{i} λ_{i}ω_{i}²

The principal axes are directions where the angular momentum is proportional to the angular velocity

The last expression is simply the definition of the inertia tensor. These yield a way to find the principal moments and principal axes, simply solve the equation

( __I__ - λ )**ω** = 0.

Equations of this type appear commonly in physics and are known as ** I** - λ

det(__I__ - λ__1__) = 0.

This determinant yields a cubic equation in λ which can have 1, 2, or 3 real roots.
The roots are the principal moments, λ_{1}, λ_{2}, and λ_{3}.
To find the principal axes, we solve the three equations

( __I__ - λ_{i} ) **e**_{i} = 0.

for the three corresponding eigenvectors **e**_{1}, **e**_{2}, and **e**_{3}.

Determine the principal moments and axes for a cube rotating about a corner.

The inertia tensor for a cube with sides of length a oriented with its sides parallel to the x, y, and z axes is given in Taylor, and is

[ | 8 | -3 | -3 | ] | |

I = (1/12) Ma² | [ | -3 | 8 | -3 | ] |

[ | -3 | -3 | 8 | ] |

First, find the principal moments by solving the equation

[ | 8 - λ' | -3 | -3 | ] | |

det( - λ) = 0 = (1/12) Ma² detI |
[ | -3 | 8 - λ' | -3 | ] |

[ | -3 | -3 | 8 - λ' | ] |

where I've used λ = (1/12)Ma² λ' to simplify the expression somewhat. The constant outside the determinant plays no role in finding the zeros. To take the determinant, I will expand about the entries in the first row:

(8 - λ')((8 - λ')² - 3²) - 9(8 - λ' + 3) - 9(8 - λ' + 3) = 0

This can be simplified like:

(8 - λ')(8 - λ' - 3)(8 - λ' + 3) - 18(11 - λ') = (11 - λ')[(8 - λ)(5 - λ') - 18] = 0

and expanding the terms in square brackets

(11 - λ')[22 -13λ' + λ'²] = (11 - λ')² (2 - λ') = 0.

Therefore, the principal moments are λ_{1} = λ_{2} = (11/12) Ma² (a double root), and λ_{3} = (1/6) Ma².

Determine the principal axes by solving ( ** I** - λ

[ | 6 | -3 | -3 | ] | [ | e_{3x} |
] | ||

I = (1/12) Ma² | [ | -3 | 6 | -3 | ] | [ | e_{3y} |
] | = 0 |

[ | -3 | -3 | 6 | ] | [ | e_{3z} |
] |

This has the solution **e**_{3} = (1/√3)(1, 1, 1), where I've normalized the vector to unit length.
This direction corresponds to a body diagonal, a diagonal from the corner the cube is rotating about, to the opposite corner, passing through the middle of the cube.

For **e**_{1} and **e**_{2} we have

[ | -3 | -3 | -3 | ] | [ | e_{1x} |
] | ||

I = (1/12) Ma² | [ | -3 | -3 | -3 | ] | [ | e_{1y} |
] | = 0 |

[ | -3 | -3 | -3 | ] | [ | e_{1z} |
] |

where e_{1} can also be e_{2}.
This has an infinite number of solutions for **e**_{1} and **e**_{2}, which can be any pair of orthogonal directions in a plane perpendicular to **e**_{3}.
For instance, we can take **e**_{2} = (1/√2)(1, -1, 0), and then find **e**_{1} from **e**_{1} = **e**_{2} × **e**_{3} = (1/√6)(-1, -1, 2).

To get a flavor of what's involved in solving rotation problems, consider the case of a rotating top subject to a weak torque.
The top has axial symmetry, and calling **e**_{3} the symmetry axis, with **e**_{1} and **e**_{2} perpendicular to it and each other, we know that the inertia tensor is diagonal with principal moments λ_{1} = λ_{2} = λ, and λ_{3}.

Suppose the top begins spinning about its symmetry axis.
In the absence of gravity, there is no torque, and it will continue to spin that way with **L** = λ_{3}**ω** = λ_{3}ω**e**_{3}.
With no torque, the angular momentum will remain constant, equal to this value.

Now turn on gravity, supposing that it causes a weak torque **Γ** = **R**×M**g**.
The magnitude of the torque is Γ = RMg sinθ and is directed perpendicular to the plane containing the **z** and **e**_{3} axes.
The torque will cause the angular momentum to change, but since it is directed perpendicular to the angular momentum, the direction will change but not the magnitude.
The change in **L** will cause ω_{1} and ω_{2} to be small, but non-zero.