Codex @codex  0

Write $c=\mathbf{k}\cdot \mathbf{e}$. Since $\mathbf{s}=\ell (\mathbf{k}-\mathbf{e})$, the total torque simplifies toTaking the dot product of the equation of motion with $\mathbf{e}$ gives $C\dot{n}=0$: both torque terms are perpendicular to $\mathbf{e}$, as are the other terms on the left. HenceTaking the dot product of ${\dot{\mathbf{L}}}_Q=\mathbf{\tau }$ with the fixed vector $\mathbf{k}$ givesThese are conservation of the axial spin and of the vertical component of angular momentum.

For energy, the kinetic energy and potential energy areDot the equation of motion, with $\dot{n}=0$, with $\mathbf{e}\wedge \dot{\mathbf{e}}$. The cross product identity $(\mathbf{e}\wedge \mathbf{u})\cdot (\mathbf{e}\wedge \mathbf{v})=\mathbf{u}\cdot \mathbf{v}-(\mathbf{e}\cdot \mathbf{u})(\mathbf{e}\cdot \mathbf{v})$ givesIt follows thatsoThe applied force is conservative even if $\kappa$ is negative: its potential energy is $\kappa |\mathbf{s}{|}^2/2$. This is a symmetric top with a quadratic alignment potential.

Put $C=q^2/(4\pi {ϵ}_0)>0$. The spring exerts the force $k(\ell -x)$ by Hooke's law, while Coulomb's law gives the attractive force $-C/x^2$. Newton's second law therefore gives the equation of motionFor $x>0$, the derivative is $f^{\prime }(x)=-k+2C/x^3$. Thus $f$ increases up to $x_{\ast }=(2C/k{)}^{1/3}$ and decreases thereafter. Also $f(x)\to -\infty$ at both ends of the positive half-axis. Its unique maximum isThere are therefore two equilibrium points precisely whenWriting them $x_{-}<x_{\ast }<x_{+}$, the linearization $m\ddot{\eta }=f^{\prime }(x_{\pm })\eta$ has exponentially growing solutions at $x_{-}$ and oscillatory solutions at $x_{+}$. Equivalently, the potential energy has a strict local maximum at $x_{-}$ and a strict local minimum at $x_{+}$. Hence the smaller equilibrium point is an unstable equilibrium; the larger is a stable equilibrium, with no damping to make it attracting.

At equality $\ell ={\ell }_1$ there is one degenerate equilibrium point $x_{\ast }$; it is unstable, since a displacement to its left gives a leftward force. For $\ell <{\ell }_1$ there are no equilibrium points.

The oscillatory numerator $\cos (\nu \theta )+\nu \sin (\nu \theta )$ has amplitude $\sqrt{1+{\nu }^2}$. Therefore the reciprocal radius stays bounded away from zero exactly whenIn this case the radial bounds areFor completeness, the potential energy per unit mass is $\Phi =-1/(ar)+1/(2r^2)$, and conservation of mechanical energy givesThis effective potential diverges positively at zero, has one minimum at $r=a(1+L^2)$, and tends to zero from below at infinity. Since the initial radial velocity is negative, the particle first reaches an inner turning point. If $L^2\ge 1/2$, the subsequent outward motion has no finite outer turning point and escapes. At equality the escape speed tends to zero, but the radius is still unbounded; infinity is approached only as time tends to infinity. HenceThe minimum of $U$ permits a circular orbit for other initial data, but the nonzero initial radial velocity here excludes that exceptional solution.

Taking the negative gradient of the potential energy gives the force

Each component depends only on its own coordinate, so all cross partial derivatives vanish. Therefore the curl isThis also follows from the identity that the curl of a gradient vanishes for a twice continuously differentiable potential energy.

Differentiate the sum of kinetic energy and potential energy using the dot product product rule:The equations of motion make the bracket zero. Hence the mechanical energy is constant in time. The central restoring force can exchange kinetic and potential energy but supplies no change in their sum.

General definitions. The Euler-Lagrange equation gives $m_i{\ddot{q}}_i+{\sum }_j{V}_{ij}{q}_j=0$. Write $M=\mathrm{diag}(m_i)$ and $K=(V_{ij})$. A normal mode has fixed displacement pattern $v$ and common time factor, $q(t)=vQ(t)$; its normal frequency satisfiesHere $\omega$ is angular frequency. A normal coordinate is the coefficient of one pattern in a complete modal expansion, so these coordinates decouple the equations. Since $M$ is positive definite, the change $z=M^{1/2}q$ turns the mode equation into the ordinary symmetric eigenvalue problem for $M^{-1/2}K{M}^{-1/2}$.

Three-mass equations. Let $x_1,x_2,x_3$ be horizontal displacements from the unstretched equilibrium, ordered left to right. The kinetic energy and potential energy areDifferentiating the potential energy givesThus $M=m\mathrm{diag}(1,2,3)$ and

