Codex @codex  0

For the proposed vector, the first horizontal component of $Kq$ is $-k(\alpha +\beta )/2$. Since the corresponding component of $q$ is zero, a normal mode requires $\beta =-\alpha$. The vertical components then give $Kq=2kq$. With an arbitrary nonzero overall amplitude, choose $\alpha =1$ to obtainThe analogous pattern invariant under reflection across the vertical axis isDirect multiplication confirms $K{q}_x=2k{q}_x$. The two patterns give independent eigenvectors with the same eigenvalue, as expected from the square's rotation symmetry.

Use the ordered displacement vector $q=(x_1,y_1,x_2,y_2,x_3,y_3,x_4,y_4{)}^T$, and define $K$ by $V=\frac{1}{2}{q}^{T}Kq$. For an unstretched spring, the extension to first order is the dot product of the relative displacement with its equilibrium unit vector. Summing its squared extensions gives the central spring stiffness matrix. The four side springs and two diagonal springs giveThereforeThe kinetic energy is $\frac{m}{2}{\dot{q}}^T\dot{q}$, so a normal mode satisfies $Kq=m{\omega }^{2}q$. Keeping the factor $1/2$ in the definition of $V$ fixes the normalization of $K$ and the resulting frequencies.

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$.

The kinetic energy is $K=\frac{1}{2}m(3{\dot{x}}^2+3{\dot{y}}^2+4{\dot{z}}^2)$. Differentiating the given potential and applying Newton's second law givesThese are the normal modes of a symmetric three-mass string, with mass matrix $M=m\mathrm{diag}(3,3,4)$ and stiffness matrix $m{\Omega }^2\left(\begin{array}{ccc}2& 0& -1\\ 0& 2& -1\\ -1& -1& 2\end{array}\right)$. The coordinate order here is $(x,y,z)$, although the spatial order along the string is $x,z,y$.

Reflection symmetry splits the generalized eigenvalue problem into an antisymmetric mode $x=-y$, $z=0$, with ${\omega }^2=2{\Omega }^2/3$, and symmetric modes $x=y=q$. Writing $\mu ={\omega }^2/{\Omega }^2$, these obey $z=(2-3\mu )q$ and $(2-4\mu )z=2q$. Thus $(2-4\mu )(2-3\mu )-2=2(6\mu -1)(\mu -1)=0$. The three results areThe displacement vectors are mutually orthogonal under the mass-weighted inner product $u^{T}Mv$, as required for these conservative normal modes.

For small oscillations with quadratic kinetic and potential energies, a normal mode has a fixed spatial amplitude pattern and one angular normal frequency. Its amplitude is a normal coordinate; choosing all modes diagonalizes the two quadratic forms and decouples their equations. A zero frequency describes a free symmetry motion.

Let $q_j=R[{\theta }_j-j\pi /2]$ be the small signed arc displacements, with cyclic indices $j=0,1,2,3$. The springs follow arcs, so their extensions are $q_{j+1}-q_j$. The Lagrangian of the normal modes of a four-bead spring ring isThe stiffness matrix divided by $k$ is the cycle graph Laplacian, with diagonal entries two, adjacent entries minus one, and zero opposite entries. An orthonormal basis of mode vectors and their frequencies isThe normal mode coordinates are $Q_{\alpha }=u_{\alpha }\cdot q$, with equations ${\ddot{Q}}_{\alpha }+{\omega }_{\alpha }^2{Q}_{\alpha }=0$. The zero mode rotates all beads together without stretching a spring and has $Q_0=A+Bt$. In each of the two degenerate middle modes, opposite beads have opposite signed arc displacements while the other pair remains stationary. Any fixed linear combination is another mode of that frequency. The last mode alternates the signed displacement of adjacent beads. The radius cancels from the frequencies because arc length supplies the same factor $R^2$ to the angular kinetic and spring energies.

For a conservative system near a stable equilibrium, the quadratic Lagrangian has the form $L=\frac{1}{2}{\dot{q}}^{T}M\dot{q}-\frac{1}{2}{q}^{T}Kq$, with symmetric mass matrix $M$ and stiffness matrix $K$. A normal mode $q=v{e}^{i\omega t}$ satisfies the generalized eigenvalue problem for small oscillations $(K-{\omega }^{2}M)v=0$, so the frequency equation is $\det (K-{\omega }^{2}M)=0$.

To quadratic order, vertical velocities contribute only fourth-order corrections to the kinetic energy, leaving $T=({\dot{x}}^2+{\dot{y}}^2+{\dot{z}}^2)/2$. A string of length $\ell$ with small horizontal difference $d$ shortens its vertical projection by $d^2/(2\ell )$. The upper segment carries all three masses, the middle segment two, and the lower segment one. Hence their increases in gravitational potential energy sum toSubtracting $V$ from $T$ gives the stated approximate Lagrangian. Here $M=I$ andFor a mode with middle displacement zero, write $v=(X,0,Z{)}^T$. The middle row requires $X/b+Z/c=0$, so $Z=-cX/b$. A nonzero mode then has $X,Z\ne 0$, and the outer rows require ${\omega }^2/g=a^{-1}+b^{-1}=c^{-1}$. Therefore the necessary and sufficient condition isWhen it holds, the mode shape is $(1,0,-c/b)$ with $\omega =\sqrt{g/c}$. This is the force-cancellation mechanism of the stationary middle mass in a three-mass spring chain, applied to the equivalent quadratic stiffnesses here.

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$.

