With gravity omitted, the initial data contain only the first fixed-end normal mode, so write
Then
Because , this oscillator is critically damped, and
Therefore
The function is strictly convex. Subject to , the Jensen inequality gives
with equality only at
This is the minimum of the relative potential, so it is a stable relative equilibrium.
Let near an equally spaced configuration, and put
The gap perturbations are
whose sum vanishes. A Taylor expansion gives the quadratic energies
Thus the linearized equations are
The matrix in parentheses is the Graph Laplacian of . Its constant eigenvector has eigenvalue zero, while every vector whose components sum to zero has eigenvalue three. One orthonormal set of three normal modes and their angular frequencies is therefore
The zero mode is rigid rotation. The positive and equal squared frequencies on the two-dimensional relative subspace prove stability and agree with the general normal modes of three equal masses with a symmetric gap potential.
The linearized equations are homogeneous and have constant coefficients in and . Consequently, spatial Fourier modes evolve independently, and the superposition principle reconstructs arbitrary sufficiently regular initial disturbances from them. It is therefore enough to seek the normal modes
Here the real is the wavenumber, is the possibly complex temporal exponent, and taking the real part produces a physical displacement.
The steady radial equation first gives the base-pressure gradient
For the stated normal mode, put
Retaining terms linear in the disturbance in the Euler equations for an inviscid fluid gives
together with incompressible flow
The and terms are the linearized centrifugal and angular-momentum couplings of the swirling base flow.
Linearization about removes the quadratic term and gives
The three homogeneous boundary conditions select the vertical normal modes
Take a horizontal Laplacian eigenfunction satisfying
On , the full Laplacian has eigenvalue . The growth rate is consequently
Neutrality occurs at
For fixed nonzero , this increases strictly with , so the first unstable vertical mode is and
For the normal mode
the linear material derivative becomes
Substitution into the equation from part b and division by gives the Taylor–Goldstein equation
where
and the Scorer parameter is
When , set . The reduced ordinary differential equation is . Enforcing the same initial conditions gives
The corresponding overdamped normal mode can also be written
Both rates in the exponential functions are negative because , so the hyperbolic sine factor does not prevent decay.
When , the reduced ordinary differential equation is . The initial conditions , give
This is the critically damped normal mode, also obtained by taking in .
Substitute the separated solution into the damped wave equation. Its time derivatives give
whereas . Thus
The last two conditions follow from and . For , let . The simple harmonic motion equation has solution of a differential equation
This is an underdamped normal mode: its oscillations have the decaying factor .
For , let denote the time factor normalized by , . By superposition the Fourier sine series solution is
At , its initial velocity is . Multiply by , integrate on , and use the permitted interchange and Fourier orthogonality. This yields the requested integral expression
Under , is unchanged and . Hence every even-index Fourier coefficient is zero, including .
For odd , use and the product-to-sum formula:
Thus
In particular and . An explicit version of the solution of a differential equation is
Only the first, overdamped normal mode, the third, critically damped normal mode, and the odd underdamped normal modes of index at least five occur. No odd-index Fourier coefficient vanishes.
For , substitution of the separated solutions into the heat equations gives
The zero exterior boundary values are incorporated by
To obtain a eigenvalue equation for heat conduction in two joined slabs, write and . Continuity of temperature and heat flux at zero gives the homogeneous linear system
A nonzero amplitude pair exists precisely when its determinant vanishes, equivalently
If , division by gives the printed transcendental equation
The printed necessity claim needs the qualification that both tangents are defined. If just one of vanishes, . If both cosines vanish, however, and legitimate normal modes remain. Their temperature at the thermal interface need not be zero; the temperature condition fixes the amplitude ratio, while each heat flux is zero there.
For a concrete counterexample, take identical slabs: , , , and . Then
with each formula restricted to its own slab, satisfies both heat equations and all four boundary conditions. But , so neither tangent is defined. This directly disproves the unqualified printed assertion. The determinant condition retains these missing normal modes.
Finally, must be checked before dividing by it. Then , , and interface continuity gives , while heat flux continuity gives . Positivity of the thermal conductivities forces . Thus In the determinant formula, its apparent zero at is spurious because the chosen sine basis degenerates there.
Measure as horizontal displacements from the respective downward equilibrium positions; the absolute horizontal coordinates are and . Each vertical displacement is of second order in its small angle. Thus the spring is horizontal to first order, its extension is , and Hooke's law gives horizontal forces and . Combining these with the linearized simple pendulum restoring force yields
For the normal modes of two spring-coupled equal pendulums, add and subtract these ordinary differential equations. The normal mode coordinates and satisfy
Hence the two natural angular frequencies and displacement patterns are
Ordinary frequencies are and . In the first normal mode, both particles swing together with equal amplitudes and the spring is unextended to first order. In the second normal mode, the particles swing in opposite directions with equal amplitudes, so the spring adds a restoring force and increases the angular frequency.
Figure 1.
The two normal modes: equal displacements together and equal displacements in opposite directions
. Details.
Use separation of variables to seek a normal mode , with angular frequency and wavenumber . The endpoint boundary conditions are built into
Let and denote the mass's displacement amplitude by . Continuity and the mass equation become
The two possibilities must be separated before dividing by .
Antisymmetric modes. If , then , , and . The force equation gives , since . Thus and the mass stays at rest. These allowed angular frequencies satisfy
The restoring forces on the mass cancel even though the string vibrates.
Symmetric modes. If , continuity forces . Cancelling the nonzero amplitude in the mass equation and using gives
These modes have and move the mass. In dimensionless form their equation is , where . This is the spectral splitting for a string with a central point mass.
Figure 1.
Antisymmetric zeros and symmetric frequency intersections
. Details.
For the antisymmetric family, the graph of crosses zero at every positive multiple of . For the symmetric family, on each interval , , the graph of decreases continuously from to zero while increases. By the intermediate value theorem and strict monotonicity, there is exactly one intersection in each interval. On the remaining half-intervals , so there are no further positive roots. This proves graphically and analytically that both families contain infinitely many normal modes. The figure uses for illustration.
The algebraic case is not a mode: a zero-frequency profile would be linear on each half, continuous and zero at the endpoints; the zero-acceleration mass equation also makes its slopes agree, forcing the entire profile to be zero. More generally, integration by parts on the two halves and the matching condition give
for a nonzero real mode. Thus no nonzero static or exponentially growing mode has been omitted.
For a short string element of length , the small-slope approximation makes the upward components of tension at its two ends and . Its mass is . Newton's second law gives
Divide by and take the limit to obtain the wave equation on a string,
Here is linear mass density, and the derivation assumes small transverse slopes, uniform tension, and a flexible string.
At the right endpoint the string pulls the ring toward the interior. Its vertical force is therefore to first order; the frictionless wire supplies no vertical force. The ring obeys
This is the dynamic endpoint law for a string with an endpoint mass.
For a separated oscillatory solution, write with and . The fixed left end removes the spatial cosine, so the nonzero normal modes have
The endpoint law gives . Since would also force the cosine to vanish, it can be divided out for a nonzero mode. Using gives
For physical , these give all positive modal frequencies. There is no nonzero static linear mode: the left value and right slope conditions would force it to vanish. Positivity also follows by integration by parts from
For , the right endpoint is free, with the Neumann boundary condition . At , the allowed wave numbers are , , and . Zero initial displacement removes every cosine-in-time coefficient. The initial velocity contains only the modes and , each with coefficient one. Dividing those velocity coefficients by their angular frequencies gives
At each term vanishes, at each spatial derivative vanishes, and at the displacement and its time derivative reproduce the two initial conditions.
At and , the fixed-end boundary conditions require . In the oscillatory family, the first condition removes the cosine term and the second requires , so with . The linear and hyperbolic spatial families satisfy both boundary conditions only when identically zero.
Thus the nonzero normal modes and their superposition give
For initial displacement and velocity compatible with the fixed ends, orthogonality of the sine functions determines
These are all the string solutions in the usual class where the series has the regularity needed to satisfy the wave equation.
The initial displacement is already the first spatial normal mode, so sine orthogonality gives and every other . Zero initial velocity gives every . Consequently
Both fixed-end values vanish for all time, and differentiation shows that the initial velocity is zero. This is the fundamental standing wave, with angular frequency .
General definitions. The Euler-Lagrange equation gives . Write and . A normal mode has fixed displacement pattern and common time factor, ; its normal frequency satisfies
Here is angular frequency. A normal coordinate is the coefficient of one pattern in a complete modal expansion, so these coordinates decouple the equations. Since is positive definite, the change turns the mode equation into the ordinary symmetric eigenvalue problem for .
Three-mass equations. Let be horizontal displacements from the unstretched equilibrium, ordered left to right. The kinetic energy and potential energy are
Differentiating the potential energy gives
Thus and
Frequencies and patterns. Put . Expansion of the generalized eigenvalue problem for small oscillations gives
Solving at each root gives
The 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 is symmetric, two mode equations imply
Distinct eigenvalues therefore give the mass-matrix orthogonality of normal modes:
The Euclidean dot product need not vanish for unequal masses; for example , while . The modal squared norms are .
Figure 1.
The three normal modes of the unequal free mass chain
. Masses from left to right are m, 2m and 3m, joined by springs k and 3k. Arrow lengths show the modal amplitude ratios (1,1,1), (-3,0,1), and (1,-2,1). Details.
Normal coordinates. Expand . Project using the mass-weighted inner product and divide by each modal norm to obtain
For example, the inverse equations are , , . They reduce the energies to
and hence
The first coordinate is the centre of mass and has free linear motion. Multiplying by respectively gives mass-normalized normal coordinates if desired.
For the normal modes of an equilateral spring triangle, the complement of the three rigid motions has dimension three. One convenient mutually orthogonal basis is
These vectors are also orthogonal to , and each has squared length . Their geometry comes from applying symmetric linear deformations to the equilibrium triangle: is uniform expansion, stretches one axis while compressing the other, and is the second independent shear pattern. Their extensions are
Insert each row into the force equations. For example, gives and , with the same proportionality in every component. The other two rows similarly give force times their displacement vector. Thus
and the remaining normal modes and normal frequencies are
Any orthogonal pair spanning the same two-dimensional shear subspace is equally valid, because its eigenvalue is repeated. The six displayed vectors are nonzero and mutually orthogonal, so they supply all six independent normal modes without solving a characteristic polynomial.
Figure 1.
Rotation, breathing and two shear patterns of the three-mass triangle
. Arrows show infinitesimal displacements using the equilibrium coordinates in the paper. Translations supply two further zero modes. The breathing frequency squared is 3k/m; both shear frequencies squared are 3k/(2m). Details.
For explicit normal coordinates, normalize by and by , and write . Since all masses are equal,
Positive-frequency coordinates have ; the three rigid coordinates have the linear form found above.
Order coordinates as . For either uniform displacement, or , every relative displacement and every extension vanishes. These normal modes have zero normal frequency and represent two independent planar translations.
For , substitution gives
Thus too. A small rigid rotation through angle moves equilibrium position by . With the specified positions this displacement is , identifying the third vector as infinitesimal rigid rotation. Hence
A zero-frequency normal coordinate obeys , so . This permits a constant displacement or a uniform velocity at the level of the linearized equations. In particular, the rotational pattern preserves lengths to first order; finite-angle rotation follows a curved path outside a single linear displacement pattern.
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, denotes the angular frequency, so the frequency in cycles per unit time is . A normal coordinate is the amplitude of one such pattern, chosen so that the linear equations decouple. For mass matrix and stiffness matrix , the patterns solve . Zero eigenvalues describe free rigid motions rather than periodic restoring oscillations.