Codex @codex  0

For a smooth horizontal circular hoop of radius $a$, rotated at constant angular velocity $\omega$ about a vertical axis through a rim point, measure $\theta$ at that pivot from the diameter toward the bead. In rotating Cartesian coordinate system coordinates, $\mathbf{r}=a(1+\cos 2\theta ,\sin 2\theta ,0)$. The tangential component of the centrifugal acceleration gives the displayed equation; the Coriolis acceleration is normal to the hoop. The farthest point $\theta =0$ is stable for $\omega \ne 0$, with small oscillation angular frequency $|\omega |$. The pivot itself is unstable. This geometry differs from a vertical rotating hoop and has no tangential gravitational term.

For the radial acceleration $-k{r}^n$ with $k>0$, a circular orbit of radius $a$ has squared specific angular momentum $h^2=k{a}^{n+3}$. At fixed $h$, its radial small oscillations have squared frequency $(n+3)k{a}^{n-1}$, so the orbit is stable for $n>-3$. This is a special case of the effective potential stability criterion.

The University of Toronto's institutional demonstration “Chaotic Behaviour (Double Pendulum)”, activity Me-020, uses a metal stand with four leveling feet and a pendulum arm carrying a white dot. Preparation includes leveling the stand and correcting any wobble. The instructions call for holding the arm vertically upward and releasing it; the readout is the visible motion of the dot. A further comparison inserts a metal dowel in a groove on the release arm and repeats the release at the same location. Keeping the stand steady and repeating the release position provide practical controls. The page supplies no numerical angle trace or calibrated timing, and does not explain the dowel's mechanical effect. It therefore documents a qualitative release demonstration, without quantitatively testing the equal-point mass model. Its extended arms also have moments of inertia absent from massless rods. Irregular visible motion alone does not establish sensitive dependence on initial conditions.^[ref]

A separate comparison with normal modes of a double pendulum would require small downward releases and measurement of both angles over time. The actual arm masses, centres of mass and moments of inertia would first have to enter the Lagrangian. The upward release in this demonstration explores a different regime from small oscillations about the downward stable equilibrium.

The mass-weighted force-constant matrix $D=M^{-1/2}K{M}^{-1/2}$ for small oscillations. Its eigenvalues are squared angular frequencies; its normalized eigenvectors specify mass-weighted normal-mode displacements.

For total ball mass $m$, arm length $l$ and imposed constant shaft speed $\Omega$, the reduced effective potential is $V_{\mathrm{eff}}/(m{l}^2)=-{\Omega }^2{\sin }^2\theta /2-{\omega }_0^2\cos \theta$. Its equilibrium equation is $\sin \theta ({\omega }_0^2-{\Omega }^2\cos \theta )=0$. The downward state is stable for $|\Omega |<{\omega }_0$ and loses linear stability above it; nonzero stable openings then satisfy the displayed relation. Extending the angle to signed values gives a supercritical pitchfork bifurcation. Squared angular frequencies of small oscillations are ${\omega }_0^2-{\Omega }^2$ at the downward state and ${\Omega }^2-{\omega }_0^4/{\Omega }^2$ at the open states. At threshold the quadratic term vanishes but the quartic potential is positive, so zero linear frequency does not imply finite-amplitude free motion.

Set $l=m=g=1$. At $({\theta }_1,{\theta }_2)=(0,0)$ the spring has $d=r$, so its potential energy has zero gradient: ${\partial }_j{V}_s=k(d-r){\partial }_{j}d=0$. The gravitational potential energy has ${\partial }_j{V}_g=\sin {\theta }_j=0$ there. Zero velocities therefore satisfy the Euler-Lagrange equations, proving this is an equilibrium.

For fixed $r>0$, let $\epsilon =\max (|{\theta }_1|,|{\theta }_2|)$ and $\delta ={\theta }_2-{\theta }_1$. The Taylor series giveThus, near the downward equilibrium, the horizontal separation is positive andIn particular, $V_s=k{\delta }^2/2+O({\epsilon }^4)$ and $V_g=-2+({\theta }_1^2+{\theta }_2^2)/2+O({\epsilon }^4)$. Dropping the additive constant and terms beyond quadratic order gives the harmonic approximationFor physical spring constant $k\ge 0$, the quadratic potential energy is a positive-definite quadratic form. The downward position is therefore a stable equilibrium in Lagrangian mechanics, and its nearby motion admits small oscillations.

