Foucault pendulum

One of the easiest examples is the Foucault pendulum. An easy explanation in terms of geometric phases is given by Wilczek and Shapere:^[7]

How does the pendulum precess when it is taken around a general path C? For transport along the equator, the pendulum will not precess. [...] Now if C is made up of geodesic segments, the precession will all come from the angles where the segments of the geodesics meet; the total precession is equal to the net deficit angle which in turn equals the solid angle enclosed by C modulo 2π. Finally, we can approximate any loop by a sequence of geodesic segments, so the most general result (on or off the surface of the sphere) is that the net precession is equal to the enclosed solid angle.

To put it in different words, there are no inertial forces that could make the pendulum precess, so the precession (relative to the direction of motion of the path along which the pendulum is carried) is entirely due to the turning of this path. Thus the orientation of the pendulum undergoes parallel transport. For the original Foucault pendulum, the path is a circle of latitude, and by the Gauss–Bonnet theorem, the phase shift is given by the enclosed solid angle.^[8]

Derivation

Parallel transport of a vector around a closed loop on the sphere: The angle by which it twists, α, is proportional to the area inside the loop.

In a near-inertial frame moving in tandem with the Earth, but not sharing the rotation of the Earth about its own axis, the suspension point of the pendulum traces out a circular path during one sidereal day.

At the latitude of Paris, 48 degrees 51 minutes north, a full precession cycle takes just under 32 hours, so after one sidereal day, when the Earth is back in the same orientation as one sidereal day before, the oscillation plane has turned by just over 270 degrees. If the plane of swing was north–south at the outset, it is east–west one sidereal day later.

This also implies that there has been exchange of momentum; the Earth and the pendulum bob have exchanged momentum. The Earth is so much more massive than the pendulum bob that the Earth's change of momentum is unnoticeable. Nonetheless, since the pendulum bob's plane of swing has shifted, the conservation laws imply that an exchange must have occurred.

Rather than tracking the change of momentum, the precession of the oscillation plane can efficiently be described as a case of parallel transport. For that, it can be demonstrated, by composing the infinitesimal rotations, that the precession rate is proportional to the projection of the angular velocity of Earth onto the normal direction to Earth, which implies that the trace of the plane of oscillation will undergo parallel transport. After 24 hours, the difference between initial and final orientations of the trace in the Earth frame is α = −2π sin φ, which corresponds to the value given by the Gauss–Bonnet theorem. α is also called the holonomy or geometric phase of the pendulum. When analyzing earthbound motions, the Earth frame is not an inertial frame, but rotates about the local vertical at an effective rate of 2π sin φ radians per day. A simple method employing parallel transport within cones tangent to the Earth's surface can be used to describe the rotation angle of the swing plane of Foucault's pendulum.^[9][10]

From the perspective of an Earth-bound coordinate system (the measuring circle and spectator are Earth-bounded, also if terrain reaction to Coriolis force is not perceived by spectator when he moves), using a rectangular coordinate system with its x-axis pointing east and its y-axis pointing north, the precession of the pendulum is due to the Coriolis force (other fictitious forces as gravity and centrifugal force have not direct precession component, Euler's force is low because Earth's rotation speed is nearly constant). Consider a planar pendulum with constant natural frequency ω in the small angle approximation. There are two forces acting on the pendulum bob: the restoring force provided by gravity and the wire, and the Coriolis force (the centrifugal force, opposed to the gravitational restoring force, can be neglected). The Coriolis force at latitude φ is horizontal in the small angle approximation and is given by

$${\displaystyle {\begin{aligned}F_{c,x}&=2m\Omega {\dfrac {dy}{dt}}\sin \varphi \\F_{c,y}&=-2m\Omega {\dfrac {dx}{dt}}\sin \varphi \end{aligned}}}$$

where Ω is the rotational frequency of Earth, F_c,x is the component of the Coriolis force in the x-direction and F_c,y is the component of the Coriolis force in the y-direction.

The restoring force, in the small-angle approximation and neglecting centrifugal force, is given by

$${\displaystyle {\begin{aligned}F_{g,x}&=-m\omega ^{2}x\\F_{g,y}&=-m\omega ^{2}y.\end{aligned}}}$$

