arXiv:1807.03867v2 [physics.class-ph] 21 Feb 2019 eetrve)adohrtertcladcalculational paper and theoretical other for advantages. (8) and Ref. review) (see to recent properties due a objects. material recently to classical popular applications rotating their become of have are phase phases features Berry Geometric a surprising of and cation debated myster- being new revealed. still motion, are classical ies and mundane seemingly a nl n a eeoe ntemd18sb Han- Han- by work. the 1980s original is mid Berry’s by the phase inspired in others, Berry developed and of nay was and analog angle classical ana- nay The ‘classical term the to log’. version classical the relegating lsia ytmta a eidcmto.A xml of example a An consider motion. example, periodic simple has a that As system are classical (15). examples Ref. real-world in several calculated of references. phases several features geometric in The basic seen be The can calculations courses. and undergraduate in discussed iae ute nta etadtecneto othe Hannay to by work connection would original which the work the their predated and have in included text not was that angle Hannay in further (14), tigated Ref. (see position.” not original its does § to top return was the time This any that “this at result stated time. the they in has when incommensurability periodic Lifshitz and perfectly Landau not by identified was top Lifshitz’s the and asym- of Landau the In the century. of down 19th Mechanics, motion write the in the to top for first metric expressions the was analytic sub- Jacobi exact the not Evidently yet of apparent, treatments ject. is elementary angle in Hannay discussed the typically top, a of spinning aoia xml fteHna nl sal involves usually angle Hannay pendulum. Foucault’s the of example canonical 7 ae10 oee,ti bevto a o inves- not was observation this However, 120) page 37, h oino pnigojcsi bqios rma From ubiquitous. is objects spinning of motion The emti hss(anyage)aenttypically not are angles) (Hannay phases Geometric 11 eciigteqatmmcaia ffc and effect mechanical quantum the describing 2–4 14 9,10 akadfrhi eidcmne sse napplrYouTub popular a in seen as manner periodic a in forth and back ewrs emti hs,Hna nl,Brypae holo phase, Berry angle, Hannay phase, geometric Keywords: Berry of numbers: analog PACS classical sy the dynamical explore simplest to the particul example of the canonical one a with As formula phas this purposes. geometric of demonstration the consequences hence and some u phase explore explicitly total we for Here th expression of problem. an form concepts classical equations recent this However, Euler on century. the perspective a by than more described for known completely is motion its n motn bevto steidentifi- the is observation important One twsarayapeitdta h motion the that appreciated already was it nltcfruafrteGoercPaeo nAymti To Asymmetric an of Phase Geometric the for formula Analytic h oino adesinn nsaehsa d eair I behavior. odd an has space in spinning handle a of motion The xmlshv eido er’ original Berry’s on relied have Examples .INTRODUCTION I. 12,13 oee,ee nteforce-free the in even However, h ahlcUiest fAmerica of University Catholic The 7 n Berry. and ihlsA Mecholsky A. Nicholas irosSaeLaboratory State Vitreous ahntn C20064 DC Washington, 5,7,15–18 Dtd 1/21/2019) (Dated: 6,7 The 5–7 11 ie akLv,pbihdpapers published Levi, Mark One time, ‘active’ forcing. and this ‘passive’ to the Hamilto- due examples cases. two original phase these the extra call to an might returned is then there is nian, pendulum. Foucault system influence. the external the be some If could this to of due example is slowly An Hamiltonian change whose to system Alternatively,forced periodic a ‘geometric’. consider hence may we and time of dependent atta eed ntetm yaisadteohris other Angle, Hannay the the and dynamic called dynamics total the part, time geometric The is the the One on paper. parts. depends this two that of part in combination phase’ a be ‘total is will phase the phase extra as angular This an to as state. referred viewed initial the be from can are difference that this variables Typically the re- for periodic. except almost not state system original the its variables, to periodic In turns the periodic. of not cycle are one angles) Euler an- the However (like the variables vectors). (like other velocity periodic space. angular and are in momentum variables motion gular dynamic force-free the in of top Some asymmetric the is this uafrcmuigtettlpaeo ii bodies rigid of phase total the computing for mula nteps e eae sefreapeRf 2,the 22), an of Ref. phase example total the for for (see formula explicit decades and few connection past the in hand mo- right sign). angular the a (using the assign frame of to path body of rule the the fraction by in angular a vector out is mentum swept the part sphere onto unit geometric projected a The velocity vector. angular momentum the of gle indsldae wp u yteaglrmomentum the angular is part the dynamic Ω The by and vector. out (16)) swept Eq. area from it solid compute signed we (here space where n19,RcadMngmr,adaon h same the around and Montgomery, Richard 1991, In vntog uhwr a endn nti area this in done been has work much though Even oy ue o,rgdbody rigid top, Euler nomy, eteeutoso oint n closed a find to motion of equations the se T phase. reapeo h pnighnl for handle spinning the of example ar ftefrefe smerctpand top asymmetric force-free the of e h qain fmto aebeen have motion of equations the d ∗ stepro fteaglrvlct etrin vector velocity angular the of period the is tm,teaymti o hudbe should top asymmetric the stems, emti hs aealwdanew a allowed have phase geometric e video. e em ouepcel flip unexpectedly to seems t 1 ∆ sa smercltop, asymmetrical an As α = 2  M ET 2 M ET −  6,21 Ω , steitgae an- integrated the is htdrvdafor- a derived that p 19 hti in- is that 13,20 (1) 2 asymmetric top was never published. Part of the goal is given by of this paper is to take the last step and cast the exact expressions for the motion of the asymmetric top in the 5 2 3 3 l2w + l1 framework of the total phase. A closed form expression a = (3a) for the total phase will be presented and explored. We use q 2p√l1 + l2 the example of the T-shaped handle (seen flipping around 5 2 2 2 2 4 3 2l1l2 (3l2w +2l2 +2w )+ l1w + l2 on the International Space Station in a recent YouTube b = (3b) 1 2 video ) to demonstrate some of these expressions, how- q p 2 (l1 + l2) ever the formulas are valid for a general asymmetric top. 1 5 p In some ways, this is the most elementary example of c = w, (3c) total phase available in classical physics. Even in this 2r3 simple system, a wide variety of questions and lines of shown superimposed in Fig. 1. Although this is correct investigation are possible. To that end, an exploration 1 and some observations of the formula are in order. for all l1, l2, and w, here we wish l1 large enough with respect to l2 (and w small enough) so that a>b>c. The paper begins by reviewing the dynamics of the It will turn out that the only important parameters asymmetric top. The and Euler angles for the description of the body are those that determine are determined as functions of time. From these expres- the direction that the points in the sions, the total phase is determined directly and com- moment of inertia space (the polar and azimuthal angles pared to Montgomery’s formula. Finally, some observa- relative to the I1, I2, and I3 axes). Thus the absolute tions about the total phase are made and some conclu- magnitude of the vector of the eigenvalues of the mo- sions are discussed. ment of inertia is not important. We will find it useful to use a moment of inertia that is parameterized by two variables, b and c, the semi-axes of a solid ellipsoid. The third semi-axis of this ellipsoid is determined such that

II. DYNAMICS AND TOTAL PHASE

As an example for visualization, and without loss of generality, we may use a T-shaped handle to explore our dynamics such as the dynamics seen in the video.1 For any three physical moments of inertia there is a corre- sponding T-shaped handle and thus we are losing no gen- erality. Consider a handle of uniform density and unit with a cross piece of length l1, a shaft of length l2, and a square cross-sectional width of w (Fig. 1). The total moments of inertia can be computed using the parallel axis theorem with individual square prisms. This will give a moment of inertia tensor with diagonal elements

2 2 2 2 4 3 6l1l2w +2l1w + l2w (6l1 + w)+ l2 +4l1l2 I1 = 2 (2a) 12 (l1 + l2) 2 3 (l1 +2l2) w + l1 I2 = (2b) FIG. 1: Two equivalent rigid bodies. A T-shaped han- 12 (l1 + l2) dle (green) with l1 = 8, l2 = 4, and w = 1 and the I3 = equivalent solid ellipsoid that has the identical principal mo- 2 2 2 2 2 2 2 4 ments of inertia given by Eqs. (3), with semi-axes (a,b,c) = l1l2 l1 +5w +4l2 + l1 l1 + w + l2w (6l1 + w)+ l2 2 , (4.23281, 3.07318, 0.645497). The ellipsoid is shrunk by 1/4  12 (l1 + l2)  for presentation. (2c) where the coordinate axes have been chosen to have a diagonal moment of inertia tensor. An ellipsoid (with 1 Otherwise a relabeling of a, b, and c is needed so that the mo- semi-axes lengths a, b, and c) with the same moment of ments of inertia are correctly ordered. 3 the moment of inertia vector (composed from the eigen- reference (determined by the principle moments of inertia values) has unit magnitude. The moment of inertia for I1 < I2 < I3 in a right handed coordinate system), the this particular choice is given by, initial angular velocity is given by

2 2 ω I1 = b + c (4a) 0 = (ω01,ω02,ω03). (9)

1 2 4 2 2 4 2 We will show below that the total phase will not depend I2 = b + 2 3b 2b c 3c + c (4b) 2 − − − − on the magnitude of ω0. Thus we shall take ω0 = 1. 1  p  | | I = b2 + 2 3b4 2b2c2 3c4 c2 (4c) The time dynamics of the angular velocity is then de- 3 2 − − − − termined by these initial conditions. In -free mo-  p  tion, there are two constants of motion for the asymmet- with the third semi-axis of the ellipsoid given by ric top. One is the total angular momentum of the top, given by √2 3b4 2b2c2 3c4 b2 c2 a = − − − − − (5) p √2 M = I ω = I2ω2 + I2ω2 + I2ω2 , (10) | · 0| 1 01 2 02 3 03 q This particular choice of semi-axes of the ellipsoid is al- and the total of (assume the center of ways chosen to be labeled so that I1 < I2 < I3. However, mass velocity of the asymmetric top is 0), given by as long as a>b>c in the ellipsoid, this will always be the case. 1 1 E = ω I ω = I ω2 + I ω2 + I ω2 . (11) In the next section we take a small digression to iden- 2 0 · · 0 2 1 01 2 02 3 03 tify the full space of moments of inertia. The motion of the angular momentum vector in space co- ordinates is, of course, fixed since no torque is applied to A. Space of Possible Moments of Inertia 0.8 Not all triplets (I , I , I ) of positive numbers are valid 1 2 3 Prolate Spheroids moments of inertia. To have a valid (physical) moment π/2 of inertia we must satisfy the following relations, π 0.7 2π

I1 + I2 I3 (6) ≥ Sphere I2 + I3 I1 (7) ≥ 0.6 I + I I . (8) 3 1 ≥ 2 2 these may be referred to as the inertial inequalities. I Without loss of generality (we may relabel axes in the 0.5 object if necessary) we may further restrict ourselves to Degenerate the region, I1 < I2 < I3. For labeling purposes, consider 2 2 2 Ellipsoids only the moments of inertia where I1 + I2 + I3 =1. In this case (we will see that these are the only cases we 0.4 (c=0) need consider) we can show the region of possible physi- cal moments of inertia as those moments of inertia inside Oblate Spheroids the red boundary in the I I plane in Fig. 2. The 1 − 2 0.3 curves that bound this region are the curve of all pro- 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 2 1 I1 I1 late spheroids I2 = − , the curve of all oblate r 2 ! spheroids (I2 = I1), and a curve of degenerate ellipsoids FIG. 2: Full parameter space of asymmetric rigid bodies of 1 fixed moment of inertia magnitude equal to 1. I1 and I2 where c = 0, I = 2 3I2 I . Thus picking 2 2 − 1 − 1 are two of the principal moments of inertia labeled so that 0

I I −I − 3( 2 1) 2 − Π I I −I k MT 4M(I3 I1) 1( 3 2) 2 + 2π M > 2EI2 I3 I1I3 I −I M2− EI  ( 3 2)( 2 1) I I I −  1 2 3 r I M2− EI ∆α = 3 2 1 (20)  − ( ) 2  Π I EI −M2 k  1(2 3 )  4M(I3−I1) !  MT + M 2 < 2EI I3 I1I3 I −I EI −M2 2 ( 2 1)(2 3 ) I I I 1 2 3  r   where Π(n m) is the Complete Elliptic Integral of the where k is defined in Eq. (16) (b). Third Kind| defined by It is worthwhile to note that the total phase formula in Eq. (20) seems to depend on E, M, I1, I2, and I3. π/2 dθ However it can be shown that it does not depend on the Π n m = , (21) (1 n sin2 θ)(1 m sin2 θ)1/2 magnitude of the moment of inertia (I2 = I2 + I2 + I2) Z0 1 2 3  − −

6 nor on the magnitude of the initial angular velocity. angular momentum and total rotational energy. Those, To see the independence of the magnitude I, note that along with the moments of inertia are: E = 0.264805 J, 2 2 the energy, Eq. (11), is homogeneous of order 1 in I. M =0.533252 Js, I1 =0.286827 Js , I2 =0.533256 Js , 2 The angular momentum, Eq. (10), is also homogeneous and I3 = 0.795844 Js . The dynamics of ω in the body of order 1 in I. The functions k and T of Eqs. (16) are frame are shown in Fig. 3(b). In this case, the period therefore homogeneous of order 0 in I. It then follows calculated from Eq. (16) (a) is 23.8181 s. We can see that Eq. (20) is independent of the magnitude of the from the shape of the curves that this corresponds visu- moment of inertia. ally to the period of Fig. 3(b). One aspect to note from For the independence of the magnitude of ω0, let us this plot is that the dashed green curve shows the body consider the initial polar and azimuthal angles of the an- y-component of the angular velocity (the axis of the in- gular velocity (ω0) with respect to the body axes (θω and termediate moment of inertia) suddenly reverses twice φω) fixed for the moment. Note that in this case, the an- describing the flip known as the Dzhanibekov effect or 14,23,26,27 gular momentum is homogeneous of order 1 in ω0 , and the tennis racket theorem. Interestingly, most the energy is homogeneous of order 2 in ω . This| means| classic text in dynamics describe this behavior as ‘unsta- | 0| that k is homogeneous of order 0 in ω0 and T is homo- ble’, however the motion is completely determined and is geneous of order -1 in ω . Then looking| | at the total periodic, eventually coming back to the ‘unstable’ point. | 0| phase formula, Eq. (20), this means that the total phase We may now compute some of the quantities that have is independent from the magnitude of the initial angular been developed. The dynamic part of the total phase ω velocity, 0 . 2ET | | from Eq. (1) is 23.6555. If we integrate the Considering that the total phase formula is itself a con- M   stant of motion (it is composed of constants of motion), projected angular velocity (from Eq. (A.1) or numerically the initial angles, θω and φω, which set the constants of integrating Eq. (15)) along the angular momentum axis the motion (E and M) along with the direction of the moment of inertia (or simply I1 and I2 along with the 2 2 2 constraint that I1 + I2 + I3 = 1) are the only parameters needed to determine the total phase. Thus every point along the closed curve on the Binet ellipsoid has the same total phase since every point has the same quantities of E, M, I1, I2, and I3. Thus curves of the same total phase foliate the Binet ellipsoid. An example2 of this can be seen in Fig. 4 for 2 2 the case of (I1, I2, I3) = (0.286827 Js , 0.533256 Js , 0.795844 Js2). Here the color of the curve corresponds to the total phase modulo 2π. In the next section, we check Eq. (20) by evaluating Eq. (20), by calculating Eq. (1) numerically, and by solving for a rotation matrix after one period.

D. Demonstration of formula

For illustrative purposes, we choose a handle shape where l1 = 8, l2 = 4, and w = 1 as shown in Fig. 1. We further choose the total mass (from a uniform density) to be such that the moment of inertia has unit magnitude. The initial angular velocity is chosen so that it has unit magnitude, and the direction is given by a polar angle FIG. 4: Total phase (mod 2π) of the body given in Fig. 1 (measured from the body z axis), θω = π/2 0.15, and 2 − for an asymmetric top with (I1, I2, I3) = (0.286827 Js , an azimuthal angle, φω = π/2 0.2, measured from body 2 2 − 0.533256 Js , 0.795844 Js ). For each initial angle, a given x axis. This corresponds to a large initial velocity near curve on the Binet ellipsoid is represented. The Binet ellip- the middle moment of inertia, I2. This sets the total soid is the ellipsoid of constant energy in the angular momen- tum space. Given an initial angular momentum, the resulting dynamics traces out a closed curve on this surface. The final angular mismatch (mod 2 π) is the total phase and corre- 2 We use Js2 in the units of the moment of inertia as opposed to sponding colors relate to corresponding values from 0 to 2π. kg m2 to make transparent the relationship between the energy Note that there are multiple curves where the phase is near (J), angular momentum (Js), moments of inertia (Js2), and 0. angular velocity (s−1) in Eqs. (16) and (20). 7 for one period, we get 0 < θω < π/2, this covers all initial conditions (using symmetry of the angle, direction of spin, and indepen- T ω M dence from ω0 and I ). Equation (20) has been checked M· dt = 23.6555, (22) | | | | 0 exhaustively against Eq. (1) for 292 sets of moments of Z | | inertia equally spaced throughout the region in Fig. 2 which matches the dynamic part. If we compute the with 22 angles between 0 < θω < π/2 for each moment solid angle of intersection between the sphere of angu- of inertia choice. This corresponds to 6446 initial con- lar momentum and the Binet ellipsoid numerically using ditions. In all cases, both Eq. (20) and Eq. (1) give the Mathematica (seen in Fig. 3 (a)), we get the signed area same value. Ω=2.1253. III. DISCUSSION OF FORMULA This gives a calculated total phase of 21.5302. The formula from Eq. (20) gives a value of 21.5303. The dif- A. Symmetries of the total phase ference is mostly due to the numerical discretization of the area of intersection and can be made to be smaller with a finer discretization. As mentioned before, the total phase formula has no- Another way of checking both formulas is by numeri- table symmetries including its independence from I = 2 2 2 ω 2 2 2 | | cally integrating the Euler equations over one period and I1 + I2 + I3 and 0 = ω01 + ω02 + ω03. solving a system of equations for the rotation matrix af- | |I pIndependence from pfollows from the observation ter one period. The rotation matrix would give the ge- that I is directly proportional| | to the total mass of a rigid ometric rotation of the initial system compared to the body.| Changing| I amounts to changing the mass scale. final orientation after one period. In this case, the ro- If the mass scale| were| changed, the dynamics would be tation around the angular momentum vector is given by unchanged, especially the geometric quality of the rota- 2.68074 which is precisely the value of the total phase tion. formula (21.5303), modulo 2π. To qualitatively understand the independence of Eq. To further simplify the total phase formula, we may use (20) from ω0 , we see that changing the magnitude of the observation that there are related initial conditions | | ω0 changes the timescale for the dynamics, but geo- that describe the same total phase (since the total energy metrically,| | does not affect the phase. If we had made and magnitude of the angular momentum is the same a video recording of the spinning, and played it back at throughout the evolution). These lie on the same closed a slower speed, the geometrical quality of the rotation curve on the Binet ellipsoid (Fig. 4). Without loss of would also be unchanged. generality, we may choose the initial conditions to lie on Thus we only need three parameters to specify the to- the space x z plane so that the azimuthal angle of the − 3 tal phase: two moments of inertia, I1 and I2 (or two initial angular velocity, φω is equal to 0. Thus for the ellipsoidal semi-axes (b and c)), and the polar angle of remainder of the paper we take the initial angular velocity in the body frame, θω.

ω0 = (sin θω, 0, cos θω). (23) We may explicitly plug in the expressions for E and M (Eqs. (11) and (10)) along with Eq. (23) and the defini- If we consider all combinations of moments of inertia tions for T and k (Eqs. (16)) into the total phase formula in the region indicated in Fig. 2, along with all initial (Eq. (20)) to give 3 θω It might be objected that by the time the initial condition ( rescale our unit of time to give |ω0| = 1. This requires scaling and φω) makes it to the x − z plane, the angular velocity may the time by the ratio of the periods. The Binet ellipsoid and no longer be of unit magnitude (in fact, in general, it will not). angular momentum sphere are simply magnified or reduced, but ω However since Eq. (20) is independent of | 0|, we may renor- the polhode curve shape is maintained. malize and retain the same ∆α, but not necessarily the same exact dynamics. Equivalently (as discussed in §III A), we may

I I2−I2(1−sec2 θω) 2 2( 3 1 ) 2 I3 2 2 4 β 2 I1K β tan θω + (I3 I1)Π 2 β β tan θω 2π M > 2EI2 I I3(I3−I1)(I2−I1) I ∆α =  r 1 − − 1 − , (24) I I2+I2 cot2 θω  2 2  2 2   2( 1 3 ) cot θω  I3 cot θω cot θω 2 4 2 I1K + (I3 I1)Π 2 β M < 2EI2  I I3(I3−I1)(I2−I1) β I β r 1 − − 1        where and K and Π are defined in Eqs. (17) and (21), respec- tively. Here we might have eliminated I3 as well using 2 2 I (I I ) I3 = 1 I1 I2 , but that would perhaps not benefit β = 1 2 − 1 , (25) − − I (I I ) the exposition. 3 3 − 1 p 8

B. Minimum Total Phase When the total phase modulo 2π is equal to 0, trajecto- ries that are started on these curves are closed and each Figure 5 shows two plots for the above demonstrated Euler angle (mod 2π) is periodic with a period whose case. The solid red curve is the total phase for a given integer multiple is the period of the angular velocity vec- initial condition described by the polar angle (θω) of the tor. As the trajectories on the Binet ellipsoid approach 2 angular velocity with respect to the body z axis and an the curve M = 2EI2, the total phase increases with- azimuthal angle of zero (φω ) with respect to the body out bound. One question that is apparent is what is the x axis. The total phase (solid red) is the difference be- minimum total phase for a given object? For example are tween the dynamic phase (dotted green) and the geomet- there objects (I1 and I2 or b and c values) where there ric phase (dot-dashed blue). The dashed black curve is are no initial conditions that give rise to a total phase of the corresponding time period, T , for the angular veloc- exactly 2π? ity or angular momentum vectors. For a given object, In Fig. 2 we see that the minimum total phase can be as the total phase, dynamic phase, and period diverge at high as desired. Looking at all objects, the largest mini- an angle that corresponds to the separatrix for the inter- mum total phase is unbounded for objects with moments section of the Binet ellipsoid and the moment of inertia of inertia approaching I1 = I2 = I3 = √3/3. There is an sphere. The critical angle is given by the condition apparent paradox in this behavior. If all moments of in- ertia are equal and the total angular momentum is fixed, 2 M =2EI2 (26) then the body spins around whatever axis it is initial spinning about. The total phase would seem to be zero. which reduces to However the total phase formula goes to infinity. The resolution may be seen in the behavior of the period T . −1 θcrit = cot β , (27) When the moments of inertia each approach √3/3, the period is seen to diverge. Thus with T going to infinity, p  with β defined in Eq. (25). the total phase is permitted to correspondingly diverge. This critical angle is displayed as a gray vertical line on Note that when I1 = I2 = I3 = √3/3 exactly, the object Figs. 5 and 6. Note that the minimum value for the total is no longer dynamically asymmetric (an example of this phase occurs at the extremes of 0 and π. This appears is a sphere). to be a generic behavior in all of the initial conditions On the other hand, the asymmetric top with the small- investigated. It might not be too difficult to differentiate est minimum possible total phase for any initial condition Eqs. (24) and investigate the sign of the derivative, but is a total phase of 0 for objects with moments of inertia may not add much. Additionally, near the critical angle, approaching I1 = 0 and I2 = I3 = √2/2. An example of the total phase diverges. this is a thin solid rod. If we look at initial conditions near the body z or x axes (θω 0 or θω π/2, respectively), we may work out the minimum≈ total≈ phase directly from linearizing the Euler equations (Eqs. (15)). The Euler equations reduce to a simple harmonic oscillator equation where the period may be identified directly. Since in this case, the geometric part of the total phase is zero (the angular momentum does not enclose any solid angle), the total phase is given by the dynamic part, for which there is a simple formula (see Eq. (1)). In the case where θω 0 the minimum total phase reduces to ≈

I I ∆α =2π 1 2 . (28) min,z (I I )(I I ) s 3 − 2 3 − 1 FIG. 5: For the demonstrated example of Fig. 3, the solid red In the case where θω π/2 the minimum total phase is curve indicates the total phase for different initial conditions given by ≈ with varying polar angle of the angular velocity. A polar angle of 0 is an angular velocity along the body z axis. The total I3I2 phase (solid red) is the difference between the dynamic phase ∆α =2π . (29) min,z (I I )(I I ) (dotted green) and the geometric phase (dot-dashed blue). s 2 − 1 3 − 1 Here we see that regardless of the initial condition, the total phase will be larger than about 2π radians. The black dashed These formulas may also be derived by expanding Eq. curve plots the corresponding period of the angular velocity (24) near θω 0 and θω π/2. It is worth noting or angular momentum vector in seconds. that the total≈ phase near these≈ extremes is composed of entirely the dynamic part of the total phase. For initial 9 conditions away from these extremes, the total phase is a in the plot, but for the even multiples of 2 π, (4π,8π, 12π, combination of the dynamic and geometric parts as can etc.), the herpolhodes also have a horizontal symmetry. be seen in Fig. 5. Additionally, changing the magnitude of ω0 expands or The locus of points where Eq. (28) and Eq. (29) are contracts the herpolhode image, but otherwise| | keeps it equal is exactly the location of the kinks in the contours unchanged. of Fig. 2. It is unclear what, if any, connection these closed curves and the associated have to the Bohr- Sommerfeld quantization in quantum mechanics.

IV. APPLICATIONS AND POSSIBLE CONNECTIONS

The primary purpose of this paper was to use the asym- metric top as an intuitive and natural example of total phase. However we give three examples of the connec- tions that could be made with other areas of physics. Quantization of the asymmetric top28 leads to quan- tized energy levels of the quantum mechanical asymmet- ric top. How this corresponds to the closed orbits of the last section still needs to be studied. When molecules are no longer symmetric, the asym- metry leads to observable spectroscopic modes that allow FIG. 6: Total phase for moments of inertia given by (I1, I2, I3) 29–31 2 2 2 for analysis and characterization. (See Ch. 1, 4 of = (0.26 Js , 0.55 Js , 0.793662 Js ). Initial angular velocity § angle is given as a polar angle from the body z axis. When the Ref. (30) and Ch. 4 of Ref. (31)) Understanding the dy- total phase is an integer multiple of 2 π, the herpolhode (curve namics of the asymmetric top classically gives us an un- of the terminus of the angular velocity in space coordinates) derstanding of the connection with the rotational modes is a close plane curve. Various examples are plotted near the of an asymmetric molecule. Making further connections intersection of multiples of 2π and the total phase curve. with the total phase arising from simple rotation of these molecules could lead to further insight. Large nuclei such as 135Pr and 163Lu, provide exam- ples of asymmetric top systems in nuclear physics.32,33 Through X-ray spectroscopy, the so-called ‘wobble’ C. Closed Herpolhodes bands provide a way to measure angular momentum tran- sitions of spinning nuclei. The connection to closed total We may also use the total phase formula (Eq. (20)) phase orbits still needs to be made. to produce curves of closed herpolhodes. A herpolhode is the curve traced out by the terminus of the angular velocity vector in the space coordinates. In general these V. CONCLUSIONS do not close after a period. However if we use initial con- ditions that have a total phase that is an integer multiple In this paper, we have found the total phase formula of 2π, we find interesting closed curve patterns. (and thus trivially the geometric phase) explicitly in Figure 6 shows the total phase for all initial polar an- terms of special functions. This leads immediately to ob- gles (θω) of the angular velocity (with φω = 0). The servations about the total phase for the asymmetric top particular object (whose moments of inertia are given by in force-free motion. The total phase of a rigid object 2 2 2 (I1, I2, I3) = (0.26 Js , 0.55 Js , 0.793662 Js )) has ini- is independent of the magnitude of the moment of iner- tial conditions with closed herpolhode trajectories. The tia, independent of the magnitude of the angular velocity, corresponding angles that result in the smallest closed and all initial conditions along the path of the angular herpolhodes are (from left to right) 0.9562, 1.004, 1.011, velocity vector in space have the same total phase and 1.012, 1.014, 1.021, 1.061, 1.2742 radians. These are ex- geometric phase. Thus for a given , the total actly the angles where the total phase is a multiple of phase, dynamic phase, and geometric phase only depend 2π. In this case the multiples are 2, 3, 4, 5, 5, 4, 3, and on a single angle, the polar angle that the initial angu- 2 respectively. Note that closed herpolhodes are possible lar velocity vector makes when it crosses the x z plane for all integer multiples of 2π greater than 2. The critical (for instance). Additionally, in the parameter− space of angle is 1.01262. Some of the closed herpolhode figures the moments of inertia, each class of asymmetric rigid are depicted in Fig. 6. bodies has a minimum total phase; that is, a minimum Two other comments about closed herpolhodes should angle by which an object must rotate after undergoing be made. There is vertical symmetry for all herpolhodes one period of the angular velocity vector in space. For 10 some objects, the total phase can be as high as desired. for all I1 < I2 < I3 and any E and M. Closed trajectories were also highlighted and closed her- We may further express the Euler angles as explicit polhodes were illustrated and discussed briefly. functions of time. We choose the Euler angle convention of both Refs. (23) and (14) (Figure 4.7, p. 152 and Fig. 47 35, respectively). This describes a sequence of intrinsic ′ ′′ Acknowledgments elemental§ rotations about the body z-x -z axes. In this convention, two of the Euler angles (θ and ψ) are given The author would like to thank Richard Cecil for many in terms of the angular velocity functions (see Ref. 14, conversations relating to geometric phase and other top- Eq. 37.14), ics. Additionally, the author would like to thank the in- I ω put of Richard Montgomery for many observations and cos θ = 3 3 (A.4a) suggesting the line of analysis for Eqs. (28) and (29). M I ω This work was funded by the Vitreous State Laboratory. tan ψ = 1 1 . (A.4b) I2ω2 The final Euler angle, φ, is given as a first-order non- Appendix: Explicit Derivation of the Time Dependence of the Euler Angles linear differential equation as follows. Using Eq. 37.16 of Ref. 14, we may substitute in the expressions of Eq. 37.6 of Ref. 14. Simplifying and using The solution of the time dependence of the angular ve- the explicit expressions of the angular velocity in terms locity for an asymmetric top in force-free motion is given of Jacobi elliptic functions (Eq. 37.10, Ref. 14 and above in Eq. 37.10 of Ref. 14, reproduced here for convenience: Eqs. (A.1)), we arrive at the expression of the time rate of change of φ,

2EI M 2 dφ M A ω (t)= 3 cn(Dt k2) (A.1a) 1 − = + 2 2 (A.5) s I1(I3 I1) | dt I3 B + C sn (Dt κ ) − | 2 2EI3 M 2 with ω2(t)= − sn(Dt k ) (A.1b) s I2(I3 I2) | − A = M(I1 + I3)(I3 I2) < 0 (A.6a) 2 − − M 2EI1 2 B = I1I3(I3 I2) > 0 (A.6b) ω3(t)= − dn(Dt k ). (A.1c) − s I3(I3 I1) | C = I2(I I ) > 0 (A.6c) − 3 2 − 1 2 2 Here sn(u k ) is the Jacobi Elliptic Sine function defined (I3 I2)(M 2EI1) by inverting| the Jacobi Elliptic Integral of the First Kind, D = − − , (A.6d) s I1I2I3 F, through the pair of equations and κ defined in Eq. (A.3). This is valid for all initial 2 sin φ = sn(u k ) (A.2a) conditions. | φ dt When we integrate this, we get u = F(φ k2)= , (A.2b) 2 | 0 1 k2 sin t t Dt Z − dφ M A du φ = dt = t + 2 2 2 , p dt I3 BD 1 α sn (u κ ) and cn(u k2) = 1 sn2(u k2) and dn(u k2) = Z0   Z0 − | | − | | (A.7) 1 k2sn2(u k2). The definition of D is given in Eq. p with (A.6)(d).− | p 2 The above equations (Eqs. (A.1)) are valid when M > 2 C 2 α = . (A.8) 2EI2. When M < 2EI2, all indices 1 3 (including in −B the definition of k and ω). Alternatively,↔ we may use Eqs. (A.1) for all initial conditions without making any index We may use Eq. 400.01 (in the definition of the Incom- 2 plete Elliptic Integral of the Third Kind of Ref. 34) changes (extending k > 1 in the region M < 2EI2). This is the most convenient, but care must be taken. For u dx the rest of the Appendix, we use this extended definition Π α2; φ κ2 = , (A.9) 1 α2 sn2 (x κ2) of k to simplify the notation but at the end of the Ap- Z0 − | pendix produce equations that are valid in the notation  where u is defined in Eq. (A.2) (a), to do the above in- of the text. To distinguish it, we use κ. Thus tegral. The exact expressions for each region of initial 2 2 conditions (M > 2EI2 and M < 2EI2) must be done 2 (I2 I1)(2EI3 M ) with care to avoid problems with the definitions of built- κ = − 2 − , (A.3) (I3 I2)(M 2EI1) in functions. s − − 11

2 For the case where initial conditions give M > 2EI2 lus transformation for the Jacobi Elliptic Sine function (0 <κ< 1), we get (162.01 p. 39 of Ref. 34) and definition 22.16.1 of Ref. 35 to give M φ(t)= t I3 M φ(t)= t M(I3 I1) I3(I2 I1) + − Π − ; am(Dt κ2) κ2 , I3 I1I3D −I1(I3 I2) |  −  M(I3 I1) 1 Π I3(I2 I1) am −2 −2 (A.10) + − 2 − ; (κDt κ ) κ . I1I3D κ −κ I1(I3 I2) |  −  (A.11) where am u m is the Jacobi Amplitude. If u = F φ m then φ = am u m , where F is the Incomplete Elliptic   Integral of the First Kind defined in Eq. (A.2) (b). Again, As before κ is defined in Eq. (A.3), and D is given in Eq.  κ is defined in Eq. (A.3), and D is given in Eq. (A.6)(d). (A.6)(d). 2 2 For the case M < 2EI2 (κ> 1) we must use the recip- Note that in the region M < 2EI2, κ = 1/k with k rocal modulus transformation on the Incomplete Elliptic defined through Eq. (16) (b). Thus we may rewrite Eqs. Integral of the Third Kind (19.7.4 of Ref. 35, p. 492 or (A.10) and (A.11) with notation in the rest of the paper 162.02 p. 39 of Ref. 34) as well as the reciprocal modu- as

− − M M(I3 I1) I3(I2 I1) 2 2 2 t + Π − ; am(Dt k ) k M > 2EI2 I3 I1I3D − I1(I3 I2) | φ(t)= − 2 − . (A.12) M M(I3 I1) k I3(I2 I1) Dt 2 2 2  t + kΠ − ; am k  k M < 2EI2  I3 I1I3D − I1(I3 I2) k |    

∗ Electronic address: [email protected]; Corresponding 18 J. E. Marsden et al., Lectures on mechanics, vol. 174 (Cam- Author bridge University Press, 1992). 1 Plasma Ben, Dancing T-handle in zero-g, HD, URL 19 J. Robbins, Journal of Physics A: Mathematical and The- https://www.youtube.com/watch?v=1n-HMSCDYtM. oretical 49, 431002 (2016). 2 H. Moffatt, Nature 404, 833 (2000). 20 J. B. Hart, R. E. Miller, and R. L. Mills, American Journal 3 M. A. Jalali, M. S. Sarebangholi, and M.-R. Alam, Physical of Physics 55, 67 (1987). Review E 92, 032913 (2015). 21 M. Levi, Archive for rational mechanics and analysis 122, 4 M. Berry and P. Shukla, European Journal of Physics 32, 213 (1993). 115 (2010). 22 L. Bates, R. Cushman, and E. Savev, Zeitschrift f¨ur Ange- 5 M. Berry and J. Hannay, Journal of Physics A: Mathemat- wandte Mathematik und Physik (ZAMP) 56, 183 (2005). ical and General 21, L325 (1988). 23 H. Goldstein, C. Poole, and J. Safko, 6 R. Montgomery, Am. J. Phys 59, 394 (1991). (Addison Wesley, San Francisco, 2002). 7 J. H. Hannay, Journal of Physics A: Mathematical and 24 E. T. Whittaker, A treatise on the analytical dynamics General 18, 221 (1985). of particles and rigid bodies (Cambridge University Press, 8 D. Xiao, M.-C. Chang, and Q. Niu, Reviews of modern 1988). physics 82, 1959 (2010). 25 J. W. Zwanziger, M. Koenig, and A. Pines, Annual Review 9 R. Resta, Journal of Physics: Condensed Matter 12, R107 of Physical Chemistry 41, 601 (1990). (2000). 26 M. S. Ashbaugh, C. C. Chicone, and R. H. Cushman, Jour- 10 C. A. Mead, Reviews of modern physics 64, 51 (1992). nal of Dynamics and differential Equations 3, 67 (1991). 11 M. V. Berry, Proceedings of the Royal Society of London. 27 M. Levi, Classical mechanics with calculus of variations A. Mathematical and Physical Sciences 392, 45 (1984). and optimal control: an intuitive introduction, vol. 69 12 J. Anandan, Nature 360, 307 (1992). (American Mathematical Soc., 2014). 13 M. Berry, Scientific American 259, 46 (1988). 28 V. Postell and T. Uzer, Physical Review A 41, 4035 (1990). 14 L. Landau and E. Lifshitz, Mechanics, vol. 1 of Course of 29 J. W. Blaker, M. Sidran, and A. Kaercher, Tech. Rep., Theoretical Physics (1976). Grumman Aircraft Engineering Corp Bethpage NY, Re- 15 J. Lawson and M. Rave, Mathematics Magazine 89, 105 search Dept (1962). (2016). 30 G. Herzberg, Infrared and Raman spectra of polyatomic 16 S. Gil, American Journal of Physics 78, 384 (2010). molecules (D. Van Nostrand Company.; New York, 1945). 17 J. Jos´eand E. Saletan, Classical dynamics: a contemporary 31 C. H. Townes and A. L. Schawlow, Microwave spectroscopy approach (2000). (Courier Corporation, 2013). 12

32 S. Frauendorf and F. D¨onau, Phys. Rev. C 89, 014322 34 P. Byrd and M. Friedman, Berlin, Heidelberg. New York (2014). (1971). 33 S. Ødeg˚ard, G. B. Hagemann, D. R. Jensen, 35 F. W. Olver, D. W. Lozier, R. F. Boisvert, and C. W. M. Bergstroem, B. Herskind, G. Sletten, S. Toer- Clark, NIST handbook of mathematical functions hardback maenen, J. Wilson, P. Tjøm, I. Hamamoto, et al., Physical and CD-ROM (Cambridge University Press, 2010). review letters 86, 5866 (2001).