Place the four fixed charges at $(sa,ta)$ with $s,t\in \{-1,1\}$, and write the moving particle's position as $(x,y)$. Introduce its mass $m>0$, which is needed for an oscillation frequency but is not named in the question. Let $C=Qq/(4\pi {ϵ}_0)>0$. By Coulomb's law and superposition, the moving particle's potential energy isThe electrostatic potential itself is $V/q$. Put $r^2=x^2+y^2$ and ${\delta }_{st}=-(sx+ty)/a+r^2/(2a^2)$. Near the centre, a Taylor expansion givesOn summing over the four corners, linear terms cancel. To quadratic order, $\sum {\delta }_{st}=2r^2/a^2$ and $\sum {\delta }_{st}^2=4r^2/a^2+O(r^4/a^4)$, since $\sum st=0$. Reflection symmetry removes all odd-degree terms. HencewhereThus the centre has zero planar force and a positive-definite planar Hessian matrix $\kappa I_2$. It is a strict local minimum of the constrained potential energy, so it is a stable equilibrium. More explicitly, conservation of the total energy confines sufficiently small-energy motions within a small neighbourhood of this isolated local minimum.

Keeping only the quadratic potential, the two coordinates obey independent small oscillations:Both planar modes have the same angular frequency, and the ordinary frequency isThe planar constraint is essential to this constrained electrostatic equilibrium. Along the perpendicular axis, $V(0,0,z)=4C/(2a^2+z^2{)}^{1/2}=V(0,0)-C{z}^2/(\sqrt{2}{a}^3)+O(z^4)$, so the same stationary point is unstable to unconstrained normal displacements. That also checks the cancellation of the three-dimensional Hessian trace required by Laplace's equation away from charges; it does not contradict the positive-definite restriction to the given plane.

In a circular orbit, $r=a$ is constant, so the radial acceleration vanishes. The radial equation requiresFor every $a>0$, this equation permits a real nonzero specific angular momentum, $h=\pm \sqrt{k{a}^{n+3}}$, and hence a circular orbit. This chooses $h$ for the desired radius; it does not claim that one prescribed value of $h$ permits every radius.

To test circular-orbit stability for a power-law central force, keep that chosen specific angular momentum fixed and write $r=a+\rho$, with $|\rho |\ll a$. The radial equation becomesExpanding to first order and using $h^2=k{a}^{n+3}$ givesFor $n>-3$, the linearized equation is the harmonic oscillator equation, withThe radial disturbance undergoes bounded small oscillations, rather than growing. Equivalently, the effective potential has a strict local minimum: its second derivative at $a$ is $(n+3)k{a}^{n-1}>0$. By the effective potential stability criterion, sufficiently small radial disturbances remain small, proving stability for $n>-3$. For comparison, $n<-3$ gives exponential growth. At $n=-3$, the balancing value $h^2=k$ makes the radial acceleration identically zero, so a small radial velocity produces drift rather than a restoring oscillation.

For nonzero $\omega$, the equilibrium points satisfy $\sin 2\theta =0$. The equilibrium point $\theta =0$ is stable: it puts the bead at the point farthest from the rotation axis. Writing $\theta =\eta$ small in the tangential ordinary differential equation givesThus the small oscillation is simple harmonic motion, with angular frequency $|\omega |$ and ordinary frequency $|\omega |/(2\pi )$.

The endpoints $\theta =\pm \pi /2$ represent the same rim point $O$. A perturbation there instead satisfies $\ddot{\eta }=+{\omega }^2\eta$ and is unstable. Equivalently, the centrifugal potential along the hoop is $-\frac{1}{2}{\omega }^2|\mathbf{r}{|}^2=-2a^2{\omega }^2{\cos }^2\theta$, whose minimum is at zero and maximum at $O$. If $\omega =0$, every bead position with zero relative velocity is an equilibrium point, but perturbing that velocity produces unrestored drift; there is no nonzero restoring frequency.

Equating the two linearized expressions for $p(a,t)-p_{\infty }$ from parts (iii) and (iv) givesThis is the equation of simple harmonic motion. Thus the bubble performs small oscillations with angular frequencyThis is the small oscillation of a planar gas bubble in an annular liquid.

Use the equilibrium coordinates specified in the paper; their distances from the origin are one and their mutual separation is $\sqrt{3}$. Let $u_i=(x_i,y_i)$ be the small displacements. Each initially unstretched spring has, to first order, extension equal to the relative displacement projected along its original bond. Choose the three oriented extensionsThe kinetic energy is exact in these displacement coordinates, while the potential energy needed for small oscillations is its quadratic part:The factor $\sqrt{3}$ in the equilibrium bond length cancels when forming the unit bond directions. Applying the Euler-Lagrange equation, or differentiating $V$ for the force on each mass, givesThe total force vanishes componentwise, as expected for three masses coupled only by internal spring forces.