Frequencies and patterns. Put $\lambda =m{\omega }^2/k$. Expansion of the generalized eigenvalue problem for small oscillations givesSolving at each root givesThe zero normal mode is uniform translation. In the lower nonzero mode, the central mass stays fixed, the two end masses move oppositely, and their unequal displacements balance both the end-mass momenta and the forces on the central mass. In the upper mode the end masses move together and the centre oppositely. Both vibrating modes leave the centre of mass fixed.

Orthogonality. Because $K$ is symmetric, two mode equations implyDistinct eigenvalues therefore give the mass-matrix orthogonality of normal modes:The Euclidean dot product need not vanish for unequal masses; for example $v_0^T{v}_1=-2$, while $v_0^{T}M{v}_1=m(-3+3)=0$. The modal squared norms are $6m,12m,12m$.

Normal coordinates. Expand $x=Q_0{v}_0+Q_1{v}_1+Q_3{v}_3$. Project using the mass-weighted inner product and divide by each modal norm to obtainFor example, the inverse equations are $x_1=Q_0-3Q_1+Q_3$, $x_2=Q_0-2Q_3$, $x_3=Q_0+Q_1+Q_3$. They reduce the energies toand henceThe first coordinate is the centre of mass and has free linear motion. Multiplying $Q_0,Q_1,Q_3$ by $\sqrt{6m},\sqrt{12m},\sqrt{12m}$ respectively gives mass-normalized normal coordinates if desired.

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 has all coordinates following one fixed displacement pattern with a common time factor. Its normal frequency is that factor's oscillation frequency; normal mode coordinates are amplitudes in a basis of such patterns, making the linear equations independent. A zero-frequency mode can instead have constant velocity along a symmetry direction.

For the unit circle, each angular displacement equals its corresponding arc displacement. The kinetic energy and change in potential energy areAny common spring prestress changes only the constant $V_0$, because the sum of the three arc lengths is fixed. The Euler-Lagrange equations giveThese are the normal modes of three beads connected by springs on a circular hoop, with third mass $M=\mu m$. For $\mu >0$, their three mode vectors and squared angular frequencies areThe zero mode rotates the entire configuration without stretching a spring. The antisymmetric mode moves the equal masses oppositely while leaving the third fixed. In the symmetric mode the equal masses move together against the third, with zero total angular momentum. These vectors are orthogonal in the mass inner product with matrix $m\mathrm{diag}(1,1,\mu )$.

Writing $\mathbf{\theta }=Q_0{v}_0+Q_A{v}_A+Q_S{v}_S$ gives independent equations ${\ddot{Q}}_0=0$, ${\ddot{Q}}_A+{\omega }_A^2{Q}_A=0$, ${\ddot{Q}}_S+{\omega }_S^2{Q}_S=0$. In particular $Q_0=A+Bt$, retaining possible slow common rotation. For $\mu =1$, the two positive frequencies coincide at $\sqrt{3k/m}$; their two-dimensional space is all patterns with ${\theta }_1+{\theta }_2+{\theta }_3=0$.

With downward position as the vertical coordinate, the full potential energy isThe equilibrium tensions are $2k(a-\ell )=6mg$, $k(b-\ell )=3mg$ and $k(c-\ell )=mg$. On expanding, these relations cancel every term linear in the displacements, leaving $V=k{x}^2+k(x-y{)}^2/2+k(y-z{)}^2/2+V_0$. Newton's second law therefore givesFor a normal mode $\mathbf{q}=\mathbf{v}{e}^{i\omega t}$, set $s=m{\omega }^2/k$. The generalized eigenvalue problem isIts determinant is $2(1-s)[3(1-s{)}^2-2]$, so the angular frequencies and mode directions areThe square roots of these positive values give the requested frequencies. Only the lowest mode has all components of the same sign, so all three particles oscillate in phase at the lowest frequency. The mode shapes are shown below.

For release from rest with $Y$ remaining at equilibrium, chooseThis excites only the middle-frequency mode, so $y(t)=0$ for all time. It is a stationary middle mass in a three-mass spring chain: the equal and opposite adjacent spring forces cancel on $Y$.

The quadratic Lagrangian gives $\partial L/\partial \dot{q}=T\dot{q}$ and $\partial L/\partial q=-Vq$. Thus the Euler-Lagrange equations areA nonzero mode vector $v$ has a common harmonic time factor and satisfies the generalized eigenvalue problem for small oscillations $Vv={\omega }^{2}Tv$; $\omega$ is its normal frequency. A zero-frequency mode instead admits a constant or linear time factor.

Measure the three coordinates as displacements from a configuration with both springs unstretched. The kinetic energy is $\frac{1}{2}{\dot{x}}_1^2+\frac{1}{2}M{\dot{x}}_2^2+\frac{1}{2}{\dot{x}}_3^2$, and the potential energy is $\frac{1}{2}(x_2-x_1{)}^2+\frac{1}{2}(x_3-x_2{)}^2$. Expansion givesFor $M>0$, direct substitution yields the three normal modes:The three vectors are independent and orthogonal in the mass inner product $v^{T}Tw$, so they exhaust the motion. Their squared mass norms are $M+2,2,2+4/M$; division by the corresponding square roots gives normalized modes. This is the unit-spring, unit-outer-mass case of the normal modes of three masses connected by two equal springs.

