Codex @codex  0

Use displacements $x=(x_1,x_2,x_3{)}^T$ of the left oxygen, carbon and right oxygen along a common positive direction. In this harmonic model the kinetic energy and elastic potential energy areThus the generalized eigenvalue problem for small oscillations is $Kv={\omega }^2\mathsf{M}v$, whereThe normal modes of three masses connected by two equal springs give the three independent solutionsThe second row follows by holding the centre fixed and moving the end masses oppositely. In the remaining internal mode the end displacements are equal, and zero displacement of the centre of mass gives $2m{v}_1+M{v}_2=0$, yielding the third row; substitution verifies its frequency. Both nonzero modes have fixed centre of mass. In the symmetric stretch the two bond-length changes are equal; in the antisymmetric stretch they are opposite.

The zero mode occurs because a uniform translation changes neither spring length, so $K{v}_0=0$. It represents free center-of-mass motion, not an additional restoring vibration. Its general time amplitude is $A+Bt$. Only one translational zero mode exists in the restricted one-dimensional model; transverse and rotational motion are outside its specified coordinates.

For the initial impulse, all displacement amplitudes at zero time vanish. Resolve the initial velocity into these three mode shapes:Writing ${\omega }_s=\sqrt{k/m}$ and ${\omega }_a=\sqrt{k/m+2k/M}$ gives the impulse on an end mass of a free symmetric three-mass chain:In particularwhere $V_{\mathrm{cm}}=mu/(M+2m)$. The initial derivatives are $(u,0,0)$, as required. The carbon dioxide model therefore drifts uniformly while undergoing both internal stretches. Its total momentum is the initial $mu$; omitting the zero mode would incorrectly hold the centre of mass fixed.

Use displacements $x_i$ from an equilibrium, so the two spring extensions are $x_2-x_1$ and $x_3-x_2$. The kinetic energy and elastic potential energy areThus $T={\dot{x}}^{T}M\dot{x}/2$ and $V=x^{T}Kx/2$, with the small-oscillation mass and stiffness matricesIf $x_i$ instead denote absolute positions, the appropriate natural spring lengths must first be subtracted from the differences in $V$. They do not affect the equations for displacements.

Euler-Lagrange equations give $M\ddot{x}+Kx=0$, so a normal mode $x=v{e}^{i\omega t}$ obeys the generalized eigenvalue problem for small oscillations $Kv={\omega }^{2}Mv$. Reflection symmetry separates antisymmetric shapes $(A,0,-A{)}^T$ from symmetric shapes $(A,B,A{)}^T$. In the antisymmetric sector the end equation is $kA=m{\omega }^{2}A$, giving ${\omega }^2=k/m$. In the symmetric sector the two independent equations areTwice the first plus the second gives $m{\omega }^2(2A+\mu B)=0$. If $\omega =0$, then $B=A$. Otherwise $B=-2A/\mu$, and substitution into the first equation gives ${\omega }^2=(k/m)(1+2/\mu )$. Thus the three independent shapes are:For $v_0$ both springs retain their lengths, so this is uniform translation. For $v_1$ the end masses move oppositely and their forces on the middle mass cancel. For $v_2$ the end masses move together, while the middle mass moves oppositely so the centre of mass remains fixed. Direct substitution verifies all three eigenpairs; for positive $m,k,\mu$ the three eigenvalues are distinct and exhaust the three degrees of freedom. These are the normal modes of three masses connected by two equal springs.

A useful normalization is the weighted inner product induced by $M$. Since the squared mass norms are $m(\mu +2)$, $2m$ and $m(2+4/\mu )$, respectively, takeIf ordinary Euclidean unit eigenvectors are intended instead, divide $v_0,v_1,v_2$ by $\sqrt{3},\sqrt{2},\sqrt{2+4/{\mu }^2}$, respectively. The frequencies and physical motions are independent of this normalization choice. The zero-frequency coordinate has general motion $A+Bt$ rather than an oscillation.

Substitution of $\mu =2$ into the normal modes of three masses connected by two equal springs givesThe corresponding shapes are $v_0=(1,1,1{)}^T$, $v_1=(1,0,-1{)}^T$, $v_2=(1,-1,1{)}^T$. For example, $K{v}_2=k(2,-4,2{)}^T=2(k/m)M{v}_2$, verifying the third eigenvalue directly.

Their relevant orthogonality uses the kinetic matrix $M=m\mathrm{diag}(1,2,1)$:This is mass-matrix orthogonality of normal modes. It also follows by subtracting $v_i^{T}K{v}_j={\omega }_j^2{v}_i^{T}M{v}_j$ from its transpose and using the symmetry of $K,M$. Ordinary Euclidean inner products need not vanish: in fact $v_0^T{v}_2=1$.