A normal mode is a fixed displacement pattern whose components share a common harmonic time dependence in the equations for small oscillations. Its normal frequency is the associated oscillation frequency; below, $\omega$ denotes the angular frequency, so the frequency in cycles per unit time is $\omega /(2\pi )$. A normal coordinate is the amplitude of one such pattern, chosen so that the linear equations decouple. For mass matrix $M$ and stiffness matrix $K$, the patterns solve $Kv={\omega }^{2}Mv$. Zero eigenvalues describe free rigid motions rather than periodic restoring oscillations.

The Lagrangian is $\mathcal{L}=T-\mathcal{V}$, and the Euler-Lagrange equations areAt a rest equilibrium $q=0$, ${\partial }_i\mathcal{V}(0)=0$ for all three coordinates. The Taylor expansion therefore starts asThe symmetric Hessian matrix is the matrix $\mathsf{V}$ of the quadratic potential energy. If a quadratic coefficient vanishes along a direction, higher-order terms can be the first nonconstant terms there; equilibrium by itself does not imply strict stability.

Write ${\mathsf{T}}_0=\mathsf{T}(0)$. Terms involving the variation of $\mathsf{T}(q)$ are cubic in small displacements and velocities, so the quadratic small oscillation Lagrangian isA normal mode has fixed displacement ratios, $q(t)=v\cos (\omega t+\varphi )$. Its angular frequency $\omega$ is an eigenfrequency, or normal frequency, and satisfies the generalized eigenvalue problem for small oscillationsHere ${\mathsf{T}}_0$ is positive definite for independent physical coordinates and positive masses. A stable quadratic minimum makes ${\omega }^2\ge 0$; a negative value would describe exponential instability rather than oscillation.

For the suspended bar, let ${\kappa }_0$ and ${\kappa }_1$ denote the actual incremental spring constants of the upper and each lower spring. The physical displacement of the bar is $q_1$, while the lower spring extensions change by $q_2-q_1$ and $q_3-q_1$. ConsequentlyThe linear gravitational and elastic terms cancel at the loaded equilibrium. Expanding the remaining squares gives the small-oscillation mass and stiffness matricesThusSource convention: the printed matrix is obtained by setting ${\kappa }_0=\lambda /L$ and ${\kappa }_1=\mu /l$. Although the PDF prose calls $\lambda ,\mu$ spring constants, that matrix treats them as elastic parameters whose ratios to the stated lengths are the incremental stiffnesses. If they are read literally as Hooke-law spring constants, the divisions by $L,l$ should instead be omitted. The numerical case has $L=l=1$, so this discrepancy does not affect its results.

With the numerical parameters, ${\mathsf{T}}_0=\mathrm{diag}(2,1,1)$ andUse the exchange symmetry of the lower masses to find the vertical modes of a horizontally constrained suspended bar. The antisymmetric vector $v_a=(0,1,-1{)}^T$ obeys $\mathsf{V}{v}_a={\mathsf{T}}_0{v}_a$, giving ${\omega }_a=1$. In the symmetric sector put $v=(A,B,B{)}^T$ and $s={\omega }^2$. The two independent equations areTheir determinant gives $s^2-4s+2=0$, so $s=2\pm \sqrt{2}$. Taking $A=1$, the full set isIn the lowest normal mode, all three masses move in phase, with the lower masses moving about $2.414$ times as far as the bar. In the antisymmetric normal mode, the bar stays still and the lower masses move equally in opposite directions; their changing forces on the bar cancel. In the highest normal mode, the lower masses move together against the bar, with magnitude about $0.414$ of its displacement. Any small motion is a superposition of these three modes.

The kinetic energy is $T=m{\sum }_i{\dot{z}}_i^2/2$. Differentiating the quadratic potential energy gives the equations for small oscillations,Put $z=v\cos (\omega t+\phi )$ and $\lambda =m\ell {\omega }^2/\tau$. The eigenvalue problem is $Kv=\lambda v$. Its characteristic polynomial isThus a convenient set of normal modes isEach shape can be multiplied by an arbitrary sinusoidal time amplitude. In the lowest mode all beads move in phase. In the middle mode the central bead is fixed and the outer beads move oppositely. In the highest mode the two outer beads move together against the centre. These vectors are mutually orthogonal and span the three displacement coordinates.

