# Kapitza's pendulum

Drawing showing how a Kapitza pendulum can be constructed: a motor rotates a crank at a high speed, the crank vibrates a lever arm up and down, which the pendulum is attached to with a pivot.

Kapitza's pendulum is a rigid pendulum in which the pivot point vibrates in a vertical direction, up and down. It is named after Russian Nobel laureate physicist Pyotr Kapitza, who in 1951 developed a theory which successfully explains some of its unusual properties.[1] The unique feature of the Kapitza pendulum is that the vibrating suspension can cause it to balance stably in an inverted position, with the bob above the suspension point. In the usual pendulum with a fixed suspension, the only stable equilibrium position is with the bob hanging below the suspension point; the inverted position is a point of unstable equilibrium, and the smallest perturbation moves the pendulum out of equilibrium. In nonlinear control theory the Kapitza pendulum is used as an example of a parametric oscillator that demonstrates the concept of "dynamic stabilization".

The pendulum was first described by A. Stephenson in 1908, who found that the upper vertical position of the pendulum might be stable when the driving frequency is fast[2] Yet until the 1950s there was no explanation for this highly unusual and counterintuitive phenomenon. Pyotr Kapitza was the first to analyze it in 1951.[1] He carried out a number of experimental studies and as well provided an analytical insight into the reasons of stability by splitting the motion into "fast" and "slow" variables and by introducing an effective potential. This innovative work created a new subject in physics, that is vibrational mechanics. Kapitza's method is used for description of periodic processes in atomic physics, plasma physics and cybernetical physics. The effective potential which describes the "slow" component of motion is described in "Mechanics" volume of the Landau's Course of Theoretical Physics.[3]

Another interesting feature of the Kapitza pendulum system is that the bottom equilibrium position, with the pendulum hanging down below the pivot, is no longer stable. Any tiny deviation from the vertical increases in amplitude with time.[4] Parametric resonance can also occur in this position, and chaotic regimes can be realized in the system when strange attractors are present in the Poincaré section .

## Notation

Kapitza's pendulum scheme

Denote the vertical axis as ${\displaystyle y}$ and the horizontal axis as ${\displaystyle x}$ so that the motion of pendulum happens in the (${\displaystyle x}$${\displaystyle y}$) plane. The following notation will be used

Denoting the angle between pendulum and downward direction as ${\displaystyle \varphi }$ the time dependence of the position of pendulum gets written as

${\displaystyle {\begin{matrix}\left\{{\begin{matrix}x&=&l\sin \varphi \\y&=&-l\cos \varphi -a\cos \nu t\end{matrix}}\right.\end{matrix}}}$

## Energy

The potential energy of the pendulum is due to gravity and is defined by of the vertical position as

${\displaystyle {\begin{matrix}E_{POT}=-mg(l\cos \varphi +a\cos \nu t)\;.\end{matrix}}}$

The kinetic energy in addition to the standard term ${\displaystyle E_{KIN}=ml^{2}{\dot {\varphi }}^{2}/2}$, describing velocity of a mathematical pendulum, there is a contribution due to vibrations of the suspension

${\displaystyle {\begin{matrix}E_{KIN}={\frac {ml^{2}}{2}}{\dot {\varphi }}^{2}+mal\nu ~\sin(\nu t)\sin(\varphi )~{\dot {\varphi }}+{\frac {ma^{2}\nu ^{2}}{2}}\sin ^{2}(\nu t)\;.\end{matrix}}}$

The total energy is given by the sum of the kinetic and potential energies ${\displaystyle E=E_{KIN}+E_{POT}}$ and the Lagrangian by their difference ${\displaystyle L=E_{KIN}-E_{POT}}$.

The total energy is conserved in a mathematical pendulum, so time ${\displaystyle t}$ dependence of the potential ${\displaystyle E_{POT}}$ and kinetic ${\displaystyle E_{KIN}}$ energies is symmetric with respect to the horizontal line. According to the virial theorem the mean kinetic and potential energies in harmonic oscillator are equal. This means that the line of symmetry corresponds to half of the total energy.

In the case of vibrating suspension the system is no longer is a closed one and the total energy is no longer conserved. The kinetic energy is more sensitive to vibration compared to the potential one. The potential energy ${\displaystyle E_{POT}=mgy}$ is bound from below and above ${\displaystyle -mg(l+a) while the kinetic energy is bound only from below ${\displaystyle E_{KIN}\geq 0}$. For high frequency of vibrations ${\displaystyle \nu }$ the kinetic energy can be large compared to the potential energy.

## Equations of motion

Motion of pendulum satisfies Euler–Lagrange equations. The dependence of the phase ${\displaystyle \varphi }$ of the pendulum on its position satisfies the equation:[5]

${\displaystyle {\begin{matrix}{\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {\varphi }}}}={\frac {\partial L}{\partial \varphi }}\;,\end{matrix}}}$

