Distinct squared traction-free eigenfrequencies in a bounded lossless elastic body have eigenfunctions orthogonal in , by the Betti identity for elastodynamic fields. Normalize one field from each eigenspace. This is an orthonormal set in a separable Hilbert space, so it is countable: disjoint balls of radius less than around its unit vectors must contain distinct points of a countable dense set. Positive and negative frequency partners add at most a factor of two. This argument proves countability without requiring a separate compact-resolvent proof.
Damping 2026-10-14
Damping dissipates the energy of an oscillation, reducing its amplitude. A linearly damped string applies this mechanism to each normal mode.
Why can stronger damping make a return to equilibrium slower? A weak brake lets a spring exchange potential energy with kinetic energy many times; a very strong brake makes the final motion creep. Critical damping marks the boundary between these behaviours. The distinction depends on the force law: the classification below concerns a single damped harmonic oscillator with viscous damping, rather than every mechanism that removes energy.
Let be displacement from equilibrium, the moving mass, the spring constant, and the constant coefficient of viscous damping. With no driving force, Newton's second law gives
The two initial conditions specify displacement and velocity. Define the undamped eigenfrequency as an angular frequency and the damping ratio as . The ordinary differential equation becomes . A constant weight in a vertical spring merely shifts equilibrium; it is already absorbed in the definition of .[ref]
For example, the ideal parameters and give and critical damping at . Values give damping ratios , respectively. These are model choices, not measured values. For release from rest with , set and ; then , where primes now denote derivatives with respect to .
Figure 1.
Ideal viscous oscillator responses: release from rest in three damping regimes, and the effect of initial velocity at critical damping
. Details.
The left panel uses identical initial conditions , . Underdamping produces repeated crossings of equilibrium; critical damping and overdamping give positive decreasing responses. The right panel keeps critical damping but changes initial velocity: gives one crossing, although there is no oscillation. Both panels are exact model predictions.
Distinct linearly independent normal modes are degenerate if they have the same eigenfrequency. For , their eigenvectors span a common eigenspace; linear combinations in that space are also modes with the same frequency. A symmetry preserving both quadratic energies maps this space into itself. Higher-dimensional irreducible symmetry representations can enforce degeneracy, while perturbations breaking that symmetry can split it. Degeneracies can also occur accidentally without such symmetry protection.
The University of California, Santa Barbara's cord demonstration uses a function generator and mechanical driver loosely attached near a clamped cord end. The other end passes over a pulley to a mass hanger, setting tension. Its photograph shows operation with a strobe; the page reports measured and calculated lowest-mode frequencies. Varying drive frequency at fixed loading selects higher normal modes. This clamp arrangement differs from a string tied directly to a moving driver.[ref]
A controlled frequency scan holds the cord, endpoint separation and hanging load fixed, then varies frequency slowly to locate large steady amplitudes. Read stationary nodes of a standing wave and count loops; compare the inferred with . As additional checks, compare small drive settings, repeat the scan in reverse, and change only the hanging load to test the predicted square-root tension dependence. These checks test the ideal model rather than adding measurements to the demonstration.
An elastic cord's stretched linear mass density can differ from its relaxed value; the published calculation uses the relaxed value. A strobe can freeze motion at misleading settings, so identify nodes of a standing wave across several snapshots. A moving attachment need not be an exact node of a standing wave. These limits qualify any comparison with ideal fixed-end eigenfrequencies.[ref]
For the ideal uniform wave equation on a string, the wave speed is , where is constant tension and is linear mass density. Assume small transverse slopes, negligible bending stiffness and no damping. Fix both ends at zero displacement, giving the homogeneous Dirichlet boundary conditions . Our standing wave already satisfies the first. For a nonzero amplitude, the second requires , hence
The value gives zero displacement in this sine family, not a vibrating normal mode. Each positive fits half-wavelengths into the string: it has antinodes, interior nodes of a standing wave, and two endpoint nodes of a standing wave. For and , ; the third normal mode has and wavelength .
Figure 1.
Third-mode nodes and antinodes, with the discrete frequency spectrum of a one-metre ideal string at wave speed twenty metres per second
. Details.
The equally spaced eigenfrequencies follow from this string's uniform wave speed and fixed-end boundary conditions. They are not a universal spectrum for every vibrating object. Changing an endpoint constraint changes the allowed mode shapes; changing tension changes wave speed. At fixed and , quadrupling doubles every . The more general construction of arbitrary initial motion belongs to normal mode and wave equation on a string.
For an isothermal gas slab, a horizontal Fourier mode has , , , and . Multiply momentum by and integrate by parts, then use the Poisson equation for Newtonian gravity. Vanishing pressure-work and gravitational boundary terms give the displayed real quadratic identity, with . Thus every admissible squared eigenfrequency is real; a negative value gives gravitational growth rather than an oscillatory instability.
Normal mode Created 2026-09-24 Updated 2026-10-14
A normal mode, also called an eigenmode, is a pattern of free oscillation in a linear system whose relative displacements stay fixed while a single temporal factor changes their size and sign. For an undamped stable vibrational mode,
Here is the mode shape, sets the amplitude, is the eigenfrequency expressed as an angular frequency, and is the oscillation phase. Every snapshot is a scalar multiple of the same mode shape; neighbouring regions may move in opposite directions. Normal modes of a double pendulum show that normal modes also describe coupled discrete systems, beyond strings. Their shapes are eigenvectors, rather than spatial eigenfunctions. [ref]
For the ideal wave equation on a string, take a uniform string of length with constant tension and linear mass density . Assume small transverse slope, negligible bending stiffness, no damping, and fixed endpoints. Its wave speed is , and the nonzero normal modes are
The angular frequency is in radians per second; the frequency is in hertz. The wavelength is . The mode shape has loops separated by nodes of a standing wave, including both fixed endpoints. [ref]
As a worked example, choose , , and . Then and , , . Their oscillation periods are , , and . These are deliberately simple model parameters, rather than measurements of a musical instrument. The plot divides displacement by so that the three mode shapes can be compared, and uses each mode's own oscillation period for its snapshots.
Figure 1.
First three ideal fixed-string normal modes at fractions of their own periods; black dots mark permanent nodes
. Details.
For , at every point has zero displacement, but most points have nonzero velocity. Only zeros of the mode shape are permanent nodes of a standing wave. At the original mode shape has reversed sign, and at it returns. The synchronized conversion between kinetic energy and potential energy is shown in normal-mode energy exchange.
The short derivation uses separation of variables: substituting into the wave equation gives and . The fixed-end boundary conditions remove the cosine spatial term and require , hence and . Different geometries and boundary conditions give different eigenfunctions; a sine mode shape is a property of this uniform fixed string, not a universal form of an eigenmode.
The usefulness is that a complicated vibration becomes a collection of harmonic oscillators. For sufficiently regular fixed-end initial data, the Fourier sine series gives
The initial displacement and velocity determine and . This is the superposition principle in a Fourier sine basis, and modal energy distribution measures how strongly each normal mode is excited. The striking simplification is that each component keeps its identity while their sum can continually change shape. [ref]
These independent normal modes belong to the linear model. At larger slopes, changing tension introduces nonlinear effects and can couple the modes. A linearly damped string still admits a Fourier sine basis, but its temporal factors decay and need not oscillate; its energy dissipation identity for a linearly damped string replaces conservation of energy. More general damping need not leave undamped mode shapes independent.
For three equal simple pendulums in a straight open chain with springs only between neighbours, the spring stiffness matrix in displacement coordinates is . Its eigenvectors , and have respective dimensionless eigenvalues . Adding the gravitational restoring term yields the displayed squared eigenfrequencies. These are free spring ends: fixed wall springs would change the matrix and mode shapes.
Choose field to be at frequency , and field to be at . Zero body forces and zero boundary tractions make all source and surface terms vanish. Therefore
For distinct squared frequencies this proves weighted orthogonality,
Assume positive mass density bounded above and below and finite-energy fields in the ordinary bounded domain. The weighted space is a separable Hilbert space. Select one unit eigenfunction for each distinct squared frequency. They form an orthonormal set. Let be a countable dense subset of ; assign to each eigenfunction a point within distance . Two orthogonal unit eigenfunctions have distance , so the same cannot be assigned to both. This gives an injection into a countable set.
Hence there are at most countably many squared eigenfrequencies, and at most twice that many signed frequencies. Multiplicity at one frequency does not change the argument. A zero frequency associated with rigid translations or rotations is included. This proves the countability of elastic eigenfrequencies directly; proving discreteness or finite multiplicity would require additional spectral regularity arguments.
If two harmonic solutions have the same body force and prescribed boundary traction, their difference has no body force and is traction-free. A nonzero difference would therefore be a free eigenfunction at that frequency. For any frequency outside the countable exceptional set found above, the difference must vanish, so the forced boundary-value solution, when it exists, is unique. At an eigenfrequency, a free eigenfunction can be added to any existing solution, and a resonant forcing may fail to admit a steady harmonic response. This distinguishes uniqueness from existence.
For the full time-dependent problem, let be the difference between two fields with identical forces, tractions and initial displacement field and velocity. It obeys , zero boundary traction, and zero initial data. Its energy uniqueness for traction-driven elasticity follows by direct integration:
Multiply the equation by and integrate by parts. Major symmetry turns the strain term into the derivative of elastic energy, giving . Stable elasticity makes , and gives . Thus equals its initial value zero. Even rigid zero-strain motions cannot survive these initial conditions.
The PDF correctly specifies zero initial displacement field and zero initial velocity; the TeX duplicates displacement field in the second condition. The complete causal displacement field is unique, including when its Fourier representation contains resonant modal contributions fixed by the initial data.
Eigenfrequencies and mode geometry. Write and . The characteristic polynomial is . Therefore the two distinct real eigenvalues are
The sign assumptions imply and ; they do not alone imply . A positive lower frequency additionally requires , as in the usual nondegenerate physical Laplace-Lagrange secular matrix.
Distinct eigenvalues give an eigenbasis, and each secular eigenmode evolves by multiplication by . Hence
with nonnegative amplitudes and initial phases. The second row of each eigenvalue equation gives . For the upper eigenvalue this ratio is negative, whereas for the lower one it is positive. Thus
The upper secular eigenmode has anti-aligned apsides, and the lower has aligned apsides. If one amplitude vanishes, its phase is immaterial.
The initial vector sums. The planet-1 modal vectors point at and . Their resultant at requires . The planet-2 vectors point at and ; their resultant at requires . These follow by resolving components perpendicular to each resultant. In particular the specified phases constrain the matrix: . They cannot be imposed on an arbitrary matrix satisfying only the sign assumptions.
Put and let . For a matrix compatible with those phases,
At zero the resultants have lengths and and precisely the specified longitudes of periapsis. The diagram uses , as a drawing normalization, not as an additional assumption on the physical planets.
Figure 1.
Initial secular eigenmode vectors and their resultants at apsidal longitudes 60 and 150 degrees
. Details.
Each modal vector rotates counterclockwise if its eigenvalue is positive; their sum executes a beat at frequency . For these phases,
The orbital eccentricities exchange their maxima and minima, with respective ranges to and to . The longitude of periapsis is the argument of each displayed vector sum; it does not in general advance uniformly at either individual eigenfrequency. The complex eccentricity trajectories are generally quasiperiodic; they close only for commensurate modal frequencies. This example has a dominating upper mode in planet 1 and a dominating lower mode in planet 2, so their apsidal phases have different winding rates.
Amplitude product and the angular-momentum convention. Since ,
Under the ratio supplied in the paper, this is approximately , to the same accuracy as that supplied ratio. There is a physical convention issue: with the usual planet labels and , conservation of quadratic angular momentum deficit requires , giving . This follows directly by differentiating . Thus the paper's stated angular-momentum ratio is reversed for the conventional column-vector equation . The exact algebraic product above is independent of that naming discrepancy.
Tidal dissipation and damping. Assume the tidal dissipation contributes only the specified linear eccentricity damping, without changing the conservative matrix or adding apsidal precession. Then
The plus sign in is crucial: multiplication by produces the negative real damping rate. Its exact eigenvalues are
The square-root branch is chosen to approach as . With , their first-order imaginary parts are
Both are positive because , so both secular eigenmodes decay as , although only planet 1 feels direct eccentricity damping. The coupling transmits that damping to planet 2. Their precession frequencies change only at second order; their eigenvectors acquire small complex corrections, so exact alignment or anti-alignment becomes a small apsidal phase lag. The more slowly damped mode eventually dominates, unless its initial coefficient vanishes: it is the lower aligned mode for and the upper anti-aligned mode for . At both leading damping rates equal .
The stated slow-secular condition presumes a nonzero positive . It need not imply near a nearly degenerate pair, so the exact square-root formula is the appropriate answer there. Exceptionally, the discriminant vanishes at , ; then a generalized eigenvector supplies a term. There is still decay. Indeed positive diagonal weights satisfying give
Nonzero off-diagonal coupling excludes a nondecaying mode supported solely on planet 2. This checks the damping signs without relying on the weak-damping expansion.
Translate the generalized coordinates so that the equilibrium is at zero. Its potential gradient vanishes, and its quadratic potential energy is , where is the symmetric Hessian matrix. The Euler-Lagrange equations for small motion are
A normal mode has fixed component ratios, , with
The angular frequency is its eigenfrequency. A zero mode instead admits a constant or linear time factor. For modes , symmetry gives
Hence mass-matrix orthogonality of normal modes is
Within a degenerate eigenspace, a basis can be chosen orthogonal in this positive weighted inner product.
For three masses connected pairwise, the exact energies are
A fixed common translation preserves all differences and velocities. A fixed rigid rotation 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 has zero spring extensions, hence zero force and minimum potential, proving that it is an equilibrium.
For the specified displacement coordinates, write
The exact potential is . For an equilibrium edge vector and small perturbation ,
Thus only each first-order extension is needed for the quadratic potential. These extensions, with an immaterial sign choice for the horizontal edge, are
Consequently , with remainder . If , then . Expanding the squares gives
The small-oscillation mass and stiffness matrices are and the displayed . Hence
The two common translations give zero-mode vectors
An infinitesimal counterclockwise rotation sends each equilibrium vertex to , where . Therefore a third zero-mode vector is
Each vector makes all three vanish, so lies in the kernel of a linear map and therefore in the kernel of . The rows of are independent: the and entries first force the coefficients of its last two rows to vanish, then the first coefficient vanishes. Thus , and the rank-nullity theorem proves that these are exactly the three zero modes.
For radial motion take
so that to first order. The force on mass is , since . Newton's equation therefore gives . Equivalently . Hence the breathing normal mode has
All 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 ; together with the breathing mode and three rigid zero modes they account for all six degrees of freedom.
Figure 1.
Two translations, an infinitesimal rotation and the radial breathing mode of the equilateral spring triangle
. Details.
The Lagrangian is , and the Euler-Lagrange equations are
At a rest equilibrium , for all three coordinates. The Taylor expansion therefore starts as
The symmetric Hessian matrix is the matrix 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 . Terms involving the variation of are cubic in small displacements and velocities, so the quadratic small oscillation Lagrangian is
A normal mode has fixed displacement ratios, . Its angular frequency is an eigenfrequency, or normal frequency, and satisfies the generalized eigenvalue problem for small oscillations
Here is positive definite for independent physical coordinates and positive masses. A stable quadratic minimum makes ; a negative value would describe exponential instability rather than oscillation.
For the suspended bar, let and denote the actual incremental spring constants of the upper and each lower spring. The physical displacement of the bar is , while the lower spring extensions change by and . Consequently
The linear gravitational and elastic terms cancel at the loaded equilibrium. Expanding the remaining squares gives the small-oscillation mass and stiffness matrices
Thus
Source convention: the printed matrix is obtained by setting and . Although the PDF prose calls 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 should instead be omitted. The numerical case has , so this discrepancy does not affect its results.
With the numerical parameters, and
Use the exchange symmetry of the lower masses to find the vertical modes of a horizontally constrained suspended bar. The antisymmetric vector obeys , giving . In the symmetric sector put and . The two independent equations are
Their determinant gives , so . Taking , the full set is
In the lowest normal mode, all three masses move in phase, with the lower masses moving about 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 of its displacement. Any small motion is a superposition of these three modes.
Figure 1.
Relative vertical displacements in the three normal modes, with the bar constrained to remain horizontal
. Details.
Standing wave Created 2026-09-24 Updated 2026-10-18
A standing wave results from the superposition principle applied to equal-frequency, equal-amplitude waves travelling in opposite directions. Its spatial nodes of a standing wave remain fixed while its displacement oscillates in time.
For a pure standing wave, the two travelling waves also have equal amplitudes: their interference leaves a fixed mode shape, rather than a profile moving along the string. Choose the position origin and time origin so that their transverse displacements are
Here is each constituent's amplitude, is the wavenumber, is the wavelength, and is the angular frequency corresponding to frequency . Their wave speed is . The superposition principle and the sine addition formula give
The spatial factor stays put; the temporal factor changes its size and sign. Other constant wave phases shift the position or time origin. Unequal constituent amplitudes generally leave moving displacement at the would-be nodes of a standing wave, rather than this exact factorization.[ref]
For an original ideal-string example, set , , and choose the third normal mode. Then , , and by fixed-end standing wave frequencies. With each constituent amplitude ,
Its nodes of a standing wave are at and its antinodes at . These are constructed model parameters, not experimental readings.
Figure 1.
Opposite travelling waves add to a standing third mode; permanent nodes survive every phase snapshot
. Details.
At , every point has zero displacement, but most have nonzero transverse velocity: the string is crossing its straight configuration. The black dots alone mark permanent nodes of a standing wave. At an antinode the greatest transverse speed is , much smaller than the chosen wave speed. One is motion of string material; the other is propagation of the constituent profiles.
The beauty of a standing wave is that motion reveals geometry: fixed nodes of a standing wave turn a moving disturbance into a spatial ruler. Their separation measures wavelength, and boundary conditions turn the continuous choice of wavenumber into a discrete set of normal modes. This is useful both for interpreting string vibrations and for understanding how a bounded wave equation selects its eigenfrequencies.
Toroidal normal mode of the Earth Created 2026-10-20 Updated 2026-10-25
A toroidal normal mode of the Earth is a free oscillation with tangential, divergence-free displacement. It is the spherical counterpart of the transverse shear motion in a Love wave. A bounded spherical Earth has discrete angular orders and eigenfrequencies, so the arbitrarily low-frequency continuum of a flat infinite layer is not a literal statement about global nonrigid oscillations.