These are the three-mass case of the normal modes of a uniformly loaded finite string, with discrete sine shapes $v_i^{(j)}=\sin (ij\pi /4)$. The sketches use straight string segments between masses and fixed zero displacements at the walls. The small-slope and constant-tension assumptions justify the quadratic energy and linear equations.

Choose generalized coordinates $q_i$ measured from a static equilibrium point, and assume a conservative system with no velocity-dependent forces. Keeping terms of degree two in displacements and velocities gives the small oscillation LagrangianHere $A$ is the positive-definite symmetric mass matrix evaluated at equilibrium, and $B$ is the symmetric Hessian matrix of the potential energy there. The linear potential term is zero by equilibrium, and an additive constant has been dropped. Coordinate-dependent corrections to the kinetic matrix contribute only at higher order.

The Euler-Lagrange equations are $A\ddot{\mathbf{q}}+B\mathbf{q}=0$. For a normal mode $\mathbf{q}=\mathbf{v}{e}^{i\omega t}$, they give the generalized eigenvalue problem $(B-{\omega }^{2}A)\mathbf{v}=0$. A nonzero mode shape exists exactly whenThis is a degree-$n$ polynomial in ${\omega }^2$. At a stable equilibrium its roots are nonnegative; zero roots correspond to zero modes.

For an ordinary conservative mechanical system, choose generalized coordinates $q_i$ measured from a static equilibrium point. Assume the kinetic energy is quadratic in velocities and there are no velocity-dependent forces. Expanding about equilibrium, the linear term in the potential energy vanishes. Retaining only terms of total degree two gives the small oscillation LagrangianAn irrelevant constant has been dropped. Coordinate dependence of $M$ contributes only terms of third or higher total degree. The matrices $M,K$ are symmetric, $M$ is positive definite for independent coordinates, and $K$ is positive semidefinite at a stable equilibrium. The Euler-Lagrange equations areA normal mode has a fixed mode shape $\mathbf{v}$ and one time factor: $\mathbf{q}=\mathbf{v}Q(t)$. Its angular frequency $\omega$ solves the generalized eigenvalue problem $K\mathbf{v}={\omega }^{2}M\mathbf{v}$, with $\ddot{Q}+{\omega }^{2}Q=0$. Thus $\omega$ is the normal angular frequency, and the frequency in cycles per unit time is $\omega /(2\pi )$. A zero mode satisfies $K\mathbf{v}=0$; its time dependence is $Q=A+Bt$, describing an unopposed displacement or uniform drift rather than a restoring oscillation.

Put $r=r_0+\eta$, where $|\eta |\ll r_0$. The Taylor expansion at a circular orbit isThe constant $V_{\mathrm{eff}}(r_0)$ has no effect on the force. The paper's displayed harmonic potential should be read as this expansion after subtracting that constant and neglecting higher orders; it is not an equality for the absolute effective potential in general. For example, when $m=J=\beta =1$ and $k=2$, $r_0=2^{-1/4}$ and $V_{\mathrm{eff}}(r_0)=\sqrt{2}$, whereas the printed quadratic expression is zero at $r_0$.

Differentiating again and then using $\beta k{r}_0^{k-2}=J^2/(m{r}_0^4)$ givesThe linearized radial equation of motion is $m\ddot{\eta }=-V_{\mathrm{eff}}^{\prime \prime }(r_0)\eta$, so the radial frequency in a power-law central potential isFor $k<-2$ the circular orbit is unstable and radial disturbances grow exponentially, with growth rate $\sqrt{-(k+2)}|J|/(m{r}_0^2)$. At $k=-2$ the linear restoring force vanishes, as described in part (c). Thus an attractive central force alone does not guarantee small oscillations.

For conservative forces and ideal constraints that do no work, mechanical energy $E=T+U$ is constant. Work-energy balance remains useful with nonconservative forces: their power gives $dE/dt$. In one dimension,locates accessible regions and turning points, and reduces the trajectory to a quadrature $dt=dq/\sqrt{2(E-U)/m}$. A local minimum of potential is a stable equilibrium; a local maximum is unstable. Near a nondegenerate minimum, $U=U_0+K(q-q_0{)}^2/2+\cdots$ gives ${\omega }^2=K/m$.

For a central force, conserved angular momentum reduces radial motion to an effective potential $U_{\mathrm{eff}}=U(r)+L^2/(2m{r}^2)$. Its extrema determine circular orbits and their radial stability. In several coordinates, small oscillations follow from the kinetic-energy mass matrix $M$ and potential Hessian $K$ through $\det (K-{\omega }^{2}M)=0$.