where the Lagrangian ${\displaystyle L}$ reads

${\displaystyle {\begin{matrix}L={\frac {ml^{2}}{2}}{\dot {\varphi }}^{2}+ml(g+a~\nu ^{2}\cos \nu t)\cos \varphi \;,\end{matrix}}}$

up to irrelevant total time derivative terms. The differential equation

${\displaystyle {\begin{matrix}{\ddot {\varphi }}=-(g+a~\nu ^{2}\cos \nu t){\frac {\sin \varphi }{l}}\;,\end{matrix}}}$

which describes the movement of the pendulum is nonlinear due to the ${\displaystyle \sin \varphi }$ factor. The presence of the nonlinear term might lead to chaotic motion and to appearance of strange attractors.

## Equilibrium positions

Model of Kapitza's pendulum is more general than the model of mathematical pendulum. The latter is reproduced in the limit ${\displaystyle a=0}$. Its phase portrait is a simple circle ${\displaystyle x^{2}+y^{2}=l^{2}=const}$. If the energy in the initial moment was larger than the maximum of the potential energy ${\displaystyle E>mgl}$ then the trajectory will be closed and cyclic. If the initial energy is smaller ${\displaystyle E then the pendulum will oscillate close to the only stably point ${\displaystyle (x,y)=(0,-l)}$.

When the suspension is vibrating with a small amplitude ${\displaystyle a\ll l}$ and with a frequency ${\displaystyle \nu \gg \omega _{0}}$ much higher than the proper frequency ${\displaystyle \omega _{0}}$, the angle ${\displaystyle \varphi }$ may be viewed as a superposition ${\displaystyle \varphi =\varphi _{0}+\xi }$ of a "slow" component ${\displaystyle \varphi _{0}}$ and a rapid oscillation ${\displaystyle \xi }$ with small amplitude due to the small but rapid vibrations of the suspension. Technically, we perform a perturbative expansion in the "coupling constants" ${\displaystyle (a/l),(\omega _{0}/\nu )\ll 1}$ while treating the ratio ${\displaystyle (a/l)(\nu /\omega _{0})}$ as fixed. The perturbative treatment becomes exact in the double scaling limit ${\displaystyle a\to 0,\nu \to \infty }$. More precisely, the rapid oscillation ${\displaystyle \xi }$ is defined as

${\displaystyle {\begin{matrix}\xi ={\frac {a}{l}}\sin \varphi _{0}~\cos \nu t\;.\end{matrix}}}$

The equation of motion for the "slow" component ${\displaystyle \varphi _{0}}$ becomes

${\displaystyle {\begin{array}{rcl}{\ddot {\varphi }}_{0}={\ddot {\varphi }}-{\ddot {\xi }}&=&-(g+a~\nu ^{2}\cos \nu t){\frac {\sin \varphi }{l}}-{\frac {a}{l}}\left({\ddot {\varphi }}_{0}\cos \varphi _{0}~\cos \nu t-{\dot {\varphi }}_{0}^{2}\sin \varphi _{0}~\cos \nu t-2\nu {\dot {\varphi }}_{0}\cos \varphi _{0}~\sin \nu t-\nu ^{2}\sin \varphi _{0}~\cos \nu t\right)\\&=&-{\frac {g}{l}}\sin \varphi _{0}-(g+a~\nu ^{2}\cos \nu t){\frac {1}{l}}\left(\xi \cos \varphi _{0}+O(\xi ^{2})\right)-{\frac {a}{l}}\left({\ddot {\varphi }}_{0}\cos \varphi _{0}~\cos \nu t-{\dot {\varphi }}_{0}^{2}\sin \varphi _{0}~\cos \nu t-2\nu {\dot {\varphi }}_{0}\cos \varphi _{0}~\sin \nu t\right)\;.\end{array}}}$

Time averaging over the rapid ${\displaystyle \nu }$-oscillation yields to leading order

${\displaystyle {\begin{matrix}{\ddot {\varphi }}_{0}=-{\frac {g}{l}}\sin \varphi _{0}-{\frac {1}{2}}({\frac {a\nu }{l}})^{2}\sin \varphi _{0}\cos \varphi _{0}\;.\end{matrix}}}$

The "slow" equation of motion becomes

${\displaystyle {\begin{matrix}ml^{2}{\ddot {\varphi }}_{0}=-{\frac {\partial V_{\mathrm {eff} }}{\partial \varphi _{0}}}\;,\end{matrix}}}$