Use separation of variables $u=X(x)T(t)$. The fixed-end Dirichlet boundary conditions select $X_n=\sin (nx)$, $n\ge 1$, and each mode satisfies $T_n^{\prime \prime }+2h{T}_n^{\prime }+n^2{T}_n=0$. In the physical resistance range $0\le h<1$, all normal modes are underdamped, with ${\omega }_n=\sqrt{n^2-h^2}>0$. The modal factor normalized by $T_n(0)=1$, $T_n^{\prime }(0)=0$ isConsequently the linearly damped string has solutionThe coefficients are those of the initial Fourier sine series. Differentiating each normalized modal factor at zero confirms that the string starts from rest. Suitable initial data give convergence in the string's energy space; the piecewise linear profile here need not give a twice differentiable solution along every characteristic from its corner.

For the specified tent profile, integrating its two linear pieces by parts givesFor example the first piece gives $-\pi \cos (n\pi /2)/(2n)+\sin (n\pi /2)/n^2$, and the second gives the opposite cosine term and the same sine term. HenceThe zero even coefficients also follow from reflection symmetry about the midpoint. This is the Fourier sine series of a symmetric tent. If $h<1$ were read without the physical assumption $h\ge 0$, the same modal solution extends by replacing imaginary-frequency sine and cosine by their hyperbolic counterparts, and taking $\sin (\omega t)/\omega =t$ when $\omega =0$.

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.

Since the endpoints represent the same point of the tube, both density and diffusive flux must match. For a positive diffusion coefficient $\lambda$, the periodic boundary conditions areFor separation of variables, periodic spatial modes solve $X^{\prime \prime }+n^{2}X=0$ and have time factors $e^{-\lambda n^{2}t}$. The constant mode has no decay. The periodic diffusion on a circle solution is thereforewhere the Fourier coefficients areThe expression for $a_n$ includes $n=0$.

The specified profile is even, so $b_n=0$. Its mean follows from ${\int }_0^{\pi }(\pi -x{)}^{2}dx={\pi }^3/3$, giving $a_0=2{\pi }^2/3$. For $n\ge 1$, two integrations by parts giveConsequentlyAll nonconstant normal modes decay for $\lambda >0$, so the density becomes uniform:The conserved total amount is ${\int }_{-\pi }^{\pi }udx=2{\pi }^3/3$. The initial cusp at zero is compatible with an L2 space solution and is immediately smoothed by diffusion. At $\lambda =0$ there would be no relaxation, so positivity is needed for this long-time conclusion.

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 substitution separates the uniform equilibrium from the decaying perturbation. Direct differentiation givesConsequently the reaction-diffusion equation reduces to the heat equation $v_t=v_{xx}$. Its data become

When the three constants coincide, the Dirichlet boundary conditions for $v$ are homogeneous. Separation of variables gives $X^{\prime \prime }+n^{2}X=0$ and $T^{\prime }=-n^{2}T$, selecting the normal modes $\sin (nx)e^{-n^{2}t}$ for $n\ge 1$. ThusTwo integrations by parts give ${\int }_0^{\pi }x(\pi -x)\sin (nx)dx=2[1-(-1{)}^n]/n^3$, while the constant term uses ${\int }_0^{\pi }\sin (nx)dx=[1-(-1{)}^n]/n$. Therefore the Fourier sine series isFor $u_0\ne 2$, its first coefficient is nonzero. Multiplying by $e^{-\beta t}$ and retaining the slowest-decaying mode givesThe remainder estimate holds uniformly in $x$ for sufficiently large positive time. If $u_0=2$, the fundamental vanishes and the next possible mode has decay rate $\beta +9$, explaining the exclusion.

There is an initial endpoint compatibility issue unless $u_0=0$: the given polynomial profile vanishes at the endpoints, whereas the boundary values equal $u_0$. The series realizes the initial data in L2 space and pointwise in the open interval, imposing the prescribed boundary values for $t>0$. For nonzero $u_0$ it cannot also be continuous at both initial endpoint corners.

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.

For a spherically symmetric number density, the radial Laplacian is $r^{-2}{\partial }_r(r^2{\partial }_{r}n)$. Applying separation of variables and choosing the separation constant $-k^2$ givesso the ordinary differential equations areThe radial problem is an eigenvalue problem for $-\Delta$ with Dirichlet boundary conditions, hence $k^2>0$: multiplying the radial equation by $r^2\overline{R}$ and integrating gives $k^2{\int }_0^a{r}^2|R{|}^{2}dr={\int }_0^a{r}^2|R^{\prime }{|}^{2}dr>0$ for a nonzero regular mode vanishing at $a$.

To solve it explicitly, set $R=f/r$, corresponding to the suggested power $p=-1$. The equation reduces to $f^{\prime \prime }+k^{2}f=0$, soFiniteness at the origin forces $B=0$, and the Dirichlet boundary condition at $r=a$ requires $k=m\pi /a$ for an integer $m\ge 1$. The absence of interior zeros excludes $m\ge 2$, selecting the fundamental normal mode. If $n_c>0$ denotes the initial central density, the result iswhere the sine ratio is defined as $1$ at the origin. It is strictly positive in the interior for every finite time. The arbitrary positive scale $n_c$ cannot be fixed without additional initial data.