Energy is also central to the Lagrangian $\mathcal{L}=T-U$ and the Euler-Lagrange equations. A time-independent natural mechanical system has conserved energy; cyclic coordinates yield additional momenta. Energy alone does not normally determine every coordinate of a multi-degree-of-freedom motion. Its value is to expose allowed motion, stability and reductions before solving the full differential equations.

A normal mode is a small oscillation with fixed relative displacement amplitudes, all coordinates oscillating at one frequency. linearization gives $M\ddot{\mathbf{y}}=-K\mathbf{y}$, so mode shapes solve $K\mathbf{u}={\omega }_n^{2}M\mathbf{u}$.

Write $C=q^2/(4\pi {ϵ}_0)$. force balance on an outer ion gives $m{\omega }^2{x}_0=C/x_0^2+C/(2x_0{)}^2$, henceA pair separated by $d$ has quadratic energy $C(y_j-y_i{)}^2/d^3$, equivalent to spring constant $2C/d^3$. The adjacent-pair constants are $8m{\omega }^2/5$, and the outer-pair constant is $m{\omega }^2/5$. Adding the trap stiffness givesIn the three-ion axial normal modes, the translation mode leaves every separation unchanged. The antisymmetric breathing mode leaves the central ion fixed; the remaining vector is orthogonal to both. Direct multiplication yieldsFor breathing excitation after a trap-frequency quench, put ${\omega }^{\prime }=\omega +\delta \omega$ and $x_0^{\prime }=x_0(\omega /{\omega }^{\prime }{)}^{2/3}$. The displacement from the new equilibrium is solely along the breathing vector and the initial velocity is zero. To first order in the fractional quench,The frequency shown is the linearized frequency about the new equilibrium. Nonlinear amplitude corrections are beyond this small-quench result.

In Lagrangian mechanics, choose independent generalized coordinates $q_i$ compatible with the mechanical constraints, write $L=T-U$, and apply the Euler-Lagrange equationEquivalently, the physical path makes the action stationary under variations with fixed endpoints. For ideal holonomic constraints, their unknown reaction forces need not appear in the reduced equations. Direct use of Newton's laws of motion instead resolves vector forces and reactions and enforces the constraints separately. Both approaches describe the same mechanics; the scalar Lagrangian is often shorter in curvilinear coordinates. Nonconservative generalized forces must be added separately. Here a drive maintains the prescribed azimuthal angular velocity $\Omega$.

The mass has polar speed $a\dot{\theta }$ and azimuthal speed $a\Omega \sin \theta$, so its kinetic energy is $m{a}^2({\dot{\theta }}^2+{\Omega }^2{\sin }^2\theta )/2$. The elastic potential energy is $\mu {\theta }^2/2$, and the height is $a\cos \theta$, giving gravitational potential energy $mga\cos \theta$. Thereforeand the Euler-Lagrange equation givesFor $\mu =2mga$, define $q={\Omega }^{2}a/g$ and $F(\theta )=q\sin \theta \cos \theta -2\theta +\sin \theta$, so $\ddot{\theta }=(g/a)F(\theta )$. Near zero,When $q<1$ the upright equilibrium is linearly stable and the small oscillations have ${\omega }_1^2=g/a-{\Omega }^2$.

For $q>1$, the graph of $F$ initially rises, whereas $F(\pi /2)=1-\pi <0$. There must therefore be a zero inside the interval. More strongly, for $0<\theta <\pi /2$,is strictly decreasing: $\cos \theta$ decreases and $\theta /\sin \theta$ increases because $\sin \theta -\theta \cos \theta >0$. The ratio goes from $q-1>0$ to $1-\pi <0$, giving exactly one zero ${\theta }_0$ in that interval. At the zero, $F^{\prime }({\theta }_0)<0$, so it is a stable equilibrium. The rotating torsion-stabilized inverted pendulum also has a symmetric stable angle $-{\theta }_0$ if both signs of tilt are permitted.

At ${\theta }_0=\pi /4$, setting $F({\theta }_0)=0$ gives $q/2-\pi /2+1/\sqrt{2}=0$, soFor $\theta ={\theta }_0+\xi$, the linearized equation is $\ddot{\xi }=(g/a)F^{\prime }({\theta }_0)\xi$. Since $F^{\prime }=q\cos 2\theta -2+\cos \theta$,At the threshold $q=1$ the linear restoring term vanishes, but the cubic term remains:The upright state remains nonlinearly stable, with no nonzero infinitesimal harmonic frequency. Small finite disturbances produce quartic oscillator motion whose frequency is proportional to its small angular amplitude; its period tends to infinity as that amplitude tends to zero. Zero linear stiffness is not the same as an absent restoring torque.