by introducing an effective potential

${\displaystyle {\begin{matrix}V_{\mathrm {eff} }=-mgl\cos \varphi _{0}+m({\frac {a\nu }{2}}\sin \varphi _{0})^{2}\;.\end{matrix}}}$

It turns out[1] that the effective potential ${\displaystyle V_{\mathrm {eff} }}$ has two minima if ${\displaystyle (a\nu )^{2}>2gl}$, or equivalently, ${\displaystyle (a/l)(\nu /\omega _{0})>{\sqrt {2}}}$. The first minimum is in the same position ${\displaystyle (x,y)=(0,-l)}$ as the mathematical pendulum and the other minimum is in the upper vertical position ${\displaystyle (x,y)=(0,l)}$. As a result the upper vertical position, which is unstable in a mathematical pendulum, can become stable in Kapitza's pendulum.

## Rotating solutions

The rotating solutions of the Kapitza's pendulum occur when the pendulum rotates around the pivot point at the same frequency that the pivot point is driven. There are two rotating solutions, one for a rotation in each direction. We shift to the rotating reference frame using ${\displaystyle \varphi \rightarrow \varphi ^{\prime }\pm \nu t}$ and the equation for ${\displaystyle \varphi }$ becomes:

${\displaystyle {\begin{matrix}{\ddot {\varphi }}^{\prime }=-{\frac {1}{l}}\left[{\frac {1}{2}}a\nu ^{2}\sin(\varphi ^{\prime })+g\sin(\varphi ^{\prime }\pm \nu t)+{\frac {1}{2}}a\nu ^{2}\sin(\varphi ^{\prime }\pm 2\nu t)\right]\;.\end{matrix}}}$

Again considering the limit in which ${\displaystyle \nu }$ is much higher than the proper frequency ${\displaystyle \omega _{0}}$, we find that the rapid-${\displaystyle \nu }$ slow-${\displaystyle \varphi _{0}^{\prime }}$ limit leads to the equation:

${\displaystyle {\begin{matrix}{\ddot {\varphi }}_{0}^{\prime }=-{\frac {1}{2l}}a\nu ^{2}\sin \varphi _{0}^{\prime }\;.\end{matrix}}}$

The effective potential is just that of a simple pendulum equation. There is a stable equilibrium at ${\displaystyle \varphi _{0}^{\prime }=0}$ and an unstable equilibrium at ${\displaystyle \varphi _{0}^{\prime }=\pi }$.

## Phase portrait

Interesting phase portraits might be obtained in regimes which are not accessible within analytic descriptions, for example in the case of large amplitude of the suspension ${\displaystyle a\approx l}$.[6][7] Increasing the amplitude of driving oscillations to half of the pendulum length ${\displaystyle a=l/2}$ leads to the phase portrait shown in the figure.

Further increase of the amplitude to ${\displaystyle a\approx l}$) leads to full filling of the internal points of the phase space, if before some points of the phase space were not accessible, now system can reach any of the internal points. This situation holds also for larger values of ${\displaystyle a}$.

## Interesting facts

• Kapitza noted that a pendulum clock with a vibrating pendulum suspension always goes faster than a clock with a fixed suspension.
• Walking is defined by an 'inverted pendulum' gait in which the body vaults over the stiff limb or limbs with each step. Increased stability during walking might be related to stability of Kapitza's pendulum. This applies regardless of the number of limbs - even arthropods with six, eight or more limbs.

## Literature

1. P. L. Kapitza, "Dynamic stability of a pendulum when its point of suspension vibrates", Soviet Phys. JETP 21, 588–592 (1951); P. L. Kapitza, "Pendulum with a vibrating suspension" Usp. Fiz. Nauk, 44, 7-15 (1951).
2. A. Stephenson "On an induced stability" Phil. Mag. 15, 233 (1908)
3. {{#invoke:citation/CS1|citation |CitationClass=book }}
4. Бутиков Е. И. «Маятник с осциллирующим подвесом (к 60-летию маятника Капицы»), учебное пособие.
5. {{#invoke:citation/CS1|citation |CitationClass=book }}
6. G.E. Astrakharchik, N.A. Astrakharchik «Numerical study of Kapitza pendulum» arXiv:1103.5981 (2011)
7. Time motion of Kapitza's pendulum can be modeled in online java applets he sites http://www.myphysicslab.com/beta/Inverted-pendulum.html http://faculty.ifmo.ru/butikov/Nonlinear/index.html Arbitrary parameters of the system can be used and can be inserted manually.