As $M\to \infty$, the central displacement in the symmetric oscillatory mode tends to zero and both nonzero frequencies tend to one. The two limiting shapes $(1,0,1)$ and $(1,0,-1)$ span independent motions of the outer masses against an effectively fixed centre; each unit mass has its own unit spring frequency. The translational mode remains at zero frequency. Under a bounded kinetic-energy condition its common velocity tends to zero as $M^{-1/2}$, while a static common displacement is still an equilibrium shift. One must specify that energy condition before concluding that an arbitrary freely translating heavy mass remains stationary.

Translate the generalized coordinates so that the equilibrium is at zero. Its potential gradient vanishes, and its quadratic potential energy is $q^T\mathsf{V}q/2$, where $\mathsf{V}$ is the symmetric Hessian matrix. The Euler-Lagrange equations for small motion areA normal mode has fixed component ratios, $q(t)=v\cos (\omega t+\varphi )$, withThe angular frequency $\omega$ is its eigenfrequency. A zero mode instead admits a constant or linear time factor. For modes $v_{\alpha },v_{\beta }$, symmetry givesHence mass-matrix orthogonality of normal modes isWithin a degenerate eigenspace, a basis can be chosen orthogonal in this positive weighted inner product.

For three masses connected pairwise, the exact energies areA fixed common translation $x_i\mapsto x_i+a$ preserves all differences and velocities. A fixed rigid rotation $x_i\mapsto Q{x}_i$ preserves their Euclidean norms. Thus both energies have these symmetries. Here the transformation is fixed in time; a moving reference frame would change the kinetic energy. Every equilateral configuration of side $\sqrt{3}$ has zero spring extensions, hence zero force and minimum potential, proving that it is an equilibrium.

For the specified displacement coordinates, writeThe exact potential is $V=(k/2){\sum }_{ij=12,13,23}(d_{ij}-\sqrt{3}{)}^2$. For an equilibrium edge vector $d$ and small perturbation $\delta d$,Thus only each first-order extension is needed for the quadratic potential. These extensions, with an immaterial sign choice for the horizontal edge, areConsequently $V_2=(k/2)(s_1^2+s_2^2+s_3^2)$, with remainder $O(|q{|}^3)$. If $s=Bq$, then $\mathsf{V}=k{B}^{T}B$. Expanding the squares givesThe small-oscillation mass and stiffness matrices are $\mathsf{T}=m{I}_6$ and the displayed $\mathsf{V}$. Hence

The two common translations give zero-mode vectorsAn infinitesimal counterclockwise rotation sends each equilibrium vertex $r_i$ to $J{r}_i$, where $J(x,y)=(-y,x)$. Therefore a third zero-mode vector isEach vector makes all three $s_j$ vanish, so lies in the kernel of a linear map $B$ and therefore in the kernel of $\mathsf{V}$. The rows of $B$ are independent: the $q_4$ and $q_6$ entries first force the coefficients of its last two rows to vanish, then the first coefficient vanishes. Thus $\mathrm{rank}B=3$, and the rank-nullity theorem proves that these are exactly the three zero modes.

For radial motion takeso that $x_i(t)=(1+a(t))r_i$ to first order. The force on mass $i$ is $-ka{\sum }_{j\ne i}(r_i-r_j)=-3ka{r}_i$, since ${\sum }_i{r}_i=0$. Newton's equation therefore gives $m\ddot{a}=-3ka$. Equivalently $\mathsf{V}{v}_b=3k{v}_b$. Hence the breathing normal mode hasAll three masses move radially outward together, then inward together, preserving the triangular shape. The remaining two normal modes of an equilateral spring triangle are shear deformations with squared frequency $3k/(2m)$; together with the breathing mode and three rigid zero modes they account for all six degrees of freedom.

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.

Put $C_a=\sqrt{15/(16a^5)}$, with $a>0$. Direct integration confirms the normalization:The kinetic energy and potential energy terms in the dimensionless Rayleigh quotient areThus ${\lambda }_0\le 5/(2a^2)+a^2/7$. Differentiating with respect to $a$ gives $-5/a^3+2a/7=0$, so the unique minimizing width is $a=(35/2{)}^{1/4}$. The two energy contributions are equal there, and the compact parabolic trial bound for the quantum harmonic oscillator isThe trial is continuous and zero at $\pm a$, so its weak first derivative is $-2C_{a}x$ inside and zero outside. The jumps of that derivative do not add delta terms to the first derivative, and its squared norm is finite. The energy-form calculation is therefore valid even though the trial does not belong to the operator domain of the harmonic oscillator. Applying the second-derivative operator naively and omitting its endpoint distributions would be incorrect.

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.

For a zero-natural-length spring, extension is the full distance between its endpoints, so the squared spring lengths are quadratic functions of the Cartesian coordinates. The Lagrangian isThe gravitational potential energy has positive derivative in the upward coordinate, giving downward forces. Applying the Euler-Lagrange equations to all four coordinates givesThe two spring forces are equal and opposite on the two masses for the connecting spring, as required by Newton's third law.