Graphs of precession period and precession per sidereal day vs latitude. The sign changes as a Foucault pendulum rotates anticlockwise in the Southern Hemisphere and clockwise in the Northern Hemisphere. The example shows that one in Paris precesses 271° each sidereal day, taking 31.8 hours per rotation.

Using Newton's laws of motion this leads to the system of equations

$${\displaystyle {\begin{aligned}{\dfrac {d^{2}x}{dt^{2}}}&=-\omega ^{2}x+2\Omega {\dfrac {dy}{dt}}\sin \varphi \\{\dfrac {d^{2}y}{dt^{2}}}&=-\omega ^{2}y-2\Omega {\dfrac {dx}{dt}}\sin \varphi .\end{aligned}}}$$

Switching to complex coordinates z = x + iy, the equations read

$${\displaystyle {\frac {d^{2}z}{dt^{2}}}+2i\Omega {\frac {dz}{dt}}\sin \varphi +\omega ^{2}z=0\,.}$$

To first order in Ω/ω this equation has the solution

$${\displaystyle z=e^{-i\Omega \sin \varphi t}\left(c_{1}e^{i\omega t}+c_{2}e^{-i\omega t}\right)\,.}$$

If time is measured in days, then Ω = 2π and the pendulum rotates by an angle of −2π sin φ during one day.

Polarized light in an optical fiber

A second example is linearly polarized light entering a single-mode optical fiber. Suppose the fiber traces out some path in space, and the light exits the fiber in the same direction as it entered. Then compare the initial and final polarizations. In semiclassical approximation the fiber functions as a waveguide, and the momentum of the light is at all times tangent to the fiber. The polarization can be thought of as an orientation perpendicular to the momentum. As the fiber traces out its path, the momentum vector of the light traces out a path on the sphere in momentum space. The path is closed, since initial and final directions of the light coincide, and the polarization is a vector tangent to the sphere. Going to momentum space is equivalent to taking the Gauss map. There are no forces that could make the polarization turn, just the constraint to remain tangent to the sphere. Thus the polarization undergoes parallel transport, and the phase shift is given by the enclosed solid angle (times the spin, which in case of light is 1).

Stochastic pump effect

A stochastic pump is a classical stochastic system that responds with nonzero, on average, currents to periodic changes of parameters. The stochastic pump effect can be interpreted in terms of a geometric phase in evolution of the moment generating function of stochastic currents.^[11]

Spin 1⁄2

The geometric phase can be evaluated exactly for a spin-1⁄2 particle in a magnetic field.^[1]

Geometric phase defined on attractors

While Berry's formulation was originally defined for linear Hamiltonian systems, it was soon realized by Ning and Haken^[12] that similar geometric phase can be defined for entirely different systems such as nonlinear dissipative systems that possess certain cyclic attractors. They showed that such cyclic attractors exist in a class of nonlinear dissipative systems with certain symmetries.^[13] There are several important aspects of this generalization of Berry's phase: 1) Instead of the parameter space for the original Berry phase, this Ning-Haken generalization is defined in phase space; 2) Instead of the adiabatic evolution in quantum mechanical system, the evolution of the system in phase space needs not to be adiabatic. There is no restriction on the time scale of the temporal evolution; 3) Instead of a Hermitian system or non-hermitian system with linear damping, systems can be generally nonlinear and non-hermitian.

Exposure in molecular adiabatic potential surface intersections

There are several ways to compute the geometric phase in molecules within the Born–Oppenheimer framework. One way is through the "non-adiabatic coupling ${\displaystyle M\times M}$ matrix" defined by

$${\displaystyle \tau _{ij}^{\mu }=\langle \psi _{i}|\partial ^{\mu }\psi _{j}\rangle ,}$$

where ${\displaystyle \psi _{i}}$ is the adiabatic electronic wave function, depending on the nuclear parameters ${\displaystyle R_{\mu }}$. The nonadiabatic coupling can be used to define a loop integral, analogous to a Wilson loop (1974) in field theory, developed independently for molecular framework by M. Baer (1975, 1980, 2000). Given a closed loop ${\displaystyle \Gamma }$, parameterized by ${\displaystyle R_{\mu }(t),}$ where ${\displaystyle t\in [0,1]}$ is a parameter, and ${\displaystyle R_{\mu }(t+1)=R_{\mu }(t)}$. The D-matrix is given by

$${\displaystyle D[\Gamma ]={\hat {P}}e^{\oint _{\Gamma }\tau ^{\mu }\,dR_{\mu }}}$$

(here ${\displaystyle {\hat {P}}}$ is a path-ordering symbol). It can be shown that once ${\displaystyle M}$ is large enough (i.e. a sufficient number of electronic states is considered), this matrix is diagonal, with the diagonal elements equal to ${\displaystyle e^{i\beta _{j}},}$ where ${\displaystyle \beta _{j}}$ are the geometric phases associated with the loop for the ${\displaystyle j}$-th adiabatic electronic state.

For time-reversal symmetrical electronic Hamiltonians the geometric phase reflects the number of conical intersections encircled by the loop. More accurately,

$${\displaystyle e^{i\beta _{j}}=(-1)^{N_{j}},}$$

where ${\displaystyle N_{j}}$ is the number of conical intersections involving the adiabatic state ${\displaystyle \psi _{j}}$ encircled by the loop ${\displaystyle \Gamma .}$

An alternative to the D-matrix approach would be a direct calculation of the Pancharatnam phase. This is especially useful if one is interested only in the geometric phases of a single adiabatic state. In this approach, one takes a number ${\displaystyle N+1}$ of points ${\displaystyle (n=0,\dots ,N)}$ along the loop ${\displaystyle R(t_{n})}$ with ${\displaystyle t_{0}=0}$ and ${\displaystyle t_{N}=1,}$ then using only the j-th adiabatic states ${\displaystyle \psi _{j}[R(t_{n})]}$ computes the Pancharatnam product of overlaps:

$${\displaystyle I_{j}(\Gamma ,N)=\prod \limits _{n=0}^{N-1}\langle \psi _{j}[R(t_{n})]|\psi _{j}[R(t_{n+1})]\rangle .}$$

In the limit ${\displaystyle N\to \infty }$ one has (see Ryb & Baer 2004 for explanation and some applications)

$${\displaystyle I_{j}(\Gamma ,N)\to e^{i\beta _{j}}.}$$

Geometric phase and quantization of cyclotron motion

An electron subjected to magnetic field ${\displaystyle B}$ moves on a circular (cyclotron) orbit. Classically, any cyclotron radius ${\displaystyle R_{c}}$ is acceptable. Quantum-mechanically, only discrete energy levels (Landau levels) are allowed, and since ${\displaystyle R_{c}}$ is related to electron's energy, this corresponds to quantized values of ${\displaystyle R_{c}}$. The energy quantization condition obtained by solving Schrödinger's equation reads, for example, ${\displaystyle E=(n+\alpha )\hbar \omega _{c},}$ ${\displaystyle \alpha =1/2}$ for free electrons (in vacuum) or ${\textstyle E=v{\sqrt {2(n+\alpha )eB\hbar }},\quad \alpha =0}$ for electrons in graphene, where ${\displaystyle n=0,1,2,\ldots }$. Although the derivation of these results is not difficult, there is an alternative way of deriving them, which offers in some respect better physical insight into the Landau level quantization. This alternative way is based on the semiclassical Bohr–Sommerfeld quantization condition

$${\displaystyle \hbar \oint d\mathbf {r} \cdot \mathbf {k} -e\oint d\mathbf {r} \cdot \mathbf {A} +\hbar \gamma =2\pi \hbar (n+1/2),}$$

which includes the geometric phase ${\displaystyle \gamma }$ picked up by the electron while it executes its (real-space) motion along the closed loop of the cyclotron orbit.^[14] For free electrons, ${\displaystyle \gamma =0,}$ while ${\displaystyle \gamma =\pi }$ for electrons in graphene. It turns out that the geometric phase is directly linked to ${\displaystyle \alpha =1/2}$ of free electrons and ${\displaystyle \alpha =0}$ of electrons in graphene.